International Journal of Biomedical Imaging

International Journal of Biomedical Imaging / 2010 / Article
Special Issue

Mathematical Methods for Images and Surfaces

View this Special Issue

Research Article | Open Access

Volume 2010 |Article ID 291874 | https://doi.org/10.1155/2010/291874

Xiaowei He, Jimin Liang, Xiaochao Qu, Heyu Huang, Yanbin Hou, Jie Tian, "Truncated Total Least Squares Method with a Practical Truncation Parameter Choice Scheme for Bioluminescence Tomography Inverse Problem", International Journal of Biomedical Imaging, vol. 2010, Article ID 291874, 11 pages, 2010. https://doi.org/10.1155/2010/291874

Truncated Total Least Squares Method with a Practical Truncation Parameter Choice Scheme for Bioluminescence Tomography Inverse Problem

Academic Editor: Guo Wei
Received29 Sep 2009
Revised04 Mar 2010
Accepted08 Mar 2010
Published19 May 2010

Abstract

In bioluminescence tomography (BLT), reconstruction of internal bioluminescent source distribution from the surface optical signals is an ill-posed inverse problem. In real BLT experiment, apart from the measurement noise, the system errors caused by geometry mismatch, numerical discretization, and optical modeling approximations are also inevitable, which may lead to large errors in the reconstruction results. Most regularization techniques such as Tikhonov method only consider measurement noise, whereas the influences of system errors have not been investigated. In this paper, the truncated total least squares method (TTLS) is introduced into BLT reconstruction, in which both system errors and measurement noise are taken into account. Based on the modified generalized cross validation (MGCV) criterion and residual error minimization, a practical parameter-choice scheme referred to as improved GCV (IGCV) is proposed for TTLS. Numerical simulations with different noise levels and physical experiments demonstrate the effectiveness and potential of TTLS combined with IGCV for solving the BLT inverse problem.

1. Introduction

In recent years, molecular imaging has emerged as a promising tool in basic, preclinical and clinical research for monitoring a variety of molecular and cellular processes in living organisms [14]. As one of molecular imaging modality, bioluminescence tomography (BLT) has attracted much attention due to its exquisite sensitivity and cost effectiveness.

The key problem of BLT is to reconstruct the bioluminescent source distribution inside a biological tissue from the optical signals detected on the body surface, which is a highly ill-posed inverse problem. By using numerical method such as finite element method (FEM), the inverse problem of BLT can be formulated into a nonsquare matrix equation, where the coefficient matrix is typically ill-conditioned [5]. Hence overcoming the ill-posedness and seeking a stable solution of the matrix equation are the major issues of BLT inverse problem. For this purpose, the inverse problem is often transformed to a least squares problem incorporated with the regularization technique. Tikhonov regularization is the most widely-used method in BLT reconstruction [68]. It is aiming to stabilize the inverse of an ill-conditioned operator and minimize the effects of the inevitable error by minimizing a trade-off between the loss function and the -norm [8]. However, previous studies based on Tikhonov regularization only consider noise in the measurement. In fact, some system errors also exist in the computed coefficient matrix of the system equation. These errors may take place in such aspects as FEM discretization, geometrical mismatch, optical parameters inaccuracy and model approximation, and so forth. System errors as well as measurement noise are inevitable in real BLT experiment, which may lead to large errors for the reconstruction results.

The total least squares (TLS) method is a generalization of the least squares approximation method when the data in both sides of the matrix equation are perturbed [9, 10]. Based on the TLS method, the truncated total least squares (TTLS) method is proposed for regularization of ill-conditioned linear systems [11]. It is inspired by truncated singular value decomposition (TSVD) which aims at limiting the contribution of noise by cutting off a certain number of terms in the singular value decomposition of coefficient matrix [12]. Truncation level plays the role of regularization parameter in truncation methods, which has great influence on the quality of the solution. As a result, determining an appropriate truncation level for TTLS is a critical step in the inverse procedure. Most existing parameter-choice schemes such as L-curve, discrepancy principle, generalized cross-validation (GCV), and zero crossing methods assume that the coefficient matrix is exactly known, that is, it is not contaminated by noises or errors [1316]. In [17], a truncation level choice criterion named modified GCV (MGCV) is proposed for TTLS method; theoretical analysis and simulation tests show its potential for solving ill-posed linear system. However, it has been recognized that choice schemes of regularization parameter are mostly problem-dependent and practical parameter-choice scheme for BLT reconstruction deserves further study.

In this paper, the aim of our study is to extend the BLT reconstruction to the case including both the measurement noise and the system errors. For this purpose, TTLS method combined with a practical scheme termed as improved GCV (IGCV) is proposed to solve the BLT inverse problem. In the next section, our methodology of solving the inverse problem in BLT is described. In Section 3, we demonstrate the performance of the TTLS method combined with IGCV scheme in BLT reconstruction using numerical simulation and physical experiments in various source and noise level settings. Finally, we draw a conclusion and discuss the relevant issues.

2. Methodology

2.1. Diffusion Approximation and Boundary Condition

In general, light propagation in living subjects is mainly hindered by both tissue scattering and absorption [7, 8]. Considering that bioluminescent photons belong to the near-infrared region where scattering predominates over absorption [3], the propagation of photon can be well modeled by the following steady-state diffusion equation [18]: where is the bounded domain, represents the photon flux density, and denotes the energy density distribution of an internal bioluminescence source, is the optical diffusion coefficient with being the optical absorption coefficient and the reduced scattering coefficient, respectively.

Assuming that the BLT experiment is performed in a totally dark environment, the equation is subject to a Robin boundary condition [18]: where denotes the boundary, represents the unit outer normal on and , n is the ratio of optical reflective index of the inner tissue to that outside the boundary, and is close to 1.0 when the subject is in air [8]. In a bioluminescent imaging experiment, the measurable photon flux density on can be calculated by the following outgoing radiation [18]:

2.2. The Model of BLT Reconstruction

Based on (1), (2), and (3), the essence of the BLT reconstruction is to estimate the light source distribution inside the biological tissues from the measured flux on the surface, given the corresponding optical parameters of the tissues. In order to solve the BLT inverse problem, FEM was introduced to solve the diffusion equation in [8, 1820] because of its capability to process volume with arbitrary geometries. After the discretization using FEM, the linear relation between the bioluminescence source intensity and the photon flux density can be expressed as the following matrix form: where and are the collection of all the nodal values of the photon flux density and source density, is a positive-definite matrix, and K, C, and B are called the mass, stiff, and boundary matrix, respectively. The photon density can be obtained from . In fact, only partial photon on the boundary can be acquired in the BLT experiment, therefore, can be partitioned into the measurable boundary data and other immeasurable values and thus the reconstruction of the bioluminescent source is to identify the unknown vector from the photon flux density According to the uniqueness theorem, the BLT solution is not unique in the general case [21]. Some prior information or constraints such as permissible area of source should be imposed on the unknown variables to obtain a meaningful reconstruction result. Considering the source permissible region, we can obtain the linear relation between the photon flux density and the source energy density distribution in the light source permissible region, that is, where the coefficient matrix is ill-conditioned and can cause severe numerical instabilities in the solution. Therefore, it cannot be directly solved using a simple least squares method.

2.3. Regularization

In order to obtain a stable solution, regularization methods are typically used for solving inverse problems [8, 22, 23]. The commonly used Tikhonov regularization method approximately solves (5) by converting it into the following minimization problem: where is a properly chosen regularization parameter. As a function of the regularization parameter, the solution of (5) is given by However, the reconstructions with Tikhonov regularization method assume that the coefficient matrix is exactly known and noises only exist in the measurement. The regularization solutions computed by (7) do not take system errors into account.

As mentioned in the introduction, TLS method is designed for the case that both sides of the matrix equation are subject to errors. BLT inverse problem can be stated with TLS formulation as follows: where denotes the Frobenius norm, and are the error versions of and respectively, and is the augmented matrix that combines matrix and vector by using as the last column of the new matrix. Based on the TLS method, TTLS method is proposed by Fierro in [11] for regularization of ill-conditioned linear systems. In TTLS, the redundant information in associated with the small singular values, is discarded and the original ill-conditioned problem is replaced with another appropriate and more well-conditioned problem.

The TTLS algorithm used in this paper can be summarized as follows.(1)Compute the SVD of the augmented matrix where is m×(n+1), is and and is an diagonal matrix with the singular values on the diagonal.(2)Select a truncation parameter (3)Partition the matrix such that (4)Then the TTLS solution is given by

In fact, the aim of TTLS regularization is to appropriately identify an optimal truncation level, and then to construct a truncated solution that can capture the essential features of the unknown true solution, without explicit knowledge about the true solution and even without a priori knowledge about the magnitude of the noise in the data. For this purpose, truncation level must be carefully determined.

2.4. Choice of the Truncation Parameter

MGCV criterion proposed by Sima in [17] makes use of the filter factor formulation of the TTLS solution proved in [11]: where the filter factor values

The property used for choosing the truncation parameter is that when the parameter is greater than a certain crucial value, the TTLS solution is very sensitive to the noise or errors. Specifically, for small truncation level the filter factors with indices stay close to 1 and the filter factors with indices stay close to 0; when the truncation level gradually increases to a certain critical value, the filter factors with indices nearby increase dramatically. It implies a way to identify the value of where the filter factors change their steady behavior into erratic growth behavior.

As for the regularization problem in BLT, the choice of regularization parameter with classical GCV is by means of minimizing the GCV function: where presents the pseudoinverse of . With filter factors, the denominator can be computed by means of the following expression: where p is the rank of matrix with the singular values on the diagonal. We denote the sum of the filter factors of TTLS solution by as the effective number of parameters:

According to the properties of filter factors mentioned above, for a above a certain critical value, the filter factors for TTLS solutions with indices nearby are larger than 1. A fact can be derived that is greater than when reaches this critical value, which is used to modify the above classic GCV function to suit the TTLS case. And then the MGCV criterion for TTLS is obtained

However, the regularization parameter directly identified by (17) may be not optimal for the specific BLT reconstruction problem. Inspired by L-curve method, we propose a hybrid scheme that combines MGCV with the minimization of the corresponding residual norm for regularization parameter choice. The IGCV scheme for TTLS is summarized in the following steps.

Step 1. Use the MGCV criterion to get an initial truncation parameter and compute where is the maximum such that at the same time an array of GCV function values is obtained,

Step 2. For find the local minimum points that satisfy the conditions: and

Step 3. For all the local minimum points, compute the residual error where is an approximation of the TTLS solution for a given truncation level that is, only the top 70% of the nodal values are kept for computation convenience, and then the final truncation parameter

Thus, a proper truncation parameter for TTLS is sought according the above IGCV scheme.

3. Experiments and Results

The experiments implemented in this section are to test the performance of TTLS combined with IGCV for BLT inverse problem. To demonstrate the effectiveness of the proposed scheme, we compare the following reconstruction algorithms: Tikhonov method with classical GCV (Tik-GCV), TTLS method with MGCV (TTLS-M), and TTLS method with the proposed IGCV (TTLS-I). The parameter-choice scheme of Tikhonov method is different from that of TTLS method because MGCV and IGCV are specially designed for TTLS. A similar scheme, namely, classical GCV, is adopted in Tikhonov method for comparison convenience. The qualities of the reconstruction are assessed by the following quantitative indices: relative residual error (RRE), reconstructed location error, and reconstructed source power. Here, RRE is used to depict the extent of the solution fitting the measured data and is defined as Absolute error (AE) of the reconstructed source location is used to describe the accuracy of the reconstruction, which is defined by, where is the reconstructed center of each source and ( the actual center. Considering the ill-posedness of the BLT inverse problem, it is difficult to discriminate the influence of small source of high density and large one of low density [24]. So we prefer reconstructed source power compared with the actual value to source density for evaluating the quality of the reconstruction results. And the source power is estimated by computing the integral of the source intensity over its support [25].

3.1. Numerical Simulation Verifications

In the numerical simulation, a 30 mm diameter and 30 mm high cylindrical mouse chest phantom is designed to evaluate the performance of the reconstruction method. The structure of the phantom is shown in Figure 1(a). The phantom is heterogeneous and the corresponding optical parameters are set as in Table 1 [25]. Two sphere sources of 0.5 mm diameter with 1 energy density are located in the left lung and the centers are and respectively. The power of each source is 0.5236 nW. In the following single source case, only the source centered at is considered.


Material (mm-1) (mm-1)

Tissue0.0071.031
Lung0.0232.000
Heart0.0111.096
Bone0.0010.060

In order to reduce the ill-posedness of the inverse problem, a priori information of the source permissible region (PR) is incorporated to our method, which is shown in Figure 1(b) as [25], where is the coordinates of the corresponding FEM mesh vertices.

Generally speaking, simulated data used in reconstruction algorithms for inverse problems often come from the numerical solution of the forward problem. To avoid the typical issue of inverse crime, we use different FEM discretization for the forward process and reconstruction algorithms. Specifically, the forward model contains 11997 mesh vertices corresponding to 66334 tetrahedral elements, whereas the reconstruction model is consisting of 5277 vertices and 27465 tetrahedral elements. In addition, we employ Lagrange-Quadratic interpolation function in the forward process owing to the observation that high-order interpolation function can improve the numerical accuracy of the forward solution [26, 27].

To comprehensively simulate the noise and system errors involved in real BLT experiment, the photon flux density is added with Gaussian white noise, and the coefficient matrix A is added with a system errors matrix. Due to the complexity of error sources, it is difficult to have an exact mathematic model to describe the system errors accurately. Hence we adopted the commonly used Gaussian white noise [2830] and exponential noise to simulate the errors in matrix respectively.

As discussed in Section 2, regularization parameter is the crucial factor that affects the quality of regularization solution to inverse problem. Figure 2 illustrates the determination of regularization parameters in single source case with measurement noise level of 10% and Gaussian system error level of 1%. Among them, Figure 2(c) shows the residual error values of all the local minimum points described in our improved scheme IGCV, which are used for the selection of an optimal truncation parameter for TTLS. It should be noticed that the parameter k identified by MGCV is 64, whereas the optimal parameter obtained by IGCV is 78. It is because is 0.0038 and is 0.0046, which indicate that 64 is not the optimal parameter value according to IGCV criterion. The determination of regularization parameter in double sources case is similar to that of single source case. For space limitation, we just provide the final regularization parameter obtained in various noise settings in Tables 2 and 3.


Sys. errorMeas. noiseRecon. methodSys. error patternRegular. param.RRERecons. power (nW)

Without errorsWithout noisesTik-GCVN/A0.009460.02630.4284
TTLS-MN/A750.02430.4661
TTLS-IN/A780.02360.4797
10%Tik-GCVN/A0.014370.05850.4096
TTLS-MN/A690.05170.4535
TTLS-IN/A780.04190.5232
20%Tik-GCVN/A0.022660.10830.3825
TTLS-MN/A650.08090.4641
TTLS-IN/A690.07030.5214

1%10%Tik-GCVGaus.0.018470.05960.3502
Exp.0.022880.07050.3074
TTLS-MGaus.640.05900.3982
Exp.710.05620.3843
TTLS-IGaus.780.05240.4389
Exp.820.05570.4113
20%Tik-GCVGaus.0.028440.11500.4295
Exp.0.022810.08970.3432
TTLS-MGaus.520.11420.4410
Exp.730.05990.4755
TTLS-IGaus.640.09150.5427
Exp.870.05970.4829

5%10%Tik-GCVGaus.0.037710.08500.3069
Exp.0.050160.10230.2760
TTLS-MGaus.650.08230.3498
Exp.530.10180.2960
TTLS-IGaus.940.07490.3553
Exp.850.09780.2995
20%Tik-GCVGaus.0.054140.16260.2722
Exp.0.065860.18710.2618
TTLS-MGaus.520.15360.3088
Exp.450.17580.2927
TTLS-IGaus.610.14440.3359
Exp.1090.16730.3044


Sys. errorMeas. noiseRecon. methodSys. error patternRegular. param.RRERecon. power (nW)

Without errorsWithout noisesTik-GCVN/A0.008580.02520.46710.2677
TTLS-MN/A760.02280.48240.2911
TTLS-IN/A760.02280.48240.2911
10%Tik-GCVN/A0.015330.05860.43340.2148
TTLS-MN/A750.04180.50470.3560
TTLS-IN/A750.04180.50470.3560
20%Tik-GCVN/A0.030010.13630.36330.1787
TTLS-MN/A590.10740.44250.2069
TTLS-IN/A710.07390.50160.3004

1%10%Tik-GCVGaus.0.017110.05600.41640.2448
Exp.0.021480.06180.41740.1974
TTLS-MGaus.740.04500.45810.3670
Exp.710.05340.51700.2440
TTLS-IGaus.740.04500.45810.3670
Exp.710.05340.51700.2440
20%Tik-GCVGaus.0.031120.12480.35060.1891
Exp.0.030340.12130.36880.1862
TTLS-MGaus.590.10230.45290.2042
Exp.500.13470.40850.2086
TTLS-IGaus.720.07800.47680.2782
Exp.580.11590.42090.2253

5%10%Tik-GCVGaus.0.041790.09690.26130.2347
Exp.0.051540.11510.35280.2427
TTLS-MGaus.630.08900.31940.2647
Exp.500.11920.34950.2576
TTLS-IGaus.810.08800.33990.2660
Exp.650.11410.35270.2906
20%Tik-GCVGaus.0.050770.14670.28590.1574
Exp.0.067300.18580.31220.2480
TTLS-MGaus.510.13760.37130.1574
Exp.480.17410.35570.3304
TTLS-IGaus.870.13390.36920.2147
Exp.540.16980.35620.3306

In single source test, we found that all the methods under consideration can detect the source with the same center location in different noise levels, but the reconstructed source power varies with different reconstruction methods. Although the absolute error of the source location is the reconstructed source center is the nearest node to the original location in the aforementioned FEM discretization. Figure 3 only shows the reconstruction results by our proposed method with measurement noise level of 10% and Gaussian system error level of 1%. The detailed quantitative reconstruction results for the single-source model in various noise settings are listed in Table 2. The optimal results are listed in bold. Based on the simulation results in single-source case, it is clear that all the reconstruction methods can estimate the source location with no matter Gaussian or exponential noise in matrix but TTLS combined with IGCV performs best in all quantitative indices under different noise or error levels.

In the double sources case, both of the two sphere sources located in the left lung are tested. The final reconstruction results are listed in Table 3. Under all the noise conditions considered in this paper, the three methods can reconstruct the two sources at and = which are 0.911mm and 0.467mm away from the actual ones, respectively. In fact, they are the nearest nodes to the original source locations under the FEM mesh used in our tests. However, with the increase of noise or error level, besides the optimal nodes and some artifacts appear in the reconstruction results, which are illustrated in Figure 4. Simulation results in double sources case further show that although there are differences between the results of different noise pattern in matrix similar conclusions can be obtained. As shown in Table 3, the reconstruction results of TTLS combined with MGCV are comparable to that of TTLS combined with IGCV when noise level is low; whereas with the increase of noise or error, TTLS combined with IGCV outperforms the other methods in all quantitative indices.

For BLT inverse problem, permission region is an effective way to regularize the solution by restricting the source distribution within a proper permissible region. In order to further test the proposed method, a ball shape permissible region of 10 mm in diameter is utilized, which is expressed as . The sources settings in this section are the same as the aforementioned double sources case. The source distribution in the ball permission region was reconstructed, and the results are summarized in Table 4. Considering that the different system error pattern has little effect on the reconstruction results in the foregoing simulations, we only add Gaussian noise to the system matrix in this section.


Sys. errorMeas. noiseRecon. methodRegular. param.RRERecon. position (mm) and power (nW)

1%10%Tik-GCV0.027280.1720 0.2102 0.2509
TTLS-M190.1535 0.2308 0.2373
TTLS-I 300.1523 0.2978 0.3062
20%Tik-GCV0.093970.4383 0.3633 0.3501
TTLS-M90.4288 0.3782 0.3651
TTLS-I280.4226 0.5808 0.4632

5%10%Tik-GCV0.096600.2350 0.3693 0.3784
TTLS-M90.2350 0.3752 0.3853
TTLS-I180.2171 0.3899 0.4234
20%Tik-GCV0.128040.4367 0.3549 0.3578
TTLS-M150.3865 0.3900 0.3653
TTLS-I190.3532 0.4538 0.3807

It is shown in Table 4 that TTLS combined with IGCV still performs best under all the noise levels in terms of RRE, reconstructed power and source location. Compared with the results in Table 3, the location accuracy for ball shape permission region is lower. For example, the largest deviation of the reconstructed position of is up to 1.2 mm. It is clear that all reconstruction methods under consideration suffer from performance degradation with the relaxation of the permission region. However; the proposed method outperforms the other two methods and produces acceptable reconstruction results in our tests.

3.2. Physical Experiment Verifications

A physical experiment was carried out to further investigate the performance of the proposed method. A cylindrical phantom of 45 mm height and 22.5 mm radius was designed to evaluate different methods. The phantom shown in Figure 5(a) was made from nylon, and one small hole of 2.95 mm radius and 21 mm depth was drilled in the phantom to inject luminescent mixed solution used as the light source. In our physical experiment, the total volume of the mixed solution injected into the hole is 0.15 mL, thus a cylindrical source with a 2.95 mm radius and 5.4 mm height is centered at as shown in Figure 5(b). The optical parameters of the phantom were determined by a time-correlated single photon counting (TCSPC) system specifically constructed for the optical properties of the turbid medium [31]. The measured values of absorption and reduced scattering coefficients at the wavelength around 660 nm are 0.91 mm-1 and 0.0138 mm-1, respectively.

A scientific cooled back-illuminated CCD camera (PIXIS 2048B) is used to collect the outgoing photons from the phantom surface. The photon flux density from different angles can be acquired by rotating the stage under the phantom, as illustrated in Figure 5(c). Figures 6(a)6(d) exhibits the four views of the cylindrical phantom obtained by the CCD camera, respectively. Because the data captured by CCD camera is planar, mapping it onto 3D surface of the cylindrical phantom must be accomplished before reconstruction, which will also bring some inevitable errors to the measured data [32]. The mapping result was shown in Figure 6(e).

According to the photon flux density distribution on the phantom surface, the source permissible region is set as In the reconstruction process, the phantom model consists of 2734 vertices corresponding to13551 tetrahedral elements. The schemes for the selection of regularization parameters are identical to those in numerical simulations. The final reconstruction results and the corresponding regularization parameter are listed in Table 5. The 3D views of the reconstructed results using different methods are presented in Figure 7, which verified the feasibility and effectiveness of the proposed method. As is evident from the images in Figure 7 and the data in Table 5, TTLS combined with IGCV successfully reconstructed the luminescent source with the minimum distance of 1.76 mm away from the actual source center.


Recon. methodRegular. param.RRERecon. source position (mm)AE (mm)

Tik-GCV0.000030.9513(7.64,4.42,27.18)3.71
TTLS-M400.9015 2.97
TTLS-I450.8318(9.86,1.00,28.39)1.76

4. Discussion and Conclusion

BLT reconstruction is a highly ill-posed inverse problem where small measurement noise and system errors in the input data can produce large changes in the results. In addition, bioluminescence signals are generally very weak, thus the noise or errors will significantly affect the reconstruction quality. Regularization technique has played an important role in solving BLT inverse problem. And most of the previous works assume that there is only measurement noise, which affects the right-hand side of the system equations. However, the computed coefficient matrix in the model also has some errors, which may be caused by the calculation errors, the geometrical approximation, optical parameter inaccuracy, as well as the assumption of diffusion equation model itself. For example, the FEM discretization typically adds some errors to the matrix Hence, there is a need for seeking methods that can deal with the errors in both sides of the system equation. TTLS is a truncation regularization method that can take account of both system errors and measurement noise in the reconstruction process. This method depends on a parameter called truncation level; this single parameter has a significant influence on the regularization solutions. In this paper, IGCV, a practical scheme for determining the truncation parameter, is proposed to be combined with TTLS method for solving BLT inverse problem

Simulations considering both system errors and measurement noise are conducted to investigate the performance of the proposed reconstruction method. Due to the lack of an accurate model to describe the system errors arising from multiple sources, commonly used Gaussian white noise and exponential noise are adopted to simulate the errors in matrix respectively. In addition, physical phantom experiments further test the proposed method.

Both the numerical simulations and physical experiments demonstrate the effectiveness of the proposed method. Tests with different noise levels show that TTLS with combined IGCV is able to produce much better reconstruction results than Tikhonov method, and TTLS combined with IGCV performs better than TTLS combined with MGCV, especially when both sides of the system equation are contaminated by measurement noise and system errors. Based on the experiments in this paper, we can draw a preliminary conclusion that TTLS combined with IGCV criterion is a potential reconstruction method for BLT inverse problem. Further investigation of the performance of the proposed method on animal experiments will be conducted in our future work.

Acknowledgments

This work is supported by the Program of the National Basic Research and Development Program of China (973) under Grant no. 2006CB705700, the Cheung Kong Scholars and Innovative Research Team in University (PCSIRT) under Grant no. IRT0645, the Chair Professors of Cheung Kong Scholars Program of Ministry of Education of China, CAS Hundred Talents Program, the National Natural Science Foundation of China under Grant no. 30873462, 60532050, 30900334, the Beijing Municipal Natural Science Foundation of China under Grants no. 4071003, the CAS Scientific Research Equipment Develop Program (YZ0642, YZ200766), the Natural Science Basic Research Plan in Shaanxi Province of China under Grant no. 2009JQ8018, and the Science Foundation of Northwest University under Grant no. 09NW34.

References

  1. J. K. Willmann, N. van Bruggen, L. M. Dinkelborg, and S. S. Gambhir, “Molecular imaging in drug development,” Nature Reviews Drug Discovery, vol. 7, pp. 591–607, 2008. View at: Google Scholar
  2. V. Ntziachristos, J. Ripoll, L. V. Wang, and R. Weissleder, “Looking and listening to light: the evolution of whole-body photonic imaging,” Nature Biotechnology, vol. 23, no. 3, pp. 313–320, 2005. View at: Publisher Site | Google Scholar
  3. J. Tian, J. Bai, X.-P. Yan et al., “Multimodality molecular imaging: improving image quality,” IEEE Engineering in Medicine and Biology Magazine, vol. 27, no. 5, pp. 48–57, 2008. View at: Publisher Site | Google Scholar
  4. C. H. Contag and M. H. Bachmann, “Advances in in vivo bioluminescence imaging of gene expression,” Annual Review of Biomedical Engineering, vol. 4, pp. 235–260, 2002. View at: Publisher Site | Google Scholar
  5. M. Schweiger, S. R. Arridge, M. Hiraoka, and D. T. Delpy, “The finite element method for the propagation of light in scattering media: boundary and source conditions,” Medical Physics, vol. 22, no. 11, pp. 1779–1792, 1995. View at: Publisher Site | Google Scholar
  6. A. N. Tikhonov and V. Y. Aresenin, Solutions of Ill-Posed Problems, V. H. Winston & Sons, Washington, DC, USA, 1977.
  7. G. Wang, W. Cong, K. Durairaj et al., “In vivo mouse studies with bioluminescence tomography,” Optics Express, vol. 14, no. 17, pp. 7801–7809, 2006. View at: Publisher Site | Google Scholar
  8. Y. Lu, X. Zhang, A. Douraghy et al., “Source reconstruction for spectrally-resolved bioluminescence tomography with sparse A priori information,” Optics Express, vol. 17, no. 10, pp. 8062–8080, 2009. View at: Publisher Site | Google Scholar
  9. S. V. Huffel and J. Vandewalle, The Total Least Squares Problem: Computational Aspects and Analysis, Frontiers in Applied Mathematics, Society for Industrial Mathematics, Philadelphia, Pa, USA, 1991.
  10. G. H. Golub and C. F. Van Loan, “An analysis of the total least squares problem,” SIAM Journal on Numerical Analysis, vol. 17, no. 6, pp. 883–893, 1997. View at: Google Scholar
  11. R. D. Fierro, G. H. Golub, P. C. Hansen, and D. P. O'Leary, “Regularization by truncated total least squares,” SIAM Journal of Scientific Computing, vol. 18, no. 4, pp. 1223–1241, 1997. View at: Google Scholar
  12. P. C. Hansen, “Truncated singular value decomposition solutions to discrete ill-posed problems with ill-determined numerical rank,” SIAM Journal on Scientific and Statistical Computing, vol. 11, no. 3, pp. 503–518, 1990. View at: Google Scholar
  13. P. C. Hansen, “Analysis of discrete ill-posed problems by means of the L-curve,” SIAM Review, vol. 34, no. 4, pp. 561–580, 1992. View at: Google Scholar
  14. P. C. Hansen, M. E. Kilmer, and R. H. Kjeldsen, “Exploiting residual information in the parameter choice for discrete ill-posed problems,” BIT Numerical Mathematics, vol. 46, no. 1, pp. 41–59, 2006. View at: Publisher Site | Google Scholar
  15. G. H. Golub, M. Heath, and G. Wahba, “Generalized cross-validation as a method for choosing a good ridge parameter,” Technometrics, vol. 21, pp. 215–223, 1979. View at: Google Scholar
  16. J. E. Moody, “The effective number of parameters: an analysis of generalization and regularization in nonlinear learning systems,” in Advances in Neural Information Processing Systems, Moody, Hanson, and Lippmann, Eds., vol. 4, pp. 847–854, Morgan Kaufmann, Palo Alto, Calif, USA, 1992. View at: Google Scholar
  17. D. M. Sima and S. V. Huffel, “Level choice in truncated total least squares,” Computational Statistics and Data Analysis, vol. 52, no. 2, pp. 1103–1118, 2007. View at: Publisher Site | Google Scholar
  18. W. Cong, G. Wang, D. Kumar et al., “Practical reconstruction method for bioluminescence tomography,” Optics Express, vol. 13, no. 18, pp. 6756–6771, 2005. View at: Publisher Site | Google Scholar
  19. X. Song, D. Wang, N. Chen, J. Bai, and H. Wang, “Reconstruction for free-space fluorescence tomography using a novel hybrid adaptive finite element algorithm,” Optics Express, vol. 15, no. 26, pp. 18300–18317, 2007. View at: Publisher Site | Google Scholar
  20. X. Song, J. Yi, and J. Bai, “A parallel reconstruction scheme in fluorescence tomography basedon contrast of independent inversed absorption properties,” International Journal of Biomedical Imaging, vol. 2006, no. 26, Article ID 70839, 7 pages, 2006. View at: Google Scholar
  21. G. Wang, Y. Li, and M. Jiang, “Uniqueness theorems in bioluminescence tomography,” Medical Physics, vol. 31, no. 8, pp. 2289–2299, 2004. View at: Publisher Site | Google Scholar
  22. A. X. Cong and G. Wang, “Multispectral bioluminescence tomography: methodology and simulation,” International Journal of Biomedical Imaging, vol. 2006, Article ID 57614, 7 pages, 2006. View at: Publisher Site | Google Scholar
  23. X. He, J. Tian, Y. Wu, Y. Hou, N. Ren, and K. Penga, “Study of four regularization methods for the inverse problem in bioluminescence tomography,” in Medical Imaging 2009: Biomedical Applications in Molecular, Structural, and Functional Imaging, vol. 7262 of Proceedings of SPIE, p. 10, Lake Buena Vista, Fla, USA, 2009, 7262B. View at: Publisher Site | Google Scholar
  24. H. Dehghani, S. C. Davis, and B. W. Pogue, “Spectrally resolved bioluminescence tomography using the reciprocity approach,” Medical Physics, vol. 35, no. 11, pp. 4863–4871, 2008. View at: Publisher Site | Google Scholar
  25. M. Jiang, T. Zhou, J. Cheng, W. Cong, and G. Wang, “Image reconstruction for bioluminescence tomography from partial measurement,” Optics Express, vol. 15, no. 18, pp. 11095–11116, 2007. View at: Publisher Site | Google Scholar
  26. Y. Hou, J. Tian, Y. Wu, and J. Liang, “A new numerical method for BLT forward problem based on high-order finite elements,” Communications in Numerical Methods in Engineering, vol. 25, no. 6, pp. 667–681, 2009. View at: Publisher Site | Google Scholar
  27. R. Han, J. Liang, X. Qu et al., “A source reconstruction algorithm based on adaptive hp-FEM for bioluminescence tomography,” Optics Express, vol. 15, no. 18, pp. 11095–11116, 2009. View at: Google Scholar
  28. W. Zhu, Y. Wang, Y. Yao, J. Chang, H. L. Graber, and R. L. Harbour, “Iterative total least-squares image reconstruction algorithm for optical tomography by the conjugate gradient method,” Journal of the Optical Society of America A, vol. 14, no. 4, pp. 799–807, 1997. View at: Google Scholar
  29. W. Zhu, Y. Wang, and J. Zhang, “Total least-squares reconstruction with wavelets for optical tomography,” Journal of the Optical Society of America A, vol. 15, no. 10, pp. 2639–2650, 1998. View at: Google Scholar
  30. G. Shou, L. Xia, M. Jiang, Q. Wei, F. Liu, and S. Crozier, “Truncated total least squares: a new regularization method for the solution of ECG inverse problems,” IEEE Transactions on Biomedical Engineering, vol. 55, no. 4, pp. 1327–1335, 2008. View at: Publisher Site | Google Scholar
  31. D. Qin, H. Zhao, Y. Tanikawa, and F. Gao, “Experimental determination of optical properties in turbid medium by TCSPC technique,” in Progress in Biomedical Optics and Imaging, vol. 6434 of Proceedings of SPIE, p. 10, San Jose, Calif, USA, 2007, 64342E. View at: Publisher Site | Google Scholar
  32. J. Ripoll, D. Yessayan, G. Zacharakis, and V. Ntziachristos, “Experimental determination of photon propagation in highly absorbing and scattering media,” Journal of the Optical Society of America A, vol. 22, no. 3, pp. 546–551, 2005. View at: Publisher Site | Google Scholar

Copyright © 2010 Xiaowei He et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views1609
Downloads526
Citations

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.