<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 12 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
        {font-family:Tahoma;
        panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
        {font-family:Consolas;
        panose-1:2 11 6 9 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
p
        {mso-style-priority:99;
        mso-margin-top-alt:auto;
        margin-right:0cm;
        mso-margin-bottom-alt:auto;
        margin-left:0cm;
        font-size:12.0pt;
        font-family:"Times New Roman","serif";}
pre
        {mso-style-priority:99;
        mso-style-link:"HTML Preformatted Char";
        margin:0cm;
        margin-bottom:.0001pt;
        font-size:10.0pt;
        font-family:"Courier New";}
span.HTMLPreformattedChar
        {mso-style-name:"HTML Preformatted Char";
        mso-style-priority:99;
        mso-style-link:"HTML Preformatted";
        font-family:Consolas;}
span.EmailStyle20
        {mso-style-type:personal;
        font-family:"Calibri","sans-serif";
        color:#1F497D;}
span.EmailStyle21
        {mso-style-type:personal-reply;
        font-family:"Calibri","sans-serif";
        color:#1F497D;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-size:10.0pt;}
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
        {page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-GB" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D">After discussion with Donald, it appears that the sequence in use on our GE scanner (a research sequence) does not write the gradient directions to the DICOM
 file so there is no way to extract them from the DICOM.&nbsp; I have always previously used a manual gradient table for Camino / FSL.&nbsp; Since Camino uses its own voxel order format, the rotation does not matter as it assumes the gradients are relative to the acquisition
 plane.&nbsp; Similarly from my experiments, FSL seems to apply the same logic as it does not make any difference what is in the qform/sform header.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D"><o:p>&nbsp;</o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D">However MRtrix does use the transforms in the header so when you tilt the acquisition, the fitting is incorrect unless you change the transform to an identity
 as I suggested or rotate the gradient vectors.&nbsp; Using the identity transform is a quick fix but is not ideal as some of the other processing I do with the same relies on these being present.&nbsp; So does anybody have a ready-made script which takes a set of gradient
 directions and rotates them according to the qform/sform?&nbsp; I seem to remember seeing something like this recently but cannot find in on the MRtrix list so it must have been elsewhere.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D"><o:p>&nbsp;</o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D">Many thanks<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D"><o:p>&nbsp;</o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D">Gavin<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D"><o:p>&nbsp;</o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:&quot;Calibri&quot;,&quot;sans-serif&quot;;color:#1F497D"><o:p>&nbsp;</o:p></span></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:10.0pt;font-family:&quot;Tahoma&quot;,&quot;sans-serif&quot;">From:</span></b><span lang="EN-US" style="font-size:10.0pt;font-family:&quot;Tahoma&quot;,&quot;sans-serif&quot;"> Donald Tournier [mailto:d.tournier@brain.org.au]
<br>
<b>Sent:</b> 16 October 2012 00:18<br>
<b>To:</b> Winston, Gavin<br>
<b>Cc:</b> mrtrix-discussion@www.nitrc.org<br>
<b>Subject:</b> Re: [Mrtrix-discussion] Specifying gradient directions on a GE Signa HDx scanner<o:p></o:p></span></p>
</div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
<p class="MsoNormal">Hi Gavin,<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">Yes, the fun never ends with DICOM...<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">MRtrix was actually coded up based on the same wiki information, but I only implemented support for the latest version, as per the last few entries in the wiki:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:13.5pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;;background:white">In GE DWI images (software version 12.0)</span><o:p></o:p></p>
<div style="border:dashed #2F6FAB 1.0pt;padding:12.0pt 12.0pt 12.0pt 12.0pt">
<pre style="line-height:13.2pt;background:#F9F9F9"><o:p>&nbsp;</o:p></pre>
<pre style="line-height:13.2pt;background:#F9F9F9"><o:p>&nbsp;</o:p></pre>
<pre style="line-height:13.2pt;background:#F9F9F9">(0008,1090) LO [SIGNA EXCITE]<o:p></o:p></pre>
<pre style="line-height:13.2pt;background:#F9F9F9">(0018,1020) LO [12\LX\MR Software release:12.0_M4_0520.a]<o:p></o:p></pre>
</div>
<p style="mso-margin-top-alt:4.8pt;margin-right:0cm;margin-bottom:6.0pt;margin-left:0cm;line-height:20.25pt;background:white">
<span style="font-size:13.5pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;">diffusion directions are stored under the following tags:<o:p></o:p></span></p>
<div style="border:dashed #2F6FAB 1.0pt;padding:12.0pt 12.0pt 12.0pt 12.0pt">
<pre style="line-height:13.2pt;background:#F9F9F9"><o:p>&nbsp;</o:p></pre>
<pre style="line-height:13.2pt;background:#F9F9F9"><o:p>&nbsp;</o:p></pre>
<pre style="line-height:13.2pt;background:#F9F9F9">(0019,10bb) DS [0.430617]<o:p></o:p></pre>
<pre style="line-height:13.2pt;background:#F9F9F9">(0019,10bc) DS [-0.804161]<o:p></o:p></pre>
<pre style="line-height:13.2pt;background:#F9F9F9">(0019,10bd) DS [-0.420008]<o:p></o:p></pre>
</div>
<div>
<p class="MsoNormal"><a name="Private_vendor:_Siemens" id="Private_vendor:_Siemens"></a><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">This has worked fine for all the GE data I've had access to (admittedly not a lot), and it does take care of re-orienting the directions into scanner coordinates. I guess it would be good to include support for these other data, especially
 if this is actually the GE 'standard'. I'd be happy to have a go if you can send me some data to try it out. It sounds like a relatively trivial change, I just need to parse the (0021,105A) tag... Actually, just a dump of the output of &quot;read_dicom -a &lt;DICOM
 file&gt;&quot; should be sufficient for me to figure out how it's stored. If you can send me that, I'll add support for it.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">Otherwise, the directions you get from the headers sound like they would be similar to the FSL bvecs/bvals format. There's been a bit of discussion about how to convert between FSL &amp; MRtrix formats recently, with a few people contributing
 their scripts. A quick google search for MRtrix &amp; bvecs should give you the info you need.&nbsp;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">Finally, I see no reason why you couldn't use the trick of resetting the transform matrix to identity, unless you actually need to preserve the original orientation.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">Hope this helps.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Cheers,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><br>
Donald.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal">On 15 October 2012 23:08, Gavin Winston &lt;<a href="mailto:g.winston@ucl.ac.uk" target="_blank">g.winston@ucl.ac.uk</a>&gt; wrote:<o:p></o:p></p>
<p class="MsoNormal">Dear all,<br>
<br>
I would like to use MRtrix on DTI data from a GE Signa HDx 3.0T scanner (software release 14.0). &nbsp;My query relates to how to specify the gradient directions in the correct reference frame for MRtrix.<br>
<br>
The gradient table file on the scanner specifies the three components of the diffusion gradients relative to the acquisition (i.e. if you tilt the acquisition plane by 10 degrees, the applied gradients also tilt by 10 degrees). &nbsp;I think this behaviour differs
 from Siemens scanners. &nbsp;If I use the gradient table directly in MRtrix and the acquisition is tilted the fibre orientations are also incorrectly tilted. &nbsp;There seem to be several ways to approach this so which do people advise?<br>
<br>
Firstly, extracting the gradient directions from the DICOM header. &nbsp;I understand from previous posts that GE scanners are a bit funny about how they do this. &nbsp;Using read_dicom on our data seems to suggest that it is stored in the fields 0019,10E0 (DTI diffusion
 directions), 0021,105A (diffusion direction), 0043,1039 (b value) as per <a href="http://wiki.na-mic.org/Wiki/index.php/NAMIC_Wiki:DTI:DICOM_for_DWI_and_DTI#Private_vendor:_GE" target="_blank">
http://wiki.na-mic.org/Wiki/index.php/NAMIC_Wiki:DTI:DICOM_for_DWI_and_DTI#Private_vendor:_GE</a>. &nbsp;However mrinfo does not seem to support this as it just produces a gradient table which is entirely zeros.<br>
<br>
Secondly, I could reorientate the entries of the gradient table according to the transforms from the DICOM file (or the qform/sform once converted to Nifti). &nbsp;Are there are ready-made tools to do this?<br>
<br>
Finally, it is possibly to implement a quick fix by modifying the qform/sform headers. &nbsp;If you convert DICOM -&gt; NIFTI, using fslorient as below then convert NIFTI -&gt; MIF with mrconvert things work fine. &nbsp;For example on a 1.875x1.875x2.4mm dataset:<br>
<br>
fslorient -setqform -1.875 0 0 0 0 1.875 0 0 0 0 2.4 0 0 0 0 1 &lt;file&gt;<br>
fslorient -setsform -1.875 0 0 0 0 1.875 0 0 0 0 2.4 0 0 0 0 1 &lt;file&gt;<br>
<br>
Normally our acquisitions are virtually aligned with the scanner axes so the error is small, but I acquired a dataset this morning where the acquisition was tilted by 45 degrees to check the above worked.<br>
<br>
Many thanks,<br>
<br>
Gavin<br>
_______________________________________________<br>
Mrtrix-discussion mailing list<br>
<a href="mailto:Mrtrix-discussion@www.nitrc.org" target="_blank">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><o:p></o:p></p>
</div>
<p class="MsoNormal"><br>
<br clear="all">
<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<p class="MsoNormal">-- <br>
<b><span style="font-size:7.5pt;color:#FF6600">Dr Jacques-Donald Tournier</span></b><o:p></o:p></p>
<div>
<p class="MsoNormal"><span style="font-size:7.5pt;color:#FF6600">Research Fellow</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p>&nbsp;</o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:7.5pt">The Florey Institute of Neuroscience and Mental Health</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:7.5pt">Melbourne Brain Centre - Austin Campus</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:7.5pt">245 Burgundy Street</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:7.5pt">Heidelberg &nbsp;Vic &nbsp;3084</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:7.5pt">Ph: &nbsp;&#43;61 3 9035&nbsp;7033</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:7.5pt">Fax: &nbsp;&#43;61 3 9035 7307</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:7.5pt"><a href="http://www.florey.edu.au" target="_blank">www.florey.edu.au</a></span><o:p></o:p></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><o:p>&nbsp;</o:p></p>
</div>
</div>
</body>
</html>