รับละติจูดและลองจิจูดของจุดที่ฉายโดยใช้ ArcPy? [ปิด]

ฉันมีคุณสมบัติจุดในคลาสคุณลักษณะที่ ArcPy เข้าถึงได้ คาดว่าจะมีจุด แต่ฉันต้องหาวิธีที่มีประสิทธิภาพเพื่อให้ได้ละติจูดและลองจิจูดที่ไม่มีโครงการสำหรับจุดนี้

มีวิธีอื่นนอกเหนือจากการปฏิเสธ (ไม่ได้คาดการณ์) รับเคอร์เซอร์ค้นหาในคลาสคุณลักษณะใหม่การค้นหาสถานที่จากนั้นนำ lat / lon ออกจากรูปร่างของสถานที่หรือไม่

SearchCursor สนับสนุนการระบุ reference- เชิงพื้นที่ในกรณีนี้คุณจะต้องการภูมิศาสตร์ระบบพิกัดเช่น WGS 1984 แล้วคุณย้ำผ่านเคอร์เซอร์และคว้า x & Y จากรูปร่างดูที่นี่

คำตอบอื่น ๆ ส่วนใหญ่ถูกโพสต์เมื่อ ArcGIS 10.0 เป็นซอฟต์แวร์ล่าสุด ที่ ArcGIS 10.1 ฟังก์ชั่นใหม่ของ ArcPy นั้นมีวางจำหน่ายแล้ว คำตอบนี้ใช้ประโยชน์จากฟังก์ชั่นใหม่นั้น มันจะไม่เหมาะสำหรับ 10.0 แต่มีประสิทธิภาพและการทำงานที่เพิ่มขึ้นสำหรับ 10.1 และใหม่กว่า

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

ข้อมูลโค้ดนี้ใช้wkidที่จะสร้างวัตถุอ้างอิงเชิงพื้นที่มากกว่าการพิมพ์ออกมาเป็นตัวแทนสตริงใช้ที่ทันสมัยมากขึ้นเคอร์เซอร์การเข้าถึงข้อมูลและโครงการเรขาคณิตแต่ละวัตถุที่ใช้ projectAs วิธี

เพื่ออธิบายรายละเอียดเกี่ยวกับข้อเสนอแนะของ James ต่อไปนี้เป็นตัวอย่างโค้ดขั้นต่ำโดยใช้ 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()

ไม่ว่าคุณจะเรียกมันว่าฉายภาพหรือไม่ก็ตามผมค่อนข้างแน่ใจว่าตามคำนิยามเมื่อคุณกำลังแปลค่าพิกัดจากระบบอ้างอิงเชิงพื้นที่หนึ่งไปยังอีกระบบหนึ่งคุณก็กำลังฉายอีกครั้ง

ฉันไม่คุ้นเคยกับ ArcPy แต่ใน arcgisscripting ที่ 9.3 คุณจะต้องฉายคลาสฟีเจอร์ทั้งหมด

ขึ้นอยู่กับความซับซ้อนของอัลกอริธึมการฉายภาพ / การเปลี่ยนรูปแบบที่คุณต้องการคุณสามารถหมุนการฉายภาพของคุณเองเพื่อหาพิกัดในคณิตศาสตร์ไพ ธ อนขั้นพื้นฐาน สิ่งนี้จะช่วยให้คุณสามารถประสานการประมาณการค่าที่ระดับคุณสมบัติ

หากคุณเปิดให้ใช้การผูก Python ของ OGR คุณสามารถฉายภาพที่ระดับสถานที่ภายในบางสิ่งเช่น 'เคอร์เซอร์ค้นหา'

ที่ ArcPy 10.0 ไม่มีความสามารถในการแสดงรูปทรงเรขาคณิตของแต่ละบุคคล อย่างไรก็ตามคุณสามารถสร้างชุดคุณลักษณะ (หรือคลาสคุณลักษณะในหน่วยความจำ) และโครงการที่แทนคลาสฟีเจอร์เต็มรูปแบบในพื้นที่ทำงานบนดิสก์หรือในฐานข้อมูลบางแห่ง

เหตุผลหลักที่ฉันไม่เห็นว่าต้องการสร้างคลาสคุณลักษณะเป็นเพราะ arcpy.CreateFeatureclass_management อาจช้า คุณยังสามารถใช้ arcpy.da.NumPyArrayTofeatureClass ซึ่งเป็นอินสแตนซ์ที่เพิ่มมากขึ้นหรือน้อยลงสำหรับคลาสคุณสมบัติ in_memory:

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