我有两个具有相同尺寸的二进制文件(栅格)。第一个文件名为over
,第二个文件名为corr
。每当over
大于0.5时,我都希望NA
替换corr
中的值。
阅读第一个文件:
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
阅读第二个文件:
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
写下结果:
to.write = file(paste("/orcomplete.bin", sep=""), "wb")
writeBin(as.double(results), to.write, size = 4)
答案 0 :(得分:2)
确切地看到你想要的东西是很棘手的,但是:
over[corr > 0.5] = NA
或者更复杂的东西:
over[corr < 0.2 | corr > 0.4] = NA