基于MOD02/MYD02获得亮度温度再转冰温

就是一只白

已于 2024-01-17 18:56:19 修改

阅读量875

点赞数 9

分类专栏： python RS 文章标签： python 遥感

于 2024-01-15 05:10:58 首次发布

本文链接：https://blog.csdn.net/xinemachen/article/details/135592020

版权

python 同时被 2 个专栏收录

50 篇文章 4 订阅

订阅专栏

3 篇文章 0 订阅

订阅专栏

用HEG处理MOD02/MYD02，提取里面的EV_1KM_Emissive波段，band为11和12（其实就是band 31和32）。注意这里的band和output dile type

1. 获得之后，转辐射亮度。

参考：https://www.cnblogs.com/enviidl/p/16539422.html

radiance_scales，和radiance_offset这两项参数代表波段的增益和偏移量，是辐射定标的系数。比如要计算31波段的辐射亮度，读取到scales为0.00084002，offsets为1577.33972168，带入MODIS辐射定标的通用公式：Radiance=scales*(DN-offsets)，即可以得到该波段的辐射亮度。

2. 辐射亮度转亮温

参考：将图像DN值定标维热辐射强度之后,可用Planck函数求解出星上亮度温度_亮度温度计算公式-CSDN博客

用普朗克方程可以求出

使用bandmath工具计算31和32的亮温。

31波段亮温=1304.413871/alog(1+729.541636/b31)

B31: 31波段辐射亮度值

32波段亮温=1196.978785 /alog(1+474.684780/b32)

就可以获得亮温了

3. 亮温转冰温

https://nsidc.org/sites/default/files/myd29-v005-userguide_0.pdf

这里的角度是传感器扫描角度

_,_,_,_,zenith_data=read_tif(sensorzenithpath)
zenith_data=zenith_data.astype(float) * 0.01
R0=6378 #km
height=705 #km
temp=(R0*1.0/(R0+height))*np.sin(np.deg2rad(zenith_data))
angle=np.arcsin(temp)

可以从MOD02产品中获得sensor zeith可以转为角度

参考：MODIS传感器：扫描角、传感器天顶角（sensor zenith），MODTRAN里的sensor zenith - 知乎