Just estimate the transformation without writing, then write without estimating using the calculated transform to warp an image with the coordinates you want normalized set to 1, and all others 0. Then just look for all non-zero voxels after the transformation.

Cheers,
Michael

On Sat, Jun 4, 2011 at 9:49 PM, Yonatan Dinai <[log in to unmask]> wrote:
Hello all,

I've used the "normalize" function of SPM, and was able to warp two
volumes succesfully.

I would now like to use the output parameters of the normalization
process - the 4x4 matrix of the affine transformation & the 4D matrix
of the DCT coefficients to transform an input point (X,Y,Z in the
source volume) into it's mapped point in the reference coordinates?

Does anybody has a code the computes that transformation? Is there
such a function on SPM? If not, I would be happy to hear of any fast
way to compute such a transformation.

Thanks a lot,

Yonatan



--
Research Associate
Gazzaley Lab
Department of Neurology
University of California, San Francisco