Некоторые данные:
require(sp)
require(rgeos)
pts1 = cbind(seq(50,850,100), 50)
pts2 = cbind(seq(50,850,100), 150)
df = data.frame(1:18)
ptss <- rbind(pts1,pts2)
Spts <- SpatialPointsDataFrame(ptss, df)
sPts.B <- gBuffer(Spts, byid=TRUE, width=48, capStyle="SQUARE")
plot(sPts.B)
Результат - это то, что я ищу, своего рода сетка.Построение sPts.B с lwd=4 приводит к квадратам с закругленными углами, но я бы хотел получить квадратные углы.используя параметр построения lend, вы получите один квадратный угол.
Как сделать все углы квадратными?