Je travaille sur un DEM qui dans une zone où il y a peu de soulagement. Je voudrais «brûler» le réseau fluvial dans le DEM pour pouvoir calculer avec précision l'accumulation et la longueur du débit. J'utilise ArcGIS Desktop
est-il possible de diviser un fichier de formes par fonctionnalité en python? (Le mieux serait une solution où je peux enregistrer temporairement les objets vectoriels résultants dans la mémoire à la place sur le disque). La raison: je veux utiliser la fonction gdal rasterizeLayer avec plusieurs...
J'ai un shapefile de ligne représentant un réseau routier. Je souhaite pixelliser ces données, les valeurs résultantes dans le raster indiquant la longueur totale des lignes qui se trouvent dans la cellule raster. Les données sont dans la projection du British National Grid, les unités seront donc...
J'essaie de graver un fichier de formes dans un raster en utilisant RasterizeLayer de GDAL. Je pré-crée une raster de zone d'intérêt à partir d'un fichier de formes différent, étant donné une taille de pixel spécifique. Cet AOI sert ensuite de base à toutes les rastérisations suivantes (même nombre...
Je veux créer un raster avec une résolution de 25 mètres × 25 mètres, où chaque cellule contient la distance jusqu'au littoral le plus proche, calculée à partir du centre de la cellule. Pour ce faire, je n'ai qu'un fichier de formes des côtes de la Nouvelle-Zélande . J'ai essayé de suivre le...
Il y a 100 SpatialPointsDataFrames dans mon espace de travail. Je les ai chargés de cette façon: filenames <- list.files(path="", pattern="XYhectareTravelTimes_ez+.*shp") for(i in filenames){ filepath <- file.path("/",i) assign(i, readShapePoints(filepath)) Ils sont appelés comme ceci:...
J'ai une couche de points qui montre la distribution d'une espèce. Est-il possible dans QGIS de convertir cette couche de points en une grille raster dans laquelle la valeur de chaque cellule correspond à la quantité de points dans cette cellule? Jusqu'à présent, je n'arrivais qu'à attacher des...
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...