12-专题：实用技术.pptxVIP

技术专题：实用技术;Landsat8TIRS反演地表温度 基于CART的自动决策树分类 国产资源三号数据处理 水深反演和船只提取 图像数据类型转换 图像批量/分幅裁剪技术;1. Landsat8TIRS反演地表温度;热红外遥感（Infrared Remote Sensing）是指传感器工作波段限于红外波段范围之内的遥感。即利用星载或机载传感器收集、记录地物的热红外信息，并利用这种热红外信息来识别地物和反演地表参数如温度、湿度和热惯量等。目前有很多的卫星携带了热红外传感器，包括ASTER、AVHRR、MODIS、TM/ETM+/ TIRS等。 练习数据:\ 3-专题：基于像元二分模型的植被覆盖度反演\1-Landsat8 ;名词;目前，地表温度反演算法主要有以下三种 大气校正法（也称为辐射传输方程：Radiative Transfer Equation——RTE） 单通道算法 分裂窗算法 本实例是基于大气校正法，利用Landsat8 TIRS反演地表温度。 基本原理：首先估计大气对地表热辐射的影响, 然后把这部分大气影响从卫星传感器所观测到的热辐射总量中减去, 从而得到地表热辐射强度, 再把这一热辐射强度转化为相应的地表温度。;卫星传感器接收到的热红外辐射亮度值Lλ由三部分组成：大气向上辐射亮度L↑，地面的真实辐射亮度经过大气层之后到达卫星传感器的能量；大气向下辐射到达地面后反射的能量。卫星传感器接收到的热红外辐射亮度值Lλ的表达式可写为（辐射传输方程）： Lλ = [εB(TS) + (1-ε)L↓]τ + L↑ （1.1） 式中，ε为地表比辐射率，TS为地表真实温度(K)，B(TS)为黑体热辐射亮度，τ为大气在热红外波段的透过率。则温度为T的黑体在热红外波段的辐射亮度B(TS)为： B(TS) = [Lλ - L↑- τ(1-ε)L↓]/τε （1.2） Ts可以用普朗克公式的函数获取。 TS = K2/ln(K1/ B(TS)+ 1) （1.3） 对于TIRS Band10，K1= 774.89 W/(m2*μm*sr)，K2 = 1321.08K。;1. 处理流程;Band10 热红外数据辐射定标，得到B10辐射亮度图像。 工具：Radiometric Correction/Radiometric Calibration OLI辐射定标和大气校正为可选操作。;TIRS的Band10热红外波段与TM/ETM+6热红外波段具有近似的波谱范围，本例采用TM/ETM+6相同的地表比辐射率计算方法。使用Sobrino提出的NDVI阈值法计算地???比辐射率。 ε=0.004Pv+0.986 其中， Pv是植被覆盖度，用以下公式计算： Pv = [(NDVI- NDVISoil)／(NDVIVeg - NDVISoil)] NDVI为归一化植被指数，NDVISoil为完全是裸土或无植被覆盖区域的NDVI值，NDVIVeg则代表完全被植被所覆盖的像元的NDVI值。取经验值NDVIVeg = 0.70和NDVISoil = 0.05。;第一步：计算NDVI Spectral/Vegetation/NDVI 第二步：计算植被覆盖度 (b1 gt 0.7)*1+(b1 lt 0.05)*0+(b1 ge 0.05 and b1 le 0.7)*((b1-0.05)/(0.7-0.05)) b1：NDVI计算结果 第三步：计算地表比辐射率 0.004*b1+0.986 b1：植被覆盖度图像 ;在NASA公布的网站查询（），输入成影时间：2013-10-03 02:55和中心经纬度（Lat：40Lon：116，以及其他相应的参数，得到大气剖面信息为： 大气在热红外波段的透过率τ：0.90 大气向上辐射亮度L↑：0.75 W/(m2·sr·μm) 大气向下辐射亮辐射亮度L↓：1.29W/(m2·sr·μm) 提示：由于缺少地表相关参数（气压、温度、相对湿度等信息），得到的结果是基于模型计算的结果。;计算得到同温度下的黑体辐射亮度图像 (b2-0.75-0.9*(1-b1)*1.29)/(0.9*b1) b1：地表比辐射率图像 b2：Band10辐射亮度图像 计算得到地表温度图像（单位为℃） (1321.08)/alog(774.89/b1+1)-273 b