Skip to main content

A sparse-projection computed tomography reconstruction method for in vivo application of in-line phase-contrast imaging

Abstract

Background

In recent years, X-ray phase-contrast imaging techniques have been extensively studied to visualize weakly absorbing objects. One of the most popular methods for phase-contrast imaging is in-line phase-contrast imaging (ILPCI). Combined with computed tomography (CT), phase-contrast CT can produce 3D volumetric images of samples. To date, the most common reconstruction method for phase-contrast X-ray CT imaging has been filtered back projection (FBP). However, because of the impact of respiration, lung slices cannot be reconstructed in vivo for a mouse using this method. Methods for reducing the radiation dose and the sampling time must also be considered.

Methods

This paper proposes a novel method of in vivo mouse lung in-line phase-contrast imaging that has two primary improvements compared with recent methods: 1) using a compressed sensing (CS) theory-based CT reconstruction method for the in vivo in-line phase-contrast imaging application and 2) using the breathing phase extraction method to address the lung and rib cage movement caused by a live mouse’s breathing.

Results

Experiments were performed to test the breathing phase extraction method as applied to the lung and rib cage movement of a live mouse. Results with a live mouse specimen demonstrate that our method can reconstruct images of in vivo mouse lung.

Conclusions

The results demonstrate that our method could deal with vivo mouse’s breathing and movements, meanwhile, using less sampling data than FBP while maintaining the same high quality.

Background

Lung disease has become increasingly common. Lung diseases range from the common cold to life-threatening conditions such as bacterial pneumonia, pulmonary embolisms, and lung cancer. The lungs are the essential respiration organs in many air-breathing animals. Their principal function is to transport oxygen from the atmosphere into the bloodstream and to release carbon dioxide from the bloodstream into the atmosphere. Lung tissue is a typical soft tissue that has received extensive attention from many research groups to reveal its detail structures.

X-rays [1] have been widely used to reveal the detailed structures of bodies in medical applications. However, it is difficult to image soft tissue with X-rays. Imaging the lung tissue results in a low absorption contrast because the lung’s anatomical structure has multiple air-tissue alveoli. Phase-contrast imaging (PCI) is a new imaging mechanism intended to solve the problem of soft tissue imaging. PCI is based on using X-rays’ phase change through objects to obtain a phase contrast image. For hard X-rays (10–100 keV), the X-ray absorption of soft tissues is low, but the phase change caused by soft tissue is approximately 1000 times. Therefore, PCI performs better than absorption imaging and substantially increases the X-ray contrast resolution of soft tissue imaging. Currently there are three primary types of techniques that have been used for PCI: X-ray interferometry [2], diffraction-enhanced imaging (DEI) [3] and in-line holography [4].

The in-line phase contrast X-ray imaging method was developed by Snigirev et al. at European Synchrotron Radiation Facility (ESRF) [5] and by Wilkins et al. at Commonwealth Scientific and Industrial Research Organization (CSIRO) [6]. In-line holography of PCI is a powerful phase-sensitive technique that generates high spatial resolution and super contrast of weakly absorbing objects compared with conventional radiography by using a propagation-based technique. In this method, X-rays are transmitted through the sample object at various angles and then received by the detector within a certain distance. When the X-rays propagate through the sample, the shape of the wave-front changes because the varying thickness of the sample causes differences in the X-ray refractive index, then, an edge-enhanced contrast is observed.

Phase contrast X-ray imaging methods have been applied to mouse’s lung [79] and propagation-based methods have been applied to mouse’s lung and exhibit a speckled intensity pattern that is attributed to the air-filled alveoli [10]. For three-dimensional observation of lung tissues using this method, phase contrast X-ray computed tomography (PCX-CT) [11] has been developed. PCX-CT was achieved by introducing the technique of X-ray phase contrast imaging into X-ray CT. It has the characteristics of high sensitivity and three-dimensional imaging of lung tissues [12]. A 3D lung model is built, and the experiment is implemented in a dead mouse. It remains challenging to implement for a live mouse. During normal breathing, expiration is passive and no muscles are contracted (the diaphragm relaxes). The rib cage itself is also able to expand and contract to some degree through the action of other respiratory and accessory respiratory muscles. Consequently, air is transported into or expelled out of the lungs. This type of lung is known as a bellows lung, and it causes difficulty in lung tissue registration and alignment at various angles.

In this paper, we present a compressed sensing-based (CS-based) optimization method for CT reconstruction and apply this method to an in vivo in-line phase contrast imaging experiment. The proposed method has two primary contributions. 1) It addresses the problems caused by a live mouse’s respiration. Respiratory phases are classified and extracted into different stages. For example, the rib cage’s two stages are expansion to maximum and contraction to minimum. 2) It addresses the sparse 3D reconstruction problem. Respiratory phase extraction causes the number of images to be reduced at each angle. A compressed sensing-based (CS-based) optimization method for CT reconstruction is applied to address this problem. This research is also important for performing CT image reconstruction from a small number of projections, which can shorten the scan time and reduce the radiation dose.

Methods

Optical system

The basic theory of in-line phase-contrast X-ray imaging is Fresnel diffraction. X-rays propagate through the sample object and are received by the detector within a certain distance. During the propagation, the shape of the X-ray wave front changes because of variations in the thickness and the X-ray refractive index of the sample [13].

Figure 1 shows a schematic of an in-line phase contrast imaging configuration. In Figure 1, R1 is the source-to-object distance and R2 is the object-to-detector distance (In the experiment, R1 is set to 34 m, R2 is set to 1.2 m). To achieve phase-contrast imaging based on diffraction effects, the X-ray source should provide a sufficient degree of spatial coherence. In the X-ray region, a wave field propagating in the z direction is expressed under the paraxial approximation by

2 π λ z I x , y , z = I x , y , z Φ x , y , z ,
(1)

which is known as the transport of intensity equation (TIE) [14, 15], where is the two-dimensional gradient operator acting in the x-y plane. In the near-field Fresnel region, TIE can address the change of the wave front caused by propagation as Φ(x, y, z) in eq. (1).

Figure 1
figure 1

Schematic of an in-line phase contrast imaging configuration.

In the case of weak absorption and unit-amplitude plane-wave illumination, eq. (1) is simplified to

I x , y , z 1 + λz 2 π 2 Φ x , y , 0
(2)

Where z is the distance of the detector to the object, λ is the wavelength, I(x, y, z) is the detected intensity on the detector plane, and Φ(x, y, 0) is the object’s phase function. Eq. (2) establishes a linear relation between the Laplace transform of the phase function and the measured intensity data. It enables generation of contrast-outlining surfaces and structural boundaries, where the refractive index changes abruptly. The process of phase retrieval is to reformulate Eq. (2) to obtain the decrement of the real part of the object’s refractive index δ(x, y, 0) from knowledge of I(x, y, z).

The phase retrieval method we used is the error reduction algorithm [16], which consists of the following four steps: 1) Forward Fresnel transform an estimate of the complex amplitude of the object plane. 2) Replace the modulus of the resulted Fresnel transform with the measured intensity to form an estimate of the complex amplitude of the detector plane. 3) Backward Fresnel transform the estimated complex amplitude of the detector plane. 4) Constrain the amplitude to be zero out of the object zone and then form a new estimate of the complex amplitude of the object plane.

Define the complex amplitude as U(x, y, z) = ρ(x, y, z)exp[(x, y, z)]. The relation between the intensity I(x, y, z) and the complex amplitude U(x, y, z) is I(x, y, z) = |U(x, y, z)|2.

Respiratory phase extraction and alignment

Lung and rib cage movements are caused by a live mouse’s breathing, there are two ways to deal with this problem: one is the use of animal ventilator to control rat’s breathing rhythm; the other is anesthetized spontaneous breathing, and respiratory rhythm parameters are applied to extract the same phase of breathing images, and then rebuilt the 3D model. In this paper, we utilize the second way. The proposed alignment method is fully automatic and it includes two steps: 1) design a respiratory phase extraction algorithm that can divide the images into different phases and extract two stages. 2) aligning the stable region between each two images within the same respiration stage. So the gating information is not relying on threshold but our 2-step alignment algorithm.

In a live mouse, respiration is passive and the rib cageexpands and contracts to some degree through the action of other respiratory and accessory respiratory muscles. We took 20 images at each angle, which sampled several respiratory cycles of the live mouse.

As shown in Figure 2, there are two respiratory cycles in one angle. The aim is to design a respiratory phase extraction algorithm that can divide the images into different phases and extract two stages, expanded to maximum and contracted to minimum. Two parameters, AP (average parameter) and MP (mean parameter), are calculated as

AP = mask I x 1 I x 2 N mask
MP = max mask I x 1 I x 2
(3)
Figure 2
figure 2

Images of two respiratory cycles viewed from one angle.

x1 and x2 are in a silhouette mask, which has nonzero pixels where the motion occurs. Respiratory cycle curves can be drawn after the calculation of the two parameters, as shown in Figure 3. At each angle, AP and MP are consistent, and respiratory phases are divided by AP and MP. Two stages, expand to maximum and contract to minimum, are extracted.

Figure 3
figure 3

Respiratory cycles curve for one angle.

Alignment of images taken at different angles is the next step. As shown in Figure 4, images taken at different angles have similar covariant characteristics if the pixels in the two images that correspond to the partial image region satisfying the following mapping:

x = H s x = sR θ t 0 T 1 x
(4)
Figure 4
figure 4

Respiratory phase images extracted at different angles. These images are called projection image which are taken at different angle from 1° to 180° and there are 180 images to be aligned.

Where x = (x, y, 1)T and x′ = (x′, y′, 1)T denote the corresponding pixels in the two images’ homogeneous coordinates,t denotes the translation vector, R(θ) indicates the parameter that is a rotation matrix of the θ, and the scalar s indicates a scaling factor. The similarity transformation has four degrees of freedom (two translation parameters, a rotation parameter and a scale parameter); therefore, the shape remains after the similarity transformation.

Aligning the stable region between each two images separated by steps of 1° is important for aligning all the images together. The scheme is implemented as follows:

  1. 1)

    Use a 2D Gaussian kernel to continuously smooth the given gray scale image (if the input is a color image, it is necessary to pre-convert to gray scale images). Use down-sampling to construct the multi-scale representation, namely, the Gaussian pyramid. In the Gaussian pyramid for each image, by smoothing using the Gaussian kernel function, the parameter is k = 2. The adjacent images within each layer of the Gaussian pyramid are subtracted to obtain the DoG pyramids.

  2. 2)

    Solve for local extrema in the DoG scale space. Then, according to principal curvatures, extreme points in the instability of the position of the image edge are removed to obtain a stable (x, y, s). Here, (x, y) represents the position of the stability of extreme points and s represents the scale of extreme points.

  3. 3)

    Utilize the Scale Invariant Feature Transform (SIFT) feature extraction [17, 18] method for each of the stable extreme points.

  4. 4)

    Use an exhaustive search method to solve the nearest neighbour matching relationship between two images. Then, the geometric constraints between the two image models and the RANdom SAmple Consensus (RANSAC) algorithm [19] are combined to exclude incorrect candidate matching pairs, and the geometric mapping relationship between two images is estimated.

Compressed sensing-based CT reconstruction method

After the respiratory phase extraction and alignment scheme, images taken at different angles are aligned and prepared for reconstruction. The number is sparse: in this paper, each angle has two images (expand to maximum and contract to minimum). The reconstruction methods share the idea of compressed sensing (CS). The main premise of compressed sensing (CS) [2022] is that although the signal is not necessarily sparse in real space or Fourier space, it is sparse or compressible in some basis. If we consider a real signal X = {x i }1≤iN and define a real basis Ψ = {ψi,j}1 ≤ i ≤ N;1 ≤ j ≤ T, then we can say that the decomposition α = {α j }1≤jT is the column vector of weighting coefficients α j . The signal X can be expressed as

X = Ψα
(5)

The signal X is sparse or compressible if there are only K < < N weighting coefficients that are nonzero or significant.

Suppose that each pixel is denoted by p t 1 , t 2 , 0 t 1 , t 2 < N , and at the same position (t1, t2), the pixel in images is denoted by q(t1, t2). Then,

q t 1 , t 2 = p t 1 + 1 , t 2 p t 1 , t 2 2 + p t 1 , t 2 + 1 p t 1 , t 2 2 .
(6)

With any N × N digital image p, define a real number

p TV = t 1 = 0 N 2 t 2 = 0 N 2 p t 1 + 1 , t 2 p t 1 , t 2 2 + p t 1 , t 2 + 1 p t 1 , t 2 2 .
(7)

This is known as the image’s total variation (TV) [23]. We define that the position of (t1 + 1, t2) is ‘below’ and (t1, t2 + 1) is ‘to the right’ of (t1, t2).

We then probe the signal using M real linear measurements Y = {y r }1≤rM in some sensing basis Φ = {φr,i}1 ≤ r ≤ M;1 ≤ i ≤ N such that

Y = Φ X = ΦΨ α = Θ α ,
(8)

where Θ = ΦΨ RM × T. The measurement process is not adaptive, which means that Ф is fixed and does not depend on the signal X. The image processing procedure consists of designing 1) a stable measurement matrix Ф and 2) a reconstruction algorithm to recover X from only M measurements Y.

The measurement matrix Ф must allow the reconstruction of the length-N signal X from M < N measurements (the vector Y). If X is K-sparse and the K locations of the nonzero coefficients in αare known, the problem can be solved with the provided M ≥ K. A necessary and sufficient condition for this simplified problem is that for any vector α K that shares the same K nonzero entries as αand for some 0 < δ K  < 1, 

1 δ K α K 2 2 Θ α K 2 2 1 + δ K α K 2 2 .
(9)

Defining the l p -norm of vector u, u p  = ( ∑ |u i |p)1/p. Usually, the locations of the K nonzero entries in αare unknown. However, a sufficient condition for a stable solution for both K-sparse and compressible signals is that Θ satisfies Eq. (7) for an arbitrary 3 K-sparse vector α K . This condition is referred to as the restricted isometry property (RIP) [24]. A related condition, referred to as incoherence, requires that the rows {φ r } of Ф cannot sparsely represent the columns {ψ j } of Ѱ, and vice versa.

The reconstruction task is as follows: suppose that X is an unknown image and that we are given M measurements Y. We must find an image X * that is a ‘good’ approximation of X. This X * has the minimal TV, defined in Eq. (6), of all those that satisfy the constraint of Eq. (7).

min X * TV s . t . Y ΦX * 2 ϵ .
(10)

The method used to minimize X *  TV is the gradient descent algorithm. The general form of the gradient descent algorithm is

x next = x current + β current υ current ,
(11)

where β current is a positive real number that denotes the current step length. For a digital image X and any kN, let sk ∂ X TV be a sub-gradient of X TV at xk and

υ k = s k s k , s k 0 0 , s k = 0 .
(12)

We now describe the iterative steps of this reconstruction algorithm for phase sensitive CT images. Denote the digital image estimate by Xk, its ith pixel by x i k, the measurement by Y, and the measurement matrix by Ф. The algorithm is performed according to the following pseudo code:

  1. 1.

    X 0 = 0; β = 1; 

  2. 2.

    x i 1 = x i 0 + λ ϕ r , i Y r i = 1 N ϕ r , i x i 0 i = 1 N ϕ r , i 2

  3. 3.

    for k = 1 , , K

  4. 4.

    if (X k TV  < X k−1 TV )

  5. 5.

    x i k = x i k + βυk; 

  6. 6.

    x i k + 1 = x i k + λ ϕ r , i Y r i = 1 N ϕ r , i x i k i = 1 N ϕ r , i 2 ;

  7. 7.

    if i x i k + 1 x i k < ϵ

  8. 8.

    break;

  9. 9.

    else

  10. 10.

    β = 0.8β; 

  11. 11.

    end if

  12. 12.

    else

  13. 13.

    β = 0.8β; 

  14. 14.

    x i k = x i k + βυk

  15. 15.

    end if

ϵ is a small positive value to control the iterative process. To avoid infinite loops, we also introduced an additional stopping criterion based on the maximum loop iterations, K. When the iterative process is stopped, the phase retrieval process is begun. These parameter values are experimentally set, β = 0.8β is experimentally tested. In our experiment, gradient descent can start the process of decline more rapidly, as extreme approach, the rate of decline to slow down to make sure it converges to the optimum. λ is in the range [0, 2]. It is called relaxation factor, which is affecting convergence speed and reconstruction quality. There is a large number of specialized research focus on this topic. In this paper, relaxation factor values is in [0, 2] range. The bigger of Relaxation factor, the faster iteration convergence, but the result is more unstable. The smaller of Relaxation factor, the slower iteration convergence, but can guarantee convergence to a better optimum.

Results

Experiment implementation

The experiments were performed at the BL13W1 beamline at the Shanghai Synchrotron Radiation Facility (SSRF) in China. The light source was a hybrid-type wiggler with periodic length of 14 cm and period number 8. The radiating power was varied from 8.0 to 72.5 keV by tuning the gap from 17 mm to 35 mm. A fixed-exit double-crystal cryogenic-cooling monochromator was placed 28 m away from the source. The monochromator crystal was a combination of an Si(111) orientation crystal and an Si(311) orientation crystal. A high-precision sample platform was used to position the sample and rotate the sample axis perpendicular to the beam. The space resolution of the sample platform was better than 1 micron. Figure 5 shows a schematic of the BL13W1 beamline at the SSRF [24].

Figure 5
figure 5

Schematic of the BL13W1 beamline at the SSRF.

The specimen was a 4-week-old live mouse. It was fixed to the sample rotation platform. The edge-enhancement images were measured at 180 different angular settings of the specimen. The white beam X-ray (wide band photoelectron spectroscopy X-ray) was monochromatized at 25 keV using a silicon double-crystal monochromator. The distance from the sample to the detector could be as great as 1.2 m, which would achieve the greatest phase shift and the best image in this experiment. The sample was placed on a high-precision sample platform and was scanned rotationally at 1° per step over 180°. The exposure time was 24 ms for each scanned step. The transmitted radiation was detected by the X-ray CCD (13 μm/pixel). All experiments and procedures carried out on the animals were approved by the animal welfare committee of Capital Medical University and the approval ID is SCXK-(Army) 2007–004.

Respiratory phase extraction and alignment

Respiratory cycle curves can be drawn after the calculation of the two parameters. Two stages, expand to maximum and contract to minimum, are extracted. Alignment of images taken at different angles is the next step (Results shown in Figure 6). Images taken at different angles have similar covariant characteristics. In this paper, 180 images from 1° to 180° are aligned. The performance of alignment algorithm is listed in Table 1. Here we take 1° to 5° for example to illustrate the numerical results of time and number of matching points.

Figure 6
figure 6

Alignment results for two images at different angles, 0° and 1.

Table 1 Alignment results

Reconstruction results of in vivo mouse lung

Figure 7 shows three in-line phase-contrast projection images of the mouse’s chest taken at different angles. Figure 8 shows the CT reconstructed images for one slice reconstructed using two algorithms under three sampling conditions. All the images come from the same slice, and the total number of slices is 105. The images are (a) the CT image reconstructed via the CS-based algorithm from 60 views, (b) the CT image reconstructed via the CS-based algorithm from 180 views, (c) the CT image reconstructed via the CS-based algorithm from 30 views, and (d) the CT image reconstructed via the FBP algorithm from 180 views. From the results, we can see that the CT images reconstructed via the CS-based algorithm have better quality than the CT images from the FBP algorithm. Even the CS-based method with less views, such as 30, has better performance than the FBP method with more views, such as 180. This result implies that the proposed novel CS-based CT reconstruction method could be applied with a lower exposure time and dose. Our algorithm could address the respiration of a live mouse and reconstruct the 3D rib cage using very sparse image data.

Figure 7
figure 7

Three in-line phase-contrast experimental projection images of the mouse’s chest taken from several different angles.

Figure 8
figure 8

Alignment results for two images at different angles, 0° and 1. The CT reconstructed image of the 70th slice using (a) the CS-based method from 60 projection images, (b) the CS-based method from 180 projection images, (c) the CS-based method from 30 projection images, and (d) the FBP method from 180 projection images.

Conclusions

A novel compressed sensing (CS) theory-based CT reconstruction method was proposed, and an experiment was performed to test the breathing phase extraction method as applied to the lung and rib cage movement of a live mouse. The results demonstrate that our method can reconstruct the images using less sampling data than FBP while maintaining the same high quality.

Discussion

Several topics must be further studied and discussed:

  1. 1)

    The alignment of images taken at different angles requires further study. This paper could address the respiration of a live mouse by respiratory phase extraction. However, lung tissue reconstruction requires further study on the alignment topic because a more accurate alignment algorithm would result in more robust 3D reconstruction of lung tissue.

  2. 2)

    The algorithm remains time-consuming, which presents a difficulty for clinical applications. Methods for designing a fast and real-time in-line phase-contrast X-ray CT system are important and require further research.

Abbreviations

ILPCI:

In-line phase-contrast imaging

CT:

Computed tomography

FBP:

Filtered back projection

CS:

Compressed sensing

PCI:

Phase-Contrast imaging

DEI:

Diffraction-enhanced imaging

ESRF:

European synchrotron radiation facility

CSIRO:

Commonwealth scientific and industrial research organization

PCX-CT:

Phase contrast x-ray computed tomography

TIE:

Transport of intensity equation

SIFT:

Scale invariant feature transform

RANSAC:

RANdom sample consensus

SSRF:

Shanghai synchrotron radiation facility.

References

  1. Wall BF, Hart D: Revised Radiation Doses for Typical X-Ray Examinations. Br J Radiol 1997, 70: 437–439.

    Article  Google Scholar 

  2. Momose A, Takeda T, Hirano K: Phase-contrast x-ray computed tomography for observing biological soft tissues. Nat Med 1996, 2: 473–475. 10.1038/nm0496-473

    Article  Google Scholar 

  3. Davis TJ, Gao D, Gureyev TE, et al.: Phase-contrast imaging of weakly absorbing materials using hard x-rays. Nature 1995, 373: 595–597. 10.1038/373595a0

    Article  Google Scholar 

  4. Wilkins SW, Gureyev TE: Phase-contrast imaging using polychromatic hard X-rays. Nature 1996, 384: 335–338. 10.1038/384335a0

    Article  Google Scholar 

  5. Snigirev A, Snigireva I, Kohn V, Kuznetsov S, Schelokov I: On the possibilities of X-ray phase contrast microimaging by coherent high-energy synchrotron radiation. Rev Sci Instrum 1995, 66(12):5486–5492. 10.1063/1.1146073

    Article  Google Scholar 

  6. Wilkins SW, Gureyev TE, Gao D, Pogany A, Stevenson AW: Phase-contrast imaging using polychromatic hard X-rays. Nature 1996, 384: 335–338. 10.1038/384335a0

    Article  Google Scholar 

  7. Hwan Seok Y, Eun-Young K, Yoon Kyung K, et al.: Phase Contrast Microradiography of Mouse Lung Using Synchrotron X-ray: Correlation with Optical Microscopy. Yonsei Med J 2009, 50(3):422–426. 10.3349/ymj.2009.50.3.422

    Article  Google Scholar 

  8. Schleede S, Meinel FG, Bech M, et al.: Emphysema diagnosis using X-ray dark-field imaging at a laser-driven compact synchrotron light source. Proc Natl Acad Sci USA 2012, 109(44):17880–17885. 10.1073/pnas.1206684109

    Article  Google Scholar 

  9. Leong AFT, Andreas Fouras M, Sirajul I, et al.: High spatiotemporal resolution measurement of regional lung air volumes from 2D phase contrast x-ray images. Med Phys. 2013, 40(4):041909. 10.1118/1.4794926

    Article  Google Scholar 

  10. Kitchen MJ, et al.: On the origin of speckle in x-ray phase contrast images of lung tissue. Phys Med Biol 2004, 49(18):4335–4348. 10.1088/0031-9155/49/18/010

    Article  Google Scholar 

  11. Takeda T, Momose A, Ueno E, Itai Y: Phase-contrast X-ray CT image of breast tumor. J Synchrotron Radia 1998, 5: 1133–1135. 10.1107/S0909049597020116

    Article  Google Scholar 

  12. Zhang L, Li D, Luo S: Non-invasive microstructure and morphology investigation of the mouse lung: Qualitative description and quantitative measurement. PLoS One 2011, 6(2):e1 7400.

    Article  Google Scholar 

  13. Chen B, Zhu PP, Chen CC, Shu H, Liu YJ, Yuan QX, Wang JY, Huang WX, Ming H, Wu ZY: Theory and experiment of in-line phase contrast imaging on non-uniformly distributed source. Spectrochimica Acta 2007, 62(6–7):6,636–6,641.

    Google Scholar 

  14. Teague MR: Deterministic phase retrieval: A Green's function solution. J Opt Soc Am 1983, 73(11):1434–1441. 10.1364/JOSA.73.001434

    Article  Google Scholar 

  15. Siu KKW, Morgan KS, Paganin DM, Boucher R, Uesugi K, Yagi N, Parson DW: Phase contrast X-ray imaging for the non-invasive detection of airway surfaces and lumen characteristics in mouse models of airway disease. Eur J Radiol 2008, 68S: S22-S26.

    Article  Google Scholar 

  16. Li YH, Li Z, Liu Z, Zhang L: Phase retrieval of X-ray in-line phase con trast imaging based on iterative algorithm. CT Theory and Applications 2010, 19(2):53–60.

    Google Scholar 

  17. Mikolajczyk K, Tuytelaars T, Schmid C, et al.: A comparison of affine region detectors. Int J Comput Vis 2005, 65: 43–72. 10.1007/s11263-005-3848-x

    Article  Google Scholar 

  18. Mikolajczyk K, Schmid C: A performance evaluation of local descriptors. IEEE Trans Pattern Anal Mach Intell 2005, 27: 1615–1630.

    Article  Google Scholar 

  19. Fischler MA, Bolles RC: Random sample consensus: a paradigm for model fitting with application to image analysis and automated cartography. Commun ACM 1981, 24: 381–395. 10.1145/358669.358692

    Article  MathSciNet  Google Scholar 

  20. Donoho DL: Compressed sensing. IEEE Trans. on Information Theory 2006, 52: 1289–1306.

    Article  MathSciNet  Google Scholar 

  21. Candès EJ: Compressive sampling. Spain: In Proceedings of International Congress of Mathematicians. Madrid; 2006:1–20.

    Google Scholar 

  22. Baraniuk RG: Compressive Sensing. IEEE Signal Proc Mag 2007, 24(4):118–121.

    Article  Google Scholar 

  23. Rudin LI, Osher S, Fatemi E: Nonlinear total variation based noise removal algorithms. Physica D 1992, 60: 259–268. 10.1016/0167-2789(92)90242-F

    Article  Google Scholar 

  24. Candès EJ, Romberg J, Tao T: Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information. IEEE Trans. on Information Theory 2006, 52: 489–509.

    Article  Google Scholar 

Download references

Acknowledgements

This study was supported by the National Natural Science Foundation of China, Grant Nos. 61227802, 60532090 and 30770593, and by the 7th Framework Programme of the European Community, Grant Agreement No. PIRSES-GA-2009-269124. The authors thank the staff of the BL13W1 beamline of the Shanghai Synchrotron Radiation Facility in China for their kind assistance with our experiments.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Shuqian Luo.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

LW carried out the respiratory phase extraction and alignment studies, participated in the testing of algorithm performance and drafted the manuscript. XL carried out the CS reconstruction algorithm and implementation. MW performed the statistical analysis. SL conceived of the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.

Liting Wang, Xueli Li contributed equally to this work.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Wang, L., Li, X., Wu, M. et al. A sparse-projection computed tomography reconstruction method for in vivo application of in-line phase-contrast imaging. BioMed Eng OnLine 12, 75 (2013). https://doi.org/10.1186/1475-925X-12-75

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1475-925X-12-75

Keywords