python使用cartopy库绘制台风路径代码
作者:oceanography-Rookie 发布时间:2023-08-22 15:50:18
标签:python,cartopy
使用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)
来源:https://blog.csdn.net/weixin_44237337/article/details/122889219
0
投稿
猜你喜欢
- 1、简介QListWidget 是一个升级版本的QListView , 它已经建立一个基于数据存储模型(QListWidgetItem),直
- 本文记录了mysql 8.0.17 安装配置的图文教程,供大家参考,具体内容如下1.下载安装包下载地址:安装包2.解压安装包到目录这里解压到
- 这篇文章主要介绍了python中的Elasticsearch操作汇总,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习
- 目录1. 理解进程2. 进程的语法3. join自定义进程类4. 守护进程总结1. 理解进程进程的概念:(process)进程就是正在运行的
- 误区 #1:在服务器故障转移后,正在运行的事务继续执行 这当然是错误的! 每次故障转移都伴随着某种形式的恢复。但是如果当正在执行的事务没有C
- Python工程师面试必备25条Python知识点:1.到底什么是Python?你可以在回答中与其他技术进行对比下面是一些关键点:Pytho
- 本文为大家分享的Python工程师面试题主要与Python Web相关,供大家参考,具体内容如下1、解释一下 WSGI 和 FastCGI
- 服务器重新启动的情况当mysqld启动时,所有的授权表内容被读进存储器并且从那时开始生效。被服务器立即应用的情况用GRANT、REVOKE或
- 很神奇的一个晚上,居然在以前老同事的群里跟同事讨论起CSS的东西来了,不过很意外的还是有收获。在IE中常常会碰到如果将容器定位后,出现容器内
- 使用正则库爬取淘宝商品的商品信息,首先我们需要确定想要爬取的对象我们在淘宝里搜索“python”,出来的结果从url连接中可以得到搜索商品的
- 在 python 中除了用 opencv,也可以用 matplotlib 和 PIL 这两个库操作图片。本人偏爱 matpoltlib,因为
- 最近为了测试网速情况怎么样,由于部分业务服务器需要关闭icmp,这样的话采用普通的ping就无法适应我的需求,于是自己简单的写了一个基于tc
- Python之Selenium(自动化浏览器测试)1.安装seleniumpip install selenium -i https://p
- Python中打开文本使用的是with语句,比如打开一个文件并读取每一行with open(filename) as fp: f
- 本人 python新手,使用的环境是python2.7,勿喷# -*- coding:utf8 -*-import random
- 功能与实现点击小图片可以查看大图实现就是把大图放置在顶层(z-index大于当前页面的),并且还可以加一些额外的比如透明度什么的。大图以动画
- 在使用操作XML文件时,我们可以使用Load方法直接加载xml文件即可,在ie和ff下通用。但是是XML字符串,则在两种浏览器下就会有所不同
- 一、sqlserver优点:易用性、适合分布式组织的可伸缩性、用于决策支持的数据仓库功能、与许多其他服务器软件紧密关联的集成性、良好的性价比
- 所谓的列表推导式,就是指的轻量级循环创建列表。基本使用方式# 创建一个0-10的列表a = [x for x in range(11)]pr
- 本文实例讲述了PHP常用函数之获取汉字首字母功能。分享给大家供大家参考,具体如下://获取汉字的首字母function getFirstCh