«r» etiketlenmiş sorular

Yığılmış bir taramadan maksimum hücre değerini nasıl bulabilirim. Rmax <- maxValue(RAD1998.all[[1]]) iyi çalışıyor ama Rmax <- maxValue(RAD1998.all[[2]]) NA verir. Tabii ki yığılmış bir taramada değil. İşte benim kod: RAD1998 <- raster(paste(getwd(), "/1998bil/1998ASC5min_ppt_spas1214_0001_19980202_0810_UTC.asc.bil", sep = "")) list.ras <- mixedsort(list.files(paste(getwd(), "/1998bil/", sep = ""), full.names = T, pattern = ".asc.bil")) RAD1998.all <- …

9 raster  r 

Tür dağılımlarının yüzey alanını hesapladım (şekil dosyalarından çokgenleri birleştirmek), ancak bu alan oldukça uzak çokgenlerden oluşabileceğinden, bir miktar dağılım ölçüsünü hesaplamak istiyorum. Şimdiye kadar yaptığım, her çokgenin centroidlerini geri almak, aralarındaki mesafeyi hesaplamak ve bunları aşağıdaki kukla örnekte olduğu gibi varyasyon katsayısını hesaplamak için kullanmaktır; require(sp) require(ggplot2) require(mapdata) require(gridExtra) require(scales) …

9 polygon  r  distance  centroids 

Herkes, rahatsızlık için özür dilerim, ama çok önemli bir zorlukla karşı karşıya kaldım: Farklı yıllar boyunca işsizlik değişiklikleri ile birlikte hareketli bir Russin haritası oluşturmak istiyorum. Resimde bir yıllık verileri görebilirsiniz require(sp) require(maptools) require(RColorBrewer) require(rgdal) rus<-url("http://www.filefactory.com/file/4h1hb5c1cw7r/n/RUS_adm1_RData") print(load(rus)) unempl <- read.delim2(file="C:\\unempl1.txt", header = TRUE, sep = ";",quote = "", dec=",", stringsAsFactors=F) …

9 cartography  r  animation 

Kesişen Poly2 Poly1 üzerine kaldırmaya çalışıyorum. Aşağıda 2 çokgenin koordinatları bulunmaktadır. > coords1 [,1] [,2] [1,] 992.0161 7462.531 [2,] 950.1962 7413.532 [3,] 902.3632 7421.207 [4,] 706.1985 7472.378 [5,] 654.0139 7711.760 [6,] 657.5960 7726.950 [7,] 786.2667 7844.380 [8,] 789.5935 7842.925 [9,] 981.7046 7498.659 [10,] 983.4246 7493.271 [11,] 990.6680 7469.888 [12,] 992.0161 …

9 r  polygon  intersection  symmetricdifference 

Küçük bir bahçe arsası için birim alan başına ortalama güneş yoğunluğunu gösteren bir tarama yüzeyi oluşturmakla ilgileniyorum. Modelin girdilerinin şu olacağını hayal ediyorum: Arsa alanı - basitlik uğruna, bahçe arsasının dairesel olduğunu varsayalım. Çevredeki orman kenarının ortalama ağaç yüksekliği Çeşitli güneş açılarında güneş yoğunluğu Zamanla bahçe arsasına göre ortalama güneş …

9 arcgis-desktop  python  arcpy  raster  r 

Ekteki görüntü, daire olarak temsil edilen kızıl çam ve haçlar olarak temsil edilen beyaz çam ile bir orman boşluğunu göstermektedir. İki çam ağacı türü arasında pozitif veya negatif bir ilişki olup olmadığını belirlemekle ilgileniyorum (yani aynı alanlarda büyüyorlar mı yoksa büyüyorlar mı). R spatstat paketindeki Kcross ve Kmulti'nin farkındayım. Bununla …

9 arcgis-desktop  python  arcpy  r  spatial-statistics 

Dokümantasyonunda fonksiyon argümanının sadece 2 sütun kabul spettiğini gördüm . 3 boyutlu çokgenler oluşturmak mümkün değil mi? Bir üçüncü mekansal boyut olarak zaman ekleyerek paketi kullanarak uzamsal-zamansal kriging yapıyorum . Bunu çokgenler ile kriging'i blok olarak engellemek için uzatmak istiyorum. Bunu yapmanın bir yolu var mı?coordsSpatialPolygonsgstatSpatialPointsDataFrame

9 r  interpolation  spatial-statistics  kriging  gstat 

Farklı tematik haritalar için ana ızgara olarak bir Mekansal Izgaraya ihtiyacım var. Tüm NA piksellerini atarak bir taramadan nasıl bir Mekansal Izgara üretebilirim?

9 r 

Çalışma Alanımda 100 SpatialPointsDataFrames var. Onları bu şekilde yükledim: filenames <- list.files(path="", pattern="XYhectareTravelTimes_ez+.*shp") for(i in filenames){ filepath <- file.path("/",i) assign(i, readShapePoints(filepath)) Bunlara şöyle denir: XYhectareTravelTimes_ez10.*shp XYhectareTravelTimes_ez11.*shp XYhectareTravelTimes_ez12.*shp vb. Bunları Çalışma Alanı'nda döngü yapan rasterlere nasıl dönüştürebilirim? Çok yeni bir R kullanıcısıyım ve yardım bulmayı umuyorum. çok teşekkürler.

9 r  rasterization 

Benim daha önceki değinen soruya kullanarak R - -, nasıl / dönüştürme "coğrafi kod" do harita üzerinde bir katman olarak gösterilebilir bir biçime oluşturulmuş bir raster görüntü? Kısacası, ticari SpatialKey ( ekran görüntüsü ) yazılımının sonuçlarını yeniden oluşturmayı hedefliyorum . Haritalar için verileri barındırmak için bir GeoServer örneği kullanıyorum, ancak …

9 raster  r  georeferencing 

Muhtemelen çok basit bir soru. Yaklaşık bin potansiyel coğrafi konumun (lat lon) ve "en çok yayılmış" 200 konumu seçmem gereken bir listem var. Bence bu toplam ortalama mesafenin en yüksek olduğu 200 puan. Bir şehirdeki mağazaları düşünün. Bunu yapmak için tanımlanmış bir yöntem var mı? Belki bir R paketinde? Bunu …

9 r  spatial-statistics  business 

Bu yığın akışından bir crosspost . GIS yazılımı hakkında çok az şey biliyorum ve çok basit bir soru ise, R'deki özür dileme işlemimi yapıyorum. Diyelim ki farklı kaynaklardan ancak farklı özelliklere sahip iki şekil dosyası var. Diyelim ki biri Texas'ın ( boundaries.shp) idari sınırları , diğeri ise Texas nehirleri ( …

9 r  overlay 

Csv biçiminde bir veri kümesi var, gibi: latitude, longitude, value -45, 45, 10 .... ArcGIS veya R + çimen kullanarak bu verilerden ilçe düzey araçlarını hesaplamak istiyorum. ArcGIS USA Counties katmanım var, ancak uzamsal bir birleşmeyi başarıyla tamamlayamadım. Bunu nasıl yapabilirim?

9 arcgis-desktop  data  geoprocessing  r