afni.nimh.nih.gov

You could use p2dsetstat to get the desired threshold for volume AAA in the original DSET_WITH_INFO, for a given pvalue PVAL, say for a bisided stat:

p2dsetstat \
  -bisided  \
   -inset DSET_WITH_INFO"[AAA]" \
   -pval PVAL

... and then you could use that threshold value directly within your warped dset.

--pt