如何读取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)

相关推荐

  • 2021-2022-1 20211410周意凯《信息安全专业导论》第六周学习总结

    技术2021-2022-1 20211410周意凯《信息安全专业导论》第六周学习总结 2021-2022-1 20211410周意凯《信息安全专业导论》第六周学习总结学习目标Polya如何解决问题
    简单

    礼包 2021年11月1日
  • 怎样浅谈Java Servlets与CGI程序

    技术怎样浅谈Java Servlets与CGI程序怎样浅谈Java Servlets与CGI程序,相信很多没有经验的人对此束手无策,为此本文总结了问题出现的原因和解决方法,通过这篇文章希望你能解决这个问题。从网络的早期阶

    攻略 2021年11月20日
  • 抖音粉丝在哪里买,抖音活粉多少钱一个?

    技术抖音粉丝在哪里买,抖音活粉多少钱一个?抖音快速增长粉料的方法抖音无疑是目前新媒体中增长粉料最简单、增长最快的平台。从前,成都小甜甜一夜涨粉五百万,后来,灵魂当铺一天涨粉七十五万。这样的涨粉速度在其他平台是难以想象的。

    测评 2021年10月20日
  • JavaScript中的useRef 和 useState有什么用

    技术JavaScript中的useRef 和 useState有什么用小编给大家分享一下JavaScript中的useRef 和 useState有什么用,相信大部分人都还不怎么了解,因此分享这篇文章给大家参考一下,希望

    攻略 2021年11月25日
  • Hibernate在HQL中查询的方法有哪些

    技术Hibernate在HQL中查询的方法有哪些本篇内容介绍了“Hibernate在HQL中查询的方法有哪些”的有关知识,在实际案例的操作过程中,不少人都会遇到这样的困境,接下来就让小编带领大家学习一下如何处理这些情况吧

    攻略 2021年12月4日
  • 【9036期】JUC多线程—AQS抽象队列同步器原理

    技术【9036期】JUC多线程—AQS抽象队列同步器原理 【9036期】JUC多线程—AQS抽象队列同步器原理AQS 的工作原理
    什么是 AQS
    AQS,Abstract Queued Sync

    礼包 2021年12月23日