r - Trazar UTM / Coordenadas

CorePress2024-01-24  9

Soy nuevo en la programación R, tengo un archivo csv/excel de 20 ciudades en un país que contiene el siguiente formato,

Pueblos UTM Cordinates UTM Cordinaetes xxxxxxx 1377777 249514 yyyyyyyy 142145 228942

No puedo trazarlos en un mapa. ¿Alguien tiene alguna idea de cómo trazar estas coordenadas UTM? ¿Es posible trazar ciudades en programación R con UTM? Si es así, ¿alguien puede ayudarme aquí? También tengo conmigo el archivo de forma del país. Pero no estoy seguro de cómo procesarlo.

myfilepath <- file.choose()
Cordinates<-read.csv(myfilepath,header=TRUE)
Cordinates
str(Cordinates)
library(rgdal)
library(sp)
library(data.table)
library(ggplot2)
library(tibble)
myfilepath <- file.choose()
Shapefile<-readOGR(myfilepath)
plot(Shapefile)
ggmap(Shapefile)+geom_point(aes(x=Easting,y=Northing,col=Cordinates))

Cualquier ayuda será apreciada.



------------------------------------

Aquí hay un sf-solsolución, aprovechando todo el arduo trabajo de @Dave2e para encontrar el sistema de coordenadas correcto utilizado...

#convert to simple feature
library( sf )
mysf <- sf::st_as_sf( mydata, coords = c("Easting", "Northing"), crs = 29902)

#plot for visual inspection
mapview::mapview(mysf)



------------------------------------

El truco consiste en determinar el sistema de red utilizado. Después de mucha búsqueda, el código para la cuadrícula estándar de la República de Irlanda es epsg:29902

El primer paso es transformar las coordenadas de la cuadrícula irlandesa en latitud y longitud estándar. Esto se logra con el "rgdal" biblioteca.

library(rgdal)    

points <- read.table(header=TRUE, text = "Towns Easting Northing 
                      Belclare 137777 249514
                      Carnmore 142145 228942")

#Pull out the location columns and rename
point <- points[,2:3]
names(point) <-c('x', 'y')

#convert to corrdinates and define initial coordinates systems
coordinates(point) <- c('x', 'y')
proj4string(point)=CRS("+init=epsg:29902")  #29903 is the grid for Ireland

#Transform the Ireland's grid to longitude & latitude
coord<-spTransform(point,CRS("+init=epsg:4326"))
coord

Esto transformará su lista de coordenadas; busque en este sitio cómo trazar en un mapa. Hay muchas opciones disponibles.

1

Muy tarde, lo sé, pero gracias por la increíble explicación.

Usuario15473160

10 de mayo de 2023 a las 18:33

Su guía para un futuro mejor - libreflare
Su guía para un futuro mejor - libreflare