首页 > 解决方案 > Fortify and tidy (from broom) 导致多边形不正确

问题描述

我正在尝试使用莱索托地图并在其上绘制各种其他数据。我提取的 shapefile 是来自此链接的 ESRI 文件:http: //www.mapmaker.com/maplibrary/library/stacks/Africa/Lesotho/index.htm

当我通过 readOGR 导入 shapefile,然后使用基本绘图功能绘制它时,多边形处于适当的方向。一旦我使用 fortify 或 tidy,它们就不再处于正确的方向。我已经包含了下面的代码。请注意,我只包含了下面的整洁代码,但也尝试使用 fortify 解决同样的问题。

我已经尝试指定 proj4string(在这个形状文件中很好)并且我尝试了 spTransform,但是当我使用 ggplot 绘图时,没有一个正确的区域方向。我不确定要包括哪些其他信息可能会有所帮助......我在 R 3.4.4 中。任何和所有的帮助表示赞赏!

s1 <- readOGR(".","LSO_adm1")
summary(s1)
plot(s1) #This is actually what the districts in Lesotho look like
s2 <- tidy(s1)
ggplot()+
  geom_polygon(data = s2, aes(x=lat, y=long, group=group), col="black") 
#This is then some sort of weird alternative Lesotho

标签: rggplot2shapefilemap-projectionsbroom

解决方案


推荐阅读