*** empty log message ***

fslmeants.cc

@@ -137,6 +137,7 @@ int main(int argc,char *argv[])
     nvox = MISCMATHS::round(mask.sum()+0.5);
     nt += 3;
   }
+
   Matrix meants(nt,nvox);
   meants = 0;
   long int num=0;
@@ -159,10 +160,16 @@ int main(int argc,char *argv[])
     }
   }
 
+  
   if (verbose.value()) {
     cout << "Number of voxels used = " << num << endl;
   }
 
+  if(nvox!=num)
+    {
+      meants=meants.Columns(1,num);
+    }
+
   // normalise for number of valid entries if averaging
   if (!showall.value() && !showallasavw.value()) {
     if (num>0) meants /= (float) num;