netcdf4相关内容

如何使用CDO在NetCDF文件上应用海掩码

假设我创建了一个netCDF文件,其中有很多零。我希望对数据应用海洋掩码,以便只对海洋中的区域进行真正的掩码,而保留陆地上的数据。 我的数据在陆地上有很多零(这是正确的),但在海中也有很多零(这是不正确的)。 我本可以使用cdo setmissval,nan input.nc output.nc,但这也会将土地上的值更改为NaN。 有人有什么(好的)解决方案吗? 推荐答案 ..
发布时间:2022-08-20 23:30:40 其他开发

R在打开netcdf文件时崩溃

我下载了一个netcdf文件,并尝试在R中打开它。以下是我的代码 download.file("https://data.giss.nasa.gov/impacts/agmipcf/agmerra/AgMERRA_1980_prate.nc4",destfile = "AgMERRA_1980_prate.nc4", method="libcurl") 我要使用R 打开netcdf文件 ..
发布时间:2022-05-13 19:21:26 其他开发

使用 Python 将 NetCDF 文件转换为 CSV 或文本

我正在尝试使用 Python 将 netCDF 文件转换为 CSV 或文本文件.我已阅读这篇文章 但我仍然缺少一步(我是 Python 新手).它是一个数据集,包括纬度、经度、时间和降水数据. 这是我目前的代码: 导入 netCDF4将熊猫导入为 pdprecip_nc_file = 'file_path'nc = netCDF4.Dataset(precip_nc_file, mode= ..
发布时间:2021-12-22 20:08:53 Python

子集 xarray.Dataset 关于多个坐标

假设我使用 xarray.open_dataset(..., decode_times=False) 加载了一个 xarray.Dataset 对象,打印时看起来像这样: 尺寸:(边界:2,纬度:15,经度:34,列维:8,时间:3650)坐标:* 时间(时间) float64 3.322e+04 3.322e+04 3.322e+04 3.322e+04 ...* plev (plev) fl ..
发布时间:2021-09-01 18:32:17 Python

迭代一个函数以将许多栅格堆栈合并为一个

在这个问题上坚持了一段时间.到处寻找答案,但我似乎在 Stack 上找不到任何东西.大家可以提供的任何帮助将不胜感激. 我的主要问题是我需要导入很多很多 netcdf4 文件,为每个文件创建光栅砖,然后将许多砖组合起来为每个变量制作一个“主砖".给你一个更清楚的例子,我有 40 年(netcdf = 40)的许多气候变量(n = 15),它们每天都在解析.目标是每月汇总,但首先我必须得到这个 ..
发布时间:2021-07-02 20:30:55 其他开发

使用气候数据运营商 (CDO) 的每日数据的月雨天总和

我有具有每日时间分辨率的气候数据,并且想要按月和按年计算有降水(例如,大于 1 毫米/天)的天数. 我尝试过 eca_pd,1 和 eca_rr1,但这些命令返回所有年份的雨天总数. 例如,cdo eca_pd,1 infile outfile 是否有返回每个月和/或年雨天的命令? 解决方案 您可以使用 CDO 的屏蔽功能完成此任务. 第一步是制作一个等效文件,如果 ..
发布时间:2021-06-07 19:17:17 其他开发

使用 r 中的栅格包聚合季节性均值

我正在尝试将每日数据(35 年)汇总到每月,然后使用 R 中的栅格包计算季节性平均值(我知道如何使用 CDO 进行计算).下面是我的代码,它输出所有年份的 4 个季节性平均值(140 层).如何循环输出仅 4 层(4 个季节)?.我感谢您的帮助. dailydata ..
发布时间:2021-06-07 19:17:11 其他开发

将 netcdf 时间变量转换为 R 日期对象

我有一个带有时间序列的 netcdf 文件,时间变量具有以下典型元数据: 双倍时间(时间);时间:标准名称=“时间";时间:边界=“time_bnds";time:units = "自 1979-1-1 00:00:00 以来的天数";时间:日历=“标准";时间:轴=“T"; 在 R 中,我想将时间转换为 R 日期对象.我目前通过读取单位属性并拆分字符串并使用第三个条目作为我的原点(因此假设 ..
发布时间:2021-06-07 19:16:52 其他开发

Jetified-aws-android-sdk和jetified-aws-java-sdk中的重复类

我有一个Android应用程序,该应用程序明确使用最新版本(v2.15.+)的Amazon AWS库,但它也依赖于一个似乎庞大的库(netcdf4),该库似乎是自己的(旧的))复制(v1.11.236)亚马逊的非Android Java库,并在构建时导致“重复类"错误. 我知道我需要以某种方式修改build.gradle,以便它仅将旧版本(v1.11.236)用于实现实现组:'edu.uca ..

使用并行NetCDF保存分布式3D复杂阵列

我有一个用Fortran编写的基于MPI的程序,该程序在每个节点(2D时间序列的各个部分)上生成复杂数据的3D数组.我想使用并行I/O将这些数组写入单个文件,可以在python中相对容易地打开该文件以进行进一步的分析/可视化.理想情况下,我希望该解决方案能够提高内存效率(即避免创建中间临时数组). 使用NetCDF,我设法调整了子例程可以实现3D实数数组.但是,对于复杂数组,我遇到了绊脚石. ..
发布时间:2021-05-09 19:17:45 Python

如何使用虹膜模块绘制大气的垂直剖面以及地形?

我有一个netcdf,其风速处于模型级别.在同一netcdf上,我具有每个模型级别的高度.我将netcdf转换为一个立方体,因此每个级别的高度都变成了辅助坐标.我想绘制一个横截面(经度x经度),并希望模型级别遵循地形.我尝试使用Iris模块文档示例( https://scitools.org.uk/iris/docs/latest/examples/General/cross_section.ht ..

更改netCDF文件的网格大小

让我们假设我有2个netCDF数据文件,这些文件具有相同区域(例如,南美,非洲等)的数据,但不同网格的大小分别为0.5度x 0.5度和1.0度x 1.0度.我想将其网格大小增加或减少为其他值,例如0.25 x 0.25或1.0 x 1.0,以便可以轻松地将其用于栅格计算和比较等. 是否有一种使用任何bash脚本,CDO等进行此操作的方法. 可以从此处下载示例数据. https://ww ..
发布时间:2021-04-14 20:05:20 其他开发

H5PY不坚持分块规范?

问题: 我有现有的netCDF4文件(其中约有5000个)(通常形状为96x3712x3712)数据点(float32).这些文件的第一个维度是时间(每天1个文件),第二个和第三个空间维度是. 目前,由于以下原因,在第一维上切片(甚至是局部切片)将花费大量时间: netCDF文件以1x3712x3712的块大小进行了分块.在时间维度上进行切片基本上会读取整个文件. 在所有较小的文件中循环( ..
发布时间:2020-11-22 01:18:45 Python

使用R下载NetCDF文件:手动工作,download.file产生错误

我正在尝试从以下网站下载一组NetCDF文件:当我手动下载文件时,连接没有问题,但是当我使用download.file并尝试连接时,出现以下错误: 断言失败! 程序:C:\ Program Files \ Rstudio \ bin \ rsession.exe 文件:nc4file.c,第2771行 表达式:0 此应用程序已请求运行时以异常方式终止它. 请与应用程序的支持团队联 ..
发布时间:2020-11-11 01:00:56 其他开发

取得一段时间内每个季节的最大降雨值(xarray)

我试图找到10年内每个季节(DJF, MAM, JJA, SON)的最大降雨量.我正在使用netcdf数据和xarray尝试执行此操作.数据包括降雨(每3小时记录一次),纬度和经度数据.现在,我有以下代码: ds.groupby('time.season).max('time') 但是,当我这样做时,输出的形状为(4,145,192),表示在整个期间内每个季节都取最大值.我想要每年每个季 ..
发布时间:2020-07-28 05:23:20 Python