我的目标是从netcdf文件访问数据并以下列格式写入CSV文件.

Latitude Longitude Date1 Date2 Date3

100 200

到目前为止,我已经访问了变量,将标题写入文件并填充了lat / lons.

如何访问指定的lon,lat坐标和日期的MIN_SFC值,然后写入CSV文件.

我是一个蟒蛇新手,如果有更好的方法来解决这个问题请告诉我.

NetCDF文件信息：

Dimensions:

time = 7

latitude = 292

longitude =341

Variables:

float MIN_SFC (time=7, latitude = 292, longitude = 341)

这是我尝试过的：

from netCDF4 import Dataset, num2date

filename = "C:/filename.nc"

nc = Dataset(filename, 'r', Format='NETCDF4')

print nc.variables

print 'Variable List'

for var in nc.variables:

print var, var.units, var.shape

# get coordinates variables

lats = nc.variables['latitude'][:]

lons = nc.variables['longitude'][:]

sfc= nc.variables['Min_SFC'][:]

times = nc.variables['time'][:]

# convert date, how to store date only strip away time?

print "Converting Dates"

units = nc.variables['time'].units

dates = num2date (times[:], units=units, calendar='365_day')

#print [dates.strftime('%Y%m%d%H') for date in dates]

header = ['Latitude', 'Longitude']

# append dates to header string

for d in dates:

print d

header.append(d)

# write to file

import csv

with open('Output.csv', 'wb') as csvFile:

outputwriter = csv.writer(csvFile, delimiter=',')

outputwriter.writerow(header)

for lat, lon in zip(lats, lons):

outputwriter.writerow( [lat, lon] )

# close the output file

csvFile.close()

# close netcdf

nc.close()

更新：

我已经更新了写入CSV文件的代码,有一个属性错误,因为lat / lon是双倍的.

AttributeError：’numpy.float32’对象没有属性’append’

有什么方法可以在python中转换为字符串？你认为它会起作用吗？

当我向控制台打印值时,我注意到一些值返回为“ – ”.我想知道这是否代表fillValue或missingValue定义为-32767.0.

我也想知道是否应该通过lats = nc.variables [‘latitude’] [：] [：]或lats = nc.variables [‘latitude’] [：] [：,]来访问3d数据集的变量]？

# the csv file is closed when you leave the block

with open('output.csv', 'wb') as csvFile:

outputwriter = csv.writer(csvFile, delimiter=',')

for time_index, time in enumerate(times): # pull the dates out for the header

t = num2date(time, units = units, calendar='365_day')

header.append(t)

outputwriter.writerow(header)

for lat_index, lat in enumerate(lats):

content = lat

print lat_index

for lon_index, lon in enumerate(lons):

content.append(lon)

print lon_index

for time_index, time in enumerate(times): # for a date

# pull out the data

data = sfc[time_index,lat_index,lon_index]

content.append(data)

outputwriter.writerow(content)