Auswählen eines Prozentsatzes (75%) einer Gruppe von Punkten basierend auf der Entfernung von einem separaten Punkt in ArcGIS?

Dies ist ArcGIS-spezifisch.

Ich habe 2 Punkt-Shapefiles Aund Bdas erste ( A) ist ein einzelner Punkt, der eine Lat-Länge enthält, das zweite ( B) ist eine Vielzahl von Punkten (über 12 KB), die jeweils ihre Lat- und Long-Punkte enthalten. Ich versuche, die Auswahl von 75% der Shapefile B-Punkte basierend auf dem Abstand zum Shapefile zu automatisieren A. Mit anderen Worten, ich möchte die 75% der Shapefile- BPunkte auswählen, die dem Aeinen Punkt des Shapefiles am nächsten kommen .

Sie können einen Mehrfachringpuffer für Shapefile A erstellen und dann eine räumliche Verknüpfung des Puffers mit Shapefile B durchführen. Wenn Sie eine räumliche Verknüpfung von Polygonen und Punkten durchführen, erhalten Sie eine Anzahl der Punkte in jedem Polygon im Attribut Tabelle des Joins. Wenn Sie dann die Gesamtzahl der Punkte in den Puffern untersuchen, können Sie 75% der Punkte in Shapefile B erreichen.

Ein etwas anderer Ansatz wäre, dies in Python zu skripten und in einer Schleife nach 75% zu suchen. Wenn es sich jedoch um eine einmalige Berechnung handelt, benötigen Sie diese möglicherweise nicht.

Für 1200 Punkte (oder sogar bis zu 12 Millionen Punkte?) Würde ich sie einfach als generische Sammlung speichern - in diesem Fall als sortierte Liste von Listen . Dies könnte vereinfacht werden, indem nur Punkte übersprungen werden, wenn Sie auf eine Situation mit mehreren Punkten stoßen, die denselben Abstand vom Ursprungspunkt haben. Aus Gründen der Leistung sollten Sie auch eine Hashtabelle anstelle einer SortedList verwenden und nach dem Einfügen aller Entfernungen einmal sortieren. Das würde allerdings noch ein paar Zeilen Code erfordern (?).

Ich hatte keine Zeit, dies zu testen, aber dieses c # könnte Ihnen den Einstieg erleichtern:

private void SelectNTile(string layer1, string layer2, double nTile)
{
    var fLayer1 = FindLayer(ArcMap.Document.FocusMap, "LayerWithLotsofPoints");
    var fLayer2 = FindLayer(ArcMap.Document.FocusMap, "LayerWithOneSelectedPoint");
    IFeature feat = GetSingleFeature(fLayer2);
    var distList = MakeDistList(fLayer1.FeatureClass,(IPoint)feat.ShapeCopy);
    // assume not many points exactly same distance
    var nRecs = (int)(distList.Count * nTile); // nTile would be 0.75 for 75%
    var Oids = new List<int>();
    foreach (KeyValuePair<double, List<int>> kvp in distList)
    {
        Oids.AddRange(kvp.Value);
        if (Oids.Count > nRecs)
            break;
    }
    var fSel = fLayer1 as IFeatureSelection;
    var OidArray = Oids.ToArray();
    fSel.SelectionSet.AddList(Oids.Count, ref OidArray[0]);                
}

private SortedList<double, List<int>> MakeDistList(IFeatureClass fc, IPoint pnt)
{
    var outList = new SortedList<double, List<int>>();
    var proxOp = pnt as IProximityOperator;
    IFeatureCursor fCur = null;
    try
    {
        fCur = fc.Search(null, true); // recycling is faster, we just need OIDs
        IFeature feat;
        while ((feat = fCur.NextFeature()) != null)
        {
            double dist = proxOp.ReturnDistance(feat.Shape);
            if (!outList.ContainsKey(dist))
                outList.Add(dist, new List<int> { feat.OID });
            else
                outList[dist].Add(feat.OID);  // this should rarely happen
        }
    }
    catch
    {
        throw;
    }
    finally
    {
        if (fCur != null)
            System.Runtime.InteropServices.Marshal.FinalReleaseComObject(fCur);
    }
    return outList;
}
private IFeature GetSingleFeature(IFeatureLayer fLayer)
{
    var fSel = fLayer as IFeatureSelection;
    if (fSel.SelectionSet.Count != 1)
        throw new Exception("select one feature in " + fLayer.Name + " first");
    var enumIDs = fSel.SelectionSet.IDs;
    enumIDs.Reset();
    IFeature feat = fLayer.FeatureClass.GetFeature(enumIDs.Next());
    return feat;
}
private IFeatureLayer FindLayer(IMap map, string name)
{
    throw new NotImplementedException();
}

Ein Python-Geoverarbeitungsskript ist eine naheliegende Wahl:

1. Verwenden Sie die Punktentfernung Werkzeug , um die Entfernung von Ihren Funktionen in Feature - Class B (die „Eingabe - Features“ Parameter des Werkzeugs) bis zu dem Punkt in Feature - Class A zu berechnen (die „Near - Features“ Parameter des Werkzeugs).
2. Sortieren Sie die Tabelle nach der berechneten Entfernung.
3. Wählen Sie die ersten 75% der Objekt-IDs in der Ausgabetabelle (Spalte "Input_FID") aus und verwenden Sie diese, um Ihre Auswahl aus den ursprünglichen Features in Feature-Class B zu treffen.

Ich hatte dieses Problem vor ein paar Jahren. Ich fand es einfacher, die Daten als "flache Daten" zu halten, alle Daten zu durchlaufen, die Entfernung manuell zu berechnen und dann die besten 75% zu nehmen (ich habe tatsächlich die oberen 10% behalten). Ich habe dann dasselbe in ArcIMS mit ihren Entfernungsberechnungen gemacht und es hat viel länger gedauert.

Das Puffern ist ein enormer Aufwand, aber mathematische Berechnungen sind eine große Stärke. Wenn Sie 12.000 Punkte puffern, werden Sie wahrscheinlich Probleme mit der Leistung haben.