У меня есть большое количество шейп-файлов, представляющих интересующие области для анализа, который будет проводиться с использованием различных источников спутниковых изображений (IKONOS, RapidEye и т. Д.). К сожалению, на снимках не используется система траекторий, как, например, Landsat,...
Новичок здесь борется с ГИС. Я пытаюсь наметить приходы для города Милуаки, используя шейп-файлы, найденные на их веб-сайте округа . Я следую за темой здесь с некоторым успехом. Мой код дает: from pyproj import Proj, transform # wisconsing EPSG:32054 # epsg:4326 is for the entire world, wgs...
Я пытаюсь узнать, как использовать OGR в Python, используя наборы данных о стране и населенных пунктах из http://www.naturalearthdata.com/downloads/50m-cultural-vectors/, Я пытаюсь использовать фильтры и буферы для поиска точек (ne_50m_populated_places.shp) в указанном буфере названной страны...
Я использую алгоритм Берча из пакета Python scipy-learn для кластеризации набора точек в одном маленьком городе по 10. Я использую следующий код: no = len(list_of_points)/10 brc = Birch(branching_factor=50, n_clusters=no, threshold=0.05,compute_labels=True) По моей идее, я бы всегда заканчивал...
Я пытаюсь создать полигоны для спутниковых орбитальных полос. Пока у меня есть метод для генерации двух линий, которые представляют край каждой полосы в [lat, long]. Некоторые полосы пересекли международную линию дат и так обернулись: Я смог решить это с ogr2ogr -wrapdateline: ogr2ogr -wrapdateline...
Я использую следующую строку, чтобы сообщить пользователю о статусе: iface.mainWindow().statusBar().showMessage("Status:" + str(i)) Плагину требуется около 2 минут для запуска моего набора данных, но Windows обнаруживает его как «не отвечающий» и перестает показывать обновления статуса. Для нового...
Я хочу создать функцию, которая выбирает функцию и увеличивает ее (аналогично QGIS). Поэтому есть следующая функция: QgsMapLayerRegistry.instance().addMapLayer(self.vlayer) def zoomTo(self): layer = self.vlayer atable = self.ui.table selectList=[] for i in atable.selectionModel().selectedRows(): ID...
Я пытаюсь создать шейп-файл полигона из списка координат, используя инструменты Python с открытым исходным кодом. Следующее - это то, что я до сих пор получил, взломанный вместе с кулинарной книгой Python GDAL / OGR, и этот ответ GIS SE . Есть похожий вопрос Python: Как создать шейп-файл полигона...
Можно ли изменить цвет при добавлении векторного слоя? У меня есть векторный слой с областями воды и я хочу загрузить его в стиле синего цвета. Слой добавлен с: QgsMapLayerRegistry.instance().addMapLayer(self.vlayer) Можно ли как-то изменить цвет после или во время...
Я нашел данные высотных отметок (1/3 - 1/9 угловых секунд) из набора данных национальных высот, предоставленного USGS. Он поставляется в форматах IMG, Arcgrid и GridFloat. Я знаю, что ArcGrid соответствует платному программному пакету, но я стараюсь использовать свободно доступные инструменты. У...
Я пытаюсь создать плагин для загрузки композитора печати из файла, создания атласа и экспорта в изображение. До сих пор мне удалось загрузить шаблон и экспортировать его в изображение. Мне не удалось добавить ни один из слоев в легенде (которые также находятся в toc) на экспортированную карту, что...
Мне было поручено создать модель с использованием ArcGIS ModelBuilder, которая сначала выполняет несколько задач геообработки, а затем соединяет таблицу Excel с результатами этих геообработок. Однако модели необходимо сначала открыть файл Excel, чтобы пользователь мог отредактировать его вручную, а...
Есть ли способ получить список установленных форматов gdal в системе с python? В командной строке это легко сделать с помощью: gdalinfo --formats ... но можно ли это сделать на...
Я хочу сделать простую функцию if / else внутри моей Модели, и я хотел знать, возможно ли это или нет. Я хочу, чтобы это делалось так, если выбран Вход для «Сегментного слоя», чтобы он не выполнял «Сегментацию (означает смещение)» и использует вместо этого «Сегментный слой» в качестве Входных...
Есть ли способ проверить, имеют ли два заданных растровых слоя одинаковое содержимое ? У нас есть проблема с нашим корпоративным общим хранилищем: теперь он настолько большой, что для полного резервного копирования требуется более 3 дней. Предварительное расследование показало, что одним из самых...
Я хочу привязать растр с помощью pythonи GDAL. Мой текущий подход состоит в том, чтобы вызвать gdal_translateи gdalwarpиспользовать os.systemи ужасный список наземных контрольных точек. Я действительно хотел бы сделать это изначально внутри python. Это текущий процесс, который я использую: import...
Безусловно, Доступ к ArcObjects из Python? мои наиболее часто задаваемые вопросы и ответы на GIS Stack Exchange. Несмотря на этот успех, это, вероятно, одна из моих самых слабых областей, когда дело доходит до реального использования. Большая часть этого плохого показа проистекает из моей плохой...
Как я могу получить проекционные координаты, а также фактические значения пикселей в этих координатах из файла GeoTiff, а затем сохранить их в массив массивов? У меня есть файл arsenci020l.tif, и его координаты указаны в метрах. Ниже приведен сокращенный вывод gdalinfo, на котором я работал. ~$...
Интересно, есть ли способ открыть консоль Python непосредственно при запуске программы. Существует возможность установить ярлык для консоли Python, но я не могу найти такой вариант для запуска программы...
У меня есть набор линий, представляющих автобусные линии. Некоторые линии перекрываются и проходят по одним и тем же дорогам. Я могу извлечь узлы. Однако я заинтересован в извлечении только таких переходов: Как я могу это сделать? Я ищу способы с QGIS или Python. Я попробовал метод пересечения из...