¿Encuentra la distancia mínima de borde a borde de polígonos usando ArcGIS Desktop?

Tengo un mapa de aproximadamente 3000 polígonos en ArcGIS 10. Estoy buscando encontrar la distancia entre cada uno de ellos. Sé cómo hacerlo usando las coordenadas lat y long del centroide, pero estoy buscando la distancia más corta en línea recta desde el borde más cercano de un polígono hasta el borde más cercano del otro polígono. ¿Algunas ideas?

Es un código agradable, pero no tan bueno como (suponiendo que su tabla esté en coordenadas geográficas, si no solo elimine los modelos a la geografía)

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;

¿He mencionado que las bases de datos espaciales son geniales? Ellas hacen. Oh, ellos lo hacen.

La distancia de A a B es igual a B a A, y la distancia de A a A es cero, por lo tanto, una media matriz le ahorrará algo de trabajo.

IProximityOperator devuelve la distancia desde el borde. El siguiente código utiliza una proyección azimutal centrada en el centroide de cada polígono (también debería funcionar con líneas). Si los polígonos no son demasiado complejos (o si tiene mucha memoria) cargando todas las geometrías en la memoria, proyectarlas sería más rápido. (Esto no se ha probado 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;
    }
}

El cálculo no es algo que ArcGIS proporciona de forma nativa, pero di una solución al problema en Shapely que puede ayudar, o este hilo en los viejos foros de ArcGIS.

Creo que la herramienta de tabla cercana funcionaría para lo que quieres:

Determina las distancias desde cada entidad en las entidades de entrada a una o más entidades cercanas en las entidades cercanas, dentro del radio de búsqueda. Los resultados se registran en la tabla de salida.

Simplemente deje el radio de búsqueda en blanco.