Advertisement
Guest User

Untitled

a guest
Jul 28th, 2016
76
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.45 KB | None | 0 0
  1. # Define the coordinate systems
  2.  
  3. wgs_co <- data.frame(lon=x, lat=y)
  4.  
  5. coordinates(wgs_co) <- c("lon", "lat")
  6.  
  7. proj4string(wgs_co) <- CRS("+init=epsg:4326") # WGS 84
  8.  
  9.  
  10. #convert coordinates
  11.  
  12. UTM60H <- "+proj=utm +zone=60 +datum=NAD83 +units=m +no_defs"
  13.  
  14. UTM60H <- paste(UTM60H, "+ellps=GRS80 +towgs84=0,0,0")
  15.  
  16. UTM60H <- CRS(UTM60H)
  17.  
  18. [4003,] 352172.5 -4445766
  19.  
  20. [4004,] 352174.1 -4445767
  21. ...
  22. ...
  23.  
  24. [5000,] 352657.5 -4445754
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement