¥ The result:
HThresholds the dataset on subrick #1 (-1tindex 1) at value 0.4 (-1thresh 0.4)
HClusters together all the surviving nonzero voxels using
a contiguity test of 1mm and keeping only
clusters at least 500 mm3 in volume (-1custorder 1 500)
HAll voxels in the largest cluster are assigned value 1
(say, 1200 voxels in one cluster), the second
largest are assigned value 2 (say, 920 voxels in one cluster), etc., and the result is written to disk in
dataset ROI_func+tlrc
¥You can use this dataset as a
mask, edit it with the drawing plugin, etcÉ