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

27
Чтение класса объектов в файловой базе геоданных с использованием R?

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

21
Чтение таблицы из файловой базы геоданных ESRI (.gdb) с использованием R

Я пытаюсь прочитать таблицу непосредственно из файловой базы геоданных ESRI в R. Пример файла данных можно скачать здесь . База данных содержит класс точечных объектов (Zone9_2014_01_Broadcast) и две связанные таблицы (Zone9_2014_01_Vessel и Zone9_2014_01_Voyage). Вы можете прочитать шейп readOGR-...

19
Проверка, попадают ли точки в полигональный шейп-файл

У Zillow есть набор шейп-файлов для разных районов крупных городов США. Я хотел проверить, присутствовали ли определенные здания в определенных районах, используя R: library(rgeos) library(sp) library(rgdal) df <- data.frame(Latitude =c(47.591351, 47.62212,47.595152), Longitude =...

19
Извлечение областей пересечения в R

У меня есть два полигона. Одно содержит поля (X, Y, Z), а другое содержит типы почвы (A, B, C, D). Я хочу знать, какая область каждого поля содержит, какой тип почвы. Я попробовал следующее: library(rgdal) library(rgeos) Field<-readOGR("./","Field") Soil<-readOGR("./","Soil")...

13
Создать полигоны экстента и границы данного растра в R

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

12
написать OGR с пространственным полигоном, упрощенным с помощью gSimplify

Я использую gSimplify(пакет rgeos) для упрощения геометрии шейп-файла. Функция работает хорошо, но теперь я не могу записать вывод в новый шейп-файл. Я попробовал несколько способов: writeOGR(simplyshape, file, driver="ESRI Shapefile", layer='test') я получил объект должен быть...

11
R - Создать ограничивающий прямоугольник, преобразовать в класс Polygon и Plot

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

11
Запись нескольких слоев в GeoPackage с использованием writeOGR () в R?

Я пытаюсь записать несколько слоев в один и тот же GeoPackage в R, но получаю ошибку Creation of output file failed. Я попытался найти документацию для чтения и записи в файлы .gpkg с помощью RGDAL, например, чтобы выяснить, writeOGR()поддерживает ли на самом деле несколько слоев, с небольшим...

10
Читать только часть шейп-файла ESRI в R

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

9
Обработка вектора в растре быстрее с помощью R

Я конвертирую вектор в растр в R. Однако процесс был слишком долгим. Есть ли возможность поместить скрипт в многопоточную или графическую обработку, чтобы сделать это быстрее? Мой скрипт в растеризованный вектор. r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do raster...