Pertanyaan yang diberi tag «r»

Ketika saya ingin melakukan kriging itu hanya berfungsi kadang-kadang, tergantung pada nilai apa yang saya gunakan di datatable saya. Sebagai hasil dari fungsi krige yang saya dapatkan var1.pred: NA NA NA ...dan var1.var: NA NA NA ...(tetapi hanya ketika saya menggunakan nilai "salah" dalam data saya.) Sebagai contoh: ini berfungsi …

10 r  interpolation  kriging  gstat 

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 mengalami masalah dalam memplot data spasial saya menggunakan ggplot2. Peta terlihat baik-baik saja ketika diplot menggunakan spplot, jadi saya berasumsi bahwa robekan terjadi pada tahap fortifikasi. Kode tersebut adalah sebagai berikut: #install the packages library(rgdal) library(mapproj) library(raster) library(rgeos) library(ggplot2) library(plyr) if (!require(gpclib)) install.packages("gpclib", type="source") gpclibPermit() setwd("C:/Users/My Documents") #read in …

10 r  ggplot2  rgeos 

Saya memiliki serangkaian gambar satelit (5 band) dan ingin mengklasifikasikannya berdasarkan kmeans dalam R. Skrip saya berfungsi dengan baik (loop melalui gambar saya, mengonversi gambar ke data.frame, mengelompokkannya, dan mengubahnya kembali menjadi raster): for (n in files) { image <- stack(n) image <- clip(image,subset) ###classify raster image.df <- as.data.frame(image) cluster.image …

10 r  classification 

Saya membaca posting tentang peta interaktif dengan R menggunakan leafletpaket. Pada artikel ini, penulis membuat peta panas seperti ini: X=cbind(lng,lat) kde2d <- bkde2D(X, bandwidth=c(bw.ucv(X[,1]),bw.ucv(X[,2]))) x=kde2d$x1 y=kde2d$x2 z=kde2d$fhat CL=contourLines(x , y , z) m = leaflet() %>% addTiles() m %>% addPolygons(CL[[5]]$x,CL[[5]]$y,fillColor = "red", stroke = FALSE) Saya tidak terbiasa dengan bkde2Dfungsinya, …

10 r  leaflet  heat-map 

Saya menggunakan levelplotfungsi rasterVispaket luar biasa untuk membuat peta yang nilainya berbeda sekitar nol. Saya merencanakan file raster menggunakan palet merah ke biru, tetapi saya mencoba untuk menetapkan nilai abu-abu ke nol di peta. Secara khusus, saya ingin mereproduksi warna gambar ini: Perhatikan bahwa di mana skala merah ke biru …

9 raster  r  legend  visualisation  color-ramp 

Tutup. Pertanyaan ini di luar topik . Saat ini tidak menerima jawaban. Ingin meningkatkan pertanyaan ini? Perbarui pertanyaan sehingga sesuai topik untuk Sistem Informasi Geografis. Ditutup 9 bulan lalu . Saya memiliki satu set data dengan 1000 lintang-bujur yang berbeda. Saya ingin mengekstraksi suhu tahunan rata-rata dan curah hujan tahunan …

9 r  climate  extract 

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 

Menggunakan ggmap, saya ingin memasukkan batas kota (poligon) dari shapefile pada peta yang berisi beberapa titik lokasi. Script ini melakukan segalanya kecuali merencanakan poligon: library(rgdal) library(ggmap) # Get shapefile with Drammen municipality borders tmpzip<-tempfile() tmpdir<-tempfile() dir.create(tmpdir) download.file("http://www.kartverket.no/Documents/Kart/N50-N5000%20Kartdata/33_N5000_shape.zip",tmpzip) unzip(tmpzip, exdir=tmpdir) kommune <- readOGR(dsn=tmpdir, layer="NO_AdminOmrader_pol") kommune<-kommune[kommune$NAVN=="Drammen",] kommune<-spTransform(kommune, CRS("+init=epsg:4326")) # Get location point …

9 shapefile  r  ggmap 

Saya mencoba membuat choropleth()peta yang benar di R menggunakan GISToolspaket. Saya gunakan choro.legend()untuk menampilkan legenda saya. Namun saya tidak dapat membuat penempatan legenda yang benar, DAN / ATAU Saya tidak mengerti peran fmtparameter dalam choro.legend()fungsi saya . Sepertinya saya yang fmtbisa mengurangi jarak antara warna dan deskripsi legenda saya. Saya …

9 r  legend  choropleth 

Saya telah membuat data.frame di R. Saya ingin mendapatkan SpatialLinesDataFrame sebagai hasilnya, jadi saya mengerti bahwa saya harus mengubah data.frame menjadi garis, garis menjadi SpatialLines dan SpatialLines menjadi SpatialLinesDataFrame. Saya melakukan ini filedata.frame=data.frame(matrix(file),ncol=14, byrow=T)) file=Lines(filedata.frame) Saya mendapatkan data.frame yang saya inginkan tetapi tidak garis. Saya yakin dasar tetapi mengapa saya …

9 r  line 

Saya ingin membuat dataset vektor titik 10.000 poin (atau lebih besar) dalam raster biner, di mana titik harus dibatasi ke area di mana nilai raster adalah 1. Saya mencoba langkah-langkah berikut. Polygonize raster QGIS: Vektor -> Alat Penelitian -> Poin Acak Ini berfungsi dengan baik hingga 2000 poin tetapi apa …

9 qgis  arcgis-desktop  python  r  point-in-polygon 

Saya mencoba memproyeksikan Raster. Di R ada projectRaster()fungsi untuk ini (di bawah contoh yang sepenuhnya direproduksi): # example Raster require(raster) r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40) r <- setValues(r, 1:ncell(r)) projection(r) # project to newproj <- "+init=epsg:4714" # using raster package to reproject pr1 <- projectRaster(r, crs = …

9 raster  gdal  r  coordinate-system  gdalwarp 

Saya ingin memuat file .adf ke dalam R. Data berasal dari halaman ini: http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057 Saya mencoba kode berikut yang saya temukan setelah beberapa penelitian di internet. Masalahnya adalah, bahwa di kelas RasterLayer saya mendapatkan nilai negatif yang seharusnya tidak ada di sana. Aku tidak tahu mengapa ini terjadi, jadi semoga …

9 raster  data  r  esri-grid-format 

Saya ingin melakukan overlay spasial untuk mengidentifikasi poligon di mana satu set titik jatuh. Namun saya pertama-tama ingin buffer dan melarutkan poligon, sehingga poin jatuh di mana saja di dalam poligon gabungan (tetapi tidak dalam lubang) akan diberi label yang sama dengan prosedur overlay. Sayangnya, proses buffering dan / atau …

9 r  sp