如何读取netcdf数据并在matplotlib Basemap上绘图

技术如何读取netcdf数据并在matplotlib Basemap上绘图这篇文章主要为大家展示了“如何读取netcdf数据并在matplotlib Basemap上绘图”,内容简而易懂,条理清晰,希望能够帮助大家解决疑

本文主要向您展示如何读取netcdf数据并在matplotlib底图上绘制。内容简单易懂,条理清晰。希望能帮你解开疑惑。让边肖带领您学习如何阅读netcdf数据并绘制matplotlib底图。

Netcdf是一种数据格式;

Matplotlib是python中的一个绘图插件。

使用python进行地图效果是工作的需要。气象数据一般存储在netcdf中,python语言还没有深入研究,但已经按照要求达到了实现效果。

下面的代码详细列出了如何从netcdf读取数据以及在底图地图上绘制数据的效果。

Python的开发环境和各种轮子在下载代码,请百度搜索部署,我不太熟悉。下面的代码是基于python2.7开发的,netcdf的样例数据也可以从网上下载。

#设置编码格式,参考各种需要的类库。

#coding=utf-8

将matplotlib.pyplot导入为plt

将numpy作为np导入

导入scipy.ndimage

来自mpl_toolkits .底图导入底图,cm

将netCDF4导入为nc

将matplotlib.colors作为颜色导入

使用反斜杠的#netcdf文件的路径设置不同于windows常用的路径设置。

filename=' e :/git/grib2json/20171011/20171011-3KM _ west . NC '

#打开nc文件,从中根据变量获取数据,详细学习nc数据的格式。

数据=nc。数据集(文件名)

lat=data . variables[' XLAT '][:]。挤压()

lon=data . variables[' XLONG '][:]。挤压()

pm25=data . variables[' pm25 '][:]。挤压()

cutpm=pm25[20,4,]

#初始化底图对象,将坐标系设置为纬度和经度,然后是范围。

m=底图(epsg=4326,分辨率='i ',llcrnrlat=23,urcrnrlat=37,llcrnrlon=97,urcrnrlon=112)

图1=PLT . fig()

x,y=m(lon,lat)

#初始化颜色级别,根据不同的值渲染不同的颜色。

级别=[0,35,75,115,150,250,350,500]

myRBB=np.array([[153,219,231],[26,161,206],[0,123,217],[103,255,136],[10,227,50],[50,244,62],[228,243,62],[226,151,62]])

cmap=颜色。ListedColormap(myRBB/255.0)

norm=颜色。边界范数(边界=级别,ncolors=len(级别))

#根据上面的x\y\值数组在地图上绘制。

CS2=m.contourf(x,y,cutpm,cmap=cmap,alpha=0.75)

#画海岸线和国界

m . draw海岸线(线宽=0.5)

m.drawcounties(线宽=0.3)

plt.title('pm25 ',大小=20)

#加载底部附带的地图图像

m.etopo()

#显示绘图窗口

plt.show()

注意:使用时请去掉所有注释,很容易报错。以下是效果图:

如何读取netcdf数据并在matplotlib  Basemap上绘图

以上就是“如何读取netcdf数据并在matplotlib底图上绘制”一文的全部内容。感谢您的阅读!相信大家都有一定的了解,希望分享的内容对大家有所帮助。想了解更多知识,请关注行业资讯频道!

内容来源网络,如有侵权,联系删除,本文地址:https://www.230890.com/zhan/136913.html

(0)

相关推荐

  • oracle18c和12c的区别(oracle20c新特性)

    技术Oracle 数据库12c新特性有哪些这篇文章将为大家详细讲解有关Oracle 数据库12c新特性有哪些,小编觉得挺实用的,因此分享给大家做个参考,希望大家阅读完这篇文章后可以有所收获。  1. 在线重命名和重新定位

    攻略 2021年12月13日
  • HBase原理、基本概念、基本架构-3

    技术HBase原理、基本概念、基本架构-3 HBase原理、基本概念、基本架构-3HBase是Apache Hadoop的数据库,能够对大型数据提供随机、实时的读写访问。HBase的目标是存储并处理大型

    礼包 2021年12月22日
  • 经典Javascript正则表达式有哪些

    技术经典Javascript正则表达式有哪些本篇内容介绍了“经典Javascript正则表达式有哪些”的有关知识,在实际案例的操作过程中,不少人都会遇到这样的困境,接下来就让小编带领大家学习一下如何处理这些情况吧!希望大

    攻略 2021年11月26日
  • 喵喵喵是什么意思,猫往人身上蹭是什么意思

    技术喵喵喵是什么意思,猫往人身上蹭是什么意思猫咪蹭人可是有着很丰富的含义呢,各位铲屎官快点跟着喵大侠来了解一下吧喵喵喵是什么意思。1、如果猫咪用小脑袋的一侧、身体的一侧,或者是用尾巴对着你蹭来蹭去,这是给你打个招呼呢!一

    生活 2021年10月26日
  • AspNetCore应用注意哪些点

    技术AspNetCore应用注意哪些点今天就跟大家聊聊有关AspNetCore应用注意哪些点,可能很多人都不太了解,为了让大家更加了解,小编给大家总结了以下内容,希望大家根据这篇文章可以有所收获。背景已经有很多文章记录了

    攻略 2021年11月18日
  • Python使用消息框的语法显示相关消息方法是什么

    技术Python使用消息框的语法显示相关消息方法是什么这篇文章主要介绍“Python使用消息框的语法显示相关消息方法是什么”,在日常操作中,相信很多人在Python使用消息框的语法显示相关消息方法是什么问题上存在疑惑,小

    攻略 2021年11月25日