<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">Hi Jonas,</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">what you probably need is <a href="https://itk.org/Doxygen/html/classitk_1_1ImageBase.html#af4a7c9c3787e9fdafbaaade2e02efa25">TransformPhysicalPointToIndex</a> method or one of its friends.</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">Regards,</div><div class="gmail_default"><font face="verdana, sans-serif">Dženan Zukić, PhD, Senior R&D Engineer, Kitware (Carrboro, N.C.)</font></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, May 27, 2017 at 3:34 AM, Jonas Teuwen <span dir="ltr"><<a href="mailto:jonasteuwen@gmail.com" target="_blank">jonasteuwen@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Hi all,</div><div><br></div>I want to extract patches around points using CropImageFilter. To do this, I need to calculate the correct bounds.<div>The point itself is an itk::Point object, and then I calculate the lower_crop and upper_crop. It can be possible these are negative in some components, which is checked later on and compensated with ConstantPadImageFilter.</div><div><br></div><div>So I need to switch between a Size and Point, which apparently cannot be implicitly casted. </div><div>I am using a for loop to do this, but this is probably not the right way to do it. So what is the proper way to do this in ITK?</div><div>I understand this question is probably very elementary, but your help would be greatly appreciated.</div><div><br></div><div>  itk::Size<3> image_size = image->GetBufferedRegion().<wbr>GetSize();<br></div><div><div><div>  std::cout << "Input size: " << image_size << std::endl;</div><div><br></div><div><b>  itk::Vector<double, 3> image_size2;</b></div><div><b>  for(int i=0; i<3; ++i)</b></div><div><b>    {</b></div><div><b>      image_size2[i] = image_size[i];</b></div><div><b>    }</b></div><div><br></div><div><br></div><div>  itk::Vector<double, 3> lower_crop = point - crop_size;</div><div>  itk::Vector<double, 3> upper_crop = image_size2 - (point + crop_size) - 1;</div></div><div><br></div></div><div><br></div><div>Best,</div><div>Jonas Teuwen</div></div>
<br>______________________________<wbr>_______<br>
Powered by <a href="http://www.kitware.com" rel="noreferrer" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at<br>
<a href="http://www.kitware.com/opensource/opensource.html" rel="noreferrer" target="_blank">http://www.kitware.com/<wbr>opensource/opensource.html</a><br>
<br>
Kitware offers ITK Training Courses, for more information visit:<br>
<a href="http://www.kitware.com/products/protraining.php" rel="noreferrer" target="_blank">http://www.kitware.com/<wbr>products/protraining.php</a><br>
<br>
Please keep messages on-topic and check the ITK FAQ at:<br>
<a href="http://www.itk.org/Wiki/ITK_FAQ" rel="noreferrer" target="_blank">http://www.itk.org/Wiki/ITK_<wbr>FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://public.kitware.com/mailman/listinfo/insight-users" rel="noreferrer" target="_blank">http://public.kitware.com/<wbr>mailman/listinfo/insight-users</a><br>
<br></blockquote></div><br></div>