Estoy tratando de mapa de accidentes de carretera en Australia, que viene con coordenadas x e y (por ejemplo, x = 1328459.1 y y = 1670878 de Adelaida), así como la siguiente información:
PROJECTION LAMBERT
UNITS METERS
DATUM GDA94 SEVEN /* GDA94 SPHEROID GRS80 PARAMETERS
-28 00 00 /* 1st standard parallel
-36 00 00 /* 2nd standard parallel
135 00 00 /* Central meridian
-32 00 00 /* Latitude of projections origin
1000000 /* False easting (meters)
2000000 /* False northin (meters)
END
El uso de R y esta lista, he tratado de traducir de la siguiente manera
# Projection for accident coordinates (EPSG:3112)
proj <- "+proj=lcc +lat_1=-28 +lat_2=-36 +lat_0=-32 +lon_0=135 +init=epsg:3112 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs"
# Code spatial points
points <- proj4::project(accidents[, c("ACCLOC_X", "ACCLOC_Y")], proj = proj, inverse = TRUE)
accidents$longitude <- points$x
accidents$latitude <- points$y
names(accidents)[names(accidents) == "longitude"] <- "LONG"
names(accidents)[names(accidents) == "latitude"] <- "LAT
antes de la asignación, pero mi enfoque parece ser incorrecta dada la resultante ggplot2-mapa (shapefiles se obtuvieron de los Naturales de la Tierra y crea un subconjunto de a Australia):
¿Alguien tiene una idea de dónde está el error?
