Skip to content
Prev 5571 / 29559 Next

Find incircle of a polygon

On Wed, 22 Apr 2009, Etienne Bellemare Racine wrote:

            
You can read (polygon) shapefiles into SpatialPolygon* objects using 
readShapeSpatial in maptools and readOGR in rgdal. Once there, coerce to 
an owin object with methods in maptools, something like:

x <- readOGR(".", "myshape")
ox <- as(as(x, "SpatialPolygons"), "owin")

should do it, but watch for topology checking in spatstat. Spatstat ships 
with a vignette on using shapefiles, see for an online link:

http://cran.r-project.org/web/packages/spatstat/vignettes/shapefiles.pdf

which is up to date, humourous, and should tell you what you need to know. 
Using the advice there should resolve any topology problems. If you are 
stuck after trying things from there, please follow up on this list.

Hope this helps,

Roger