[Insight-users] Problem implementing GradientAnisotropicDiffusionImageFilter

[Bill Lorensen]

In itk, the K divisor is average_gradient_squared *
conductance_squared. I suspect that accounts for some of the
difference.

On Tue, Dec 1, 2009 at 9:01 PM, O'Connor, Michael
<Michael.OConnor at umassmed.edu> wrote:
> Hello Luis,
>      The pixel spacing is .2 mm.   I did not SetSpacing in ITK code.
> Thanks for the insight.   I will check that out.
>
> Regards,
>    Michael
>
>
> -----Original Message-----
> From: Luis Ibanez [mailto:luis.ibanez at kitware.com]
> Sent: Tue 12/1/2009 7:56 PM
> To: O'Connor, Michael
> Cc: Bill Lorensen; insight-users at itk.org
> Subject: Re: [Insight-users] Problem implementing
> GradientAnisotropicDiffusionImageFilter
>
> Hi Michael,
>
> What is the pixel spacing of your image ?
>
> I believe that the ITK implementation may
> be taking the pixel spacing into account
> when using the value of conductance.
>
>
>    Please let us know,
>
>
>              Thanks
>
>
>                        Luis
>
>
> ------------------------------------------------------------------
> On Tue, Dec 1, 2009 at 5:17 PM, O'Connor, Michael
> <Michael.OConnor at umassmed.edu> wrote:
>> Hi Bill,
>>     I use the exact same file with ITK (the one attached to the previous
>> email with MATLAB code).
>>
>> Regards,
>>    Michael
>>
>>
>> -----Original Message-----
>> From: Bill Lorensen [mailto:bill.lorensen at gmail.com]
>> Sent: Tue 12/1/2009 4:17 PM
>> To: O'Connor, Michael
>> Cc: insight-users at itk.org
>> Subject: Re: [Insight-users] Problem implementing
>> GradientAnisotropicDiffusionImageFilter
>>
>> Can you provide a png file or meta file, or whatever file you tried with
>> itk?
>>
>> On Tue, Dec 1, 2009 at 4:08 PM, O'Connor, Michael
>> <Michael.OConnor at umassmed.edu> wrote:
>>> Hi Bill,
>>>     Thanks.   Attached is simple MATLAB test code  + the MATLAB filter
>>> function used (Dr. Guy Gilboa's code).   I've also take the liberty of
>>> attaching a test image (specifically opened in MATLAB code).   The test
>>> image is 400x400 float with background = 0.024 + additive Gaussian Noise
>>> and
>>> foreground (circle) = 0.031 + additive Gaussian Noise.   As said
>>> previously
>>> MATLAB code seems to work as expected as shown in plot.   In ITK there is
>>> no
>>> similar intraregion filtering (i.e. filtered image ~ input )
>>>
>>> Regards,
>>>    Michael
>>>
>
>
>