This command would take a mask of an ROI_MASK and a correlation value dataset CORR_DSET, and make a new dset NEW that is the intersection of
+ where the correlation values in CORR_DSET lie within [0.5, 0.6], and
+ where the ROI is:
3dcalc \
-a ROI_MASK \
-b CORR_DSET \
-expr "step(a)*within(b,0.5,0.6)" \
-prefix NEW
If your ROI is located in a map of lots of ROIs, you could select it by its numerical index IDX or attached label LAB (if it has one) with:
3dcalc \
-a ROI_MAP"<IDX>" \
-b CORR_DSET \
-expr "step(a)*within(b,0.5,0.6)" \
-prefix NEW
or
3dcalc \
-a ROI_MAP"<LAB>" \
-b CORR_DSET \
-expr "step(a)*within(b,0.5,0.6)" \
-prefix NEW
You can make multiple new dsets for different bands.
You can overlay the NEW dset on your other dsets to make sure that has gone right.
--pt