Grads处理风云FY2静止卫星AWX数据总结

arctichare

工作中正好遇到，其实我之前习惯用matlab和ncl。AWX数据占用空间小，grads批处理快都是值得应用的。grads接触时间不长，有不那么灵巧的地方还请指教首先，FY2系列AWX数据的获取
http://satellite.nsmc.org.cn/PortalSite/Data/DataView.aspx?SatelliteType=0&DataCategoryCode=Atmosphere&DataTypeCode=ASL
可以从中挑选想要的产品，连续性还是不错的

9210格式云分类产品	AWX	--	2015-06-03	2016-07-06	11844	11.1	[url=]查看[/url]
	标称格式云分类产品	HDF	--

下载后犯了拿来主义错误，总是什么都想要现成的，后来发现次级产品中经纬度不是统一的，grads做ctl要注意。
http://satellite.nsmc.org.cn/PortalSite/StaticContent/DocumentDownload.aspx?TypeID=3
辅助文档都有写，针对AWX的有卫星产品分发格式规范2.1版，明明白白写着如
射出长波辐 射出长波辐 射(OLR) (OLR)
50 °N-50 °S， 55 °E-155 °

相当黑体亮 相当黑体亮 度温（TBBTBBTBB）
60°N-60°S， 45°E-165°E

之类的选取范围和精度，估计到这里很多帖子中的问题已经解决了。

再说CTT，PRE等一些因存储格式导致grads不能顺利读出的产品，原因有老师在之前帖子中已分析。
由于强迫症同类产品想出一致的图，我的做法比较简单粗暴，该写AWX成grads可读即可。
给出matlab处理的部分脚本

file='D:FY2G_CTT_MLT_OTG_20150727_2100.AWX';

fid = fopen(file,'rb');
line1 = fread(fid,1001);%跳过文件头
line2 = fread(fid,1001);%再次扩展段 check 是否漂移
data0 = fread(fid,[1001 1001],'ubit16');% 读取TBB数据
data0 = data0';%二进制读取读取之后需要进行一个转置操作，才能是正常经纬度放置方式data0 = data0-200;
fclose(fid); %关闭文件
caxis([0 400])
colorbar
[X Y]=meshgrid(1:1001,1:1001);
pcolor(X,Y,rot90(data0',1));
shading flat
caxis([0 400])
colorbar
% already checked OK
% rewrite for grads

fid = fopen('test.AWX', 'w');
fwrite(fid, line1);
fwrite(fid, line2);
fwrite(fid, data0,'ubit8');
fclose(fid);
% one more test

可以看到很多文段也借鉴了本坛帖子。相应的处理方法fortran等估计更简单，个人习惯而已。

还有，文档说明较老，有一些新产品我是用matlab对应下载时的样图试出范围的，很笨，不知道有没有更好的方法。

以及，IR1,VIS等数据精度显示的是5km，坛子中也有说非等经纬度的grads处理不了。
我也试了下。发现5km就是0.05出的图和样图一致。

嗯想到别的再说，先这样了。

对了，色标用的是岛国小哥的分享，因为我[url=http://bbs.06climate.com/forum.php?mod=viewthread&tid=975]点我查看给积分不够的朋友的建议[/url]无法下载（自动回复：请不要使用迅雷等下载工具，点我查看下载帮助）坛子里的脚本::>_<::http://kodama.fubuki.info/wiki/wiki.cgi/GrADS/script/color.gs?lang=en

雪颜

各位吧友注意不要照搬范围，具体看卫星星下点的位置

60°N-60°S， 45°E-165°E

只是是FY-2C用的

wuwei2163

之前的太零碎了，这个真心不错，赞一个~楼主要是能够把总结出的代码打包上传，分享给大家就更好了

贫道敬孔

多谢楼主分享{:5_213:}

墨墨卡托

感谢楼主分享

tiantianyixin

感谢楼主分享，对于卫星数据处理IDL是个很好的软件

雪颜

并且你说的
射出长波辐 射出长波辐 射(OLR) (OLR)
50 °N-50 °S， 55 °E-155 °

相当黑体亮 相当黑体亮 度温（TBBTBBTBB）
60°N-60°S， 45°E-165°E

这个在2.1规范里是FY-2C的，FY-2C是105°E，正负60所以才有这个范围
其他卫星要看具体星下点位置

雨恬魅影

非常有用，感谢分享!

香山木叶

大神有没有具体操作神马的？求拯救呀

arctichare

香山木叶 发表于 2017-1-5 16:35
大神有没有具体操作神马的？求拯救呀

你要的亮温grads脚本论坛里都有的。我就是拿来用的呀，现在哪一步走不通？说的具体点吧