«r» etiketlenmiş sorular

Ggplot2 kullanarak uzamsal verilerimi çizerken sorun yaşıyorum. Harita spplot kullanılarak çizildiğinde iyi görünüyor, bu yüzden yırtılmanın takviye aşamasında gerçekleştiğini varsayıyorum. Kod aşağıdaki gibidir: #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 laa to regional mapping #must aggregate to higher level …

10 r  ggplot2  rgeos 

Bir dizi uydu görüntüsü (5 bant) var ve R'de kmeans tarafından sınıflandırmak istiyorum. Betiğim iyi çalışıyor (resimlerim arasında döngü, görüntüleri data.frame dönüştürmek, küme ve geri dönüştürmek raster): for (n in files) { image <- stack(n) image <- clip(image,subset) ###classify raster image.df <- as.data.frame(image) cluster.image <- kmeans(na.omit(image.df), 10, iter.max = 10, …

10 r  classification 

Ben değerleri sıfıra yakın bir harita oluşturmak levelplotiçin harika rasterVispaket işlevini kullanıyorum . Kırmızıdan maviye palet kullanarak bir tarama dosyası çiziyorum, ancak haritadaki sıfıra değerlere gri atamaya çalışıyorum. Özellikle, bu rakamın renklerini çoğaltmak istiyorum: Kırmızıdan maviye bir ölçek kullanıldığında, ancak sıfır değerlerinin gri ile renklendirildiğine dikkat edin. Şu anda haritalarım …

9 raster  r  legend  visualisation  color-ramp 

Kapalı. Bu soru konu dışı . Şu anda cevapları kabul etmiyor. Bu soruyu geliştirmek ister misiniz? Sorunuzu güncelleyin o yüzden konuyla Coğrafi Bilgi Sistemleri Stack Exchange için. 9 ay önce kapalı . 1000 farklı enlem-boylam içeren bir veri setim var. Bu koordinatların her biri için ortalama yıllık sıcaklık ve yıllık …

9 r  climate  extract 

R'de raster vektörü dönüştürüyorum. Ancak süreç çok uzundu. Daha hızlı yapmak için betiği çok iş parçacığına veya GPU işlemeye koyma olasılığı var mı? Rasterleştirilmiş vektöre benim senaryom. r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do raster res(r.raster) = 10 #definindo o tamanho do pixel setor.r = rasterize(setor, r.raster, …

9 raster  r  rasterization  rgdal  parallel 

Ggmap kullanarak, bazı konum noktaları içeren bir harita üzerinde bir şekil dosyasından belediye sınırlarını (çokgen) dahil etmek istiyorum. Bu komut dosyası, çokgeni çizmek dışında her şeyi yapar: 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")) # …

9 shapefile  r  ggmap 

Paketi choropleth()kullanarak R doğru harita oluşturmaya çalışıyorum GISTools. choro.legend()Efsanemi göstermek için kullanıyorum . Ancak her iki doğru efsane yerleşimini oluşturamıyorum VE / VEYA İşlevimdeki fmtparametrenin rolünü anlamıyorum choro.legend(). Bana öyle geliyor ki fmtefsane renklerim ve tanımım arasındaki boşluğu azaltabilir. Bulduğum choro.legend yardım sayfasını devletler bu fmtolduğunu "choropleth sınıf sınırların üzerinde …

9 r  legend  choropleth 

Noktaların raster değerinin 1 olduğu alanlarla sınırlandırılması gereken ikili bir raster içinde 10000 nokta (veya daha büyük) bir nokta vektör veri kümesi oluşturmak istiyorum. Aşağıdaki adımları denedim. Poligonize tarama QGIS: Vektör -> Araştırma Araçları -> Rastgele Noktalar Bu 2000 noktaya kadar iyi çalışır, ancak yukarıdaki herhangi bir şey QGIS'in çökmesine …

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

R'de bir data.frame oluşturdum. Sonuç olarak SpatialLinesDataFrame almak istiyorum, böylece satırlar, SpatialLines ve SpatialLinesDataFrame içine SpatialLines içine satırları veri dönüştürmek zorunda anlamak. bunu ben yaptım filedata.frame=data.frame(matrix(file),ncol=14, byrow=T)) file=Lines(filedata.frame) İstediğim data.frame'i alıyorum ama satırları değil. Eminim basit ama neden hatları alamıyorum?

9 r  line 

Bir Raster yansıtmaya çalışıyorum. R'de bunun projectRaster()işlevi vardır (tamamen tekrarlanabilir bir örneğin altında): # 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 = CRS(newproj), method = …

9 raster  gdal  r  coordinate-system  gdalwarp 

.Adf dosyalarını R'ye yüklemek istiyorum. Veriler bu sayfadan: http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057 İnternette bazı araştırmalardan sonra bulduğum aşağıdaki kodu denedim. Sorun şu ki, RasterLayer sınıfında olmaması gereken negatif değerler alıyorum. Bunun neden olduğunu bilmiyorum, umarım biri bana yardım edebilir !? Kod: library(rgdal) library(RColorBrewer) dpath<- path... x <- new("GDALReadOnlyDataset", dpath) getDriver(x) getDriverLongName(getDriver(x)) xx<-asSGDF_GROD(x) r …

9 raster  data  r  esri-grid-format 

İçinde bir dizi noktanın düştüğü çokgeni tanımlamak için uzamsal bir bindirme yapmak istiyorum. Bununla birlikte, ilk önce çokgenleri tamponlamak ve çözmek istiyorum, böylece birleştirilmiş çokgenler içinde herhangi bir yere düşen noktalar (ancak herhangi bir delik içinde değil) üst üste bindirme prosedürü ile benzer şekilde etiketlenecektir. Ne yazık ki, kullandığım tamponlama …

9 r  sp 

Belirtilen klasördeki tüm çoklu mozaik rasterlerini (her mozaik 7 katmandan oluşur) R'ye almak istiyorum. Sonra bunlara tek katmanlı rasterler olarak erişin. Ben de yaptım: # read all mosaics named "mos....img" in R raster_data<-list.files(path=getwd(), pattern="mos.*.img$") # read files as rasters s <- stack(raster_data) # check my imported rasters p.ex. raster n°8 …

9 raster  r  remote-sensing  landsat 

SpatialLinesR'de iki nesne yarattım . Bu nesneler şu şekilde yaratıldı: library(sp) xy <- cbind(x,y) xy.sp = sp::SpatialPoints(xy) spl1 <- sp::SpatialLines(list(Lines(Line(xy.sp), ID="a"))) Şimdi bir şekilde bunun aynı çizgi döndürülmüş ve çevrilmiş olduğu ve aralarındaki farkın 0'a eşit olduğu (yani şekil eşit olduğu) sonucuna varmak istiyorum. Bunu yapmak için, maptoolspaket kullanabilir ve …

9 r  geometry  similarity 

Danimarkalı belediyeler için bir şekil dosyasının bazı centroidlerini bulmaya çalışıyorum ve daha sonra aralarındaki sürüş süresini buluyorum. Kullandığım R'ın readShapePolydan fonksiyonu maptoolsile kombine gCentroidgelen fonksiyonu rgeosve her şey eserler. Ancak, SpatialPoints: x y 1 571860.7 6225016 Coordinate Reference System (CRS) arguments: NA Bu, seyahat sürelerini yakalamak için Google'da kullanabileceğim bir …

9 shapefile  r  coordinate-system