1. 本文原创链接:https:////www.cnblogs.com/zhanling/p/12192978.html
  2.  
  3. 1 import numpy as np
  4. import xarray as xr
  5. import cartopy.crs as ccrs
  6. import cartopy.feature as cfeat
  7. from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER
  8. import matplotlib.pyplot as plt
  9.  
  10. ds = xr.open_dataset('2039071310.003.nc')
  11. t = ds['value']
  12. lons = ds.lon.data
  13. lats = ds.lat.data
  14. temp = xr.DataArray(t.data.T, coords=[lats,lons], dims=['latitude','longitude'])
  15.  
  16. # 创建画图空间
  17. proj = ccrs.PlateCarree() #创建投影
  18. fig = plt.figure(figsize=(16,9)) #创建页面
  19. ax = fig.subplots(1, 1, subplot_kw={'projection': proj}) #子图
  20. # 设置地图属性:加载国界、海岸线、河流、湖泊
  21. ax.add_feature(cfeat.BORDERS.with_scale('50m'), linewidth=0.8, zorder=1)
  22. ax.add_feature(cfeat.COASTLINE.with_scale('50m'), linewidth=0.6, zorder=1)
  23. ax.add_feature(cfeat.RIVERS.with_scale('50m'), zorder=1)
  24. ax.add_feature(cfeat.LAKES.with_scale('50m'), zorder=1)
  25. # 设置网格点属性
  26. gl = ax.gridlines(crs=ccrs.PlateCarree(), draw_labels=True,
  27. linewidth=1.2, color='k', alpha=0.5, linestyle='--')
  28. gl.xlabels_top = False #关闭顶端标签
  29. gl.ylabels_right = False #关闭右侧标签
  30. gl.xformatter = LONGITUDE_FORMATTER #x轴设为经度格式
  31. gl.yformatter = LATITUDE_FORMATTER #y轴设为纬度格式
  32. # 设置colorbar
  33. cbar_kwargs = {
  34. 'orientation': 'horizontal',
  35. 'label': 'Potential',
  36. 'shrink': 0.8,
  37. }
  38. # 画图
  39. levels = np.arange(0,1,0.1)
  40. temp.plot.contourf(ax=ax, levels=levels, cmap='Spectral_r',
  41. cbar_kwargs=cbar_kwargs, transform=ccrs.PlateCarree())
  42. plt.savefig('test.jpg')

示例效果(强对流概率预报结果):

 

Python-气象-大气科学-可视化绘图系列(一)——利用xarray读取netCDF文件并画图(代码+示例)的更多相关文章

  1. Python-气象-大气科学-可视化绘图系列(二)——利用basemap叠加地图,并添加白化效果(代码+示例)

    本文为原创链接: https:////www.cnblogs.com/zhanling/p/12193031.html 白化单图代码: import numpy as np import xarray ...

  2. Python-气象-大气科学-可视化绘图系列(三)—— 地图上自动标注省会名称(demo调整中)(代码+示例)

    本文为原创文章 本文链接:https://www.cnblogs.com/zhanling/p/12606990.html # -*- coding: utf-8 -*- ''' Author: He ...

  3. 使用Python处理Excel文件的一些代码示例

    笔记:使用Python处理Excel文件的一些代码示例,以下代码来自于<Python数据分析基础>一书,有删改 #!/usr/bin/env python3 # 导入读取Excel文件的库 ...

  4. 使用Python处理CSV文件的一些代码示例

    笔记:使用Python处理CSV文件的一些代码示例,来自于<Python数据分析基础>一书,有删改 # 读写CSV文件,不使用CSV模块,仅使用基础Python # 20181110 wa ...

  5. python利用xlrd读取excel文件始终报错原因

    1.代码按照网上百度的格式进行书写如下: 但运行后,始终报错如下: 百度了xlrd网页: 分明支持xls和xlsx两种格式的文件,但运行始终报错. 最后找到原因是因为我所读取的文件虽然是以.xls命名 ...

  6. 【 D3.js 进阶系列 — 1.2 】 读取 CSV 文件时乱码的解决方法

    在 D3 中使用 d3.csv 读取 CSV 文件时,有时会出现乱码问题. 怎么解决呢? 1. 乱码问题 使用 d3.csv 读取 xxx.csv 文件时.假设 xxx.csv 文件使用的是 UTF- ...

  7. perl 利用管道读取压缩文件内容

    perl的文件句柄不仅支持普通文件, 还支持管道,今天需要统计一个fastq文件中的序列数和碱基数,而NGS的fastq文件一般都是gzip压缩的,所以 需要读取压缩文件中的内容,代码如下: my ( ...

  8. python常识系列07-->python利用xlwt写入excel文件

    前言 读书之法,在循序而渐进,熟读而精思.--朱熹 抽空又来写一篇,毕竟知识在于分享! 一.xlwt模块是什么 python第三方工具包,用于往excel中写入数据:(ps:只能创建新表格,不能修改表 ...

  9. 利用XPath读取Xml文件

    之所以要引入XPath的概念,目的就是为了在匹配XML文档结构树时能够准确地找到某一个节点元素.可以把XPath比作文件管理路径:通过文件管理路 径,可以按照一定的规则查找到所需要的文件:同样,依据X ...

随机推荐

  1. 参加Folding@Home(FAH)项目,为战胜新冠肺炎贡献出自己的一份力量

    鉴于新冠病毒(COVID-19)在全球范围内的大规模传播,PCMR和NVIDIA呼吁全球PC用户加入Folding@home项目贡献自己闲置的GPU计算力,协助抗击新冠状病毒疫情. 目前全球有超过40 ...

  2. C# 获取基类或者接口的所有继承类方法

    static class ReflectionHelper { public static IEnumerable<T> CreateAllInstancesOf<T>() { ...

  3. 强化学习之七:Visualizing an Agent’s Thoughts and Actions

    本文是对Arthur Juliani在Medium平台发布的强化学习系列教程的个人中文翻译,该翻译是基于个人分享知识的目的进行的,欢迎交流!(This article is my personal t ...

  4. TensorFlow 中文资源精选,官方网站,安装教程,入门教程,实战项目,学习路径。

    Awesome-TensorFlow-Chinese TensorFlow 中文资源全集,学习路径推荐: 官方网站,初步了解. 安装教程,安装之后跑起来. 入门教程,简单的模型学习和运行. 实战项目, ...

  5. if-else、switch、while、for

    文章主要会涉及如下几个问题: if-else 和 switch-case 两者相比谁的效率会高些?在日常开发中该如何抉择? 如何基于赫夫曼树结构减少 if-else 分支判断次数? 如何巧妙的应用 d ...

  6. sql MySQL5.7 安装 centos docker

    MySQL5.7 安装 1 . 普通安装 # 准备工作 停止以前的数据库并删除配置文件 systemctl stop mysqld rpm -e mysql-community-server-5.7. ...

  7. CSS样式的4种写法 | 以及选择器的几种用法

    CSS样式: 1.内部样式表 <style type="text/css">    样式表写法   </style> 2.使用link标签,在文档中声明使用 ...

  8. 各种杂记关于Linux

    修改Linux 日期 修改Linux时间

  9. eclipse导入项目乱码问题及快键键收集

    1.eclipse导入的项目出现乱码,选择该项目,右键选择最下方的Properties,第一个出现的Resource页面: 在Text file encoding项中,选择other. 目前中文编码主 ...

  10. 手动搭建I/O网络通信框架2:Socket和ServerSocket入门实战,实现单聊

    第一章:手动搭建I/O网络通信框架1:Socket和ServerSocket入门实战,实现单聊 在第一章中运用Socket和ServerSocket简单的实现了网络通信.这一章,利用BIO编程模型进行 ...