如何從netcdf文件中減去連續的每日數據及其后繼日期？

Question

我想用第二個日期減去一個月的第一個每日日期，而不是用第三個日期減去第二個日期數據，依此類推，在R.213-Lon，78-Lat，30天數據中使用DIff。 我的輸出為零。

我的代碼開始：

cdf <-nc_open（file.choose（））//打開我的netcdf文件

rain1 <-ncvar_get（cdf，“ RAINNC”）// extrct變量

暗（rain1）

[1] 213 78 30

obs1 <-array（0，c（213,78,30））//用於輸出的虛擬數組

否<-array（0，c（213,78,30））

for（i in 1：dim（rain1）[1]）{

     for( j in 1 :dim(rain1)[2]){

      for( k in 1 :dim(rain1)[3]){

            if(!is.na(rain1)[3]){

              obs[i,j] <- diff(k ,lag = 1)}

             else { no[i,j] <- NA }
}}}

Answer 1

你可以用這樣的東西

# generate data
my_dates <- as.Date(c("2019-01-01", "2019-01-02", "2019-01-05"), "%Y-%m-%d")
# substract each date data from the preceding one
my_dates[1:(length(my_dates)-1)] - my_dates[2:length(my_dates)]
# see output
Time differences in days
[1] -1 -3

如果您需要其他方法（從以下數據中減去日期數據），您當然可以簡單地將順序更改為my_dates[2:length(my_dates)] - my_dates[1:(length(my_dates)-1)]