У меня есть список нескольких сотен городов с их широтой / долготой. Учитывая другое местоположение (также в лат / лонг) мне нужно найти ближайший город. Поскольку я не использую ГИС, к настоящему времени очевидным алгоритмом является создание цикла для всех городов, вычисляющего расстояние между...
Я пытаюсь выяснить, как загрузить растр в базу данных PostGIS2.0 (я задавал предыдущие вопросы по этой теме здесь и здесь ). Я пытаюсь использовать raster2pgsql.exeпрограмму, поставляемую с PostGIS2.0. После того, как выяснилось, что командную строку в Windows нужно запускать с правами...
Можно ли создавать представления PostgreSQL над таблицей, содержащей растровые данные? Мне удалось создать представление, но запись raster_columnsне правильная (отсутствует вся растровая информация). Кроме того, я не мог открыть представление SQL с помощью QGIS (скорее всего, из-за неправильного...
У меня небольшие проблемы с Python OGR API. То, что я пытаюсь сделать, это получить все координаты каждой вершины внешнего кольца многоугольника. Это то, что я до сих пор: import osgeo.ogr import glob path = "/home/woo/maps/" out = path + 'output.txt' file = open(out,'w') for filename in...
Под обычаем я не подразумеваю другую структуру файла KML, но в <description></description>разделе у меня есть данные с тегами HTML такого типа: <?xml version="1.0" encoding="UTF-8"?> <kml xmlns="http://www.opengis.net/kml/2.2"> <Document> <Placemark>...
С помощью gis.stackexchange.com я создал лидерные линии для точечных символов. Линии основаны на символике генератора геометрии: Сокращение геометрии линии с помощью диалогового выражения в QGIS? Теперь я хочу реализовать выносные линии для перемещенных надписей. Для меня очень важно, чтобы линии...
У меня есть Contour Featureclass, с интервалом 0,2 метра. Я хочу выбрать только те контуры, которые с интервалом 1 м? Я попытался использовать следующий синтаксис в диалоговом окне «Выбор по атрибуту» (в ArcMap), но он выбирает все контуры: Mod("ELEVATION", 1)=0 Какой запрос я должен...
Я пытаюсь получить число путем подсчета количества точек в шейп-файле. И это работает, за исключением того, что у меня тогда возникают проблемы с использованием этого числа где-то еще. В конце концов, я буду использовать этот счет в некоторой математике (полевом калькуляторе), но во время отладки я...
Этот пост является сообществом Wiki . Отредактируйте существующие ответы, чтобы улучшить этот пост. В настоящее время он не принимает новые ответы. Я занимаюсь разработкой некоторого кода на Python и хочу использовать привязки Python GDAL / OGR. Какой минимум мне нужно установить, чтобы...
Я хочу преобразовать GeoPDF с векторной сущностью в векторный формат и сохранить реальные координаты этих векторов. GeoPDF был экспортирован из ArcMap и содержит слои. Я попробовал GDAL, но он растеризовал PDF перед экспортом. Что не хорошо для меня. Также я не хочу решения, которое конвертировало...
У меня есть голая земная матрица, сделанная из данных лидара (A). Он хорошо работает и прошел все тесты на точность. Используется в моделях штормовых нагонов. Проблема, которую мы имеем, находится на подземных переходах. Мне нужно позволить воде течь через подземные переходы, когда одна дорога...
Я пытаюсь выполнить объединение на общем поле после объединения двух смежных шейп-файлов. Шейп-файлы заканчиваются как минимум одним узким промежутком между ними. Когда я пытаюсь объединиться, я получаю следующую ошибку: Ошибка в createPolygonsComment (p): rgeos_PolyCreateComment: потерянная дыра,...
Я надеюсь, что кто-то уже понял это. Я только что установил Geoserver 2.9 на ванильный дистрибутив Ubuntu 16.04. Метод Geoserver 2.8 для включения CORS с классом shanbe.hezoun больше не работает с Jetty 9.2.13. Есть упоминания о том, что поддержка CORS уже включена в Jetty 9.2.13 в файле...
Я хотел бы обобщить SpatialPolygonsDataFrame в R, чтобы сделать построение графиков более эффективным. Полигоны, которые я имею, приводят к огромным pdf-графикам, которые медленно загружаются и трудно интегрируются в документ. Я попытался использовать gSimplifyфункцию из rgeosпакета, но, к...
У меня есть растр (на самом деле USGS DEM), и мне нужно разделить его на более мелкие фрагменты, как показано на рисунке ниже. Это было достигнуто в ArcGIS 10.0 с использованием инструмента Split Raster. Я хотел бы метод FOSS, чтобы сделать это. Я посмотрел на GDAL, думая, что он обязательно это...
Я рассчитываю рассчитать фокальную статистику для каждой ячейки растра в окрестности указанного критерия. Справочная информация - у меня есть три бинарных растра, каждый из которых представляет один интересующий тип растительности. Я хотел бы рассчитать процент покрытия каждого типа растительности...
В настоящее время я использую следующий сценарий для добавления некоторых данных атрибутов из таблицы во множество отдельных шейп-файлов: library(rgdal) specieslist <- read.csv("SpeciesList1.txt", header=F) attdata <- read.table("TestAtt.csv", sep = ",", header=T) for (n in...
Глядя на эту информационную страницу для набора данных, размещенного на esri.com: http://fema-services2.esri.com/arcgis/rest/services/2012_Sandy/ImageCat_NLT/MapServer/layers Попытка выяснить, как я могу получить доступ к необработанным данным (либо необработанные координаты широты / долготы, либо...
У меня есть шейп-файл с полигонами. И у меня есть глобальный растровый файл. Я хочу наложить полигоны шейп-файла на растровую сетку и рассчитать среднее значение растра для каждого полигона. Как я могу сделать это с помощью GDAL, записав результаты в...
Я использую ArcGIS 10.2 для рабочего стола, и я просматривал этот сайт, и пока не могу найти свой ответ. Очень плохо знаком с VBA и Python, но много лет с ArcGIS. Я знаю, что могу сделать это медленно с помощью Select By Attributes, но это отнимает много времени. Я пытаюсь выполнить...