[发明专利]稀土矿的地表温度空间降尺度方法、装置、设备及介质

权利要求书

1.一种稀土矿的地表温度空间降尺度方法，其特征在于，包括：

获取指定时段内的目的稀土矿区域的原始遥感影像及对应的高分影像；

对所述原始遥感影像进行图像融合，得到预设空间分辨率的多光谱数据，并对融合后的多光谱数据进行图像裁剪；

根据所述目的稀土矿区域的实际地表情况及所述高分影像确定所述目的稀土矿区域的主体端元；

根据所述多光谱数据、所述主体端元及全约束的线性光谱混合模型确定所述目的稀土矿区域的地表组分丰度值；

对所述原始遥感影像依次进行辐射定标、大气校正、图像裁剪，得到预处理后的遥感影像；

根据所述预处理后的遥感影像确定典型端元地表温度；

利用所述地表组分丰度值及所述典型端元地表温度，实现所述目的稀土矿区域的地表温度空间降尺度；

所述主体端元包括下列中的至少一项：植被、裸土、尾砂；

所述地表组分丰度值的表达式为：

其中，R_pi为i波段上p像元的光谱反射率；N为像元包含的基本组分数目；f_jp为p像元的j端元组分的丰度值，且，0≤f_jp≤1，ε_pi为i波段上p像元的剩余残差；

所述利用所述地表组分丰度值及典型端元地表温度，实现所述目的稀土矿区域的地表温度空间降尺度，包括：

根据所述地表组分丰度值及典型端元地表温度，利用普朗克黑体辐射定律得到所述目的稀土矿区域的高分辨率的地表温度；

所述目的稀土矿区域的高分辨率的地表温度的表达式为：

R(T_L)＝ε₁₅R_B(T_L)

T_L＝K₂/ln[K₁/R_B(T_L)+1]

其中，R(T_L)为辐射亮度值；ε₁₅为15m空间分辨率的比辐射率；R_B(T_L)为同温度的黑体辐射亮度值；K₁、K₂为Landsat卫星发射前的预设数值，K₁＝774.89W/(m²·μm·sr)，K₂＝1321.08K；

其中，

其中，ε₁₅为15m空间分辨率的比辐射率；ε_i为地表组分i的比辐射率；f_i为地表组分i的丰度值。

2.根据权利要求1所述的稀土矿的地表温度空间降尺度方法，其特征在于，所述原始遥感影像为陆地卫星Landsat 8卫星影像，所述高分影像为谷歌地球Google Earth高分影像。

3.根据权利要求1所述的稀土矿的地表温度空间降尺度方法，其特征在于，所述对所述原始遥感影像进行图像融合，得到预设空间分辨率的多光谱数据，包括：

基于格拉姆-施密特正交化Gram-Schmidt算法，利用ENVI的图像融合模块对所述原始遥感影像进行图像融合，得到15m空间分辨率的多光谱数据。

4.根据权利要求1所述的稀土矿的地表温度空间降尺度方法，其特征在于，所述对所述原始遥感影像依次进行辐射定标、大气校正、图像裁剪，得到预处理后的遥感影像，包括：

利用遥感图像处理平台ENVI的辐射定标Radiometric Calibration工具对所述原始遥感影像进行辐射定标，得到第一预处理遥感影像；

利用所述ENVI的大气校正模块FLAASH对所述第一预处理遥感影像进行大气校正，得到第二预处理遥感影像；

根据所述目的稀土矿区域的矢量边界对所述第二预处理遥感影像进行图像裁剪，得到预处理后的遥感影像。