[发明专利]一种利用GDAL和SHAPELIB在遥感图像上快速叠加显示海岸线方法

摘要

本发明公开了一种利用GDAL和SHAPELIB在遥感图像上快速叠加显示海岸线方法。包括如下步骤：1)利用SHAPELIB开源库对海岸线数据建立Q树索引，同时获取海岸线数据的坐标系；2)利用GDAL读取遥感图像中的5个参数；3)计算显示范围；4)判断海岸线文件和遥感图像文件的投影坐标系是否一致；5)如果投影坐标系不一致，则需要将显示范围进行坐标转换；6)利用得到的显示范围坐标，通过SHAPELIB开源库得到显示对象；7)如果海岸线文件和遥感图像文件的投影坐标系不一致，则将对象进行投影转换。本发明通过GDAL和SHAPELIB的联合应用，实现了在遥感图像上快速叠加显示海岸线，具有处理速度快、转换时间短的优点，对遥感图像处理与分析领域有着积极意义。

主权项

1.一种利用GDAL和SHAPELIB在遥感图像上快速叠加显示海岸线方法，其特征在于，包括如下步骤：1)利用SHAPELIB开源库对海岸线数据建立Q树索引，同时获取海岸线数据的坐标系；2)利用GDAL读取遥感图像中的投影坐标系、左上角坐标、分辨率、遥感图像宽度、遥感图像高度5个参数；3)根据步骤2)得到的参数计算显示范围；4)判断海岸线文件和遥感图像文件的投影坐标系是否一致；5)如果海岸线文件和遥感图像文件的投影坐标系不一致，则需要将步骤3)得到的显示范围进行坐标转换，如果一致，则进行步骤6)；6)利用得到的显示范围坐标，通过SHAPELIB开源库得到显示对象；7)如果海岸线文件和遥感图像文件的投影坐标系不一致，则将步骤5)中获取到的对象进行投影转换，否则进行步骤8)；8)将获取到的对象进行绘制；所述步骤3)中计算公式为：x_min＝x_lt‑width*x_res*ratio  (1)y_min＝y_lt‑height*y_res*ratio  (2)x_max＝x_lt+width*x_res*ratio  (3)y_max＝y_lt+height*y_res*ratio  (4)其中，(x_lt，y_lt)是遥感图像左上角的坐标，width是遥感图像的宽度，height代表遥感图像的高度，x_res是遥感图像的横向分辨率，y_res是遥感图像的纵向分辨率，ratio代表当前遥感图像的缩放因子；(x_min，y_min)代表当前显示范围的左上角坐标，(x_max，y_max)代表当前显示范围的右上角坐标；所述步骤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)时，则认为两者坐标系一致，否则判定为不一致；所述的步骤6)具体为：根据获取到的显示范围，利用SHAPELIB开源库，对显示对象进行提取。