[Insight-users] how the itk pipeline get the update fied for all the voxels--self answer and further question

[Baoyun Li]

Dear All:
Seems itkDenseFiniteDifferenceImageFilter iterate all the pixel to get the updated fied. But I am still not clear about the pipeline.

In the case, if I want to use itk to write fluid registration algorithm proposed by Christensen, what should I do.

Should I derive a class from itkPDEDeformationRegistrationFilter?

Thanks

Baoyun




________________________________
From: Baoyun Li <baoyun_li123 at yahoo.com>
To: luis.ibanez at kitware.com; insight-users at itk.org
Sent: Friday, June 5, 2009 3:58:49 PM
Subject: in itk Demons Registration---how the itk pipeline get the update fied for all the voxels


Dear Luis and All:

I am learning the how the itk registration pipeline works using Demons as example. After reading the source code, I found in DemonsRegistrationFunction::ComputerUpdate, calculate the update for one voxel, but I could not find where to use this function to iterate the whole image to get the update field over the whole image.

Can somebody teach me.

Baoyun

00175 template <class TFixedImage, class TMovingImage, class TDeformationField>
00176 typename DemonsRegistrationFunction<TFixedImage,TMovingImage,TDeformationField>
00177 ::PixelType
00178 DemonsRegistrationFunction<TFixedImage,TMovingImage,TDeformationField>
00179 ::ComputeUpdate(const NeighborhoodType &it, void * gd,
00180                 const FloatOffsetType& itkNotUsed(offset))



      
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.itk.org/pipermail/insight-users/attachments/20090605/8e0dc737/attachment-0001.htm>