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

**|**HHS Author Manuscripts**|**PMC2958053

Formats

Article sections

- Abstract
- 1. Introduction
- 2. SVD analysis
- 3. Numerical results
- 4. Discussion and future work
- 5. Conclusions
- References and links

Authors

Related links

Opt Express. Author manuscript; available in PMC 2010 October 20.

Published in final edited form as:

PMCID: PMC2958053

NIHMSID: NIHMS240862

The publisher's final edited version of this article is available at Opt Express

See other articles in PMC that cite the published article.

Tomosynthesis is an emerging technique with potential to replace mammography, since it gives 3D information at a relatively small increase in dose and cost. We present an analytical singular-value decomposition of a tomosynthesis system, which provides the measurement component of any given object. The method is demonstrated on an example object. The measurement component can be used as a reconstruction of the object, and can also be utilized in future observer studies of tomosynthesis image quality.

The tomosynthesis method [1, 2] is currently considered to be of great potential, especially in detection of breast cancer [3]. The idea is to replace mammography, which gives two-dimensional (2D) projection, with tomosynthesis which gives three-dimensional (3D) information e.g. as slice images. The ultimate purpose is to improve the detection of breast cancer, at a relatively small increase in dose and cost. Prototype tomosynthesis systems have already been developed and pilot clinical trials performed [4–7]. These trials have so far not demonstrated distinct advantages in detection compared to mammography, but indicate that tomosynthesis images have better visualization of structures without masking by overlapping tissue. Improved detection has been clearly demonstrated in studies on mastectomy specimens [8].

Tomosynthesis can be said to be half-way between a single projection as used in mammography, and the full-angle scan performed in computed tomography (CT). The first produces a 2D representation and the other a 3D representation. Tomosynthesis, where a number of projections are taken at a limited angle, provides partial 3D information — the resolution in one direction must be limited. The hope is that tomosynthesis will provide the advantage of 3D information over conventional mammography, without the increased dose and cost of a full CT.

The 3D object, for example, the breast tissue, must be reconstructed from the 2D projection images. Several techniques are used [9,10], mainly back-projection methods, iterative methods, or statistical methods. The reconstruction task is more difficult than in tomography, since the information is incomplete and a true inverse transform cannot be found. Especially in analytical methods such as back-projection, this has a tendency to produce artifacts. We propose to use singular-value decomposition (SVD) [11–13] for the reconstruction, in a manner very similar to the technique we developed for through-focus imaging systems [14]. The result produced is the measurement component of the object, which could be regarded as its pseudoinverse [11] reconstruction and is as close to the inverse as we can get. The measurement component is the part of the object that is actually transmitted into the image. An advantage of the method is its mathematical stringency: what does not show up in the measurement component simply cannot be reconstructed, and the measurement component is unique for each set of measurements. A disadvantage is that the measurement component is different from a conventional reconstruction since it contains negative values. The measurement component could be used for further studies on image quality and detectability, using a channelized Hotelling observer [11, 15, 16].

At present, we are analyzing a parallel-path setup where the x-ray source moves in a plane parallel to the detector plane. Other geometries [2] could be analyzed but it would require some work to adjust the analysis. We consider the geometrical approximation, disregarding any effects of diffraction. The object is seen as purely absorbing, with no scattering or phase contribution. Scattering cannot be included, since it breaks the lateral shift-invariance. As we analyze the ideal case, we do not yet consider noise.

We consider a geometry where all source positions are in one transverse plane, and images captured in a second plane parallel to the first and placed behind the object. This geometry is illustrated in Fig. 1. In the figure the object is shown as limited in both *z*, *x* and *y* direction. We present this case of limited object and image support as the most general, but will mainly consider the simplified case of infinite transverse extent. The support in the *z* direction is per definition limited to at most the distance between the source and image planes. The variables used are **r** = (*x*, *y*) for transverse coordinates in object space and *z* for longitudinal coordinate in object space. The image plane is placed at the origin (*z* = 0), and has transverse coordinates **r*** _{d}* = (

Geometry of the considered tomographic system. The different source positions are *χ*_{m}, and the screen is placed along the *z* axis.

A relation between object and image can be found from the steady-state solution of the Boltzmann transport equation in absorbing media, assuming a point source [11]:

$${I}_{m}({\mathbf{r}}_{d})=\frac{Acos\theta}{{\mid {\mathbf{r}}_{1}-{\mathbf{r}}_{0}\mid}^{2}}exp\left[-{\int}_{0}^{\mid {\mathbf{r}}_{1}-{\mathbf{r}}_{0}\mid}\text{d}{l}^{\prime}\phantom{\rule{0.16667em}{0ex}}f\left({\mathbf{r}}_{1}-\frac{{\mathbf{r}}_{1}-{\mathbf{r}}_{0}}{\mid {\mathbf{r}}_{1}-{\mathbf{r}}_{0}\mid}{l}^{\prime}\right)\right],$$

where *I _{m}*(

$$log\left(\frac{{I}_{m}({\mathbf{r}}_{d}){\mid {\mathbf{r}}_{1}-{\mathbf{r}}_{0}\mid}^{2}}{Acos\theta}\right)=\frac{\mid {\mathbf{r}}_{1}-{\mathbf{r}}_{0}\mid}{s}{\int}_{0}^{-s}\text{d}z\phantom{\rule{0.16667em}{0ex}}f\left({\mathbf{r}}_{0}+({\mathbf{r}}_{0}-{\mathbf{r}}_{1})\frac{z}{s}\right).$$

Using cos *θ* = *s*/|**r**_{1} − **r**_{0}| and |**r**_{1} − **r**_{0}| = (|**r*** _{d}* −

$$log\left[\frac{1}{As}{({\mid {\mathbf{r}}_{d}-{\mathbf{r}}_{m}\mid}^{2}+{s}^{2})}^{3/2}{I}_{m}({\mathbf{r}}_{d})\right]={\int}_{0}^{-s}\text{d}z\phantom{\rule{0.16667em}{0ex}}f\left({\mathbf{r}}_{d}+({\mathbf{r}}_{d}-{\mathbf{r}}_{m})\frac{z}{s},z\right)$$

(1)

We introduce a relation between the irradiance *I _{m}*(

$${g}_{m}({\mathbf{r}}_{d})=\frac{s}{{({s}^{2}+{\mid {\mathbf{r}}_{d}-{\mathbf{r}}_{m}\mid}^{2})}^{1/2}}log\left[\frac{1}{As}{({s}^{2}+{\mid {\mathbf{r}}_{d}-{\mathbf{r}}_{m}\mid}^{2})}^{3/2}{I}_{m}({\mathbf{r}}_{d})\right].$$

(2)

This definition of *g _{m}*(

(3)

and its adjoint (or backprojection) operator

(4)

where *z*_{1} and *z*_{2} limit the object extension in the *z* direction, −*s* < *z*_{1} < *z*_{2} < 0, M is the number of source positions, *b _{f}* (

$${h}_{m}(\mathbf{r},{\mathbf{r}}_{d},z)=\delta \phantom{\rule{0.16667em}{0ex}}\left(\mathbf{r}-{\mathbf{r}}_{d}-({\mathbf{r}}_{d}-{\mathbf{r}}_{m})\frac{z}{s}\right).$$

(5)

The object and image support contained in *b _{f}* (

$$\stackrel{\sim}{\mathbf{r}}=(\stackrel{\sim}{x},\stackrel{\sim}{y})=\left(\frac{s}{s+z}x,\frac{s}{s+z}y\right),\stackrel{\sim}{z}=\frac{s}{s+z}z$$

(6)

with the inverse relations

$$\mathbf{r}=(x,y)=\left(\frac{s}{s+\stackrel{\sim}{z}}\stackrel{\sim}{x},\frac{s}{s-\stackrel{\sim}{z}}\stackrel{\sim}{y}\right),z=\frac{s}{s-\stackrel{\sim}{z}}\stackrel{\sim}{z}.$$

(7)

Intuitively, the transformation in Eq. (6) means the source plane is moved to minus infinity and the rays from one source point become parallel as illustrated in Fig. 2. The position **r*** _{m}* of the source determines the inclination of the rays. The object is distorted by this transformation, but the image remains the same. In these coordinates, the impulse-response function becomes

$${\stackrel{\sim}{h}}_{m}(\stackrel{\sim}{\mathbf{r}}-{\mathbf{r}}_{d},\stackrel{\sim}{z})={\left(\frac{s-\stackrel{\sim}{z}}{s}\right)}^{2}\delta \left(\stackrel{\sim}{\mathbf{r}}-{\mathbf{r}}_{d}+\frac{\stackrel{\sim}{z}}{s}{\mathbf{r}}_{m}\right),$$

(8)

where * _{m}*( −

$$\left|\frac{\partial (x,y,z)}{\partial (\stackrel{\sim}{x},\stackrel{\sim}{y},\stackrel{\sim}{z})}\right|=\frac{{s}^{4}}{{(s-\stackrel{\sim}{z})}^{4}}$$

the integral relations now become

(9)

and

(10)

where (, ) = *f* (**r**, *z*) is the distorted object function and * _{f}* (, ) =

$${\langle {\stackrel{\sim}{f}}_{1}(\stackrel{\sim}{\mathbf{r}},\stackrel{\sim}{z}),{\stackrel{\sim}{f}}_{2}(\stackrel{\sim}{\mathbf{r}},\stackrel{\sim}{z})\rangle}_{\text{obj}}=\int {\int}_{-\infty}^{\infty}{\text{d}}^{2}\stackrel{\sim}{r}{\int}_{{\stackrel{\sim}{z}}_{1}}^{{\stackrel{\sim}{z}}_{2}}\text{d}\stackrel{\sim}{z}\phantom{\rule{0.16667em}{0ex}}{\stackrel{\sim}{b}}_{f}(\stackrel{\sim}{\mathbf{r}},\stackrel{\sim}{z}){\stackrel{\sim}{f}}_{1}(\stackrel{\sim}{\mathbf{r}},\stackrel{\sim}{z}){\stackrel{\sim}{f}}_{2}^{\ast}(\stackrel{\sim}{\mathbf{r}},\stackrel{\sim}{z})\frac{{s}^{4}}{{(s-\stackrel{\sim}{z})}^{4}}$$

(11)

and

$${\langle {\mathbf{g}}_{1}({\mathbf{r}}_{d}),{\mathbf{g}}_{2}({\mathbf{r}}_{d})\rangle}_{\text{im}}=\sum _{m=1}^{M}\int {\int}_{-\infty}^{\infty}{\text{d}}^{2}{r}_{d}{b}_{g}({\mathbf{r}}_{d}){g}_{1m}({\mathbf{r}}_{d}){g}_{2m}^{\ast}({\mathbf{r}}_{d}).$$

The weighted inner product of Eq. (11) is required in order to fulfill the adjoint relation

Combining Eqs. (9) and (10) finally yields the image-space Hermitian operator as

(12)

where

$${k}_{m{m}^{\prime}}({\mathbf{r}}_{d},{\mathbf{r}}_{d}^{\prime})={\int}_{{\stackrel{\sim}{z}}_{1}}^{{\stackrel{\sim}{z}}_{2}}\text{d}\stackrel{\sim}{z}\int {\int}_{-\infty}^{\infty}{\text{d}}^{2}\stackrel{\sim}{r}{\stackrel{\sim}{b}}_{f}(\stackrel{\sim}{\mathbf{r}},\stackrel{\sim}{z})\frac{{s}^{4}}{{(s-\stackrel{\sim}{z})}^{4}}{\stackrel{\sim}{h}}_{m}(\stackrel{\sim}{\mathbf{r}}-{\mathbf{r}}_{d},\stackrel{\sim}{z}){\stackrel{\sim}{h}}_{{m}^{\prime}}^{\ast}(\stackrel{\sim}{\mathbf{r}}-{\mathbf{r}}_{d}^{\prime},\stackrel{\sim}{z}).$$

(13)

Similarly, the object-space Hermitian operator will be

where

$$p(\stackrel{\sim}{\mathbf{r}},{\stackrel{\sim}{\mathbf{r}}}^{\prime},\stackrel{\sim}{z},{\stackrel{\sim}{z}}^{\prime})=\frac{{s}^{4}}{{(s-\stackrel{\sim}{z})}^{4}}\sum _{m=1}^{M}\int {\int}_{-\infty}^{\infty}{\text{d}}^{2}{r}_{d}{b}_{g}({\mathbf{r}}_{d}){\stackrel{\sim}{h}}_{m}^{\ast}(\stackrel{\sim}{\mathbf{r}}-{\mathbf{r}}_{d},\stackrel{\sim}{z}){\stackrel{\sim}{h}}_{m}({\stackrel{\sim}{\mathbf{r}}}^{\prime}-{\mathbf{r}}_{d},{\stackrel{\sim}{z}}^{\prime}).$$

If both object and image are considered infinite in the transverse direction, and the object confined in the longitudinal direction by *z*_{1} < *z* < *z*_{2}, both support functions *b _{f}* (

$${k}_{m{m}^{\prime}}({\mathbf{r}}_{d}-{\mathbf{r}}_{d}^{\prime})={\int}_{{\stackrel{\sim}{z}}_{1}}^{{\stackrel{\sim}{z}}_{2}}\text{d}\stackrel{\sim}{z}\delta \left({\mathbf{r}}_{d}-{\mathbf{r}}_{d}^{\prime}-\frac{\stackrel{\sim}{z}}{s}({\mathbf{r}}_{m}-{\mathbf{r}}_{m}^{\prime})\right).$$

(14)

This integral could easily be evaluated, but the subsequent analysis gets easier if it is left as it is. Assuming singular functions of the form

$${\mathbf{v}}_{\rho ,j}({\mathbf{r}}_{d})={\mathbf{V}}_{j}(\rho )exp(2\pi \text{i}\rho \xb7{\mathbf{r}}_{d})$$

(15)

yields, after insertion into Eq. (12), the eigenequation [14]

(16)

where (*ρ*) is an *M* × *M* matrix with elements

$${K}_{m{m}^{\prime}}=\int {\int}_{-\infty}^{\infty}{\text{d}}^{2}{r}_{d}{k}_{m{m}^{\prime}}({\mathbf{r}}_{d})exp(-2\pi \text{i}\rho \xb7{\mathbf{r}}_{d}).$$

(17)

Inserting Eq. (14) into Eq. (17) and performing the integral in **r*** _{d}* yields

$${K}_{m{m}^{\prime}}={\int}_{{\stackrel{\sim}{z}}_{1}}^{{\stackrel{\sim}{z}}_{2}}\text{d}\stackrel{\sim}{z}exp\phantom{\rule{0.16667em}{0ex}}\left(-2\pi \text{i}\rho \xb7({\mathbf{r}}_{m}-{\mathbf{r}}_{{m}^{\prime}})\frac{\stackrel{\sim}{z}}{s}\right),$$

(18)

which is integrated to yield

$${K}_{m{m}^{\prime}}={\stackrel{\sim}{z}}_{2}-{\stackrel{\sim}{z}}_{1}$$

for *m* = *m*′ or *ρ* = (0, 0), and

$${K}_{m{m}^{\prime}}=\frac{\text{i}s}{2\pi \rho \xb7({\mathbf{r}}_{m}-{\mathbf{r}}_{{m}^{\prime}})}\left\{exp\left[-2\pi \text{i}\rho \xb7({\mathbf{r}}_{m}-{\mathbf{r}}_{{m}^{\prime}})\frac{{\stackrel{\sim}{z}}_{2}}{s}\right]-exp\left[-2\pi i\rho \xb7({\mathbf{r}}_{m}-{\mathbf{r}}_{{m}^{\prime}})\frac{{\stackrel{\sim}{z}}_{1}}{s}\right]\right\}$$

for *m*≠ *m*′. Thus the elements of have been found, and Eq. (16) can be solved numerically to yield the vectors **V*** _{j}*(

The solution can also be obtained in object space. The results will be identical to those obtained from the image-space analysis, but we still outline the analysis required. The kernel of the Hermitian operator is

$$p(\stackrel{\sim}{\mathbf{r}}-{\stackrel{\sim}{\mathbf{r}}}^{\prime},\stackrel{\sim}{z},{\stackrel{\sim}{z}}^{\prime})={\left(\frac{s-\stackrel{\sim}{z}}{s-{\stackrel{\sim}{z}}^{\prime}}\right)}^{2}\sum _{m=1}^{M}\delta \left(\stackrel{\sim}{\mathbf{r}}-{\stackrel{\sim}{\mathbf{r}}}^{\prime}+{\mathbf{r}}_{m}\frac{\stackrel{\sim}{z}-{\stackrel{\sim}{z}}^{\prime}}{s}\right),$$

(19)

leading to the object-space eigenequation

$${\int}_{{\stackrel{\sim}{z}}_{1}}^{{\stackrel{\sim}{z}}_{2}}\text{d}{\stackrel{\sim}{z}}^{\prime}{\left(\frac{s-\stackrel{\sim}{z}}{s-{\stackrel{\sim}{z}}^{\prime}}\right)}^{2}{U}_{j}(\rho ,{\stackrel{\sim}{z}}^{\prime})\sum _{m=1}^{M}exp\left(2\pi \text{i}\rho \xb7{\mathbf{r}}_{m}(\stackrel{\sim}{z}-{\stackrel{\sim}{z}}^{\prime})\frac{1}{s}\right)={\mu}_{\rho ,j}{U}_{j}(\rho ,\stackrel{\sim}{z})$$

where it has been assumed the singular functions can be written as

$${u}_{j,\rho}(\stackrel{\sim}{\mathbf{r}},\stackrel{\sim}{z})={U}_{j}(\rho ,\stackrel{\sim}{z})exp(2\pi \text{i}\rho \xb7\stackrel{\sim}{\mathbf{r}}).$$

Introducing *Û* (*ρ*, )(*s* − )^{2} = *U*(*ρ*, ) we find the integral equation

$${\widehat{U}}_{j}(\rho ,\stackrel{\sim}{z})={\int}_{{\stackrel{\sim}{z}}_{1}}^{{\stackrel{\sim}{z}}_{2}}d{\stackrel{\sim}{z}}^{\prime}\widehat{U}(\rho ,{\stackrel{\sim}{z}}^{\prime})\sum _{m=1}^{M}exp\left(2\pi \text{i}\frac{\stackrel{\sim}{z}}{s}\rho \xb7{\mathbf{r}}_{m}\right)exp\left(-2\pi \text{i}\frac{{\stackrel{\sim}{z}}^{\prime}}{s}\rho \xb7{\mathbf{r}}_{m}\right).$$

This integral equation has a degenerate kernel, and such equations can be solved using standard methods [19].

Equation (19) is the kernel of the object-space Hermitian operator, but can also be seen the impulse-response function of the system. A delta function at (′, ′) in the object will cause a distribution *p*( − ′, , ′) in the reconstruction. The impulse response function is a sum of *M* lines whose intensity varies with , and that all meet at = ′, = ′.

In section 2, the singular functions have been found for the simplified case of infinite object and image. These results can be illustrated numerically. Under certain circumstances, the SVD technique is not very time-consuming and can be run on a normal desktop computer. The calculation times are on the order of minutes or tens of minutes. The most time-consuming part is actually the simulation of the object (see section 3.2). The technique for numerical SVD calculations is very similar to Burvall et al. [14].

The geometry is the same as in Fig. 1. We have assumed a distance of 30 cm from source to image plane, and that the object is located between 10 and 1 cm from the image plane, i.e. *s* = 0.3 m, *z*_{1} = −0.1 m, and *z*_{2} = −0.01 m. Furthermore, it is assumed that all the sources are placed along a line or arc rather than in a 2D array. While it would be interesting to investigate other arrangements, this option simplifies the analysis significantly and is often employed in tomosynthesis. Assuming the source distribution is along the *x* axis, we have **r*** _{m}* = (

The assumption that **r*** _{m}* = (

Once the image-space singular functions are known, the object-space singular functions *U _{ρ}*

(20)

This equation is also known as part of the shifted-eigenvalue equations [20]. Insertion of Eqs. (9), (8) and (15) into Eq. (20) yields

$${u}_{\rho ,j}(\stackrel{\sim}{\mathbf{r}},\stackrel{\sim}{z})=exp(2\pi \text{i}\rho \xb7\stackrel{\sim}{\mathbf{r}}){U}_{j}(\rho ,\stackrel{\sim}{z})$$

where

$${U}_{j}(\rho ,\stackrel{\sim}{z})=\frac{1}{\sqrt{{\mu}_{\rho ,j}}}{\left(\frac{s-\stackrel{\sim}{z}}{s}\right)}^{2}\sum _{m=1}^{M}{[{\mathbf{V}}_{j}(\rho )]}_{m}exp\left(2\pi \text{i}\frac{\stackrel{\sim}{z}}{s}\rho \xb7{\mathbf{r}}_{m}\right).$$

So *U _{j}*(

Absolute values of the singular functions (a) *U*_{1}(*ρ*_{x}, ), (b) *U*_{2}(*ρ*_{x}, ), (c) *U*_{6}(*ρ*_{x}, ), and (d) *U*_{11}(*ρ*_{x}, ).

Unlike the 3D imaging case [14], where the singular functions are real, these functions are complex. In Fig. 4 their absolute value is shown. Interpreting the results intuitively is fairly difficult, as there are no obvious image planes and the coordinates are a combination of spatial () and frequency (*ρ*). It is easier to interpret the expansion of an example object into these eigenfunctions, i.e., its measurement component. This requires a relevant sample object.

We illustrate the method using an example of at least some relevance to the future task, namely tomosynthesis of breast tissue. The object consists of two parts: a designer nodule and a clustered lumpy background (CLB). This approach is chosen since the expected advantage of tomosynthesis is the ability to detect low-contrast masses, such as the nodule, that are masked by overlapping tissue, such as the clustered lumpy background. The designer nodule has a projection onto the image plane given by given by [21]

$${g}_{m}({\mathbf{r}}_{d})=A\text{rect}\left(\frac{2\mid {\mathbf{r}}_{d}\mid}{R}\right)\phantom{\rule{0.16667em}{0ex}}{\left(1-\frac{{\mid {\mathbf{r}}_{d}\mid}^{2}}{{R}^{2}}\right)}^{v}$$

(21)

where *A* is the amplitude of the nodule, *R* its radius, and *v* a real and non-negative shape factor that determines the sharpness of the edges. We have chosen *R* = 3mm, *A* = 0.03, and *v* = 1.5. This value of *v* is reported to give the best fit to average tumor profiles [21]. The 3D object that gives this projection is [22]

$$f(\mathbf{r})=A\frac{3}{4R}\text{rect}\left(\frac{2\mid \mathbf{r}\mid}{R}\right)\phantom{\rule{0.16667em}{0ex}}\left(1-\frac{{\mid \mathbf{r}\mid}^{2}}{{R}^{2}}\right),$$

(22)

as can be shown by finding its projection for a far-away source, i.e., by integrating Eq. (22) with respect to z.

The CLB simulates the 2D projection of the background by distributing lumps according to a statistical model [23]. The projection is given by

$$g({\mathbf{r}}_{d})=\sum _{k=1}^{K}\sum _{n=1}^{{N}_{k}}b\left(\frac{1}{{a}_{kn}}({\mathbf{r}}_{d}-{\mathbf{r}}_{k}-{\mathbf{r}}_{kn}),{\mathbf{R}}_{{\theta}_{k}n}\right)$$

where *K* is the numbers of clusters, each centered on **r*** _{k}*, and

followed the classic CLB as introduced by Bochud [23]. The number *K* was a Poisson random variable with expectation value 150 and *N _{k}* a Poisson random variable with expectation value 20. The scale

An example object (, ) is shown in Fig. 5, in the distorted coordinates and . It is shown in two slices: one through the middle of the designer nodule in the *ỹ* plane, and one through the middle of the designer nodule in the plane. The object is not a proper 3D breast model, since it is not connected enough for a breast CT slice [24], but it will give reasonable results once projected and will serve as an illustration.

Once an object has been generated, and the singular functions of the system found, they can be combined to generate the measurement component. It is an expansion of the object into the object-space singular functions according to

$${\stackrel{\sim}{f}}_{\text{meas}}(\stackrel{\sim}{\mathbf{r}},\stackrel{\sim}{z})=\int {\int}_{-\infty}^{\infty}{\text{d}}^{2}\rho \sum _{j=1}^{M}{A}_{j}(\rho ){U}_{j}(\rho ,\stackrel{\sim}{z})exp2\pi \text{i}\rho \xb7\stackrel{\sim}{\mathbf{r}},$$

where the integral is numerically found as a Fast Fourier Transform. The expansion coefficients *A _{j}*(

$${A}_{j}(\rho )={\int}_{{\stackrel{\sim}{z}}_{1}}^{{\stackrel{\sim}{z}}_{2}}\text{d}\stackrel{\sim}{z}\stackrel{\sim}{F}(\rho ,\stackrel{\sim}{z}){U}_{j}^{\ast}(\rho ,\stackrel{\sim}{z})\frac{{s}^{4}}{{(s-\stackrel{\sim}{z})}^{4}}$$

(23)

where (*ρ*, ) is the 2D Fourier transform of the object (, ) with respect to . Two things should be noted, as they differ from the expression used in Ref. [14]. First, the complex conjugate in Eq. 23 was omitted in [14] since the singular functions were real, but must be included here. Second, the weighting function *w*(, ) = *s*^{4}/(*s* − )^{4} follows from the earlier change of variables. Slices through the measurement component are shown in Fig. 6. The slices are taken at the same positions as those in Fig. 5. In part (b), the low axial resolution is showing up as elongation of both nodule and background.

The zero-angle projection of the object is shown in Fig. 7(a), along with the same projection of its measurement component in part (b). The images are generated by direct numerical propagation of the quantities involved, using a sampling of 400 points in the transverse directions and 300 points in the longitudinal direction. This serves as a test of both the analysis and the numerical accuracy: the image of an object should be identical to the image of its measurement component. As can be seen in Fig. 8, the difference between them is below 1% except at the edges were edge effects could cause it to rise to around 10%.

An interesting point is how the slice through the measurement component is a hybrid between the projection and the similar object slice. This illustrates the very principle of tomosynthesis. Unlike tomography, which provides a full reconstruction of every plane of the object, tomosynthesis provides only partial 3D information. The longitudinal resolution will always be low, and contributions from close-by transverse planes will be present in every slice of the measurement component.

The measurement component of an object provides information on how well this object is transmitted through the imaging system. Comparing the measurement component in Fig. 6 to its object in Fig. 5, we can see that the information received in the image plane is not complete. This illustrates the fundamental limit of the acquisition process, and serves as an estimate of the best reconstruction that could be achieved. We can also see that the measurement component looks a lot like the object, but that it also borrows traits from the projection in Fig. 7. This reflects the small angle used in tomosynthesis and its effect on the results. The reconstruction, while given in 3D, is in many ways an improved projection or a projection with some separation of information into different layers. It is not the full 3D reconstruction that could be achieved from e.g. tomography.

It is also becoming evident that there is a lack of three-dimensional models of the breast tissue. Statistical models for the normal breast, such as the clustered lumpy background [23], are made for 2D and not for 3D. Since mammography has been the dominating investigation method, it has simply been sufficient to work only on the 2D projection through the breast. Extensions to 3D are either intuitive and lack the statistical back-up required for a systematic study, like the example presented in this paper, or advanced models designed to produce one really good phantom and not the large number of samples required for a statistical study. There is need for a model of the breast tissue that has the same statistical properties as experimental images in 3D but also in its 2D projections, that looks similar to experimental images in 3D but also in its 2D projections, and that allows for generation of numerous realizations. Strong candidates at the moment are a more sophisticated extension of the clustered lumpy background [23], the threshold model by Abbey and Boone [25], or the breast phantoms by Zhang et al. [26]. The example object in Fig. 5 is sufficient as an illustration of the SVD method, but is not interconnected enough for a real slice through breast tissue.

The next step would be to use the SVD of a more realistic breast tissue model in an image-quality study, including the effect of noise. The SVD model is particularly suitable for this since the statistics of the object background are relatively easily transferred to the image through the use of characteristic functionals [11, Sec. 8.5.3]. The SVD singular functions can also be used as channels for a channelized observer [27]. This allows for many comparisons, for example how or if the detection seems to improve with tomosynthesis compared to mammography, or how the SVD method compares to other reconstruction algorithms.

The mixed continuous and discrete operators of a tomosynthesis system have been derived from the Boltzmann equation, and used to develop an analytical singular-value decomposition of a tomosynthesis system. In combination with numerical calculations, it provides the measurement component of any object, where object refers to the breast tissue. The measurement component of an example object has been presented. Although the example object is not truly realistic, interesting information on the tomosynthesis process and its effects on reconstruction have been found. Future prospects involve more realistic 3D breast tissue models, and an image-quality study involving a channelized Hotelling observer.

C. Dainty acknowledges support from Science Foundation Ireland Grant No 07/IN.1/1906, H. Barrett from the Center for Gamma-Ray Imaging Grant 5 P41EB002035-11 and the SPECT Imaging and Parallel Computing Grant 5 R37 EB000803-19, and A. Burvall from the Swedish Research Council.

1. Dobbins JT, III, Godfrey DJ. Digital x-ray tomosynthesis: current state of the art and clinical potential. Phys Med Biol. 2003;48:R65–R106. [PubMed]

2. Dobbins JT., III Tomosynthesis: at translational crossroads. Med Phys. 2009;36:1956–1967. [PubMed]

3. Niklason LT. Digital tomosynthesis in breast imaging. Radiology. 1997;1997:399–406. [PubMed]

4. Gennaro G, Toledano A, di Maggio C, Baldan E, Bezzon E, La Grassa M, Pescarini L, Polico I, Proietti A, Toffoli A, Muzzio PC. Digital breast tomosynthesis versus digital mammography: a clinical performance study. Eur Radiol. 2009 doi: 10.1007/s00330-009-1699-5. [PubMed] [Cross Ref]

5. Andersson I, Ikeda DM, Zackrisson S, Ruschin M, Svahn T, Timberg P, Tingberg A. Breast tomosynthesis and digital mammography: a comparison of breast cancer visibility and BIRADS classification in a population of cancers with subtle mammographic findings. Eur Radiol. 2008;18:2817–2825. [PubMed]

6. Good WF, Abrams GS, Catullo VJ, Chough DM, Ganott MA, Hakim CM, Gur D. Digital breast tomosynthesis: a pilot observer study. Am J Radiology. 2008;190:865–869. [PubMed]

7. Poplack SP, Tosteson TD, Kogel CA, Nagy HM. Digital breast tomosyntheis: Initial experiance in 98 women with abnormal digital screening mammography. Am J Radiology. 2007;189:616–623. [PubMed]

8. Chawla AS, Lo JY, Baker JA, Samei E. Optimized image acquisition for breast tomosynthesis in projection and reconstruction space. Med Phys. 2009;36:4859–4869. [PubMed]

9. Wu T, Moore RH, Rafferty EA, Kopans DB. A comparison of reconstruction algorithms for breast tomosynthesis. Med Phys. 2004;31:2636–2647. [PubMed]

10. Zhang Y, Chan HP, Sahiner B, Wei J, Goodsitt MM, Hadjiiiski LM, Ge J, Zhou C. A comparative study of limited-angle cone-beam reconstruction methods for breast tomosynthesis. Med Phys. 2006;33:3781–3795. [PMC free article] [PubMed]

11. Barrett HH, Myers KJ. Foundations of Image Science. John Wiley; Hoboken, New Jersey: 2004.

12. Bertero M, Boccacci P. Inverse Problems in Imaging. Institute of Physics Publishing; Bristol, UK: 1998.

13. Barrett HH, Aarsvold JN, Roney TJ. Null functions and eigenfunctions: tools for the analysis of imaging systems. Lect Notes Comput Sci. 1991;11:211–226. [PubMed]

14. Burvall A, Barrett HH, Dainty C, Myers KJ. Singular-value decomposition for through-focus imaging systems. J Opt Soc Am A. 2006;23:2440–2448. [PubMed]

15. Yao J, Barrett HH. Predicting human performance by a channelized Hotelling observer. Proc SPIE. 1992;1768:161–168.

16. Barrett HH, Yao J, Rolland JP, Myers KJ. Model observers for assessment of image quality. Proc Natl Acad Sci USA. 1993;90:9758–9765. [PubMed]

17. Chiu MY, Barrett HH, Simpson RG, Chou C, Arendt JW, Gindi GR. Three-dimensional radiographic imaging with a restricted view angle. J Opt Soc Am. 1979;69:1323–1333.

18. Pierri R, Liseno A, Soldovieri F, Solimene R. In-depth resolution for a strip source in the Fresnel zone. J Opt Soc Am A. 2001;18:352–359. [PubMed]

19. Polianin AD, Manzhirov AV. Handbook of integral equations. chapter 11.2 CRC Press; Florida: 1998.

20. Lanczos C. Linear differential operators. Van Nostrand; London: 1961.

21. Burgess AE. Visual signal detection with two-component noise: low-pass spectrum effects. J Opt Soc Am A. 1999;16:694–704. [PubMed]

22. Reiser I, Nishikawa RM. Task-based assessment of breast tomosynthesis: Effects of acquisition parameters and quantum noise. Med Phys. 2010;37:1591–1600. [PubMed]

23. Bochud FO, Abbey CK, Eckstein MP. Statistical texture synthesis of mammographic images with clustered lumpy backgrounds. Opt Express. 1998;4:33–43. [PubMed]

24. Metheany KG, Abbey CK, Packard N, Boone JM. Characterizing anatomical variability in breast CT images. Med Phys. 2008;35:4685–4694. [PubMed]

25. Abbey CK, Boone JM. An ideal observer for a model of x-ray imaging in breast parenchymal tissue. In: Krupinski EA, editor. IWDM 2008, LNCS. Vol. 5116. 2008. pp. 393–400.

26. Zhang C, Bakic PR, Maidment ADA. Development of an anthropomorphic breast software phantom based on region growing algorithm. Proc SPIE. 2008;6918:69180V.

27. Park S, Witten JM, Myers KJ. Singular vectors of a linear imaging system as efficient channels for the bayesian ideal observer. IEEE Trans Med Imaging. 2009;28:657–668. [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. |