CN112082655B - 一种基于横向剪切干涉信号测量体温的方法 - Google Patents
一种基于横向剪切干涉信号测量体温的方法 Download PDFInfo
- Publication number
- CN112082655B CN112082655B CN202010805349.8A CN202010805349A CN112082655B CN 112082655 B CN112082655 B CN 112082655B CN 202010805349 A CN202010805349 A CN 202010805349A CN 112082655 B CN112082655 B CN 112082655B
- Authority
- CN
- China
- Prior art keywords
- temperature
- arm
- phase difference
- measuring
- body 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
Links
- 230000036760 body temperature Effects 0.000 title claims abstract description 27
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000010008 shearing Methods 0.000 title claims abstract description 12
- 238000009529 body temperature measurement Methods 0.000 claims abstract description 13
- 239000013078 crystal Substances 0.000 description 5
- 230000005855 radiation Effects 0.000 description 5
- 238000005259 measurement Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000002277 temperature effect Effects 0.000 description 2
- 208000025721 COVID-19 Diseases 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/0022—Radiation pyrometry, e.g. infrared or optical thermometry for sensing the radiation of moving bodies
- G01J5/0025—Living bodies
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/02—Constructional details
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/02—Constructional details
- G01J5/08—Optical arrangements
- G01J5/0803—Arrangements for time-dependent attenuation of radiation signals
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Radiation Pyrometers (AREA)
- Measuring Temperature Or Quantity Of Heat (AREA)
Abstract
一种基于横向剪切干涉信号测量体温的方法,依据横向剪切干涉信号非接触测量体温。该方法首先采集无测量对象时的干涉信号,识别背景条纹,采集存在手臂温度影响时的干涉信号,识别测量条纹;依据背景条纹和测量条纹的相位分布,确定手臂温度影响区域的宽度以及造成的相位差分布;采用线性拟合方法确定靠近手臂皮肤位置处的相位差;依据手臂横截面的宽度和高度,结合手臂温度影响区域的宽度和实时测量的环境温度,设定可能的体温范围和步长,遍历设定的体温范围,计算每个温度值对应的相位差;选择计算相位差与已确定的靠近手臂皮肤处相位差最接近的温度值,并依据相邻温度值计算得到的相位差插值确定体温测量结果。
Description
技术领域
本发明涉及体温测量方法,特别是涉及一种基于横向剪切干涉信号测量体温的方法,属于干涉测温技术领域。
背景技术
体温测量是发现COVID-19疾病的一种有效方式。为了实现防疫的目标,测温设备需要不接触人体皮肤,因此只能采用非接触式测温方法。常用的红外测温枪采集人体发射的热辐射信号,处理产生的电信号,推算被测的人体温度,依据人体自身发射的热辐射信号进行体温测量。但红外测温枪的测量精度有限,具体原因包括:1)操作时测温距离不稳定,环境衰减作用不同,影响热辐射信号的接收;2)人为操作,无法完全垂直于待测皮肤表面,与理想运行工况存在偏差;3)人体体表状况不同,辐射能力存在差异。上述问题涵盖了测量环境、操作过程以及测温原理等多个方面,是红外测温枪无法避免的问题。所以发展一种可靠的非接触式体温测量方法具有重大的实际意义。
发明内容
为了避免基于人体热辐射信号进行体温测量所遇到的问题,本发明提供了一种基于横向剪切干涉信号测量体温的方法,该方法能够基于横向剪切干涉信号可靠地非接触测量体表温度。
本发明所采用的技术方案如下:
一种基于横向剪切干涉信号测量体温的方法,该方法包括如下步骤:
1)采集无测量对象时的干涉信号,识别背景条纹,确定相位分布;
2)采集存在手臂温度影响时的干涉信号,识别测量条纹,确定相位分布;
3)测量并记录环境温度;
4)依据背景条纹和测量条纹的相位分布,确定手臂温度影响区域的相位差分布,确定手臂温度影响区域的宽度;
5)依据相位差分布,采用线性拟合方法,确定靠近手臂皮肤处的相位差;
6)设定可能的体温范围为30℃至42℃,以0.1℃为步长,遍历设定的体温范围,共121个温度值。依据手臂横截面的宽度和高度,结合手臂温度影响区域的宽度和环境温度,通过公式(Ⅰ)计算每个温度值对应的相位差:
式中,T为温度,T0为环境温度,λ为波长,p为大气压强,M为气体分子量,K为气体折射系数,R为理想气体常数,L为手臂温度影响区域的宽度,w为手臂横截面宽度,h为手臂横截面高度,xb和yb分别为手臂面上距离离散点(xi,h/2)最近点的横坐标和纵坐标,Δx为离散微元的长度;
7)结合已确定的靠近手臂皮肤处相位差,选择使公式(Ⅱ)最小的温度值:
f(T)=|θmea-θcal,T| (Ⅱ)
式中,θmea为测量并拟合得到的靠近手臂皮肤处相位差,θcal,T为体温为T时通过公式(Ⅰ)计算得到的相位差。选择计算相位差与靠近手臂皮肤处相位差最接近的温度值,并依据相邻温度值计算得到的相位差插值确定体温测量结果。
本发明提供的方法能够基于横向剪切干涉信号测量体温,不接触待测皮肤,结果可靠。
附图说明
图1为本发明方法的总体流程图。
图2为手臂温度影响区域及光束传播过程示意图。
图3为实测背景条纹和存在手臂温度影响时的干涉条纹。
图4为中间高度位置的干涉信号。
图5为整理得到的干涉条纹。
图6为波峰位置的相位差及线性拟合结果。
图7为30℃至42℃不同温度值计算得到的相位差及靠近手臂皮肤处相位差。
具体实施方式
图1为本发明方法的总体流程图,该方法首先采集无测量对象时的干涉信号,识别并用作背景条纹,采集存在手臂温度影响时的干涉信号,识别测量条纹;依据背景条纹和测量条纹的相位分布,确定手臂温度影响区域的宽度以及造成的相位差分布;采用线性拟合方法确定靠近手臂皮肤位置处的相位差;依据手臂横截面的宽度和高度,结合手臂温度影响区域的宽度和实时测量的环境温度,设定可能的体温范围和步长,遍历设定的体温范围,计算每个温度值对应的相位差;结合已确定的靠近手臂皮肤处相位差,选择计算相位差与其最接近的温度值,并依据相邻温度值计算得到的相位差插值确定体温测量结果。具体实施过程如下:
一、横向剪切干涉条纹的产生
单色线光源传播经过待测手臂区域,部分被手臂挡住,剩余部分穿过手臂周围空气环境,一部分光由平行平晶前表面反射,一部分光穿透平行平晶前表面被后表面反射并穿透前表面,两部分光因相位差而产生干涉条纹。依据理想气体状态方程,密度ρ可以表示为:
式中,p为气体压强,M为气体分子量,R为理想气体常数,T为温度。依据格拉斯顿-戴尔方程,折射率n可以表示为密度的函数:
n=1+ρK (Ⅳ)
其中K为气体折射系数,也称作格拉斯顿-戴尔系数。联立可得:
波长为λ的光束穿过折射率不均匀的空气区域,波前会发生扭曲,与环境的温差会造成波面畸变。光束的相函数可以表示为:
因此,环境温度为T0的条件下,相对相函数可以表示为:
在光程为l的条件下,测量相位差θT可以表示为:
平行平晶前后表面反射的波前在横向方向会有位移s,两束波前相互叠加,产生干涉条纹。相位差Θ可以表示为:
式中,ΔL为光在平晶内部的光程,对于确定的系统,ΔL为常数。背景条件下的相位差Θ0和测量条件下的相位差Θ1分别表示为:
当剪切量s大于手臂温度影响区域长度L时,在剪切区域任意位置,θT(x)和θT(x-s)至少有一个为0。因此,测量相位差可以表示为:
二、体温的计算
分别依据测量条纹和背景条纹识别相位分布,进而可以计算得到相位差的分布状况,并确定手臂温度影响区域的宽度,采用线性拟合方法确定靠近手臂皮肤处的相位差。手臂温度影响区域及光束传播过程如图2的示意图所示,在仅考虑热传导作用的情况下,手臂附近的温度随与手臂的距离变化呈近似线性分布。假定手臂横截面为椭圆,沿着与光束平行的方向建立x轴,与光束垂直的方向建立y轴,公式(Ⅷ)可以写成公式(Ⅰ)的离散形式:
式中,λ为波长,p为大气压强,M为气体分子量,K为气体折射系数,R为理想气体常数,L为手臂温度影响区域的宽度,w为手臂横截面宽度,h为手臂横截面高度,xb和yb分别为手臂表面上距离离散点(xi,h/2)最近点的横坐标和纵坐标,Δx为离散微元的长度。
设定可能的体温范围为30℃至42℃,以0.1℃为步长,遍历设定的体温范围,共121个温度值,通过公式(Ⅰ)计算每个温度值对应的相位差。结合已确定的靠近手臂皮肤处的相位差,选择使公式(Ⅱ)最小的温度值,为计算相位差与其最接近的温度值:
f(T)=|θmea-θcal,T| (Ⅱ)
式中,θmea为已确定的靠近手臂皮肤处的相位差,θcal,T为体温为T时通过公式(Ⅰ)计算得到的相位差。最后,依据选择得到的温度值和相邻温度值计算得到的相位差插值确定体温测量结果。
三、实施例
图3为实测背景条纹和存在手臂温度影响时的条纹。提取中间高度位置的干涉信号,如图4所示,图中干涉信号的空间分辨率为0.0052mm,可以确定手臂温度影响区域的宽度为3.49mm。将图4中的干涉信号均值滤波两次,滤去毛刺,将小于等于30的数值赋值为0,大于30的数值赋值为1,完成二值化操作。依据每个全1区域的中心横坐标确定干涉信号波峰位置,依据每个全0区域的中心横坐标确定干涉信号波谷位置,整理可得图5所示的干涉条纹。依据每个波峰位置的偏移以及右侧相邻周期的坐标范围,确定每个波峰位置的相位差,如图6所示。手臂皮肤位置不位于波峰或波谷,采用线性拟合方法确定靠近手臂皮肤处的相位差,如图6所示,为-0.6038rad。
测量手臂部位靠近肘部,宽度为64mm,高度为88mm,波长为632.8nm,环境温度为32.0℃。图7给出了不同温度下依据公式(Ⅰ)计算得到的相位差以及已确定的靠近手臂皮肤处相位差。由图7可知,交点位于36.0℃和36.1℃之间,36.0℃下计算的相位差为-0.5900rad,36.1℃下计算的相位差为-0.6045rad,插值可得体温为36.0952℃。
对于上述的本发明的实施例,方案中公知的具体结构及特性等常识未作过多描述;各实施例采用递进的方式描述,各实施例中所涉及到的技术特征在彼此之间不构成冲突的前提下可以相互组合,各实施例之间相同相似部分互相参见即可。以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明所述原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为落入本发明的保护范围。
Claims (2)
1.一种基于横向剪切干涉信号测量体温的方法,其特征在于,所述方法包括如下步骤:
1)采集无测量对象时的干涉信号,识别背景条纹;
2)采集存在手臂温度影响时的干涉信号,识别测量条纹;
3)测量并记录环境温度;
4)依据背景条纹和测量条纹的相位分布,确定手臂温度影响区域的相位差分布,确定手臂温度影响区域的宽度;
5)依据相位差分布,采用线性拟合方法,确定靠近手臂皮肤处的相位差;
6)依据手臂横截面的宽度和高度,结合手臂温度影响区域的宽度和环境温度,设定可能的体温范围和步长,遍历设定的体温范围,计算每个温度值对应的相位差;所述每个温度值对应的相位差依据公式(Ⅰ)计算:
式中,T为温度,T0为环境温度,λ为波长,p为大气压强,M为气体分子量,K为气体折射系数,R为理想气体常数,L为手臂温度影响区域的宽度,w为手臂横截面宽度,h为手臂横截面高度,xb和yb分别为手臂面上距离离散点(xi,h/2)最近点的横坐标和纵坐标,Δx为离散微元的长度;
7)选择计算相位差与靠近手臂皮肤处相位差最接近的温度值,并依据相邻温度值计算得到的相位差插值确定体温测量结果。
2.根据权利要求1所述的基于横向剪切干涉信号测量体温的方法,其特征在于:所述步骤7)中与测量相位差最接近的温度值依据使公式(II)最小进行选择:
f(T)=|θmea-θcal,T| (II)
式中,θmea为测量并拟合得到的靠近手臂皮肤处相位差,θcal,T为体温为T时通过公式(Ⅰ)计算得到的相位差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010805349.8A CN112082655B (zh) | 2020-08-12 | 2020-08-12 | 一种基于横向剪切干涉信号测量体温的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010805349.8A CN112082655B (zh) | 2020-08-12 | 2020-08-12 | 一种基于横向剪切干涉信号测量体温的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112082655A CN112082655A (zh) | 2020-12-15 |
CN112082655B true CN112082655B (zh) | 2022-08-12 |
Family
ID=73729156
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010805349.8A Active CN112082655B (zh) | 2020-08-12 | 2020-08-12 | 一种基于横向剪切干涉信号测量体温的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112082655B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103698022A (zh) * | 2013-12-09 | 2014-04-02 | 西南技术物理研究所 | 横向剪切干涉仪波前测量方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5904665B2 (ja) * | 2012-04-27 | 2016-04-13 | 株式会社 エフケー光学研究所 | 位相計測装置、位相計測プログラム及び位相計測方法 |
CN103776559A (zh) * | 2014-01-13 | 2014-05-07 | 南京航空航天大学 | 断层成像激光剪切干涉三维温度测量装置及温度测量方法 |
CN104484577B (zh) * | 2014-12-30 | 2017-06-16 | 华中科技大学 | 一种基于山脊能量校正的山地中带状地下目标的探测方法 |
JP6718160B2 (ja) * | 2017-02-23 | 2020-07-08 | 国立研究開発法人産業技術総合研究所 | 残留熱ひずみ測定方法、残留熱ひずみ測定装置、及びそのプログラム |
CN110243786B (zh) * | 2019-06-10 | 2024-08-13 | 佛山科学技术学院 | 基于二维剪切干涉测量待测气体参数的装置 |
-
2020
- 2020-08-12 CN CN202010805349.8A patent/CN112082655B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103698022A (zh) * | 2013-12-09 | 2014-04-02 | 西南技术物理研究所 | 横向剪切干涉仪波前测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112082655A (zh) | 2020-12-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Swojak et al. | Assessment of selected metrological properties of laser triangulation sensors | |
CN108956534B (zh) | 一种基于开放腔法珀干涉仪的折射率测量方法 | |
CN112082655B (zh) | 一种基于横向剪切干涉信号测量体温的方法 | |
CN103528499A (zh) | 形貌补偿式双光轴线位移激光干涉仪校准方法与装置 | |
JP3935469B2 (ja) | 光沢金属のストリップ上の距離を測定するための方法と装置 | |
Franco et al. | Theoretical analysis of straightness errors in coordinate measuring machines (CMM) with three linear axes | |
US5333053A (en) | Apparatus for measuring straightness | |
CN108267743B (zh) | 基于拟合插值的快速迭代水下定位方法 | |
US6501550B1 (en) | Method for improving the accuracy of measurements made using a laser interferometer | |
JP2885422B2 (ja) | 形状評価装置および形状評価方法 | |
JP3799453B2 (ja) | 内部温度分布測定装置及び内部温度分布測定方法 | |
JP3913519B2 (ja) | 走査型間隙量検出による真直度測定方法 | |
JPH0634647A (ja) | 速度計測装置 | |
JP2008286598A (ja) | 追尾式レーザ干渉計の波長推定方法 | |
CN116045824B (zh) | 一种基于白光共焦原理的高精度检测装置及方法 | |
CN113686241B (zh) | 一种高温表面线激光几何测量误差分析方法 | |
CN204988174U (zh) | 一种对比式抗干扰微动角反射镜激光干涉仪 | |
JP3078733B2 (ja) | 三次元曲率の測定方法 | |
Schmalzried et al. | A flexible and, robust 3D coordinate measurement system based on white light interferometry for calibration of industrial systems | |
JPH06331447A (ja) | 放射温度計の精度評価方法および装置 | |
Vorburger et al. | Nist surface roughness and step height calibrations | |
CN116861528A (zh) | 一种建筑工程结构变形重构方法 | |
CN117034571A (zh) | 激光干涉条纹间距三维数学模型的构建方法 | |
JPS648762B2 (zh) | ||
Brosed Dueso et al. | Thermal Error Compensation Strategy in Simultaneous Laser Multilateration for Volumetric Verification |
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 |