CN103799984A - Method for obtaining depth of internal heat source of biological tissue on basis of infrared thermal imaging - Google Patents
Method for obtaining depth of internal heat source of biological tissue on basis of infrared thermal imaging Download PDFInfo
- Publication number
- CN103799984A CN103799984A CN201410067769.5A CN201410067769A CN103799984A CN 103799984 A CN103799984 A CN 103799984A CN 201410067769 A CN201410067769 A CN 201410067769A CN 103799984 A CN103799984 A CN 103799984A
- Authority
- CN
- China
- Prior art keywords
- temperature
- gray
- pixel
- depth
- heat source
- 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
Images
Abstract
The invention provides a method for obtaining the depth of an internal heat source of biological tissue on the basis of infrared thermal imaging. According to the method, a relational expression among the surface temperature difference, the depth of the internal heat source and the strength of the internal heat source is obtained on the basis of a classic biological heat transfer theory, a grayscale image is extracted through an infrared video camera, the grayscale of each pixel is extracted, grayscale segmentation is carried out, distribution of temperature points is detected, data are substituted into the relational expression among the surface temperature difference, the depth of the internal heat source and the strength of the internal heat source, and the depth of an abnormal temperature point is obtained. According to the method for obtaining the depth of the internal heat source of the biological tissue on the basis of infrared thermal imaging, the depth of the internal heat source of the biological tissue can be easily obtained only according to the environment object temperature, the black body temperature, grayscale distribution and the like.
Description
Technical field
The invention belongs to infrared thermal imaging technique field, be specifically related to a kind of biological tissue's endogenous pyrogen degree of depth acquisition methods based on infrared thermal imaging.
Background technology
The infrared thermal imagery that existing Medical Infrared Thermal Imaging Technology distributes by human surface temperature can be diagnosed its inside pathological changes qualitatively, but can not obtain the variation of organization internal temperature and the relevant information of focus, cannot position focus.
1997, develop the first in the world platform Thermal texture maps (Thermal Texture Maps System), claim again TTM system.The principle of this system is to utilize infrared emanation receiver-scanner to receive human infrared radiation, reconstruct the relative metabolism intensity distribution of cell corresponding to human body institute check point based on specific rule and algorithm, and tomography in addition, measure the degree of depth of infrared source, according to the poor health status of assessing human body of heat radiation in normal and abnormal structure region.But TTM technology is that tissue is considered as homogeneous structure, adopt that heat--the method for electrical analogue finds the relation of human body surface heat distribution and heat depth, the relation of so-called Gaussian curve half power points obtains the information such as the degree of depth of internal heat resource.Can find out, this theoretical foundation has too been simplified with respect to complicated tissue, does not consider the impact of some physical parameters of organizing self, and use procedure complexity, produces error large, awaits further development and improvement.
Summary of the invention
The present invention proposes a kind of biological tissue's endogenous pyrogen degree of depth acquisition methods based on infrared thermal imaging, can only utilize environment object temperature, blackbody temperature, and intensity profile etc. are obtained the degree of depth of biological tissue's endogenous pyrogen simplely.
In order to solve the problems of the technologies described above, the invention provides a kind of biological tissue's endogenous pyrogen degree of depth acquisition methods based on infrared thermal imaging, comprise the following steps:
The first step, at floor, a room temperature object close with biological in vitro detecting object color is set, and with the black matrix of assigned temperature as room temperature object background, use infrared video camera that detecting object and room temperature object are taken in same gray level image, obtain gray level image.Natural temperature and the blackbody temperature at records photographing scene simultaneously;
The gray value of each pixel in second step, the gray level image that obtains according to step 1, obtain the coordinate of the pixel that in gray-scale map, gray value is the highest, then obtain in gray-scale map the pixel the highest with gray value and have the coordinate of each pixel of same abscissa according to the coordinate of the highest pixel of gray value;
The 3rd step, detecting object gray level image from gray-scale map intercept one and cross full scale image being detected on object labelling in advance.Physical length size corresponding to each pixel in simultaneously calculating gray level image according to this object quantity of pixel and proportionate relationship of its length and width on the longer wide-edge of gray-scale map;
The gray value differences of the 4th step, calculating background blackbody temperature and room temperature object, and calculating utilizes the temperature difference of blackbody temperature and natural temperature, calculate according to temperature difference the corresponding proportionate relationship with gray value differences, each gray scale is carried out to temperature division, obtain the temperature of the highest gray scale point and its all gray scale points on same cross shaft surface;
The 5th step, the temperature substitution formula (1) of arbitrfary point in truncated picture in step 3 is solved to the endogenous pyrogen degree of depth of detecting object,
T
ρ＝T
0+kd(T
m-T
0)/[hρ(ρ-d)+kd] （1）
In formula (1), the endogenous pyrogen degree of depth that d is detecting object, T
ρfor arbitrfary point, detecting object surface temperature, T
mfor detecting object hot face temperature, ρ is the distance of arbitrfary point apart from initial point O, ρ=(d
2+ a
2)
1/2, k is utilizing object conducting heat coefficient, h is the convection transfer rate of testee and air.
Compared with prior art, its remarkable advantage is in the present invention, and (1) calculates effectively simple, on the physical model of setting up, in the surface temperature of deriving and heat depth relational expression, only have a degree of depth unknown quantity, make to calculate simply, directly can in substitution relational expression, calculate.(2) concrete operations are simple, for the Scaling Problem of infrared video camera, utilize pixel, gray scale, and the relation of indoor temperature can easily solve, and has saved and has calibrated this loaded down with trivial details problem, can directly use fast.
Accompanying drawing explanation
Fig. 1 is the inventive method flow chart.
Fig. 2 detecting object Hepar Sus domestica and background gray-scale map.
Fig. 3 is the in vitro detecting object gray-scale map of biology that during the present invention tests, step 1 obtains.
Fig. 4 is step 1 environment room temperature object gray-scale map during the present invention tests.
Fig. 5 is that during the present invention tests, 40 degrees Celsius of backgrounds of step 1 black matrix manufacture refer to gray-scale map.
The specific embodiment
As shown in Figure 1, the present invention is based on biological tissue's endogenous pyrogen degree of depth acquisition methods of infrared thermal imaging, step is as follows:
The first step, at floor, a room temperature object close with biological in vitro detecting object color is set, and with the black matrix of assigned temperature as room temperature object background, use infrared video camera that detecting object and room temperature object are taken in same gray level image, obtain gray level image.Natural temperature and the blackbody temperature at records photographing scene simultaneously.
The gray value of each pixel in second step, the gray level image that obtains according to step 1, obtain the coordinate of the pixel that in gray-scale map, gray value is the highest, then obtain in gray-scale map the pixel the highest with gray value and have the coordinate of each pixel of same abscissa according to the coordinate of the highest pixel of gray value.
The 3rd step, detecting object gray level image from gray-scale map intercept one and cross full scale image being detected on object labelling in advance.Physical length size corresponding to each pixel in simultaneously calculating gray level image according to this object quantity of pixel and proportionate relationship of its length and width on the longer wide-edge of gray-scale map.
The gray value differences of the 4th step, calculating background blackbody temperature and room temperature object, and calculating utilizes the temperature difference of blackbody temperature and natural temperature, calculate according to temperature difference the corresponding proportionate relationship with gray value differences, each gray scale is carried out to temperature division, obtain the temperature of the highest gray scale point and its all gray scale points on same cross shaft surface.
The 5th step, the temperature substitution formula (1) of arbitrfary point in truncated picture in step 3 is solved to the endogenous pyrogen degree of depth of detecting object.
T
ρ＝T
0+kd(T
m-T
0)/[hρ(ρ-d)+kd] （1）
In formula (1), the endogenous pyrogen degree of depth that d is detecting object, T
ρfor arbitrfary point, detecting object surface temperature, T
mfor detecting object hot face temperature, ρ is the distance of arbitrfary point apart from initial point O, ρ=(d
2+ a
2)
1/2, k is utilizing object conducting heat coefficient, h is the convection transfer rate of testee and air.
The derivation of formula (1) refers to " based on the Medical Infrared Thermal Imaging Technology of biology heat transfer effect " (Zhou Minhua, Institutes Of Technology Of Nanjing's doctorate paper, 2010)
The present invention can further illustrate by following experiment:
At temperature stabilization in the room of 27 ℃, get a fresh pig liver, use the suitable resistance of size as detecting object endogenous pyrogen, resistance temperature maintained and stablely use thermal camera to take after forming a constant temperature field, obtain Hepar Sus domestica hot face temperature with and position.Endogenous pyrogen degree of depth 0.75CM is set, and endogenous pyrogen temperature is 40 degrees Celsius, and setting blackbody temperature is 40 degrees Celsius, the region that size of labelling is 12cm × 10cm on Hepar Sus domestica.
At floor, a plastic casing close with Hepar Sus domestica color is set according to the inventive method step 1, and with specifying the black matrix of 40 ℃ as room temperature object background, use infrared video camera that Hepar Sus domestica and plastic casing are taken in same gray level image, obtain gray level image, record as shown in Figure 2 40 ℃ of 27 ℃ of on-the-spot natural temperatures and blackbody temperatures simultaneously.
The gray value of each pixel in the gray level image obtaining according to step 1, obtains the coordinate of the pixel that in Fig. 2, gray value is the highest for (85,1), as shown in table 1; Then obtain in gray-scale map the pixel the highest with gray value and have the coordinate of each pixel of same abscissa according to the coordinate of the highest pixel of gray value.
Table 1 gray scale and pixel distribution
Use MATLAB software to obtain on gray level image the gray scale of other points on the transverse axis coordinate of the highest gray scale point place, as shown in table 2
Gray value corresponding to pixel on the transverse axis coordinate that table 2 is sitting at the highest gray scale point
| Pixel | 87 | 88 | 89 | 90 | 91 | ......... | 123 | 124 | 125 |
| Gray value | 65 | 64 | 63 | 63 | 63 | ......... | 59 | 59 | 59 |
From gray-scale map 2, obtain the gray-scale map of Hepar Sus domestica marked area according to step 3 as shown in Figure 3, this marked area is the image for pixel 125 × 100.The length scale that just can obtain each pixel is 0.1CM.
The gray scale that obtains 40 ℃ of background black matrixes according to step 4 is 199, and the gray scale that room temperature plastic casing is 27 ℃ is 45.The temperature difference that marks off the interval representative of each gray value is 0.0844 degree Celsius, can obtain the temperature of the highest gray scale point and its all gray scale points on same cross shaft surface.
Can calculate the each point of tissue and the distance of maximum temperature point and the temperature level of each point according to gray value corresponding to pixel.Apart from d: each point is poor with the pixel of the highest gray value.Each temperature T
ρ: the long-pending indoor temperature that adds of the difference of the gray value of gray value and indoor temperature and gray value corresponding temperature value.
According to depth calculation formula T
ρ=T
0+ kd (T
m-T
0)/[h
cρ (ρ-d)+kd], substitution data, obtain depth d.
As shown in table 3:
The corresponding relation of table 3 tissue surface temperature and distance
Experiment finishes the depth d=0.5833cm calculating.Error is 7.4%, and error amount is less than 10%, has illustrated that relative gray scale draws theoretical feasibility and the correctness of obtaining internal heat source information method at infrared thermal imaging.
Claims (1)
1. the biological tissue's endogenous pyrogen degree of depth acquisition methods based on infrared thermal imaging, is characterized in that, comprises the following steps:
The first step, at floor, a room temperature object close with biological in vitro detecting object color is set, and with the black matrix of assigned temperature as room temperature object background, use infrared video camera that detecting object and room temperature object are taken in same gray level image, obtain gray level image.Natural temperature and the blackbody temperature at records photographing scene simultaneously;
The gray value of each pixel in second step, the gray level image that obtains according to step 1, obtain the coordinate of the pixel that in gray-scale map, gray value is the highest, then obtain in gray-scale map the pixel the highest with gray value and have the coordinate of each pixel of same abscissa according to the coordinate of the highest pixel of gray value;
The 3rd step, detecting object gray level image from gray-scale map intercept one and cross full scale image being detected on object labelling in advance.Physical length size corresponding to each pixel in simultaneously calculating gray level image according to this object quantity of pixel and proportionate relationship of its length and width on the longer wide-edge of gray-scale map;
The gray value differences of the 4th step, calculating background blackbody temperature and room temperature object, and calculating utilizes the temperature difference of blackbody temperature and natural temperature, calculate according to temperature difference the corresponding proportionate relationship with gray value differences, each gray scale is carried out to temperature division, obtain the temperature of the highest gray scale point and its all gray scale points on same cross shaft surface;
The 5th step, the temperature substitution formula (1) of arbitrfary point in truncated picture in step 3 is solved to the endogenous pyrogen degree of depth of detecting object,
T
ρ＝T
0+kd(T
m-T
0)/[hρ(ρ-d)+kd] （1）
In formula (1), the endogenous pyrogen degree of depth that d is detecting object, T
ρfor arbitrfary point, detecting object surface temperature, T
mfor detecting object hot face temperature, ρ is the distance of arbitrfary point apart from initial point O, ρ=(d
2+ a
2)
1/2, k is utilizing object conducting heat coefficient, h is the convection transfer rate of testee and air.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410067769.5A CN103799984B (en) | 2014-02-26 | 2014-02-26 | Based on biological tissue's endogenous pyrogen degree of depth acquisition methods of infrared thermal imaging |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410067769.5A CN103799984B (en) | 2014-02-26 | 2014-02-26 | Based on biological tissue's endogenous pyrogen degree of depth acquisition methods of infrared thermal imaging |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103799984A true CN103799984A (en) | 2014-05-21 |
| CN103799984B CN103799984B (en) | 2016-04-20 |
Family
ID=50697727
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410067769.5A Expired - Fee Related CN103799984B (en) | 2014-02-26 | 2014-02-26 | Based on biological tissue's endogenous pyrogen degree of depth acquisition methods of infrared thermal imaging |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103799984B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105758400A (en) * | 2016-02-15 | 2016-07-13 | 上海卫星工程研究所 | Extracting method for fixed star sensitive east and west parameters of stationary satellite imaging navigation and registration |
| CN107802243A (en) * | 2017-09-29 | 2018-03-16 | 武汉昊博科技有限公司 | The analysis method and device of organism internal heat resource intensity and Temperature Distribution based on heat point source model |
| CN108606782A (en) * | 2018-04-28 | 2018-10-02 | 泰州市榕兴医疗用品股份有限公司 | A kind of surface of a wound imaging system |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000166950A (en) * | 1998-12-02 | 2000-06-20 | Lion Corp | Periodontal pocket detector |
| US6381488B1 (en) * | 1999-06-15 | 2002-04-30 | Sandia Corporation | Method and apparatus to measure the depth of skin burns |
| CN1771882A (en) * | 2005-09-02 | 2006-05-17 | 武汉市昊博科技有限公司 | Method and device for obtaining internal heat source information from the surface temperature distribution of living body |
| CN101088454A (en) * | 2007-07-11 | 2007-12-19 | 武汉昊博科技有限公司 | Medical thermal chromatographic imaging system |
-
2014
- 2014-02-26 CN CN201410067769.5A patent/CN103799984B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000166950A (en) * | 1998-12-02 | 2000-06-20 | Lion Corp | Periodontal pocket detector |
| US6381488B1 (en) * | 1999-06-15 | 2002-04-30 | Sandia Corporation | Method and apparatus to measure the depth of skin burns |
| CN1771882A (en) * | 2005-09-02 | 2006-05-17 | 武汉市昊博科技有限公司 | Method and device for obtaining internal heat source information from the surface temperature distribution of living body |
| CN101088454A (en) * | 2007-07-11 | 2007-12-19 | 武汉昊博科技有限公司 | Medical thermal chromatographic imaging system |
Non-Patent Citations (2)
| Title |
|---|
| 周敏华: "生物组织中异常热源信息获取方法", 《南京理工大学学报(自然科学版)》, vol. 33, no. 3, 30 June 2009 (2009-06-30), pages 379 - 387 * |
| 王春燕: "体内异常热源参数与体表温度关系的热像研究", 《激光与红外》, vol. 42, no. 1, 31 January 2012 (2012-01-31), pages 31 - 35 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105758400A (en) * | 2016-02-15 | 2016-07-13 | 上海卫星工程研究所 | Extracting method for fixed star sensitive east and west parameters of stationary satellite imaging navigation and registration |
| CN105758400B (en) * | 2016-02-15 | 2018-08-03 | 上海卫星工程研究所 | Fixed statellite imaging based navigation be registrated fixed star sensitivity thing parameter extracting method |
| CN107802243A (en) * | 2017-09-29 | 2018-03-16 | 武汉昊博科技有限公司 | The analysis method and device of organism internal heat resource intensity and Temperature Distribution based on heat point source model |
| CN107802243B (en) * | 2017-09-29 | 2020-11-03 | 武汉昊博科技有限公司 | Method and device for analyzing heat source intensity and temperature distribution in organism based on point heat source model |
| CN108606782A (en) * | 2018-04-28 | 2018-10-02 | 泰州市榕兴医疗用品股份有限公司 | A kind of surface of a wound imaging system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103799984B (en) | 2016-04-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10401156B2 (en) | System and method for quantifying deformation, disruption, and development in a sample | |
| CN107507195B (en) | The multi-modal nasopharyngeal carcinoma image partition method of PET-CT based on hypergraph model | |
| CN105608737B (en) | A kind of human foot three-dimensional rebuilding method based on machine learning | |
| CN107677372A (en) | A kind of tunnel detection method based on binocular vision | |
| CN103799984B (en) | Based on biological tissue's endogenous pyrogen degree of depth acquisition methods of infrared thermal imaging | |
| CN105628208B (en) | A kind of thermometry based on infrared imaging system | |
| GB2551594A (en) | Surface simulation | |
| CN102779354A (en) | Three-dimensional reconstruction method for traditional Chinese medicine inspection information surface based on photometric stereo technology | |
| CN108805939A (en) | The caliberating device and method of trinocular vision system based on statistics feature | |
| CN110060304B (en) | Method for acquiring three-dimensional information of organism | |
| CN109480788A (en) | A kind of human skin epidermis and subdermal muscle layer three-dimensional infrared thermal imaging testing method | |
| Pesce et al. | A 12-camera body scanning system based on close-range photogrammetry for precise applications | |
| CN105521562B (en) | A kind of temperature field indication of tumor thermotherapy and control device and method | |
| CN104981709B (en) | It is imaged thermal measurement | |
| CN115523847A (en) | Monocular camera ranging method and system | |
| CN103955687B (en) | A kind of method for rapidly positioning of the light spot image center based on centroid method | |
| Krefer et al. | A method for generating 3D thermal models with decoupled acquisition | |
| IL258134A (en) | Method and system for correcting image data | |
| CN112633113A (en) | Cross-camera human face living body detection method and system | |
| CN111486961A (en) | Efficient forehead temperature estimation method based on wide-spectrum human forehead imaging and distance sensing | |
| CN103335727A (en) | Thermal imaging image processing method based on setting of multiple emissivities for visible light divided area | |
| CN107644679B (en) | Information pushing method and device | |
| Nasehi Tehrani et al. | Lung surface deformation prediction from spirometry measurement and chest wall surface motion | |
| Hess et al. | Multimodal registration of high-resolution thermal image mosaics for the non-destructive evaluation of structures | |
| CN205850007U (en) | The temperature field indication of a kind of tumor thermotherapy and control device |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160420 Termination date: 20180226 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |