<br>Hi Raheleh,<br><br><br>One of the powerful features of Open Source<br>software is that you can study the software <br>in order to fully understand what it does.<br><br><br> "Glimpsing at the Source<br>
leaves no doubt"<br><br><br><br>Please look at the files:<br><br>Insight/Code/Algorithms/<br> itkMultiResolutionPyramidImageFilter.h<br> itkMultiResolutionPyramidImageFilter.txx<br> itkRecursiveMultiResolutionPyramidImageFilter.h<br>
itkRecursiveMultiResolutionPyramidImageFilter.txx<br><br><br>You will find that these classes use the<br><br> DiscreteGaussianImageFilter<br><br>to smooth the images, an then uses the<br><br><br> ShrinkImageFilter <br>
<br>for doing subsampling.<br><br><br><br> Regards,<br><br><br> Luis<br><div class="gmail_quote"><br><br>-------------------------------------------------------------------------<br>On Mon, Jun 21, 2010 at 1:34 AM, raheleh m <span dir="ltr"><<a href="mailto:meshgolok@yahoo.com">meshgolok@yahoo.com</a>></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;"><div><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;">
<span style="font-weight: bold;">Dear Mr. Ibanez</span><br><br>Thank you so much for your attention. I am really grateful for your help.<br>I am using 3 Resolution levels and the subsampling schedule is<br>[8 8 1], [4 4 1], [2 2 1] for Moving Image and [4 4 1],[2 2 1], [1 1 1] for Fixed Image. <br>
It would be really helpful if you could please help me to understand what prepare pyramids do? I mean for down sampling, what kind of filters does it use? Can I improve the procedure by removing some of these filters or it would affect the results?<br>
<br><span style="font-style: italic;">Best Regards,</span><br>Sincerely Yours<br><br>Rahil Meshgolok<br><br><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><font face="Tahoma" size="2"><div class="im">
<hr size="1"><b><span style="font-weight: bold;">From:</span></b> Luis Ibanez <<a href="mailto:luis.ibanez@kitware.com" target="_blank">luis.ibanez@kitware.com</a>><br><b><span style="font-weight: bold;">To:</span></b> raheleh m <<a href="mailto:meshgolok@yahoo.com" target="_blank">meshgolok@yahoo.com</a>><br>
<b><span style="font-weight: bold;">Cc:</span></b> itk <<a href="mailto:Insight-users@itk.org" target="_blank">Insight-users@itk.org</a>><br></div><b><span style="font-weight: bold;">Sent:</span></b> Sun, June 20, 2010 5:30:20 PM<br>
<b><span style="font-weight: bold;">Subject:</span></b> Re: Thanks, Another Question<br></font><div><div></div><div class="h5"><br>
Hi Rahil<br><br>1) Before you start measuring computation time,<br> please make sure that you configured ITK to<br> build for "Release". This is done by running<br> CMake and setting the variable<br><br>
CMAKE_BUILD_TYPE<br><br>2) Regarding the time for building the pyramid,<br> 170 seconds looks a bit too long...<br><br> How many resolution levels are you using ?<br> What is the subsampling schedule ?<br>
<br><br><br> Please let us know,<br><br><br> Thanks<br><br><br> Luis<br><br>-------------------------------------------<br><div class="gmail_quote">On Sun, Jun 20, 2010 at 8:28 AM, raheleh m <span dir="ltr"><<a rel="nofollow" href="mailto:meshgolok@yahoo.com" target="_blank">meshgolok@yahoo.com</a>></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;"><div><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><div>
<span style="font-weight: bold;">Dear Mr. Ibanez</span><br>
<br>Thanks a lot for your attention and help. Actually, I have another question and I would be really grateful<br>if you could please help me with this one too. I am using Multi Resolution Registration method of itk, and my data size is 512*512*200. The prepare pyramids function takes about 170 seconds, but the optimization loop is fast, about 12 seconds. The total registration time is about 202 seconds. Is this normal, or I miss something? How can I reduce the time of prepare pyramids? I know it is a lot to ask and I appreciate it a lot. <br>
<br><span style="font-style: italic;">Thanks in advance, </span><br>Sincerely Yours<br><br>Rahil Meshgolok<br><font face="Tahoma" size="2"> </font></div><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;">
<br><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><font face="Tahoma" size="2"><hr size="1"><b><span style="font-weight: bold;">From:</span></b> Luis Ibanez <<a rel="nofollow" href="mailto:luis.ibanez@kitware.com" target="_blank">luis.ibanez@kitware.com</a>><br>
<b><span style="font-weight: bold;">To:</span></b> raheleh m <<a rel="nofollow" href="mailto:meshgolok@yahoo.com" target="_blank">meshgolok@yahoo.com</a>><br><b><span style="font-weight: bold;">Cc:</span></b> itk <<a rel="nofollow" href="mailto:Insight-users@itk.org" target="_blank">Insight-users@itk.org</a>><br>
<b><span style="font-weight: bold;">Sent:</span></b> Thu, June 17, 2010 8:40:41 PM<br><b><span style="font-weight: bold;">Subject:</span></b> Re: Permute Order in MultiResMIRegistration<br></font><br>
<br>Hi Rahil,<br><br><br>Please read the documentation of the class at:<br><span><span><a href="http://public.kitware.com/Insight/Doxygen/html/classitk_1_1FlipImageFilter.html" target="_blank">http://public.kitware.com/Insight/Doxygen/html/classitk_1_1FlipImageFilter.html</a></span></span><br>
<br> "The image is flipped across axes for which array[i] is true."<br><br><span><span><a href="http://public.kitware.com/Insight/Doxygen/html/classitk_1_1PermuteAxesImageFilter.html" target="_blank">http://public.kitware.com/Insight/Doxygen/html/classitk_1_1PermuteAxesImageFilter.html</a></span></span><br>
<br> "The elements of the array must be a rearrangement <br> of the numbers from 0
to ImageDimension - 1.
<p> "The i-th axis of the output image corresponds with <br> the order[i]-th
axis of the input image."</p><br>If you pass [0 1 2] as parameter, then the output image will be identical<br>to the input image. If you pass [1 0 2] then the output image will use as<br>X axis, what used to be the Y axis of the input image, and will use as <br>
Y axis, what used to be the X axis of the input image.<br><br><br>Note that in general, you should NOT use the flip / permute<br>filters. Particularly not if you are using medical data.<br><br><br> Regards,<br><br><br>
Luis<br><br><br>-------------------------------------------------------------------------<br><div class="gmail_quote">On Mon, Jun 14, 2010 at 8:12 AM, raheleh m <span dir="ltr"><<a rel="nofollow" href="mailto:meshgolok@yahoo.com" target="_blank">meshgolok@yahoo.com</a>></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;"><div><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><div>
<i>Greetings!<br>
</i></div>
<div>I am using ITK Registration. I would be really grateful if you could please let me know the meaning of permute order and flip axes arrays. I mean when permute order is [0 1 2] what is the order of axes? what about [2 0 1]?</div>
<div>Does the first array mean (x,y,z) and the second (z,x,y)? I have the same question about flip axes.</div>
<div> </div>
<div>Thanks in advance</div>
<div>Rahil Meshgolok</div></div><br>
</div></blockquote></div><br>
</div></div>
</div><br>
</div></blockquote></div><br>
</div></div></div></div>
</div><br>
</div></blockquote></div><br>