@inproceedings{Paetsch2019, author = {Paetsch, Olaf}, title = {Possibilities and Limitations of Automatic Feature Extraction shown by the Example of Crack Detection in 3D-CT Images of Concrete Specimen}, booktitle = {iCT 2019}, year = {2019}, abstract = {To assess the influence of the alkali-silica reaction (ASR) on pavement concrete 3D-CT imaging has been applied to concrete samples. Prior to imaging these samples have been drilled out of a concrete beam pre-damaged by fatigue loading. The resulting high resolution 3D-CT images consist of several gigabytes of voxels. Current desktop computers can visualize such big datasets without problems but a visual inspection or manual segmentation of features such as cracks by experts can only be carried out on a few slices. A quantitative analysis of cracks requires a segmentation of the whole specimen which could only be done by an automatic feature detection. This arises the question of the reliability of an automatic crack detection algorithm, its certainty and limitations. Does the algorithm find all cracks? Does it find too many cracks? Can parameters of that algorithm, once identified as good, be applied to other samples as well? Can ensemble computing with many crack parameters overcome the difficulties with parameter finding? By means of a crack detection algorithm based on shape recognition (template matching) these questions will be discussed. Since the author has no access to reliable ground truth data of cracks the assessment of the certainty of the automatic crack is restricted to visual inspection by experts. Therefore, an artificial dataset based on a combination of manually segmented cracks processed together with simple image processing algorithms is used to quantify the accuracy of the crack detection algorithm. Part of the evaluation of cracks in concrete samples is the knowledge of the surrounding material. The surrounding material can be used to assess the detected cracks, e.g. micro-cracks within the aggregate-matrix interface may be starting points for cracks on a macro scale. Furthermore, the knowledge of the surrounding material can help to find better parameter sets for the crack detection itself because crack characteristics may vary depending on their surrounding material. Therefore, in addition to a crack detection a complete segmentation of the sample into the components of concrete, such as aggregates, cement matrix and pores is needed. Since such a segmentation task cannot be done manually due to the amount of data, an approach utilizing convolutional neuronal networks stemming from a medical application has been applied. The learning phase requires a ground truth i.e. a segmentation of the components. This has to be created manually in a time-consuming task. However, this segmentation can be used for a quantitative evaluation of the automatic segmentation afterwards. Even though that work has been performed as a short term subtask of a bigger project funded by the German Research Foundation (DFG) this paper discusses problems which may arise in similar projects, too. [1.2MB | id=23664 ] iCT 2019 Session: Short talks Thu 13:50 Auditorium 2019-03 M{\"o}glichkeiten und Grenzen automatischer Merkmalserkennung am Beispiel von Risserkennungen in 3D-CT-Aufnahmen von Betonproben O. Paetsch11 Visualisation and Data Analysis; Konrad-Zuse-Institut Berlin (ZIB)18, Berlin, Germany Abstract [1MB | id=23104 ] DE DGZfP 2018 Session: Bauwesen 2018-09 Quantitative Rissanalyse im Fahrbahndeckenbeton mit der 3D-Computertomographie D. Meinel125, K. Ehrig128, F. Weise16, O. Paetsch211 1Division 8.5; BAM Federal Institute for Materials Research and Testing1277, Berlin, Germany 2Visualisation and Data Analysis; Konrad-Zuse-Institut Berlin (ZIB)18, Berlin, Germany concrete, ROI tomography, in-situ-CT, 3D-CT, Beton, AKR, Feuchtetransport, automatic crack detection Abstract [0.7MB | id=18980 ] DE DGZfP 2015 Session: CT Algorithmen 2016-04 3D Corrosion Detection in Time-dependent CT Images of Concrete O. Paetsch111, D. Baum15, S. Prohaska17, K. Ehrig228, D. Meinel225, G. Ebell24 1Visualisation and Data Analysis; Konrad-Zuse-Institut Berlin (ZIB)18, Berlin, Germany 2Division 8.5; BAM Federal Institute for Materials Research and Testing1277, Berlin, Germany CT, multi-angle radiography, defect detection, Feature Extraction, image processing, concrete, corrosion Abstract [0.5MB | id=18043 ] DIR 2015 Session: Quantitative imaging and image processing 2015-08 Korrosionsverfolgung in 3D-computertomographischen Aufnahmen von Stahlbetonproben O. Paetsch111, D. Baum15, G. Ebell24, K. Ehrig228, A. Heyn2, D. Meinel225, S. Prohaska17 1Konrad-Zuse-Institut Berlin (ZIB)18, Berlin, Germany 2Division VIII.3; BAM Federal Institute for Materials Research and Testing1277, Berlin, Germany Computertomographie [0.4MB | id=17375 ] DE DGZfP 2014 Session: Bauwesen 2015-03 Examination of Damage Processes in Concrete with CT D. Meinel125, K. Ehrig128, V. L'Hostis2, B. Muzeau2, O. Paetsch311 1BAM Federal Institute for Materials Research and Testing1277, Berlin, Germany 2Laboratoire d'Etude du Comportement des B{\´e}tons et des Argiles; Commissariat Energie Atomique (CEA)287, Gif-Sur-Yvette, France 3Konrad-Zuse-Institut Berlin (ZIB)18, Berlin, Germany X-ray computed tomography, concrete, corrosion, crack detection, 3D visualization Abstract [4.9MB | id=15692 ] iCT 2014 Session: Non-destructive Testing and 3D Materials Characterisation of... 2014-06 3-D-Visualisierung und statistische Analyse von Rissen in mit Computer-Tomographie untersuchten Betonproben O. Paetsch111, D. Baum15, D. Breßler1, K. Ehrig228, D. Meinel225, S. Prohaska1,17 1Konrad-Zuse-Institut Berlin (ZIB)18, Berlin, Germany 2Division VIII.3; BAM Federal Institute for Materials Research and Testing1277, Berlin, Germany Radiographic Testing (RT), statistical analysis, 3D Computed Tomography, visualization, concrete structural damage, automated crack detection [1MB | id=15343 ] DE DGZfP 2013 Session: Computertomographie 2014-03 Vergleich automatischer 3D-Risserkennungsmethoden f{\"u}r die quantitative Analyse der Schadensentwicklung in Betonproben mit Computertomographie O. Paetsch111, K. Ehrig228, D. Meinel225, D. Baum15, S. Prohaska1,1,17 1Konrad-Zuse-Institut Berlin (ZIB)18, Berlin, Germany 2Division VIII.3; BAM Federal Institute for Materials Research and Testing1277, Berlin, Germany Radiographic Testing (RT), visualization, crack detection, Visualisierung, computer tomography, template matching, Hessian eigenvalues, ZIBAmira, automated crack detection, percolation [0.9MB | id=14269 ] DE DGZfP 2012 Session: Computertomographie 2013-05 Automated 3D Crack Detection for Analyzing Damage Processes in Concrete with Computed Tomography O. Paetsch111, D. Baum15, K. Ehrig228, D. Meinel225, S. Prohaska1,1,17 1Konrad-Zuse-Institut Berlin (ZIB)18, Berlin, Germany 2Division VIII.3; BAM Federal Institute for Materials Research and Testing1277, Berlin, Germany computed tomography, template matching, Hessian eigenvalues, crack statistics, visualization, crack surface, ZIBAmira [0.6MB | id=13736 ] iCT 2012 Session: Poster - Analysis and Algorithms 2012-12 3-D-Visualisierung von Radar- und Ultraschallecho-Daten mit ZIBAmira D. Streicher112, O. Paetsch211, R. Seiler2, S. Prohaska27, M. Krause360 [Profile of Krause] , C. Boller178 1Saarland University74, Saarbr{\"u}cken, Germany 2Konrad-Zuse-Institut Berlin (ZIB)18, Berlin, Germany 3BAM Federal Institute for Materials Research and Testing1277, Berlin, Germany [0.4MB | id=12284 ] DE DGZfP 2011 Session: Bauwesen 2012-05 Comparison of Crack Detection Methods for Analyzing Damage Processes in Concrete with Computed Tomography K. Ehrig128, J. Goebbels153, D. Meinel125, O. Paetsch211, S. Prohaska27, V. Zobel2 1Division VIII.3; BAM Federal Institute for Materials Research and Testing1277, Berlin, Germany 2Konrad-Zuse-Institut Berlin (ZIB)18, Berlin, Germany [0.7MB | id=11150 ] DIR 2011 Session: Poster 2011-11 Actual Cooperations 10th International Workshop NDT in Progress 2019 2019 Oct 7-9 11th International Symposium on NDT in Aerospace 2019 2019 Nov 13-15 3rd Singapore International NDT Conference \& Exhibition, SINCE 2019 2019 Dec 4-5 10th Conference on Industrial Computed Tomography (iCT) 2020 2020 Feb 4-7 34th European Conference on Acoustic Emission Testing (EWGAE 2020) 2020 Sep 9-11 Contribute Papers and Proceedings to NDT.net Share... Home Exhibition Archive Forum Jobs Members Events Directory NDT A-Z Advertise Privacy Policy Contact About © NDT.net - Where expertise comes together. The Largest Open Access Portal of Nondestructive Testing (NDT)- since 1996}, language = {en} } @article{AlchikhConradMaetal.2019, author = {Alchikh, Maren and Conrad, Tim and Ma, Xiaolin and Broberg, Eeva K. and Penttinen, P. and Reiche, J. and Biere, Barbara and Schweiger, Brunhilde and Rath, Barbara and Hoppe, Christian}, title = {Are we missing respiratory viral infections in infants and children? Comparison of a hospital-based quality management system with standard of care}, volume = {25}, journal = {Clinical Microbiology and Infection}, number = {3}, issn = {1469-0691}, doi = {10.1016/j.cmi.2018.05.023}, pages = {380.e9 -- 380.e16}, year = {2019}, language = {en} } @inproceedings{IravaniConrad2019, author = {Iravani, Sahar and Conrad, Tim}, title = {Deep Learning for Proteomics Data for Feature Selection and Classification}, volume = {11713}, booktitle = {Machine Learning and Knowledge Extraction. CD-MAKE 2019}, editor = {Holzinger, A. and Kieseberg, P. and Tjoa, A. and Weippl, E.}, publisher = {Springer, Cham}, doi = {10.1007/978-3-030-29726-8_19}, year = {2019}, language = {en} } @inproceedings{NavaYazdaniHegevonTycowicz2019, author = {Nava-Yazdani, Esfandiar and Hege, Hans-Christian and von Tycowicz, Christoph}, title = {A Geodesic Mixed Effects Model in Kendall's Shape Space}, volume = {11846}, booktitle = {Proc. 7th MICCAI workshop on Mathematical Foundations of Computational Anatomy (MFCA)}, doi = {10.1007/978-3-030-33226-6_22}, pages = {209 -- 218}, year = {2019}, abstract = {In many applications, geodesic hierarchical models are adequate for the study of temporal observations. We employ such a model derived for manifold-valued data to Kendall's shape space. In particular, instead of the Sasaki metric, we adapt a functional-based metric, which increases the computational efficiency and does not require the implementation of the curvature tensor. We propose the corresponding variational time discretization of geodesics and apply the approach for the estimation of group trends and statistical testing of 3D shapes derived from an open access longitudinal imaging study on osteoarthritis.}, language = {en} } @article{BanyassadyChiuKormanetal.2019, author = {Banyassady, Bahareh and Chiu, Man-Kwun and Korman, Matias and Mulzer, Wolfgang and van Renssen, Andr{\´e} and Roeloffzen, Marcel and Seiferth, Paul and Stein, Yannik and Vogtenhuber, Birgit and Willert, Max}, title = {Routing in polygonal domains}, volume = {87}, journal = {Computational Geometry, Theory and Applications}, publisher = {Elsevier}, doi = {10.1016/j.comgeo.2019.101593}, year = {2019}, abstract = {We consider the problem of routing a data packet through the visibility graph of a polygonal domain P with n vertices and h holes. We may preprocess P to obtain a "label" and a "routing table" for each vertex of P. Then, we must be able to route a data packet between any two vertices p and q of P, where each step must use only the label of the target node q and the routing table of the current node. For any fixed epsilon > 0, we present a routing scheme that always achieves a routing path whose length exceeds the shortest path by a factor of at most 1 + epsilon. The labels have O(log n) bits, and the routing tables are of size O(((epsilon^-1)+h)log n). The preprocessing time is O((n^2)log n). It can be improved to O(n^2) for simple polygons.}, language = {en} } @misc{BuchmannKaplanPowelletal.2019, author = {Buchmann, Jens and Kaplan, Bernhard and Powell, Samuel and Prohaska, Steffen and Laufer, Jan}, title = {Quantitative PA tomography of high resolution 3-D images: experimental validation in tissue phantoms}, issn = {1438-0064}, doi = {10.1016/j.pacs.2019.100157}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-75416}, year = {2019}, abstract = {Quantitative photoacoustic tomography aims recover the spatial distribution of absolute chromophore concentrations and their ratios from deep tissue, high-resolution images. In this study, a model-based inversion scheme based on a Monte-Carlo light transport model is experimentally validated on 3-D multispectral images of a tissue phantom acquired using an all-optical scanner with a planar detection geometry. A calibrated absorber allowed scaling of the measured data during the inversion, while an acoustic correction method was employed to compensate the effects of limited view detection. Chromophore- and fluence-dependent step sizes and Adam optimization were implemented to achieve rapid convergence. High resolution 3-D maps of absolute concentrations and their ratios were recovered with high accuracy. Potential applications of this method include quantitative functional and molecular photoacoustic tomography of deep tissue in preclinical and clinical studies.}, language = {en} } @misc{LindowBruenigDercksenetal.2019, author = {Lindow, Norbert and Br{\"u}nig, Florian and Dercksen, Vincent J. and Fabig, Gunar and Kiewisz, Robert and Redemann, Stefanie and M{\"u}ller-Reichert, Thomas and Prohaska, Steffen}, title = {Semi-automatic Stitching of Serial Section Image Stacks with Filamentous Structures}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-73739}, year = {2019}, abstract = {In this paper, we present a software-assisted workflow for the alignment and matching of filamentous structures across a stack of 3D serial image sections. This is achieved by a combination of automatic methods, visual validation, and interactive correction. After an initial alignment, the user can continuously improve the result by interactively correcting landmarks or matches of filaments. This is supported by a quality assessment that visualizes regions that have been already inspected and, thus, allows a trade-off between quality and manual labor. The software tool was developed in collaboration with biologists who investigate microtubule-based spindles during cell division. To quantitatively understand the structural organization of such spindles, a 3D reconstruction of the numerous microtubules is essential. Each spindle is cut into a series of semi-thick physical sections, of which electron tomograms are acquired. The sections then need to be stitched, i.e. non-rigidly aligned; and the microtubules need to be traced in each section and connected across section boundaries. Experiments led to the conclusion that automatic methods for stitching alone provide only an incomplete solution to practical analysis needs. Automatic methods may fail due to large physical distortions, a low signal-to-noise ratio of the images, or other unexpected experimental difficulties. In such situations, semi-automatic validation and correction is required to rescue as much information as possible to derive biologically meaningful results despite of some errors related to data collection. Since the correct stitching is visually not obvious due to the number of microtubules (up to 30k) and their dense spatial arrangement, these are difficult tasks. Furthermore, a naive inspection of each microtubule is too time consuming. In addition, interactive visualization is hampered by the size of the image data (up to 100 GB). Based on the requirements of our collaborators, we present a practical solution for the semi-automatic stitching of serial section image stacks with filamentous structures.}, language = {en} } @article{NyakaturaBaumgartenBaumetal.2019, author = {Nyakatura, John and Baumgarten, Roxane and Baum, Daniel and Stark, Heiko and Youlatos, Dionisios}, title = {Muscle internal structure revealed by contrast-enhanced μCT and fibre recognition: The hindlimb extensors of an arboreal and a fossorial squirrel}, volume = {99}, journal = {Mammalian Biology}, doi = {10.1016/j.mambio.2019.10.007}, pages = {71 -- 80}, year = {2019}, abstract = {In individuals of similar body mass representing closely related species with different lifestyles, muscle architectural properties can be assumed to reflect adaptation to differing, lifestyle-related functional demands. We here employ a fiber recognition algorithm on contrast-enhanced micro-computed tomography (μCT) scans of one specimen each of an arboreal (Sciurus vulgaris) and a fossorial (Spermophilus citellus) sciuromorph rodent. The automated approach accounts for potential heterogeneity of architectural properties within a muscle by analyzing all fascicles that compose a muscle. Muscle architectural properties (volume, fascicle length, and orientation, and force-generating capacity) were quantified in 14 hindlimb (hip, knee, and ankle) extensor muscles and compared between specimens. We expected the arboreal squirrel to exhibit greater force-generating capacity and a greater capacity for length change allowing more powerful hindlimb extension. Generally and mostly matching our expectations, the S. vulgaris specimen had absolutely and relatively larger extensor muscles than the S. citellus specimen which were thus metabolically more expensive and demonstrate the relatively larger investment into powerful hindlimb extension necessary in the arboreal context. We conclude that detailed quantitative data on hindlimb muscle internal structure as was gathered here for a very limited sample further lends support to the notion that muscle architecture reflects adaptation to differential functional demands in closely related species with different locomotor behaviors and lifestyles.}, language = {en} } @misc{NavaYazdaniHegevonTycowicz2019, author = {Nava-Yazdani, Esfandiar and Hege, Hans-Christian and von Tycowicz, Christoph}, title = {A Geodesic Mixed Effects Model in Kendall's Shape Space}, journal = {Proc. 7th MICCAI workshop on Mathematical Foundations of Computational Anatomy (MFCA)}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-74621}, year = {2019}, abstract = {In many applications, geodesic hierarchical models are adequate for the study of temporal observations. We employ such a model derived for manifold-valued data to Kendall's shape space. In particular, instead of the Sasaki metric, we adapt a functional-based metric, which increases the computational efficiency and does not require the implementation of the curvature tensor. We propose the corresponding variational time discretization of geodesics and apply the approach for the estimation of group trends and statistical testing of 3D shapes derived from an open access longitudinal imaging study on osteoarthritis.}, language = {en} } @article{FroehlerElberfeldMoelleretal.2019, author = {Fr{\"o}hler, Bernhard and Elberfeld, Tim and M{\"o}ller, Torsten and Hege, Hans-Christian and Weissenb{\"o}ck, Johannes and De Beenhouwer, Jan and Sijbers, Jan and Kastner, Johann and Heinzl, Christoph}, title = {A Visual Tool for the Analysis of Algorithms for Tomographic Fiber Reconstruction in Materials Science}, volume = {38}, journal = {Computer Graphics Forum}, number = {3}, doi = {10.1111/cgf.13688}, pages = {273 -- 283}, year = {2019}, abstract = {We present visual analysis methods for the evaluation of tomographic fiber reconstruction algorithms by means of analysis, visual debugging and comparison of reconstructed fibers in materials science. The methods are integrated in a tool (FIAKER) that supports the entire workflow. It enables the analysis of various fiber reconstruction algorithms, of differently parameterized fiber reconstruction algorithms and of individual steps in iterative fiber reconstruction algorithms. Insight into the performance of fiber reconstruction algorithms is obtained by a list-based ranking interface. A 3D view offers interactive visualization techniques to gain deeper insight, e.g., into the aggregated quality of the examined fiber reconstruction algorithms and parameterizations. The tool was designed in close collaboration with researchers who work with fiber-reinforced polymers on a daily basis and develop algorithms for tomographic reconstruction and characterization of such materials. We evaluate the tool using synthetic datasets as well as tomograms of real materials. Five case studies certify the usefulness of the tool, showing that it significantly accelerates the analysis and provides valuable insights that make it possible to improve the fiber reconstruction algorithms. The main contribution of the paper is the well-considered combination of methods and their seamless integration into a visual tool that supports the entire workflow. Further findings result from the analysis of (dis-)similarity measures for fibers as well as from the discussion of design decisions. It is also shown that the generality of the analytical methods allows a wider range of applications, such as the application in pore space analysis.}, language = {en} } @article{GoubergritsHellmeierBrueningetal.2019, author = {Goubergrits, Leonid and Hellmeier, Florian and Bruening, Jan Joris and Spuler, Andreas and Hege, Hans-Christian and Voss, Samuel and Janiga, G{\´a}bor and Saalfeld, Sylvia and Beuing, Oliver and Berg, Philipp}, title = {Multiple Aneurysms AnaTomy CHallenge 2018 (MATCH): Uncertainty Quantification of Geometric Rupture Risk Parameters}, volume = {18}, journal = {BioMedical Engineering OnLine}, number = {35}, doi = {10.1186/s12938-019-0657-y}, year = {2019}, abstract = {Background Geometric parameters have been proposed for prediction of cerebral aneurysm rupture risk. Predicting the rupture risk for incidentally detected unruptured aneurysms could help clinicians in their treatment decision. However, assessment of geometric parameters depends on several factors, including the spatial resolution of the imaging modality used and the chosen reconstruction procedure. The aim of this study was to investigate the uncertainty of a variety of previously proposed geometric parameters for rupture risk assessment, caused by variability of reconstruction procedures. Materials 26 research groups provided segmentations and surface reconstructions of five cerebral aneurysms as part of the Multiple Aneurysms AnaTomy CHallenge (MATCH) 2018. 40 dimensional and non-dimensional geometric parameters, describing aneurysm size, neck size, and irregularity of aneurysm shape, were computed. The medians as well as the absolute and relative uncertainties of the parameters were calculated. Additionally, linear regression analysis was performed on the absolute uncertainties and the median parameter values. Results A large variability of relative uncertainties in the range between 3.9 and 179.8\% was found. Linear regression analysis indicates that some parameters capture similar geometric aspects. The lowest uncertainties < 6\% were found for the non-dimensional parameters isoperimetric ratio, convexity ratio, and ellipticity index. Uncertainty of 2D and 3D size parameters was significantly higher than uncertainty of 1D parameters. The most extreme uncertainties > 80\% were found for some curvature parameters. Conclusions Uncertainty analysis is essential on the road to clinical translation and use of rupture risk prediction models. Uncertainty quantification of geometric rupture risk parameters provided by this study may help support development of future rupture risk prediction models.}, language = {en} } @misc{AmbellanZachowvonTycowicz2019, author = {Ambellan, Felix and Zachow, Stefan and von Tycowicz, Christoph}, title = {An as-invariant-as-possible GL+(3)-based Statistical Shape Model}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-74566}, year = {2019}, abstract = {We describe a novel nonlinear statistical shape model basedon differential coordinates viewed as elements of GL+(3). We adopt an as-invariant-as possible framework comprising a bi-invariant Lie group mean and a tangent principal component analysis based on a unique GL+(3)-left-invariant, O(3)-right-invariant metric. Contrary to earlier work that equips the coordinates with a specifically constructed group structure, our method employs the inherent geometric structure of the group-valued data and therefore features an improved statistical power in identifying shape differences. We demonstrate this in experiments on two anatomical datasets including comparison to the standard Euclidean as well as recent state-of-the-art nonlinear approaches to statistical shape modeling.}, language = {en} } @incollection{Baum2019, author = {Baum, Daniel}, title = {An Evaluation of Color Maps for Visual Data Exploration}, booktitle = {Science in Color: Visualizing Achromatic Knowledge}, editor = {Bock von W{\"u}lfingen, Bettina}, publisher = {De Gruyter}, address = {Berlin}, pages = {147 -- 161}, year = {2019}, language = {en} } @misc{AmbellanZachowvonTycowicz2019, author = {Ambellan, Felix and Zachow, Stefan and von Tycowicz, Christoph}, title = {A Surface-Theoretic Approach for Statistical Shape Modeling}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-74497}, year = {2019}, abstract = {We present a novel approach for nonlinear statistical shape modeling that is invariant under Euclidean motion and thus alignment-free. By analyzing metric distortion and curvature of shapes as elements of Lie groups in a consistent Riemannian setting, we construct a framework that reliably handles large deformations. Due to the explicit character of Lie group operations, our non-Euclidean method is very efficient allowing for fast and numerically robust processing. This facilitates Riemannian analysis of large shape populations accessible through longitudinal and multi-site imaging studies providing increased statistical power. We evaluate the performance of our model w.r.t. shape-based classification of pathological malformations of the human knee and show that it outperforms the standard Euclidean as well as a recent nonlinear approach especially in presence of sparse training data. To provide insight into the model's ability of capturing natural biological shape variability, we carry out an analysis of specificity and generalization ability.}, language = {en} } @inproceedings{AmbellanZachowvonTycowicz2019, author = {Ambellan, Felix and Zachow, Stefan and von Tycowicz, Christoph}, title = {An as-invariant-as-possible GL+(3)-based Statistical Shape Model}, volume = {11846}, booktitle = {Proc. 7th MICCAI workshop on Mathematical Foundations of Computational Anatomy (MFCA)}, publisher = {Springer}, doi = {10.1007/978-3-030-33226-6_23}, pages = {219 -- 228}, year = {2019}, abstract = {We describe a novel nonlinear statistical shape model basedon differential coordinates viewed as elements of GL+(3). We adopt an as-invariant-as possible framework comprising a bi-invariant Lie group mean and a tangent principal component analysis based on a unique GL+(3)-left-invariant, O(3)-right-invariant metric. Contrary to earlier work that equips the coordinates with a specifically constructed group structure, our method employs the inherent geometric structure of the group-valued data and therefore features an improved statistical power in identifying shape differences. We demonstrate this in experiments on two anatomical datasets including comparison to the standard Euclidean as well as recent state-of-the-art nonlinear approaches to statistical shape modeling.}, language = {en} } @misc{OezelKulkarniHasanetal.2019, author = {{\"O}zel, M. Neset and Kulkarni, Abhishek and Hasan, Amr and Brummer, Josephine and Moldenhauer, Marian and Daumann, Ilsa-Maria and Wolfenberg, Heike and Dercksen, Vincent J. and Kiral, F. Ridvan and Weiser, Martin and Prohaska, Steffen and von Kleist, Max and Hiesinger, Peter Robin}, title = {Serial synapse formation through filopodial competition for synaptic seeding factors}, issn = {1438-0064}, doi = {10.1016/j.devcel.2019.06.014}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-74397}, year = {2019}, abstract = {Following axon pathfinding, growth cones transition from stochastic filopodial exploration to the formation of a limited number of synapses. How the interplay of filopodia and synapse assembly ensures robust connectivity in the brain has remained a challenging problem. Here, we developed a new 4D analysis method for filopodial dynamics and a data-driven computational model of synapse formation for R7 photoreceptor axons in developing Drosophila brains. Our live data support a 'serial synapse formation' model, where at any time point only a single 'synaptogenic' filopodium suppresses the synaptic competence of other filopodia through competition for synaptic seeding factors. Loss of the synaptic seeding factors Syd-1 and Liprin-α leads to a loss of this suppression, filopodial destabilization and reduced synapse formation, which is sufficient to cause the destabilization of entire axon terminals. Our model provides a filopodial 'winner-takes-all' mechanism that ensures the formation of an appropriate number of synapses.}, language = {en} } @article{KraemerMaggioniBrissonetal.2019, author = {Kr{\"a}mer, Martin and Maggioni, Marta and Brisson, Nicholas and Zachow, Stefan and Teichgr{\"a}ber, Ulf and Duda, Georg and Reichenbach, J{\"u}rgen}, title = {T1 and T2* mapping of the human quadriceps and patellar tendons using ultra-short echo-time (UTE) imaging and bivariate relaxation parameter-based volumetric visualization}, volume = {63}, journal = {Magnetic Resonance Imaging}, number = {11}, doi = {10.1016/j.mri.2019.07.015}, pages = {29 -- 36}, year = {2019}, abstract = {Quantification of magnetic resonance (MR)-based relaxation parameters of tendons and ligaments is challenging due to their very short transverse relaxation times, requiring application of ultra-short echo-time (UTE) imaging sequences. We quantify both T1 and T2⁎ in the quadriceps and patellar tendons of healthy volunteers at a field strength of 3 T and visualize the results based on 3D segmentation by using bivariate histogram analysis. We applied a 3D ultra-short echo-time imaging sequence with either variable repetition times (VTR) or variable flip angles (VFA) for T1 quantification in combination with multi-echo acquisition for extracting T2⁎. The values of both relaxation parameters were subsequently binned for bivariate histogram analysis and corresponding cluster identification, which were subsequently visualized. Based on manually-drawn regions of interest in the tendons on the relaxation parameter maps, T1 and T2⁎ boundaries were selected in the bivariate histogram to segment the quadriceps and patellar tendons and visualize the relaxation times by 3D volumetric rendering. Segmentation of bone marrow, fat, muscle and tendons was successfully performed based on the bivariate histogram analysis. Based on the segmentation results mean T2⁎ relaxation times, over the entire tendon volumes averaged over all subjects, were 1.8 ms ± 0.1 ms and 1.4 ms ± 0.2 ms for the patellar and quadriceps tendons, respectively. The mean T1 value of the patellar tendon, averaged over all subjects, was 527 ms ± 42 ms and 476 ms ± 40 ms for the VFA and VTR acquisitions, respectively. The quadriceps tendon had higher mean T1 values of 662 ms ± 97 ms (VFA method) and 637 ms ± 40 ms (VTR method) compared to the patellar tendon. 3D volumetric visualization of the relaxation times revealed that T1 values are not constant over the volume of both tendons, but vary locally. This work provided additional data to build upon the scarce literature available on relaxation times in the quadriceps and patellar tendons. We were able to segment both tendons and to visualize the relaxation parameter distributions over the entire tendon volumes.}, language = {en} } @inproceedings{NeumannHellwichZachow2019, author = {Neumann, Mario and Hellwich, Olaf and Zachow, Stefan}, title = {Localization and Classification of Teeth in Cone Beam CT using Convolutional Neural Networks}, booktitle = {Proc. of the 18th annual conference on Computer- and Robot-assisted Surgery (CURAC)}, isbn = {978-3-00-063717-9}, pages = {182 -- 188}, year = {2019}, abstract = {In dentistry, software-based medical image analysis and visualization provide efficient and accurate diagnostic and therapy planning capabilities. We present an approach for the automatic recognition of tooth types and positions in digital volume tomography (DVT). By using deep learning techniques in combination with dimensionality reduction through non-planar reformatting of the jaw anatomy, DVT data can be efficiently processed and teeth reliably recognized and classified, even in the presence of imaging artefacts, missing or dislocated teeth. We evaluated our approach, which is based on 2D Convolutional Neural Networks (CNNs), on 118 manually annotated cases of clinical DVT datasets. Our proposed method correctly classifies teeth with an accuracy of 94\% within a limit of 2mm distance to ground truth labels.}, language = {en} } @inproceedings{JoachimskyMaIckingetal.2019, author = {Joachimsky, Robert and Ma, Lihong and Icking, Christian and Zachow, Stefan}, title = {A Collision-Aware Articulated Statistical Shape Model of the Human Spine}, booktitle = {Proc. of the 18th annual conference on Computer- and Robot-assisted Surgery (CURAC)}, pages = {58 -- 64}, year = {2019}, abstract = {Statistical Shape Models (SSMs) are a proven means for model-based 3D anatomy reconstruction from medical image data. In orthopaedics and biomechanics, SSMs are increasingly employed to individualize measurement data or to create individualized anatomical models to which implants can be adapted to or functional tests can be performed on. For modeling and analysis of articulated structures, so called articulated SSMs (aSSMs) have been developed. However, a missing feature of aSSMs is the consideration of collisions in the course of individual fitting and articulation. The aim of our work was to develop aSSMs that handle collisions between components correctly. That way it becomes possible to adjust shape and articulation in view of a physically and geometrically plausible individualization. To be able to apply collision-aware aSSMs in simulation and optimisation, our approach is based on an e� cient collision detection method employing Graphics Processing Units (GPUs).}, language = {en} } @misc{Joachimsky2019, type = {Master Thesis}, author = {Joachimsky, Robert}, title = {Approaching Spinal Kinematics using a Collision-Aware Articulated Deformable Model}, pages = {84}, year = {2019}, abstract = {Statistical Shape Models (SSMs) allow for a compact representation of shape and shape variation and they are a proven means for model-based 3D anatomy reconstruction from medical image data. In orthopaedics and biomechanics, SSMs are increasingly employed to individualize measurement data or to create individualized anatomical models. The human spine is a versatile and complex articulated structure and thus is an interesting candidate to be modeled using an advanced type of SSMs. For modeling and analysis of articulated structures, so called articulated SSMs (aSSMs) have been developed. However, a missing feature of aSSMs is the consideration of collisions in the course of individual fitting and articulation. The aim of this thesis is to develop an aSSM of two adjacent vertebrae that handles collisions between components correctly. The model will incorporate the two major aspects of variability: Shape of a single vertebra and the relative positioning of neighboring vertebrae. That way it becomes possible to adjust shape and articulation in view of a physically and geometrically plausible individualization. To be able to apply collision-aware aSSMs in simulation and optimisation in future work, the approach is based on a parallelized collision detection method employing Graphics Processing Units (GPUs).}, language = {en} } @misc{MahnkeArltBaumetal.2019, author = {Mahnke, Heinz-Eberhard and Arlt, Tobias and Baum, Daniel and Hege, Hans-Christian and Herter, Felix and Lindow, Norbert and Manke, Ingo and Siopi, Tzulia and Menei, Eve and Etienne, Marc and Lepper, Verena}, title = {Virtual unfolding of folded papyri}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-74338}, year = {2019}, abstract = {The historical importance of ancient manuscripts is unique since they provide information about the heritage of ancient cultures. Often texts are hidden in rolled or folded documents. Due to recent impro- vements in sensitivity and resolution, spectacular disclosures of rolled hidden texts were possible by X-ray tomography. However, revealing text on folded manuscripts is even more challenging. Manual unfolding is often too risky in view of the fragile condition of fragments, as it can lead to the total loss of the document. X-ray tomography allows for virtual unfolding and enables non-destructive access to hid- den texts. We have recently demonstrated the procedure and tested unfolding algorithms on a mockup sample. Here, we present results on unfolding ancient papyrus packages from the papyrus collection of the Mus{\´e}e du Louvre, among them objects folded along approximately orthogonal folding lines. In one of the packages, the first identification of a word was achieved, the Coptic word for "Lord".}, language = {en} } @misc{Baum2019, author = {Baum, Daniel}, title = {An Evaluation of Color Maps for Visual Data Exploration}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-74259}, year = {2019}, language = {en} } @article{BaumWeaverZlotnikovetal.2019, author = {Baum, Daniel and Weaver, James C. and Zlotnikov, Igor and Kn{\"o}tel, David and Tomholt, Lara and Dean, Mason N.}, title = {High-Throughput Segmentation of Tiled Biological Structures using Random-Walk Distance Transforms}, journal = {Integrative And Comparative Biology}, doi = {10.1093/icb/icz117}, year = {2019}, abstract = {Various 3D imaging techniques are routinely used to examine biological materials, the results of which are usually a stack of grayscale images. In order to quantify structural aspects of the biological materials, however, they must first be extracted from the dataset in a process called segmentation. If the individual structures to be extracted are in contact or very close to each other, distance-based segmentation methods utilizing the Euclidean distance transform are commonly employed. Major disadvantages of the Euclidean distance transform, however, are its susceptibility to noise (very common in biological data), which often leads to incorrect segmentations (i.e. poor separation of objects of interest), and its limitation of being only effective for roundish objects. In the present work, we propose an alternative distance transform method, the random-walk distance transform, and demonstrate its effectiveness in high-throughput segmentation of three microCT datasets of biological tilings (i.e. structures composed of a large number of similar repeating units). In contrast to the Euclidean distance transform, this random-walk approach represents the global, rather than the local, geometric character of the objects to be segmented and, thus, is less susceptible to noise. In addition, it is directly applicable to structures with anisotropic shape characteristics. Using three case studies—stingray tessellated cartilage, starfish dermal endoskeleton, and the prismatic layer of bivalve mollusc shell—we provide a typical workflow for the segmentation of tiled structures, describe core image processing concepts that are underused in biological research, and show that for each study system, large amounts of biologically-relevant data can be rapidly segmented, visualized and analyzed.}, language = {en} } @inproceedings{Paetsch2019, author = {Paetsch, Olaf}, title = {Possibilities and Limitations of Automatic Feature Extraction shown by the Example of Crack Detection in 3D-CT Images of Concrete Specimen}, booktitle = {iCT 2019}, year = {2019}, abstract = {To assess the influence of the alkali-silica reaction (ASR) on pavement concrete 3D-CT imaging has been applied to concrete samples. Prior to imaging these samples have been drilled out of a concrete beam pre-damaged by fatigue loading. The resulting high resolution 3D-CT images consist of several gigabytes of voxels. Current desktop computers can visualize such big datasets without problems but a visual inspection or manual segmentation of features such as cracks by experts can only be carried out on a few slices. A quantitative analysis of cracks requires a segmentation of the whole specimen which could only be done by an automatic feature detection. This arises the question of the reliability of an automatic crack detection algorithm, its certainty and limitations. Does the algorithm find all cracks? Does it find too many cracks? Can parameters of that algorithm, once identified as good, be applied to other samples as well? Can ensemble computing with many crack parameters overcome the difficulties with parameter finding? By means of a crack detection algorithm based on shape recognition (template matching) these questions will be discussed. Since the author has no access to reliable ground truth data of cracks the assessment of the certainty of the automatic crack is restricted to visual inspection by experts. Therefore, an artificial dataset based on a combination of manually segmented cracks processed together with simple image processing algorithms is used to quantify the accuracy of the crack detection algorithm. Part of the evaluation of cracks in concrete samples is the knowledge of the surrounding material. The surrounding material can be used to assess the detected cracks, e.g. micro-cracks within the aggregate-matrix interface may be starting points for cracks on a macro scale. Furthermore, the knowledge of the surrounding material can help to find better parameter sets for the crack detection itself because crack characteristics may vary depending on their surrounding material. Therefore, in addition to a crack detection a complete segmentation of the sample into the components of concrete, such as aggregates, cement matrix and pores is needed. Since such a segmentation task cannot be done manually due to the amount of data, an approach utilizing convolutional neuronal networks stemming from a medical application has been applied. The learning phase requires a ground truth i.e. a segmentation of the components. This has to be created manually in a time-consuming task. However, this segmentation can be used for a quantitative evaluation of the automatic segmentation afterwards. Even though that work has been performed as a short term subtask of a bigger project funded by the German Research Foundation (DFG) this paper discusses problems which may arise in similar projects, too.}, language = {en} } @misc{Neumann2019, type = {Master Thesis}, author = {Neumann, Mario}, title = {Localization and Classification of Teeth in Cone Beam Computed Tomography using 2D CNNs}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-74045}, pages = {77}, year = {2019}, abstract = {In dentistry, software-based medical image analysis and visualization provide effcient and accurate diagnostic and therapy planning capabilities. We present an approach for the automatic recognition of tooth types and positions in digital volume tomography (DVT). By using deep learning techniques in combination with dimension reduction through non-planar reformatting of the jaw anatomy, DVT data can be effciently processed and teeth reliably recognized and classified, even in the presence of imaging artefacts, missing or dislocated teeth. We evaluated our approach, which is based on 2D Convolutional Neural Networks (CNNs), on 118 manually annotated cases of clinical DVT datasets. Our proposed method correctly classifies teeth with an accuracy of 94\% within a limit of 2mm distancr to ground truth landmarks.}, language = {en} } @inproceedings{FroehlerdaCunhaMeloWeissenboecketal.2019, author = {Fr{\"o}hler, Bernhard and da Cunha Melo, Lucas and Weissenb{\"o}ck, Johannes and Kastner, Johann and M{\"o}ller, Torsten and Hege, Hans-Christian and Gr{\"o}ller, Eduard M. and Sanctorum, Jonathan and De Beenhouwer, Jan and Sijbers, Jan and Heinzl, Christoph}, title = {Tools for the analysis of datasets from X-ray computed tomography based on Talbot-Lau grating interferometry}, booktitle = {Proceedings of iCT 2019, (9th Conference on Industrial Computed Tomography, Padova, Italy - iCT 2019, February 13-15, 2019)}, number = {paper 52}, pages = {8}, year = {2019}, abstract = {This work introduces methods for analyzing the three imaging modalities delivered by Talbot-Lau grating interferometry X-ray computed tomography (TLGI-XCT). The first problem we address is providing a quick way to show a fusion of all three modal- ities. For this purpose the tri-modal transfer function widget is introduced. The widget controls a mixing function that uses the output of the transfer functions of all three modalities, allowing the user to create one customized fused image. A second problem prevalent in processing TLGI-XCT data is a lack of tools for analyzing the segmentation process of such multimodal data. We address this by providing methods for computing three types of uncertainty: From probabilistic segmentation algorithms, from the voxel neighborhoods as well as from a collection of results. We furthermore introduce a linked views interface to explore this data. The techniques are evaluated on a TLGI-XCT scan of a carbon-fiber reinforced dataset with impact damage. We show that the transfer function widget accelerates and facilitates the exploration of this dataset, while the uncertainty analysis methods give insights into how to tweak and improve segmentation algorithms for more suitable results.}, language = {en} } @misc{BaumWeaverZlotnikovetal.2019, author = {Baum, Daniel and Weaver, James C. and Zlotnikov, Igor and Kn{\"o}tel, David and Tomholt, Lara and Dean, Mason N.}, title = {High-Throughput Segmentation of Tiled Biological Structures using Random-Walk Distance Transforms}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-73841}, year = {2019}, abstract = {Various 3D imaging techniques are routinely used to examine biological materials, the results of which are usually a stack of grayscale images. In order to quantify structural aspects of the biological materials, however, they must first be extracted from the dataset in a process called segmentation. If the individual structures to be extracted are in contact or very close to each other, distance-based segmentation methods utilizing the Euclidean distance transform are commonly employed. Major disadvantages of the Euclidean distance transform, however, are its susceptibility to noise (very common in biological data), which often leads to incorrect segmentations (i.e. poor separation of objects of interest), and its limitation of being only effective for roundish objects. In the present work, we propose an alternative distance transform method, the random-walk distance transform, and demonstrate its effectiveness in high-throughput segmentation of three microCT datasets of biological tilings (i.e. structures composed of a large number of similar repeating units). In contrast to the Euclidean distance transform, this random-walk approach represents the global, rather than the local, geometric character of the objects to be segmented and, thus, is less susceptible to noise. In addition, it is directly applicable to structures with anisotropic shape characteristics. Using three case studies—stingray tessellated cartilage, starfish dermal endoskeleton, and the prismatic layer of bivalve mollusc shell—we provide a typical workflow for the segmentation of tiled structures, describe core image processing concepts that are underused in biological research, and show that for each study system, large amounts of biologically-relevant data can be rapidly segmented, visualized and analyzed.}, language = {en} } @inproceedings{AmbellanZachowvonTycowicz2019, author = {Ambellan, Felix and Zachow, Stefan and von Tycowicz, Christoph}, title = {A Surface-Theoretic Approach for Statistical Shape Modeling}, volume = {11767}, booktitle = {Proc. Medical Image Computing and Computer Assisted Intervention (MICCAI), Part IV}, publisher = {Springer}, doi = {10.1007/978-3-030-32251-9_3}, pages = {21 -- 29}, year = {2019}, abstract = {We present a novel approach for nonlinear statistical shape modeling that is invariant under Euclidean motion and thus alignment-free. By analyzing metric distortion and curvature of shapes as elements of Lie groups in a consistent Riemannian setting, we construct a framework that reliably handles large deformations. Due to the explicit character of Lie group operations, our non-Euclidean method is very efficient allowing for fast and numerically robust processing. This facilitates Riemannian analysis of large shape populations accessible through longitudinal and multi-site imaging studies providing increased statistical power. We evaluate the performance of our model w.r.t. shape-based classification of pathological malformations of the human knee and show that it outperforms the standard Euclidean as well as a recent nonlinear approach especially in presence of sparse training data. To provide insight into the model's ability of capturing natural biological shape variability, we carry out an analysis of specificity and generalization ability.}, language = {en} } @article{ZhangKlusConradetal.2019, author = {Zhang, Wei and Klus, Stefan and Conrad, Tim and Sch{\"u}tte, Christof}, title = {Learning chemical reaction networks from trajectory data}, volume = {18}, journal = {SIAM Journal on Applied Dynamical Systems (SIADS)}, number = {4}, arxiv = {http://arxiv.org/abs/1902.04920}, doi = {10.1137/19M1265880}, pages = {2000 -- 2046}, year = {2019}, abstract = {We develop a data-driven method to learn chemical reaction networks from trajectory data. Modeling the reaction system as a continuous-time Markov chain and assuming the system is fully observed,our method learns the propensity functions of the system with predetermined basis functions by maximizing the likelihood function of the trajectory data under l^1 sparse regularization. We demonstrate our method with numerical examples using synthetic data and carry out an asymptotic analysis of the proposed learning procedure in the infinite-data limit.}, language = {en} } @article{ThielDjurdjevacConradNtinietal.2019, author = {Thiel, Denise and Djurdjevac Conrad, Natasa and Ntini, Evgenia and Peschutter, Ria and Siebert, Heike and Marsico, Annalisa}, title = {Identifying lncRNA-mediated regulatory modules via ChIA-PET network analysis}, volume = {20}, journal = {BMC Bioinformatics}, number = {1471-2105}, doi = {10.1186/s12859-019-2900-8}, year = {2019}, abstract = {Background: Although several studies have provided insights into the role of long non-coding RNAs (lncRNAs), the majority of them have unknown function. Recent evidence has shown the importance of both lncRNAs and chromatin interactions in transcriptional regulation. Although network-based methods, mainly exploiting gene-lncRNA co-expression, have been applied to characterize lncRNA of unknown function by means of 'guilt-by-association', no strategy exists so far which identifies mRNA-lncRNA functional modules based on the 3D chromatin interaction graph. Results: To better understand the function of chromatin interactions in the context of lncRNA-mediated gene regulation, we have developed a multi-step graph analysis approach to examine the RNA polymerase II ChIA-PET chromatin interaction network in the K562 human cell line. We have annotated the network with gene and lncRNA coordinates, and chromatin states from the ENCODE project. We used centrality measures, as well as an adaptation of our previously developed Markov State Models (MSM) clustering method, to gain a better understanding of lncRNAs in transcriptional regulation. The novelty of our approach resides in the detection of fuzzy regulatory modules based on network properties and their optimization based on co-expression analysis between genes and gene-lncRNA pairs. This results in our method returning more bona fide regulatory modules than other state-of-the art approaches for clustering on graphs. Conclusions: Interestingly, we find that lncRNA network hubs tend to be significantly enriched in evolutionary conserved lncRNAs and enhancer-like functions. We validated regulatory functions for well known lncRNAs, such as MALAT1 and the enhancer-like lncRNA FALEC. In addition, by investigating the modular structure of bigger components we mine putative regulatory functions for uncharacterized lncRNAs.}, language = {en} } @incollection{AmbellanLameckervonTycowiczetal.2019, author = {Ambellan, Felix and Lamecker, Hans and von Tycowicz, Christoph and Zachow, Stefan}, title = {Statistical Shape Models - Understanding and Mastering Variation in Anatomy}, volume = {3}, booktitle = {Biomedical Visualisation}, number = {1156}, editor = {Rea, Paul M.}, edition = {1}, publisher = {Springer Nature Switzerland AG}, isbn = {978-3-030-19384-3}, doi = {10.1007/978-3-030-19385-0_5}, pages = {67 -- 84}, year = {2019}, abstract = {In our chapter we are describing how to reconstruct three-dimensional anatomy from medical image data and how to build Statistical 3D Shape Models out of many such reconstructions yielding a new kind of anatomy that not only allows quantitative analysis of anatomical variation but also a visual exploration and educational visualization. Future digital anatomy atlases will not only show a static (average) anatomy but also its normal or pathological variation in three or even four dimensions, hence, illustrating growth and/or disease progression. Statistical Shape Models (SSMs) are geometric models that describe a collection of semantically similar objects in a very compact way. SSMs represent an average shape of many three-dimensional objects as well as their variation in shape. The creation of SSMs requires a correspondence mapping, which can be achieved e.g. by parameterization with a respective sampling. If a corresponding parameterization over all shapes can be established, variation between individual shape characteristics can be mathematically investigated. We will explain what Statistical Shape Models are and how they are constructed. Extensions of Statistical Shape Models will be motivated for articulated coupled structures. In addition to shape also the appearance of objects will be integrated into the concept. Appearance is a visual feature independent of shape that depends on observers or imaging techniques. Typical appearances are for instance the color and intensity of a visual surface of an object under particular lighting conditions, or measurements of material properties with computed tomography (CT) or magnetic resonance imaging (MRI). A combination of (articulated) statistical shape models with statistical models of appearance lead to articulated Statistical Shape and Appearance Models (a-SSAMs).After giving various examples of SSMs for human organs, skeletal structures, faces, and bodies, we will shortly describe clinical applications where such models have been successfully employed. Statistical Shape Models are the foundation for the analysis of anatomical cohort data, where characteristic shapes are correlated to demographic or epidemiologic data. SSMs consisting of several thousands of objects offer, in combination with statistical methods ormachine learning techniques, the possibility to identify characteristic clusters, thus being the foundation for advanced diagnostic disease scoring.}, language = {en} } @misc{BuchmannKaplanPowelletal.2019, author = {Buchmann, Jens and Kaplan, Bernhard and Powell, Samuel and Prohaska, Steffen and Laufer, Jan}, title = {3D quantitative photoacoustic tomography using an adjoint radiance Monte Carlo model and gradient descent}, issn = {1438-0064}, doi = {10.1117/1.JBO.24.6.066001}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-72995}, year = {2019}, abstract = {Quantitative photoacoustic tomography aims to recover maps of the local concentrations of tissue chromophores from multispectral images. While model-based inversion schemes are promising approaches, major challenges to their practical implementation include the unknown fluence distribution and the scale of the inverse problem. This paper describes an inversion scheme based on a radiance Monte Carlo model and an adjoint-assisted gradient optimization that incorporates fluence-dependent step sizes and adaptive moment estimation. The inversion is shown to recover absolute chromophore concentrations, blood oxygen saturation and the Gr{\"u}neisen parameter from in silico 3D phantom images for different radiance approximations. The scattering coefficient was assumed to be homogeneous and known a priori.}, language = {en} } @article{BuchmannKaplanPowelletal.2019, author = {Buchmann, Jens and Kaplan, Bernhard and Powell, Samuel and Prohaska, Steffen and Laufer, Jan}, title = {3D quantitative photoacoustic tomography using an adjoint radiance Monte Carlo model and gradient descent}, volume = {24}, journal = {Journal of Biomedical Optics}, number = {6}, doi = {10.1117/1.JBO.24.6.066001}, pages = {066001}, year = {2019}, abstract = {Quantitative photoacoustic tomography aims to recover maps of the local concentrations of tissue chromophores from multispectral images. While model-based inversion schemes are promising approaches, major challenges to their practical implementation include the unknown fluence distribution and the scale of the inverse problem. This paper describes an inversion scheme based on a radiance Monte Carlo model and an adjoint-assisted gradient optimization that incorporates fluence-dependent step sizes and adaptive moment estimation. The inversion is shown to recover absolute chromophore concentrations, blood oxygen saturation and the Gr{\"u}neisen parameter from in silico 3D phantom images for different radiance approximations. The scattering coefficient was assumed to be homogeneous and known a priori.}, language = {en} } @misc{AmbellanTackEhlkeetal.2019, author = {Ambellan, Felix and Tack, Alexander and Ehlke, Moritz and Zachow, Stefan}, title = {Automated Segmentation of Knee Bone and Cartilage combining Statistical Shape Knowledge and Convolutional Neural Networks: Data from the Osteoarthritis Initiative}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-72704}, year = {2019}, abstract = {We present a method for the automated segmentation of knee bones and cartilage from magnetic resonance imaging (MRI) 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 shape models and neural networks employed are trained using data from the Osteoarthritis Initiative (OAI) and the MICCAI grand challenge "Segmentation of Knee Images 2010" (SKI10), respectively. We evaluate our method on 40 validation and 50 submission datasets from the SKI10 challenge.For the first time, an accuracy equivalent to the inter-observer variability of human readers is achieved in this challenge.Moreover, the quality of the proposed method is thoroughly assessed using various measures for data from the OAI, i.e. 507 manual segmentations of bone and cartilage, and 88 additional manual segmentations of cartilage. Our method yields sub-voxel accuracy for both OAI datasets. We make the 507 manual segmentations as well as our experimental setup publicly available to further aid research in the field of medical image segmentation.In conclusion, combining localized classification via CNNs with statistical anatomical knowledge via SSMs results in a state-of-the-art segmentation method for knee bones and cartilage from MRI data.}, language = {en} } @misc{AmbellanLameckervonTycowiczetal.2019, author = {Ambellan, Felix and Lamecker, Hans and von Tycowicz, Christoph and Zachow, Stefan}, title = {Statistical Shape Models - Understanding and Mastering Variation in Anatomy}, issn = {1438-0064}, doi = {10.1007/978-3-030-19385-0_5}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-72699}, year = {2019}, abstract = {In our chapter we are describing how to reconstruct three-dimensional anatomy from medical image data and how to build Statistical 3D Shape Models out of many such reconstructions yielding a new kind of anatomy that not only allows quantitative analysis of anatomical variation but also a visual exploration and educational visualization. Future digital anatomy atlases will not only show a static (average) anatomy but also its normal or pathological variation in three or even four dimensions, hence, illustrating growth and/or disease progression. Statistical Shape Models (SSMs) are geometric models that describe a collection of semantically similar objects in a very compact way. SSMs represent an average shape of many three-dimensional objects as well as their variation in shape. The creation of SSMs requires a correspondence mapping, which can be achieved e.g. by parameterization with a respective sampling. If a corresponding parameterization over all shapes can be established, variation between individual shape characteristics can be mathematically investigated. We will explain what Statistical Shape Models are and how they are constructed. Extensions of Statistical Shape Models will be motivated for articulated coupled structures. In addition to shape also the appearance of objects will be integrated into the concept. Appearance is a visual feature independent of shape that depends on observers or imaging techniques. Typical appearances are for instance the color and intensity of a visual surface of an object under particular lighting conditions, or measurements of material properties with computed tomography (CT) or magnetic resonance imaging (MRI). A combination of (articulated) statistical shape models with statistical models of appearance lead to articulated Statistical Shape and Appearance Models (a-SSAMs).After giving various examples of SSMs for human organs, skeletal structures, faces, and bodies, we will shortly describe clinical applications where such models have been successfully employed. Statistical Shape Models are the foundation for the analysis of anatomical cohort data, where characteristic shapes are correlated to demographic or epidemiologic data. SSMs consisting of several thousands of objects offer, in combination with statistical methods ormachine learning techniques, the possibility to identify characteristic clusters, thus being the foundation for advanced diagnostic disease scoring.}, language = {en} } @misc{SakuraiOnoCarretal.2019, author = {Sakurai, Daisuke and Ono, Kenji and Carr, Hamish and Nonaka, Jorji and Kawanabe, Tomohiro}, title = {Flexible Fiber Surfaces: A Reeb-Free Approach}, journal = {Topological Methods in Data Analysis and Visualization V}, editor = {Carr, Hamish and Fujishiro, Issei and Sadlo, Filip and Takahashi, Shigeo}, publisher = {Springer}, pages = {14}, year = {2019}, abstract = {The fiber surface generalizes the popular isosurface to multi-fields, so that pre-images can be visualized as surfaces. As with the isosurface, however, the fiber surface suffers from visual occlusion. We propose to avoid such occlusion by restricting the components to only the relevant ones with a new component-wise flexing algorithm. The approach, flexible fiber surface, generalizes the manipulation idea found in the flexible isosurface for the fiber surface. The flexible isosurface in the original form, however, relies on the contour tree. For the fiber surface, this corresponds to the Reeb space, which is challenging for both the computation and user interaction. We thus take a Reeb-free approach, in which one does not compute the Reeb space. Under this constraint, we generalize a few selected interactions in the flexible isosurface and discuss the implication of the restriction.}, language = {en} } @misc{SakuraiOnoCarretal.2019, author = {Sakurai, Daisuke and Ono, Kenji and Carr, Hamish and Nonaka, Jorji and Kawanabe, Tomohiro}, title = {Flexible Fiber Surfaces: A Reeb-Free Approach}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-72599}, year = {2019}, abstract = {The fiber surface generalizes the popular isosurface to multi-fields, so that pre-images can be visualized as surfaces. As with the isosurface, however, the fiber surface suffers from visual occlusion. We propose to avoid such occlusion by restricting the components to only the relevant ones with a new component-wise flexing algorithm. The approach, flexible fiber surface, generalizes the manipulation idea found in the flexible isosurface for the fiber surface. The flexible isosurface in the original form, however, relies on the contour tree. For the fiber surface, this corresponds to the Reeb space, which is challenging for both the computation and user interaction. We thus take a Reeb-free approach, in which one does not compute the Reeb space. Under this constraint, we generalize a few selected interactions in the flexible isosurface and discuss the implication of the restriction.}, language = {en} } @masterthesis{Prendke2019, type = {Bachelor Thesis}, author = {Prendke, Mona}, title = {Comparison of 2D and 3D CNNs for Classification of Knee MRI}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-72439}, pages = {53}, year = {2019}, language = {en} } @article{HildebrandtBrueningSchmidtetal.2019, author = {Hildebrandt, Thomas and Bruening, Jan Joris and Schmidt, Nora Laura and Lamecker, Hans and Heppt, Werner and Zachow, Stefan and Goubergrits, Leonid}, title = {The Healthy Nasal Cavity - Characteristics of Morphology and Related Airflow Based on a Statistical Shape Model Viewed from a Surgeon's Perspective}, volume = {35}, journal = {Facial Plastic Surgery}, number = {1}, doi = {10.1055/s-0039-1677721}, pages = {9 -- 13}, year = {2019}, abstract = {Functional surgery on the nasal framework requires referential criteria to objectively assess nasal breathing for indication and follow-up. Thismotivated us to generate amean geometry of the nasal cavity based on a statistical shape model. In this study, the authors could demonstrate that the introduced nasal cavity's mean geometry features characteristics of the inner shape and airflow, which are commonly observed in symptom-free subjects. Therefore, the mean geometry might serve as a reference-like model when one considers qualitative aspects. However, to facilitate quantitative considerations and statistical inference, further research is necessary. Additionally, the authorswere able to obtain details about the importance of the isthmus nasi and the inferior turbinate for the intranasal airstream.}, language = {en} } @article{HildebrandtBrueningLameckeretal.2019, author = {Hildebrandt, Thomas and Bruening, Jan Joris and Lamecker, Hans and Zachow, Stefan and Heppt, Werner and Schmidt, Nora and Goubergrits, Leonid}, title = {Digital Analysis of Nasal Airflow Facilitating Decision Support in Rhinosurgery}, volume = {35}, journal = {Facial Plastic Surgery}, number = {1}, doi = {10.1055/s-0039-1677720}, pages = {1 -- 8}, year = {2019}, abstract = {Successful functional surgery on the nasal framework requires reliable and comprehensive diagnosis. In this regard, the authors introduce a new methodology: Digital Analysis of Nasal Airflow (diANA). It is based on computational fluid dynamics, a statistical shape model of the healthy nasal cavity and rhinologic expertise. diANA necessitates an anonymized tomographic dataset of the paranasal sinuses including the complete nasal cavity and, when available, clinical information. The principle of diANA is to compare the morphology and the respective airflow of an individual nose with those of a reference. This enablesmorphometric aberrations and consecutive flow field anomalies to localize and quantify within a patient's nasal cavity. Finally, an elaborated expert opinion with instructive visualizations is provided. Using diANA might support surgeons in decision-making, avoiding unnecessary surgery, gaining more precision, and target-orientation for indicated operations.}, language = {en} } @article{LelievreZhang2019, author = {Leli{\`e}vre, Tony and Zhang, Wei}, title = {Pathwise estimates for effective dynamics: the case of nonlinear vectorial reaction coordinates}, journal = {Multiscale Modeling and Simulation}, number = {17}, arxiv = {http://arxiv.org/abs/1805.01928}, doi = {10.1137/18M1186034}, pages = {1019 -- 1051}, year = {2019}, abstract = {Effective dynamics using conditional expectation was proposed in [F. Legoll and T. Leli{\`e}vre, Nonlinearity, 2010] to approximate the essential dynamics of high-dimensional diffusion processes along a given reaction coordinate. The approximation error of the effective dynamics when it is used to approximate the behavior of the original dynamics has been considered in recent years. As a continuation of the previous work [F. Legoll, T. Leli{\`e}vre, and S. Olla, Stoch. Process. Appl, 2017], in this paper we obtain pathwise estimates for effective dynamics when the reaction coordinate function is either nonlinear or vector-valued.}, language = {en} } @article{Zhang2019, author = {Zhang, Wei}, title = {Ergodic SDEs on submanifolds and related numerical sampling schemes}, journal = {ESAIM: Mathematical Modelling and Numerical Analysis}, arxiv = {http://arxiv.org/abs/1702.08064}, year = {2019}, abstract = {In many applications, it is often necessary to sample the mean value of certain quantity with respect to a probability measure \$\mu\$ on the level set of a smooth function ξ:R^d→R^k, 1≤k