Question

基本上，我想打开一个netcdf文件，读出单个像素的时间戳，然后将时间戳写入一个新文件。这是我的伪代码：

f10 = Dataset(nc_f10, 'r')
Time_UTC_10 = np.transpose(f10.variables['TIME_UTC'][:]) #shape is [92,104]
radiance_10 = f10.variables['RADIANCE'][:] #shape is [92,104]
f10.close()

#Manipulate Radiance Information

#python separates the characters in the timestamp, so join it back up:
for i in np.arange(92):
    for j in np.arange(104):
        joined_16 = ''.join(Time_UTC_16[:,i,j])
        datetime_16[i,j] = datetime.datetime.strptime(joined_16, '%Y-%m-%dT%H:%M:%S.%fZ')

#Create and fill the netcdf 
nc_out = Dataset(output_directory+nc_out_file, 'w', format='NETCDF4')

y = nc_out.createDimension('y',104)
x = nc_out.createDimension('x',92)

times = nc_out.createVariable('time', np.unicode_, ('x','y'))
O5s = nc_out.createVariable('O5s', np.float32, ('x', 'y'))

times[:] = datetime_16
O5s[:] = radiance_10

nc_out.close()

但是当我尝试运行它时，出现以下错误： TypeError：只能将numpy字符串，unicode或对象数组分配给VLEN str var slices

我觉得我可能在这里误解了一些重要的东西。关于如何纠正此代码以将时间戳写入netcdf中的变量的任何想法？

Answer 1

我真的不知道，为什么要将时间变量保留为字符串（这是错误消息所显示的：值可以是字符串，unicode或对象），但是一个示例如下：

#!/usr/bin/env ipython
# ----------------------
import numpy as np
from netCDF4 import Dataset,num2date,date2num
# ----------------------
ny=104;
nx=92
# ----------------------
radiance_10=np.random.random((ny,nx));
datetime_16=np.ones((ny,nx))
# ----------------------
nc_out = Dataset('test.nc', 'w', format='NETCDF4')

y = nc_out.createDimension('y',ny)
x = nc_out.createDimension('x',nx)

times = nc_out.createVariable('time', np.unicode_, ('x','y'))
O5s = nc_out.createVariable('O5s', np.float32, ('x', 'y'))

O5s[:] = radiance_10
for ii in range(ny):
    for jj in range(nx):
        times[jj,ii] = "2011-01-01 00:00:00"

nc_out.close()

基本上，写入时间变量的值现在是在每个网格点“ 2011-01-01 00:00:00”具有值的字符串。

尽管如此，我还是将时间值用作从任意选择的时机开始经过的时间。这是如何在netCDF文件中保留时间的最常用方法。让我们假设我们每时每刻的数据都是2014年4月11日23:59。然后我可以将其保存为2014年4月1日以来的秒数。这是我要使用的代码：

import numpy as np
from netCDF4 import Dataset,num2date,date2num
import datetime
# ----------------------
ny=104;
nx=92
# ----------------------
radiance_10=np.random.random((ny,nx));
# ---------------------------------------------------
timevalue = datetime.datetime(2014,4,11,23,59)
time_unit_out= "seconds since 2014-04-01 00:00:00"
# ---------------------------------------------------
nc_out = Dataset('test_b.nc', 'w', format='NETCDF4')
y = nc_out.createDimension('y',ny)
x = nc_out.createDimension('x',nx)
times = nc_out.createVariable('time', np.float64, ('x','y'))
times.setncattr('unit',time_unit_out);

O5s = nc_out.createVariable('O5s', np.float32, ('x', 'y'))

O5s[:] = radiance_10
times[:] = date2num(timevalue,time_unit_out);
nc_out.close()

如果您检查时间变量中现在的值，则为950340，这是从2014-04-01 00:00到2014-04-11 23:59的秒数。

如何将日期时间信息写入netcdf4？

1 个答案: