Encontrando distância mínima de ponta a ponta dos polígonos usando o ArcGIS Desktop?

Eu tenho um mapa de cerca de 3000 polígonos no ArcGIS 10. Estou procurando encontrar a distância entre cada um deles. Eu sei como fazê-lo usando as coordenadas lat e long do centróide, mas estou procurando a menor distância em linha reta da borda mais próxima de um polígono à borda mais próxima do outro polígono. Alguma ideia?

Esse é um bom código, mas não tão bom quanto (supondo que sua tabela esteja em coordenadas geográficas, se não apenas remova as transmissões para a 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;

Eu mencionei que os bancos de dados espaciais balançam? Eles fazem. Oh eles fazem.

A distância de A a B é igual a B a A e a distância de A a A é zero, portanto, uma meia matriz poupará algum trabalho.

IProximityOperator retorna a distância da borda. O código abaixo usa uma projeção azimutal centralizada no centróide de cada polígono (também deve funcionar com linhas). Se os polígonos não forem muito complexos (ou se você tiver muita memória), carregar todas as geometrias na memória e projetá-las seria mais rápida. (Isso não foi totalmente testado).

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;
    }
}

O cálculo não é algo que o ArcGIS fornece nativamente, mas dei uma solução para o problema no Shapely que pode ajudar, ou esse tópico nos antigos fóruns do ArcGIS.

Eu acho que a ferramenta Near Table funcionaria para o que você deseja:

Determina as distâncias de cada recurso nos recursos de entrada para um ou mais recursos próximos nos recursos próximos, dentro do raio de pesquisa. Os resultados são registrados na tabela de saída.

Apenas deixe o raio de pesquisa em branco.