[發明專利]一種基于Google Earth Engine云平臺的潮間帶灘涂快速遙感提取方法

在審

• 申請號: 201910494413.2
• 申請日: 2019-06-06
• 公開（公告）號: CN110296690A
• 公開（公告）日: 2019-10-01
• 發明（設計）人: 賈明明;王宗明;任春穎;毛德華
• 申請（專利權）人: 中國科學院東北地理與農業生態研究所
• 主分類號: G01C11/04
• 分類號: G01C11/04;G06K9/00;G06K9/32
• 代理公司: 哈爾濱華夏松花江知識產權代理有限公司 23213
• 代理人: 孟憲會

1.一種基于Google Earth Engine云平臺的潮間帶灘涂快速遙感提取方法，其特征在于，它包括如下步驟：

步驟一：利用Google Earth Engine云平臺，篩選研究區一年內云量較少的Sentinel-2MSI影像，形成一個影像集；

步驟二：計算影像集中每一個影像的水體指數NDWI，得到空間上的NDWI圖像，并將NDWI圖像當做一個波段，合并到原始影像中，NDWI計算公式如下；

NDWI＝(P_Green-P_NIR)/(P_Green+P_NIR)

式中，P_Green和P_NIR代表綠波段和近紅外波段的反射率；

步驟三：利用Google Earth Engine中qualityMosaic函數選取影像集中NDWI值最大的像元，形成最大NDWI像元組成的圖像；

步驟四：設置NDWI閾值，選出一年中至少有一次水體覆蓋的范圍，利用mask函數去掉沒有水體的區域；

步驟五：根據灘涂在Sentinel-2 MSI影像的8A波段反射率大于水體反射率的原理，利用Google Earth Engine中qualityMosaic函數選取影像集中8A波段值最大的像元，形成最大8A反射率的像元組成的圖像；

步驟六：設置8A閾值，選出最大8A反射率圖像中的非水面，利用mask函數去掉水面的區域；

步驟七：步驟四等到的是高潮能夠覆蓋的所有水面區域，步驟七得到低潮時的非水面區域，對這兩個區域進行空間疊加，就得到了潮間帶灘涂的分布圖。

2.根據權利要求1所述的一種基于Google Earth Engine云平臺的潮間帶灘涂快速遙感提取方法，其特征在于，其特征在于，在步驟一中，所述一年內云量較少的Sentinel-2MSI影像是指一年內云量少于20％的影像。

3.根據權利要求2所述的一種基于Google Earth Engine云平臺的潮間帶灘涂快速遙感提取方法，其特征在于，在步驟三中，qualityMosaic函數是一種基于像素的影像鑲嵌方法，用于將一系列時間序列上的像素中的值最大的像素選取出來組成一幅新的影像。

4.根據權利要求1、2或3所述的一種基于Google Earth Engine云平臺的潮間帶灘涂快速遙感提取方法，在步驟四中，mask函數可直接將閾值在一定范圍內的像素進行掩膜，使這部分像素區域被擦除，所述一定范圍是指NDWI小于或等于0.1。

5.根據權利要求4所述的一種基于Google Earth Engine云平臺的潮間帶灘涂快速遙感提取方法，其特征在于，在步驟五中，利用Google Earth Engine中qualityMosaic函數選取影像集中8A波段值最大的像元的過程中將8A波段設置為qualityMosaic的基準圖層。

6.根據權利要求5所述的一種基于Google Earth Engine云平臺的潮間帶灘涂快速遙感提取方法，其特征在于，在步驟六中，8A閾值的范圍為大于或等于1500。

7.根據權利要求6所述的一種基于Google Earth Engine云平臺的潮間帶灘涂快速遙感提取方法，其特征在于，在步驟六中，設置8A閾值，選出最大8A反射率圖像中的非水面，利用mask函數去掉非水面的區域的過程中，利用mask函數對最大8A反射率圖像中的每個像素進行掩膜，將像素值小于閾值的所有像素擦除。

8.根據權利要求7所述的一種基于Google Earth Engine云平臺的潮間帶灘涂快速遙感提取方法，其特征在于，在步驟七中，兩個區域進行空間疊加后，進行空間相交分析，得到潮間帶灘涂的分布圖。