[发明专利]一种基于矩阵计算的卫星遥感图像网格共定位方法

权利要求书

1.一种基于矩阵计算的卫星遥感图像网格共定位方法，其特征在于，包括以下步骤：

获取待共定位卫星遥感源图像；

构建共定位网格坐标文件；

将待共定位卫星遥感源图像输入至共定位网格坐标文件，得到共定位后的目标图像；

所述构建共定位网格坐标文件，包括：

获取相同网格和分辨率的卫星遥感源图像；

根据卫星遥感源图像得到源图像经纬度矩阵；

定义目标图像经纬度矩阵并将其进行重塑以及切分，得到多个切片；

将多个网格点进行合并，得到多个批；

基于批对批的方法对多个切片进行索引循环，得到多个索引矩阵，并将多个索引矩阵导出，得到共定位网格坐标文件。

2.如权利要求1所述的一种基于矩阵计算的卫星遥感图像网格共定位方法，其特征在于，所述卫星遥感图像的形状为(m,n)，源图像经纬度矩阵形状为(m,n,2)；目标图像经纬度矩阵的形状为(s,t,2)，重塑后的目标图像经纬度矩阵的形状(s×t,2)，m、s均为宽度，n、t均为高度。

3.如权利要求2所述的一种基于矩阵计算的卫星遥感图像网格共定位方法，其特征在于，根据内存或显存大小选择批的大小，根据批的大小将重塑后的目标图像经纬度矩阵沿第一个轴进行切分，得到多个形状为(batch,2)的切片。

4.如权利要求3所述的一种基于矩阵计算的卫星遥感图像网格共定位方法，其特征在于，根据下式选择批的大小：

式中，Memory_req表示程序需要的显存或内存大小，batch表示批的大小，且为正整数，其取值范围为[1,s×t]。

5.如权利要求3所述的一种基于矩阵计算的卫星遥感图像网格共定位方法，其特征在于，通过将s×t的值与batch的值做除法，得到切片的数量，若batch不能被s×t整除，则最后一次循环的切片形状为(s×t％batch,2)，其中％表示取余运算。

6.如权利要求3所述的一种基于矩阵计算的卫星遥感图像网格共定位方法，其特征在于，基于批对批的方法对多个切片进行索引循环，得到多个索引矩阵，并将多个索引矩阵导出，得到共定位网格坐标文件，包括以下步骤：

定义一个元素个数为s×t全零列向量Indices0；

定义循环变量i，0≤i≤N，且N＝s×t；

输入循环变量i，当i＝0时，从源图像经纬度矩阵中查找与第一个切片中batch个点距离最近的batch个点，获得batch个点在源图像中的索引矩阵；

将索引矩阵存放在Indices0中索引位置为i～i+batch处，得到Indices1；

当循环变量i大于N时，将得到的Indices1导出为共定位网格坐标文件。

7.如权利要求6所述的一种基于矩阵计算的卫星遥感图像网格共定位方法，其特征在于，当循环变量i小于N时，将循环变量的值增加batch，从源图像经纬度矩阵中查找与下一个切片中batch个点距离最近的batch个点，获得新的索引矩阵，将Indices1进行更新。

8.如权利要求7所述的一种基于矩阵计算的卫星遥感图像网格共定位方法，其特征在于，将待共定位卫星遥感源图像输入至共定位网格坐标文件，得到共定位后的目标图像，包括：

将源图像展平为包括m×n个元素的列向量，根据Indices1得到一个元素个数为s×t的目标列向量，将其重塑成形状为(s,t)的二维矩阵，得到共定位后的目标图像(s,t)。