Hi Cesar,
Daniel's good suggestion aside, what decimals are
you seeing? Are they just things like 2.0000? If
so, that just suggests the input is (stored as)
float.
Also, there is no need for more than 1 3dcalc
command here. A zero volume is not necessary.
Also, there does not seem to be an index in the
name of the output. Consider:
3dcalc -LPI -a MNI152_T1_2009c+tlrc -prefix ROI_parcels_MNI.$iROI \
-expr "a+${iROI}*(step($sphereROI-(x-${xcor})*(x-${xcor})-(y-${ycor})*(y-${ycor})-(z-${zcor})*(z-${zcor})))"
Of course, Daniel's suggestion is a nice way to go.
- rick
