0

私は一種の新しいRなので、愚かな質問かもしれませんが、自分で理解できません。

これは私の問題です。同じグリッドサイズの複数の asc ファイルがあり、同じ領域をカバーしています。すべての asc ファイルから各グリッドの最大値を取得したいと考えています。私は複数のことを試しました：

   for (i in 1:144){
    asc0<-rasterToPoints(raster(asc0))
        asc1<-rasterToPoints(raster(asc[i]))
        asc0[,3] <-pmax(asc0[,3], asc1[,3])
    }

これは、ファイルをループしてスローすると失敗します。これは、NA が除外されるためです。そのため、asc0 (ベース ファイル) は次のファイル asc1[2] とは異なるサイズになります。

誰もこれへの方法を知っていますか? すべてのファイルを通過する準備ができているループがあり、13x144 ファイルがあります。しかし、最大値を取得して保存し、次のファイルと比較する方法がわかりません。

ありがとう、本当にあなたの助けに感謝します!!!

4

1 に答える 1