Menemukan jarak minimum ujung ke ujung poligon menggunakan ArcGIS Desktop?

Saya memiliki peta sekitar 3000 poligon di ArcGIS 10. Saya mencari untuk menemukan jarak di antara mereka. Saya tahu bagaimana melakukannya dengan menggunakan koordinat lat dan panjang dari centroid, tetapi saya sedang mencari jarak garis lurus terpendek dari tepi terdekat satu poligon ke tepi terdekat dari poligon lain. Ada ide?

Itu adalah sepotong kode yang bagus, tetapi tidak sebagus (dengan asumsi tabel Anda berada dalam koordinat geografis, jika tidak hanya membuang gips ke geografi)

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;

Sudahkah saya menyebutkan bahwa basis data spasial menjadi batu? Mereka melakukannya. Oh, benar.

Jarak dari A ke B sama dengan B ke A, dan jarak dari A ke A adalah nol, oleh karena itu setengah matriks akan menghemat pekerjaan Anda.

IProximityOperator mengembalikan jarak dari tepi. Kode di bawah ini menggunakan proyeksi azimut yang berpusat pada centroid setiap poligon (harus bekerja dengan garis juga). Jika poligon tidak terlalu rumit (atau jika Anda memiliki banyak memori) memuat semua geometri ke dalam memori memproyeksikannya akan lebih cepat. (Ini belum diuji secara menyeluruh).

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

Perhitungannya bukan sesuatu yang secara asli disediakan oleh ArcGIS, tetapi saya memberikan solusi untuk masalah di Shapely yang mungkin membantu, atau utas ini di forum ArcGIS lama.

Saya pikir alat tabel dekat akan bekerja untuk apa yang Anda inginkan:

Menentukan jarak dari setiap fitur dalam fitur input ke satu atau lebih fitur terdekat di fitur dekat, dalam radius pencarian. Hasilnya dicatat dalam tabel output.

Biarkan radius pencarian kosong.