基于椭球分层大气折射的光学卫星影像精确对地定位方法技术

本发明专利技术提供一种基于椭球分层大气折射的光学卫星影像精确对地定位方法，包括建立全球任意点的高度与大气折射指数之间的对应关系；根据光线初始点，确定入射光线矢量到下一分层的参考椭球面的交点，确定下一分层上的折射光线矢量；通过对大气精细分层逐渐降低定位点的高度，迭代解算光线与参考椭球面的交点，直到与真实地形表面相交，得到基于椭球分层大气折射的光学卫星影像精确对地定位点。本发明专利技术提出对大气折射的处理，使得地面点坐标更加精确，并且提出通过计算折射点处参考椭球面的法向量实现准确提取大气折射影响下的分层折射光线矢量，利用大气精细分层、逐渐迭代解算的方式实现，因此实施简便，定位精度高。

Accurate positioning of optical satellite images based on ellipsoid stratified atmospheric refraction

【技术实现步骤摘要】
基于椭球分层大气折射的光学卫星影像精确对地定位方法
本专利技术属于光学卫星影像
，包括一种基于椭球分层大气折射的光学卫星影像精确对地定位方法。
技术介绍
现有的光学卫星直接对地定位模型大多数采用共线方程模型，没有考虑大气折射的影响[1]；有的虽然考虑了大气折射的影响，但是基于圆地球模型而不是椭球模型，或者仅仅计算了大气折射造成的地心张角，并且大气分层也非常粗糙、不符合实际[2，3]。本专利技术提出的基于椭球分层大气折射的定位模型克服了以上不足，可以用于光学卫星影像精确对地定位。参考文献：[1]单杰.光束法平差简史与概要[J].武汉大学学报(信息科学版),2018,43(12):1797-1810.[2]严明,王智勇,汪承义,于冰洋.大气折射对光学卫星遥感影像几何定位的影响分析[J].测绘学报,2015,44(09):995-1002.[3]田原.敏捷光学卫星全球无地面控制高精度对地定位方法研究[D].武汉大学.2018
技术实现思路
为了克服以上技术问题，本专利技术提供一种基于椭球分层大气折射的新的光学卫星影像精确对地定位方法。本专利技术技术方案提供一种基于椭球分层大气折射的光学卫星影像精确对地定位方法，包括以下步骤，步骤1，建立全球任意点的高度与大气折射指数之间的对应关系；步骤2，根据光线初始点，确定入射光线矢量到下一分层的参考椭球面的交点，实现方式如下，设成像时刻卫星位置(L0,B0)，其中L0为地理经度，B0为地理纬度，入射光线矢量为P0，光线所在初始高度为卫星高度H0，点(L0,B0,H0)处大气折射指数N0通过步骤1所得高度与大气折射指数之间的对应关系内插得到；设大气分层厚度为δH，计算入射光线矢量P0与高度为H1＝H0-δH的参考椭球面的交点(L1,B1)，得到计算经过点(L0,B0,H0)的光线矢量P0与高为H1的椭球面相交的点的坐标(L1,B1,H1)；步骤3，确定下一分层上的折射光线矢量，实现方式如下，点(L1,B1,H1)对应的大气折射指数N1通过步骤1所得高度与大气折射指数之间的对应关系内插得到；根据折射定律得到N0×sin(a0)＝N1×sin(a1)，其中a0和a1分别是高度H1处两个大气层的光线入射角和折射角；设折射点(L1,B1,H1)处参考椭球面的法向量为Q0，记折射光线矢量为P1，由P1与Q0夹角为a1，以及P1与P0夹角为a0-a1，实现光线矢量从P0到P1的折射方向的改变，计算得到折射光线矢量P1；步骤4，通过对大气精细分层逐渐降低定位点的高度，迭代解算光线与参考椭球面的交点，直到与真实地形表面相交，得到基于椭球分层大气折射的光学卫星影像精确对地定位点。而且，步骤1中，大气折射指数是基于多层大气物理量计算得到。而且，步骤2中，计算入射光线矢量P0与高度为H1＝H0-δH的参考椭球面的交点(L1,B1)，基于共线方程模型实现。或者，步骤2中，计算入射光线矢量P0与高度为H1＝H0-δH的参考椭球面的交点(L1,B1)，基于考虑光行差的几何成像模型实现。而且，步骤3中，计算得到折射光线矢量P1的实现方式如下，1)假设参考椭球面方程为，(X2+Y2)/A2+Z2/B2＝1A＝Ra+H1，B＝Rb+H1，式中，Ra、Rb分别为参考椭球的长半轴和短半轴，A、B为椭球面的长半轴和短半轴，H1为发生折射处的高度，(X,Y,Z)为椭球面上的点的地心直角坐标；根据曲面方程偏导数，得到该点(L1,B1,H1)处曲面法向量Q0＝[X0B2,Y0B2,Z0A2]，其中坐标(X0,Y0,Z0)为该点处的地心直角坐标，由地理坐标(L1,B1,H1)计算得到；2)计算入射光线P0与法向量Q0夹角a0＝acos(P0.*Q0)/|P0|/|Q0|，式中.*表示点积；折射角a1＝asin(N0×sin(a0)/N1)；其中，|P0|为P0的模，|Q0|为Q0的模；3)记折射光线矢量为P1，设P1＝[Px/Pz,Py/Pz,1]，由P1与Q0夹角为a1得：cos(a1)＝(P1.*Q0)/|P1|/|Q0|，由P1与P0夹角为a0-a1得：cos(a0-a1)＝(P1.*P0)/|P1|/|P0|，联立两式求得Px/Pz和Py/Pz，从而得到折射光线矢量P1；其中，|P1|为P1的模，Px、Py、Pz分别是是矢量P1的分量。本专利技术与现有技术相比的多项优点在于：(1)大气折射指数的计算是基于实际观测或计算所得多层大气物理量，和现有技术中采用通过理论模型计算得到的物理量相比，精度更高。(2)本专利技术在现有的共线方程模型基础上进行拓展，增加对大气折射的处理，使得地面点坐标更加精确。并且步骤2中，通过现有技术，计算经过点(L0,B0,H0)的光线矢量P0与高为H1的椭球面相交的点的坐标(L1,B1,H1)，是基于共线方程模型，但不限于共线方程模型(如，考虑光行差的几何成像模型，等)。(3)步骤3中，通过计算折射点处参考椭球面的法向量、实现光线矢量从P0到P1的折射方向的改变，从而准确提取大气折射影响下的分层折射光线矢量，提高了对地定位的精度，是本技术专利技术的主要特色之一。(4)步骤4中，通过对大气精细分层、逐渐降低定位点的高度，迭代解算光线与参考椭球面的交点，直到与真实地形表面相交，是本技术专利技术的主要特色之一。综上所述，本专利技术实施简便，高效可靠，可以推广适用于各种卫星影像产品应用中。附图说明图1为本专利技术实施例的原理示意图。具体实施方式为了便于本领域普通技术人员理解和实施本专利技术，下面结合附图及实施例对本专利技术作进一步的详细描述，应当理解，此处所描述的实施例仅用于说明和解释本专利技术，并不用于限定本专利技术。本专利技术实施例提出的基于椭球分层大气折射的光学卫星影像精确对地定位方法是在现有的共线方程模型基础上拓展而成的，通过进一步提出在计算像点对应的地面点坐标时，增加对大气折射的处理，使得地面点坐标更加精确，提供了一种新的光学卫星影像精确对地定位模型。实施例的方法实现包括如下步骤：(1)建立全球任意点的高度与大气折射指数之间的对应关系：首先获取成像时刻对应的多层大气物理量，包括：气压、气温、相对湿度等，利用现有的技术计算成像时刻对应的全球多层大气折射指数，从而建立了该成像时刻、全球任意点的高度与大气折射指数之间的对应关系(即任意点处，高度与大气折射指数一一对应)。具体实施时，计算全球多层大气折射指数实现可参见郑高峰.靶场测量大气折射率误差修正系统研究与设计[D].电子科技大学.2013.本专利技术不予赘述。本步骤大气折射指数的计算是基于多层大气物理量，这些值是实际观测或计算结果，是随时间动态变化的物理量。使用前本专利技术技术方案前，可以从网上下载。例如：https://rda.ucar.edu/datasets/ds083.2/#metadata/detailed.html？_do＝y&view＝level相比较而言，参考文献[2]中则是通过理论模型计算得到的物理量，比如大气温度是随高度变化的物理量，而不是真正观测得到的物理量。(2)根据光线初始点，确定入射光线矢量到下一分层的参考椭球面的交点：设成像时刻卫星位置(L0,B0)，其中L0为地理经度，B0为地理纬度，入射光线矢量为P0，(备注：根据光路可逆，入射光线和折射光线实际是反本文档来自技高网...

【技术保护点】
1.一种基于椭球分层大气折射的光学卫星影像精确对地定位方法，其特征在于：包括以下步骤，步骤1，建立全球任意点的高度与大气折射指数之间的对应关系；步骤2，根据光线初始点，确定入射光线矢量到下一分层的参考椭球面的交点，实现方式如下，设成像时刻卫星位置(L0,B0)，其中L0为地理经度，B0为地理纬度，入射光线矢量为P0，光线所在初始高度为卫星高度H0，点(L0,B0,H0)处大气折射指数N0通过步骤1所得高度与大气折射指数之间的对应关系内插得到；设大气分层厚度为δH，计算入射光线矢量P0与高度为H1＝H0‑δH的参考椭球面的交点(L1,B1)，得到计算经过点(L0,B0,H0)的光线矢量P0与高为H1的椭球面相交的点的坐标(L1,B1,H1)；步骤3，确定下一分层上的折射光线矢量，实现方式如下，点(L1,B1,H1)对应的大气折射指数N1通过步骤1所得高度与大气折射指数之间的对应关系内插得到；根据折射定律得到N0×sin(a0)＝N1×sin(a1)，其中a0和a1分别是高度H1处两个大气层的光线入射角和折射角；设折射点(L1,B1,H1)处参考椭球面的法向量为Q0，记折射光线矢量为P1，由P1与Q0夹角为a1，以及P1与P0夹角为a0‑a1，实现光线矢量从P0到P1的折射方向的改变，计算得到折射光线矢量P1；步骤4，通过对大气精细分层逐渐降低定位点的高度，迭代解算光线与参考椭球面的交点，直到与真实地形表面相交，得到基于椭球分层大气折射的光学卫星影像精确对地定位点。...

【技术特征摘要】
1.一种基于椭球分层大气折射的光学卫星影像精确对地定位方法，其特征在于：包括以下步骤，步骤1，建立全球任意点的高度与大气折射指数之间的对应关系；步骤2，根据光线初始点，确定入射光线矢量到下一分层的参考椭球面的交点，实现方式如下，设成像时刻卫星位置(L0,B0)，其中L0为地理经度，B0为地理纬度，入射光线矢量为P0，光线所在初始高度为卫星高度H0，点(L0,B0,H0)处大气折射指数N0通过步骤1所得高度与大气折射指数之间的对应关系内插得到；设大气分层厚度为δH，计算入射光线矢量P0与高度为H1＝H0-δH的参考椭球面的交点(L1,B1)，得到计算经过点(L0,B0,H0)的光线矢量P0与高为H1的椭球面相交的点的坐标(L1,B1,H1)；步骤3，确定下一分层上的折射光线矢量，实现方式如下，点(L1,B1,H1)对应的大气折射指数N1通过步骤1所得高度与大气折射指数之间的对应关系内插得到；根据折射定律得到N0×sin(a0)＝N1×sin(a1)，其中a0和a1分别是高度H1处两个大气层的光线入射角和折射角；设折射点(L1,B1,H1)处参考椭球面的法向量为Q0，记折射光线矢量为P1，由P1与Q0夹角为a1，以及P1与P0夹角为a0-a1，实现光线矢量从P0到P1的折射方向的改变，计算得到折射光线矢量P1；步骤4，通过对大气精细分层逐渐降低定位点的高度，迭代解算光线与参考椭球面的交点，直到与真实地形表面相交，得到基于椭球分层大气折射的光学卫星影像精确对地定位点。2.根据权利要求1所述保持原始分辨率的遥感虚拟线阵参数提取方法，其特征在于：步骤1中，大气折射指数是基于多层大气物理量计算得到。3.根据权利要求1所述保持原始分辨率的遥感虚拟线阵参数提取...

【专利技术属性】
技术研发人员：金淑英，王艳丽，王密，程宇峰，
申请(专利权)人：武汉大学，
类型：发明
国别省市：湖北,42