首先感谢ENVI-IDL技术殿堂在ENVI平台下的指导,原文链接,现将其在GEE平台实现。不同的是,该博文称使用原始DN值计算NDVI对反演影响不大,但通过实践发现还是有差异,所以代码采用了经过大气校正的地表反射率产品进行NDVI的计算。
注意的地方:代码中44和45行需要自行替换,也可以通过filter过滤功能进行查找得到。另一个需要注意的是原博文中提到的三个变量参数还是需要自己去NASA获取,获取地址,中心经纬度可以通过自己在地图上添加geometry得到。
体谅的地方:接触GEE才开始接触Javascript,所以格式可能不规范,尽量多学习。GEE里面好像不能通过1.divide(image)这样来对图像求倒数,所以很笨的用image.divide(image.multiply(image))这样的方式来实现。
代码基于辐射传输方程实现(如果后期有时间学习分裂窗算法,将再次更新):
/*
* Calculate Land Surface Temperature
* Based on Radiative Transfer Equation
* from Landsat-8 Surface Reflectance Product
* and Landsat-8 Raw Scenes