如何在文件2中用NA替换值对应于文件1中的指定值?

时间:2013-03-04 08:57:24

标签: r binary raster

我有两个具有相同尺寸的二进制文件(栅格)。第一个文件名为over,第二个文件名为corr。每当over大于0.5时,我都希望NA替换corr中的值。

  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. 阅读第二个文件:

    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)
    

1 个答案:

答案 0 :(得分:2)

确切地看到你想要的东西是很棘手的,但是:

over[corr > 0.5] = NA

或者更复杂的东西:

over[corr < 0.2 | corr > 0.4] = NA