点特征的距离R中最近的多边形(Distance of point feature to neares

我工作的一个项目的那一刻，我在那里有一个点要素 - 点要素包括142点 - 和多个多边形（约10）。我想计算每一个点和R.最近的多边形要素之间的距离

我目前的做法是乏味的，有点长篇大论。我目前打算计算每一个点和每一个多边形之间的距离。例如，我将计算142点和多边形A之间的距离，所述142点和多边形B，142点和多边形C，等。这里之间的距离之间的距离是这些距离计算中的一个的示例代码：

dist_cen_polya <- dist2Line(centroids_coor, polygonA_shp)

做这些计算后，我会写代码来选择每一个点和最近的多边形之间的最小/最近距离。问题是，这个过程是单调乏味的。

有谁知道一个软件包/代码这将减少计算的工作量/计算时间？我真的想用一个包，比较单一指向最近的面要素或计算点和所有感兴趣的多边形之间的距离？

谢谢。

Answer 1:

在这里，我现在用的是gDistance在rgeos拓扑库函数。我用的是蛮力双环却是出奇的快。它只需不到2秒142分10周的多边形。我相信，有执行循环更elgant方式。

   require(rgeos)

    # CREATE SOME DATA USING meuse DATASET
    data(meuse)
      coordinates(meuse) <- ~x+y
        pts <- meuse[sample(1:dim(meuse)[1],142),]  
    data(meuse.grid) 
      coordinates(meuse.grid) = c("x", "y") 
        gridded(meuse.grid) <- TRUE 
          meuse.grid[["idist"]] = 1 - meuse.grid[["dist"]]    
        polys <- as(meuse.grid, "SpatialPolygonsDataFrame")
          polys <- polys[sample(1:dim(polys)[1],10),]   
    plot(polys)
      plot(pts,pch=19,cex=1.25,add=TRUE)      

    # LOOP USING gDistance, DISTANCES STORED IN LIST OBJECT
    Fdist <- list()
      for(i in 1:dim(pts)[1]) {
        pDist <- vector()
          for(j in 1:dim(polys)[1]) { 
            pDist <- append(pDist, gDistance(pts[i,],polys[j,])) 
          }
        Fdist[[i]] <- pDist
      } 

    # RETURN POLYGON (NUMBER) WITH THE SMALLEST DISTANCE FOR EACH POINT  
    ( min.dist <- unlist(lapply(Fdist, FUN=function(x) which(x == min(x))[1])) ) 

    # RETURN DISTANCE TO NEAREST POLYGON
    ( PolyDist <- unlist(lapply(Fdist, FUN=function(x) min(x)[1])) ) 

    # CREATE POLYGON-ID AND MINIMUM DISTANCE COLUMNS IN POINT FEATURE CLASS
    pts@data <- data.frame(pts@data, PolyID=min.dist, PDist=PolyDist)

    # PLOT RESULTS
    require(classInt)
    ( cuts <- classIntervals(pts@data$PDist, 10, style="quantile") )
       plotclr <- colorRampPalette(c("cyan", "yellow", "red"))( 20 )
         colcode <- findColours(cuts, plotclr)
    plot(polys,col="black")
      plot(pts, col=colcode, pch=19, add=TRUE)

所述min.dist矢量表示多边形的行号。例如，你可以使用此向量作为这样的子集最近的多边形。