[發明專利]一種基于監督分類和DEM的露天開采區邊界提取方法

在審

• 申請號: 202010995020.2
• 申請日: 2020-09-21
• 公開（公告）號: CN112232135A
• 公開（公告）日: 2021-01-15
• 發明（設計）人: 包勤躍;汪宏;趙守恒;汪兵兵
• 申請（專利權）人: 銅陵有色金屬集團股份有限公司
• 主分類號: G06K9/00
• 分類號: G06K9/00;G06K9/46;G06K9/62
• 代理公司: 安徽匯樸律師事務所 34116
• 代理人: 朱克儉

1.一種基于監督分類和DEM的露天開采區邊界提取方法，其特征在于，該方法包括以下步驟：

S1：選取指定區域開采后對應的質量較優的Landsat 8OLI影像，并對其進行預處理，得到Image1；

S2：以與步驟1選取的Landsat 8OLI影像同一時期的高分辨率Google Earth影像為參照建立訓練樣本，并對Image1進行監督分類，獲得含分類信息的Image2；

S3：將Image2中分類信息為露天開采區的柵格轉為矢量，得到矢量數據Vector1；

S4：在ArcGIS中對矢量數據Vector1中的矢量要素進行空洞和碎斑處理，得到只含露天開采區邊界1信息的矢量數據A；

S5：以ArcGIS為平臺，對指定區域開采前后的兩期DEM分別作投影處理，得到影像DEM1和影像DEM2；

S6：對DEM1和DEM2進行配準并作相減運算，得到含高程差值信息的影像Elevation1；

S7：對影像Elevation1作二值化處理過濾掉低于閾值的高程差信息，獲得含露天開采區高程差信息的影像Elevation2；

S8：提取影像Elevation2中有露天開采區高程差信息的區域對應的柵格并轉為矢量，得矢量數據Vector2；

S9：對矢量數據Vector2進行空洞和碎斑處理，得到含露天開采區邊界2信息的矢量數據B；

S10：在ArcGIS中對矢量數據A和矢量數據B作相交運算，得到含真實露天開采區邊界信息的矢量數據C。

2.根據權利要求1所述的方法，其特征在于，S1中所述質量較優的Landsat 8OLI影像，是以云量低于10％且指定區域內無云為條件篩選獲得的。

3.根據權利要求1所述的方法，其特征在于，所述S1中的預處理為輻射定標和大氣校正。

4.根據權利要求1所述的方法，其特征在于，所述S2中的監督分類為支持向量機法，且分類后進行聚類處理。

5.根據權利要求1所述的方法，其特征在于，所述S3和S8中柵格轉矢量的邊界均為原邊界，未進行光滑處理。

6.根據權利要求1所述的方法，其特征在于，所述S4和S9中的碎斑處理，具體為篩選面積小于0.02km²的多邊形并刪除。

7.根據權利要求1所述的方法，其特征在于，S5所述的投影處理為WGS1984地理坐標系、UTM投影坐標系。

8.根據權利要求1所述的方法，其特征在于，S6所述的配準是在ArcGIS中以開采前的DEM1為基準對開采后的DEM2作配準。

9.根據權利要求1所述的方法，其特征在于，S7所述的二值化閾值為5m。