如何为各国绘制沿海边界

时间:2017-03-08 00:57:32

标签: r gis

我想画国家,突出他们的沿海边界。我有以下技术,几乎有效。

Table1.column1                 Table1.column2             Duplicate?
=====================================================================
ABC DEF                        DEF ABC                    Y
ABC DEF                        GHI ABC                    N
ABCD EFGH IJKL MNOP            IJKL MNOP ABCD EFGH        Y
ABCD EFGH IJKL MNOP            IJKL QRST EFGH ABCD        N
ABC ABC DEF                    DEF ABC DEF                N 

输出:

enter image description here

我想改进两个方面。

首先,我如何摆脱海岸的不连续性(见阿根廷的图像)。是否有一个函数或方法(在rgeos或wherever中)允许模糊确定重叠?

其次,如何使用更高分辨率的地图多边形做类似的事情?

1 个答案:

答案 0 :(得分:0)

您的方法非常好,您只需要使用setScale gIntersection来缩小比例,以便更加宽容:

library(rworldmap)
library(rgeos)
library(maps)

data("countriesCoarseLessIslands")
data("coastsCoarse")

country <- "Argentina"

countryOutline <- countriesCoarseLessIslands[countriesCoarseLessIslands$NAME %in% country,]
maps::map(countryOutline,col="light grey",fill=TRUE,border=0)

setScale(1e+04) # Change Tolerance to avoid dicontinuities

coastalBorder <- gIntersection(coastsCoarse,countryOutline)
plot(coastalBorder, col="blue",lwd=2,add=TRUE)

enter image description here