1

月間製品(バイナリファイル) https://echange-fichiers.inra.fr/get?k=Jr3tvgeKnWUiC1B0MMYの年間合計を計算しようとしています。したがって、12 個のファイルがあり、合計を計算する必要があります。以下のコードは正常に機能しましたが、必要なのは、各ファイルにその年のその月に対応する日数を掛けてから、合計を計算することです。私の場合、これらのデータは 2000 年 (閏年) のものなので、ファイル番号 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)

ジェームズの答えから、これは日数を計算します：

    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)
4

1 に答える 1

4

ここで以前に使用されたトリックは次のとおりです。翌月の最初を計算し、1日さかのぼって現在の月の最後の日を取得します。

R> seq(as.Date("2000-02-01"), as.Date("2001-01-01"), by="1 month")
[1]  "2000-02-01" "2000-03-01" "2000-04-01" "2000-05-01" "2000-06-01" 
[6]  "2000-07-01" "2000-08-01" "2000-09-01" "2000-10-01" "2000-11-01" 
[11] "2000-12-01" "2001-01-01"
R> seq(as.Date("2000-02-01"), as.Date("2001-01-01"), by="1 month") - 1
[1]  "2000-01-31" "2000-02-29" "2000-03-31" "2000-04-30" "2000-05-31" 
[6]  "2000-06-30" "2000-07-31" "2000-08-31" "2000-09-30" "2000-10-31" 
[11] "2000-11-30" "2000-12-31"

次に、POSIXlt「月の日」フィールドを持つ表現を使用します。

R> as.POSIXlt(seq(as.Date("2000-02-01"), as.Date("2001-01-01"), 
+             by="1 month") - 1)$mday
[1] 31 29 31 30 31 30 31 31 30 31 30 31
R> sum(as.POSIXlt(seq(as.Date("2000-02-01"), as.Date("2001-01-01"), 
+                     by="1 month") - 1)$mday)
[1] 366
R>
于 2012-10-24T17:47:34.647 に答える