<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:arial, helvetica, sans-serif;font-size:10pt"><div>Hi Luis,</div><div style="font-family:arial, helvetica, sans-serif;font-size:10pt"><div style="font-family:times new roman, new york, times, serif;font-size:12pt"><div style="font-family:arial, helvetica, sans-serif;font-size:10pt;"><div><br></div><div>I am already using the&nbsp;</div><div><br></div><div><span class="Apple-style-span" style="font-family:'times new roman';font-size:16px;">&nbsp;registration-&gt;SetFixedImageRegion( region )&nbsp;</span></div><div><font class="Apple-style-span" face="'times new roman'" size="4"><span class="Apple-style-span" style="font-size:16px;"><br></span></font></div><div><font class="Apple-style-span" face="'times new roman'" size="4"><span class="Apple-style-span" style="font-size:16px;">My scaling parameters are around 0.6 and&nbsp;rotation of 6-7
 degrees,&nbsp;I have to use a big region&nbsp;&nbsp;in the fixed image&nbsp;(I am initialiazing with scale =0.5, and hence using a region of size 1280*960). &nbsp;The performance has not improved
 much.&nbsp;</span></font></div><div><font class="Apple-style-span" face="'times new roman'" size="4"><span class="Apple-style-span" style="font-size:16px;"><br></span></font></div><div><font class="Apple-style-span" face="'times new roman'" size="4"><span class="Apple-style-span" style="font-size:16px;">Hence I want to <b>sample pixels from the moving image,</b> and do a <b>backward mapping to pixels in the fixed image</b>( by the inverse transform parameters, that are estimated in each iteration of the optimizer). In this case, I think that almost all the points will be mapped, and hence the metric will be better. Also, it will take care of the scale and translation, as in the previous case, a wrong initial scale will make me miss out the the pixel values in the fixed image.</span></font></div><div><font class="Apple-style-span" face="'times new roman'" size="4"><span class="Apple-style-span" style="
font-size:16px;"><br></span></font></div><div><font class="Apple-style-span" face="'times new roman'" size="4"><span class="Apple-style-span" style="font-size:16px;"><br></span></font></div><div><font class="Apple-style-span" face="'times new roman'" size="4"><span class="Apple-style-span" style="font-size:16px;">Do I need to change the functions in the respective files, or is there a simpler way to do this, or is already implemented?<br><br>Thanks,<br><br>Sharath Venkatesha</span></font></div><div style="font-family:arial, helvetica, sans-serif;font-size:10pt;"><br><div style="font-family:times new roman, new york, times, serif;font-size:12pt;"><font size="2" face="Tahoma"><hr size="1"><b><span style="font-weight:bold;">From:</span></b> Luis Ibanez &lt;luis.ibanez@kitware.com&gt;<br><b><span style="font-weight:bold;">To:</span></b> Sharath Venkatesha &lt;sharath20284@yahoo.com&gt;<br><b><span style="font-weight:bold;">Cc:</span></b> Insight users
 &lt;insight-users@itk.org&gt;<br><b><span style="font-weight:bold;">Sent:</span></b> Saturday, July 11, 2009 8:49:51 AM<br><b><span style="font-weight:bold;">Subject:</span></b> Re: [Insight-users] Registration - Mutual Information Metric -  Sampling from Moving Image<br></font><br>
<br>Hi Sharath,<br><br>You are right in that it is not a good idea to take samples<br>from the entire 2000 x 2000 pixels image if we expect that<br>only a regions of 640 x 480 is going to match the Moving<br>image.<br><br>
The solution is to simply call <br><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; registration-&gt;SetFixedImageRegion( region )<br><br>&nbsp;and set "region" to the region of interest from the <br>2000 x 2000 image, where you anticipate to find matches<br>
for the moving image.<br><br><br>You will find examples on the use of <br><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; SetFixedImageRegion()<br><br>in the ITK software Guide<br><br><span><span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <a target="_blank" href="http://www.itk.org/ItkSoftwareGuide.pdf">http://www.itk.org/ItkSoftwareGuide.pdf</a></span></span><br>
<br>and most (if not all) the corresponding Image Registration<br>examples in <br><br>&nbsp; &nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Insight/Examples/Registration<br><br><br><br><br>&nbsp;&nbsp; Regards,<br><br><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Luis<br><br><br><br>---------------------------------------------------------------<br>
<div class="gmail_quote">On Fri, Jul 10, 2009 at 9:15 PM, Sharath Venkatesha <span dir="ltr">&lt;<a rel="nofollow" ymailto="mailto:sharath20284@yahoo.com" target="_blank" href="mailto:sharath20284@yahoo.com">sharath20284@yahoo.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="border-left:1px solid rgb(204, 204, 204);margin:0pt 0pt 0pt 0.8ex;padding-left:1ex;">
<br>
Hi,<br>
<br>
I am working with images where the Fixed image is 2000*2000 pixels and the moving image is 640*480 (sub region of the fixed image). In all the Mutual Information metrics, the sampling is done on the Fixed image, and the transformed points are mapped on the moving image.<br>

<br>
I figured out that it is incorrect to sample the points randomly in the full fixed image ( as in my case, the moving image is a sub-part of the fixed image, and only 30-40% points were mapped), I used an approximate region of the fixed image from which sampling should be done. This has not improved the performance/accuracy of registration.<br>

<br>
I am thinking of sampling pixels from the moving image, and do a backward mapping to pixels in the fixed image. In this case, I think that almost all the points will be mapped, and hence the metric will be better.<br>
<br>
Do I need to change the functions in the respective files, or is there a simpler way to do this, or is already implemented?<br>
<br>
Thanks,<br>
<br>
Sharath Venkatesha<br>
<br>
<br>
<br>
<br>
<br>
_____________________________________<br>
Powered by <a rel="nofollow" target="_blank" href="http://www.kitware.com">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at<br><span><span>
<a target="_blank" href="http://www.kitware.com/opensource/opensource.html">http://www.kitware.com/opensource/opensource.html</a></span></span><br>
<br><span><span>
Please keep messages on-topic and check the ITK FAQ at: <a target="_blank" href="http://www.itk.org/Wiki/ITK_FAQ">http://www.itk.org/Wiki/ITK_FAQ</a></span></span><br>
<br>
Follow this link to subscribe/unsubscribe:<br><span><span>
<a target="_blank" href="http://www.itk.org/mailman/listinfo/insight-users">http://www.itk.org/mailman/listinfo/insight-users</a></span></span><br>
</blockquote></div><br>
</div></div><div style=""></div></div><br>

      </div></div><div style="position:fixed"></div></div><br>

      </body></html>