MOD08_M3产品处理求助交流贴

tobewealthy

大四生正在准备毕业论文，需要利用MOD08_M3月数据气溶胶产品计算吸收性气溶胶季度均值与年均值（AAOD = (1 - SSA) * AOD），之前没接触过meteoinfo和python，编程小白，希望可以和大家互相交流数据处理心得。
结合王老师给的脚本：http://bbs.06climate.com/forum.p ... B%C6%F8%C8%DC%BD%BA
做了一点修改，改正空间还有很大，希望大家可以批评指正，代码如下：

import glob
year = 2008
sday = 032
eday = 365
aods = []
ssas = []
while sday <= eday:
fn = glob.glob('D:\modis\xi\MOD08_M3.A%i%i*.hdf' % (year, sday))[0]
print fn
f = addfile(fn)
v_aod = f['Aerosol_Optical_Depth_Land_Ocean_Mean_Mean']
v_ssa = f['Deep_Blue_Single_Scattering_Albedo_Land_Mean_Mean']

#Get data
aod = v_aod['20:40', '90:115']
ssa = v_ssa['20:40', '90:115']
aods.append(aod)
ssas.append(ssa)
aaod = (1 - ssa)* aod
for i in range(1,12):
  if (i%2 == 1):
   sady += 31
  if (i%2 == 0 and i != 2):
   sday += 30
  else:
   sday += 29
#Mean
aaod = mean(aaods)

#Plot
axesm()
geoshow('cn_province', linecolor='gray')
geoshow('country')
levels = arange(0, 1.0, 0.02)
layer = imshowm(result, levels)
title('Aerosol Optical Depth year_mean at 550 nm')
colorbar(layer)
savefig('D:/modis/year_mean08.pdf', 500, 600)
        

CHL1998

为你点赞，一起努力，我来啦

tobewealthy

CHL1998 发表于 2020-3-13 21:04
为你点赞，一起努力，我来啦

这小帖子我都想删了，不知道咋删

CHL1998

一步一步来，也许，她会记录下我们的整个过程，以后想起来也会觉得很美妙

张张要努力学习

你好，能否问一下，aaod = (1 - ssa)* aod  这一步是为什么呀？

tobewealthy

张张要努力学习 发表于 2020-5-19 20:57
你好，能否问一下，aaod = (1 - ssa)* aod  这一步是为什么呀？

我的论文是研究aaod的，这是一个计算式，不过后来数据有所更改，不用modis了

zhuyidong

tobewealthy 发表于 2020-5-21 14:04
我的论文是研究aaod的，这是一个计算式，不过后来数据有所更改，不用modis了

您好，这个for i in range(1,12):这循环是干啥的

tobewealthy

zhuyidong 发表于 2022-4-12 19:39
您好，这个for i in range(1,12):这循环是干啥的

计算月平均