중심과 다각형의 가장 먼 지점 사이의 거리

6,00,000 개 이상의 레코드가있는 마을 다각형 레이어가 있습니다. 각 마을의 중심을 계산했습니다. 중심과 각 다각형의 가장 먼 노드 사이의 거리를 찾고 싶습니다. 아래 이미지를 참조하십시오. 검은 선은 다각형 경계입니다.

MapInfo로 작업하는 것 같습니다. "거리 계산기"(도구 관리자에서 찾기)에는 많은 옵션이 있으며이 작업도 처리 할 수 ​​있다고 생각합니다. 여기에 대한 기사가 있습니다 : http://web.pb.com/mapinfopro-jul-2013/Toolbox-Distance-Calculator

PostGIS를 사용하여 ST_ConvexHull을 사용하여 다각형을 단순화하여 더 빠른 결과를 얻었습니다.

가장 먼 점을 얻으십시오 :

SELECT Villages_v4_Trial_region.geom as FarPoint from (
SELECT ST_PointN(ST_ExteriorRing(ST_ConvexHull(Villages_v4_Trial_region.geom)),
generate_series(1, ST_NPoints(ST_ExteriorRing(ST_ConvexHull(Villages_v4_Trial_region.geom))))) as points, 
geom
FROM Villages_v4_Trial_region
ORDER BY ST_MaxDistance(points,ST_Centroid(Villages_v4_Trial_region.geom)) DESC
LIMIT 1;

그리고 중심에서 원을 만드는 데 관심이 있다면 :

SELECT ST_Buffer(Center,ST_Distance(Center,FarPoint)) as Circle
FROM (
SELECT Villages_v4_Trial_region.geom as FarPoint, Center from (
    SELECT ST_PointN(ST_ExteriorRing(ST_ConvexHull(Villages_v4_Trial_region.geom)),
    generate_series(1, ST_NPoints(ST_ExteriorRing(ST_ConvexHull(Villages_v4_Trial_region.geom))))) as points,
    ST_Centroid(Villages_v4_Trial_region.geom) as Center, 
    geom
    FROM Villages_v4_Trial_region
    ) as Villages_v4_Trial_region
    ORDER BY ST_MaxDistance(points,Center) DESC
    LIMIT 1) as foo;

다음 PyQGIS 코드 사용 :

from math import sqrt

layer = iface.activeLayer()

feats = [ feat for feat in layer.getFeatures() ]

n = len(feats)

centroids = [ feat.geometry().centroid().asPoint() for feat in feats ]
polygons = [ feat.geometry().asPolygon()[0] for feat in feats ]

lengths = []

for i, pol in enumerate(polygons):
    max_dist = 0
    idx_j = 0
    for j, point in enumerate(pol):
        dist = sqrt(centroids[i].sqrDist(point))
        if dist > max_dist:
            max_dist = dist
            idx_j = j
    print i, idx_j, max_dist
    lengths.append([centroids[i], pol[idx_j]])

crs = layer.crs()
epsg = crs.postgisSrid()

uri = "LineString?crs=epsg:" + str(epsg) + "&field=id:integer""&index=yes"

mem_layer = QgsVectorLayer(uri,
                           'max_distance',
                           'memory')

prov = mem_layer.dataProvider()

feats = [ QgsFeature() for i in range(n) ]

for i, feat in enumerate(feats):
    feat.setAttributes([i])
    feat.setGeometry(QgsGeometry.fromPolyline(lengths[i]))

prov.addFeatures(feats)

QgsMapLayerRegistry.instance().addMapLayer(mem_layer)

이 shapefile (11 개의 기능이 있음) :

폴리 라인이 중심과 각 다각형의 가장 먼 지점 사이의 거리 (기능) 인 메모리 레이어가 있습니다. 다음 이미지에서 볼 수 있습니다.

QGIS의 Python Console에서는 중심점으로부터의 거리가 최대이고 마지막으로 최대 거리 인 피쳐의 인덱스 인 피쳐 인덱스도 인쇄되었습니다.

MapInfo를 사용하고있는 것처럼 여기에 내가 작업했던 사내 도구에 대해 잠시 전에 쓴 MapBasic 함수가 있습니다. 소스 노드 (중심점)와 영역 객체 (다각형)를 인수로 사용하여 소스 점에서 다각형의 가장 먼 노드에있는 포인트 객체를 반환합니다.

Function GetFurthest(ByVal oNode1 as Object, ByVal oObj as Object) as Object

Dim sourceE,sourceN,East,North,Longest,Dist as Float,
    nNodes,nPolys,i,j as SmallInt,
    oNode2 as Object

    sourceE = CentroidX(oNode1)
    sourceN = CentroidY(oNode1)
    Longest = 0

    nPolys = ObjectInfo(oObj,OBJ_INFO_NPOLYGONS)
    For i = 1 to nPolys
        nNodes = ObjectInfo(oObj,OBJ_INFO_NPOLYGONS+nPolys)
        For j = 1 to nNodes
            East = ObjectNodeX(oObj,i,j)
            North = ObjectNodeY(oObj,i,j)
            Dist = Distance(sourceE,sourceN,East,North,"m")
            If Dist > Longest then
                Longest = Dist
                oNode2 = CreatePoint(East,North)
            End if
        Next
    Next

    GetFurthest = oNode2

End Function