Sto riscontrando problemi nel tracciare i miei dati spaziali usando ggplot2. La mappa ha un bell'aspetto quando viene tracciata usando spplot, quindi presumo che lo strappo avvenga nella fase fortificata. Il codice è il seguente: #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 …
Ho una serie temporale di immagini satellitari (5 bande) e desidero classificarle in base ai kmean in R. La mia sceneggiatura funziona bene (passa in rassegna le mie immagini, converte le immagini in data.frame, raggruppale e riconvertile in un raster): for (n in files) { image <- stack(n) image <- …
Sto usando la levelplotfunzione del fantastico rasterVispacchetto per creare una mappa i cui valori divergono intorno allo zero. Sto disegnando un file raster usando una tavolozza da rossa a blu, ma sto cercando di assegnare valori da grigio a zero nella mappa. In particolare, voglio riprodurre i colori di questa …
Chiuso. Questa domanda è fuori tema . Al momento non accetta risposte. Vuoi migliorare questa domanda? Aggiorna la domanda in modo che sia in argomento per lo scambio di stack dei sistemi di informazione geografica. Chiuso 9 mesi fa . Ho un set di dati con 1000 diverse latitudini-longitudini. Desidero …
Sto convertendo il vettore in raster in R. Tuttavia il processo è stato troppo lungo. Esiste la possibilità di mettere lo script in elaborazione multithread o GPU per farlo più velocemente? La mia sceneggiatura al vettore rasterizzato. r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do raster res(r.raster) = …
Usando ggmap, voglio includere i bordi del comune (poligono) da un file di forma su una mappa contenente alcuni punti di posizione. Questo script fa tutto tranne che stampare il poligono: 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") …
Sto cercando di creare una choropleth()mappa corretta in R usando il GISToolspacchetto. Uso choro.legend()per visualizzare la mia leggenda. Tuttavia, non riesco a creare il corretto posizionamento della legenda, E / O Non capisco il ruolo del fmtparametro nella mia choro.legend()funzione. Mi sembra che fmtpossa diminuire lo spazio tra i colori …
Voglio creare un set di dati vettoriali di punti di 10000 punti (o superiore) all'interno di un raster binario, in cui i punti dovrebbero essere vincolati ad aree in cui il valore raster è 1. Ho provato i seguenti passaggi. Poligono raster QGIS: Vector -> Strumenti di ricerca -> Punti …
Ho creato un data.frame in R. Voglio ottenere SpatialLinesDataFrame come risultato, quindi capisco che devo convertire data.frame in linee, le linee in SpatialLines e SpatialLines in SpatialLinesDataFrame. l'ho fatto filedata.frame=data.frame(matrix(file),ncol=14, byrow=T)) file=Lines(filedata.frame) Ottengo il data.frame che voglio ma non le linee. Sono sicuro che è di base, ma perché non …
Sto cercando di proiettare un raster. In R c'è la projectRaster()funzione per fare ciò (sotto un esempio completamente riproducibile): # 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 <- …
Voglio caricare i file .adf in R. I dati provengono da questa pagina: http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057 Ho provato il seguente codice che ho trovato dopo alcune ricerche su Internet. Il problema è che nella classe RasterLayer ottengo valori negativi che non dovrebbero esserci. Non so perché questo accada, quindi spero che qualcuno …
Vorrei eseguire una sovrapposizione spaziale per identificare il poligono all'interno del quale rientrano una serie di punti. Tuttavia, per prima cosa voglio bufferizzare e dissolvere i poligoni, in modo tale che i punti che ricadono ovunque all'interno dei poligoni uniti (ma non all'interno di alcun foro) saranno etichettati in modo …
Vorrei importare tutti i raster di mosaico multipli (ogni mosaico è composto da 7 strati) dalla cartella indicata nel mio R. Quindi accedervi come singoli raster multistrato. Così ho fatto: # 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 …
Ho creato due SpatialLinesoggetti in R: . Questi oggetti sono stati creati in questo modo: library(sp) xy <- cbind(x,y) xy.sp = sp::SpatialPoints(xy) spl1 <- sp::SpatialLines(list(Lines(Line(xy.sp), ID="a"))) Ora voglio in qualche modo concludere che questa è la stessa linea ruotata e capovolta e che la differenza tra loro è uguale a …
Sto cercando di trovare alcuni centroidi di uno shapefile per i comuni danesi e successivamente trovare il tempo di guida tra di loro. Uso Rla readShapePolyfunzione di maptoolscombinata con la gCentroidfunzione di rgeose tutto funziona. Tuttavia, ottengo punti spaziali come SpatialPoints: x y 1 571860.7 6225016 Coordinate Reference System (CRS) …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.