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)

到此这篇关于python使用cartopy库绘制台风路径代码的文章就介绍到这了,更多相关python cartopy绘制台风路径内容请搜索我们以前的文章或继续浏览下面的相关文章希望大家以后多多支持我们!

(0)

相关推荐

  • python Cartopy的基础使用详解

    前言 常用地图底图的绘制一般由Basemap或者cartopy模块完成,由于Basemap库是基于python2开发的一个模块,目前已经不开发维护.故简单介绍cartopy模块的一些基础操作. 一.基础介绍 首先导入相关模块. import numpy as np import matplotlib.pyplot as plt import cartopy.crs as ccrs import cartopy.feature as cfeature from cartopy.mpl.ticker

  • 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) #得到文件夹下的

  • Python利用turtle库绘制彩虹代码示例

    语言:Python IDE:Python.IDE 需求 做出彩虹效果 颜色空间 RGB模型:光的三原色,共同决定色相 HSB/HSV模型:H色彩,S深浅,B饱和度,H决定色相 需要将HSB模型转换为RGB模型 代码示例: #-*- coding:utf-8 –*- from turtle import * def HSB2RGB(hues): hues = hues * 3.59 #100转成359范围 rgb=[0.0,0.0,0.0] i = int(hues/60)%6 f = hues/

  • Python使用turtle库绘制小猪佩奇(实例代码)

    turtle(海龟)是Python重要的标准库之一,它能够进行基本的图形绘制.turtle图形绘制的概念诞生于1969年,成功应用于LOGO编程语言. turtle库绘制图形有一个基本框架:一个小海龟在坐标系中爬行,其爬行轨迹形成了绘制图形.刚开始绘制时,小海龟位于画布正中央,此处坐标为(0,0),前进方向为水平右方. 在Python3系列版本安装目录的Lib文件夹下可以找到turtle.py文件. 下面通过代码给大家介绍Python使用turtle库绘制小猪佩奇, 具体代码如下所示: # -*

  • python+matplotlib实现动态绘制图片实例代码(交互式绘图)

    本文研究的主要是python+matplotlib实现动态绘制图片(交互式绘图)的相关内容,具体介绍和实现代码如下所示. 最近在研究动态障碍物避障算法,在Python语言进行算法仿真时需要实时显示障碍物和运动物的当前位置和轨迹,利用Anaconda的Python打包集合,在Spyder中使用Python3.5语言和matplotlib实现路径的动态显示和交互式绘图(和Matlab功能类似). Anaconda是一个用于科学计算的Python发行版,支持 Linux, Mac, Windows系统

  • python使用turtle库绘制树

    本文实例为大家分享了python使用turtle库绘制树的具体代码,供大家参考,具体内容如下 # -*- coding: utf-8 -*- """ Spyder Editor This is a temporary script file. """ import turtle, datetime def drawGap(): #绘制数码管间隔 turtle.penup() turtle.fd(5) def drawLine(draw): #绘制

  • Python使用Turtle库绘制一棵西兰花

    Turtle库是Python中一个强大的绘制图像的函数库,灵活使用Turtle库可以绘制各种好看的图像. 下面介绍使用Turtle库绘制一棵西兰花. 绘制一棵西兰花,从主干出发以一定的角度向左向右生成对称的枝干,再从每个枝干出发向左向右生成对称的枝干,循环此动作,并最终绘制出一棵漂亮的西兰花. 首先导入Turtle库,并设置画笔大小.画笔速度及颜色,并隐藏画笔 from turtle import Turtle p=Turtle() p.pensize(5) p.color(clr) p.hid

  • python 利用turtle库绘制笑脸和哭脸的例子

    我就废话不多说了,直接上代码吧! import turtle turtle.pensize(5) turtle.pencolor("yellow") turtle.fillcolor("red") turtle.penup() turtle.goto(0,-200) turtle.pendown() turtle.circle(200) turtle.penup() turtle.goto(-100,50) turtle.pendown() turtle.begin

  • 使用Python的Turtle库绘制森林的实例

    这是由一个小作业引发的对Python的Turtle库的学习 下面是官方手册: Turtle官方手册 1.配置编程环境 由于现在的笔记本是临时借的,编程环境不是熟悉的环境,又由于种种原因没有安装成功Anaconda,就尝试了下其他的IDE: 最早接触的Enthought Canopy跑示例程序时各种报错无法解决(Python Kernal Crashed): 最著名的Python IDE是JetBeans的Pycharm,装好以后啥都没跑就占了1G内存(虽然舍友电脑上的没问题): 好在之前装了No

  • python使用turtle库绘制奥运五环

    Turtle库是Python语言中一个很流行的绘制图像的函数库,想象一个小乌龟,在一个横轴为x.纵轴为y的坐标系原点,(0,0)位置开始,它根据一组函数指令的控制,在这个平面坐标系中移动,从而在它爬行的路径上绘制了图形. 效果图: #奥运五环 import turtle turtle.setup(1.0,1.0) #设置窗口大小 turtle.title("奥运五环") #蓝圆 turtle.penup() turtle.right(90) turtle.forward(-50) tu

  • Python利用Turtle库绘制一个小老虎

    目录 导语 1.定义库以及初始化界面 2.画出左右两只耳朵 3.画出小老虎头部轮廓 4. 画出老虎的两只眼睛 5.画出老虎的鼻子和嘴巴 6.画出小老虎的左右肢体和脚趾 7.在需要的位置写上我们的新年祝福 8. 显示倒数3,2,1 9.显示我们需要的文字 10.设定代码运行入口,调用目标函数 成果展示 导语 哈喽铁汁们好久不见吖~小编已经复工了于是马不停蹄赶来给大家准备新年礼物算开工礼物吧! 海龟来作图 虎年就是要画老虎 2022不用纸和笔~ 今晚画老虎~ 一二三四五 老虎宝宝示意图 虎年怎么能少

随机推荐