Hi Gavin,<div><br></div><div>Yes, the fun never ends with DICOM...</div><div><br></div><div>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:</div>
<div><br></div><div><span style="background-color:rgb(255,255,255);font-family:Arial,sans-serif;font-size:18px;line-height:27px">In GE DWI images (software version 12.0)</span><pre style="padding:1em;border:1px dashed rgb(47,111,171);background-color:rgb(249,249,249);line-height:1.1em">
(0008,1090) LO [SIGNA EXCITE]
(0018,1020) LO [12\LX\MR Software release:12.0_M4_0520.a]
</pre><p style="margin:0.4em 0px 0.5em;line-height:27px;font-family:Arial,sans-serif;font-size:18px;background-color:rgb(255,255,255)">diffusion directions are stored under the following tags:</p><pre style="padding:1em;border:1px dashed rgb(47,111,171);background-color:rgb(249,249,249);line-height:1.1em">
(0019,10bb) DS [0.430617]
(0019,10bc) DS [-0.804161]
(0019,10bd) DS [-0.420008]
</pre><a name="Private_vendor:_Siemens" id="Private_vendor:_Siemens" style="color:rgb(0,43,184);background-image:none;background-color:rgb(255,255,255);font-family:Arial,sans-serif;font-size:18px;line-height:27px"></a><div>
<br></div><div>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 "read_dicom -a <DICOM file>" should be sufficient for me to figure out how it's stored. If you can send me that, I'll add support for it.</div>
<div><br></div><div>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 & MRtrix formats recently, with a few people contributing their scripts. A quick google search for MRtrix & bvecs should give you the info you need. </div>
<div><br></div><div>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.</div><div><br></div><div>Hope this helps.</div>
<div>Cheers,</div><div><br>Donald.</div><div><br></div><div><br></div><div class="gmail_quote">On 15 October 2012 23:08, Gavin Winston <span dir="ltr"><<a href="mailto:g.winston@ucl.ac.uk" target="_blank">g.winston@ucl.ac.uk</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">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). 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). I think this behaviour differs from Siemens scanners. If I use the gradient table directly in MRtrix and the acquisition is tilted the fibre orientations are also incorrectly tilted. 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. I understand from previous posts that GE scanners are a bit funny about how they do this. 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/<u></u>index.php/NAMIC_Wiki:DTI:<u></u>DICOM_for_DWI_and_DTI#Private_<u></u>vendor:_GE</a>. 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). 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. If you convert DICOM -> NIFTI, using fslorient as below then convert NIFTI -> MIF with mrconvert things work fine. 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 <file><br>
fslorient -setsform -1.875 0 0 0 0 1.875 0 0 0 0 2.4 0 0 0 0 1 <file><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>
______________________________<u></u>_________________<br>
Mrtrix-discussion mailing list<br>
<a href="mailto:Mrtrix-discussion@www.nitrc.org" target="_blank">Mrtrix-discussion@www.nitrc.<u></u>org</a><br>
<a href="http://www.nitrc.org/mailman/listinfo/mrtrix-discussion" target="_blank">http://www.nitrc.org/mailman/<u></u>listinfo/mrtrix-discussion</a><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><br>
</div>