# Estimating view parameters from random projections for Tomography using spherical MDS

- Yi Fang
^{1}Email author, - Sundar Murugappan
^{1}and - Karthik Ramani
^{1, 2}

**10**:12

**DOI: **10.1186/1471-2342-10-12

© Fang et al; licensee BioMed Central Ltd. 2010

**Received: **26 December 2009

**Accepted: **18 June 2010

**Published: **18 June 2010

## Abstract

### Background

During the past decade, the computed tomography has been successfully applied to various fields especially in medicine. The estimation of view angles for projections is necessary in some special applications of tomography, for example, the structuring of viruses using electron microscopy and the compensation of the patient's motion over long scanning period.

### Methods

This work introduces a novel approach, based on the spherical *multidimensional scaling (sMDS)*, which transforms the problem of the angle estimation to a sphere constrained embedding problem. The proposed approach views each projection as a high dimensional vector with dimensionality equal to the number of sampling points on the projection. By using SMDS, then each projection vector is embedded onto a 1D sphere which parameterizes the projection with respect to view angles in a globally consistent manner. The parameterized projections are used for the final reconstruction of the image through the inverse radon transform. The entire reconstruction process is non-iterative and computationally efficient.

### Results

The effectiveness of the sMDS is verified with various experiments, including the evaluation of the reconstruction quality from different number of projections and resistance to different noise levels. The experimental results demonstrate the efficiency of the proposed method.

### Conclusion

Our study provides an effective technique for the solution of 2D tomography with unknown acquisition view angles. The proposed method will be extended to three dimensional reconstructions in our future work. All materials, including source code and demos, are available on https://engineering.purdue.edu/PRECISE/SMDS.

## Background

This work studies the problem of 2D tomography with unknown view angles and discusses the potential applications of our work. We give the background of our work, reviews of the existing methods and a brief introduction of our proposed method in the following subsections.

### Tomography with unknown view angles

The computed tomography has been successfully applied to various fields over the past decades, for example, medical imaging, synthetic aperture radar (SAR) and Cryo-electron microscopy (cryoEM) for structuring viruses [1–4]. The traditional tomography is defined as a process of recovering the object from the measurements that are line integrals of that object at some set of known orientations (view angles). However, in some special situations, obtaining the view angles is difficult or suffers from accurate measurement. For example, 1) the patient's motion owing to long scanning period can result in uncertainty of view angles, 2) the data acquisition of single particle cryoEM are the line integrals of many identical copies of virus molecules at random orientations. Therefore, the research of a more generalized tomography independent of known acquisition view angles is worthwhile to study for the reconstruction of the objects under various circumstances. Our overarching goal is the reconstruction of 3D virus from 2D cryoEM images. As one step towards the goal, we address the problem of image alignment (orientation determination). In the application of cryoEM, the projections for the macromolecules at a preset angle are captured for a large number of identical macromolecules at different unknown and random orientations. If we assume there is no overlapping o the projection, the imaging process is the same as having projections of the macromolecule at multiple but unknown angles while fixing the position of a single homogeneous macromolecule [5]. Hence, like [5–7], we propose a computational model which is equivalent to the real-world scenario, as having projections at different, unknown and random angles while fixing the position of a single homogeneous macromolecule.

Recently, the uniqueness and feasibility of tomography with unknown view angles has been discussed in [3, 8]. An object can be uniquely determined using the projection data with a certain number of unknown and distinct view angles. The reconstructed object is subjected to a global arbitrary spatial rotation, which doesn't influence the study of the investigated objects. A review of recently proposed methods dealing with this issue is described below.

### Review of existing methods

Several approaches for analyzing projections measured from unknown view angles have been developed over last two decades [6, 8–18]. Those methods can be roughly categorized into two classes, iterative and non-iterative methods. The iterative methods described in [8, 18] use the moment characterization of the range space of the Radon Transform, known as the Helgasson-Ludwig (HL) consistency conditions to reconstruct the image. The authors proposed a Bayesian approach for the view angle estimation for tomography in [16]. An integrated statistical technique for volume reconstruction with unordered sequential slices is presented in [17]. The limitations of the iterative methods cited above arise because of a huge computation complexity owing to the solution of a large nonlinear problem at each iteration [5]. Recently, non-iterative approaches have been developed in [5, 6, 13] to achieve a fast way of structuring the object from its projections with unknown view angles. Yagle introduced a simple non-iterative algorithm based on circular harmonic expansion in [5]. The work decouples the view angles estimation problem from the image estimation problem and thereby largely reduces the computational expenses. Two manifold learning based techniques showed great performance at solving the view angles uncertainty [6, 13]. Georg etc. [13] applied the manifold learning to sort the time-ordered slab data for automatic estimation of lung volume without any external breath measurements. The goal is to piece together the local slab data on a proper position in a globally consistent manner, with regards to breath phase circle. Coifman et al. [6] presented a Laplacian graph based manifold learning method to enforce the view angles embedded on the circle. All of the non-iterative approaches have a common attractive property: computation is really efficient as the angle estimation problem is transformed into a matrix eigenvalue problem with the size equal to the number of angles.

### sMDS for view angles uncertainty

The manifold learning methods show attractive properties in terms of estimating view angles from the acquisition data [6, 13, 19]. As mentioned in [19], the popular manifold learning methods such as, Multidimensional scaling (MDS) [20], ISOMAP [21], and locally linear embedding (LLE) [22], cannot handle the view angle uncertainty problem directly because these methods can only embed points in a flat space while the view angles are intrinsically distributed on a sphere. To solve this problem, our study introduces a reconstruction scheme based on spherical MDS (sMDS) proposed in [23], which is able to embed points on a spherical manifold. Our work is inspired by the projection-slice theorem, which states that the Fourier transform of the projection of a 2D function is equal to a slice through the origin of the 2D Fourier transform of that function. The algorithm for estimation of view angle for each projection consists of three steps. First, Fourier transform is applied to all of the projection data (e.g.
denotes the Fourier transform of the *ith* projection). Then the distance between the pairwise Fourier transform (e.g.
and
) is measured to build the distance matrix. Last, sMDS is applied on this distance matrix to estimate the view angle for each projection data. The algorithm details are presented in Section 2. Essentially the embedding process assigns each projection a point in a low dimensional intrinsic parameter space. This process assembles the similar projections that are close to each other without using any prior knowledge. There are two reasons to build the distance matrix in the Fourier domain 1) it is more flexible to perform the computations on the Fourier domain as it helps us to handle noise by allowing us to choose a proper range of frequency instead of the complete range. This is even useful when the noise level is high as the signal in low frequency range is comparable to noise while the signal in high frequency range is totally buried by noise. Therefore, the distance computed in the Fourier domain is more robust to noise by controlling the range of frequency in Fourier domain. The algorithm for choosing an optimal range in Fourier domain is one of our future work. Currently, we use either half of the range or the full range itself if noise is low. 2) the magnitude of the Fourier vector is invariant to the center shift of the image. Thus, even when the image shifts during projecting, the distance would not be affected if the computation is only based on the magnitude of the Fourier values.

## Methods

### Spherical MDS

- 1.
Build the pairwise distance matrix

**M**, and - 2.
Compute the dot-product form of the distance matrix , where Γ denotes a operator applied on the distance matrix

**M**and*r*denotes the radius of the sphere calculated using , and - 3.
Choose the first the

*k*+ 1 eigenvectors of the distance matrix Γ(**M**). The*ith*eigenvector is denoted as = (*V*_{ i }(1),*V*_{ i }(2)), ...,*V*_{ i }(*n*) where n is the number of the points. The coordinate of the*jth*embedded point is = (*V*_{1}(*j*),*V*_{2}(*j*), ...,*V*_{ k+1}(*j*)) where*k*is the dimensionality of the sphere.

This procedure briefly introduces the general ideas about the sMDS. The details about how to apply sMDS for 2D tomography application is provided in the following sections.

### Problem definition

#### Fourier slice theorem

*P*

_{ θ }), from a view angle, mathematically, an integral of the object density function along the sight parallel line. The figure on the right shows the two dimensional fourier transform of the object. The green color points represent a slice of 2D Fourier transform and the projection

*P*

_{ θ }is a one-dimensional Fourier transform pair according to the theorem.

#### Derivation of the Model

##### Projection and its Fourier transform

Each projection data is viewed as a high dimensional vector (PVector)
= (*P*(*t*
_{1}), *P*(*t*
_{2}), ..., *P*(*t*
_{
n
})) where *t*
_{1}, *t*
_{2}, ..., *t*
_{
n
}are equally spaced sampling points on the projection data and *n* is the dimensionality of the vector. The Fourier transform of the PVector is also viewed as a high dimensional vector (FTVector), represented as
(*F*(*ρ*
_{1}), *F*(*ρ*
_{2}), ..., *F*(*ρ*
_{
n
})) where *ρ*
_{1}, *ρ*
_{2}, ..., *ρ*
_{
n
}are represented as equally spaced sampling points on the Fourier transform of the projection data and *n* is the dimensionality of the vector. According to the theorem, each projection collected from unknown view angle has a unique Fourier transform pair, which is a slice of the 2D Fourier transform of the object. Therefore, the problem of sorting projections in space domain is equivalent to orienting their corresponding Fourier transforms in frequency domain.

##### Orienting Fourier transforms

The rearrangement of Fourier transforms is essentially a dimensionality reduction problem with the internal structure constrained by spherical manifold. As we can see in Figure 2, the slices are distributed along the radial direction. We can reason analytically and imagine that a FTVector (a slice),
is intrinsically restricted on a circle since a line passing through the origin could be uniquely parameterized by its orientation. We provide an observation in the appendix section to explain why the underlying structure is a 1D sphere (a circle). Mathematically, the FTVector
= (*F*
_{
i
}(*ρ*
_{1}), *F*
_{
i
}(*ρ*
_{2}), ..., *F*
_{
i
}(*ρ*
_{
n
})) could be intrinsically reduced to a two dimensional point
= (*X*
_{
i
}, *Y*
_{
i
}) on a circle, where i denotes the *ith* projection and *X* and *Y* denote the principal axes in cartesian coordinate. The FTVectors can be oriented based on the corresponding two dimensional point set {
, *i* ∈ (1, 2, ..., *n*)} where n is the number of points.

##### Orienting the projections

As we discussed above, one projection is uniquely associated with a slice in the 2D Fourier transform space. The low dimensional intrinsic parameter for each FTVector could be assigned to the corresponding PVector directly. The embedded point set reveals the orientation of the projection in a coherent global manner, which clearly organize the projections with unknown view angles. Note that the estimated organization is subjected to a global arbitrary rotation, which doesn't influence the analysis of the reconstruction result.

This model basically converts the view angle uncertainty problem to a dimensionality reduction problem and the reduced low dimensional structure is used to recover the view angles.

### View angles estimation algorithm

#### Pairwise distance matrix

In this work, the pairwise distance matrix is built by following the steps described in [21]. The computation of the pairwise distance matrix consists of two steps. First, given a distance threshold **T**, the high dimensional points (e.g. the FTVector in this work) which are neighbors, are determined based on distance *d*(*A*, *B*) between pairs of points *A*, *B* ( Eq.1). A particular point connects to other points if the pairwise distance d(A, B) is less than the predefined threshold **T**. A weighted graph **G** with weight *d*(*A*, *B*) between neighboring points can be used to describe these neighborhood relations. Second, the Dijkstra's algorithm [24] is applied to the weighted graph **G** to compute the shortest path distance in the graph, which well approximates the geodesic distance between all pairs of the points.

*L*1 norm,

*L*2 norm, λ

^{2}measure, and Bhattacharyya distance. In this work, we adopt a widely used distance measurement,

*L*2 norm.

where and denote differently oriented FTVector respectively.

#### Estimation of view angles

##### Determination of coordinates

**M**), intrinsically reflect the relative positions of the embedded points in a globally consistent way. The coordinates of the embedded points are solved by the following procedure.

- 1.
Rank the eigenvalues of the Γ(

**M**) in decreasing order and choose first two eigenvectors, and , corresponding to the first two eigenvalues, and - 2.
Normalize each pair (

*V*_{1}(*i*),*V*_{2}(*i*)) to a unit vector, where*i*is the index of the embedded points. The normalized vectors will be used as the coordinates of the embedded point on the circle.

##### Determination of view angles

- 1.Apply the inverse triagonometric function, arctangent, to the coordinates of embedded points for calculating the initial set of view angles,(2)

*ith*point,

*Y*

_{ i }and

*X*

_{ i }denote the coordinate of the

*ith*embedded point and

*φ*

_{ i }the estimated view angle for the

*ith*projection.

- 2.
Sort the initial set of view angles

*φ*_{1},*φ*_{2}, ...,*φ*_{ n }) in ascending order, uniformly rearrange the view angles along the circle, and associate the refined view angle set to the original projection data.

## Results

We have implemented sMDS for 2D tomography with unknown view angles and assessed its performance from the experimental results. The algorithms presented in the paper are implemented on a Pentium D 3.2 GHz computer with 1 G RAM running Windows XP. The images for the experiments have been chosen from the Whole Brain Atlas http://www.med.harvard.edu/AANLIB/home.html. The Whole Brain Atlas provides a large set of MR images for both normal and diseased brain. There are three parameters in the experimental setting, the number of the projections, **N**, the threshold, **T**, which is used to determine the neighborhood relation between a pair of points, and the number of the sampling points along each projection, **S**.

### Orienting projection data

**N**, is set at 360 and the number of the sampling points of each projection is set at 299. We used a normal brain MR image as the experimental object, shown in Figure 3. Figure 4 shows the visualization of the projection data. The Figure 4(A) illustrates the projections with random view angles and Figure 4(B) shows the sorted projection data. The horizontal axis in the figure is labelled with the index of projection data, and the vertical axis of the histogram is the number of the sampling points of each projection. The colorbar on the right side indicates the value of the projection data at each sampled position. The value range of horizontal axis is between 1 and 360, indicating that there are total 360 projections acquired from different view angles. The value range of vertical axis is from 1 to 299, indicating that number of the sampling points,

**S**, set at 299 for each projection data.

A visual comparison between the two figures demonstrates that sMDS sorts the projection data well. We can see from the Figure 4(A) that the projection data are shuffled owing to the random projections, while in the Figure 4(B) the projection data are clearly ordered with the global relative orientation.

### Verification of reconstruction quality

**N**, and threshold value,

**T**. Figure 5 shows reconstructed results of the brain MR image (see Figure 3) using different number of projections. The threshold,

**T**, used for the neighbors detection are slightly varied according to the number of the projections. With the increase in the number,

**N**, the distance metric between pairwise points decrease. Therefore, we need to lower the value of the threshold to assure that a proper neighbor relation for each projection is preserved. The reconstruction quality is enhanced with the increase in the number of the projections(see Figure 3). In addition, we observe that the reconstructed images are subjected to an arbitrary rotation of the original image. We further verify quantitatively the reconstruction performance by comparing the original image with the reconstructed image. In our test, we cannot directly subtract the original from the reconstructed image as there is an arbitrary rotation of the reconstructed image. However, the registration of the two images can remove the effect of the arbitrary global rotation. In this experiment, the registration between original and reconstructed images is trivial as the correspondence between the randomly shuffled projections and the re-organized projections can be easily tracked. The global rotation angle can be easily retrieved by finding the difference between any corresponding pairs. For example, we can track the first projection of the original image, find its relative position in the re-organized projections sequence, then compute the rotation angle. We provide two quantitative measures for the reconstruction performance: one is the peak signal-to-noise ratio (PSNR), and the other is mean squared error (MSE), which are calculated using Eq.3 and Eq.4.

### Resistance to noise

## Discussion

### Analysis of experimental results

The experiments conducted in this work verify the performance of our proposed method in dealing with the 2D Tomography with unknown view angles. The number of the projection data **(N)** is one of the most important parameters in our experiments. The manifold learning methods generally require a sufficient number of samples. To meet this requirement, **N** should be relatively big enough, e.g. 256. Secondly, the threshold for neighbors detections is extremely important for our method. Since the sMDS is essentially a method piecing together the local information in a global manner, the local or neighbor relations between points ultimately determine the global embedding. In our experiment, the setting of the threshold is mainly based on **N**. As we understand, larger the value **N**, less the distance between them, and thus a smaller threshold is chosen to build a proper set of neighbors.

### Applications

The estimation of view angles for projections is necessary in some special applications of tomography. We discuss the potential applications of our proposed method in this section. Since the main contribution of this work is the tomography with unknown view angles, this work would be crucial to address the needs of patients in designing next generation tomography equipment. Various studies have shown that it could make the analysis less complicated if the patients are at ease and less anxious. Nowadays, the patients are required to remain motionless during a long scanning period. The difficulty of being motionless and the increasing discomfort and anxiety of the patients lead to the measurement error and analysis complication. Our method shows great potential to minimize these constraints. In addition, there are some situations where the view angles or the acquisition positions of projections are really hard or not possible to be known to us, for instance, the macromolecule structure determination by electron microcopy. Our proposed method offers a very promising approach to reconstruct the object efficiently.

### Limitations

The estimation of view angles or acquisition positions could be categorized as the inverse problem, which is defined as the inference of model parameters from the observed data. There is an unavoidable limitation of methods for the estimation view angles, including the method presented in this work. A perfectly symmetrical object would lead to the fact that a certain number of projections from different view angles are identical to each other. In this case, it is not possible to associate the projections with different view angles because the projection data themselves are not distinguishable, even though they are obtained from different view angles. To overcome this problem, the pre-estimation of the symmetry of the investigated object could be combined with the view angles estimation procedure. There exist some works for detection of the symmetry of an object in computer vision, image processing and computer graphics areas [26–28]. Our future work is to extend our method to three dimensional tomography with unknown view angles and combine the symmetry detection for the symmetrical object reconstruction.

## Conclusions

The uniqueness and feasibility of tomography with unknown view angles have been proved in the earlier works, which offer the theoretical fundamentals for our method. We have introduced an efficient reconstruction procedure for 2D tomography with unknown view angles by means of sMDS. The experimental results indicate our method performs well with high quality image reconstruction, even in the case of highly corrupted noisy data. Our method would potentially provide an alternative approach in dealing with special cases of tomography with unknown acquisition parameters.

## Observation

## Declarations

### Acknowledgements

We would like to thank the Institute for Pure & Applied Mathematics for providing fellowship to support Karthik Ramani's visit. We would also like to thank Amit Singer, Yoel Shkolnisky, Ronald Coifman and Peter W. Jones for providing valuable insights on graph laplacian. The experimental images are provided by the whole brain database from Harvard university. This material is partly based upon work supported by the National Institute of Health (GM-075004). Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the National Institute of Health.

## Authors’ Affiliations

## References

- Hiriyannaiah H: X-ray computed tomography for medical imaging. Signal Processing Magazine, IEEE. 1997, 14 (2): 42-59. 10.1109/79.581370.View ArticleGoogle Scholar
- Smith BD: Cone-beam tomography: recent advances and a tutorial review. Optical Engineering. 1990, 29:Google Scholar
- Basu S, Bresler Y: Feasibility of Tomography with Unknown View Angles. IEEE Transactions on Image Processing. 2000, 9: 1107-1122. 10.1109/83.846252.View ArticlePubMedGoogle Scholar
- Lewitt RM, Matej S: Overview of Methods for Image Reconstruction From Projections in Emission Computed Tomography. Proceedings Of The IEEE. 91: 1588-1611. 10.1109/JPROC.2003.817882.
- Yagle AE: A simple non-iterative algorithm for 2-D tomography with unknown view angles. Appeared as a preprint copy. [http://www.eecs.umich.edu/~aey/recent/angle.pdf]
- Coifman RR, Shkolnisky Y, Sigworth FJ, Singer A: Graph Laplacian Tomography From Unknown Random Projections. IEEE Transactions on Image Processing. 2008, 17: 1891-1899. 10.1109/TIP.2008.2002305.View ArticlePubMedGoogle Scholar
- Borg I, Groene P: Introduction to 3D reconstruction of macromolecules using single particle electron microscopy. Acta Pharmacologica Sinica. 2005, 26: 1153-1164. 10.1111/j.1745-7254.2005.00203.x.View ArticleGoogle Scholar
- Basu S, Bresler Y: Uniqueness of Tomography with Unknown View Angles. IEEE Transactions on Image Processing. 2000, 9: 1094-1106. 10.1109/83.846251.View ArticlePubMedGoogle Scholar
- Zeng R, Fessler JA, Balter JM, Balter PA: Iterative sorting for four-dimensional CT images based on internal anatomy motion. Medical Physics. 2008, 35: 917-926. 10.1118/1.2837286.View ArticlePubMedPubMed CentralGoogle Scholar
- Zeng R, Fessler J, Balter J: Estimating 3-D respiratory motion from orbiting views by tomographic image registration. IEEE Transcations on Medical Imaging. 2007, 26: 153-163. 10.1109/TMI.2006.889719.View ArticlePubMedGoogle Scholar
- Lauren PD, Nandhakumar N: Estimating the Viewing Parameters of Random, Noisy Projections of Asymmetric Objects for Tomographic Reconstruction. IEEE Transactions on Pattern Analysis and Machine Intelligence. 1997, 19: 417-430. 10.1109/34.589202.View ArticleGoogle Scholar
- Radermacher M: Three Dimensional Reconstruction From Random Projections: Oriented Alignment Via Radon Transforms. Ultramicroscopy. 1994, 53: 121-136. 10.1016/0304-3991(94)90003-5.View ArticlePubMedGoogle Scholar
- Georg M, Souvenir R, Hope A, Pless R: Manifold learning for 4D CT reconstruction of the lung. Computer Vision and Pattern Recognition Workshops, CVPRW '08. 2008, 1-8. full_text.Google Scholar
- Heel MV: Angular Reconstitution: A Posteriori Assignment of Projection Directions For 3D Reconstruction. Ultramicroscopy. 1987, 21: 111-123. 10.1016/0304-3991(87)90078-7.View ArticlePubMedGoogle Scholar
- Vainshtein B, Goncharov A: Determination of the Spatial Orientation of Arbitrarily Arranged Identical Particles of an Unknown Structure From Their Projections. Dokl Akad Nauk SSSR. 1986, 287: 1131-1134.Google Scholar
- Mallick S, Agarwal S, Kriegman D, Belongie S, Carragher B, Potter C: Structure and View Estimation for Tomographic Reconstruction: A Bayesian Approach. Computer Vision and Pattern Recognition, 2006 IEEE Computer Society Conference on. 2006, 2: 2253-2260.Google Scholar
- Yu Y, Trouve A, Wang J, Chalmond B: An integrated statistical approach for volume reconstruction from unregistered sequential slices. Inverse Problems. 2008, 24 (4): 10.1088/0266-5611/24/4/045004.
- Salzman D: A method of general moments for orienting 2D projections of unknown 3D objects. Computer Vision, Graphics, and Image Processing. 1990, 50: 129-156. 10.1016/0734-189X(90)90038-W.View ArticleGoogle Scholar
- Coifman RR, Shkolnisky Y, Sigworth FJ, Singer A: Cryo-EM Structure Determination through Eigenvectors of Sparse Matrices. Technical Report 1389, Yale University, Department of Computer Science. 2007Google Scholar
- Borg I, Groene P: Modern Multidimensional Scaling: Theory and Application. 1997, Springer-VerlagView ArticleGoogle Scholar
- Tenenbaum JB, de Silva V, Langford JC: A Global Geometric Framework for Nonlinear Dimensionality Reduction. Science. 2000, 290: 2319-2322. 10.1126/science.290.5500.2319.View ArticlePubMedGoogle Scholar
- Roweis S, Saul L: Nonlinear dimensionality reduction by locally linear embedding. Science. 2000, 290: 2323-2326. 10.1126/science.290.5500.2323.View ArticlePubMedGoogle Scholar
- Pless R, Simon I: Embedding images in non-flat spaces. Proc. of the International Conference on Imaging Science, Systems, and Technology. 2001Google Scholar
- Dijkstra EW: A note on two problems in connexion with graphs. In Numerische Mathematik. 1959, 1: 269-271. 10.1007/BF01386390.View ArticleGoogle Scholar
- Osada R, Funkhouser T, Chazelle B, Dokin D: Shape distributions. ACM Transactions on Grphics. 2002, 33: 133-154.Google Scholar
- Zabrodsky H, Peleg S, Avnir D: Symmetry as a continuous feature. IEEE Transactions on Pattern Analysis and Machine Intelligence. 1995, 17: 1154-1166. 10.1109/34.476508.View ArticleGoogle Scholar
- Keller Y, Shkolnisky Y: A signal processing approach to symmetry detection. IEEE Transactions on Image Processing. 2006, 15: 2198-2207. 10.1109/TIP.2006.875227.View ArticlePubMedGoogle Scholar
- Mitra NJ, Guibas L, Pauly M: Partial and Approximate Symmetry Detection for 3D Geometry. 2006, 25 (3): 560-568.Google Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2342/10/12/prepub

### Pre-publication history

## Copyright

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