Research Article
On the Usage of GPUs for Efficient Motion Estimation in Medical Image Sequences
Algorithm 1
The estimation algorithm without the volume preserving term (see Section 3.3. The meanings of the symbols remain as in the text. represents a static image and represents a moving image . First, the Image is initialized and partial derivatives are computed. denotes the partial derivative operator. Following this, the algorithm is applied repeatedly times. Each time of the iteration, as a precursor to solve (6), values of and are computed using Jacobi method. The routine Jacobi performs this operation. One this is available, all the unknowns, given by , are solved iteratively using the routine eqnSolve, which solves a system of linear equations. Then the solutions are used to resample the image to estimate the moving image .
| ||||||||||||||||||||||||||||