How to plot rectangle/bounding box in R?
up vote
0
down vote
favorite
How do I generate a box rather than the hourglass shape this code is returning? I've tried changing the order of the coordinates but always have that "x" forming rather than a box, sometimes as a sideways hourglass.
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)
Returns:
r sp
asked Nov 19 at 18:20
Emily
11810
add a comment |
up vote
0
down vote
favorite
How do I generate a box rather than the hourglass shape this code is returning? I've tried changing the order of the coordinates but always have that "x" forming rather than a box, sometimes as a sideways hourglass.
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)
Returns:
r sp
asked Nov 19 at 18:20
Emily
11810
Please start the script with a call tolibrary(pkgname)to load the package where functionsCRS,SpatialPolygons, etc are to be found.
– Rui Barradas
Nov 19 at 18:26
Sorry, I should know that
– Emily
Nov 19 at 18:33
add a comment |
up vote
0
down vote
favorite
up vote
0
down vote
favorite
How do I generate a box rather than the hourglass shape this code is returning? I've tried changing the order of the coordinates but always have that "x" forming rather than a box, sometimes as a sideways hourglass.
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)
Returns:
r sp
asked Nov 19 at 18:20
Emily
11810
How do I generate a box rather than the hourglass shape this code is returning? I've tried changing the order of the coordinates but always have that "x" forming rather than a box, sometimes as a sideways hourglass.
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)
Returns:
r sp
r sp
asked Nov 19 at 18:20
Emily
11810
asked Nov 19 at 18:20
Emily
11810
edited Nov 19 at 18:32
asked Nov 19 at 18:20
Emily
11810
asked Nov 19 at 18:20
Emily
11810
asked Nov 19 at 18:20
Emily
11810
11810
Please start the script with a call tolibrary(pkgname)to load the package where functionsCRS,SpatialPolygons, etc are to be found.
– Rui Barradas
Nov 19 at 18:26
Sorry, I should know that
– Emily
Nov 19 at 18:33
add a comment |
Please start the script with a call tolibrary(pkgname)to load the package where functionsCRS,SpatialPolygons, etc are to be found.
– Rui Barradas
Nov 19 at 18:26
Sorry, I should know that
– Emily
Nov 19 at 18:33
Please start the script with a call to
library(pkgname) to load the package where functions CRS, SpatialPolygons, etc are to be found.– Rui Barradas
Nov 19 at 18:26
Please start the script with a call to
library(pkgname) to load the package where functions CRS, SpatialPolygons, etc are to be found.– Rui Barradas
Nov 19 at 18:26
Sorry, I should know that
– Emily
Nov 19 at 18:33
Sorry, I should know that
– Emily
Nov 19 at 18:33
add a comment |
1 Answer
1
active
oldest
votes
up vote
1
down vote
accepted
Not sure why rearranging the coordinates hasn't worked so far. This worked for me.
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)
answered Nov 19 at 18:58
Sarah
656
Thanks! For others, this is the order coordinate pairs should be written, where TL = top left, RB = right bottom etc. referring to vertices of a rectangle: 1) TL, 2) TR, 3) BR, 4) BL, optionally 5) TL to close (clockwise)
– Emily
Nov 19 at 19:07
add a comment |
1 Answer
1
active
oldest
votes
1 Answer
1
active
oldest
votes
active
oldest
votes
active
oldest
votes
up vote
1
down vote
accepted
Not sure why rearranging the coordinates hasn't worked so far. This worked for me.
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)
answered Nov 19 at 18:58
Sarah
656
Thanks! For others, this is the order coordinate pairs should be written, where TL = top left, RB = right bottom etc. referring to vertices of a rectangle: 1) TL, 2) TR, 3) BR, 4) BL, optionally 5) TL to close (clockwise)
– Emily
Nov 19 at 19:07
add a comment |
up vote
1
down vote
accepted
Not sure why rearranging the coordinates hasn't worked so far. This worked for me.
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)
answered Nov 19 at 18:58
Sarah
656
Thanks! For others, this is the order coordinate pairs should be written, where TL = top left, RB = right bottom etc. referring to vertices of a rectangle: 1) TL, 2) TR, 3) BR, 4) BL, optionally 5) TL to close (clockwise)
– Emily
Nov 19 at 19:07
add a comment |
up vote
1
down vote
accepted
up vote
1
down vote
accepted
Not sure why rearranging the coordinates hasn't worked so far. This worked for me.
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)
answered Nov 19 at 18:58
Sarah
656
Not sure why rearranging the coordinates hasn't worked so far. This worked for me.
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)
answered Nov 19 at 18:58
Sarah
656
answered Nov 19 at 18:58
Sarah
656
answered Nov 19 at 18:58
Sarah
656
answered Nov 19 at 18:58
Sarah
656
656
Thanks! For others, this is the order coordinate pairs should be written, where TL = top left, RB = right bottom etc. referring to vertices of a rectangle: 1) TL, 2) TR, 3) BR, 4) BL, optionally 5) TL to close (clockwise)
– Emily
Nov 19 at 19:07
add a comment |
Thanks! For others, this is the order coordinate pairs should be written, where TL = top left, RB = right bottom etc. referring to vertices of a rectangle: 1) TL, 2) TR, 3) BR, 4) BL, optionally 5) TL to close (clockwise)
– Emily
Nov 19 at 19:07
Thanks! For others, this is the order coordinate pairs should be written, where TL = top left, RB = right bottom etc. referring to vertices of a rectangle: 1) TL, 2) TR, 3) BR, 4) BL, optionally 5) TL to close (clockwise)
– Emily
Nov 19 at 19:07
Thanks! For others, this is the order coordinate pairs should be written, where TL = top left, RB = right bottom etc. referring to vertices of a rectangle: 1) TL, 2) TR, 3) BR, 4) BL, optionally 5) TL to close (clockwise)
– Emily
Nov 19 at 19:07
add a comment |
draft saved
draft discarded
Thanks for contributing an answer to Stack Overflow!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.
Some of your past answers have not been well-received, and you're in danger of being blocked from answering.
Please pay close attention to the following guidance:
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.
draft saved
draft discarded
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f53380528%2fhow-to-plot-rectangle-bounding-box-in-r%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Please start the script with a call to
library(pkgname)to load the package where functionsCRS,SpatialPolygons, etc are to be found.– Rui Barradas
Nov 19 at 18:26
Sorry, I should know that
– Emily
Nov 19 at 18:33