[发明专利]大尺度地距星载SAR图像镶嵌方法

权利要求书

1.大尺度地距星载SAR图像镶嵌方法，包括多幅已编号地距星载SAR图像，其特征在于：该方法包括以下步骤：

步骤1、图像预处理；

对每个已编号地距星载SAR图像进行参数修正后获得含地理定位信息的待处理SAR图像；

步骤2、计算每个待处理SAR图像中每个像素的经、纬度坐标；

步骤3、判定所有待处理SAR图像两两之间是否重叠；

根据步骤2所得到的每个待处理SAR图像中每个像素的经、纬度坐标，判断所有待处理SAR图像两两之间是否重叠；

存在重叠的，记录所有两两之间存在重叠的待处理SAR图像对应的编号以及两两之间重叠区域范围，并生成含有重叠区域范围信息的编号组进行存储；

步骤4、计算每个编号组对应重叠区域范围之间的匹配偏移量；

对于峰均比低的区域，采用互相关法计算每个编号组对应重叠区域范围之间的匹配偏移量；

对于峰均比高的区域，采用基于过零点的二值匹配统计法来计算每个编号组对应重叠区域范围之间的匹配偏移量；

步骤5、创建全零像素图像和全局结构体变量

根据步骤2中得到的待处理SAR图像中每个像素的经、纬度坐标和步骤4得到的每个编号组对应重叠区域范围之间的匹配偏移量获得所有待处理SAR图像修正后的地理坐标，得到修正后每个待处理SAR图像对应的待拼接图像；

创建与最终要生成的镶嵌图像同样大小的全零像素图像，同时建立一个全局结构体变量；

步骤6、采用Wallis滤波器对步骤5获得的每个编号组对应的重叠区域范围进行匀色；

步骤7、对匀色后的所有待拼接图像进行拼接；

拼接方法如下：

a.在全局结构体变量中为所有待拼接图像设置拼接匀色标记，并复位为0；

b.输入一幅待拼接图像，基于编号组找出所有与该幅待拼接图像重叠且匀色标记为0的待拼接图像；

c.假设与当前输入待拼接图像重叠且匀色标记为0的待拼接图像总数为Num，以当前输入待拼接图像为参考图像，依次按照步骤6对这Num幅待拼接图像进行Wallis匀色，且每一次匀色均更新Wallis滤波器中相应的参数，将匀色标记置为1并保存匀色后的图像；

d.Num次Wallis匀色后，将该输入待拼接图像写入步骤5中创建的全零像素图像中，并将拼接匀色标记设置为1；

e.对所有待拼接图像重复上述b、c、d步骤，处理完所有待拼接图像后，得到最终的镶嵌结果。

2.根据权利要求1所述的大尺度地距星载SAR图像镶嵌方法，其特征在于，步骤1中所述参数修正具体为如下过程：

对每个已编号地距星载SAR图像进行相干斑滤波，以去除相干斑噪声，获得无相干斑噪声的SAR图像；

对无相干斑噪声的SAR图像进行辐射定标，以消除系统辐射误差，并将定标后的SAR图像用后向散射系数表示，获得用后向散射系数表示的SAR图像；

对用后向散射系数表示的SAR图像进行椭球体校正，并将辅助信息中的地理定位信息写入SAR图像，得到待处理SAR图像。

3.根据权利要求2所述的大尺度地距星载SAR图像镶嵌方法，其特征在于，步骤3中每个待处理SAR图像中每个像素的经、纬度坐标的计算公式如下：

E_geo＝A(0)+X_pix×A(1)+Y_pix×A(2) (1)

N_geo＝A(3)+X_pix×A(4)+Y_pix×A(5) (2)

其中：E_geo和N_geo分别表示所述待处理SAR图像中像素的经度和纬度坐标；

X_pix和Y_pix分别表示所述待处理SAR图像中像素的列坐标和行坐标；

所述地理定位信息具体由地理定位数组构成，所述地理定位数组具体由六个参数A(0),A(1),...A(5)构成，其中：

A(0)和A(3)表示分别表示所述待处理图像左上角的经、纬度坐标；

A(1)和A(5)分别表示所述待处理图像的横向与纵向分辨率；

A(2)和A(4)分别表示所述待处理图像中像素在纬、经度方向上的旋转系数。