如何创建多边形并使其与R中的空间数据协调

Question

I Would like to know if there is a way to set a polygon interactively in a plot spatial data image and get It's coordinates. I know there is the locator() function, but I would like to set a circles or square and get it coordinates.

Answer 1

The question has been updated so maybe I better understand the aim now. You want to plot a region with a lot of boundaries and then click a polygon on top? With the shape file provided in the comment you can do:

library(spatstat)
library(sp)
library(maptools)

obj <- readShapeSpatial("mun_lin.shp")
region <- unmark(as(obj, "psp"))
x11() # Necessary when using RStudio (at least on my machine)
plot(region, main = "")
x <- clickpoly(add = TRUE, col = 2, lwd = 2) # Start locator overlayed on region

Then you just left click the points you need and middle/right click when you are done. The points are stored as a window ( owin ) in x . You can get a list with the x,y coordinates:

vertices(x)

Answer 2

So I had used drawPoly() from the raster library, here is the code:

library(spatstat)
library(raster)
#load the shapefile
plot(shape)#plot it
points(data.coords,col="green",pch=4)#add data coordenates points
polygon<- drawPoly()#draw a polygon in the window
#take it limits
xlim<-polygon@bbox[1,]
ylim<-polygon@bbox[2,]
new_points<-ppp(data.coords[,1], data.coords[,2], xlim, ylim)#take points inside the polygon

new_points<- data.frame(new_points)#convert to data.frame

plot(shape)
points(new_points ,col="red",pch=4)#plot the selected points in the shape