Как я могу извлечь NetCDF в CSV для определенного места? - PullRequest
Новая
       16

Как я могу извлечь NetCDF в CSV для определенного места?

0 голосов
/ 08 апреля 2019

Я хочу извлечь временные ряды из NetCDF в csv для определенного местоположения.Я получил это далеко с кодом, но он дает TypeError: требуется байтоподобный объект, а не 'str'

Как я могу преодолеть эту проблему?Кроме того, будет ли этот разработанный код иметь в качестве вывода: время / mwp для определенного местоположения? 

import netCDF4 
import pandas as pd
import matplotlib.pyplot as plt
import csv
import numpy as np
from netCDF4 import Dataset, num2date
from pylab import *
import xarray

f = netCDF4.Dataset('Wave_period_global.nc')
f.variables.keys()
print (f)

lat = f.variables['latitude'][:]
lon = f.variables['longitude'][:]
time_var = f.variables['time']
dtime = netCDF4.num2date(time_var[:],time_var.units)
mwp = f.variables['mwp'][:]

print(lon.min(), lon.max())
print(lat.min(), lat.max())

longitude = 172
latitude = 50

def near(array,value):
    idx=(abs(array-value)).argmin()
    return idx

ix = near(lat, latitude)
iy = near(lon, longitude)

print ('Latitude =',ix)
print ('Longitude =',iy)

with open ('Wave_period_global.csv', 'wb') as csvfile:
    filewriter = csv.writer(csvfile, delimiter=',',
                            quotechar='|', quoting=csv.QUOTE_MINIMAL)
    for ln in range(len(lon)):
        for lt in range(len(lat)):
            value=f.variables['mwp'][0][lt][ln]
            dtime = netCDF4.num2date(time_var[0],time_var.units)
            print(dtime,lat[lt],lon[ln,],value)
            filewriter.writerow([dtime,lat[lt],lon[ln,],value])

1 Ответ

0 голосов
/ 08 апреля 2019

Смотрите мой комментарий выше.Я оставляю вывод как csv из своего ответа, чтобы сосредоточиться на вашем реальном вопросе. 

from netCDF4 import Dataset, num2date

f = netCDF4.Dataset('Wave_period_global.nc')

lat = f.variables['latitude'][:]
lon = f.variables['longitude'][:]
time_var = f.variables['time']
# These are all datetime object.
dtime = netCDF4.num2date(time_var[:],time_var.units)
mwp = f.variables['mwp'][:]

print(lon.min(), lon.max())
print(lat.min(), lat.max())

longitude = 172
latitude = 50

def near(array,value):
    idx=(abs(array-value)).argmin()
    return idx

ix = near(lat, latitude)
iy = near(lon, longitude)

print ('Latitude =',ix)
print ('lat[] =', lat[ix])
print ('Longitude =',iy)
print ('lon[] =', lon[iy])

for i in range(len(dtime)):
  value=f.variables['mwp'][i][ix][iy]
  print(dtime[i],lat[ix],lon[iy,],value)
...