Desperate with coordinates != WGS84

[hip]

Hi everyone,

I'm freaking out about plotting coordinates which are not in WGS84 standard, can anybody help me how to transform:
Latitude Longitude
1 6459155 1845760
2 6459155 1845827
3 6430265 1853418
4 6427214 1850734
5 6472200 1850005
6 6472165 1850477
7 6472165 1850427
8 ....
into useful GPS positions?

Thanks in advance!

Cheers
Hip

[mfherman]

If you're working with spatial data in R, I'd highly recommend getting to know the sf package.

As for your specific question, the function sf::st_transform() will change your coordinates into WGS84. Below, I show an example of how you would do this. I'm not sure what projection your original data is in, but you can just replace the crs argument of st_as_sf with the EPSG code of your data.

library(sf)
#> Linking to GEOS 3.6.1, GDAL 2.2.3, PROJ 4.9.3

xy_coord <- tibble::tribble(
  ~Latitude, ~Longitude,
   6459155, 1845760,
   6459155, 1845827,
   6430265, 1853418,
   6427214, 1850734,
   6472200, 1850005,
   6472165, 1850477,
   6472165, 1850427
  ) %>% 
  st_as_sf(
    coords = c("Latitude", "Longitude"),
    crs = 2263
    )

xy_coord %>% 
  st_transform(4263)
#> Simple feature collection with 7 features and 0 fields
#> geometry type:  POINT
#> dimension:      XY
#> bbox:           xmin: -53.32455 ymin: 43.44243 xmax: -53.16039 ymax: 43.4814
#> epsg (SRID):    4263
#> proj4string:    +proj=longlat +ellps=clrk80 +towgs84=-92,-93,122,0,0,0,0 +no_defs
#>                     geometry
#> 1 POINT (-53.21212 43.44243)
#> 2  POINT (-53.21207 43.4426)
#> 3  POINT (-53.31103 43.4814)
#> 4  POINT (-53.32455 43.4762)
#> 5 POINT (-53.16068 43.44532)
#> 6  POINT (-53.16039 43.4466)
#> 7 POINT (-53.16043 43.44646)

^{Created on 2019-02-15 by the reprex package (v0.2.1)}

[hip]

Thanks man, I'll try that tomorrow, sounds very good