绘制UTM /坐标

Question

我是R编程新手，我有一个csv/excel文件，其中包含了一个国家的20个城镇，其中包含以下格式，

城镇UTM协调UTM坐标

xxxxxxxx 1377777 249514

yyyyyyyy 142145 228942

我无法将它们绘制成地图，是否有人知道如何绘制这些UTM坐标。是否有可能用UTM在R编程中规划城镇？如果是的话有人能帮我一下吗。我也有这个国家的形状文件。但我不知道该怎么处理。

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))

任何帮助都将不胜感激。

Answer 1

这是一个sf-solution，它利用@Dave2e中所有的辛勤工作来找到正确的坐标系统。

#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)

​

​

Answer 2

诀窍是确定所使用的网格系统。经过多次搜索，标准爱尔兰共和国网格的代码是epsg:29902

第一步是将爱尔兰的网格坐标转换为标准的经纬度。这是通过"rgdal“库完成的。

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

这将转换您的坐标列表，请搜索该网站如何绘制成一张地图。有许多可供选择的方案。