Как построить прямоугольник / ограничивающий прямоугольник в R? - PullRequest
Новая
       14

Как построить прямоугольник / ограничивающий прямоугольник в R?

0 голосов
/ 19 ноября 2018

Как мне сгенерировать ящик, а не форму песочных часов, которую возвращает этот код? Я пытался изменить порядок координат, но у меня всегда была форма «х», а не прямоугольник, иногда в виде боковых песочных часов. 

library(sp)
mybb <- cbind(x=c(363498.5,480497.5,363498.5, 480497.5,363498.5), y=c(5894630,5894630,5806524,5806524,5894630))
crs <-CRS("+proj=utm +zone=12 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0")
mybb <- SpatialPolygons(list(Polygons(list(Polygon(mybb)),"1")), proj4string=crs)
plot(mybb)

Возвращает:

Coordinates not plotted as a rectangle/box

1 Ответ

0 голосов
/ 19 ноября 2018

Не уверен, почему перестановка координат не сработала до сих пор. Это сработало для меня. 

library(rgdal)
library(sp)

mybb <- cbind(x=c(363498.5, 480497.5, 480497.5, 363498.5), y=c(5894630, 5894630, 5806524, 5806524))
crs <-CRS("+proj=utm +zone=12 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0")
mybb <- SpatialPolygons(list(Polygons(list(Polygon(mybb)),"1")), proj4string=crs)
plot(mybb)
...