Guest User

Untitled

a guest
Oct 23rd, 2017
72
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.20 KB | None | 0 0
  1. library(rgdal)
  2. library(rgeos)
  3.  
  4. states = readOGR(dsn=path.expand(),layer = "states")
  5. hawaii = states[states$STATE_NAME == "Hawaii", ]
  6. hawaii = spTransform(hawaii, CRS("+proj=utm +zone=4 ellps=WGS84"))
Add Comment
Please, Sign In to add comment