Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- library(rgdal)
- library(rgeos)
- states = readOGR(dsn=path.expand(),layer = "states")
- hawaii = states[states$STATE_NAME == "Hawaii", ]
- hawaii = spTransform(hawaii, CRS("+proj=utm +zone=4 ellps=WGS84"))
Add Comment
Please, Sign In to add comment