Запитання з тегом «rgdal»

У мене є функція, що міститься в базі даних геоданих, яка перевищує 2 ГБ як експортований файл форми. Мені потрібно запустити функцію вилучення в R, щоб віднести багатокутники до даних із растрового файлу. Експорт функції у вигляді таблиці не є рішенням. Як я можу читати класи функцій, що містяться у …

27 r  features  file-geodatabase  feature-class  rgdal 

Я намагаюся прочитати таблицю безпосередньо з бази даних геоданих ESRI у Р. Приклад файлу даних можна завантажити тут . База даних містить точковий клас функцій (Zone9_2014_01_Broadcast) та дві пов'язані таблиці (Zone9_2014_01_Vessel та Zone9_2014_01_Voyage). Ви можете прочитати файл форми у R за readOGRдопомогою rgeosпакету: 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 

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 

Я хочу створити два багатокутники. Один з прямокутних протяжностей мого растра. Я знаю, що я можу отримати масштаб растра, використовуючи r <- raster("band5.tif") e <- extent(r) plot(e)мені це: Як я можу створити ОДИН багатокутник межі растра, як показано нижче? Ось як виглядає мій растр:

13 raster  r  rgdal 

Я використовую gSimplify(пакет rgeos) для спрощення геометрії форми файлів. Функція працює добре, але тепер я не можу записати вихід у новому форматі. Я спробував кілька способів: writeOGR(simplyshape, file, driver="ESRI Shapefile", layer='test') я отримав Obj повинен бути SpatialPointsDataFrame, SpatialLinesDataFrame або SpatialPolygonsDataFrame і з: writePolyShape(simplyshape, file) Я отримав: Помилка: is (x, "SpatialPolygonsDataFrame") …

12 shapefile  r  simplify  rgdal 

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

11 r  rgdal 

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

11 gdal  r  rgdal  geopackage  sf 

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

10 shapefile  r  rgdal 

Я перетворюю вектор у растр у Р. Однак процес був занадто довгим. Чи є можливість поставити скрипт в обробку багатопотокових чи графічних процесорів, щоб зробити це швидше? Мій сценарій для растрованих вектор. r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do raster res(r.raster) = 10 #definindo o tamanho do …

9 raster  r  rasterization  rgdal  parallel