0

同じ次元の 2 つのバイナリ ファイル (ラスター) があります。最初のファイルは呼び出されover 、2 番目のファイルはcorr. が 0.5 より大きいときはいつでも値overを置き換えたいです。NAcorr

  1. 最初のファイルを読み取ります。

    conn <- file("C:\\corr.bin", "rb")
    over <- readBin(conn, numeric(), size=4,  n=1440*720, signed=TRUE)
    y <- t(matrix((data=corr), ncol=720, nrow=1440))## binary
    r1 <- raster(t(matrix((data=over), ncol=720, nrow=1440)))## raster
    
  2. 2 番目のファイルを読み取ります。

    conne <- file("C:\\cor06.bin", "rb")
    corr <- readBin(conne, numeric(), size=4,  n=1440*720, signed=TRUE)
    y1 <- t(matrix((data=corr), ncol=720, nrow=1440))## binary
    r <- raster(t(matrix((data=corr), ncol=720, nrow=1440)))##raster
    
  3. 結果を書きます:

    to.write = file(paste("/orcomplete.bin", sep=""), "wb")
    writeBin(as.double(results), to.write, size = 4)
    
4

1 に答える 1

2

あなたが望むものを正確に見るのは難しいですが、どうですか:

over[corr > 0.5] = NA

またはもう少し複雑なもの:

over[corr < 0.2 | corr > 0.4] = NA
于 2013-03-04T08:59:30.733 に答える