基于先验DEM数据的星载SAR图像几何精校正方法技术

本发明专利技术公开了一种基于先验DEM数据的星载SAR图像几何精校正方法，其主要思路为：获取待校正的实际星载SAR图像，并建立待校正的实际星载SAR图像的距离多普勒模型，记为成像模型，然后利用所述成像模型对待校正的实际星载SAR图像中的P×N个像素点分别进行定位，分别得到待校正的实际星载SAR图像P×N个像素点各自的定位坐标；根据待校正的实际星载SAR图像P×N个像素点各自的定位坐标，进而计算待校正的实际星载SAR图像P×N个像素点各自的电磁波入射角，以及最终反演的理想星载SAR图像Sima；分别计算待校正的实际星载SAR图像相对于Sima在方位向的像素偏差量na，以及待校正的实际星载SAR图像相对于Sima在距离向的像素偏差量nr，进而计算得到几何精校正后的实际星载SAR图像。

Geometric precision correction method for Spaceborne SAR images based on prior DEM data

The invention discloses a SAR DEM data based on prior image geometric correction method, the main idea is: access to the actual star to be corrected SAR image, range Doppler model and the establishment of an actual star to be corrected SAR images, recorded as imaging model, the actual location and then treated respectively by star correction the imaging model of the load in the SAR image of P * N pixel image, P * N pixel coordinates of their actual star to be corrected SAR respectively; according to the image of P * N pixel coordinates of their actual star to be corrected loaded SAR * N pixel electromagnetic the actual incident wave and then calculate the star borne SAR image correction of P angle, and the ideal star of the inversion SAR Sima image; image relative to the Sima in the azimuth pixel deviation Na actual star to be corrected were calculated to load SAR And the corrected deviation of the actual spaceborne SAR image relative to the Sima in the range direction NR, and then calculate the actual spaceborne SAR images after geometric correction.

【技术实现步骤摘要】
基于先验DEM数据的星载SAR图像几何精校正方法
本专利技术属于信号处理
，特别涉及一种基于先验DEM数据的星载SAR图像几何精校正方法，适用于减小星载SAR因为系统误差造成的SAR图像场景中的像素点成像位置偏离真实位置点的问题，从而改善星载SAR的成像质量，有助于用户对SAR图像的判读和解译。
技术介绍
星载SAR是从空间对地观测的一种有效手段，具有全天时、全天候、多波段及高分辨率的优点，能够详细、准确地测绘地形、地貌，获取地球表面的信息，最终生成目标场景的高分辨地图。但是星载SAR在从雷达脉冲的发射一直到最后通过成像算法生成星载SAR图像的过程中，会有很多非理想因素误差引入其中，这样会使得实际得到的星载SAR图像中每个像素点对应的场景中位置和实际场景中的位置都会有一定的偏移，使得实际得到的星载SAR图像与真实场景并不完全对应，进而无法详细、准确地测绘地形、地貌，获取地球表面的信息。
技术实现思路
针对上述现有技术存在的不足，本专利技术的目的在于提出一种基于先验DEM数据的星载SAR图像几何精校正方法，该种基于先验DEM数据的星载SAR图像几何精校正方法使用星载SAR回波数据进行本文档来自技高网...

【技术保护点】
一种基于先验DEM数据的星载SAR图像几何精校正方法，其特征在于，包括以下步骤：步骤1，获取待校正的实际星载SAR图像，所述待校正的实际星载SAR图像为P×N维矩阵，P表示待校正的实际星载SAR图像的方位向包含的像素点个数，N表示待校正的实际星载SAR图像的距离向包含的像素点个数，P、N分别为大于0的正整数；步骤2，建立待校正的实际星载SAR图像的距离多普勒模型，记为成像模型，然后利用所述成像模型对待校正的实际星载SAR图像中的P×N个像素点分别进行定位，分别得到待校正的实际星载SAR图像中P×N个像素点各自的定位坐标；步骤3，根据待校正的实际星载SAR图像中P×N个像素点各自的定位坐标，计算...

【技术特征摘要】
1.一种基于先验DEM数据的星载SAR图像几何精校正方法，其特征在于，包括以下步骤：步骤1，获取待校正的实际星载SAR图像，所述待校正的实际星载SAR图像为P×N维矩阵，P表示待校正的实际星载SAR图像的方位向包含的像素点个数，N表示待校正的实际星载SAR图像的距离向包含的像素点个数，P、N分别为大于0的正整数；步骤2，建立待校正的实际星载SAR图像的距离多普勒模型，记为成像模型，然后利用所述成像模型对待校正的实际星载SAR图像中的P×N个像素点分别进行定位，分别得到待校正的实际星载SAR图像中P×N个像素点各自的定位坐标；步骤3，根据待校正的实际星载SAR图像中P×N个像素点各自的定位坐标，计算得到待校正的实际星载SAR图像中P×N个像素点各自的电磁波入射角；步骤4，根据待校正的实际星载SAR图像中P×N个像素点各自的电磁波入射角对待校正的实际星载SAR图像进行反演，得到最终反演的理想星载SAR图像；步骤5，分别计算待校正的实际星载SAR图像相对于最终反演的理想星载SAR图像在方位向的像素偏差量na，以及待校正的实际星载SAR图像相对于最终反演的理想星载SAR图像在距离向的像素偏差量nr；步骤6，根据待校正的实际星载SAR图像相对于最终反演的理想星载SAR图像在方位向的像素偏差量na，以及待校正的实际星载SAR图像相对于最终反演的理想星载SAR图像在距离向的像素偏差量nr，对待校正的实际星载SAR图像进行几何精校正，得到几何精校正后的实际星载SAR图像。2.如权利要求1所述的一种基于先验DEM数据的星载SAR图像几何精校正方法，其特征在于，还包括:步骤7，对几何精校正后的实际星载SAR图像进行地理编码，得到地理编码后的几何精校正星载SAR图像。3.如权利要求1所述的一种基于先验DEM数据的星载SAR图像几何精校正方法，其特征在于，步骤2的子步骤为：2.1建立待校正的实际星载SAR图像的距离多普勒模型，在所述距离多普勒模型中，卫星沿着自身轨道飞行，且卫星自带的GPS会对星载雷达的空间位置坐标进行记录，星载雷达位于卫星上，星载雷达的空间位置包括M个点，M个点各自的坐标都是WGS84坐标系中的坐标，WGS84坐标系的原点是地球的中心；赤道和0度经线的交点，与地球中心的连线构成WGS84坐标系的X轴；北极点和地球的中心连线构成WGS84坐标系的Z轴，再由右手坐标系原则确定出WGS84坐标系的Y轴；将星载雷达的空间位置包括的M个点中第i个点记为Si，i＝1,2,3,...,M，第i个点的坐标为(Xi,Yi,Zi)，Xi表示第i个点Si在WGS84坐标系的X轴坐标值，Yi表示第i个点Si在WGS84坐标系的Y轴坐标值，Zi表示第i个点Si在WGS84坐标系的Z轴坐标值；将M个点各自的坐标记为星载雷达POS数据，星载雷达POS数据对应包括M个POS数据；将M个点各自的坐标进行逐点连接后形成星载雷达航迹；根据星载雷达POS数据，并由速度和位移关系，计算得到星载雷达的空间位置中第l个点处的瞬时速度l＝1,2,3,...,M-1，PRF表示为脉冲重复频率,Sl+1表示第l+1个点,Sl表示第l个点,M表示星载雷达的空间位置包含的点个数,与POS数据个数取值相等；星载雷达在卫星轨道上飞行的过程中，会向地面发射电磁波，星载雷达发射的电磁波照射到地面上时也会形成一个椭圆的照射区域，随着星载雷达从第1个POS数据到第M个POS数据运动，该椭圆的照射区域会在地面上移动形成一个长方形条带，将所述长方形条带记为SAR图像成像区域；将与星载雷达航迹平行的SAR图像成像区域的对应方向,记作SAR图像成像区域的方位向；将与星载雷达航迹垂直的SAR图像成像区域的对应方向,记作SAR图像成像区域的距离向，且SAR图像成像区域对应待校正的实际星载SAR图像，待校正的实际星载SAR图像的方位向包含P个像素点，待校正的实际星载SAR图像的距离向包含N个像素点，则待校正的实际星载SAR图像为P×N维矩阵；2.2将待校正的实际星载SAR图像中第m行、第n列处的像素点记作Tmn，Tmn的坐标记作(Xmn,Ymn,Zmn),Xmn表示待校正的实际星载SAR图像中第m行、第n列处的像素点在WGS84坐标系的X轴坐标值，Ymn表示待校正的实际星载SAR图像中第m行、第n列处的像素点在WGS84坐标系的Y轴坐标值，Zmn表示待校正的实际星载SAR图像中第m行、第n列处的像素点在WGS84坐标系的Z轴坐标值；星载雷达空间位置包含的点个数与待校正的实际星载SAR图像的方位向包含的像素数个数取值相同，且一一对应；2.3使用迭代法对待校正的实际星载SAR图像中第m行、第n列处的像素点Tmn进行定位，Tmn的坐标为(Xmn,Ymn,Zmn)，n＝1,2,3,…,N，m＝1,2,3,…,P，P表示待校正的实际星载SAR图像的方位向包含的像素点个数，N表示待校正的实际星载SAR图像的距离向包含的像素点个数，m和n的初始值分别为1；初始化：令f表示第f次迭代，f的初始值为1，令h0表示待校正的实际星载SAR图像中每个像素点的海拔高度预设值，其中每个像素点分别为SAR图像场景中每个像素点；设定待校正的实际星载SAR图像中包含J个像素点，令j表示第j个像素点，j∈{1,2,…,J}，j的初始值为1，J＝P×N；2.4根据第f-1次迭代后第j个像素点的海拔高度预设值hf-1，得到如下方程组：其中，Xi表示第i个点Si在WGS84坐标系的X轴坐标值，Xjf表示第f次迭代后第i个点Si在WGS84坐标系的X轴坐标值，Yi表示第i个点Si在WGS84坐标系的Y轴坐标值，Yjf表示第f次迭代后第i个点Si在WGS84坐标系的Y轴坐标值，Zi表示第i个点Si在WGS84坐标系的Z轴坐标值，Zjf表示第f次迭代后第i个点Si在WGS84坐标系的Z轴坐标值，表示第f次迭代后第i个点Si和待校正的实际星载SAR图像中第j个像素点之间的斜距，fdjf表示第f次迭代后第i个点Si和待校正的实际星载SAR图像中第j个像素点之间的斜距对应的多普勒频率，表示星载雷达的空间位置中第i个点处的瞬时速度矢量，表示第f次迭代后第i个点Si和待校正的实际星载SAR图像中第j个像素点之间的斜距矢量，→表示矢量，λ表示星载雷达载频对应的波长，||表示求绝对值操作，Re表示地球的赤道半径，hj(f-1)表示第f-1次迭代后第j个像素点的海拔高度预设值，Rp表示地球的极地半径；求解以上方程组，计算得到第f次迭代后第j个像素点的坐标(Xjf,Yjf,Zjf)，所述第f次迭代后第j个像素点的坐标(Xjf,Yjf,Zjf)是在WGS84坐标系中的，将(Xjf,Yjf,Zjf)转换到经纬度坐标系中，得到第f次迭代后第j个像素点的经纬度坐标(Ljf,Bjf,Hjf)，其中Ljf为第f次迭代后第j个像素点的经度，Bjf为第f次迭代后第j个像素点的纬度，Hjf为第f次迭代后第j个像素点的海拔高度；然后利用SRTM库中的DEM数据确定地球上经纬度(Ljf,Bjf)对应的海拔高度Hjf'，Hjf'为经纬度(Ljf,Bjf)对应的标准海拔高度,且此时的第f-1次迭代后第j个像素点的海拔高度预设值hf-1和Hjf'不相等；2.5如果|hf-1-Hjf'|≥ε0,ε0表示设定的门限阈值，则令f加1，将经纬度(Ljf,Bjf)对应的海拔高度Hjf'作为第f-1次迭代后第j个像素点的海拔高度预设值hf-1，返回2.4；如果|hf-1-Hjf'|<ε0，则迭代停止，并将迭代停止时得到的第f次迭代后第j个像素点的坐标(Xjf,Yjf,Zjf)，作为待校正的实际星载SAR图像中第j个像素点的定位坐标。2.6令j加1，返回2.4，直到得到待校正的实际星载SAR图像中第1个像素点的定位坐标至待校正的实际星载SAR图像中第J个像素点的定位坐标，此时求得了待校正的实际星载SAR图像...

【专利技术属性】
技术研发人员：李真芳，田锋，邢超，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西,61