<div dir="ltr">Hi Stefan,<div><br></div><div>You where right, there was a problem with the voxel spacing. Once again, I made a mistake assuming the aquisition of the data had gone according to the protocol. In this case, the voxel size was a bit off (1.17 instead of 0.93 mm in-plane), which messed up the euler based registration. After I made this parameter dynamic as well, the registration went as expected, with acceptable SSD.</div>
<div><br></div><div>Thanks for your feedback, and for the time and effort you guys put into this great program!</div><div><br></div><div>Best,</div><div><br></div><div>Tom</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On Fri, May 9, 2014 at 1:18 PM, S. Klein <span dir="ltr"><<a href="mailto:s.klein@erasmusmc.nl" target="_blank">s.klein@erasmusmc.nl</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div bgcolor="#FFFFFF" text="#000000">
Hi Tom,<br>
Interesting. Are some of your CT scans perhaps made with gantry angle? If it is not taken into account in the direction cosines, it could explain a shearing transformation. Or perhaps a mistake is made somewhere with the voxel spacing in one of the images?
 That could explain why you need a global scaling.<br>
<br>
To investigate why affine gives so much better results than rigid, it could be good to inspect the parameters of the resulting affine transformation. Does the resulting affine transform matrix induce scaling/shearing? Or is it equivalent to a rigid transformation?<br>

<br>
Best,<br>
Stefan<div><div class="h5"><br>
<br>
<div>On 08/05/2014 08:50, Tom Knoop wrote:<br>
</div>
<blockquote type="cite">
<div dir="ltr"><span style="font-family:arial,sans-serif;font-size:13px">Dear Kasper, all,</span>
<div style="font-family:arial,sans-serif;font-size:13px"><br>
</div>
<div style="font-family:arial,sans-serif;font-size:13px">Thank you for your suggestion. I am a bit hesitant to use affine transformations since I try to do intra-patient registrations, which should work with rigid transformations only. In the end (after I cleared
 another problem which had to do with me assuming the patient position was alwas feet first, thus messing up the dicom import into matlab I made), I build a workflow that performs a rigid transformation first, and if the final result is higher than a certain
 threshold, continues with a affine transformation. I am quite happy with the results, so thank you! I have included a visualisation for the registration process for both a failed and working rigid transformation, and the consecutive working affine transformation.
 I included the affine parameterfile I used as well. The files can be found here:</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br>
</div>
<div style="font-family:arial,sans-serif;font-size:13px"><a href="https://filesender.surfnet.nl/?vid=14ca4c6e-3339-6728-72d1-000055597207" target="_blank">https://filesender.surfnet.nl/?vid=14ca4c6e-3339-6728-72d1-000055597207</a><br>

</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br>
</div>
<div style="font-family:arial,sans-serif;font-size:13px">This leaves me with the question why the rigid transformation does not work in all intra-patient cases. By the way, I tried to manually register some points, and do a multi metric transformation as well,
 this had no succes (although my makeshift landmark selection might have been a bit off as wel) If anyone has additional insights in this, I am very interested.</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br>
</div>
<div style="font-family:arial,sans-serif;font-size:13px">Best,</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br>
</div>
<div style="font-family:arial,sans-serif;font-size:13px">Tom</div>
</div>
<div class="gmail_extra"><br>
<br>
<div class="gmail_quote">On Sat, May 3, 2014 at 2:08 PM, Kasper Marstal <span dir="ltr">
<<a href="mailto:kaspermarstal@gmail.com" target="_blank">kaspermarstal@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">Dear Tom,
<div><br>
</div>
<div>I had the same problem when initializing registration of MRIs of heads using a rigid transformation (I can see you have specified "EulerTransform" in your parameter file). In my case, some heads would align nicely and other would seem to get stuck in random
 places yet close to the solution. This typically happens when you apply rigid transformations (scaling, rotation, translation) to registrations where they cannot account for anatomical differences between scans that violates this rigid body assumption (fx
 in inter-patient registrations). </div>
<div><br>
</div>
<div>I solved my problem using an affine transform which also allows for stretching and shearing of the deformation field. This typically works well for bones. Use this transform by changing the line in your parameter file from <span style="color:rgb(0,0,0);font-family:monospace;font-size:13.63636302947998px">(Transform
 "EulerTransform")</span> to <span style="color:rgb(0,0,0);font-family:monospace;font-size:14px">(Transform "AffineTransform")</span>. </div>
<div><br>
</div>
<div>If you are feeling adventurous you might even consider subsequent application of a non-rigid transform, such as the  <span style="color:rgb(0,0,0);font-family:monospace;font-size:13.63636302947998px">(Transform "BSplineTransform")</span>, to capture the
 remaining differences near bone edges that the affine transformation is unable to recover. This really depends on the assumptions about your problem that you are willing to make and/or sacrifice, however.</div>
<div><br>
</div>
<div>Let us know how it goes!</div>
<div><br>
</div>
<div>Cheers,</div>
<div>Kasper</div>
</div>
<div class="gmail_extra"><br>
<br>
<div class="gmail_quote">
<div>
<div>On 30 April 2014 12:10, Tom Knoop <span dir="ltr"><<a href="mailto:t.h.knoop@student.utwente.nl" target="_blank">t.h.knoop@student.utwente.nl</a>></span> wrote:<br>
</div>
</div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div>
<div dir="ltr"><span style="font-family:arial,sans-serif;font-size:13px">Dear All,</span>
<div style="font-family:arial,sans-serif;font-size:13px"><br>
</div>
<div style="font-family:arial,sans-serif;font-size:13px">I am trying to register the long bones in the leg (femur) in multiple CT scans of the same patients. There are a number of reasons to do this, but one of them is that I would like to get a segmentation
 of the consecutive femurs by transforming the femur segmentation I made for the first image to the other images in the series. Because the rotation of the legs are different in every image, I do this per femur</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br>
</div>
<div style="font-family:arial,sans-serif;font-size:13px">To do this, I use the segmentation of the first image as a fixed mask, and mask out the controlateral side of the body in my moving mask. I have build a parameter file that uses AdvancedMeanSquares metric
 and five smoothing steps with a high number of iterations to make the process more robust. See also parameterfile rigid_bone_forward.txt in the zip-file which you can find at the download-link here:
<a href="https://filesender.surfnet.nl/?vid=1dbe6839-ef04-b228-03e8-00001a643b33" target="_blank">
https://filesender.surfnet.nl/?vid=1dbe6839-ef04-b228-03e8-00001a643b33</a></div>
<div style="font-family:arial,sans-serif;font-size:13px"><br>
</div>
<div style="font-family:arial,sans-serif;font-size:13px">This works fine for most of my images, with final SSD of around 5000. However, on some images, the registration just doesn't seem to work. It does get in the right direction (the femurs definitely have
 more overlap after the registration than before) but the voxel-accuracy I normally get is never reached. I have tried fooling around with the sp_A parameters and the resolution steps and maximum stepsizes, but to no avail. Also in the zip file are the logfiles
 of a successfull and failed registration.</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br>
</div>
<div style="font-family:arial,sans-serif;font-size:13px">Unfortunately, I cannot share the scans themselves. However, it are CT scans of the lower body (containing complete femurs, and part of the knee joint and pelvis) made on an philips brilliance big bore
 scanner, with a resolution of 512x512x±210 voxels with a voxel size of .9375x.9375x3mm. The images are saved as mhd double datafiles, with their anatomical position taken from the original dicom data. masks are saved as binary mhd files.<br>

</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br>
</div>
<div style="font-family:arial,sans-serif;font-size:13px">Has anyone have any idea about where this goes wrong? Visually, I do not see a difference between the failed and succesfull scans. It is possible that a single scan fails, while another scan for the same
 patient does register succesfully. What can I do to make this more robust?</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br>
</div>
<div style="font-family:arial,sans-serif;font-size:13px">Many thanks for your help!</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br>
</div>
<div style="font-family:arial,sans-serif;font-size:13px">Best,</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br>
</div>
<div style="font-family:arial,sans-serif;font-size:13px">Tom</div>
</div>
<br>
</div>
</div>
_______________________________________________<br>
Elastix mailing list<br>
<a href="mailto:Elastix@bigr.nl" target="_blank">Elastix@bigr.nl</a><br>
<a href="http://lists.bigr.nl/cgi-bin/mailman/listinfo/elastix" target="_blank">http://lists.bigr.nl/cgi-bin/mailman/listinfo/elastix</a><br>
<br>
</blockquote>
</div>
<br>
</div>
<br>
_______________________________________________<br>
Elastix mailing list<br>
<a href="mailto:Elastix@bigr.nl" target="_blank">Elastix@bigr.nl</a><br>
<a href="http://lists.bigr.nl/cgi-bin/mailman/listinfo/elastix" target="_blank">http://lists.bigr.nl/cgi-bin/mailman/listinfo/elastix</a><br>
<br>
</blockquote>
</div>
<br>
</div>
<br>
<fieldset></fieldset> <br>
<pre>_______________________________________________
Elastix mailing list
<a href="mailto:Elastix@bigr.nl" target="_blank">Elastix@bigr.nl</a>
<a href="http://lists.bigr.nl/cgi-bin/mailman/listinfo/elastix" target="_blank">http://lists.bigr.nl/cgi-bin/mailman/listinfo/elastix</a>
</pre>
</blockquote>
<br>
</div></div><span class="HOEnZb"><font color="#888888"><pre cols="72">-- 
Stefan Klein
<a href="tel:%2B31%2010%207043442" value="+31107043442" target="_blank">+31 10 7043442</a>
<a href="http://www.bigr.nl/people/StefanKlein" target="_blank">http://www.bigr.nl/people/StefanKlein</a></pre>
</font></span></div>

</blockquote></div><br></div>