PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2017; 12(9): e0183396.
Published online 2017 September 8. doi:  10.1371/journal.pone.0183396
PMCID: PMC5590738

Robust phase-retrieval-based X-ray tomography for morphological assessment of early hepatic echinococcosis infection in rats

Huiqiang Liu, Conceptualization, Data curation, Formal analysis, Funding acquisition, Methodology, Writing – original draft, Writing – review & editing,1,* Chuanshan Zhang, Resources, Validation,2 Xiaoxi Fan, Visualization,1 Yingni Duan, Visualization,1 Tiqiao Xiao, Supervision,3 Guohao Du, Software, Visualization,3 Yanan Fu, Resources, Validation,3 Haigang Liu, Data curation, Validation,3,* and Hao Wen, Project administration2,*
Thomas Abraham, Editor

Abstract

Propagation-based phase-contrast computed micro-tomography (PPCT) dominates the non-destructive, three-dimensional inner-structure measurement in synchrotron-based biomedical research due to its simple experimental setup. To quantitatively visualize tiny density variations in soft tissues and organs closely related to early pathological morphology, an experimental study of synchrotron-based X-ray PPCT combined with generalized phase and attenuation duality (PAD) phase retrieval was implemented with the hepatic echinococcosis (HE) infection rat model at different stages. We quantitatively analyzed and evaluated the different pathological characterizations of hepatic echinococcosis during the development of this disease via our PAD-based PPCT and especially provided evidence that hepatic alveolar echinococcosis invades the liver tissue and spreads through blood flow systems with abundant blood supply in the early stage. Additionally, the infiltration of tiny vesicles in HE lesions can be clearly observed by our PAD-PPCT technique due to the striking contrast-to-noise ratio (CNR) and mass density resolution, which cannot be found by the medical imaging techniques, such as magnetic resonance imaging (MRI), computed tomography (CT), and ultrasound, in hospitals. The results demonstrated that our PAD-PPCT technique has a great potential for indicating the subtle structural information of pathological changes in soft biomedical specimens, especially helpful for the research of early micro-morphology of diseases.

Introduction

Hepatic echinococcosis (HE) arises from the intrahepatic growth of Echinococcus multilocularis larvae and generally spreads to other tissue and organs by the infiltration of small vesicles. Similar to a slowly developing malignant tumor, it has a high mortality rate and a high incidence in countries of the northern hemisphere that feature the industry of animal husbandry. The micro-morphological characterization of HE is essential to assessments of disease development, metastases and therapeutic effects at different pathological stages, especially for the prolonged drug therapy of albendazole (ABDZ) inpreventing HE parasite growth [13]. A variety of imaging techniques have been employed in the investigation of HE lesion severities in recent decade, including morphological and functional imaging methods. However, an effective early detection of HE is usually hard to achieve based on conventional imaging techniques due to atypical characteristics of HE, and a definitive diagnosis may further rely on serological assessment [46]. Therefore, it is necessary to present a new imaging method to assess the morphological HE lesions at an early stage. Although the typical imaging characteristics of HE by ultrasound (US), computed tomography (CT) and magnetic resonance imaging (MRI) have been reported with their own clinical features and advantages [712], further three-dimensional, high-contrast, and micron-scale morphological research on ex vivo HE pathology is essential and required to assess the HE lesion micro-changes and developing mechanism. With respect to optical or electron microscopy, the resolution can achieve a quasi-nanoscale observation on the basis of micron-level slices and staining that cannot meet the needs of nondestructive three-dimensional (3D) measurement. In our work, we presented the employment of novel synchrotron-based X-ray computed tomography to study the high-resolution 3D characteristic structure and pathological morphology of HE lesions at an early and developed stage.

In many interactions with X-rays and matter, an attenuation effect has been widely and successfully used for medical imaging applications for several decades. Conventional attenuation-based X-ray imaging, due to the penetration power of X-rays, is the central imaging technique for observing internal structures of opaque objects nondestructively [1316]. However, it is difficult to reveal weakly absorbing structures of an object consisting of low-Z materials, such as soft tissues and organs or light materials, because these structures are almost transparent to hard X-rays. Thus, attenuation-based imaging typically produces poor image contrast for low-Z samples. When X-rays traverse an object, the X-ray undergoes a considerably large phase shift as well, described by the X-ray complex refractive index n(x, y, z) = 1 − δ(x, y, z) + (x, y, z), where the decrement of the refractive index δ(x, y, z) is approximately 1000 times higher than the absorption coefficient β(x, y, z) for low-Z elements using hard X-ray imaging. Thus, acquisition of the phase contrast, rather than attenuation contrast, is anticipated to provide a much higher image contrast and a lower radiation dose for low-Z biomedical samples. Therefore, X-ray phase sensitive imaging techniques enable us to nondestructively observe inner weak-absorption structures with high phase sensitivity, which is particularly suitable for visualizing and quantifying important information of structures with fine mass density differences in early pathological micro-morphology research or material sciences [1723]. Currently, hard X-ray propagation-based phase-contrast computed micro-tomography (PPCT) techniques, particularly based on synchrotron radiation facilities, are widely and successfully applied to acquire non-destructive three-dimensional datasets of soft biomedical specimens due to the simpler experimental setup compared to other hard X-ray imaging techniques. The PPCT technique is the most direct method to achieve phase contrast by placing a detector at a distance downstream of an object without any additional X-ray optics required, which renders phase variations visible through translating the phase shift into the accessibly recorded intensity variations. The reconstructed intensity distribution contains the 3D map of the linear attenuation coefficients and the 3D map of Laplacian refraction index decrements of a sample, scaled by the object-detector distance [2428]. To obtain the quantitative extraction of pure phase maps, phase retrieval algorithms were proposed on the basis of the general X-ray intensity propagation equation or the simplified transport of intensity equation (TIE) [2936]. There are several studies on different propagation-based phase-contrast micro-tomography methods for biological samples. Their adopted phase retrieval methods of the propagation-based techniques used in these studies were limited in their performance. For example, in these studies, the single-distance phase retrieval method based on the homogeneous material assumption was used. According to this simplifying assumption, the ratio δ(x, y, z) / β(x, y, z) of the decrement of the refractive index over the absorption coefficient should be assumed to be constant across the specimen volume. To perform phase retrieval with this method, one has to assume a prioria value of δ(x, y, z) / β(x, y, z), which was usually set to the δ(x, y, z) / β(x, y, z) value of water at the given X-ray energy. However, for the biological soft tissue samples imaged in these studies, the tissues are inhomogeneous; hence, the actual δ(x, y, z) / β(x, y, z) is not a constant, and it varied by approximately 300% as reported in a study. The breaking down of the constant δ(x, y, z) / β(x, y, z) assumption in these comparative studies adversely affects the validity of the contrast and density resolution analyses of these studies. Therefore, a single-distance phase retrieval method valid for inhomogeneous low-Z material such as soft tissues is desirable. For better use of the robust phase retrieval tool for pathological micro-morphological research of hepatic echinococcosis, a robust single-distance phase retrieval method applicable to inhomogeneous soft tissue samples was presented and demonstrated to be applicable to our HE disease research. In brief, we combined our different propagation-based phase retrieval method, namely, the so-called phase-attenuation duality (PAD) method [3743], with quantitative measurement of biomedical specimens by using hard X-rays from synchrotron radiation. This PAD method of phase retrieval is based on the fact that for low-Z material such as soft tissues, regardless of homogeneous or inhomogeneous compositions, the electron density completely determines both the decrement of the refractive index and the absorption coefficient for hard X-rays, which will be explained in detail in Section II.

Materials and methods

HE specimens

The biomedical specimens used for this study were two groups of HE rat models, totaling 40 SD male rats (200 ± 20g) equally divided into a two-week feeding group and a two-month feeding group after HE inoculation, provided by the breeding center of the First Affiliated Hospital of Xinjiang Medical University (XJMU), China. They were euthanized immediately by cervical dislocation, and the ex vivo hepatic echinococcosis HE specimens without staining were kept in formalin solution (10% formalin neutral buffer solution) at room temperature. The operative procedures were carried out by strictly conforming to the guide for the care and use of laboratory animals, and the experimental protocol was approved by the animal experiment ethics committee of the First Affiliated Hospital of XJMU, China. Both HE groups were screened by using an MRI scanner (GE USA), while only the two-month group can be diagnosed as the HE rat models definitively, as shown in Fig 1. The two-week group cannot be confirmed due to the limited resolution of MRI for the early HE lesions. A process of serial graded dehydration with ethanol solution (short: graded ethanol concentration fixation-GECF) was employed in the 48 hours before our experiments, and then, samples were placed in an upright position in a plastic tube to facilitate the synchrotron-based PPCT experimental setup and to avoid motion artifacts associated with deformation or degeneration of the fresh tissues during the measurement with a highly brilliant synchrotron X-ray source.

Fig 1
Hepatic echinococcosis (HE) experimental specimen from the two-month disease rat model (a) and the corresponding MRI examination results: (b) Coronal section image; (c) Transverse section image (T1WI); (d) Transverse section image (T2WI).

Image acquisition system

The experiments were performed at the X-ray imaging and medical application beamline (BL13W1)of Shanghai Synchrotron Radiation Facility (SSRF), depicted in Fig 2. The synchrotron beam, produced from the 3.5GeV storage ring (with ring current I = 250mA, top-up mode), is monochromated with a double crystal monochromator (DCM) Si(111) to select an X-ray energy from 8 to 72keV. The optical stage is located atapproximately 35m downstream of the source, where the acquisition of raw experimental data was implemented with the PPCT technique, which is a free-space propagation method without any optical element between the object and detector. The detector (Hamamatsu, Japan) in our experiments consists of a scientific CMOS camera with a 6.5μm pixel size, coupled with a P43 powder (Gd2O2S: Tb+) phosphor screen with a 10μm thickness for converting the incident hard X-rays into visible light imaging. Propagation over an object-to-detector distance (ODD = Zp), as shown in Fig 2, results in a Fresnel diffraction pattern that contains the contribution from the phase distortion induced by the object. To yield highly detailed three-dimensional views of the non-destructive inner structure of HE rat models based on this PPCT system, a HE sample is rotated around the rotation axis, and a series of raw images are recorded for each projection angle. The tomographic reconstruction can then be achieved from the phase-retrieved projection dataset by using the standard filtered back-projection (FBP) algorithm. The PPCT experimental conditions were as follows: the number of raw imageswas1080 projections over 180° of rotation, acquired at Zp = 200mmwith aflat-field correction interval of 36°and an exposure time of 500ms per image at 60keV.

Fig 2
Schematic representation of the PPCTexperimental setup at SSRF.

Phase retrieval

The phase retrieval (PR) algorithm based on generalized phase-attenuation duality (PAD)proposed in our previous works [4143], suitable for meeting the requirement of limited radiation exposure in biomedical imaging, was adopted to recover phase maps induced by specimens from only a single free-space propagation radiograph. The PAD phase retrieval holds in good approximation in our biomedical specimen, as a soft liver tissue is composed mainly of low-Z elements (Z<10) and the hard X-ray energy (≥ 60keV) is high enough to make the Compton scattering dominate the soft tissue attenuation [37, 38]. We can apply a generalized PAD relationship to solve for the phase map [var phi]θ of a given view angle θ, which can be retrieved from a single projection image Iθ in this angle as follows [41, 42]:

φθ(r)=λreσKNln{[cos(λR24πM2)(2λreσKNλR24πM2)sin(λR24πM2)]1(M2Iθ(rD)IIN)}
(1)

where λ is the X-ray wavelength, re is the classic atomic radius, and σKN is the Klein-Nishina total crosssection for X-ray photon Compton scattering from a single free electron. M = (R1+R2)/R1, where R1 is the source-to-sample distance and R2 is the sample-to-detector distance (R2 = Zp); here, for synchrotron radiation source, M≈1 because R1>>R2. In addition, [cos(λR24πM2)(2λreσKNλR24πM2)sin(λR24πM2)]1 is a convolved pseudo-differential operator. Note that the large value of the ratio λre / σKN reflects the high sensitivity of phase-contrast imaging, which is more than 1000 on the basis of our experimental condition. Prior to the process of phase retrieval for each projection, background correction was performed through the normalization of dark signals and flat-field images.

Results and discussion

Calibration of PAD-PPCT with standard specimens

To quantitatively investigate HE specimens, three standard polymer test samples, such as polypropylene (PP, 0.92 g/cm3), polystyrene (PS, 1.05 g/cm3), and acrylic (PMMA, 1.19 g/cm3) tubes, were measured to demonstrate the strikingly enhanced phase contrast of our robust phase retrieval algorithm and calibrate the PAD-PPCT system. Fig 3(a) was directly reconstructed from the PPCT technique without the PAD algorithm and shows strong edge-enhanced material boundaries as a result of the refraction and diffraction of the phase-shifted X-rays. The reconstructed abnormal linear attenuation coefficients at interfaces and boundaries enable the outlines of the three standard materials to be seen clearly, but their bulk attenuation contrasts are very weak because of the weakly absorbing materials. This results in complete overlapping of the three different materials' peaks on a histogram, displayed in Fig 3(c), which is not in accordance with their actual material density distribution. However, the reconstructed section image based on our PAD-PPCT technique, shown in Fig 3(b), exhibited good conformity on both the real material structure and density distributions of the standard specimens, with totally separated material peaks in Fig 3(d). The comparison of both experimental results demonstrated that the PAD-PPCT technique has great potential for discriminating and visualizing the tiny density and structural differences of low-Z materials or soft biological tissues.

Fig 3
Comparison of reconstructed section imagesof standard material specimens in terms of contrast and histogram.

For quantitative analysis of experimental results and estimation of the HE specimen's density distributions, it is necessary to calibrate this imaging system through the measured pixel values obtained from the PAD-PPCT experiments, which were compared to their theoretical phase shift values based on the fitting line by using the least square method. The calibration result is described in Fig 4, and a linear relationship was found between the experimental and theoretical values for the quantitative measurement. The calibration factor of the PAD-PPCT system is 1.5, mainly depending on the X-ray detector.

Fig 4
Relationship between the theoretical phase-shift values and the measuredvalues based on PAD-PPCT.

Measurement of HE specimens based on PAD-PPCT

Without using a contrast agent or staining methods, the phase contrast sectional images of a two-week HE rat model are displayed in Fig 5(a), 5(b) and 5(c), respectively representing transversal, sagittal and coronal slices reconstructed only by the PAD-PPCT technique. It can be obviously seen that the experimental results show clear contrast between the HE lesions and periphrastic normal tissues with the gray scale of mass density, such as the early tiny HE cystic / necrosis lesions due to alveolar Echinococcus multilocularis infection indicated by the red arrows in Fig 5(b) and 5(c). However, the lesions cannot be found by MRI examination, and such a small change could not be clearly visualized by conventional absorption contrast imaging. These small changes in the density of tissues were depicted as changes in the apparent image contrast via the presented phase contrast imaging technique and the three-dimensional renderings of the HE specimen, shown in purple color in Fig 5(d) and 5(e) with different scale levels. This technique clearly revealed the entire volumetric endovascular HE lesion microstructure and distribution at an early stage due to the advantage of density resolution as depicted in Fig 3, demonstrating that the early tumor-like HE disease propagates along the microvascular blood networks. With careful visual inspection, there are some light dark areas, such as those indicated with a red rectangle in Fig 5(b), and they were examined by H&E staining microscopy for further verification, shown in Fig 6(a) and 6(b). It can be obviously found that there are some bright HE immune response rings, as indicated by red arrows in Fig 6(a). The enlarged inset of Fig 6(b), corresponding to the yellow rectangle in Fig 6(a), shows the early hepatic echinococcosis and the triggered cluster of microvesicles smaller than 2μmsurrounded by macrophages, epithelioid cells and lymphocyte infiltration.

Fig 5
Reconstructed results ofthe early HE specimen of the two-week rat model based on PAD-PPCT.
Fig 6
Histopathological examinations of the early HE specimens of the two-week rat modelwith thehematoxylin and eosin (H&E) staining process, corresponding to the light dark area indicated by the red rectangle in Fig 5(b), using 10× and 40× ...

Another experiment was performed with an HE specimen selected from the two-month rat model group to investigate the micro-morphological characteristics of developed HE based on our PAD-PPCT technique. Fig 7(a), 7(b) and 7(c) are reconstructed sectional images through the different orientations displayed on a mass density gray scale, showing the diverse HE lesions, such as microcysts (green arrows), micro-calcifications (white arrows), and inflammation areas (yellow arrows). Moreover, the 3D rendering of Fig 7(d) enables the visualization of different HE lesions with different colors vividly, such as micro-calcifications (white), fibrosis tissues (yellow), HE parenchyma (blue), E.multilocularis metacestodes and entire HE micro-cavities (green holes). Due to the high density sensitivity and striking contrast-to-noise ratio of our PAD-PPCT technique, 3Drenderings of quantitative segmentation, as shown in Fig 7(e) and 7(f), exhibited the volumetric distribution and structural characteristics of different HE lesions with small density differences, such as HE parenchyma rendering (e) and inflammation area rendering (f). In addition, Fig 8(a) and its enlarged inset (b) indicated with the yellow rectangle in Fig 8(a) show the results of histopathological examination using H&E staining of the HE-induced inflammation area, corresponding to the red rectangle in Fig 7(a). This figure clearly demonstrated that the developed cysts induced by hepatic echinococcosis triggered a strong HE immune response, agreeing with the inflammation areas in Fig 7(a)–7(c).

Fig 7
Reconstructed results ofthe developed HE specimen of the two-monthrat model based on PAD-PPCT.
Fig 8
Histopathological examinations of the early HE specimen of the two-monthrat model with thehematoxylin and eosin (H&E) staining process, corresponding to the area indicated by the red rectangle in Fig 7(a), using 10× and 40× objective ...

Analysis of contrast-to-noise ratio (CNR)

To quantitatively analyze the quality of the tomographic images obtained from the PAD-PPCT technique, we examined the density resolution (in view of the contrast-to-noise ratio: CNR). For the calculation of CNR based on these selected pairs of ROIs, as seen in Fig 3(a) and 3(b), the following formula was used:

CNR=|SobjSbg|σobj2+σbg2
(2)

where S represents the average of a homogeneous material, and σ represents its standard deviation. The manually defined object (red solid rectangles: I/PP, II/PS, III/PMMA) and background region (green dashed rectangle), indicated in Fig 3(a) and 3(b), are denoted with the subscripts obj and bg, respectively. The calculated results are the following: CNR_PPCT_I, II, and III equal to 4.35, 5.21, and 5.98 and CNR_PAD-PPCT_I, II, and III equal to45.61, 55.36, and 61.23, respectively. These results indicated that the bulk CNRs with the PAD-PPCT technique were more than 10 times higher than those with the direct PPCT method. This significant enhancement in the bulk CNR implies that the PAD-PPCT method has great potential for suppressing noise, enhancing image contrast, and reducing the radiation dose. Actually, our experimental data demonstrated that the PAD-PPCT method plays an important role in quantitative measurement and segmentation of HE specimens in this disease research.

Estimation of tissue's mass density

For the PAD-PPCT method, the reconstructed results generally describe the distribution of the object’s electron density, in which a pixel value denotes a specimen’s refraction index decrement δ. To obtain the mass density distribution of our HE specimens, the relationship of both mass density ρ and refraction index decrement δ can be calculated from the formula below:

δ=λ2re2πNaρwiZiAi
(3)

where λ, re, and NA, respectively, represent the X-ray wavelength, the classical electron radius (2.82×10-15m), and Avogadro’s constant (6.022×1023), wi is the weight fraction of the i-th element of a molecule, and Zi and Ai are the atomic number and atomic weight of the corresponding i-th element, respectively. For our estimation, the ratio value of Zi/Ai was assumed as 0.55; based on the average for soft tissues (referring to the NIST database),∑wi is actually equal to 1. According to the Eq (3), the density variance Δρ is expediently estimated from the reconstructed Δδ in X-ray PAD-PPCT. A measured Δδ value can be calibrated using standard specimens under a given experimental condition. Our experimentally measured Δδ values were multiplied by the calibration factor 1.5. Thus, the mass density distributions of the specimens from both HE rat model groups are displayed with mass density scalebars in Figs Figs55 and and77.

Quantitative segmentation and statistical analysis

Without the staining process and micron-level slice cutting of our experimental HE specimens, the quantitative pathological measurement of HE specimens based on the single-distance synchrotron PAD-PPCT technique with a simple experimental setup was investigated and demonstrated the potential of our presented method, especially for the quantitative segmentation of different types of HE lesions due to the striking CNR advantage and high density resolution of the PAD-PPCT method, as shown in Fig 7 (e) and 7(f). Therefore, it is feasible to obtain improved image contrast and higher spatial resolution to extract a type of HE lesions with tiny density changes from weak-absorption biomedical soft samples by using the robust phase retrieval based segmentation, which benefits visualization and quantitative measurement of inner pathological structures and lesion characterization in HE tissues. For example, it can visually supply are liable morphological indicator of disease activity, and its absence could indicate degeneration, revealing the parasitic metabolic activity at the early stage of pathology.

In addition, the statistical analysis of different HE lesion volumes based on the effective 3D segmentation can be obtained and is presented for three HE rat disease models (specimens I, II, and III), as shown in Fig 9, providing the quantitative estimations of the damage level of different types of HE lesions.

Fig 9
Statistical volumetric results of three types of HE lesions based on quantitative PAD-PPCT segmentation.

Conclusion

The PAD-PPCT imaging technique showed excellent density sensitivity and spatial resolution in 3D visualization of hepatic echinococcosis (HE) pathology based on our experimental results with rat HE disease models, effectively displaying enhanced phase contrast of early HE lesions of blood vascular invasion and different types of micro-morphological characteristics of HE lesions in mature rat models, which will play an important role in the pathological examination and assessment of therapeutic effects. Generally, it is helpful to quantitatively extract and analyze the morphological characteristics of soft tissues or organs closely related to the occurrence and development of HE or other diseases.

Acknowledgments

The authors thank the SSRF for the allocated beamtime at the X-ray imaging and biomedical application beamline (BL13W1) and appreciate the technical support from the duty staff of BL13W1for our experiment. This work was performed within the cooperative protocol between XJMU and SSRF.

Funding Statement

This work was supported by the National Natural Science Foundation of China (Grant No. 81660306) to HQL, Postdoctoral Science Foundation of Xinjiang Medical University, and the Technology Foundation for Selected Overseas Chinese Scholar. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Data Availability

Data Availability

All relevant data are within the paper.

References

1. Donald PM, Li Z, Yang S. Case studies emphasising the difficulties in the dagnosis and management of alveolar echinococcosis in rural China. Parasites & Vectors. 2011;4:196. [PMC free article] [PubMed]
2. Czermak BV, Unsinn KM, Gotwald T, Waldenberger P, Freund MC, Bale RJ, et al. Echinococcus multilocularris revisited. AJR Am J Roentgenol. 2001;176, 1207–1212. doi: 10.2214/ajr.176.5.1761207 [PubMed]
3. Laura C, Elissondo C, Moreno L, Dopchiz M, Sanchez BS. Albendazole treatment in cystic echinococcosis: pharmacokinetics and clinical efficacy of two different aqueous formulations.Parasitol Res.2008;103: 355–362. doi: 10.1007/s00436-008-0980-x [PubMed]
4. Kern P, Bardonnet K, Renner E, Auer H, Pawlowski Z.European echinococcosis registry: human alveolar echinococcosis Europe 1982-2000.Emerg Infect Dis. 2003;3:343–349. [PMC free article] [PubMed]
5. Yang YR, Craig PS, Ito A, Vuitton DA, Giraudoux P. A correlative study of ultrasound with serology in an area in China co-enemic for human alveolar and cystic echinococcosis. Trop Med Int Health. 2007; 12:637–646. doi: 10.1111/j.1365-3156.2007.01834.x [PubMed]
6. Bresson-Hadni S, Delabrousse E, Blagosklonov O, Bartholomot B, Koch S. Imaging aspects and non-surgical interventional treatment in human alveolar echinococcoisis. ParasitolInt. 2006; 55:S267–S272. [PubMed]
7. Caoduro C, Porot C, Vuitton DA, Bressonhadni S, Grenouillet F. The role of delayed 18F-FDG PET imaging in the follow-up of patients with alveolar echinococcosis. J. Nucl. Med.2013;54:1–6. [PubMed]
8. Amel A, Oleg B, Ahmed L, Louis B, Vuitton DA, Solange BH, et al. Alveolar echinococcosis: correlation between hepatic MRI findings and FDG-PET/CT metabolic activity.Abdominal Imaging. 2015;40:56–63. doi: 10.1007/s00261-014-0183-0 [PubMed]
9. Zeng H, Wang J, Xie WD, Liu WY, Wen H.Assessment of early hepatic Echinococcus Multilocularis Infection in Rats with real-time contrast-enhanced ultrasonography.Ultrasound in Med. & Biol.2012;38(11):1982–1988. [PubMed]
10. Piscaglia F, Lencioni R, Sagrini E, Pina CD, Cioni D. Characterization of focal liver lesions with contrast-enhanced ultrasound, Ultrosound in Med. & Biol. 2010; 36: 531–550. [PubMed]
11. Reuter S, Gruner B, Buck AK, Blumstein N, Kerm P.Long-term follow-up of metabolic activity in human alveolar echinococcosis using FDG-PET.Nuklearmedizin. 2008;4: 147–152. [PubMed]
12. Coskun A, Ozturk M, Karahan OI, Erdogan N, Isin S, Gulec M. (2004) Alveolar echinococcosis of the liver: correlative color Doppler US, CT, and MRI study.ActaRadiol. 2004;45: 492–498. [PubMed]
13. Groso A, Stampanoni M, Abela R, Schneider P, Linga S et al. (2006)Phase contrast tomography: An alternative approach. APPLIED PHYSICS LETTERS 88, 214104.
14. Nesterets YI, Wilkins SW, Gureyev TE, Pogany A, Stevenson AW. On the optimization of experimental parameters for X-ray in-line phase contrast imaging. REVIEW OF SCIENTIFIC INSTRUMENTS. 2005;76: 093706.
15. Lewis RA.Medical phase contrast X-ray imaging: current status and future prospects. PHYS. MED. BIOL. 2004; 49: 3573–3583. [PubMed]
16. Cloetens P, Mache R, Schlenker M,Lerbs MS. Quantitative phase tomography of Arabidopsis seeds reveals intercellular void network. PNAS. 2006; 103: 14626–14630. doi: 10.1073/pnas.0603490103 [PubMed]
17. Lovric G, Peter O, Mohacsi I, Stampanoni M, Mokso R. A robust tool for photon source geometry measurements using the fractional Talbot effect. OPT. EXPRESS. 2014; 22: 2745 doi: 10.1364/OE.22.002745 [PubMed]
18. Revol V, Kottler C, Kaufmann R, Straumann U, Urban C.Noise analysis of grating-based X-ray differential phase contrast imaging. REVIEW OF SCIENTIFIC INSTRUMENTS. 2010; 81: 073709 doi: 10.1063/1.3465334 [PubMed]
19. Pfeiffer F, David C, Bunk O, Donath T, Bech M.Region-of-Interest Tomography for Grating-based X-ray Differential Phase-Contrast Imaging. PHYSICAL REVIEW LETTERS. 2008;101: 168101 doi: 10.1103/PhysRevLett.101.168101 [PubMed]
20. Pfeiffer F, Bech M, Bunk O, Kraft P, Eikenberry EF.Hard-X-ray dark-field imaging using a grating interferometer. NATURE materials. 2008;7: 134 doi: 10.1038/nmat2096 [PubMed]
21. Bech M, Jensen TH, Bunk O, Donath T,David C. Advanced contrast modalities for X-ray radiology: Phase-contrast and dark-field imaging using a grating interferometer. Z. Med. Phys. 2010; 20: 7–16. doi: 10.1016/j.zemedi.2009.11.003 [PubMed]
22. Zanette I, Lang S, Rack A, Dominietto M.Holotomography versus X-ray grating interferometry: A quantitative study. Appl. Phys. Lett. 2013; 103: 244105.
23. Lang S, Zanette I, Dominietto M, Langer M, Rack A, David C,et al. Experimental comparison of grating- and propagation-based hard X-ray phase tomography of soft tissue. JOURNAL OF APPLIED PHYSCIS. 2014; 116: 154903(1–12).
24. Zhou T, Larsson DH, Stampanoni M, Hertz HM. Comparison of propagation-based X-ray phase-contrast imaging techniques with a liquid-metal-jet source. Proceedings of SPIE. 2014; 9033(1): 131–135.
25. Liu y, Nelson J, Holzner C, Andrews JC, Pianetta P.Recent advances in Synchrotron-based hard X-ray phase contrast imaging. J. Phys. D: Appl. Phys. 2013; 46: 494001.
26. David MP, Myers GR, Gureyev TE. Stability of phase-contrast tomography. J. Opt. Soc. Am. A. 2007;24: 2516. [PubMed]
27. Groso A, Abela R, and Stampanoni M.Implementation of a fast method for high resolution phase contrast tomography. Opt. Express. 2006; 14: 8103–8110.; [PubMed]
28. Paganin D, Mayo SC, Gureyev TE, Miller PR, Wilkins SW. Simultaneous phase and amplitude extraction from a single defocused image of a homogeneous object. J. Microscopy.2002; 206: 33–40. [PubMed]
29. Langer M, Cloetens P, and Peyrin F.Regularization of phase retrieval with phase-attenuation duality prior for 3-D Holotomography. IEEE Trans. Image-Process, 2010;19(9): 2428–2436. doi: 10.1109/TIP.2010.2048608 [PubMed]
30. Gureyev TE, Paganin DM, Myers GR,Nesterets YI, Wilkins SW.Phase-and-amplitude computer tomography. Appl. Phys. Lett.2006;89: 034102.
31. Guigay JP, Langer M, Boistel R, Cloetens P.Mixed transfer function and transport of intensity approach for phase retrieval in the Fresnel region. Opt. Lett. 2007; 32: 1617–1619. [PubMed]
32. Gureyev TE, Nesterets YI, Paganin DM, Pogany A, Wilkins SW.Linear algorithms for phase retrieval in the Fresnel region. Optics Communications 2004; 231: 53–70.
33. Cloetens P, Ludwig W, Baruchel J, Dyck DV, Landuyt JV.Holotomography: Quantitative phase tomography with micrometer resolution using hard synchrotron radiation x rays. Appl. Phys. Lett. 1999; 75: 2912.
34. Thuering T, Modregger P, Grund T, Kenntner J, David C.High resolution, large field of view X-ray differential phase contrast imaging on a compact setup. APPLIED PHYSICS LETTERS. 2011;99: 04111.
35. Engelhardt M, Baumann J, Schuster M, Kottler C, Pfeiffer F.High-resolution differential phase contrast imaging using a magnifying projection geometry with a microfocus X-ray source. APPLIED PHYSICS LETTERS. 2007;90: 224101.
36. Satoshi T, Kazumasa H, Taro K, Yamamoto Y, Yoneyama A.Imaging and density mapping of tetrahydrofuran clathrate hydrates by phase-contrast X-ray computed tomography. APPLIED PHYSICS LETTERS. 2007;90: 081920.
37. Wu XZ, and Liu H. A general theoretical formalism for X-ray phase contrast imaging. Journal of X-ray Science and technology. 2003;11: 33–42. [PubMed]
38. Wu XZ, Liu H, and Yan AM.X-ray phase-attenuation duality and phase retrieval. Opt. Lett. 2005; 30: 379–381. [PubMed]
39. Wu XZ, Liu H, and Yan AM.Phase-contrast X-ray tomography: Contrast mechanism and roles of phase retrieval.EJR. 2008; 68S: S8–S12. [PubMed]
40. Wu XZ, and YAN AM.Phase retrieval from one single phase contrast x-ray image. Opt. Express. 2009; 17: 11187–11196. [PubMed]
41. Liu HQ, Wu XZ, andXiao TQ. Optimization of reconstructed quality of hard X-ray phase micro-tomography. APPLIED OPTICS. 2015;54:5610 doi: 10.1364/AO.54.005610 [PMC free article] [PubMed]
42. Liu HQ, Wu XZ, and Xiao TQ. Technical Note: Synchrotron-based High Energy X-ray Phase Sensitive Micro-tomography for Biomedical Research. Med. Phys. 2015; 42(10): 5595–5603. doi: 10.1118/1.4929551 [PubMed]
43. Liu HQ, Ren YQ, Guo H,Xue YL, Xie HL, Wu XZ,et al. Phase retrieval for hard X-ray computed tomography of samples with hybrid compositions. Chin. Opt. Lett. 2012; 10: 121101 doi: 10.3788/COL201210.121101 [PMC free article] [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science