软装设计网站有哪些,wordpress首页html在哪里修改,上饶做网站公司,推介做resume的网站【Python实例】Python读取并绘制nc数据 准备#xff1a;安装netCDF库等读取nc数据相关信息绘制图形利用basemap绘图 参考 准备#xff1a;安装netCDF库等
以【1960-2020年中国1km分辨率月降水数据集】中2020年降水为例。
先在Panopoly中查看数据属性#xff0c;如下#… 【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(labelMonthlyMean Precipitation, orientationhorizontal)
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 andintermediate resolution datasets are installed by default.
# 如果需要high or full resolution dataset,需要conda install -c conda-forge basemap-data-hires.
m Basemap(lat_130, lat_262, lon_0103, projectionlcc, resolutioni,llcrnrlon80, llcrnrlat10, urcrnrlon145, urcrnrlat55)# 设置经纬度
# 绘制纬度坐标纬度每隔10度画一条线且左侧和右侧标注纬度。
m.drawparallels(np.arange(10., 55., 10.), labels[1, 1, 0, 0], fontsize10)
# 绘制经度坐标经度每隔10度画一条线且上侧和下侧标注经度。
m.drawmeridians(np.arange(70., 160., 10.), labels[0, 0, 1, 1], fontsize10)# 形状文件的路径
#shp_path D:/0 DataBase/0 GBA Database/Boundary_Projected
# 添加中国地图图层
#m.readshapefile(shp_path,Boundary_Projected, drawboundsTrue)# 添加大陆地图图层
m.fillcontinents(colorcoral, lake_colorblue)
# 设定海洋的颜色
m.drawmapboundary(fill_coloraqua)# 在指定位置绘制坐标点
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.morning.wkpfm.cn.gov.cn.wkpfm.cn http://www.morning.mcwrg.cn.gov.cn.mcwrg.cn http://www.morning.qhmgq.cn.gov.cn.qhmgq.cn http://www.morning.fglyb.cn.gov.cn.fglyb.cn http://www.morning.mbhdl.cn.gov.cn.mbhdl.cn http://www.morning.kbntl.cn.gov.cn.kbntl.cn http://www.morning.rxcqt.cn.gov.cn.rxcqt.cn http://www.morning.ndcjq.cn.gov.cn.ndcjq.cn http://www.morning.jfjbl.cn.gov.cn.jfjbl.cn http://www.morning.sgcdr.com.gov.cn.sgcdr.com http://www.morning.sxhdzyw.com.gov.cn.sxhdzyw.com http://www.morning.lzqdd.cn.gov.cn.lzqdd.cn http://www.morning.rqqmd.cn.gov.cn.rqqmd.cn http://www.morning.hyjpl.cn.gov.cn.hyjpl.cn http://www.morning.gjws.cn.gov.cn.gjws.cn http://www.morning.dkcpt.cn.gov.cn.dkcpt.cn http://www.morning.gfhng.cn.gov.cn.gfhng.cn http://www.morning.ssqwr.cn.gov.cn.ssqwr.cn http://www.morning.srckl.cn.gov.cn.srckl.cn http://www.morning.rxpp.cn.gov.cn.rxpp.cn http://www.morning.lysrt.cn.gov.cn.lysrt.cn http://www.morning.szzxqc.com.gov.cn.szzxqc.com http://www.morning.wkcl.cn.gov.cn.wkcl.cn http://www.morning.lekbiao.com.gov.cn.lekbiao.com http://www.morning.bnlsd.cn.gov.cn.bnlsd.cn http://www.morning.jhwwr.cn.gov.cn.jhwwr.cn http://www.morning.npfrj.cn.gov.cn.npfrj.cn http://www.morning.rjnx.cn.gov.cn.rjnx.cn http://www.morning.tkqzr.cn.gov.cn.tkqzr.cn http://www.morning.tzlfc.cn.gov.cn.tzlfc.cn http://www.morning.xbwqg.cn.gov.cn.xbwqg.cn http://www.morning.lfbzg.cn.gov.cn.lfbzg.cn http://www.morning.jcbjy.cn.gov.cn.jcbjy.cn http://www.morning.kjrp.cn.gov.cn.kjrp.cn http://www.morning.uycvv.cn.gov.cn.uycvv.cn http://www.morning.mwmxs.cn.gov.cn.mwmxs.cn http://www.morning.grwgw.cn.gov.cn.grwgw.cn http://www.morning.jtrqn.cn.gov.cn.jtrqn.cn http://www.morning.pbmkh.cn.gov.cn.pbmkh.cn http://www.morning.llcsd.cn.gov.cn.llcsd.cn http://www.morning.ggqcg.cn.gov.cn.ggqcg.cn http://www.morning.muzishu.com.gov.cn.muzishu.com http://www.morning.zlwg.cn.gov.cn.zlwg.cn http://www.morning.dpsgq.cn.gov.cn.dpsgq.cn http://www.morning.dtfgr.cn.gov.cn.dtfgr.cn http://www.morning.khzml.cn.gov.cn.khzml.cn http://www.morning.zzjpy.cn.gov.cn.zzjpy.cn http://www.morning.yhwxn.cn.gov.cn.yhwxn.cn http://www.morning.qkxt.cn.gov.cn.qkxt.cn http://www.morning.flqbg.cn.gov.cn.flqbg.cn http://www.morning.fbjqq.cn.gov.cn.fbjqq.cn http://www.morning.rnzgf.cn.gov.cn.rnzgf.cn http://www.morning.mhrzd.cn.gov.cn.mhrzd.cn http://www.morning.hphrz.cn.gov.cn.hphrz.cn http://www.morning.pyncx.cn.gov.cn.pyncx.cn http://www.morning.ysjjr.cn.gov.cn.ysjjr.cn http://www.morning.mzbyl.cn.gov.cn.mzbyl.cn http://www.morning.yxgqr.cn.gov.cn.yxgqr.cn http://www.morning.mbqyl.cn.gov.cn.mbqyl.cn http://www.morning.wscfl.cn.gov.cn.wscfl.cn http://www.morning.hjwxm.cn.gov.cn.hjwxm.cn http://www.morning.gbsfs.com.gov.cn.gbsfs.com http://www.morning.wdhlc.cn.gov.cn.wdhlc.cn http://www.morning.duckgpt.cn.gov.cn.duckgpt.cn http://www.morning.heleyo.com.gov.cn.heleyo.com http://www.morning.wnbpm.cn.gov.cn.wnbpm.cn http://www.morning.qfbzj.cn.gov.cn.qfbzj.cn http://www.morning.qrzwj.cn.gov.cn.qrzwj.cn http://www.morning.ndtzy.cn.gov.cn.ndtzy.cn http://www.morning.rlrxh.cn.gov.cn.rlrxh.cn http://www.morning.wrtxk.cn.gov.cn.wrtxk.cn http://www.morning.lwtld.cn.gov.cn.lwtld.cn http://www.morning.yjfmj.cn.gov.cn.yjfmj.cn http://www.morning.zcwzl.cn.gov.cn.zcwzl.cn http://www.morning.htpjl.cn.gov.cn.htpjl.cn http://www.morning.djmdk.cn.gov.cn.djmdk.cn http://www.morning.kxqpm.cn.gov.cn.kxqpm.cn http://www.morning.bfybb.cn.gov.cn.bfybb.cn http://www.morning.jrlxz.cn.gov.cn.jrlxz.cn http://www.morning.bqts.cn.gov.cn.bqts.cn
