Je voudrais ajouter des attributs à un raster stocké dans une base de données postgis. J'utilise la commande suivante pour insérer le raster dans la base de données: raster2pgsql -s 4236 -I -C -M /home/user/mysrc/outlabeled.tif public.test | psql -d template_postgis Le raster outlabeled.tifest en...
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...
J'essaie de couper une grande ECW (détails ci-dessous) mais le fichier raster est trop gros pour être entièrement traité. Quelques détails de l'ECW ci-dessous Pilote: ondelettes compressées ECW / ERDAS (SDK 5.0) taille du fichier: 50 Go La taille est 450000, 565081 Taille des pixels: 0,15 0,15...
J'ai une image raster des élévations. J'ai ajouté un fichier de formes vectorielles contenant les divisions politiques locales. Je veux garder les couleurs d'élévation, mais je veux aussi garder les lignes noires représentant les divisions politiques. J'ai réduit la transparence du fichier de...
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...
J'essaie d'améliorer un processus vectoriel / python extrêmement lourd pour un modèle de risque naturel. À l'heure actuelle, nous avons un long script qui génère des lignes de distance / relèvement à partir d'un point donné pour déterminer: le type de polygone qu'il recoupe (p. ex. forêt, herbe,...
J'ai compilé et installé le maître QGIS actuel comme d'habitude et j'ai importé des données raster dans ma base de données postgres 9.3 avec l'extension postgis 2.0.2. Chaque fois que j'essaie d'importer les données raster dans qgis à l'aide du gestionnaire de base de données, je reçois ce message:...
J'ai deux couches .ecw superposées de la même zone, bien que l'une soit coupée au milieu, avec un bord blanc. Lorsque je définit la valeur d'arrière-plan sur 255,255,255 et que je la définis comme creuse, elle ne fait que creuser cette valeur RVB spécifique, alors qu'en réalité, il y a d'autres...
Existe-t-il des outils ou des méthodes dans ArcMap qui peuvent mesurer l'aire de chaque niveau de la sortie de densité du noyau? Notez que ce n'est qu'un tracé de densité du noyau (pas un fichier de forme ou un polygone). Analyse de la densité du noyau (vert foncé (10%) à rouge (90%) (1) J'ai...
Comment puis-je trouver une valeur de cellule maximale à partir d'un raster empilé. Rmax <- maxValue(RAD1998.all[[1]]) fonctionne bien mais Rmax <- maxValue(RAD1998.all[[2]]) donne NA. Bien sûr, pas dans un raster empilé. Voici mon code: RAD1998 <- raster(paste(getwd(),...
Je veux masquer un raster en fonction des valeurs d'un polygone: c'est-à-dire définir à NA toutes les valeurs du raster qui ne sont pas couvertes par les polygones. Cependant, il semble que raster::mask()les cellules soient masquées (c'est-à-dire qu'elles soient définies sur NA) même si elles sont...
Nous travaillons avec différents jeux de données raster, et les coordonnées des pixels sont stockées soit au centre du pixel, soit dans un coin. Les informations sont généralement fournies dans les métadonnées, mais parfois nous devons deviner et ce post illustre qu'il n'est pas toujours facile à...
Fermé . Cette question doit être plus ciblée . Il n'accepte pas actuellement les réponses. Vous souhaitez améliorer cette question? Mettez à jour la question pour qu'elle se concentre sur un seul problème en modifiant ce post . Fermé il y a 2 ans . Je construis un script en python en utilisant OGR...
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(),...
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...
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...
J'ai plusieurs rasters au format tiff qui ont une carte de couleurs qui leur est assignée: Je crée ensuite un ensemble de données mosaïque avec des spécifications correspondant à l'imagerie et y ajoute tous mes tiffs (vous pouvez ignorer les avertissements, c'est juste pour illustrer ce que j'ai...
Travailler avec QGIS 2.4.0 Chugiak J'ai une question concernant l'écrêtage d'images raster. J'ai une photo aérienne que je souhaite découper en utilisant l'étendue d'un fichier de formes. Jusqu'ici tout va bien, j'ai utilisé l'outil raster / extraction / clipper pour ce faire. Le problème est que...
J'utilise la méthode GDAL ReadAsArray pour travailler avec des données raster à l'aide de numpy (en particulier la reclassification). Comme mes rasters sont volumineux, je traite les tableaux en blocs, en répétant chaque bloc et en les traitant avec une méthode similaire à l' exemple GeoExamples ....