Переезд из списка рассылки gdal-dev: В понедельник, 2 сентября 2013 года, в 19:09 Дэвид Шин написал: Привет список, я пытаюсь упаковать временные ряды растров GTiff с идентичной проекцией / экстентом / разрешением как один файл NetCDF для распространения. Я провел последний час, консультируясь с...
Я читаю шейп-файл с помощью Fiona: with fiona.open('data/boroughs/boroughs_n.shp') as source: mpolys = [p for p in source] candidate = polys[0]['geometry'] это дает мне диктат, содержащий список списков координат, и типа «MultiPolygon». Однако я не уверен, как использовать данные координат для...
У меня есть буферный слой (зеленый многоугольник), который я хочу разделить на два многоугольника всякий раз, когда он пересекает барьер (синяя линия). Я пытался использовать метод "splitGeometry", но я просто не могу заставить его работать. Мой код пока таков: while...
Каков наилучший метод для расчета расстояния от ближайшей дороги для ряда точек адреса в QGIS, имеющих очень большой набор данных? (QGIS 2.0.1, Win XP) Я читаю ответы на предыдущие подобные вопросы: Измерение расстояния между линиями и точками в QGIS Рисование перпендикулярных линий в PyQGIS?...
Я делаю мелкомасштабный рабочий процесс в блокноте ipython с геопандами и фигурными данными, извлекая кучу геопространственных данных, иногда из шейп-файлов, иногда из postgis (где выполняется более дорогая обработка). Теперь я перетаскиваю таблицы postgis в python с помощью...
У меня есть очки в WGS84 широта / длина, и я хотел бы измерить «малые» (скажем, 5 км) расстояния между ними. Я могу использовать формулу haversine из http://www.movable-type.co.uk/scripts/latlong.html, и она работает очень хорошо. Я бы хотел использовать библиотеки Python Shapely, чтобы я мог...
Я пытаюсь создать повернутый прямоугольник в координатах широта / долгота. В принципе, у меня есть координаты центра и угол поворота вокруг центра. Можно ли сделать с помощью, shapely или gdal, или что-то...
Закрыто. Этот вопрос не по теме . В настоящее время он не принимает ответы. Хотите улучшить этот вопрос? Обновите вопрос так это на тему для географической информационной системы Stack Exchange. Закрыто 3 года назад . Я довольно новичок в Python и ищу руководство для вопроса, который может...
Я создаю несколько растров с плавающей точкой в GDAL с драйвером GeoTIFF. Когда я загружаю результирующие изображения в QGIS или Arc, символы по умолчанию устанавливают min на -3.40282e38, а max на 3.40282e38, поэтому растр выглядит точно серым. Есть ли способ, которым я могу записать фактический...
Я использую Ubuntu 14.4, и я установил виртуальную среду Python. Затем я установил несколько библиотек Python (gdal, lxml, psycopg2 и т. Д.). Также в моей системе у меня есть модуль для QGIS. Поэтому, когда я набираю python, я могу импортировать этот модуль как: import qgis.core Но когда я пытаюсь...
Я пишу метод Python, который добавляет два растра и генерирует один выходной растр. По независящим от меня причинам экстенты входных растров отличаются, но они перекрываются. Можно ли работать исключительно с входными растровыми пикселями, которые перекрываются в 2 перекрывающихся областях, чтобы...
У нас есть большой (более 10 000) репозиторий пространственных данных (векторы shp / gdb / tab и растры ecw / tif, а также данные las, xyz и DEM), где у некоторых есть записи метаданных в xml с использованием редактора метаданных anzlic. Я хочу проверить, имеет ли каждый набор пространственных...
Я пытаюсь изучить канвы обработки изображений с помощью дистанционного зондирования с использованием привязок Python GDAL и numpy. В качестве первой попытки я читаю файл геотекса Landsat8, делаю простую манипуляцию и записываю результат в новый файл. Приведенный ниже код работает нормально, за...
На ftp-сайте находится файловая база геоданных, которую я хотел бы скачать с помощью скрипта Python. Сейчас я думаю, что один из способов сделать это - скопировать базу геоданных ftp в базу геоданных на моем компьютере. Ниже приведен скрипт, который я начал. Кто-нибудь знает, как я могу изменить...
Я реализую инструмент, который позволяет пользователю создавать профиль объекта (в моем случае городские участки). Можно ли автоматически сделать снимок экрана этой функции (отображается в окне QGIS) и отобразить его в профиле PDF? На рисунке должны отображаться только размеры функции и, возможно,...
Учитывая ограничивающий прямоугольник, как я могу подсчитать количество определенных значений (скажем, меня интересует количество value == 1) в растре в 1) arcpy или 2) растровом пакете в...
С помощью следующего кода я могу установить 100% прозрачность для растрового значения 0: map=None for lyr in QgsMapLayerRegistry.instance().mapLayers().values(): if lyr.name() == "some_layer": map= lyr tr=None tr=QgsRasterTransparency() tr.initializeTransparentPixelList(0)...
Как только я создал слой, как я могу его скрыть / показать? Я могу включить / отключить рендеринг определенного слоя, установив флажок через QGIS, но мне нужно сделать это программно из кода Python. Как я могу показать / скрыть (не удалить) метку из кода Python? Я ищу что-то вроде: aLayer =...
Как правильно вставить Point в базу данных PostGIS, используя
Я хочу изучить математическую схему, чтобы понять, как можно преобразовать данную широту и долготу в десятичных градусах в EPSG: 3857. Кто-нибудь может указать на хорошую ссылку или возможный исходный код библиотеки Ppython / C / C ++ / Java для проверки...