Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- # Define the coordinate systems
- wgs_co <- data.frame(lon=x, lat=y)
- coordinates(wgs_co) <- c("lon", "lat")
- proj4string(wgs_co) <- CRS("+init=epsg:4326") # WGS 84
- #convert coordinates
- UTM60H <- "+proj=utm +zone=60 +datum=NAD83 +units=m +no_defs"
- UTM60H <- paste(UTM60H, "+ellps=GRS80 +towgs84=0,0,0")
- UTM60H <- CRS(UTM60H)
- [4003,] 352172.5 -4445766
- [4004,] 352174.1 -4445767
- ...
- ...
- [5000,] 352657.5 -4445754
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement