Я обслуживаю векторные плитки с помощью TileStache , у меня все настроено так, как я хочу. Мои данные хранятся в Postgres, и я использую провайдер VecTiles для обслуживания плиток GeoJSON . Я хочу кэшировать все свои тайлы, чтобы они быстрее работали. Я использую tilestache-seed.py для заполнения...
Как выбрать объекты с помощью PyQGIS с помощью выражения? Я пытался использовать, QgsExpressionно метод выбора не принимает его: exp = QgsExpression("'ogc_fid' = 482") cLayer = canvas.currentLayer() cLayer.select(exp) Возможно ли это, и если да, то как мне это...
Мне удалось добавить QgsVectorLayerв QGISчерез питон плагин, например: vl = self.iface.addVectorLayer(uri.uri(), layerName, self.dbConn.getProviderName()) Но как я могу удалить
У меня есть геопанда, GeoDataFrame содержащая сотни стройныхPolygon и MultiPolygonгеометрических фигур . Полигоны перекрываются во многих местах. Я хотел бы сделать новую геометрию, которая содержит подсчет того, сколько из них перекрывается. Что-то вроде этого: У кого-нибудь есть идеи, как к этому...
Смотрите эту ссылку для более подробной информации. Проблема: Я хочу перебрать непрерывный растр (тот, у которого нет таблицы атрибутов), ячейка за ячейкой и получить значение ячейки. Я хочу взять эти значения и выполнить для них условия, эмулируя шаги алгебры карты, подробно описанные ниже, без...
Есть ли способ, при котором пользовательский плагин, созданный для QGIS, может появляться вместе с оболочкой Python при запуске плагина, как экран отладки? Я спрашиваю это, чтобы узнать, есть ли способ увидеть команды печати через Python Script при программировании плагина для...
В документации по gdal_calc указан растровый калькулятор командной строки с пустым синтаксисом . Позже есть несколько примеров, где в одном из них: gdal_calc.py -A input.tif --outfile = result.tif --calc = "A * (A> 0)" --NoDataValue = 0 - означает установить значения от нуля и ниже до нуля К...
По этому вопросу было множество вопросов, но ни один из них не решил мою проблему. У меня есть машина с: Windows 7 x64 Python 3.4.3 x64 GDAL установлено: gdalinfo --version --> GDAL 1.11.4, released 2016/01/25; GDAL была установлена с колесом GDAL-1.11.4-cp34-none-win_amd64.whlот сюда GDAL...
У меня проблемы с пониманием использования пространственных индексов с RTree. Пример: у меня есть 300 буферизованных точек, и мне нужно знать область пересечения каждого буфера с помощью полигонального шейп-файла. Шейп-файл полигонов имеет> 20000 полигонов. Было предложено использовать...
у меня есть шейп-файл страновых подразделений, и я хотел бы извлечь массив ограничивающих лат и долг для каждого подразделения ... это можно
У меня есть растр, с которым я бы хотел провести точечную интерполяцию. Вот где я нахожусь: from osgeo import gdal from numpy import array # Read raster source = gdal.Open('my_raster.tif') nx, ny = source.RasterXSize, source.RasterYSize gt = source.GetGeoTransform() band_array =...
Я хотел бы создать шейп-файл из фрейма данных Pandas, используя строки lon & lat. У меня есть файл CSV, и я обрабатываю его с помощью панд, чтобы сделать кадр данных, который легче обрабатывать Возможно ли сделать это без создания цикла строка за...
Некоторое время назад я проводил исследование ориентации крыши по соседству для размещения солнечных панелей, визуально проверяя каждую крышу и отмечая по отдельности, ориентированы ли они на юг или почти на юг. Поскольку это занимает много времени, я хотел бы найти способ автоматизировать это, но...
Изучал геообработку с Python, используя курс MSU GEOG 825 (python геообработка). К сожалению, ссылка сейчас не работает. Я знаю, что это длинный путь, но есть ли у кого-нибудь копии материалов этого курса? Я действительно хотел бы закончить это. Уже ознакомились с другими онлайн-курсами...
У меня возникли проблемы с использованием функции Polygonize в Python. Пример поваренной книги для этого можно найти здесь . Соответствующая часть моего кода: sourceRaster = gdal.Open('myraster.tif') band = sourceRaster.GetRasterBand(1) bandArray = band.ReadAsArray() outShapefile = "polygonized"...
Я хотел бы обновить один атрибут функции. Тем не менее, я не знаю, чтобы обновить его с помощью объекта функции. Я должен использовать поставщика данных, чтобы обновить его. layers = QgsMapLayerRegistry.instance().mapLayersByName('my_line') layer = layers[0] dp = layer.dataProvider() it =...
При попытке установить Shapely на мой 64-битный компьютер с Windows я не могу заставить работать библиотеку GEOS. До сих пор я запустил установщик OSGeo4W, из которого я установил GDAL (я считаю, что библиотека geos включена в этот пакет). После этого я проверил и у меня есть geos_c.dll в моем...
У меня есть двоичная сетка Arc / Info - в частности, растр накопления потока ArcGIS - и я хотел бы идентифицировать все ячейки, имеющие определенное значение (или диапазон значений). В конечном счете, мне нужен шейп-файл точек, представляющих эти ячейки. Я могу использовать QGIS, чтобы открыть...
Я ищу метод для вычисления координат X, Y внутренней точки (также известный как labelPoint в arcmap) (неправильного) многоугольника. Я пытаюсь найти решение без использования модуля ArcGIS Разница между «центроидом» и меткой PointPoint заключается в том, что центроид может находиться за пределами...
Есть ли способ выполнить ту же задачу, что и утилита gdalbuildvrt, используя привязки GDAL Python? До сих пор я не нашел никакого другого способа сделать это, кроме как создать vrt одного набора данных и вручную отредактировать xml. Я хотел бы создать VRT из нескольких растров (по сути, выполняя...