当前位置: 首页 > news >正文

sem推广培训百度seo发帖推广

sem推广培训,百度seo发帖推广,搬瓦工搭建wordpress,专门做搜索种子的网站有哪些【Python实例】Python读取并绘制nc数据 准备:安装netCDF库等读取nc数据相关信息绘制图形利用basemap绘图 参考 准备:安装netCDF库等 以【1960-2020年中国1km分辨率月降水数据集】中2020年降水为例。 先在Panopoly中查看数据属性,如下&#…

【Python实例】Python读取并绘制nc数据

  • 准备:安装netCDF库等
  • 读取nc数据相关信息
  • 绘制图形
    • 利用basemap绘图
  • 参考

准备:安装netCDF库等

以【1960-2020年中国1km分辨率月降水数据集】中2020年降水为例。

先在Panopoly中查看数据属性,如下:
在这里插入图片描述
右键绘制图形,如下:
在这里插入图片描述

读取nc数据相关信息

Python代码如下:

import netCDF4
from netCDF4 import Dataset
import matplotlib.pyplot as plt
import numpy as np# 打开下载的.nc文件
ncfile = Dataset("D:/0 DataBase/CN_Prec_MonthlyMean_1km_2020.nc/CN_Prec_MonthlyMean_1km_2020.nc")
print(ncfile.variables.keys())# 读取变量
lat = ncfile.variables['latitude'][:]
long = ncfile.variables['longitude'][:]
[X, Y] = np.meshgrid(lat, long)    # 合成矩阵,将向量latitude和longitude定义的区域转换成矩阵X和Y
pre = ncfile.variables['pre'][:]
time = ncfile.variables['time'][:]# 绘制图形
plt.contourf(Y, X,pre[0])
plt.colorbar(label="MonthlyMean Precipitation", orientation="horizontal")
plt.show()# 关闭nc文件
ncfile.close()

输出属性如下:
在这里插入图片描述

绘制图形

利用basemap绘图

Python代码如下:

#利用basemap绘图
import netCDF4
from mpl_toolkits import basemap
from netCDF4 import Dataset
from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import numpy as np# 打开下载的.nc文件
ncfile = Dataset("D:/0 DataBase/CN_Prec_MonthlyMean_1km_2020.nc/CN_Prec_MonthlyMean_1km_2020.nc")# 读取变量
lat = ncfile.variables['latitude'][:]
long = ncfile.variables['longitude'][:]
[X, Y] = np.meshgrid(lat, long)    # 合成矩阵,将向量latitude和longitude定义的区域转换成矩阵X和Y
pre = ncfile.variables['pre'][:]
time = ncfile.variables['time'][:]# 创建一个图形和一个轴
fig = plt.figure(figsize=(10, 8))
ax = fig.add_subplot(111)# 设定标题
titleName = "MonthlyMean Precipitation"# 初始化Basemap对象,设置地图参数
# 设置投影方式:cyl为圆柱投影、还可设置merc为mercator投影,lcc兰伯特投影等。
# 参数 "resolution" 用于控制地图面积边缘的精细程度,有 'crude', 'low' and'intermediate' resolution datasets are installed by default.
# 如果需要'high' or 'full' resolution dataset,需要`conda install -c conda-forge basemap-data-hires`.
m = Basemap(lat_1=30, lat_2=62, lon_0=103, projection='lcc', resolution='i',llcrnrlon=80, llcrnrlat=10, urcrnrlon=145, urcrnrlat=55)# 设置经纬度
# 绘制纬度坐标,纬度每隔10度画一条线,且左侧和右侧标注纬度。
m.drawparallels(np.arange(10., 55., 10.), labels=[1, 1, 0, 0], fontsize=10)
# 绘制经度坐标,经度每隔10度画一条线,且上侧和下侧标注经度。
m.drawmeridians(np.arange(70., 160., 10.), labels=[0, 0, 1, 1], fontsize=10)# 形状文件的路径
#shp_path = "D:/0 DataBase/0 GBA Database/Boundary_Projected"
# 添加中国地图图层
#m.readshapefile(shp_path,'Boundary_Projected', drawbounds=True)# 添加大陆地图图层
m.fillcontinents(color="coral", lake_color="blue")
# 设定海洋的颜色
m.drawmapboundary(fill_color="aqua")# 在指定位置绘制坐标点
lon_me = long.mean()
lat_me = lat.mean()
lon_do, lat_do = np.meshgrid(lon_me, lat_me)  # 生成网格点坐标矩阵# 转换经纬度到Basemap坐标系
xi, yi = m(lon_do, lat_do)
# print(xi,yi)
m.plot(xi, yi, "ro")  # 绘制红色点
# 在指定位置标注文字
plt.text(xi, yi, "A1(%5.1fW, %3.1fN)" % (lon_me, lat_me))# 关闭.nc文件
ncfile.close()# 显示图形
plt.show()

图形绘制如下:
在这里插入图片描述

参考

http://www.shuangfujiaoyu.com/news/42572.html

相关文章:

  • 完整的网站后台权限编码知乎软文推广
  • 免费一级做网站搜索引擎下载入口
  • 网站开启伪静态成都今天重大新闻事件
  • 滨海网站建设公司优化大师官方免费
  • 女装网站建设规划书自动提取关键词的软件
  • 网站设计一个月多少钱手机网站模板下载
  • 做网站有用nodejsseo站内优化教程
  • 做网站注册哪类商标关键词优化搜索引擎
  • 站内优化怎么做网络营销的特点有哪些特点
  • 南京网站制作招聘网地推扫码平台
  • 页面设计工作内容自述网站优化联系
  • 淘宝客服推销做网站的技巧seo是免费的吗
  • 专题探索网站开发教学模式的结构郑州网络优化实力乐云seo
  • 旅游网站需求分析百度推广登陆首页
  • 中国建设银行官方网站纪念币免费网站创建
  • 百度站点管理优化大师下载安装app
  • 企业网站图片渐变效果怎么做正安县网站seo优化排名
  • 公司做公司网站宣传宁波网站建设的公司
  • 做网站如何引流seo是什么职业做什么的
  • 百度地图网站后台更新竞价推广套户渠道商
  • 地方网站还有得做吗中国局势最新消息今天
  • 分类网站怎么做seo网络营销推广主要做什么
  • 北京做网站的公司拟杭州网络
  • 那可以做网站武汉网络关键词排名
  • wordpress 数学主题win10优化软件哪个好
  • 网站开发和运营合同分开签么北京网站建设公司哪家好
  • 纯div css做网站简洁版杭州网站设计公司
  • 重庆网站建设索q479185700公司官网开发制作
  • 泰安市人才服务平台芜湖seo
  • wordpress 邮件文本优化设计答案六年级