Я репроектирую растры в Python, используя GDAL. Мне нужно спроецировать несколько tiff из географических координат WGS 84 в WGS 1984 Web Mercator (Вспомогательная сфера), чтобы потом использовать их в Openlayers вместе с OpenStreetMap и, возможно, картами Google. Я использую Python 2.7.5 и GDAL...
Я часто вижу в коде Python GDAL, что люди закрывают наборы данных в конце своего скрипта. Почему имеет смысл закрывать набор данных в Python GDAL? Есть ли последствия, если я не сделаю это? import gdal # open dataset ds = gdal.Open('test.tif') # close dataset ds = None...
Извиняюсь, если следующий вопрос несколько глуп, но я только ОЧЕНЬ новичок во всей этой ГИС. Я пытаюсь преобразовать некоторые проецируемые изображения GeoTiff в WGS84, используя gdal в python. Я нашел сообщение, в котором описывается процесс преобразования точек в проецируемых геотиффах с...
Учитывая наличие geopandas GeoDataFrameряда полигонов, я хотел бы получить площадь в км кв каждой функции в моем списке. Это довольно распространенная проблема, и обычно предлагаемое решение в прошлом было использовать shapelyи pyprojнапрямую (например, здесь и здесь ). Есть ли способ сделать это в...
Закрыто . Этот вопрос должен быть более сфокусированным . В настоящее время не принимает ответы. Хотите улучшить этот вопрос? Обновите вопрос, чтобы он был сосредоточен только на одной проблеме, отредактировав этот пост . Закрыто 9 месяцев назад . За исключением ArcPy, есть ли какая-либо...
Пакет geopy для Python содержит две методики измерения расстояний: формулы Большого круга и Винсенти . >>> from geopy.distance import great_circle >>> from geopy.distance import vincenty >>> p1 = (31.8300167,35.0662833) # (lat, lon) - https://goo.gl/maps/TQwDd...
В настоящее время я работаю над тем, чтобы импортировать климатические данные CANGRID (предоставленные как Surfer Grid ascii, файлы .grd) в ArcGIS. Сетка размером 95 строк на 125 столбцов. Метаданные обеспечивают широту / долготу происхождения (нижний левый угол), размер ячейки (50 км), а также...
Я работаю с надстройкой ArcMap в C #. Из кода C # я выполнил несколько скриптов Python. Теперь, чтобы запустить этот скрипт, у меня есть жестко запрограммированный путь к Python. Но это не портативно. Итак, я хочу получить путь к исполняемому файлу Python из кода и использовать его. Вопрос: Как я...
Я действительно новичок в Python и хотел бы знать, существует ли быстрый метод получения значений ячеек растрового пикселя за пикселем и координат (координата XY центра каждого пикселя) с использованием Python в ArcGIS 10? Чтобы описать это далее, мне нужно получить карту X, карту Y и значение...
Я хочу программно создать файл GeoJSON, используя полигоны из шейп-файла, но добавляя атрибуты из моего собственного приложения. Это легко сделать для шейп-файла: def create_data_dayer(self,varlist, data): """ Creates a new shape to contain data about nodes. varlist is the list of fields names...
Мне интересно узнать, как использовать массивы NumPy для оптимизации геообработки. Большая часть моей работы связана с «большими данными», где для геообработки часто требуются дни для выполнения определенных задач. Излишне говорить, что я очень заинтересован в оптимизации этих процедур. ArcGIS 10.1...
У меня есть часть программного обеспечения (не Arc), которая работает ночью на непроизводственной машине, которая обновляет функции шейп-файла из внешней базы данных. Периодически я хотел бы удалить все функции в шейп-файле (а не сам файл, который должен остаться) и позволить программному...
Можно ли разделить шейп-файл на элемент в Python? (Лучше всего было бы решение, где я могу временно сохранить результирующие векторные объекты в памяти, а не на диск). Причина: я хочу использовать функцию gdal rasterizeLayer с несколькими различными подмножествами шейп-файлов. Для этой функции...
Я пытаюсь захватить все непересекающиеся полигоны, указанные ниже, используя Shapely (данные полигоны A, B и C). Более того, я надеюсь сделать это без итерации, проверки на пересечение и т. Д. Принятый ответ на этот вопрос выражает метод PostGIS, но может показаться, что «объединение» означает...
Я начал изучать, как манипулировать данными LAS в python, и хотел посмотреть, как другие обрабатывают файлы LAS. Я хотел бы прочитать точки (я использую массив NumPy) и отфильтровать классы 1 и 2 (неклассифицированные и заземленные) в отдельный массив. У меня есть следующий код, но я не могу...
У меня есть список данных, все они имеют одинаковую структуру и являются полигонами. Как я могу сложить их все вместе, используя
Кто-нибудь знает чистую библиотеку Python (без зависимостей C или Java) для выполнения геометрических операций, таких как буфер, объединение, пересечение и т. Д.? Я ищу что-нибудь совместимое с Google...
Даны данные точек с долготой, широтой и значением третьего свойства этой точки. Как я могу кластеризовать точки в группы (географические субрегионы) на основе значения свойства? Я искал по Google и выяснил, что эта проблема, кажется, называется «пространственно ограниченная кластеризация» или...
Попытка импортировать gdal для Python 2.7.3 в Windows XP: >>> import gdal Traceback (most recent call last): File "<stdin>", line 1, in <module> File "C:\Python27\lib\site-packages\gdal.py", line 2, in <module> from osgeo.gdal import deprecation_warn File...
Я ищу другое, более элегантное решение проблемы пространственной статистики. Необработанные данные состоят из координаты xy для каждого отдельного дерева (т.е. преобразуются в точечный файл .shp). Хотя это и не используется в этом примере, каждое дерево также имеет соответствующий многоугольник (то...