一种地球同步轨道合成孔径雷达干涉的高程反演方法
【专利摘要】本发明提供一种地球同步轨道合成孔径雷达干涉的高程反演方法，其包括：步骤1，选取GEO?SAR获取干涉数据的轨道，采集轨道上的干涉数据；步骤2，根据步骤1获得的干涉数据利用BP算法进行GEO?SAR成像处理；步骤3，根据经步骤2处理后的GEO?SAR图像建立GEO干涉模型，当相位矢量与成像平面分离时，根据GEO干涉模型进行GEO?SAR干涉高程反演。本发明的高程反演方法通过建立合理的GEO?SAR干涉模型，实现有效相位矢量与成像平面分离情况下的高程反演处理，解决了GEO?SAR干涉处理的高程反演的核心问题——有效相位矢量与成像平面分离问题，实现了任意位置处利用GEO?SAR干涉处理的高程反演。
【专利说明】一种地球同步轨道合成孔径雷达干涉的高程反演方法
【技术领域】
[0001]本发明属于合成孔径雷达【技术领域】，尤其涉及一种地球同步轨道合成孔径雷达干涉的高程反演方法。
【背景技术】
[0002]合成孔径雷达(SAR)是一种全天候、全天时的高分辨率微波遥感成像雷达，可安装在飞机、卫星、导弹等飞行平台上。自上世纪50年代发明以来，已经在很多领域取得了越来越广泛的应用，例如灾害控制、植被分析、微波遥感等领域。地球同步轨道合成孔径雷达(GEO SAR)是运行在36000km高度地球同步椭圆轨道上的SAR卫星。相比于低轨SAR(LE0SAR，轨道高度低于IOOOKm)而言，GEO SAR具有成像范围大、重访时间短、抗打击与抗摧毁能力强等特点，目前已成为国内外的研究热点。
[0003]其中干涉处理是GEO SAR研究的一个重要方面。GEO SAR干涉处理的最大难点在于GEO SAR高轨道高度、较大的轨道空变形、严重的摄动影响以及较大的等效前斜角度，使得GEO SAR的有效相位矢量不在成像平面中。对于传统成像算法近直线运动轨迹而言，带高程的地面点在进行二维成像时会投向地面上的某一点。这两个点与轨道上垂直处的点构成了成像平面。然而对于目标点而言，它的有效相位中心位于孔径中心，因此孔径中心与目标连线产生的有效相位矢量并不位于成像平面内。这种分离会导致传统的低轨高程反演的干涉模型严重失效(有时可以相差10倍)；同时由于GEO干涉具有大场景干涉处理的特点，但对于大场景的边缘的目标点采用低轨高程反演的干涉模型进行处理时，同样会使计算得到的模糊高度与真实值相差较大，这会导致一个多周期的目标场景的高程误差在几十米甚至上百米以上。因此，GEO SAR干涉处理的核心问题就是如何解决有效相位矢量与成像平面分离的问题，并且达到较好的精度，这在现有的各种GEO SAR干涉处理算法中并未有提及。
【发明内容】
[0004]为解决上述问题，本发明提供一种地球同步轨道合成孔径雷达干涉的高程反演方法，该高程反演方法能够通过建立合理的GEO SAR干涉模型，实现有效相位矢量与成像平面分离情况下的高程反演处理，实现任意位置处利用GEOSAR干涉处理的高程反演。
[0005]本发明的地球同步轨道合成孔径雷达干涉的高程反演方法包括:
[0006]步骤1，选取GEO SAR获取干涉数据的轨道，采集所述轨道上的干涉数据:
[0007]步骤2，根据步骤I获得的干涉数据利用BP算法进行GEO SAR成像处理；
[0008]步骤3，根据经步骤2处理后的GEO SAR图像建立GEO干涉模型，当相位矢量与成像平面分离时，根据所述GEO干涉模型进行GEO SAR干涉高程反演；
[0009]进一步的，所述步骤I包括:
[0010]步骤11，固定卫星的波束指向和下视角，在卫星成像工作时间对应的轨道区间上选取距轨道起始工作时刻和终止工作时刻大于半个合成孔径时间的任意一个位置做为合成孔径中心位置，确定指定视角对应的地面点的经纬度，并将该地面点设定为场景中心[0011]利用合成孔径时间Ts和脉冲重复时间PRT确定所述第一轨轨道上每个脉冲发射时刻的卫星位置；
[0012]步骤12，根据两轨之间的时间间隔确定第二轨的位置；
[0013]进一步的，所述步骤2包括:
[0014]步骤21，根据所述步骤I采集的干涉数据建立雷达回波信号模型，该雷达回波信号模型如式(I)所示:.(/4irl? [ ( IR Y),,、
【权利要求】
1.一种地球同步轨道合成孔径雷达干涉的高程反演方法，其特征在于，包括: 步骤1，选取GEO SAR获取干涉数据的轨道，采集所述轨道上的干涉数据； 步骤2，根据步骤I获得的干涉数据利用BP算法进行GEO SAR成像处理； 步骤3，根据经步骤2处理后的GEO SAR图像建立GEO干涉模型，当相位矢量与成像平面分离时，根据所述GEO干涉模型进行GEO SAR干涉高程反演； 进一步的，所述步骤I包括: 步骤11，固定卫星的波束指向和下视角，在卫星成像工作时间对应的轨道区间上选取距轨道起始工作时刻和终止工作时刻大于半个合成孔径时间的任意一个位置做为合成孔径中心位置，确定指定视角对应的地面点的经纬度，并将该地面点设定为场景中心点； 利用合成孔径时间Ts和脉冲重复时间PRT确定所述第一轨轨道上每个脉冲发射时刻的卫星位置； 步骤12，根据两轨之间的时间间隔确定第二轨的位置； 进一步的，所述步骤2包括: 步骤21，根据所述步骤I采集的干涉数据建立雷达回波信号模型，该雷达回波信号模型如式⑴所示:
2.如权利要求1所述的地球同步轨道合成孔径雷达干涉的高程反演方法，其特征在于，所述步骤3包括: 步骤31，根据步骤2获得进行GEO SAR成像处理后的两幅GEO SAR图像，通过所述两幅GEO SAR图像的互相关法计算偏移量插值，然后根据该偏移量插值进行图像配准； 步骤32，配准完成后，将所述两幅GEO SAR图像的主图像和辅图像共轭相乘得到干涉相位图，并对所述干涉相位图进行去平地处理得到反映地形高度变化的稀疏干涉条纹； 步骤33，利用所述干涉相位图及稀疏干涉条纹生成相关系数图； 步骤34，对所述相关系数图采用Goldstein相位滤波方法进行滤除干涉相位噪声处理，采用相位解缠方法进行相位解缠处理； 步骤35，根据经步骤34处理后的相关系数图，将相位矢量与成像平面分离获得GEOSAR干涉模型，然后将主辅卫星位置投影到成像平面进行分析，根据所述GEO SAR干涉模型推导获得高程反演表达式进行高程反演: 步骤36，通过地理编码将所述GEO SAR干涉模型及高程反演表达式中的斜距方位坐标转化到地距方位坐标。
3.如权利要求1所述的地球同步轨道合成孔径雷达干涉的高程反演方法，其特征在于，所述步骤12包括: 步骤121，利用卫星通过所述第一轨轨道的孔径中心位置、所述场景中心点及该卫星下点位置所组成平面的时刻确定第二轨的孔径中心时刻； 步骤122，以第二轨的孔径中心对应的卫星时刻为基准时刻向卫星轨道的前后时间方向每一个脉冲重复时间PRT在卫星轨道上确定一个第二轨轨道上每个脉冲发射时刻的卫星位置，直到延伸至基准时刻的前后半合成孔径时间长度。
【文档编号】G01S13/90GK103543453SQ201310517927
【公开日】2014年1月29日 申请日期:2013年10月28日 优先权日:2013年10月28日
【发明者】胡程, 龙腾, 李元昊, 李延, 曾涛, 丁泽刚 申请人:北京理工大学, 北京空间飞行器总体设计部