R: Convert SpatialPolygons object for Stata tmap command

sp2tmap {maptools} | R Documentation

## Convert SpatialPolygons object for Stata tmap command

### Description

The function converts a SpatialPolygons object for use with the Stata tmap command, by creating a data frame with the required columns.

### Usage

    sp2tmap(SP)

### Arguments

`SP` | a SpatialPolygons object

### Value

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

### Author(s)

Roger Bivand

### References

http://www.stata.com/search.cgi?query=tmap

### See Also

`write.dta`

### Examples

    ## 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)) {
    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)

[Package maptools version 0.9-3 Index]