python 生成nc文件_如何python写nc文件

今天试了一下用PyNIO写nc文件，难用的一腿

算了，还是用netcdf模块写了。

代码如下：

# -*- coding: utf-8 -*-

import numpy as np
import sys
import os
from netCDF4 import Dataset

def creatspinc(value, filename):

    gridspi = Dataset(filename, 'w', format='NETCDF4')

   # dimensions

    gridspi.createDimension('time', None)

    gridspi.createDimension('lat', 1000)   #len(lat)

    gridspi.createDimension('lon', 2000)

   # Create coordinate variables for dimensions

    times = gridspi.createVariable('time', np.float64, ('time',))

    latitudes = gridspi.createVariable('lat', np.float32, ('lat',))

    longitudes = gridspi.createVariable('lon', np.float32, ('lon',))

    # Create the actual variable

    var = gridspi.createVariable('var', np.float32, ('time', 'lat', 'lon',))

    # Global Attributes

    import time

    gridspi.description = 'var'

    gridspi.history = 'Created ' + time.ctime(time.time())

    gridspi.source = 'netCDF4 python module tutorial'

    # Variable Attributes

    latitudes.units = 'degree_north'

    longitudes.units = 'degree_east'

    times.units = 'days since 2019-01-01 00:00:00'

    times.calendar = 'gregorian'

    # data

    lats = np.linspace(25.,35.,1000) #notice: the last numb is not included

    lons = np.linspace(95.,115.,2000) #notice: the last numb is not included

    latitudes[:] = lats

    longitudes[:] = lons

    #--Fill in values  

    var[0:value.shape[0],:,:] = value

    # Fill in times

    from datetime import datetime, timedelta

    from netCDF4 import num2date, date2num

    dates = []

    for n in range(var.shape[0]):

        dates.append(datetime(2019, 1, 1) + n * timedelta(days=1))

    times[:] = date2num(dates, units = times.units,calendar = times.calendar)

    print 'time values (in units %s): ' % times.units +'\n', times[:]

    dates = num2date(times[:], units=times.units, calendar=times.calendar)

      gridspi.close()

 
    return

out_dir='~/'
out_file_nc = out_dir+'out.nc'
data_nc = np.zeros(shape=(61,1000,2000)) 

creatspinc(data_nc, out_file_nc)

最后这一段是我随便写了个变量调用写.nc文件的函数来测试的。

最近有小伙伴问我为什么公众号文章不能留言，我百度了一下，说这几个月开通的公众号账户文章都没有留言功能，也不知道以后会不会给留言。不过可以赞赏的哟！