Preguntas etiquetadas con r

Tengo un archivo de forma ESRI muy grande (> 8000 polígonos, el archivo .shp es> 32 MB), que quiero leer en R en Windows 7. Estoy familiarizado con rgdal y readOGR (). Sin embargo, con este tamaño de shapefile, leer todo el shapefile en la memoria es un proceso muy …

10 shapefile  r  rgdal 

Tengo problemas para trazar mis datos espaciales usando ggplot2. El mapa se ve bien cuando se traza usando spplot, por lo que supongo que el desgarro se produce en la etapa de fortificación. El código es el siguiente: #install the packages library(rgdal) library(mapproj) library(raster) library(rgeos) library(ggplot2) library(plyr) if (!require(gpclib)) install.packages("gpclib", …

10 r  ggplot2  rgeos 

Tengo una serie temporal de imágenes de satélite (5 bandas) y quiero clasificarlas por kmeans en R. Mi script está funcionando bien (recorrer mis imágenes, convertir las imágenes en data.frame, agruparlas y convertirlas nuevamente en un trama): for (n in files) { image <- stack(n) image <- clip(image,subset) ###classify raster …

10 r  classification 

Leí una publicación sobre mapas interactivos con R usando el leafletpaquete. En este artículo, el autor crea un 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) No …

10 r  leaflet  heat-map 

Soy nuevo en R. ¿Cómo puedo reemplazar los NA con 0 para mis datos ráster que tienen información espacial?

10 raster  r  missing-data  replace 

Estoy usando la levelplotfunción del rasterVispaquete impresionante para crear un mapa cuyos valores divergen alrededor de cero. Estoy trazando un archivo ráster con una paleta de rojo a azul, pero estoy tratando de asignar valores grises a cero en el mapa. Específicamente, quiero reproducir los colores de esta figura: Observe …

9 raster  r  legend  visualisation  color-ramp 

Cerrado. Esta pregunta está fuera de tema . Actualmente no está aceptando respuestas. ¿Quieres mejorar esta pregunta? Actualice la pregunta para que sea sobre el tema para el intercambio de pila de sistemas de información geográfica. Cerrado hace 9 meses . Tengo un conjunto de datos con 1000 latitudes y …

9 r  climate  extract 

Estoy convirtiendo vector a ráster en R. Sin embargo, el proceso fue demasiado largo. ¿Existe la posibilidad de poner el script en procesamiento multiproceso o GPU para hacerlo más rápido? Mi guión al vector rasterizado. r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do raster res(r.raster) = 10 #definindo …

9 raster  r  rasterization  rgdal  parallel 

Usando ggmap, quiero incluir los bordes del municipio (polígono) de un archivo de forma en un mapa que contiene algunos puntos de ubicación. Este script hace todo excepto trazar el 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") …

9 shapefile  r  ggmap 

Estoy tratando de crear el choropleth()mapa correcto en R usando el GISToolspaquete. Yo uso choro.legend()para mostrar mi leyenda. Sin embargo, no puedo crear la ubicación correcta de la leyenda, Y / O NO entiendo el papel del fmtparámetro en mi choro.legend()función. Me parece que fmtpuede disminuir el espacio entre los …

9 r  legend  choropleth 

He creado un data.frame en R. Como resultado, quiero obtener SpatialLinesDataFrame, por lo que entiendo que tengo que convertir data.frame en líneas, las líneas en SpatialLines y SpatialLines en SpatialLinesDataFrame. hice esto filedata.frame=data.frame(matrix(file),ncol=14, byrow=T)) file=Lines(filedata.frame) Obtengo el data.frame que quiero pero no las líneas. Estoy seguro de que es básico, …

9 r  line 

Quiero crear un conjunto de datos vectoriales de puntos de 10000 puntos (o más grandes) dentro de un ráster binario, donde los puntos deben estar restringidos a áreas donde el valor del ráster es 1. Intenté los siguientes pasos. Poligonalizar ráster QGIS: Vector -> Herramientas de investigación -> Puntos aleatorios …

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

Estoy tratando de proyectar una trama. En R existe la projectRaster()función para esto (debajo de un ejemplo completamente reproducible): # 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 <- …

9 raster  gdal  r  coordinate-system  gdalwarp 

Quiero cargar archivos .adf en R. Los datos provienen de esta página: http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057 Intenté el siguiente código que encontré después de algunas investigaciones en Internet. El problema es que en la clase RasterLayer obtengo valores negativos que no deberían estar allí. ¡No sé por qué sucede esto, así que espero …

9 raster  data  r  esri-grid-format 

Me gustaría importar todos los rásteres de mosaico múltiples (cada mosaico está compuesto por 7 capas) de la carpeta indicada en mi R. Luego, acceda a ellos como rásteres de múltiples capas individuales. Así que lo hice: # read all mosaics named "mos....img" in R raster_data<-list.files(path=getwd(), pattern="mos.*.img$") # read files …

9 raster  r  remote-sensing  landsat