将数据集保存为 .netCDF 并将属性保存在 Python (xarray)
[英]save dataset as netCDF and also save attributes in Python (xarray)
我生成了新数据并首先保存为数据集 (.to_dataset)。 然后我为每个变量创建了属性 然后将数据合并为： 最后另存为.netCDF： 然而，当我打开并查看每个变量时，属性不可用。 有没有其他方法可以分配新属性并让它们生成 new.netCDF？ 提前致谢。 ...
使用 Python 将 hdf5 文件转换为 nc 文件
[英]Converting hdf5 file to nc file using Python
伙计们，我的目录中有一个 hdf5 文件列表，我想将 hdf5 文件转换为 nc 文件。 Python有没有办法做到这一点？ 我尝试使用 bash 但在 wsl 中遇到了很多问题。 所以我想在 Python 中执行此操作。我想有一种方法可以使用 .netCDF 库来执行此操作。 你能帮忙吗？ ...
大型 RasterStack 上的错误计算分位数
[英]Error computing quantiles on large RasterStack
我想在 R 上的一组堆叠的 RasterLayers 上跨层获取每个单元格的分位数。 我的 RasterStack 看起来像这样尺寸：721、1440、1038240、18628（nrow、ncol、ncell、nlayers） 我运行以下行以获得每个单元格的值为第 25 个和第 75 个百分位数的 ...
Python：用日期时间值替换 dataframe 中的列
[英]Python : Replace a column in a dataframe by datetime values
我正在尝试用我处理的日期时间值替换一列 4 列的数组。 问题是很难在 dataframe、数组、.... 的不同格式之间保持相同的形式。 如果我直接将转换后的日期放入数组中，结果如下：1.32542e+18 我试过 并将其包含在数组中，然后使用 我收到错误： 我试过也直接放： ...
如何用.netCDF土地利用数据和Key Biodiversity Areas (.shp)数据进行空间数据分析
[英]How to conduct spatial data analysis with netCDF land use data and Key Biodiversity Areas (.shp) data
我有一个 .netCDF 文件（“SSP119.nc”），其中包含 2015 年到 2100 年的土地利用变化信息，包括在内（来源： https://luh.umd.edu/data.shtml ）。 我还有一个包含全球关键生物多样性区域 (KBA) 多边形的 shapefile。 我想知道在某些年 ...
如何使用CDO统计满足设定条件的长度为5<=N<10天的事件数？
[英]How to count number of events of length 5<=N<10 days meeting a set condition using CDO?
我需要有关对 .netcdf 文件进行 CDO 操作的帮助。 我在网格区域从 ERA5 下载了 40 年的数据集，并使用 cdo 将范围 (30-50) 的变量值屏蔽为 1，将其他值屏蔽为 0。 现在我想计算每个网格单元在过去 40 年中连续 5 天但少于 10 天记录 var2= 1 的次数。 ...
将栅格文件（4 维）转换为允许进行随机森林分类的结构
[英]Turn raster files (4-dimensional) into structure that allows to conduct a random forest classification
我的目标是对农业用地 forms（作物分类）进行随机森林分类。 我对所有课程都有几个基本事实点。 此外，我有 37 个栅格文件 (.tif)，每个文件都具有相同的 12 个波段和相同的范围，一个文件代表时间序列中的一个日期。 时间序列不是恒定的。 下面显示了文件、日期和乐队名称以及用 terra 读 ...
使用xarray的open_mfdataset打开一系列.nc文件
[英]Using xarray's open_mfdataset to open a series of .nc files
我有一系列名为 0.nc、1.nc、2.nc 的文件，我希望使用 open_mfdataset 按文件名顺序一次打开它们。 但是，当我运行命令时： 我得到错误 我尝试通过以下方式添加坐标： 然后加载这些文件，但无济于事（同样的错误）。 有什么方法可以将这些文件作为一个数据集加载，而不会遇到这些错误？ ...
如何将 3 维数组平均化为 2 维数组 Python
[英]How to average a 3-D Array Into 2-D Array Python
我想从 python 中的 .netcdf 文件中获取一个温度变量，并对所有卫星扫描进行平均。 温度变量给出如下： tdisk 数组的形状是 68,52,46。 卫星每天进行 68 次扫描。 经度和纬度变量给出如下： 其大小为 52,46。 我想将每个 nscan 的温度平均在一起以获得每日平均 ...
如何提高光栅plot的分辨率？
[英]How to improve the resolution of a raster plot?
这是我的数据： https://github.com/dschmidt42/Benguela_Chla_Danielle_Schmidt 。 到目前为止我做了什么： 这是生成的 plot。分辨率非常小，所以非常块状。 我可以更改分辨率还是可以调到最高？ ...
将刻度颜色更改为 R 中的特定数字
[英]Changing scale colour to specific number in R
我已将我的数据上传到 github： https://github.com/dschmidt42/Benguela_Chla_Danielle_Schmidt 。 设置我的数据： 以下是结果图： 正如您所看到的，尽管它们具有不同的值，但比例都是相同的。 我怎样才能改变它，使一张图表上的 4.0 与 ...
如何交换 x 轴和 y 轴上的纬度和经度？
[英]How to swap latitude and longitude on x and y axes?
我已将我的数据上传到 github： https://github.com/dschmidt42/Benguela_Chla_Danielle_Schmidt 。 它涉及非洲东海岸的 chl a 数据，时间跨度为 26 年。 我试图将数据分成每个月，然后比较几个月内的差异。 这是我第一次使用 .n ...
如何打开多个 nc 文件，过滤所需的时间段并将结果写入单个 txt 文件？
[英]how to open several nc files, filter a desired period and write the result in a single txt file?
我有几个 nc 文件，我需要将它们全部打开。 过滤所需的时间段并将所有结果写入单个 txt 文件。 nc 文件对应一个月（一月、二月、三月...），并有四个变量（温度、露点、u 和 v）。 我需要组装一个表，在特定时期内并排显示所有变量。 例如，从一月到十月。 第一列是温度，第二列是露点，第三列 ...
无法从 .netcdf 文件中绘制箭袋
[英]Trouble plotting quivers from netcdf file
我正在尝试从 NetCDF 文件 plot 风箭袋，但是，我只能 plot 一行箭袋，直接穿过底部 map，如图所示。 代码如下。 非常感谢您的帮助:) 数据在这里，请替换为onedrive文件夹，谢谢https://drive.google.com/file/d/160121aFx0Ys6G1jd ...
前日12UTC至12UC之间的累计降水量
[英]Accumulated precipitation between 12UTC to 12UC of the other day
我有一个 .nc 文件，其中包含一整年每 6 小时降水量的数据，我的兴趣是计算每日降水量并与观测数据进行比较，为此我必须使它们在时间上重合。 为此，降水应在一天的 12 点到第二天的 12 点之间累积。 有人对如何使用 CDO 实现这一目标有建议吗？ 谢谢！ ...
适用于非直线网格的 CDO mergegrid？
[英]CDO mergegrid applicable on a non-rectilinear grid?
我在全球范围内有 16.netcdf 文件（没有时间维度），其中包含相同数量的经度网格点（即 432）和相同数量的纬度网格点（即 324）。 每个.netcdf 文件包含一个数据片段，例如在文件 1 中有一个纬度范围为 90 到 80 的数据片段，在文件 2 中有一个纬度范围为 80 到 70 的 ...
访问边界多边形内的 NetCDF 值
[英]Accessing NetCDF values within a bounding polygon
我正在尝试从 NetCDF 访问每日温度值以进行分析，但想在不同的行政单位内创建临时工摘要（即温度范围内的总天数）。 我还有一个全局 nc 文件和一个带有管理单元的 shapefile。 我的计划是通过循环访问纬度、经度和时间（三个临时参数）来读取临时数据，并将所需数据保存到列表中，但是我无法概念化 ...
无法在 xarray 中打开 Netcdf 变量
[英]Unable to open Netcdf variable in xarray
我正在使用 xarray 打开一个数据集，这条线是： DS = xr.open_dataset(dataDIR) 。 当我使用DS.variables时，我看到有一个名为“u-vel”的变量，但是，当我尝试使用DS.u-vel打开这个变量时，它会返回以下内容： “数据集”object 没有属性“u ...
如何使用单层 shapefile 屏蔽 multi-layer.netCDF 或栅格？
[英]How do I mask a multi-layer netCDF or raster with a single-layer shapefile?
我目前正在处理 .netCDF 格式的每日降水数据。 数据分辨率为 4 公里，覆盖美国。 但是，我想为特定地理区域（大约一个县的大小）使用更高分辨率的 shapefile 来屏蔽/剪辑数据。 最终，我希望 output 成为更小区域的每日降水数据，无论是高分辨率还是原始 4km 分辨率。 我尝试了 ...
无法反转 .netcdf 文件中的轴
[英]Inability to reverse axes in a netcdf file
我真的希望有人能帮助我。 我正在处理一个 NetCDF 文件，我设法沿深度维度提取加权平均值。 下一步是将文件排列成 ASCII 文件。 然而，当我将它上传到QGIS时，它出现顺时针和倒转（在图片下方找到）。 我发现，如果我将 NetCDF 作为网格文件上传到 QGIS 中，它会像图片中那样显示相反 ...