基于查找表的遥感影像快速重投影方法与相关技术

图片简介:

本技术介绍了一种基于查找表的遥感影像快速重投影方法，涉及遥感影像处理技术领域，包括以下步骤：预先生成四个空查找表LUX1、LUY1、LUX2和LUY2，并记录每个像素重投影后坐标信息；确定正变换Transform将遥感影像Raster的四个顶点坐标变换到目标投影坐标系

DST中；确定四个顶点原有顺序，将重投影后的四个坐标点连结成四边形，并获取该四边形

的外接矩形OutRect。本技术实现遥感影像重投影的精度，可根据需要选择重投影方法，不仅采用查找表的方式进行遥感影像重投影，减少了重投影耗时的计算过程，保证了重投影效

率，另外只需生成一次查找表，不会明显增加数据量。

技术要求

1.一种基于查找表的遥感影像快速重投影方法，其特征在于，包括以下步骤：

预先生成四个空查找表LUX1、LUY1、LUX2和LUY2，并记录每个像素重投影后坐标信息，其中，每个查找表的列数与行数分别为RasterXSize和RasterYSize；

确定正变换Transform将遥感影像Raster的四个顶点坐标变换到目标投影坐标系DST中，表示为：

XREvertex,YREvertex＝Transform(XEvertex,YEvertex)，

其中，XEvertex,YEvertex为顶点坐标，XREvertex,YREvertex为顶点坐标在目标投影坐标系

DST中的坐标；

确定四个顶点原有顺序，将重投影后的四个坐标点连结成四边形，并获取该四边形的外接矩形OutRect；

扩展外接矩形OutRect，获取扩展后的矩形ExtendRect，生成空重投影影像Grid，表示为：

对获取的Grid中像素按照行列号进行遍历，确定每个像素的中心点行列号，将行列号变换成目标投影坐标系DST中坐标，表示为：

XE*＝XE+(COL+0.5)*Xres，YE*＝YE―(ROW+0.5)*YRes，

其中，XE*和YE*为Grid中第COL列第ROW行像素的中心点在目标投影坐标系DST中的坐标；

将获取的像素中心点DST坐标变换到SRC中，表示为：

XRE*,YRE*＝Transform_Inv(XE*,YE*)，

其中，XRE*,YRE*为目标投影坐标系DST中坐标[XE*,YE*]在源投影坐标系SRC中对应的坐标；

将获取的坐标点[XRE*,YRE*]转换成原始卫星遥感影像Raster的图像坐标，表示为：

COL*＝(XRE*―XRE)/Xres，ROW*＝(YRE―YRE*)/YRes，

其中，COL*和ROW*为坐标点[XRE*，YRE*]在遥感影像Raster中的列号与行号；提取Grid中OutRect范围内所有空像素，并将这些像素行列号写入查找表LUX2和LUY2中；

获取一景遥感影像Image，生成空重投影影像RepImage，从查找表获取非空像素处所记录的行列号，将Image中对应于非空像素位置处的像素值写入RepImage中对应于非空像素所记录行列位置处。

2.根据权利要求1所述的基于查找表的遥感影像快速重投影方法，其特征在于，所述查找

表其一景遥感影像Raster的宽高分别为RasterXSize和RasterYSize，分辨率为XRes和YRes，左上角坐标为[XRE,YRE]。

3.根据权利要求1所述的基于查找表的遥感影像快速重投影方法，其特征在于，所述确定

正变换Transform将遥感影像Raster的四个顶点坐标变换到目标投影坐标系DST中，包括以下步骤：

标定遥感影像Raster的投影坐标系为源投影坐标系SRC；标定目标投影坐标系为DST；

标定其SRC坐标系坐标转换为DST坐标系坐标的变换为正变换Transform；标定其DST坐标系坐标转换为SRC坐标系坐标的变换称为逆变换Transform_Inv。

4.根据权利要求2所述的基于查找表的遥感影像快速重投影方法，其特征在于，所述生成

空重投影影像Grid，进一步包括以下步骤：

确定ExtendRect的宽与高分别能被XRes与YRes整除；以XRes、YRes为分辨率，以ExtendRect范围为空间范围；标定Grid左上角投影坐标为[XE,YE]。

5.根据权利要求1所述的基于查找表的遥感影像快速重投影方法，其特征在于，所述提取Grid中OutRect范围内所有空像素，包括以下步骤：

判断其COL*和ROW*为坐标点[XRE*，YRE*]在遥感影像Raster中的列号与行号是否落入遥感影像Raster范围内；

其落入在Raster范围内，对行列号向下取整数得到整型行列号R和C；将获取的COL写入LUX1中第C列第R行，将ROW写入LUXY1中第C列第R行。

技术说明书

一种基于查找表的遥感影像快速重投影方法技术领域

本技术涉及遥感影像处理技术领域，具体来说，涉及一种基于查找表的遥感影像快速重投影方法。背景技术

随着卫星遥感技术的成熟、商业化的推进，卫星遥感得到快速发展，在越来越多的领域得到了应用。遥感影像预处理是卫星遥感应用的前提，而遥感影像重投影则是其中一个重要环节。遥感影像所采用的坐标系种类繁多，遥感影像与使用者需要的坐标系也未必相同，因此，为了在所需坐标系下对卫星遥感影像进行处理与分析，通常需要对遥感影像进行重投影。

遥感影像重投影计算过程极度耗时，为了提高遥感影像重投影效率。因此亟需一种基于查找表的遥感影像快速重投影方法。针对相关技术中的问题，目前尚未提出有效的解决方案。技术内容

针对相关技术中的问题，本技术提出一种基于查找表的遥感影像快速重投影方法，通过具有同一空间范围和空间分辨率的遥感影像配置共同的查找表，该查找表记录了遥感影像空间内所有像素所对应的重投影后的坐标信息。当需要对新的遥感影像进行重投影操作时，只需从查找表中对应像素处读取所记录的坐标信息，即可以将该像素值写入目标投影坐标系下的重投影的影像中。以克服现有相关技术所存在的上述技术问题。本技术的技术方案是这样实现的：

一种基于查找表的遥感影像快速重投影方法，包括以下步骤：

步骤S1，预先生成四个空查找表LUX1、LUY1、LUX2和LUY2，并记录每个像素重投影后坐标信息，其中，每个查找表的列数与行数分别为RasterXSize和RasterYSize；步骤S2，确定正变换Transform将遥感影像Raster的四个顶点坐标变换到目标投影坐标系

DST中，表示为：

XREvertex,YREvertex＝Transform(XEvertex,YEvertex)，

其中，XEvertex,YEvertex为顶点坐标，XREvertex,YREvertex为顶点坐标在目标投影坐标系

DST中的坐标；

步骤S3，确定四个顶点原有顺序，将重投影后的四个坐标点连结成四边形，并获取该四边形的外接矩形OutRect；

步骤S4，扩展外接矩形OutRect，获取扩展后的矩形ExtendRect，生成空重投影影像Grid，表示为：