CN101103905A - Optical system for detecting turbidity medium and biological tissue optical parameter and detecting method - Google Patents
Optical system for detecting turbidity medium and biological tissue optical parameter and detecting method Download PDFInfo
- Publication number
- CN101103905A CN101103905A CNA2007100585937A CN200710058593A CN101103905A CN 101103905 A CN101103905 A CN 101103905A CN A2007100585937 A CNA2007100585937 A CN A2007100585937A CN 200710058593 A CN200710058593 A CN 200710058593A CN 101103905 A CN101103905 A CN 101103905A
- Authority
- CN
- China
- Prior art keywords
- optical
- optical fiber
- biological tissue
- measuring
- parameter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000003287 optical effect Effects 0.000 title claims abstract description 149
- 238000000034 method Methods 0.000 title abstract description 19
- 239000013307 optical fiber Substances 0.000 claims abstract description 121
- 238000004364 calculation method Methods 0.000 claims abstract description 22
- 238000012545 processing Methods 0.000 claims abstract description 17
- 238000001228 spectrum Methods 0.000 claims abstract description 15
- 239000000523 sample Substances 0.000 claims description 53
- 238000005259 measurement Methods 0.000 claims description 42
- 238000012360 testing method Methods 0.000 claims description 42
- 239000000835 fiber Substances 0.000 claims description 31
- 230000005540 biological transmission Effects 0.000 claims description 16
- 238000003556 assay Methods 0.000 claims description 15
- 238000005286 illumination Methods 0.000 claims description 12
- 230000011664 signaling Effects 0.000 claims description 11
- 230000003321 amplification Effects 0.000 claims description 5
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 5
- 230000005611 electricity Effects 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 3
- 238000003780 insertion Methods 0.000 claims 1
- 230000037431 insertion Effects 0.000 claims 1
- 238000012546 transfer Methods 0.000 abstract description 9
- 210000001519 tissue Anatomy 0.000 description 37
- 230000006870 function Effects 0.000 description 14
- 230000005855 radiation Effects 0.000 description 9
- 238000002834 transmittance Methods 0.000 description 7
- 238000009792 diffusion process Methods 0.000 description 6
- 238000002835 absorbance Methods 0.000 description 4
- 230000031700 light absorption Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 239000004065 semiconductor Substances 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- 238000000149 argon plasma sintering Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000008676 import Effects 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000005055 memory storage Effects 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 230000006854 communication Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 210000004087 cornea Anatomy 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 230000002500 effect on skin Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 210000002615 epidermis Anatomy 0.000 description 1
- 210000003238 esophagus Anatomy 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 239000008187 granular material Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 229940028435 intralipid Drugs 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000000049 pigment Substances 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 210000003437 trachea Anatomy 0.000 description 1
Images
Abstract
Disclosed are an optical fiber system and a measuring method for measuring a turbid medium and a biological tissue optical parameter, comprising light source inputting, optical fiber measuring, controlling, data processing and calculating. The light source inputting is connected with the optical fiber measuring which is connected with controlling, data processing and calculating which are connected with the light source, and the optical fiber measuring is connected with the turbid medium or biological issue. The method is that real-time spectrum signals output by a measuring apparatus stores optical signals with various wavelengths according to the wavelengths; in a light transfer theoretical calculation sub-procedure, the stored real-time optical signals with various wavelengths are demodulated by mathematic treatment and then the actual measured optical signals are calculated; the actual measured optical signals are compared with the calculated optical signals, and then the optical parameter of the measured sample is adjusted to make the differences between the calculated optical signals and the actual measured optical signals less than the preset value, the optical parameter of the measured sample is output finally. The measuring system of the invention is simple and can accurately measure the optical parameter of an original body or living body in no need of leaving the living body sample.
Description
Technical field
The present invention relates to a kind of fiber-optic probe system.Particularly relate to a kind of turbid medium and the mensuration turbid medium of biological tissue optical parameter and the fibre system and assay method of biological tissue optical parameter that can be used for measuring exactly body or live body.
Background technology
The quantitative Diagnosis and the treatment that are used for disease along with optical means are more and more monitored, and research or medical personnel usually need to know the irradiates light energy is how to distribute and how to produce to reflect and transmission signal after entering patient's organ biological tissue.Except that minority tissues such as eye's cornea, the biological tissue of human body or animal is the optical opacity medium, also both light absorption and light scattering and be stored in light and the interaction of biological tissue, based on scattering, interactional room and time distributes and has at random character at visible light and near infrared band.The energy transmission of light in tissue and be distributed as a complex conditions problem in this case need be described and find the solution with mathematics physics model accurately.The turbid medium and the biological organism optical model of extensive use at present are radiation transfer theory, and this theory is defined as absorptance with the optical parametric of biological tissue, scattering coefficient and scattering phase function.Absorptance is represented photon absorbed average time of unit propagation distance in medium, usually with in-house heterogeneity such as biomacromolecule kind and concentration, blood, the size of pigment granule with what etc. relevant.Scattering coefficient and scattering phase function are then represented the Probability Distribution of photon average time that the unit propagation distance is scattered in medium and scattering angle.Under the situation of known scattering phase functional form, as Han Ni-Greenstein (Henyey-Greenstein) scattering phase function etc., the scattering phase function can be by one or more scalar parameters decisions.Han Ni-Greenstein scattering phase function is determined by an anisotropic parameters, anisotropic parameters is defined as the meansigma methods of scattering angle cosine, go through visible list of references (Z.Song for example, K.Dong, X.H.Hu, and J.Q.Lu, " MonteCarlo Simulation of Converging Laser Beams Propagating in Biological Tissues ", Applied Optics, vol.38, pp.2944-2949 (1999)).
Many research documents show that the optical parametric of turbid medium is made up of the microstructure of medium and macroscopic view distribution decision; For the turbid medium of biological tissue's class, then by decisions such as its cell category and distributed architectures.Between different tissues in the human organ such as cancer focus and normal structure, the optical parametric difference between inhomogeneity normal structure such as epidermis and the dermal tissue.The energy transmission of light in turbid medium or biologic-organ-tissue and distribution can be determined according to radiation transfer theory by these optical parametrics.Therefore, the optical parametric of turbid medium or biological tissue is that the optical property of this type of material is carried out the requisite critical data of quantitative analysis.In addition, these optical parametrics change with optical wavelength usually, are the functions of wavelength therefore.The optical parametric of different media is specific and different function of wavelength.
Measure turbid medium or biological tissue optical parameter the irradiation incident illumination that medium is used, the different directions scattered light signal identical with position measurement and lambda1-wavelength in medium or outside the medium then need be provided.Light in communication process because and the interaction between scattering particles or the scattering center cause light or photon to relay the change of direction.The space of these scattering centers and time distribute and are generally at random in turbid medium or biological tissue, and its statistical nature is by the optical property decision of medium.The random nature of light scattering causes emergent light randomness on direction and the position distribution in three dimensions, and this makes the scattered light signal of turbid medium or biological tissue measure relatively difficulty.The measuring method of reporting at present comprises the method for utilizing integrating sphere measurement diffuse-reflectance and diffuse transmission signal and utilizes discrete detector or the method for detector array measurement reflected signal distribution.The shortcoming of these methods or be the measuring system complexity, or, all can't be used for measuring exactly body or live body optical parametric for leaving the sample of live body.Go through visible list of references (C.Chen for example, J.Q.Lu, H.Ding, K.M.Jacobs, Y.Du, and X.H.Hu, " A primary method for determination of opticalparameters of turbid samples and application to intralipid between 550 and 1630nm ", Optics Express, vol.14, pp.7420-7435 (2006)).
Summary of the invention
Technical problem to be solved by this invention is that a kind of turbid medium and the mensuration turbid medium of biological tissue optical parameter and the fibre system and assay method of biological tissue optical parameter that can be used for measuring exactly body or live body is provided.
The technical solution adopted in the present invention is: a kind of fibre system and assay method of measuring turbid medium and biological tissue optical parameter, and wherein, fibre system includes: the light source importation; The optical fiber measurement part; Handle and calculating section with control data, wherein, the light source importation partly is connected with optical fiber measurement, the optical fiber measurement part is connected with control, date processing and calculating section, control, date processing and calculating section link to each other with the light source importation again, and the optical fiber measurement part also connects tested turbid medium or biological tissue.
Described light source importation includes light source electricity consumption source temperature controller and light source, wherein, the input of power supply temperature controller links to each other with control, date processing and calculating section by cable, outfan links to each other with light source by power supply temperature control cable, and the output coupled end of light source partly is connected with optical fiber measurement.
Described optical fiber measurement partly includes and excites guiding optical fiber, measuring instrument and measuring optical fiber bundle, wherein, excite an end of guiding optical fiber to be connected with the light source of light source importation, measuring instrument connects control, date processing and calculating section and measuring optical fiber bundle respectively, excites the common tested turbid medium or the biological tissue of inserting of the other end of the other end and the measuring optical fiber bundle of guiding optical fiber.
Described control, date processing and calculating section are to be finished by computer.
Described light source can be obtained by one or more laser coherence light source.
Described light source also can be obtained by the incoherent light source of continuous spectrum, and its irradiates light wavelength is continuous distribution in the spectral domain of setting.
The irradiates light intensity modulated of described light source can by source current or mechanical chopper both one of, at 0.1 hertz to 100 mhz frequencies scope internal modulations.
The described diameter of guiding optical fiber that excites is between 10 microns～10 millimeters.
The described planar normal of guiding optic fibre end that excites becomes the interior a certain angle of 0 degree～90 degree scopes with fiber axis.
Described measuring optical fiber bundle is made up of 1～10 measuring optical fiber, and the diameter of the optical fiber that measuring optical fiber is intrafascicular is between 10 microns～10 millimeters.
The planar normal of described measuring optical fiber Shu Duantou is parallel with fiber axis or become 0 to spend the interior a certain angle of～90 degree scopes.
Described excite guiding optical fiber and measuring optical fiber bundle to be fixed in the same syringe needle after, insert in turbid medium or the biological tissue.
The described intracavity passage that excites guiding optical fiber and the fixing back of measuring optical fiber Shu Caiyong cable overcoat to enter human body by the service aisle of endoscope.
Described measuring instrument is made up of beam split part, photodetection part and analog-digital converter, and the optical signal that its measuring optical fiber is exported is at first launched its wave spectrum by the beam split part and become by photodetector to deliver to analog-digital converter after the signal of telecommunication amplification and become the real time spectrum signal.
Described measuring instrument photodetection part and analog-digital converter are formed, and the optical signals photodetector that measuring optical fiber is exported becomes to be delivered to analog-digital converter after the signal of telecommunication amplification and become the real-time optical signal.
Wherein, the optical fiber assay method is by finishing as the next stage:
Phase I: the optical signalling on a plurality of wavelength is stored respectively by its wavelength from the real time spectrum signal of measuring instrument output;
Second stage: the real-time optical signal of different wave length of storage is being undertaken calculating optical signal after the demodulation by Mathematical treatment on the different wave length;
Phase III: enter the subprogram of optical transmission Theoretical Calculation, according to the incident illumination parameter of input, test sample probe form parameter and test sample optical parametric initial value obtain to calculate optical signal;
Quadravalence section: compare measured light signal and the difference of calculating optical signal, and by adjust repeatedly the test sample optical parametric until the difference of calculating optical signal and measured light signal less than till the setting value in advance, and export test sample optical parametric numerical value on the irradiates light wavelength.
Optical signalling on described a plurality of wavelength includes the incident illumination parameter, test sample probe form parameter and test sample optical parametric initial value.
The subprogram of described optical transmission Theoretical Calculation is by finishing as the next stage:
Phase I: the input optical parameter, determine tracked total number of light photons N
0, and determine photon total distance of advancing at random according to absorptance;
Second stage: judge that tracked accumulative total number of photons N is greater than 1? also promptly advance direction and determine the photon free path of advancing at random of decision photon scattering angle, and follow the trail of photon to next scattering point;
Phase III: judge that the accumulative total traveling distance is greater than total distance? whether is measured optical fiber accepted? overflow test sample? Contact Boundary whether?
Quadravalence section: calculate optical signal, and judge that accumulative total number of photons N is whether greater than the total number of light photons N of incident beam
0
Five-stage: optical signal is calculated in EP (end of program) and output.
The fibre system and the assay method of mensuration turbid medium of the present invention and biological tissue optical parameter, measuring system is simple, and is easy to use, and the sample that does not need to leave live body just can be measured body or live body optical parametric exactly.
Description of drawings
Fig. 1 is a system structure sketch map of the present invention;
Fig. 2 is the sketch map that optical fiber is provided with in the embodiments of the invention;
Fig. 3 is the flow chart of assay method of the present invention;
Fig. 4 is the optical transmission Theoretical Calculation subroutine flow chart in the assay method of the present invention.
Wherein:
1: power supply temperature controller 2: power supply temperature control cable
3: light source 4: excite guiding optical fiber
5: protection needle guard 6: biological tissue
7: measuring optical fiber bundle 8: measuring instrument
9: cable 10: computer
11: cable
The specific embodiment
Below, describe the fibre system and the assay method of mensuration turbid medium of the present invention and biological tissue optical parameter in conjunction with the accompanying drawings and embodiments in detail.
Accurately measure the absorptance of turbid medium or biological tissue, scattering coefficient and anisotropic parameters need be measured the diffusion light on the different directions.Scattered light signal is defined as under the constant condition of optical wavelength its direction of propagation and departs from the optical signal that detects behind the former direction of propagation, the diffusion light signal definition is that light is after light source or the outgoing of radiation delivery system terminal, the optical signal that after the single or multiple scattering, detects, the diffusion light signal of forward direction diffusion light signal definition for measuring in extremely positive and negative 90 degree of polar angle 0 degree of emergent light direction, the back is to the diffusion light signal of diffusion light signal definition for measuring in positive and negative 90 degree of the polar angle of emergent light direction are spent to 180.
As shown in Figure 1, the fibre system of measurement turbid medium of the present invention and biological tissue optical parameter includes: the light source importation; The optical fiber measurement part; Handle and calculating section with control data, wherein, the light source importation partly is connected with optical fiber measurement, the optical fiber measurement part is connected with control, date processing and calculating section, control, date processing and calculating section link to each other with the light source importation again, and the optical fiber measurement part also connects tested turbid medium or biological tissue.
Described light source importation includes light source electricity consumption source temperature controller 1 and light source 3, and wherein, power supply temperature controller 1 provides light source required electric current, the temperature of control light source and output light modulation frequency.The input of power supply temperature controller 1 links to each other with control, date processing and calculating section by cable 11, and outfan links to each other with light source 3 by power supply temperature control cable 2, and the output coupled end of light source 3 partly is connected with optical fiber measurement.
A kind of implementation of power supply temperature controller 1 can comprise a temperature control current source and a laser current source of carrying out negative feedback control by temp-sensing element (as critesistor) signal.The electric current in temperature control current source exports the semiconductor cooler with power supply to, and reaches the purpose of control light-source temperature.Temperature control current source and laser current source belong to common circuit, have how tame photoelectric device company to sell, and the model of selling as the Thorlabs company of the U.S. is the power supply temperature controller of ITC102.Power supply temperature control cable 2 is made up of many cables, comprises the laser diode current line, the semiconductor cooler power line, and temp-sensing element holding wires etc. are common cable.
Described light source 3 can be obtained by one or more laser coherence light source, is single or discrete multi-wavelength.Also can be obtained by the incoherent light source of continuous spectrum, its irradiates light wavelength is continuous distribution in the spectral domain of setting.The irradiates light intensity modulated of light source 3 can by source current or mechanical chopper both one of, at 0.1 hertz to 100 mhz frequencies scope internal modulations.
Described optical fiber measurement partly includes and excites guiding optical fiber 4, measuring instrument 8 and measuring optical fiber bundle 7, wherein, excite an end of guiding optical fiber 4 to be connected with the light source 3 of optical signal importation, measuring instrument 8 connects control, date processing and calculating section and measuring optical fiber bundle 7 respectively, excites the common tested turbid medium or the biological tissue of inserting of the other end of the other end and the measuring optical fiber bundle 7 of guiding optical fiber 4.
The described diameter of guiding optical fiber 4 that excites is between 10 microns～10 millimeters.And excite the planar normal in guiding optical fiber 4 terminations to become the interior a certain angle of 0 degree～90 degree scopes with fiber axis.
Described measuring optical fiber bundle 7 is made up of 1～10 measuring optical fiber, and the diameter of the optical fiber in the measuring optical fiber bundle 7 is between 10 microns～10 millimeters.Be the diameter of the intrafascicular optical fiber of measuring optical fiber can be identical also can be inequality.Described measuring optical fiber bundle 7 planar normals are parallel with fiber axis or become the interior a certain angle of 0 degree～90 degree scopes, to accept the forward direction transmission, diffuse transmission or irreflexive scattered light signal.
Described guiding optical fiber 4 and the measuring optical fiber bundle 7 of exciting be by with after different methods of attachment or directly being fixed in the same syringe needle (as protection needle guard 5), inserts the optical parametric of mensuration turbid medium in turbid medium or the biological tissue or biological tissue 6.
Described guiding optical fiber 4 and the measuring optical fiber bundle 7 of exciting can also adopt the fixing back of cable overcoat to enter the intracavity passage of human body such as the optical parametric of trachea and esophagus mensuration tube wall tissue by the service aisle of endoscope.
Guiding optical fiber 4 and the measuring optical fiber bundle 7 bonded a kind of methods for designing of exciting of the present invention can be by shown in Fig. 2 a and Fig. 2 b, and Fig. 2 a is a front view, and Fig. 2 b is a lateral plan.Excite guiding optical fiber 4 to be used to import the irradiates light of exciting light scattered signal among the figure.Optical fiber a, b, c, d are measuring optical fiber bundle 7 among the figure, and wherein, optical fiber a and optical fiber b are for measuring the optical fiber of forward direction and diffuse transmission optical signal, and optical fiber c and optical fiber d are for measuring the optical fiber of the signal that diffuses.
Described measuring instrument 8 is made up of beam split part, photodetection part and analog-digital converter, and the optical signal that its measuring optical fiber is exported is at first launched its wave spectrum by the beam split part and become by photodetector to deliver to analog-digital converter after the signal of telecommunication amplification and become the real time spectrum signal.
Because the output light of light source 3 exportable multi-wavelengths guides optical fiber 4 to exciting, so the scattered light signal that measuring optical fiber 7 is received can be multi-wavelength.Therefore, under the irradiates light of single wavelength and the insignificant condition of fluorescence part in the scattered light signal, the beam splitting system of measuring instrument 8 can be made up of grating or prism, the scattered light signal of multi-wavelength is separated in different angles by its wavelength, become respectively by photodetector (being generally a line style photoelectronic detecting array) and to deliver to analog-digital converter after the signal of telecommunication amplifies and become the real time spectrum signal, be sent to computer 10 by cable 9 then.
Described control, date processing and calculating section are to be finished by computer.
As shown in Figure 3, the optical fiber assay method of mensuration turbid medium of the present invention and biological tissue optical parameter is by finishing as the next stage:
Phase I: the optical signalling on a plurality of wavelength is stored respectively by its wavelength from the real time spectrum signal of measuring instrument output;
Second stage: the real-time optical signal of the different wave length of storage is calculated the measured light signal after carrying out demodulation by Mathematical treatment such as Fourier transform etc. on the different wave length;
Phase III: enter the subprogram of optical transmission Theoretical Calculation, according to the incident illumination parameter of input, test sample probe form parameter and test sample optical parametric initial value obtain to calculate optical signal;
Quadravalence section: compare measured light signal and the difference of calculating optical signal, and by adjust repeatedly the test sample optical parametric until the difference of calculating optical signal and measured light signal less than till the setting value in advance, and export test sample optical parametric numerical value on the irradiates light wavelength.
Optical signalling on described a plurality of wavelength includes the incident illumination parameter, test sample probe form parameter and test sample optical parametric initial value.
Insert under the condition of test sample inside termination in measuring optical fiber 7, and phase III described optical transmission Theoretical Calculation subprogram can be finished by following steps:
S1: input incident illumination parameter and tracked total number of light photons N
0
S2: input test sample optics and boundary geometrical parameter;
S3: the incident direction by photon determines its initial direct of travel;
S4: set tracked number of photons N=1;
S5: determine the total distance of photon at random according to absorptance;
S7: judge that tracked number of photons N is greater than 1?, N enters S8 greater than 1, otherwise enters S9 according to S3;
S8: determine the photon scattering angle direction of also promptly advancing at random according to the scattering phase function, and enter S9;
S9: determine the free traveling distance of photon at random according to scattering coefficient;
S10: follow the trail of photon to next scattering point;
S11: judge that the accumulative total traveling distance is greater than total distance? be to enter S12, otherwise enter S13;
S12: judge after photon is absorbed to enter S6;
S6: will enter S5 after the tracked number of photons N increase by 1;
S13: judge whether Contact Boundary? be to enter S14, otherwise enter S8;
S14: judge whether to overflow test sample? be to enter S15, otherwise enter S16;
S15: judge after photon is escaped to enter S6;
S16: judge whether measured optical fiber acceptance? be to enter S17, otherwise enter S8;
S17: calculate optical signal;
S18: judge that tracked number of photons N is whether greater than the total number of light photons N of incident beam
0Be to enter S19, otherwise enter S6;
S19: optical signal is calculated in EP (end of program) and output.
Optical fiber assay method to mensuration turbid medium of the present invention and biological tissue optical parameter provides further instruction below.
The calculating of fibre system of the present invention and assay method determines that the core of program part is the optical signalling computational methods based on radiation transfer theory in the turbid media.Radiation transfer theory can be expressed as the integro-differential equation of a radiation transfer equation, and the form that itself and time have nothing to do can be expressed as follows
s·L(r，s)＝-(μ
a+μ
s)L(r，s)+μ
s∫
4πp(s，s′)L(r，s′)dΩ′。
S is the unit vector of optical propagation direction in the following formula, representative vector dot product operator, and representative vector gradient operator, r are the coordinate vector in the three dimensions, (r s) is luminous flux (luminous power in the unit are unit solid angle), μ to L
aBe absorptance, μ
sBe scattering coefficient, p (s, s ') is a scattering phase function (being proportional to light scatters to the s direction from s ' direction probability),
Representative is the three-dimensional viewpoin integration of the s ' direction of 4 π to the total solid angle of three bit spaces.Boundary value problem based on radiation transfer equation has two kinds of solutions usually: numerical solution method and be the statistical method of representative with the DSMC.The numerical solution method is found the solution according to boundary condition after aforesaid radiation transfer equation is converted into DIFFERENCE EQUATIONS.DSMC then is according to the described optical delivery process of aforesaid radiation transfer equation, represents incident beam with many photons, calculates the track of advancing of each photon transmission course in three dimensions.The track of advancing of photon determines by a plurality of stochastic variables, and the distribution function of these stochastic variables is respectively by absorptance, scattering coefficient and the decision of scattering phase function.Near the advance track of the photon zone boundary of being considered is usually according to the reflection coefficient formula manipulation of collimated light beam on boundary face.Calculate at the track of advancing to all photons (hundreds of thousands or more) and to carry out statistical analysis again after finishing, what those were detected photon that optical fiber collects and incident photon sum represents the calculating optical signal than promptly.
Fig. 4 is a kind of logic flow that realizes above-mentioned DSMC.This method is equivalent to a medium that comprises the light absorption center and the light scattering center of random distribution with turbid medium, the concentration at light absorption center and light scattering center and the absorptance of turbid medium, scattering coefficient is relevant respectively, and the random distribution at light absorption center and light scattering center then embodies by the random distribution to total distance of photon and free path.Before DSMC calculates beginning, the number of photons N that needs to import incident illumination parameter such as beam energy distribution and incident direction and represent incident beam
0, and optical parametric and the boundary geometry parameter of representing test sample.Because DSMC is a statistical method, its result can comprise statistical error, so need carry out the number of photons N of following calculation
0Must be enough big, could drop to the statistical error in the result of calculation enough little.But N
0Crossing conference causes computation time long.N generally speaking
0Between 9 powers of 4 powers to 10 10.
As shown in Figure 4, the Monte Carlo Calculation method need be to N
0It also is that photon or detected sample absorb or overflow test sample (also being that photon is escaped) or the measuring optical fiber that is inserted in the test sample is accepted until the end of advancing of this photon that individual incident photon carries out its traveling distance following calculation in test sample one by one.Before the traveling distance following calculation to each photon began, the Monte Carlo Calculation program will be according to total distance of being determined photon by the random distribution of test sample absorptance decision with according to the free traveling distance length of being determined photon by the random distribution of test sample scattering coefficient decision.The first step of photon following calculation supposes that for following the tracks of the initial direct of travel in photon edge to the position that its free traveling distance determined photon is scattered in this position.Begin next free path at photon and will whether be absorbed this photon before advancing, or overflow, or condition such as measured optical fiber acceptance is tested.Be satisfied as one of above-mentioned condition, then begin the traveling distance following calculation of next photon.All be not satisfied as above-mentioned condition, after the Monte Carlo Calculation program will also be the direction of next free path of advancing according to the scattering phase function of test sample (or at anisotropic parameters of determining under phase function form such as the Han Ni-Greenstein scattering phase function according to test sample) decision scattering angle, again according to determining the free path length of photon, thereby begin repeated calculation that the traveling distance of this photon is followed the tracks of until the photon end of advancing by the random distribution of test sample scattering coefficient decision.Accept i.e. conduct of this photon and the relevant calculating data record of calculating optical signal as the measured optical fiber of tracked photon (the measuring optical fiber bundle 7 among Fig. 1 or the optical fiber a to d among Fig. 2).After certain photon traveling distance following calculation is finished, the accumulative total N of the more tracked calculating photon of Monte Carlo Calculation program, as N greater than N
0, Monte Carlo Calculation finishes, otherwise N is increased by 1 back next incident photon is begun following calculation.When to N
0After the following calculation of individual incident photon is all finished, add up number of photons and N that all measured optical fiber are accepted
0Ratio promptly export from the Monte Carlo Calculation program as calculating optical signal.
Program is determined in the described calculating of realization Fig. 3, is to be finished by optical transmission Theoretical Calculation subprogram shown in Figure 4.Calculate and determine that the needed input data of program are: distribute power and area in incident illumination parameter such as the optical wavelength, beam cross section; The structure of test sample probe form parameter such as turbid media or biological tissue's test sample, longitudinal thickness, breadth wise dimension and refractive index and the diameter that excites guiding optical fiber and measuring optical fiber, paragraph header shape, numerical aperture, distance and refractive index between optical fiber; Test sample optical parametric initial value rule of thumb determines.Above-mentioned input data are imported by user interface by user usually.According to these input data, calculate the light Conveying Theory counting subroutine part of determining in the program and represent the photon of incident illumination energy collecting by measuring optical fiber through in the test sample model, being scattered the back with DSMC calculating by after exciting the outgoing of guiding optical fiber.Use the definition identical with the measured light signal, the aforementioned calculation subprogram is according to the number of photons of measuring optical fiber collection and the specific output and the measured light signal corresponding calculated optical signal of incident illumination subnumber.Determined calculating to determine whether program finishes or iterative computation and calculate with the difference of surveying optical signalling.As calculating the setting value that determines less than the experimental error according to the measured light signal with the difference of surveying optical signalling, then test sample optical parametric initial value is correct test sample optical parametric, calculates and determines EP (end of program) and export the test sample optical parametric.As calculating difference with the actual measurement optical signalling greater than setting value, calculate to determine that then program enters iteration cycle process, also promptly adjust the test sample optical parametric repeatedly and reenter light Conveying Theory counting subroutine and partly calculate optical signal until calculating difference with the actual measurement optical signalling less than setting value.
Test sample optical parametric modulation in the above-mentioned iteration cycle process can be based on following principle design.At first determine the modulation direction of attenuation quotient (for absorptance and scattering coefficient sum): as surveying the forward direction absorbance greater than calculating the forward direction absorbance according to the difference of the forward direction absorbance in actual measurement and the calculating optical signal, then reduce attenuation quotient, otherwise then increase attenuation quotient.According to the difference of diffuse transmittance in actual measurement and calculating optical signal and diffuse reflectance sum determine the modulation direction of absorptance: as surveying diffuse transmittance and diffuse reflectance sum greater than calculating diffuse transmittance and diffuse reflectance sum thereafter, then reduce absorptance, otherwise then increase absorptance.Determine the modulation direction of anisotropy coefficient according to the diffuse transmittance in actual measurement and the calculating optical signal and the difference of the ratio of diffuse reflectance then: as the ratio of actual measurement diffuse transmittance and diffuse reflectance greater than the ratio of calculating diffuse transmittance with diffuse reflectance, then increase anisotropy coefficient, otherwise then reduce anisotropy coefficient.
Claims (18)
1. a fibre system of measuring turbid medium and biological tissue optical parameter is characterized in that, includes: the light source importation; The photo measure part; Handle and calculating section with control data, wherein, the light source importation partly is connected with optical fiber measurement, the optical fiber measurement part is connected with control, date processing and calculating section, control, date processing and calculating section link to each other with the light source importation again, and the optical fiber measurement part also connects tested turbid medium or biological tissue.
2. the fibre system of measurement turbid medium according to claim 1 and biological tissue optical parameter, it is characterized in that, described light source importation includes light source electricity consumption source temperature controller (1) and light source (3), wherein, the input of power supply temperature controller (1) links to each other with control, date processing and calculating section by cable (11), outfan links to each other with light source (3) by power supply temperature control cable (2), and the output coupled end of light source (3) partly is connected with optical fiber measurement.
3. the fibre system of measurement turbid medium according to claim 1 and biological tissue optical parameter, it is characterized in that, described optical fiber measurement partly includes and excites guiding optical fiber (4), measuring instrument (8) and measuring optical fiber bundle (7), wherein, excite an end of guiding optical fiber (4) to be connected with the light source (3) of light source importation, measuring instrument (8) connects control, date processing and calculating section and measuring optical fiber bundle (7) respectively, excites the common tested turbid medium or the biological tissue of inserting of the other end of the other end and the measuring optical fiber bundle (7) of guiding optical fiber (4).
4. the fibre system of measurement turbid medium according to claim 1 and biological tissue optical parameter is characterized in that described control, date processing and calculating section are to be finished by computer.
5. the fibre system of measurement turbid medium according to claim 2 and biological tissue optical parameter is characterized in that, described light source (3) can be obtained by one or more laser coherence light source.
6. the fibre system of measurement turbid medium according to claim 2 and biological tissue optical parameter is characterized in that, described light source (3) also can be obtained by the incoherent light source of continuous spectrum, and its irradiates light wavelength is continuous distribution in the spectral domain of setting.
7. the fibre system of measurement turbid medium according to claim 2 and biological tissue optical parameter, it is characterized in that, the irradiates light intensity modulated of described light source (3) can by source current or mechanical chopper both one of, at 0.1 hertz to 100 mhz frequencies scope internal modulations.
8. the fibre system of measurement turbid medium according to claim 3 and biological tissue optical parameter is characterized in that, the described diameter of guiding optical fiber (4) that excites is between 10 microns～10 millimeters.
9. the fibre system of measurement turbid medium according to claim 3 and biological tissue optical parameter is characterized in that, the described planar normal in guiding optical fiber (4) termination that excites becomes the interior a certain angle of 0 degree～90 degree scopes with fiber axis.
10. the fibre system of measurement turbid medium according to claim 3 and biological tissue optical parameter, it is characterized in that, described measuring optical fiber bundle (7) is made up of 1～10 measuring optical fiber, and the diameter of the optical fiber in the measuring optical fiber bundle (7) is between 10 microns～10 millimeters.
11. the fibre system of measurement turbid medium according to claim 3 and biological tissue optical parameter is characterized in that, the planar normal in described measuring optical fiber bundle (7) termination is parallel with fiber axis or become 0 to spend the interior a certain angle of～90 degree scopes.
12. the fibre system of measurement turbid medium according to claim 3 and biological tissue optical parameter is characterized in that, described excite guiding optical fiber (4) and measuring optical fiber bundle (7) to be fixed in the same syringe needle after, in insertion turbid medium or the biological tissue.
13. the fibre system of measurement turbid medium according to claim 3 and biological tissue optical parameter, it is characterized in that described guiding optical fiber (4) and the measuring optical fiber bundle (7) of exciting adopts the fixing back of cable overcoat to enter the intracavity passage of human body by the service aisle of endoscope.
14. the fibre system of measurement turbid medium according to claim 3 and biological tissue optical parameter, it is characterized in that, described measuring instrument (8) is made up of beam split part, photodetection part and analog-digital converter, and the optical signal that its measuring optical fiber is exported is at first launched its wave spectrum by the beam split part and become by photodetector to deliver to analog-digital converter after the signal of telecommunication amplification and become the real time spectrum signal.
15. the fibre system of measurement turbid medium according to claim 3 and biological tissue optical parameter, it is characterized in that, described measuring instrument (8) photodetection part and analog-digital converter are formed, and the optical signals photodetector that measuring optical fiber is exported becomes to be delivered to analog-digital converter after the signal of telecommunication amplification and become the real-time optical signal.
16. an optical fiber assay method of measuring turbid medium and biological tissue optical parameter is characterized in that, is by finishing as the next stage:
Phase I: the optical signalling on a plurality of wavelength is stored respectively by its wavelength from the real time spectrum signal of measuring instrument output;
Second stage: the real-time optical signal of different wave length of storage is being undertaken calculating optical signal after the demodulation by Mathematical treatment on the different wave length;
Phase III: enter the subprogram of optical transmission Theoretical Calculation, according to the incident illumination parameter of input, test sample probe form parameter and test sample optical parametric initial value obtain to calculate optical signal;
Quadravalence section: compare measured light signal and the difference of calculating optical signal, and by adjust repeatedly the test sample optical parametric until the difference of calculating optical signal and measured light signal less than till the setting value in advance, and export test sample optical parametric numerical value on the irradiates light wavelength.
17. the optical fiber assay method of mensuration turbid medium according to claim 16 and biological tissue optical parameter, it is characterized in that, optical signalling on described a plurality of wavelength includes the incident illumination parameter, test sample probe form parameter and test sample optical parametric initial value.
18. the optical fiber assay method of mensuration turbid medium according to claim 16 and biological tissue optical parameter is characterized in that, the subprogram of described optical transmission Theoretical Calculation is by finishing as the next stage:
Phase I: the input optical parameter, determine tracked total number of light photons N
0, and determine photon total distance of advancing at random according to absorptance;
Second stage: judge that tracked accumulative total number of photons N is greater than 1? also promptly advance direction and determine the photon free path of advancing at random of decision photon scattering angle, and follow the trail of photon to next scattering point;
Phase III: judge that the accumulative total traveling distance is greater than total distance? whether is measured optical fiber accepted? overflow test sample? Contact Boundary whether?
Quadravalence section: calculate optical signal, and judge that accumulative total number of photons N is whether greater than the total number of light photons N of incident beam
0
Five-stage: optical signal is calculated in EP (end of program) and output.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2007100585937A CN100531669C (en) | 2007-08-06 | 2007-08-06 | Optical fiber system for detecting turbidity medium and biological tissue optical parameter |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2007100585937A CN100531669C (en) | 2007-08-06 | 2007-08-06 | Optical fiber system for detecting turbidity medium and biological tissue optical parameter |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101103905A true CN101103905A (en) | 2008-01-16 |
| CN100531669C CN100531669C (en) | 2009-08-26 |
Family
ID=38998135
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB2007100585937A Expired - Fee Related CN100531669C (en) | 2007-08-06 | 2007-08-06 | Optical fiber system for detecting turbidity medium and biological tissue optical parameter |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN100531669C (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102058393A (en) * | 2010-10-30 | 2011-05-18 | 华中科技大学 | Method and system for measuring kin physiology parameters and optical property parameters based on reflective spectral measurement |
| CN102292019A (en) * | 2009-01-23 | 2011-12-21 | 皇家飞利浦电子股份有限公司 | Optical examination device adapted to be at least partially inserted into a turbid medium |
| CN104603601A (en) * | 2012-08-16 | 2015-05-06 | 新加坡国立大学 | Diagnostic instrument and methods relating to raman spectroscopy |
| CN106018286A (en) * | 2016-05-18 | 2016-10-12 | 哈尔滨工业大学 | Reconstruction detector based on light-field camera and frequency modulation laser and applied to optical parameter distribution of dispersion media and reconstruction method |
| CN106706585A (en) * | 2016-12-31 | 2017-05-24 | 中国科学院昆明动物研究所 | Device and method for recording living body fluorescence signals |
| CN112384786A (en) * | 2018-05-09 | 2021-02-19 | 国立研究开发法人量子科学技术研究开发机构 | Tissue recognition device, tissue recognition system and method, examination method, endoscope system, tissue recognition program, and recording medium |
-
2007
- 2007-08-06 CN CNB2007100585937A patent/CN100531669C/en not_active Expired - Fee Related
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102292019A (en) * | 2009-01-23 | 2011-12-21 | 皇家飞利浦电子股份有限公司 | Optical examination device adapted to be at least partially inserted into a turbid medium |
| CN102292019B (en) * | 2009-01-23 | 2014-03-26 | 皇家飞利浦电子股份有限公司 | Optical examination device adapted to be at least partially inserted into a turbid medium |
| CN102058393A (en) * | 2010-10-30 | 2011-05-18 | 华中科技大学 | Method and system for measuring kin physiology parameters and optical property parameters based on reflective spectral measurement |
| CN102058393B (en) * | 2010-10-30 | 2012-10-31 | 华中科技大学 | Method for measuring kin physiology parameters and optical property parameters based on reflective spectral measurement |
| CN104603601A (en) * | 2012-08-16 | 2015-05-06 | 新加坡国立大学 | Diagnostic instrument and methods relating to raman spectroscopy |
| CN104603601B (en) * | 2012-08-16 | 2017-08-22 | 新加坡国立大学 | The diagnostic instrments relevant with Raman spectroscopy and method |
| CN106018286A (en) * | 2016-05-18 | 2016-10-12 | 哈尔滨工业大学 | Reconstruction detector based on light-field camera and frequency modulation laser and applied to optical parameter distribution of dispersion media and reconstruction method |
| CN106018286B (en) * | 2016-05-18 | 2018-12-11 | 哈尔滨工业大学 | Reconstruction detection device and method for reconstructing based on the distribution of the dispersive medium optical parameter of light-field camera and laser of frequency modulation |
| CN106706585A (en) * | 2016-12-31 | 2017-05-24 | 中国科学院昆明动物研究所 | Device and method for recording living body fluorescence signals |
| CN112384786A (en) * | 2018-05-09 | 2021-02-19 | 国立研究开发法人量子科学技术研究开发机构 | Tissue recognition device, tissue recognition system and method, examination method, endoscope system, tissue recognition program, and recording medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN100531669C (en) | 2009-08-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhu et al. | Review of Monte Carlo modeling of light transport in tissues | |
| CN100531669C (en) | Optical fiber system for detecting turbidity medium and biological tissue optical parameter | |
| Patterson et al. | Time resolved reflectance and transmittance for the noninvasive measurement of tissue optical properties | |
| Bonner et al. | Model for photon migration in turbid biological media | |
| CN102144154B (en) | Methods and systems for optically characterizing a turbid material using a structured incident beam | |
| US7440659B2 (en) | Depth-resolved reflectance instrument and method for its use | |
| JP4038179B2 (en) | Method and apparatus for the determination of light transport parameters and analytes in biological matrices | |
| CN101612034B (en) | Time-resolved measurement system and method for reconstructing optical parameters of turbid media | |
| US8263947B2 (en) | Method and device for localising fluorophores or absorbers in a surrounding medium | |
| JPH07209177A (en) | Method and apparatus for measuring scattering absorbing member | |
| CN101140222B (en) | Spectrometer system and method for measuring whole optical parameter including turbidity dielectric materials | |
| US9983127B2 (en) | Optical detection device and optical detection method | |
| CN111358473A (en) | Tissue blood flow blood oxygen imaging device and method based on near infrared spectrum | |
| US20020165456A1 (en) | Estimation of the average size of white light scatterers in normal and cancerous tissue using light scattering spectrum | |
| CN102495027B (en) | Method and device for measuring optical parameter of chaotic medium by multi-wavelength space resolution phase lock | |
| CN107300538A (en) | A kind of unicellular accurate test experience device based on optical fiber | |
| Sung et al. | Enhancing the sensitivity to scattering coefficient of the epithelium in a two-layered tissue model by oblique optical fibers: Monte Carlo study | |
| Wang et al. | Measurement of internal tissue optical properties at ultraviolet and visible wavelengths: Development and implementation of a fiberoptic-based system | |
| DePaoli et al. | Anisotropic light scattering from myelinated axons in the spinal cord | |
| JP5652599B2 (en) | Concentration determination apparatus, concentration determination method and program | |
| Gomes et al. | Algorithm for automated selection of application-specific fiber-optic reflectance probes | |
| Boros et al. | Optical and thermal modeling of an optrode microdevice for infrared neural stimulation | |
| WO2002040971A1 (en) | Method and apparatus for determination of optical parameters of turbid media and use thereof | |
| Lu et al. | Determination optical properties of tissue-like phantoms using diffuse reflectance and transmittance spectroscopy | |
| US20070232931A1 (en) | Optical Tomography System Using Short-Pulse Laser for Early Cancer Diagnostics |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| ASS | Succession or assignment of patent right |
Owner name: TIANJIN TIANYI-KANGDA CO. LTD. Free format text: FORMER OWNER: TIANJIN WEIFU MEDICAL TECHNOLOGY CO., LTD. Effective date: 20110218 |
|
| C41 | Transfer of patent application or patent right or utility model | ||
| COR | Change of bibliographic data |
Free format text: CORRECT: ADDRESS; FROM: 300204 2-6-507-510, RONGHUALI, GUANGDONG ROAD, HEXI DISTRICT, TIANJIN CITY TO: 300400 BESIDE QINJIANQIAO BRIDGE, BEICHEN DISTRICT, TIANJIN CITY |
|
| TR01 | Transfer of patent right |
Effective date of registration: 20110218 Address after: 300400, Beichen District, Tianjin province beside the bridge of diligence and frugality Patentee after: Tianjin Tianyi Kangda Photoelectric Technology Co., Ltd. Address before: 300204, Ronghua Road, Guangdong Road, Tianjin, Hexi District 2-6-507-510 Patentee before: Tianjin Weifu Medical Technology Co., Ltd. |
|
| DD01 | Delivery of document by public notice |
Addressee: Shang Hong Document name: Approval notice for cost mitigation |
|
| DD01 | Delivery of document by public notice |
Addressee: Tianjin Tianyi Kangda Photoelectric Technology Co., Ltd. Document name: Notification of Passing Examination on Formalities |
|
| DD01 | Delivery of document by public notice |
Addressee: Shang Hong Document name: Notification of Decision on Request for Restoration of Right |
|
| ASS | Succession or assignment of patent right |
Owner name: TIANJIN WEIFU MEDICAL TECHNOLOGY CO., LTD. Free format text: FORMER OWNER: TIANJIN TIANYI-KANGDA CO. LTD. Effective date: 20130416 |
|
| C41 | Transfer of patent application or patent right or utility model | ||
| COR | Change of bibliographic data |
Free format text: CORRECT: ADDRESS; FROM: 300400 BEICHEN, TIANJIN TO: 300457 HANGU, TIANJIN |
|
| TR01 | Transfer of patent right |
Effective date of registration: 20130416 Address after: 300457 experimental building 220, Dongting Road, Tianjin economic and Technological Development Zone, S615 Patentee after: Tianjin Wei Wei medical technology Co., Ltd. Address before: 300400, Beichen District, Tianjin province beside the bridge of diligence and frugality Patentee before: Tianjin Tianyi Kangda Photoelectric Technology Co., Ltd. |
|
| CI01 | Correction of invention patent gazette |
Correction item: Patentee Correct: Tianjin Weifu Medical Technology Co., Ltd.|300457 experimental building 220, Dongting Road, Tianjin economic and Technological Development Zone, S615 False: Tianjin Wei Wei medical technology Co., Ltd.|300457 experimental building 220, Dongting Road, Tianjin economic and Technological Development Zone, S615 Number: 19 Volume: 29 |
|
| ERR | Gazette correction | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090826 Termination date: 20190806 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |