Я хотел бы использовать скрипт Python, не основанный на arcpy, для выполнения таких задач, как запрос шейп-файла по атрибутам, создание нового слоя из выделения, вычисление областей многоугольника и преобразование многоугольников в точки.
У кого-нибудь есть примеры кода использования других модулей или библиотек Python? Я могу сделать это легко, используя arcpy, но я хотел изучить другие варианты.
python
geoprocessing
шерпы
источник
источник
Ответы:
Это странно, как если бы люди внезапно обнаружили мощь Python (без ArcPy, который является всего лишь одним модулем Python среди других), см., Например, вопрос Визуализация шейп-файла в Python :
Вы можете объединить все (Pysal с shapely, ...) и смешать их с другими научными модулями.
Таким образом, для примеров сценариев Python найдите Pyshp Fiona, ogr, gdal или shapely в gis.stackexchange или в Интернете (много примеров, не только на английском).)
Один из них на французском языке (сценарии и рисунки универсальны!):
- Python: использование векторных и растровых слоев в геологической перспективе, без программного обеспечения ГИС,
другое на английском языке:
- ГИС с Python, Shapely и Fiona
и на испанском языке
- Определение областей неправильных многоугольников с использованием координат вершин
в gis.stackexchange
- Профиль высоты 10 км с каждой стороны линии
- Обновление атрибутов с помощью Pyshp
- Как создать 3D-файл формы из растра?
- Python Script для получения разницы высот между двумя точками и
т. Д.
Сценарий, представленный Аароном, можно написать проще с помощью Fiona, которая использует только словари Python:
и если вы используете стройный в дополнение:
Есть также две книги:
Python Geospatial Development Эрика Вестра.
Изучение геопространственного анализа с помощью Python of Joel Lawhead
Python также используется в качестве языка сценариев в других ГИС-приложениях, таких как QGIS (Quantum GIS), GRASS GIS, gvSIG или OpenJump, или в 3D-моделлерах, таких как Paraview (и Blender !). И вы можете использовать большинство геопространственных модулей во всех этих приложениях (см. Визуализация данных QGIS с помощью Blender )
источник
Я настоятельно рекомендую начать геообработку сайта USU с использованием Python с использованием ГИС с открытым исходным кодом . В основном они используют библиотеку GDAL / OGR на протяжении всего упражнения. Установка GDAL / OGR может быть сложной задачей, поэтому эта запись в блоге может быть вам полезна: Установка GDAL (и OGR) для Python в Windows . Также ознакомьтесь с альтернативами использованию Arcpy на GIS.SE.
Следующий пример сценария геообработки с открытым исходным кодом (с сайта USU) используется для извлечения данных атрибутов и записи их в текстовый файл:
источник
.Destroy
это удивительное имя метода: pВы можете быть заинтересованы в GDAL / OGR .
GDAL используется для обработки растров, а OGR - для векторов. Обе библиотеки с открытым исходным кодом.
Если вы хотите удалить некоторую зависимость от ArcPy, вы можете имитировать некоторые возможности, читая информацию в массив и выполняя ваши собственные вычисления на чистом Python.
Я недавно сделал это с помощью выбора точек в многоугольнике, как показано здесь . Он использует алгоритм наведения лучей, чтобы определить, находится ли точка внутри многоугольника, учитывая координаты вершин многоугольника.
источник
Я никогда не использовал это лично, но другие в офисе любят использовать аккуратно: https://pypi.python.org/pypi/Shapely
источник