我有以下数据框:
crs Lai_500m FparLai_QC FparExtra_QC Lai_ref
lat lon month
42.152083 13.356250 1 -127.0 0.20 8.0 160.0 0.50
2 -127.0 0.10 8.0 160.0 0.50
3 -127.0 0.20 10.0 160.0 0.70
4 -127.0 0.40 8.0 160.0 1.40
5 -127.0 0.70 2.0 160.0 2.10
6 -127.0 1.10 0.0 160.0 1.70
7 -127.0 1.10 0.0 128.0 1.00
8 -127.0 1.00 0.0 128.0 0.90
9 -127.0 0.80 2.0 136.0 0.70
10 -127.0 0.70 2.0 136.0 0.95
11 -127.0 0.50 2.0 136.0 0.70
12 -127.0 0.50 5.0 160.0 0.70
我想将此数据帧转换为netcdf文件。为此,我在使用以下代码:
from netCDF4 import Dataset
import xarray as xr
xri = xr.Dataset.from_dataframe(df)
xri.to_netcdf('test.nc')
但是当我想检查netCDF文件时,变量Lai_500m不正确。
dataset = Dataset("test.nc")
print dataset.variables['Lai_500m'][42.152083, 13.356250,:]
[0.5 0.5 0.6 1. 3.25 4.9 4.95 5.1 4.6 1.4 0.7 0.6]
有人知道我在哪里做错吗?