如何叠加地图并添加从多边形到点的最短距离线?

时间:2014-09-15 02:07:42

标签: r

我有一张海岸线地图,并在地图上添加了点数。我试图在海岸线外添加一个包含多边形的Shapefile,然后测量从多边形到地图上各个点的最短距离。地图和示例点如下:

map(database="worldHires", xlim = c(-97.1, -80), ylim = c(25, 32.5), col = "grey", fill=TRUE,     lwd = 2, bg = "lightblue2", mar=rep(0, 4))
#add point
points(-90.0122, 29.2278, pch = 21, col = "blue", bg = "blue", cex=2)

然后我想将Shapefile添加到地图中,但是当我尝试这样做时,我不会收到错误消息,但它也不会显示为地图的附加内容。也许我使用的是错误的套餐?我试过了:

may <- readShapePoly("directory/satellite maps.shp")
plot(may, add=TRUE, col = "red")

我是R的新手,所以如果这是初级的,我会道歉。解决了这个问题之后,我想测量从地图中的点到多边形边缘的最短距离。任何建议将不胜感激!

补充:我尝试了以下内容:

maptest <- readOGR("directory", "satellite maps")#Shapefile without extension
maptest1 <- transform(maptest, CRS("+proj=longlat +datum=WGS84"))
plot(maptest1)

其中&#34; maptest1&#34;是转换的&#34;可能&#34;从上面的问题。这不起作用,情节只显示两点,而不是多边形。如果我试图覆盖&#34; maptest1&#34;在原始的海岸线地图上,我收到了几条关于&#34;添加&#34;的错误信息。不是图形参数。

0 个答案:

没有答案