CN117112966B - 一种校正红外遥感器非线性响应系数的方法、系统和设备 - Google Patents
一种校正红外遥感器非线性响应系数的方法、系统和设备 Download PDFInfo
- Publication number
- CN117112966B CN117112966B CN202311264840.4A CN202311264840A CN117112966B CN 117112966 B CN117112966 B CN 117112966B CN 202311264840 A CN202311264840 A CN 202311264840A CN 117112966 B CN117112966 B CN 117112966B
- Authority
- CN
- China
- Prior art keywords
- response coefficient
- nonlinear
- remote sensor
- infrared remote
- linear response
- 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
- 230000004044 response Effects 0.000 title claims abstract description 284
- 238000000034 method Methods 0.000 title claims abstract description 69
- 230000005855 radiation Effects 0.000 claims abstract description 35
- 238000012937 correction Methods 0.000 claims abstract description 23
- 238000012935 Averaging Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 abstract description 4
- 238000004364 calculation method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 238000007430 reference method Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000011088 calibration curve Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000001747 pupil Anatomy 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- 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/90—Testing, inspecting or checking operation of radiation pyrometers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Operations Research (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
本发明涉及遥感器技术领域，具体涉及一种校正红外遥感器非线性响应系数的方法、系统和设备。方法包括：构建红外遥感器的红外辐射定标方程组，求解得到线性响应系数和非线性响应系数，并根据线性响应系数和非线性响应系数计算得到非线性形状参数；获取预设温度下红外遥感器的初始观测值，利用非线性形状参数对非线性响应系数红外遥感器的非线性响应系数进行迭代，得到收敛后的线性响应系数和收敛后的非线性响应系数。本发明解决了传统校正方法中将非线性响应系数视为保持不变的常量、忽略线性响应系数与非线性响应系数之间的动态变化关系从而导致红外遥感器定标精度低的问题。
Description
技术领域
本发明涉及遥感器技术领域，具体涉及一种校正红外遥感器非线性响应系数的方法、系统和设备。
背景技术
随着定量遥感技术的发展，需要利用气象卫星进行持续的观测，研究人员对卫星遥感器的数据质量也提出了越来越高的要求，因此直接影响遥感数据准确性的定标精度就显得尤为重要，其中，影响定标精度的主要因素之一是遥感器系统的非线性校正。
热红外通道辐射定标采用非线性二次项形式，基于星上黑体和冷空观测，建立仪器响应计数值与入瞳辐射间的关系。在红外遥感器中，非线性响应系数一般是在发射前的实验室定标中计算得到的，并且非线性响应系数被视为不受其他条件影响的、在轨期间也保持不变的常量，也就是在遥感器的在轨定标中通常只校正线性响应系数，而非线性响应系数是直接应用实验室得到的数据。
由于红外遥感器的线性响应系数和非线性响应系数都会受到工作环境的热变化和其他因素的影响，即非线性响应系数与线性响应系数之间存在动态变化关系，而非保持不变，导致红外遥感器在非线性响应的计算方面也将受到影响，因此对红外遥感器非线性校正响应有必要进行优化。
发明内容
针对现有技术的以上缺陷或改进需求，本发明提供了一种校正红外遥感器非线性响应系数的方法、系统和设备，其目的在于优化红外遥感器非线性校正得到的非线性响应系数和线性响应系数的计算，由此解决传统校正方法中将非线性响应系数粗略视为保持不变的常量、忽略线性响应系数与非线性响应系数之间的动态变化关系从而导致红外遥感器定标精度低的技术问题。
为实现上述目的，按照本发明的一个方面，提供了一种校正红外遥感器非线性响应系数的方法，所述方法包括：
构建红外遥感器的红外辐射定标方程组，求解得到线性响应系数和非线性响应系数/>，并根据所述线性响应系数/>和所述非线性响应系数/>计算得到非线性形状参数/>；
所述红外遥感器通过预设温度黑体得到初始观测值，利用所述非线性形状参数对所述红外遥感器的非线性响应系数进行迭代，得到收敛后的线性响应系数和收敛后的非线性响应系数。
作为对上述方案进一步的完善和补充，本发明还包括以下附加技术特征。
优选地，构建所述红外辐射定标方程组具体包括：
利用实验室预设温度的定标黑体以及所述红外遥感器接收所述定标黑体的有效辐射亮度和所述红外遥感器输出的响应信号计数值/>，构建用于求解线性响应系数和非线性响应系数/>的超定方程组：
；
其中：表示定标黑体的预设温度，单位为开尔文（K），/>表示观测/>温度的黑体时，去除冷黑体辐射后红外遥感器接收的有效净辐射亮度，/>表示观测/>温度的黑体时，去除冷黑体观测计数值后红外遥感器的输出信号计数值，/>为线性响应系数，/>为非线性响应系数。
优选地，所述红外辐射定标方程组的求解方法包括：
使用至少3个所述超定方程组通过最小二乘法求解非线性响应系数和线性响应系数；
得到组结果，对非线性响应系数和线性响应系数结果求平均值后，分别作为所述线性响应系数/>和所述非线性响应系数/>。
优选地，所述非线性形状参数与所述线性响应系数/>和所述非线性响应系数/>的关系为：/>。
优选地，所述初始观测值包括：所述红外遥感器观测所述预设温度黑体接收得到的信号计数值及所述红外遥感器输出的所述预设温度黑体的辐射亮度/>。
优选地，所述利用所述非线性形状参数对所述非线性响应系数红外遥感器的非线性响应系数进行迭代的具体方法包括：
根据所述预设温度黑体得到的信号计数值及所述预设温度黑体的辐射亮度计算得到线性响应系数初始值/>非线性响应系数初始值/>；
修正初始辐射亮度中的非线性响应系数，获得修正后的辐射亮度/>；
根据修正后的辐射亮度和所述预设温度黑体得到的信号计数值/>计算出线性响应系数修正值/>；
由所述线性响应系数初始值和所述线性响应系数修正值/>求得线性响应系数更新值/>；
判断修正后的线性响应系数是否符合收敛性；
若满足收敛条件，得到收敛后的线性响应系数和收敛后的非线性响应系数；若不满足收敛条件，将所述线性响应系数更新值赋值给所述线性响应系数初始值/>，使用更新后的线性响应系数初始值参与迭代，继续进行修正，直至收敛。
优选地，所述判断修正后的线性响应系数是否符合收敛性的方法包括：
；
其中：；/>为预设收敛阈值。
优选地，所述收敛后的非线性响应系数的计算方法包括：
通过收敛后的线性响应系数和所述非线性形状参数求出收敛后的非线性响应系数。
优选地，所述预设收敛阈值与所述红外遥感器的定标精度负相关。
第二方面，本发明还提供了一种校正星载红外遥感器非线性的迭代设备，设备包括：一个或多个处理器；存储装置，用于存储一个或多个程序，当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如第一方面所述的校正红外遥感器非线性响应系数的方法。
总体而言，通过本发明所构思的以上技术方案与现有技术相比，具有如下有益效果：
本发明提出一种能够对星载红外遥感器非线性进行精确动态校正的迭代算法，提取非线性形状参数作为表征遥感器非线性特性的物理参数，结合两点定标所得线性响应系数和非线性形状参数/>计算非线性响应系数，通过多次迭代修正非线性响应的影响后将线性响应系数和非线性响应系数均收敛到准确值，为各类型遥感器在轨定标中的动态非线性校正过程提供新的参考方法。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对本发明实施例中所需要使用的附图作简单地介绍。显而易见地，下面所描述的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本实施例一提供的校正红外遥感器非线性响应系数的方法流程图；
图2是本实施例一提供的校正红外遥感器非线性响应系数的方法示意图；
图3是本实施例二提供的校正红外遥感器非线性的迭代系统的结构示意图；
图4是本实施例三中提供的校正红外遥感器非线性的迭代设备示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
实施例一
本实施例一提供一种校正红外遥感器非线性响应系数的方法，所述方法包括以下步骤，如图1所示：
S101：构建红外遥感器的红外辐射定标方程组，求解得到线性响应系数和非线性响应系数/>，并根据所述线性响应系数/>和所述非线性响应系数/>计算得到非线性形状参数/>。
构建基于实验室观测结果的辐射定标方程组，在最小二乘准则下求解线性和非线性响应系数，并根据线性和非线性响应系数计算非线性形状参数，具体过程为：
S201：如图2所示，利用实验室定标中红外遥感器观测不同温度（）定标黑体时红外遥感器接收的有效辐射亮度/>和其响应信号计数值，构造用于求解线性和非线性响应系数的超定方程组，具体如下：
（1）；
式（1）中，表示定标黑体的预设温度，单位为开尔文（K），/>表示观测/>温度的黑体时，去除冷黑体辐射后红外遥感器接收的有效净辐射亮度，/>表示观测/>温度的黑体时，去除冷黑体观测计数值后红外遥感器的输出信号计数值，/>为线性响应系数，/>为非线性响应系数。
S202：由于式（1）给出的方程组中包含两个未知量，因此可用3个及以上方程式通过最小二乘法求解出和/>，共可得到/>组/>和/>，求出所有组数据的平均值作为线性响应系数/>和非线性响应系数/>的实验室定标结果。
S203：由于微波遥感器非线性校正方法中对定标增益的定义为：
（2）；
式（1）中，分别表示热黑体和冷黑体的亮温，/>则分别是热黑体和冷黑体对应的遥感器输出信号计数值。
红外遥感器定标一般使用辐射亮度进行计算，因此将式（2）中的亮温替换成辐射亮度可以得到适用于红外遥感器的定标增益：
（3）；
在不含非线性响应系数的理想情况下，线性响应系数/>可以用两点定标方法计算得到：
（4）；
因此在不含非线性响应系数的理想情况下可以用线性响应系数/>代替定标增益/>。
S204：提取非线性形状参数，将非线性形状参数/>用非线性响应系数/>表示为：
（5）；
非线性形状参数代表了遥感器定标曲线的形状，是一个不随环境温度等外部因素而变化的参数。
S102：所述红外遥感器通过预设温度黑体得到初始观测值，利用所述非线性形状参数对所述非线性响应系数红外遥感器的非线性响应系数进行迭代，得到收敛后的线性响应系数和收敛后的非线性响应系数。
迭代算法修正系数包括：获取预设温度下的初始观测值，根据S101中得到的非线性形状参数，利用所述非线性形状参数/>对所述红外遥感器的非线性响应系数进行迭代后，得到收敛后的线性响应系数和收敛后的非线性响应系数，最终得到精确的定标系数，具体过程如图2所示：
S301：根据红外遥感器观测预设温度黑体得到的信号计数值及预设温度黑体的辐射亮度/>，用两点定标法得到线性响应系数初始值/>：
（6）；
S302：根据式（5），在不含非线性响应系数的理想情况下，通过线性响应系数初始值/>和步骤S101中得到的非线性形状参数/>求出非线性响应系数初始值/>：
（7）；
S303：修正初始辐射亮度中的非线性响应，即从/>中减去根据非线性响应系数初始值/>和初始/>计算得到的非线性响应，以获得修正后的辐射亮度/>：
（8）；
S304：根据修正后的辐射亮度和初始/>计算出线性响应系数修正值/>：
（9）；
S305：为了避免步骤S304对线性响应系数的修正过度，由所述线性响应系数初始值和所述线性响应系数修正值/>求得线性响应系数更新值/>，本实施例一中将线性响应系数初始值/>和线性响应系数修正值/>之和平均后求得线性响应系数更新值/>：
（10）；
S306：根据步骤S304得到的线性响应系数修正值和步骤S305得到的线性响应系数更新值/>，判断修正前后线性响应系数的偏差与预设收敛阈值/>的大小关系：
（11）；
式（11）中的为用于判断修正后的线性响应系数是否收敛的预设收敛阈值，预设收敛阈值/>的设置取决于红外遥感器定标系数所需精度。
如果判断式（11）不成立，则认为修正后的线性响应系数仍存在较大偏差，不满足收敛条件，将所述线性响应系数更新值赋值给所述线性响应系数初始值/>：
（12）；
使用根据式（12）更新后的线性响应系数初始值重新参与迭代，继续从步骤S302开始进行修正，直至收敛。
S307：如果判断步骤S306中的式（11）成立，则可以认为修正前后的线性响应系数差异足够小，满足了收敛性判断条件，可以得到收敛后的线性响应系数和收敛后的非线性响应系数，并且收敛后的线性响应系数和收敛后的非线性响应系数较为准确。且步骤S303得到的修正后的辐射不含非线性响应系数，也就是此时得到的线性响应系数更新值/>比最开始参与迭代计算的初始值更为准确，根据步骤S306中的式（12）将线性响应系数初始值更新为此时的线性响应系数更新值/>，再结合步骤S302中的式（7）算出更新后的非线性响应系数初始值/>，最终得到线性响应系数和非线性响应系数的准确值，分别将线性响应系数初始值/>作为收敛后的线性响应系数，将非线性响应系数初始值/>作为收敛后的线性响应系数，用收敛后的线性响应系数和收敛后的线性响应系数对红外遥感器定标，得到的结果更为精确。
收敛后的线性响应系数和收敛后的非线性响应系数是辐射定标系数，红外遥感器的定标方程一般用非线性（二次项）形式，如公式（1），表示辐射亮度I与信号计数值DN值之间的关系。
S102中提到的利用非线性形状参数进行迭代用于修正在轨定标阶段的非线性响应系数的算法，可以同样适用于红外遥感器和微波遥感器。
为了得到实验室定标阶段的线性响应系数和非线性响应系数/>，结合本发明实施例，还存在一种优选的实现方案，具体的，构建所述红外辐射定标方程组具体包括：
利用实验室预设温度的定标黑体以及所述红外遥感器接收所述定标黑体的有效辐射亮度和所述红外遥感器输出的响应信号计数值/>，构建用于求解线性响应系数和非线性响应系数/>的超定方程组：
；
其中：表示定标黑体的预设温度，单位为开尔文（K），/>表示观测/>温度的黑体时，去除冷黑体辐射后红外遥感器接收的有效净辐射亮度，/>表示观测/>温度的黑体时，去除冷黑体观测计数值后红外遥感器的输出信号计数值，/>为线性响应系数，/>为非线性响应系数。
为了得到实验室定标阶段的线性响应系数和非线性响应系数/>，结合本发明实施例，还存在一种优选的实现方案，具体的，所述红外辐射定标方程组的求解方法包括：
使用至少3个所述超定方程组通过最小二乘法求解非线性响应系数和线性响应系数；
得到组结果，对非线性响应系数和线性响应系数结果求平均值后，分别作为所述线性响应系数/>和所述非线性响应系数/>。
在实验室定标阶段获得线性响应系数和非线性响应系数/>，是为了通过线性响应系数/>和非线性响应系数/>求出代表遥感器特性的非线性形状参数/>，在卫星轨迹运行中线性响应系数/>和非线性响应系数/>会发生变化，不能直接用实验室定标结果，而线性响应系数/>和非线性响应系数/>可看作常数能够参与到在轨定标的计算中。所以在红外遥感器在轨运行阶段，由于环境发生变化，在轨定标的线性响应系数/>和非线性响应系数与实验室定标得到的线性响应系数/>和非线性响应系数/>不同，因此需要利用非线性形状参数/>参与迭代计算并得到在轨定标的线性响应系数/>和非线性响应系数/>。
为了提取得到非线性形状参数，结合本发明实施例，还存在一种优选的实现方案，具体的，如图2所示，所述非线性形状参数/>与所述线性响应系数/>和所述非线性响应系数/>的关系为：/>。
结合本发明实施例，还存在一种优选的实现方案，具体的，如图2所示，所述初始观测值包括：所述红外遥感器观测所述预设温度黑体接收得到的信号计数值及所述红外遥感器输出的所述预设温度黑体的辐射亮度/>。
为了对星载红外遥感器非线性进行精确动态校正，结合本发明实施例，还存在一种优选的实现方案，具体的，如图2所示，所述利用所述非线性形状参数对所述非线性响应系数红外遥感器的非线性响应系数进行迭代的具体方法包括：
根据所述预设温度黑体得到的信号计数值及所述预设温度黑体的辐射亮度计算得到线性响应系数初始值/>非线性响应系数初始值/>；
修正初始辐射亮度中的非线性响应系数，获得修正后的辐射亮度/>；
根据修正后的辐射亮度和所述预设温度黑体得到的信号计数值/>计算出线性响应系数修正值/>；
由所述线性响应系数初始值和所述线性响应系数修正值/>求得线性响应系数更新值/>；
判断修正后线性响应系数是否符合收敛性；
若满足收敛条件，得到收敛后的线性响应系数和收敛后的非线性响应系数；若不满足收敛条件，将所述线性响应系数更新值赋值给所述线性响应系数初始值/>，使用替换后的线性响应系数初始值参与迭代，继续进行修正，直至收敛。
为了判断每次更新得到的线性响应系数是否满足收敛条件，结合本发明实施例，还存在一种优选的实现方案，具体的，如图2所示，所述判断修正后的线性响应系数是否符合收敛性的方法包括：
；
其中：；/>为预设收敛阈值。
为了获得收敛后的非线性响应系数，结合本发明实施例，还存在一种优选的实现方案，具体的，如图2所示，所述收敛后的非线性响应系数的计算方法包括：
通过收敛后的线性响应系数和所述非线性形状参数求出收敛后的非线性响应系数。
具体计算公式为：；
为定标增益。
为了动态调整红外遥感器的定标精度，预设收敛阈值也可以调整，结合本发明实施例，还存在一种优选的实现方案，具体的，所述预设收敛阈值与所述红外遥感器的定标精度负相关。红外遥感器的定标精度越高，预设收敛阈值越小。
综上所述，本实施例一提出一种能够对星载红外遥感器非线性进行精确动态校正的迭代算法，提取非线性形状参数作为表征遥感器非线性特性的物理参数，结合两点定标所得线性响应系数和非线性形状参数计算非线性响应系数，通过多次迭代修正非线性响应的影响后将线性响应系数和非线性响应系数均收敛到准确值，为各类型遥感器在轨定标中的动态非线性校正过程提供新的参考方法。
实施例二：
本实施例二中提供一种校正红外遥感器非线性响应系数的方法的系统，如图3所示，系统包括：
提取模块，所述提取模块用于提取非线性形状参数；
计算模块，所述计算模块用于迭代计算所述红外遥感器的非线性响应系数，得到收敛后的线性响应系数和收敛后的非线性响应系数；
判断模块，所述判断模块用于判断线性响应系数是否满足收敛条件。
实施例三：
本实施例三中提供一种校正星载红外遥感器非线性的迭代的设备，用于执行实施例一中提供的校正红外遥感器非线性响应系数的方法，设备包括：
一个或多个处理器；
存储装置，用于存储一个或多个程序，当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如实施例一所述的校正红外遥感器非线性响应系数的方法。
图4为本发明实施例三提供的校正红外遥感器非线性响应系数的方法的设备结构示意图。图4示出了适于用来实现本发明实施方式的示例性校正红外遥感器非线性响应系数的方法的设备框图。图4显示的校正红外遥感器非线性响应系数的方法的设备仅仅是一个示例，不应对本发明实施例的功能和使用范围带来任何限制。
如图4所示，校正红外遥感器非线性响应系数的方法的设备以通用设备的形式表现。校正红外遥感器非线性响应系数的方法的设备的组件可以包括但不限于：一个或者多个处理器或者处理单元，存储器，连接不同系统组件(包括存储器和处理单元)的总线。
总线表示几类总线结构中的一种或多种，包括存储器总线或者存储器控制器，外围总线，图形加速端口，处理器或者使用多种总线结构中的任意总线结构的局域总线。举例来说，这些体系结构包括但不限于工业标准体系结构(ISA)总线，微通道体系结构(MAC)总线，增强型ISA总线、视频电子标准协会(VESA)局域总线以及外围组件互连(PCI)总线。
校正红外遥感器非线性响应系数的方法的设备典型地包括多种计算机系统可读介质。这些介质可以是任何能够被校正红外遥感器非线性响应系数的方法的设备访问的可用介质，包括易失性和非易失性介质，可移动的和不可移动的介质。
存储器可以包括易失性存储器形式的计算机系统可读介质，例如随机存取存储器(RAM)30和/或高速缓存存储器。校正红外遥感器非线性响应系数的方法的设备可以进一步包括其它可移动/不可移动的、易失性/非易失性计算机系统存储介质。仅作为举例，存储系统可以用于读写不可移动的、非易失性磁介质(图4未显示，通常称为“硬盘驱动器”)。尽管图4中未示出，可以提供用于对可移动非易失性磁盘(例如“软盘”)读写的磁盘驱动器，以及对可移动非易失性光盘(例如CD-ROM，DVD-ROM或者其它光介质)读写的光盘驱动器。在这些情况下，每个驱动器可以通过一个或者多个数据介质接口与总线相连。存储器可以包括至少一个程序产品，该程序产品具有一组(例如至少一个)程序模块，这些程序模块被配置以执行本发明各实施例的功能。
具有一组(至少一个)程序模块的程序/实用工具，可以存储在例如存储器中，这样的程序模块包括但不限于操作系统、一个或者多个应用程序、其它程序模块以及程序数据，这些示例中的每一个或某种组合中可能包括网络环境的实现。程序模块通常执行本发明所描述的实施例中的功能和/或方法。
校正红外遥感器非线性响应系数的方法的设备也可以与一个或多个外部设备(例如键盘、指向设备、显示器等)通信，还可与一个或者多个使得用户能与该校正红外遥感器非线性响应系数的方法的设备交互的设备通信，和/或与使得该校正红外遥感器非线性响应系数的方法的设备能与一个或多个其它设备进行通信的任何设备(例如网卡，调制解调器等等)通信。这种通信可以通过输入/输出(I/O)接口进行。并且，校正红外遥感器非线性响应系数的方法的设备还可以通过网络适配器与一个或者多个网络(例如局域网(LAN)，广域网(WAN)和/或公共网络，例如因特网)通信。如图4所示，网络适配器通过总线与校正红外遥感器非线性响应系数的方法的设备的其它模块通信。应当明白，尽管图中未示出，可以结合校正红外遥感器非线性响应系数的方法的设备使用其它硬件和/或软件模块，包括但不限于：微代码、设备驱动器、冗余处理单元、外部磁盘驱动阵列、RAID系统、磁带驱动器以及数据备份存储系统等。
处理单元通过运行存储在存储器中的程序，从而执行各种功能应用以及数据处理，例如实现本发明任意实施例提供的校正红外遥感器非线性响应系数的方法的方法。也即：构建红外遥感器的红外辐射定标方程组，求解得到线性响应系数和非线性响应系数，并根据线性响应系数和非线性响应系数计算得到非线性形状参数；所述红外遥感器通过预设温度黑体得到初始观测值，利用所述非线性形状参数对所述非线性响应系数红外遥感器的非线性响应系数进行迭代，得到收敛后的线性响应系数和收敛后的非线性响应系数。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (7)
1.一种校正红外遥感器非线性响应系数的方法，其特征在于，所述方法包括：
构建红外遥感器的红外辐射定标方程组，求解得到线性响应系数和非线性响应系数，并根据所述线性响应系数/>和所述非线性响应系数/>计算得到非线性形状参数/>；所述非线性形状参数/>与所述线性响应系数/>和所述非线性响应系数/>的关系为：；
所述红外遥感器通过预设温度黑体得到初始观测值，利用所述非线性形状参数对所述红外遥感器的非线性响应系数进行迭代，得到收敛后的线性响应系数和收敛后的非线性响应系数；所述利用所述非线性形状参数/>对所述红外遥感器的非线性响应系数进行迭代的具体方法包括：
根据所述预设温度黑体得到的信号计数值及所述预设温度黑体的辐射亮度/>计算得到线性响应系数初始值/>非线性响应系数初始值/>；
修正初始辐射亮度中的非线性响应系数，获得修正后的辐射亮度/>；
根据修正后的辐射亮度和所述预设温度黑体得到的信号计数值/>计算出线性响应系数修正值/>；
由所述线性响应系数初始值和所述线性响应系数修正值/>求得线性响应系数更新值/>；
判断修正后的线性响应系数是否符合收敛性；
若满足收敛条件，得到收敛后的线性响应系数和收敛后的非线性响应系数；若不满足收敛条件，将所述线性响应系数更新值赋值给所述线性响应系数初始值/>，使用更新后的线性响应系数初始值参与迭代，继续进行修正，直至收敛。
2.如权利要求1所述的校正红外遥感器非线性响应系数的方法，其特征在于，构建所述红外辐射定标方程组具体包括：
利用实验室预设温度的定标黑体以及所述红外遥感器接收所述定标黑体的有效辐射亮度和所述红外遥感器输出的响应信号计数值/>，构建用于求解线性响应系数/>和非线性响应系数/>的超定方程组：
；
其中：表示定标黑体的预设温度，单位为开尔文（K），/>表示观测/>温度的黑体时，去除冷黑体辐射后红外遥感器接收的有效净辐射亮度，/>表示观测温度的黑体时，去除冷黑体观测计数值后红外遥感器的输出信号计数值，/>为线性响应系数，/>为非线性响应系数。
3.如权利要求2所述的校正红外遥感器非线性响应系数的方法，其特征在于，所述红外辐射定标方程组的求解方法包括：
使用至少3个所述超定方程组通过最小二乘法求解非线性响应系数和线性响应系数；
得到组结果，对非线性响应系数和线性响应系数结果求平均值后，分别作为所述线性响应系数/>和所述非线性响应系数/>。
4.如权利要求1所述的校正红外遥感器非线性响应系数的方法，其特征在于，所述初始观测值包括：所述红外遥感器观测所述预设温度黑体接收得到的信号计数值及所述红外遥感器输出的所述预设温度黑体的辐射亮度/>。
5.如权利要求1所述的校正红外遥感器非线性响应系数的方法，其特征在于，所述判断修正后的线性响应系数是否符合收敛性的方法包括：
；
其中：为线性响应系数更新值；/>为线性响应系数修正值；/>为预设收敛阈值。
6.如权利要求1所述的校正红外遥感器非线性响应系数的方法，其特征在于，所述收敛后的非线性响应系数的计算方法包括：
通过收敛后的线性响应系数和所述非线性形状参数求出收敛后的非线性响应系数。
7.一种校正红外遥感器非线性响应系数的方法的设备，其特征在于，设备包括：
一个或多个处理器；
存储装置，用于存储一个或多个程序，当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如权利要求1-6中任一项所述的校正红外遥感器非线性响应系数的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311264840.4A CN117112966B (zh) | 2023-09-28 | 2023-09-28 | 一种校正红外遥感器非线性响应系数的方法、系统和设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311264840.4A CN117112966B (zh) | 2023-09-28 | 2023-09-28 | 一种校正红外遥感器非线性响应系数的方法、系统和设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117112966A CN117112966A (zh) | 2023-11-24 |
| CN117112966B true CN117112966B (zh) | 2023-12-22 |
Family
ID=88805585
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311264840.4A Active CN117112966B (zh) | 2023-09-28 | 2023-09-28 | 一种校正红外遥感器非线性响应系数的方法、系统和设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117112966B (zh) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103335716A (zh) * | 2013-06-21 | 2013-10-02 | 中国科学院西安光学精密机械研究所 | 基于变积分时间的面阵红外相机定标与非均匀性校正方法 |
| CN110132416A (zh) * | 2019-05-31 | 2019-08-16 | 国家卫星气象中心(国家空间天气监测预警中心) | 宽波段遥感器在轨光谱定标方法及装置 |
| CN110411585A (zh) * | 2019-08-05 | 2019-11-05 | 中国科学院光电技术研究所 | 一种高精度红外辐射测量方法 |
| CN111044153A (zh) * | 2019-12-31 | 2020-04-21 | 华中科技大学 | 一种图谱关联系统红外光谱的非线性定标方法及装置 |
| CN113341434A (zh) * | 2020-02-18 | 2021-09-03 | 国家卫星气象中心(国家空间天气监测预警中心) | 一种红外波段在轨传递定标方法及装置 |
| CN113670191A (zh) * | 2021-08-18 | 2021-11-19 | 国家卫星气象中心(国家空间天气监测预警中心) | 一种红外高光谱干涉仪关键定标参数在轨优化方法 |
| CN115752757A (zh) * | 2022-10-28 | 2023-03-07 | 北京空间机电研究所 | 探测器和内黑体分离的红外扫描相机在轨非均匀校正方法 |
| CN116183036A (zh) * | 2023-02-09 | 2023-05-30 | 中国科学院合肥物质科学研究院 | 一种偏振遥感器短波红外波段背景辐射响应校正方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9143703B2 (en) * | 2011-06-10 | 2015-09-22 | Flir Systems, Inc. | Infrared camera calibration techniques |
-
2023
- 2023-09-28 CN CN202311264840.4A patent/CN117112966B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103335716A (zh) * | 2013-06-21 | 2013-10-02 | 中国科学院西安光学精密机械研究所 | 基于变积分时间的面阵红外相机定标与非均匀性校正方法 |
| CN110132416A (zh) * | 2019-05-31 | 2019-08-16 | 国家卫星气象中心(国家空间天气监测预警中心) | 宽波段遥感器在轨光谱定标方法及装置 |
| CN110411585A (zh) * | 2019-08-05 | 2019-11-05 | 中国科学院光电技术研究所 | 一种高精度红外辐射测量方法 |
| CN111044153A (zh) * | 2019-12-31 | 2020-04-21 | 华中科技大学 | 一种图谱关联系统红外光谱的非线性定标方法及装置 |
| CN113341434A (zh) * | 2020-02-18 | 2021-09-03 | 国家卫星气象中心(国家空间天气监测预警中心) | 一种红外波段在轨传递定标方法及装置 |
| CN113670191A (zh) * | 2021-08-18 | 2021-11-19 | 国家卫星气象中心(国家空间天气监测预警中心) | 一种红外高光谱干涉仪关键定标参数在轨优化方法 |
| CN115752757A (zh) * | 2022-10-28 | 2023-03-07 | 北京空间机电研究所 | 探测器和内黑体分离的红外扫描相机在轨非均匀校正方法 |
| CN116183036A (zh) * | 2023-02-09 | 2023-05-30 | 中国科学院合肥物质科学研究院 | 一种偏振遥感器短波红外波段背景辐射响应校正方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117112966A (zh) | 2023-11-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110793635B (zh) | 热成像测温方法、装置及终端设备 | |
| CN107677375B (zh) | 一种红外辐射测量系统定标装置及定标方法 | |
| Delabie et al. | An accurate and efficient Gaussian fit centroiding algorithm for star trackers | |
| WO2021249437A1 (zh) | 一种热像仪温度修正系统、方法、热像仪及计算机可读取记录介质 | |
| CN111595462B (zh) | 红外成像测温系统标定方法、装置、计算设备及存储介质 | |
| CN112525359B (zh) | 一种基于红外图像的红外测温校正方法及终端 | |
| CN113049109B (zh) | 基于参考黑体的红外测温方法及计算机可读存储介质 | |
| CN109870240A (zh) | 具有环境自适应性的目标红外图像实时定量处理系统 | |
| CN117112966B (zh) | 一种校正红外遥感器非线性响应系数的方法、系统和设备 | |
| CN110726700A (zh) | 烟雾透过率分布测量获取方法和装置 | |
| CN113866768B (zh) | 一种时序干涉雷达相位优化估计方法 | |
| CN114235167B (zh) | 一种温度补偿方法、热成像设备和计算机可读存储介质 | |
| CN111696093B (zh) | 恒星辐射照度相对测量方法和装置 | |
| Liu et al. | An improved Richardson-Lucy algorithm for star image deblurring | |
| CN109932059A (zh) | 一种基于贝叶斯推断的红外热成像温度高精度标定方法 | |
| US20160070025A1 (en) | Weather forecasting apparatus and weather forecasting method | |
| CN110132416B (zh) | 宽波段遥感器在轨光谱定标方法及装置 | |
| KR20210131115A (ko) | 구름의 광학 두께 및 미세물리량 산출 방법 및 장치 | |
| CN114066782B (zh) | 一种气动光学效应校正方法及装置 | |
| Lei et al. | A Maximum Likelihood approach to determine sensor radiometric response coefficients for NPP VIIRS reflective solar bands | |
| Wen et al. | Image sensor correction algorithm for photon transfer curve based on neural network | |
| CN118392324B (zh) | 一种基于背景板扣除内部杂散辐射的红外辐射测量方法 | |
| KR100423061B1 (ko) | 적외선 열상 장비의 결점을 검출하는 방법 | |
| CN115272108A (zh) | 红外成像系统的多点校正方法、装置、设备及存储介质 | |
| CN116429273A (zh) | 基于tof技术的远距离温度补偿方法、装置、终端及介质 |
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 |