1.实习目的
了解和熟悉植被指数/地表温度特征空间与土壤湿度的关系。
2.实习内容
利用给定植被区域的地表温度图像和植被指数图像，分析二者之间的关系，并能给出简单解释；
了解利用温度植被干旱指数反映土壤湿度的方法。
3.实习原理与步骤
3.1 植被区域的地表温度图像和植被指数图像关系分析
(1)利用“分类图像(60m)ok”数据，提取植被分布ROI区域；
具体步骤：
打开加载”分类图像(60m)ok“, "NDVI", "lst"文件，从图中可以看出，原始"分类图像(60m)ok"图像已经对地物进行了分类，其中类型7为植被覆盖区。
提取ROI区域：Toolbox-ROI-Band Threshold to ROI-选择“分类图像(60m)ok"。并且阈值参数只提取7。
加载选定的ROI区域，下图中其中红色部分为植被覆盖区域。
接下来进行掩膜提取ndvi图像中的植被覆盖区域：1）首先需要将刚才的ROI区域保存为”vegROI“。2）Toolbox-搜索”mask“-Build Mask-选择图像”分类图像(60m)ok“。在Mask Definition窗口下-Options-import ROI-并将新生成的mask保存输出，命名为”mask-7“
下图为生成的mask-7掩膜图像，红十字表示当前指定的目标区域。可以看出亮白色区域的DN值为1，黑色区域的DN值为0。即生成了一张植被覆盖区域的值为1，其他区域值为0的二值化图。
(2)绘制研究区域和植被区域的植被指数/地表温度分布散点图，分析散点图的特点，并给出解释。
具体步骤：在上方工具栏中找到Scatter Plot Tool工具（用于绘制散点图）。横坐标设置为NDVI, 纵坐标设置为地表温度lst。
然后再import ROI将之前生成的感兴趣区加载进去。这样既可以查看研究区域的植被指数/地表温度散点图，又能查看植被区域的散点图。由于软件显示问题，纵坐标的显示分别是”-10 0 10 20 30 40 50“，横坐标的显示分别是” 0.0 0.2 0.4 0.6“。
分析：从图像中可以看出植被区域的ndvi值基本都处在0.0到0.6之间，lst温度基本都在20以上，不高于45，且分布相对集中。与分类图中其他地物相比较，其ndvi值基本都在0以上，且随着ndvi值的提高，植被区域所占比例越高，至0.4以后上时基本全部是植被区域。这可以用于区分植被区域的参考标准。
3.2 温度植被干旱指数计算
(1)利用tvdi_main.sav工具，绘制研究区域的温度和NDVI的特征空间，并计算干湿边。
具体步骤：首先将文件“lst”中的摄氏温度转化为华氏温度。利用波段运算来处理。Band Math-新建公式“b1+273.15”, b1为lst中的数据。保存名“LST+273.15”.
打开扩展包目录下的tvdi_main.sav工具（前置操作为需要先将数据包中的该扩展包自行放置在“ENVI53\extensions”目录下），设置好对应的edvi和lst参数。从lst图像的数据中可以看出，温度基本分布在20-40℃之间，于是需要缩小有效LST范围，改为“290”“330”.然后即可选择输出路径。
(2)绘制研究区域的温度植被干旱指数分布图。
下图为生成的干湿边散点图和干湿边的系数拟合图像。
在envi中加载tvdi分布图发现边缘区域仍有数值，所以还需要再进行一次掩膜处理，筛选出研究区域，剔除边缘多余信息。
步骤：1)提取ROI区域：Toolbox-ROI-Band Threshold to ROI-选择“分类图像(60m)ok"。并且阈值参数选取为1-8。
2）Toolbox-搜索”mask“-Build Mask-选择图像”new_tvdi“。在Mask Definition窗口下-Options-import ROI-并将新生成的mask保存输出，命名为”mask_tvdi“。
3)Apply Mask。选择处理图像“new_tvdi”, 下面的“Select Band Mask”选择“mask_tvdi”, 默认参数保存命名为“pro_tvdi”。可以发现边缘区域的值均为0。
3.3制作干旱专题图
（1）在ArcGIS中绘制专题图。参考以下表格，根据温度植被干旱指数TVDI与干旱等级的对应关系，绘制植被区域干旱专题图。要求：干旱情况分为湿润、正常、微旱、干旱、重旱五个等级，栅格颜色由蓝色渐变为红色，可参考以下分类图。
（2）根据以下表格统计研究区域的土壤湿度情况：
