python使用cartopy库绘制台风路径代码
更新时间:2022年2月13日 17:06 点击:277 作者:oceanography-Rookie
使用python基于cartopy库绘制台风路径
使用python 绘制西太平洋进入我国的台风路径,文件为.dat
格式,内容如下所示:
代码如下:
import netCDF4 as nc4 import matplotlib.pyplot as plt import numpy as np import datetime import os import cartopy.crs as ccrs path='E://' #文件路径 files= os.listdir(path) #得到文件夹下的所有文件名称 fig=plt.figure(figsize=(20,12)) #设置画布大小 parallels = np.arange(0.,90.,3.) meridians = np.arange(0.0,360.,3.) ax = plt.axes(projection=ccrs.Robinson()) #设置投影方式 # Set figure extent & ticks ax.set_extent([100, 150, 5, 50]) #设置纬度范围 # plt.grid(linestyle=':',color='y')# for file in files: #按照顺序在 files 里面进行每一个文件的 数据名称 循环读取 f = open(path+file,'r') # 打开第一个 dat 文件 records = f.readlines() # 读取这个文件里面的所有数据 f.close() # 关闭这个dat文件 date_t = '' # 设置一个用来表示空格 btk_lat = [] # 设置一个空的 待传入数据的纬度 btk_lon = [] # 设置一个空的 待传入数据的经度 btk_vmax = [] # 风速最大值 Maximum sustained wind speed in knots: 0 - 300 kts. btk_time = [] # 时间 btk_rmw = [] # 最大风速半径 radius of max winds, 0 - 999 n mi. btk_name = [] # 台风名称 for rcd in records: # 对这个dat文件里面,已经读取的每一行数据进行循环处理 strs = rcd.split(',') #将每一个','分隔开 if(len(strs)<21): #判断语句,如果这个被分割开的字符 长度<21 ,继续进行处理 continue date_str = strs[2].strip(' ') #将strs这个list的索引为2的值赋给data_str,既年月日时 if date_str == date_t:#判读如果是一个空格值,赋给data——str continue dt = datetime.datetime(int(date_str[0:4]),int(date_str[4:6]),int(date_str[6:8]),\ int(date_str[8:]),0,0,0) btk_time.append(nc4.date2num(dt,units='second since 1970-1-1 00:00:00'))#计算距离给的时间有多少秒,并从后往前排列 #处理纬度 lat_str = strs[6].strip() #判断南北纬 if lat_str[-1] == 'N': lat_t = float(lat_str[0:-1])*0.1 else: lat_t = float(lat_str[0:-1])*-0.1 btk_lat.append(lat_t) #处理经度 lon_str = strs[7].strip() #判断 东西经 if lon_str[-1] == 'E': lon_t = float(lon_str[0:-1])*0.1 else: lon_t = float(lon_str[0:-1])*-0.1 btk_lon.append(lon_t) #处理最大风速 vmax = strs[8].strip() btk_vmax.append(float(vmax))#转换为单浮点型,(带小数点) #时间 date_t = date_str #最大风速半径 rmw = strs[19].strip() btk_rmw.append(float(rmw)) #处理台风名称 if(len(strs) < 27): btk_name.append('noname') else: name = strs[27].strip() btk_name.append(name) #============================================================================== btk_lat = np.array(btk_lat) #将得到的list 值转换为数组型的值,为了便于绘图。因为绘图的横纵坐标都是数组排列 btk_lon = np.array(btk_lon)%360 #因为原始经度为-180 - 0 -180 ,出现断隔,为解决问题,化为 0-360 btk_time = np.array(btk_time) #时间转换 btk_vmax = np.array(btk_vmax)*0.5144 #风速换算公式 btk_rmw = np.array(btk_rmw)*1.852 # #判断,如果数组纬度的值是0,则为nan值,既无法计算的值(无穷大,,),否则即为台风的名称 if(len(btk_lat) == 0): tc_name = 'noname' else: index = btk_vmax.argmax() tc_name = btk_name[index] #进行绘图,经度、纬度曲线 ax.plot(btk_lon,btk_lat,color='k',linewidth=0.5,transform=ccrs.PlateCarree()) #散点图绘制,经度、纬度、最大风速, cb = ax.scatter(btk_lon,btk_lat,c=btk_vmax,s=10.0,transform=ccrs.PlateCarree() ,vmin=10,vmax=60) ax.coastlines() ax.gridlines(draw_labels=True, dms=True, x_inline=False, y_inline=False) plt.colorbar(cb,label='Vmax (m/s)',pad=0.07,orientation='vertical',shrink=1) plt.title(' path') # 保存绘制图片 ,注意保存路径不能放在dat文件夹中 #fig.savefig(path2+'tester.tiff',format='tiff',dpi=100)
到此这篇关于python使用cartopy库绘制台风路径代码的文章就介绍到这了,更多相关python cartopy绘制台风路径内容请搜索猪先飞以前的文章或继续浏览下面的相关文章希望大家以后多多支持猪先飞!
原文出处:https://blog.csdn.net/weixin_44237337/article/details/122889
上一篇: Python装饰器实现函数运行时间的计算
下一篇: Python的语法基础你真的了解吗
相关文章
- 这篇文章主要介绍了python-opencv-画外接矩形框的实例代码,代码简单易懂,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友可以参考下...2021-09-04
Python astype(np.float)函数使用方法解析
这篇文章主要介绍了Python astype(np.float)函数使用方法解析,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下...2020-06-08- 2022虎年新年即将来临,小编为大家带来了一个利用Python编写的虎年烟花特效,堪称全网最绚烂,文中的示例代码简洁易懂,感兴趣的同学可以动手试一试...2022-02-14
- 在本篇文章里小编给大家分享的是一篇关于python中numpy.empty()函数实例讲解内容,对此有兴趣的朋友们可以学习下。...2021-02-06
python-for x in range的用法(注意要点、细节)
这篇文章主要介绍了python-for x in range的用法,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧...2021-05-10- 这篇文章主要介绍了Python 图片转数组,二进制互转操作,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧...2021-03-09
- 这篇文章主要介绍了Python中的imread()函数用法说明,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧...2021-03-16
- 这篇文章主要介绍了python如何实现b站直播自动发送弹幕,帮助大家更好的理解和学习使用python,感兴趣的朋友可以了解下...2021-02-20
python Matplotlib基础--如何添加文本和标注
这篇文章主要介绍了python Matplotlib基础--如何添加文本和标注,帮助大家更好的利用Matplotlib绘制图表,感兴趣的朋友可以了解下...2021-01-26- 这篇文章主要介绍了解决python 使用openpyxl读写大文件的坑,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧...2021-03-13
- 今天小编就为大家分享一篇python 计算方位角实例(根据两点的坐标计算),具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧...2020-04-27
- 这篇文章主要为大家详细介绍了python实现双色球随机选号,文中示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下...2020-05-02
- 在本篇文章里小编给大家整理的是一篇关于python中使用np.delete()的实例方法,对此有兴趣的朋友们可以学习参考下。...2021-02-01
- 这篇文章主要介绍了使用Python的pencolor函数实现渐变色功能,本文通过实例代码给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友可以参考下...2021-03-09
- 这篇文章主要介绍了python自动化办公操作PPT的实现,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧...2021-02-05
Python getsizeof()和getsize()区分详解
这篇文章主要介绍了Python getsizeof()和getsize()区分详解,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧...2020-11-20- 这篇文章主要为大家详细介绍了python实现学生通讯录管理系统,文中示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下...2021-02-25
- 这篇文章主要介绍了PyTorch一小时掌握之迁移学习篇,本文给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友可以参考下...2021-09-08
- 这篇文章主要介绍了解决python 两个时间戳相减出现结果错误的问题,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧...2021-03-12
- 这篇文章主要介绍了Python绘制的爱心树与表白代码,本文通过实例代码给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友可以参考下...2021-04-06