X-ray phase contrast imaging holds great promise for improving the visibility of light-element materials such as soft tissues and tumors. The single-mask differential phase contrast imaging method stands out as a simple and effective approach to yield differential phase contrast. In this work, we introduce a model for a single-mask phase imaging system based on the transport-of-intensity equation. Our model provides an accessible understanding of signal and contrast formation in single-mask x-ray phase imaging, offering a clear perspective on the image formation process, for example, the origin of alternate bright and dark fringes in phase contrast intensity images. Aided by our model, we present an efficient retrieval method that yields differential phase contrast imagery in a single acquisition step. Our model gives insight into the contrast generation and its dependence on the system geometry and imaging parameters in both the initial intensity image as well as retrieved images. The model validity as well as the proposed retrieval method are demonstrated via both experimental results on a system developed in house as well as Monte Carlo simulations. In conclusion, our work not only provides a model for an intuitive visualization of image formation but also offers a method to optimize differential phase imaging setups, holding tremendous promise for advancing medical diagnostics and other applications.
X-ray phase contrast imaging (PCI) combined with phase retrieval has the potential to improve soft-material visibility and discrimination. This work examined the accuracy, image quality gains, and robustness of a spectral phase retrieval method proposed by our group. Spectroscopic PCI measurements of a physical phantom were obtained using state-of-the-art photon-counting detectors in combination with a polychromatic x-ray source. The phantom consisted of four poorly attenuating materials. Excellent accuracy was demonstrated in simultaneously retrieving the complete refractive properties (photoelectric absorption, attenuation, and phase) of these materials. Approximately 10 times higher SNR was achieved in retrieved images compared to the original PCI intensity image. These gains are also shown to be robust against increasing quantum noise, even for acquisition times as low as 1 s with a low-flux microfocus x-ray tube (average counts of 250 photons/pixels). We expect that this spectral phase retrieval method, adaptable to several PCI geometries, will allow significant dose reduction and improved material discrimination in clinical and industrial x-ray imaging applications.
more » « less- Award ID(s):
- 1652892
- NSF-PAR ID:
- 10205283
- Publisher / Repository:
- Optical Society of America
- Date Published:
- Journal Name:
- Journal of the Optical Society of America A
- Volume:
- 38
- Issue:
- 1
- ISSN:
- 1084-7529; JOAOD6
- Page Range / eLocation ID:
- Article No. 71
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
more » « less
Abstract Purpose The constrained one‐step spectral CT image reconstruction (cOSSCIR) algorithm with a nonconvex alternating direction method of multipliers optimizer is proposed for addressing computed tomography (CT) metal artifacts caused by beam hardening, noise, and photon starvation. The quantitative performance of cOSSCIR is investigated through a series of photon‐counting CT simulations.
Methods cOSSCIR directly estimates basis material maps from photon‐counting data using a physics‐based forward model that accounts for beam hardening. The cOSSCIR optimization framework places constraints on the basis maps, which we hypothesize will stabilize the decomposition and reduce streaks caused by noise and photon starvation. Another advantage of cOSSCIR is that the spectral data need not be registered, so that a ray can be used even if some energy window measurements are unavailable. Photon‐counting CT acquisitions of a virtual pelvic phantom with low‐contrast soft tissue texture and bilateral hip prostheses were simulated. Bone and water basis maps were estimated using the cOSSCIR algorithm and combined to form a virtual monoenergetic image for the evaluation of metal artifacts. The cOSSCIR images were compared to a “two‐step” decomposition approach that first estimated basis sinograms using a maximum likelihood algorithm and then reconstructed basis maps using an iterative total variation constrained least‐squares optimization (MLE+TV). Images were also compared to a nonspectral TV reconstruction of the total number of counts detected for each ray with and without normalized metal artifact reduction (NMAR) applied. The simulated metal density was increased to investigate the effects of increasing photon starvation. The quantitative error and standard deviation in regions of the phantom were compared across the investigated algorithms. The ability of cOSSCIR to reproduce the soft‐tissue texture, while reducing metal artifacts, was quantitatively evaluated.
Results Noiseless simulations demonstrated the convergence of the cOSSCIR and MLE+TV algorithms to the correct basis maps in the presence of beam‐hardening effects. When noise was simulated, cOSSCIR demonstrated a quantitative error of −1 HU, compared to 2 HU error for the MLE+TV algorithm and −154 HU error for the nonspectral TV+NMAR algorithm. For the cOSSCIR algorithm, the standard deviation in the central iodine region of interest was 20 HU, compared to 299 HU for the MLE+TV algorithm, 41 HU for the MLE+TV+Mask algorithm that excluded rays through metal, and 55 HU for the nonspectral TV+NMAR algorithm. Increasing levels of photon starvation did not impact the bias or standard deviation of the cOSSCIR images. cOSSCIR was able to reproduce the soft‐tissue texture when an appropriate regularization constraint value was selected.
Conclusions By directly inverting photon‐counting CT data into basis maps using an accurate physics‐based forward model and a constrained optimization algorithm, cOSSCIR avoids metal artifacts due to beam hardening, noise, and photon starvation. The cOSSCIR algorithm demonstrated improved stability and accuracy compared to a two‐step method of decomposition followed by reconstruction.
-
Inertial confinement fusion (ICF) holds increasing promise as a potential source of abundant, clean energy, but has been impeded by defects such as micro-voids in the ablator layer of the fuel capsules. It is critical to understand how these micro-voids interact with the laser-driven shock waves that compress the fuel pellet. At the Matter in Extreme Conditions (MEC) instrument at the Linac Coherent Light Source (LCLS), we utilized an x-ray pulse train with ns separation, an x-ray microscope, and an ultrafast x-ray imaging (UXI) detector to image shock wave interactions with micro-voids. To minimize the high- and low-frequency variations of the captured images, we incorporated principal component analysis (PCA) and image alignment for flat-field correction. After applying these techniques we generated phase and attenuation maps from a 2D hydrodynamic radiation code (xRAGE), which were used to simulate XPCI images that we qualitatively compare with experimental images, providing a one-to-one comparison for benchmarking material performance. Moreover, we implement a transport-of-intensity (TIE) based method to obtain the average projected mass density (areal density) of our experimental images, yielding insight into how defect-bearing ablator materials alter microstructural feature evolution, material compression, and shock wave propagation on ICF-relevant time scales.more » « less
-
Telecystoscopy can lower the barrier to access critical urologic diagnostics for patients around the world. A major challenge for robotic control of flexible cystoscopes and intuitive teleoperation is the pose estimation of the scope tip. We propose a novel real-time camera localization method using video recordings from a prior cystoscopy and 3D bladder reconstruction to estimate cystoscope pose within the bladder during follow-up telecystoscopy. We map prior video frames into a low-dimensional space as a dictionary so that a new image can be likewise mapped to efficiently retrieve its nearest neighbor among the dictionary images. The cystoscope pose is then estimated by the correspondence among the new image, its nearest dictionary image, and the prior model from 3D reconstruction. We demonstrate performance of our methods using bladder phantoms with varying fidelity and a servo-controlled cystoscope to simulate the use case of bladder surveillance through telecystoscopy. The servo-controlled cystoscope with 3 degrees of freedom (angulation, roll, and insertion axes) was developed for collecting cystoscope videos from bladder phantoms. Cystoscope videos were acquired in a 2.5D bladder phantom (bladder-shape cross-section plus height) with a panorama of a urothelium attached to the inner surface. Scans of the 2.5D phantom were performed in separate arc trajectories each of which is generated by actuation on the angulation with a fixed roll and insertion length. We further included variance in moving speed, imaging distance and existence of bladder tumors. Cystoscope videos were also acquired in a water-filled 3D silicone bladder phantom with hand-painted vasculature. Scans of the 3D phantom were performed in separate circle trajectories each of which is generated by actuation on the roll axis under a fixed angulation and insertion length. These videos were used to create 3D reconstructions, dictionary sets, and test data sets for evaluating the computational efficiency and accuracy of our proposed method in comparison with a method based on global Scale-Invariant Feature Transform (SIFT) features, named SIFT-only. Our method can retrieve the nearest dictionary image for 94–100% of test frames in under 55[Formula: see text]ms per image, whereas the SIFT-only method can only find the image match for 56–100% of test frames in 6000–40000[Formula: see text]ms per image depending on size of the dictionary set and richness of SIFT features in the images. Our method, with a speed of around 20 Hz for the retrieval stage, is a promising tool for real-time image-based scope localization in robotic cystoscopy when prior cystoscopy images are available.more » « less
-
Abstract Acquiring detailed 3D images of samples is needed for conducting thorough investigations in a wide range of applications. Doing so using nondestructive methods such as X-ray computed tomography (X-ray CT) has resolution limitations. Destructive methods, which work based on consecutive delayering and imaging of the sample, face a tradeoff between throughput and resolution. Using focused ion beam (FIB) for delayering, although high precision, is low throughput. On the other hand, mechanical methods that can offer fast delayering, are low precision and may put the sample integrity at risk. Herein, we propose to use femtosecond laser ablation as a delayering method in combination with optical and confocal microscopy as the imaging technique for performing rapid 3D imaging. The use of confocal microscopy provides several advantages. First, it eliminates the 3D image distortion resulting from non-flat layers, caused by the difference in laser ablation rate of different materials. It further allows layer height variations to be maintained within a small range. Finally, it enables material characterization based on the processing of material ablation rate at different locations. The proposed method is applied on a printed circuit board (PCB), and the results are validated and compared with the X-ray CT image of the PCB part.more » « less
