Dear all,
I need to compute the geodesic distance (i.,e. the distance on the cortical
surface) between a series of vertices. Specifically, I want to compute the
distance between all the vertices in the right hemisphere (>100000) and a
series of vertices of interest (around 20000), always located in the right
hemisphere. Ideally, my output should be a matrix
right_verticesXvertices_of_interest where in each entry there is the geodesic
distance between the pair of vertices. And I need to do this for hundreds of
subjects!
I tried to use mris_pmake:
for i = right_vertices; for j=vertices_of_interest; command1 =
['/Applications/freesurfer_v5.3/bin/mris_pmake --subject 001 --hemi rh
--surface0 white --curv0 sulc --curv1 sulc --mpmOverlay euclidean --mpmProg
pathFind --mpmArgs startVertex:'num2str(i) ',' 'endVertex:' num2str(j)]
T(i)=evalc('system([command1])') endend
But the problem is that it is way too much time consuming and it is difficult
to store the results in a matrix.
Is there a way to adapt this command for such high computations? I Are there
other ways to obtain the geodesic distance between vertices, either than
mris_pmake (for example a matlab function)?
Thank you in advance,Marica_______________________________________________
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.
