Python-ArcMAP Точки InsertCursor ArcPy отображаются неправильноPython

Программы на Python
Anonymous
 Python-ArcMAP Точки InsertCursor ArcPy отображаются неправильно

Сообщение Anonymous »

Конечной целью является обновление постоянного ShapeFile с отсутствующими данными скважины. Пока полевые агенты находят отсутствующие скважины, они могут обновить электронную таблицу с помощью GPS -координат (десятичные градусы). Примерно раз в месяц я могу запустить код Python в ArcMap и обновить ShapeFile с новыми точками. Я пытаюсь вставить координаты GPS в качестве точек в файл Shape в ArcMap 10.7, используя Python 2.7. < /P>
Все, кажется, работает правильно, таблица данных и атрибуты верны, вновь добавленные точки кажутся правильными по сравнению друг с другом. Проблема в том, когда я пытаюсь обновить существующий файл. Недавно добавленные очки находятся в другом месте. В примере здесь я попытался добавить те же самые GPS -координаторы (в десятичных градусах), что и существующая точка. Вы можете увидеть первую запись, Y_COORD и X_COORD для WOUL_1 такие же, как 1.
Однако новая точка находится ниже других точек. Это какая-то проблема координат, и поиск по проблеме также указывает на проблему координаты, но у меня есть миллион раз, чтобы убедиться, что ShapeFile, пограничный слой (только визуальная проверка) и DataFrame-все в правильном уровне Географические и проекционные системы координат. Используется GCS IS GCS_WGS_1984 (WKID: 4326), а PCS IS WGS_1984_UTM_ZONE_16N (WKID: 32616). < /P>
Я вывел тройной, тройной, миллион раз. правильный. Я пытался специально назначить его на питоне (хотя я не уверен, что называю это правильно. Как будто я просто называю объектом, но на самом деле не инициирую его.
Если нет, есть ли способ конкретно сказать Python применить это в правильных GCS /PCS? < /p>

Код: Выделить всё

import arcpy

#one coordinate only for testing
coords = [(1, (43.20779, -89.78513))]

# link to shapefile
shapefile = r'M:\ARC\gis_users\Missing_Wells\Missing_Wells.shp'

# Define the spatial reference (GCS_WGS_1984 = 4326, PCS_WGS_1984_UTM16N = 32616)
arcpy.SpatialReference(32616)

#use SHAPE@XY syntax to add point features to a point feature class
with arcpy.da.InsertCursor(shapefile, ['Well_ID', 'Y_Coord', 'X_Coord', 'SHAPE@XY']) as cursor:
#Insert new rows that include the well ID and a x,y coordinate pair
for Well_ID, (lat, lon) in coords:
cursor.insertRow((Well_ID, lat, lon, (lat, lon)))
print('added')
...

sys.exit()
p.s. да, я знаю, что в ArcGIS Online есть инструменты, которые могут достичь этой цели намного проще и быстрее. Наша организация перейдет на ArcGIS Pro в следующем году. На данный момент у нас просто нет на это бюджета.


Подробнее здесь: https://stackoverflow.com/questions/793 ... -correctly

Вернуться в «Python»