Hi Luis,<br>     thanks for the quick reply. Nope, I do not believe I have enabled this flag. I am confused, I did take a look at the code snippet you suggested and figured that time step is &lt;= 1/2^(N+1) not 1/2^N (Judging from line 75 and/through 81). In reality I believe this has to be 1/2^N as you suggested, at least if the implementation is same as Perona and Malik.. <br>
<br>just in case my metadata looks as follows :<br><br>NDims = 3<br>DimSize = 256 256 256<br>ElementType = MET_SHORT<br>ElementDataFile = raw256.raw<br><br><br>Again, thanks for the help!<br><br>Cheers,<br><br>C.S.N<br><br>
<br><br><div class="gmail_quote">On Mon, Feb 15, 2010 at 2:00 PM, Luis Ibanez <span dir="ltr">&lt;<a href="mailto:luis.ibanez@kitware.com">luis.ibanez@kitware.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;">
Hi Natarajan,<br>
<br>
<br>
<br>
                      &quot;Glimpsing at the Source, leaves no doubt&quot;<br>
<br>
<br>
This is what the file:<br>
<br>
           Insight/Code/BasicFilters/itkAnisotropicDiffusionImageFilter.txx<br>
<br>
has in lines 75-82:<br>
<br>
<br>
  if ( m_TimeStep &gt;  (minSpacing / vcl_pow(2.0,<br>
static_cast&lt;double&gt;(ImageDimension) + 1))  )<br>
    {<br>
    //    f-&gt;SetTimeStep(1.0 / vcl_pow(2.0,<br>
static_cast&lt;double&gt;(ImageDimension)));<br>
    itkWarningMacro( &lt;&lt; &quot;Anisotropic diffusion unstable time step: &quot;<br>
                     &lt;&lt; m_TimeStep &lt;&lt; std::endl<br>
<div class="im">                     &lt;&lt; &quot;Stable time step for this image must be smaller than &quot;<br>
</div>                     &lt;&lt; minSpacing / vcl_pow(2.0,<br>
static_cast&lt;double&gt;(ImageDimension+1)));<br>
    }<br>
<br>
<br>
In summary:<br>
<br>
                    Timestep    &lt;     pixelSpacing  /  ( 2^N )<br>
<br>
where N is the image dimension (e.g. 3 for 3D images, and 2 for 2D images).<br>
<br>
<br>
<br>
<br>
Did you enable the flag:<br>
<br>
               UseImageSpacing  ?<br>
<br>
If so, what is the pixel spacing of your image ?<br>
<br>
<br>
      Regards,<br>
<br>
<br>
<br>
            Luis<br>
<br>
<br>
-------------------------------------------<br>
<div><div></div><div class="h5">On Mon, Feb 15, 2010 at 12:27 PM, Natarajan CS &lt;<a href="mailto:csnataraj@gmail.com">csnataraj@gmail.com</a>&gt; wrote:<br>
&gt; Hello all,<br>
&gt;       Not sure if this is a bug or a implementation detail I am overlooking,<br>
&gt; I am tending towards the latter, so any info would be helpful !<br>
&gt;  I get  the following warning from the<br>
&gt; GradientAnisotropicDiffusionImageFilter (The filter is being applied to a<br>
&gt; 256^3 dataset of type short)<br>
&gt;<br>
&gt; WARNING: In<br>
&gt; /hpc/soft/natac0/itk/include/InsightToolkit/BasicFilters/itkAnisotropicDiffusionImageFilter.txx,<br>
&gt; line 78<br>
&gt; GradientAnisotropicDiffusionImageFilter (0xbc4120): Anisotropic diffusion<br>
&gt; unstable time step: 0.07<br>
&gt; Stable time step for this image must be smaller than 0.0625<br>
&gt;<br>
&gt; i was under the impression that the ITK implimentation was a explicit method<br>
&gt; with the time step requirement for 3D images is &lt;= 1/8 and &lt;= 1/4 for 2D. Am<br>
&gt; i completely off-track here?<br>
&gt;<br>
&gt; Appreciate any help!<br>
&gt;<br>
&gt; Thanks,<br>
&gt;<br>
&gt; C.S.N<br>
&gt;<br>
</div></div>&gt; _____________________________________<br>
&gt; Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
&gt;<br>
&gt; Visit other Kitware open-source projects at<br>
&gt; <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
&gt;<br>
&gt; Kitware offers ITK Training Courses, for more information visit:<br>
&gt; <a href="http://www.kitware.com/products/protraining.html" target="_blank">http://www.kitware.com/products/protraining.html</a><br>
&gt;<br>
&gt; Please keep messages on-topic and check the ITK FAQ at:<br>
&gt; <a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
&gt;<br>
&gt; Follow this link to subscribe/unsubscribe:<br>
&gt; <a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
&gt;<br>
&gt;<br>
</blockquote></div><br>