[发明专利]基于遥感数据的矿山和固废填埋场自动提取方法和装置

权利要求书

1.一种基于遥感数据的矿山和固废填埋场自动提取方法，其特征在于，所述方法包括：

S1：获取待检测区的亚米级遥感影像和热红外遥感影像；

S2：对所述亚米级遥感影像和热红外遥感影像进行预处理；

S3：根据归一化植被指数方法，从预处理后的亚米级遥感影像上提取非植被覆盖区；

S4：根据K-T变换方法，从预处理后的亚米级遥感影像上提取非建筑用地区；

S5：根据预处理后的热红外遥感影像计算地表温度，并从预处理后的热红外遥感影像上提取高温异常区；

S6：将高温异常区与非植被覆盖区、非建筑用地区进行空间相交分析，其交集为矿山和固废填埋场图斑；

S7：根据矿山和固废填埋场图斑的亚米级遥感影像的光谱特征对所述矿山和固废填埋场图斑进行分类，得到矿山区域和固废填埋场区域。

2.根据权利要求1所述的基于遥感数据的矿山和固废填埋场自动提取方法，其特征在于，所述S3包括：

S31：基于预处理后的亚米级遥感影像，计算每个像元的归一化植被指数NDVI；

其中，NDVI为归一化植被指数，ρ_NIR、ρ_RED分别为预处理后的亚米级遥感影像的近红外波段、红波段的遥感反射率；

S32：将NDVI＜0的像元提取为非植被覆盖区。

3.根据权利要求2所述的基于遥感数据的矿山和固废填埋场自动提取方法，其特征在于，所述S4包括：

S41：基于预处理后的亚米级遥感影像，根据K-T变换模型Y＝R^Tx+r进行主成分分析，提取每个像元的第一主成分分量Y1和第二主成分分量Y2；

其中，Y为变换后的主成分分量，Y＝{Y1、Y2......Yn}；R为K-T变换系数，R＝{R1、R2......Rn}，R^T为R的转置，x为预处理后的亚米级遥感影像各波段的灰度值，x＝{x1、x2......xn}，r为常数偏移量，r＝{r1、r2......rn}；

S42：根据每个像元的第一主成分分量Y1和第二主成分分量Y2计算每个像元的比值建筑指数RBI；

RBI＝Y₁/Y₂

S43：将RBI＜-0.05的像元提取为非建筑用地区。

4.根据权利要求3所述的基于遥感数据的矿山和固废填埋场自动提取方法，其特征在于，所述S5包括：

S51：根据预处理后的热红外遥感影像计算每个像元的地表温度T_s；

其中，C₁、C₂为Plank函数参数，λ为有效波长，B(T_s)＝[L_λ-L_up-τ*(1-ε)*L_down]/(ε*τ)，L_λ为预处理后的热红外遥感影像λ波段的辐射亮度，L_up、L_down、τ分别为大气上行辐射、下行辐射和大气透过率，ε为比辐射率；

S52：对每个像元的地表温度T_s通过如下公式进行正态分布拟合，并计算拟合优度p-value，得到p-value∈(0.75，1)时对应的标准差σ；

其中，σ为标准差，μ为T_s的均值；

S53：将T_s＞1.96σ的像元提取为高温异常区。