首页 理论教育大沽河流域水文要素监测体系建设流程

大沽河流域水文要素监测体系建设流程

【摘要】:使用Matlab编程,将2个时相8d合成的地表温度数据合成为1个时相16d合成的地表温度数据。利用上一步骤中提取的数据,在Excel中对NDVI和最大和最小陆地表面温度进行线性拟合,获得干边和湿边方程的系数a1、b1、a2和b2。使用Matlab编程,根据式计算图像上每个像元的温度植被干旱指数值,获取流域TVDI的分布图,根据TVDI等级划分形成大沽河流域土壤湿度分布图。图5.4显示了部分Matlab程序文件列表。

图像处理过程由ENVI图像处理软件和Matlab软件完成,处理过程如下:

(1)首先利用MODIS数据产品投影变换软件MRT 进行投影变换,采用WGS-1984椭球体,UTM 投影。

(2)在ENVI中使用大沽河流域界线矢量图形剪切出流域的地表温度图像和植被指数图像(NDVI)。

(3)使用Matlab编程,将2个时相8d合成的地表温度数据合成为1个时相16d合成的地表温度数据(单位为℃)。合成时,如果某一像元2个时相都有值,则取其平均作为合成后的像元值;如果只有1个时相有值,就使用这个时相的值作为合成后的像元值。

(4)使用Matlab编程,利用16d合成的地表温度数据和植被指数数据,提取某一NDVI对应的所有地表温度中的最大值和最小值,将不同NDVI下的最大和最小陆地表面温度保存于Excel文件中。

(5)利用上一步骤中提取的数据,在Excel中对NDVI和最大和最小陆地表面温度进行线性拟合,获得干边和湿边方程的系数a1、b1、a2和b2

(6)使用Matlab编程,根据式(5.11)计算图像上每个像元的温度植被干旱指数(TVDI)值,获取流域TVDI的分布图,根据TVDI等级划分形成大沽河流域土壤湿度分布图。

图5.4显示了部分Matlab程序文件列表。

图5.4 研究中编写的部分Matlab程序文件