Pytania otagowane jako r

Mam problem z wykreśleniem moich danych przestrzennych za pomocą ggplot2. Mapa wygląda dobrze po wykreśleniu za pomocą spplotu, więc zakładam, że łzawienie występuje na etapie wzmocnienia. Kod jest następujący: #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 …

10 r  ggplot2  rgeos 

Mam szereg zdjęć satelitarnych (5 pasm) i chcę je sklasyfikować według kmeanów w R. Mój skrypt działa dobrze (przeglądaj moje obrazy, konwertuj obrazy na data.frame, grupuj je i konwertuj z powrotem na 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 

Korzystam z levelplotfunkcji niesamowitego rasterVispakietu, aby stworzyć mapę, której wartości różnią się od zera. Planuję plik rastrowy przy użyciu palety od czerwonej do niebieskiej, ale próbuję przypisać wartości szare do zera na mapie. W szczególności chcę odtworzyć kolory tej figury: Zauważ, że tam, gdzie używana jest skala od czerwonej do …

9 raster  r  legend  visualisation  color-ramp 

Zamknięte. To pytanie jest nie na temat . Obecnie nie przyjmuje odpowiedzi. Chcesz poprawić to pytanie? Zaktualizuj pytanie, aby było na temat wymiany stosów systemów informacji geograficznej. Zamknięte 9 miesięcy temu . Mam zestaw danych zawierający 1000 różnych szerokości i długości geograficznych. Chcę wyliczyć średnią roczną temperaturę i roczne opady …

9 r  climate  extract 

Przekształcam wektor na raster w R. Jednak proces ten był zbyt długi. Czy istnieje możliwość włączenia skryptu w proces wielowątkowy lub przetwarzanie GPU, aby zrobić to szybciej? Mój skrypt do zrasteryzowanego wektora. 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 

Korzystając z ggmap, chcę dołączyć granice gminy (wielokąt) z pliku kształtu na mapie zawierającej pewne punkty lokalizacji. Ten skrypt robi wszystko poza wykreślaniem wielokąta: 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 

Próbuję utworzyć poprawną choropleth()mapę w R za pomocą GISToolspakietu. Używam choro.legend()do wyświetlania mojej legendy. Nie mogę jednak utworzyć poprawnego położenia legendy ORAZ / LUB Nie rozumiem roli fmtparametru w mojej choro.legend()funkcji. Wydaje mi się, że fmtmogę zmniejszyć odstęp między kolorami mojej legendy a opisem. Znalazłem choro.legend stronę pomocy , by …

9 r  legend  choropleth 

Chcę utworzyć zestaw danych wektora punktowego o wartości 10000 punktów (lub więcej) w binarnym rastrze, w którym punkty powinny być ograniczone do obszarów, w których wartość rastra wynosi 1. Próbowałem następujące kroki. Poligonizuj raster QGIS: Wektor -> Narzędzia badawcze -> Losowe punkty Działa to dobrze do 2000 punktów, ale wszystko …

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

Utworzyłem ramkę danych w języku R. Chcę w rezultacie uzyskać SpatialLinesDataFrame, więc rozumiem, że muszę przekonwertować data.frame na linie, linie na SpatialLines, a SpatialLines na SpatialLinesDataFrame. ja to zrobiłem filedata.frame=data.frame(matrix(file),ncol=14, byrow=T)) file=Lines(filedata.frame) Dostaję data.frame chcę, ale nie linie. Jestem pewien, że to podstawa, ale dlaczego nie mogę uzyskać linii?

9 r  line 

Próbuję wyświetlić raster. W R jest do tego projectRaster()funkcja (poniżej w pełni odtwarzalnego przykładu): # 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 

Chcę załadować pliki .adf do R. Dane pochodzą z tej strony: http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057 Wypróbowałem następujący kod, który znalazłem po kilku badaniach w Internecie. Problem polega na tym, że w klasie RasterLayer otrzymuję wartości ujemne, których nie powinno tam być. Nie wiem, dlaczego tak się dzieje, więc mam nadzieję, że ktoś może …

9 raster  data  r  esri-grid-format 

Chciałbym zaimportować wszystkie wielokrotne rastry mozaikowe (każda mozaika składa się z 7 warstw) ze wskazanego folderu do mojego R. Następnie uzyskać do nich dostęp jako pojedyncze rastry wielowarstwowe. Więc zrobiłem: # 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 …

9 raster  r  remote-sensing  landsat 

Chciałbym wykonać nakładkę przestrzenną, aby zidentyfikować wielokąt, w którym mieści się zbiór punktów. Najpierw jednak chcę buforować i rozpuszczać wielokąty, tak aby punkty mieszczące się w dowolnym miejscu w połączonych wielokątach (ale nie w żadnych dziurach) zostaną podobnie oznaczone za pomocą procedury nakładania. Niestety proces buforowania i / lub rozpuszczania, …

9 r  sp 

Stworzyłem dwa SpatialLinesobiekty w R: . Te obiekty zostały utworzone w ten sposób: library(sp) xy <- cbind(x,y) xy.sp = sp::SpatialPoints(xy) spl1 <- sp::SpatialLines(list(Lines(Line(xy.sp), ID="a"))) Teraz chcę jakoś dojść do wniosku, że jest to ta sama linia obrócona i odwrócona, a różnica między nimi jest równa 0 (tzn. Kształt jest równy). …

9 r  geometry  similarity 

Próbuję znaleźć centroidy pliku kształtu dla duńskich gmin, a następnie znaleźć czas między nimi. Używam R„s readShapePolyfunkcję z maptoolspołączeniu z gCentroidfunkcją z rgeos, i wszystko działa. Jednak dostaję punkty przestrzenne takie jak SpatialPoints: x y 1 571860.7 6225016 Coordinate Reference System (CRS) arguments: NA Które oczywiście nie jest czymś, czego …

9 shapefile  r  coordinate-system