请教：关于GrADs绘制awx格式的卫星图像偏移问题

quesera

请教坛子里的各位大神：
  使用卫星中心下载awx格式的数据 FY2E_SEC_IR3_MLS_20140725_0000.AWX ，使用fortran读取后用grads作图。仔细对比同一文件在swap平台上的图和grads作出的图后发现，两者在东西向没有偏移，而南北向则偏移了接近5°。对比一下两图中，swap图上位于湖南中部的小块云团，在grads图上南偏到了广西东部；swap图上位于台湾岛西南侧的季风云团，在grads图上也明显南偏。两张图对比发现，在最南端偏移几乎没有，越往北偏移越多，似乎偏差是纬向叠加的。

swap和micaps显示的云图一致，说明这应该是自己编程作图上出现了问题。在fortran读取awx格式数据头文件时，显示是数据是麦卡托投影，而grads中似乎没有这种投影方式。以下将ctl文件贴出，请大侠指导这到底是哪里出了问题？


dset E:\awx\FY2E_SEC_IR3_MLS_201507%d2_%h2%n2.AWX.bin
options template
options yrev
undef 999.0
xdef 512 linear 100.57 0.1344531
ydef 512 linear -10.38 0.1185352
zdef 1 linear 1000 1
tdef 9 linear 00:15z11JUL2015 6hr
vars  2
t 0 99 grey
endvars

quesera

ctl文件里的起始经纬度、步长都是由读取的头文件信息给出的，经反复检查，没有错误

MeteoInfo

可以考虑用MeteoInfoLab来做，参考此贴：MeteoInfoLab脚本示例：AWX数据
http://bbs.06climate.com/forum.p ... 312&fromuid=106
(出处: 气象家园)

Pilgrim

请问楼主能否将Fortran程序贴出来分享一下吗

quesera

MeteoInfo 发表于 2016-5-30 16:36
可以考虑用MeteoInfoLab来做，参考此贴：MeteoInfoLab脚本示例：AWX数据
http://bbs.06climate.com/forum. ...

尝试下载了等经纬度投影的云图数据，出图无误。可能是这个投影格式的云图数据，ctl得用另一种方法定义呢。正要学习meteroinfo，谢谢指教！

安达充

quesera 发表于 2016-5-30 23:29
尝试下载了等经纬度投影的云图数据，出图无误。可能是这个投影格式的云图数据，ctl得用另一种方法定义呢 ...

求分享一下fortran读取数据并出图的ctl。gs