sp2tmap | R Documentation |
The function converts a SpatialPolygons object for use with the Stata tmap command, by creating a data frame with the required columns.
sp2tmap(SP)
SP |
a SpatialPolygons object |
a data frame with three columns:
_ID |
an integer vector of polygon identifiers in numeric order |
_X |
numeric x coordinate |
_Y |
numeric y coordinate |
and an ID_n
attribute with the named polygon identifiers
Roger Bivand
https://www.stata.com/search.cgi?query=tmap
write.dta
## Not run:
xx <- readShapePoly(system.file("shapes/sids.shp", package="maptools")[1],
IDvar="FIPSNO", proj4string=CRS("+proj=longlat +ellps=clrk66"))
plot(xx, border="blue", axes=TRUE, las=1)
tmapdf <- sp2tmap(as(xx, "SpatialPolygons"))
if (require(foreign)) {
td <- tempdir()
write.dta(tmapdf, file=file.path(td, "NCmap.dta"), version=7)
NCdf <- as(xx, "data.frame")
NCdf$ID_n <- attr(tmapdf, "ID_names")
write.dta(NCdf, file=file.path(td, "NC.dta"), version=7)
}
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.