R语言中使用readOGR读取shp文件坐标转换|经纬度值很大

我在绘制比利时海岸线的研究地点地图时,下载了geo.be的shp文件,使用如下代码导入的时候出现了问题:经纬度值非常大

be_map <- readOGR("SiteMap/resources/AD_3_District.shp")

这是由于使用的坐标系不同导致的,因此需要做转换:

be_map <- spTransform(be_map, CRS("+proj=longlat +datum=WGS84"))
be_map <- fortify(be_map)

搞定:

本文采用 CC BY-NC-SA 4.0 创作共享协议,转载请标注以下信息:
【转载自LiuJason's Blog:https://www.liujason.com/article/219.html

说点什么

Please Login to comment
avatar
  Subscribe  
提醒