亚星游戏官网-www.yaxin868.com

专利名称：：一种基于星地坐标转换的星载sar图像地理编码方法
技术领域：
：本发明涉及一种图像地理编码方法，特别涉及一种基于星地坐标转换的星载SAR图像地理编码方法。属于图像处理
技术领域：
。
背景技术：
：合成孔径雷达(SAR)卫星近些年来发展迅速，由于SAR卫星不受天气、地理、时间等因素的限制，能够对地进行全天时的观测，且具有一定的穿透力，因而被广泛的应用于军事侦察、资源探测、海洋观测、生态监测、快速救援等方面。随着应用需求的不断增加，星载SAR图像产品的等级也不断提高，其中对星载SAR图像产品进行地理编码是一个重要的星载SAR图像处理内容，具有广泛的应用前景，比如快速营救、精确打击、资源探测、图像拼接等。对于星载SAR图像，目前常用的星载SAR图像地理编码基本可以分为两类1、基于距离-多普勒(R-D)定位模型的方法基于距离-多普勒(R-D)定位模型的方法是一种精确的SAR图像地理编码方法，其核心是利用了3个方程求解SAR图像中的像素在转动地心坐标系下的坐标，进而获取该像素位置所对应的经纬度，3个方程分别为地球椭球模型方程、斜距方程及多普勒方程。由于这3个方程表示形式较为复杂，故最后的问题归结到求解非线性方程组的问题。一般情况下，非线性方程组不具有显示解，需要利用数值分析的方法求解近似值，而当初值选取不合适时，结果可能不收敛，同时迭代处理需要消耗过多的时间，严重影响了SAR图像的地理编码效率。2、基于近似模型的方法正由于R-D定位模型需要迭代运算，处理效率低，在实际应用中往往对处理模型进行修改，在不影响应用的前提下对模型进行各种近似，进而获取某种条件下的显示解，这种方法虽然提高了处理效率，伹处理精度受到了很大的影响，难以实现星载SAR图像的高精度地理编码。本方法解决了上述的问题，通过建立星地坐标之间的6组转换关系，利用雷达天线视角这一中间变量，将问题转换为求解一元四次方程的问题，在不损失精度的前提下显著地提高了处理的效率，是一种全新的星载SAR图像地理编码方法，对实现星载SAR图像高精度实时地理编码具有重要意义。
发明内容本发明的目的是提供一种基于星地坐标转换的星载SAR图像地理编码方法，它是通过星地坐标转换矩阵，将雷达天线视角作为未知量带入星地坐标转换矩阵，将其转换到转动地心坐标系中，并带入地球椭球模型方程进行求解，通过求解一元四次方程，求解出雷达天线视角大�。�再将雷达天线视角带入星地坐标转换矩阵即可求解SAR图像中某像素所对应的经讳度，依次对星载SAR图像中的像素点重复上述操作，即可快速、精度的完成星载SAR图像的地理编码。为了更好对本方法进行介绍，首先需要介绍星地6个坐标系之间的转换关系。首先说明星地转换6个坐标系的定义1.不转动的地心坐标系￡。坐标原点地球球心Z轴沿地球的自转轴指向正北极X轴在赤道平面内，指向春分点Y轴在赤道平面内，使该坐标系构成右手直角坐标系不转动的地心坐标系是在不转动的地心赤道参考系(惯性参考系)中建立的直角坐标系。2.转动的地心坐标系^坐标原点地球球心Z轴沿地球的自转轴指向正北极X轴在赤道平面内，通过格林威治子午线上半分支Y轴在赤道平面内，使该坐标系构成右手坐标系3.卫星轨道平面坐标系^坐标原点卫星椭圆轨道的一个焦点(即地球球心)Z轴垂直于卫星轨道平面，正向指向卫星的角动量矢量方向Y轴在卫星平面内，正向指向近心点X轴在卫星轨道平面内，使该坐标系构成右手直角坐标系4.卫星平台坐标系A坐标原点卫星质心Z轴垂直于卫星轨道平面，正向指向卫星的角动量矢量方向X轴在卫星轨道平面内，陀螺平台住纵轴方向(卫星的设计分性方向)Y轴在卫星轨道平面内，使该坐标系构成右手直角坐标系5.卫星星体坐标系i^,坐标原点卫星质心；X轴沿卫星星体纵轴方向(卫星的真实飞行方向)Y轴、Z轴沿卫星星体的另外两个惯性主轴方向6.天线坐标系E。坐标原点天线相位中心点；X轴正向指向卫星的真实飞行方向；Y轴沿天线瞄准线，指向地球方向为正向；Z轴右手准则给出，使该坐标系构成右手直角坐标系。其次说明星地六个坐标系之间的转换关系4、4、4、4、4、4。44。坐标系转换示意图如上图所示，4^为4。的逆矩阵(4。^。g=/)，以下类1.转动的地心坐标系^/不转动的地心坐标系A不转动的地心坐标系绕轴逆时针转过一个春分点的格林威治时角就得到转动的地心坐标系￡g。<formula>formulaseeoriginaldocumentpage7</formula>2.不转动的地心坐标系A/轨道平面坐标系A不转动的地心坐标系五。经三次旋转得到轨到平面坐标系A。第一次，将不转动的地心坐标系绕Z轴逆时针旋转一个角Q;第二次，将得到的坐标系再绕X同<轴逆时针旋转一个角度"第三次，再将所得坐标系绕z轴逆时针旋转一个角度最后得到轨道平面坐标系《。<formula>formulaseeoriginaldocumentpage8</formula>3.轨道平面坐标系A/卫星平台坐标系A卫星平台坐标系&绕Z轴逆时针旋转一个角度90°+^-r得到卫星轨道平面坐标系A。其中，S是卫星的真近心角，y是卫星的航迹角。真近心角0由开普勒方程及真近心角与偏心角E的对应关系求得-"tan—=1+e五—、-tan—22航迹角/HS90。sin(6>——cos(^—y)0人00l乂其中，A为引力场常数，取3.986013e14，r为星上某一时刻。4.卫星平台坐标系A/卫星形体坐标系A卫星星体坐标系A经三次旋转得到卫星平台坐标系A。第一次，将卫星星体坐标系《绕X轴顺时针旋转一个横滚角度《；第二次，将得到的坐标系绕Z轴顺时针旋转一个俯仰角度^;第三次，再将所得到的坐标系绕Y轴逆时针旋转一个角度^，最后得到卫星平台坐标系A。<formula>formulaseeoriginaldocumentpage8</formula>5.卫星星体坐标系&/天线坐标系^。天线坐标系￡。绕X轴逆时针旋转一个角度^得到卫星星体坐标系&<formula>formulaseeoriginaldocumentpage9</formula>注意这些变换矩阵的转置矩阵就是逆矩阵(综上所述，本发明一种基于星地坐标转换的星载SAR图像地理编码方法，其具体操作步骤如下步骤一依照星载SAR图像数据产品说明书，对相关信息副产品进行读�。�获取所需的星历参数，包括偏心率"半长轴"、升交点赤径Q、近地点俯角、轨道倾角Z、波前斜距i，、信号釆样率/、、脉冲重复频率J^/、过近地点时刻r、光速c、格林威治时角//c;本步骤的实施条件在于需要提供SAR图像相关信息的副产品，卫星星历参数都将包含在该副产品内，目前国际上标准SAR图像数据都提供相关信息的副产品。步骤二依照星载SAR图像数据产品说明书，读入星载SAR图像数据，计算图像每一列所对应的斜距大�。�相关计算公式如下，其中，i」代表第j个距离门对应的斜距，A，代表第一个距离门对应的斜距，乂为采样率，c为光速，j为距离门编号<formula>formulaseeoriginaldocumentpage9</formula>本步骤的实施条件在于获取相关SAR图像信息，尤其是波前斜距/，、采样率乂.和SAR图像大小。步骤三对星载SAR图像中某一像素点(i，j)，其中i为行数、j为列数，依照星载SAR图像数据产品说明书，对相关信息副产品进行读�。�获取其对应的星上时间，计算该时刻下卫星位置矢量在转动地心坐标系下的坐标，相关计算公式如下，其中，(、、,々、,，、)代表卫星位置矢量在转动地心坐标系下的坐标，《，。分别代表了不转动地心坐标系到转动地心坐标系的转换矩阵和轨道坐标系到不转动地心坐标系的转换矩阵，^代表真近心角，M代表平均真近心角，r代表极矢径，"代表半长轴，//代表地球引力常数，r代表过近心点时刻，f代表星上时间，e代表偏心率:<formula>formulaseeoriginaldocumentpage10</formula>本步骤的实施条件在于获取相关SAR图像信息，尤其是偏心率"半长轴"、升交点赤径Q、近地点俯角w、轨道倾角/、脉冲重复频率、过近地点时刻r、光速c、格林威治时角//6。步骤四对星载SAR图像中某一像素点(i,j)，其中i为行数、j为列数，依照星载SAR图像数据产品说明书，对相关信息副产品进行读�。�获取其对应的星上时间、三轴姿态控制角(A，^,《)、偏心率e、半长轴"、升交点赤径Q、近地点俯角w、轨道倾角/、脉冲重复频率^/、过近地点时刻r、光速c、格林威治时角&，计算该时刻下卫星天线相位中心位置矢量(&,XA)在转动地心坐标系下的坐标(x併,;^，zj,相关计算公式如下，其中，(&，K,z。)代表卫星天线相位中心在卫星星体坐标系下的坐标，(x^尺w,^)代表卫星天线相位中心在转动地心坐标系下的坐标，《,、<v、4、々'分别代表不转动地心坐标系到转动地心坐标系的转换矩阵、轨道坐标系到不转动地心坐标系的转换矩阵、卫星平台坐标系到轨道坐标系的转换矩阵、卫星星体坐标系到卫星平台坐标系的转换矩阵<formula>formulaseeoriginaldocumentpage10</formula>本步骤的实施条件在于获取相关SAR图像信息，尤其是偏心率e、半长轴"、升交点赤径Q、近地点俯角w、轨道倾角/、脉冲重复频率、过近地点时刻f、光速c、格林威治时角/^，同时还需要获取卫星天线相位中心坐标。步骤五对星载SAR图像中某一像素点(i，j)，其中i为行数、j为列数，依照星载SAR图像数据产品说明书，对相关信息副产品进行读�。�获取其对应的星上时间、三轴姿态控制角(^,,《)、偏心率e、半长轴"、升交点赤径Q、近地点俯角w、轨道倾角"脉冲重复频率^/、过近地点时刻r、光速c、格林威治时角，计算该时刻下卫星星体坐标系到转动地心坐标系的转换矩阵4^,，相关计算公式如下，其中，^w代表由卫星星体坐标系到转动地心坐标系的转换矩阵，《。、4一4,、4。分别代表不转动地心坐标系到转动地心坐标系的转换矩阵、轨道坐标系到不转动地心坐标系的转换矩阵、卫星平台坐标系到轨道坐标系的转换矩阵、卫星星体坐标系到卫星平台坐标系的转换矩阵足"=」,'".,^4""A。=V仏"2"U"32"33乂本步骤的实施条件在于获取相关SAR图像信息，尤其是偏心率e、半长轴"、升交点赤径Q、近地点俯角w、轨道倾角/、过近地点时刻r、光速c、格林威治时角Z/e。步骤六对星载SAR图像中某一像素点(i,j)，其中i为行数、j为列数，计算其对应的斜距，求解一元四次方程组，解出此时的雷达天线视角，相关计算公式如下，其中，i」代表第j个距离门对应的斜距，i，代表第一个距离门对应的斜距，,为采样率，c为光速，j为距离门编号，^代表地球长半轴，^代表地球短半轴，^代表雷达天线视角,附，+w42)x4十2(/^m2+/w4m5)x3+(2^/773+/n22—m42+/w52)x2+2(m2;w3—/w4m5)x3+(附32—m52)=072如122+222+A:."322-"132-232_&'a332)m2=2.,(a12.6,+a22.62+&.a32.63)m3=7y2.(a132+a232+A:'a"2)+(6卩+622+632-^2)附<》25JC=COS6,本步骤的实施条件在于需要考虑边界条件，避免求解值超出(0,;r/2)范围<步骤七将求解出的雷达天线视角带入星地坐标转换矩阵，求解星载SAR图像某一像素点(i,j)在转动地心坐标系下的坐标，进而求出该像素点所对应的经纬度(八,O)，相关计箅公式如下，其中，《,、<v、4、4、4。分别代表不转动地心坐标系到转动地心坐标系的转换矩阵、轨道坐标系到不转动地心坐标系的转换矩阵、卫星平台坐标系到轨道坐标系的转换矩阵、卫星星体坐标系到卫星平台坐标系的转换矩阵、天线坐标系到卫星星体坐标系的转换矩阵，i,代表第j个距离门对应的斜距，,&、,zft、)代表卫星位置矢量在转动地心坐标系下的坐标，,)代表天线指向点在转动地心坐标系下的坐标:<formula>formulaseeoriginaldocumentpage12</formula>本步骤的实施条件在于获取相关SAR图像信息，尤其是偏心率e、半长轴"、升交点赤径Q、近地点俯角"、轨道倾角/、过近地点时刻r、光速c、格林威治时角Z/e。步骤八每完成一个像素点的地理编码，计算i-M和j-N的值，其中M和N分别代表SAR图像的行数和列数，判断此时j-N是否小于O，如果小于0，则从步骤六进行重复操作；如果大等于O,则判断此时i-M是否小于0，如果小于0，则转入图像下一行处理，从步骤2进行重复操作，否则星载SAR图像编码完成。本步骤的实施条件在于硬件设施所能开辟的内存大小不能小于SAR图像大小。本发明一种基于星地坐标转换的星载SAR图像地理编码方法的优点在于(1)本方法能够利用星地转换关系，通过求解雷达天线视角进一步求解图像像素所对应的经纬度，避免了直接求解所需的迭代运算，显著提高了处理效率。(2)本方法能够利用星地转换关系，通过求解一元四次方程来得到精确的地理编码结果，显著提高了处理结果的精度。(3)本方法可避免迭代运算，因而避免了传统方法中由于采用迭代运算而导致结果不收敛的缺陷。(4)本方法计算结果精度不受卫星速度影响。传统地理编码方法中，卫星速度是影响定位精度的最大因素，而采用本方法将不受卫星速度误差的影响，即当星历参数不理想时，采用本方法所得到的地理编码误差将远小于传统方法的地理编码误差。(5)本方法计算效率不受初始迭代值影响。传统地理编码方法中，其编码效率受初始迭代值影响，当初始迭代值选取不合适时，其编码效率将很低，甚至出现结果不收敛、效率为0的状况，而本方法可避免这种状况，在任何初始条件下都具有更高的处理效率。(6)本方法对星载SAR图像采用逐行处理方式进行图像地理编码，由于星载SAR图像行与行之间的处理过程完全独立，故可对星载SAR图像进行分块处理，进一步提高处理效率。图1基于星地坐标转换的星载SAR图像地理编码方法的流程图图2SAR图像说明示意图图3计算卫星位置矢量在转动地心坐标系下坐标的流程图图4计箅卫星天线相位中心位置矢量在转动地心坐标系的流程图图5计算转动地心坐标系到卫星星体坐标系的转换矩阵的流程图图6计算求解雷达天线视角的流程图图7仿真场景摆放示意图图中符号说明如下4。不转动地心坐标系到转动地心坐标系的转换矩阵^。v轨道坐标系到不转动地心坐标系的转换矩阵Ar卫星平台坐标系到轨道坐标系的转换矩阵Ae卫星星体坐标系到卫星平台坐标系的转换矩阵A。天线坐标系到卫星星体坐标系的转换矩阵及，波前距离，即第一个距离门所对应的斜距w,第j个距离门所对应的斜距^第1行SAR图像数据对应的星上时间,,第i行SAR图像数据对应的星上时间(&。,力,^)目标像素点在转动地心坐标系下的坐标具体实施例方式见图1、图2、图3、图4、图5、图6、图7所示，本发明一种基于星地坐标转换的星载SAR图像地理编码方法，以一幅大小为2048x2048的SAR图像为例，该方法具体步骤如下步骤一读入卫星星历参数，包括偏心率e=0.003、半长轴a=7071140.0m、升交点赤径Q-69.517572°、近地点俯角"=0.0°、轨道倾角!'=99.038314°、波前斜距^,=787796.306附、信号采样率,.=75.0M他、脉冲重复频率p/=2300.0/fe、过近地点时刻r--9,6465s、光速c二3.0e8w/"格林威治时角//G=O.Os。本步骤的实施条件在于需要提供SAR图像相关信息的副产品，卫星星历参数都将包含在该副产品内，目前国际上标准SAR图像数据都提供相关信息的副产品。步骤二读入星载SAR图像，计算图像每一列所对应的斜距大小。星载SAR成像机制在距离向是一种距离分辨机制，其每一距离门(即星载SAR图像每一列，如图2所示)中的像素对应相同的斜距，即像素所代表的目标距天线相位中心具有相同的距离。星载SAR图像第j列所代表的斜距可按下式进行计算，其中，/，代表第一个距离门所对应的斜距(即图像的第一列)，其中，,=75,0MJ/z,imin=787796.306m，_2-乂本步骤的实施条件在于获取相关SAR图像信息，尤其是波前斜距i■=787796.306m、采样率,=75.0MZ/z禾[JSAR图像大小。步骤三对星载SAR图像中某一像素点(i，j)，其中i为行数、j为列数，计算该时刻下卫星位置矢量在转动地心坐标系下的坐标，流程图见图3，以1=1024，j-1024为例，转动地心坐标系下的坐标(、、，，、)求取方法如下式<formula>formulaseeoriginaldocumentpage15</formula>其中，M代表真近心角，M^(r-r)々(一3)，〃为引力常数，取3.986013e14。r=a'G—e2)/(l+e.cose)=7049927.6994m(4)本步骤的实施条件在于获取相关SAR图像信息，尤其是偏心率e-0.003、半长轴fl二7071140.0m、升交点赤径0=69.517572°、近地点俯角"=0.0°、轨道倾角/=99.038314°、波前斜距i,-787796.306m、信号采样率乂.=75.0M/fe、脉冲重复频率=2300.0//z、过近地点时刻r--9.6465s、光速c^3.0e8w/"格林威治时角/^;=0.0"步骤四对星载SAR图像中某一像素点(i，j)，其中i为行数、j为列数，计算该时刻下卫星天线相位中心在转动地心坐标系下的坐标，流程图见图4。首先读取卫星天线相位中心在卫星星体坐标系下的坐标(xt,,k,^)及三轴姿态控制角(《.，，《)，一般情况下，该坐标取(O,O,O)，并利用坐标转换矩阵4、之、4、4将其转换到转动地心坐标系下的坐标(~,;^，、)，以1=1024，)=1024为例，则(^,凡,^)=(0.0,0.0，0.0)，<formula>formulaseeoriginaldocumentpage15</formula>(5)本步骤的实施条件在于获取相关SAR图像信息，尤其是偏心率e-0.003、半长轴a=7071140.0m、升交点赤径Q=69.517572°、近地点俯角《=0.0°、轨道倾角!'=99.038314°、波前斜距/画=787796.306w、信号采样率/、.=75.0M//Z、脉冲重复频率jw/=2300.0//z、过近地点时刻r--9.6465s、光速c二3.0e8m/"格林威治时角/^-0.(Xy，同时还需要获取卫星天线相位中心坐标(^^4)=(0.0,0.0,0.0)及三轴姿态控制角^,^,《)=(3.79°,0.0°，0.0°)。步骤五对星载SAR图像中某一像素点(i,j)，其中i为行数、j为列数，计算该时刻下卫星星体坐标系到转动地心坐标系的转换矩阵，流程图见图5，以i=1024，=1024为例，为了便于表示，将卫星星体坐标系到转动地心坐标系的转换矩阵用4e表示，则<formula>formulaseeoriginaldocumentpage16</formula>(6)本步骤的实施条件在于获取相关SAR图像信息，尤其是偏心率e0.003、半长轴"=7071140.0m、升交点赤径Q=69.517572°、近地点俯角"=0.0。、轨道倾角"99.038314°、波前斜距i,=787796.306^1、信号采样率乂=75.0M历、脉冲重复频率jw/=2300.07/z、过近地点时刻7=-9.6465s、光速c=3.0e8w〃、格林威治时角//(;=0.0"步骤六对星载SAR图像中某一像素点(i,j),其中i为行数、j为列数，根据公式1计算其对应的斜距，求解一元四次方程组，解出此时的雷达天线视角，流程图见图6，以i-1024，j-1024为例，在转动地球坐标系下，卫星位置矢量、天线相位中心位置矢量、目标到天线相位中心的距离矢量、目标点矢量形成了闭合的四边形，故其矢量和为0，故有，<formula>formulaseeoriginaldocumentpage16</formula>(7)其中，对于SAR图像每一行而言，后两项为固定值，为了便于分析，作如下进一步化简，<formula>formulaseeoriginaldocumentpage16</formula>(8)其中，《即为雷达天线视角，(x^^。,^,)为目标点在转动地心坐标系下的坐标。由于(。,斗，~。)为转动地心坐标系下的坐标，故其满足地球椭球方程，<formula>formulaseeoriginaldocumentpage16</formula>(9)」、万分别代表了地球的长、短半轴，分别取^=6378140.0mS=6356755.0m，将^g。,^。,z伊)用《表示，并代入表达式(8)，则可化简为关于cos^的一元四次方程，其中，(mj+m42)x4十2(mm2+m4m5)x3+(2/^/^+m22—m42+m52)x2,、(10)附,=7,2.(a122+a222+A:■"322—a132—q232—A:'a332)m4=27—/.("12"13+"22"23+A."32a33)=2.7,("13W2+"3363)j、，附5■X=cos《L则问题最终转换为求解一元四次方程，其求解过程可参考各种数学手册即可，最终解出A=30.0°。本步骤的实施条件在于需要考虑边界条件，避免求解值超出(0，;r/2)范围。步骤七将求解出的雷达天线视角带入星地坐标转换矩阵，求解星载SAR图像某一像素点(i,j)在转动地心坐标系下的坐标。此时由于cos&和sin《已经求出，将其代入公式(7)即可求解出(xg。,;v^)，则其对应的经度为A=atan&=73.071556°(11)、'纬度为C>=atango。+力。+zs。乂=0.205607°(12)本步骤的实施条件在于获取相关SAR图像信息，尤其是偏心率6=0.003、半长轴fl=7071140.0m、升交点赤径Q=69.517572°、近地点俯角^;=0.0°、轨道倾角,'=99.038314°、波前斜距imin=787796.306m、信号采样率/、.=75.0M他、脉冲重复频率=、过近地点时刻7=-9.6465s、光速c=3.0e8w/s、格林威治时角/^-0.0s。步骤八判断此时j是否小于图像列数，如果小于图像列数，则从步骤六进行重复操作；如果大于图像列数，则判断此时i是否小于图像行数，如果小于图像行数，则转入图像下一行处理，从步骤二进行重复操作，否则星载SAR图像编码完成。本步骤的实施条件在于硬件设施所能开辟的内存大小不能小于SAR图像大小。为了说明该方法的有效性，进行如下仿真试验，仿真参数如表l。点阵摆放如图7所示，在地球表面摆放3x3的点阵共9个点，每一个点所代表的经炜度如表2:表1仿真参数<table>tableseeoriginaldocumentpage18</column></row><table><table>tableseeoriginaldocumentpage19</column></row><table>表5两种地理编码方法对比结果<table>tableseeoriginaldocumentpage19</column></row><table>从表3中可以看出，用基于星地坐标转换的方法对目标进行地理编码，其结果可以精确到小数点第6位，具有非常高的精度；从表4中可以看出，每一个点的地理编码误差不超过0.4m，也验证本方法的精确性；从表5中可以看出，对一个2048x2048的SAR图像进行地理编码，达到同样数量极的地理编码精度，基于R-D模型的地理编码方法消耗的时间约为基于星地坐标转换的地理编码方法的67倍，且SAR图像越大，两种方法消耗的时间之比也将越来越大，这也验证了本方法的快速性。权利要求1、一种基于星地坐标转换的星载SAR图像地理编码方法，其特征在于，它包括如下步骤步骤一依照星载SAR图像数据产品说明书，对相关信息副产品进行读�。�获取所需的星历参数，包括偏心率e、半长轴a、升交点赤径Ω、近地点俯角ω、轨道倾角i、波前斜距R_min、信号采样率fs、脉冲重复频率prf、过近地点时刻τ、光速c、格林威治时角HG；步骤二依照星载SAR图像数据产品说明书，读入星载SAR图像数据，计算图像每一列所对应的斜距大�。�相关计算公式如下，其中，R_j代表第j个距离门对应的斜距，R_min代表第一个距离门对应的斜距，fs为采样率，c为光速，j为距离门编号步骤三对星载SAR图像中某一像素点(i，j)，其中i为行数、j为列数，依照星载SAR图像数据产品说明书，对相关信息副产品进行读�。�获取其对应的星上时间，计算该时刻下卫星位置矢量在转动地心坐标系下的坐标，相关计算公式如下，其中，(xgs，ygs，zgs)代表卫星位置矢量在转动地心坐标系下的坐标，Ago·Aov分别代表了不转动地心坐标系到转动地心坐标系的转换矩阵和轨道坐标系到不转动地心坐标系的转换矩阵，θ代表真近心角，M代表平均真近心角，r代表极矢径，a代表半长轴，μ代表地球引力常数，τ代表过近心点时刻，t代表星上时间，e代表偏心率r＝a·(1-e2)/(1+e·cosθ)步骤四对星载SAR图像中某一像素点(i，j)，其中i为行数、j为列数，依照星载SAR图像数据产品说明书，对相关信息副产品进行读�。�获取其对应的星上时间、三轴姿态控制角(θy，θp，θr)、偏心率e、半长轴a、升交点赤径Ω、近地点俯角ω、轨道倾角i、脉冲重复频率prf、过近地点时刻τ、光速c、格林威治时角HG，计算该时刻下卫星天线相位中心位置矢量(xe，ye，ze)在转动地心坐标系下的坐标(xge，yge，zge)，相关计算公式如下，其中，(xe，ye，ze)代表卫星天线相位中心在卫星星体坐标系下的坐标，(xge，yge，zge)代表卫星天线相位中心在转动地心坐标系下的坐标，Ago、Aov、Avr、Are分别代表不转动地心坐标系到转动地心坐标系的转换矩阵、轨道坐标系到不转动地心坐标系的转换矩阵、卫星平台坐标系到轨道坐标系的转换矩阵、卫星星体坐标系到卫星平台坐标系的转换矩阵步骤五对星载SAR图像中某一像素点(i，j)，其中i为行数、j为列数，依照星载SAR图像数据产品说明书，对相关信息副产品进行读�。�获取其对应的星上时间、三轴姿态控制角(θy，θp，θr)、偏心率e、半长轴a、升交点赤径Ω、近地点俯角ω、轨道倾角i、脉冲重复频率prf、过近地点时刻τ、光速c、格林威治时角HG，计算该时刻下卫星星体坐标系到转动地心坐标系的转换矩阵Age，相关计算公式如下，其中，Age代表由卫星星体坐标系到转动地心坐标系的转换矩阵，Ago、Aov、Avr、Are分别代表不转动地心坐标系到转动地心坐标系的转换矩阵、轨道坐标系到不转动地心坐标系的转换矩阵、卫星平台坐标系到轨道坐标系的转换矩阵、卫星星体坐标系到卫星平台坐标系的转换矩阵步骤六对星载SAR图像中某一像素点(i，j)，其中i为行数、j为列数，计算其对应的斜距，求解一元四次方程组，解出此时的雷达天线视角，相关计算公式如下，其中，R_j代表第j个距离门对应的斜距，R_min代表第一个距离门对应的斜距，fs为采样率，c为光速，j为距离门编号，A代表地球长半轴，B代表地球短半轴，θL代表雷达天线视角m2＝2·R_j·(a12·b1+a22·b2+k·a32·b3)m5＝2·R_j·(a13b1+a23b2+k·a33b3)x＝cosθL步骤七将求解出的雷达天线视角带入星地坐标转换矩阵，求解星载SAR图像某一像素点(i，j)在转动地心坐标系下的坐标，进而求出该像素点所对应的经纬度(Λ，Φ)，相关计算公式如下，其中，Ago、Aov、Avr、Are、Aea分别代表不转动地心坐标系到转动地心坐标系的转换矩阵、轨道坐标系到不转动地心坐标系的转换矩阵、卫星平台坐标系到轨道坐标系的转换矩阵、卫星星体坐标系到卫星平台坐标系的转换矩阵、天线坐标系到卫星星体坐标系的转换矩阵，R_j代表第j个距离门对应的斜距，(xgs，ygs，zgs)代表卫星位置矢量在转动地心坐标系下的坐标，(xgo，ygo，zgo)代表天线指向点在转动地心坐标系下的坐标步骤八每完成一个像素点的地理编码，计算i-M和j-N的值，其中M和N分别代表SAR图像的行数和列数，判断此时j-N是否小于0，如果小于0，则从步骤六进行重复操作；如果大等于0，则判断此时i-M是否小于0，如果小于0，则转入图像下一行处理，从步骤2进行重复操作，否则星载SAR图像编码完成。全文摘要本发明提出了一种基于星地坐标转换的星载SAR图像地理编码方法，它是利用星地坐标之间的6个坐标系关系，巧妙的借用观测视角这一中间变量，在不需要迭代运算的前提下完成对星载SAR图像的高精度快速地理编码。传统的星载SAR图像地理编码方法大多需要进行迭代运算，处理效率低，在计算中是采用近似处理，影响定位精度，特别是卫星速度误差，它是影响传统星载SAR图像地理编码结果精度的最大因素，而本方法却不受卫星速度误差的影响，在星历参数非理想时，本方法的定位精度远远优于传统定位的方法，因而更适合于实际工程应用。本方法是一种全新的星载SAR图像定位方法，对实现星载SAR图像快速高精度地理编码具有重要意义。文档编号G01S7/295GK101414003SQ200810226999公开日2009年4月22日申请日期2008年11月28日优先权日2008年11月28日发明者周荫清,李春升,威杨,王鹏波,杰陈申请人:北京航空航天大学