Hi Dan and all, <br><br>We use here your nice <a href="http://www.midasjournal.org/browse/publication/213">contribution</a> for the centerline extraction, and we also had some trouble with centerlines near the border, with a crash and that kind of backtrace :<div>
<span style="background-color:rgb(255,255,221);font-family:Verdana,Arial,'Bitstream Vera Sans',Helvetica,sans-serif;font-size:13px"><br></span></div>1 - 0x00000001001a0d28: itk::InterpolateImageFunction<itk::Image<float, 3u>, double>::Evaluate(itk::Point<double, 3u> const ) const in libThEndosize.dylib<br>
2 - 0x0000000100194eb4: itk::PhysicalCentralDifferenceImageFunction<itk::Image<float, 3u>, double>::Evaluate(itk::Point<double, 3u> const ) const in libThEndosize.dylib<br>3 - 0x00000001001951c4: itk::SingleImageCostFunction<itk::Image<float, 3u> >::GetDerivative(itk::Array<double> const , itk::Array<double> ) const in libThEndosize.dylib<br>
4 - 0x0000000100620ae1: itk::RegularStepGradientDescentBaseOptimizer::ResumeOptimization() in libThItk.dylib<br>5 - 0x0000000100620f62: itk::RegularStepGradientDescentBaseOptimizer::StartOptimization() in libThItk.dylib<br>
6 - 0x00000001001ac692: itk::ArrivalFunctionToPathFilter<itk::Image<float, 3u>, itk::PolyLineParametricPath<3u> >::GenerateData() in libThEndosize.dylib<br>7 - 0x00000001001b0651: itk::SpeedFunctionToPathFilter<itk::Image<float, 3u>, itk::PolyLineParametricPath<3u> >::GenerateData() in libThEndosize.dylib<div>
<br>Our assumption was also that in some case when the centerline was near the border the gradient descent optimizer ask for a point outside of the image, leading to that crash, but we disdn't use the "padding" hack. So my question is : how much should be the padding to be sure that a crash like that won't occurs ?<br>
<span style="background-color:rgb(255,255,221);font-family:Verdana,Arial,'Bitstream Vera Sans',Helvetica,sans-serif;font-size:13px"><br></span>Thanks again for your contribution.</div><div><br></div><div>-Simon<span style="background-color:rgb(255,255,221);font-family:Verdana,Arial,'Bitstream Vera Sans',Helvetica,sans-serif;font-size:13px"><br>
</span><br><div class="gmail_quote">On Sun, Apr 29, 2012 at 00:28, Dan Mueller <span dir="ltr"><<a href="mailto:dan.muel@gmail.com" target="_blank">dan.muel@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Charu,<br>
<br>
(I was just re-reading my email, I don't think I was clear enough.)<br>
<br>
I used Fast Marching to extract the paths in your 3d dataset. Then I<br>
projected both the dataset and the paths onto a 2d image, which was<br>
attached to my previous email. In short: the method works for 3d data.<br>
<br>
Cheers, Dan<br>
<div class="HOEnZb"><div class="h5"><br>
On 29 April 2012 07:54, Dan Mueller <<a href="mailto:dan.muel@gmail.com">dan.muel@gmail.com</a>> wrote:<br>
> Hi Charu,<br>
><br>
> The image I attached to my previous email was a MIP of your 3d<br>
> dataset. This method works with both 2d or 3d data.<br>
><br>
> Cheers, Dan<br>
><br>
> On 29 April 2012 05:17, <a href="mailto:chans@cs.uh.edu">chans@cs.uh.edu</a> <<a href="mailto:chans@cs.uh.edu">chans@cs.uh.edu</a>> wrote:<br>
>> Hi Dan,<br>
>><br>
>> Thank You for your reply. I looked at Fast Marching minimization. It works<br>
>> for 2d images. What, I have here is 3d data.<br>
>><br>
>> Thank You<br>
>> Charu<br>
>><br>
>> --<br>
>> View this message in context: <a href="http://itk-insight-users.2283740.n2.nabble.com/Centerline-tp7507501p7509617.html" target="_blank">http://itk-insight-users.2283740.n2.nabble.com/Centerline-tp7507501p7509617.html</a><br>
>> Sent from the ITK Insight Users mailing list archive at Nabble.com.<br>
>> _____________________________________<br>
>> Powered by <a href="http://www.kitware.com" 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" target="_blank">http://www.kitware.com/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" target="_blank">http://www.kitware.com/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" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
>><br>
>> Follow this link to subscribe/unsubscribe:<br>
>> <a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
_____________________________________<br>
Powered by <a href="http://www.kitware.com" 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" target="_blank">http://www.kitware.com/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" target="_blank">http://www.kitware.com/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" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
</div></div></blockquote></div><br><br clear="all"><div>-- <br>------------------------------------------------------------------<br>Simon Esneault - Therenva<br>Centre d'Innovation Technologique<br>Centre Cardio-Pneumologique<br>
CHU Pontchaillou<br>Rennes, France<br>Tel : +33 (0)6 64 61 30 94<br>Mail : <a href="mailto:simon.esneault@therenva.com" target="_blank">simon.esneault@therenva.com</a><br>------------------------------------------------------------------</div>
<br>
</div>