Геоинформационные системы

На картографической службе армии США (AMS C762) 1: 50 000 карт, сделанных из Исландии 1948–1951 гг., Вверху есть шкала с обозначениями градусов, я полагаю, углов. На изображении ниже показано 25–29 °. Оно увеличивается влево, что предполагает, что оно может относиться к градусам долготы на запад, но цифры не верны; карта …

14 maps  topography  angles  iceland 

Как мне упростить sfмногоугольник без введения пробелов и осколков? С шейп-файлом, например, я бы использовал rmapshaper::ms_simplify(): library("pryr") library("rgdal") library("rmapshaper") download.file("https://borders.ukdataservice.ac.uk/ukborders/easy_download/prebuilt/shape/England_gor_2011.zip", destfile = "regions.zip") unzip("regions.zip") regions <- readOGR(".", "england_gor_2011") object_size(regions) # ~13MB regions <- ms_simplify(regions) object_size(regions) # < 1MB Я пробовал, sf::st_cast()который, из man-страниц, заявляет: Приведите геометрию к другому типу: либо …

14 r  simplify  sf 

То, что я хочу сделать, это поместить 100-миллиметровый буфер на многоугольник, но вместо того, чтобы буфер находился на фиксированной ширине «по прямой линии», я хочу, чтобы он следовал склонам моей лидаровой ЦМР. У меня есть ArcMap (с 3D Analyst). К сожалению, нет пространственного или сетевого аналитического расширения с ArcMap. Я …

14 arcgis-desktop  distance  elevation  3d-analyst  cost-path 

USDA (через Национальную карту) предоставляет 4-полосный NAIP ArcGIS «ImageServer» ( сервисы ) для большинства штатов США. Я могу успешно подключиться к сервису WMS в QGIS, как описано здесь , который загружает изображения в QGIS с использованием диапазонов RGB (1,2,3) по умолчанию. Тем не менее, я хотел бы изменить изображения для …

14 qgis  configuration  image-service 

Я пытаюсь разделить около 4000 полигонов в их средней точке, перпендикулярно их самой длинной оси (то есть поперек ширины в средней точке), как показано на диаграмме ниже. В идеале я хотел бы сделать это автоматически и избегать разбиения каждого полигона вручную. Я извлек среднюю точку полигонов, преобразовав самые длинные линии, …

14 arcpy  polygon  splitting 

Я изо всех сил пытаюсь четко понять различие между классификацией на основе пикселей и объектами в области дистанционного зондирования, и я надеюсь, что кто-то из этого сообщества сможет обеспечить понимание. Основываясь на информации, которую я имею до сих пор, мое текущее понимание имеет следующие черты: Классификация на основе пикселей : …

14 remote-sensing  classification  machine-learning 

В пакете HistData для R ( https://r-forge.r-project.org/R/?group_id=574 ) у меня есть наборы данных, относящиеся к карте Джона Сноу о вспышке холеры в Лондоне в 1854 году. Я считаю, что они являются авторитетными, будучи тщательно оцифрованы под наблюдением Уолтера Тоблера. Некоторые подробности об этих наборах данных описаны Джоном Маккензи по адресу …

14 r  coordinate-system  scale 

Не могли бы вы объяснить, что означает «буфер» в диалоге Вороного, то есть что означает число?

14 qgis  buffer  voronoi-thiessen 

На мой взгляд, если два объекта, подлежащих маркировке, расположены слишком близко друг к другу, то по сути есть три варианта. По умолчанию QGIS пропускает одну из меток. Если я включу свойство «показать все метки», либо для рассматриваемого слоя, либо для глобальных настроек надписей, то оно будет печатать перекрывающиеся метки. Если …

14 qgis  labeling 

В QGIS 2.14.5 LTR (и недавно в QGIS 3.0) я создал красивый дизайн карты с градуированным визуализатором и маркерами векторного поля для каждого диапазона значений (длина стрелки прямо пропорциональна скорости потока, которая колеблется в пределах 0,001 - ~ 5 м / с). Это сработало великолепно, но мне интересно, пропустил ли …

14 qgis  legend  markers  qgis-3  vector-field 

Я пытаюсь привязать две линии друг к другу, используя Shapely / Geopandas, но результат привязки очень странный. Я старался : import geopandas as gpd from shapely.geometry import * from shapely.ops import snap lines1 = gpd.GeoDataFrame.from_file('lines1.shp') lines1 = lines1.to_crs({'init': 'epsg:2227'}) lines2 = gpd.GeoDataFrame.from_file('lines2.shp') lines2 = lines2.to_crs({'init': 'epsg:2227'}) res = lines1 lines2_union …

14 shapely  snapping  geopandas 

Чтобы получить значение одного пикселя в точке растра с помощью растерио, приведите пример здесь: https://github.com/mapbox/rasterio/pull/275 Однако существует ли прямой API в растерио (а не в cli), который можно использовать для извлечения значения в одной точке растра? -- РЕДАКТИРОВАТЬ with rasterio.drivers(): # Read raster bands directly to Numpy arrays. # with …

14 python  numpy  rasterio 

Есть ли способ обновить таблицу атрибутов слоя после редактирования таблицы (например, добавить столбец) в DB-Manager?

14 qgis  attribute-table  database 

Я новичок в области дистанционного зондирования. Иногда я путался, когда читал о предварительной обработке изображений. Может ли кто-нибудь рассказать мне о разнице и примерах методов радиометрической коррекции и радиометрической калибровки? Является ли преобразование значений DN в значения отражательной способности Top Of Atmospheric (TOA) из 2 наборов изображений (например, Landsat 7) …

14 remote-sensing 

У меня есть приложение, которое должно выполнить HTTP-запрос (который будет обрабатывать некоторые преобразования карты), прежде чем получать плитки WMS с сервера. Есть ли способ приостановить автоматическое обновление TileWMS после перетаскивания карты? Установка невидимого или удаление и добавление слоя будет работать, но я не хочу, чтобы слой исчезал при выполнении HTTP-запроса. …

14 wms  openlayers