<div class="gmail_quote"><font color="#888888">Sorry. I forgot to change the subject of my previous email. Here I re-send it again.<br><br>====================================================================<br><br></font></div>
Dear ITK Users and Developers,<br><br>        In deformable image
registration, it is very common to find the image gradient at non-grid
positions. For example, to obtain the driving force by taking the
derivative of SSD similarity measure. I am wondering how this procedure
is implemented in ITK. I used to construct a continuous image based on
Bspline model and then find the derivitive at non-grid positions.
However, I recently found a MATLAB code which resamples the image based
on the current deformation field and then find the gradient at the
resampled image. In this way, there is no need to find the gradient at
non-grid positions. Afterwards, the gradient information is used to
construct the update of the deformation field and the new deformation
field is obtained by ADDING the update to the original deformation
field. Is this a common practice to implement deformable image
registration algorithms like demons, fluid, and FFD in ITK? Any
response is appreciated.<br><font color="#888888">
<br>Chen      </font>