<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><DIV>&nbsp;</DIV>
<DIV>Hi Luis,</DIV>
<DIV>&nbsp;</DIV>
<DIV>The&nbsp; landmarks are&nbsp;recorded in terms&nbsp;of pixel coordinates. &nbsp;The origin&nbsp;is (0,0,0) and orientation is trans-axial.</DIV>
<DIV>&nbsp;</DIV>
<DIV>The code mentioned in my mail&nbsp;was correct&nbsp;.I found that the issue was with the <SPAN>paramter values</SPAN>&nbsp;that I had used for demons. I adjusted&nbsp;certain parameters of demons (iterations, histogram matching criteria). Now the qualitative results match with the error vectors.</DIV>
<DIV>&nbsp;</DIV>
<DIV>Best Wishes</DIV>
<DIV>Yogish</DIV>
<DIV>&nbsp;</DIV>
<DIV>&nbsp;<BR><BR>--- On <B>Thu, 21/5/09, Luis Ibanez <I>&lt;luis.ibanez@kitware.com&gt;</I></B> wrote:<BR></DIV>
<BLOCKQUOTE style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: rgb(16,16,255) 2px solid"><BR>From: Luis Ibanez &lt;luis.ibanez@kitware.com&gt;<BR>Subject: Re: [Insight-users] Accuracy of demons registration<BR>To: "Yogish Mallya" &lt;mallyayogish@yahoo.co.in&gt;<BR>Cc: insight-users@itk.org<BR>Date: Thursday, 21 May, 2009, 8:34 AM<BR><BR>
<DIV class=plainMail><BR>Hi Yogish,<BR><BR><BR>&nbsp; Why are you assigning the landmarks to Indexes ?<BR><BR>&nbsp; Are they really defined in pixel coordinates ?<BR><BR>&nbsp; Usually landmarks are recorded in Physical Space<BR>&nbsp; (taking origin,spacing and orientation into account).<BR><BR><BR><BR>&nbsp;&nbsp;&nbsp;Luis<BR><BR><BR><BR>-------------------<BR>Yogish Mallya wrote:<BR>&gt; Dear All,<BR>&gt;&nbsp; Here is a piece of code that I have written to estimate accuracy of demons registration.&nbsp; The code reads index of the landmark points on fixed and moving image from a file and calculates reference vectors. Accuracy is calculated by taking difference between reference vector field and vector field generated by demons at fixed landmark position. Qualitativly or visually the registration results look good (deformed moving image overlaid on to fixed image). However the components of error vector are very high. Is there any coding errors
 ?<BR>&gt;&nbsp; FILE *read_fixed_landmarks = fopen(argv[7],"rt");<BR>&gt; FILE *read_moving_landmarks = fopen(argv[8],"rt");<BR>&gt; FILE *error_file = fopen(argv[9],"w");<BR>&gt;&nbsp; VectorType est_displacement;<BR>&gt; VectorType ref_displacement;<BR>&gt;&nbsp; double error[3];<BR>&gt; ;<BR>&gt;&nbsp; while (!feof(read_fixed_landmarks) &amp;&amp; !feof(read_moving_landmarks))<BR>&gt; {<BR>&gt; float fixed_x, fixed_y, fixed_z;<BR>&gt; float moving_x, moving_y, moving_z;<BR>&gt;&nbsp; fscanf(read_fixed_landmarks,"%f %f %f\n", &amp;fixed_x, &amp;fixed_y, &amp;fixed_z);<BR>&gt; fscanf(read_moving_landmarks,"%f %f %f\n", &amp;moving_x, &amp;moving_y, &amp;moving_z);<BR>&gt;&nbsp; DeformationFieldType::IndexType index;<BR>&gt; index[0] = fixed_x;<BR>&gt; index[1] = fixed_y;<BR>&gt; index[2] = fixed_z;<BR>&gt; est_displacement = field-&gt;GetPixel(index);<BR>&gt; ref_displacement[0] = moving_x*moving_spacing[0]- fixed_x*fixed_spacing[0];<BR>&gt;
 ref_displacement[1] = moving_y*moving_spacing[1]- fixed_y*fixed_spacing[1];<BR>&gt; ref_displacement[2] = moving_z*moving_spacing[2]- fixed_z*fixed_spacing[2];<BR>&gt;&nbsp; error[0] = est_displacement[0]- ref_displacement[0];<BR>&gt; error[1] = est_displacement[1]- ref_displacement[1];<BR>&gt; error[2] = est_displacement[2]- ref_displacement[2];<BR>&gt;&nbsp; }<BR>&gt;&nbsp; Thank you in advance<BR>&gt; Yogish Mallya<BR>&gt; <BR>&gt; <BR>&gt; ------------------------------------------------------------------------<BR>&gt; Explore and discover exciting holidays and getaways with Yahoo! India Travel Click here! &lt;<A href="http://in.rd.yahoo.com/tagline_Travel_1/*http://in.travel.yahoo.com/" target=_blank>http://in.rd.yahoo.com/tagline_Travel_1/*http://in.travel.yahoo.com/</A>&gt;<BR>&gt; <BR>&gt; <BR>&gt; ------------------------------------------------------------------------<BR>&gt; <BR>&gt; _____________________________________<BR>&gt; Powered by
 www.kitware.com<BR>&gt; <BR>&gt; Visit other Kitware open-source projects at<BR>&gt; <A href="http://www.kitware.com/opensource/opensource.html" target=_blank>http://www.kitware.com/opensource/opensource.html</A><BR>&gt; <BR>&gt; 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>&gt; <BR>&gt; Follow this link to subscribe/unsubscribe:<BR>&gt; <A href="http://www.itk.org/mailman/listinfo/insight-users" target=_blank>http://www.itk.org/mailman/listinfo/insight-users</A><BR></DIV></BLOCKQUOTE></td></tr></table><br>


      <!--3--><hr size=1></hr> Explore and discover exciting holidays and getaways with Yahoo! India Travel <a href="http://in.rd.yahoo.com/tagline_Travel_1/*http://in.travel.yahoo.com/"> Click here!</a>