请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 11765|回复: 19

modis 多年平均气溶胶光学厚度

[复制链接]

新浪微博达人勋

发表于 2017-10-30 19:22:04 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
@MeteoInfo 老师您好,前段时间看了您写的用脚本文件处理MOD08和MOD04的气溶胶光学厚度 是针对某一天或者某一个月的,但是最近准备做一个多年尺度的,2000年到2010年这十年期间的平均气溶胶光学厚度的时空分布特征,老师您和我说用循环读出所有文件中的aod数组,放入一个列表中,用mean函数来求平均,但自己尝试真的写不出来(刚接触强大的MeteoInfo,还不是很熟悉),非常希望老师您能写个脚本文件,论坛里的朋友应该也有需要这个脚本文件的,感激涕零!!在线等!!劳烦老师了!!非常感谢~~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-1-2 16:04:20 | 显示全部楼层
sfz13851581045 发表于 2017-12-29 16:09
好的,我再发个帖子请教一下王老师

你好,可以加你的QQ吧?相互讨论一下,我的QQ 834630747
密码修改失败请联系微信:mofangbao
回复 支持 0 反对 1

使用道具 举报

新浪微博达人勋

发表于 2017-10-31 00:41:47 | 显示全部楼层
这两天很忙,周末吧
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-10-31 01:03:41 | 显示全部楼层
最好,最简便的方法就是把这些文件按照时间维合并成一个,然后就行循环读取,求平均,应该很容易,你的是MOD08月产品?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-10-31 15:12:15 | 显示全部楼层
同等同关注中。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-11-2 16:13:55 | 显示全部楼层
MeteoInfo 发表于 2017-10-31 00:41
这两天很忙,周末吧

麻烦老师了!!感激!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-11-3 17:54:26 | 显示全部楼层
参考此脚本:
  1. import glob
  2. year = 2015
  3. sday = 121
  4. eday = 122
  5. aods = []
  6. while sday <= eday:
  7.     fn = glob.glob('D:/Temp/hdf/MOD08_D3.A%i%i*.hdf' % (year, sday))[0]
  8.     print fn
  9.     f = addfile(fn)
  10.     v_aod = f['Optical_Depth_Land_And_Ocean_Mean']
  11.     #Get data
  12.     aod = v_aod[:,:]
  13.     aods.append(aod)
  14.     sday += 1

  15. #Mean
  16. aod = mean(aods)

  17. #Set x/y
  18. x = linspace(-180.0,180.0, 360)
  19. y = linspace(-90.0,90.0, 180)

  20. #Plot
  21. axesm()
  22. geoshow('country', linecolor='gray')
  23. levels = arange(0, 1.0, 0.02)
  24. layer = imshowm(x, y, aod, levels)
  25. title('Aerosol Optical Depth at 550 nm')
  26. colorbar(layer)


密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-11-3 23:14:15 | 显示全部楼层

老师您好,非常感谢您百忙之中能写脚本,运行的时候出现了错误您看看>>> run script...
Traceback (most recent call last):
  File "<iostream>", line 7, in <module>
IndexError: index out of range: 0


想请问一下您:1)sday = 121 eday = 122 这个开始的和结束的天 121 和122是什么意思呢?2)我的数据格式是这样的,MOD08_D3.A2016133.006.2016184004724.hdf  这个应该怎么写呢? 谢谢您~~
报错.jpg
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-11-6 17:33:37 | 显示全部楼层
我觉得http://www.ncl.ucar.edu/Applications/Scripts/hdf4sds_3.ncl这个脚本,可能比较方便一点,,,
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-12-27 22:37:10 来自手机 | 显示全部楼层
我试了一下脚本,好像也是第七行报错
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-12-27 22:37:15 来自手机 | 显示全部楼层
我试了一下脚本,好像也是第七行报错
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表