Вопросы с тегом «r»

Я использовал интеграцию R с QGIS и нашел ее очень полезным инструментом. Взглянув на расширение набора инструментов R для 9 и расширение на 10, кажется, что они оба предлагают только точечную кластеризацию, чтобы продемонстрировать подтверждение концепции. Могут ли другие методы анализа, доступные в R, быть реализованы в ArcGIS Desktop с …

22 arcgis-desktop  r  spatstat 

Я создал 2 полигона в QGIS. Используя их в R, полигоны автоматически становятся SpatialPolygonsDataFrame (SPDF). Я хотел бы объединить их в один SPDF (как это легко сделать в ArcGis, используя инструмент Merge ). Я уверен, что должен быть простой способ, как это сделать в R, но я не могу найти …

22 r  spatial-analyst  merge  maptools 

Я пытаюсь выполнить пространственное соединение между данными точек и данными многоугольников. У меня есть данные, которые указывают пространственные координаты события в моем CSV-файле A, и у меня есть другой файл, шейп-файл B, который содержит границы области в виде полигонов. head(A) month longitude latitude lsoa_code crime_type 1 2014-09 -1.550626 53.59740 E01007359 …

21 r  spatial-join  point-in-polygon 

Я пытаюсь прочитать таблицу непосредственно из файловой базы геоданных ESRI в R. Пример файла данных можно скачать здесь . База данных содержит класс точечных объектов (Zone9_2014_01_Broadcast) и две связанные таблицы (Zone9_2014_01_Vessel и Zone9_2014_01_Voyage). Вы можете прочитать шейп readOGR- rgeosфайл в R используя пакет: library(rgeos) library(downloader) download("https://coast.noaa.gov/htdata/CMSP/AISDataHandler/2014/01/Zone9_2014_01.zip", dest="Zone9_2014_01.zip", mode="wb") unzip("Zone9_2014_01.zip", exdir …

21 r  file-geodatabase  table  rgdal 

Есть ли функция обрезки sf карты, похожая на maptools::pruneMap(lines, xlim= c(4, 10), ylim= c(10, 15))используемую для SpatialPolygon или SpatialLine? Я рассматриваю, st_intersection()но может быть правильный путь.

20 r  clip  sf 

Учитывая большую (~ 1 миллион) выборку неравномерно распределенных точек - возможно ли создать нерегулярную сетку (по размеру, но также может быть неправильной формы, если это возможно?), Которая будет содержать указанное минимальное количество из n точек? Для меня менее важно, если генерируемые «ячейки» такой сетки содержат ровно n точек или хотя …

20 r  clustering  quadtree 

Как я могу преобразовать точки lon-lat в простые объекты (sfg), а затем поместить их в простой набор объектов (sfc)? Вот MWE, который не работает, но самый близкий, который я получил. library(data.table) library(sf) # The DT data.table is the data I have (but 10,000s of rows, each row is a point) …

20 r  sf 

Я строю некоторые геолокации (lon, lat)твитов, собранных с помощью Stream API с 7mi. радиус. На рисунке ниже показаны две карты с одинаковыми данными, но с разной альфа / прозрачностью для каждого наблюдения. Глядя на карту слева, я заметил, что геолокации следуют в основном за некой воображаемой «сеткой» с конкретными узлами. …

19 r  geolocation  twitter 

coast<-readShapeSpatial("coastline.shp") landc<-readShapeSpatial("landcover.shp") ro<-readShapeSpatial("roads.shp") bc<-gBuffer(ro,width=100) landc$ratings=1 landc$ratings[landc$LANDUSE_ID==4]=0 Выше я беру любую категорию, которая имеет 4, и в новом столбце ставлю ее как 0. На данный момент, я хочу, чтобы столбец также был назван ratingsдля bc, где он будет принимать 0, если он внутри буфера, и 1, если он снаружи. Проблема в …

19 r  cartography 

У меня есть набор растровых изображений, которые представляют определенный месяц по годам, и я хочу создать временную шкалу средних значений области, используя шейп-файл. Как извлечь значения из растров и импортировать их в R так, чтобы я мог их использовать?

19 r  spatial-statistics 

Морана I , мера пространственной автокорреляции, не является особенно надежной статистикой (она может быть чувствительной к искаженному распределению атрибутов пространственных данных). Какие существуют более надежные методы измерения пространственной автокорреляции? Меня особенно интересуют решения, которые легко доступны / могут быть реализованы на языке сценариев, например R. Если решения применимы к уникальным …

19 algorithm  r  spatial-statistics 

У Zillow есть набор шейп-файлов для разных районов крупных городов США. Я хотел проверить, присутствовали ли определенные здания в определенных районах, используя R: library(rgeos) library(sp) library(rgdal) df <- data.frame(Latitude =c(47.591351, 47.62212,47.595152), Longitude = c(-122.332271,-122.353985,-122.331639), names = c("Safeco Field", "Key Arena", "Century Link")) coordinates(df) <- ~ Latitude + Longitude wa.map <- …

19 r  rgdal  zillow 

У меня есть два полигона. Одно содержит поля (X, Y, Z), а другое содержит типы почвы (A, B, C, D). Я хочу знать, какая область каждого поля содержит, какой тип почвы. Я попробовал следующее: library(rgdal) library(rgeos) Field<-readOGR("./","Field") Soil<-readOGR("./","Soil") Results<-gIntersects(Soil,Field,byid=TRUE) rownames(Results)<-Field@data$FieldName colnames(Results)<-Soil@data$SoilType > Results A B C D Z TRUE FALSE …

19 shapefile  polygon  r  intersection  rgdal 

Я пытаюсь выполнить объединение на общем поле после объединения двух смежных шейп-файлов. Шейп-файлы заканчиваются как минимум одним узким промежутком между ними. Когда я пытаюсь объединиться, я получаю следующую ошибку: Ошибка в createPolygonsComment (p): rgeos_PolyCreateComment: потерянная дыра, не удается найти содержащий полигон для дыры с индексом 17 Я загрузил воспроизводимый пример …

18 polygon  r  merge 

Моя цель - изменить существующий шейп-файл путем объединения определенных полигонов. После импорта шейп-файла и использования команды UnionSpatialPolygons я получаю желаемый контур многоугольника. Однако теперь это объект SpatialPolygons, а не SpatialPolygonsDataFrame, поэтому я не могу экспортировать его в шейп-файл с помощью writeOGR. Как я могу обойти эту проблему?

18 shapefile  cartography  r  ogr