18 Dec
2012
18 Dec
'12
12:06 p.m.
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.
- 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
greve@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/