请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5026|回复: 7

MeteoInfoLab脚本示例:Himawari standard data

[复制链接]

新浪微博达人勋

发表于 2017-9-21 06:21:03 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 MeteoInfo 于 2018-9-12 09:25 编辑

Himawari卫星数据提供两种格式:netCDF和HSD数据(Himawari standard data),数据下载需要先在此网站注册:http://www.eorc.jaxa.jp/ptree/index.html,然后用提供的FTP下载数据。HSD数据是二进制数据,数据格式的详细描述见此文档:http://www.data.jma.go.jp/mscweb ... rs_guide_en_v12.pdf。这里示例读取HSD数据并绘图:


  1. import struct

  2. def read_h8(fn):
  3.     #Read data header   
  4.     f = open(fn, 'rb')
  5.     hlen = 0
  6.     #1 Basic information block
  7.     f.read(282)
  8.     hlen += 282
  9.     #2 Data information block
  10.     f.read(5)
  11.     ncol, = struct.unpack('<h', f.read(2))
  12.     nrow, = struct.unpack('<h', f.read(2))
  13.     f.read(41)
  14.     hlen += 50
  15.     #3 Projection information block
  16.     #f.read(127)
  17.     f.read(19)
  18.     sx, = struct.unpack('<f', f.read(4))
  19.     sy, = struct.unpack('<f', f.read(4))
  20.     f.read(127 - 27)
  21.     hlen += 127
  22.     #4 Navigation information block
  23.     f.read(139)
  24.     hlen += 139
  25.     #5 Calibration information block
  26.     f.read(147)
  27.     hlen += 147
  28.     #6 Inter-calibration information block
  29.     f.read(259)
  30.     hlen += 259
  31.     #7 Segment information block
  32.     #f.read(47)
  33.     f.read(3)
  34.     tns, = struct.unpack('b', f.read(1))
  35.     ssn, = struct.unpack('b', f.read(1))
  36.     fln, = struct.unpack('<h', f.read(2))
  37.     f.read(40)
  38.     hlen += 47
  39.     #8 Navigation correction information block
  40.     f.read(1)
  41.     blen, = struct.unpack('<h', f.read(2))
  42.     f.read(blen - 3)
  43.     hlen += blen
  44.     #9 Observation time information block
  45.     f.read(1)
  46.     blen, = struct.unpack('<h', f.read(2))
  47.     f.read(blen - 3)
  48.     hlen += blen
  49.     #10 Error information block
  50.     f.read(1)
  51.     blen, = struct.unpack('<h', f.read(2))
  52.     f.read(blen - 3)
  53.     hlen += blen
  54.     #11 Spare block
  55.     f.read(259)
  56.     hlen += 259
  57.    
  58.     f.close()

  59.     #Read data
  60.     data = binread(fn, [nrow, ncol], 'short', skip=hlen)
  61.     data = data.astype('float')
  62.     data[data<0] = nan
  63.     return data, ncol, nrow, fln

  64. #Read data files
  65. segments = range(1, 11)
  66. for segment in segments:
  67.     fn = 'E:/Temp/himawari8/HS_H08_20170921_0410_B16_FLDK_R20_S%02i10.DAT' % segment
  68.     print fn
  69.     data1,ncol,nrow1,fln1 = read_h8(fn)
  70.     if segment == segments[0]:
  71.         data = data1
  72.         fln = fln1
  73.         nrow = nrow1
  74.     else:
  75.         data = concatenate([data, data1], axis=0)
  76.         nrow += nrow1
  77. data = data[::-1,:]

  78. #Plot
  79. sx = -5500000
  80. sy = 5500000 - segments[-1] * 550 * 2000
  81. x = arange1(sx, ncol, 2000)
  82. y = arange1(sy, nrow, 2000)
  83. ax = axesm(proj='geos', lon_0=140.7, h=35785863, gridlabel=True, gridline=True, frameon=False)
  84. geoshow('country', edgecolor='b')
  85. cmap = 'MPL_gist_gray'
  86. levs = arange(800, 2001, 50)
  87. layer = imshowm(x, y, data, levs, cmap=cmap, proj=ax.proj)
  88. colorbar(layer, shrink=0.8)


himawari8_hsd.png



密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-9-21 06:28:49 | 显示全部楼层
紧跟业务需求
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-9-21 07:55:49 | 显示全部楼层
王老师真给力!!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-9-21 10:01:30 | 显示全部楼层
王老师威武!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-9-22 22:08:11 | 显示全部楼层
多谢大神分享
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-9-25 06:10:06 | 显示全部楼层
膜拜{:eb502:}
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2022-3-7 15:43:26 | 显示全部楼层
王老师威武!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-3-7 16:53:15 | 显示全部楼层
本帖最后由 小渔007 于 2022-3-7 16:57 编辑

王老师,请问下载了FileZilla软件但是为什么连不上ftp呢,是什么原因?
微信图片_20220307165409.png
微信图片_20220307165528.png
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表