CN105894537A - Emission type optical projection tomography attenuation and scattering correction method - Google Patents

Abstract

The invention discloses an emission type optical projection tomography (OPT) attenuation and scattering correction method. An OPT model under the condition of weak scattering and an OPT fluorescence reconstruction algorithm are combined. OPT fluorescence signal attenuation correction is conducted using a normalized Born ratio method; the scattering component in measured data is simulated using the Monte Carlo technology, and the scattering component is removed from detection data; and through an iterative framework effectively combining attenuation correction and scattering component estimation, quantitative 3D reconstruction of OPT fluorescent dye concentration is realized. The influence of both absorption and scattering is considered comprehensively in quantitative reconstruction of fluorescent dye concentration, the scattering component in signals is estimated using the Monte Carlo technology and removed, and quantitative reconstruction of fluorescence dye concentration is realized.

Description

A kind of emission-type optical projection fault imaging decay and scatter correction method

Technical field

The invention belongs to medical image processing technical field, particularly relate to a kind of emission-type optical projection fault imaging decay and scatter correction method.

Background technology

Optical projection fault imaging OPT is a kind of novel high-resolution three-dimension molecular image imaging technique, and its image-forming principle is similar with the principle of X ray computer fault imaging.OPT uses visible ray as irradiation source, the imaging of the transparent or semitransparent sample being highly suitable in the 1-10mm range scales such as toy embryo, organ, fruit bat, nematicide.OPT technology makes people can obtain the three dimensional structure picture of high-resolution without the integrity losing biological tissue's organ, and can realize molecular specificity imaging by fluorescent labelling techniques.OPT has the advantages that bimodulus merges, and this just meets field of biomedical research and utilizes multimodality fusion more fully to understand the demand of Biont information.Therefore, OPT has had been subjected to the extensive concern of each research unit, it is applied in the research of gene expression, protein interaction, Preclinical Drug research and development etc., imaging object covers toy organ, embryo, fruit bat, nematicide, Brachydanio rerio, locust, arabidopsis isotype biology, provide strong research means for small sample entirety imaging, a series of biological basis progress of research will be promoted.

At body bio-imaging, emission-type OPT introduced some new problems, the fluorescent photon produced that is excited in emission-type OPT imaging is understood some photon in biological tissue's communication process and is absorbed by biological tissue, and traditional emission-type OPT method for reconstructing does not accounts for the impact of the decay that fluorescence signal is caused by this absorption effect, therefore its fluorescence signal Quantitative Reconstruction is inaccurate.Related experiment proves, even if sample interior places the fluorescent dye of same concentrations, owing to fluorescence information transmits out signal attenuation degree difference from sample, the fluorescence signal intensity also resulting in reconstruction exists notable difference；Being simultaneous for weak scattering sample, the scattering effect of sample also can affect the accuracy that fluorescence signal is quantitative.

Summary of the invention

It is an object of the invention to provide a kind of emission-type optical projection fault imaging for weak scattering sample to decay and scatter correction method, aim to solve the problem that traditional emission-type OPT method for reconstructing does not accounts for the impact of the decay that fluorescence signal is caused by absorption effect, and the inaccurate problem of fluorescence signal Quantitative Reconstruction caused because of weak scattering.

nullThe present invention is realized in，A kind of emission-type optical projection fault imaging for weak scattering sample is decayed and scatter correction method，OPT imaging model under weak scattering situation is combined with OPT fluorescence algorithm for reconstructing by the decay of described emission-type optical projection fault imaging with scatter correction method，First scattering impact is ignored，Directly method based on normalization Born ratio is utilized to carry out the preliminary correction for attenuation of OPT fluorescence signal detection data，Then filter back-projection reconstruction algorithm is utilized to calculate the concentration of fluorescent dye in sample，Again based on the concentration rebuilding the fluorescent dye obtained，The Monte Carlo technology emulation of employing generates corresponding scattering component，And from detection data, scattering component is removed，And then further the detection data removing scattering component is carried out correction for attenuation based on normalization Born ratio、The sequence of operations such as the concentration reestablishing of fluorescent dye，Constitute correction for attenuation and scattering component estimates effective iteration framework combined，Until the fluorescent dye concentration of twice reconstruction front and back is less than the threshold value set，Finally realize the quantitative three-dimensional reconstruction of OPT fluorescent dye concentration.

Further, the decay of described emission-type optical projection fault imaging comprises the following steps with scatter correction method:

Step one, uses normalized Born ratio method that measurement data is carried out correction for attenuation；

Step 2, utilizes filtered back projection (FBP) algorithm for reconstructing to calculate the concentration of fluorescent dye in sample；

Step 3, by the scattering component in the fluorescence data that Monte Carlo simulation estimate is measured；

Step 4, removes scattering component, again carries out correction for attenuation and the calculating of fluorescent dye concentration；

Step 5, the difference of twice reconstructed results and threshold comparison, iteration carries out simulation calculation.

Further, described employing normalized Born ratio method carries out correction for attenuation to measurement data and specifically includes:

According to the imaging model of emission-type OPT fluorescence imaging, exciting light fromPoint incides sample match liquid pool,Sample is excited by point, and fromPlace's injection sample cell, the fluorescence signal that CCD receives is expressed as:

$g_{out}^f={\∫}_{\stackrel{\→}{r_s}}^{\stackrel{\→}{r_d}}S\left(\stackrel{\→}{r_s}\right){\∫}_{\stackrel{\→}{r_s}}^{\stackrel{\→}{r_d}}\exp (-u_t(\stackrel{\→}{r}\left)\right){\mathrm{drdr}}_s;$

$S\left(\stackrel{\→}{r_s}\right)=g_{in}{\∫}_{\stackrel{\→}{r_0}}^{\stackrel{\→}{r_s}}\exp (-u_t(\stackrel{\→}{r}\left)\right)dr\·\mathrm{\ϵ}\mathrm{\η}\mathrm{\σ}\left(\stackrel{\→}{r_s}\right);$

WhereinRepresent the signal that when excitation wavelength does not mates, CCD collects with filter plate,Represent shot point,Represent transmitted wave and filter plate contact point,Representing the fluorescence signal intensity exciting generation, ε and η represents extinction coefficient and conversion quantum efficiency respectively,Represent the concentration of fluorescent dye, thus know:

$g_{out}^f=g_{in}{\∫}_{\stackrel{\→}{r_0}}^{\stackrel{\→}{r_d}}\exp (-u_t(\stackrel{\→}{r}\left)\right)dr\·{\∫}_{\stackrel{\→}{r_0}}^{\stackrel{\→}{r_d}}\mathrm{\ϵ}\mathrm{\η}\mathrm{\σ}\left(\stackrel{\→}{r_s}\right){\mathrm{dr}}_s;$

The parameter meaning repeated in the most several formula is just the same, context-sensitive part seen from concrete meaning.When signals collecting is to use the filter plate mated with excitation wavelength to carry out signals collecting, the signal that CCD detects is:

$g_{out}^0=g_{in}{\∫}_{\stackrel{\→}{r_0}}^{\stackrel{\→}{r_d}}\exp (-u_t(\stackrel{\→}{r}\left)\right)dr;$

HaveWherein G_fFor exciting light along the direction of propagationRadon conversion,Being respectively two kinds of signals exciting CCD camera under situation to collect, other parameter has been introduced above.

Further, during the described FBP of utilization algorithm for reconstructing calculates sample, the concentration of fluorescent dye specifically includes:

Converted, i.e. by the Radon obtaining fluorescent dye concentrationG_fFor exciting light along the direction of propagationRadon conversion,Being respectively two kinds of signals exciting CCD camera under situation to collect, ε and η represents extinction coefficient and conversion quantum efficiency respectively,Concentration for fluorescent dye.The Radon conversion of multiple directions can reconstruct fluorescent dye concentrationCollect measurement data G of 360 degree_fAfter, i.e. calculated the concentration of fluorescent dye by inverse Radon transform, i.e. calculate σ=FBP (G_f), obtain fluorescent dye concentration σ.

Further, the scattering component in the described fluorescence data measured by Monte Carlo simulation estimate is specifically included:

Utilize fluorescent dye concentration σ reconstructed, the information such as the absorptance and the scattering coefficient that simultaneously utilize sample, set up model by Monte Carlo emulation, by adjustment simulation parameter to simulate actual experiment situation, therefrom can estimate the scattering component in fluorescence data.

Further, described removal scattering component, the calculating again carrying out correction for attenuation and fluorescent dye concentration specifically includes:

By the scattering component obtained, measurement data is corrected further, from measurement data, i.e. rejects scattering component, then carry out correction for attenuation again and fluorescent dye concentration and calculate, i.e. repeat following steps:

Utilize FBP algorithm for reconstructing to calculate the concentration of fluorescent dye in sample, specifically include:

Converted, i.e. by the Radon obtaining fluorescent dye concentrationG_fFor exciting light along the direction of propagationRadon conversion, multiple directions Radon conversion can reconstruct fluorescent dye concentrationCollect measurement data G of 360 degree_fAfter, i.e. calculated the concentration of fluorescent dye by inverse Radon transform, i.e. calculate σ=FBP (G_f), obtain fluorescent dye concentration σ；

Further, the difference of twice reconstructed results and threshold comparison, iteration carries out simulation calculation, the fluorescent dye concentration reconstructed and last reconstructed results contrast, if the two difference is less than given threshold value, then terminator, output result, complete decay and scatter correction, if being unsatisfactory for end condition, then proceeding MonteCarlo emulation, the detection data of removal scattering component carries out the sequence of operations such as the concentration reestablishing of correction for attenuation based on normalization Born ratio, fluorescent dye, until meeting end condition.

The emission-type optical projection fault imaging that the present invention provides is decayed and scatter correction method, is combined with weak scattering OPT imaging model, it is achieved rebuild at body OPT fluorescent dye concentration quantitative；OPT imaging model under weak scattering situation is combined with OPT fluorescence algorithm for reconstructing, and utilize the absorptance and scattering coefficient rebuilding out in the imaging of transmission-type OPT, method based on normalization Born ratio is utilized to carry out OPT fluorescent signal decay correction, use the scattering component in Monte Carlo technology simulated measurement data, and remove it from detection data, to eliminate the impact that fluorescent quantitation is rebuild by scattering, correction for attenuation and scattering component are estimated the iteration framework effectively combined, fluorescent signal decay correction and the removal of scattering component is completed by successive ignition, thus realize the quantitative three-dimensional reconstruction of OPT fluorescent dye concentration.

The present invention absorbs and the impact of scattering realizing having considered during fluorescent dye concentration quantitative is rebuild, and uses scattering component that Monte Carlo technology estimated in signal and removed, and then realizes the Quantitative Reconstruction of fluorescent dye concentration.

Accompanying drawing explanation

Fig. 1 is that the emission-type optical projection fault imaging that the embodiment of the present invention provides is decayed and scatter correction method flow chart.

Fig. 2 is the flow chart of the embodiment that the embodiment of the present invention provides.

Fig. 3 is the emission-type OPT image acquisitions schematic diagram that the embodiment of the present invention provides.

Fig. 4 is fluorescence excitation and the signal attenuation schematic diagram of embodiment of the present invention offer.

Detailed description of the invention

In order to make the purpose of the present invention, technical scheme and advantage clearer, below in conjunction with embodiment, the present invention is further elaborated.Should be appreciated that specific embodiment described herein, only in order to explain the present invention, is not intended to limit the present invention.

OPT imaging model under weak scattering situation is combined by the present invention with OPT fluorescence algorithm for reconstructing, method based on normalization Born ratio is utilized to carry out OPT fluorescent signal decay correction, use the scattering component in Monte Carlo technology simulated measurement data, and remove it from detection data, study the iteration framework that correction for attenuation and scattering component are estimated effectively to be combined, thus realize the quantitative three-dimensional reconstruction of OPT fluorescent dye concentration.

Below in conjunction with the accompanying drawings the application principle of the present invention is explained in detail.

As it is shown in figure 1, the emission-type optical projection fault imaging decay of the embodiment of the present invention comprises the following steps with scatter correction method:

S101: use normalized Born ratio method that measurement data is carried out correction for attenuation；

S102: utilize FBP algorithm for reconstructing to calculate the concentration of fluorescent dye in sample；

S103: by the scattering component in the fluorescence data that Monte Carlo simulation estimate is measured；

S104: remove scattering component, again carry out correction for attenuation and the calculating of fluorescent dye concentration；

The difference of S105: twice reconstructed results and threshold comparison, iteration carries out simulation calculation.

Below in conjunction with specific embodiment, the application principle of the present invention is further described.

As in figure 2 it is shown, the emission-type optical projection fault imaging decay of the embodiment of the present invention comprises the following steps with scatter correction method:

Step 1, uses normalized Born ratio method that measurement data carries out correction for attenuation:

First, according to the imaging model of emission-type OPT fluorescence imaging, sample is through exciting the fluorescence signal portion of generation directly to be received by CCD through ballistic propagation, and part is received by CCD after being scattered, remaining photon direct projection or scatter to other directions and cannot be detected by CCD.Assume exciting light fromPoint incides sample match liquid pool,Sample is excited by point, and fromPlace's injection sample cell.Wavelength ratio in view of exciting light and transmitting light is closer to, and ignores the two difference between optical coefficient, simultaneously takes account of fluorescent dye and is likely distributed in the optional position that light passes, then the fluorescence signal that CCD receives can be expressed as:

$g_{out}^f={\∫}_{\stackrel{\→}{r_s}}^{\stackrel{\→}{r_d}}S\left(\stackrel{\→}{r_s}\right){\∫}_{\stackrel{\→}{r_s}}^{\stackrel{\→}{r_d}}\exp (-u_t(\stackrel{\→}{r}\left)\right){\mathrm{drdr}}_s;$

$S\left(\stackrel{\→}{r_s}\right)=g_{in}{\∫}_{\stackrel{\→}{r_0}}^{\stackrel{\→}{r_s}}\exp (-u_t(\stackrel{\→}{r}\left)\right)dr\·\mathrm{\ϵ}\mathrm{\η}\mathrm{\σ}\left(\stackrel{\→}{r_s}\right);$

WhereinRepresenting the fluorescence signal intensity exciting generation, ε and η represents extinction coefficient and conversion quantum efficiency respectively,Represent the concentration of fluorescent dye, thus understand:

$g_{out}^f=g_{in}{\∫}_{\stackrel{\→}{r_0}}^{\stackrel{\→}{r_d}}\exp (-u_t(\stackrel{\→}{r}\left)\right)dr\·{\∫}_{\stackrel{\→}{r_0}}^{\stackrel{\→}{r_d}}\mathrm{\ϵ}\mathrm{\η}\mathrm{\σ}\left(\stackrel{\→}{r_s}\right){\mathrm{dr}}_s;$

Notice that the signal that CCD detects is when using the filter plate mated with excitation wavelength to carry out signals collecting:

$g_{out}^0=g_{in}{\∫}_{\stackrel{\→}{r_0}}^{\stackrel{\→}{r_d}}\exp (-u_t(\stackrel{\→}{r}\left)\right)dr;$

Thus have

The method of fluorescent signal decay carried out above correction is referred to as normalized Born and compares method.

Step 2, utilizes FBP algorithm for reconstructing to calculate the concentration of fluorescent dye in sample:

The Radon conversion of fluorescent dye concentration is obtained, i.e. by step 1From above formula, G_fFor exciting light along the direction of propagationRadon conversion, Radon convert and inverse transformation theory, measure multiple directions Radon convert can reconstruct fluorescent dye concentrationThe realization in OPT imaging of multidirectional DATA REASONING can utilize automatically controlled rotary apparatus to drive, and rotary sample utilizes light source to irradiate simultaneously, CCD carries out data acquisition and both may be used.When measurement data G collecting 360 degree_fAfter, the concentration of fluorescent dye can be calculated by inverse Radon transform.

I.e. calculate σ=FBP (G_f), thus obtain fluorescent dye concentration σ.

Step 3, the scattering component by the fluorescence data that Monte Carlo simulation estimate is measured:

Concrete methods of realizing is to utilize fluorescent dye concentration σ reconstructed in step 2, the information such as the absorptance and the scattering coefficient that simultaneously utilize sample, model is set up by Monte Carlo emulation, by adjustment simulation parameter to simulate actual experiment situation, the scattering component in fluorescence data therefrom can be estimated.

Amount of calculation based on Monte Carlo simulation calculating is the biggest, also there is the feature of highly-parallel simultaneously, it is also admirably suitable for GPU to accelerate, therefore the photon based on Monte Carlo method of exploitation propagates emulation platform in biological tissues before utilizing this laboratory, Molecular Optical Simulation Environment (MOSE), research uses GPU technology to carry out algorithm acceleration, it is achieved the acceleration of scattering component Monte Carlo simulation calculation.

Step 4, removes scattering component, again carries out correction for attenuation and the calculating of fluorescent dye concentration:

By the scattering component obtained in step 3, measurement data is corrected further, from measurement data, i.e. reject scattering component, then carry out correction for attenuation again and fluorescent dye concentration calculates, i.e. repeat step 2 and 3.

Step 5, the difference of twice reconstructed results and threshold comparison, iteration carries out simulation calculation.

Reconstructed results and last reconstructed results contrast, if the two difference is less than given threshold value, then terminator, output result, completes decay and scatter correction, if being unsatisfactory for end condition, then proceed the sequence of operations such as Monte Carlo emulation, until meeting end condition.

Below in conjunction with the accompanying drawings 3, the reconstructed results of the present invention is described in detail by accompanying drawing 4.

Accompanying drawing 3 is the emission-type OPT image acquisitions schematic diagram of the present invention.Sample is through exciting the fluorescence signal portion of generation directly to be received by CCD through ballistic propagation, and part is received by CCD after being scattered, remaining photon direct projection or scatter to other directions and cannot be detected by CCD.Exciting light and launch light and propagate all to be decayed by biological tissue in biological tissue and affected, this is irrespective in traditional OPT imaging, is the most just the core content that considered of the present invention.

Accompanying drawing 4 is the decay schematic diagram of fluorescence excitation and signal.Wherein figure (a) be the filter plate that mates with wavelength of transmitted light of employing schematic diagram when carrying out signals collecting, filter plate now does not mates with excitation wavelength, and figure (b) be the filter plate that mates with excitation wavelength of employing schematic diagram when carrying out signals collecting.Two width schematic diagrams all being assumed, exciting light is from pointIncide sample match liquid pool,Sample is excited by point, and fromPlace's injection sample cell.

The foregoing is only presently preferred embodiments of the present invention, not in order to limit the present invention, all any amendment, equivalent and improvement etc. made within the spirit and principles in the present invention, should be included within the scope of the present invention.

Claims (7)

1. an emission-type optical projection fault imaging decay and scatter correction method, it is characterised in that described Emission-type optical projection fault imaging decay with scatter correction method by OPT imaging model under weak scattering situation with OPT fluorescence algorithm for reconstructing combines, and utilizes method based on normalization Born ratio to carry out OPT detection data The preliminary correction for attenuation of fluorescence signal；Then filter back-projection reconstruction algorithm is utilized to calculate fluorescent dye in sample Concentration, then based on rebuilding the concentration of fluorescent dye obtained；It is right that the Monte Carlo technology emulation of employing generates The scattering component answered, and from detection data, scattering component is removed, to the detection data removing scattering component Carry out the concentration reestablishing operation of correction for attenuation based on normalization Born ratio, fluorescent dye, constitute correction for attenuation Effective iteration framework combined is estimated, until the fluorescent dye concentration of twice reconstruction front and back is less than with scattering component The threshold value set；Finally realize the quantitative three-dimensional reconstruction of OPT fluorescent dye concentration.

2. emission-type optical projection fault imaging decay as claimed in claim 1 and scatter correction method, its Being characterised by, the decay of described emission-type optical projection fault imaging comprises the following steps with scatter correction method:

Step one, uses normalized Born ratio method that measurement data is carried out correction for attenuation；

Step 2, utilizes FBP algorithm for reconstructing to calculate the concentration of fluorescent dye in sample；

Step 3, by the scattering component in the fluorescence data that Monte Carlo simulation estimate is measured；

Step 4, removes scattering component, again carries out correction for attenuation and the calculating of fluorescent dye concentration；

Step 5, the difference of twice reconstructed results and threshold comparison, iteration carries out simulation calculation.

3. emission-type optical projection fault imaging decay as claimed in claim 2 and scatter correction method, its Being characterised by, described employing normalized Born ratio method carries out correction for attenuation to measurement data and specifically includes:

According to the imaging model of emission-type OPT fluorescence imaging, exciting light fromPoint incides sample match liquid Pond,Sample is excited by point, and fromPlace's injection sample cell, the fluorescence signal that CCD receives represents For:

$g_{out}^f={\∫}_{\stackrel{\→}{r_s}}^{\stackrel{\→}{r_d}}S\left(\stackrel{\→}{r_s}\right){\∫}_{\stackrel{\→}{r_s}}^{\stackrel{\→}{r_d}}\exp (-u_t(\stackrel{\→}{r}\left)\right){\mathrm{drdr}}_s;$

$S\left(\stackrel{\→}{r_s}\right)=g_{in}{\∫}_{\stackrel{\→}{r_0}}^{\stackrel{\→}{r_s}}\exp (-u_t(\stackrel{\→}{r}\left)\right)dr\·\mathrm{\ϵ}\mathrm{\η}\mathrm{\σ}\left(\stackrel{\→}{r_s}\right);$

WhereinRepresent the signal that when excitation wavelength does not mates, CCD collects with filter plate,Expression excites Point,Represent transmitted wave and filter plate contact point,Represent the fluorescence signal intensity exciting generation, ε and η Represent extinction coefficient and conversion quantum efficiency respectively,Represent the concentration of fluorescent dye, thus know:

$g_{out}^f=g_{in}{\∫}_{\stackrel{\→}{r_0}}^{\stackrel{\→}{r_d}}\exp (-u_t(\stackrel{\→}{r}\left)\right)dr\·{\∫}_{\stackrel{\→}{r_0}}^{\stackrel{\→}{r_d}}\mathrm{\ϵ}\mathrm{\η}\mathrm{\σ}\left(\stackrel{\→}{r_s}\right){\mathrm{dr}}_s;$

When signals collecting is to use the filter plate mated with excitation wavelength to carry out signals collecting, CCD detects To signal be:

$g_{out}^0=g_{in}{\∫}_{\stackrel{\→}{r_0}}^{\stackrel{\→}{r_d}}\exp (-u_t(\stackrel{\→}{r}\left)\right)dr;$

HaveWherein G_fFor exciting light along the direction of propagationRadon conversion, It is respectively two kinds of signals exciting CCD camera under situation to collect.

4. emission-type optical projection fault imaging decay as claimed in claim 2 and scatter correction method, its Being characterised by, the described FBP of utilization algorithm for reconstructing calculates the concentration of fluorescent dye in sample and specifically includes:

Converted, i.e. by the Radon obtaining fluorescent dye concentrationG_fFor exciting light edge The direction of propagationRadon conversion, Being respectively two kinds excites CCD camera under situation to collect Signal, ε and η represents extinction coefficient and conversion quantum efficiency respectively,Concentration for fluorescent dye；Many The Radon conversion in individual direction can reconstruct fluorescent dye concentrationCollect the measurement data of 360 degree G_fAfter, i.e. calculated the concentration of fluorescent dye by inverse Radon transform, i.e. calculate σ=FBP (G_f), obtain Fluorescent dye concentration σ.

5. emission-type optical projection fault imaging decay as claimed in claim 2 and scatter correction method, its Being characterised by, the scattering component in the described fluorescence data measured by Monte Carlo simulation estimate is concrete Including:

Utilize fluorescent dye concentration σ reconstructed, utilize absorptance and the scattering coefficient information of sample simultaneously, Model is set up, by adjusting simulation parameter to simulate actual experiment situation, therefrom by Monte Carlo emulation The scattering component in fluorescence data can be estimated.

6. emission-type optical projection fault imaging decay as claimed in claim 2 and scatter correction method, its Being characterised by, described removal scattering component, the calculating again carrying out correction for attenuation and fluorescent dye concentration is concrete Including:

By the scattering component obtained, measurement data is corrected further, from measurement data, i.e. reject scattering Component, then carries out correction for attenuation again and fluorescent dye concentration and calculates, and i.e. repeats following steps:

Utilize FBP algorithm for reconstructing to calculate the concentration of fluorescent dye in sample, specifically include:

Converted, i.e. by the Radon obtaining fluorescent dye concentrationG_fFor exciting light edge The direction of propagationRadon conversion, can to reconstruct fluorescent dye dense in the Radon conversion of multiple directions DegreeCollect measurement data G of 360 degree_fAfter, fluorescent dye is i.e. calculated by inverse Radon transform Concentration, i.e. calculate σ=FBP (G_f), obtain fluorescent dye concentration σ；

By the scattering component in the fluorescence data that Monte Carlo simulation estimate is measured, specifically include:

Utilize fluorescent dye concentration σ reconstructed, utilize absorptance and the scattering coefficient information of sample simultaneously, Model is set up, by adjusting simulation parameter to simulate actual experiment situation, therefrom by Monte Carlo emulation I.e. estimate the scattering component in fluorescence data.

7. emission-type optical projection fault imaging decay as claimed in claim 2 and scatter correction method, its Being characterised by, the difference of twice reconstructed results and threshold comparison, iteration carries out simulation calculation, and reconstructed results is with upper Reconstructed results contrast, if the two difference is less than given threshold value, then terminator, exports result, Completing decay and scatter correction, if being unsatisfactory for end condition, then proceeding Monte Carlo emulation, Until meeting end condition.