From 09cf5f3bd245309c9c816bbab5ffca9a8613c1d2 Mon Sep 17 00:00:00 2001
From: Mark Jenkinson <mark@fmrib.ox.ac.uk>
Date: Wed, 11 Apr 2012 09:39:46 +0000
Subject: [PATCH] Fixed mm scalings

---
 calc_grad_perc_dev.cc | 12 +++++++++---
 1 file changed, 9 insertions(+), 3 deletions(-)

diff --git a/calc_grad_perc_dev.cc b/calc_grad_perc_dev.cc
index f0d00ef..92f80eb 100644
--- a/calc_grad_perc_dev.cc
+++ b/calc_grad_perc_dev.cc
@@ -57,11 +57,17 @@ int do_work(int argc, char* argv[])
   volume4D<float> fw, xpd, ypd, zpd;
   read_volume4D(fw,fullwarp.value());
   gradient(fw[0],xpd);
-  xpd/=fw.xdim();
+  xpd[0]/=fw.xdim();
+  xpd[1]/=fw.ydim();
+  xpd[2]/=fw.zdim();
   gradient(fw[1],ypd);
-  ypd/=fw.ydim();
+  ypd[0]/=fw.xdim();
+  ypd[1]/=fw.ydim();
+  ypd[2]/=fw.zdim();
   gradient(fw[2],zpd);
-  zpd/=fw.zdim();
+  zpd[0]/=fw.xdim();
+  zpd[1]/=fw.ydim();
+  zpd[2]/=fw.zdim();
   string bname=fslbasename(outname.value());
   save_volume4D(xpd,bname+"_x");
   save_volume4D(ypd,bname+"_y");
-- 
GitLab