If you load the volume in with MRIread
a = MRIread('yourvolume.mgz');
a.vox2ras1 is the one-based vox2ras matrix which takes a 1-based column,
row, slice and converts it to RAS. So you just need to convert your
indices to col, row, slice and apply this matrix.
doug
On 12/17/2012 06:53 PM, octavian lie wrote:
> Dear Douglas,
>
> Thank you very much for your input.
>
> 1. I generated the ind matrix as per your script, but what I wanted is
> to get a matrix with the RAS coordinates for the pixels indexed in
> ind. Is there a way to do this?
>
> Thank you,
>
> Octavian
--
Douglas N. Greve, Ph.D.
MGH-NMR Center
gr...@nmr.mgh.harvard.edu
Phone Number: 617-724-2358
Fax: 617-726-7422
Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting
FileDrop: www.nmr.mgh.harvard.edu/facility/filedrop/index.html
Outgoing: ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/
_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
The information in this e-mail is intended only for the person to whom it is
addressed. If you believe this e-mail was sent to you in error and the e-mail
contains patient information, please contact the Partners Compliance HelpLine at
http://www.partners.org/complianceline . If the e-mail was sent to you in error
but does not contain patient information, please contact the sender and properly
dispose of the e-mail.
