Distanza tra centroide e punto più lontano del poligono

Ho uno strato poligonale di villaggio che ha più di 6.000.000 di record. Ho calcolato il centroide di ogni villaggio. Voglio trovare la distanza tra il centroide e il nodo più lontano di ciascun poligono. Controlla l'immagine qui sotto per riferimento. Le linee nere sono contorni poligonali.

Sembra che tu stia lavorando con MapInfo. Il "Distance Calculator" (lo trovi nel gestore degli strumenti) ha molte opzioni e penso che possa gestire anche questo compito. C'è un articolo a riguardo qui: http://web.pb.com/mapinfopro-jul-2013/Toolbox-Distance-Calculator

Usando PostGIS, ho usato ST_ConvexHull per semplificare il poligono per un risultato più veloce:

Ottieni il punto più lontano:

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;

E se sei interessato a creare un cerchio dal centroide:

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;

Utilizzando il prossimo codice 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)

e questo shapefile (con 11 funzionalità):

Ho ottenuto uno strato di memoria in cui le polilinee erano la distanza tra il centroide e il punto più lontano di ciascun poligono (caratteristica); come si può vedere nella prossima immagine:

Alla Python Console di QGIS, è stato anche stampato l'indice della caratteristica, l'indice del punto nella caratteristica in cui la distanza dal centroide è una distanza massima e, infine, massima.

Sembra che tu stia usando MapInfo, ecco una funzione MapBasic che ho scritto qualche tempo fa per uno strumento interno su cui stavo lavorando. Prende un nodo sorgente (il punto centroide) e un oggetto regione (il poligono) come argomenti e restituisce un oggetto punto nel nodo più lontano nel poligono dal punto sorgente.

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