一种高频次重复凝视成像空基遥感载荷在轨测试方法与流程

1.本发明涉及定量综合地评价遥感影像质量工程技术领域，尤其涉及一种高频次重复凝视成像空基遥感载荷在轨测试方法。


背景技术：

2.gf-4卫星是我国
″
高分专项
″
的重要组成部分，旨在获取极高的时间分辨率、和中等空间分辨率全国范围及周边地区的多光谱、中红外图像信息，能够满足减灾、林业、气象等用户广泛需求，具有极高的社会经济效益，并且具有技术引领性和先进性。该星配置地球同步轨道50米分辨率的凝视相机组(可见光波段与红外波段)，成像区域7000km
×
7000km，单景覆盖区域400km
×
400km，可实现分钟级高时间分辨率的遥感监测任务，设计寿命8年。
3.本发明为gf-4卫星图像质量评价技术的研究。卫星遥感影像质量优劣直接关系到遥感影像信息获取的可靠性和准确性，是遥感成像系统设计、遥感信息处理与应用的重要基础。由于光学遥感成像过程受到大气、辐照、温度、传感器以及地物特性等多种因素影响，定量综合地评价遥感影像质量是一项复杂的系统工程。就遥感影像本身而言，光学遥感影像质量总体表现为辐射质量和几何质量，其中，辐射质量主要包括系统调制传递函数(modulation transfer function，mtf)、信噪比(signal noise ratio，snr)，也是光学遥感成像质量评价的重要依据。根据总体方案安排，本研究基于gf-4卫星成像的特点，开展在轨调制传递函数(mtf)和信噪比(snr)测试工作。


技术实现要素：

4.为解决上述技术问题，本发明的目的是提供一种高频次重复凝视成像空基遥感载荷在轨测试方法，该方法可在数分钟之内完成任务机动响应，并对用户指定的目标区域成像，就遥感影像本身而言，光学遥感影像质量总体表现为辐射质量和几何质量，其中，辐射质量主要包括系统调制传递函数(modulation transfer function，mtf)、信噪比(signal noise ratio，snr)，也是光学遥感成像质量评价的重要依据，本研究为mtf、sar的获取提供了一种更高效、更准确的方法。
5.本发明的目的通过以下的技术方案来实现：
6.一种高频次重复凝视成像空基遥感载荷在轨测试方法，包括：
7.a、选择满足mtf和snr测试需求的场地，布设辐射特性靶标；
8.b、监测太阳光谱信息和地球辐射收支信息，使用光谱仪对靶标及自然地物反射率进行测量；
9.c、基于同步大气气溶胶观测数据，对gf-4连续观测的影像数据进行筛选，剔除大气扰动时刻获取的卫星影像数据；
10.d、通过对线扩散函数lsf的计算与处理得到得到截止频率内不同频率的mtf值；
11.e、对均匀场景地物的连续多帧影像数据在时间维度上进行列向量统计计算，获取均方根误差，基于配准后的多帧影像均匀地物计算snr。
12.与现有技术相比，本发明的一个或多个实施例可以具有如下优点：
13.可以定量综合、高效地评价遥感影像质量，可在数分钟之内完成任务机动响应，并对用户指定的目标区域成像，包括对目标区域进行高频重复凝视观测、快速对某一区域进行拼接成像、对多个区域进行巡查成像。
附图说明
14.图1是高频次重复凝视成像空基遥感载荷在轨测试方法流程图；
15.图2是研究区概况图；
16.图3是基于gf-4高频重复凝视观测数据的脉冲法mtf测试流程图；
17.图4是辐射靶标布设示意图；
18.图5是高分四号卫星高频重复凝视观测(微错位)数据配准示意图；
19.图6是基于配准后的多帧gf-4卫星数据恢复lsf示意图；
20.图7是基于gf-4高频重复凝视观测数据的snr测试的流程图。
具体实施方式
21.为使本发明的目的、技术方案和优点更加清楚，下面将结合实施例及附图对本发明作进一步详细的描述。
22.本次试验主要区域在农场，地理坐标为东经125
°
45-126
°
30
′
，北纬46
°
12
′‑
46
°
22
′
。
23.如图1所示，为高频次重复凝视成像空基遥感载荷在轨测试方法流程包括：
24.步骤10、选择满足mtf和snr测试需求的场地，布设辐射特性靶标；
25.步骤20、监测太阳光谱信息和地球辐射收支信息，使用光谱仪对靶标及自然地物反射率进行测量；
26.步骤30、基于同步大气气溶胶观测数据，对gf-4连续观测的影像数据进行筛选，剔除大气扰动时刻获取的卫星影像数据；
27.步骤40、通过对线扩散函数lsf的计算与处理得到得到截止频率内不同频率的mtf值；
28.步骤50、对均匀场景地物的连续多帧影像数据在时间维度上进行列向量统计计算，获取均方根误差，基于配准后的多帧影像均匀地物计算snr。
29.本实施例以某研究区(图2)进行测试的结果为例，详细说明利用本发明高频次重复凝视成像空基遥感载荷在轨测试的过程。
30.如图3所示，基于gf-4高频重复凝视观测数据的脉冲法mtf测试技术包括：
31.1)布设辐射特性靶标
32.如图4所示，选取低反射率区域作为背景地物，布设高反射率辐射特性靶标(60％)。考虑到高分四号的地面分辨率(50米)及靶标成本、人工等因素，布设的高反射率辐射靶标至少覆盖3
×
4个像素，预计布设靶标约为180米
×
270米(48000平方米)。
33.2)地面同步反射率观测
34.高分四号卫星gf-4对试验区连续观测成像时候，地面开展同步光谱反射率试验，测量的靶标及背景地物反射率作为输入数据。
35.3)gf-4卫星高频次重复凝视观测数据高精度几何配准
36.如图5，gf-4卫星可在数分钟之内完成任务机动响应，并对用户指定的目标区域成像，包括对目标区域进行高频重复凝视观测、快速对某一区域进行拼接成像、对多个区域进行巡查成像。本研究主要基于gf-4卫星对同一目标区域进行高频重复凝视观测成像模式。对于单谱段成像时存在5s，不同谱段影像成像时存在50s的差异，被摄目标场景差异基本可以忽略，而卫星位置和姿态可能存在较小差异的特点，对高频重复凝视成像的数据进行高精度几何配准，经过处理后的单谱段多帧影像数据的配准精度(1σ)优于0.1～0.2个像素。
37.4)线扩散函数lsf获取
38.基于地面连续大气参数观测，去除大气扰动时刻的卫星影像数据，基于配准后的多帧影像恢复lsf线扩散函数。基于配准后的多帧gf-4卫星数据恢复lsf示意如图6。
39.5)计算mtf
40.计算恢复得到线扩散函数lsf；其次，对lsf函数经归一化处理；最后，对归一化的lsf进行离散傅里叶变换，并取变换后截止频率以内各分量的模值，并以0频率处的mtf值为基准进行归一化，得到截止频率内不同频率的mtf值，形成全频段mtf曲线。如下列公式所示：
41.mtf
x
＝normal(fft(normal(lsf
x
)))
42.其中，lsf
x
是x方向的线扩散函数；mtf
x
为x方向的调制传递函数；fft代表傅里叶变换处理；normal代表归一化处理。
43.2、基于gf-4高频重复凝视观测数据的snr测试
44.图7是基于gf-4高频重复凝视观测数据的snr测试的流程图，其方法步骤包括：
45.1)选取合适试验区域
46.选取不同反射率地物(低、中、高)的试验区，其中低反射率对应深海区，中反射率对应敦煌卫星辐射定标场的戈壁滩，高反射率对应敦煌石膏矿区的亮沙地，分别测试不同地物反射率条件下的信噪比。
47.2)制定gf-4卫星成像计划
48.针对选取的试验区域，制定成像计划，对均匀场景进行高频重复凝视观测。
49.3)gf-4卫星高频次重复凝视观测数据几何配准
50.本研究主要基于gf-4卫星对同一目标区域进行高频重复凝视观测成像模式。选择卫星影像上的同名点地物进行高精度的几何配准，经过处理后的单谱段多帧影像数据的配准精度(1σ)优于1个像素即可。
51.4)均匀地物在多帧影像时间维上列向量统计
52.对均匀场景地物进行长时间序列的变化趋势分析，对经过几何配准后的均匀场景的连续多帧影像数据在时间维度上进行列向量的统计计算，获取均方根误差，计算该成像条件下的snr。
53.本研究中基于gf-4高频重复凝视观测数据的脉冲法mtf测试技术和snr测试技术都取得了一定的结果；但是配准后多帧影像恢复的线扩散函数lsf计算mtf、基于配准后的多帧影像均匀地物计算snr的前提是要求很高的几何配准精度，下一步拟对多帧影像配准方法进行改进优化，利用地面及高分影像控制点数据库对gf4的多帧影像进行几何配准。
54.虽然本发明所揭露的实施方式如上，但所述的内容只是为了便于理解本发明而采
用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所揭露的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的专利保护范围，仍须以所附的权利要求书所界定的范围为准。