python - 读取 Netcdf 文件并用 Python 计算 RMSE

Question

我有一个 netcdf 文件。我在这个文件中有两个变量:wspd_wrf_m 和 wspd_sodar_o。我想读取 netcdf 文件并计算 wspd_wrf_m 和 wspd_sodar_o 之间的 RMSE 值。

变量的维度为(天、时间)，即 (1094, 24)我想计算最近 365 天的文件的 RMSE。你能帮我解决这个问题吗？

我知道我需要使用:

from netCDF4 import Dataset
import numpy as np

g = Dataset('station_test_new.nc','r',format='NETCDF3_64BIT')
wspd_wrf = g.variables["wspd_wrf_m"][:,:]
wspd_sodar = g.variables["wspd_sodar_o"][:,:]

但是如何选择我需要的过去 365 天的每小时数据并从中计算 RMSE？

Answer 1

选择过去 365 天是将数组切片为正确大小的问题。例如:

import numpy as np
var = np.zeros((1094, 24))
print(var.shape, var[729:,:].shape, var[-365:,:].shape)

打印:

(1094, 24) (365, 24) (365, 24)

因此 var[729:,:] 和 var[-365:,:] 都将过去 365 天(包含所有每小时值)从 1094 天中分割出来日大小的数组。

Numpy manual中有更多信息/更多示例.

有很多关于如何在 Python 中计算 RMSE 的示例(例如 this 一个)。请尝试一下，如果您无法使其工作，请根据您的尝试更新您的问题。