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

2
Чтение класса объектов в файловой базе геоданных с использованием R?
У меня есть функция, содержащаяся в базе геоданных размером более 2 ГБ в качестве экспортируемого шейп-файла. Мне нужно запустить функцию извлечения в R, чтобы приписать полигоны с данными из растрового файла. Экспорт функции в виде таблицы не является решением. Как я могу прочитать классы объектов, содержащиеся в файловой базе геоданных …

5
Чтение таблицы из файловой базы геоданных ESRI (.gdb) с использованием R
Я пытаюсь прочитать таблицу непосредственно из файловой базы геоданных 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 …

3
Проверка, попадают ли точки в полигональный шейп-файл
У 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 

2
Извлечение областей пересечения в R
У меня есть два полигона. Одно содержит поля (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 …

1
Создать полигоны экстента и границы данного растра в R
Я хочу создать два полигона. Один из прямоугольных размеров моего растра. Я знаю, что могу получить степень использования растра, r <- raster("band5.tif") e <- extent(r) plot(e)дает мне это: Как я могу создать один полигон границы растра, как показано ниже? Вот как выглядит мой растр:
13 raster  r  rgdal 

2
написать OGR с пространственным полигоном, упрощенным с помощью gSimplify
Я использую gSimplify(пакет rgeos) для упрощения геометрии шейп-файла. Функция работает хорошо, но теперь я не могу записать вывод в новый шейп-файл. Я попробовал несколько способов: writeOGR(simplyshape, file, driver="ESRI Shapefile", layer='test') я получил объект должен быть SpatialPointsDataFrame, SpatialLinesDataFrame или SpatialPolygonsDataFrame и с: writePolyShape(simplyshape, file) Я получил: Ошибка: is (x, "SpatialPolygonsDataFrame") не …
12 shapefile  r  simplify  rgdal 

2
R - Создать ограничивающий прямоугольник, преобразовать в класс Polygon и Plot
У меня есть NE lat lngи SW lat lng. Моя цель проста, я использую RStudio, и я хочу создать ограничивающий прямоугольник из двух вышеупомянутых латов и длин, а затем сгенерировать Polygon, и я хочу добавить случайные точки в квадратный многоугольник. У меня есть NELat / Lng и SW-Lat / Lng …
11 r  rgdal 

1
Запись нескольких слоев в GeoPackage с использованием writeOGR () в R?
Я пытаюсь записать несколько слоев в один и тот же GeoPackage в R, но получаю ошибку Creation of output file failed. Я попытался найти документацию для чтения и записи в файлы .gpkg с помощью RGDAL, например, чтобы выяснить, writeOGR()поддерживает ли на самом деле несколько слоев, с небольшим успехом. Это вообще …
11 gdal  r  rgdal  geopackage  sf 

3
Читать только часть шейп-файла ESRI в R
У меня есть очень большой шейп-файл ESRI (> 8000 полигонов, файл .shp> 32 МБ), который я хочу прочитать в R в Windows 7. Я знаком с rgdal и readOGR (). Однако при таком размере шейп-файла считывание всего шейп-файла в память является очень медленным процессом. Кроме того, меня интересуют только несколько …
10 shapefile  r  rgdal 

1
Обработка вектора в растре быстрее с помощью R
Я конвертирую вектор в растр в R. Однако процесс был слишком долгим. Есть ли возможность поместить скрипт в многопоточную или графическую обработку, чтобы сделать это быстрее? Мой скрипт в растеризованный вектор. r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do raster res(r.raster) = 10 #definindo o tamanho do pixel …
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.