Hi,<br><br>I'm now using a smoothed distance map (output from SignedDanielssonDistanceMapImageFilter + RecursiveGaussianImageFilter) as an input to the AverageOutwardFluxFilter.<br><br>The inputs for MedialCurveImageFilter are the smoothed distance map and the output from the AverageOutwardFluxFilter.<br>
<br>The method still fails at the MedialCurveFilter (after the computation of the AverageOutwardFluxFilter).<br><br>I removed the comment from line 617 of itkMedialCurveImageFilter.txx (std::cout<<q<<" Is deleted"<<std::endl;) - some time after entering the MedialCurveFilter it gives me this message in a loop, deleting all the voxels of the image (?). I don't understand this.<br>
<br>Grateful for your help,<br><br><div class="gmail_quote">2011/5/25 Pedro Névoa <span dir="ltr"><<a href="mailto:pedronevoa@gmail.com">pedronevoa@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
OK, I'll do that.<br><br>Thank you,<div><div></div><div class="h5"><br><br><div class="gmail_quote">2011/5/25 Sergio Vera <span dir="ltr"><<a href="mailto:sergio.vera@alma3d.com" target="_blank">sergio.vera@alma3d.com</a>></span><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div>Hi</div>The combination that worked for me is to use a Smoothed version of the distance map (dist map + gaussian filter) however, check the following things: distance map is correctly applied (ie: values are bigger at the center of your mask).<div>
<br></div><div>Also check for any other wring input image. Input to the AverageOutwardFluxFilter should be the smoothed distance map. Input for the MedialSurfaceImageFilter should be the smoothed dmap + the output of the AverageOutwardFluxFilter...</div>
<div><br></div><div>If this fails, maybe having a look at the code might help. </div><div><br></div><div>Regards.<div><div></div><div><br><br><div class="gmail_quote">On Tue, May 24, 2011 at 7:01 PM, Pedro Névoa <span dir="ltr"><<a href="mailto:pedronevoa@gmail.com" target="_blank">pedronevoa@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div>Yes, the image contains a binary mask of the bone... What should I do then? <br><br>I used the DanielssonDistanceMapImageFilt<div>
er to obtain the distance map and tried to run MedialCurve with the output of this algorithm, but it still fails.<br>
<br>Thank you,</div><br><br></div><div class="gmail_quote">No dia 24 de Maio de 2011 17:31, Sergio Vera <span dir="ltr"><<a href="mailto:sergio.vera@alma3d.com" target="_blank">sergio.vera@alma3d.com</a>></span> escreveu:<div>
<div></div><div><br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div>Hi</div>Yes, vtk file, but what the image contains? it is a binary mask of the bone? If this is true, medialAxis will fail. <div><br></div><div>regards<div><div></div><div><br><br><div class="gmail_quote">
On Tue, May 24, 2011 at 5:49 PM, Pedro Névoa <span dir="ltr"><<a href="mailto:pedronevoa@gmail.com" target="_blank">pedronevoa@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Thank you for the answers.<br><br>I am using a .vtk image as input to the filter.<div><div>
</div><div><br><br>
<div class="gmail_quote">2011/5/24 Sergio Vera <span dir="ltr"><<a href="mailto:sergio.vera@alma3d.com" target="_blank">sergio.vera@alma3d.com</a>></span><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Pedro, <div>what image are you passing as input to the filter?</div><div><br></div><div>MedialCurveImageFilter requires a distance map in order to work properly. Even more, you need to supply a smoothed distance map in. Although not explicitly specified in the insight journal article, in the code example you can see how it is done. My experience is that unsmoothed distance map does not work, and that a gaussian blurred dist map is required.</div>
<div><br></div><div>Hope it helps<br><br><div class="gmail_quote"><div><div></div><div>On Tue, May 24, 2011 at 12:16 AM, Daniel Blezek <span dir="ltr"><<a href="mailto:Blezek.Daniel@mayo.edu" target="_blank">Blezek.Daniel@mayo.edu</a>></span> wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div><div></div><div>
<div>
<font face="Calibri, Verdana, Helvetica, Arial"><span style="font-size: 11pt;">Pedro,<br>
<br>
We’ve found the MedialCurve method to be very sensitive, and difficult to use correctly. In contrast, the BinaryThinningImageFilter3D (also from an Insight Journal article), does a very nice job.<br>
<br>
-dan<div><div></div><div><br>
<br>
<br>
On 5/23/11 3:26 PM, "Pedro Névoa" <<a href="http://pedronevoa@gmail.com" target="_blank">pedronevoa@gmail.com</a>> wrote:<br>
<br>
</div></div></span></font><div><div></div><div><blockquote><font face="Calibri, Verdana, Helvetica, Arial"><span style="font-size: 11pt;">Hello,<br>
<br>
I'm trying to reduce a femur to a simple line, in order to achieve an image of the femur's medial axis. I tried to use the MedialCurve method to accomplish this but no results were obtained.<br>
<br>
The pixel type of the input image is signed short. I realize that the input pixel type accepted by Medial Curve is float so I changed the pixel type of my image to float using ImageReadCastWrite but it kept on not giving any results - the algorithm is processed endlessly.<br>
<br>
Does anyone know what may be the problem?<br>
<br>
Grateful for your attention,<br>
</span></font></blockquote></div></div><font face="Calibri, Verdana, Helvetica, Arial"><span style="font-size: 11pt;"><br>
-- <br>
<b>Daniel Blezek, PhD<br>
</b>Medical Imaging Informatics Innovation Center<br>
<br>
P 127 or (77) 8 8886<br>
T 507 538 8886<br>
E <a href="http://blezek.daniel@mayo.edu" target="_blank">blezek.daniel@mayo.edu</a><br>
<br>
Mayo Clinic<br>
200 First St. S.W.<br>
Harwick SL-44<br>
Rochester, MN 55905<br>
<a href="http://mayoclinic.org" target="_blank">mayoclinic.org</a><br>
"It is more complicated than you think." -- RFC 1925<br>
</span></font>
</div>
<br></div></div>_____________________________________<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.html" target="_blank">http://www.kitware.com/products/protraining.html</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></blockquote></div><br><br clear="all"><br>-- <br>Sergio Vera<br><br> Alma IT Systems<br> C/ Vilana, 4B, 4º 1ª<br> 08022 Barcelona<br> T. (+34) 932 380 592<br> <a href="http://www.alma3d.com" target="_blank">www.alma3d.com</a><br>
</div>
</blockquote></div><br><br clear="all"><br></div></div>-- <br><font color="#888888">Pedro Névoa<br><br>
</font></blockquote></div><br><br clear="all"><br>-- <br>Sergio Vera<br><br> Alma IT Systems<br> C/ Vilana, 4B, 4º 1ª<br> 08022 Barcelona<br> T. (+34) 932 380 592<br> <a href="http://www.alma3d.com" target="_blank">www.alma3d.com</a><br>
</div></div></div>
</blockquote></div></div></div><br><br clear="all"><br>-- <br><font color="#888888">Pedro Névoa<br><br>
</font></blockquote></div><br><br clear="all"><br>-- <br>Sergio Vera<br><br> Alma IT Systems<br> C/ Vilana, 4B, 4º 1ª<br> 08022 Barcelona<br> T. (+34) 932 380 592<br> <a href="http://www.alma3d.com" target="_blank">www.alma3d.com</a><br>
</div></div></div>
</blockquote></div><br><br clear="all"><br></div></div>-- <br><font color="#888888">Pedro Névoa<br><br>
</font></blockquote></div><br><br clear="all"><br>-- <br>Pedro Névoa<br><br>