2017-04-01 8 views
1

은 위도 경도를 위해 특정 시간 범위에 대한 netCDF의 파일로부터 값을 선택하는 xarray 함수의 netCDF4의 equaivalent 거기 CDO 또는 NCO를 사용하지 않으추출 값 형태의 netCDF는

+1

이미 수행 한 것처럼 xarray에서 작업하지 마십시오. 마지막에 .values를 사용하여 데이터를 np.array로 직접 추출 하시겠습니까? – dreab

답변

1

가 특정 데이터 세트에 따라 차이가있을 것이다, 그러나이 같은이 기능은 특정 좌표와 타임 슬라이스에서 값을 꺼내해야합니다

def variable_pull(site_lat, site_lon, lat, lon, var, time_slice): 
    lat_idx = (np.abs(lat-site_lat)).argmin() 
    lon_idx = (np.abs(lon-site_lon)).argmin() 
    return var[time_slice, lat_idx, lon_idx] 

NetCDF4에서 다음과 비슷한 데이터 세트를 열었습니다.

f = ncdf.Dataset("filepath","r") 
lat = f.variables['lat'][:] 
lon = f.variables['lon'][:] 
var = f.variables['var'][:] 
f.close() 
site_lat = your coordinates of interest 
site_lon = your coordinates of interest 
values = variable_pull(site_lat, site_lon, lat, lon, var, time_slice)