Towards retrospective motion correction and reconstruction for clinical 3D brain MRI protocols with a reference contrast.

OBJECT: In a typical MR session, several contrasts are acquired. Due to the sequential nature of the data acquisition process, the patient may experience some discomfort at some point during the session, and start moving. Hence, it is quite common to have MR sessions where some contrasts are well-resolved, while other contrasts exhibit motion artifacts. Instead of repeating the scans that are corrupted by motion, we introduce a reference-guided retrospective motion correction scheme that takes advantage of the motion-free scans, based on a generalized rigid registration routine. MATERIALS AND METHODS: We focus on various existing clinical 3D brain protocols at 1.5 Tesla MRI based on Cartesian sampling. Controlled experiments with three healthy volunteers and three levels of motion are performed. RESULTS: Radiological inspection confirms that the proposed method consistently ameliorates the corrupted scans. Furthermore, for the set of specific motion tests performed in this study, the quality indexes based on PSNR and SSIM shows only a modest decrease in correction quality as a function of motion complexity. DISCUSSION: While the results on controlled experiments are positive, future applications to patient data will ultimately clarify whether the proposed correction scheme satisfies the radiological requirements.

Quantifying the effect of X-ray scattering for data generation in real-time defect detection.

BACKGROUND: X-ray imaging is widely used for the non-destructive detection of defects in industrial products on a conveyor belt. In-line detection requires highly accurate, robust, and fast algorithms. Deep Convolutional Neural Networks (DCNNs) satisfy these requirements when a large amount of labeled data is available. To overcome the challenge of collecting these data, different methods of X-ray image generation are considered. OBJECTIVE: Depending on the desired degree of similarity to real data, different physical effects should either be simulated or can be ignored. X-ray scattering is known to be computationally expensive to simulate, and this effect can greatly affect the accuracy of a generated X-ray image. We aim to quantitatively evaluate the effect of scattering on defect detection. METHODS: Monte-Carlo simulation is used to generate X-ray scattering distribution. DCNNs are trained on the data with and without scattering and applied to the same test datasets. Probability of Detection (POD) curves are computed to compare their performance, characterized by the size of the smallest detectable defect. RESULTS: We apply the methodology to a model problem of defect detection in cylinders. When trained on data without scattering, DCNNs reliably detect defects larger than 1.3 mm, and using data with scattering improves performance by less than 5%. If the analysis is performed on the cases with large scattering-to-primary ratio (1 < SPR < 5), the difference in performance could reach 15% (approx. 0.4 mm). CONCLUSION: Excluding the scattering signal from the training data has the largest effect on the smallest detectable defects, and the difference decreases for larger defects. The scattering-to-primary ratio has a significant effect on detection performance and the required accuracy of data generation.

Maximum-likelihood estimation in ptychography in the presence of Poisson-Gaussian noise statistics.

Optical measurements often exhibit mixed Poisson-Gaussian noise statistics, which hampers the image quality, particularly under low signal-to-noise ratio (SNR) conditions. Computational imaging falls short in such situations when solely Poissonian noise statistics are assumed. In response to this challenge, we define a loss function that explicitly incorporates this mixed noise nature. By using a maximum-likelihood estimation, we devise a practical method to account for a camera readout noise in gradient-based ptychography optimization. Our results, based on both experimental and numerical data, demonstrate that this approach outperforms the conventional one, enabling enhanced image reconstruction quality under challenging noise conditions through a straightforward methodological adjustment.

Maximum-likelihood estimation in ptychography in the presence of Poisson-Gaussian noise statistics: publisher's note.

This publisher's note contains a correction to Opt. Lett.48, 6027 (2023)10.1364/OL.502344.

Tomosipo: fast, flexible, and convenient 3D tomography for complex scanning geometries in Python.

Tomography is a powerful tool for reconstructing the interior of an object from a series of projection images. Typically, the source and detector traverse a standard path (e.g., circular, helical). Recently, various techniques have emerged that use more complex acquisition geometries. Current software packages require significant handwork, or lack the flexibility to handle such geometries. Therefore, software is needed that can concisely represent, visualize, and compute reconstructions of complex acquisition geometries. We present tomosipo, a Python package that provides these capabilities in a concise and intuitive way. Case studies demonstrate the power and flexibility of tomosipo.

ForametCeTera, a novel CT scan dataset to expedite classification research of (non-)foraminifera.

This paper introduces ForametCeTera, a pioneering dataset designed to address the challenges associated with automating the analysis of benthic foraminifera in sediment cores. Foraminifera are sensitive sentinels of environmental change and are a crucial component of carbonate-denominated ecosystems, such as coral reefs. Studying their prevalence and characteristics is imperative in understanding climate change. However, analysis of foraminifera contained in core samples currently requires washing, sieving and manual quantification. These methods are thus time-consuming and require trained experts. To overcome these limitations, we propose an alternative workflow utilizing 3D X-ray computational tomography (CT) for fully automated analysis, saving time and resources. Despite recent advancements in automation, a crucial lack of methods persists for segmenting and classifying individual foraminifera from 3D scans. In response, we present ForametCeTera, a diverse dataset featuring 436 3D CT scans of individual foraminifera and non-foraminiferan material following a high-throughput scanning workflow. ForametCeTera serves as a foundational resource for generating synthetic digital core samples, facilitating the development of segmentation and classification methods of entire core sample CT scans.

Task-Adaptive Angle Selection for Computed Tomography-Based Defect Detection.

Sparse-angle X-ray Computed Tomography (CT) plays a vital role in industrial quality control but leads to an inherent trade-off between scan time and reconstruction quality. Adaptive angle selection strategies try to improve upon this based on the idea that the geometry of the object under investigation leads to an uneven distribution of the information content over the projection angles. Deep Reinforcement Learning (DRL) has emerged as an effective approach for adaptive angle selection in X-ray CT. While previous studies focused on optimizing generic image quality measures using a fixed number of angles, our work extends them by considering a specific downstream task, namely image-based defect detection, and introducing flexibility in the number of angles used. By leveraging prior knowledge about typical defect characteristics, our task-adaptive angle selection method, adaptable in terms of angle count, enables easy detection of defects in the reconstructed images.

Enabling 3D CT-scanning of cultural heritage objects using only in-house 2D X-ray equipment in museums.

Visualizing the internal structure of museum objects is a crucial step in acquiring knowledge about the origin, state, and composition of cultural heritage artifacts. Among the most powerful techniques for exposing the interior of museum objects is computed tomography (CT), a technique that computationally forms a 3D image using hundreds of radiographs acquired in a full circular range. However, the lack of affordable and versatile CT equipment in museums, combined with the challenge of transporting precious collection objects, currently keeps this technique out of reach for most cultural heritage applications. We propose an approach for creating accurate CT reconstructions using only standard 2D radiography equipment already available in most larger museums. Specifically, we demonstrate that a combination of basic X-ray imaging equipment, a tailored marker-based image acquisition protocol, and sophisticated data-processing algorithms, can achieve 3D imaging of collection objects without the need for a costly CT imaging system. We implemented this approach in the British Museum (London), the J. Paul Getty Museum (Los Angeles), and the Rijksmuseum (Amsterdam). Our work paves the way for broad facilitation and adoption of CT technology across museums worldwide.

CT-based data generation for foreign object detection on a single X-ray projection.

Although X-ray imaging is used routinely in industry for high-throughput product quality control, its capability to detect internal defects has strong limitations. The main challenge stems from the superposition of multiple object features within a single X-ray view. Deep Convolutional neural networks can be trained by annotated datasets of X-ray images to detect foreign objects in real-time. However, this approach depends heavily on the availability of a large amount of data, strongly hampering the viability of industrial use with high variability between batches of products. We present a computationally efficient, CT-based approach for creating artificial single-view X-ray data based on just a few physically CT-scanned objects. By algorithmically modifying the CT-volume, a large variety of training examples is obtained. Our results show that applying the generative model to a single CT-scanned object results in image analysis accuracy that would otherwise be achieved with scans of tens of real-world samples. Our methodology leads to a strong reduction in training data needed, improved coverage of the combinations of base and foreign objects, and extensive generalizability to additional features. Once trained on just a single CT-scanned object, the resulting deep neural network can detect foreign objects in real-time with high accuracy.

2DeteCT - A large 2D expandable, trainable, experimental Computed Tomography dataset for machine learning.

Recent research in computational imaging largely focuses on developing machine learning (ML) techniques for image reconstruction, which requires large-scale training datasets consisting of measurement data and ground-truth images. However, suitable experimental datasets for X-ray Computed Tomography (CT) are scarce, and methods are often developed and evaluated only on simulated data. We fill this gap by providing the community with a versatile, open 2D fan-beam CT dataset suitable for developing ML techniques for a range of image reconstruction tasks. To acquire it, we designed a sophisticated, semi-automatic scan procedure that utilizes a highly-flexible laboratory X-ray CT setup. A diverse mix of samples with high natural variability in shape and density was scanned slice-by-slice (5,000 slices in total) with high angular and spatial resolution and three different beam characteristics: A high-fidelity, a low-dose and a beam-hardening-inflicted mode. In addition, 750 out-of-distribution slices were scanned with sample and beam variations to accommodate robustness and segmentation tasks. We provide raw projection data, reference reconstructions and segmentations based on an open-source data processing pipeline.

Unsupervised Foreign Object Detection Based on Dual-Energy Absorptiometry in the Food Industry.

X-ray imaging is a widely used technique for non-destructive inspection of agricultural food products. One application of X-ray imaging is the autonomous, in-line detection of foreign objects in food samples. Examples of such inclusions are bone fragments in meat products, plastic and metal debris in fish, and fruit infestations. This article presents a processing methodology for unsupervised foreign object detection based on dual-energy X-ray absorptiometry (DEXA). A novel thickness correction model is introduced as a pre-processing technique for DEXA data. The aim of the model is to homogenize regions in the image that belong to the food product and to enhance contrast where the foreign object is present. In this way, the segmentation of the foreign object is more robust to noise and lack of contrast. The proposed methodology was applied to a dataset of 488 samples of meat products acquired from a conveyor belt. Approximately 60% of the samples contain foreign objects of different types and sizes, while the rest of the samples are void of foreign objects. The results show that samples without foreign objects are correctly identified in 97% of cases and that the overall accuracy of foreign object detection reaches 95%.

Micromagnetic Tomography for Paleomagnetism and Rock-Magnetism.

Our understanding of the past behavior of the geomagnetic field arises from magnetic signals stored in geological materials, e.g., (volcanic) rocks. Bulk rock samples, however, often contain magnetic grains that differ in chemistry, size, and shape; some of them record the Earth's magnetic field well, others are unreliable. The presence of a small amount of adverse behaved magnetic grains in a sample may already obscure important information on the past state of the geomagnetic field. Recently it was shown that it is possible to determine magnetizations of individual grains in a sample by combining X-ray computed tomography and magnetic surface scanning measurements. Here we establish this new Micromagnetic Tomography (MMT) technique and make it suitable for use with different magnetic scanning techniques, and for both synthetic and natural samples. We acquired reliable magnetic directions by selecting subsets of grains in a synthetic sample, and we obtained rock-magnetic information of individual grains in a volcanic sample. This illustrates that MMT opens up entirely new venues of paleomagnetic and rock-magnetic research. MMT's unique ability to determine the magnetization of individual grains in a nondestructive way allows for a systematic analysis of how geological materials record and retain information on the past state of the Earth's magnetic field. Moreover, by interpreting only the contributions of known magnetically well-behaved grains in a sample, MMT has the potential to unlock paleomagnetic information from even the most complex, crucial, or valuable recorders that current methods are unable to recover.

Task-Driven Learned Hyperspectral Data Reduction Using End-to-End Supervised Deep Learning.

An important challenge in hyperspectral imaging tasks is to cope with the large number of spectral bins. Common spectral data reduction methods do not take prior knowledge about the task into account. Consequently, sparsely occurring features that may be essential for the imaging task may not be preserved in the data reduction step. Convolutional neural network (CNN) approaches are capable of learning the specific features relevant to the particular imaging task, but applying them directly to the spectral input data is constrained by the computational efficiency. We propose a novel supervised deep learning approach for combining data reduction and image analysis in an end-to-end architecture. In our approach, the neural network component that performs the reduction is trained such that image features most relevant for the task are preserved in the reduction step. Results for two convolutional neural network architectures and two types of generated datasets show that the proposed Data Reduction CNN (DRCNN) approach can produce more accurate results than existing popular data reduction methods, and can be used in a wide range of problem settings. The integration of knowledge about the task allows for more image compression and higher accuracies compared to standard data reduction methods.

Point spread function based image reconstruction in optical projection tomography.

As a result of the shallow depth of focus of the optical imaging system, the use of standard filtered back projection in optical projection tomography causes space-variant tangential blurring that increases with the distance to the rotation axis. We present a novel optical tomographic image reconstruction technique that incorporates the point spread function of the imaging lens in an iterative reconstruction. The technique is demonstrated using numerical simulations, tested on experimental optical projection tomography data of single fluorescent beads, and applied to high-resolution emission optical projection tomography imaging of an entire zebrafish larva. Compared to filtered back projection our results show greatly reduced radial and tangential blurring over the entire [Formula: see text] mm2 field of view, and a significantly improved signal to noise ratio.