静止气象卫星成像仪的地球固定网格映射方法技术

本发明专利技术公开了一种静止气象卫星成像仪的地球固定网格映射方法，包括如下步骤：步骤一、将地球表面特征物体映射到标称理想情况下静止气象卫星成像仪所成的图像中，将映射点作为地球固定网格；步骤二、通过比对实际遥感图像与固定网格的差异，作为L1级图像输出的参考基准，评价成像仪图像定位精度，分析成像仪内部光路变形特性和外部安装误差，分析卫星当前轨道位置；所述的标称理想情况包括如下描述：卫星位于赤道平面与定点经度半圈平面交线上、与地心距离为42164.172公里；成像仪光学基准坐标系的Z轴指向地心，X轴位于瞬时真赤道面内指向当地东方；成像仪内部光路准直，不存在畸变和失配。本发明专利技术可以提高静止气象卫星地面应用系统的数据处理效率和精度。

【技术实现步骤摘要】
静止气象卫星成像仪的地球固定网格映射方法
本专利技术涉及静止遥感卫星成像仪地球固定网格映射方法，具体涉及一种静止气象卫星成像仪的地球固定网格映射方法。
技术介绍
静止遥感卫星利用扫描式成像仪感知地球表面物体的辐射信息，可以用于地图测绘、天气诊断、灾害监测，以及生成其他定性、定量遥感产品。这需要成像仪遥感图像对应的绝对地理位置定位准确，相邻图像的相对位置关系配准精确。但卫星平台和遥感仪器受到天体力学、空间力学环境、热环境变化等因素的影响，卫星会存在轨道漂移和姿态指向偏差，成像仪会产生几何变形和失配，影响遥感图像的定位与配准精度。因此，静止遥感卫星的定量产品需要有一个固定基准，作为补偿校正成像仪视线、定量对比评估图像结果的参考依据。地球固定网格是静止气象卫星成像仪L1级数据的基准，美国、欧洲等均制定了关于固定网格的规范文件。但均由于历史的原因存在一定局限性，不能适应新一代静止遥感卫星成像仪生成图像的需求。气象卫星协调组(CoordinationGroupofMeteorologicalSatellites，缩写为CGMS)在其2013年发布的《通信传输数据全球规范》(LRIT/HRITGlobalSpecification)[1]中描述了地理坐标与卫星扫描角之间的投影关系。其固定网格描述为成像仪扫描角坐标；地表特征目标选用了海岸线；计算地表目标在地心地固系中的位置时采用WGS84椭球模型代表地球表面。欧洲气象卫星组织2015年公布了正在研制的MTG卫星主载荷FCI的数据集用户指南[2]。其中，对地理坐标投影的描述与文献[1]一致。美国Harris公司2015年公布了美国新一代静止气象卫星GOES-R的主载荷ABI的产品定义与用户手册[3]。文中描述了地理坐标到固定网格南北/东西角之间的映射关系。其固定网格描述为成像仪扫描角坐标；地表特征目标选用了海岸线；计算地表目标在地心地固系中的位置时采用WGS84椭球模型代表地球表面。针对中国风云二号静止气象卫星成像仪图像定位的论文[4]中描述了地球表面物体与自旋卫星扫描成像的对应关系，但未给出固定网格的定义和推导过程。文献[5]中采用定义卫星、地球和成像仪相关坐标系的方法描述地球表面物体与三轴稳定静止气象卫星成像仪坐标之间的关系，并在这种关系的基础上分析了各种姿态误差对成像结果的影响，以及基于卡尔曼滤波对误差进行处理的方法。综上，当前各国对静止气象卫星固定网格存在的缺陷主要包括以下几点：1)固定网格的地表特征目标仅选用海岸线，未选择内陆湖泊轮廓线。以至于在内陆地区缺少遥感图像参考或评价的特征目标。2)在计算投影关系时未考虑地表目标海拔高度，而是采用WGS84椭球模型代表地球表面。这将使某些高海拔地区目标的图像与其对应的网格投影存在偏差，对固定网格精度的影响最大可达4个像元。3)按照早期自旋卫星成像仪工作方式建立成像仪视线模型，不能更好的适应新一代静止气象卫星的带有两面扫描镜成像仪的成像特性。4)未建立从成像仪扫描镜转角到图像坐标的投影关系，也未建立从地理坐标到图像坐标的投影关系。以至于L1级数据产品仍需要进行角度坐标到图像坐标的转换处理。
技术实现思路
针对现有静止遥感卫星成像仪地球固定网格的缺陷，本专利技术提出了一种静止气象卫星成像仪的地球固定网格映射方法，可以提高静止气象卫星地面应用系统的数据处理效率和精度，为实现上述目的，本专利技术具体通过以下技术方案实现：静止气象卫星成像仪的地球固定网格映射方法，其特征在于，包括如下步骤：步骤一、将地球表面特征物体映射到标称理想情况下静止气象卫星成像仪所成的图像中，将映射点作为地球固定网格；步骤二、通过比对实际遥感图像与固定网格的差异，作为L1级图像输出的参考基准，评价成像仪图像定位精度，分析成像仪内部光路变形特性和外部安装误差，分析卫星当前轨道位置；所述的标称理想情况包括如下描述：卫星位于赤道平面与定点经度半圈平面交线上、与地心距离为42164.172公里；成像仪光学基准坐标系的Z轴指向地心，X轴位于瞬时真赤道面内指向当地东方；成像仪内部光路准直，不存在畸变和失配。所述的地球表面特征物体包括海岸线、内陆湖泊轮廓线，以及地球表面其他遥感辐射梯度较大的固定特征轮廓线。所述的静止气象卫星成像仪安装在静止轨道三轴稳定卫星平台上，通过二维扫描机构的旋转运动来往复移动视线，实现对地球的成像覆盖。本专利技术优选的静止气象卫星成像仪二维扫描机构，带有两面反射镜21,23和两根相互垂直的旋转轴20,22。从地球发出的光线首先入射在南北镜23上，经反射后再入射到东西镜21上，再次反射后进入内部成像系统。所述成像仪的光学基准坐标系的原点位于成像仪上的某个特征位置，Z轴为成像仪视场中心视线方向，X轴垂直于中心视线指向星下点图像正东方对应的方向，Y轴按照右手法则确定。所述步骤一具体包括如下步骤：首先将地球表面特征物体地理坐标映射到成像仪转角坐标，再将成像仪转角坐标映射到成像仪图像坐标；所述的地理坐标包括地理纬度φ、经度λ，海拔高度h。地理经纬度的单位为度，海拔高度的单位为公里。所述的成像仪转角坐标包括东西镜扫描角度ε和南北镜扫描角度η，单位为度。本专利技术优选的定义坐标零点和极性如下：在成像仪视线运动到视场中心视线时，成像仪转角坐标位于(0,0)点；定义南北镜绕成像仪光学基准坐标系X轴正方向按照右手方向旋转为ε的正方向，定义东西镜绕成像仪光学基准坐标系Z轴正方向按照右手方向旋转为η的正方向。所述的成像仪图像坐标包括南北方向像素号m和东西方向像素号n，单位为像素。本专利技术优选的定义图像左上角的像素为(1,1)点14，从上至下的方向为南北方向像素号递增方向，从左至右的方向为东西方向像素号递增方向。其中，地球表面特征物体地理坐标到成像仪转角坐标的映射，计算过程包括以下步骤：步骤1、根据地球表面特征物体8的地理坐标，按照地球的WGS84椭球模型计算其在地心地固坐标系中的位置。步骤2、计算卫星指向地球表面特征物体的视线矢量在地心地固坐标系中的分量。步骤3、计算卫星指向地球表面特征物体的成像仪视线向量。步骤4、根据成像仪视线向量与成像仪转角坐标的转换关系，计算地球表面特征物体对应的成像仪扫描镜转角坐标。所述的成像仪视线向量与成像仪转角坐标的转换关系包括两个转换：成像仪转角坐标(ε,η)转换为成像仪视线13向量(Lx,Ly,Lz)的表达式为：成像仪视线13向量(Lx,Ly,Lz)转换为成像仪转角坐标(ε,η)的表达式为：其中，成像仪转角坐标到成像仪图像坐标的映射，计算过程包括以下步骤：步骤1、根据成像仪的光学特性确定成像仪像素的南北方向张角κ和东西方向张角τ，并根据图像方位确定κ和τ的符号以调整成像仪图像与方位的对应关系，确定图像坐标(1,1)点对应的成像仪转角坐标点(ε1,η1)。结合本专利技术优选的成像仪图像坐标极性定义，本专利技术优选的将κ和τ取为负值。为适应静止轨道成像仪星下点分辨率为1公里的技术指标，本专利技术优选的将κ和τ同取为1/35786.035弧度。根据静止轨道成像仪视场和扫描机构运动行程，本专利技术优选的设置成像仪转角坐标(6,6)点对应图像坐标(1,1)点。步骤2、对于任意成像仪转角坐标[ε,η]，成像仪图像坐标的南北方向像素号m可表示为：步骤3、对于任意成像仪转角坐本文档来自技高网...

【技术保护点】
1.静止气象卫星成像仪的地球固定网格映射方法，其特征在于，包括如下步骤：步骤一、将地球表面特征物体映射到标称理想情况下静止气象卫星成像仪所成的图像中，将映射点作为地球固定网格；步骤二、通过比对实际遥感图像与固定网格的差异，作为L1级图像输出的参考基准，评价成像仪图像定位精度，分析成像仪内部光路变形特性和外部安装误差，分析卫星当前轨道位置；所述的标称理想情况包括如下描述：卫星位于赤道平面与定点经度半圈平面交线上、与地心距离为42164.172公里；成像仪光学基准坐标系的Z轴指向地心，X轴位于瞬时真赤道面内指向当地东方；成像仪内部光路准直，不存在畸变和失配。

【技术特征摘要】
1.静止气象卫星成像仪的地球固定网格映射方法，其特征在于，包括如下步骤：步骤一、将地球表面特征物体映射到标称理想情况下静止气象卫星成像仪所成的图像中，将映射点作为地球固定网格；步骤二、通过比对实际遥感图像与固定网格的差异，作为L1级图像输出的参考基准，评价成像仪图像定位精度，分析成像仪内部光路变形特性和外部安装误差，分析卫星当前轨道位置；所述的标称理想情况包括如下描述：卫星位于赤道平面与定点经度半圈平面交线上、与地心距离为42164.172公里；成像仪光学基准坐标系的Z轴指向地心，X轴位于瞬时真赤道面内指向当地东方；成像仪内部光路准直，不存在畸变和失配。2.如权利要求1所述的静止气象卫星成像仪的地球固定网格映射方法，其特征在于，所述的地球表面特征物体包括海岸线、内陆湖泊轮廓线，以及地球表面其他遥感辐射梯度较大的固定特征轮廓线。3.如权利要求1所述的静止气象卫星成像仪的地球固定网格映射方法，其特征在于，所述的静止气象卫星成像仪安装在静止轨道三轴稳定卫星平台上，通过二维扫描机构的旋转运动来往复移动视线，实现对地球的成像覆盖，该静止气象卫星成像仪二维扫描机构带有两面反射镜和两根相互垂直的旋转轴，从地球发出的光线首先入射在南北镜上，经反射后再入射到东西镜上，再次反射后进入内部成像系统。4.如权利要求1所述的静止气象卫星成像仪的地球固定网格映射方法，其特征在于，所述成像仪的光学基准坐标系的原点位于成像仪上的某个特征位置，Z轴为成像仪视场中心视线方向，X轴垂直于中心视线指向星下点图像正东方对应的方向，Y轴按照右手法则确定。5.如权利要求1所述的静止气象卫星成像仪的地球固定网格映射方法，其特征在于，所述步骤一具体包括如下步骤：首先将地球表面特征物体地理坐标映射到成像仪转角坐标，再将成像仪转角坐标映射到成像仪图像坐标；地理坐标包括地理纬度φ、经度λ，海拔高度h。地理经纬度的单位为度，海拔高度的单位为公里。成像仪转角坐标包含东西镜扫描角度ε和南北镜扫描角度η，单位为度，零点和极性如下...

【专利技术属性】
技术研发人员：吕旺，董瑶海，宋效正，沈毅力，周徐斌，陈文强，
申请(专利权)人：上海卫星工程研究所，
类型：发明
国别省市：上海,31