Abrufen von Breiten- und Längengrad des projizierten Punkts mit ArcPy? [geschlossen]

Ich habe ein Punkt-Feature in einer Feature-Class, auf die ArcPy zugreift. Der Punkt ist projiziert, aber ich muss ein effizientes Mittel finden, um den nicht projizierten Breiten- und Längengrad für diesen Punkt zu ermitteln.

Gibt es eine andere Methode als die Neuprojektion (Projektion aufheben), einen Suchcursor für die neue Feature-Class abzurufen, das Feature zu finden und das Lat / Lon von der Form des Features zu entfernen?

Der SearchCursor unterstützt die Angabe eines Raumbezugs. In diesem Fall möchten Sie ein geografisches Koordinatensystem, z. B. WGS 1984, verwenden. Durchlaufen Sie dann den Cursor und nehmen Sie das x & y aus der Form, siehe hier .

Die meisten anderen Antworten wurden veröffentlicht, als ArcGIS 10.0 die neueste Software war. Ab ArcGIS 10.1 sind viele neue ArcPy-Funktionen verfügbar. Diese Antwort nutzt diese neue Funktionalität. Es ist nicht für 10.0 geeignet, bietet jedoch eine höhere Leistung und Funktionalität für 10.1 und höher.

import arcpy

input_feature_class = 'C:\your_feature_class.shp'
wkid = 4326 # wkid code for wgs84
spatial_reference = arcpy.SpatialReference(wkid)

fields_to_work_with = ['SHAPE@']

with arcpy.da.SearchCursor(input_feature_class,
                           fields_to_work_with) as s_cur:
    for row in s_cur:
        point_in_wgs84 = row[0].projectAs(spatial_reference)
        print point_in_wgs84.firstPoint.X, point_in_wgs84.firstPoint.Y

Dieser Code - Schnipsel , die verwendet wkid um ein räumliches Referenzobjekt zu erstellen , anstatt die Eingabe aus einer String - Darstellung, verwendet die modernere Datenzugriff Cursor und Projekte der einzelnen Geometrieobjekte mit Hilfe der projectAs () -Methode .

Um auf James 'Vorschlag näher einzugehen, hier ein minimales Codebeispiel mit Python / arcpy:

import arcpy

def main():
    projectedPointFC = r'c:\point_test.shp'
    desc = arcpy.Describe(projectedPointFC)
    shapefieldname = desc.ShapeFieldName

    rows = arcpy.SearchCursor(projectedPointFC, r'', \
                              r'GEOGCS["GCS_WGS_1984",' + \
                              'DATUM["D_WGS_1984",' + \
                              'SPHEROID["WGS_1984",6378137,298.257223563]],' + \
                              'PRIMEM["Greenwich",0],' + \
                              'UNIT["Degree",0.017453292519943295]]')

    for row in rows:
        feat = row.getValue(shapefieldname)
        pnt = feat.getPart()
        print pnt.X, pnt.Y

if __name__ == '__main__':
    main()

Ganz gleich, ob Sie es Projektion nennen oder nicht, ich bin mir ziemlich sicher, dass Sie per Definition beim Übersetzen der Koordinatenwerte von einem Raumbezugssystem in ein anderes eine Neu- / Neuprojektion durchführen.

Ich bin nicht so vertraut mit ArcPy, aber in ArcGisscripting in Version 9.3 müssten Sie die gesamte Feature-Class projizieren.

Je nachdem, wie komplex ein Projektions- / Transformationsalgorithmus ist, können Sie in der einfachen Python-Mathematik jederzeit eine eigene Projektion für die Koordinaten ausführen. Auf diese Weise können Sie die Wertprojektion auf Feature-Ebene koordinieren.

Wenn Sie offen für die Verwendung der OGR-Python-Bindungen waren, können Sie auf Feature-Ebene in so etwas wie einem "Suchcursor" projizieren.

In ArcPy 10.0 können keine einzelnen Geometrien projiziert werden. Sie können jedoch ein Feature-Set (oder eine speicherinterne Feature-Class) erstellen und dieses anstelle einer vollständigen Feature-Class in einem Arbeitsbereich auf der Festplatte oder in einer Datenbank irgendwo projizieren.

Der Hauptgrund, warum ich sehe, dass ich keine Feature-Class erstellen möchte, ist, dass arcpy.CreateFeatureclass_management langsam sein kann. Sie können auch arcpy.da.NumPyArrayTofeatureClass verwenden, was für in_memory-Feature-Classes mehr oder weniger augenblicklich ist:

In [1]: import arcpy

In [2]: import numpy as np

In [3]: geosr = arcpy.SpatialReference('Geographic Coordinate Systems/Spheroid-based/WGS 1984 Major Auxiliary Sphere')

In [4]: tosr = arcpy.SpatialReference('Projected Coordinate Systems/World/WGS 1984 Web Mercator (auxiliary sphere)')

In [5]: npai=list(enumerate(((-115.12799999956881, 36.11419999969922), (-117, 38.1141))))

In [6]: npai
Out[6]: [(0, (-115.12799999956881, 36.11419999969922)), (1, (-117, 38.1141))]

In [7]: npa=np.array(npai, np.dtype(([('idfield', np.int32), ('XY', np.float, 2)])))

In [8]: npa
Out[8]: 
array([(0, [-115.12799999956881, 36.11419999969922]),
       (1, [-117.0, 38.1141])], 
      dtype=[('idfield', '<i4'), ('XY', '<f8', (2,))])

In [9]: fcName = arcpy.CreateScratchName(workspace='in_memory', data_type='FeatureClass')

In [10]: arcpy.da.NumPyArrayToFeatureClass(npa, fcName, ['XY'], geosr)

In [11]: with arcpy.da.SearchCursor(fcName, 'SHAPE@XY', spatial_reference=tosr) as cur:
    ...:     print list(cur)
    ...:     
[((-12815990.336048, 4316346.515041453),), ((-13024380.422813002, 4595556.878958654),)]

import arcpy
dsc = arcpy.Describe(FC)
cursor = arcpy.UpdateCursor(FC, "", "Coordinate Systems\Geographic Coordinate   Systems\World\WGS 1984.prj")
for row in cursor:
  shape=row.getValue(dsc.shapeFieldName)
  geom = shape.getPart(0)
  x = geom.X
  y = geom.Y
  row.setValue('LONG_DD', x)
  row.setValue('LAT_DD', y)
  cursor.updateRow(row)

del cursor, row