Abstand zwischen Schwerpunkt und weitestem Punkt des Polygons

Ich habe einen Dorfpolygon-Layer mit mehr als 6,00.000 Datensätzen. Ich habe den Schwerpunkt jedes Dorfes berechnet. Ich möchte den Abstand zwischen dem Schwerpunkt und dem entferntesten Knoten jedes Polygons finden. Überprüfen Sie das Bild unten als Referenz. Schwarze Linien sind Polygongrenzen.

Anscheinend arbeiten Sie mit MapInfo. Der "Entfernungsrechner" (finden Sie im Tool-Manager) hat viele Optionen und ich denke, er kann diese Aufgabe auch bewältigen. Hier gibt es einen Artikel dazu: http://web.pb.com/mapinfopro-jul-2013/Toolbox-Distance-Calculator

Mit PostGIS habe ich ST_ConvexHull verwendet, um das Polygon für ein schnelleres Ergebnis zu vereinfachen:

Holen Sie sich den weitesten Punkt:

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;

Und wenn Sie daran interessiert sind, einen Kreis aus dem Schwerpunkt zu erstellen:

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;

Verwenden des nächsten PyQGIS- Codes:

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)

und dieses Shapefile (mit 11 Features):

Ich habe eine Memory-Ebene erhalten, in der die Polylinien den Abstand zwischen dem Schwerpunkt und dem entferntesten Punkt jedes Polygons (Feature) darstellen. wie es auf dem nächsten Bild zu sehen ist:

In der Python-Konsole von QGIS wurde auch der Merkmalsindex gedruckt, der Index des Merkmalspunkts, bei dem der Abstand zum Schwerpunkt ein Maximum und schließlich der maximale Abstand ist.

Anscheinend verwenden Sie MapInfo. Hier ist eine MapBasic-Funktion, die ich vor einiger Zeit für ein internes Tool geschrieben habe, an dem ich gearbeitet habe. Es verwendet einen Quellknoten (Ihren Schwerpunkt) und ein Bereichsobjekt (das Polygon) als Argumente und gibt ein Punktobjekt am entferntesten Knoten des Polygons vom Quellpunkt zurück.

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