![](/img/trans.png)
[英]How can I filter netCDF variables with specific conditions in Python?
[英]How can I extract netCDF to csv for a specific location?
我想從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])
請參閱上面的評論。 我將輸出作為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)
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.