CN110298902B - Medical infrared image reconstruction method, system, computer and storage medium - Google Patents

Medical infrared image reconstruction method, system, computer and storage medium Download PDF

Info

Publication number
CN110298902B
CN110298902B CN201910689274.9A CN201910689274A CN110298902B CN 110298902 B CN110298902 B CN 110298902B CN 201910689274 A CN201910689274 A CN 201910689274A CN 110298902 B CN110298902 B CN 110298902B
Authority
CN
China
Prior art keywords
image
pixels
pixel
activity ratio
temperature
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.)
Active
Application number
CN201910689274.9A
Other languages
Chinese (zh)
Other versions
CN110298902A (en
Inventor
陈耀弘
武力
王华伟
谢庆胜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xi'an Zhongke Feitu Photoelectric Technology Co ltd
Original Assignee
Xi'an Zhongke Feitu Photoelectric Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xi'an Zhongke Feitu Photoelectric Technology Co ltd filed Critical Xi'an Zhongke Feitu Photoelectric Technology Co ltd
Priority to CN201910689274.9A priority Critical patent/CN110298902B/en
Publication of CN110298902A publication Critical patent/CN110298902A/en
Application granted granted Critical
Publication of CN110298902B publication Critical patent/CN110298902B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/001Texturing; Colouring; Generation of texture or colour
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10048Infrared image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30096Tumor; Lesion
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30101Blood vessel; Artery; Vein; Vascular

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

The invention belongs to a medical infrared image processing method, and particularly relates to a medical infrared image reconstruction method, a system, a computer and a storage medium. The method solves the problems of insufficient accuracy caused by insufficient contrast of the infrared image and relatively blurred edge profile after reconstruction in the existing infrared image processing. The reconstruction method comprises the steps of collecting a stable preprocessed image, dividing the preprocessed image into a plurality of pixels, obtaining the activity ratio of each pixel according to the temperature recovery rate and the refrigeration rate of each pixel, and obtaining a corresponding activity ratio image through color mapping to finish reconstruction; the reconstruction system comprises a thermal infrared imager, an image processing module and a display module, and is used for realizing the method. In addition, the reconstruction method of the present invention can also be stored on a computer readable storage medium or on a processor of a computer for execution. An image with better contrast and definition can be obtained.

Description

Medical infrared image reconstruction method, system, computer and storage medium
Technical Field
The invention belongs to a medical infrared image processing method, and particularly relates to a medical infrared image reconstruction method, a system, a computer and a storage medium.
Background
According to clinical studies, when some parts of the human body are diseased, the change of the heat energy distribution is generated at the same time, so that the distribution change of the body surface temperature is caused, therefore, the temperature is one of the most commonly used indexes for observing and measuring whether the human body is normal or not, and the diagnosis by acquiring the distribution change of the whole body or local temperature is also a conventional medical diagnosis means.
The medical infrared thermal imaging system can acquire real-time temperature distribution of the body surface of the human body, compares the infrared influence under the condition of human illness with the infrared influence under the normal physiological state, judges the pathological state through the difference, has poor contrast ratio in the traditional medical infrared thermal imaging, is difficult to accurately carry out qualitative analysis on lesions such as subcutaneous blood vessels and tumors, and limits the clinical application and popularization of the lesions.
It has been confirmed that the skin surface is subjected to external contact stimulation (most commonly cold stimulation with an ice pack), and the recovery rate of the subcutaneous vascular and tumor tissue is significantly faster than that of the adjacent areas during the recovery of the skin temperature to the normal state after the stimulus source is removed. Based on the method, a dynamic continuous image sequence with gradual change of the skin temperature field is obtained by using an infrared thermal imager after cold stimulation, the state of the skin in the temperature recovery process can be represented, an infrared image can be obtained after image reconstruction is carried out on the images, and the morphological characteristics of the target can be clearly reflected. The dynamic continuous image sequence is analyzed by taking the pixels as units, the most commonly used reconstruction method at present is to generate a tau image, and define the time required for the temperature recovery to reach a stable state (the numerical value of the pixels in the dynamic continuous image sequence is not changed any more) after the cold stimulus of the ith row and the jth row of pixels is removed as t ij Each pixel value tau in the tau image ij The calculation mode of (a) is as follows: τ ij =0.63t ij
However, although the τ image can more clearly represent the difference between the affected part and the healthy part than the conventional infrared image, the edge contour of the τ image is blurred, and the accuracy thereof is still to be improved.
Disclosure of Invention
The invention mainly aims to solve the problems of insufficient accuracy caused by the fact that the reconstructed edge contour is still blurred due to insufficient contrast of an infrared image in the existing infrared image processing. A medical infrared image reconstruction method, system, computer and storage medium are provided.
In order to achieve the above purpose, the present invention provides the following technical solutions:
the medical infrared image reconstruction method is used for non-disease diagnosis and is characterized by comprising the following steps of:
step 1, collecting infrared images
Continuously acquiring the image of the part to be detected after cold treatment by using a thermal infrared imager until the image is no longer changed, stopping acquisition, and recording as a preprocessed image;
step 2, dividing the pixels
Dividing the preprocessed image into i x j pixels, wherein i is the number of rows of the pixels, and j is the number of columns of the pixels;
step 3, reconstructing an image
Calculating the temperature recovery rate of each pixel:
Figure BDA0002147368610000021
the refrigeration rate of each pixel is calculated:
Figure BDA0002147368610000022
wherein t is 1 For cold treatment time, t 2 For the duration of the acquired image; t (T) 1ij T is the temperature before cold treatment of the ith row and jth column pixels 2ij The temperature after cold treatment of the ith row and the jth column of pixels; the highest temperature value of each pixel in the process of acquiring the infrared image in the step 1 is recorded as T maxij Reference temperature T R For each picture element T maxij Is the minimum of (2);
obtaining the activity ratio A of each pixel ij
Figure BDA0002147368610000023
And (3) mapping the activity ratio of each pixel to obtain an activity ratio image through color mapping, and completing reconstruction.
In step 1, the cold treatment is specifically that the part to be measured is exposed for at least 2min under the constant temperature and humidity environment with the temperature range of 24+/-2 ℃ and the humidity range of 40-60%, and then is cold-rolled by an ice bag, and the ice bag is filled with an ice-water mixture.
Further, in the step 1, the continuous acquisition of the images of the part to be detected which is subjected to the cold treatment by using the thermal infrared imager is performed under the constant temperature and humidity environment with the temperature range of 24+/-2 ℃ and the humidity range of 40-60%.
Further, in step 1, the thermal infrared imager is used for continuously collecting images of the part to be measured after the cold treatment, and the frequency of continuous collection is 0.5 s/sheet.
Further, in step 1 and step 2, the t 1 And t 2 Is accurate in seconds.
Further, in step 2, the preprocessed image is divided into i×j pixels, where i is 240 and j is 320.
The reconstruction system for realizing the medical infrared image reconstruction method is characterized by comprising a thermal infrared imager, an image processing module and a display module;
the infrared thermal imager continuously collects the images of the part to be detected after the cold treatment until the images are no longer changed and recorded as preprocessed images, and the collection is stopped;
the image processing module is used for dividing the preprocessed image into i x j pixels, wherein i is the number of rows of the pixels, and j is the number of columns of the pixels; according to the refrigerating rate C of each pixel cold treatment ij And a temperature recovery rate R during continuous image acquisition ij Obtaining the activity ratio of each pixel
Figure BDA0002147368610000031
Performing color mapping on the activity ratio of each pixel to obtain an activity ratio image, and sending the activity ratio image to a display module;
and the display module displays the activity ratio image.
A computer readable storage medium having stored thereon a computer program, characterized in that the program when executed by a processor performs the steps of:
step 1, dividing pixels
Dividing the preprocessed image into i x j pixels, wherein i is the number of rows of the pixels, and j is the number of columns of the pixels; the preprocessing image is obtained by continuously collecting the image of the part to be detected after cold processing by using a thermal infrared imager until the image is not changed any more, stopping collecting, and recording the image as the preprocessing image;
step 2, reconstructing an image
Calculating the temperature recovery rate of each pixel:
Figure BDA0002147368610000032
the refrigeration rate of each pixel is calculated:
Figure BDA0002147368610000041
wherein t is 1 For cold treatment time, t 2 For the duration of the acquired image; t (T) 1ij T is the temperature before cold treatment of the ith row and jth column pixels 2ij The temperature after cold treatment of the ith row and the jth column of pixels; the highest temperature value of each pixel in the process of acquiring the infrared image in the step 1 is recorded as T maxij Reference temperature T R For each picture element T maxij Is the minimum of (2);
obtaining the activity ratio A of each pixel ij
Figure BDA0002147368610000042
And (3) mapping the activity ratio of each pixel to obtain an activity ratio image through color mapping, and completing reconstruction.
A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the following steps when executing the program:
step 1, dividing pixels
Dividing the preprocessed image into i x j pixels, wherein i is the number of rows of the pixels, and j is the number of columns of the pixels; the preprocessing image is obtained by continuously collecting the image of the part to be detected after cold processing by using a thermal infrared imager until the image is not changed any more, stopping collecting, and recording the image as the preprocessing image;
step 2, reconstructing an image
Calculating the temperature recovery rate of each pixel:
Figure BDA0002147368610000043
the refrigeration rate of each pixel is calculated:
Figure BDA0002147368610000044
wherein t is 1 For cold treatment time, t 2 For the duration of the acquired image; t (T) 1ij T is the temperature before cold treatment of the ith row and jth column pixels 2ij The temperature after cold treatment of the ith row and the jth column of pixels; the highest temperature value of each pixel in the process of acquiring the infrared image in the step 1 is recorded as T maxij Reference temperature T R For each picture element T maxij Is the minimum of (2);
obtaining the activity ratio A of each pixel ij
Figure BDA0002147368610000051
And (3) mapping the activity ratio of each pixel to obtain an activity ratio image through color mapping, and completing reconstruction.
Compared with the prior art, the invention has the beneficial effects that:
1. according to the medical infrared image reconstruction method, after cold treatment of the part to be detected, a dynamic continuous image sequence with gradual change of a skin temperature field is obtained by using an infrared thermal imager, the state of skin in the temperature recovery process is represented, a corresponding activity ratio image can be obtained after the image sequence is reconstructed, the activity ratio image is better than the original infrared image and tau image in contrast ratio by introducing the temperature recovery rate and the refrigeration rate, morphological characteristics of the part to be detected can be reflected more clearly, the image definition and contrast pair are better, and details of target characteristics are more clear. The method of the invention can be widely applied to simulation or teaching.
2. The invention fully exposes the part to be measured under the condition of constant temperature and constant humidity, thereby eliminating the environmental interference.
3. The invention uses the thermal infrared imager to collect images under the condition of constant temperature and humidity, the accuracy is higher, and the thermal infrared imager is less influenced by environmental change.
4. The frequency of the acquired image is 0.5 s/sheet, especially when the part to be measured is about to recover after cold treatmentWhen in use, the preprocessing image can be more accurately captured, and the acquisition is stopped in time, so that the duration t of the acquired image is further prolonged 2 Is more accurate.
5. The time accuracy recorded in the invention is s, and the reconstruction structure is more accurate.
6. The invention divides the preprocessed image into 240 x 320 pixels, thereby ensuring the diagnosis precision and reasonably saving the reconstruction time due to the division number.
7. According to the medical infrared image reconstruction system, images are continuously acquired through the thermal infrared imager until a preprocessed image is obtained in a stable state, an image processing module is used for dividing, calculating and mapping colors of the preprocessed image, and finally an activity ratio image is obtained on a display module. The reconstruction method can be effectively realized by the reconstruction system of the invention, and the active ratio image with higher contrast and contour definition is obtained.
8. The computer readable storage medium of the present invention stores the method of image reconstruction, converts the method of reconstruction into a program executable by a processor, and completes the reconstruction in a more convenient method.
9. The computer equipment stores the reconstruction method on the memory in a mode of an executable program, and can automatically complete the reconstruction by operating the computer equipment through the processor and only inputting a preprocessed image in advance.
Drawings
FIG. 1 is a flow chart of an infrared image reconstruction method according to the present invention;
FIG. 2 is a diagram showing images of a part to be measured after being subjected to cold treatment continuously acquired by a thermal infrared imager in accordance with a first embodiment of the present invention;
FIG. 3 is an activity ratio image obtained in the first embodiment of the present invention;
FIG. 4 is a τ image obtained in the first embodiment of the present invention
Fig. 5 is a τ' image obtained in the first embodiment of the present invention.
Detailed Description
The technical solutions of the present invention will be clearly and completely described below with reference to the embodiments of the present invention and the accompanying drawings, and it is apparent that the described embodiments do not limit the present invention.
An infrared image reconstruction method as shown in fig. 1, comprising the steps of:
step 1, collecting infrared images
Continuously acquiring the image of the part to be detected after cold treatment by using a thermal infrared imager until the image is no longer changed, stopping acquisition, and recording as a preprocessed image;
step 2, dividing the pixels
Dividing the preprocessed image into i x j pixels, wherein i is the number of rows of the pixels, and j is the number of columns of the pixels;
step 3, reconstructing an image
Calculating the temperature recovery rate of each pixel:
Figure BDA0002147368610000061
the refrigeration rate of each pixel is calculated:
Figure BDA0002147368610000062
wherein t is 1 For cold treatment time, t 2 For the duration of the acquired image; t (T) 1ij T is the temperature before cold treatment of the ith row and jth column pixels 2ij The temperature after cold treatment of the ith row and the jth column of pixels; the highest temperature value of each pixel in the process of acquiring the infrared image in the step 1 is recorded as T maxij Reference temperature T R For each picture element T maxij Is the minimum of (2);
obtaining the activity ratio A of each pixel ij
Figure BDA0002147368610000063
And (3) mapping the activity ratio of each pixel to obtain an activity ratio image through color mapping, and completing reconstruction.
The invention uses the temperature change characteristic of each part of human body when the part is diseased, namely the part can correspondingly change the heat energy distribution, in additionAfter the bits are cold treated, the temperature recovers faster than the adjacent parts where no lesions occur. The method comprises the steps of obtaining a continuous image sequence in a recovery process by adopting a thermal infrared imager, obtaining a preprocessed image after the image is stable after no change, dividing the preprocessed image into pixels, reconstructing the pixels, and defining a temperature recovery rate R of each pixel ij And refrigerating rate C ij The ratio of (2) is the activity ratio, and the corresponding activity ratio image is obtained through color mapping of the activity ratio and is used for diagnosing the part to be tested. The temperature of the part with lesions can be recovered faster, so that diagnosis can be completed with higher efficiency by adopting the method; meanwhile, the activity obtained through reconstruction is clearer than that of an image, the contrast is higher, and the outline of the image is clearer.
Example 1
(1) Exposing thyroid position at 24+ -2deg.C and humidity of 40-60% under constant temperature and humidity for 2min, and cold compressing with ice bag containing ice-water mixture for 30s;
(2) Fixing the thermal infrared imager above the thyroid gland part in the same constant temperature and humidity environment, starting to collect images after the thermal infrared imager is started for 5min, continuously collecting the images of the thyroid gland part after cold compress by using the thermal infrared imager, continuously collecting the images at the frequency of 0.5 s/sheet until the images are not changed any more and recorded as preprocessed images, stopping collecting the images, and collecting the images for 150s;
(3) Dividing the preprocessed image into 240 x 340 pixels, and marking the highest temperature value of each pixel in the process of collecting the image as T maxij Taking all pixels T maxij The minimum value of (2) is recorded as the reference temperature T R
(4) Calculating the temperature recovery rate R of each pixel ij And refrigeration rate C ij Obtaining the activity ratio of each pixel
Figure BDA0002147368610000071
And then the activity ratio image is obtained through color mapping to complete reconstruction.
Example two
(1) Exposing the forearm part at 24+ -2deg.C and 40-60% humidity for 2min, and cold compressing with ice bag containing ice-water mixture for 30s;
(2) Fixing the thermal infrared imager above the forearm under the same constant temperature and humidity environment, starting to collect images after the thermal infrared imager is started for 5min, continuously collecting the images of the forearm after cold compress by using the thermal infrared imager, continuously collecting the images at the frequency of 0.5 s/sheet until the images are no longer changed and recorded as preprocessed images, stopping collecting the images, and collecting the images for 140s. Obtaining a series of images as shown in fig. 2, it can be seen that no blood vessel is visible at 0.5s after removal of the cold treatment, the blood vessel region can be clearly seen at 20s, the blood vessel region is basically stable at 130s, and the blood vessel region is completely stable at 140 s;
(3) Dividing the preprocessed image into 240 x 340 pixels, and marking the highest temperature value of each pixel in the process of collecting the image as T maxij Taking all pixels T maxij The minimum value of (2) is recorded as the reference temperature T R
(4) Calculating the temperature recovery rate R of each pixel ij And refrigeration rate C ij Obtaining the activity ratio of each pixel
Figure BDA0002147368610000081
And obtaining an activity ratio image shown in figure 3 through color mapping, and completing reconstruction.
Fig. 4 and 5 show τ and τ' images obtained by the same cold treatment process and the same thermal infrared imager under the same environment. As can be seen from fig. 3, fig. 4 and fig. 5, the morphological detail features of subcutaneous blood vessels can be clearly represented by using the activity ratio reconstructed image, and the contrast ratio is better. Wherein the conventional tau image is a time-based reconstructed image, and the temperature recovery time tau after cold treatment of each pixel in the image is due to the shorter recovery time required for the subcutaneous region of vascular or tumor tissue ij The numerical value is smaller, in order to be convenient for comparison with the activity ratio image of the invention, the tau image is inverted to obtain a tau 'image, and the numerical value of each pixel in the tau' image is as follows:
Figure BDA0002147368610000082
wherein τ max Is the maximum value of the temperature recovery time of all pixels in the tau image. The values of the vascular regions in the τ' images are larger and are easier to compare with the activity ratio images.
In order to realize the infrared image reconstruction method, the invention also comprises an infrared image reconstruction system which comprises a thermal infrared imager, an image processing module and a display module;
the infrared thermal imager continuously collects the images of the part to be detected after the cold treatment until the images are no longer changed and recorded as preprocessed images, and the collection is stopped;
the image processing module is used for dividing the preprocessed image into i x j pixels, wherein i is the number of rows of the pixels, and j is the number of columns of the pixels; according to the refrigerating rate C of each pixel cold treatment ij And a temperature recovery rate R during continuous image acquisition ij Obtaining the activity ratio of each pixel
Figure BDA0002147368610000091
Performing color mapping on the activity ratio of each pixel to obtain an activity ratio image, and sending the activity ratio image to a display module;
and the display module displays the activity ratio image.
The reconstruction system can complete the acquisition and reconstruction of the infrared image, and finally an activity ratio image with higher contrast ratio is obtained.
In addition, the infrared image reconstruction method of the invention can be applied to a computer readable storage medium or a computer device, and the acquired preprocessed image is only required to be stored or input in advance.
If applied to a computer readable storage medium, the computer readable storage medium stores a computer program which when executed by a processor performs the above-described method steps; if applied to a computer device, the computer device includes a memory, a processor, and a computer program stored on the memory and executable on the processor, the processor implementing the steps described above when executing the program. The reconstruction method can be more conveniently executed by taking the readable storage medium and the computer as media, and the result is more accurate.
In addition, the method related by the invention is not applied to the purpose of disease diagnosis in the meaning of patent law, but can be used for teaching, research experiments, simulation, establishment of enterprise simulation databases and the like, and is beneficial to wider research by obtaining an activity ratio image with better contrast and contour definition.
The foregoing description is only illustrative of the present invention and is not intended to limit the scope of the invention, and all equivalent structural changes made by the present invention and the accompanying drawings, or direct or indirect application in other related technical fields, are included in the scope of the present invention.

Claims (9)

1. A medical infrared image reconstruction method for non-disease diagnosis purposes, comprising the steps of:
step 1, collecting infrared images
Continuously acquiring the image of the part to be detected after cold treatment by using a thermal infrared imager until the image is no longer changed, stopping acquisition, and recording as a preprocessed image;
step 2, dividing the pixels
Dividing the preprocessed image into i x j pixels, wherein i is the number of rows of the pixels, and j is the number of columns of the pixels;
step 3, reconstructing an image
Calculating the temperature recovery rate of each pixel:
Figure FDA0002147368600000011
the refrigeration rate of each pixel is calculated:
Figure FDA0002147368600000012
wherein t is 1 For cold treatment time, t 2 For the duration of the acquired image; t (T) 1ij Is the ith row and the jth columnTemperature before cold treatment of pixel, T 2ij The temperature after cold treatment of the ith row and the jth column of pixels; the highest temperature value of each pixel in the process of acquiring the infrared image in the step 1 is recorded as T maxij Reference temperature T R For each picture element T maxij Is the minimum of (2);
obtaining the activity ratio A of each pixel ij
Figure FDA0002147368600000013
And (3) mapping the activity ratio of each pixel to obtain an activity ratio image through color mapping, and completing reconstruction.
2. A medical infrared image reconstruction method as set forth in claim 1, wherein: in the step 1, the cold treatment is specifically that the part to be measured is exposed for at least 2min under the constant temperature and humidity environment with the temperature range of 24+/-2 ℃ and the humidity range of 40-60%, and then is cold-compressed by an ice bag, and the ice bag is filled with ice-water mixture.
3. A medical infrared image reconstruction method as set forth in claim 2, wherein: in the step 1, the thermal infrared imager is used for continuously collecting images of the part to be detected after the cold treatment, and the images are carried out under the constant temperature and humidity environment with the temperature range of 24+/-2 ℃ and the humidity range of 40-60%.
4. A medical infrared image reconstruction method as set forth in claim 1, wherein: in the step 1, continuously acquiring images of the part to be detected subjected to cold treatment by using a thermal infrared imager, wherein the continuous acquisition frequency is 0.5 s/sheet.
5. A medical infrared image reconstruction method as set forth in claim 1, wherein: in step 1 and step 2, the t 1 And t 2 Is accurate in seconds.
6. A medical infrared image reconstruction method as set forth in claim 1, wherein: in step 2, the preprocessed image is divided into i×j pixels, where i is 240 and j is 320.
7. A reconstruction system for implementing the medical infrared image reconstruction method according to any one of claims 1 to 6, characterized in that: the system comprises a thermal infrared imager, an image processing module and a display module;
the infrared thermal imager continuously collects the images of the part to be detected after the cold treatment until the images are no longer changed and recorded as preprocessed images, and the collection is stopped;
the image processing module is used for dividing the preprocessed image into i x j pixels, wherein i is the number of rows of the pixels, and j is the number of columns of the pixels; according to the refrigerating rate C of each pixel cold treatment ij And a temperature recovery rate R during continuous image acquisition ij Obtaining the activity ratio of each pixel
Figure FDA0002147368600000021
Performing color mapping on the activity ratio of each pixel to obtain an activity ratio image, and sending the activity ratio image to a display module;
and the display module displays the activity ratio image.
8. A computer readable storage medium having stored thereon a computer program, characterized in that the program when executed by a processor performs the steps of:
step 1, dividing pixels
Dividing the preprocessed image into i x j pixels, wherein i is the number of rows of the pixels, and j is the number of columns of the pixels; the preprocessing image is obtained by continuously collecting the image of the part to be detected after cold processing by using a thermal infrared imager until the image is not changed any more, stopping collecting, and recording the image as the preprocessing image;
step 2, reconstructing an image
Calculating the temperature recovery rate of each pixel:
Figure FDA0002147368600000031
the refrigeration rate of each pixel is calculated:
Figure FDA0002147368600000032
wherein t is 1 For cold treatment time, t 2 For the duration of the acquired image; t (T) 1ij T is the temperature before cold treatment of the ith row and jth column pixels 2ij The temperature after cold treatment of the ith row and the jth column of pixels; the highest temperature value of each pixel in the process of acquiring the infrared image in the step 1 is recorded as T maxij Reference temperature T R For each picture element T maxij Is the minimum of (2);
obtaining the activity ratio A of each pixel ij
Figure FDA0002147368600000033
And (3) mapping the activity ratio of each pixel to obtain an activity ratio image through color mapping, and completing reconstruction.
9. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the following steps when executing the program:
step 1, dividing pixels
Dividing the preprocessed image into i x j pixels, wherein i is the number of rows of the pixels, and j is the number of columns of the pixels; the preprocessing image is obtained by continuously collecting the image of the part to be detected after cold processing by using a thermal infrared imager until the image is not changed any more, stopping collecting, and recording the image as the preprocessing image;
step 2, reconstructing an image
Calculating the temperature recovery rate of each pixel:
Figure FDA0002147368600000034
the refrigeration rate of each pixel is calculated:
Figure FDA0002147368600000041
wherein t is 1 For cold treatment time, t 2 For the duration of the acquired image; t (T) 1ij T is the temperature before cold treatment of the ith row and jth column pixels 2ij The temperature after cold treatment of the ith row and the jth column of pixels; the highest temperature value of each pixel in the process of acquiring the infrared image in the step 1 is recorded as T maxij Reference temperature T R For each picture element T maxij Is the minimum of (2);
obtaining the activity ratio A of each pixel ij
Figure FDA0002147368600000042
And (3) mapping the activity ratio of each pixel to obtain an activity ratio image through color mapping, and completing reconstruction.
CN201910689274.9A 2019-07-29 2019-07-29 Medical infrared image reconstruction method, system, computer and storage medium Active CN110298902B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910689274.9A CN110298902B (en) 2019-07-29 2019-07-29 Medical infrared image reconstruction method, system, computer and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910689274.9A CN110298902B (en) 2019-07-29 2019-07-29 Medical infrared image reconstruction method, system, computer and storage medium

Publications (2)

Publication Number Publication Date
CN110298902A CN110298902A (en) 2019-10-01
CN110298902B true CN110298902B (en) 2023-04-21

Family

ID=68032165

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910689274.9A Active CN110298902B (en) 2019-07-29 2019-07-29 Medical infrared image reconstruction method, system, computer and storage medium

Country Status (1)

Country Link
CN (1) CN110298902B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3005540B1 (en) * 1998-06-25 2000-01-31 勲 内田 Calculation image creation method
CA2333406A1 (en) * 1998-05-27 2000-10-22 Omnicorder Technologies, Inc. Method of detecting cancerous lesions
CN1823683A (en) * 2006-04-03 2006-08-30 何宗彦 Method of dynamic detecting organism parameter and its medical detecting analyzer

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013172963A1 (en) * 2012-05-18 2013-11-21 Sanjay Krishna A method and system for feature extraction and decision making from series of images
WO2015061588A1 (en) * 2013-10-23 2015-04-30 The Trustees Of Dartmouth College Surgical vision augmentation system
US9865052B2 (en) * 2016-03-18 2018-01-09 Niramai Health Analytix Pvt Ltd Contour-based determination of malignant tissue in a thermal image

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2333406A1 (en) * 1998-05-27 2000-10-22 Omnicorder Technologies, Inc. Method of detecting cancerous lesions
JP3005540B1 (en) * 1998-06-25 2000-01-31 勲 内田 Calculation image creation method
CN1823683A (en) * 2006-04-03 2006-08-30 何宗彦 Method of dynamic detecting organism parameter and its medical detecting analyzer

Also Published As

Publication number Publication date
CN110298902A (en) 2019-10-01

Similar Documents

Publication Publication Date Title
Li et al. The obf database: A large face video database for remote physiological signal measurement and atrial fibrillation detection
Adam et al. Automated characterization of diabetic foot using nonlinear features extracted from thermograms
US8360986B2 (en) Non-contact and passive measurement of arterial pulse through thermal IR imaging, and analysis of thermal IR imagery
US20040076316A1 (en) Method and apparatus for measuring physiology by means of infrared detector
CN105342565A (en) Traditional Chinese medicine channels and collaterals analysis method based on human body surface infrared medical science image
WO2006060727A9 (en) System and method for eeg imaging of cerebral activity using electrode sets
Bergquist et al. Body surface potential mapping: contemporary applications and future perspectives
CN112294282A (en) Self-calibration method of emotion detection device based on RPPG
Li et al. An 8-layer residual U-Net with deep supervision for segmentation of the left ventricle in cardiac CT angiography
CA2524807C (en) System and method for identifying and classifying dynamic thermodynamic processes in mammals and discriminating between and among such processes
CN110298902B (en) Medical infrared image reconstruction method, system, computer and storage medium
US8582852B2 (en) Dual-spectrum heat pattern separation algorithm for assessing chemotherapy treatment response and early detection
US7408156B2 (en) System and method for identifying and classifying dynamic thermodynamic processes in mammals and discriminating between and among such processes
US8233968B1 (en) Method and apparatus for high resolution dynamic digital infrared imaging
CN115474930A (en) Hyperspectral image reconstruction-based noninvasive hemoglobin detection method
JP3005540B1 (en) Calculation image creation method
CN114916916A (en) Integrated medical equipment, scanning control method and device thereof, and computer equipment
Arabi et al. Early Diagnosis of Diabetic Vasculopathy by Thermo-Regulatory Vascular Impairment
Cipresso et al. Big data in preclinical ECG alterations research
Guryleva et al. Application Software for Visualization and Quantification of Capillary Networks
Pakarinen et al. IRlab-Platform for thermal video analysis in evaluation of peripheral thermal behavior and blood perfusion
CN114209299B (en) IPPG technology-based human physiological parameter detection channel selection method
CN113143264B (en) Blood glucose detection area selection method based on blood perfusion imaging
Czajkowska et al. An automatic method of chronic wounds segmentation in multimodal images
CN117357105A (en) Optical noninvasive blood glucose concentration detection method based on time-spectrum-space multiparameter fusion

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant