<html>
<body>
<br><br>
Hi all,<br><br>
I'm trying to modify the <b>MultiResolutionImageRegistration Method
</b>in order to improve the results of a global algorithm containing
it.<br>
What I wanted to do is to <b>replace the DiscreteGaussianFiltering of the
smoothing step of every level in the pyramid by the
AnisotropicDiffusionImageFiltering</b>, so that the sharp boundaries in
the image that help to distinguish the anatomical structures would be
preserved...<br><br>
For doing so, I have tried to modify:<br><br>
itkMYMultiResolutionImageRegistrationMethod.txx<x-tab> </x-tab><br>
itkMYMultiResolutionPyramidImageFilter.h <br><br>
itkMYRecursiveMultiResolutionPyramidImageFilter.txx <br>
itkMYMultiResolutionImageRegistrationMethod.h<br><br>
itkMYMultiResolutionPyramidImageFilter.txx <br>
itkMYRecursiveMultiResolutionPyramidImageFilter.h<br><br>
The problem is that I could change the smoothing filter both in <br>
void MYRecursiveMultiResolutionPyramidImageFilter<TInputImage,
TOutputImage> ::GenerateData() <br>
<x-tab> </x-tab>(in
itkMYRecursiveMultiResolutionPyramidImageFilter.txx )<br>
void MYMultiResolutionPyramidImageFilter<TInputImage, TOutputImage>
::GenerateData()<br>
<x-tab> </x-tab>(in
itkMYMultiResolutionPyramidImageFilter.txx )<br><br>
so that I included properly the <br>
<x-tab> </x-tab>typedef
GradientAnisotropicDiffusionImageFilter<TOutputImage, TOutputImage>
SmootherType;<br>
and changed the different parameters.<br><br>
But then I found that the methods<br>
void MYMultiResolutionPyramidImageFilter<TInputImage,
TOutputImage>::GenerateInputRequestedRegion()<br>
void
MYRecursiveMultiResolutionPyramidImageFilter<TInputImage,
TOutputImage>::GenerateOutputRequestedRegion(DataObject * ptr )<br>
void
MYRecursiveMultiResolutionPyramidImageFilter<TInputImage,
TOutputImage>::GenerateInputRequestedRegion()<br>
are using, <b>to compute requirements for the smoothing part, the kernel
implementation:<br>
typedef
GaussianOperator<OutputPixelType,ImageDimension>
OperatorType;<br><br>
How could I also replace this operator by an anisotropic one? </b>I
couldn't find that implementation...<br>
Maybe I'm not getting how does it work, and the kernel implemented
in the <b>GaussianOperator </b>for the smoothing step is also valid for
the anisotropic case or just I shouldn't need it, but I don't think
so...<br><br>
Would somebody help me?<br>
Thank you very much in advance.<br>
Rosario<br><br>
<br>
<x-sigsep><p></x-sigsep>
<font face="Arial TUR, Helvetica">Rosario Sance Garzón<br>
Dpto. Ingeniería Electrónica<br>
E.T.S.I. Telecomunicación<br>
Avda. Ciudad Universitaria s/n<br>
E-28040 - Madrid (SPAIN)<br><br>
</font><i>Teléfono:</i> +34 91 549 57 00 (Ext. 4220) <br>
<i>Fax: </i>+34<i> </i>91 336
73 23<br>
<font face="Eras Medium ITC">
___________________________________________________</font> <br>
</body>
</html>