Python-合并具有时间重复的netcdf文件

时间:2019-03-20 17:19:59

标签: python file merge netcdf python-xarray

我有来自MODIS卫星数据的几个netCDF文件。我想合并那些文件。为此,我使用以下命令(在python中):

data = data_1.merge(data_2)

或者我只是使用以下命令打开文件:

data = xr.open_mfdataset('MCD43A3*.nc')

我的问题是1个文件的结尾与下一个文件的开头重复。所以我最终得到的文件是重复的:

lat;lon;time;var1;var2
10;10;2000-01-01;22;55
      2000-01-02;12;87
      2000-01-03;57;65
      2000-01-04;45;67
      2000-01-04;45;67
      2000-01-05;78;15

如何删除时间序列的最后一个时间戳,或者如何合并而不重复?

注意:输入文件具有儒略历,因此,不幸的是CDO / NCO命令不起作用

2 个答案:

答案 0 :(得分:2)

您是否使用CDO命令“ seldate”进行选择?

您也许仍然可以使用seltimestep来做到这一点:

cdo seltimestep,first/last in.nc out.nc 

另一种选择是尝试先使用转换日历

cdo setcalendar,standard in.nc out.nc 

然后使用seldate ...

答案 1 :(得分:1)

我认为您应该可以使用nco。例如如果您有4个时间戳,并且要删除最后一个时间戳,

ncks -d time,0,2 in.nc out.nc

或删除第一个:

ncks -d time,1,3 in.nc out.nc