Учитывая ограничивающий прямоугольник, как я могу подсчитать количество определенных значений (скажем, меня интересует количество value == 1) в растре в 1) arcpy или 2) растровом пакете в...
Я использую gdalwarp для манипулирования плитками SRTM рядом с датой (т.е. 180 °, иначе антимеридиан). Плитки SRTM имеют очень небольшое (1/2 пикселя) перекрытие с меридианом. Вы можете увидеть это с помощью gdalinfo: gdalinfo S16W180.hgt Driver: SRTMHGT/SRTMHGT File Format Files: S16W180.hgt Size...
У меня есть данные округа, в которых FIPs (на самом деле FIPS штата и округа FIPS вместе взятые). У меня есть данные, импортированные как CSV с теми же данными FIPS под названием STCounty. Однако шейп-файл FIPS указан как Real, а CSV - как Integer, и функция Join, похоже, не работает. Откидные...
Я использую OpenLayers с mapserver для отображения большого количества растровых изображений и запрашиваемых слоев wms. Все работает нормально, слои можно включать / выключать в переключателе слоев OpenLayers ... Но теперь я хочу включить / выключить слои с помощью внешних флажков и не знаю, как...
У меня есть пара точек координат (начальная и конечная точки), которые я должен преобразовать в линии. До сих пор я использовал добавление обеих координат в a pippo.Point(), a, pippo.CalculateGeometry()чтобы определить геометрию каждого пионта и pippo.append(defined geometry)идентифицировать пару...
Я пытаюсь пересечь два SpatialPolygonsDataFrames и получить SpatialPolygonsDataFrame в результате. К сожалению, используя gIntersectionфункцию from rgeos(которая работает впечатляюще быстро для пересечения полигонов), я не могу извлечь соответствующие кадры данных. Рассмотрим следующий пример: >...
У меня есть линейные функции, которые должны использоваться для этикеток. Эти линии правильно ориентированы (повернуты) и расположены таким образом, чтобы текст, помещенный на них, давал хороший картографический вывод. Сравните два изображения ниже и их ярлыки. В частности слово «Waverton» и его...
Я использую osm2pgsql для заполнения базы данных OSM в Северной Америке из planet.osm. Мне нравится, насколько проста полученная схема PostGIS, но меня немного смущает ее содержание. У меня есть две таблицы planet_osm_roadsи planet_osm_line. Я предположил, что planet_osm_roads содержит все данные о...
Я ищу способ реализовать видимость (он же видимость). Я нашел предыдущий вопрос GIS SE, в котором предлагалось несколько приложений ГИС (не то, что мне нужно - я пытаюсь встроить вычисления), а также SAGA. SAGA больше похож на то, что я ищу (классы C ++, Java, C #), но я не могу понять достаточно...
Я нахожусь в процессе ggplotсоставления картограммы населения в административных районах Уэльса. Я скачал данные Boundary Line из Ordnance Survey и извлек что-то, что кажется правильным шейп-файлом (community_ward_region.shp). Используя R, я дошел до чтения в шейп-файле. require(maptools) shape...
Я только начал работать с пространственными базами данных и хочу написать запрос SQL (PostGIS) для автоматического обобщения необработанных GPS-треков (с фиксированной частотой отслеживания). Первое, на чем я работаю, - это запрос, который идентифицирует точки остановки в форме запроса, такого как...
Я много искал, но не нашел ответа: при редактировании строк с помощью QGIS мне часто нужно разделить линейную строку точно по одной из его существующих внутренних вершин, но я не могу найти способ; Я пробовал инструмент «Разделить функции», думая, что он может привязываться к вершинам, но я ошибся,...
Рассмотрим простую ситуацию, когда три ребра соединяются в узле: Я хотел бы написать краткое и четкое описание отношений между A и B таким образом, чтобы отличать их от отношений между A и C. Что-то вроде «при прохождении узла по часовой стрелке A соседствует? к Б, но А не соседствует? в C. » Но...
Я хотел бы узнать, как выполнить неконтролируемую классификацию сцены Landsat, используя i.cluster> i.maxlikв GRASS, используя заостренные изображения с разрешением 15 м (пример, приведенный в его вики, использует обычные разрешения с 30 м). i.pansharpenСначала я попытался с помощью создания...
Я хочу перепроектировать файл растровых изображений из проекции A (ETRS89 / UTM Zone 32) в проекцию B (DHDN / Gauss-Kruger zone 2). Две проекции имеют разные данные, поэтому я должен использовать преобразование данных (NTv2 BETa2007 в моем случае). Я хочу создать новый растровый файл с проекцией B,...
Я хочу импортировать шейп-файл в PostGIS. В нем есть одна база данных, postgresно таблица еще не добавлена. Я использую эту команду: shp2pgsql -I -s 4326 districts.shp districts | psql -d postgres -U postgres и это дает мне ошибку: Невозможно преобразовать значение данных в UTF-8 (iconv сообщает...
У нас есть некоторое программное обеспечение как часть более крупного проекта, написанного с использованием ArcObjects через .NET. Это программное обеспечение переносит и объединяет данные из нескольких файловых баз геоданных в одной схеме в одну файловую базу геоданных в аналогичной, но другой...
Я планирую кластеризовать GeoServer. В Интернете я могу найти несколько результатов, говорящих о том, что часто кластеризация настраивается путем настройки нескольких узлов GeoServer, которые все используют один и тот же каталог данных. В то же время есть некоторые результаты (...
Позволяет ли ArcGIS Desktop или какое-либо из его собственных расширений создавать ЦМР из стереопар (приобретенных через спутник GeoEye 1 )? Если ответом будет «нет», «коричневые очки» перейдут к ответу, предлагающему жизнеспособную альтернативу с использованием...
У меня есть два набора полигонов в двух таблицах. Наборы перекрывают друг друга. Для каждого многоугольника в наборе A я хотел бы получить идентификатор многоугольника в наборе B, который он перекрывает больше всего. Я использую PostgreSQL с расширением PostGIS. Я достаточно разбираюсь в SQL, чтобы...