Python netCDF4中group的使用

netCDF4中group的使用

group可以当作一个dataset一样处理,两个类拥有大多数相同的函数,在一个group中也能继续创建子group。group中的维度可以与全局维度一起使用,如果维度名称相同,则在创建变量的时候会使用group中的维度。

创建group

import netCDF4 as nc
# 创建数据
dataset = nc.Dataset(r'H:\testgroup.nc',mode='w')
# 创建group
group = dataset.createGroup('1')
dataset.close()

创建group

在group中创建变量

import netCDF4 as nc
import numpy as np
dataset = nc.Dataset(r'H:\testgroup.nc',mode='w')
group = dataset.createGroup('1')
# 创建维度和变量
group.createDimension('lat',10)
group.createDimension('lon',10)
group.createVariable('lat',np.float32,('lat'))
group.createVariable('lon',np.float32,('lon'))
group['lat'][:] = np.arange(10)
group['lon'][:] = np.arange(10)
dataset.close()

在group中写入数据

group维度与全局维度

import netCDF4 as nc
import numpy as np
dataset = nc.Dataset(r'H:\testgroup.nc',mode='w')
# 创建全局维度
dataset.createDimension('lat',11)
dataset.createDimension('lon',11)
dataset.createVariable('lat',np.float32,('lat'))
dataset.createVariable('lon',np.float32,('lon'))
dataset.variables['lat'][:] = np.arange(11)
dataset.variables['lon'][:] = np.arange(11)
# 创建group
group = dataset.createGroup('1')
group.createDimension('lon',10)
group.createVariable('lat',np.float32,('lat'))
group.createVariable('lon',np.float32,('lon'))
group['lat'][:] = np.arange(11)
# 将10改成11会报错,说明这里的lon使用的是group中的维度
group['lon'][:] = np.arange(10)
dataset.close()

全局维度与局部维度

创建与读取多组group

根据我的需求,我需要创建上千个group,它们使用同一组维度,因此将维度信息存储在全局中,只存储一次,每一组中存放不同时期的数据。

写入数据

import netCDF4 as nc
import numpy as np

dataset = nc.Dataset(r'H:\testgroup.nc',mode='w')

dataset.createDimension('lat',10)
dataset.createDimension('lon',10)
dataset.createVariable('lat',np.float32,('lat'))
dataset.createVariable('lon',np.float32,('lon'))
dataset.variables['lat'][:] = np.arange(10)
dataset.variables['lon'][:] = np.arange(10)

for i in range(10):
    group = dataset.createGroup(str(i))
    group.createVariable('water_temp',np.float32,('lat','lon'))
    group['water_temp'][:] = np.random.rand(10,10)
    
dataset.close()

groups

读取数据

import netCDF4 as nc
import numpy as np

dataset = nc.Dataset(r'H:\testgroup.nc')
for group_name in dataset.groups:
    group = dataset.groups[group_name]
    for var_name in group.variables.keys():
        print(group.name,group[var_name][:])
dataset.close()

output

  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值