使用Python绘制台风轨迹图的示例代码


Posted in Python onSeptember 21, 2020

参考:

2.Basemap生成的图中绘制轨迹

使用CMA热带气旋最佳路径数据集,对我国周边的台风进行绘制

import re
import os
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap




path=r"E:\Computer Science\数学建模\第二次模拟赛题\附件"
files= os.listdir(path) #得到文件夹下的所有文件名称
data=[]
all=[]
for file in files: #遍历文件夹
  if not os.path.isdir(file): # 判断是否是文件夹,不是文件夹才打开
    f = open(path + "/" + file) # 打开文件
    tmp=f.readlines()
    for i in tmp:
      line=i.split()
      if(line[0]=='66666'):
        if(len(all)>0):
          data.append(all)
        # print(line)
        all=[]
      else:
        we=(int(line[2])*0.1,int(line[3])*0.1)
        all.append(we)
# print(data)
print(len(data))
CHN='E:\Computer Science\数学建模\python_basemap'
plt.figure(figsize=(20,12))
map=Basemap(llcrnrlon=70,llcrnrlat=2,urcrnrlon=170,urcrnrlat=58)
map.drawcoastlines()
map.drawcountries()
#添加河流
# map.drawrivers(color='blue',linewidth=0.3)
#添加大陆
map.readshapefile(CHN+'\gadm36_CHN_shp\gadm36_CHN_1',
         'states',color='blue',drawbounds=True)
map.readshapefile(CHN+'\gadm36_TWN_shp\gadm36_TWN_1',
         'taiwan',color='blue',drawbounds=True)
#添加经纬线
parallels = np.linspace(3,55,5)
# print(parallels)
map.drawparallels(parallels,labels=[False,True,False,False],fontsize=5)
meridians = np.linspace(70,170,5)
# print(meridians)
map.drawmeridians(meridians,labels=[False,False,False,True],fontsize=5)
plt.rcParams['savefig.dpi'] = 300 #图片像素
plt.rcParams['figure.dpi'] = 300 #分辨率

ans=1
x=[]
y=[]
for typhoon in data:
  length=len(typhoon)
  print("%d is process!" % ans)
  ans += 1
  for i in range(length):
    x.append(typhoon[i][1])
    y.append(typhoon[i][0])
  # print(x,y)
  # map.plot(x, y, 'c*-', linewidth=2)
  map.plot(x, y, color='r',linewidth=1.5)
  x = []
  y = []
map.fillcontinents()
plt.title(r'$China\ Typhoon$',fontsize=24)
# plt.ylim(70, 170)
# plt.xlim(2, 58)
plt.show()

效果图:

使用Python绘制台风轨迹图的示例代码

以上就是使用Python绘制台风轨迹图的示例代码的详细内容,更多关于Python绘制轨迹图的资料请关注三水点靠木其它相关文章!

Python 相关文章推荐
在Python中处理字符串之isdigit()方法的使用
May 18 Python
Python中set与frozenset方法和区别详解
May 23 Python
Python编程实现二分法和牛顿迭代法求平方根代码
Dec 04 Python
python实现隐马尔科夫模型HMM
Mar 25 Python
python远程调用rpc模块xmlrpclib的方法
Jan 11 Python
简单介绍python封装的基本知识
Aug 10 Python
PyTorch中 tensor.detach() 和 tensor.data 的区别详解
Jan 06 Python
python scatter函数用法实例详解
Feb 11 Python
树莓派升级python的具体步骤
Jul 05 Python
Python常用数据分析模块原理解析
Jul 20 Python
Python列表推导式实现代码实例
Sep 09 Python
正确的理解和使用Django信号(Signals)
Apr 14 Python
利用python绘制中国地图(含省界、河流等)
Sep 21 #Python
python利用google翻译方法实例(翻译字幕文件)
Sep 21 #Python
python统计mysql数据量变化并调用接口告警的示例代码
Sep 21 #Python
python3.8动态人脸识别的实现示例
Sep 21 #Python
如何在Anaconda中打开python自带idle
Sep 21 #Python
Python 实现国产SM3加密算法的示例代码
Sep 21 #Python
python如何实现DES加密
Sep 21 #Python
You might like
DC最新动画电影:《战争之子》为何内容偏激,毁了一个不错的漫画
2020/04/09 欧美动漫
YII模块实现绑定二级域名的方法
2014/07/09 PHP
PHP登录验证码的实现与使用方法
2016/07/07 PHP
PHP中危险的file_put_contents函数详解
2017/11/04 PHP
PHP-FPM和Nginx的通信机制详解
2019/02/01 PHP
PHP lcfirst()函数定义与用法
2019/03/08 PHP
Javascript 获取字符串字节数的多种方法
2009/06/02 Javascript
使用jQuery Ajax功能时需要注意的一个问题(内存溢出)
2012/05/30 Javascript
jQuery根据元素值删除数组元素的方法
2015/06/24 Javascript
JS实现的通用表单验证插件完整实例
2015/08/20 Javascript
jQuery点击输入框显示验证码图片
2016/05/19 Javascript
JavaScript数组复制详解
2017/02/02 Javascript
Angular使用过滤器uppercase/lowercase实现字母大小写转换功能示例
2018/03/27 Javascript
JS实现用特殊符号替换字符串的中间部分区域的实例代码
2018/07/24 Javascript
WebGL three.js学习笔记之阴影与实现物体的动画效果
2019/04/25 Javascript
详解vue-property-decorator使用手册
2019/07/29 Javascript
vue基于better-scroll实现左右联动滑动页面
2020/06/30 Javascript
python自动zip压缩目录的方法
2015/06/28 Python
Python利用matplotlib生成图片背景及图例透明的效果
2017/04/27 Python
Python3中简单的文件操作及两个简单小实例分享
2017/06/18 Python
Python基于多线程实现ping扫描功能示例
2018/07/23 Python
Python函数式编程指南:对生成器全面讲解
2019/11/19 Python
使用python处理题库表格并转化为word形式的实现
2020/04/14 Python
西班牙电子产品购物网站:Electronicamente
2018/07/26 全球购物
Crocs波兰官方商店:女鞋、男鞋、童鞋、洞洞鞋
2019/10/08 全球购物
临床医学大学生求职信
2013/09/28 职场文书
工业设计专业个人求职信范文
2013/12/28 职场文书
生日派对邀请函
2014/01/13 职场文书
幼儿园春游活动方案
2014/01/19 职场文书
团日活动策划书
2014/02/01 职场文书
2014年端午节活动方案
2014/03/11 职场文书
2014年乡镇卫生院工作总结
2014/11/24 职场文书
离婚协议书范本(通用篇)
2014/11/30 职场文书
汽车转让协议书
2015/01/29 职场文书
python实战之用emoji表情生成文字
2021/05/08 Python
Python办公自动化之Excel(中)
2021/05/24 Python