MOD09A1数据下载与预处理-地表干湿度指数的计算

最新推荐文章于 2025-09-09 17:10:57 发布

空中旋转篮球

最新推荐文章于 2025-09-09 17:10:57 发布

阅读量6.3k

点赞数 7

CC 4.0 BY-SA版权

分类专栏： Python遥感图像处理文章标签： MOD09A1 干湿度计算

本文链接：https://blog.youkuaiyun.com/soderayer/article/details/120950350

Python遥感图像处理专栏收录该内容

52 篇文章 ¥29.90 ¥99.00

订阅专栏

超级会员免费看

本文介绍了MOD09A1数据，包括数据来源、波段信息及下载处理步骤。通过MRT工具进行数据预处理，转换为TIFF格式。接着详细讲述了如何计算湿度指数Wet和干度指数NDBSI，提供了计算公式，并说明了在ENVI Band Math中进行计算的方法。

1.MOD09A1数据介绍

1.1 数据介绍网址：LP DAAC - Search

中分辨率成像光谱仪 (MODIS) Terra MOD09A1 6.1 版产品提供了对 Terra MODIS 波段 1 至 7 的表面光谱反射率的估计，并针对大气条件（如气体、气溶胶和瑞利散射）进行了校正。除了七个 500 米 (m) 反射带之外，还有两个质量层和四个观测带。对于每个像素，从 8 天合成周期内的所有采集中选择一个值。像素选择的标准包括云天顶和太阳天顶。当多次采集满足标准时，使用具有最小通道 3（蓝色）值的像素。（原文：he Moderate Resolution Imaging Spectroradiometer (MODIS) Terra MOD09A1 Version 6.1 product provides an estimate of the surface spectral reflectance of Terra MODIS Bands 1 through 7 corrected for atmospheric conditions such as gasses, aerosols, and Rayleigh scattering. Along with the seven 500 meter (m) reflectance bands are two quality layers and four observation bands.

了解本专栏

超级会员免费看

37 条评论

12223zdgg 2022.05.04
您好，请问mod09a1计算干度的时候，没有中红外波段吧，那是怎么计算的呀，还有近红外波段是指近红外1还是近红外2啊
- 空中旋转篮球回复12223zdgg 2022.05.05
  是的
- 12223zdgg回复空中旋转篮球 2022.05.05
  是近红外用NIR-1,中红外用SWIR-1吗？我好像是这样用的，结果似乎有问题[face]emoji:029.png[/face]
- 空中旋转篮球回复12223zdgg 2022.05.04
  这里近红外用的是短波那个，即第2波段，中红外用的短波红外，波长较短的那个短波红外

空中旋转篮球 2022.04.27
理论上都应该要乘以系数转化为真实反射率，有除法的实际上抵消系数了，结果应该不影响，湿度指数没乘以数就会很大的。

灵凌零 2022.04.27
您好，请问是转换成tif格式，直接用这个公式计算提取指标么？需不需要×0.0001，我出来的值非常非常大，接下来需要将它控制到-1到1吗？跪求😭
- 空中旋转篮球回复灵凌零 2022.04.27
  理论上都应该要乘以系数转化为真实反射率，有除法的实际上抵消系数了，结果应该不影响，湿度指数没乘以数就会很大的。

六月十三见 2022.04.21
您好，是计算指标再矢量裁减还是先将一个一个的波段裁剪再计算呢
- 六月十三见回复空中旋转篮球 2022.04.21
  不过我今天按照您的wet公式计算数值很大不在-1～1，是要波段计算除以10000是嘛。我除以了10000之后出现了大于1的数，我又用去除异常值公式去掉了。可是统计结果显示概率从40%跳到了99%，因为有一个值数量过大导致，我该怎么办呢
- 六月十三见回复空中旋转篮球 2022.04.21
  解释的这么详细非常感谢
- 六月十三见回复空中旋转篮球 2022.04.21
  多谢大佬
- 空中旋转篮球回复六月十三见 2022.04.21
  不需要，按照这个公式你乘以0.0001结果不变（分子分母都乘以这个数结果一样），可能还增加了计算量
- 空中旋转篮球回复六月十三见 2022.04.21
  这个都可以，各有各的好处，先计算再裁剪，边界地带后期好处理一点（比如不同数据边界不一致问题等），先裁剪再计算数据量降低，效率会高一点，这个数据不是很大，可以先计算在最后裁剪得到自己研究区结果。
- 六月十三见回复六月十三见 2022.04.21
  波段用不用先乘以0·0001呢

立志学会编程 2022.04.01
你好，请问批量处理三个月份的数据，利用IDL的代码的能不能请教一下

weixin_45813650 2022.03.24
您好！请问可以用这个处理方法做NDVI嘛
- weixin_45813650回复空中旋转篮球 2022.03.25
  好的！谢谢您！
- 空中旋转篮球回复weixin_45813650 2022.03.25
  换一下公式和波段，一样可以计算的

空中旋转篮球 2022.03.18
这里的公式我记得是徐涵秋教授文章中定义的，用于计算RSEI指数的，公式是固定的，如果只用裸土指数就不用除以2了。

记忆~ 2022.03.18
你好研究区是个山区如果计算干度只用裸土指数以后还要/2吗

楠楠难难 2022.01.20
您好，这个数据mod09A1不需要进行大气以及辐射校核了嘛？
- 空中旋转篮球回复m0_49155171 2022.03.10
  处理好了的，一般可以不用做。
- m0_49155171回复空中旋转篮球 2022.03.10
  你好，那mod09A1需要进行几何校正吗
- dede_hotpot回复空中旋转篮球 2022.02.22
  好滴明白啦！谢谢你！
- 空中旋转篮球回复dede_hotpot 2022.02.21
  可以查看数据说明https://lpdaac.usgs.gov/products/mod09a1v061/，应该要除以10000，Scale Factor是0.0001，数值范围是-100-16000.
- dede_hotpot回复空中旋转篮球 2022.02.21
  您好，请问mod09a1提取第三波段得到地表反射率数据，使用前是否需要真值转换呢，因为从envi中查看的这个数值比较大，不在地表反射率数值的范围内
- 楠楠难难回复空中旋转篮球 2022.01.21
  哦哦，好的！
- 空中旋转篮球回复楠楠难难 2022.01.20
  这个已经是地表反射率产品了，不需要另外做

mramir 2022.01.08
SI = (float((b3+b5)-(b1+b4)))/((b3+b5)+(b1+b4)) 你好在上面湿度计算当中B1代表红波段来着但在这里说b1是率波段，是这样没问题吗
- 空中旋转篮球回复mramir 2022.01.08
  按照当前部分对应波段，前后不同的，分开看！