如何使用R指定特定年份的每个月中的天数？

Question

我正在尝试计算月度产品（二进制文件）的年度总数https://echange-fichiers.inra.fr/get?k=Jr3tvgeKnWUiC1B0MMY 。 所以我有12个文件，我必须计算总和。 下面给出的代码工作正常，但我需要的是将每个文件乘以与该年中该月份相对应的天数，然后计算总和。 就我而言，这些数据是针对2000年（a年）的，所以我想将文件编号1的结果乘以31，将文件编号2的结果乘以29，依此类推。

files<- list.files("C:\\jung file_2000_img", "*.img", full.names = TRUE)
x<- do.call(rbind,(lapply(files, readBin  , double() , size = 4 , n=360 * 720 , signed =TRUE)))
results <- colSums(x)
fileName <- sprintf("C:\\annual_ET2000_without999_1.img")
writeBin(as.double(results), fileName, size = 4)

根据Jame的答案，这将计算出天数：

numDays <- function(month,year){
    as.numeric(strftime(as.Date(paste(year+month%/%12,month%%12+1,"01",sep="-"))-1,"%d"))
}

那么现在如何将其应用于上面的代码？ 我的意思是说我想获取每个文件的结果，并将其乘以天数，例如：

file1*numDays(1,2000)+file2*numDays(2,2000)+file3*numDays(3,2000)
   .............file12*numDays(12,2000) 

Answer 1

您可以从下个月的第一天退后一天来获得该月的最后一天。 因此，请尝试以下操作：

numDays <- function(month,year){
  as.numeric(strftime(as.Date(paste(year+month%/%12,month%%12+1,"01",sep="-"))-1,"%d"))
}

测试：

numDays(2,2012)
 [1] 29
numDays(2,2011)
 [1] 28
numDays(1,2011)
 [1] 31
numDays(12,2011)
 [1] 31
# vector of all month days in year 2000
sapply(1:12,numDays,2000)
 [1] 31 29 31 30 31 30 31 31 30 31 30 31

Answer 2

我建议您使用lubridate软件包，该软件包对于此类事情具有许多有用的功能。 例如：

> lubridate::days_in_month(as.Date("2000-02-01"))
Feb 
 29 

Answer 3

@James为您提供了一个函数来查找一个月中的天数（已将其合并到问题中）。 下一步是按天数加权每月值。 问题是您把所有的几个月都放在一起了，无法得知之后是哪个。 因此，不要将它们组合得太多：

files<- list.files("C:\\jung file_2000_img", "*.img", full.names = TRUE)
x<- lapply(files, readBin  , double() , size = 4 , n=360 * 720 , signed =TRUE)
results <- sapply(x, colSums)
ann.results <- sum(sapply(1:12, function(m) {results[[m]] * numDays(m, 2000)}))

这是一个巨大的假设： files[[1]]对应于一月， files[[12]]对应于12月，并且所有匹配项也是如此。

这是做一个数据帧列表，每个月一个。 然后，它使用整个列表中的colSums来创建一个向量，该向量包含每个月的总数（与files相同的顺序）。 最后，将月份迭代为具有天数的相应results值的倍数（再次假定是正确的相应月份）。