CN111257238B - 一种基于星载太阳漫射板的探测元件间相对定标方法 - Google Patents
一种基于星载太阳漫射板的探测元件间相对定标方法 Download PDFInfo
- Publication number
- CN111257238B CN111257238B CN202010080321.2A CN202010080321A CN111257238B CN 111257238 B CN111257238 B CN 111257238B CN 202010080321 A CN202010080321 A CN 202010080321A CN 111257238 B CN111257238 B CN 111257238B
- Authority
- CN
- China
- Prior art keywords
- diffusion plate
- relative
- solar
- probe
- satellite
- 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
- 238000009792 diffusion process Methods 0.000 title claims abstract description 114
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000001514 detection method Methods 0.000 title description 6
- 239000000523 sample Substances 0.000 claims abstract description 124
- 238000002310 reflectometry Methods 0.000 claims abstract description 59
- 238000012937 correction Methods 0.000 claims abstract description 36
- 230000005855 radiation Effects 0.000 claims abstract description 21
- 230000004044 response Effects 0.000 claims abstract description 20
- 210000001747 pupil Anatomy 0.000 claims description 14
- 230000010287 polarization Effects 0.000 abstract description 13
- 239000000443 aerosol Substances 0.000 abstract description 4
- 238000004364 calculation method Methods 0.000 description 10
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 6
- 230000008859 change Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 229910002092 carbon dioxide Inorganic materials 0.000 description 3
- 239000001569 carbon dioxide Substances 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000009533 lab test Methods 0.000 description 2
- 238000011545 laboratory measurement Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 239000005427 atmospheric aerosol Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000013028 emission testing Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/21—Polarisation-affecting properties
-
- 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
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/447—Polarisation spectrometry
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
本发明公开了一种基于星载太阳漫射板的探测元件间相对定标方法，该方法包括以下步骤：计算遥感器探元相对定标系数初始值；确定遥感器探元间的相对辐射响应初始状态；计算漫射板不同部位反射率的分布订正系数；分析并评估太阳漫射板不同部位反射率的分布情况；计算遥感器在轨运行期间的探元间相对定标系数；更新探元间相对定标系数。通过该方法，解决了线阵推扫式偏振遥感器在轨工作寿命期内，探元间相对辐射响应随时间发生变化问题，解决了线阵推扫式遥感成像仪图像纵向条纹现象不断恶化的订正问题；同时针对线阵推扫式偏振遥感器的探元间相对定标问题，提供了基于在轨太阳漫射板的相对定标方法，成功应用于云与气溶胶偏振成像仪。
Description
技术领域
本发明涉及探测元件定标技术领域，具体来说，涉及一种基于星载太阳漫射板的探测元件间相对定标方法。
背景技术
现有技术方案在处理遥感器多探元的辐射响应不一致问题时，一般采用发射前测试法和对地观测数据统计法。
发射前测试法是基于卫星发射前，遥感器在实验室的测试结果，用于遥感器在轨工作期间的探元间相对定标。遥感器在实验室测试时，使用积分球作为理想辐射定标源，积分球出射光的空间均匀性、角度均匀性和时间稳定性方面都很理想。遥感器对积分球出口观测的计数值，直接可以获取探元间的相对辐射定标系数。存在的问题是遥感器发射后，在轨工作环境与实验室测试环境存在较大差异，遥感器长期在轨工作后，各个探元的辐射响应会发生不同的变化。实验室测试得到的相对辐射定标系数不适用于遥感器全生命期内观测数据的订正，需要在轨更新相对辐射定标系数。
另一种相对定标方法是基于对地观测数据统计特征而计算相对定标系数。统计各个探元的对地观测数据，提取统计量(如平均值，输出码值的累积概率等)，在不同探元间建立联系，计算出相对定标系数。统计法假设了各个探元所观测到的数据在从统计学角度是一致的，但是对于线阵推扫式成像仪，地面图像左侧和右侧的像素在光照上存在系统的差异，因此统计法在针对线阵推扫式成像仪的相对定标问题时，其前提假设不成立，因此需要做出另外的改进。
发明内容
针对相关技术中的上述技术问题，本发明提出一种基于星载太阳漫射板的探测元件间相对定标方法，能够克服现有技术的上述不足。
为实现上述技术目的，本发明的技术方案是这样实现的：
一种基于星载太阳漫射板的探测元件间相对定标方法，该方法包括以下步骤：
S1：基于发射前实验室积分球观测数据，计算遥感器探元相对定标系数初始值；
S2：确定遥感器探元间的相对辐射响应初始状态；
S3：基于卫星在轨后的首次对太阳漫射板观测数据，计算漫射板不同部位反射率的订正系数；
S4：分析并评估太阳漫射板不同部位反射率的分布情况；
S5：基于遥感器在轨工作期间对太阳漫射板的观测数据，计算遥感器在轨运行期间的探元间相对定标系数；
S6：更新探元间相对定标系数。
进一步的，所述步骤S1包括以下步骤：
S11：遥感器输出数据，其中，数据表示为
DN(t,i)＝A·RPrelaunch(i)·E(t,i)+n(i,t)，
式中，DN(t,i)表示t时刻第i个探元的输出计数值，A为绝对辐射定标系数，R(i)为第i个探元的相对定标系数，E(t,i)为t时刻第i个探元入瞳处能量，n(t,i)表示t时刻第i个探元的随机误差；
S12：获取并分析发射前实验室积分球数据；
S13：计算探元间相对定标系数；
S14：生成数据初始值。
进一步的，所述步骤S3包括以下步骤：
S31：在轨后卫星首次采集太阳漫射板数据；
S32：分析并输出数据，其中，输出数据表示为
DN(t,i)＝A·R(i)·ρ·k(i)Es(t)+n(i)，
E(t,i)＝ρ·k(i)·Es(t)，
式中t为帧时间，对应到图像上即为行号，i为第i个探元的编号，入瞳处能量可以表示为太阳漫射板的反射率与照射到漫射板上的太阳能量的乘积，ρ为标准探元对应太阳漫射板部分的反射率，k(i)为反射率空间订正系数随着i连续缓慢变化，Es随着时间t连续缓慢变化；
S33：计算分布订正系数。
进一步的，所述步骤S5包括以下步骤：
S51：获取在轨工作期间太阳漫射板的数据；
S52：分析太阳漫射板的数据；
S53：计算运行期间的探元间相对定标系数。
进一步的，所述步骤S11中，所有探元的随机误差同分布于均值为0，标准差为σ的正态分布；所有E(t,i)数值相同；标准探元的相对定标系数等于1。
进一步的，所述步骤S3和步骤S5中，标准探元对应的太阳漫射板反射率订正系数等于1。
本发明的有益效果：通过该方法，解决了线阵推扫式偏振遥感器在轨工作寿命期内，探元间相对辐射响应随时间发生变化问题，解决了线阵推扫式遥感成像仪图像纵向条纹现象不断恶化的订正问题；同时针对线阵推扫式偏振遥感器的探元间相对定标问题，提供了基于在轨太阳漫射板的相对定标方法，成功应用于云与气溶胶偏振成像仪(CAPI)。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是根据本发明实施例所述的一种基于星载太阳漫射板的探测元件间相对定标方法的流程框图；
图2是根据本发明实施例所述的一种基于星载太阳漫射板的探测元件间相对定标方法的积分球原始观测数据图；
图3是根据本发明实施例所述的一种基于星载太阳漫射板的探测元件间相对定标方法的初始化的相对定标系数图；
图4是根据本发明实施例所述的一种基于星载太阳漫射板的探测元件间相对定标方法的首次太阳漫射板观测数据图；
图5是根据本发明实施例所述的一种基于星载太阳漫射板的探测元件间相对定标方法的全部探元单帧漫射板观测计数值图；
图6是根据本发明实施例所述的一种基于星载太阳漫射板的探测元件间相对定标方法的漫射板反射率订正系数图；
图7是根据本发明实施例所述的一种基于星载太阳漫射板的探测元件间相对定标方法的在轨工作后太阳漫射板观测数据图；
图8是根据本发明实施例所述的一种基于星载太阳漫射板的探测元件间相对定标方法的全部探元单帧漫射板观测计数值图；
图9是根据本发明实施例所述的一种基于星载太阳漫射板的探测元件间相对定标方法的相对定标系数更新值图；
图10是根据本发明实施例所述的一种基于星载太阳漫射板的探测元件间相对定标方法的在轨相对定标系数与发射前测量值的绝对差异图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，根据本发明实施例所述的种基于星载太阳漫射板的探测元件间相对定标方法，包括以下步骤：
S1：基于发射前实验室积分球观测数据，计算遥感器探元相对定标系数初始值；
S2：确定遥感器探元间的相对辐射响应初始状态；
S3：基于卫星在轨后的首次对太阳漫射板观测数据，计算漫射板不同部位反射率的订正系数；
S4：分析并评估太阳漫射板不同部位反射率的分布情况；
S5：基于遥感器在轨工作期间对太阳漫射板的观测数据，计算遥感器在轨运行期间的探元间相对定标系数；
S6：更新探元间相对定标系数。
步骤S1包括以下步骤：
S11：遥感器输出数据，其中，数据表示为
DN(t,i)＝A·RPrelaunch(i)·E(t,i)+n(i,t)，
式中，DN(t,i)表示t时刻第i个探元的输出计数值，A为绝对辐射定标系数，R(i)为第i个探元的相对定标系数，E(t,i)为t时刻第i个探元入瞳处能量，n(t,i)表示t时刻第i个探元的随机误差；
S12：获取并分析发射前实验室积分球数据；
S13：计算探元间相对定标系数；
S14：生成数据初始值。
步骤S3包括以下步骤：
S31：在轨后卫星首次采集太阳漫射板数据；
S32：分析并输出数据，其中，输出数据表示为
DN(t,i)＝A·R(i)·ρ·k(i)Es(t)+n(i)，
E(t,i)＝ρ·k(i)·Es(t)，
式中t为帧时间，对应到图像上即为行号，i为第i个探元的编号，入瞳处能量可以表示为太阳漫射板的反射率与照射到漫射板上的太阳能量的乘积，ρ为标准探元对应太阳漫射板部分的反射率，k(i)为反射率空间订正系数随着i连续缓慢变化，Es随着时间t连续缓慢变化；
S33：计算分布订正系数。
步骤S5包括以下步骤：
S51：获取在轨工作期间太阳漫射板的数据；
S52：分析太阳漫射板的数据；
S53：计算运行期间的探元间相对定标系数。
在本发明的一个具体实施例中，所述步骤S11中，所有探元的随机误差同分布于均值为0，标准差为σ的正态分布；所有E(t,i)数值相同；标准探元的相对定标系数等于1。
在本发明的一个具体实施例中，所述步骤S3和步骤S5中，标准探元对应的太阳漫射板反射率订正系数等于1。
为了方便理解本发明的上述技术方案，以下通过具体使用方式上对本发明的上述技术方案进行详细说明。
本发明基于遥感器对太阳漫射板的观测数据，实现寿命期内探元间辐射响应的相对定标。
本发明主要流程为遥感器探元间相对定标系数的初始值计算，太阳漫射板反射率空间分布模型计算，相对定标系数在轨更新计算。
遥感器输出图像可用如下模型表示：
DN＝coef·E+n
式中DN为输出计数值，E为入瞳处能量，coef为探元的辐射响应系数，其描述了遥感器探元对入射能量的辐射响应情况。n为随机噪声，一般为均值为0，标准差为σ的正态分布。
对于线阵推扫式遥感成像仪，其单帧输出数据可表示为：
DN(i)＝A·R(i)·E(i)+n(i)
coef(i)＝A·R(i)
R(standard)＝1
式中i为线阵探测器的第i个探元的编号，coef可以表示为绝对定标系数A和相对定标系数R的乘积，绝对定标系数A描述了标准探元对入射能量的输出响应情况，相对定标系数R描述了线阵探元相对于标准探元的相对辐射响应情况，由定义可知标准探元的R(i)等于1。
1.遥感器探元间相对定标系数的初始值计算模块
用于确定遥感器探元间的相对辐射响应初始状态。基于发射前实验室积分球观测数据，计算探元间的相对定标系数。观测积分球时，遥感器输出数据可以表示为：
DN(t,i)＝A·RPrelaunch(i)·E(t,i)+n(i,t)
式中，DN(t,i)表示t时刻第i个探元的输出计数值，A为绝对辐射定标系数，R(i)为第i个探元的相对定标系数，E(t,i)为t时刻第i个探元入瞳处能量，n(t,i)表示t时刻第i个探元的随机误差，所有探元的随机误差同分布于均值为0，标准差为σ的正态分布。考虑到积分球的时间稳定性和空间均匀性，所有E(t,i)数值相同。另外，标准探元的相对定标系数等于1，因此相对定标系数的初始值的最优估计，可以利用下式计算：
考虑到约束条件Rprelaunch(standard)＝1，则相对定标系数初始值计算式如下：
2.太阳漫射板反射率空间分布订正系数计算模块
用于评估太阳漫射板不同部位反射率的分布情况。基于卫星在轨后的首次对太阳漫射板观测数据，计算漫射板不同部位反射率的订正系数。标准探元对应的太阳漫射板反射率订正系数等于1。在轨对太阳漫射板观测时，其输出数据可以表示为：
DN(t,i)＝A·Rprelaunch(i)·ρ·k(i)Es(t)+n(i)
E(t,i)＝ρ·k(i)·Es(t)
式中t为帧时间，对应到图像上即为行号，i为第i个探元的编号，入瞳处能量可以表示为太阳漫射板的反射率与照射到漫射板上的太阳能量的乘积。推扫式成像仪观测漫射板时，不同探元对应于漫射板不同的部位，因此漫射板反射率的空间分布应该考虑到模型中，ρ为标准探元对应太阳漫射板部分的反射率，k(i)为反射率空间订正系数，随着i连续缓慢变化。由于数据为卫星入轨工作初期获取，这里认为遥感器的探元间相对定标系数与发射前实验室测量结果一致，没有产生变化。由于太阳照到漫射板的角度随着时间变化，因此Es随着时间t连续缓慢变化。另外太阳照射到漫射板的能量是均匀地，因此漫射板不同部位受到的太阳照射能量是一致的。
太阳漫射板反射率的空间分布计算可从下式：
考虑到约束条件k(standard)＝1，则太阳漫射板反射率空间订正系数的最优估计如下式：
3.相对定标系数在轨更新计算模块
用于计算在轨工作后，探元间相对定标系数。基于遥感器在轨工作期间对太阳漫射板的观测数据，漫射板反射率分布订正系数不随时间变化，计算遥感器在轨运行期间的探元间相对定标系数。标准探元对应的太阳漫射板反射率订正系数等于1。在轨对太阳漫射板观测时，其输出数据可以表示为：
DN(t,i)＝A·R(i)·ρ·k(i)Es(t)+n(i)
E(t,i)＝ρ·k(i)·Es(t)
式中t为帧时间，对应到图像上即为行号，i为第i个探元的编号，入瞳处能量可以表示为太阳漫射板的反射率与照射到漫射板上的太阳能量的乘积。ρ为标准探元对应太阳漫射板部分的反射率，k(i)为反射率空间订正系数，随着i连续缓慢变化。漫射板的反射率一般不随时间变化。由于太阳照到漫射板的角度随着时间变化，因此Es随着时间t连续缓慢变化。另外太阳照射到漫射板的能量是均匀地，因此漫射板不同部位受到的太阳照射能量是一致的。
遥感器探元间相对定标系数计算如下式：
考虑到约束条件k(standard)＝1，则相对定标系数的在轨更新值的最优估计如下式：
在一个实施例中，云与气溶胶偏振成像仪(Cloud and Aerosol PolarizationImager，CAPI)是全球二氧化碳监测科学试验卫星(TanSat)的载荷之一。其主要功能是为卫星主载荷高光谱二氧化碳探测仪反演二氧化碳柱总量(XCO2)提供大气和地表信息。遥感数据可以用于大气气溶胶、云量和地表反射率等参数反演。
CAPI是线阵推扫式偏振成像仪器，获取紫外到近红外波段的太阳反射能量，在670nm和1640nm波段设置了偏振探测通道，实现0°、60°、120°三个方向的偏振测量功能。在380nm、870nm和1375nm等波段设置非偏振通道。CAPI幅宽为375km，可以得到250m和1000m的两种空间分辨率对地观测数据。
表1通道技术指标：
以870nm波段，即通道5的数据为实施例详述本发明方案。
1.遥感器探元间相对定标系数的初始值计算模块
用于确定遥感器探元间的相对辐射响应初始状态。基于发射前实验室积分球观测数据，计算探元间的相对定标系数。观测积分球时，遥感器输出数据可以表示为：
DN(t,i)＝A·RPrelaunch(i)·E(t,i)+n(i,t)
式中，DN(t,i)表示t时刻第i个探元的输出计数值，A为绝对辐射定标系数，R(i)为第i个探元的相对定标系数，E(t,i)为t时刻第i个探元入瞳处能量，n(t,i)表示t时刻第i个探元的随机误差，所有探元的随机误差同分布于均值为0，标准差为σ的正态分布。考虑到积分球的时间稳定性和空间均匀性，所有E(t,i)数值相同。另外，标准探元的相对定标系数等于1，因此相对定标系数的初始值的最优估计，可以利用下式计算：
考虑到约束条件Rprelaunch(standard)＝1，则相对定标系数初始值计算式如下：
实验室测试数据如图2所示。图中横轴为探元编号，纵轴为积分球观测输出码值，以第16、17和18帧数据为例作图。图中可知，积分球输出的能量是均匀的，但各个探元的输出码值存在差异，这种差异就是探元间辐射响应差异引起的。以探元间相对定标系数描述这种探元间相对辐射响应差异。按照相对定标系数初始值计算式，计算得出的发射前相对定标系数如图3所示。
2.太阳漫射板反射率空间分布订正系数计算模块
用于评估太阳漫射板不同部位反射率的分布情况。基于卫星在轨后的首次对太阳漫射板观测数据，计算漫射板不同部位反射率的订正系数。
标准探元对应的太阳漫射板反射率订正系数等于1。在轨对太阳漫射板观测时，其输出数据可以表示为：
DN(t,i)＝A·Rprelaunch(i)·ρ·k(i)Es(t)+n(i)
E(t,i)＝ρ·k(i)·Es(t)
式中t为帧时间，对应到图像上即为行号，i为第i个探元的编号，入瞳处能量可以表示为太阳漫射板的反射率与照射到漫射板上的太阳能量的乘积。推扫式成像仪观测漫射板时，不同探元对应于漫射板不同的部位，因此漫射板反射率的空间分布应该考虑到模型中，ρ为标准探元对应太阳漫射板部分的反射率，k(i)为反射率空间订正系数，随着i连续缓慢变化。由于数据为卫星入轨工作初期获取，这里认为遥感器的探元间相对定标系数与发射前实验室测量结果一致，没有产生变化。由于太阳照到漫射板的角度随着时间变化，因此Es随着时间t连续缓慢变化。另外太阳照射到漫射板的能量是均匀的，因此漫射板不同部位受到的太阳照射能量是一致的。
太阳漫射板反射率的空间分布计算可从下式：
考虑到约束条件k(standard)＝1，则反射率空间订正系数初始值的最优估计如下式：
CAPI的首次在轨太阳漫射板观测于2017年5月24日完成。太阳漫射板观测数据的分布存在以下特点。以通道5的标准探元(第800号探元)所观测到的数据为例。随着卫星的姿态机动，太阳照射太阳漫射板的角度发射变化，导致进入遥感器的能量随着时间发生变化，如图4所示。
另外，在同一帧内，各个探元对太阳漫射板的观测数据分布如图5所示。输出码值存在差异，这一方面是探元间的相对辐射响应差异的影响，另外，太阳漫射不同部位的反射率也存在差异。利用首次对太阳漫射板的观测数据，提取出不同部位反射率的订正系数。
按照太阳漫射板反射率空间订正系数的计算式处理以上数据，结果如图6所示。
3.相对定标系数在轨更新计算模块
用于计算在轨工作后，探元间相对定标系数。基于遥感器在轨工作期间对太阳漫射板的观测数据，漫射板反射率分布订正系数不随时间变化，计算遥感器在轨运行期间的探元间相对定标系数。标准探元对应的太阳漫射板反射率订正系数等于1。在轨对太阳漫射板观测时，其输出数据可以表示为：
DN(t,i)＝A·R(i)·ρ·k(i)Es(t)+n(i)
E(t,i)＝ρ·k(i)·Es(t)
式中t为帧时间，对应到图像上即为行号，i为第i个探元的编号，入瞳处能量可以表示为太阳漫射板的反射率与照射到漫射板上的太阳能量的乘积。ρ为标准探元对应太阳漫射板部分的反射率，k(i)为反射率空间订正系数，随着i连续缓慢变化。漫射板的反射率一般不随时间变化。由于太阳照到漫射板的角度随着时间变化，因此Es随着时间t连续缓慢变化。另外太阳照射到漫射板的能量是均匀的，因此漫射板不同部位受到的太阳照射能量是一致的。
遥感器探元间相对定标系数计算如下式：
考虑到约束条件k(standard)＝1，则相对定标系数的在轨更新值的最优估计如下式：
CAPI于2017年9月7日再次实现太阳漫射板观测。与首次太阳漫射板观测数据相比，此次数据分布相似。以通道5的标准探元(第800号探元)所观测到的数据为例。随着卫星的姿态机动，太阳照射太阳漫射板的角度发射变化，导致进入遥感器的能量随着时间发生变化，如图7所示。
另外，在同一帧内，各个探元对太阳漫射板的观测数据分布如8图所示。输出码值存在差异，这一方面是探元间的相对辐射响应差异的影响，另外，太阳漫射不同部位的反射率也存在差异。利用首次对太阳漫射板的观测数据，提取出不同部位反射率的订正系数。
按照太阳漫射板反射率空间订正系数的计算式处理以上数据，结果如图9所示。
与发射前相对定标系数的相对差异如图10所示。相对定标系数变化集中于-0.2到0.4之间。
综上所述，借助于本发明的上述技术方案，通过该方法，解决了线阵推扫式偏振遥感器在轨工作寿命期内，探元间相对辐射响应随时间发生变化问题，解决了线阵推扫式遥感成像仪图像纵向条纹现象不断恶化的订正问题；同时针对线阵推扫式偏振遥感器的探元间相对定标问题，提供了基于在轨太阳漫射板的相对定标方法，成功应用于云与气溶胶偏振成像仪(CAPI)。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (6)
1.一种基于星载太阳漫射板的探测元件间相对定标方法，其特征在于，包括以下步骤：
S1：基于发射前实验室积分球观测数据，计算遥感器探元相对定标系数初始值；
S2：确定遥感器探元间的相对辐射响应初始状态；
S3：基于卫星在轨后的首次对太阳漫射板观测数据，计算漫射板不同部位反射率的反射率分布订正系数k(i)，公式为：
S4：分析并评估太阳漫射板不同部位反射率的分布情况；
S5：基于遥感器在轨工作期间对太阳漫射板的观测数据，计算遥感器在轨运行期间的探元间相对定标系数；
S6：更新探元间相对定标系数。
2.根据权利要求1所述的一种基于星载太阳漫射板的探测元件间相对定标方法，其特征在于，所述步骤S1包括以下步骤：
S11：遥感器输出数据，其中，数据表示为
DN(t,i)＝A·RPrelaunch(i)·E(t,i)+n(i,t)，
式中，DN(t，i)表示t时刻第i个探元的输出计数值，A为绝对辐射定标系数，R(i)为第i个探元的相对定标系数，下脚标Prelaunch表示的是卫星发射前地面实验室测量值，用以区分在轨计算值，E(t，i)为t时刻第i个探元入瞳处能量，n(i，t)表示t时刻第i个探元的随机误差；
S12：获取并分析发射前实验室积分球数据；
S13：计算探元间相对定标系数；
S14：生成数据初始值。
3.根据权利要求1所述的一种基于星载太阳漫射板的探测元件间相对定标方法，其特征在于，所述步骤S3包括以下步骤：
S31：在轨后卫星首次采集太阳漫射板数据；
S32：分析并输出数据，其中，输出数据表示为
DN(t，i)＝A·R(i)·ρ·k(i)Es(t)+n(i)，
E(t，i)＝ρ·k(i)·Es(t)，
式中t为帧时间，对应到图像上即为行号，DN(t，i)表示t时刻第i个探元的输出计数值，A为绝对辐射定标系数，i为第i个探元的编号，R(i)为第i个探元的相对定标系数，入瞳处能量表示为太阳漫射板的反射率与照射到漫射板上的太阳能量的乘积，ρ为标准探元对应太阳漫射板部分的反射率，k(i)为反射率分布订正系数随着i连续缓慢变化，Es(t)为t时刻的太阳常数，Es随着时间t连续缓慢变化，n(i)表示第i个探元的随机误差，E(t，i)为t时刻第i个探元入瞳处能量；
S33：计算反射率分布订正系数。
4.根据权利要求1所述的一种基于星载太阳漫射板的探测元件间相对定标方法，其特征在于，所述步骤S5包括以下步骤：
S51：获取在轨工作期间太阳漫射板的数据；
S52：分析太阳漫射板的数据；
S53：计算运行期间的探元间相对定标系数。
5.根据权利要求2所述的一种基于星载太阳漫射板的探测元件间相对定标方法，其特征在于，所述步骤S11中，所有探元的随机误差同分布于均值为0，标准差为σ的正态分布；所有E(t，i)数值相同；标准探元的相对定标系数等于1。
6.根据权利要求1所述的一种基于星载太阳漫射板的探测元件间相对定标方法，其特征在于，所述步骤S3和步骤S5中，标准探元对应的太阳漫射板反射率分布订正系数等于1。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010080321.2A CN111257238B (zh) | 2020-02-05 | 2020-02-05 | 一种基于星载太阳漫射板的探测元件间相对定标方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010080321.2A CN111257238B (zh) | 2020-02-05 | 2020-02-05 | 一种基于星载太阳漫射板的探测元件间相对定标方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111257238A CN111257238A (zh) | 2020-06-09 |
| CN111257238B true CN111257238B (zh) | 2023-03-14 |
Family
ID=70945644
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010080321.2A Active CN111257238B (zh) | 2020-02-05 | 2020-02-05 | 一种基于星载太阳漫射板的探测元件间相对定标方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111257238B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112070414B (zh) * | 2020-09-17 | 2023-11-07 | 国核电力规划设计研究院重庆有限公司 | 一种基于多站点观测数据的太阳能资源评估方法、系统及介质 |
| CN112161943B (zh) * | 2020-09-29 | 2021-04-16 | 中国科学院地理科学与资源研究所 | TanSat卫星XCO2反演数据的偏差校正方法和系统 |
| CN113534158B (zh) * | 2021-06-10 | 2024-06-11 | 国家卫星气象中心(国家空间天气监测预警中心) | 基于深度学习的静止轨道气象卫星雷达反射率因子反演方法 |
| CN114218786B (zh) * | 2021-12-13 | 2024-07-05 | 北京理工大学 | 非偏振卫星传感器在轨偏振辐射特性反演方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102324098A (zh) * | 2011-08-23 | 2012-01-18 | 中国资源卫星应用中心 | 一种结合实验室定标和均匀景统计的相对辐射定标方法 |
| CN104065892A (zh) * | 2014-06-24 | 2014-09-24 | 中国资源卫星应用中心 | 一种凝视卫星面阵ccd相机的相对辐射校正方法 |
| CN106485683A (zh) * | 2016-10-20 | 2017-03-08 | 中国科学院上海技术物理研究所启东光电遥感中心 | 一种基于场景的图像自适应非均匀校正方法 |
| CN106600646A (zh) * | 2016-11-25 | 2017-04-26 | 北京空间机电研究所 | 一种红外摆扫相机在轨图像均匀性校正方法 |
| CN107036629A (zh) * | 2017-04-20 | 2017-08-11 | 武汉大学 | 视频卫星在轨相对辐射定标方法及系统 |
| CN107633487A (zh) * | 2017-09-01 | 2018-01-26 | 天津津航技术物理研究所 | 一种航空摆扫式多光谱扫描仪图像的系统级相对辐射校正方法 |
| CN109712089A (zh) * | 2018-12-14 | 2019-05-03 | 航天恒星科技有限公司 | 适用于出口型遥感卫星红外短波载荷相对辐射校正的方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7688438B2 (en) * | 2005-12-20 | 2010-03-30 | Raytheon Company | Scanning solar diffuser relative reflectance monitor |
| CN110120018A (zh) * | 2019-04-10 | 2019-08-13 | 武汉大学 | 一种面阵高动态范围夜光成像卫星在轨相对辐射定标方法 |
-
2020
- 2020-02-05 CN CN202010080321.2A patent/CN111257238B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102324098A (zh) * | 2011-08-23 | 2012-01-18 | 中国资源卫星应用中心 | 一种结合实验室定标和均匀景统计的相对辐射定标方法 |
| CN104065892A (zh) * | 2014-06-24 | 2014-09-24 | 中国资源卫星应用中心 | 一种凝视卫星面阵ccd相机的相对辐射校正方法 |
| CN106485683A (zh) * | 2016-10-20 | 2017-03-08 | 中国科学院上海技术物理研究所启东光电遥感中心 | 一种基于场景的图像自适应非均匀校正方法 |
| CN106600646A (zh) * | 2016-11-25 | 2017-04-26 | 北京空间机电研究所 | 一种红外摆扫相机在轨图像均匀性校正方法 |
| CN107036629A (zh) * | 2017-04-20 | 2017-08-11 | 武汉大学 | 视频卫星在轨相对辐射定标方法及系统 |
| CN107633487A (zh) * | 2017-09-01 | 2018-01-26 | 天津津航技术物理研究所 | 一种航空摆扫式多光谱扫描仪图像的系统级相对辐射校正方法 |
| CN109712089A (zh) * | 2018-12-14 | 2019-05-03 | 航天恒星科技有限公司 | 适用于出口型遥感卫星红外短波载荷相对辐射校正的方法 |
Non-Patent Citations (4)
| Title |
|---|
| Prelaunch radiometric calibration of the TanSat atmospheric carbon dioxide grating spectrometer;Zhongdong Yang,Yuquan Zhen;《IEEE TRANSACTION ON GEOSCIENCE AND REMOTE SENSING》;20180731;第56卷(第7期);第4225-4233页 * |
| 基于太阳漫反射板线阵CCD相对辐射定标方法研究;许和鱼，张黎明，李鑫;《光学学报》;20191211;第40卷(第6期);第179-187页 * |
| 基于太阳漫射板的星上定标方法研究;黄文薪，张黎明;《光谱学与光谱分析》;20170325;第37卷(第3期);第952-957页 * |
| 基于月球反射的遥感器定标跟踪监测;吴荣华，张鹏，杨忠东;《遥感学报》;20160325;第20卷(第2期);第278-289页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111257238A (zh) | 2020-06-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111257238B (zh) | 一种基于星载太阳漫射板的探测元件间相对定标方法 | |
| Liu et al. | Keck Laser Guide Star Adaptive Optics Monitoring of 2MASS J15344984–2952274AB: First Dynamical Mass Determination of a Binary T Dwarf | |
| Tubiana et al. | Scientific assessment of the quality of OSIRIS images | |
| Xiong et al. | Terra MODIS on-orbit spatial characterization and performance | |
| Yoshida et al. | Improvement of the retrieval algorithm for GOSAT SWIR XCO 2 and XCH 4 and their validation using TCCON data | |
| Hagolle et al. | Results of POLDER in-flight calibration | |
| US8067738B1 (en) | Space based calibration transfer spectroradiometer | |
| CN105092055B (zh) | 基于冷云目标的气象卫星太阳反射波段辐射定标方法 | |
| Lenhard et al. | Independent laboratory characterization of NEO HySpex imaging spectrometers VNIR-1600 and SWIR-320m-e | |
| Kern et al. | Intercomparison of SO2 camera systems for imaging volcanic gas plumes | |
| US11711429B2 (en) | Method for calibrating networks of environmental sensors | |
| CN112798013B (zh) | 一种对光学载荷在轨绝对辐射定标结果进行验证的方法 | |
| CN108226059B (zh) | 一种卫星高光谱co2探测仪在轨辐射定标方法 | |
| US6294785B1 (en) | Infrared sensor response calibration using atmospheric limb emission measurements | |
| CN116735008B (zh) | 红外交叉辐射的定标方法、装置、电子设备及介质 | |
| CN112665829B (zh) | 一种光学遥感卫星的波段间定标方法 | |
| CN117848977A (zh) | 一种星载大气主要温室气体探测仪在轨星上定标装置与方法 | |
| Pan et al. | Solar observation of Ozone Mapping and Profiler Suite nadir system during the first 3 years of on-orbit operation | |
| Tukiainen et al. | Retrieval of ozone profiles from GOMOS limb scattered measurements | |
| CN109655415A (zh) | 一种波长偏移校正方法及装置及计算机设备 | |
| CN111256818B (zh) | 一种线阵推扫式遥感成像仪月球辐照度观测值的计算方法 | |
| Lichtenberg et al. | SCIAMACHY Level1 data: Calibration concept and in-flight calibration | |
| Xiong et al. | On-orbit spatial characterizations for Terra MODIS | |
| US6987566B1 (en) | Methods and apparatus for analyzing mirror reflectance | |
| CN118129905B (zh) | 一种星载高光谱成像仪曲线型辐射定标方法及系统 |
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 |