Publications of the computing group for ultrafast imaging
van de Kamp T., Schwermann A.H., dos Santos Rolo T., Losel P.D., Engler T., Etter W., Farago T., Gottlicher J., Heuveline V., Kopmann A., Mahler B., Mors T., Odar J., Rust J., Tan Jerome N., Vogelgesang M., Baumbach T., Krogmann L.
in Nature Communications, 9 (2018), 3325. DOI:10.1038/s41467-018-05654-y
© 2018, The Author(s). About 50% of all animal species are considered parasites. The linkage of species diversity to a parasitic lifestyle is especially evident in the insect order Hymenoptera. However, fossil evidence for host–parasitoid interactions is extremely rare, rendering hypotheses on the evolution of parasitism assumptive. Here, using high-throughput synchrotron X-ray microtomography, we examine 1510 phosphatized fly pupae from the Paleogene of France and identify 55 parasitation events by four wasp species, providing morphological and ecological data. All species developed as solitary endoparasitoids inside their hosts and exhibit different morphological adaptations for exploiting the same hosts in one habitat. Our results allow systematic and ecological placement of four distinct endoparasitoids in the Paleogene and highlight the need to investigate ecological data preserved in the fossil record.
Buzmakov A.V., Asadchikov V.E., Zolotov D.A., Roshchin B.S., Dymshits Y.M., Shishkov V.A., Chukalina M.V., Ingacheva A.S., Ichalova D.E., Krivonosov Y.S., Dyachkova I.G., Balzer M., Castele M., Chilingaryan S., Kopmann A.
in Crystallography Reports, 63 (2018) 1057-1061. DOI:10.1134/S106377451806007X
© 2018, Pleiades Publishing, Inc. Abstract: The design of a new automatic X-ray microtomograph is described. The parameters of optical schemes and X-ray detectors in use are presented. Methods for automating experiments, processing tomographic data, and getting access to them are reported.
Ametova E., Ferrucci M., Chilingaryan S., Dewulf W.
in Measurement Science and Technology, 29 (2018), 065007. DOI:10.1088/1361-6501/aab1a1
© 2018 IOP Publishing Ltd. The recent emergence of advanced manufacturing techniques such as additive manufacturing and an increased demand on the integrity of components have motivated research on the application of x-ray computed tomography (CT) for dimensional quality control. While CT has shown significant empirical potential for this purpose, there is a need for metrological research to accelerate the acceptance of CT as a measuring instrument. The accuracy in CT-based measurements is vulnerable to the instrument geometrical configuration during data acquisition, namely the relative position and orientation of x-ray source, rotation stage, and detector. Consistency between the actual instrument geometry and the corresponding parameters used in the reconstruction algorithm is critical. Currently available procedures provide users with only estimates of geometrical parameters. Quantification and propagation of uncertainty in the measured geometrical parameters must be considered to provide a complete uncertainty analysis and to establish confidence intervals for CT dimensional measurements. In this paper, we propose a computationally inexpensive model to approximate the influence of errors in CT geometrical parameters on dimensional measurement results. We use surface points extracted from a computer-aided design (CAD) model to model discrepancies in the radiographic image coordinates assigned to the projected edges between an aligned system and a system with misalignments. The efficacy of the proposed method was confirmed on simulated and experimental data in the presence of various geometrical uncertainty contributors.
Cavadini P., Weinhold H., Tonsmann M., Chilingaryan S., Kopmann A., Lewkowicz A., Miao C., Scharfer P., Schabel W.
in Experiments in Fluids, 59 (2018), 61. DOI:10.1007/s00348-017-2482-z
© 2018, Springer-Verlag GmbH Germany, part of Springer Nature. To understand the effects of inhomogeneous drying on the quality of polymer coatings, an experimental setup to resolve the occurring flow field throughout the drying film has been developed. Deconvolution microscopy is used to analyze the flow field in 3D and time. Since the dimension of the spatial component in the direction of the line-of-sight is limited compared to the lateral components, a multi-focal approach is used. Here, the beam of light is equally distributed on up to five cameras using cubic beam splitters. Adding a meniscus lens between each pair of camera and beam splitter and setting different distances between each camera and its meniscus lens creates multi-focality and allows one to increase the depth of the observed volume. Resolving the spatial component in the line-of-sight direction is based on analyzing the point spread function. The analysis of the PSF is computational expensive and introduces a high complexity compared to traditional particle image velocimetry approaches. A new algorithm tailored to the parallel computing architecture of recent graphics processing units has been developed. The algorithm is able to process typical images in less than a second and has further potential to realize online analysis in the future. As a prove of principle, the flow fields occurring in thin polymer solutions drying at ambient conditions and at boundary conditions that force inhomogeneous drying are presented.
Jerome N.T., Kopmann A.
in VISIGRAPP 2018 – Proceedings of the 13th International Joint Conference on Computer Vision, Imaging and Computer Graphics Theory and Applications, 3 (2018) 341-348.
© 2018 by SCITEPRESS – Science and Technology Publications, Lda. All rights reserved. With the advancement of instrument precision, research facilities are generating data at an unprecedented rate. These experimental results are stored in a digital library platform which the contents are later accessible from within the facility or the public. However, the sheer volume of collected data is overwhelming the capacity of researchers and impedes the process of browsing for the desired data. In this paper, we present a concept of Digital Visual Exploration Library (DVEL) based on the confluence of two major research domains-digital library and visualisation-that enables efficient browsing of the growing data within a digital library. We complement the current state-of-the-art textual metadata description by integrating visual exploration to address big complex data, i.e., data of large size, multimodal data and multivariate data. We describe our concept based on use cases from three unique domains: climate research with Doppler wind lidar, X-ray-imaging for entomology research, and medical imaging with ultrasound computer tomography.
Ametova E., Ferrucci M., Chilingaryan S., Dewulf W.
in Precision Engineering (2018). DOI:10.1016/j.precisioneng.2018.05.016
© 2018 Elsevier Inc. X-ray computed tomography (CT) is an imaging technique that allows the reconstruction of an imaged part in the form of a three-dimensional attenuation map. The CT data acquisition process consists of acquiring X-ray transmission images from multiple perspectives. Analysis of the reconstructed attenuation map can provide dimensional and material information about the measured part(s). Therefore, CT is recognized as a solution for quality control tasks, for example dimensional inspection of complex objects with intricate inner geometries. CT measurements can suffer from various sources of error in the measurement procedure. One such influence is the geometrical alignment of the CT instrument components. Typical tomographic reconstruction algorithms impose strict requirements on the relative position and orientation of the three main components: X-ray source, rotation axis of the sample stage, and X-ray detector. Any discrepancy in the actual CT geometry from the geometry assumed by the reconstruction algorithm will contribute to errors in measurements performed on the reconstructed data. There is currently no standardized or easily implementable method for users to compensate geometrical misalignments of the CT instrument. In many cases, the procedure of mechanical adjustment of CT instrument is time consuming and impractical. In this paper, we show that software-based compensation of deviations in CT instrument geometry is an effective alternative to mechanical adjustment of CT instrument. Through computer simulations, we compare qualitatively and quantitatively two methods to compensate CT instrument misalignment: radiographic re-binning (interpolation) and a modified conventional reconstruction algorithm with embedded misalignment compensation.
Farago T., Mikulik P., Ershov A., Vogelgesang M., Hanschke D., Baumbach T.
in Journal of Synchrotron Radiation, 24 (2017) 1283-1295. DOI:10.1107/S1600577517012255
© International Union of Crystallography, 2017. An open-source framework for conducting a broad range of virtual X-ray imaging experiments, syris, is presented. The simulated wavefield created by a source propagates through an arbitrary number of objects until it reaches a detector. The objects in the light path and the source are time-dependent, which enables simulations of dynamic experiments, e.g. four-dimensional time-resolved tomography and laminography. The high-level interface of syris is written in Python and its modularity makes the framework very flexible. The computationally demanding parts behind this interface are implemented in OpenCL, which enables fast calculations on modern graphics processing units. The combination of flexibility and speed opens new possibilities for studying novel imaging methods and systematic search of optimal combinations of measurement conditions and data processing parameters. This can help to increase the success rates and efficiency of valuable synchrotron beam time. To demonstrate the capabilities of the framework, various experiments have been simulated and compared with real data. To show the use case of measurement and data processing parameter optimization based on simulation, a virtual counterpart of a high-speed radiography experiment was created and the simulated data were used to select a suitable motion estimation algorithm; one of its parameters was optimized in order to achieve the best motion estimation accuracy when applied on the real data. syris was also used to simulate tomographic data sets under various imaging conditions which impact the tomographic reconstruction accuracy, and it is shown how the accuracy may guide the selection of imaging conditions for particular use cases.The flexible and efficient framework syris is presented and its capabilities for the simulation of four-dimensional X-ray imaging experiments are demonstrated by two exemplary applications.
Kopmann A., Chilingaryan S., Vogelgesang M., Dritschler T., Shkarin A., Shkarin R., Dos Santos Rolo T., Farago T., Van De Kamp T., Balzer M., Caselle M., Weber M., Baumbach T.
in 2016 IEEE Nuclear Science Symposium, Medical Imaging Conference and Room-Temperature Semiconductor Detector Workshop, NSS/MIC/RTSD 2016, 2017-January (2017), 8069895. DOI:10.1109/NSSMIC.2016.8069895
© 2016 IEEE. New imaging stations aim for high spatial and temporal resolution and are characterized by ever increasing sampling rates and demanding data processing workflows. Key to successful imaging experiments is to open up high-performance computing resources. This includes carefully selected components for computing hardware and development of advanced imaging algorithms optimized for efficient use of parallel processor architectures. We present the novel UFO computing platform for online data processing for imaging experiments and image-based feedback. The platform handles the full data life cycle from the X-ray detector to long-term data archives. Core components of this system are an FPGA platform for ultra-fast data acquisition, the GPU-based UFO image processing framework, and the fast control system “Concert”. Reconstruction algorithms implemented in the UFO framework are optimized for the latest GPU architectures and provide a reconstruction throughput in the GB/s-range. The control system “Concert” integrates high-speed computing nodes and fast beamline devices and thus enables image-based control loops and advanced workflow automation for efficient beam time usage. Low latencies are ensured by direct communication between FPGA and GPUs using AMDs DirectGMA technology. Time resolved tomography is supported by cutting edge regularization methods for high quality reconstructions with a reduced number of projections. The new infrastructure at ANKA has dramatically accelerated tomography from hours to second and resulted in new application fields, like high-throughput tomography, pump-probe radiography and stroboscopic tomography. Ultra-fast X-ray cine-tomography for the first time allows one to observe internal dynamics of moving millimeter-sized objects in real-time.
Schmelzle S., Heethoff M., Heuveline V., Losel P., Becker J., Beckmann F., Schluenzen F., Hammel J.U., Kopmann A., Mexner W., Vogelgesang M., Jerome N.T., Betz O., Beutel R., Wipfler B., Blanke A., Harzsch S., Hornig M., Baumbach T., Van De Kamp T.
in Proceedings of SPIE – The International Society for Optical Engineering, 10391 (2017), 103910P. DOI:10.1117/12.2275959
© 2017 SPIE. Beamtime and resulting SRμCT data are a valuable resource for researchers of a broad scientific community in life sciences. Most research groups, however, are only interested in a specific organ and use only a fraction of their data. The rest of the data usually remains untapped. By using a new collaborative approach, the NOVA project (Network for Online Visualization and synergistic Analysis of tomographic data) aims to demonstrate, that more efficient use of the valuable beam time is possible by coordinated research on different organ systems. The biological partners in the project cover different scientific aspects and thus serve as model community for the collaborative approach. As proof of principle, different aspects of insect head morphology will be investigated (e.g., biomechanics of the mouthparts, and neurobiology with the topology of sensory areas). This effort is accomplished by development of advanced analysis tools for the ever-increasing quantity of tomographic datasets. In the preceding project ASTOR, we already successfully demonstrated considerable progress in semi-automatic segmentation and classification of internal structures. Further improvement of these methods is essential for an efficient use of beam time and will be refined in the current NOVAproject. Significant enhancements are also planned at PETRA III beamline p05 to provide all possible contrast modalities in x-ray imaging optimized to biological samples, on the reconstruction algorithms, and the tools for subsequent analyses and management of the data. All improvements made on key technologies within this project will in the long-term be equally beneficial for all users of tomography instrumentations.
Mohr H., Dritschler T., Ardila L.E., Balzer M., Caselle M., Chilingaryan S., Kopmann A., Rota L., Schuh T., Vogelgesang M., Weber M.
in Journal of Instrumentation, 12 (2017), C04019. DOI:10.1088/1748-0221/12/04/C04019
© 2017 IOP Publishing Ltd and Sissa Medialab srl. In this work, we investigate the use of GPUs as a way of realizing a low-latency, high-throughput track trigger, using CMS as a showcase example. The CMS detector at the Large Hadron Collider (LHC) will undergo a major upgrade after the long shutdown from 2024 to 2026 when it will enter the high luminosity era. During this upgrade, the silicon tracker will have to be completely replaced. In the High Luminosity operation mode, luminosities of 5-7 × 1034 cm-2s-1 and pileups averaging at 140 events, with a maximum of up to 200 events, will be reached. These changes will require a major update of the triggering system. The demonstrated systems rely on dedicated hardware such as associative memory ASICs and FPGAs. We investigate the use of GPUs as an alternative way of realizing the requirements of the L1 track trigger. To this end we implemeted a Hough transformation track finding step on GPUs and established a low-latency RDMA connection using the PCIe bus. To showcase the benefits of floating point operations, made possible by the use of GPUs, we present a modified algorithm. It uses hexagonal bins for the parameter space and leads to a more truthful representation of the possible track parameters of the individual hits in Hough space. This leads to fewer duplicate candidates and reduces fake track candidates compared to the regular approach. With data-transfer latencies of 2 μs and processing times for the Hough transformation as low as 3.6 μs, we can show that latencies are not as critical as expected. However, computing throughput proves to be challenging due to hardware limitations.