Refine
Document Type
- Article (10)
- In Proceedings (5)
- ZIB-Report (5)
- Research data (3)
- Doctoral Thesis (1)
- Master's Thesis (1)
- Poster (1)
Is part of the Bibliography
- no (26)
Keywords
- Sparse Geometry Reconstruction (2)
- Statistical Shape Models (2)
- 4D-CT (1)
- Deep Learning (1)
- Deep learning (1)
- Gruppenweise Registrierung (1)
- Knee Arthroplasty (1)
- Total Knee Arthoplasty (1)
- cartilage (1)
- cartilage morphometry (1)
Institute
In this study we investigate methods for fitting a Statistical Shape Model (SSM) to intraoperatively acquired point cloud data from a surgical navigation system. We validate the fitted models against the pre-operatively acquired Magnetic Resonance Imaging (MRI) data from the same patients.
We consider a cohort of 10 patients who underwent navigated total knee arthroplasty. As part of the surgical protocol the patients’ distal femurs were partially digitized. All patients had an MRI scan two months pre-operatively. The MRI data were manually segmented and the reconstructed bone surfaces used as ground truth against which the fit was compared. Two methods were used to fit the SSM to the data, based on (1) Iterative Closest Points (ICP) and (2) Gaussian Mixture Models (GMM).
For both approaches, the difference between model fit and ground truth surface averaged less than 1.7 mm and excellent correspondence with the distal femoral morphology can be demonstrated.
The reconstruction of an object's shape or surface from a set of 3D points is a common topic in materials and life sciences, computationally handled in computer graphics. Such points usually stem from optical or tactile 3D coordinate measuring equipment. Surface reconstruction also appears in medical image analysis, e.g. in anatomy reconstruction from tomographic measurements or the alignment of intra-operative navigation and preoperative planning data. In contrast to mere 3D point clouds, medical imaging yields contextual information on the 3D point data that can be used to adopt prior information on the shape that is to be reconstructed from the measurements. In this work we propose to use a statistical shape model (SSM) as a prior for surface reconstruction. The prior knowledge is represented by a point distribution model (PDM) that is associated with a surface mesh. Using the shape distribution that is modelled by the PDM, we reformulate the problem of surface reconstruction from a probabilistic perspective based on a Gaussian Mixture Model (GMM). In order to do so, the given measurements are interpreted as samples of the GMM. By using mixture components with anisotropic covariances that are oriented according to the surface normals at the PDM points, a surface-based tting is accomplished. By estimating the parameters of the GMM in a maximum a posteriori manner, the reconstruction of the surface from the given measurements is achieved. Extensive experiments suggest that our proposed approach leads to superior surface reconstructions compared to Iterative Closest Point (ICP) methods.
In this study we investigate methods for fitting a Statistical Shape Model (SSM) to intraoperatively acquired point cloud data from a surgical navigation system. We validate the fitted models against the pre-operatively acquired Magnetic Resonance Imaging (MRI) data from the same patients.
We consider a cohort of 10 patients who underwent navigated total knee arthroplasty. As part of the surgical protocol the patients’ distal femurs were partially digitized. All patients had an MRI scan two months pre-operatively. The MRI data were manually segmented and the reconstructed bone surfaces used as ground truth against which the fit was compared. Two methods were used to fit the SSM to the data, based on (1) Iterative Closest Points (ICP) and (2) Gaussian Mixture Models (GMM).
For both approaches, the difference between model fit and ground truth surface averaged less than 1.7 mm and excellent correspondence with the distal femoral morphology can be demonstrated.
The reconstruction of an object’s shape or surface from a set of 3D points plays an important role in medical image analysis, e.g. in anatomy reconstruction from tomographic measurements or in the process of aligning intra-operative navigation and preoperative planning data. In such scenarios, one usually has to deal with sparse data, which significantly aggravates the problem of reconstruction. However, medical applications often provide contextual information about the 3D point data that allow to incorporate prior knowledge about the shape that is to be reconstructed. To this end, we propose the use of a statistical shape model (SSM) as a prior for surface reconstruction. The SSM is represented by a point distribution model (PDM), which is associated with a surface mesh. Using the shape distribution that is modelled by the PDM, we formulate the problem of surface reconstruction from a probabilistic perspective based on a Gaussian Mixture Model (GMM). In order to do so, the given points are interpreted as samples of the GMM. By using mixture components with anisotropic covariances that are “oriented” according to the surface normals at the PDM points, a surface-based fitting is accomplished. Estimating the parameters of the GMM in a maximum a posteriori manner yields the reconstruction of the surface from the given data points. We compare our method to the extensively used Iterative Closest Points method on several different anatomical datasets/SSMs (brain, femur, tibia, hip, liver) and demonstrate superior accuracy and robustness on sparse data.
Gruppenweise Registrierung zur robusten Bewegungsfeldschätzung in artefaktbehafteten 4D-CT-Bilddaten
(2015)
Das Ziel der Strahlentherapie ist, eine möglichst hohe Dosis in den Tumor zu applizieren und zeitgleich die Strahlenexposition des Normalgewebes zu minimieren. Insbesondere bei thorakalen und abdominalen Tumoren treten aufgrund der Atmung während der Bestrahlung große, komplexe und patientenspezifisch unterschiedliche Bewegungen der Gewebe auf. Um den Einfluss dieser Bewegung auf die i.d.R. statisch geplante Dosisverteilung abzuschätzen, können unter Verwendung der nicht-linearen Bildregistrierung anhand von 3D-CT-Aufnahmen eines Atmungszyklus - also 4D-CT-Daten - zunächst die Bewegungsfelder für die strahlentherapeutisch relevanten Strukturen, beispielsweise für die Lunge, berechnet werden. Diese Informationen bilden die Grundlage für sogenannte 4D-Dosisberechnungs- oder Dosisakkumulationsverfahren. Deren Genauigkeit hängt aber wesentlich von der Genauigkeit der Bewegungsfeldschätzung ab.
Klassisch erfolgt die Berechnung der Bewegungsfelder mittels paarweiser Bildregistrierung, womit für die Berechnung des Bewegungsfeldes zwischen zwei Bildern im Allgemeinen eine sehr hohe Genauigkeit erreicht wird. Auch für CT-Bilder, die Bewegungsartefakte, wie beispielsweise doppelte oder unvollständige Strukturen, enthalten, wird unter Verwendung der paarweisen Bildregistrierung im Kontext der Registrierung eine exakte Abbildung der anatomischen Strukturen zwischen den beiden Bildern erreicht. Dabei erfolgt aber eine physiologisch unplausible Anpassung der Felder an die Artefakte. Bei Verwendung der paarweisen Bildregistrierung müssen weiterhin für einen Atemzyklus die Voxel-Trajektorien aus Bewegungsfeldern zwischen mehreren dreidimensionalen Bildern zusammengesetzt werden. Durch Bewegungsartefakte entsprechen diese Trajektorien dann teilweise keiner natürlichen Bewegung der anatomischen Strukturen. Diese Ungenauigkeit stellt in der klinischen Anwendung ein Problem dar; dies gilt umso mehr, wenn Bewegungsartefakte im Bereich eines Tumors vorliegen.
Im Gegensatz zu der paarweisen Registrierung kann mit der gruppenweisen Registrierung das Problem der durch Bewegungsartefakte hervorgerufenen ungenauen Abbildung der physiologischen Gegebenheiten dadurch reduziert werden, dass im Registrierungsprozess Bildinformationen aller Bilder, also in diesem Kontext der CT-Daten zu unterschiedlichen Atemphasen, gleichzeitig genutzt werden. Es kann bereits im Registrierungsprozess eine zeitliche Glattheit der Voxel-Trajektorien gefordert werden.
In dieser Arbeit wird eine Methode zur B-Spline-basierten zeitlich regularisierten gruppenweisen Registrierung entwickelt. Die Genauigkeit der entwickelten Methode wird für frei zugängliche klinische Datensätze landmarkenbasiert evaluiert. Dabei wird mit dem Target Registration Error (TRE) die durchschnittliche dreidimensionale euklidische Distanz zwischen den korrespondierenden Landmarken nach Transformation der Landmarken bezeichnet. Eine Genauigkeit in der Größenordnung von aktuellen paarweisen Registrierungen verdeutlicht die Qualität des vorgestellten Registrierungs-Algorithmus. Anschließend werden die Vorteile der gruppenweisen Registrierung durch Experimente an einem Lungenphantom und an manipulierten, artefaktbehafteten klinischen 4D-CT-Bilddaten demonstriert. Dabei werden unter Verwendung der gruppenweisen Registrierung im Vergleich zu der paarweisen Registrierung glattere Trajektorien berechnet, die der realen Bewegung der anatomischen Strukturen stärker entsprechen. Für die Patientendaten wird außerdem anhand von automatisch detektierten Landmarken der TRE ausgewertet. Der TRE verschlechterte sich für die paarweise Bildregistrierung unter Vorliegen von Bewegungsartefakten von durchschnittlich 1,30 mm auf 3,94 mm. Auch hier zeigte sich für die gruppenweise Registrierung die Robustheit gegenüber Bewegungsartefakten und der TRE verschlechterte sich nur geringfügig von 1,45 mm auf 1,71 mm.
In der Strahlentherapie von Lungentumoren kann mittels Dosisakkumulation der Einfluss von Atembewegungen auf statisch geplante Dosisverteilungen abgeschätzt werden. Grundlage sind 4D-CT-Daten des Patienten, aus denen mittels nicht-linearer Bildregistrierung eine Sequenz von Bewegungsfeldern berechnet wird. Typischerweise werden Methoden der paarweisen Bildregistrierung eingesetzt, d.h. konsekutiv zwei Atemphasen aufeinander registriert. Hierbei erfolgt i.d.R. eine physiologisch nicht plausible Anpassung der Felder an CT-Bewegungsartefakte. Gruppenweise Registrierungsansätze berücksichtigen hingegen gleichzeitig sämtliche Bilddaten des 4D-CT-Scans und ermöglichen die Integration von zeitlichen Konsistenzbetrachtungen. In diesem Beitrag wird der potentielle Vorteil der gruppen- im Vergleich zur paarweisen Registrierung in artefaktbehafteten 4D-CT-Daten untersucht.
Precise voxel trajectory estimation in 4D CT images is a prerequisite for reliable dose accumulation during 4D treatment planning. 4D CT image data is, however, often affected by motion artifacts and applying standard pairwise registration to such data sets bears the risk of aligning anatomical structures to artifacts – with physiologically unrealistic trajectories being the consequence. In this work, the potential of a novel non-linear hybrid intensity- and feature-based groupwise registration method for robust motion field estimation in artifact-affected 4D CT image data is investigated. The overall registration performance is evaluated on the DIR-lab datasets; Its robustness if applied to artifact-affected data sets is analyzed using clinically acquired data sets with and without artifacts. The proposed registration approach achieves an accuracy comparable to the state-of-the-art (subvoxel accuracy), but smoother voxel trajectories compared to pairwise registration. Even more important: it maintained accuracy and trajectory smoothness in the presence of image artifacts – in contrast to standard pairwise registration, which yields higher landmark-based registration errors and a loss of trajectory smoothness when applied to artifact-affected data sets.
Volumetry of cartilage of the knee is needed for knee osteoarthritis (KOA) assessment. It is typically performed manually in a tedious and subjective process. We developed a method for an automated, segmentation-based quantification of cartilage volume by employing 3D Convolutional Neural Networks (CNNs). CNNs were trained in a supervised manner using magnetic resonance imaging data and cartilage volumetry readings performed by clinical experts for 1378 subjects provided by the Osteoarthritis Initiative. It was shown that 3D CNNs are able to achieve volume measures comparable to the magnitude of variation between expert readings and
the real in vivo situation. In the future, accurate automated cartilage volumetry might support both, diagnosis of KOA as well as longitudinal analysis of KOA progression.
We present a method based on a generative model for detection of disturbances such as prosthesis, screws, zippers, and metals in 2D radiographs. The generative model is trained in an unsupervised fashion using clinical radiographs as well as simulated data, none of which contain disturbances. Our approach employs a latent space consistency loss which has the benefit of identifying similarities, and is enforced to reconstruct X-rays without disturbances. In order to detect images with disturbances, an anomaly score is computed also employing the Frechet distance between the input X-ray and the reconstructed one using our generative model. Validation was performed using clinical pelvis radiographs. We achieved an AUC of 0.77 and 0.83 with clinical and synthetic data, respectively. The results demonstrated a good accuracy of our method for detecting outliers as well as the advantage of utilizing synthetic data.
We present a method for the quantification of knee alignment from full-leg X-Rays. A state-of-the-art object detector, YOLOv4, was trained to locate regions of interests (ROIs) in full-leg X-Ray images for the hip joint, the knee, and the ankle. Residual neural networks (ResNets) were trained to regress landmark coordinates for each ROI.Based on the detected landmarks the knee alignment, i.e., the hip-knee-ankle (HKA) angle, was computed. The accuracy of landmark detection was evaluated by a comparison to manually placed landmarks for 360 legs in 180 X-Rays. The accuracy of HKA angle computations was assessed on the basis of 2,943 X-Rays. Results of YARLA were compared to the results of two independent image reading studies(Cooke; Duryea) both publicly accessible via the Osteoarthritis Initiative. The agreement was evaluated using Spearman's Rho, and weighted kappa as well as regarding the correspondence of the class assignment (varus/neutral/valgus). The average difference between YARLA and manually placed landmarks was less than 2.0+- 1.5 mm for all structures (hip, knee, ankle). The average mismatch between HKA angle determinations of Cooke and Duryea was 0.09 +- 0.63°; YARLA resulted in a mismatch of 0.10 +- 0.74° compared to Cooke and of 0.18 +- 0.64° compared to Duryea. Cooke and Duryea agreed almost perfectly with respect to a weighted kappa value of 0.86, and showed an excellent reliability as measured by a Spearman's Rho value of 0.99. Similar values were achieved by YARLA, i.e., a weighted kappa value of0.83 and 0.87 and a Spearman's Rho value of 0.98 and 0.99 to Cooke and Duryea,respectively. Cooke and Duryea agreed in 92% of all class assignments and YARLA did so in 90% against Cooke and 92% against Duryea. In conclusion, YARLA achieved results comparable to those of human experts and thus provides a basis for an automated assessment of knee alignment in full-leg X-Rays.