擅长:python、mysql、java
<p><a href="https://docs.scipy.org/doc/numpy-1.10.1/reference/generated/numpy.diff.html" rel="nofollow noreferrer">numpy's diff</a>计算连续条目的差异</p>
<p>我怀疑您的文件中有一个多维变量,因此下面是一个通用示例:</p>
<pre><code>import netCDF4
import numpy as np
ncfile = netCDF4.Dataset('./myfile.nc', 'r')
var = ncfile.variables['variable'][:,:,:] # [time x lat x lon]
# Differences with a step of 1 along the 'time' axis (0)
var_diff = np.diff(var, n=1, axis=0)
ncfile.close()
# Write out the new variable to a new file
ntim, nlat, nlon = np.shape(var_diff)
ncfile_out = netCDF4.Dataset('./outfile.nc', 'w')
ncfile_out.createDimension('time', ntim)
ncfile_out.createDimension('lat', nlat)
ncfile_out.createDimension('lon', nlon)
var_out = ncfile_out.createVariable('variable', 'f4', ('time', 'lat', 'lon',))
var_out[:,:,:] = var_diff[:,:,:]
ncfile_out.close()
</code></pre>