## Abstract

The purpose of this study is to propose a strategy for organ reconstruction in fluorescence molecular tomography (FMT) without prior information from other imaging modalities, and to overcome the high cost and ionizing radiation caused by the traditional *structural prior* strategy. The proposed strategy is designed as an iterative architecture to solve the inverse problem of FMT. In each iteration, a short time Fourier transform (STFT) based algorithm is used to extract the *self-prior* information in the space-frequency energy spectrum with the assumption that the regions with higher fluorescence concentration have larger energy intensity, then the cost function of the inverse problem is modified by the self-prior information, and lastly an iterative Laplacian regularization algorithm is conducted to solve the updated inverse problem and obtains the reconstruction results. Simulations and *in vivo* experiments on liver reconstruction are carried out to test the performance of the self-prior strategy on organ reconstruction. The organ reconstruction results obtained by the proposed self-prior strategy are closer to the ground truth than those obtained by the iterative Tikhonov regularization (ITKR) method (traditional *non-prior* strategy). Significant improvements are shown in the evaluation indexes of relative locational error (RLE), relative error (RE) and contrast-to-noise ratio (CNR). The self-prior strategy improves the organ reconstruction results compared with the non-prior strategy and also overcomes the shortcomings of the traditional structural prior strategy. Various applications such as metabolic imaging and pharmacokinetic study can be aided by this strategy.

© 2017 Optical Society of America

## 1. Introduction

Fluorescence molecular tomography (FMT) visualizes the physiological and pathological process in biological tissues at cellular and molecular levels, through noninvasively monitoring the three-dimensional (3-D) distribution of fluorescent probes in the tissues *in vivo*. Fluorescent probes are emitted by an excitation source and the fluorescence signals are detected in multiple angles. Spatial and quantitative distribution of fluorescent probes in the tissues can be obtained by tomographic reconstruction.

Tumor detection and organic functional imaging are two typical applications of FMT [1–5]. In the studies of tumor detection, early-stage tumors are generally selected as the studied targets and labeled with specific fluorescent probes. In this situation the imaged targets are usually much smaller than the background in the image domain. The reconstruction task can be treated as a sparse problem. A series of L1-based reconstruction algorithms such as greedy reconstruction, restarted nonlinear conjugate gradient based L1 regularization, adaptive support driven reweighted L1-minimization and discrete cosine transform-based regularization are proposed to solve this problem [6–10].

In the situations of organ functional imaging, the assumption of sparsity of the imaging targets becomes invalid, due to the large-scale accumulation of fluorescent probes in the organ. The non-sparse and complex distribution of fluorescent probes also introduce the mixture of fluorescent signals from the observed organ with signals from non-targeted regions and/or noises. The reconstruction quality suffers from the above issues, and therefore the accuracy and visibility of the observed organ are decreased. Generally, structural priors obtained from other imaging modalities are applied to constrain the reconstruction of the investigated organs. Structural prior information from magnetic resonance imaging (MRI) or X-ray computed tomography (XCT) has been used to constrain the reconstruction of organ or tumor [11, 12] and to assist the localization of the metabolic region in small animals [4]. Various multi-modality imaging systems such as computed tomography (CT)/FMT system, MRI/FMT system and positron emission computed tomography (PET)/FMT system, can provide different kinds of prior information for FMT reconstruction [13–15]. Although the structural prior strategy is currently the mainstream solution and can provide the fluorescence reconstruction results with significantly higher quality than other reconstruction methods, the introduction of other imaging modalities increases the complexity and cost of the imaging system. More noise interference and ionizing radiation also decrease the stability and safety of FMT.

In this study, a self-prior reconstruction strategy is proposed to realize the reconstruction of organ in FMT. Unlike the structural prior methods, this strategy only uses the data collected by FMT. The proposed strategy is designed as an iterative architecture and in each iteration, the self-prior information is extracted from the reconstruction results obtained at the previous iteration. The reconstructed distribution of fluorescent probes in the observed organ is improved significantly by updating the self-prior information and cost function of the inverse problem. Therefore, the self-prior strategy can be applied to the reconstruction of organ in the functional imaging, and also metabolic analysis such as dynamic fluorescence molecular tomography (DFMT).

The paper is organized as follows. Section 2 describes the details of the proposed strategy. Sections 3 and 4 describe the setups of simulations and *in vivo* experiments. Section 5 shows the reconstruction results and quantitative analysis. The discussion and conclusion are presented in Section 6.

## 2. Methods

#### 2.1 Forward model

In FMT, within the near-infrared spectral window, the photon migration through highly scattering media can be modeled using the coupled diffusion equations (DE), which are known as the low-order approximation of the radiative transfer equation (RTE), as follows [16]

Here the subscripts *e* and *m* denote the excitation and emission, respectively. *D* is the diffusion coefficient given by $D=1/\left(3\left({\mu}_{a}+{\mu}_{s}^{\text{'}}\right)\right)$, where ${\mu}_{a}$ is the absorption coefficient and ${\mu}_{s}^{\text{'}}$ is the reduced scattering coefficient. $\Phi $ is the optical flux density and $S\left(r\right)$ is the excitation source term. $\eta $ and ${\mu}_{af}\left(r\right)$ are the quantum yield and absorption coefficient of the fluorescent probe. The product of the two variables is the fluorescence yield $x\left(r\right)$ which reflects the distribution of fluorescent probes in the media.

If the optical properties are known, the fluorescence measurement $\Phi \left({r}_{d},{r}_{s}\right)$ can be obtained according to the first-order Born-type approximation [17].

*V*is the image domain, ${G}_{e}\left(r,{r}_{s}\right)$ and ${G}_{m}\left({r}_{d},r\right)$ denote the functions describing the photon propagation from the source position to an arbitrary position and from the arbitrary position to the detector position, respectively. $\Theta $ is a calibration parameter depending on the system characteristics.

The finite element method is used to solve ${G}_{e}\left(r,{r}_{s}\right)$ and ${G}_{m}\left({r}_{d},r\right)$ [18]. After discretization the forward model can be formulated as a linear matrix equation:

where ${W}_{\left(M\times N\right)}$ is a weight matrix and each element represents the contribution of each voxel in the discrete volume of interest to the fluorescence measurement ${\Phi}_{m}{}_{\left(M\times 1\right)}$. Vector ${x}_{\left(N\times 1\right)}$ is the fluorescence yield to be reconstructed, which represents the fluorescence distribution.*M*and

*N*stand for the numbers of discrete voxels in the detector domain and image domain, respectively.

#### 2.2 Tikhonov regularization

The goal of the inverse problem of FMT is to obtain the fluorescence yield ${x}_{\left(N\times 1\right)}$ by solving Eq. (3). To overcome the ill-posed nature of the inverse problem, Eq. (3) can be modified by employing the Tikhonov regularization method, which is a reconstruction strategy without any prior information. The problem is transferred into an optimization problem presented as [19, 20]

*k_itkr*denotes the index of iteration and ${W}^{H}{}_{\left(N\times M\right)}$ is the conjugate transpose matrix of ${W}_{\left(M\times N\right)}$.

#### 2.3 Self-prior strategy and Laplacian regularization

The self-prior strategy aims at extracting the locational information of fluorescence yield from the initial reconstruction results as a self-prior information, and then modifying the inverse problem of the FMT. The strategy is designed as an iterative architecture and each iteration includes three steps: extraction of the self-prior information, modification of the cost function of the inverse problem and reconstruction of the distribution of fluorescence yield. The fluorescence distribution ${x}_{0\left(N\times 1\right)}$ reconstructed by the ITKR algorithm is used as the initial input of the iterative architecture. Figure 1 shows the flowchart of the iterative architecture of the self-prior strategy.

### 2.3.1 Extraction

The self-prior spatial information is extracted from the reconstructed fluorescence distribution ${x}_{t\left(N\times 1\right)}$ (*t* stands for the index of iteration in the self-prior strategy; for the first iteration, the fluorescence distribution is *x*_{0}). In general, the spatial distribution of fluorescence tends to have large variations, leading to the difficulties in extracting the self-prior information directly from the spatial domain with threshold based methods. Because the self-prior information is more pronounced and easier to clarify in the spatial-frequency domain than in the spatial domain due to the energy accumulation phenomenon, the self-prior information will be extracted in the spatial-frequency domain in this study, in order to obtain the positions of self-prior information in both the spatial and frequency domains. In detail, a short time Fourier transform (STFT) is performed on the one-dimensional (1-D) spatial sequence *x _{t}* to obtain the space-frequency spectrum,

*H*stands for the width of the window function.

*H*can influence both the spatial and frequency resolutions. In order to extract relatively accurate information from the frequency domain, the value of

*H*is set to be

*N*/1000 after multiple tests with different widths of window function, where

*N*denotes the length of the 1-D spatial sequence

*x*. The space-frequency energy spectrum ${F}_{t}\left(n,{\omega}_{k}\right)$ of

_{t}*x*can be presented as

_{t}By assuming that the regions with fluorescence distribution contain higher energy than the background regions, the self-prior information can be extracted in the space-frequency energy domain by selecting the spectral components with relatively high energy. An energy threshold ${T}_{F}$ is selected and the components higher than this threshold are chosen to be the self-prior information of fluorescence distribution, and the regions corresponding to these components are defined as the self-prior regions. The discrete voxels within the self-prior regions are located in the space-frequency energy spectrum, corresponding to the 1-D spatial sequence. A self-prior vector ${K}_{t}{}_{\left(N\times 1\right)}$ corresponding to the current reconstructed *x _{t}* is then formed as follows:

### 2.3.2 Modification

The cost function of the inverse problem is modified using the self-prior information obtained from the previous iteration. In detail, a Laplacian-type constraint matrix ${L}_{t\left(N\times N\right)}$ is generated based on the 1-D self-prior vector ${K}_{t\left(N\times 1\right)}$:

The cost function of the inverse problem is updated by replacing the identity matrix with ${L}_{t\left(N\times N\right)}$ in the penalty item and the new reconstructed fluorescence distribution *x _{t + 1}* can be updated by solving:

### 2.3.3 Reconstruction

After the modification of the inverse problem, the updated cost function can be optimized by the iterative algorithm as follows:

*k_self*denotes the index of iteration in the optimization process. The reconstructed fluorescence distribution ${x}_{t+1}$ is updated by the solution of Eq. (11) after

*k_self*iterations and will be imported into the next iteration until the iterative process is terminated.

## 3. Simulations

Numerical simulations are performed to evaluate the performance of the proposed strategy. Digimouse, a 3-D mouse atlas, is used to construct a 3-D biological simulation model [22]. Figure 2(a) shows the anatomic structure of the atlas, and four kinds of organs including the heart, lungs, liver and kidneys are contained in this model. As shown in Table 1, different optical properties are assigned to these organs to constitute a heterogeneous model.

Indocyanine green (ICG) is chosen as the fluorescent probe. The metabolism of ICG in realistic environment is a dynamic process and different organs have different metabolic characteristics. Figure 2(b) shows the ICG concentration curves which simulate the metabolic process of ICG in different organs. In the simulations, the distribution of ICG concentration is uniform in a single organ at a specific time point. The ICG is injected at the time of 0 minute. The liver is selected to be the investigated organ and the distribution of ICG at the time of 9th minute is chosen to generate the fluorescent signal used for reconstruction because it reaches the peak value at this time point, according to Fig. 2(b). The detailed ICG concentrations of each organ at this time point are also listed in Table 2.

The mouse torso from the neck to the base of the abdomen is chosen as the image domain, with a length of 3.2 cm. 24 fluorescence projections are generated equidistantly in a full view of 360° with an angular increment of 15°. The rotation axis of the mouse is defined as the z axis with the bottom plane set as z = 0.0 cm. In the reconstruction, the image domain is discretized into 70 × 42 × 32 mesh voxels using COMSOL Multiphysics 3.5 (COMSOL Inc., Stockholm, Sweden) is used for the discretization process. Since the structural information is not obtainable for the reconstruction process, the discretization of the image domain is homogeneous and not constrained by any prior information such as tissue boundary.

Three sets of simulations are conducted to test the influence of three different factors, i.e., the energy threshold, number of iterations and noise level, on the proposed strategy. The studies on the influence of the energy threshold and number of iterations, which are two important parameters in the self-prior reconstruction process, provide a set of reference values for the studies on the influence of different noise levels.

## 4. *In vivo* experiments

*In vivo* animal experiments are conducted to further test the performance of the proposed strategy, under realistic physiological environment. In the *in vivo* experiments, the liver is selected as the objective organ to be reconstructed, which is consistent with the simulations. The data acquisition is based on a hybrid FMT/XCT system previously developed by our laboratory [13]. A Xenon lamp (300 W, MAX-302, Asahi Spectra, Torrance, CA, USA) is employed as the excitation source. A fiber is attached to the lamp to generate a line-shaped excitation source with a length of 4 cm. A 780 ± 6 nm band-pass filter (XBPA780, Asahi Spectra, Torrance, CA, USA) is used to select the excitation light from the fiber. The emitted fluorescence is filtered by an 840 ± 6 nm band-pass filter (FF01-840/12-25, Semrock, Rochester, NY, USA) and detected by a 512 × 512 pixel, −70 °C cooled charge coupled device (CCD) (iXon DU-897, Andor Technologies, Belfast, Northern Ireland, UK). The XCT module is used for acquiring the XCT images to make a comparison with the FMT reconstruction results.

All animal studies are conducted under the protocol approved by the Ethical Committee of Tsinghua University. A female nude mouse (BALB/c, 8 weeks old) is used in the experiments. A bolus of ICG (0.1 mL, 50 μg/mL) is injected via the tail vain. With similar setups to the simulations, the mouse is fixed on the rotation stage, and the fluorescence image frames are collected by the CCD with an angular increment of 15°. The fluorescence data collection lasts for 20 minutes and every 24 frames are arranged as a single data set, corresponding to a full view of 360°. Similar to the simulation, the data set at the time point when the corresponding fluorescent signal of the liver reaches the highest intensity, is extracted for reconstruction. After the fluorescence data acquisition, the XCT signals are collected by a complementary metal oxide semiconductor flat-panel detector (C7921-02, Hamamatsu, Japan).

A set of 62 × 62 × 26 mesh voxels are generated to discretize the image domain, without any constraint of structural prior information. All the parameters of the self-prior reconstruction strategy are set to be the same as those in the simulations. Figure 3 shows the white light image of the nude mouse fixed on the rotation stage and its meshed volume of the image domain.

## 5. Results

In this section, the reconstruction results of the simulations and *in vivo* experiments are presented. In order to compare different reconstruction methods and evaluate the performance of the proposed strategy, the results of the non-prior strategy (i.e., the ITKR method) and structural prior strategy are also presented.

#### 5.1 Evaluation indices

The RLE demonstrates the reconstruction accuracy of spatial structure (or locational accuracy). The normalized reconstruction result (${x}_{recon}$) is binarized to obtain the locational reconstruction result (${x}_{loc-recon}$), as follows,

*i*stands for the

*i*th element of ${x}_{loc-recon}$. The true locational result (${x}_{loc-true}$) of the objective region (${R}_{obj}$) can be described as

The RE evaluates the reconstruction accuracy of the fluorescence yields. It is defined as the relative difference between the normalized reconstruction result (${x}_{recon}$) and the true distribution of the fluorescence yields (${x}_{true}$), as follows,

The CNR describes how much the reconstructed fluorescence distribution is recovered in the objective region or lost in the background, and it is defined as

_{VOI}= S

_{VOI}/ (S

_{VOI}+ S

_{BG}) and ω

_{BG}= S

_{BG}/ (S

_{VOI}+ S

_{BG}) are the weights of the VOI and BG, where S

_{VOI}and S

_{BG}are the volumes of the VOI and BG, respectively. ${\sigma}_{VOI}^{\text{2}}$ and ${\sigma}_{BG}^{\text{2}}$ are the variances of the reconstructed fluorescence yields in the VOI and BG, respectively. In this work, the region which contains the main distribution of the fluorescent yields is selected as the volume of interest, and then the rest of the image domain is set as the background.

#### 5.2 Results of simulations

### 5.2.1 Influence of energy threshold

The energy thresholds are set to be 25%, 45%, 65% and 85% of the maximum value in the space-frequency energy spectrum of the reconstruction results in each iteration, respectively. To eliminate the influence of other factors, the number of iterations is set to be 5 and zero-mean Gaussian noise with the signal-to-noise ratio (SNR) of 26 dB is added to the fluorescence data.

Figure 4 shows the reconstruction results at different energy thresholds in the simulations. The energy threshold significantly influences the reconstruction results. A higher energy threshold (i.e., 65% and 85%) will generate a relatively smaller self-prior region and may lead to the insufficient reconstruction results. On the contrary, a lower threshold (i.e., 25%) results in more non-target or background information in the self-prior region and the reconstruction results are more likely to contain the true distribution of fluorescent probes in the liver. Table 3 shows the results of REs. The energy threshold of 45% obtains the best reconstruction accuracy with the minimum RE, which is in agreement of the qualitative results in Fig. 4.

### 5.2.2 Influence of number of iterations

The number of iterations is set to vary from 1 to 10 for a single reconstruction task, with an SNR of 26 dB and an energy threshold of 45% of the maximum energy in each iteration.

Figure 5 presents the variations of REs and CNRs along with the number of iterations. The REs decrease dramatically after the second iteration and then reach a plateau afterwards. The CNRs show an opposite tendency of variation. This phenomenon demonstrates that the second iteration shrinks the reconstruction region and the following iterations search for the final results gradually in this region. Because no self-prior information can be used in the first iteration (i.e. the ITKR method), insufficient performance is obtained, with high RE and low CNR. After the self-prior information is extracted iteratively, the reconstruction results turn to be much closer to the ground truth.

### 5.2.3 Influence of noise level

Three sets of simulations with additive zero-mean Gaussian noise are conducted, with an SNR of 26 dB, 20 dB and 14 dB, respectively. According to the reference values provided by the previous studies in 1) and 2), in the self-prior reconstruction strategy, the number of iterations is configured to be 5. In each iteration, the energy threshold in the extraction step is set to be 45% of the maximum value in the space-frequency energy spectrum, and 100 iterations are executed in the optimization process.

Figure 6 shows the reconstruction results at different noise levels in the simulations. Figure 6(a) shows the true distribution of ICG in the liver in the 3-D and transverse views. It can be seen that the distribution of ICG is uniform in the liver. Figure 6(b) shows the space-frequency energy spectrum of the reconstruction result of the non-prior strategy with an SNR of 26 dB. Figures 6(c)-6(k) display the reconstruction results by the non-prior, self-prior and structural prior strategies at different noise levels, respectively. The reconstruction quality of the self-prior strategy is observed to be better than that of the non-prior strategy at all noise levels and is closer to that of the structural prior strategy on the aspects of the size, shape and location of the reconstructed region, by comparing the reconstruction results with the liver in Fig. 6(a). As the noise level increases, the self-prior strategy also shows greater noise resistance than the non-prior strategy.

Quantitative evaluation of different reconstruction strategies is shown in Table 4. Compared with the non-prior strategy at SNR levels of 26 dB, 20 dB and 14 dB, the self-prior strategy reduces the RLEs by 70.9%, 70.0% and 78.4%, reduces the REs by 67.6%, 72.8% and 67.8%, and increases the CNRs by 322.0%, 416.1% and 479.1%, respectively. This means that the proposed strategy obtains higher reconstruction accuracy and improves the detectability of the fluorescent probes. The RLEs and REs of the self-prior strategy maintain lower value than the non-prior strategy when the SNR decreases from 26 dB to 14 dB, which indicates that the self-prior strategy obtains higher reconstruction accuracy without sacrificing the stability when the noise level increases. Additionally, the locational accuracies (calculated by the RLEs) of the self-prior strategy reach 92.6%, 89.4% and 88.0% of those of the structural prior strategy with SNRs of 26 dB, 20 dB and 14 dB, respectively. In contrast, the RLEs of the non-prior strategy are 74.3%, 63.7% and 42.0% of those of the structural prior strategy, respectively. Similar results are found for the REs and CNRs. These results indicate that the performance of the self-prior strategy is much better than that of the non-prior strategy and it approaches the performance of the structural prior strategy in the simulations.

#### 5.3 Results of in vivo experiments

Figure 7 presents the tomographic reconstruction results of different methods in the *in vivo* experiments. Figure 7(g) shows the XCT image in the coronal view and Figs. 7(a)-7(f) shows the tomographic images of different methods in the transverse view at the position indicated by the red dotted line shown in Fig. 7(g). The FMT result of the self-prior strategy has higher quality and more uniform distribution of fluorescence yield than the non-prior strategy, as can be seen in Figs. 7(a) and 7(d). The FMT images are merged with XCT images to present the locational performance of each method. Figs. (c) and (f) show the merged images of FMT/XCT reconstructed by the non-prior strategy and the self-prior strategy, respectively. The regions of the liver are covered by the reconstructed fluorescent signals in both images, which means that the accumulation of ICG in the liver is detected and located by both methods. However, the self-prior strategy shows higher locational quality of the liver in the merged images according to the coverage of the FMT results and XCT results.

Figure 8 shows the reconstruction results of the *in vivo* experiments in the 3-D view. Figures 8(a) and 8(c) shows the FMT results in the 3-D view. Then they are merged with XCT results in the 3-D view to exhibit the relative relation between reconstructed fluorescent signals and anatomic structures, as shown in Figs. 8(b) and 8(d). Figs. 8(e) and 8(f) are the XCT image as the reference for location in the transverse and 3-D views, respectively. The self-prior strategy obtains higher accuracy in locating the liver and higher smoothness of the fluorescence yield within the liver than the non-prior strategy, as shown in Fig. 8. Quantitative evaluation indexes of the RLE and CNR are demonstrated in Table 5. The REs are not calculated because the true distribution of fluorescence yield in the mouse cannot be obtained. The XCT results are set to be the ground truth in the calculation of the RLEs. Compared with the non-prior strategy, the self-prior strategy reduces the RLE by 81.1% and increases the CNR by 952.6%.

## 6. Discussion and conclusion

Organ reconstruction is an intractable task when applying FMT to organ functional imaging and metabolic research. The large-scale accumulation of fluorescent probes in the organ makes the problem non-sparse and complex. Structural prior strategies based reconstruction methods have been developed, but the introduction of the anatomical imaging modalities increase the cost and system complexity. Since the structural prior strategies considerably rely on the information provided by the introduced anatomical imaging modalities, its performance would directly influence the accuracy and stability of the reconstruction results. Additionally, the ionizing radiation from the anatomical imaging modalities may decrease the safety of FMT.

In this study, a self-prior reconstruction strategy is proposed to realize the organ reconstruction in FMT. The strategy is designed as an iterative architecture and the self-prior information is extracted to constrain the inverse problem solving. In the simulations conducted on the Digimouse model, the distribution of fluorescent probes in the liver is reconstructed by the proposed strategy. The effects of three influence factors (i.e., the energy threshold, number of iterations and noise levels) on the performance of the proposed strategy are investigated. The energy threshold mainly influences the locational accuracy of the reconstruction, because a small or large threshold leads to insufficient or excessive self-prior region. The number of iterations influences both the locational accuracy and CNR of the reconstruction results. The lack of the self-prior information causes significant deviations in the first iteration and the reconstruction results approach the ground truth gradually as the number of iterations increases.

To investigate the influence of the noise level, comparison is conducted among the non-prior strategy, self-prior strategy and structural prior strategy. The results show that the self-prior strategy reduces the REs by 70.0-78.4%, reduces the RLEs by 67.6%-72.8% and increases the CNRs by 322.0-479.1% compared with the non-prior strategy. Additionally, the self-prior strategy can reach 88.0-92.6% of the locational accuracy evaluated by the RLE, 79.0-90.9% of the reconstruction accuracy of fluorescence yield evaluated by the RE and 73.5-85.9% of the CNR of the structural prior strategy. However, the corresponding values are 42.0-74.3% (RLE), 32.6-71.4% (RE) and 12.7%-19.8% (CNR) for the non-prior strategy. *In vivo* experiments are also conducted to further test the performance of the methods. Similar results are found that the self-prior strategy achieves higher performance than the traditional non-prior strategy in the liver reconstruction of a nude mouse. By using the self-prior strategy, the RLE is reduced by 81.1% while the CNR is increased by 952.6%, compared with the non-prior strategy. The proposed strategy is demonstrated to significantly improve the reconstruction performance without the need of the prior information provided by other imaging modalities. By using the proposed strategy, the size, shape and location of the reconstructed organ is observed to approach the ground truth. However, the reconstruction results also show some disparities between the proposed strategy and the structural prior strategy, especially at high level of noise.

The extracted self-prior regions are affected significantly by the energy threshold. Selection of this parameter is a non-trivial task. In this study, the energy threshold is determined empirically. In addition, a single threshold used for the whole domain may also affect the accuracy of the reconstruction results. In the situation of multiple reconstruction targets, it will be difficult to classify all the targets using a single threshold. In the future, an adaptive method will be studied to select the energy threshold automatically.

In the reconstruction step of the self-prior iterative architecture, only the regularized item (i.e., $\lambda {\Vert Lx\Vert}_{2}^{2}$) is updated by reassembling the Laplacian matrix according to the self-prior information, while the forward model remains the same in each iteration. In fact, the biological tissue is assumed to be homogeneous in the construction of the forward model, leading to an inaccurate cost function of the inverse problem. In order to solve this issue, a heterogeneous forward model can be constructed by modifying the optical properties of the updated self-prior regions in the iterative process.

From the methodological aspect, more reliable algorithms will be developed for accurate and stable self-prior information extraction. For biomedical applications, the proposed strategy can be combined with DFMT and multispectral FMT [5, 23, 24] to achieve dynamic and real-time metabolic imaging of organs *in vivo*. In the present study, only liver has been selected as the objective organ to be reconstructed, which has limited utility in terms of biological applications. To generalize this method into more scenarios of organ imaging, tests on different organs should be conducted in the future. The technique of FMT can also be used in pharmacokinetic study and drug development, without the need of other imaging modalities.

## Disclosures

The authors declare that there are no conflicts of interest related to this article.

## Funding

This work is supported by the National Natural Science Foundation of China under Grant Nos. 81227901, 81271617, 61322101, and 61361160418 and the National Major Scientific Instrument and Equipment Development Project under Grant No. 2011YQ030114.

## Acknowledgments

The authors would like to thank Dr. Fei Liu in their lab for her help on experiments. The authors declare that there are no conflicts of interest related to this article.

## References and links

**1. **N. C. Deliolanis and V. Ntziachristos, “Fluorescence molecular tomography of brain tumors in mice,” Cold Spring Harb. Protoc. **2013**(5), 438–443 (2013). [CrossRef] [PubMed]

**2. **Y. Zhang, B. Zhang, F. Liu, J. Luo, and J. Bai, “In vivo tomographic imaging with fluorescence and MRI using tumor-targeted dual-labeled nanoparticles,” Int. J. Nanomedicine **9**, 33–41 (2014). [PubMed]

**3. **F. Liu, X. Cao, W. He, J. Song, Z. Dai, B. Zhang, J. Luo, Y. Li, and J. Bai, “Monitoring of tumor response to cisplatin by subsurface fluorescence molecular tomography,” J. Biomed. Opt. **17**(4), 040504 (2012). [CrossRef] [PubMed]

**4. **G. Zhang, F. Liu, B. Zhang, Y. He, J. Luo, and J. Bai, “Imaging of pharmacokinetic rates of indocyanine green in mouse liver with a hybrid fluorescence molecular tomography/x-ray computed tomography system,” J. Biomed. Opt. **18**(4), 040505 (2013). [CrossRef] [PubMed]

**5. **G. Zhang, F. Liu, H. Pu, W. He, J. Luo, and J. Bai, “A Direct Method With Structural Priors for Imaging Pharmacokinetic Parameters in Dynamic Fluorescence Molecular Tomography,” IEEE Trans. Biomed. Eng. **61**(3), 986–990 (2014). [CrossRef] [PubMed]

**6. **J. Shi, X. Cao, F. Liu, B. Zhang, J. Luo, and J. Bai, “Greedy reconstruction algorithm for fluorescence molecular tomography by means of truncated singular value decomposition conversion,” J. Opt. Soc. Am. A **30**(3), 437–447 (2013). [CrossRef] [PubMed]

**7. **J. Shi, B. Zhang, F. Liu, J. Luo, and J. Bai, “Efficient L1 regularization-based reconstruction for fluorescent molecular tomography using restarted nonlinear conjugate gradient,” Opt. Lett. **38**(18), 3696–3699 (2013). [CrossRef] [PubMed]

**8. **J. Shi, F. Liu, J. Luo, and J. Bai, “Depth compensation in fluorescence molecular tomography using an adaptive support driven reweighted L1-minimization algorithm,” Proc. SPIE **9216**, 921603 (2014).

**9. **J. Shi, F. Liu, G. Zhang, J. Luo, and J. Bai, “Enhanced spatial resolution in fluorescence molecular tomography using restarted L1-regularized nonlinear conjugate gradient algorithm,” J. Biomed. Opt. **19**(4), 046018 (2014). [CrossRef] [PubMed]

**10. **J. Shi, F. Liu, J. Zhang, J. Luo, and J. Bai, “Fluorescence molecular tomography reconstruction via discrete cosine transform-based regularization,” J. Biomed. Opt. **20**(5), 055004 (2015). [CrossRef] [PubMed]

**11. **S. C. Davis, H. Dehghani, J. Wang, S. Jiang, B. W. Pogue, and K. D. Paulsen, “Image-guided diffuse optical fluorescence tomography implemented with Laplacian-type regularization,” Opt. Express **15**(7), 4066–4082 (2007). [CrossRef] [PubMed]

**12. **P. K. Yalavarthy, B. W. Pogue, H. Dehghani, C. M. Carpenter, S. Jiang, and K. D. Paulsen, “Structural information within regularization matrices improves near infrared diffuse optical tomography,” Opt. Express **15**(13), 8043–8058 (2007). [CrossRef] [PubMed]

**13. **X. Guo, X. Liu, X. Wang, F. Tian, F. Liu, B. Zhang, G. Hu, and J. Bai, “A Combined Fluorescence and Microcomputed Tomography System for Small Animal Imaging,” IEEE Trans. Biomed. Eng. **57**(12), 2876–2883 (2010). [CrossRef] [PubMed]

**14. **S. C. Davis, B. W. Pogue, R. Springett, C. Leussler, P. Mazurkewitz, S. B. Tuttle, S. L. Gibbs-Strauss, S. S. Jiang, H. Dehghani, and K. D. Paulsen, “Magnetic resonance-coupled fluorescence tomography scanner for molecular imaging of tissue,” Rev. Sci. Instrum. **79**(6), 064302 (2008). [CrossRef] [PubMed]

**15. **L. Cao and J. Peter, “Bayesian reconstruction strategy of fluorescence-mediated tomography using an integrated SPECT-CT-OT system,” Phys. Med. Biol. **55**(9), 2693–2708 (2010). [CrossRef] [PubMed]

**16. **S. R. Arridge, “Optical tomography in medical imaging,” Inverse Probl. **15**(2), R41–R93 (1999). [CrossRef]

**17. **M. A. O’Leary, D. A. Boas, X. D. Li, B. Chance, and A. G. Yodh, “Fluorescence lifetime imaging in turbid media,” Opt. Lett. **21**(2), 158–160 (1996). [CrossRef] [PubMed]

**18. **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,” Med. Phys. **22**(11), 1779–1792 (1995). [CrossRef] [PubMed]

**19. **W. He, G. Zhang, H. Pu, F. Liu, X. Cao, J. Luo, and J. Bai, “Modified forward model for eliminating the time-varying impact in fluorescence molecular tomography,” J. Biomed. Opt. **19**(5), 056012 (2014). [CrossRef] [PubMed]

**20. **A. Tikhonov, “Solving ill-conditioned and singular linear systems: A tutorial on regularization,” Soviet Math. Dokl. **5**, 1035–1038 (1963).

**21. **M. Chen, H. Su, Y. Zhou, C. Cai, D. Zhang, and J. Luo, “Automatic selection of regularization parameters for dynamic fluorescence molecular tomography: A comparison of L-curve and U-curve methods,” Biomed. Opt. Express **7**, 5021–5041 (2016).

**22. **B. Dogdas, D. Stout, A. F. Chatziioannou, and R. M. Leahy, “Digimouse: a 3D whole body mouse atlas from CT and cryosection data,” Phys. Med. Biol. **52**(3), 577–587 (2007). [CrossRef] [PubMed]

**23. **S. C. Davis, K. S. Samkoe, K. M. Tichauer, K. J. Sexton, J. R. Gunn, S. J. Deharvengt, T. Hasan, and B. W. Pogue, “Dynamic dual-tracer MRI-guided fluorescence tomography to quantify receptor density in vivo,” Proc. Natl. Acad. Sci. U.S.A. **110**(22), 9025–9030 (2013). [CrossRef] [PubMed]

**24. **Y. Zhou, H. Guang, H. Pu, J. Zhang, and J. Luo, “Unmixing multiple adjacent fluorescent targets with multispectral excited fluorescence molecular tomography,” Appl. Opt. **55**(18), 4843–4849 (2016). [CrossRef] [PubMed]