Вопросы по теме 'sp'
Неправильные значения, возвращаемые функцией sp::over?
Я извлекаю данные о высоте маршрута из цифровой модели высоты, используя
my.elev <- over(new.points, mygrid)
new.points — это объект SpatialPoints с координатами (долгота/широта) около 7000 точек, преобразованными в CRS mygrid.
mygrid...
117 просмотров
schedule
05.12.2023
R удалить повторяющиеся пространственные точки в соответствии с атрибутом
В R у меня есть дублированная точка SpatialPointsDataFrame (координаты и атрибуты), я хотел бы удалить все точки с теми же данными ...
Я нашел в пакете sp функцию remove.duplicates() , но, похоже, она удаляется только на месте ... Есть ли другой...
5610 просмотров
schedule
09.12.2021
Как создать объект SpatialLine
Я использую пакет sp для создания объекта SpatialLines и сохранения его в списке объектов allLines . Позже мне нужно будет сравнить SpatialLines друг с другом, но это выходит за рамки текущего вопроса.
Пока мне нужно построить только...
6060 просмотров
schedule
31.05.2022
Использование простого цикла for для пространственных данных
Мне очень жаль, что это будет вопрос для цикла 101. Я изо всех сил пытаюсь написать простой цикл для создания таблицы расстояний между городами на основе данных долготы и широты.
locations <-read.csv("distances.csv")
Location возвращает...
444 просмотров
schedule
13.04.2023
Совокупный кадр данных для более грубого пространственного разрешения
У меня есть набор данных с разрешением сетки 0,25 * 0,25 градуса. У меня есть другой с разрешением 1 * 1 градус, и я хочу сделать их сопоставимыми (оба изменены на разрешение 1 * 1). Область составляет от 28,5 до 36,5 долготы и от -4,5 до 4,5 широты...
1157 просмотров
schedule
11.07.2023
Использование grid.newpage для создания карты ggplot с картой-врезкой
Следуя очень полезному примеру здесь . Я собрал код для карты, созданной в ggplot, со вставкой.
# Clear WD
rm(list=ls())
library(sp)
library(ggplot2)
library(raster)
library(rgeos)
library(gridExtra)
# Download the Level 0 Map for Inset
ind0...
518 просмотров
schedule
10.06.2023
В чем разница между z~1 и z~LON+LAT при построении вариограммы в R?
library(sp)
data(meuse)
coordinates(meuse)<-~x+y
v<-variogram(log(zinc)~1,meuse)
v
np dist gamma dir.hor dir.ver id
1 57 79.29244 0.1234479 0 0 var1
2 299 163.97367 0.2162185 0 0 var1...
128 просмотров
schedule
07.12.2022
Как создать карту в R из пространственного многоугольника, заданного как тип данных списка
У меня есть список (импортированный из файла Stata dta), содержащий записи, определяющие пространственные полигоны для карты, чтобы дать представление:
> typeof(aux3)
[1] "list"
> mode(aux3)
[1] "list"
> head(aux3)
_ID _X...
588 просмотров
schedule
27.07.2022
finalPolygonCRS + пакет функций over () в R
Я создаю SpatialPolygon таким образом
#make spatial points
#assign Original CRS WGS84 EPSG:4326 (DATUM --> on sphere)
cornersEPSG4326 <- SpatialPoints(coords=cbind(x,y), proj4string = CRS("+proj=longlat +ellps=WGS84 +datum=WGS84...
1016 просмотров
schedule
22.09.2022
Недопустимое предупреждение `.internal.selfref`, столбец не обновляется при работе с `SpatialPolygonsDataFrame`
Я пытаюсь выполнить некоторый геопространственный анализ в R , который будет включать добавление атрибутов в SpatialPolygonsDataFrame для раскрашивания и т. д. во время построения графика.
Для организации я хотел бы добавить эти атрибуты в свои...
520 просмотров
schedule
24.06.2022
Искажение многоугольника, взвешенное по населению (картограммы)
Я пытаюсь создать карту в R, которая передает как форму базовой геометрии (т. е. физические границы), и относительную важность объекта с точки зрения связанного значения.
Для конкретики я хотел бы сосредоточиться на воспроизведении (версии)...
1383 просмотров
schedule
03.04.2022
точки наложения на многоугольники в R
Мне нужно наложить 40 000 точек на 100 000 полигонов. Полигоны были созданы с использованием:
polygon <- gBuffer(pc, width=500, byid=TRUE)
Оверлей был создан с использованием:
test <- over(pts, polygon)
В результате у меня...
1932 просмотров
schedule
15.09.2022
Объединяйте пространственные точки и используйте интерполяцию, чтобы заполнить пробелы
У меня есть набор фреймов данных, и я хотел бы использовать их для создания одной растровой сетки. Вот два примера файлов: dat1 dat2
Я использую пакет sp для построения пространственных точек:
library(raster)
library(sp)...
545 просмотров
schedule
13.12.2022
Ошибка: не удалось найти оверлей функций после переустановки
Я переустановил свой компьютер, а также переустановил R, но почему-то получаю сообщение об ошибке при запуске моей программы R:
Ошибка: не удалось найти функцию "оверлей"
что было хорошо на предыдущей установке.
Новая версия R — 3.2.4.,...
1381 просмотров
schedule
04.12.2022
sp R пакет SpatialPolygons-Class график прозрачный
Я хочу построить объекты класса SpatialPolygons с прозрачным цветом заливки. Немного изменен пример со страниц справки SpatialPolygons-Class:
# simple example, from vignette("sp"):
Sr1 = Polygon(cbind(c(2,4,4,1,2),c(2,3,5,4,2)))
Sr2 =...
221 просмотров
schedule
03.08.2022
Как эффективно рассчитать расстояние между парой координат с помощью data.table: =
Я хочу найти наиболее эффективный (самый быстрый) метод вычисления расстояний между парами координат широта и долгота.
Не очень эффективное решение было представлено (здесь) с использованием sapply и spDistsN1{sp} . Я считаю, что это можно...
7201 просмотров
schedule
01.02.2023
Как изменить класс с фрейма данных на пространственный полигон?
Я нашел то же самое здесь Как преобразовать фрейм данных в пространственные координаты . Но в моем случае я получил очень большие данные.
exchange longitude latitude
AB 103.3281386 1.594218196
AB 103.3285929 1.593990735
AB 103.3312494...
4169 просмотров
schedule
01.12.2022
Разделить полилинию полилинией
Я хочу разделить строку на строку. (Я хочу сделать как алгоритм QGIS "разделить строки линиями")
функция "gsection" в пакете stplanr делает это.
library(stplanr)
data(routes_fast)
result <- gsection(routes_fast)
class(result)
Но...
441 просмотров
schedule
14.03.2023
Задача построения (простых) полигонов с буклетом в r
Я хотел бы использовать брошюру, чтобы помочь мне построить многоугольники тепловой карты по Великобритании, но листовка ()%>% addPolygons (...) дает мне странный результат даже для очень простого случая.
Я использую следующий код для построения...
689 просмотров
schedule
06.06.2022
Количество сегментов в SpatialLinesDataFrame
После импорта кадра данных пространственных линий с помощью readOGR я хотел бы определить количество сегментов линии (узлов) в каждой строке. Я не могу найти практический способ просто экспортировать длины в виде вектора.
routes@lines приведет...
198 просмотров
schedule
03.10.2022