ArcGIS Desktop을 사용하여 다각형의 최소 가장자리 거리를 찾으십니까?

ArcGIS 10에는 약 3000 개의 다각형으로 구성된지도가 있습니다. 각 다각형 사이의 거리를 찾고 싶습니다. 중심의 위도 및 경도 좌표를 사용하는 방법을 알고 있지만 한 다각형의 가장 가까운 가장자리에서 다른 다각형의 가장 가까운 가장자리까지의 최단 직선 거리를 찾고 있습니다. 어떤 아이디어?

그것은 좋은 코드 조각이지만 (거의 지리적으로 캐스트를 제거하지 않는 경우 테이블이 지리적 좌표에 있다고 가정)만큼 좋지는 않습니다.

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;

공간 데이터베이스가 흔들리는 것을 언급 했습니까? 그들이하다. 아, 그래

A에서 B까지의 거리는 B에서 A와 동일하고 A에서 A까지의 거리는 0이므로 반 행렬은 작업을 절약합니다.

IProximityOperator는 가장자리에서 거리를 반환합니다. 아래 코드는 각 다각형의 중심을 중심으로 한 방위각 투영을 사용합니다 (선과 함께 작동해야 함). 다각형이 너무 복잡하지 않거나 메모리가 많은 경우 모든 도형을 메모리에로드하면 투영 속도가 빨라집니다. (이것은 철저히 테스트되지 않았습니다).

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

계산은 ArcGIS가 기본적으로 제공하는 것이 아니지만 Shapely 의 문제 에 대한 해결책을 주었 거나 오래된 ArcGIS 포럼 의이 스레드 에 도움 이 되었습니다.

니어 테이블 도구 는 원하는대로 작동 한다고 생각합니다 .

검색 반경 내에서 입력 형상의 각 형상에서 근거리 형상의 하나 이상의 주변 형상까지의 거리를 결정합니다. 결과는 출력 테이블에 기록됩니다.

검색 반경을 비워 두십시오.