[Insight-users] Deformation field smoothing

[Luis Ibanez]

Hi Miki,

Yes, this could be easily implemented.

Simply create a class that derives from the current
itkPDEDeformableRegistrationFilter and overload
the method:

template <class TFixedImage, class TMovingImage, class TDeformationField>
void
PDEDeformableRegistrationFilter<TFixedImage,TMovingImage,TDeformationField>
::SmoothDeformationField()


with one that implements the smoothing approach
that you want to use.


Note that you get the deformation field pointer with:

  DeformationFieldPointer field = this->GetOutput();


as shown in line 417 of  the file

   Insight/Code/Algorithms/
        itkPDEDeformableRegistrationFilter.txx


while you get access to the deformation update
by calling:

  DeformationFieldPointer field = this->GetUpdateBuffer();

as shown in line 490 of the same file.



     Regards,


           Luis


---------------------------------------------------------------
On Mon, Mar 29, 2010 at 4:55 AM, Miki Maric <mikimaric1975 at gmail.com> wrote:
> Hello to all,
>
> I'm interesting to change Gaussian smoothing to another similar approach in
>
> PDEDeformableRegistrationFilter<TFixedImage,TMovingImage,TDeformationField>
> ::SmoothUpdateField() and SmoothDeformationField()
>
> from the itkPDEDeformableRegistrationFilter.txx
>
> 1) Can this be easily implemented?
>
> 2) How to access raw values for deformation field in SmoothUpdateField()?
>
> Best regards,
>
> Miki
>
> _____________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> Kitware offers ITK Training Courses, for more information visit:
> http://www.kitware.com/products/protraining.html
>
> Please keep messages on-topic and check the ITK FAQ at:
> http://www.itk.org/Wiki/ITK_FAQ
>
> Follow this link to subscribe/unsubscribe:
> http://www.itk.org/mailman/listinfo/insight-users
>
>