问题标签 [netcdf4]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
1086 浏览

python - 理解 pyresample 将不规则网格数据重新网格化为规则网格

我需要将不规则网格(朗伯圆锥形)上的数据重新网格化为规则网格。我认为 pyresample 是我最好的选择。事实上,我原来的纬度,经度不是一维的(似乎需要使用 basemap.interp 或 scipy.interpolate.griddata)。

我发现这个 SO 的答案很有帮助。但是我得到空的插值数据。我认为这与我的影响半径的选择以及我的数据被包装的事实有关(??)。

这是我的代码:

lat_n 和 lon_n 是不规则的,经纬度对应投影坐标 x,y。

由于 lon_n 的方式,我补充说:

所以现在如果我绘制它们,它们看起来不错并且还可以:

在此处输入图像描述

然后我创建了一组新的常规坐标:

按照上面的答案,我编写了以下代码:

从 pyresample.geometry 导入 SwathDefinition 从 pyresample.kd_tree 导入 resample_nearest

数据的分辨率大约是30km,所以我放了70km,我放的fill_value是数据中的那个,当然我也可以放0或者nan。

但是我得到一个空数组。

我做错了什么?另外-如果有另一种方法,我有兴趣知道它。Pyresample 文档有点薄,我需要更多帮助。

我确实发现这个答案建议使用另一个 griddata 函数:

似乎没问题:

在此处输入图像描述

但我想更多地了解 pyresample,因为它看起来很强大。

0 投票
1 回答
1565 浏览

python - 在 python xarray netcdf4 中创建输出变量和复制属性

我可以像这样在 netcdf4 中创建变量并复制属性:

xarray 对应的版本是什么?

0 投票
1 回答
1615 浏览

r - 将栅格名称从栅格堆栈导出到 R 中的 NetCDF 文件

我正在尝试从 netcdf 文件中提取栅格图层名称,就像之前从栅格堆栈中写入的那样。将栅格堆栈导出到 ncdf 可以正常工作。例如:

但是,当我将 ncdf 文件读回 R 时,不会保留 Z 维度(即年份)。例如:

所以我们得到的是波段编号,即1,2,3。但我需要的是由 Year(例如 2014,2015,2016)定义的文本,定义如下:

是否有可能做到这一点??这个问题并不新鲜,参考这里: https ://gis.stackexchange.com/questions/122167/export-band-names-with-netcdf-file-in-r

有一些复杂的解决方法可以得到所需的东西,但它们似乎效率很低(即将堆栈转换为矩阵,然后从这里操作它)。只是想知道是否有更优雅的方法,而无需编写大量额外代码并占用不必要的 RAM。

0 投票
2 回答
3541 浏览

python - 使用 xarray open_mfdataset 函数时出错

我正在尝试组合多个具有相同尺寸的 netCDF 文件,它们的尺寸如下:

但是,在使用 open_mfdataset 时,我收到此错误:

如何修复此错误?我的尺寸在所有文件中都相同

0 投票
1 回答
475 浏览

python - 尝试将 xarray 写入 netcdf 时出现“ValueError:chunksize 不能超过维度大小”

尝试将 xarray 对象写入 netcdf 文件时出现以下错误:

数据对我的记忆来说太大了,需要分块。
套路基本如下:

所以基本上,我只是操纵内容并重写。然而,大块似乎很糟糕。与重新分块到一个数组相同。我也不能重写ds。知道如何追踪错误或解决这个问题吗?

netCDF4 版本是 1.2.4
xarray(以前的 xray)版本是 0.8.2
dask 版本是 0.10.1

0 投票
3 回答
4506 浏览

python - 更改 netCDF 文件中的 _FillValue

是否有 python netCDF4 命令/示例来更改 netCDF 文件中的全局元数据 _FillValue?我尝试替换 netCDF 文件中的所有 -ve 值,但在设置 _FillValue 属性之前,这不起作用

0 投票
1 回答
98 浏览

python - 从数据中提取点的更简单方法?

我正在分析天气模型的输出,使用 python 模块 netCDF4 并使用 mpl_toolkits Basemap 制作图像。

我做了一些气压的等高线图,这很容易。现在我想为风添加倒钩,最后得到一个只有海风的情节。经过一番摆弄,我发现底图函数 fillcontinents() 实际上覆盖了这些风倒钩,即使稍后调用了 barbs 函数。因此,删除 fillcontinents() 函数后,所有的风倒钩都得到了想要的结果。

有没有办法填满大陆,但情节中仍然有我所有的倒钩?

我的代码看起来有点像这样(我遗漏了一些图形参数,它们对我的问题并不重要):

0 投票
0 回答
392 浏览

python - 安装 netcdf4-python 时出错

我正在尝试按照https://code.google.com/archive/p/netcdf4-python/wikis/UbuntuInstall.wiki安装 netcdf4-python 我已经成功安装了 HDF5,但是当我安装 netCDF4 时,在“ make ”命令中,收到以下错误

0 投票
3 回答
1379 浏览

python - 在 R 中创建/连接多维 NetCDF

将此用作来源:如何在 Windows 7 上使用 NCO 或 R 将每月 TRMM netCDF 文件连接到单个 netCDF 文件中?

我能够创建一个 rds 文件。但是,我真的很想将它保存为 nc 文件。我尝试通过以下方式创建一个具有 12 个步骤(每个月一个)时间维度的新 netCDF 文件:

现在的挑战是将降水数据添加到每个月。
在第一个脚本之后,我尝试使用该ncvar_put函数但没有成功。

我得到的

ncvar_put 错误(precip.nccreate, precip.ncvar, data2d, start = c(1), : object 'precip.nccreate' not found

nc_create 中的错误(“precip.nccreate.nc”,precip.ncvar,force_v4 = FALSE,:找不到对象 'precip.ncvar'

无论如何,我想我只是想找到一种简单的方法将多个 netcdf 文件连接成一个 netcdf。

谢谢

0 投票
2 回答
2570 浏览

r - 读取 R 中的 .nc4 文件时出错(ncdf4 包)

我正在尝试使用从NASA.nc4 files下载的数据集。

NCDF4 格式由该来源确认。我在 R 中使用下载 .file 来获取数据库,然后使用简单的 nc_open(ncdf4 包)来测试文件。不幸的是,结果是“未知文件格式”错误。

这是我的复制文件和我的脚本:

到目前为止一切都很好,但是在测试文件时:

我错过了什么吗?
谢谢你。