Arcpy – Как изменить пространственную привязку объектов геометрии в пространственно включенном фрейме данных?

Я использую следующий код для экспорта данных из класса пространственных объектов arcgis pro в пространственный фрейм данных. Пространственная привязка экспортируемых данных должна быть 4326.

out_coordinate_system = arcpy.SpatialReference(4326)
    arcpy.env.addOutputsToMap = False
    arcpy.Project_management(fc, fc.name + "_wgs84", out_coordinate_system)
    arcpy.env.addOutputsToMap = True

    #Create a spatially enabled dataframe from featureclass 
    sdf = pd.DataFrame.spatial.from_featureclass(fc.name + "_wgs84")

    #Cleanup temporary layer created for reprojection
    arcpy.Delete_management(fc.name + "_wgs84")

Есть ли способ избежать создания временного класса объектов с помощью инструмента «Проект» и перепроецировать объекты геометрии в SR 4326 в столбце SHAPE объекта sdf с помощью arcpy.geometry?

arcgis arcpy

22.03.2022 08:33

Стоит ли изучать PHP в 2023-2024 годах?

Привет всем, сегодня я хочу высказать свои соображения по поводу вопроса, который я уже много раз получал в своем сообществе: "Стоит ли изучать PHP в...

Поведение ключевого слова "this" в стрелочной функции в сравнении с нормальной функцией

В JavaScript одним из самых запутанных понятий является поведение ключевого слова "this" в стрелочной и обычной функциях.

Приемы CSS-макетирования - floats и Flexbox

Здравствуйте, друзья-студенты! Готовы совершенствовать свои навыки веб-дизайна? Сегодня в нашем путешествии мы рассмотрим приемы CSS-верстки - в...

Тестирование функциональных ngrx-эффектов в Angular 16 с помощью Jest

В системе управления состояниями ngrx, совместимой с Angular 16, появились функциональные эффекты. Это здорово и делает код определенно легче для...

Концепция локализации и ее применение в приложениях React ⚡️

Локализация - это процесс адаптации приложения к различным языкам и культурным требованиям. Это позволяет пользователям получить опыт, соответствующий...

Пользовательский скаляр GraphQL

Листовые узлы системы типов GraphQL называются скалярами. Достигнув скалярного типа, невозможно спуститься дальше по иерархии типов. Скалярный тип...

Перейти к ответу Данный вопрос помечен как решенный

Ответы 1

Ответ принят как подходящий

Удалось заставить его работать с функцией проект_как(), поскольку объекты, хранящиеся в столбце SHAPE файла sdf, имеют тип данных arcgis.geometry objects.

out_coordinate_system = arcpy.SpatialReference(4326)
sdf = pd.DataFrame.spatial.from_featureclass(fc.name)
sdf.dropna(inplace=True)
sdf['SHAPE'] = sdf['SHAPE'].apply(lambda x: x.project_as(out_coordinate_system))

sdf['SHAPE'].apply(lambda x: type(x))
0    <class 'arcgis.geometry._types.Polygon'>
1    <class 'arcgis.geometry._types.Polygon'>
2    <class 'arcgis.geometry._types.Polygon'>
Name: SHAPE, dtype: object

dir(sdf['SHAPE'][0])
['EWKT', 'JSON', 'WKB', 'WKT', '_HASARCPY', '_HASSHAPELY', '__class__', '__contains__', '__delattr__', '__delitem__', '__dict__', '__dir__', '__doc__', '__eq__', '__format__', '__ge__', '__geo_interface__', '__getattr__', '__getattribute__', '__getitem__', '__getstate__', '__gt__', '__hash__', '__init__', '__init_subclass__', '__iter__', '__le__', '__len__', '__lt__', '__module__', '__ne__', '__new__', '__reduce__', '__reduce_ex__', '__repr__', '__setattr__', '__setattribute__', '__setitem__', '__setstate__', '__sizeof__', '__str__', '__subclasshook__', '__weakref__', '_ao', '_check_geometry_engine', '_class_attributes', '_from_geojson', '_repr_svg_', '_shoelace_area', '_type', '_wkt', 'angle_distance_to', 'area', 'as_arcpy', 'as_shapely', 'boundary', 'buffer', 'centroid', 'clear', 'clip', 'contains', 'convex_hull', 'coordinates', 'copy', 'crosses', 'cut', 'densify', 'difference', 'disjoint', 'distance_to', 'envelope', 'equals', 'extent', 'first_point', 'from_shapely', 'fromkeys', 'generalize', 'geoextent', 'geometry_type', 'get', 'get_area', 'get_length', 'get_part', 'has_m', 'has_z', 'hull_rectangle', 'intersect', 'is_empty', 'is_multipart', 'is_valid', 'items', 'keys', 'label_point', 'last_point', 'length', 'length3D', 'measure_on_line', 'overlaps', 'part_count', 'point_count', 'point_from_angle_and_distance', 'pop', 'popitem', 'position_along_line', 'project_as', 'query_point_and_distance', 'rotate', 'scale', 'segment_along_line', 'setdefault', 'skew', 'snap_to_line', 'spatial_reference', 'svg', 'symmetric_difference', 'touches', 'translate', 'true_centroid', 'type', 'union', 'update', 'values', 'within']

22.03.2022 22:34

Другие вопросы по теме

Как отсортировать список имен файлов на основе предопределенного порядка из другого списка после сопоставления шаблонов?

TypeError: ufunc add не может использовать операнды с типами dtype('<M8[ns]') и dtype('<M8[ns]')

Существует ли функция R, эквивалентная FuzzyOverlay( ,"AND") в ArcPy?

Перекрытия всех объектов в одном наборе данных в arcpy (st_overlaps)

Ошибка при использовании geopandas to_crs epsg «столбцы не отображаются»

AttributeError: объект «DataFrame» не имеет атрибута «аллах1__27»

«datetime.datetime» не имеет атрибута «datetime» в python

Не все строки добавляются в список в SearchCursor

Добавить файлы dbf в один файл dbf

Установка Arcpy в среде Anaconda

Arcpy – Как изменить пространственную привязку объектов геометрии в пространственно включенном фрейме данных?

Ответы 1

Другие вопросы по теме

Похожие вопросы