<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><span class="Apple-style-span" style="font-family: 'Times New Roman'; font-size: medium; "><pre>Hi Sibin,</pre><pre><br></pre><pre>You can use the ITK to read series of DICOM images, compute the </pre><pre><br></pre><pre>histogram or gradient image, perform the registration and evaluate</pre><pre><br></pre><pre>the registration quality.</pre><pre><br></pre><pre>The "IO" folder in the ITK examples directory has some source code</pre><pre><br></pre><pre>samples for reading DICOM series, also chapter 7 of the software</pre><pre><br></pre><pre>guide contains more information about the samples and the methods used.</pre><pre><br></pre><pre>Computing the histogram can be performed using the histogram generator</pre><pre><br></pre><pre>in ITK. There are code samples included in the "statistics" folder</pre><pre><br></pre><pre>of the examples directory and chapter 10
 of the software guide contains</pre><pre><br></pre><pre>the documentation regarding the samples. Once you compute the histogram&nbsp;</pre><pre><br></pre><pre>you can access the histogram information using one of the following,&nbsp;</pre><pre><br></pre><pre>instance identifiers, n-dimensional indices or iterators&nbsp;(I would&nbsp;</pre><pre><br></pre><pre>recommend the iterators). Gradient filter is also supported by the ITK.&nbsp;</pre><pre><br></pre><pre>Examples are included in the "filtering" folder and chapter 6 contains&nbsp;</pre><pre><br></pre><pre>the documentation for this filter.</pre><pre><br></pre><pre>Evaluating the registration accuracy is a bit more tricky. If you want</pre><pre><br></pre><pre>to use mean squares you have two options, first is to evaluate the mean</pre><pre><br></pre><pre>square metric after the registration by setting the transformation of</pre><pre><br></pre><pre>a registration to identity which as I recall already
 has been discussed&nbsp;</pre><pre><br></pre><pre>in the&nbsp;in the mailing list, the other option is to use the "GetRMSChange()"&nbsp;</pre><pre><br></pre><pre>method supported by some registration functions which gives you the Root&nbsp;</pre><pre><br></pre><pre>Mean Square (RMS) in the deformation filed. The first approach will be&nbsp;</pre><pre><br></pre><pre>something like the "MeanSquareImageMetric" example in the "registration"&nbsp;</pre><pre><br></pre><pre>folder, as for the second approach you can extract the deformation field&nbsp;</pre><pre><br></pre><pre>after the registration and pass it to the registration function as a&nbsp;</pre><pre><br></pre><pre>template parameter (along with the fixed and moving image) then get&nbsp;</pre><pre><br></pre><pre>the RMS value.</pre><pre><br></pre><pre>Hope this helps,</pre><pre><br></pre><pre>Dawood
 Masslawi</pre><pre><br></pre><pre><br></pre><pre>&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;</pre><pre><br></pre><pre><span class="Apple-style-span" style="font-family: 'Times New Roman'; white-space: normal; "><pre>&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;</pre><pre><br></pre></span></pre><pre><br></pre><pre>Hi,
I am new to itk, I am currently doing project on Demons registration for
multi modality using preprocessing approaches such as Histogram matching and
gradient images.

Now let me tell you in detail.......

first,
As we know in CT, values are normalized(calibrated to housenfield unit) but
not in MR.
So to normalize MR image we go for preprocessing.....
Here i need help

two approaches:
Histogram Matching
Use of Gradient images

Please note my input would be reading a *set* of *DICOM* images...CT and MR
from directory not single images.


after normalize i need to do Demons registration and then need to measure
accuracy of demons...................How can i Do that using ITK

i think mean square will show me accuracy of Demons but can i do that using
ITK


Please help....
i know there is code in ITK DeformableRegistation2.cxx in which demons are
used but how to implement my problem....in that code....Please help...!!!</pre></span></td></tr></table><br>