Robust Cardiac Motion Estimation using Ultrafast Ultrasound Data: A Low-Rank-Topology-Preserving Approach

Cardiac motion estimation is an important diagnostic tool to detect heart diseases and it has been explored with modalities such as MRI and conventional ultrasound (US) sequences. US cardiac motion estimation still presents challenges because of the complex motion patterns and the presence of noise. In this work, we propose a novel approach to estimate the cardiac motion using ultrafast ultrasound data. -- Our solution is based on a variational formulation characterized by the L2-regularized class. The displacement is represented by a lattice of b-splines and we ensure robustness by applying a maximum likelihood type estimator. While this is an important part of our solution, the main highlight of this paper is to combine a low-rank data representation with topology preservation. Low-rank data representation (achieved by finding the k-dominant singular values of a Casorati Matrix arranged from the data sequence) speeds up the global solution and achieves noise reduction. On the other hand, topology preservation (achieved by monitoring the Jacobian determinant) allows to radically rule out distortions while carefully controlling the size of allowed expansions and contractions. Our variational approach is carried out on a realistic dataset as well as on a simulated one. We demonstrate how our proposed variational solution deals with complex deformations through careful numerical experiments. While maintaining the accuracy of the solution, the low-rank preprocessing is shown to speed up the convergence of the variational problem. Beyond cardiac motion estimation, our approach is promising for the analysis of other organs that experience motion.Comment: 15 pages, 10 figures, Physics in Medicine and Biology, 201

A Heterogeneous and Multi-Range Soft-Tissue Deformation Model for Applications in Adaptive Radiotherapy

During fractionated radiotherapy, anatomical changes result in uncertainties in the applied dose distribution. With increasing steepness of applied dose gradients, the relevance of patient deformations increases. Especially in proton therapy, small anatomical changes in the order of millimeters can result in large range uncertainties and therefore in substantial deviations from the planned dose. To quantify the anatomical changes, deformation models are required. With upcoming MR-guidance, the soft-tissue deformations gain visibility, but so far only few soft-tissue models meeting the requirements of high-precision radiotherapy exist. Most state-of-the-art models either lack anatomical detail or exhibit long computation times. In this work, a fast soft-tissue deformation model is developed which is capable of considering tissue properties of heterogeneous tissue. The model is based on the chainmail (CM)-concept, which is improved by three basic features. For the first time, rotational degrees of freedom are introduced into the CM-concept to improve the characteristic deformation behavior. A novel concept for handling multiple deformation initiators is developed to cope with global deformation input. And finally, a concept for handling various shapes of deformation input is proposed to provide a high flexibility concerning the design of deformation input. To demonstrate the model flexibility, it was coupled to a kinematic skeleton model for the head and neck region, which provides anatomically correct deformation input for the bones. For exemplary patient CTs, the combined model was shown to be capable of generating artificially deformed CT images with realistic appearance. This was achieved for small-range deformations in the order of interfractional deformations, as well as for large-range deformations like an arms-up to arms-down deformation, as can occur between images of different modalities. The deformation results showed a strong improvement in biofidelity, compared to the original chainmail-concept, as well as compared to clinically used image-based deformation methods. The computation times for the model are in the order of 30 min for single-threaded calculations, by simple code parallelization times in the order of 1 min can be achieved. Applications that require realistic forward deformations of CT images will benefit from the improved biofidelity of the developed model. Envisioned applications are the generation of plan libraries and virtual phantoms, as well as data augmentation for deep learning approaches. Due to the low computation times, the model is also well suited for image registration applications. In this context, it will contribute to an improved calculation of accumulated dose, as is required in high-precision adaptive radiotherapy

Dynamic Deformation and Mechanical Properties of Brain Tissue

Traumatic brain injury is an important medical problem affecting millions of people. Mathematical models of brain biomechanics are being developed to simulate the mechanics of brain injury and to design protective devices. However, because of a lack of quantitative data on brain-skull boundary conditions and deformations, the predictions of mathematical models remain uncertain. The objectives of this dissertation are to develop methods and obtain experimental data that will be used to parameterize and validate models of traumatic brain injury. To that end, this dissertation first addresses the brain-skull boundary conditions by measuring human brain motion using tagged magnetic resonance imaging. Magnetic resonance elastography was performed in the ferret brain to measure its mechanical properties in vivo. Brain tissue is not only heterogeneous, but may also be anisotropic. To characterize tissue anisotropy, an experimental procedure combining both shear testing and indentation was developed and applied to white matter and gray matter. These measurements of brain-skull interactions and mechanical properties of the brain will be valuable in the development and validation of finite element simulations of brain biomechanics

Analysis of Venous Blood Flow and Deformation in the Calf under External Compression

Deep vein thrombosis (DVT) is a common post-operative complication, and a serious threat to the patient’s general recovery. In recent years, there has been increasing awareness of the risk of DVT in healthy individuals after prolonged immobility, such as people taking long-period flights or sitting at a computer. Mechanical methods of DVT prophylaxis, such as compression stockings, have gained widespread acceptance, but the haemodynamic mechanism of their action is still not well understood. In this study, computational modelling approaches based on magnetic resonance (MR) images are used to (i) predict the deformation of calf and deep veins under external compression, (ii) determine blood flow and wall shear stress in the deep veins of the calf, and (iii) quantify the effect of external compression on flow and wall shear stress in the deep veins. As a first step, MR images of the calf obtained with and without external compression were analysed, which indicated different levels of compressibility for different calf muscle compartments. A 2D finite element model (FEM) with specifically tailored boundary conditions for different muscle components was developed to simulate the deformation of the calf under compression. The calf tissues were described by a linear elastic model. The simulation results showed a good qualitative agreement with the measurements in terms of deep vein deformation, but the area reduction predicted by the FEM was much larger than that obtained from the MR images. In an attempt to improve the 2D FEM, a hyperelastic material model was employed and a finite element based non-rigid registration algorithm was developed to calculate the bulk modulus of the calf tissues. Using subject-specific bulk modulus derived with this method together with a hyperelastic material model, the numerical results showed better quantitative agreement with MR measured deformations of deep veins and calf tissues. In order to understand the effect of external compression on flow in the deep veins, MR imaging and real-time flow mapping were performed on 10 healthy volunteers before and after compression. Computational fluid dynamics was then employed to calculate the haemodynamic wall shear stress (WSS), based on the measured changes in vessel geometry and flow waveforms. The overall results indicated that application of the compression stocking led to a reduction in both blood flow rate and cross sectional area of the peroneal veins in the calf, which resulted in an increase in WSS, but the individual effects were highly variable. Finally, a 3D fluid-structure interactions (FSI) model was developed for a segment of the calf with realistic geometry for the calf muscle and bones but idealised geometry for the deep vein. The hyperelastic material properties evaluated previously were employed to describe the solid behaviours. Some predictive ability of the FSI model was demonstrated, but further improvement and validation are still needed

Left Ventricle Myocardium Segmentation from 3D Cardiac MR Images using Combined Probabilistic Atlas and Graph Cut-based Approaches

Medical imaging modalities, including Computed Tomography (CT) Magnetic Resonance Imaging (MRI) and Ultrasound (US) are critical for the diagnosis and progress monitoring of many cardiac conditions, planning, visualization and delivery of therapy via minimally invasive intervention procedures, as well as for teaching, training and simulation applications. Image segmentation is a processing technique that allows the user to extract the necessary information from an image dataset, in the form of a surface model of the region of interest from the anatomy. A wide variety of segmentation techniques have been developed and implemented for cardiac MR images. Despite their complexity and performance, many of them are intended for specific image datasets or are too specific to be employed for segmenting classical clinical quality Magnetic Resonance (MR) images. Graph Cut based segmentation algorithms have been shown to work well in regards to medical image segmentation. In addition, they are computationally efficient, which scales well to real time applications. While the basic graph cuts algorithms use lower-order statistics, combining this segmentation approach with atlas-based methods may help improve segmentation accuracy at a lower computational cost. The proposed technique will be tested at each step during the development by assessing the segmentation results against the available ground truth segmentation. Several metrics will be used to quantify the performance of the proposed technique, including computational performance, segmentation accuracy and fidelity assessed via the Sørensen-Dice Coefficient (DSC), Mean Absolute Distance (MAD) and Hausdorff Distance (HD) metrics

Alignment of contrast enhanced medical images

The re-alignment of series of medical images in which there are multiple contrast variations is difficult. The reason for this is that the popularmeasures of image similarity used to drive the alignment procedure do not separate the influence of intensity variation due to image feature motion and intensity variation due to feature enhancement. In particular, the appearance of new structure poses problems when it has no representation in the original image. The acquisition of many images over time, such as in dynamic contrast enhanced MRI, requires that many images with different contrast be registered to the same coordinate system, compounding the problem. This thesis addresses these issues, beginning by presenting conditions under which conventional registration fails and proposing a solution in the form of a ’progressive principal component registration’. The algorithm uses a statistical analysis of a series of contrast varying images in order to reduce the influence of contrast-enhancement that would otherwise distort the calculation of the image similarity measures used in image registration. The algorithm is shown to be versatile in that it may be applied to series of images in which contrast variation is due to either temporal contrast enhancement changes, as in dynamic contrast-enhanced MRI or intrinsically in the image selection procedure as in diffusion weighted MRI

Biomechanical modelling of the pelvic system: improving the accuracy of the location of neoplasms in MRI-TRUS fusion prostate biopsy

Background An accurate knowledge of the relocation of prostate neoplasms during biopsy is of great importance to reduce the number of false negative results. Prostate neoplasms are visible in magnetic resonance images (MRI) but it is difficult for the practitioner to locate them at the time of performing a transrectal ultrasound (TRUS) guided biopsy. In this study, we present a new methodology, based on simulation, that predicts both prostate deformation and lesion migration during the biopsy. Methods A three-dimensional (3-D) anatomy model of the pelvic region, based on medical images, is constructed. A finite element (FE) numerical simulation of the organs motion and deformation as a result of the pressure exerted by the TRUS probe is carried out using the Code-Aster open-source computer software. Initial positions of potential prostate lesions prior to biopsy are taken into consideration and the final location of each lesion is targeted in the FE simulation output. Results Our 3-D FE simulations show that the effect of the pressure exerted by the TRUS probe is twofold as the prostate experiences both a motion and a deformation of its original shape. We targeted the relocation of five small prostate lesions when the TRUS probe exerts a force of 30 N on the rectum inner wall. The distance travelled by these lesions ranged between 5.6 and 13.9 mm. Conclusions Our new methodology can help to predict the location of neoplasms during a prostate biopsy but further studies are needed to validate our results. Moreover, the new methodology is completely developed on open-source software, which means that its implementation would be affordable to all healthcare providers