Perguntas com a marcação «r»

Estou com problemas para plotar meus dados espaciais usando o ggplot2. O mapa parece bom quando plotado usando spplot, então estou assumindo que o rasgo ocorre no estágio de fortificação. O código é o seguinte: #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 …

10 r  ggplot2  rgeos 

Eu tenho uma série temporal de imagens de satélite (5 bandas) e quero classificá-las por kmeans em R. Meu script está funcionando bem (faça um loop pelas minhas imagens, converta as imagens em data.frame, agrupe-as e converta-as novamente em varredura): for (n in files) { image <- stack(n) image <- …

10 r  classification 

Eu li um post sobre mapas interativos com R usando o leafletpacote. Neste artigo, o autor cria um mapa de calor como este: 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) Eu …

10 r  leaflet  heat-map 

Sou novo em R. Como posso substituir NAs por 0 para meus dados rasterizados que possuem informações espaciais?

10 raster  r  missing-data  replace 

Estou usando a levelplotfunção do rasterVispacote incrível para criar um mapa cujos valores divergem em torno de zero. Estou plotando um arquivo raster usando uma paleta de vermelho para azul, mas estou tentando atribuir valores de cinza a zero no mapa. Especificamente, quero reproduzir as cores desta figura: Observe que …

9 raster  r  legend  visualisation  color-ramp 

Fechadas. Esta questão está fora de tópico . No momento, não está aceitando respostas. Deseja melhorar esta pergunta? Atualize a pergunta para que ela esteja no tópico do Geographic Information Systems Stack Exchange. Fechado há 8 meses . Eu tenho um conjunto de dados com 1000 latitudes e longitudes diferentes. …

9 r  climate  extract 

Eu tenho um arquivo de formato ESRI muito grande (> 8000 polígonos, o arquivo .shp é> 32 MB), que eu quero ler no R no Windows 7. Estou familiarizado com rgdal e readOGR (). No entanto, com esse tamanho de arquivo de forma, a leitura de todo o arquivo de …

9 shapefile  r  rgdal 

Estou convertendo vetor para varredura em R. No entanto, o processo foi muito longo. Existe a possibilidade de colocar o script no processamento multithread ou GPU para fazer isso mais rapidamente? Meu script para vetor rasterizado. r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do raster res(r.raster) = 10 …

9 raster  r  rasterization  rgdal  parallel 

Usando o ggmap, quero incluir bordas de município (polígono) de um shapefile em um mapa que contém alguns pontos de localização. Este script faz tudo, exceto a plotagem do polígono: 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",] …

9 shapefile  r  ggmap 

Estou tentando criar o choropleth()mapa correto em R usando o GISToolspacote. Eu uso choro.legend()para exibir minha legenda. No entanto, não consigo criar o posicionamento correto da legenda, E / OU Não entendo o papel do fmtparâmetro na minha choro.legend()função. Parece-me que fmtpode diminuir o espaço entre as cores e a …

9 r  legend  choropleth 

Quero criar um conjunto de dados de vetores de pontos de 10000 pontos (ou maiores) em uma varredura binária, em que os pontos devem ser restritos a áreas onde o valor da varredura é 1. Eu tentei os seguintes passos. Polygonize raster QGIS: Vetor -> Ferramentas de pesquisa -> Pontos …

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

Eu criei um data.frame no R. Quero obter SpatialLinesDataFrame como resultado, para que eu entenda que preciso converter o data.frame em linhas, as linhas em SpatialLines e as SpatialLines em SpatialLinesDataFrame. eu fiz isso filedata.frame=data.frame(matrix(file),ncol=14, byrow=T)) file=Lines(filedata.frame) Eu recebo o data.frame que eu quero, mas não as linhas. Tenho certeza …

9 r  line 

Estou tentando projetar uma varredura. Em R existe a projectRaster()função para isso (abaixo de um exemplo totalmente reprodutível): # 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, …

9 raster  gdal  r  coordinate-system  gdalwarp 

Desejo carregar arquivos .adf no R. Os dados são desta página: http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057 Eu tentei o seguinte código que encontrei depois de algumas pesquisas na internet. O problema é que, na classe RasterLayer, recebo valores negativos que não deveriam estar lá. Não sei por que isso acontece, então espero que alguém …

9 raster  data  r  esri-grid-format 

Gostaria de importar todos os múltiplos rasters de mosaico (cada mosaico é composto por 7 camadas) da pasta indicada para o meu R. Em seguida, acesse-os como rasters de multicamadas únicas. Então eu fiz: # read all mosaics named "mos....img" in R raster_data<-list.files(path=getwd(), pattern="mos.*.img$") # read files as rasters s …

9 raster  r  remote-sensing  landsat