sp2tmap {maptools} | 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
http://www.stata.com/search.cgi?query=tmap
## 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")) write.dta(tmapdf, file="NCmap.dta", version=7) NCdf <- as(xx, "data.frame") NCdf$ID_n <- attr(tmapdf, "ID_names") write.dta(NCdf, file="NC.dta", version=7) ## End(Not run)