An attenuation-partition based iterative phase retrieval algorithm for in-line phase-contrast imaging

Phase retrieval-based phase-contrast CT for vascular imaging with microbubble contrast agent

PURPOSE

The introduction of microbubble contrast agent into tissues can create significant phase shifts. Phase retrieval (PR)-based phase-contrast computed tomography (PCCT) is an imaging method for retrieving and reconstructing the phase shifts within an object. This study aimed to evaluate the feasibility of PR-based PCCT with microbubble contrast agent for vascular imaging.

METHODS

Projection phase-contrast images of individual microbubbles and a cluster of microbubbles were captured and compared. Contrast enhancement from microbubbles was evaluated by comparing to the gold standard iodine-based contrast agent in vitro. The arterial systems of 14 Sprague-Dawley rats were perfused with microbubbles or saline. The rat hearts and the arterial systems were excised and imaged ex vivo. CT imaging was performed at the energy of 22 keV. PR was performed using the phase-attenuation duality (PAD) method with different δ/β values (PAD property). The contrast-to-noise ratio (CNR) was used for quantitatively assessing the contrast enhancement.

RESULTS

Individual microbubbles functioned as a lens to focus the x rays, whereas, a cluster of microbubbles scattered the x rays. In the in vitro experiment, the contrast enhancement from iodine was significantly greater than that from microbubbles (P < 0.05). In the heart samples, the CNRs for microbubbles on PR-based PCCT were significantly greater than those on absorption-contrast CT (ACCT) and PR-free PCCT (both P < 0.001). The CNRs for microbubbles were also significantly greater than those for saline on PR-based PCCT in the samples (P < 0.001). Although they provided weaker contrast enhancement than that from iodine, microbubbles could still provide sufficient contrast enhancement to clearly show the 3D architecture of rat aortas and their main branches.

CONCLUSION

The imaging modality can currently be used as a complement or alternative to absorption-based microCT for imaging vessels in biological samples.

Dose and detectability improvements with high energy phase sensitive x-ray imaging in comparison to low energy conventional imaging

The objective of this study was to demonstrate the potential benefits of using high energy x-rays for phase sensitive breast imaging through a comparison with conventional mammography imaging. We compared images of a contrast-detail phantom acquired on a prototype phase sensitive x-ray imaging system with images acquired on a commercial flat panel digital mammography unit. The phase contrast images were acquired using a micro-focus x-ray source with a 50 µm focal spot at 120 kVp and 4.5 mAs, with a magnification factor of 2.46 and a 50 µm pixel pitch. A phase attenuation duality-based phase retrieval algorithm that requires only a single phase contrast image was applied. Conventional digital mammography images were acquired at 27 kVp, 131 mAs and 28 kVp, 54 mAs. For the same radiation dose, both the observer study and signal-to-noise ratio (SNR)/figure of merit comparisons indicated a large improvement by the phase retrieved image as compared to the clinical system for the larger disc sizes, but the improvement was not enough to detect the smallest discs. Compared to the double dose image acquired with the clinical system, the observer study also indicated that the phase retrieved image provided improved detection capabilities for all disc sizes except the smallest discs. Thus the SNR improvement provided by phase contrast imaging is not yet enough to offset the noise reduction provided by the clinical system at the doubled dose level. However, the potential demonstrated by this study for high energy phase sensitive x-ray imaging to improve lesion detection and reduce radiation dose in mammography warrants further investigation of this technique.

In-line phase shift tomosynthesis

PURPOSE

The purpose of this work is to (1) demonstrate laboratory measurements of phase shift images derived from in-line phase-contrast radiographs using the attenuation-partition based algorithm (APBA) of Yan et al. [Opt. Express 18(15), 16074-16089 (2010)], (2) verify that the APBA reconstructed images obey the linearity principle, and (3) reconstruct tomosynthesis phase shift images from a collection of angularly sampled planar phase shift images.

METHODS

An unmodified, commercially available cabinet x-ray system (Faxitron LX-60) was used in this experiment. This system contains a tungsten anode x-ray tube with a nominal focal spot size of 10 μm. The digital detector uses CsI∕CMOS with a pixel size of 50×50 μm. The phantoms used consisted of one acrylic plate, two polystyrene plates, and a habanero pepper. Tomosynthesis images were reconstructed from 51 images acquired over a ±25° arc. All phase shift images were reconstructed using the APBA.

RESULTS

Image contrast derived from the planar phase shift image of an acrylic plate of uniform thickness exceeded the contrast of the traditional attenuation image by an approximate factor of two. Comparison of the planar phase shift images from a single, uniform thickness polystyrene plate with two polystyrene plates demonstrated an approximate linearity of the estimated phase shift with plate thickness (-1600 rad vs -2970 rad). Tomographic phase shift images of the habanero pepper exhibited acceptable spatial resolution and contrast comparable to the corresponding attenuation image.

CONCLUSIONS

This work demonstrated the feasibility of laboratory-based phase shift tomosynthesis and suggests that phase shift imaging could potentially provide a new imaging biomarker. Further investigation will be needed to determine if phase shift contrast will be able to provide new tissue contrast information or improved clinical performance.

A robust general phase retrieval method for medical applications

From medical imaging perspective the robustness of a phase retrieval method is of critical importance. In this presentation we compare the robustness of two general phase retrieval methods, namely the transport of intensity equation inversion (TIE-inversion) method and the attenuation partition based (AP-based) method. We showed that the TIE-inversion method, regardless if being assisted with the Tikhonov regularization, failed to retrieve the phase maps in two experimental studies. The failure exposes this method's weakness as being unstable against the noise. In contrast, the sample phase maps are retrieved successfully by using the AP-based method. The stark performance differences of the two methods are rooted in their different techniques dealing with the singularity problem. This comparison shows that the robust AP-based phase retrieval method will be superior to the TIE-inversion method for medical imaging applications where radiation doses are stringently limited.

Robustness of phase retrieval methods in X-ray phase contrast imaging: a comparison

PURPOSE

The robustness of the phase retrieval methods is of critical importance for limiting and reducing radiation doses involved in x-ray phase contrast imaging. This work is to compare the robustness of two phase retrieval methods by analyzing the phase maps retrieved from the experimental images of a phantom.

METHODS

Two phase retrieval methods were compared. One method is based on the transport of intensity equation (TIE) for phase contrast projections, and the TIE-based method is the most commonly used method for phase retrieval in the literature. The other is the recently developed attenuation-partition based (AP-based) phase retrieval method. The authors applied these two methods to experimental projection images of an air-bubble wrap phantom for retrieving the phase map of the bubble wrap. The retrieved phase maps obtained by using the two methods are compared.

RESULTS

In the wrap's phase map retrieved by using the TIE-based method, no bubble is recognizable, hence, this method failed completely for phase retrieval from these bubble wrap images. Even with the help of the Tikhonov regularization, the bubbles are still hardly visible and buried in the cluttered background in the retrieved phase map. The retrieved phase values with this method are grossly erroneous. In contrast, in the wrap's phase map retrieved by using the AP-based method, the bubbles are clearly recovered. The retrieved phase values with the AP-based method are reasonably close to the estimate based on the thickness-based measurement. The authors traced these stark performance differences of the two methods to their different techniques employed to deal with the singularity problem involved in the phase retrievals.

CONCLUSIONS

This comparison shows that the conventional TIE-based phase retrieval method, regardless if Tikhonov regularization is used or not, is unstable against the noise in the wrap's projection images, while the AP-based phase retrieval method is shown in these experiments to be superior to the TIE-based method for the robustness in performing the phase retrieval.

Potential for imaging engineered tissues with X-ray phase contrast

As the field of tissue engineering advances, it is crucial to develop imaging methods capable of providing detailed three-dimensional information on tissue structure. X-ray imaging techniques based on phase-contrast (PC) have great potential for a number of biomedical applications due to their ability to provide information about soft tissue structure without exogenous contrast agents. X-ray PC techniques retain the excellent spatial resolution, tissue penetration, and calcified tissue contrast of conventional X-ray techniques while providing drastically improved imaging of soft tissue and biomaterials. This suggests that X-ray PC techniques are very promising for evaluation of engineered tissues. In this review, four different implementations of X-ray PC imaging are described and applications to tissues of relevance to tissue engineering reviewed. In addition, recent applications of X-ray PC to the evaluation of biomaterial scaffolds and engineered tissues are presented and areas for further development and application of these techniques are discussed. Imaging techniques based on X-ray PC have significant potential for improving our ability to image and characterize engineered tissues, and their continued development and optimization could have significant impact on the field of tissue engineering.

Performance analysis of the attenuation-partition based iterative phase retrieval algorithm for in-line phase-contrast imaging

The phase retrieval is an important task in x-ray phase contrast imaging. The robustness of phase retrieval is especially important for potential medical imaging applications such as phase contrast mammography. Recently the authors developed an iterative phase retrieval algorithm, the attenuation-partition based algorithm, for the phase retrieval in inline phase-contrast imaging [1]. Applied to experimental images, the algorithm was proven to be fast and robust. However, a quantitative analysis of the performance of this new algorithm is desirable. In this work, we systematically compared the performance of this algorithm with other two widely used phase retrieval algorithms, namely the Gerchberg-Saxton (GS) algorithm and the Transport of Intensity Equation (TIE) algorithm. The systematical comparison is conducted by analyzing phase retrieval performances with a digital breast specimen model. We show that the proposed algorithm converges faster than the GS algorithm in the Fresnel diffraction regime, and is more robust against image noise than the TIE algorithm. These results suggest the significance of the proposed algorithm for future medical applications with the x-ray phase contrast imaging technique.

Phase and amplitude retrieval of objects embedded in a sinusoidal background from its diffraction pattern

Efforts of phase and amplitude retrieval from diffraction patterns have almost exclusively been applied for nonperiodic objects. We investigated the quality of retrieval of nonperiodic objects embedded in a sinusoidal background, using the approach of iterative hybrid input-output with oversampling. Two strategies were employed; one by filtering in the frequency domain prior to phase retrieval, and the other by filtering the phase or amplitude image after retrieval. Results obtained indicate better outcomes with the latter approach provided detector noise is not excessive.

Phase retrieval from one single phase contrast x-ray image

Phase retrieval is required for achieving artifact-free x-ray phase-sensitive 3D imaging. A phase-retrieval approach based on the phase-attenuation duality with high energy x-rays can greatly facilitate for phase sensitive imaging by allowing phase retrieval from only one single projection image. The previously derived phase retrieval formula is valid only for small Fresnel propagator phases corresponding to common clinical imaging tasks. In this work we presented a new duality-based phase retrieval formula that can be applied for cases with large Fresnel propagator phases corresponding to high spatial resolution imaging. The computer simulation demonstrated superiority of this new formula over the previous phase retrieval formula in reconstructing the high frequency components of imaged objects. A modified Tikhonov regularization technique has been devised for phase retrieval in cases of very high resolution and large object-detector distance such that some Fresnel propagator phases may be close or greater than pi. This new phase retrieval formula lays the foundation for implementing high-resolution phase-sensitive 3D imaging of soft tissue objects.

Closed-form algorithms for phase retrieval with an additive point signal

The phase retrieval problem is a problem of a reconstruction of a signal from the magnitude of its Fourier transform. In this paper, we consider the problem of reconstructing an unknown one-dimensional signal from the magnitude of its Fourier transform and the magnitude of the Fourier transform of another signal that is made by the addition of a point signal at the center of the original signal. After we show that there are exactly two solution signals that satisfy the given condition, we develop a closed-form algorithm that may reconstruct the two solution signals.