Non-diffeomorphic registration of brain tumor images by simulating tissue loss and tumor growth

Although a variety of diffeomorphic deformable registration methods exist in the literature, application of these methods in the presence of space-occupying lesions is not straightforward. The motivation of this work is spatial normalization of MR images from patients with brain tumors in a common stereotaxic space, aiming to pool data from different patients into a common space in order to perform group analyses. Additionally, transfer of structural and functional information from neuroanatomical brain atlases into the individual patient's space can be achieved via the inverse mapping, for the purpose of segmenting brains and facilitating surgical or radiotherapy treatment planning. A method that estimates the brain tissue loss and replacement by tumor is applied for achieving equivalent image content between an atlas and a patient's scan, based on a biomechanical model of tumor growth. Automated estimation of the parameters modeling brain tissue loss and displacement is performed via optimization of an objective function reflecting feature-based similarity and elastic stretching energy, which is optimized in parallel via APPSPACK (Asynchronous Parallel Pattern Search). The results of the method, applied to 21 brain tumor patients, indicate that the registration accuracy is relatively high in areas around the tumor, as well as in the healthy portion of the brain. Also, the calculated deformation in the vicinity of the tumor is shown to correlate highly with expert-defined visual scores indicating the tumor mass effect, thereby potentially leading to an objective approach to quantification of mass effect, which is commonly used in diagnosis.