遥感卫星地球临边背景特性的红外成像仿真方法
【专利摘要】本发明公开了一种遥感卫星地球临边背景特性的红外成像仿真方法，包括如下步骤：根据遥感卫星轨道和成像传感器参数，通过模拟遥感卫星遥测场景，导出遥感卫星不同时刻在地球固联坐标系下的位置和速度信息；由位置和速度信息，根据地球固联坐标系、遥感卫星坐标系、成像传感器坐标系、图像坐标系之间的映射关系，得到不同时刻成像面与物面的相互转换关系；根据成像传感器坐标系到地球固联坐标系的旋转变换矩阵M和图像坐标系到成像传感器坐标系的转换关系矩阵M′，生成不同时刻、不同成像谱段下，遥感卫星地球临边背景红外能量图。本方法能对遥感卫星轨道的任意时刻进行地球临边背景特性的红外成像仿真，覆盖红外多波段，且算法简单，适应动平台。
【专利说明】遥感卫星地球临边背景特性的红外成像仿真方法
【技术领域】
[0001]本发明属于遥感图像处理【技术领域】，更具体地，涉及一种遥感卫星地球临边背景特性的红外成像仿真方法。
【背景技术】
[0002]在地球背景和深空背景边界存在渐变大气层，通常认为临边大气高度O到150km范围为临边背景。地球临边探测是遥感卫星重要的探测方式之一，因此，地球临边背景是遥感图像重要的组成部分。地球临边背景的研究越来越受关注，其辐射特性比地球背景更复杂，它包括大气、大气重力波、气辉、极光、云层等辐射，各种成分辐射特性的复杂性，导致地球临边背景的辐射特性更为复杂。由于地球临边背景的真实红外图像难以获得，同时，尚不存在全面的地球临边背景建模和仿真方法，因此，对地球临边背景特性的红外成像仿真尤为重要。
【发明内容】
[0003]针对现有技术的以上缺陷或改进需求，本发明提供了一种遥感卫星地球临边背景特性的红外成像仿真方法，其目的在于建立遥感卫星临边背景成像仿真模型，通过分析遥感卫星地球临边背景辐射特性，生成地球临边背景仿真图像，由此解决当今缺乏地球临边背景红外图像仿真方法的技术问题。该方法使临边背景特性的红外仿真更加准确，且模型简单，适用于动平台，可覆盖红外多波段仿真。
[0004]为实现上述目的，本发明提供了一种遥感卫星地球临边背景特性的红外成像仿真方法，其特征在于，包括如下步骤:
[0005]( I)根据遥感卫星轨道和成像传感器参数，通过模拟遥感卫星遥测场景，导出遥感卫星不同时刻在地球固联坐标系下的位置和速度信息；
[0006](2)由步骤(I)得到的遥感卫星不同时刻在地球固联坐标系下的位置和速度信息，根据地球固联坐标系、遥感卫星坐标系、成像传感器坐标系、图像坐标系之间的映射关系，得到不同时刻成像面与物面的相互转换关系，所述转换关系包括成像传感器坐标系到地球固联坐标系的旋转变换矩阵M和图像坐标系到成像传感器坐标系的转换关系矩阵M,；
[0007](3)根据M和M'，生成不同时刻、不同成像谱段下，遥感卫星地球临边背景红外能量图。
[0008]优选地，所述步骤(3)进一步包括如下步骤:
[0009](3-1)对于遥感卫星成像面上的任意一点P(u，V)，根据M'得到P点所在视线在传感器坐标系下的坐标(Xl，Y1, Z1)，再根据M得到P点所在视线在地球固联坐标系下的坐标(χ, y, z)；
[0010](3-2)在地球固联坐标系下，根据P点所在视线与地球的位置关系，判定P点对应的背景类型；
[0011](3-3)根据P点对应的背景类型，计算成像面上P点的辐照度E。[0012]优选地，所述步骤(3-2)的具体实现方式为:(Al)若P点所在视线与地球无交点，且切线高度大于150km，则认为P点对应深空背景；(A2)若P点所在视线与地球无交点，且切线高度小于或等于150km，则认为P点对应临边背景；(A3)若P点所在视线与地球有两个交点，则认为P点对应地球背景；(M)若P点所在视线与地球只有一个交点，则认为P点对应地球背景和临边背景的分界线。
[0013]优选地，所述步骤(3-3)中，物面背景区域对应的辐射强度按如下方法确定:(BI)地球背景:根据P点所在视线与地球的交点位置及地球背景能量图，计算其所在区域对应的辐射强度；(B2)临边背景:根据P点所在视线的切线高度，得到在该切线高度的包括大气、云层和极光在内的辐射强度；(B3)深空背景:其所在区域对应的辐射强度为O。
[0014]优选地，物面背景为临边背景时，云层的辐射强度按如下方法计算得到:(Cl)计算P点的太阳高度角h，判断h是否大于0，是则执行步骤(C2)，否则过程结束；(C2)计算P点所在视线与地球的切点所在位置的经纬度；(C3)根据经纬度选择云的覆盖率；(C4)根据云的覆盖率、厚度和高度，计算选定波段的辐射强度。
[0015]优选地，还包括步骤(4)，所述步骤(4)为:对获得的遥感背景能量图按如下公式进行量化处理，得到遥感卫星地球临边背景红外灰度图:
[0016]嫌/"她][./]=-D，
log(max)- log(min)
[0017]其中，inlmg[i] [j]为输入图像行列坐标为(i, j)的点的能量值，outImg[i] [j]为经过量化处理后输出的灰度值，max.min分别为输入图像中辐射能量的最大值和最小值，N为灰度图位数。
[0018]总体而言，通过本发明所构思的以上技术方案与现有技术相比，具有以下有益效果:
[0019]1、对临边状态背景复杂的红外辐射特性进行分析和研究，建立了地球临边背景辐射模型，解决了现有技术没有全面的地球临边背景建模和仿真方法的难题。
[0020]2、本方法考虑了包括大气、极光、重力波、云层在内的因素对地球临边背景辐射的影响，使临边背景特性的红外仿真更加准确。
[0021]3、本方法模型简单，适用于动平台，可覆盖红外多波段仿真。
【专利附图】
【附图说明】
[0022]图1是本发明实施例的遥感卫星地球临边背景特性的红外成像仿真方法流程图；
[0023]图2是各坐标系的示意图；
[0024]图3是遥感卫星成像面与物面映射示意图；
[0025]图4是地球面源到像元的辐照度计算模型；
[0026]图5是地球临边背景云层辐射强度计算方法流程图；
[0027]图6是白天2.7um地球临边背景红外灰度图；
[0028]图7是夜间2.7um地球临边背景红外灰度图；
[0029]图8是白天4.3um地球临边背景红外灰度图；
[0030]图9是夜间4.3um地球临边背景红外灰度图。【具体实施方式】
[0031]为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
[0032]如图1所示，本发明实施例的遥感卫星地球临边背景特性的红外成像仿真方法包括如下步骤:
[0033](1)根据遥感卫星轨道和成像传感器参数，通过模拟遥感卫星遥测场景，调节成像传感器光轴指向，使其能探测到临边背景，导出遥感卫星不同时刻在地球固联坐标系下的位置和速度信息。
[0034]如某轨道遥感卫星:轨道高度1600公里，倾角102.51 °，离心率0，视场
13。X13°。
[0035](2)由步骤(1)得到的遥感卫星不同时刻在地球固联坐标系下的位置和速度信息，根据地球固联坐标系、遥感卫星坐标系、成像传感器坐标系、图像坐标系之间的映射关系，得到不同时刻成像面与物面的相互转换关系。
[0036]各坐标系的示意图如图2所示，其中，地球固联坐标系:X轴指向本初子午线，Z轴指向正北，按照右手螺旋法则确定Y轴方向；遥感卫星坐标系:以卫星位置Os为原点，Zs始终指向地心，Xs指向轨道切线方向，Ys由右手定则确定；成像传感器坐标系:该坐标系原点
O。与卫星位置Os重合，成像传感器指向Z。的方向由在卫星坐标系下定义的坐标(久识)来确定，其中Θ为矢量Z。在XsOsYs平面上的投影与Xs正半轴的夹角，供定义为矢量Z。与Zs正半轴的夹角，X。轴方向由指向Z。与卫星坐标系XsOsYs平面的法向量叉乘得到，按照右手螺旋法则确定Y。轴；图像坐标系:卫星相机所成的图像以像素为单位，每个像素的坐标是该像素所在的行数和列数，像素坐标以(U，V)表示。
[0037]已知地心在地球固联坐标系的坐标为0(0，O, O)，卫星在地球固联坐标系下的坐标为Os (xs, ys, zs),卫星的速度矢量为ve (vex, vey, vez)。设卫星坐标系各坐标轴OXs, OYs, OZs
的方向向量为
【权利要求】
1.一种遥感卫星地球临边背景特性的红外成像仿真方法，其特征在于，包括如下步骤: (1)根据遥感卫星轨道和成像传感器参数，通过模拟遥感卫星遥测场景，导出遥感卫星不同时刻在地球固联坐标系下的位置和速度信息； (2)由步骤(1)得到的遥感卫星不同时刻在地球固联坐标系下的位置和速度信息，根据地球固联坐标系、遥感卫星坐标系、成像传感器坐标系、图像坐标系之间的映射关系，得到不同时刻成像面与物面的相互转换关系，所述转换关系包括成像传感器坐标系到地球固联坐标系的旋转变换矩阵M和图像坐标系到成像传感器坐标系的转换关系矩阵M,； (3)根据M和M'，生成不同时刻、不同成像谱段下，遥感卫星地球临边背景红外能量图。
2.如权利要求1所述的遥感卫星地球临边背景特性的红外成像仿真方法，其特征在于，所述步骤(3)进一步包括如下步骤: (3-1)对于遥感卫星成像面上的任意一点P(u，v)，根据M'得到P点所在视线在传感器坐标系下的坐标(Xl，Y1, Z1)，再根据M得到P点所在视线在地球固联坐标系下的坐标(χ, y, z)； (3-2)在地球固联坐标系下，根据P点所在视线与地球的位置关系，判定P点对应的背景类型； (3-3)根据P点对应的背景类型，计算成像面上P点的辐照度E。
3.如权利要求2所述的遥感卫星地球临边背景特性的红外成像仿真方法，其特征在于，所述步骤(3-2)的具体实现方式为:(Al)若P点所在视线与地球无交点，且切线高度大于150km，则认为P点对应深空背景；(A2)若P点所在视线与地球无交点，且切线高度小于或等于150km，则认为P点对应临边背景；(A3)若P点所在视线与地球有两个交点，则认为P点对应地球背景；(M)若P点所在视线与地球只有一个交点，则认为P点对应地球背景和临边背景的分界线。
4.如权利要求2所述的遥感卫星地球临边背景特性的红外成像仿真方法，其特征在于，所述步骤(3-3)中，物面背景区域对应的辐射强度按如下方法确定: (BI)地球背景:根据P点所在视线与地球的交点位置及地球背景能量图，计算其所在区域对应的辐射强度； (B2)临边背景:根据P点所在视线的切线高度，得到在该切线高度的包括大气、云层和极光在内的辐射强度； (B3)深空背景:其所在区域对应的辐射强度为O。
5.如权利要求4所述的遥感卫星地球临边背景特性的红外成像仿真方法，其特征在于，物面背景为临边背景时，云层的辐射强度按如下方法计算得到:(Cl)计算P点的太阳高度角h，判断h是否大于0，是则执行步骤(C2)，否则过程结束；(C2)计算P点所在视线与地球的切点所在位置的经纬度；(C3)根据经纬度选择云的覆盖率；(C4)根据云的覆盖率、厚度和高度，计算选定波段的辐射强度。
6.如权利要求1至5中任一项所述的遥感卫星地球临边背景特性的红外成像仿真方法，其特征在于，还包括步骤(4)，所述步骤(4)为:对获得的遥感背景能量图按如下公式进行量化处理，得到遥感卫星地球临边背景红外灰度图:
【文档编号】G01J5/00GK103743488SQ201310740551
【公开日】2014年4月23日 申请日期:2013年12月28日 优先权日:2013年12月28日
【发明者】张天序, 肖聪, 张俊青, 陈昌胜, 姚守悝 申请人:华中科技大学