<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><DIV>Dear All,</DIV>
<DIV>&nbsp;</DIV>
<DIV>Here is a piece of code that I have written to estimate accuracy of demons registration.&nbsp; The code reads&nbsp;index of the landmark points on fixed and moving image from a file and calculates reference vectors.&nbsp;Accuracy is&nbsp;calculated by&nbsp;taking difference between&nbsp;reference vector field and vector field generated by demons at fixed landmark position.&nbsp;Qualitativly or visually the&nbsp;registration&nbsp;results look good (deformed moving&nbsp;image overlaid on to fixed image). However the&nbsp;components of error vector&nbsp;are very high. Is there any&nbsp;coding errors&nbsp;?</DIV>
<DIV>&nbsp;</DIV><FONT size=2>
<DIV>FILE *read_fixed_landmarks = fopen(argv[7],</FONT><FONT color=#800000 size=2><FONT color=#800000 size=2>"rt"</FONT></FONT><FONT size=2>);</DIV>
<DIV>FILE *read_moving_landmarks = fopen(argv[8],</FONT><FONT color=#800000 size=2><FONT color=#800000 size=2>"rt"</FONT></FONT><FONT size=2>);</DIV>
<DIV>FILE *error_file = fopen(argv[9],</FONT><FONT color=#800000 size=2><FONT color=#800000 size=2>"w"</FONT></FONT><FONT size=2>);</DIV>
<DIV></DIV>
<DIV>&nbsp;</DIV>
<DIV>VectorType est_displacement;</DIV>
<DIV>VectorType ref_displacement;</DIV>
<DIV></FONT><FONT color=#0000ff size=2><FONT color=#0000ff size=2></FONT></FONT>&nbsp;</DIV>
<DIV><FONT color=#0000ff size=2><FONT color=#0000ff size=2>double</FONT></FONT><FONT size=2> error[3];</DIV>
<DIV></FONT><FONT size=2>;</DIV>
<DIV></FONT><FONT color=#008000 size=2><FONT color=#008000 size=2>&nbsp;</DIV></FONT></FONT><FONT size=2>
<DIV></FONT><FONT color=#0000ff size=2><FONT color=#0000ff size=2>while</FONT></FONT><FONT size=2> (!feof(read_fixed_landmarks) &amp;&amp; !feof(read_moving_landmarks))</DIV>
<DIV>{</DIV>
<DIV></DIV>
<DIV></FONT><FONT color=#0000ff size=2><FONT color=#0000ff size=2>float</FONT></FONT><FONT size=2> fixed_x, fixed_y, fixed_z;</DIV>
<DIV></FONT><FONT color=#0000ff size=2><FONT color=#0000ff size=2>float</FONT></FONT><FONT size=2> moving_x, moving_y, moving_z;</FONT></DIV>
<DIV><FONT size=2>&nbsp;</DIV>
<DIV>fscanf(read_fixed_landmarks,</FONT><FONT color=#800000 size=2><FONT color=#800000 size=2>"%f %f %f\n"</FONT></FONT><FONT size=2>, &amp;fixed_x, &amp;fixed_y, &amp;fixed_z);</DIV>
<DIV>fscanf(read_moving_landmarks,</FONT><FONT color=#800000 size=2><FONT color=#800000 size=2>"%f %f %f\n"</FONT></FONT><FONT size=2>, &amp;moving_x, &amp;moving_y, &amp;moving_z);</DIV>
<DIV></DIV>
<DIV>&nbsp;</DIV>
<DIV>DeformationFieldType::IndexType index;</DIV>
<DIV>index[0] = fixed_x;</DIV>
<DIV>index[1] = fixed_y;</DIV>
<DIV>index[2] = fixed_z;</DIV>
<DIV></FONT>
<DIV><FONT color=#008000 size=2><FONT color=#008000 size=2></DIV></FONT></FONT><FONT size=2>est_displacement = field-&gt;GetPixel(index);</DIV>
<DIV>ref_displacement[0] = moving_x*moving_spacing[0]- fixed_x*fixed_spacing[0];</DIV>
<DIV>ref_displacement[1] = moving_y*moving_spacing[1]- fixed_y*fixed_spacing[1];</DIV>
<DIV>ref_displacement[2] = moving_z*moving_spacing[2]- fixed_z*fixed_spacing[2];</DIV>
<DIV>&nbsp;</DIV>
<DIV>error[0] = est_displacement[0]- ref_displacement[0];</DIV>
<DIV>error[1] = est_displacement[1]- ref_displacement[1];</DIV>
<DIV>error[2] = est_displacement[2]- ref_displacement[2];</DIV>
<DIV></DIV>
<DIV></DIV>
<DIV></DIV>
<DIV>&nbsp;</DIV>
<DIV>}</DIV>
<DIV>&nbsp;</DIV>
<DIV>Thank you in advance</DIV>
<DIV>Yogish Mallya</DIV></FONT></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>