<div>Hi Sangma,<br><br>OK, there's quite a few issues at play here. First off, the convention used in MRtrix is that the gradients are specified with respect to the DICOM patient coordinate system. This is essentially the same convention as for DICOM images produced on Siemens and Philips scanners, and the newer standard DICOM tags recently introduced, with one notable exception: the DICOM patient-centered coordinate system has its x-axis running right-to-left and its y-axis running anterior-posterior, whereas MRtrix assumes these axes run in the opposite direction (as per the NIfTI standard). So MRtrix will flip the x & y components to account for this. </div>
<div><br></div><div>For Siemens & Philips images, that's essentially all that happens to the DW gradients. For GE images, it's a little different, since in this case the gradients are (or at least, used to be) stored with respect to the images axes. This means that to convert to MRtrix convention, they must be re-oriented back into the patient-centered coordinate system. This is what the rotate_DW_scheme flag specifies: when the DW gradient information was read from a GE-specific tag, this is set to true and the gradients are then rotated according to the direction of the images axes. In this case, the z component is inverted, since flipping the x & y is essentially the same as flipping z, given the symmetry of diffusion.</div>
<div><br>As far as I know, dcm2nii provides the DW directions with respect to the image axes (in the NIfTI coordinate system). To convert these back to MRtrix convention, the same transformation as you highlighted is needed, but <i>without</i> the negative sign. The transform you highlighted also accounts for differences in the directions of these axes due to the different conventions used for DICOM versus NifTI/MRtrix, so is not appropriate as-is to do convert a dcm2nii-supplied gradient table into one suitable for MRtrix. This is because the dcm2nii gradient table is provided with respect to the image axes in NIfTI coordinate space, and in MRtrix they are provided with respect to the original axes, also in NIfTI coordinate space. <br>
<br>To get back to your question, what do you mean when you say 'the original gradient table'? The one supplied by dcm2nii (in which case converting using the highlighted section of code would not work, as explained above), or the actual ones as read directly from the DICOM headers? If the latter, then unless your images were acquired on a GE scanner, then that also would not work, since that section of code would only be used for GE images. If you are trying to convert non-GE DW directions read directly from the DICOM headers, then all you need to do is flip the sign of the x & y components. </div>
<div><br></div><div>Basically, it's really hard to keep track of all the possible conventions and what each step in the processing pipeline might have done to the DW gradients. Hopefully the above will help you figure out where the problem might be...</div>
<div><br></div><div>Hope this helps.</div><div>Cheers,</div><div><br></div><div>Donald.</div><div><br></div><div><br><div class="gmail_quote">On 3 April 2013 12:57, smxie_nlpr <span dir="ltr"><<a href="mailto:smxie@nlpr.ia.ac.cn" target="_blank">smxie@nlpr.ia.ac.cn</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><u></u>
<div style="MARGIN:10px">
<div>
<div><span style="TEXT-TRANSFORM:none;TEXT-INDENT:0px;LETTER-SPACING:normal;DISPLAY:inline!important;FONT:12pt Times New Roman;WHITE-SPACE:normal;FLOAT:none;WORD-SPACING:0px">Hi,</span></div><div class="im">
<div><span style="TEXT-TRANSFORM:none;TEXT-INDENT:0px;LETTER-SPACING:normal;DISPLAY:inline!important;FONT:12pt Times New Roman;WHITE-SPACE:normal;FLOAT:none;WORD-SPACING:0px"></span> </div>
<div><span style="TEXT-TRANSFORM:none;TEXT-INDENT:0px;LETTER-SPACING:normal;DISPLAY:inline!important;FONT:12pt Times New Roman;WHITE-SPACE:normal;FLOAT:none;WORD-SPACING:0px">I
have a question about the gradient direction table used in Mrtrix. The gradient
direction table achieved from the DICOM images with dcm2nii is different from
the table achieved with the command 'mrinfo DICOMDIR -grad encoding.b'. I
checked the source code of Mrtrix and found the following codes '</span></div>
<div><span style="TEXT-TRANSFORM:none;TEXT-INDENT:0px;LETTER-SPACING:normal;DISPLAY:inline!important;FONT:12pt Times New Roman;WHITE-SPACE:normal;FLOAT:none;WORD-SPACING:0px"></span><span style="TEXT-TRANSFORM:none;TEXT-INDENT:0px;LETTER-SPACING:normal;DISPLAY:inline!important;FONT:12pt Times New Roman;WHITE-SPACE:normal;FLOAT:none;WORD-SPACING:0px">if(rotate_DW_scheme)</span></div>
<div style="TEXT-TRANSFORM:none;TEXT-INDENT:0px;LETTER-SPACING:normal;FONT:medium Tahoma;WHITE-SPACE:normal;WORD-SPACING:0px">
<div style="FONT-FAMILY:Times New Roman;FONT-SIZE:12pt">{</div>
<div style="FONT-FAMILY:Times New Roman;FONT-SIZE:12pt">G(n,0) =
image_transform(0,0) * d[0] + image_transform(0,1) * d[1] - image_transform(0,2)
* d[2];</div>
<div style="FONT-FAMILY:Times New Roman;FONT-SIZE:12pt">G(n,1) =
image_transform(1,0) * d[0] + image_transform(1,1) * d[1] - image_transform(1,2)
* d[2];</div>
<div style="FONT-FAMILY:Times New Roman;FONT-SIZE:12pt">G(n,2) =
image_transform(2,0) * d[0] + image_transform(2,1) * d[1] - image_transform(2,2)
* d[2];</div></div>
<div><span style="TEXT-TRANSFORM:none;TEXT-INDENT:0px;LETTER-SPACING:normal;DISPLAY:inline!important;FONT:12pt Times New Roman;WHITE-SPACE:normal;FLOAT:none;WORD-SPACING:0px">}</span><span style="TEXT-TRANSFORM:none;TEXT-INDENT:0px;LETTER-SPACING:normal;DISPLAY:inline!important;FONT:12pt Times New Roman;WHITE-SPACE:normal;FLOAT:none;WORD-SPACING:0px">'
in image.cpp.</span></div>
<div><span style="TEXT-TRANSFORM:none;TEXT-INDENT:0px;LETTER-SPACING:normal;DISPLAY:inline!important;FONT:12pt Times New Roman;WHITE-SPACE:normal;FLOAT:none;WORD-SPACING:0px">I
transformed the original gradient table with similar code in Matlab, but the
result is not equal to the table achieved with mrinfo. So I want to know the
procedure of transforming the original gradients to the gradients achieved with
mrinfo. And why is it '</span><span style="TEXT-TRANSFORM:none;TEXT-INDENT:0px;LETTER-SPACING:normal;DISPLAY:inline!important;FONT:12pt Times New Roman;WHITE-SPACE:normal;FLOAT:none;WORD-SPACING:0px">G(n,0)
= image_transform(0,0) * d[0] + image_transform(0,1) * d[1] -
image_transform(0,2) * d[2];</span><span style="TEXT-TRANSFORM:none;TEXT-INDENT:0px;LETTER-SPACING:normal;DISPLAY:inline!important;FONT:12pt Times New Roman;WHITE-SPACE:normal;FLOAT:none;WORD-SPACING:0px">'
rather than '</span><span style="TEXT-TRANSFORM:none;TEXT-INDENT:0px;LETTER-SPACING:normal;DISPLAY:inline!important;FONT:medium Tahoma;WHITE-SPACE:normal;FLOAT:none;WORD-SPACING:0px"><span style="FONT-FAMILY:Times New Roman;FONT-SIZE:12pt">G(n,0) =
image_transform(0,0) * d[0] + image_transform(0,1) * d[1] + image_transform(0,2)
* d[2];' ?</span></span></div></div></div>
<div style="FONT-FAMILY:Times New Roman;FONT-SIZE:12pt"> </div>
<div style="FONT-FAMILY:Times New Roman;FONT-SIZE:12pt">Thanks.</div><div class="im">
<div style="FONT-FAMILY:Times New Roman;FONT-SIZE:12pt"> </div>
<div style="FONT-FAMILY:Times New Roman;FONT-SIZE:12pt">Sangma Xie</div>
<hr style="WIDTH:210px;min-height:1px" align="left" color="#b5c4df" size="1">
<div><span>
<div><span style="font-size:10.5pt;font-family:宋体"><span style="FONT-FAMILY:Times New Roman"><span style="FONT-SIZE:10.5pt"><span style="FONT-SIZE:10.5pt">Sangma Xie , Master </span> </span>
<div><span style="FONT-FAMILY:Times New Roman">Brainnetome Center
<div style="FONT-FAMILY:Times New Roman">National Laboratory of Pattern
Recognition (NLPR)</div>Institute of Automation, Chinese Academy of Sciences
(CASIA)<br>95 Zhong Guan Cun East Road, Hai Dian District, Beijing 100190,
P.R.China</span></div></span><br><br></span></div></span></div></div></div>
<br>_______________________________________________<br>
Mrtrix-discussion mailing list<br>
<a href="mailto:Mrtrix-discussion@www.nitrc.org">Mrtrix-discussion@www.nitrc.org</a><br>
<a href="http://www.nitrc.org/mailman/listinfo/mrtrix-discussion" target="_blank">http://www.nitrc.org/mailman/listinfo/mrtrix-discussion</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><font color="#ff6600" size="1"><b>Dr Jacques-Donald Tournier<br></b></font><div><font color="#ff6600" size="1">Research Fellow</font></div><div><font size="1"><br>
</font></div><div><font size="1">The Florey Institute of Neuroscience and Mental Health</font></div><div><font size="1">Melbourne Brain Centre - Austin Campus</font></div><div><font size="1">245 Burgundy Street</font></div>
<div><font size="1">Heidelberg Vic 3084</font></div><div><font size="1">Ph: +61 3 9035 7033</font></div><div><font size="1">Fax: +61 3 9035 7307</font></div><div><font size="1"><a href="http://www.florey.edu.au" target="_blank">www.florey.edu.au</a></font></div>
<br>
</div>