<div dir="ltr">Hi Philip,<div>many thanks for the answers.</div><div><br></div><div>Best,</div><div>Alessandro</div></div><div class="gmail_extra"><br><div class="gmail_quote">2014-09-10 18:07 GMT+02:00 Philip A Cook <span dir="ltr"><<a href="mailto:cookpa@mail.med.upenn.edu" target="_blank">cookpa@mail.med.upenn.edu</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
The offset issue is a legacy thing, internally the coordinate (0.0, 0.0, 0.0) is at the edge of the voxel and not the center. Since the shift is there, FSL must place (0.0, 0.0, 0.0) at the voxel center, which is how most image processing software would do it.<br>
<br>
The Jacobian thing is standard, it's done so that an identity transform (zero displacement vectors) has a determinant of 1. When you take the log, this becomes 0. The log determinant is usually used to express contraction (logjac < 0) or expansion (logjac > 0) or no change (logjac = 0).<br>
<div><div class="h5"><br>
<br>
<br>
On Sep 1, 2014, at 9:34 AM, Alessandro Calamuneri <<a href="mailto:alecalamuneri@gmail.com">alecalamuneri@gmail.com</a>> wrote:<br>
<br>
> I camino experts,<br>
> I am writing you because I would like to better understand how reorient works on tensors. I have checked the source code, and two points are not clear to me: first one is that, when using displacement field coming from FSL fnirt command, there is a sistematic shifting on voxel position of half unit; indeed you typed down<br>
><br>
> double posX = (voxX+0.5)*CL_Initializer.voxelDims[0];<br>
> double posY = (voxY+0.5)*CL_Initializer.voxelDims[1];<br>
> double posZ = (voxZ+0.5)*CL_Initializer.voxelDims[2];<br>
><br>
> Why does this shifting occur? Is it an fsl related issue?<br>
> The other question is coming from jacobian construction. Why do you add an identity matrix to the Jacobian after getting directional derivatives of displacement field?<br>
><br>
> Best Regards,<br>
><br>
> Alessandro<br>
</div></div>> _______________________________________________<br>
> Camino-users mailing list<br>
> <a href="mailto:Camino-users@www.nitrc.org">Camino-users@www.nitrc.org</a><br>
> <a href="http://www.nitrc.org/mailman/listinfo/camino-users" target="_blank">http://www.nitrc.org/mailman/listinfo/camino-users</a><br>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div dir="ltr">Alessandro Calamuneri<br>
<div>M.Sc., PhD-student in Ricerca Clinica e Traslazionale in Neuroscienze ed Oncologia</div>
<div>Department of Neurosciences, University of Messina, Italy</div>
<div>Via Consolare Valeria, Policlinico Universitario, 98125, Messina</div>
<div><br>mail: <a href="mailto:alecalamuneri@gmail.com" target="_blank">alecalamuneri@gmail.com</a><br></div><div> <a href="mailto:alessandro.calamuneri@unime.it" target="_blank">alessandro.calamuneri@unime.it</a></div><div> <a href="mailto:acalamuneri@unime.it" target="_blank">acalamuneri@unime.it</a><br><br></div></div>
</div>