TY - JOUR A1 - von Tycowicz, Christoph A1 - Ambellan, Felix A1 - Mukhopadhyay, Anirban A1 - Zachow, Stefan T1 - An Efficient Riemannian Statistical Shape Model using Differential Coordinates JF - Medical Image Analysis N2 - We propose a novel Riemannian framework for statistical analysis of shapes that is able to account for the nonlinearity in shape variation. By adopting a physical perspective, we introduce a differential representation that puts the local geometric variability into focus. We model these differential coordinates as elements of a Lie group thereby endowing our shape space with a non-Euclidean structure. A key advantage of our framework is that statistics in a manifold shape space becomes numerically tractable improving performance by several orders of magnitude over state-of-the-art. We show that our Riemannian model is well suited for the identification of intra-population variability as well as inter-population differences. In particular, we demonstrate the superiority of the proposed model in experiments on specificity and generalization ability. We further derive a statistical shape descriptor that outperforms the standard Euclidean approach in terms of shape-based classification of morphological disorders. Y1 - 2018 U6 - https://doi.org/10.1016/j.media.2017.09.004 VL - 43 IS - 1 SP - 1 EP - 9 ER - TY - JOUR A1 - Weiser, Martin A1 - Erdmann, Bodo A1 - Schenkl, Sebastian A1 - Muggenthaler, Holger A1 - Hubig, Michael A1 - Mall, Gita A1 - Zachow, Stefan T1 - Uncertainty in Temperature-Based Determination of Time of Death JF - Heat and Mass Transfer N2 - Temperature-based estimation of time of death (ToD) can be per- formed either with the help of simple phenomenological models of corpse cooling or with detailed mechanistic (thermodynamic) heat transfer mod- els. The latter are much more complex, but allow a higher accuracy of ToD estimation as in principle all relevant cooling mechanisms can be taken into account. The potentially higher accuracy depends on the accuracy of tissue and environmental parameters as well as on the geometric resolution. We in- vestigate the impact of parameter variations and geometry representation on the estimated ToD based on a highly detailed 3D corpse model, that has been segmented and geometrically reconstructed from a computed to- mography (CT) data set, differentiating various organs and tissue types. Y1 - 2018 U6 - https://doi.org/10.1007/s00231-018-2324-4 VL - 54 IS - 9 SP - 2815 EP - 2826 PB - Springer ER - TY - CHAP A1 - Ambellan, Felix A1 - Tack, Alexander A1 - Ehlke, Moritz A1 - Zachow, Stefan T1 - Automated Segmentation of Knee Bone and Cartilage combining Statistical Shape Knowledge and Convolutional Neural Networks: Data from the Osteoarthritis Initiative T2 - Medical Imaging with Deep Learning N2 - We present a method for the automated segmentation of knee bones and cartilage from magnetic resonance imaging, that combines a priori knowledge of anatomical shape with Convolutional Neural Networks (CNNs). The proposed approach incorporates 3D Statistical Shape Models (SSMs) as well as 2D and 3D CNNs to achieve a robust and accurate segmentation of even highly pathological knee structures. The method is evaluated on data of the MICCAI grand challenge "Segmentation of Knee Images 2010". For the first time an accuracy equivalent to the inter-observer variability of human readers has been achieved in this challenge. Moreover, the quality of the proposed method is thoroughly assessed using various measures for 507 manual segmentations of bone and cartilage, and 88 additional manual segmentations of cartilage. Our method yields sub-voxel accuracy. In conclusion, combining of anatomical knowledge using SSMs with localized classification via CNNs results in a state-of-the-art segmentation method. Y1 - 2018 ER - TY - JOUR A1 - Tack, Alexander A1 - Mukhopadhyay, Anirban A1 - Zachow, Stefan T1 - Knee Menisci Segmentation using Convolutional Neural Networks: Data from the Osteoarthritis Initiative JF - Osteoarthritis and Cartilage N2 - Abstract: Objective: To present a novel method for automated segmentation of knee menisci from MRIs. To evaluate quantitative meniscal biomarkers for osteoarthritis (OA) estimated thereof. Method: A segmentation method employing convolutional neural networks in combination with statistical shape models was developed. Accuracy was evaluated on 88 manual segmentations. Meniscal volume, tibial coverage, and meniscal extrusion were computed and tested for differences between groups of OA, joint space narrowing (JSN), and WOMAC pain. Correlation between computed meniscal extrusion and MOAKS experts' readings was evaluated for 600 subjects. Suitability of biomarkers for predicting incident radiographic OA from baseline to 24 months was tested on a group of 552 patients (184 incident OA, 386 controls) by performing conditional logistic regression. Results: Segmentation accuracy measured as Dice Similarity Coefficient was 83.8% for medial menisci (MM) and 88.9% for lateral menisci (LM) at baseline, and 83.1% and 88.3% at 12-month follow-up. Medial tibial coverage was significantly lower for arthritic cases compared to non-arthritic ones. Medial meniscal extrusion was significantly higher for arthritic knees. A moderate correlation between automatically computed medial meniscal extrusion and experts' readings was found (ρ=0.44). Mean medial meniscal extrusion was significantly greater for incident OA cases compared to controls (1.16±0.93 mm vs. 0.83±0.92 mm; p<0.05). Conclusion: Especially for medial menisci an excellent segmentation accuracy was achieved. Our meniscal biomarkers were validated by comparison to experts' readings as well as analysis of differences w.r.t groups of OA, JSN, and WOMAC pain. It was confirmed that medial meniscal extrusion is a predictor for incident OA. Y1 - 2018 U6 - https://doi.org/10.1016/j.joca.2018.02.907 VL - 26 IS - 5 SP - 680 EP - 688 ER - TY - GEN A1 - Tack, Alexander A1 - Mukhopadhyay, Anirban A1 - Zachow, Stefan T1 - Knee Menisci Segmentation using Convolutional Neural Networks: Data from the Osteoarthritis Initiative (Supplementary Material) N2 - Abstract: Objective: To present a novel method for automated segmentation of knee menisci from MRIs. To evaluate quantitative meniscal biomarkers for osteoarthritis (OA) estimated thereof. Method: A segmentation method employing convolutional neural networks in combination with statistical shape models was developed. Accuracy was evaluated on 88 manual segmentations. Meniscal volume, tibial coverage, and meniscal extrusion were computed and tested for differences between groups of OA, joint space narrowing (JSN), and WOMAC pain. Correlation between computed meniscal extrusion and MOAKS experts' readings was evaluated for 600 subjects. Suitability of biomarkers for predicting incident radiographic OA from baseline to 24 months was tested on a group of 552 patients (184 incident OA, 386 controls) by performing conditional logistic regression. Results: Segmentation accuracy measured as Dice Similarity Coefficient was 83.8% for medial menisci (MM) and 88.9% for lateral menisci (LM) at baseline, and 83.1% and 88.3% at 12-month follow-up. Medial tibial coverage was significantly lower for arthritic cases compared to non-arthritic ones. Medial meniscal extrusion was significantly higher for arthritic knees. A moderate correlation between automatically computed medial meniscal extrusion and experts' readings was found (ρ=0.44). Mean medial meniscal extrusion was significantly greater for incident OA cases compared to controls (1.16±0.93 mm vs. 0.83±0.92 mm; p<0.05). Conclusion: Especially for medial menisci an excellent segmentation accuracy was achieved. Our meniscal biomarkers were validated by comparison to experts' readings as well as analysis of differences w.r.t groups of OA, JSN, and WOMAC pain. It was confirmed that medial meniscal extrusion is a predictor for incident OA. Y1 - 2018 U6 - https://doi.org/10.12752/4.TMZ.1.0 N1 - Supplementary data to reproduce results from the related publication, including convolutional neural networks' weights. ER - TY - JOUR A1 - Hoffmann, Rene A1 - Lemanis, Robert A1 - Wulff, Lena A1 - Zachow, Stefan A1 - Lukeneder, Alexander A1 - Klug, Christian A1 - Keupp, Helmut T1 - Traumatic events in the life of the deep-sea cephalopod mollusc, the coleoid Spirula spirula JF - ScienceDirect: Deep Sea Research Part I - Oceanographic Research N2 - Here, we report on different types of shell pathologies of the enigmatic deep-sea (mesopelagic) cephalopod Spirula spirula. For the first time, we apply non-invasive imaging methods to: document trauma-induced changes in shell shapes, reconstruct the different causes and effects of these pathologies, unravel the etiology, and attempt to quantify the efficiency of the buoyancy apparatus. We have analysed 2D and 3D shell parameters from eleven shells collected as beach findings from the Canary Islands (Gran Canaria and Fuerteventura), West-Australia, and the Maldives. All shells were scanned with a nanotom-m computer tomograph. Seven shells were likely injured by predator attacks: fishes, cephalopods or crustaceans, one specimen was infested by an endoparasite (potentially Digenea) and one shell shows signs of inflammation and one shell shows large fluctuations of chamber volumes without any signs of pathology. These fluctuations are potential indicators of a stressed environment. Pathological shells represent the most deviant morphologies of a single species and can therefore be regarded as morphological end-members. The changes in the shell volume / chamber volume ratio were assessed in order to evaluate the functional tolerance of the buoyancy apparatus showing that these had little effect. Y1 - 2018 U6 - https://doi.org/10.1016/j.dsr.2018.10.007 VL - 142 IS - 12 SP - 127 EP - 144 ER - TY - GEN A1 - Tack, Alexander A1 - Mukhopadhyay, Anirban A1 - Zachow, Stefan T1 - Knee Menisci Segmentation using Convolutional Neural Networks: Data from the Osteoarthritis Initiative N2 - Abstract: Objective: To present a novel method for automated segmentation of knee menisci from MRIs. To evaluate quantitative meniscal biomarkers for osteoarthritis (OA) estimated thereof. Method: A segmentation method employing convolutional neural networks in combination with statistical shape models was developed. Accuracy was evaluated on 88 manual segmentations. Meniscal volume, tibial coverage, and meniscal extrusion were computed and tested for differences between groups of OA, joint space narrowing (JSN), and WOMAC pain. Correlation between computed meniscal extrusion and MOAKS experts' readings was evaluated for 600 subjects. Suitability of biomarkers for predicting incident radiographic OA from baseline to 24 months was tested on a group of 552 patients (184 incident OA, 386 controls) by performing conditional logistic regression. Results: Segmentation accuracy measured as Dice Similarity Coefficient was 83.8% for medial menisci (MM) and 88.9% for lateral menisci (LM) at baseline, and 83.1% and 88.3% at 12-month follow-up. Medial tibial coverage was significantly lower for arthritic cases compared to non-arthritic ones. Medial meniscal extrusion was significantly higher for arthritic knees. A moderate correlation between automatically computed medial meniscal extrusion and experts' readings was found (ρ=0.44). Mean medial meniscal extrusion was significantly greater for incident OA cases compared to controls (1.16±0.93 mm vs. 0.83±0.92 mm; p<0.05). Conclusion: Especially for medial menisci an excellent segmentation accuracy was achieved. Our meniscal biomarkers were validated by comparison to experts' readings as well as analysis of differences w.r.t groups of OA, JSN, and WOMAC pain. It was confirmed that medial meniscal extrusion is a predictor for incident OA. T3 - ZIB-Report - 18-15 Y1 - 2018 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-68038 SN - 1438-0064 VL - 26 IS - 5 SP - 680 EP - 688 ER - TY - GEN A1 - Grewe, Carl Martin A1 - Le Roux, Gabriel A1 - Pilz, Sven-Kristofer A1 - Zachow, Stefan T1 - Spotting the Details: The Various Facets of Facial Expressions N2 - 3D Morphable Models (MM) are a popular tool for analysis and synthesis of facial expressions. They represent plausible variations in facial shape and appearance within a low-dimensional parameter space. Fitted to a face scan, the model's parameters compactly encode its expression patterns. This expression code can be used, for instance, as a feature in automatic facial expression recognition. For accurate classification, an MM that can adequately represent the various characteristic facets and variants of each expression is necessary. Currently available MMs are limited in the diversity of expression patterns. We present a novel high-quality Facial Expression Morphable Model built from a large-scale face database as a tool for expression analysis and synthesis. Establishment of accurate dense correspondence, up to finest skin features, enables a detailed statistical analysis of facial expressions. Various characteristic shape patterns are identified for each expression. The results of our analysis give rise to a new facial expression code. We demonstrate the advantages of such a code for the automatic recognition of expressions, and compare the accuracy of our classifier to state-of-the-art. T3 - ZIB-Report - 18-06 Y1 - 2018 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-67696 SN - 1438-0064 ER - TY - CHAP A1 - Grewe, Carl Martin A1 - le Roux, Gabriel A1 - Pilz, Sven-Kristofer A1 - Zachow, Stefan T1 - Spotting the Details: The Various Facets of Facial Expressions T2 - IEEE International Conference on Automatic Face and Gesture Recognition Y1 - 2018 U6 - https://doi.org/10.1109/FG.2018.00049 SP - 286 EP - 293 ER - TY - CHAP A1 - Krämer, Martin A1 - Maggioni, Marta A1 - Tycowicz, Christoph von A1 - Brisson, Nick A1 - Zachow, Stefan A1 - Duda, Georg A1 - Reichenbach, Jürgen T1 - Ultra-short echo-time (UTE) imaging of the knee with curved surface reconstruction-based extraction of the patellar tendon T2 - ISMRM (International Society for Magnetic Resonance in Medicine), 26th Annual Meeting 2018, Paris, France N2 - Due to very short T2 relaxation times, imaging of tendons is typically performed using ultra-short echo-time (UTE) acquisition techniques. In this work, we combined an echo-train shifted multi-echo 3D UTE imaging sequence with a 3D curved surface reconstruction to virtually extract the patellar tendon from an acquired 3D UTE dataset. Based on the analysis of the acquired multi-echo data, a T2* relaxation time parameter map was calculated and interpolated to the curved surface of the patellar tendon. Y1 - 2018 ER - TY - CHAP A1 - Joachimsky, Robert A1 - Ambellan, Felix A1 - Zachow, Stefan T1 - Computerassistierte Auswahl und Platzierung von interpositionalen Spacern zur Behandlung früher Gonarthrose T2 - Proceedings of the Jahrestagung der Deutschen Gesellschaft für Computer- und Roboterassistierte Chirurgie (CURAC) N2 - Degenerative Gelenkerkrankungen, wie die Osteoarthrose, sind ein häufiges Krankheitsbild unter älteren Erwachsenen. Hierbei verringert sich u.a. der Gelenkspalt aufgrund degenerierten Knorpels oder geschädigter Menisci. Ein in den Gelenkspalt eingebrachter interpositionaler Spacer soll die mit der Osteoarthrose einhergehende verringerte Gelenkkontaktfläche erhöhen und so der teilweise oder vollständige Gelenkersatz hinausgezögert oder vermieden werden. In dieser Arbeit präsentieren wir eine Planungssoftware für die Auswahl und Positionierung eines interpositionalen Spacers am Patientenmodell. Auf einer MRT-basierten Bildsegmentierung aufbauend erfolgt eine geometrische Rekonstruktion der 3D-Anatomie des Kniegelenks. Anhand dieser wird der Gelenkspalt bestimmt, sowie ein Spacer ausgewählt und algorithmisch vorpositioniert. Die Positionierung des Spacers ist durch den Benutzer jederzeit interaktiv anpassbar. Für jede Positionierung eines Spacers wird ein Fitness-Wert zur Knieanatomie des jeweiligen Patienten berechnet und den Nutzern Rückmeldung hinsichtlich Passgenauigkeit gegeben. Die Software unterstützt somit als Entscheidungshilfe die behandelnden Ärzte bei der patientenspezifischen Spacerauswahl. KW - Osteoarthrose KW - Kniegelenk KW - Medizinische Planung KW - Visualisierung Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-65321 VL - 16 SP - 106 EP - 111 ER - TY - JOUR A1 - Hoffmann, René A1 - Lemanis, Robert A1 - Falkenberg, Janina A1 - Schneider, Steffen A1 - Wesendonk, Hendrik A1 - Zachow, Stefan T1 - Integrating 2D and 3D shell morphology to disentangle the palaeobiology of ammonoids: A virtual approach JF - Palaeontology N2 - We demonstrate, based on data derived from computed tomography, that integrating 2D and 3D morphological data of ammonoid shells represent an important new approach to disentangle the palaeobiology of ammonoids. Characterization of ammonite morphology has long been constrained to 2D data and only a few studies collect ontogenetic data in 180° steps. We combine this traditional approach with 3D data collected from high-resolution nano-computed tomography. For this approach, ontogenetic morphological data of the hollow shell of a juvenile ammonite Kosmoceras (Jurassic, Callovian) is collected. 2D data is collected in 10° steps and show significant changes in shell morphology. Preserved hollow spines show multiple mineralized membranes never reported before, representing temporal changes of the ammonoid mantle tissue. 3D data show that chamber volumes do not always increase exponentially, as was generally assumed, but may represent a proxy for life events, e.g. stress phases. Furthermore, chamber volume cannot be simply derived from septal spacing in forms comparable to Kosmoceras. Vogel numbers, a 3D parameter for chamber shape, of Kosmoceras are similar to other ammonoids (Arnsbergites, Amauroceras) and modern cephalopods (Nautilus, Spirula). Two methods to virtually document the suture line ontogeny, used to document phylogenetic relationships of larger taxonomic entities, were applied for the first time and seem a promising alternative to hand drawings. The curvature of the chamber surfaces increases during ontogeny due to increasing strength of ornamentation and septal complexity. As increasing curvature may allow for faster handling of cameral liquid, it could compensate for decreasing SA/V ratios through ontogeny. Y1 - 2017 UR - http://dx.doi.org/10.1111/pala.12328 U6 - https://doi.org/10.1111/pala.12328 VL - 61 IS - 1 SP - 89 EP - 104 ER - TY - JOUR A1 - Moldenhauer, Marian A1 - Weiser, Martin A1 - Zachow, Stefan T1 - Adaptive Algorithms for Optimal Hip Implant Positioning JF - PAMM N2 - In an aging society where the number of joint replacements rises, it is important to also increase the longevity of implants. In particular hip implants have a lifetime of at most 15 years. This derives primarily from pain due to implant migration, wear, inflammation, and dislocation, which is affected by the positioning of the implant during the surgery. Current joint replacement practice uses 2D software tools and relies on the experience of surgeons. Especially the 2D tools fail to take the patients’ natural range of motion as well as stress distribution in the 3D joint induced by different daily motions into account. Optimizing the hip joint implant position for all possible parametrized motions under the constraint of a contact problem is prohibitively expensive as there are too many motions and every position change demands a recalculation of the contact problem. For the reduction of the computational effort, we use adaptive refinement on the parameter domain coupled with the interpolation method of Kriging. A coarse initial grid is to be locally refined using goal-oriented error estimation, reducing locally high variances. This approach will be combined with multi-grid optimization such that numerical errors are reduced. Y1 - 2017 U6 - https://doi.org/10.1002/pamm.201710071 VL - 17 IS - 1 SP - 203 EP - 204 ER - TY - CHAP A1 - Ambellan, Felix A1 - Tack, Alexander A1 - Wilson, Dave A1 - Anglin, Carolyn A1 - Lamecker, Hans A1 - Zachow, Stefan T1 - Evaluating two methods for Geometry Reconstruction from Sparse Surgical Navigation Data T2 - Proceedings of the Jahrestagung der Deutschen Gesellschaft für Computer- und Roboterassistierte Chirurgie (CURAC) N2 - 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. KW - Total Knee Arthoplasty KW - Sparse Geometry Reconstruction KW - Statistical Shape Models Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-65339 VL - 16 SP - 24 EP - 30 ER - TY - JOUR A1 - Brüning, Jan A1 - Goubergrits, Leonid A1 - Heppt, Werner A1 - Zachow, Stefan A1 - Hildebrandt, Thomas T1 - Numerical Analysis of Nasal Breathing - A Pilot Study JF - Facial Plastic Surgery N2 - Background: Currently, there is no fully sufficient way to differentiate between symptomatic and normal nasal breathing. Using the nose’s total resistance is disputed as a valid means to objectify nasal airflow, and the need for a more comprehensive diagnostic method is increasing. This work’s aim was to test a novel approach considering intranasal wall shear stress as well as static pressure maps obtained by computational fluid dynamics (CFD). Methods: X-ray computed tomography (CT) scan data of six symptom-free subjects and seven symptomatic patients were used. Patient-specific geometries of the nasal cavity were segmented from these data sets. Inspiratory and expiratory steady airflow simulations were performed using CFD. Calculated static pressures and wall shear stresses (WSS) were mapped onto a common template of the nasal septum, allowing for comparison of these parameters between the two patient groups. Results: Significant differences in wall shear stress distributions during the inspiratory phase could be identified between the two groups, whereas no differences were found for the expiratory phase. It is assumed that one essential feature of normal nasal breathing probably consists in distinctively different intranasal flow fields for inspiration and expiration. This is in accordance with previous investigations. Conclusion: The proposed method seems to be a promising tool for developing a new kind of patient-specific assessment of nasal breathing. However, more studies and a greater case number of data with an expanded focus, would be ideal. Y1 - 2017 U6 - https://doi.org/doi:10.1055/s-0037-1603789 VL - 33 IS - 4 SP - 388 EP - 395 ER - TY - GEN A1 - Weiser, Martin A1 - Erdmann, Bodo A1 - Schenkl, Sebastian A1 - Muggenthaler, Holger A1 - Hubig, Michael A1 - Mall, Gita A1 - Zachow, Stefan T1 - Uncertainty in Temperature-Based Determination of Time of Death N2 - Temperature-based estimation of time of death (ToD) can be per- formed either with the help of simple phenomenological models of corpse cooling or with detailed mechanistic (thermodynamic) heat transfer mod- els. The latter are much more complex, but allow a higher accuracy of ToD estimation as in principle all relevant cooling mechanisms can be taken into account. The potentially higher accuracy depends on the accuracy of tissue and environmental parameters as well as on the geometric resolution. We in- vestigate the impact of parameter variations and geometry representation on the estimated ToD based on a highly detailed 3D corpse model, that has been segmented and geometrically reconstructed from a computed to- mography (CT) data set, differentiating various organs and tissue types. From that we identify the most crucial parameters to measure or estimate, and obtain a local uncertainty quantifcation for the ToD. T3 - ZIB-Report - 17-18 KW - forensic medicine KW - determination of time of death KW - heat transfer equation KW - sensitivity i.r.t. thermal parameters KW - sensitivity i.r.t. geometric resolution Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-63818 SN - 1438-0064 ER - TY - GEN A1 - Ambellan, Felix A1 - Tack, Alexander A1 - Wilson, Dave A1 - Anglin, Carolyn A1 - Lamecker, Hans A1 - Zachow, Stefan T1 - Evaluating two methods for Geometry Reconstruction from Sparse Surgical Navigation Data N2 - 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. T3 - ZIB-Report - 17-71 KW - Knee Arthroplasty KW - Sparse Geometry Reconstruction KW - Statistical Shape Models Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-66052 SN - 1438-0064 ER - TY - GEN A1 - Joachimsky, Robert A1 - Ambellan, Felix A1 - Zachow, Stefan T1 - Computerassistierte Auswahl und Platzierung von interpositionalen Spacern zur Behandlung früher Gonarthrose N2 - Degenerative Gelenkerkrankungen, wie die Osteoarthrose, sind ein häufiges Krankheitsbild unter älteren Erwachsenen. Hierbei verringert sich u.a. der Gelenkspalt aufgrund degenerierten Knorpels oder geschädigter Menisci. Ein in den Gelenkspalt eingebrachter interpositionaler Spacer soll die mit der Osteoarthrose einhergehende verringerte Gelenkkontaktfläche erhöhen und so der teilweise oder vollständige Gelenkersatz hinausgezögert oder vermieden werden. In dieser Arbeit präsentieren wir eine Planungssoftware für die Auswahl und Positionierung eines interpositionalen Spacers am Patientenmodell. Auf einer MRT-basierten Bildsegmentierung aufbauend erfolgt eine geometrische Rekonstruktion der 3D-Anatomie des Kniegelenks. Anhand dieser wird der Gelenkspalt bestimmt, sowie ein Spacer ausgewählt und algorithmisch vorpositioniert. Die Positionierung des Spacers ist durch den Benutzer jederzeit interaktiv anpassbar. Für jede Positionierung eines Spacers wird ein Fitness-Wert zur Knieanatomie des jeweiligen Patienten berechnet und den Nutzern Rückmeldung hinsichtlich Passgenauigkeit gegeben. Die Software unterstützt somit als Entscheidungshilfe die behandelnden Ärzte bei der patientenspezifischen Spacerauswahl. T3 - ZIB-Report - 17-72 KW - Osteoarthrose KW - Kniegelenk KW - Medizinische Planung KW - Visualisierung Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-66064 SN - 1438-0064 ER - TY - JOUR A1 - Bernard, Florian A1 - Salamanca, Luis A1 - Thunberg, Johan A1 - Tack, Alexander A1 - Jentsch, Dennis A1 - Lamecker, Hans A1 - Zachow, Stefan A1 - Hertel, Frank A1 - Goncalves, Jorge A1 - Gemmar, Peter T1 - Shape-aware Surface Reconstruction from Sparse 3D Point-Clouds JF - Medical Image Analysis N2 - 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. Y1 - 2017 UR - http://www.sciencedirect.com/science/article/pii/S1361841517300233 U6 - https://doi.org/10.1016/j.media.2017.02.005 VL - 38 SP - 77 EP - 89 ER - TY - JOUR A1 - Lube, Juliane A1 - Flack, Natasha A1 - Cotofana, Sebastian A1 - Özkurtul, Orkun A1 - Woodley, Stephanie A1 - Zachow, Stefan A1 - Hammer, Niels T1 - Pelvic and lower extremity physiological cross-sectional areas: An MRI study of the living young and comparison to published research literature JF - Surgical and Radiologic Anatomy N2 - Purpose: Morphological data pertaining to the pelvis and lower extremity muscles are increasingly being used in biomechanical modeling to compare healthy and pathological conditions. Very few data sets exist that encompass all of the muscles of the lower limb, allowing for comparisons between regions. The aims of this study were to (a) provide physiological cross-sectional area (PCSA) data for the pelvic, thigh, and leg muscles in young, healthy participants, using magnetic resonance imaging (MRI), and (b) to compare these data with summarized PCSAs obtained from the literature. Materials and Methods: Six young and healthy volunteers participated and were scanned using 3 T MRI. PCSAs were calculated from volumetric segmentations obtained bilaterally of 28 muscles/muscle groups of the pelvis, thigh, and leg. These data were compared to published, summarized PCSA data derived from cadaveric, computed tomography, MRI and ultrasound studies. Results: The PCSA of the pelvis, thigh, and leg muscles tended to be 20–130% larger in males than in females, except for the gemelli which were 34% smaller in males, and semitendinosus and triceps surae which did not differ (<20% different). The dominant and the non-dominant sides showed similar and minutely different PCSA with less than 18% difference between sides. Comparison to other studies revealed wide ranges within, and large differences between, the cadaveric and imaging PCSA data. Comparison of the PCSA of this study and published literature revealed major differences in the iliopsoas, gluteus minimus, tensor fasciae latae, gemelli, obturator internus, biceps femoris, quadriceps femoris, and the deep leg flexor muscles. Conclusions: These volume-derived PCSAs of the pelvic and lower limb muscles alongside the data synthesised from the literature may serve as a basis for comparative and biomechanical studies of the living and healthy young, and enable calculation of muscle forces. Comparison of the literature revealed large variations in PCSA from each of the different investigative modalities, hampering omparability between studies. Sample size, age, post-mortem changes of muscle tone, chemical fixation of cadaveric tissues, and the underlying physics of the imaging techniques may potentially influence PCSA calculations. Y1 - 2017 U6 - https://doi.org/10.1007/s00276-016-1807-6 VL - 39 IS - 8 SP - 849 EP - 857 ER - TY - JOUR A1 - Sahu, Manish A1 - Mukhopadhyay, Anirban A1 - Szengel, Angelika A1 - Zachow, Stefan T1 - Addressing multi-label imbalance problem of Surgical Tool Detection using CNN JF - International Journal of Computer Assisted Radiology and Surgery N2 - Purpose: A fully automated surgical tool detection framework is proposed for endoscopic video streams. State-of-the-art surgical tool detection methods rely on supervised one-vs-all or multi-class classification techniques, completely ignoring the co-occurrence relationship of the tools and the associated class imbalance. Methods: In this paper, we formulate tool detection as a multi-label classification task where tool co-occurrences are treated as separate classes. In addition, imbalance on tool co-occurrences is analyzed and stratification techniques are employed to address the imbalance during Convolutional Neural Network (CNN) training. Moreover, temporal smoothing is introduced as an online post-processing step to enhance run time prediction. Results: Quantitative analysis is performed on the M2CAI16 tool detection dataset to highlight the importance of stratification, temporal smoothing and the overall framework for tool detection. Conclusion: The analysis on tool imbalance, backed by the empirical results indicates the need and superiority of the proposed framework over state-of-the-art techniques. Y1 - 2017 UR - https://link.springer.com/article/10.1007/s11548-017-1565-x U6 - https://doi.org/10.1007/s11548-017-1565-x N1 - Selected for final oral presentation VL - 12 IS - 6 SP - 1013 EP - 1020 PB - Springer ER - TY - JOUR A1 - Schenkl, Sebastian A1 - Muggenthaler, Holger A1 - Hubig, Michael A1 - Erdmann, Bodo A1 - Weiser, Martin A1 - Zachow, Stefan A1 - Heinrich, Andreas A1 - Güttler, Felix Victor A1 - Teichgräber, Ulf A1 - Mall, Gita T1 - Automatic CT-based finite element model generation for temperature-based death time estimation: feasibility study and sensitivity analysis JF - International Journal of Legal Medicine N2 - Temperature based death time estimation is based either on simple phenomenological models of corpse cooling or on detailed physical heat transfer models. The latter are much more complex, but allow a higher accuracy of death time estimation as in principle all relevant cooling mechanisms can be taken into account. Here, a complete work flow for finite element based cooling simulation models is presented. The following steps are demonstrated on CT-phantoms: • CT-scan • Segmentation of the CT images for thermodynamically relevant features of individual geometries • Conversion of the segmentation result into a Finite Element (FE) simulation model • Computation of the model cooling curve • Calculation of the cooling time For the first time in FE-based cooling time estimation the steps from the CT image over segmentation to FE model generation are semi-automatically performed. The cooling time calculation results are compared to cooling measurements performed on the phantoms under controlled conditions. In this context, the method is validated using different CTphantoms. Some of the CT phantoms thermodynamic material parameters had to be experimentally determined via independent experiments. Moreover the impact of geometry and material parameter uncertainties on the estimated cooling time is investigated by a sensitivity analysis. KW - temperature based death time estimation KW - finite element method KW - CT segmentation KW - sensitivity analysis Y1 - 2017 U6 - https://doi.org/doi:10.1007/s00414-016-1523-0 VL - 131 IS - 3 SP - 699 EP - 712 ER - TY - GEN A1 - Sahu, Manish A1 - Dill, Sabrina A1 - Mukhopadyay, Anirban A1 - Zachow, Stefan T1 - Surgical Tool Presence Detection for Cataract Procedures N2 - This article outlines the submission to the CATARACTS challenge for automatic tool presence detection [1]. Our approach for this multi-label classification problem comprises labelset-based sampling, a CNN architecture and temporal smothing as described in [3], which we call ZIB-Res-TS. T3 - ZIB-Report - 18-28 Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-69110 SN - 1438-0064 ER - TY - JOUR A1 - Wilson, David A1 - Anglin, Carolyn A1 - Ambellan, Felix A1 - Grewe, Carl Martin A1 - Tack, Alexander A1 - Lamecker, Hans A1 - Dunbar, Michael A1 - Zachow, Stefan T1 - Validation of three-dimensional models of the distal femur created from surgical navigation point cloud data for intraoperative and postoperative analysis of total knee arthroplasty JF - International Journal of Computer Assisted Radiology and Surgery N2 - Purpose: Despite the success of total knee arthroplasty there continues to be a significant proportion of patients who are dissatisfied. One explanation may be a shape mismatch between pre and post-operative distal femurs. The purpose of this study was to investigate a method to match a statistical shape model (SSM) to intra-operatively acquired point cloud data from a surgical navigation system, and to validate it against the pre-operative magnetic resonance imaging (MRI) data from the same patients. Methods: A total of 10 patients who underwent navigated total knee arthroplasty also had an MRI scan less than 2 months pre-operatively. The standard surgical protocol was followed which included partial digitization of the distal femur. Two different methods were employed to fit the SSM to the digitized point cloud data, based on (1) Iterative Closest Points (ICP) and (2) Gaussian Mixture Models (GMM). The available MRI data were manually segmented and the reconstructed three-dimensional surfaces used as ground truth against which the statistical shape model fit was compared. Results: For both approaches, the difference between the statistical shape model-generated femur and the surface generated from MRI segmentation averaged less than 1.7 mm, with maximum errors occurring in less clinically important areas. Conclusion: The results demonstrated good correspondence with the distal femoral morphology even in cases of sparse data sets. Application of this technique will allow for measurement of mismatch between pre and post-operative femurs retrospectively on any case done using the surgical navigation system and could be integrated into the surgical navigation unit to provide real-time feedback. Y1 - 2017 UR - https://link.springer.com/content/pdf/10.1007%2Fs11548-017-1630-5.pdf U6 - https://doi.org/10.1007/s11548-017-1630-5 VL - 12 IS - 12 SP - 2097 EP - 2105 PB - Springer ER - TY - GEN A1 - Grewe, Carl Martin A1 - Zachow, Stefan ED - Doll, Nikola ED - Bredekamp, Horst ED - Schäffner, Wolfgang T1 - Face to Face-Interface T2 - +ultra. Knowledge & Gestaltung Y1 - 2017 SP - 320 EP - 321 PB - Seemann Henschel ER - TY - JOUR A1 - Zahn, Robert A1 - Grotjohann, Sarah A1 - Ramm, Heiko A1 - Zachow, Stefan A1 - Pumberger, Matthias A1 - Putzier, Michael A1 - Perka, Carsten A1 - Tohtz, Stephan T1 - Influence of pelvic tilt on functional acetabular orientation JF - Technology and Health Care Y1 - 2016 U6 - https://doi.org/10.3233/THC-161281 VL - 25 IS - 3 SP - 557 EP - 565 PB - IOS Press ER - TY - JOUR A1 - Lemanis, Robert A1 - Korn, Dieter A1 - Zachow, Stefan A1 - Rybacki, Erik A1 - Hoffmann, René T1 - The Evolution and Development of Cephalopod Chambers and Their Shape JF - PLOS ONE N2 - The Ammonoidea is a group of extinct cephalopods ideal to study evolution through deep time. The evolution of the planispiral shell and complexly folded septa in ammonoids has been thought to have increased the functional surface area of the chambers permitting enhanced metabolic functions such as: chamber emptying, rate of mineralization and increased growth rates throughout ontogeny. Using nano-computed tomography and synchrotron radiation based micro-computed tomography, we present the first study of ontogenetic changes in surface area to volume ratios in the phragmocone chambers of several phylogenetically distant ammonoids and extant cephalopods. Contrary to the initial hypothesis, ammonoids do not possess a persistently high relative chamber surface area. Instead, the functional surface area of the chambers is higher in earliest ontogeny when compared to Spirula spirula. The higher the functional surface area the quicker the potential emptying rate of the chamber; quicker chamber emptying rates would theoretically permit faster growth. This is supported by the persistently higher siphuncular surface area to chamber volume ratio we collected for the ammonite Amauroceras sp. compared to either S. spirula or nautilids. We demonstrate that the curvature of the surface of the chamber increases with greater septal complexity increasing the potential refilling rates. We further show a unique relationship between ammonoid chamber shape and size that does not exist in S. spirula or nautilids. This view of chamber function also has implications for the evolution of the internal shell of coleoids, relating this event to the decoupling of soft-body growth and shell growth. Y1 - 2016 U6 - https://doi.org/10.1371/journal.pone.0151404 VL - 11 IS - 3 ER - TY - JOUR A1 - Dunlop, Jason A1 - Apanaskevich, Dmitry A1 - Lehmann, Jens A1 - Hoffmann, Rene A1 - Fusseis, Florian A1 - Ehlke, Moritz A1 - Zachow, Stefan A1 - Xiao, Xianghui T1 - Microtomography of the Baltic amber tick Ixodes succineus reveals affinities with the modern Asian disease vector Ixodes ovatus JF - BMC Evolutionary Biology N2 - Background: Fossil ticks are extremely rare, whereby Ixodes succineus Weidner, 1964 from Eocene (ca. 44-49 Ma) Baltic amber is one of the oldest examples of a living hard tick genus (Ixodida: Ixodidae). Previous work suggested it was most closely related to the modern and widespread European sheep tick Ixodes ricinus (Linneaus, 1758). Results: Restudy using phase contrast synchrotron x-ray tomography yielded images of exceptional quality. These confirm the fossil's referral to Ixodes Latreille, 1795, but the characters resolved here suggest instead affinities with the Asian subgenus Partipalpiger Hoogstraal et al., 1973 and its single living (and medically significant) species Ixodes ovatus Neumann, 1899. We redescribe the amber fossil here as Ixodes (Partipalpiger) succineus. Conclusions: Our data suggest that Ixodes ricinus is unlikely to be directly derived from Weidner's amber species, but instead reveals that the Partipalpiger lineage was originally more widely distributed across the northern hemisphere. The closeness of Ixodes (P.) succineus to a living vector of a wide range of pathogens offers the potential to correlate its spatial and temporal position (northern Europe, nearly 50 million years ago) with the estimated origination dates of various tick-borne diseases. Y1 - 2016 U6 - https://doi.org/10.1186/s12862-016-0777-y VL - 16 IS - 1 ER - TY - JOUR A1 - Bernard, Florian A1 - Salamanca, Luis A1 - Thunberg, Johan A1 - Tack, Alexander A1 - Jentsch, Dennis A1 - Lamecker, Hans A1 - Zachow, Stefan A1 - Hertel, Frank A1 - Goncalves, Jorge A1 - Gemmar, Peter T1 - Shape-aware Surface Reconstruction from Sparse Data JF - arXiv N2 - 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. Y1 - 2016 SP - 1602.08425v1 ER - TY - JOUR A1 - Lemanis, Robert A1 - Zachow, Stefan A1 - Hoffmann, René T1 - Comparative cephalopod shell strength and the role of septum morphology on stress distribution JF - PeerJ N2 - The evolution of complexly folded septa in ammonoids has long been a controversial topic. Explanations of the function of these folded septa can be divided into physiological and mechanical hypotheses with the mechanical functions tending to find widespread support. The complexity of the cephalopod shell has made it difficult to directly test the mechanical properties of these structures without oversimplification of the septal morphology or extraction of a small sub-domain. However, the power of modern finite element analysis now permits direct testing of mechanical hypothesis on complete, empirical models of the shells taken from computed tomographic data. Here we compare, for the first time using empirical models, the capability of the shells of extant Nautilus pompilius, Spirula spirula, and the extinct ammonite Cadoceras sp. to withstand hydrostatic pressure and point loads. Results show hydrostatic pressure imparts highest stress on the final septum with the rest of the shell showing minimal compression. S. spirula shows the lowest stress under hydrostatic pressure while N. pompilius shows the highest stress. Cadoceras sp. shows the development of high stress along the attachment of the septal saddles with the shell wall. Stress due to point loads decreases when the point force is directed along the suture as opposed to the unsupported chamber wall. Cadoceras sp. shows the greatest decrease in stress between the point loads compared to all other models. Greater amplitude of septal flutes corresponds with greater stress due to hydrostatic pressure; however, greater amplitude decreases the stress magnitude of point loads directed along the suture. In our models, sutural complexity does not predict greater resistance to hydrostatic pressure but it does seem to increase resistance to point loads, such as would be from predators. This result permits discussion of palaeoecological reconstructions on the basis of septal morphology. We further suggest that the ratio used to characterize septal morphology in the septal strength index and in calculations of tensile strength of nacre are likely insufficient. A better understanding of the material properties of cephalopod nacre may allow the estimation of maximum depth limits of shelled cephalopods through finite element analysis. Y1 - 2016 U6 - https://doi.org/10.7717/peerj.2434 VL - 4 SP - e2434 ER - TY - CHAP A1 - Mukhopadhyay, Anirban A1 - Morillo, Oscar A1 - Zachow, Stefan A1 - Lamecker, Hans T1 - Robust and Accurate Appearance Models Based on Joint Dictionary Learning Data from the Osteoarthritis Initiative T2 - Lecture Notes in Computer Science, Patch-Based Techniques in Medical Imaging. Patch-MI 2016 N2 - Deformable model-based approaches to 3D image segmentation have been shown to be highly successful. Such methodology requires an appearance model that drives the deformation of a geometric model to the image data. Appearance models are usually either created heuristically or through supervised learning. Heuristic methods have been shown to work effectively in many applications but are hard to transfer from one application (imaging modality/anatomical structure) to another. On the contrary, supervised learning approaches can learn patterns from a collection of annotated training data. In this work, we show that the supervised joint dictionary learning technique is capable of overcoming the traditional drawbacks of the heuristic approaches. Our evaluation based on two different applications (liver/CT and knee/MR) reveals that our approach generates appearance models, which can be used effectively and efficiently in a deformable model-based segmentation framework. Y1 - 2016 U6 - https://doi.org/10.1007/978-3-319-47118-1_4 VL - 9993 SP - 25 EP - 33 ER - TY - CHAP A1 - Lamecker, Hans A1 - Zachow, Stefan T1 - Statistical Shape Modeling of Musculoskeletal Structures and Its Applications T2 - Computational Radiology for Orthopaedic Interventions N2 - Statistical shape models (SSM) describe the shape variability contained in a given population. They are able to describe large populations of complex shapes with few degrees of freedom. This makes them a useful tool for a variety of tasks that arise in computer-aided madicine. In this chapter we are going to explain the basic methodology of SSMs and present a variety of examples, where SSMs have been successfully applied. Y1 - 2016 SN - 978-3-319-23481-6 U6 - https://doi.org/10.1007/978-3-319-23482-3 VL - 23 SP - 1 EP - 23 PB - Springer ER - TY - GEN A1 - Tycowicz, Christoph von A1 - Ambellan, Felix A1 - Mukhopadhyay, Anirban A1 - Zachow, Stefan T1 - A Riemannian Statistical Shape Model using Differential Coordinates N2 - We propose a novel Riemannian framework for statistical analysis of shapes that is able to account for the nonlinearity in shape variation. By adopting a physical perspective, we introduce a differential representation that puts the local geometric variability into focus. We model these differential coordinates as elements of a Lie group thereby endowing our shape space with a non-Euclidian structure. A key advantage of our framework is that statistics in a manifold shape space become numerically tractable improving performance by several orders of magnitude over state-of-the-art. We show that our Riemannian model is well suited for the identification of intra-population variability as well as inter-population differences. In particular, we demonstrate the superiority of the proposed model in experiments on specificity and generalization ability. We further derive a statistical shape descriptor that outperforms the standard Euclidian approach in terms of shape-based classification of morphological disorders. T3 - ZIB-Report - 16-69 Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-61175 UR - https://opus4.kobv.de/opus4-zib/frontdoor/index/index/docId/6485 SN - 1438-0064 ER - TY - CHAP A1 - Grewe, Carl Martin A1 - Zachow, Stefan T1 - Fully Automated and Highly Accurate Dense Correspondence for Facial Surfaces T2 - Computer Vision – ECCV 2016 Workshops N2 - We present a novel framework for fully automated and highly accurate determination of facial landmarks and dense correspondence, e.g. a topologically identical mesh of arbitrary resolution, across the entire surface of 3D face models. For robustness and reliability of the proposed approach, we are combining 2D landmark detectors and 3D statistical shape priors with a variational matching method. Instead of matching faces in the spatial domain only, we employ image registration to align the 2D parametrization of the facial surface to a planar template we call the Unified Facial Parameter Domain (ufpd). This allows us to simultaneously match salient photometric and geometric facial features using robust image similarity measures while reasonably constraining geometric distortion in regions with less significant features. We demonstrate the accuracy of the dense correspondence established by our framework on the BU3DFE database with 2500 facial surfaces and show, that our framework outperforms current state-of-the-art methods with respect to the fully automated location of facial landmarks. Y1 - 2016 U6 - https://doi.org/10.1007/978-3-319-48881-3_38 VL - 9914 SP - 552 EP - 568 PB - Springer International Publishing ER - TY - JOUR A1 - Zachow, Stefan T1 - Computational Planning in Facial Surgery JF - Facial Plastic Surgery N2 - This article reflects the research of the last two decades in computational planning for cranio-maxillofacial surgery. Model-guided and computer-assisted surgery planning has tremendously developed due to ever increasing computational capabilities. Simulators for education, planning, and training of surgery are often compared with flight simulators, where maneuvers are also trained to reduce a possible risk of failure. Meanwhile, digital patient models can be derived from medical image data with astonishing accuracy and thus can serve for model surgery to derive a surgical template model that represents the envisaged result. Computerized surgical planning approaches, however, are often still explorative, meaning that a surgeon tries to find a therapeutic concept based on his or her expertise using computational tools that are mimicking real procedures. Future perspectives of an improved computerized planning may be that surgical objectives will be generated algorithmically by employing mathematical modeling, simulation, and optimization techniques. Planning systems thus act as intelligent decision support systems. However, surgeons can still use the existing tools to vary the proposed approach, but they mainly focus on how to transfer objectives into reality. Such a development may result in a paradigm shift for future surgery planning. Y1 - 2015 U6 - https://doi.org/10.1055/s-0035-1564717 VL - 31 IS - 5 SP - 446 EP - 462 ER - TY - GEN A1 - Lamas-Rodríguez, Julián A1 - Ehlke, Moritz A1 - Hoffmann, René A1 - Zachow, Stefan T1 - GPU-accelerated denoising of large tomographic data sets with low SNR BT - Application for non-invasive analysis of paleontological data N2 - Enhancements in tomographic imaging techniques facilitate non-destructive methods for visualizing fossil structures. However, to penetrate dense materials such as sediments or pyrites, image acquisition is typically performed with high beam energy and very sensitive image intensifiers, leading to artifacts and noise in the acquired data. The analysis of delicate fossil structures requires the images to be captured in maximum resolution, resulting in large data sets of several giga bytes (GB) in size. Since the structural information of interest is often almost in the same spatial range as artifacts and noise, image processing and segmentation algorithms have to cope with a very low signal-to-noise ratio (SNR). Within this report we present a study on the performance of a collection of denoising algorithms applied to a very noisy fossil dataset. The study shows that a non-local means (NLM) filter, in case it is properly configured, is able to remove a considerable amount of noise while preserving most of the structural information of interest. Based on the results of this study, we developed a software tool within ZIBAmira that denoises large tomographic datasets using an adaptive, GPU-accelerated NLM filter. With the help of our implementation a user can interactively configure the filter's parameters and thus its effectiveness with respect to the data of interest, while the filtering response is instantly visualized for a preselected region of interest (ROI). Our implementation efficiently denoises even large fossil datasets in a reasonable amount of time. T3 - ZIB-Report - 15-14 Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-56339 SN - 1438-0064 ER - TY - JOUR A1 - Zahn, Robert A1 - Grotjohann, Sarah A1 - Ramm, Heiko A1 - Zachow, Stefan A1 - Putzier, Michael A1 - Perka, Carsten A1 - Tohtz, Stephan T1 - Pelvic tilt compensates for increased acetabular anteversion JF - International Orthopaedics N2 - Pelvic tilt determines functional orientation of the acetabulum. In this study, we investigated the interaction of pelvic tilt and functional acetabular anteversion (AA) in supine position. Y1 - 2015 U6 - https://doi.org/10.1007/s00264-015-2949-6 VL - 40 IS - 8 SP - 1571 EP - 1575 ER - TY - JOUR A1 - Zachow, Stefan A1 - Heppt, Werner T1 - The Facial Profile JF - Facial Plastic Surgery N2 - Facial appearance in our societies is often associated with notions of attractiveness, juvenileness, beauty, success, and so forth. Hence, the role of facial plastic surgery is highly interrelated to a patient's desire to feature many of these positively connoted attributes, which of course, are subject of different cultural perceptions or social trends. To judge about somebody's facial appearance, appropriate quantitative measures as well as methods to obtain and compare individual facial features are required. This special issue on facial profile is intended to provide an overview on how facial characteristics are surgically managed in an interdisciplinary way based on experience, instrumentation, and modern technology to obtain an aesthetic facial appearance with harmonious facial proportions. The facial profile will be discussed within the context of facial aesthetics. Latest concepts for capturing facial morphology in high speed and impressive detail are presented for quantitative analysis of even subtle changes, aging effects, or facial expressions. In addition, the perception of facial profiles is evaluated based on eye tracking technology. Y1 - 2015 U6 - https://doi.org/10.1055/s-0035-1566132 VL - 31 IS - 5 SP - 419 EP - 420 ER - TY - CHAP A1 - Ehlke, Moritz A1 - Frenzel, Thomas A1 - Ramm, Heiko A1 - Shandiz, Mohsen Akbari A1 - Anglin, Carolyn A1 - Zachow, Stefan T1 - Towards Robust Measurement Of Pelvic Parameters From AP Radiographs Using Articulated 3D Models T2 - Computer Assisted Radiology and Surgery (CARS) N2 - Patient-specific parameters such as the orientation of the acetabulum or pelvic tilt are useful for custom planning for total hip arthroplasty (THA) and for evaluating the outcome of surgical interventions. The gold standard in obtaining pelvic parameters is from three-dimensional (3D) computed tomography (CT) imaging. However, this adds time and cost, exposes the patient to a substantial radiation dose, and does not allow for imaging under load (e.g. while the patient is standing). If pelvic parameters could be reliably derived from the standard anteroposterior (AP) radiograph, preoperative planning would be more widespread, and research analyses could be applied to retrospective data, after a postoperative issue is discovered. The goal of this work is to enable robust measurement of two surgical parameters of interest: the tilt of the anterior pelvic plane (APP) and the orientation of the natural acetabulum. We present a computer-aided reconstruction method to determine the APP and natural acetabular orientation from a single, preoperative X-ray. It can easily be extended to obtain other important preoperative and postoperative parameters solely based on a single AP radiograph. KW - Articulated Models, Statistical Shape And Intensity Models, 2D/3D Anatomy Reconstruction, Pelvic Parameters Measurement, Total Hip Arthroplasty Y1 - 2015 ER - TY - GEN A1 - Ehlke, Moritz A1 - Frenzel, Thomas A1 - Ramm, Heiko A1 - Shandiz, Mohsen Akbari A1 - Anglin, Carolyn A1 - Zachow, Stefan T1 - Towards Robust Measurement of Pelvic Parameters from AP Radiographs using Articulated 3D Models N2 - Patient-specific parameters such as the orientation of the acetabulum or pelvic tilt are useful for custom planning for total hip arthroplasty (THA) and for evaluating the outcome of surgical interventions. The gold standard in obtaining pelvic parameters is from three-dimensional (3D) computed tomography (CT) imaging. However, this adds time and cost, exposes the patient to a substantial radiation dose, and does not allow for imaging under load (e.g. while the patient is standing). If pelvic parameters could be reliably derived from the standard anteroposterior (AP) radiograph, preoperative planning would be more widespread, and research analyses could be applied to retrospective data, after a postoperative issue is discovered. The goal of this work is to enable robust measurement of two surgical parameters of interest: the tilt of the anterior pelvic plane (APP) and the orientation of the natural acetabulum. We present a computer-aided reconstruction method to determine the APP and natural acetabular orientation from a single, preoperative X-ray. It can easily be extended to obtain other important preoperative and postoperative parameters solely based on a single AP radiograph. T3 - ZIB-Report - 15-11 KW - Articulated Models, Statistical Shape And Intensity Models, 2D/3D Anatomy Reconstruction, Pelvic Parameters Measurement, Total Hip Arthroplasty Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-53707 SN - 1438-0064 ER - TY - JOUR A1 - Lemanis, Robert A1 - Zachow, Stefan A1 - Fusseis, Florian A1 - Hoffmann, René T1 - A new approach using high-resolution computed tomography to test the buoyant properties of chambered cephalopod shells JF - Paleobiology N2 - The chambered shell of modern cephalopods functions as a buoyancy apparatus, allowing the animal to enter the water column without expending a large amount of energy to overcome its own weight. Indeed, the chambered shell is largely considered a key adaptation that allowed the earliest cephalopods to leave the ocean floor and enter the water column. It has been argued by some, however, that the iconic chambered shell of Paleozoic and Mesozoic ammonoids did not provide a sufficiently buoyant force to compensate for the weight of the entire animal, thus restricting ammonoids to a largely benthic lifestyle reminiscent of some octopods. Here we develop a technique using high-resolution computed tomography to quantify the buoyant properties of chambered shells without reducing the shell to ideal spirals or eliminating inherent biological variability by using mathematical models that characterize past work in this area. This technique has been tested on Nautilus pompilius and is now extended to the extant deep-sea squid Spirula spirula and the Jurassic ammonite Cadoceras sp. hatchling. Cadoceras is found to have possessed near-neutral to positive buoyancy if hatched when the shell possessed between three and five chambers. However, we show that the animal could also overcome degrees of negative buoyancy through swimming, similar to the paralarvae of modern squids. These calculations challenge past inferences of benthic life habits based solely on calculations of negative buoyancy. The calculated buoyancy of Cadoceras supports the possibility of planktonic dispersal of ammonite hatchlings. This information is essential to understanding ammonoid ecology as well as biotic interactions and has implications for the interpretation of geochemical data gained from the isotopic analysis of the shell. Y1 - 2015 U6 - https://doi.org/10.1017/pab.2014.17 VL - 41 IS - 2 SP - 313 EP - 329 PB - Cambridge University Press CY - Cambridge ER - TY - GEN A1 - Wilson, David A1 - Bücher, Pia A1 - Grewe, Carl Martin A1 - Anglin, Carolyn A1 - Zachow, Stefan A1 - Michael, Dunbar T1 - Validation of Three Dimensional Models of the Distal Femur Created from Surgical Navigation Point Cloud Data T2 - 15th Annual Meeting of the International Society for Computer Assisted Orthopaedic Surgery (CAOS) Y1 - 2015 ER - TY - JOUR A1 - Grewe, Carl Martin A1 - Schreiber, Lisa A1 - Zachow, Stefan T1 - Fast and Accurate Digital Morphometry of Facial Expressions JF - Facial Plastic Surgery Y1 - 2015 U6 - https://doi.org/10.1055/s-0035-1564720 VL - 31 IS - 05 SP - 431 EP - 438 PB - Thieme Medical Publishers CY - New York ER - TY - GEN A1 - Wilson, David A1 - Bücher, Pia A1 - Grewe, Carl Martin A1 - Mocanu, Valentin A1 - Anglin, Carolyn A1 - Zachow, Stefan A1 - Dunbar, Michael T1 - Validation of Three Dimensional Models of the Distal Femur Created from Surgical Navigation Data T2 - Orthopedic Research Society Annual Meeting Y1 - 2015 CY - Las Vegas, Nevada ER - TY - GEN A1 - Ehlke, Moritz A1 - Heyland, Mark A1 - Märdian, Sven A1 - Duda, Georg A1 - Zachow, Stefan T1 - Assessing the Relative Positioning of an Osteosynthesis Plate to the Patient-Specific Femoral Shape from Plain 2D Radiographs N2 - We present a novel method to derive the surface distance of an osteosynthesis plate w.r.t. the patient­specific surface of the distal femur based on 2D X­ray images. Our goal is to study from clinical data, how the plate­to­bone distance affects bone healing. The patient­specific 3D shape of the femur is, however, seldom recorded for cases of femoral osteosynthesis since this typically requires Computed Tomography (CT), which comes at high cost and radiation dose. Our method instead utilizes two postoperative X­ray images to derive the femoral shape and thus can be applied on radiographs that are taken in clinical routine for follow­up. First, the implant geometry is used as a calibration object to relate the implant and the individual X­ray images spatially in a virtual X­ray setup. In a second step, the patient­specific femoral shape and pose are reconstructed in the virtual setup by fitting a deformable statistical shape and intensity model (SSIM) to the images. The relative positioning between femur and implant is then assessed in terms of displacement between the reconstructed 3D shape of the femur and the plate. A preliminary evaluation based on 4 cadaver datasets shows that the method derives the plate­to­bone distance with a mean absolute error of less than 1mm and a maximum error of 4.7 mm compared to ground truth from CT. We believe that the approach presented in this paper constitutes a meaningful tool to elucidate the effect of implant positioning on fracture healing. T3 - ZIB-Report - 15-21 KW - 3d-­reconstruction from 2d X­rays KW - statistical shape and intensity models KW - fracture fixation of the distal femur KW - pose estimation Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-54268 SN - 1438-0064 ER - TY - GEN A1 - Ehlke, Moritz A1 - Heyland, Mark A1 - Märdian, Sven A1 - Duda, Georg A1 - Zachow, Stefan T1 - 3D Assessment of Osteosynthesis based on 2D Radiographs N2 - We present a novel method to derive the surface distance of an osteosynthesis plate w.r.t. the patient-specific surface of the distal femur based on postoperative 2D radiographs. In a first step, the implant geometry is used as a calibration object to relate the implant and the individual X-ray images spatially in a virtual X-ray setup. Second, the patient-specific femoral shape and pose are reconstructed by fitting a deformable statistical shape and intensity model (SSIM) to the X-rays. The relative positioning between femur and implant is then assessed in terms of the displacement between the reconstructed 3D shape of the femur and the plate. We believe that the approach presented in this paper constitutes a meaningful tool to elucidate the effect of implant positioning on fracture healing and, ultimately, to derive load recommendations after surgery. T3 - ZIB-Report - 15-47 KW - 3d-reconstruction from 2d X-rays KW - statistical shape and intensity models KW - osteosynthesis follow-up Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-56217 SN - 1438-0064 ER - TY - CHAP A1 - Krämer, Martin A1 - Herrmann, Karl-Heinz A1 - Boeth, Heide A1 - Tycowicz, Christoph von A1 - König, Christian A1 - Zachow, Stefan A1 - Ehrig, Rainald A1 - Hege, Hans-Christian A1 - Duda, Georg A1 - Reichenbach, Jürgen T1 - Measuring 3D knee dynamics using center out radial ultra-short echo time trajectories with a low cost experimental setup T2 - ISMRM (International Society for Magnetic Resonance in Medicine), 23rd Annual Meeting 2015, Toronto, Canada Y1 - 2015 ER - TY - CHAP A1 - Ehlke, Moritz A1 - Heyland, Mark A1 - Märdian, Sven A1 - Duda, Georg A1 - Zachow, Stefan T1 - Assessing the relative positioning of an osteosynthesis plate to the patient-specific femoral shape from plain 2D radiographs T2 - Proceedings of the 15th Annual Meeting of CAOS-International (CAOS) N2 - We present a novel method to derive the surface distance of an osteosynthesis plate w.r.t. the patient­specific surface of the distal femur based on 2D X­ray images. Our goal is to study from clinical data, how the plate­to­bone distance affects bone healing. The patient­specific 3D shape of the femur is, however, seldom recorded for cases of femoral osteosynthesis since this typically requires Computed Tomography (CT), which comes at high cost and radiation dose. Our method instead utilizes two postoperative X­ray images to derive the femoral shape and thus can be applied on radiographs that are taken in clinical routine for follow­up. First, the implant geometry is used as a calibration object to relate the implant and the individual X­ray images spatially in a virtual X­ray setup. In a second step, the patient­specific femoral shape and pose are reconstructed in the virtual setup by fitting a deformable statistical shape and intensity model (SSIM) to the images. The relative positioning between femur and implant is then assessed in terms of displacement between the reconstructed 3D shape of the femur and the plate. A preliminary evaluation based on 4 cadaver datasets shows that the method derives the plate­to­bone distance with a mean absolute error of less than 1mm and a maximum error of 4.7 mm compared to ground truth from CT. We believe that the approach presented in this paper constitutes a meaningful tool to elucidate the effect of implant positioning on fracture healing. KW - 3d-­reconstruction from 2d X­rays KW - statistical shape and intensity models KW - fracture fixation of the distal femur KW - pose estimation Y1 - 2015 ER - TY - CHAP A1 - Ehlke, Moritz A1 - Heyland, Mark A1 - Märdian, Sven A1 - Duda, Georg A1 - Zachow, Stefan T1 - 3D Assessment of Osteosynthesis based on 2D Radiographs T2 - Proceedings of the Jahrestagung der Deutschen Gesellschaft für Computer- und Roboterassistierte Chirurgie (CURAC) N2 - We present a novel method to derive the surface distance of an osteosynthesis plate w.r.t. the patient-specific surface of the distal femur based on postoperative 2D radiographs. In a first step, the implant geometry is used as a calibration object to relate the implant and the individual X-ray images spatially in a virtual X-ray setup. Second, the patient- specific femoral shape and pose are reconstructed by fitting a deformable statistical shape and intensity model (SSIM) to the X-rays. The relative positioning between femur and implant is then assessed in terms of the displacement between the reconstructed 3D shape of the femur and the plate. We believe that the approach presented in this paper constitutes a meaningful tool to elucidate the effect of implant positioning on fracture healing and, ultimately, to derive load recommendations after surgery. KW - 3d-reconstruction from 2d X-rays KW - osteosynthesis follow-up KW - statistical shape and intensity models Y1 - 2015 SP - 317 EP - 321 ER - TY - JOUR A1 - Campoli, Gianni A1 - Baka, Nora A1 - Kaptein, Bart A1 - Valstar, Edward A1 - Zachow, Stefan A1 - Weinans, Harrie A1 - Zadpoor, Amir Abbas T1 - Relationship between the shape and density distribution of the femur and its natural frequencies of vibration JF - Journal of Biomechanics N2 - It has been recently suggested that mechanical loads applied at frequencies close to the natural frequencies of bone could enhance bone apposition due to the resonance phenomenon. Other applications of bone modal analysis are also suggested. For the above-mentioned applications, it is important to understand how patient-specific bone shape and density distribution influence the natural frequencies of bones. We used finite element models to study the effects of bone shape and density distribution on the natural frequencies of the femur in free boundary conditions. A statistical shape and appearance model that describes shape and density distribution independently was created, based on a training set of 27 femora. The natural frequencies were then calculated for different shape modes varied around the mean shape while keeping the mean density distribution, for different appearance modes around the mean density distribution while keeping the mean bone shape, and for the 27 training femora. Single shape or appearance modes could cause up to 15% variations in the natural frequencies with certain modes having the greatest impact. For the actual femora, shape and density distribution changed the natural frequencies by up to 38%. First appearance mode that describes the general cortical bone thickness and trabecular bone density had one of the strongest impacts. The first appearance mode could therefore provide a sensitive measure of general bone health and disease progression. Since shape and density could cause large variations in the calculated natural frequencies, patient-specific FE models are needed for accurate estimation of bone natural frequencies. Y1 - 2014 U6 - https://doi.org/10.1016/j.jbiomech.2014.08.008 VL - 47 SP - 3334 EP - 3343 PB - Elsevier ER -