Pertanyaan yang diberi tag «rgdal»

Saya memiliki fitur yang terkandung dalam geodatabase yang lebih besar dari 2GB sebagai shapefile yang diekspor. Saya perlu menjalankan fungsi ekstrak dalam R untuk atribut poligon dengan data dari file raster. Mengekspor fitur sebagai tabel bukanlah solusi. Bagaimana saya bisa membaca kelas fitur yang terkandung dalam geodatabase file Esri?

27 r  features  file-geodatabase  feature-class  rgdal 

Saya mencoba membaca tabel langsung dari file ESRI geodatabase ke dalam R. Contoh file data dapat diunduh di sini . Basis data berisi kelas fitur titik (Zone9_2014_01_Broadcast) dan dua tabel tertaut (Zone9_2014_01_Vessel dan Zone9_2014_01_Voyage). Anda dapat membaca shapefile di R menggunakan readOGRdari rgeospaket: 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 memiliki seperangkat shapefile untuk berbagai lingkungan di kota-kota besar AS. Saya ingin memeriksa apakah ada bangunan tertentu di lingkungan tertentu menggunakan 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 

Saya punya dua poligon. Satu berisi bidang (X, Y, Z) dan yang lainnya berisi jenis tanah (A, B, C, D). Saya ingin tahu area mana dari setiap bidang yang berisi jenis tanah apa. Saya mencoba yang berikut ini: 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 …

19 shapefile  polygon  r  intersection  rgdal 

Saya ingin membuat dua poligon. Salah satunya adalah bentuk persegi panjang raster saya. Saya tahu saya bisa mendapatkan tingkat raster menggunakan r <- raster("band5.tif") e <- extent(r) plot(e)memberi saya ini: Bagaimana saya bisa membuat SATU poligon dari batas raster seperti yang ditunjukkan di bawah ini? Seperti inilah rupa raster saya:

13 raster  r  rgdal 

Saya menggunakan gSimplify(paket rgeos) untuk menyederhanakan geometri sebuah shapefile. Funcion berfungsi dengan baik, tapi sekarang saya tidak bisa menulis output di shapefile baru. Saya mencoba beberapa cara: writeOGR(simplyshape, file, driver="ESRI Shapefile", layer='test') saya mendapat obj harus berupa SpatialPointsDataFrame, SpatialLinesDataFrame, atau SpatialPolygonsDataFrame dan dengan: writePolyShape(simplyshape, file) Saya mendapat: Kesalahan: is (x, …

12 shapefile  r  simplify  rgdal 

Saya punya NE lat lngdan SW lat lng. Tujuan saya sederhana, saya menggunakan RStudio dan saya ingin membuat kotak pembatas dari dua lat dan long di atas dan kemudian menghasilkan Polygon dan saya ingin menambahkan titik acak ke dalam poligon persegi. Saya memiliki NELat / Lng dan SW-Lat / Lng …

11 r  rgdal 

Saya mencoba untuk menulis beberapa layer ke GeoPackage yang sama di R, tapi saya mendapatkan kesalahan Creation of output file failed,. Saya telah mencoba mencari dokumentasi membaca dan menulis ke file .gpkg dengan RGDAL, misalnya untuk mengetahui apakah writeOGR()sebenarnya mendukung banyak lapisan, dengan sedikit keberhasilan. Apakah ini bahkan mungkin, jika …

11 gdal  r  rgdal  geopackage  sf 

Saya memiliki shapefile ESRI yang sangat besar (> 8000 poligon, file .shp> 32MB), yang ingin saya baca di R pada Windows 7. Saya kenal dengan rgdal dan readOGR (). Namun, dengan ukuran shapefile ini, membaca seluruh shapefile ke dalam memori adalah proses yang sangat lambat. Selain itu, saya hanya tertarik …

10 shapefile  r  rgdal 

Saya mengubah vektor menjadi raster di R. Namun prosesnya terlalu lama. Apakah ada kemungkinan untuk menempatkan skrip ke dalam pemrosesan multithread atau GPU untuk melakukannya lebih cepat? Script saya untuk vektor raster. 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 

1. Pertanyaan Saya telah menemukan peringatan menggunakan fungsi projectRaster () dalam paket raster di R. Contoh lengkap yang dapat direproduksi ditempelkan di bawah ini. Warning message: In rgdal::rawTransform(projto_int, projfrom, nrow(xy), xy[, 1], : 33940 projected point(s) not finite Pertanyaan saya adalah: Apakah peringatan ini masalah yang harus saya perbaiki jika …

8 raster  r  coordinate-system  sp  rgdal