MODIS AOD 数据出图很奇怪

侠骨柔情

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

begin

fn="SATE_L3_EOT_MODIS_MWB_MOD08M3_GLB_C5-201001.nc"

fi=addfile(fn,"r")

v=fi->Optical_Depth_Land_And_Ocean_Mean_Mean

v1=v(:,:)

wks= gsn_open_wks("x11","AOD-2010-01")

res=True

res@cnFillOn = True

res@cnLinesOn=False

res@cnLevelSelectionMode = "ManualLevels" ; set manual contour levels

res@cnMinLevelValF   =0              ; set min contour level

res@cnMaxLevelValF   =200

res@cnLevelSpacingF  =20               ; set contour spacing

plot=gsn_csm_contour_map(wks,v1,res)

end

能否解释一下这是怎么回事？
第一次用MODIS数据。谢谢各位大神了。

尽头的尽头

数据是short，转换成float

侠骨柔情

尽头的尽头 发表于 2015-11-26 16:13
数据是short，转换成float

转过了，但还是不太对，感觉南半球画到了北半球。

侠骨柔情

尽头的尽头 发表于 2015-11-26 16:13
数据是short，转换成float

如果加这么一段：

v1=new((/180,360/),"float")

do i=0,179

v1(179-i,:)=v(i,:)

end do

就是颠倒了南北顺序，然后出图是这样的。

尽头的尽头

侠骨柔情 发表于 2015-11-26 16:44
如果加这么一段：
v1=new((/180,360/),"float")do i=0,179v1(179-i,:)=v(i,:)[ ...

modis数据我也没用过哦。。。具体也不清楚

侠骨柔情

尽头的尽头 发表于 2015-11-26 20:17
modis数据我也没用过哦。。。具体也不清楚

不过还是谢谢你热心的帮助！

fan

你需要给读出来的变量添加坐标变量，然后设置坐标变量的正确范围区间，这样才能和底图对上

天下唯我

我现在还无法下载（自动回复：请不要使用迅雷等下载工具，点我查看下载帮助）下来modis数据，求指教~我现在也是需要用modis反演AOD

爱高斯拉普拉斯

表示本人是来学习如何批量提取AOD，结果花了5金钱下了pdf，，，，

爱高斯拉普拉斯

不过还是学习了