MeteoInfoLab脚本示例：FY2 gpf格式数据

MeteoInfo

FY2C-E的gpf格式是二进制数据，由数据头、定标表、通道数据三部分组成，这里给出一个FY2E gpf文件（经纬度投影）读取、绘图的示例程序。先利用Python读取二进制数据文件的功能读取文件头中相关信息，然后利用MeteoInfoLab的binread函数读取通道数据至二维数组，再绘图。

1. import struct
   
2. 
   
3. fn = 'D:/Temp/binary/FY2E_2011_06_25_00_01_E_PJ3.gpf'
   
4. 
   
5. #Read data header parameters
   
6. f = open(fn, 'rb')
   
7. fileid, = struct.unpack('2s', f.read(2))
   
8. version, = struct.unpack('<h', f.read(2))
   
9. satid, = struct.unpack('<h', f.read(2))
   
10. year, = struct.unpack('<h', f.read(2))
    
11. month, = struct.unpack('<h', f.read(2))
    
12. day, = struct.unpack('<h', f.read(2))
    
13. hour, = struct.unpack('<h', f.read(2))
    
14. minute, = struct.unpack('<h', f.read(2))
    
15. chnums, = struct.unpack('<h', f.read(2))
    
16. pjtype, = struct.unpack('<h', f.read(2))
    
17. width, = struct.unpack('<h', f.read(2))
    
18. height, = struct.unpack('<h', f.read(2))
    
19. clonres, = struct.unpack('<f', f.read(4))
    
20. clatres, = struct.unpack('<f', f.read(4))
    
21. stdlat1, = struct.unpack('<f', f.read(4))
    
22. stdlat2, = struct.unpack('<f', f.read(4))
    
23. earthr, = struct.unpack('<f', f.read(4))
    
24. minlat, = struct.unpack('<f', f.read(4))
    
25. maxlat, = struct.unpack('<f', f.read(4))
    
26. minlon, = struct.unpack('<f', f.read(4))
    
27. maxlon, = struct.unpack('<f', f.read(4))
    
28. ltlat, = struct.unpack('<f', f.read(4))
    
29. ltlon, = struct.unpack('<f', f.read(4))
    
30. rtlat, = struct.unpack('<f', f.read(4))
    
31. rtlon, = struct.unpack('<f', f.read(4))
    
32. lblat, = struct.unpack('<f', f.read(4))
    
33. lblon, = struct.unpack('<f', f.read(4))
    
34. rblat, = struct.unpack('<f', f.read(4))
    
35. rblon, = struct.unpack('<f', f.read(4))
    
36. stdlon, = struct.unpack('<f', f.read(4))
    
37. centerlon, = struct.unpack('<f', f.read(4))
    
38. centerlat, = struct.unpack('<f', f.read(4))
    
39. chindex = []
    
40. for i in range(chnums):
    
41.     chindex.append(struct.unpack('b', f.read(1))[0])
    
42. f.read(128 - chnums)
    
43. plonres, = struct.unpack('<f', f.read(4))
    
44. platres, = struct.unpack('<f', f.read(4))
    
45. f.read(1808)
    
46. 
    
47. #Read calibration table
    
48. f.read(32768)
    
49. 
    
50. f.close()
    
51. 
    
52. #Read one channel data
    
53. cn = 1    #Infrared channel 1
    
54. skipn = 2048 + 32768
    
55. for i in range(1, cn):
    
56.     if chindex[i-1] < 5:    #Infrared channel
    
57.         byten = 2
    
58.     else:
    
59.         byten = 1    #Visible light channel
    
60.     skipn += width * height * byten
    
61. if cn < 5:
    
62.     data = binread(fn, [height, width], 'short', skip=skipn)
    
63. else:
    
64.     data = binread(fn, [height, width], 'byte', skip=skipn)
    
65. data = data[::-1,:]
    
66. lon = linspace(ltlon, rtlon, width)
    
67. lat = linspace(lblat, ltlat, height)
    
68. 
    
69. #Plot
    
70. axesm()
    
71. lworld = shaperead('D:/Temp/Map/country1.shp')
    
72. lchina = shaperead('D:/Temp/Map/bou2_4p.shp')
    
73. geoshow(lchina, edgecolor='b')
    
74. geoshow(lworld, edgecolor='b')
    
75. cmap = 'MPL_gist_gray'
    
76. layer = imshowm(lon, lat, data, 20, cmap=cmap)
    
77. colorbar(layer)
    
78. t = datetime.datetime(year, month, day, hour, minute)
    
79. title('FY2E (' + t.strftime('%Y-%m-%d') + ')')

import struct<br /> <br /> fn = 'D:/Temp/binary/FY2E_2011_06_25_00_01_E_PJ3.gpf'<br /> <br /> #Read data header parameters<br /> f = open(fn, 'rb')<br /> fileid, = struct.unpack('2s', f.read(2))<br /> version, = struct.unpack('<h', f.read(2))<br /> satid, = struct.unpack('<h', f.read(2))<br /> year, = struct.unpack('<h', f.read(2))<br /> month, = struct.unpack('<h', f.read(2))<br /> day, = struct.unpack('<h', f.read(2))<br /> hour, = struct.unpack('<h', f.read(2))<br /> minute, = struct.unpack('<h', f.read(2))<br /> chnums, = struct.unpack('<h', f.read(2))<br /> pjtype, = struct.unpack('<h', f.read(2))<br /> width, = struct.unpack('<h', f.read(2))<br /> height, = struct.unpack('<h', f.read(2))<br /> clonres, = struct.unpack('<f', f.read(4))<br /> clatres, = struct.unpack('<f', f.read(4))<br /> stdlat1, = struct.unpack('<f', f.read(4))<br /> stdlat2, = struct.unpack('<f', f.read(4))<br /> earthr, = struct.unpack('<f', f.read(4))<br /> minlat, = struct.unpack('<f', f.read(4))<br /> maxlat, = struct.unpack('<f', f.read(4))<br /> minlon, = struct.unpack('<f', f.read(4))<br /> maxlon, = struct.unpack('<f', f.read(4))<br /> ltlat, = struct.unpack('<f', f.read(4))<br /> ltlon, = struct.unpack('<f', f.read(4))<br /> rtlat, = struct.unpack('<f', f.read(4))<br /> rtlon, = struct.unpack('<f', f.read(4))<br /> lblat, = struct.unpack('<f', f.read(4))<br /> lblon, = struct.unpack('<f', f.read(4))<br /> rblat, = struct.unpack('<f', f.read(4))<br /> rblon, = struct.unpack('<f', f.read(4))<br /> stdlon, = struct.unpack('<f', f.read(4))<br /> centerlon, = struct.unpack('<f', f.read(4))<br /> centerlat, = struct.unpack('<f', f.read(4))<br /> chindex = []<br /> for i in range(chnums):<br />     chindex.append(struct.unpack('b', f.read(1))[0])<br /> f.read(128 - chnums)<br /> plonres, = struct.unpack('<f', f.read(4))<br /> platres, = struct.unpack('<f', f.read(4))<br /> f.read(1808)<br /> <br /> #Read calibration table<br /> f.read(32768)<br /> <br /> f.close()<br /> <br /> #Read one channel data<br /> cn = 1    #Infrared channel 1<br /> skipn = 2048 + 32768<br /> for i in range(1, cn):<br />     if chindex[i-1] < 5:    #Infrared channel<br />         byten = 2<br />     else:<br />         byten = 1    #Visible light channel<br />     skipn += width * height * byten<br /> if cn < 5:<br />     data = binread(fn, [height, width], 'short', skip=skipn)<br /> else:<br />     data = binread(fn, [height, width], 'byte', skip=skipn)<br /> data = data[::-1,:]<br /> lon = linspace(ltlon, rtlon, width)<br /> lat = linspace(lblat, ltlat, height)<br /> <br /> #Plot<br /> axesm()<br /> lworld = shaperead('D:/Temp/Map/country1.shp')<br /> lchina = shaperead('D:/Temp/Map/bou2_4p.shp')<br /> geoshow(lchina, edgecolor='b')<br /> geoshow(lworld, edgecolor='b')<br /> cmap = 'MPL_gist_gray'<br /> layer = imshowm(lon, lat, data, 20, cmap=cmap)<br /> colorbar(layer)<br /> t = datetime.datetime(year, month, day, hour, minute)<br /> title('FY2E (' + t.strftime('%Y-%m-%d') + ')')

wuwei2163

王老师，您好！通过您的这些帖子，收获很多，我现在也在不断的学习MeteoInfo并将其作为新的工具，现在有个问题向您请教：将两种类型的数据如何画在同一张图上，MeteInlab脚本该怎么写呢？比如卫星云图和降水量叠加图或者闪电定位和云图叠加等等，能不能给我们这些刚入门的小白一个关于上述MeteInlab脚本示例，让我们不停的去摸索？最后，祝您和家园的各位朋友节日快乐！！！

MeteoInfo

wuwei2163 发表于 2016-10-2 22:44
王老师，您好！通过您的这些帖子，收获很多，我现在也在不断的学习MeteoInfo并将其作为新的工具，现在有个 ...

参考此贴：http://bbs.06climate.com/forum.p ... &extra=page%3D1

wuwei2163

MeteoInfo 发表于 2016-10-3 11:18
参考此贴：http://bbs.06climate.com/forum.php?mod=viewthread&tid=44833&extra=page%3D1

谢谢王老师，准备好好去研究了