MeteoInfo的个人博客分享 http://blog.sciencenet.cn/u/MeteoInfo

博文

MeteoInfoLab读取绘制FY4A卫星全圆盘数据

已有 566 次阅读 2022-11-28 09:05 |个人分类:MeteoInfo|系统分类:科研笔记

由于FY4A数据中空间信息缺失,需要根据数据实际情况在脚本中补充。FY卫星圆盘数据是Geostationay satellite投影,在此投影的的数据网格是矩形,x, y坐标单位是米。这个数据的分辨率是4km (4000m),数据网格中心坐标是(0, 0),x方向格点数是2748,那么x方向的总长度是2748*4000,再除以2就可以得到x方向坐标起始和结束值为(-5496000.0,5496000.0)。

fn = 'D:/Temp/FY/FY4A-_AGRI--_N_DISK_1047E_L2-_CTT-_MULT_NOM_20190209140000_20190209141459_4000M_V0001.NC'
f = addfile(fn)
x = linspace(-5496000.0,5496000.0, 2748)
y = linspace(-5496000.0,5496000.0, 2748)
data = f['CTT'][::-1,:]
data[data>1000] = nan
data[data==-999] = nan
height = f['nominal_satellite_height'][:]
 
#Plot
lon0 = 104.7
ax = axesm(proj='geos', lon_0=lon0, h=height, gridlabelloc='all', griddx=30,
    griddy=30, gridline=True, frameon=False)
geoshow('coastline', color='k')
levs = arange(160, 311, 5)
layer = imshow(x, y, data, levs, proj=ax.proj)
colorbar(layer, shrink=0.8, xshift=15)
title('FY4A Cloud Top Temperature')
 
#savefig('D:/Temp/figure/fy4a_ctt.png', 500, 400, dpi=300)

fy4a_ctt.png



https://blog.sciencenet.cn/blog-611158-1365577.html

上一篇:MeteoInfoLab脚本绘制山魈脸谱
下一篇:台风大涡模拟三维可视化
收藏 IP: 114.247.188.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2023-2-3 18:43

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部