The alignment of the T1 to the template looks great but when I warp the set of ROIs, some are missing. I start with 66 points and end up with around 57. I need to retain all 66 ROIs so I was wondering if anyone has an idea about why this is happening and if there is a way to prevent it.
My guess is that some are being squashed into oblivion by the regridding in some cases. That is, by applying the NN interpolation when applying the warp, you are preserving the integer values in your data set. However, it looks like your Undump “ROIs” are just individual voxels: in the case that the warping process leads to a contraction of the source grid, your individual voxel might get “outweighed” by surrounding zero voxels when regridding.
One thing you could do is not apply the “-interp NN” in the 3dNwarpApply command-- that would get rid of your integer-values in the output dset, but you will probably still see some nonzero value for each ROI. You could then use 3dExtrema to make a binary map of all the local maxes-- basically, finding the peak location of your non-integer, mapped ROIs. If your ROIs aren’t too close together (so that 3dNwarpApply doesn’t warp them into each other), that could should work? This way of proceeding assumes that each ROI was just a ‘1’ value, because you would lose the identity of different integers.
Also consider warping the coordinates around which you put the original spheres. Use 3dNwarpXYZ to warp those coordinates. Then create the spheres with 3dUndump at those locations in native space instead. The 3dUndump typically needs a -srad option to create spheres.
Both are perfectly acceptable answers and work fine. For my purposes, it does seem 3dNwarpXYZ is more straightforward.
As always, your help is much appreciated!
Brady
The
National Institute of Mental Health (NIMH) is part of the National Institutes of
Health (NIH), a component of the U.S. Department of Health and Human
Services.