[發明專利]一種湖泊面積提取方法

有效

• 申請號: 201910584480.3
• 申請日: 2019-07-01
• 公開（公告）號: CN110276797B
• 公開（公告）日: 2022-02-11
• 發明（設計）人: 任黎;周悅;陰帥妮;徐偉;徐健
• 申請（專利權）人: 河海大學
• 主分類號: G06T7/62
• 分類號: G06T7/62
• 代理公司: 南京縱橫知識產權代理有限公司 32224
• 代理人: 范青青;董建林

1.一種湖泊面積提取方法，其特征在于，所述方法包括如下步驟：

采用FROM-GLC數據集獲取湖泊區TIFF遙感影像，對湖泊區TIFF遙感影像進行轉換，獲取湖泊面要素文件；

利用湖泊面要素文件獲取湖泊區raster文件；

根據湖泊區raster文件導出湖泊區mask文件；

根據湖泊區mask文件，批量處理長時序的湖泊區TIFF遙感影像，批量提取計算湖泊面積；

所述湖泊面要素文件的獲取方法包括：

在ARCGIS10.2軟件平臺下，根據所述湖泊區TIFF遙感影像制作湖泊線要素shape文件，然后利用ArcToolbox工具箱中的data management tools，選擇features下的feature topolygon，將所述湖泊線要素shape文件轉化為湖泊面要素shape文件；

所述湖泊區raster文件的獲取方法包括：

用湖泊面要素文件對原所述湖泊區TIFF遙感影像進行初步裁剪，獲取初步裁剪的湖泊區raster文件；

根據湖泊區TIFF遙感影像及對應湖泊區域水域特點，通過閾值設定修正掩模文件，使初步裁剪的湖泊區raster文件逐步調整至湖泊區TIFF遙感影像的大范圍之內，獲取最終裁剪好的湖泊區raster文件；

修正掩模文件的方法包括：

利用ARCGIS10.2軟件ArcToolbox工具箱中的spatial analyst tools，選擇mapalgebra，然后選擇raster calculator，把所有大于等于0的區域設置為300，大于等于0的區域設為湖泊區，再次輸出處理的湖泊區raster文件，通過若干次合并柵格文件，獲得在湖泊區TIFF遙感影像的大范圍之內的湖泊區raster文件；

最終裁剪好的湖泊區raster文件的獲取方法包括：

針對位于湖泊區TIFF遙感影像影響的大范圍之內的湖泊區raster文件，利用ARCGIS10.2軟件ArcToolBox工具箱的data management，選擇raster，然后選擇rasterdataset下的mosaic，將湖泊范圍內的區域賦值為1，其它區域賦值為-9999，輸出湖泊區raster文件，即為最終剪裁好的湖泊raster文件。

2.根據權利要求1所述的湖泊面積提取方法，其特征在于，初步裁剪的湖泊區raster文件的獲取方法包括：

利用ARCGIS10.2軟件ArcToolbox工具箱中的spatial analyst tools，選擇extraction，然后選擇extract by mask生成初步裁剪好的湖泊raster文件。

3.根據權利要求1所述的湖泊面積提取方法，其特征在于，所述FROM-GLC數據集的重訪周期為一天。

4.根據權利要求1所述的湖泊面積提取方法，其特征在于，所述湖泊區mask文件的導出方法包括：

將最終裁剪好的湖泊區raster文件導出為txt文件；

將所述txt文件導入MATLAB軟件，保存輸出湖泊區mask文件。

5.根據權利要求4所述的湖泊面積提取方法，其特征在于，所述txt文件的導出方法包括：

利用ARCGIS10.2軟件ArcToolBox工具箱的conversion tool，選擇from raster下的raster to ascii生成所述txt文件。

6.根據權利要求1所述的湖泊面積提取方法，其特征在于，所述湖泊面積的計算方法包括：

用MATLAB軟件編輯tiffread程序，根據所述湖泊區mask文件批量裁剪長時序的湖泊區TIFF遙感影像，得到每日湖泊范圍的水像元數，然后通過公式（1）計算出長時序每日的湖泊面積：

（1）

式（1）中S為每日的湖面面積，為水像元數。