MeteoInfoLab脚本示例:闪电位置图
2015-09-28 16:27
645 查看
这个脚本示例读取文本格式的闪电数据,读出每条闪电记录的经纬度和强度,在地图上绘制出每个闪电的位置,并用符号和颜色区分强度正负。
数据格式如下:
0 2009-06-06 00:01:16.6195722 纬度=35.0297 经度=114.8923 强度=23.4 陡度=9.1 误差=71.9 定位方式:四站算法
1 2009-06-06 01:34:58.5911760 纬度=34.4243 经度=115.6057 强度=-1.3 陡度=-0.4 误差=0.0 定位方式:二站振幅
2 2009-06-06 01:36:59.6324818 纬度=35.9314 经度=114.9388 强度=-35.2 陡度=-6.1 误差=0.0 定位方式:三站混合
3 2009-06-06 01:54:25.2727088 纬度=34.4163 经度=115.6121 强度=-1.9 陡度=-0.3 误差=0.0 定位方式:二站振幅
4 2009-06-06 01:58:35.1179896 纬度=35.3270 经度=106.5536 强度=-384.1 陡度=-80.6 误差=0.0 定位方式:二站振幅
5 2009-06-06 01:59:15.9764295 纬度=34.5339 经度=114.4466 强度=58.4 陡度=26.7 误差=0.0 定位方式:三站混合
7 2009-06-06 02:00:35.0869527 纬度=34.3934 经度=115.5976 强度=-3.3 陡度=-1.0 误差=0.0 定位方式:二站振幅
8 2009-06-06 02:00:43.0728269 纬度=34.4124 经度=115.7238 强度=-5.3 陡度=-3.1 误差=0.0 定位方式:二站混合
9 2009-06-06 02:00:43.2227403 纬度=34.4131 经度=115.7171 强度=-8.4 陡度=-4.1 误差=0.0 定位方式:二站混合
10 2009-06-06 02:02:21.5951061 纬度=34.3817 经度=115.6268 强度=-7.2 陡度=-1.6 误差=0.0 定位方式:二站振幅
脚本程序:
数据格式如下:
0 2009-06-06 00:01:16.6195722 纬度=35.0297 经度=114.8923 强度=23.4 陡度=9.1 误差=71.9 定位方式:四站算法
1 2009-06-06 01:34:58.5911760 纬度=34.4243 经度=115.6057 强度=-1.3 陡度=-0.4 误差=0.0 定位方式:二站振幅
2 2009-06-06 01:36:59.6324818 纬度=35.9314 经度=114.9388 强度=-35.2 陡度=-6.1 误差=0.0 定位方式:三站混合
3 2009-06-06 01:54:25.2727088 纬度=34.4163 经度=115.6121 强度=-1.9 陡度=-0.3 误差=0.0 定位方式:二站振幅
4 2009-06-06 01:58:35.1179896 纬度=35.3270 经度=106.5536 强度=-384.1 陡度=-80.6 误差=0.0 定位方式:二站振幅
5 2009-06-06 01:59:15.9764295 纬度=34.5339 经度=114.4466 强度=58.4 陡度=26.7 误差=0.0 定位方式:三站混合
7 2009-06-06 02:00:35.0869527 纬度=34.3934 经度=115.5976 强度=-3.3 陡度=-1.0 误差=0.0 定位方式:二站振幅
8 2009-06-06 02:00:43.0728269 纬度=34.4124 经度=115.7238 强度=-5.3 陡度=-3.1 误差=0.0 定位方式:二站混合
9 2009-06-06 02:00:43.2227403 纬度=34.4131 经度=115.7171 强度=-8.4 陡度=-4.1 误差=0.0 定位方式:二站混合
10 2009-06-06 02:02:21.5951061 纬度=34.3817 经度=115.6268 强度=-7.2 陡度=-1.6 误差=0.0 定位方式:二站振幅
脚本程序:
fn = 'D:/Temp/ascii/lighting/2009_06_06.txt' tf = open(fn) lats = [] lons = [] vs = [] for aline in tf: datalist = aline.split() lat = float(datalist[3].split('=')[1]) lon = float(datalist[4].split('=')[1]) v = float(datalist[5].split('=')[1]) lats.append(lat) lons.append(lon) vs.append(v) lon = array(lons) lat = array(lats) v = array(vs) axesm() mlayer = shaperead('D:/Temp/map/bou2_4p.shp') geoshow(mlayer) ss = makesymbolspec('point', {'value':(-10000,0), 'color':'b', 'marker':'m', 'size':6}, \ {'value':(0,10000), 'color':'r', 'marker':'+', 'size':6}) layer = scatterm(lon, lat, v, symbolspec=ss) xlim(90, 130) ylim(20, 50) title(u'闪电位置图', fontname=u'黑体', fontsize=18)
相关文章推荐
- unity3d 动画卡帧 动画合成 动画层次
- android 基础4
- JAVA中关于wait()和sleep()
- 新浪sae开发不能写入文件
- yum
- Linux下随机10字符病毒的清除
- 2015 ACM/ICPC合肥网赛&HDU5491 The Next
- 寻找字符串中最长回文——Manacher算法及其Java实现
- 数据结构和算法学习(2)-时间复杂度
- 字符串压缩
- Java使用jdbc连接和操作mysql的例子
- android学习地址和文章总结
- 晨遇
- 差分约束POJ3159Candies(spfa+stack)解题报告
- windows多线程编程简介(2)
- UVA 815 Flooded!
- SQL中获取树形结构所有子级数据
- Java基础知识强化之IO流笔记04:throw和throws的区别
- fseek函数
- Linux学习篇-SVN命令详解