Hi Don
the -rt nearest switch is for mri_convert. I'll leave the rest of the
questions for Doug :)
Bruce
On Fri, 1 Aug 2014, Krieger, Donald N. wrote:
Thanks for posting back so quickly, Bruce, and for the direction.
Just to clarify, here is the command which I am using so you can see the
types of files on which I am operating.
The working directory is ${FREESURFER_HOME}/subjects/bert .
I'm putting each switch on a separate line:
mri_label2label
--s Sub125
--srclabel label/rh.cortex.label
--trglabel testLabels/rh.cortext.label
--xfm mri/transforms/talairach.xfm
--xfm-invert
--regmethod volume
I plotted the output vectors for a single 1 mm slice using my tools.
This is shown in the left panel of the attached jpg.
The right panel shows the same vectors translated by hand to "fit."
The size and shape are close so this appears to be in the right direction.
I would like to find the freesurfer coords --> MRI coords transform and try
it in place of talairach.xfm .
I assume that the inverse of talairach.xfm is not the correct transform to
get from freesurfer --> MRI coordinates.
Is there a way to extract the correct transform? It must have existed early
in the "recon-all" run and must also be accessible to mri_convert -rl ...
I tried --nt as you suggested but my version of mri_label2label (stable5)
does not recognize it.
I tried --regmethod surface --hemi rh and found that it gives identical
results whether using --xfm-invert or not !
Regards,
Don
Don Krieger, Ph.D.
Department of Neurological Surgery
University of Pittsburgh
(412)648-9654 Office
(412)521-4431 Cell/Text
> -----Original Message-----
> From: freesurfer-boun...@nmr.mgh.harvard.edu [mailto:freesurfer-
> boun...@nmr.mgh.harvard.edu] On Behalf Of Bruce Fischl
> Sent: Friday, August 01, 2014 9:29 AM
> To: Freesurfer support list
> Subject: Re: [Freesurfer] FW: YANQ follow-up
>
> Hi Don
>
> try mri_label2label. Also, the if you are transforming label volumes
> you need to use nearest neighbor interpolation, so include the -rt nearest
switch.
>
> cheers
> Bruce
> On Fri, 1 Aug 2014, Krieger, Donald N. wrote:
>
> >
> > Hi All,
> >
> >
> >
> > I am still in need of finding the tool sequence to generate an ASCII
> > file containing (1) voxel RAS coordinates (2) for a specific label
> > (3) using the original MRI coordinates rather than freesurfer
coordinates.
> >
> >
> >
> > mri_convert -rl rawavg.mgz <input vol> <output vol>
> >
> > does generate a new output volume containing the “labelled” voxels
> > contained in <input vol> . I can tell because this volume shows up
> > more or less correctly in freeview although the color mapping
> > appears to have changed or perhaps freeview is doing some kind of “fill”
> > operation differently depending on the loaded volume. I assume
> > these are in the coordinates for rawavg.mgz but cannot tell for sure.
> >
> >
> >
> > I need to generate a list of these voxel coordinates from <output
> > vol> to check and for my work.
> >
> > (1) Do I need to generate a new annotation file from <output vol>
> > and then use mri_annotation2label to get there? I have not been
> > able to find the tool sequence which generates the .annot files.
> >
> > (2) An alternative appears to be mri_label2label with –xfm
> > xfmfile and possibly including --xfm-invert . I did get this
> > command to complete successfully using mri/transforms/talaraich.xfm
> > but I presume this is not the correct transform.
> >
> >
> >
> > Thanks and best regards,
> >
> >
> >
> > Don
> >
> >
> >
> > Don Krieger, Ph.D.
> >
> > Department of Neurological Surgery
> >
> > University of Pittsburgh
> >
> > (412)648-9654 Office
> >
> > (412)521-4431 Cell/Text
> >
> >
> >
> > From: Krieger, Donald N.
> > Sent: Thursday, July 31, 2014 9:36 AM
> > To: freesurfer@nmr.mgh.harvard.edu
> > Subject: YANQ follow-up
> >
> >
> >
> > Thanks again for the pointer, Bruce.
> >
> >
> >
> > You wrote:
> >
> > … mri_convert -rl rawavg.mgz <input vol> <output vol>
> >
> > which does produce an output volume from which I now need to extract
> > RAS voxel coordinates.
> >
> > A pointer to a keyboard command would be most welcome.
> >
> >
> >
> > Of interest is the difference between the color mappings for the
> > input and output volumes in freeview.
> >
> > When using the lookup table in freeview they are quite different
> > whereas in grayscale they are almost identical.
> >
> > I wonder if this is contributing to my difficulty in finding a tool
> > to extract the voxel coordinates.
> >
> >
> >
> > Thanks again.
> >
> >
> >
> > Regards,
> >
> >
> >
> > Don
> >
> >
> >
> > Don Krieger, Ph.D.
> >
> > Department of Neurological Surgery
> >
> > University of Pittsburgh
> >
> > (412)648-9654 Office
> >
> > (412)521-4431 Cell/Text
> >
> >
> >
> >
> >
_______________________________________________
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.
