Par0005 - elastix

Registration Description

intrapatient; rigid + B-spline transformation; multi-feature α-mutual information

Image data

Screen shot:

alt-text

Application

Nineteen patients were scanned. Each patient was scanned 5 times, one scan each week, except for two patients who were scanned 4 times. The data were divided into two sets. The first set contains the images at week 1 and 2 (19 image pairs), and is used for selection of the parameters. The second set consists of the remaining 55 images: 19 pairs of weeks 3 and 4, and 17 pairs of weeks 4 and 5; a total of 36 image pairs. The second set is used for comparing MI and α-MI with parameters tuned on the first set.

Registrations were performed to compare a multi-feature graph-based alpha-mutual information measure with a standard intensity-and-histogram-based mutual information measure.

The clinical purpose of registering these scans are in the field of radiotherapy and are multiple: to create a radiation dose targeting plan by means of atlas-based segmentation, to detect when this plan needs adaptation, and to define treatment margins. See the introduction of the paper for more details.

From the vast amount of choices, features that describe the local structure of images supply supplementary knowledge, which may improve the registration. We choose the set of Cartesian image structure features up to the second order derivatives, listed in Table I of the paper. These feature images are used as input for the alpha-mutual information metric, and are computed beforehand.

Registration Settings

For parameter files see the Elastix Model Zoo repository on GitHub.

elastix version: 4.003

Command line call:

elastix -f0 patient<j>_week<k>_ffeature0.mhd ... -fP patient<j>_week<k>_ffeatureP.mhd
        -m0 patient<j>_week<k>_mfeature0.mhd ... -mP patient<j>_week<k>_mfeatureP.mhd
        -p par0005.MI.rigid.txt -p par0005<X>.txt -out outputdir<k,l>

with: P the number of features
= one of the parameter files (see GitHub repository).

Published in

These registrations are described in the publication:

M. Staring, U.A. van der Heide, S. Klein, M.A. Viergever and J.P.W. Pluim, Registration of Cervical MRI Using Multifeature Mutual Information, IEEE Transactions on Medical Imaging, vol. 28, no. 9, pp. 1412 - 1421, September 2009.

Other comments

The features were computed outside of elastix and before the start of the registration. The source code of the feature computation software is made available through:

https://github.com/ITKTools/ITKTools

notably the tool pxgaussianimagefilter. The PCA analysis was performed using pxpca.

© 2020 Viktor van der Valk