Lorsque je veux effectuer le krigeage, cela ne fonctionne que parfois, selon les valeurs que j'utilise dans ma table de données. En raison de la fonction krige que j'obtiens pour var1.pred: NA NA NA ...et var1.var: NA NA NA ...(mais uniquement lorsque j'utilise les "mauvaises" valeurs dans ma table de …
J'ai un très grand fichier de formes ESRI (> 8000 polygones, le fichier .shp est> 32 Mo), que je veux lire dans R sur Windows 7. Je connais rgdal et readOGR (). Cependant, avec cette taille de fichier de formes, la lecture de l'ensemble du fichier de formes en mémoire …
J'ai du mal à tracer mes données spatiales à l'aide de ggplot2. La carte semble bien lorsqu'elle est tracée à l'aide de spplot, donc je suppose que la déchirure se produit au stade de fortification. Le code est comme suit: #install the packages library(rgdal) library(mapproj) library(raster) library(rgeos) library(ggplot2) library(plyr) if …
J'ai une série temporelle d'images satellite (5 bandes) et je souhaite les classer par kmeans dans R. Mon script fonctionne correctement (parcourir mes images, convertir les images en data.frame, les regrouper et les reconvertir en un raster): for (n in files) { image <- stack(n) image <- clip(image,subset) ###classify raster …
J'ai lu un article sur les cartes interactives avec R en utilisant le leafletpackage. Dans cet article, l'auteur a créé une carte thermique comme celle-ci: 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 = …
J'utilise la levelplotfonction du rasterVispackage génial pour créer une carte dont les valeurs divergent autour de zéro. Je trace un fichier raster à l'aide d'une palette rouge à bleu, mais j'essaie d'attribuer des valeurs de gris à zéro dans la carte. Plus précisément, je veux reproduire les couleurs de cette …
Fermé. Cette question est hors sujet . Il n'accepte pas actuellement les réponses. Vous souhaitez améliorer cette question? Mettez à jour la question afin qu'elle soit sur le sujet pour l'échange de piles de systèmes d'information géographique. Fermé il y a 9 mois . J'ai un ensemble de données avec …
Je convertis le vecteur en raster en R. Cependant, le processus était trop long. Est-il possible de mettre le script en traitement multithread ou GPU afin de le faire plus rapidement? Mon script au vecteur tramé. r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do raster res(r.raster) = 10 …
En utilisant ggmap, je veux inclure les frontières de la municipalité (polygone) à partir d'un fichier de formes sur une carte contenant des points d'emplacement. Ce script fait tout sauf tracer le polygone: 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 <- …
J'essaie de créer une choropleth()carte correcte dans R en utilisant le GISToolspackage. J'utilise choro.legend()pour afficher ma légende. Cependant, je ne peux pas créer de placement de légende correct, ET / OU Je ne comprends pas le rôle du fmtparamètre dans ma choro.legend()fonction. Il me semble que cela fmtpeut diminuer l'espace …
J'ai créé un data.frame dans R. Je veux obtenir SpatialLinesDataFrame en conséquence, donc je comprends que je dois convertir le data.frame en lignes, les lignes en SpatialLines et les SpatialLines en SpatialLinesDataFrame. j'ai fait ça filedata.frame=data.frame(matrix(file),ncol=14, byrow=T)) file=Lines(filedata.frame) J'obtiens le data.frame que je veux mais pas les lignes. Je suis …
Je veux créer un jeu de données vectorielles ponctuelles de 10000 points (ou plus) dans un raster binaire, où les points doivent être contraints aux zones où la valeur du raster est 1. J'ai essayé les étapes suivantes. Polygoniser un raster QGIS: Vector -> Outils de recherche -> Points aléatoires …
J'essaie de projeter un raster. Dans R, il y a la projectRaster()fonction à cela (ci-dessous un exemple entièrement reproductible): # 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 <- …
Je souhaite charger des fichiers .adf dans R. Les données proviennent de cette page: http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057 J'ai essayé le code suivant que j'ai trouvé après quelques recherches sur Internet. Le problème est que, dans la classe RasterLayer, j'obtiens des valeurs négatives qui ne devraient pas être là. Je ne sais pas …
Je voudrais importer tous les rasters de mosaïques multiples (chaque mosaïque est composée de 7 couches) du dossier indiqué dans mon R. Ensuite, accédez-y en tant que rasters multicouches uniques. J'ai donc fait: # read all mosaics named "mos....img" in R raster_data<-list.files(path=getwd(), pattern="mos.*.img$") # read files as rasters s <- …
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.