Registration of Cervical MRI Using Multifeature Mutual Information

• Published in: IEEE transactions on medical imaging 2009-09, Vol.28 (9), p.1412-1421
• Main Authors: Staring, M., van der Heide, U.A., Klein, S., Viergever, M.A., Pluim, J.
• Format: Article
• Language: English
• Subjects: Algorithms; alpha -mutual information; Biomedical applications of radiation; Biomedical imaging; Bladder; Cervical cancer; Cervix Uteri - anatomy & histology; Female; Humans; Image Processing, Computer-Assisted - methods; Image registration; k{\rm NN} graphs; local image structure; Magnetic resonance imaging; Magnetic Resonance Imaging - methods; Medical treatment; Mutual information; Neoplasms; nonrigid registration; Principal Component Analysis; Principal components analysis; radiation therapy; Rectum - anatomy & histology; Shannon mutual information; Stochastic processes; Urinary Bladder - anatomy & histology; Uterine Cervical Neoplasms - pathology; Uterus - anatomy & histology
• ISSN: 0278-0062; 1558-254X
• DOI: 10.1109/TMI.2009.2016560

Radiation therapy for cervical cancer can benefit from image registration in several ways, for example by studying the motion of organs, or by (partially) automating the delineation of the target volume and other structures of interest. In this paper, the registration of cervical data is addressed using mutual information (MI) of not only image intensity, but also features that describe local image structure. Three aspects of the registration are addressed to make this approach feasible. First, instead of relying on a histogram-based estimation of mutual information, which poses problems for a larger number of features, a graph-based implementation of alpha-mutual information (alpha-MI) is employed. Second, the analytical derivative of alpha-MI is derived. This makes it possible to use a stochastic gradient descent method to solve the registration problem, which is substantially faster than non-derivative-based methods. Third, the feature space is reduced by means of a principal component analysis, which also decreases the registration time. The proposed technique is compared to a standard approach, based on the mutual information of image intensity only. Experiments are performed on 93 T2-weighted MR clinical data sets acquired from 19 patients with cervical cancer. Several characteristics of the proposed algorithm are studied on a subset of 19 image pairs (one pair per patient). On the remaining data (36 image pairs, one or two pairs per patient) the median overlap is shown to improve significantly compared to standard MI from 0.85 to 0.86 for the clinical target volume (CTV, p = 2 ldr10 -2 ), from 0.75 to 0.81 for the bladder ( p = 8 ldr 10 -6 ), and from 0.76 to 0.77 for the rectum ( p = 2 ldr 10 -4 ). The registration error is improved at important tissue interfaces, such as that of the bladder with the CTV, and the interface of the rectum with the uterus and cervix.