卫星云图显示

shenwsshen

最近在做卫星云图的显示，从国家卫星气象中心下载的风云2C卫星，hdf5格式，VISSR数据，通过hdfvierer查看的时候，发现其中的NOMChannelVIS1KM 数据中一半是0，其他的最大才51，用Python把数据输出为图片后，效果与国家卫星中心数据的预览图比较如下：


我听说VISSR的坐标貌似需要转换，但是不知道有没有Python的函数可以调用，另外云的亮度好像差很多，不知道VISSR中的值到rgb值是不是还有其他转换，这个转换是自定义的，还是有约定呢？

burf

我就不多说废话了，卫星中心的HDF5数据是全球标称数据，2288*2288的格点，每个格点对应一个经纬度。对应的经纬度可以从卫星中心上下下来，或者从我这里附件下载。 NOM_ITG_2288_2288(0E0N)_LE.dat (39.94 MB, 下载次数: 165)

2015-8-13 20:09 上传

点击文件名下载附件 NOM_ITG_2288_2288(0E0N)_LE.hdr (243 Bytes, 下载次数: 40)

2015-8-13 20:09 上传

点击文件名下载附件 Readme for FY latlon.txt (1.64 KB, 下载次数: 43)

2015-8-13 20:09 上传

点击文件名下载附件

数据有三个，一个标称网格的经纬度，一个包含了星下点坐标信息的文件，一个描述文件（txt）。
读取方式如下：
第一步：
读取经纬度文件，这里请用numpy的fromfile函数进行读取，当然也可以用别的方法，我这里就说这个比较简单的。
loc=np.fromfile('NOM_ITG_2288_2288(0E0N)_LE\\NOM_ITG_2288_2288(0E0N)_LE.dat',N.float32,count=-1)
latlon=loc.reshape(2,2288,2288)
latlon[latlon==300.]=np.nan

最后得到的latlon数组既是最后能够用到的经纬度，需要注意的是，数据中的缺省部分用300来表示，这里我都换成nan了。
第二步：
读取星下点等信息，注意每个卫星的星下点信息不一样，下载的时候请注意。这里我给的数据是E星的。数据格式为hdf5格式。
filename=['FY2E_FDI_ALL_NOM_20150708_0345.hdf']

f=h5py.File(filename[0],'r')
info=f['NomFileInfo'][:][0]
print info[3],info[4]
latlon[1,:,:]=latlon[1,::-1,:]
latlon[0,:,:]=latlon[0,:,:]+info[4]
lon=np.array(latlon[0,:,:])
lat=np.array(latlon[1,:,:])

具体里面info有哪些内容，请楼主自己看看吧。info[4]为星下点的经度（为什么没有纬度呢？因为风云卫星在赤道上，星下点的纬度为0）。然后最后从latlon中把标称数据的每一个经纬度提取出来，就是lon和lat这两个数组了。
第三步：
画图了，这一步就简单了。从hdf5的风云数据中读取可见光啊、红外通道啊，都行。有一个需要注意的：可见光通道的缺省值为255，红外通道的缺省值为65535，在绘图室mask掉就行。其他的。。。就是选投影方案，投影就行了，有了经纬度和数值什么都好办。
最后附一张图：

shenwsshen

。。。一般直接谷歌学习，或者看文档学习，你是学习气象的么？你知道红外气象云图怎么处理么，我感觉自己处理出来的红外云图背景太亮了。。。我是直接i//4+3得到的灰度值，不知到这样是不是有问题。。。

shenwsshen

最后一张图片是将云的值加了100后的灰度显示效果。。。

shenwsshen

原来可见刚是64个显示等级，需要用原始值*4+10，坐标转换可以使用线性插值，但是找不到文件中的简化网格块。。。

woaichimian666

想问大神们是从哪里系统学习的Python入门，最近看了笨方法学Python，感觉那个偏向编程

ginanony

感觉这个。。。的确有点难度，之前的确遇到过类似问题。。。都是转存出来拿PS看了。。。同求教

burf

楼主是不是想画成这样？

burf

burf 发表于 2015-8-13 20:04
楼主是不是想画成这样？

顺便一说我这个是红外通道3，不过和可见光一个道理

burf

还有什么问题请楼主提问， 互相学习嘛。。。不过可能不能及时回复