我试图使用rect()函数在地图上绘制一些正方形。所有正方形应该是相同的大小,但由于内部舍入(从坐标中的小数位到屏幕/文件上的整数像素),一些正方形比其他正方形大,有些是垂直矩形,其他水平矩形等。我尝试使用symbols()而不是rect(),但问题仍然存在。
有没有办法确保它们具有相同的尺寸和外观?
这是我的代码:
rx1 <- vector()
rx2 <- vector()
ry1 <- vector()
ry2 <- vector()
gap <- 0.22
for (i in 1:nx) { # for each column
x1 <- xmin + (i-1)*side
x2 <- xmin + i*side
for (j in 1:ny) { # for each row
y1 <- ymin + (j-1)*side # side=0.5 degrees, the size of the cell
y2 <- ymin + j*side
nInd <- spdf$nInd[(i-1)*ny+j] # calculate the number of individuals for that cell
if (nInd >= 1 & nInd < 4) {
rx1[length(rx1)+1] <- x1+gap
ry1[length(ry1)+1] <- y1+gap
rx2[length(rx2)+1] <- x2-gap
ry2[length(ry2)+1] <- y2-gap
}
}
}
rect(rx1,ry1,rx2,ry2,col=NA,border="black")