0

我有两个具有相同尺寸的二进制文件(光栅)。第一个文件被调用over ,第二个文件是corr. 我想用over大于NA0.5corr的值替换值。

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

1 回答 1

2

确切地看到你想要什么是很棘手的,但是呢:

over[corr > 0.5] = NA

或更复杂的东西:

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