Hi Frederic,<br><br>The different intensities scales of your two images will not<br>be a problem if you use one of the following metrics:<br><br>* Mutual Information<br>* NormalizedCorrelation<br><br>Note that there are about six variations of the Mutual Information<br>
metric in ITK, any one of them will be insensitive to differences<br>in input intensity scaling.<br><br>For a full list of Image Metrics, please look at:<br><a href="http://public.kitware.com/Insight/Doxygen/html/group__RegistrationMetrics.html">http://public.kitware.com/Insight/Doxygen/html/group__RegistrationMetrics.html</a><br>
<br>You could also, of course, rescale the intensities of the input<br>images by using filters such as:<br><br>* RescaleIntensityImageFilter<br>* IntensityWindowingImageFilter<br>* NormalizeImageFilter<br>* HistogramMatchingImageFilter<br>
<br><br>Note also, that, this is independent of the fact that you <br>are using a ThinPlate transform. That is, you would have<br>the same issue, regardless of what Transform, and Optimizer<br>and Interpolator you decide to use.<br>
<br><br> Regards,<br><br><br> Luis<br><br><br><br>------------------------------------------------------<br><div class="gmail_quote">2009/9/9 frédéric salvador <span dir="ltr"><<a href="mailto:frederic.salvador@gmail.com">frederic.salvador@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Dear all,<br><br>I would like to register 3D CT images of the abdomen using the Thin Plate Splines Registration.<br>
My goal is to select a set of landmarks (6 or 7 for instance) both in a FixedImage and MovingImage, compute the transform from the MovingImage to the FixedImage, and apply this transform to the MovingImage (and finally save the RegisteredMovingImage as a ".mha" file).<br>
<br>Actually, I already tried to implement such a process based on the "InsightToolkit-3.14.0/Examples/Registration/ThinPlateSplineWarp.cxx".<br>The problem is that my output moving registered image presents intensities between 0 and 255, as it had been rescaled on 8 bits, whereas both my input images had intensities between lets say -1000 and +2000.<br>
Where does this come from?<br>Is it possible to avoid this?<br><br>Thank you very much in advance for your help.<br><br>Best regards.<br><br>Frederic Salvador<br><br><br>
<br>_____________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at<br>
<a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the ITK FAQ at: <a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
<br></blockquote></div><br>