Tìm cạnh tối thiểu đến khoảng cách cạnh của đa giác bằng ArcGIS Desktop?

Tôi có một bản đồ khoảng 3000 đa giác trong ArcGIS 10. Tôi đang tìm kiếm khoảng cách giữa mỗi người trong số họ. Tôi biết cách thực hiện bằng cách sử dụng tọa độ lat và long của centroid, nhưng tôi đang tìm khoảng cách đường thẳng ngắn nhất từ ​​cạnh gần nhất của một đa giác đến cạnh gần nhất của đa giác khác. Có ý kiến ​​gì không?

Đó là một đoạn mã đẹp, nhưng gần như không đẹp bằng (giả sử bảng của bạn nằm trong tọa độ địa lý, nếu không chỉ xóa các phôi thành địa lý)

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;

Tôi đã đề cập rằng cơ sở dữ liệu không gian đá? Họ làm. Ồ, họ làm.

Khoảng cách từ A đến B giống như B đến A và khoảng cách từ A đến A bằng 0, do đó, một nửa ma trận sẽ giúp bạn tiết kiệm một số công việc.

IProighborOperator trả về khoảng cách từ cạnh. Mã dưới đây sử dụng phép chiếu phương vị tập trung vào tâm của mỗi đa giác (cũng nên hoạt động với các dòng). Nếu các đa giác không quá phức tạp (hoặc nếu bạn có nhiều bộ nhớ) thì việc tải tất cả các hình học vào bộ nhớ, việc chiếu chúng sẽ nhanh hơn. (Điều này không được kiểm tra kỹ lưỡng).

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

Tính toán không phải là thứ mà ArcGIS thực sự cung cấp, nhưng tôi đã đưa ra một giải pháp cho vấn đề trong Shapely có thể giúp ích, hoặc chủ đề này trên các diễn đàn ArcGIS cũ.

Tôi nghĩ rằng công cụ gần bàn sẽ hoạt động cho những gì bạn muốn:

Xác định khoảng cách từ mỗi tính năng trong các tính năng đầu vào đến một hoặc nhiều tính năng gần đó trong các tính năng gần, trong bán kính tìm kiếm. Các kết quả được ghi lại trong bảng đầu ra.

Chỉ cần để trống bán kính tìm kiếm.