Method of High Precision Medical Infrared Thermography
-
摘要: 医学红外热像设备测得的红外数据及转换得到的温度数据难以直接判定其所属的人体区域，常需将其转为图像数据，利用图像处理技术得到感兴趣区域并从区域内温度数据得到生物特征，实现疾病的筛查或诊断。然而，从14位红外数据转换到8位图像数据存在严重的数据精度损失，导致处理性能欠佳。本文提出一种新的热像图表达方法，所得到的彩色热像图含原精度的温度数据信息，且含温度观察窗设定尺度下的彩色增强效果，同时载有温度数据记录和观察窗设定规则，通过对图像数据的逆变换，可以再现原始温度数据，并可改变彩色增强效果。该方法提供的热像图无需额外存取温度数据文件，在不同的红外热像系统间具有通用性，将更符合大数据和人工智能的发展趋势。Abstract: It is difficult to directly distinguish the human body area identified by the infrared data measured by medical infrared thermal imaging equipment from the temperature data obtained by conversion. It is often necessary to convert it into image data and use image processing technology to obtain the region of interest and the biological characteristics from the temperature data in the given area. Accordingly, disease screening or diagnosis can be realized. However, conversion from 14-bit infrared data to 8-bit image data incurs a serious loss of data accuracy, resulting in poor processing performance. In this paper, a new expression method for thermal images is proposed. The obtained color thermal image contains the original precision temperature data information and the color enhancement effect under the setting scale of the temperature observation window. At the same time, it contains the temperature data record and the setting rules of the observation window. Through the inverse transformation of the image data, the original temperature data can be reproduced and the color enhancement effect can be changed. The thermal image provided by this method can be used in different infrared thermal image systems without requiring additional access to temperature data files. This will be more aligned with the development trend of big data and artificial intelligence.
-
0. 引言
混凝土材料往往因为预埋钢筋的腐蚀、水的侵蚀、连续的冻融循环、荷载等多种因素共同作用，产生裂缝、孔洞和分层等缺陷[1]。这些内部缺陷随着时间的推移会逐渐积聚、扩展，严重威胁着结构的耐久性和安全性[2]。主动红外热像检测方法作为一种无损检测手段，由于具有检测面积大、非接触式和应用操作简单等优点在混凝土结构检测领域得到了广泛的关注，并逐渐被应用到了隧道、桥梁现场检测中。
主动红外热像检测技术是一种利用红外热像仪检测材料内部缺陷的无损检测方法，其检测的基本原理是：材料被热激励后其表面红外辐射量因内部缺陷的存在而表现出差异，红外热像仪通过记录该差异，以获得材料内部的缺陷信息[3]。红外图像的信息提取与分析是红外热像检测技术的核心内容。理论上，根据热图像缺陷区域和正常区域之间的红外辐射差异可获得缺陷的信息。但实际上，红外图像本质是根据缺陷与正常区域之间的温度和辐射率生成的灰度图，而且红外光辐射的能量远比可见光低，因此到达光学镜头的红外辐射能很小，导致了红外图像与可见光图像相比只有较低的对比度，分辨弱小目标和细节能力差[4]；另外，红外光子的随机性、红外成像系统的固有特性以及热激励不均匀等给红外图像引入了许多噪声，在一定程度上改变了图像的原始信息，增加了图像分割及后续图像处理的难度[5-6]。传统的基于空域的红外图像处理方法主要分为两个大的步骤：首先，是对红外图像进行降噪和增强，常用的方法有高斯滤波、中值滤波、引导滤波、小波变换、Contourlet变换、三维块匹配算法等[7-10]；其次，对红外图像进行边缘特征的提取，常用的方法有：基于边缘的检测算法、基于数学形态的检测算法和最近发展较快的基于网络模型的检测算等[11-12]。传统的基于单张空域的红外图像处理方法虽然在一定程度上可以消除噪声的影响、提高图像的对比度，但是仍存在一些问题。第一，采集的数据集包含大量的热图像，每一帧红外图像都对应着某一时刻，依靠视觉手动选择缺陷信息丰富的红外图像进行缺陷信息的提取的识别费时费力。第二，在图像降噪、增强和图像分割过程中设定阈值会引入主观成分，干扰红外图像的解释。第三，仅仅分析单张红外图像，忽略采集过程中的时序信息可能会导致忽略一些边缘的缺陷信息[13-14]。
针对上述问题，本文提出了一种基于时序信息的红外图像缺陷信息提取方法。首先，通过室内实验制作含缺陷分层的混凝土试块，其次利用主动红外热像检测技术进行红外数据的采集，提取每个像素点的温度特征曲线，然后采用基于时序信息的提取方法进行含分层缺陷混凝土的缺陷提取，并与传统的基于空域的处理方法进行对比分析。
1. 基本原理
1.1 主动红外热像检测技术原理和系统
主动红外热像检测技术中，当以一定的距离对混凝土材料表面持续加热时，物体会吸收入射波的能量并将其转化为热能，并以温度场的形式表现出来。如果试样存在缺陷（本研究采用的是隔热型缺陷，模拟混凝土内部的空气分层缺陷），到达缺陷的大部分能量将反射到试样表面，引起表面正常区域和有缺陷区域温度场的差异。
利用红外热像仪记录整个检测过程中物体表面温度场的空间和时间分布信息，通过对红外图像数据进行分析和处理进而可以提取材料近表层的分层缺陷信息。主动红外热像检测技术采集到的一系列红外图像组成了一个三维的图像数据集[15]（尺寸：Nx×Ny×Nt，见图 1(a)），在空间维度上，每一帧对应着某一时刻红外热像仪所采集到的物体表面温度场信息，其像素个数为Nx×Ny；在时间维度上，每个像素点可以看作是一条随时间变化的温度特征曲线（Nt时间点组成），如图 1(b)所示，其中标签NT11 DEFECT表示缺陷区域的温度特征曲线，标签NT11 NORMAL表示正常区域的温度特征曲线。
![]()
图 1 主动红外热成像采集的数据结构：(a) 三维红外数据；(b) 像素点温度特征曲线Figure 1. Active infrared thermal imaging acquisition of data structures (a) Three-dimensional infrared data; (b) Temperature characteristic curves of pixels1.2 基于时序信息的缺陷信息提取原理
主动红外热像技术采集到的三维红外图像数据，其正常区域和含缺陷的区域每个像素点的温度特征曲线在冷却阶段会存在差异，如图 1(b)所示，利用这个特点我们可以采用许多能够分辨这个特征曲线差异值的相关算法进行缺陷信息的提取，本文采用相对简单、容易实施的K-means方法来提取混凝土的缺陷信息。K-means聚类方法的原理是同一类内的实体是相似的，一个类是测试空间中点的集合，同一类内任意一个点到其类中心的距离小于其到其它类中心的距离，我们可以用此算法分辨正常区域和缺陷区域的温度特征曲线[16]。K-means方法的关键是核函数的选取，核函数决定了其分辨能力，本文采用常用的马氏距离进行度量。其具体算法原理如下[17]：
假设给定的数据集为Y={y1, y2, …, ym}，其中yi=[X1, X2, …, Xn]，将数据对象划分为k类C={c1, c2, …, ck}，每个类有一个类中心U={u1, u2, …, uk}。选取欧式距离作为相似性和距离判断准则，计算ck类内各点到聚类中心uk的距离平方和，见公式(1)：
$$ J\left( {{c_k}} \right) = \sum\limits_{{x_i} \in {C_k}} {{{\left\| {{x_i} - {u_k}} \right\|}^2}} $$ (1) 聚类的目标是使各类总的距离平方和最小，见公式(2)：
$$ J(C) = \sum\limits_{k = 1}^K J \left( {{c_k}} \right) = \sum\limits_{k = 1}^K {\sum\limits_{{x_i} \in {C_k}} {{{\left\| {{x_i} - {u_k}} \right\|}^2}} } $$ (2) 根据最小二乘法和拉格朗日原理，聚类中心uk应该取为类别ck类各数据点的平均值。
K-means算法是一个反复迭代过程，目的是使聚类域中所有的点到聚类中心距离的平方和最小。
2. 实验与数据分析
2.1 实验过程
本次实验混凝土试件的设计强度为C50，尺寸为50 cm×50 cm×20 cm，在混凝土中嵌入4块尺寸为10 cm×10 cm的聚苯乙烯材料（隔热型材料），模拟混凝土中的分层缺陷，混凝土试块见图 2(a)。试验采用美国FLIR公司的A655SC非制冷型红外热像仪，其热灵敏度为30 mK，图像分辨率为640×480，标准测温范围为－40℃~150℃，波长范围为7.5~14 μm，见图 2(b)。主动热激励系统采用自制的配有自动控制加热时间的碳化硅远红外加热板，共2块，每块加热板的加热功率为500 W，见图 2(c)。整个实验过程中加热时间为5 min，冷却时间为10 min；加热完成后，立即用红外热像仪采集试样表面的冷却温度场，整个采集系统见图 3所示，本文选用的是持续加热法。
![]()
图 2 实验试块和仪器：(a) 混凝土试块；(b) 红外热像仪；(c) 红外加热板Figure 2. Test blocks and instruments: (a) Concrete test block; (b) Infrared thermal imager; (c) Infrared heating plate2.2 红外图像数据
实验采集到的红外图像序列如图 4所示，根据红外图像可以发现一共存在4块分层缺陷，各个缺陷之间存在温度场相互影响的干扰区，而且各个缺陷的轮廓比较模糊，存在被隐藏的缺陷信息。整个采集的红外图像数据的第一帧（0 s）最清晰，其正常区域和缺陷区域存在较大的温差；之后随着冷却的进行，正常区域和缺陷区域的温差逐渐减小，缺陷信息逐渐模糊，至最后一帧很难用肉眼分辨。
3. 红外图像处理和对比分析
3.1 红外图像缺陷提取
主动红外热像检测技术所采集的图像每一帧都含有大量的噪声，为了测试基于时序信息红外图像缺陷信息提取的可行性以及性能，本次在图像特征提取前不对图像进行预处理。首先，将红外图像数据按像素点逐个提取时间序列，组成聚类数据集；然后将时序数据序列输入到K-means程序中，进而得到每一帧各个像素点的聚类信息。为了提高运算效率，本次在图像序列处理中采样间隔为1，其序列长度减小为原来的一半。聚类完成后得到图像中每一个像素点所属的类别，进而完成对图像缺陷信息的提取。
图 5为经过聚类后的提取结果，通过观察发现4个缺陷信息都被提取出来，与原始红外图像相比（图 4），其缺陷信息有了明显的增强，缺陷信息完全肉眼可见，每个识别出的缺陷区域呈近似正方形，完整性相当好，接近于预埋缺陷的形状，而且缺陷提取后的图像序列的时序信息得以保留。
3.2 对比分析
为了测试基于时序信息红外图像分层缺陷的信息提取效果，本文采用基于空域的二维K-means图像缺陷信息提取算法与其进行对比。由于采集的红外图像数据集第一帧（0 s）缺陷信息最清晰，因此选用第一帧红外图像进行基于空域的分层缺陷信息提取。两种方法的分层缺陷信息提取结果见图 6，由图 6(b)所示，基于时序信息的分层缺陷提取方法4个分层缺陷都被完全提取出来，图 6(a)中模糊、隐藏的信息也被提取出来，缺陷细节分辨能力较好。由图 6(c)所示，基于空域信息的分层缺陷信息提取方法4个分层缺陷被识别为一个缺陷，缺陷细节分辨能力较差。通过分析可能是受缺陷之间温度场叠加的相互影响，缺陷之间区域（图 6(a)中虚线框内部所示的无缺陷区域）的温度信息和分层缺陷的温度场信息较为接近，而和图 6(a)的无缺陷边缘温度场信息差异较大，导致基于空域信息的提取方法效果较差；根据含缺陷混凝土试块的物理特性可知，存在分层缺陷和正常区域的温度特征曲线是存在差异的，基于时序信息的提取方法主要是根据每个像素点的温度特征曲线差异来提取特征，其提取缺陷的分辨能力和效果更好。
![]()
图 6 处理结果对比分析：(a) 红外原图；(b) 基于时序K-means方法；(c) 基于空域K-means方法Figure 6. Comparative analysis of processing results: (a) Original infrared image; (b) K-means method based on temporal information; (c) K-means method based on spatial information4. 结论
主动红外热像检测技术中，传统的红外图像处理方法在一定程度上可以消除噪声、提高图像的对比度，但是仍存在一些问题。针对上述问题，本文根据主动红外热成像的数据特点提出了一种基于时序信息的红外图像缺陷信息提取方法。结果表明，基于时序信息的缺陷提取方法是可行的，其可以提取到隐藏的缺陷信息，分层缺陷信息提取效果优于基于空域的K-means方法。
-
![]()
图 1 红外热像机输出的原灰度图和伪彩色图
Figure 1. Original gray scale and pseudo color image output by infrared thermal imager
![]()
图 2 新方法得到的红外热图像及其色带
Figure 2. Infrared thermal image and its color band obtained by new method
![]()
图 3 对新的红外热像图通过温度观察窗进行增强
Figure 3. Enhance the new infrared thermal image through the temperature observation window
![]()
图 5 原灰度图像和新灰度图像的边缘检测对比
Figure 5. Comparison of edge detection between original gray image and new gray image
![]()
图 6 新图像去除背景信息的效果图
Figure 6. Effect picture of removing background information from new image
![]()
图 7 不同RGB排列方式得到的新的红外热像图的效果图
Figure 7. Effect picture of new infrared thermal image obtained by different RGB arrangement
-
[1] Bhowmik M K, Gogoi U R, Majumdar G, et al. Designing of Ground-Truth-Annotated DBT-TU-JU Breast Thermogram Database Toward Early Abnormality Prediction[J]. IEEE Journal of Biomedical and Health Informatics, 2018, 22: 1238-1249. DOI: 10.1109/JBHI.2017.2740500
[2] De Santana M A, Pereira J M S, Da Silva F L, et al. Breast cancer diagnosis based on mammary thermography and extreme learning machines[J]. Research on Biomedical Engineering, 2018, 34: 45-53. DOI: 10.1590/2446-4740.05217
[3] Dua G, Mulaveesala R. Applicability of active infrared thermography for screening of human breast: a numerical study[J]. Journal of Biomedical Optics, 2018, 23: 9. https://www.ncbi.nlm.nih.gov/pubmed/29560626
[4] Mambou S J, Maresova P, Krejcar O, et al. Breast Cancer Detection Using Infrared Thermal Imaging and a Deep Learning Model[J]. Sensors, 2018, 18: 19.
[5] Morales-Cervantes A, Kolosovas-Machuca E S, Guevara E, et al. an Automated Method for the Evaluation of Breast Cancer Using Infrared Thermography[J]. Excli Journal, 2018, 17: 989-998. https://www.ncbi.nlm.nih.gov/pubmed/30564079
[6] Santana M A d, Pereira J M S, Silva F L d, et al. Breast cancer diagnosis based on mammary thermography and extreme learning machines[J]. Research on Biomedical Engineering, 2018, 34: 45-53. DOI: 10.1590/2446-4740.05217
[7] Wahab A A, Salim M I M, Yunus J, et al. Comparative evaluation of medical thermal image enhancement techniques for breast cancer detection[J]. Journal of Engineering and Technological Sciences, 2018, 50: 40-52. DOI: 10.5614/j.eng.technol.sci.2018.50.1.3
[8] Abdel-Nasser M, Moreno A, Puig D. Breast Cancer Detection in Thermal Infrared Images Using Representation Learning and Texture Analysis Methods[J]. Electronics, 2019, 8: 18. http://www.researchgate.net/publication/330412196_Breast_Cancer_Detection_in_Thermal_Infrared_Images_Using_Representation_Learning_and_Texture_Analysis_Methods
[9] Singh D, Singh A K. Role of image thermography in early breast cancer detection- Past, present and future[J]. Computer Methods and Programs in Biomedicine, 2020, 183: 61-69. https://www.sciencedirect.com/science/article/pii/S0169260719311277
[10] Fokam D, Lehmann C. Clinical assessment of arthritic knee pain by infrared thermography[J]. Journal of basic and clinical physiology and pharmacology, 2018, 30: 21-25. http://www.ncbi.nlm.nih.gov/pubmed/30375348
[11] Pauk J, Wasilewska A, Ihnatouski M. Infrared thermography sensor for disease activity detection in Rheumatoid arthritis patients[J]. Sensors (Switzerland), 2019, 19: 34-48. https://www.mdpi.com/1424-8220/19/16/3444/pdf
[12] Pauk J, Ihnatouski M, Wasilewska A. Detection of inflammation from finger temperature profile in rheumatoid arthritis[J]. Medical & Biological Engineering & Computing, 2019, 57: 2629-2639. DOI: 10.1007/s11517-019-02055-1.pdf
[13] Gatt A, Mercieca C, Borg A, et al. A comparison of thermographic characteristics of the hands and wrists of rheumatoid arthritis patients and healthy controls[J]. Scientific Reports, 2019, 9: 172-180. DOI: 10.1038/s41598-018-36890-3
[14] Haq T, Crane J D, Kanji S, et al. Optimizing the methodology for measuring supraclavicular skin temperature using infrared thermography; implications for measuring brown adipose tissue activity in humans[J]. Scientific Reports, 2017, 7: 9. DOI: 10.1038/s41598-017-00053-7
[15] Jimenez-Pavon D, Corral-Perez J, Sanchez-Infantes D, et al. Infrared Thermography for Estimating Supraclavicular Skin Temperature and BAT Activity in Humans: A Systematic Review[J]. Obesity, 2019, 27: 1932-1949. DOI: 10.1002/oby.22635
[16] LIN P H, Echeverria A, Poi M J. Infrared thermography in the diagnosis and management of vasculitis[J]. Journal of vascular surgery cases and innovative techniques, 2017, 3: 112-114. DOI: 10.1016/j.jvscit.2016.12.002
[17] Gauci J, Falzon O, Formosa C, et al. Automated Region Extraction from Thermal Images for Peripheral Vascular Disease Monitoring[J]. Journal of Healthcare Engineering, 2018, 2018: 14. https://www.hindawi.com/journals/jhe/2018/5092064/
[18] Carriere M E, de Haas L E M, Pijpe A, et al. Validity of thermography for measuring burn wound healing potential[J]. Wound Repair and Regeneration, 2019, 10: 1-8. https://pubmed.ncbi.nlm.nih.gov/31777128/
[19] Knobel-Dail R B, Holditch-Davis D, Sloane R, et al. Body temperature in premature infants during the first week of life: Exploration using infrared thermal imaging[J]. Journal of Thermal Biology, 2017, 69: 118-123. DOI: 10.1016/j.jtherbio.2017.06.005
[20] Topalidou A, Ali N, Sekulic S, et al. Thermal imaging applications in neonatal care: a scoping review[J]. Bmc Pregnancy and Childbirth, 2019, 19: 14. DOI: 10.1186/s12884-018-2132-3
[21] Pereira T, Nogueira-Silva C, Simoes R. Normal range and lateral symmetry in the skin temperature profile of pregnant women[J]. Infrared Physics & Technology, 2016, 78: 84-91. https://www.sciencedirect.com/science/article/pii/S1350449516302067
[22] Martini G, Cappella M, Culpo R, et al. Infrared thermography in children: a reliable tool for differential diagnosis of peripheral microvascular dysfunction and Raynaud's phenomenon?[J]. Pediatric Rheumatology, 2019, 17: 9. DOI: 10.1186/s12969-019-0307-8
[23] Garcia-Porta N, Gantes-Nunez F J, Tabernero J, et al. Characterization of the ocular surface temperature dynamics in glaucoma subjects using long-wave infrared thermal imaging[J]. Journal of the Optical Society of America a-Optics Image Science and Vision, 2019, 36: 1015-1021. DOI: 10.1364/JOSAA.36.001015
[24] Debiec-Bak A, Wojtowicz D, Pawik L, et al. Analysis of body surface temperatures in people with Down syndrome after general rehabilitation exercise[J]. Journal of Thermal Analysis and Calorimetry, 2019, 135: 2399-2410. DOI: 10.1007/s10973-018-7345-1
[25] Hernandez-Contreras D A, Peregrina-Barreto H, Rangel-Magdaleno J D, et al. Plantar Thermogram Database for the Study of Diabetic Foot Complications[J]. IEEE Access, 2019, 7: 161296-161307. DOI: 10.1109/ACCESS.2019.2951356
[26] 丁德红. 16位高精度在线式红外热像仪的技术方案与实现[J].红外技术, 2017, 39(9): 841-847. http://www.cnki.com.cn/Article/CJFDTOTAL-HWJS201709011.htm Dehong D. Online Temperature Measurement Technology Solutions and Implementationof 16 bit Infrared Thermal Imager[J]. Infrared Technology, 2017, 39(9): 841-847. http://www.cnki.com.cn/Article/CJFDTOTAL-HWJS201709011.htm
[27] Tan J H, Acharya U R. Pseudocolours for thermography-Multi-segments colour scale[J]. Infrared Physics & Technology, 2015, 72: 140-147. https://www.sciencedirect.com/science/article/pii/S1350449515001875
[28] Kermani S, Samadzadehaghdam N, EtehadTavakol M. Automatic color segmentation of breast infrared images using a Gaussian mixture model[J]. Optik, 2015, 126: 3288-3294. DOI: 10.1016/j.ijleo.2015.08.007
[29] LI T J, WANG Y Y, CHANG C, et al. Color-appearance-model based fusion of gray and pseudo-color images for medical applications[J]. Information Fusion, 2014, 19: 103-114. DOI: 10.1016/j.inffus.2012.07.002
-
期刊类型引用(4)
1. 迟爽，吕琴，刘晓亭，于游. 太赫兹技术在中医学领域的应用与展望. 时珍国医国药. 2024(02): 425-428 . 
百度学术
2. 孙树祥，梅红樱，刘真，郑新艳，孙东雪，司仪寒. 温度对艾草不同部位药性影响的太赫兹光谱分析. 信阳农林学院学报. 2023(02): 101-105 . 百度学术
3. 王杰. 基于多光谱特征分区的油桃品质分析算法. 食品与机械. 2022(05): 133-137 . 百度学术
4. 李晓，阎萍. 沉水植物与水相互作用的太赫兹光谱研究. 科学技术创新. 2021(23): 46-48 . 百度学术
其他类型引用(2)
下载:
计量
- 文章访问数: 817
- HTML全文浏览量: 433
- PDF下载量: 111
- 被引次数: 6
