Home | About | Journals | Submit | Contact Us | Français |

**|**HHS Author Manuscripts**|**PMC2858336

Formats

Article sections

- Abstract
- I. Background
- II. Theoretical Analysis
- III. Numerical Analysis
- IV. Discussions and Conclusion
- References

Authors

Related links

Phys Med Biol. Author manuscript; available in PMC 2010 April 22.

Published in final edited form as:

Published online 2009 August 28. doi: 10.1088/0031-9155/54/18/N04

PMCID: PMC2858336

NIHMSID: NIHMS192195

Email: gro.eeei@uy-gnoygneh, Email: nc.ude.ukp@gnaysj, Email: gro.eeei@gnaij-gnim, Email: gro.eeei@gnaw-eg

The publisher's final edited version of this article is available at Phys Med Biol

See other articles in PMC that cite the published article.

Recently, in the compressed sensing framework we proved that an interior ROI can be exactly reconstructed via the total variation minimization if the ROI is piecewise constant. In the proofs, we implicitly utilized the property that if an artifact image assumes a constant value within the ROI then this constant must be zero. Here we prove this property in the space of square integrable functions.

It is well known that an interior region-of-interest (ROI) cannot be uniquely reconstructed from projection data only associated with lines through the ROI. Interestingly, our group and others independently established that this interior problem can be exactly and stably solved if a sub-region in the ROI is precisely known (Ye *et al.*, 2007; Yu *et al.*, 2008; Courdurier *et al.*, 2008; Kudo *et al.*, 2008). Recently, in the compressed sensing framework we proved that an interior ROI can be exactly reconstructed via the total variation minimization if the imaging object is piecewise constant (Yu and Wang, 2009; Han *et al.*, 2009). In our proofs of Theorems 2.2 and 2.3, after we made the conclusion that the minimum was achieved at *λ* = 0 in (Yu and Wang, 2009) we implicitly utilized the property that if an artifact image assumes a constant value within the ROI then this constant must be zero. Here we provide a proof for this property so that our compressed sensing based interior tomography theory is expressed rigorously.

Let us consider a 2D smooth image *f*(*ρ,θ*), *ρ*[0,1], *θ*[0,2*π*) on a compact support within a unit disk Ω . The Radon transform of *f*(*ρ,θ*) can be written as *R*(*s,*), *s* [–1,1], [0,*π*) . Suppose that we are only interested in reconstructing its interior part defined as *ρ*<*a* (0<*a*<1) from the corresponding local Radon transform data *R*(*s,*),|*s*|<*a* , which is also referred to as purely local projections. Based on the classic analysis, given these local data there is in general no unique solution (Natterer, 2001). Thus, any reconstructed image from such a local dataset can be viewed as a superposition of the exact reconstruction from a complete dataset *R*(*s,*),*s*[–1,1],[0,*π*) and an artificial reconstruction, which is an artifact image corresponding to another projection dataset (*s,*), *a* <|*s*|≤1,[0,*π*) . Although (*s,*)=0 for |*s*|<*a*, it can still produce a non-zero 2D local image *g*(*ρ,θ*), *ρ*[0,*a*), *θ*[0,2*π*) , inside the ROI. That is why we have the non-uniqueness of this interior problem (Natterer, 2001). It is well known that *g*(*ρ,θ*) is smooth and bounded inside the ROI if (*s,*) is continuous and bounded. Previously, we proved that both ${\scriptstyle \frac{\partial g(\rho ,\theta )}{\partial \rho}}$ and ${\scriptstyle \frac{\partial g(\rho ,\theta )}{\rho \partial \theta}}$ are smooth and bounded inside the ROI (Yu and Wang, 2009). In the this note we will prove that *g*(*ρ,θ*)*C* for *ρ*[0,*a*) if and only if *C*=0.

The organization of this report is as follows. In the next section, we first prove that the constant *C* must be zero in the symmetric case using Taylor's expansion technique and Weierstrass's approximation theorem, and then extend it to the general non-symmetric case. In the third section, we numerically verify that a non-zero constant cannot satisfy the so-called Picard condition in the framework of singular value decomposition. In the last section, we discuss the relationship between this report and our previous work as well as future directions

Assuming that a circularly symmetric artifact image *g*(*ρ,θ*)=*g*(*ρ*) is reconstructed from a projection dataset (*s,*) = \(*s*) (*a* <|*s*|≤1,[0,*π*)) and (*s*)0 for *s* [–*a,a*] (Figure 1) and *g*(*ρ*)*C* (*ρ*[0,*a*) ), where *C* is a non-zero constant. Then, by the definition of the line integral we must have the following relationship for *g*(*ρ*) (*ρ*[*a*,1]):

$$\stackrel{~}{R}\left(s\right)={\int}_{a}^{1}g\left(\rho \right)d\sqrt{{\rho}^{2}-{s}^{2}}+C\sqrt{{a}^{2}-{s}^{2}}={\int}_{a}^{1}\frac{\rho g\left(\rho \right)d\phantom{\rule{thinmathspace}{0ex}}\rho}{\sqrt{{\rho}^{2}-{s}^{2}}}+C\sqrt{{a}^{2}-{s}^{2}}=0,s\in [0,a].$$

(1)

Our main theoretical results are as follows.

**Theorem 2.1.** Assuming that a circularly symmetric artifact image *g*(*ρ,θ*) = *g*(*ρ*) is reconstructed from a projection dataset (*s,*)=(*s*) (|*s*|≤1,[0,*π*)) and (*s*) 0 for *s* [–*a,a*]. If *g*(*ρ*) is a square integrable function on [0,1] with *g*(*ρ*)*C* (*ρ*[0,*a*) ), then *C*=0.

**Proof**: Eq. (1) can be changed to

$$Ca\sqrt{1-{(s\u2215a)}^{2}}+{\int}_{a}^{1}\frac{g\left(\rho \right)d\phantom{\rule{thinmathspace}{0ex}}\rho}{\sqrt{1-{(s\u2215\rho )}^{2}}}=0,\phantom{\rule{thickmathspace}{0ex}}0\le s\le a.$$

(2)

Denote (2*k*–1)!!=1×3×5---× (2*k*–1) and *k*!=1×2×3---×*k*, we have the following Taylor's expansions

$$\frac{1}{\sqrt{1-{(s\u2215\rho )}^{2}}}=1+\sum _{k=1}^{\infty}\frac{(2k-1)!!}{{2}^{k}k!}{\left(\frac{s}{\rho}\right)}^{2k},$$

(3)

$$\sqrt{1-{(s\u2215a)}^{2}}=1-\frac{1}{2k-1}\sum _{k=1}^{\infty}\frac{(2k-1)!!}{{2}^{k}k!}{\left(\frac{s}{a}\right)}^{2k}.$$

(4)

Inserting Eqs. (3) and (4) into Eq. (2), we have

$${\int}_{a}^{1}g\left(\rho \right){(a\u2215\rho )}^{2k}d\phantom{\rule{thinmathspace}{0ex}}\rho =\frac{Ca}{2k-1},k\ge 0.$$

(5)

Making the variable transformation *t*=*a*/*ρ*, Eq. (5) becomes

$${\int}_{a}^{1}g(a\u2215t){t}^{2k-2}dt=\frac{C}{2k-1},k\ge 0.$$

(6)

Let

$$\stackrel{~}{g}\left(t\right)=\{\begin{array}{cc}g(a\u2215t)-C\hfill & t\in [a,1]\hfill \\ \phantom{\rule{1em}{0ex}}\phantom{\rule{thickmathspace}{0ex}}-C\hfill & t\in [0,a)\hfill \end{array}\phantom{\}},$$

(7)

Eq. (6) implies

$${\int}_{0}^{1}\stackrel{~}{g}\left(t\right){t}^{2k}dt=0,k\ge 0.$$

(8)

Considering the even function *ĝ*(*t*)=((|*t*|) on the interval [–1,1], we have

$${\int}_{-1}^{1}\widehat{g}\left(t\right){t}^{k}dt=0,k\ge 0.$$

(9)

By Weierstrass's approximation theorem and the fact that the continuous function space is dense in the square integrable function space on [–1,1], there exist a polynomial series ${\left\{{p}_{n}\left(t\right)\right\}}_{n=1}^{\infty}$ on interval [–1,1] such that

$${\int}_{-1}^{1}{\mid \widehat{g}\left(t\right)-{P}_{n}\left(t\right)\mid}^{2}dt\le \frac{1}{n}.$$

(10)

This leads to

$${\int}_{-1}^{1}{\widehat{g}}^{2}\left(t\right)dt=\underset{n\to \infty}{\text{lim}}{\int}_{-1}^{1}\widehat{g}\left(t\right){p}_{n}\left(t\right)dt=0.$$

(11)

Hence, *ĝ*(*t*)=0 almost everywhere on [–1,1], which implies *C*= 0.

The above theorem can be extended to the general non-symmetric case.

**Theorem 2.2.** Assuming that a circularly non-symmetric artifact image *g*(*ρ,θ*) is reconstructed from a projection dataset (*s,*) (|*s*|≤1,[0,*π*) ) and (*s,*)0 for *s*[–*a,a*]. If *g*(*ρ,θ*) is a square integrable function on *ρ*[0,1] and *θ*[0,2*π*] with *g*(*ρ,θ*)*C* (*ρ*[0,*a*) ), then *C*=0.

**Proof:** Let us denote *g*(*ρ,θ*) in the Cartesian coordinates as *g _{d}*(

$$\begin{array}{cc}\hfill {\stackrel{\u2012}{g}}_{d}(x,y)& =\frac{1}{2\pi}{\int}_{0}^{2\pi}{g}_{d}(\rho \phantom{\rule{thinmathspace}{0ex}}\text{cos}(\theta +\varphi ),\rho \phantom{\rule{thinmathspace}{0ex}}\text{sin}(\theta +\varphi ))d\varphi \hfill \\ \hfill & =\frac{1}{2\pi}{\int}_{0}^{2\pi}{g}_{d}(x\phantom{\rule{thinmathspace}{0ex}}\text{cos}\phantom{\rule{thinmathspace}{0ex}}\varphi -y\phantom{\rule{thinmathspace}{0ex}}\text{sin}\phantom{\rule{thinmathspace}{0ex}}\varphi ,x\phantom{\rule{thinmathspace}{0ex}}\text{sin}\phantom{\rule{thinmathspace}{0ex}}\varphi +y\phantom{\rule{thinmathspace}{0ex}}\text{cos}\phantom{\rule{thinmathspace}{0ex}}\varphi )d\varphi .\hfill \end{array}$$

(12)

Clearly, (*ρ,θ*) is a circularly symmetric square integrable function and satisfies (*ρ,θ*)*C* (*ρ*[0,*a*) ). For *s* [–*a,a*] the Radon transform (*s,*) of (*ρ,θ*) can be expressed as

$$\begin{array}{cc}\hfill \stackrel{\u2012}{R}(s,\phi )& ={\int}_{-\sqrt{1-{s}^{2}}}^{\sqrt{1-{s}^{2}}}{\stackrel{\u2012}{g}}_{d}(s(\text{cos}\phantom{\rule{thinmathspace}{0ex}}\phi ,\text{sin}\phantom{\rule{thinmathspace}{0ex}}\phi )+t(-\text{sin}\phantom{\rule{thinmathspace}{0ex}}\phi ,\text{cos}\phantom{\rule{thinmathspace}{0ex}}\phi ))dt\hfill \\ \hfill & ={\int}_{-\sqrt{1-{s}^{2}}}^{\sqrt{1-{s}^{2}}}{\stackrel{\u2012}{g}}_{d}(s\phantom{\rule{thinmathspace}{0ex}}\text{cos}\phantom{\rule{thinmathspace}{0ex}}\phi -t\phantom{\rule{thinmathspace}{0ex}}\text{sin}\phantom{\rule{thinmathspace}{0ex}}\phi ,s\phantom{\rule{thinmathspace}{0ex}}\text{sin}\phantom{\rule{thinmathspace}{0ex}}\phi +t\phantom{\rule{thinmathspace}{0ex}}\text{cos}\phantom{\rule{thinmathspace}{0ex}}\phi )dt\hfill \\ \hfill & ={\int}_{-\sqrt{1-{s}^{2}}}^{\sqrt{1-{s}^{2}}}\left(\frac{1}{2\pi}{\int}_{0}^{2\pi}{g}_{d}(s\phantom{\rule{thinmathspace}{0ex}}\text{cos}(\phi +\varphi )-t\phantom{\rule{thinmathspace}{0ex}}\text{sin}(\phi +\varphi ),s\phantom{\rule{thinmathspace}{0ex}}\text{sin}(\phi +\varphi )+t\phantom{\rule{thinmathspace}{0ex}}\text{cos}(\phi +\varphi ))d\varphi \right)dt\hfill \\ \hfill & =\frac{1}{2\pi}{\int}_{0}^{2\pi}\left({\int}_{-\sqrt{1-{s}^{2}}}^{\sqrt{1-{s}^{2}}}{g}_{d}(s\phantom{\rule{thinmathspace}{0ex}}\text{cos}(\phi +\varphi )-t\phantom{\rule{thinmathspace}{0ex}}\text{sin}(\phi +\varphi ),s\phantom{\rule{thinmathspace}{0ex}}\text{sin}(\phi +\varphi )+t\phantom{\rule{thinmathspace}{0ex}}\text{cos}(\phi +\varphi ))dt\right)d\varphi \hfill \\ \hfill & =\frac{1}{2\pi}{\int}_{0}^{2\pi}\left({\int}_{-\sqrt{1-{s}^{2}}}^{\sqrt{1-{s}^{2}}}{g}_{d}(s(\text{cos}(\phi +\varphi ),\text{sin}(\phi +\varphi ))+t(-\text{sin}(\phi +\varphi ),\text{cos}(\phi +\varphi )))dt\right)d\varphi \hfill \\ \hfill & =\frac{1}{2\pi}{\int}_{0}^{2\pi}\stackrel{~}{R}(s,\phi +\varphi )d\varphi =0,\hfill \end{array}$$

(13)

where Eq.(12) has been used. By Theorem 2.1., we have *C*=0.

**Remark**: For real-world applications, all the reconstructed CT artifact images are bounded function and hence square integrable. Therefore, our assumption of square integrable functions can cover all the practical cases.

Assuming that an integral kernel *K*(*s,ρ*) defined on *I _{s}*×

$$K(s,\rho )=\sum _{i=1}^{\infty}{\sigma}_{i}{u}_{i}\left(s\right){v}_{i}\left(\rho \right),$$

(14)

where σ_{1} ≥σ_{2} ≥---> 0 , and *u _{i}*(

$$\left({u}_{i},{u}_{j}\right)={\int}_{{I}_{s}}{u}_{i}\left(t\right){u}_{j}\left(t\right)d=\{\begin{array}{cc}1\hfill & \text{if}\phantom{\rule{thickmathspace}{0ex}}i=j\hfill \\ 0\hfill & \text{if}\phantom{\rule{thickmathspace}{0ex}}i\ne j\hfill \end{array}\phantom{\}},$$

(15)

and

$$\left({v}_{i},{v}_{j}\right)={\int}_{{I}_{\rho}}{v}_{i}\left(t\right){v}_{j}\left(t\right)dt=\{\begin{array}{cc}1\hfill & \text{if}\phantom{\rule{thickmathspace}{0ex}}i=j\hfill \\ 0\hfill & \text{if}\phantom{\rule{thickmathspace}{0ex}}i\ne j\hfill \end{array}\phantom{\}}.$$

(16)

For any continuous and square integral function *q*(*s*), *s**I _{s}*, the Fredholm integral equation of the first kind

$${\int}_{{I}_{\rho}}K(s,\rho )g\left(\rho \right)d\phantom{\rule{thinmathspace}{0ex}}\rho =q\left(s\right)$$

(17)

can be expressed as

$$\sum _{i=1}^{\infty}{\sigma}_{i}({v}_{i},g){u}_{i}\left(s\right)=\sum _{i=1}^{\infty}({u}_{i},q){u}_{i}\left(s\right),s\in {I}_{s}.$$

(18)

If there exists solutions for Eq. (17), one of them would be in the following form:

$$g\left(\rho \right)=\sum _{i=1}^{\infty}\frac{({u}_{i},q)}{{\sigma}_{i}}{v}_{i}\left(\rho \right),\rho \in {I}_{\rho}.$$

(19)

Eqs. (18) and (19) hold for uniform convergence (see Chapter 2.4 and Theorem 8.3.2 in (Smithies, 1958)). Although *K*(*s,ρ*) and *q*(*s*) are continuous and square integrable, *g*(*ρ*) expressed by Eq. (19) is not necessarily continuous or square integrable. In order that there exists a square integrable solution *g*(*ρ*) to Eq. (17), the right side of Eq. (19) must satisfy the Picard condition (see Chapter 1.2.3 in (Hansen, 1998)):

$$\sum _{i=1}^{\infty}{\left(\frac{({u}_{i},q)}{{\sigma}_{i}}\right)}^{2}<\infty .$$

(20)

This means that the (*u _{i}, q*) must decay faster than the singular value σ

Let $K(s,\rho )=\frac{\rho}{\sqrt{{\rho}^{2}-{s}^{2}}}$, *I _{ρ}* = [

In our case, let us discrete *s* as ${s}_{n}={\scriptstyle \frac{a(n-0.5)}{N}}$, *n* = 1,2,---,*N* and *ρ* as ${\rho}_{m}=a+{\scriptstyle \frac{(m-0.5)(1-a)}{M}}$, *m* = 1,2,---,*M*. Accordingly, *K*(*s,ρ*), *g*(*ρ*) and *q*(*s*) are discretized as *K*(*s _{n},ρ_{m}*),

$$\mathbf{KG}=\mathbf{Q},$$

(21)

where

$$\mathbf{K}=\frac{1-a}{M}\left(\begin{array}{ccccc}\hfill K({s}_{1},{\rho}_{1})\hfill & \hfill \cdots \hfill & \hfill K({s}_{1},{\rho}_{m})\hfill & \hfill \cdots \hfill & \hfill K({s}_{1},{\rho}_{M})\hfill \\ \hfill \vdots \hfill & \hfill \ddots \hfill & \hfill \vdots \hfill & \hfill \ddots \hfill & \hfill \vdots \hfill \\ \hfill K({s}_{n},{\rho}_{1})\hfill & \hfill \cdots \hfill & \hfill K({s}_{n},{\rho}_{m})\hfill & \hfill \cdots \hfill & \hfill K({s}_{n},{\rho}_{M})\hfill \\ \hfill \vdots \hfill & \hfill \ddots \hfill & \hfill \vdots \hfill & \hfill \ddots \hfill & \hfill \vdots \hfill \\ \hfill K({s}_{N},{\rho}_{1})\hfill & \hfill \cdots \hfill & \hfill K({s}_{N},{\rho}_{m})\hfill & \hfill \cdots \hfill & \hfill K({s}_{N},{\rho}_{M})\hfill \end{array}\right),$$

(22)

$$\mathbf{G}={(g\left({\rho}_{1}\right)\phantom{\rule{1em}{0ex}}\cdots \phantom{\rule{1em}{0ex}}g\left({\rho}_{m}\right)\phantom{\rule{1em}{0ex}}\cdots \phantom{\rule{1em}{0ex}}g\left({\rho}_{M}\right))}^{T},$$

(23)

and

$$\mathbf{Q}={(q\left({s}_{1}\right)\phantom{\rule{1em}{0ex}}\cdots \phantom{\rule{1em}{0ex}}q\left({s}_{n}\right)\phantom{\rule{1em}{0ex}}\cdots \phantom{\rule{1em}{0ex}}q\left({s}_{N}\right))}^{T}.$$

(24)

Let *M*=*N*, the matrix **K** has the following SVD form (Hansen, 1988):

$$\mathbf{K}=\mathbf{U}\Lambda {\mathbf{V}}^{T}=\sum _{m=1}^{M}{\sigma}_{m}{\mathbf{u}}_{m}{\mathbf{v}}_{m}^{T},$$

(25)

where **U**=(**u**_{1} --- **u*** _{m}* ---

$$\mathbf{G}=\sum _{i=m}^{M}\frac{{\mathbf{u}}_{m}^{T}\mathbf{Q}}{{\sigma}_{m}}{\mathbf{v}}_{m}.$$

(26)

According to the Picard condition, we can compute $\mid {\scriptstyle \frac{{\mathbf{u}}_{m}^{T}\mathbf{Q}}{{\sigma}_{m}}}\mid $ to test the existence of a solution to Eq. (1).

Without loss of generality, we set *a*=0.5 and $q\left(s\right)=\sqrt{{a}^{2}-{s}^{2}}$. For different combinations of sampling point numbers *M* and *N*, we computed and plotted $\mathrm{lg}\mid {\scriptstyle \frac{{\mathbf{u}}_{m}^{T}\mathbf{Q}}{{\sigma}_{m}}}\mid $ in Figure 2, which clearly show that the singular value σ* _{m}* decays faster than ${\mathbf{u}}_{m}^{T}\mathbf{Q}$, suggesting that there is no square integrable solution for

Although we assumed that an object to be reconstructed was piecewise constant both inside and outside an ROI, the piecewise constancy is only needed inside the ROI in the proofs in (Yu and Wang, 2009; Han *et al.*, 2009). In other words, the piecewise constancy outside the ROI is not necessary. Note that a generalized total variation minimization theorem was proved for compressed sensing based interior tomography in (Han et al., 2009), without involving the Dirac delta function. Because most objects in CT applications can be approximately modeled as piecewise constant, the piecewise constancy assumption is quite reasonable (Wang *et al.*, 2004), and our compressed sensing based interior tomography approach is practically useful. Most importantly, our finding suggests a new direction of interior tomography, and may lead to powerful methods for better interior reconstruction. Currently, we are working to relax the piecewise constancy assumption so that real-world images can be modeled with higher accuracy for new interior reconstruction results.

As discussed in (Yu and Wang, 2009), the key idea of compressed sensing based interior tomography is to define an appropriate sparsifying transform and an associated objective function. Then, the minimization of the objective function will lead to the true ROI image without any ambiguity. For that purpose, the most important step is to prove that the artifact image inside the ROI must be zero when its total variation is zero. Actually, it might appear in (Yu and Wang, 2009) that the artifact image inside the ROI can be any constant when the total variation is zero. Hence, it is the major contribution of this report that has proved the constant being necessarily zero, providing the uniqueness of our interior reconstruction. Using this methodology, we are working to extend compressed sensing based interior reconstruction under different conditions and for other modalities.

In conclusion, we have studied the characteristic of the artifact image inside an interior ROI subject to the local projection data constraint. Our theoretical and numerical analysis has demonstrated that the artifact image inside the ROI can be expressed as a constant if and only if the constant is zero. This analysis has consolidated our proofs of Theorems 2.2 and 2.3 in (Yu and Wang, 2009).

This work is partially supported by NIH/NIBIB grants (EB002667, EB004287, EB007288), a grant from Toshiba Medical Research Institute USA, Inc., NSFC (60325101, 60532080, 60628102), Key Laboratory of Machine Perception (Ministry of Education) of Peking University, and Microsoft Research of Asia. The first author thanks Dr. Per Christian Hansen for offering a hardcopy of his paper (Hansen, 1988).

- Courdurier M, Noo F, Defrise M, Kudo H. Solving the interior problem of computed tomography using a priori knowledge. Inverse Problems. 2008;24 Article ID 065001 , 27 pages. [PMC free article] [PubMed]
- Han W, Yu H, Wang G. A total variation minimization theorem for compressed sensing based tomography, under review. 2009. [PMC free article] [PubMed]
- Hansen PC. Computation of the singular value expansion. Computing. 1988;40:185–99.
- Hansen PC. Rank-Deficient and Discrete Ill-Posed Problems: Numerical Aspects of Linear Inversion. Society for Industrial and Applied Mathematics; Philadelphia: 1998.
- Kudo H, Courdurier M, Noo F, Defrise M. Tiny a priori knowledge solves the interior problem in computed tomography. Phys. Med. Biol. 2008;53:2207–31. [PubMed]
- Natterer F. The Mathematics of Computerized Tomography. Society for Industrial and Applied Mathematics; Philadelphia: 2001.
- Smithies F. Integral Equations. Cambridge University Press; 1958.
- Wang G, Li Y, Jiang M. Uniqueness theorems in bioluminescence tomography. Medical Physics. 2004;31:2289–99. [PubMed]
- Ye Y, Yu H, Wei Y, Wang G. A general local reconstruction approach based on a truncated Hilbert transform. International Journal of Biomedical Imaging. 2007;2007 Article ID: 63634, 8 pages. [PMC free article] [PubMed]
- Yu H, Wang G. Compressed sensing based Interior tomography. Phys Med Biol. 2009;54(9):2791–2805. [PMC free article] [PubMed]
- Yu H, Ye Y, Wang G. Local Reconstruction Using the Truncated Hilbert Transform via Singular Value Decomposition. Journal of X-Ray Science and Technology. 2008;16:243–51. [PMC free article] [PubMed]

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |