Python读取nc文件

最新推荐文章于 2024-08-08 20:17:02 发布

小徐爱吹晚风

最新推荐文章于 2024-08-08 20:17:02 发布

阅读量3.5k

点赞数 2

分类专栏： python学习新手笔记文章标签： python

本文链接：https://blog.csdn.net/ZmgandXhf/article/details/119255445

版权

python学习新手笔记专栏收录该内容

4 篇文章 1 订阅

订阅专栏

导入库

import xarray as xr
import matplotlib.pyplot as plt
import matplotlib

读取nc文件

file_path = r'D:XXX\XXX.nc' # nc文件路径
data = xr.open_mfdataset(file_path, parallel=True, decode_times=False) # 读取数据
data_name =(data.variables['data_name'])

查看数据类型

print(data_name)
# 输出：dask.array<getitem, shape=(720, 1440), dtype=float64, chunksize=(720, 1440), chunktype=numpy.ndarray>

二维数组切片观察部分区域

data_name2 = data_name[100:200,200:400]
print(data_name2)
# 输出：dask.array<getitem, shape=(100, 200), dtype=float64, chunksize=(100, 200), chunktype=numpy.ndarray>

数据可视化

fig = plt.figure(figsize = (20,10),dpi=100)
norm = matplotlib.colors.Normalize(vmin=-0.5, vmax=0.5)  
plt.rc('font', family='Microsoft YaHei')
h1 = plt.imshow(data_name2,cmap = 'seismic',norm = norm)  
c1 = plt.colorbar(h1) 
c1.ax.tick_params(labelsize=20) 
plt.contourf(data_name2,50,cmap = 'seismic',norm = norm) 
plt.title("XXX",fontsize = 20)
plt.xticks([0,50,100,150,200])
plt.yticks([0,25,50,75,100])
plt.tick_params(labelsize=20) 
plt.show()