Search

We propose a fully automated methodology for hexahedral meshing of patient-specific structures of the human knee obtained from magnetic resonance images, i.e. femoral/tibial cartilages and menisci. We select eight patients from the Osteoarthritis Initiative and validate our methodology using MATLAB on a laptop computer. We obtain the patient-specific meshes in an average of three minutes, while faithfully representing the geometries with well-shaped elements. We hope to provide a fundamentally different means to test hypotheses on the mechanisms of disease progression by integrating our patient-specific FE meshes with data from individual patients. Download both our meshes and software at http://im.engr.uconn.edu/downloads.php.

Accurate estimation of peak wall stress (PWS) is the crux of biomechanically motivated rupture risk assessment for abdominal aortic aneurysms aimed to improve clinical outcomes. Such assessments often use the finite element (FE) method to obtain PWS, albeit at a high computational cost, motivating simplifications in material or element formulations. These simplifications, while useful, come at a cost of reliability and accuracy. We achieve research-standard accuracy and maintain clinically applicable speeds by using novel computational technologies. We present a solution using our custom finite element code based on graphics processing unit (GPU) technology that is able to account for added complexities involved with more physiologically relevant solutions, e.g. strong anisotropy and heterogeneity. We present solutions up to 17× faster relative to an established finite element code using state-of-the-art nonlinear, anisotropic and nearly-incompressible material descriptions. We show a realistic assessment of the explicit GPU FE approach by using complex problem geometry, biofidelic material law, double-precision floating point computation and full element integration. Due to the increased solution speed without loss of accuracy, shown on five clinical cases of abdominal aortic aneurysms, the method shows promise for clinical use in determining rupture risk of abdominal aortic aneurysms. publisher: Elsevier articletitle: Rupture risk in abdominal aortic aneurysms: A realistic assessment of the explicit GPU approach journaltitle: Journal of Biomechanics articlelink: http://dx.doi.org/10.1016/j.jbiomech.2017.02.019 content_type: article copyright: © 2017 Elsevier Ltd. All rights reserved. ispartof: Journal of Biomechanics vol:56 pages:1-9 ispartof: location:United States status: published

Purpose: Current radiotherapy is progressing to the concept of adaptive radiotherapy, which implies the adaptation of planning along the treatment course. Nonrigid registration is an essential image processing tool for adaptive radiotherapy and image guided radiotherapy, and the three-dimensional (3D) nature of the current radiotherapy techniques requires a 3D quantification of the registration error that existing evaluation methods do not cover appropriately. The authors present a method for 3D evaluation of nonrigid registration algorithms’ performance, based on organ delineations, capable of working with near-spherical volumes even in the presence of concavities. Methods: The evaluation method is composed by a volume shape description stage, developed using a newad hoc volume reconstruction algorithm proposed by the authors, and an error quantification stage. The evaluation method is applied to the organ delineations of prostate and seminal vesicles, obtained by an automatic segmentation method over images of prostate cancer patients treated with intensity modulated radiation therapy. Results: The volume reconstruction algorithm proposed has been shown to accurately model complex 3D surfaces by the definition of clusters of control points. The quantification method, inspired by the Haussdorf–Chebysev distance, provides a measure of the largest registration error per control direction, defining a valid metric for concave-convex volumes. Summarizing, the proposed evaluation methodology presents accurate results with a high spatial resolution in a negligible computation time in comparison with the nonrigid registration time. Conclusions: Experimental results show that the metric selected for quantifying the registration error is of utmost importance in a quantitative evaluation based on measuring distances between volumes. The accuracy of the volume reconstruction algorithm is not so relevant as long as the reconstruction is tight enough on the actual volume of the organ. The new evaluation method provides a smooth and accurate volume reconstruction for both the reference and the registered organ, and a complete 3D description of nonrigid registration algorithms’ performance, resulting in a useful tool for study and comparison of registration algorithms for adaptive radiotherapy.

A correct patient-specific identification of the abdominal aortic aneurysm is useful for both diagnosis and treatment stages, as it locates the disease and represents its geometry. The actual thickness and shape of the arterial wall and the intraluminal thrombus are of great importance when predicting the rupture of the abdominal aortic aneurysms. The authors describe a novel method for delineating both the internal and external contours of the aortic wall, which allows distinguishing between vessel wall and intraluminal thrombus. The method is based on active shape model and texture statistical information. The method was validated with eight MR patient studies. There was high correspondence between automatic and manual measurements for the vessel wall area. Resulting segmented images presented a mean Dice coefficient with respect to manual segmentations of 0.88 and a mean modified Hausdorff distance of 1.14 mm for the internal face and 0.86 and 1.33 mm for the external face of the arterial wall. Preliminary results of the segmentation show high correspondence between automatic and manual measurements for the vessel wall and thrombus areas. However, since the dataset is small the conclusions cannot be generalized. ispartof: Computational and Mathematical Methods in Medicine vol:2015 ispartof: location:United States status: published

A novel method for generating patient-specific high quality conforming hexahedral meshes is presented. The meshes are directly obtained from the segmentation of patient magnetic resonance (MR) images of abdominal aortic aneu-rysms (AAA). The MRI permits distinguishing between struc-tures of interest in soft tissue. Being so, the contours of the lumen, the aortic wall and the intraluminal thrombus (ILT) are available and thus the meshes represent the actual anato-my of the patient?s aneurysm, including the layered morpholo-gies of these structures. Most AAAs are located in the lower part of the aorta and the upper section of the iliac arteries, where the inherent tortuosity of the anatomy and the presence of the ILT makes the generation of high-quality elements at the bifurcation is a challenging task. In this work we propose a novel approach for building quadrilateral meshes for each surface of the sectioned geometry, and generating conforming hexahedral meshes by combining the quadrilateral meshes. Conforming hexahedral meshes are created for the wall and the ILT. The resulting elements are evaluated on four patients? datasets using the Scaled Jacobian metric. Hexahedral meshes of 25,000 elements with 94.8% of elements well-suited for FE analysis are generated.

In order to perform finite element (FE) analyses of patient-specific abdominal aortic aneurysms, geometries derived from medical images must be meshed with suitable elements. We propose a semi-automatic method for generating conforming hexahedral meshes directly from contours segmented from medical images. Magnetic resonance images are generated using a protocol developed to give the abdominal aorta high contrast against the surrounding soft tissue. These data allow us to distinguish between the different structures of interest. We build novel quadrilateral meshes for each surface of the sectioned geometry and generate conforming hexahedral meshes by combining the quadrilateral meshes. The three-layered morphology of both the arterial wall and thrombus is incorporated using parameters determined from experiments. We demonstrate the quality of our patient-specific meshes using the element Scaled Jacobian. The method efficiently generates high-quality elements suitable for FE analysis, even in the bifurcation region of the aorta into the iliac arteries. For example, hexahedral meshes of up to 125,000 elements are generated in less than 130 s, with 94.8 % of elements well suited for FE analysis. We provide novel input for simulations by independently meshing both the arterial wall and intraluminal thrombus of the aneurysm, and their respective layered morphologies. QC 20140306

Two deformable registration methods, the Demons and the Morphon algorithms, have been used for registration of CT datasets to evaluate their usability in radiotherapy planning for prostate cancer. These methods were chosen because they can perform deformable registration in a fully automated way. The experiments show that for intrapatient registration both of the methods give useful results, although some differences exist in the way they deform the template. The Morphon method has, however, some advantageous compared to the Demons method. It is invariant to the image intensity and it does not distort the deformed data. The conclusion is therefore to recommend the Morphon method as a registration tool for this application. A more flexible regularization model is needed, though, in order to be able to catch the full range of deformations required to match the datasets.