meteinfolab展示MOD04_3K气溶胶光学厚度问题

wydyyh

刚开始研究气溶胶光学厚度，下载了modis数据产品MOD04_3K，选了整个中国的，下载了部分数据。编写meteoinfolab的脚本，如下：（也是参考王亚强老师的），现在遇到一些问题，程序运行完，显示的图片如下：为啥没有一点颜色变化呢？下载的MOD04_3K数据需要做什么处理吗？我看网上写的几何校正，矢量裁剪等，我是直接拿来作图显示了，会有什么问题吗？希望大家帮忙看看


f = addfile('/Users/wangyudong/waming/data/modis/0627/MOD04_3K.A2019178.1755.061.2019182203207.hdf')
v_aod = f['Optical_Depth_Land_And_Ocean']
#v_aod = f['AOD_550_Dark_Target_Deep_Blue_Combined_Mean']

aod = v_aod[[3,60],[75,135]]
x = linspace(-180.0,180.0, 360)
y = linspace(-90.0,90.0, 180)
#Set x/y dimensions
v_aod.setdim('Y', y)
v_aod.setdim('X', x)

axesm()
china = shaperead('shengjie.shp')
geoshow(china, linecolor='gray')

levels = arange(0, 1.0, 0.02)
layer = contourfm(aod, levels, interpolate=True)
title('Aerosol Optical Depth at 550 nm')
colorbar(layer)

wydyyh

应该是下载的数据 没有意义，全是－9999