首页 文章
  • 1 votes
     answers
     views

    Java GRIB-Decoder:从GRIB2文件中提取数据

    我从这里下载了一些grib数据文件:ftp://data-portal.ecmwf.int/20160721000000/ (file type is .bin) 并想在我的Java应用程序中从这个文件中提取数据(我想稍后将提取的数据加载到数据库中) . 我只是尝试使用ftp://wmo:essential@data-portal.ecmwf.int/20160721000000/A_HWXE8...
  • 3 votes
     answers
     views

    使用R打开grib文件

    我正在使用R来处理气象数据 . 我分两步进行: 使用来自ncar命令语言的命令行函数ncl_convert2nc将grib转换为netcdf 使用R中的包ncdf导入netcdf数据 . I still have one problem: 2-对于某些特定的grib文件,使用ncar工具进行转换不起作用 . read grib files in R 还有其他方式或技巧(除了转录成n...
  • 2 votes
     answers
     views

    将grib转换为netcdf文件

    有没有办法在Windows上将grib文件转换为NetCDF格式?我使用一个名为 tkdegrib 的软件,但是它通过一个参数捕获一个参数,而我希望我的所有grib _1168624都使用Linux .
  • 0 votes
     answers
     views

    从空间(纬度/经度)数据定义网格框并提取R中的平均值

    我想通过定义经度/纬度网格框来计算我定义的数据区域的空间平均值 . 我的数据是ECMWF海冰数据,因此它是整个北半球每个.75x.75 lon / lat坐标的时空数据 . 我已经将数据从NetCDF格式更改为R数据帧,因此head(var.df)看起来像这样:列,日期,经度,纬度,值 date_time lon lat ci 1 2016-01-01 18:00:00 0 87.75 ...
  • 2 votes
     answers
     views

    f2py与fortran netcdf

    我正在尝试使用f2py来编译一些使用netcdf库的fortran代码 . 编译似乎工作,但当我尝试在python中导入模块时,我收到错误消息: Traceback (most recent call last): File "<stdin>", line 1, in <module> ImportError: ./cdfio_simple.so: un...
  • 2 votes
     answers
     views

    在NetCDF中沿着无限维编写标量变量

    我正在尝试将流体动力学模型的时间变量写入netcdf文件(无限维变量) . 我在Fortran90中附加了一个简化的代码示例,突出了我的问题 . 根据用户指定的输出间隔(本例中为10次),在模拟期间多次调用写入netcdf文件的子例程 . 我可以在第一次调用子例程时创建文件并添加属性 . 在后续的子程序调用期间,我无法将启动和计数变量设置为正确,以便将时间变量写入文件 . 这是错误,在编写模型时间...
  • 0 votes
     answers
     views

    netcdf维度变量解释

    我试图了解NetCDF标准是否允许这样做 . 它并没有对我产生影响,但也许有一个原因是它在图书馆一级不被禁止 . Ncdump: netcdf tt { dimensions: one = 2 ; two = 1 ; variables: int64 one(two) ; data: one = 1 ; } 以及在python中生成此文件的代码: from net...
  • 1 votes
     answers
     views

    高分辨率再分析数据

    当我从netCDF文件重新分析(可变压力(SLP),01/01/2014)中提取数据时,数据的分辨率非常高(9km网格),这使得生成的图像非常嘈杂 . 我想将数据放入较低分辨率的网格(例如1度) . 我正在尝试使用meshgrid和gridata函数,但缺乏经验无法使其工作 . 有谁知道如何解决?谢谢 . from netCDF4 import Dataset import numpy a...
  • 0 votes
     answers
     views

    如何在netCDF中为现有变量添加新维度?

    代码是这个答案确实有效Add a new dimension但它没有做我想做的事情 . 为了进一步说明我的问题 - 我有以下定义时间(无限制),纬度,经度,水平为四维和一个变量位势高度 . 我想要做的是将维度dimz_Z添加到现有的变量位势高度 . 现在,上面的代码中的答案是netCDF文件与ncdump的相似之处 float hgt(time, level, lat, lon) 但是dimz_...
  • 0 votes
     answers
     views

    在netcdf文件中重命名维度

    我是 nco 的初学者,我希望对我的工作流程提出一些建议,并对我遇到的问题提供一些帮助 . 我有these data,它包含两个不同时间步的3D盐度值: dimensions: t = 780 ; z = 54 ; y = 450 ; x = 3600 ; variables: double time(t) ; time:units = &q...
  • 1 votes
     answers
     views

    在python底图中绘制gebco数据

    我已经下载了一些gebco测深数据作为netCDF文件 . 我想用python-basemap绘制它 . 我试过了, import netCDF4 from mpl_toolkits.basemap import Basemap # Load data dataset = netCDF4.Dataset('/home/david/Desktop/GEBCO/gebco_08_-30_45_5_...
  • 2 votes
     answers
     views

    PYTHON:在底图上覆盖netCDF数据(contourf)

    我是python和netCDF的新手,但是现在已经在以下问题上工作了几天而没有任何进展 . 我的代码基于我读过的各种教程和示例,我无法弄清楚出了什么问题! 基本上我希望能够将一些netCDF数据的图叠加到底图上 . 我的代码(完整)如下: import numpy as np from netCDF4 import Dataset # from scipy.io import netcdf imp...
  • 0 votes
     answers
     views

    Python ValueError:没有找到HDF5标头

    我正在尝试使用其setup.py在python中安装rvic . 当我在命令提示符下运行python setup.py install时,它显示 Value error: did not find HDF5 headers . 我也尝试安装 h5py ,但仍然会出现同样的错误 . 有什么建议? Stackoverflow不允许我发布这个,所以我添加了一些关于这个的更多细节... rvic是一个水...
  • 1 votes
     answers
     views

    NetCDF:开始计数超出维度限制

    我在Fortran中编写了一个代码,用于读取具有4维数据[时间,级别,经度,纬度]的NetCDF文件 . 但是,我的代码会产生错误 NetCDF: Start+count exceeds dimension bound 在我正在使用的任何4-bit NetCDF文件上 . 例如,http://people.sc.fsu.edu/~jburkardt/f_src/netcdf/pres_temp_...
  • 15 votes
     answers
     views

    如何在ggplot2中正确绘制投影网格化数据?

    多年来我一直在使用 ggplot2 绘制气候网格数据 . 这些通常是预计的NetCDF文件 . 单元格在模型坐标中是方形的,但取决于模型使用的投影,在现实世界中可能不是这样 . 我通常的方法是首先在合适的常规网格上重新映射数据,然后绘图 . 这引入了对数据的小修改,通常这是可以接受的 . 但是,我已经决定这已经不够好了:我想直接绘制投影数据,而不重新映射,因为其他程序(例如 ncl )可以,如果我...
  • 0 votes
     answers
     views

    xarray DataArray.where()在屏蔽时减少了坐标

    xarray新手在这里 . 非常简单的情况,我有一个降水型阵列(ntim x nlat x nlon)和一个总降水阵列(相同的尺寸) . 两者都在单独的netCDF文件中 . 我想掩盖降水阵列,其中A)降水量下降(> 1e-8 m / s率)和B)降水类型为雪(maskvar = 0.0) . 因此输出数组是“下雪的地方?”阵列 . 当使用xarray where()具有来自两个不同(但相同...

热门问题