多个降雨nc数据python出图，横纵坐标为经纬度

飞飞飞01

已于 2023-04-26 19:31:48 修改

阅读量686

点赞数

文章标签： python matplotlib 信息可视化

于 2023-04-26 19:25:28 首次发布

本文链接：https://blog.csdn.net/weixin_47194510/article/details/130392858

版权

import netCDF4 as nc
import matplotlib.pyplot as plt
import numpy as np
import cv2

path1 = r'F:\bs\new2'
nc_10 = path1 + '\\'+"16fen3.nc"

path2 = r'F:\bs\new2'
nc_p = path2 + '\\'+"deepsd163.nc"

path3 = r'F:\bs\new2\20230404srcnn\results'
nc_psr=path3+'\\'+'save_cai.nc'

# 打开.nc文件
nf1 = nc.Dataset(nc_10, 'r')
# 获取变量
var1 = nf1.variables['ppt']
# 获取某一天的数据原始数据
nc_10_1 = var1[9, :, :].astype(float)
nc_10_1 = np.ma.filled(nc_10_1, np.nan)
nc_10_2 = np.flipud(nc_10_1)
data = nc_10_2 

#重采样为40km bicubic
# 计算缩放因子
scale_factor1 = 10 / 40
# 计算输出大小
output_shape = (int(data.shape[1] * scale_factor1), int(data.shape[0] * scale_factor1))
# 向下采样数据
nc_10_4 = cv2.resize(nc_10_2, output_shape, interpolation=cv2.INTER_CUBIC)

# 获取数据
scale_factor2 = 40 / 10
output_shape = (int(data.shape[1] * scale_factor2),