Trovare la distanza minima da bordo a bordo dei poligoni usando ArcGIS Desktop?

Ho una mappa di circa 3000 poligoni in ArcGIS 10. Sto cercando di trovare la distanza tra ognuno di essi. So come farlo usando il lat e le coordinate lunghe del centroide, ma sto cercando la più breve distanza in linea retta dal bordo più vicino di un poligono al bordo più vicino dell'altro poligono. Qualche idea?

Questo è un bel pezzo di codice, ma non così bello come (supponendo che il tuo tavolo sia in coordinate geografiche, se non semplicemente rimuovi i cast in geografia)

CREATE TABLE mytable_distances AS
SELECT a.id, b.id, ST_Distance(a.geom::geography, b.geom::geography) as distance
FROM mytable a, mytable b;

Ho già detto che i database spaziali oscillano? Loro fanno. Oh, lo fanno.

La distanza da A a B è uguale a B ad A e la distanza da A ad A è zero, quindi una mezza matrice ti farà risparmiare un po 'di lavoro.

IProximity L'operatore restituisce la distanza dal bordo. Il codice seguente utilizza una proiezione azimutale centrata sul centroide di ciascun poligono (dovrebbe funzionare anche con le linee). Se i poligoni non sono troppo complessi (o se si ha molta memoria) caricare tutte le geometrie in memoria e proiettarle sarebbe più veloce. (Questo non è stato testato a fondo).

public class Pair
{
    public int Oid1;
    public int Oid2;
    public double Dist;
    public static void TestGetDistances()
    {
        IWorkspaceFactory wsf = new ESRI.ArcGIS.DataSourcesGDB.FileGDBWorkspaceFactoryClass();

        string path = @"C:\Program Files\ArcGIS\DeveloperKit10.0\Samples\data\Usa\USA.gdb";
        var fws = wsf.OpenFromFile(path, 0) as IFeatureWorkspace;
        IFeatureClass fc = fws.OpenFeatureClass("states");
        var halfMatrix = Pair.GetPairs(fc);

    }
    /// <summary>
    /// key is oid of each feature, value is pairs for features with smaller oids.
    /// </summary>
    /// <param name="fc"></param>
    /// <returns></returns>
    public static SortedList<int, List<Pair>> GetPairs(IFeatureClass fc)
    {
        ISpatialReferenceFactory3 srf = new SpatialReferenceEnvironmentClass();
        IProjectedCoordinateSystem pcs = 
        srf.CreateProjectedCoordinateSystem((int)esriSRProjCSType.esriSRProjCS_WGS1984N_PoleAziEqui);

        var outList = new SortedList<int, List<Pair>>();
        IFeatureCursor fCur = fc.Search(null, true);
        IFeature f;
        while ((f = fCur.NextFeature()) != null)
        {
            var pairs = GetDistances(f, pcs);
            Debug.Print("{0} has {1} pairs", f.OID, pairs.Count);
            outList.Add(f.OID, pairs);
        }
        System.Runtime.InteropServices.Marshal.FinalReleaseComObject(fCur);
        return outList;
    }

    private static IPoint GetGCSCentroid(IGeometry geom)
    {
        if (geom.SpatialReference is IProjectedCoordinateSystem)
        {
            geom.Project(((IProjectedCoordinateSystem)geom.SpatialReference).GeographicCoordinateSystem);
        }
        IArea a = geom is IArea ? geom as IArea : geom.Envelope as IArea;
        return a.Centroid;
    }

    /// <summary>
    /// return a list of all other features whose OID is lesser than f1
    /// </summary>
    /// <param name="f1"></param>
    /// <param name="pcs"></param>
    /// <returns></returns>
    private static List<Pair> GetDistances(IFeature f1, IProjectedCoordinateSystem pcs)
    {
        IPoint centroid = GetGCSCentroid(f1.ShapeCopy);

        pcs.set_CentralMeridian(true, centroid.X);
        ((IProjectedCoordinateSystem2)pcs).LatitudeOfOrigin = centroid.Y;
        var g1 = f1.ShapeCopy;
        g1.Project(pcs);

        var outList = new List<Pair>();
        var fc = f1.Class as IFeatureClass;
        var proxOp = g1 as IProximityOperator;
        IFeatureCursor fCur = fc.Search(null, true);
        IFeature f2 = null;
        while ((f2 = fCur.NextFeature()) != null)
        {
            if (f2.OID < f1.OID)
            {
                var g2 = f2.ShapeCopy;
                g2.Project(pcs);
                outList.Add(new Pair()
                {
                    Oid1 = f1.OID,
                    Oid2 = f2.OID,
                    Dist = proxOp.ReturnDistance(g2)
                });
            }
        }
        System.Runtime.InteropServices.Marshal.FinalReleaseComObject(fCur);
        return outList;
    }
}

Il calcolo non è qualcosa che ArcGIS fornisce in modo nativo, ma ho dato una soluzione al problema in Shapely che può aiutare, o questo thread sui vecchi forum ArcGIS.

Penso che lo strumento table vicino funzionerebbe per quello che vuoi:

Determina le distanze da ciascuna funzione nelle funzioni di input a una o più funzioni vicine nelle funzioni vicine, nel raggio di ricerca. I risultati sono registrati nella tabella di output.

Lascia vuoto il raggio di ricerca.