建议先看我的另一篇 Landsat计算TVDI进行干旱监测(二)。该篇文章由于劈窗算法并未完善好,只适用于干旱区,应用精度不佳。
影像下载
预先下载了研究区13景的Landsat8和9的影像,该影像有10个波段,其中1-9为可见光波段(30m),10-11为热红外波段(100m)。在存储时,Landsat的每个波段为TIFF形式。
OLI波段需要进行辐射定标和大气校正,两个热红外单波段要进行辐射定标
预处理步骤可见(Landsat8的辐射定标与大气校正)
做完预处理后存为TIFF影像
NDVI计算和热红外单波段提取
我图中的公式后面少了float