*** empty log message ***

cluster.cc

@@ -72,6 +72,9 @@ Option<string> inputname(string("-i,--in,-z,--zstat"), string(""),
 Option<string> copename(string("-c,--cope"), string(""),
 			string("filename of input cope volume"),
 			false, requires_argument);
+Option<string> outpvals(string("--opvals"), string(""),
+			string("filename for image output of log pvals"),
+			false, requires_argument);
 Option<string> outindex(string("-o,--oindex"), string(""),
 			string("filename for output of cluster index (in size order)"),
 			false, requires_argument);
@@ -548,6 +551,11 @@ int fmrib_main(int argc, char *argv[])
     relabel_image(labelim,relabeledim,meanvals);
     save_volume(relabeledim,outmean.value());
   }
+  if (!outpvals.unset()) {
+    volume<float> relabeledim;
+    relabel_image(labelim,relabeledim,logpvals);
+    save_volume(relabeledim,outpvals.value());
+  }
   if (!outthresh.unset()) {
     // Threshold the input volume st it is 0 for all non-clusters
     //   and maintains the same values otherwise
@@ -578,6 +586,7 @@ int main(int argc,char *argv[])
     options.add(outsize);
     options.add(outmax);
     options.add(outmean);
+    options.add(outpvals);
     options.add(pthresh);
     options.add(copename);
     options.add(voxvol);