SEBAL模型具体步骤

需要下载的数据：

1. Landsat8数据：https://earthexplorer.usgs.gov/
2. 瞬时气象数据：https://xihe-energy.com/#climate
3. 土地利用数据：https://zenodo.org/record/5816591#.ZEIse_xBxD8
4. 大气透射率数据：https://atmcorr.gsfc.nasa.gov/

目前网上已有的模型都不完整，发邮件找作者下载了完整的数据包，完成之后可以跑出来，结果也比较符合我的研究区的实际情况。但是，我的研究区位于东南地区，近四十年里找不到全景符合云量要求的图，区域内也仅有四张符合的，这四张里面三张都在研究区上空存在密集的小块云层，会极大地影响地表温度的计算，最后只有一张正确的出图。因此，研究区降水较多的不建议使用该模型。所以最后我放弃了这个方法，但是这几个月摸索的时候把操作步骤也写出来了，不想浪费所以分享一下。还有一些具体的计算excel文件等有空了也整分享出来。

步骤如下

1 遥感数据预处理

1. arcgis导入dem—检查是否是土地利用数据的投影坐标（Albers_Conic_Equal_Area）
2. LandsatMTL文件更改：

GROUP = L1_METADATA_FILE

END_GROUP = L1_METADATA_FILE

1. 辐射定标（红外thermal—spectral subset—10.9）
2. 辐射亮度计算【1321.08/alog(774.89/b1+1)】，b1为红外辐射定标
3. 计算大气反射率：辐射定标（多波段—reflectance—BIL）—大气校正（QUAC）—band math—【0.356*b1+0.13*b3+0.373*b4+0.085*b5+0.072*b7-0.0018】—提取为tif

Band1-7为定标结果中的波段，如下所示

1. 辐射定标（多波段—radiance—BIL—apply flaash）——大气校正（QUAC）
2. 计算NDVI【ENVI的NDVI功能red4 nir5】输入多波段大气校正结果—提取成tif
3. 计算FVC—【(b1 gt 0.7)*1+(b1 lt 0.05)*0+(b1 ge 0.05 and b1 le 0.7)*((b1-0.05)/(0.7-0.