CN110608800B - 大气吸收带太阳辐射照度的测量方法、装置和系统 - Google Patents
大气吸收带太阳辐射照度的测量方法、装置和系统 Download PDFInfo
- Publication number
- CN110608800B CN110608800B CN201911009410.1A CN201911009410A CN110608800B CN 110608800 B CN110608800 B CN 110608800B CN 201911009410 A CN201911009410 A CN 201911009410A CN 110608800 B CN110608800 B CN 110608800B
- Authority
- CN
- China
- Prior art keywords
- thermal imager
- infrared
- solar
- black body
- absorption band
- 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
- 230000005855 radiation Effects 0.000 title claims abstract description 118
- 238000010521 absorption reaction Methods 0.000 title claims abstract description 80
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000005286 illumination Methods 0.000 claims abstract description 23
- 238000011088 calibration curve Methods 0.000 claims abstract description 21
- 238000005259 measurement Methods 0.000 claims abstract description 21
- 230000003287 optical effect Effects 0.000 claims abstract description 7
- 230000036760 body temperature Effects 0.000 claims description 27
- 238000012360 testing method Methods 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000012546 transfer Methods 0.000 claims description 4
- 238000000691 measurement method Methods 0.000 abstract description 2
- 238000012545 processing Methods 0.000 description 6
- 230000004044 response Effects 0.000 description 5
- 238000013076 uncertainty analysis Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000005457 Black-body radiation Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 241001270131 Agaricus moelleri Species 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013524 data verification Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 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
- G01J1/00—Photometry, e.g. photographic exposure meter
- G01J1/42—Photometry, e.g. photographic exposure meter using electric radiation detectors
-
- 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
- G01J1/00—Photometry, e.g. photographic exposure meter
- G01J1/42—Photometry, e.g. photographic exposure meter using electric radiation detectors
- G01J2001/4266—Photometry, e.g. photographic exposure meter using electric radiation detectors for measuring solar light
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Radiation Pyrometers (AREA)
Abstract
本发明涉及一种大气吸收带太阳辐射照度的测量方法、装置和系统，其中方法包括：使用热像仪采集太阳红外辐射图像，其中热像仪的探测器口面前安装有大气吸收带波段的滤光片；根据采集的太阳红外辐射图像计算太阳像素的平均灰度值，并基于预先标定的红外热像仪标定曲线，计算目标区域平均辐射亮度；根据计算得到的目标区域平均辐射亮度及目标面积计算太阳辐射照度。本发明通过采集大气吸收带波段的太阳红外辐射图像，对太阳的大气吸收带红外辐射亮度进行定量分析，通过计算反演出大气吸收带太阳辐射照度，解决了大气吸收带太阳辐射照度的测量问题，提供了一种具有可实操性的太阳大气吸收带波段辐射特性测量方法。
Description
技术领域
本发明涉及一种太阳辐射照度的测量方法，尤其涉及一种大气吸收带太阳辐射照度的测量方法、装置和系统。
背景技术
太阳辐射照度是指太阳辐射经过大气层的吸收、散射、反射等作用后到达固体地球表面上单位面积单位时间内的辐射能量。其单位为：瓦特/平方米(W/m2)。太阳辐射照度是定量描述和研究太阳光辐射的重要参量。
由于大气的存在，真正到达地球表面的太阳辐射能的大小要受许多因素的影响，包括太阳高度、大气质量、大气透明度、地理纬度、日照时间及海拔高度等。
由于太阳辐射照度较强，在大气窗口波段，一般设备无法拍摄其太阳红外辐射图像。因此现有技术中通常使用仿真等方式来计算太阳辐射照度。太阳表面温度5770K，针对某一红外波段对大气外太阳常数进行波段能量积分，从而获得该波段内相应的大气外太阳总辐射照度，利用大气传输软件计算太阳穿越大气层到达地面的大气透过率，进而获得地面接收到的某一波段的辐射照度。仿真方法计算的辐射照度仅供参考，能否与地面真实接收的辐射照度吻合还需要实测数据验证。因此，亟待提供一种大气吸收带太阳辐射照度的测量方法。
发明内容
本发明要解决的技术问题在于，针对现有技术中缺乏有效地太阳辐射照度的测量方法的缺陷，提供一种大气吸收带太阳辐射照度的测量方法、装置和系统。
为了解决上述技术问题，本发明第一方面提供了一种大气吸收带太阳辐射照度的测量方法，该方法包括以下步骤：
S1、使用热像仪采集太阳红外辐射图像，其中热像仪的探测器口面前安装有大气吸收带波段的滤光片；
S2、根据采集的太阳红外辐射图像计算太阳像素的平均灰度值，并基于预先标定的红外热像仪标定曲线，计算目标区域平均辐射亮度；
S3、根据计算得到的目标区域平均辐射亮度及目标面积计算太阳辐射照度。
在根据本发明所述的大气吸收带太阳辐射照度的测量方法中，优选地，所述大气吸收带波段为2.7μm～2.95μm。
在根据本发明所述的大气吸收带太阳辐射照度的测量方法中，优选地，所述步骤S3通过以下公式计算太阳辐射照度：
F＝Ls*S*C；
其中F为太阳辐射照度，Ls为目标区域平均辐射亮度，S为目标面积，C为常数。
在根据本发明所述的大气吸收带太阳辐射照度的测量方法中，优选地，所述方法还包括通过以下公式计算测量不确定度：
其中，uxi为各直接测得量xi的测量不确定度，
为被测量估计影响的传递系数。
在根据本发明所述的大气吸收带太阳辐射照度的测量方法中，优选地，所述方法还包括以下标定步骤：
将热像仪光路轴线对准黑体中心位置，调整热像仪与黑体之间的相对位置以及黑体温度，采集不同相对位置时不同黑体温度的红外图像；
读取不同相对位置获取的同一黑体温度的红外图像灰度值，建立相对位置与红外图像灰度值的第一关系式，确定不同黑体温度时热像仪探测器口面处的红外图像灰度值；
根据不同黑体温度时热像仪探测器口面处的红外图像灰度值，建立黑体温度与热像仪探测器口面处的红外图像灰度值的第二关系式，并转换成红外图像灰度值与红外辐射亮度的对应关系，得到红外热像仪标定曲线。
在根据本发明所述的大气吸收带太阳辐射照度的测量方法中，优选地，所述方法通过以下步骤采集不同相对位置时不同黑体温度的红外图像：
1)将黑体固定于试验台上；
2)将热像仪放置于滑轨上，热像仪光路轴线对准黑体中心位置；
3)设定热像仪在滑轨的多个位置，定义黑体中心位置与热像仪探测器口中心位置的距离d作为相对位置；
4)在同一相对位置d，设定不同的黑体温度T，通过热像仪采集黑体红外图像；
5)调整热像仪的位置，重复执行4)直至所需相对位置采集完毕。
在根据本发明所述的大气吸收带太阳辐射照度的测量方法中，优选地，所述步骤方法通过以下方式计算不同黑体温度时热像仪探测器口面处的红外图像灰度值：读取不同相对位置d获取的同一黑体温度的红外图像灰度值，建立相对位置d与红外图像灰度值的第一关系式f(d)＝DL(d)，确定d＝0时的红外图像灰度值DL(0)作为热像仪探测器口面处的红外图像灰度值。
本发明第二方面，提供了一种大气吸收带太阳辐射照度的测量装置，包括：
数据获取单元，用于获取热像仪采集的太阳红外辐射图像，该太阳红外辐射图像经过大气吸收带波段滤光；
辐射亮度计算单元，用于根据采集的太阳红外辐射图像计算太阳像素的平均灰度值，并基于预先标定的红外热像仪标定曲线，计算目标区域平均辐射亮度；
辐射照度计算单元，用于根据计算得到的目标区域平均辐射亮度及目标面积计算太阳辐射照度。
在根据本发明所述的大气吸收带太阳辐射照度的测量装置中，优选地，所述测量装置还包括标定单元，用于：
获取热像仪在与黑体处于不同相对位置时不同黑体温度的红外图像；
读取不同相对位置获取的同一黑体温度的红外图像灰度值，建立相对位置与红外图像灰度值的第一关系式，确定不同黑体温度时热像仪探测器口面处的红外图像灰度值；
根据不同黑体温度时热像仪探测器口面处的红外图像灰度值，建立黑体温度与热像仪探测器口面处的红外图像灰度值的第二关系式，并转换成红外图像灰度值与红外辐射亮度的对应关系，得到红外热像仪标定曲线。
本发明第三方面，提供了一种大气吸收带太阳辐射照度的测量系统，包括：如前所述的大气吸收带太阳辐射照度的测量装置；以及与之连接的热像仪，所述热像仪的探测器口面前安装有大气吸收带波段的滤光片。
实施本发明的大气吸收带太阳辐射照度的测量方法、装置和系统，具有以下有益效果：本发明通过采集大气吸收带波段的太阳红外辐射图像，对太阳的大气吸收带红外辐射亮度进行定量分析，通过计算反演出大气吸收带太阳辐射照度，解决了大气吸收带太阳辐射照度的测量问题，提供了一种具有可实操性的太阳大气吸收带波段辐射特性测量方法。
附图说明
图1为根据本发明优选实施例的大气吸收带太阳辐射照度的测量方法的流程图；
图2a-2i为采集得到的不同时刻的太阳红外测试图像；
图3a和3b分别为边界提取前后的太阳红外辐射图像；
图4是红外热像仪在2.7μm～2.95μm波段对1300℃黑体辐射的响应随测试距离的拟合曲线；
图5是根据本发明得到的红外热像仪在典型大气吸收带波段的拟合标定曲线；
图6为根据本发明优选实施例的大气吸收带太阳辐射照度的测量装置的模块框图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
请参阅图1，为根据本发明优选实施例的大气吸收带太阳辐射照度的测量方法的流程图。如图1所示，该实施例提供大气吸收带太阳辐射照度的测量方法包括以下步骤：
首先，在步骤S1中，使用热像仪采集太阳红外辐射图像，其中热像仪的探测器口面前安装有大气吸收带波段的滤光片。该大气吸收带波段为2.7μm～2.95μm。优选地，利用中波红外热像仪采集太阳红外辐射图像，每次采集数据不低于10帧。
优选地，该步骤中可以，得到不同时刻的太阳红外测试图像如图2a-2i所示，依次为11:00、12:00、13:00、14:00、14:30、15:00、15:30、16:00和16:30。
随后，在步骤S2中，根据采集的太阳红外辐射图像计算太阳像素的平均灰度值，并基于预先标定的红外热像仪标定曲线，计算目标区域平均辐射亮度。
该步骤S2具体包括以下步骤：
(1)对采集的太阳红外辐射图像进行处理，对太阳红外辐射图像中的太阳进行边界提取，确定太阳在热像仪图像中的像素位置，即太阳像素，如图3a和3b分别为边界提取前后的太阳红外辐射图像。划定边界的方法有多种，可以对比这些方法计算出的像素数量，选取准确度高的处理方法，比如边缘检测法、阈值法、数学形态法等等，由于太阳与背景相比为超强目标，边界清晰，选用阈值法即可。
(2)计算太阳像素的平均灰度值。
(3)基于预先标定的红外热像仪标定曲线，计算目标区域平均辐射亮度。预先标定的红外热像仪标定曲线为红外图像灰度值与红外辐射亮度的对应关系，通过平均灰度值可以计算出对应的红外辐射亮度，即上述目标区域平均辐射亮度。最后，在步骤S3中，根据计算得到的目标区域平均辐射亮度及目标面积计算太阳辐射照度。
具体地，该步骤S3通过以下公式计算太阳辐射照度：
F＝Ls*S*C；
其中F为太阳辐射照度，单位：W/m2；Ls为目标区域平均辐射亮度，单位：W/(m2·sr；S为目标面积，在本发明中为太阳的面积，单位：m2；C为常数，取值为4.44444444444445E-23，单位：sr·m-2。
在本发明更优选的实施例中，上述方法还包括不确定度分析步骤，针对上述得到的大气吸收带太阳辐射照度进行不确定分析。具体采用标准不确定度B类评定方法，各测量之间相互独立，测量不确定度计算公式可表示为：
其中，uxi为各直接测得量xi的测量不确定度，
为被测量估计影响的传递系数；N为直接测得量的数量。辐射定标中，直接测得量为红外图像灰度值x1；(b)太阳面积计算中，直接测得量为像素数x2。
大气吸收带太阳红外辐射照度测量不确定度来源主要有以下两个方面：(a)辐射定标的测量不确定度；(b)太阳面积的不确定度。
本发明的大气吸收带太阳辐射照度的测量方法还对大气吸收带红外辐射亮度定标进行了研究。
红外波段包含三个大气窗口，分别为2μm～2.5μm、3μm～5μm、8μm～14μm，在大气窗口内，大气对红外线的吸收较少，因此大多数红外系统都选用大气窗口波段。
红外热像仪是一种用来探测目标物体的红外辐射，并通过光电转换、电信号处理等手段，将目标物体的温度分布瞬时可视化的技术设备。红外热像仪在出厂时经过严格的标定，其标定过程为将黑体设置到几个不同的温度，红外热像仪分别读出不同温度对应的电压值或灰度，不同温度对应的电压值在其标定范围内为线性，这样就建立了温度与电压值的关系。在测量时红外热像仪记录的是电压或是灰度，通过标定时建立的不同温度下黑体与热像仪对应电压灰度的关系算出目标和背景辐射温度并显示在显示器上。
通常对红外热像仪标定时，选定的是大气窗口，即中波3μm～5μm、长波8μm～14μm，热像仪响应(图像灰阶)在一定条件下与黑体温度为线性关系，热像仪与黑体的测试距离通常都在10m以内，属于近距离测试，大气传输对标定的影响一般忽略不计。但是，在大气吸收带波段，由于红外辐射亮度随着距离衰减严重，使得现有针对大气窗口波段的红外热像仪标定方式已经不适用。因此，本发明中研究了一种大气吸收带红外辐射亮度定标方法，可以去除测试距离的影响，获得红外热像仪在典型大气吸收带波段的标定曲线，将大气吸收带波段的红外热像数据转变成可以定量分析的数据。
因此，在本发明更优选的实施例中，上述方法还包括标定步骤：首先，执行数据获取步骤，将热像仪光路轴线对准黑体中心位置，调整热像仪与黑体之间的相对位置以及黑体温度，采集不同相对位置时不同黑体温度的红外图像。本发明中的热像仪，又可称为红外热像仪。
优选地，该数据获取步骤包括：
(1)将黑体固定于试验台上；
(2)将热像仪放置于滑轨上，热像仪光路轴线对准黑体中心位置；其中热像仪的探测器口面前安装有2.7μm～2.95μm波段的大气吸收带波段滤光片。
(3)设定热像仪在滑轨的多个位置，定义黑体中心位置与热像仪探测器口中心位置的距离d作为相对位置；
(4)在同一相对位置d，设定不同的黑体温度T，如T1、T2、T3、…、Tm，通过热像仪采集黑体红外图像；
(5)调整热像仪的位置，重复执行(4)直至所需相对位置采集完毕。由此可以获得不同相对位置d1、d2、d3、…、dn时不同黑体温度T1、T2、T3、…、Tm的红外图像。相对位置调整的范围优选为0～5m。
本发明选择的用于定标方法的典型大气吸收带为2.7μm～2.95μm，以研究2.7μm～2.95μm波段红外热像仪测试距离与黑体红外图像灰度值的关系，建立测试距离与红外图像灰度值的关系，通过不同距离的灰阶数据，反演出目标在2.7μm～2.95μm波段达到热像仪镜头口面的探测器响应灰阶。通过上述步骤采集不同测试距离(即相对位置)、不同黑体温度下的黑体的红外图像，形成不同测试要素下黑体红外图像灰阶表，如表1所示。
表1
上标中DLij为第j个相对位置dj处第i个黑体温度Ti时采集的黑体的红外图像的灰阶即灰度值。
随后，执行第一数据处理步骤，读取不同相对位置获取的同一黑体温度的红外图像灰度值，建立相对位置与红外图像灰度值的第一关系式，确定不同黑体温度时热像仪探测器口面处的红外图像灰度值。
具体地，该步骤中读取不同相对位置d获取的同一黑体温度的红外图像灰度值，建立相对位置d与红外图像灰度值的第一关系式f(d)＝DL(d)，确定d＝0m时的红外图像灰度值DL(0)作为热像仪探测器口面处的红外图像灰度值。例如，由DL11、DL12、DL13、DL14、…、DL1n拟合出在黑体温度为T1下红外图像灰阶随相对距离d的变化曲线，进而得出黑体放置于热像仪探测器口面时红外图像灰阶，即d＝0时的DL值。由此可以得到不同黑体温度T1、T2…、Tm时热像仪探测器口面处的红外图像灰度值。图4是红外热像仪在2.7μm～2.95μm波段对1300℃黑体辐射的响应随测试距离的拟合曲线，拟合函数为一阶衰减指数函数。
最后，执行第二数据处理步骤，根据不同黑体温度T时热像仪探测器口面处的红外图像灰度值DL(0)，又可记为DL(0，T)，建立各个黑体温度T与热像仪探测器口面处的红外图像灰度值DL(0，T)的第二关系式g(T)＝DL(0，T)，再转换成红外图像灰度值DL(0，T)与红外辐射亮度R的对应关系，即红外热像仪标定曲线。其中红外辐射亮度R与黑体温度T成四次平方关系。在红外探测器线性响应区，红外图像灰度值与红外辐射亮度R构成线性关系。
请参阅图5，是根据本发明得到的红外热像仪在典型大气吸收带波段的拟合标定曲线，拟合函数为一阶衰减指数函数。该曲线为红外图像灰度值与温度的第二关系式，后续再通过转换变成红外图像灰度值与红外辐射亮度的对应关系，即可得到红外热像仪标定曲线。
本发明利用红外热像仪在2.7μm～2.95μm波段的拟合标定曲线，能够成功应用于太阳在大气吸收带波段红外特性定量计算。
基于同一发明构思，本发明实施例还提供一种大气吸收带太阳辐射照度的测量装置。请参阅图6，为根据本发明优选实施例的大气吸收带太阳辐射照度的测量装置的模块框图。如图6所示，该实施例提供的装置600包括：数据获取单元601、辐射亮度计算单元602和辐射照度计算单元603。
数据获取单元601，用于获取热像仪采集的太阳红外辐射图像，该太阳红外辐射图像经过大气吸收带波段滤光。该装置可与热像仪连接，由该数据获取单元601获取热像仪采集的太阳红外辐射图像。
辐射亮度计算单元602，用于根据采集的太阳红外辐射图像计算太阳像素的平均灰度值，并基于预先标定的红外热像仪标定曲线，计算目标区域平均辐射亮度。该辐射亮度计算单元的处理过程与前述方法中步骤S2一致，在此不再赘述。
辐射照度计算单元603，用于根据计算得到的目标区域平均辐射亮度及目标面积计算太阳辐射照度。该辐射照度计算单元603的处理过程与前述方法中步骤S3一致，在此不再赘述。
优选地，该装置600还包括标定单元，用于执行以下操作：
a、获取热像仪在与黑体处于不同相对位置时不同黑体温度的红外图像；
b、读取不同相对位置获取的同一黑体温度的红外图像灰度值，建立相对位置与红外图像灰度值的第一关系式，确定不同黑体温度时热像仪探测器口面处的红外图像灰度值；
c、根据不同黑体温度时热像仪探测器口面处的红外图像灰度值，建立黑体温度与热像仪探测器口面处的红外图像灰度值的第二关系式，并转换成红外图像灰度值与红外辐射亮度的对应关系，得到红外热像仪标定曲线。
优选地，该装置600还包括不确定度分析单元，用于通过以下公式计算测量不确定度：
其中，uxi为各直接测得量xi的测量不确定度，
为被测量估计影响的传递系数。
本发明实施例还提供了一种大气吸收带太阳辐射照度的测量系统，包括：如前所述的大气吸收带太阳辐射照度的测量装置，以及与之连接的：热像仪，所述热像仪的探测器口面前安装有大气吸收带波段的滤光片。
本发明还对太阳在热像仪视场中的位置可能带来的影响进行了研究，结果表明，太阳匀速穿过热像仪视场过程中大气吸收带太阳红外辐射照度随时刻的存在一定的起伏，当太阳在热像仪视场中心时，起伏较小；因此本发明在采集太阳红外测试图像时，优选将太阳放置于热像仪视场中心。
综上所述，本发明利用中波红外热像仪(波段2.5μm～5μm)加装窄带滤光片(滤光片波段位于大气吸收带)对太阳进行测试，获得太阳在大气吸收带波段的红外热图像。并开展对红外热像仪在大气吸收带波段的定标，获取红外热像仪在大气吸收带波段的标定曲线，将太阳在大气吸收带波段的红外热图像转变成可以定量分析的数据，通过数据分析，给出太阳在大气吸收带波段的红外辐射特性。本发明还对测量方法进行了不确定度分析，可应用于某些领域采用吸收带工作波段的传感器的探测能力标定。
应该理解地是，本发明中大气吸收带太阳辐射照度的测量方法、装置和系统的原理相同，因此对大气吸收带太阳辐射照度的测量方法的实施例的详细阐述也适用于大气吸收带太阳辐射照度的测量装置和系统。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (4)
1.一种大气吸收带太阳辐射照度的测量方法，其特征在于，该方法包括以下步骤：
S1、使用热像仪采集太阳红外辐射图像，其中热像仪的探测器口面前安装有大气吸收带波段的滤光片；
S2、根据采集的太阳红外辐射图像计算太阳像素的平均灰度值，并基于预先标定的红外热像仪标定曲线，计算目标区域平均辐射亮度；
S3、根据计算得到的目标区域平均辐射亮度及目标面积计算太阳辐射照度；
所述大气吸收带波段为2.7μm～2.95μm；
所述步骤S3通过以下公式计算太阳辐射照度：
F＝Ls*S*C；
其中F为太阳辐射照度，Ls为目标区域平均辐射亮度，S为目标面积，C为常数；
所述方法还包括以下标定步骤：
将热像仪光路轴线对准黑体中心位置，调整热像仪与黑体之间的相对位置以及黑体温度，采集不同相对位置时不同黑体温度的红外图像；
读取不同相对位置获取的同一黑体温度的红外图像灰度值，建立相对位置与红外图像灰度值的第一关系式，确定不同黑体温度时热像仪探测器口面处的红外图像灰度值；
根据不同黑体温度时热像仪探测器口面处的红外图像灰度值，建立黑体温度与热像仪探测器口面处的红外图像灰度值的第二关系式，并转换成红外图像灰度值与红外辐射亮度的对应关系，得到红外热像仪标定曲线；
所述方法通过以下步骤采集不同相对位置时不同黑体温度的红外图像：
1)将黑体固定于试验台上；
2)将热像仪放置于滑轨上，热像仪光路轴线对准黑体中心位置；
3)设定热像仪在滑轨的多个位置，定义黑体中心位置与热像仪探测器口中心位置的距离d作为相对位置；
4)在同一相对位置d，设定不同的黑体温度T，通过热像仪采集黑体红外图像；
5)调整热像仪的位置，重复执行4)直至所需相对位置采集完毕；
所述方法中通过以下方式计算不同黑体温度时热像仪探测器口面处的红外图像灰度值：读取不同相对位置d获取的同一黑体温度的红外图像灰度值，建立相对位置d与红外图像灰度值的第一关系式f(d)＝DL(d)，确定d＝0时的红外图像灰度值DL(0)作为热像仪探测器口面处的红外图像灰度值。
2.根据权利要求1所述的大气吸收带太阳辐射照度的测量方法，其特征在于，所述方法还包括通过以下公式计算测量不确定度：
其中，uxi为各直接测得量xi的测量不确定度，
为被测量估计影响的传递系数。
3.一种大气吸收带太阳辐射照度的测量装置，其特征在于，包括：
数据获取单元，用于获取热像仪采集的太阳红外辐射图像，该太阳红外辐射图像经过大气吸收带波段滤光；所述大气吸收带波段为2.7μm～2.95μm；
辐射亮度计算单元，用于根据采集的太阳红外辐射图像计算太阳像素的平均灰度值，并基于预先标定的红外热像仪标定曲线，计算目标区域平均辐射亮度；
辐射照度计算单元，用于根据计算得到的目标区域平均辐射亮度及目标面积计算太阳辐射照度；
所述辐射照度计算单元通过以下公式计算太阳辐射照度：
F＝Ls*S*C；
其中F为太阳辐射照度，Ls为目标区域平均辐射亮度，S为目标面积，C为常数；
所述测量装置还包括标定单元，用于：
获取热像仪在与黑体处于不同相对位置时不同黑体温度的红外图像；
读取不同相对位置获取的同一黑体温度的红外图像灰度值，建立相对位置与红外图像灰度值的第一关系式，确定不同黑体温度时热像仪探测器口面处的红外图像灰度值；
根据不同黑体温度时热像仪探测器口面处的红外图像灰度值，建立黑体温度与热像仪探测器口面处的红外图像灰度值的第二关系式，并转换成红外图像灰度值与红外辐射亮度的对应关系，得到红外热像仪标定曲线；
其中不同相对位置时不同黑体温度的红外图像通过以下步骤采集：
1)将黑体固定于试验台上；
2)将热像仪放置于滑轨上，热像仪光路轴线对准黑体中心位置；
3)设定热像仪在滑轨的多个位置，定义黑体中心位置与热像仪探测器口中心位置的距离d作为相对位置；
4)在同一相对位置d，设定不同的黑体温度T，通过热像仪采集黑体红外图像；
5)调整热像仪的位置，重复执行4)直至所需相对位置采集完毕；
其中通过以下方式确定不同黑体温度时热像仪探测器口面处的红外图像灰度值：读取不同相对位置d获取的同一黑体温度的红外图像灰度值，建立相对位置d与红外图像灰度值的第一关系式f(d)＝DL(d)，确定d＝0时的红外图像灰度值DL(0)作为热像仪探测器口面处的红外图像灰度值。
4.一种大气吸收带太阳辐射照度的测量系统，其特征在于，包括：
权利要求3所述的大气吸收带太阳辐射照度的测量装置；以及与之连接的：
热像仪，所述热像仪的探测器口面前安装有大气吸收带波段的滤光片。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911009410.1A CN110608800B (zh) | 2019-10-23 | 2019-10-23 | 大气吸收带太阳辐射照度的测量方法、装置和系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911009410.1A CN110608800B (zh) | 2019-10-23 | 2019-10-23 | 大气吸收带太阳辐射照度的测量方法、装置和系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110608800A CN110608800A (zh) | 2019-12-24 |
| CN110608800B true CN110608800B (zh) | 2021-11-05 |
Family
ID=68894929
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911009410.1A Active CN110608800B (zh) | 2019-10-23 | 2019-10-23 | 大气吸收带太阳辐射照度的测量方法、装置和系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110608800B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113588095B (zh) * | 2021-08-17 | 2022-08-09 | 北京环境特性研究所 | 基于帧累积的点源目标红外辐射定量处理方法及装置 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3970848A (en) * | 1975-08-22 | 1976-07-20 | Calspan Corporation | Profile technique for calibrating infrared thermal imaging systems |
| CN102692271A (zh) * | 2012-06-14 | 2012-09-26 | 中国气象科学研究院 | 基于天空可见光图像的太阳直接辐射强度测量方法和装置 |
| CN108168710A (zh) * | 2017-12-28 | 2018-06-15 | 福建农林大学 | 一种基于遥感技术的城区热岛效应评估方法 |
| CN109655161A (zh) * | 2018-11-15 | 2019-04-19 | 上海无线电设备研究所 | 基于红外热像仪的目标红外积分辐射强度测试方法及设备 |
-
2019
- 2019-10-23 CN CN201911009410.1A patent/CN110608800B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3970848A (en) * | 1975-08-22 | 1976-07-20 | Calspan Corporation | Profile technique for calibrating infrared thermal imaging systems |
| CN102692271A (zh) * | 2012-06-14 | 2012-09-26 | 中国气象科学研究院 | 基于天空可见光图像的太阳直接辐射强度测量方法和装置 |
| CN108168710A (zh) * | 2017-12-28 | 2018-06-15 | 福建农林大学 | 一种基于遥感技术的城区热岛效应评估方法 |
| CN109655161A (zh) * | 2018-11-15 | 2019-04-19 | 上海无线电设备研究所 | 基于红外热像仪的目标红外积分辐射强度测试方法及设备 |
Non-Patent Citations (1)
| Title |
|---|
| 《红外热像仪的辐射定标和测温误差分析》;陆子凤;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20101115;正文第29-42页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110608800A (zh) | 2019-12-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102183237B (zh) | 一种地基双波段云高测量的装置和方法 | |
| CN102589705B (zh) | 热成像仪单黑体温控mrtd野外在线自动检测装置及检测方法 | |
| CN100464167C (zh) | 一种红外测温实时校准的方法及装置 | |
| KR101116462B1 (ko) | 원격탐사 데이터를 이용한 수질 모니터링 방법 및 이를 이용한 수질 모니터링 장치 | |
| CN109655162A (zh) | 红外热像仪测温校正系统及方法 | |
| CN102853916B (zh) | 一种用于对煤垛表面进行远距离红外测温的方法及系统 | |
| Lübcke et al. | On the absolute calibration of SO 2 cameras | |
| CN102081039A (zh) | 一种环境可控的作物营养水分高光谱图像检测装置 | |
| CN110095193B (zh) | 一种红外热像仪噪声等效温差测试方法及系统 | |
| CN106124062A (zh) | 一种基于历史数据的红外测温自动补偿方法 | |
| CN110595628A (zh) | 大气吸收带红外辐射亮度定标方法、装置和系统 | |
| CN107884005B (zh) | 一种模拟人眼对光环境感知的光学测量系统与测量方法 | |
| US7364355B2 (en) | Method and apparatus for obtaining a temperature measurement using an InGaAs detector | |
| CN103335717B (zh) | 一种基于变积分模式的红外热像仪高精度抗温漂测温方法 | |
| CN105352988A (zh) | 一种建筑物外墙保温性能评估系统及方法 | |
| CN106441808A (zh) | 一种热红外高光谱成像仪盲元检测装置和方法 | |
| CN109974854A (zh) | 一种框幅式fpi高光谱图像的辐射校正方法 | |
| CN109060731B (zh) | 红外光学系统光谱透过率测试装置及方法 | |
| CN112504463A (zh) | 一种测温系统及其测温方法 | |
| CN109060144A (zh) | 红外热像仪netd自动测试的方法 | |
| CN110608800B (zh) | 大气吸收带太阳辐射照度的测量方法、装置和系统 | |
| CN111829665A (zh) | 一种基于用人体作为虚拟黑体的体温测量的方法、装置及存储介质 | |
| CN111044153A (zh) | 一种图谱关联系统红外光谱的非线性定标方法及装置 | |
| CN113049118A (zh) | 一种红外焦平面探测器闪元测试装置及测试方法 | |
| RU2424496C2 (ru) | Способ дистанционного измерения температурного поля |
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 |