From 44961d6f903ead67f9cd717adf43b8ea0525fdeb Mon Sep 17 00:00:00 2001
From: Mark Jenkinson <mark@fmrib.ox.ac.uk>
Date: Wed, 1 Aug 2001 10:58:59 +0000
Subject: [PATCH] Fixed problem with no mask
---
fslmeants.cc | 10 ++++++----
1 file changed, 6 insertions(+), 4 deletions(-)
diff --git a/fslmeants.cc b/fslmeants.cc
index fd3d211..8792129 100644
--- a/fslmeants.cc
+++ b/fslmeants.cc
@@ -42,11 +42,12 @@ int main(int argc,char *argv[])
volume4D<float> vin;
read_volume4D(vin,argv[1]);
- bool use_mask = false;
volume<float> mask;
if (argc>=4) {
- use_mask = true;
read_volume(mask,argv[3]);
+ } else {
+ mask = vin[0];
+ mask = 1.0;
}
if (!samesize(vin[0],mask)) {
@@ -63,8 +64,7 @@ int main(int argc,char *argv[])
for (int z=mask.minz(); z<mask.maxz(); z++) {
for (int y=mask.miny(); y<mask.maxy(); y++) {
for (int x=mask.minx(); x<mask.maxx(); x++) {
- if ( (!use_mask) || (fabs(mask(x,y,z))>1e-8)) {
- // arbitrary "0" threshold
+ if (fabs(mask(x,y,z))>1e-8) { // arbitrary "0" threshold
meants += vin.voxelts(x,y,z);
num++;
}
@@ -72,6 +72,8 @@ int main(int argc,char *argv[])
}
}
+ cout << "Number of voxels used = " << num << endl;
+
// normalise for number of valid entries
meants /= (float) num;
--
GitLab