NetCDF重新排序维度
NetCDF（Network Common Data Form）是一种用于存储科学数据的文件格式和库。它可以用于存储、管理和共享多维数据集，具有平台无关性和自描述性的特点。
NetCDF文件由变量（Variable）、维度（Dimension）和属性（Attribute）组成。维度描述了数据的大小和形状，变量包含了实际的数据，而属性则是关于数据和变量的附加信息。
重新排序维度是指修改NetCDF文件中变量的维度顺序。这可以通过改变维度顺序来改变数据的存储方式和访问方式。在一些特定的数据处理和分析场景中，重新排序维度可以提高数据访问的效率和方便性。
NetCDF提供了一种灵活的数据访问方式，可以通过指定维度的顺序来对数据进行重新排序。这种操作在处理多维数据集时非常有用，特别是当数据集非常大时。
NetCDF没有直接提供重新排序维度的方法，但可以通过以下步骤实现：
- 使用NetCDF库加载并打开待处理的NetCDF文件。
- 读取需要修改的变量以及相关的维度信息。
- 根据需求重新排序维度，可以通过交换维度的顺序或者重命名维度来实现。
- 更新变量的维度信息。
- 将修改后的数据保存为新的NetCDF文件。
在云计算领域，NetCDF可以与其他云服务和平台集成，实现科学数据的存储和处理。腾讯云提供了一系列与科学计算和数据处理相关的产品和服务，包括云服务器、对象存储、云函数等。这些产品可以与NetCDF结合使用，满足科学研究、气象预测、环境监测等领域对数据存储和处理的需求。
作为腾讯云的相关产品，对象存储 COS（Cloud Object Storage）可以用于存储NetCDF文件，提供安全可靠、高扩展性的对象存储服务。您可以通过腾讯云COS产品页面（https://cloud.tencent.com/product/cos）了解更多关于该产品的详细信息。
需要注意的是，以上回答仅针对NetCDF重新排序维度的概念和相关推荐，不涉及其他云计算品牌商。如有其他问题或者需要更多信息，请告诉我。
相关·内容
1回答
我正在使用python中的netCDF4创建一个NC文件。该文件非常简单，包括1和2维变量(1 dim=N，2 dim=M)。但问题是，当我读出文件时，二维变量的维数被互换，即不是以维数(N，M)的形式出现，而是以维数(M，N)的形式出现，而所有的一维数据都是(N，)，正如预期的那样。
代码看起来像这样
root = nc.Dataset(filename, 'w', format='NETCDF4')
dimensions = ('N', 'M')
root.createDimension(dimensions[0], Non
浏览 2提问于2015-02-02得票数 0
回答已采纳
1回答
我正在为一个项目使用gridMET (http://www.climatologylab.org/gridmet.html)和MACA (http://thredds.northwestknowledge.net:8080/thredds/reacch_climate_CMIP5_macav2_catalog2.html) NetCDF文件，并且面临性能问题。在gridMET NetCDF文件上实现一个简单的函数(持续时间: 1979-2015)大约是0.01秒/网格单元。然而，具有与gridMET数据相同功能的MACA NetCDF文件(持续时间: 2016-2050)的处理速度约为0.3
浏览 50提问于2019-06-16得票数 1
回答已采纳
1回答
是否有NCO或netCDF4命令可用于扩展netCDF文件的lon维度。
我有一个与foll有关的netCDF文件。酒窝：
尺寸(尺寸)：时间(10)、瓦尔(1)、纬度(1674)、经度(4320)
我想把纬度维度扩展到整个地球，也就是说，不是1674年，应该是2160。有办法吗？应该为新单元分配一个用户指定的值，例如0.0。
浏览 2提问于2017-01-17得票数 1
回答已采纳
2回答
我试图将多个netCDF文件与相同的维度组合在一起，它们的维度如下：
OrderedDict([(u'lat', <type 'netCDF4._netCDF4.Dimension'>: name = 'lat', size = 720
), (u'lon', <type 'netCDF4._netCDF4.Dimension'>: name = 'lon', size = 1440
), (u'time', <type 'netCDF4._ne
浏览 3提问于2016-09-18得票数 3
回答已采纳
问题摘要
我正试图将a.grib2文件转换成一个netcdf，该文件代表了跨越美国大陆的一天的网格化雷达降雨数据。当.grib2缺少时间步骤时，我试图在运行xarray.Dataset.to_netcdf之前使用xarray.Dataset.reindex填充它们的NA值。但是，在重新编制数据集索引之后，由于内存分配错误，脚本将失败。如果我不重新索引它就成功了。一个线索可能是数据集块被设置为(70, 3500, 7000)，但是当调用ds.to_netcdf时，脚本会失败，因为它试图加载一个包含维度(210, 3500, 7000)的块。
访问完全可复制的示例
用于复制我的结果的代码和数据可以
浏览 18提问于2022-09-10得票数 0
1回答
如何对维度进行重新排序？
、、、
我正在处理降水数据。我需要在NetCDF文件中重新排列尺寸。我的文件中的尺寸是(time, lat, lon)顺序的。我需要它们按(lat, lon, time)顺序排列。
我尝试了以下几种方法：
ncpdq -a lat,lon,time infile.nc outfile.nc
现在，我得到了我想要的维度的顺序，也就是(lat, lon, time)。
但是，lat维度变成了UNLIMITED，这是错误的。时间维度应为UNLIMITED维度。这是我得到的：
dimensions:
lon = 720 ;
lat = UNLIMITED ; // (360 current
浏览 3提问于2019-04-05得票数 3
我能够弄清楚后面的维度声明，如下面的代码所示。我想我马上就能得到netCDF文件了。但是，我得到了一个错误。
import numpy as np
import netCDF4
import os
# load the data
path='C:\Users\.spyder2'
os.chdir(path)
# this load the file into a Nx3 array (three columns)
data = np.loadtxt('TRMM_1998_01_0100_newntcl.csv', delimiter=',')
浏览 2提问于2015-03-17得票数 1
我正在尝试组合netcdf文件，但它在netCDF4.MFDataset.init (netCDF4.c:3822)中的f= MFDataset(fili)文件"utils.pyx“和第274行中显示了"CBL_plot.py"，第11行，IOError: MFDataset THref_11:00.nc没有聚合维度。
因此，我只检查了一个netcdf文件，并且netcdf文件的信息如下所示：
float64 th_ref(u‘，)无限维= ()当前尺寸= (30，)
它似乎没有聚集维度。但是，我希望将这些netcdf文件组合在一起，而不是一个接一个地使用。是否有任何方
浏览 6提问于2013-12-03得票数 0
回答已采纳
1回答
我在处理NetCDF4文件时遇到了问题。有人能帮帮我吗。把它读进去似乎很有效。
import netCDF4
f = netCDF4.Dataset('mydata.nc', 'r')
当我试图调查文件中的内容(>>>打印f)时，我得到了一些信息：
<type 'netCDF4.Dataset'>
root group (NETCDF4 file format):
dimensions: soundings, levels
variables:
groups: Retrieval, Sounding
..。
pri
浏览 4提问于2013-10-22得票数 1
回答已采纳
我正在编写一个脚本，该脚本使用.nc打开多个未排序(时间维度中)的xarray.open_mfdataset文件，并根据时间维度对它们进行排序，然后将它们全部保存为一个.nc文件。合并数据集的大小为时间: 452652，级别: 2000。
import xarray as xr
ds = xr.open_mfdataset('2008_*.nc',combine='nested',concat_dim='time',
chunks={'time':10000,'level'
浏览 14提问于2022-05-05得票数 0
我正在尝试创建一个基于VTK文件的netCDF文件，其中包含.vti extension.The元数据或VTK文件的头信息，如下所示：
<?xml version="1.0"?>
<VTKFile type="ImageData" version="0.1" byte_order="LittleEndian">
<ImageData WholeExtent="0 256 0 192 0 128" Origin="0 0 0" Spacing=" 4.000
浏览 2提问于2015-05-18得票数 0
回答已采纳
2回答
我刚开始使用Python，也不熟悉NetCDF，所以如果我不太清楚，很抱歉。我有一个包含多个变量的nc文件，我需要以新的顺序从这些nc文件中提取数据。
我的nc文件有8个变量(经度，纬度，时间，u10，v10，swh，mwd，mwp)，我尝试的逻辑是“如果我输入经度和纬度，我的程序输出按时间排序的其他变量(u10，v10，swh，mwd，mwp)。”然后，我会将提取的数据放入另一个数据库中。
我测试我的nc文件，如下所示：
import netCDF4
from netCDF4 import Dataset
jan = Dataset('2016_01.nc')
print
浏览 24提问于2017-08-09得票数 4
回答已采纳
我正在创建一个NetCDF4文件，该文件目前有四个变量：
1)地表温度(三维阵列-时间、纬度、经度)
2)经度(每个像素中心的一维坐标)
( 3)纬度(每个像素中心的一维坐标)
4)时间(自1900-01-01 00:00:00:00起的图像采集时间)
我目前正在使用以下代码来执行此操作：
#==========================WRITE THE NETCDF FILE==========================#
newfile = nc.Dataset(export_filename, 'w', format = 'NETC
浏览 8提问于2017-01-24得票数 1
回答已采纳
假设我有一个netcdf文件，其中一个名为“time”，大小为500。
nc=netcdf(ncfile, 'write')
size(nc{'time'})
ans =
500
现在，我希望在同一个变量中存储一个更长的时间序列(比如750个成员)，但是我希望文件的其余部分保持不变。如何将此变量的维度更改为750，同时保持netcdf文件中其余变量的原样？
我尝试将其简单地附加到数组中，这不会产生错误，但是数组大小不会改变(这意味着它根本不写值)。
nc{'time'}(750)=1
我试图编写一个我试着将文件放入重新定义模式并更改大小的文件，但
浏览 6提问于2014-12-23得票数 2
回答已采纳
我有一个具有一个变量(前面)和四个维度(时间、高度、纬度和经度)的NetCDF。从下载
它是月度组合，即<xarray.DataArray 'time‘(Time151)> array('2001-01-16T12:00:00.000000000'，'2001-02-16T00:00:00.000000000'，...'2013-12-16T12:00:00.000000000'，dtype=’datetime64 is‘)。
我想为每个时间戳创建一个单独的文件(无论是NetCDF还是Geotif都没有关系)。
我试过了：
d
浏览 4提问于2020-11-24得票数 0
我正在处理使用NetCDF格式的ROMS模型输出生成的洋流数据。每个NetCDF文件都包含洋流速度的月平均值，即只有一个时间步长。到目前为止我已经到达这里了。
import netCDF4
import pandas as pd
import numpy as np
import xarray as xr
df1 = xr.open_dataset("ocean_avg_November.nc4") #NetCDF file for Nov 2013
df2 = xr.open_dataset("ocean_avg_December.nc4") #NetCDF
浏览 6提问于2022-06-18得票数 0
回答已采纳
我的目标是导入一个NetCDF文件并提取其中存储的底层变量。我一直在跟踪这个。
到目前为止，这是我的代码：
from netCDF4 import Dataset
import numpy as np
fileName = 'LES_cloud_fields.nc.cdf'
fdata = Dataset(fileName, 'r', format="NETCDF4")
print(fdata.dimensions)
print(fdata.variables['nz'][:])
fdata.close()
我得到错误的原因：K
浏览 3提问于2017-06-14得票数 2
回答已采纳
我有一个4维的netCDF文件。我想通过提供一个维度的名称从netCDF文件中提取一个切片。
我知道如何按位置来做这件事。例如。
from netCDF4 import Dataset
hndl_nc = Dataset(path_to_nc)
# Access by slice
hndl_nc.variables['name_variable'][:,5,:,:]
假设我知道维度的名称，比如A，B，C，D。如何按维度名称而不是位置访问？
浏览 5提问于2016-09-17得票数 0
回答已采纳
2回答
我试图理解NetCDF标准是否允许这样做。这对我没有意义，但也许有一个原因，为什么它不被禁止在图书馆级别。Ncdump：
netcdf tt {
dimensions:
one = 2 ;
two = 1 ;
variables:
int64 one(two) ;
data:
one = 1 ;
}
以及在python中生成此文件的代码：
from netCDF4 import Dataset
rr=Dataset('tt.nc','w')
rr.createDimension('one',2)
rr.createD
浏览 5提问于2015-07-07得票数 0
1回答
我有几个数据集(TEMP，DEWPOINT)预测数据的.nc文件大小约300MB。我需要将(TEMP)数据集转换为多个GEOTIFF (每个时间片一个.tif )。
下面是.nc文件的外观。
查看了此，但它似乎适用于整个数据集。
我尝试了GDAL，但不确定如何让它在每个时间片上工作。
有什么想法吗？netcdf4-python？
浏览 7提问于2018-08-28得票数 7
回答已采纳
1回答
我的代码：
data = Dataset('/scratch/uni/ifmto/janniklas/data.nc', 'w', format='NETCDF4_CLASSIC')
data.createDimension('z', len(zn))
data.createDimension('x', len(xn))
data.createDimension('y', len(yn))
nz = data.createVariable('zn', np.float32, &#
浏览 4提问于2017-04-26得票数 1
回答已采纳
我有一个数据数组(1500 X 2500)，我希望将其写入netcdf文件，但经度/经度网格的间距不是规则的。我可以创建一个包含数据的netcdf文件，但该文件似乎不可读和/或格式不正确。例如，当我试图用NOAA天气和气候工具包读取它时，我得到一个错误，指出文件被扫描为‘网格’文件，但没有找到网格。下面是我的代码： from netCDF4 import Dataset
ndvi_nc = Dataset("ndvi_out.nc", "w", format="NETCDF4")
print(ndvi_nc.data_model)
tim
浏览 32提问于2020-04-03得票数 2
回答已采纳
1回答
我有一个网格温度数据集df (时间: 2920 x: 349 y: 277)和同一网格mf的陆地海掩膜(时间:1 x: 349 y: 277)，其中陆地网格点mf.land =1，海洋点mf.land =0。我想使用陆地海掩码从我的温度数据集df中消除海洋点，也就是说，我只想要df中的网格点，其中mf.land = 1。
下面是df的样子：
下面是mf的样子：
我在尝试这个：
#import libraries
import os
import matplotlib.pyplot as plt
from netCDF4 import Dataset as netcdf_datase
浏览 9提问于2022-07-17得票数 1
回答已采纳
1回答
我有大量的NetCDF文件，我想从这些文件中提取一个位置的少量变量，并将它们合并到一个新的NetCDF文件中。这些文件的尺寸如下：
dimensions:
time = 18 ;
level = 65 ;
levelh = 66 ;
domain = 36 ;
我可以减去/合并所有domains的文件，如下所示：
cdo select,name=u,v file1.nc file2.nc out.nc
但是所有其他操作符似乎都与空间中的选择(例如sellonlatbox)或时间(例如seltimestep)有关，但我无法从NetCDF文件中选择一个domain。
浏览 0提问于2018-09-26得票数 2
回答已采纳
2回答
我下载了一个包含4-5个变量的netcdf文件，但它只有2个维度(经度和经度)。缺少时间，这不允许我合并时间步长或做任何有用的事情。
有没有办法通过使用CDO来解决这个问题呢？
有100个netcdf文件(没有时间维度)，我想使用时间作为合并的主要变量来合并它们。
浏览 1提问于2014-05-30得票数 2
<xarray.DataArray 'nc_file' (time: 20, latitude: 360, longitude: 720, N: 3)>
[34992000 values with dtype=float64]
Coordinates:
* latitude (latitude) float64 89.75 89.25 88.75 88.25 87.75 87.25 86.75 ...
* longitude (longitude) float64 -179.8 -179.2 -178.8 -178.2 -177.8 -177.2 ..
浏览 13提问于2017-03-12得票数 2
1回答
我想使用netCDF4-Python从netCDF文件中选择变量，我使用以下代码：
import numpy as np
import pandas as pd
from pylab import *
from netCDF4 import Dataset
pals = Dataset('pals_amplero_2003_2006_ecmwf_v1.nc4', "a",format='NETCDF4')
print (pals.variables.keys())
print (pals.variables['Rainf']
浏览 0提问于2017-11-16得票数 0
回答已采纳
1回答
我有netCDF数据集，x，y作为空间维(在Lambert投影中)，它只是枚举值0:495，0:309，和lat，lon变量作为二维网格的形状(309,495)。
我想把lon，lat变量赋给x，y维数来获得经度和纬度的投影数据，但是我找不到任何参考。
转换现有数据集将是很好的，但可接受的选项也是重新创建新的netCDF数据集从现有，或生成hdf-5代替。
很抱歉，我的问题很难理解。我将尝试另一个例子，它并不意味着地理空间投影。
我有三维netCDF数据集。垂直维度altitude以单位hPa表示，而dataset包含变量height，该变量将每个altitude级别映射为米。我希望更改ne
浏览 7提问于2014-08-30得票数 1
回答已采纳
1回答
Python绘图数据(概要文件)
、、、、
我的物理项目有一个问题，我必须使用Python，即使我不知道如何编程。我已经阅读了我的数据，我看了很多不同的教程，我不知道下一步该做什么。我为此挣扎了很长一段时间。
from netCDF4 import Dataset as ncf
file='atmPrf_C001.2010.227.00.03.G04_2013.3520_nc'
data=ncf(file,'r’)
print(infile.dimensions)
print(data.variables)
vertical=data.variables['MSL_alt']
latitude=
浏览 4提问于2017-09-14得票数 1
回答已采纳
我正在尝试处理来自非常大的netCDF文件(每个大约400 Gb )的数据。每个文件都有几个变量，都比系统内存大得多(例如，180 Gb与32 Gb RAM)。我正在尝试使用numpy和netCDF4-python对这些变量执行一些操作，方法是一次复制一个切片，然后在该切片上操作。不幸的是，仅读取每个切片就需要很长时间，这会扼杀性能。
例如，其中一个变量是shape (500, 500, 450, 300)的数组。我想对切片[:,:,0]进行操作，所以我执行以下操作：
import netCDF4 as nc
f = nc.Dataset('myfile.ncdf','
浏览 4提问于2012-08-22得票数 13
回答已采纳
2回答
尝试从NOAA服务器(即thredds aws、亚马逊等)访问原生NetCDF格式的原始NEXRAD数据并将其转换为网格NetCDF文件，以便自动进行后处理。
美国国家海洋和大气局的能够导出网格的NetCDF3文件(在下面的dropbox文件夹中)。
我已经训练了一个外部代码来处理导出的网格文件。但是，为了自动化这个过程并为我的研究分析几年，我必须弄清楚如何将NOAA服务器中的原生文件转换为这种网格格式。具体地说，我希望在最终的输出文件中有一个以dBz表示的'reflectivity‘的MxN数组，扫描角度最低，约0.5度，纬度和经度的1D数组大小分别为M和N。
为了帮助理解我要查找的
浏览 16提问于2018-09-05得票数 0
1回答
数组聚合
、、
我想在R中处理的netCDF4文件太大了。我想要写一个循环，它将读取一部分数据并对其进行总结。
我想读的变量有4个维度：“lat”、“lon”、“member”和“time”。时间有一个月分辨率和成员包含60集合运行从一个气候模式。
使用ncvar_get命令，我从netCDF中提取了12个时间片，给我留下了一个4维数组。
num[1:144,1:69:1:60,1:12]
我将如何汇总这些数据，这样我就有了年度数据。
浏览 1提问于2016-06-27得票数 0
回答已采纳
1回答
我今天的工作是在Numpy中获取一个数组，并使用file将其转储到一个(.nc)中。这个有一个特定的子模块，但是创建一个实际的netcdf文件的例子不多。
以下是进入文件所需的数据：
纬度(94)
经度(192)
与上述矩阵对应的温度数据
我已经保存了在Numpy中保存的所有3条数据，我现在必须使用NETCDF4 python包(我也有.nc python包，但我无法使用它)将这些数据放入.nc文件中。
到目前为止，这是我的代码：
#testing code about difference
f = netcdf.netcdf_file('2m air temper
浏览 3提问于2015-10-24得票数 0
回答已采纳
1回答
创建带范围的维度
、、
我正在尝试从一个大型数据集创建一个netCDF文件。我要这个新文件保存温度和时间。
我的时间是一个范围，这是造成问题时，我试图创造我的维度。我一直在犯这个错误：
无法在此组或父组中找到维度计时器。
filehandle = Dataset(path+indices,'r',format="NETCDF4")
lat = 0
lon = 0
timerange = range(0, 1459)
temps = get_var(filehandle, 'tair')[timerange, lat, lon]
ncfile = Datase
浏览 1提问于2015-11-10得票数 0
回答已采纳
我正在使用python中的netCDF4库，只是遇到了标题中提到的问题。起初，我将此归咎于组，但事实证明，这是NETCDF4和NETCDF3_CLASSIC格式之间的区别(编辑:它似乎与我们安装的netcdf库有关)。
在下面的程序中，我以两种不同的方式为相同数据创建了一个简单的时间序列netcdf文件: 1)作为NETCDF3_CLASSIC文件，2)作为NETCDF4平面文件(在netcdf4文件中创建组并没有多大区别)。通过简单的计时和ls命令，我发现：
1) NETCDF3 1.3483 seconds 1922704 bytes
2) NETCDF4 fl
浏览 4提问于2014-11-27得票数 3
回答已采纳
我已经写了一个ncks代码来从Hycom数据库中提取数据，但是，由于服务器问题，我一次只能下载一个timestep。我需要一个循环，它允许我下载一个月的时间段的单一时间步数据。
ncks -d time,2015-12-01T00:00:00Z -d depth,0.00,1000.00,1 -v water_temp,salinity,water_u,water_v,surf_el -d lon,161.45,171.28 -d lat,-4.75,3.43 'http://tds.hycom.org/thredds/dodsC/GLBv0.08/expt_53.X/data/201
浏览 6提问于2021-05-26得票数 0
我正在尝试将数据写入netcdf文件，使用此示例代码时，我得到以下错误： Traceback (most recent call last):
File "netCDF4\_netCDF4.pyx", line 4866, in netCDF4._netCDF4.Variable.__setitem__
ValueError: cannot reshape array of size 5 into shape (5,1,1,3) 下面是一个你可以自己尝试的示例测试代码； from netCDF4 import Dataset
latitude = [21.1877]
l
浏览 0提问于2020-07-14得票数 0
2回答
我有许多坐标(大约20000个)，我需要从许多NetCDF文件中提取数据，每个文件大约有30000个时间步长(未来的气候情景)。使用解决方案here效率不高，原因是在每个i，j处花费时间将"dsloc“转换为"dataframe”(请查看下面的代码)。**可以从here下载示例NetCDF文件** import pandas as pd
import xarray as xr
import time
#Generate some coordinates
coords_data = [{'lat': 68.04, 'lon': 15.20, &
浏览 88提问于2021-09-25得票数 6
回答已采纳
我正在以netcdf格式下载气候数据。对于每个变量(例如“降水量”)，我需要合并9个netcdfs，每个netcdfs属于一个独特的气候模式。每个netcdf都有相同的大小(时间、lat、lon)。如何将9个3D netcdfs合并成一个4Dnetcdf？最后，我想计算一个月的累积降水量。以下是我的当前代码：
variables = ['pr']
scenarios = ['historical', 'ssp245'] #options ['historical', 'ssp126'
浏览 18提问于2022-09-06得票数 0
在计算循环的每次迭代中，我都会将计算结果写入NetCDF变量，因为在对变量执行一次写入操作之前，我无法将所有结果数据都保存在内存中。我的假设是，每次写入NetCDF变量都只会写入我正在写入的数据片段，因此每次写入都会很快，但是每次写入似乎都需要相当长的时间，我担心的是，在每个变量赋值语句中，整个NetCDF文件都会被重写，而不仅仅是被赋值的数据片段，这会导致相当不必要的I/O开销。
结果数据的切片被写入输出NetCDF变量，如下所示：
outputSpi01MonthVariable[totalNumberOfOutputMonths, lat:lat + 1:1, len(lonDimens
浏览 1提问于2013-06-12得票数 2
1回答
使用下面的代码，我从本地保存的netCDF文件中读取数据。X和y变量是经度和纬度元素的一维列表。相反，sal_values是一个掩码数组；一个列表列表。
这绝对没问题；但我想将索引列表排序到一个2D数组中，以便根据x和y绘制我的数据。例如，当我运行sal_values 3时，我收到错误消息“sal_values 1 is out of bound for axis 0 with size 1”，这意味着sal_values不是一个实际的2D数组。
import cartopy.crs as ccrs[enter image description here][1]
import matplot
浏览 20提问于2020-02-11得票数 0
回答已采纳
1回答
我试图从netCDF文件中提取一个变量并打印出来。这是我的密码
import netCDF4
import netCDF4_utils
from netCDF4 import Dataset
from numpy.random import uniform
import csv
B = []
rootgrp = Dataset('test.cdf', 'r', format = 'NETCDF4')
f = open('testoutput.csv','wb')
b = (rootgrp.variables[&
浏览 5提问于2013-07-23得票数 1
回答已采纳
1回答
我正在使用netCDF函数读取多个xarray.open_mfdataset()格式的数据文件(WRF模型输出文件)。我正在阅读和引导各种变量，如变量QVAPOR、U和V resp。我使用以下代码读取netCDF变量并计算结果变量UQ。
import xarray as xr
def desta_var(pp,var):
flnm = xr.open_mfdataset(pp)
if var=="U":
U1 = (flnm.variables[var])
U = 0.5*(U1[:,:,:,0:-1] + U1[:,:,:,1
浏览 0提问于2018-03-09得票数 4
2回答
我已经记录了每日温度数据，但我只对冬季的月份感兴趣。
from netCDF4 import Dataset as netcdf_dataset
import numpy as np
import xarray as xr
#open NASA GISS gridded temperature netcdf file
df = xr.open_dataset('BerkeleyEarth.nc')
#pull out temperature variable
air=df.temperature
#select only winter months
WinterAir =
浏览 20提问于2022-02-28得票数 1
1回答
我想将“切片”映射从dataset保存到一个新的netcdf文件中。有关守则如下：
library(ncdf)
ncfname <- “cru10min30_tmp.nc”
crutmp <- open.ncdf(ncfname)
tmp.array <- get.var.ncdf(crutmp,"tmp")
close.ncdf(crutmp)
month <- 1
tmp.slice <- tmp.array[,,month]
数据集在上是可用的，这个代码片段是从上的一个示例中获取的。如何将"tmp.slice“二维数组保存为新的ne
浏览 5提问于2013-09-17得票数 1
回答已采纳
我有几个Python语言的列表，我正在尝试将它们写入NetCDF文件，但我不确定如何让它工作。 下面是列表结构的一个示例： 每个列表由多个列表本身组成，例如17个，每个列表有3个元素。第一个元素是一年，所以只有一个整数或字符串，第二个元素由36个数组组成，第三个元素由5个数组组成。 基于其他指南，我尝试以这种方式编写代码： fn = 'mean_yearly_list.nc'
ds = nc.Dataset(fn, 'w', format='NETCDF4')
year = ds.createDimension('date',
浏览 14提问于2021-06-21得票数 1
1回答
NcML可以用来聚合3D和4D网格吗？我不确定，因为它们在维度的数量上不同。例如，海面高度( ssh )和水温，其中ssh具有三维时间，纬度，经度，温度具有四个维度时间，深度，纬度，经度？我的测试没有成功，所以我的直觉是我必须将3D和4D变量分解成单独的目录。但我希望其他人能有其他的建议？
虽然我尝试使用以下代码片段的“联合”聚合，但时间维度没有得到适当的映射，因为3D变量从2008-12-28开始，4D变量从2008-05-08开始：
<?xml version="1.0" encoding="UTF-8"?>
<netcdf xmlns=&
浏览 20提问于2013-04-24得票数 1
回答已采纳
我正在使用xarray.DataArray向应用一个函数。它适用于某些NetCDF，而与其他在维度、坐标等方面似乎相当的NetCDF失败。但是，在代码工作的NetCDF与代码失败的NetCDF之间肯定有一些不同，希望有人在看到代码和关于下面列出的文件的一些元数据之后，能够评论问题所在。
我正在运行的执行计算的代码如下：
# open the precipitation NetCDF as an xarray DataSet object
dataset = xr.open_dataset(kwrgs['netcdf_precip'])
# get the precipitat
浏览 0提问于2018-11-01得票数 4
回答已采纳
我使用以下代码向fortran中的现有netCDF文件添加了一个新的维度-
retval = nf_open(cfn,NF_WRITE,ncid)
if (retval .ne. nf_noerr) call handle_err(retval)
retval = nf_redef(ncid)
if (retval .ne. nf_noerr) call handle_err(retval)
retval = nf_def_dim(ncid,"xyz",len,dimid_xyz)
if (retval .ne. nf_noerr) call handle_err(retval
浏览 3提问于2016-08-23得票数 0
回答已采纳
我正在尝试从IRI/LDEO气候数据库(dust_pm25_sconc10_mon)读取NetCDF文件，但我在读取此文件时遇到问题。当我选择组成数据库的变量(经度(X)、纬度(Y)和时间(T))时，X和Y的输出是一个具有观察值数量的序列(例如1，2，...，139 )。也就是说，经度和纬度的值不会正确导出。
有人能帮我解决这个问题吗？我已经尝试过用R、Python和Qgis读取这个文件，在这三个文件中，X和Y的输出是相同的。
我的代码如下(Python)。
非常感谢大家。
from netCDF4 import Dataset as dt
filestr = 'dust_pm25_
浏览 1提问于2019-04-21得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
实时音视频活动推荐
腾讯云开发者
