我有两个具有相同尺寸的二进制文件(光栅)。第一个文件被调用over
,第二个文件是corr
. 我想用over
大于NA
0.5corr
的值替换值。
读取第一个文件:
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)