[发明专利]一种适用于多投影遥感影像的海岸线数据快速自动叠加方法

摘要

本发明公开了一种适用于多投影遥感影像的海岸线数据快速自动叠加方法。包括如下步骤：1)读取遥感图像中5个参数；2)计算遥感影像范围矩阵的左上角和右下角坐标；3)判断海岸线文件和遥感图像文件的投影坐标系是否一致；4)如果投影坐标系不一致，则需要将步骤2)得到矩阵坐标进行投影转换；5)对步骤4)得到的范围矩阵进行扩充，对扩充后的范围建立内存格网索引，并获取每个格网中的海岸线对象；6)如果海岸线文件和遥感图像文件的投影坐标系不一致，则将步骤5)中获取到的对象也进行投影转换；7)从内存索引中实时获取显示对象，并将获取到的对象进行绘制。本发明能够满足以最小的内存耗费满足遥感影像范围内的海岸线快速叠加显示。

主权项

1.一种适用于多投影遥感影像的海岸线数据快速自动叠加方法，其特征在于，包括如下步骤：1)读取遥感图像中的投影坐标系、左上角坐标、分辨率、遥感图像宽度、遥感图像高度5个参数；2)根据步骤1)得到的参数计算遥感影像范围矩阵的左上角和右下角坐标；3)对步骤2)得到的范围矩阵进行扩充，使其宽度、高度均为1000的整数倍，接着对扩充后的范围建立内存格网索引；4)判断海岸线文件和遥感图像文件的投影坐标系是否一致；5)如果海岸线文件和遥感图像文件的投影坐标系不一致，则需要将步骤3)得到格网矩阵坐标进行投影转换，得到与海岸线坐标一致的矩阵坐标，如果一致，则进行步骤6)；6)利用SHAPELIB获取每个格网中的海岸线对象；7)如果海岸线文件和遥感图像文件的投影坐标系不一致，则将步骤5)中获取到的对象也进行投影转换，使对象的坐标与遥感图像文件中的坐标系一致，否则进行步骤8)；8)根据遥感影像当前显示大小，从内存索引中实时获取显示对象，并将获取到的对象进行绘制；所述步骤2)中计算公式为：x_min＝x_lt‑width*x_res  (1)y_min＝y_lt‑height*y_res  (2)x_max＝x_lt+width*x_res  (3)y_max＝y_lt+height*y_res  (4)其中，(x_lt，y_lt)是遥感图像左上角的坐标，width是遥感图像的宽度，height代表遥感图像的高度，x_res是遥感图像的横向分辨率，y_res是遥感图像的纵向分辨率；(x_min，y_min)代表当前显示范围的左上角坐标，(x_max，y_max)代表当前显示范围的右下角坐标；所述的步骤3)中内存格网索引为一个包含数据对象、格网范围、是否有数据的标识的一系列对象，其具体建立方式为：将显示范围矩阵划分为多个1000x1000的格网，并利用SHAPELIB开源库获取每个格网中的海岸线对象，并在海岸线对象填充的过程中标识该格网是否包含海岸线；所述步骤4)中的判断遥感图像和海岸线文件的投影是否一致，其判断方法为：abs(f_{false_easting}‑f'_{false_easting})＜0.0001  (5)abs(f_{false_northing}‑f'_{false_northing})＜0.0001  (6)abs(f_{scale_factor}‑f'_{scale_factor})＜0.0001  (7)abs(f_{central_meridian}‑f'_{central_meridian})＜0.0001  (8)abs(f_{latitude_origin}‑f'_{latitude_origin})＜0.0001  (9)其中，abs是绝对值函数，f_{false_easting}和f'_{false_easting}分别代表遥感图像投影坐标系和海岸线投影坐标的东伪偏移，f_{false_northing}和f'_{false_northing}分别代表遥感图像投影坐标系和海岸线投影坐标的北伪偏移，f_{scale_factor}和f'_{scale_factor}分别代表遥感图像投影坐标系和海岸线投影坐标的比例因子，f_{central_meridian}和f'_{central_meridian}分别代表遥感图像投影坐标系和海岸线投影坐标的中央经线，f_{latitude_origin}和f'_{latitude_origin}分别代表遥感图像投影坐标系和海岸线投影坐标的起始原点；当遥感图像坐标系和海岸线坐标系均满足式5)至式9)时，则认为两者坐标系一致，否则判定为不一致；所述的步骤5)的实现方式为：建立坐标原始坐标系统、建立目标坐标系统，并利用ProjNet开源库，实现坐标从原始坐标系统到目标坐标系统的转换；所述的步骤8)中实时获取显示对象的方法为：将当前遥感影像显示范围矩阵与格网求交，绘制与其有交集的格网数据；获取遥感影像当前显示范围的公式如下：x_min＝x_lt‑width*x_res*ratio  (11)y_min＝y_lt‑height*y_res*ratio  (12)x_max＝x_lt+width*x_res*ratio  (13)y_max＝y_lt+height*y_res*ratio  (14)其中，(x_lt，y_lt)是遥感图像左上角的坐标，width是遥感图像的宽度，height代表遥感图像的高度，x_res是遥感图像的横向分辨率，y_res是遥感图像的纵向分辨率，ratio代表当前遥感图像的缩放因子；(x_min，y_min)代表当前显示范围的左上角坐标，(x_max，y_max)代表当前显示范围的右下角坐标。