定量遥感：地表温度反演.docVIP

遥感数字影像处理 作品名称：黄河三角洲地表温度反演 姓名+学号： 小组成绩： 一、概述 作业背景： 地表温度是很多环境模型的一个重要参数，在大气与地表的能量与物质交换，天气预报，全球洋流循环，气候变化等研究领域有重要的应用。利用热红外遥感可以得到大范围的地表温度面状信息，与传统的地表温度测量方式相比，具有快速、便捷、测量范围大、信息连续等特点，因此利用热红外遥感数据反演地面温度得到了广泛的应用 作业意义： 黄河三角洲是黄河携带大量泥沙在渤海凹陷处沉积形成的冲积平原，位处黄河入海口处的黄河三角洲自然保护区正是以保护河口湿地生态系统和珍稀、濒危鸟类为主的湿地类型保护区。以利津为顶点，北到徒骇河口，南到小清河口，呈扇状三角形，面积5,450平方公里。地面平坦，在海拔10公尺以下向东撒开的扇状地形，海拔高程低于15米，面积达5450平方公里。三角洲属，温带季风性气候。四季分明，光照充足，区内自然资源丰富。 黄河口湿地生态旅游区占地23万亩，都处在黄河三角洲之内，地貌以芦苇沼泽，湿地为主，其次为河口滩地，带翅碱蓬盐滩湿地，灌丛疏林湿地以及人工槐林湿地等。集自然景观与人文景观为一体，既有沧海桑田的神奇与壮阔，又有黄龙入海的壮观和长河落日的静美，是人们休闲、度假、观光科普的最佳场所遥感传感器的数字量化输出值DN与其所对应视场中辐射亮度值之间的定量关系 导入后截图： 加载后数据截图 4.2、数据辐射定标： 处理后结果： 4.3、裁剪数据： 定标后的数据需要进行影像裁剪处理，选择需要分析的区域。在谷歌地图上确定研究范围的经纬度信息，在其他软件上将平面的经纬度信息转至相同投影面上的数据，作为裁剪的标准。准备好shipfile格式数据： 选择辐射定标好的影像 Shipfile格式数据转换为ROI格式： 裁剪： 裁剪完毕后对比： 对于红外波段数据的裁剪方法与上述方法相同。 4.4、可见光和近红外波段大气校正 采用大气校正来消除遥感图像中由大气散射引起的辐射误差的处理过程。 4.5、地表比辐射率值计算 NDVI指数计算： 利用TM3、4波段的象元DN值利用公式： NDVI=（NIR - R）/（NIR + R） 处理后结果： 由上统计可以看出，NDVI计算结果最小值为-2.2，最大值为1.4. （2）植被覆盖度 计算植被覆盖度采用的是混合像元分解法，将整幅影像的地类大致分为水体、植被和建筑。 利用Band Math计，根据上步操作中计算出来的NDVI指数 输入公式：(b1 lt 0.0)*0+(b1 ge 0 and b1 le 1.4)*((b1 - 0.0)/(1.4-0.0))。b1：表示获取的NDVI值。 计算结果如下图所示： 处理结果如下： （2）地表比辐射率： 根据前人研究，将遥感影像分为水体、城镇和自然表面3种类型。本次实习采用：水体像元的比辐射率赋值为0.995，自然表面和城镇像元的比辐射吕估算则根据下公式中计算 利用Band Math计算 输入公式：(b1 le 0)*0.995+(b1 gt 0 and b1 lt 0.7)*(0.9589+0.086*b2-0.0671*b2^2)+(b1 ge 0.7)*(0.9625+0.0614*b2-0.0461*b2^2)。b1：表示获取的NDVI值，b2：植被覆盖度值。 处理后结果： 4.6、大气参数获取 卫星传感器接收到的热红外辐射亮度值由三部分：大气向上辐射亮度，以及地面的真实辐射亮度经过大气层之后到达卫星传感器的能量。地面的真实辐射亮度为同温度黑体的辐射亮度值与地物发射率的乘积。 NASA官网（）中获取大气参数获取，在上网址中输入成像时间以及中心经纬度，查找下式中需要的参数。 根据上图看出，大气在热红外波段的透过率为0.68，大气向上辐射亮度为 2.55W/(m2·sr·μm)，大气向下辐射亮度为3.97W/(m2·sr·μm) 利用Band Math计算 输入公式：(b2-2.55-0.68*(1-b1)*3.97)/(0.68*b1)。b1：地表比辐射率值，b2：热红外波段辐射定标志。 获取的黑体在热红外波段的辐射亮度结果如下： 4.7地表温度反演 在获取热红外波段辐射亮度值以后根据普朗克公式的反函数，求得地表真实温度T。对于TM，K1?=607.76W/(m2·sr·μm)，K2?=1260.56K。 利用ENVI中的Band Math，在公式栏中输入：(1260.56)/alog (607.76/b1+1) 4.8反演结果分析 利用ENVI软件，采用单波段彩色变换的方法对地表真实温度的灰度图进行密度分割得到地表真实温度的分