Questions marquées «gdal»

Je veux développer un portail Web qui montrera les images produites dynamiquement par QGIS, mais QGIS a fourni des images tiff sont trop grandes et il n'est pas possible de les montrer sur le site Web en utilisant des visionneuses d'images. C'est pourquoi j'ai décidé de convertir ces images tiff en...

J'utilise GDAL en Python pour importer des fichiers de formes. Je veux connaître les noms de champs pour le fichier, ma façon actuelle est: fields = [] for i in range(1, layer.GetFeature(0).GetFieldCount()): field = layer.GetFeature(0).GetDefnRef().GetFieldDefn(i).GetName() fields.append(field)...

Je viens d'exporter une table Postgis vers shp en utilisant ces conseils, mais je ne suis pas en mesure d'importer un shp vers Postgis en utilisant la même bibliothèque (ogr). Une idée? Merci beaucoup

Je fusionne des fichiers de formes et j'ai eu des problèmes à le faire dans QGIS, donc j'utilise directement ogr2ogr. Je fais ça (en lot): ogr2ogr -overwrite %destination% %n1% ogr2ogr -update -append %destination% %n2% -nln all_new ogr2ogr -update -append %destination% %n3% -nln all_new ogr2ogr...

J'ai un fichier image tif et un fichier mondial .tfw (créé à partir de données OSM avec Mapertive). Comment puis-je (en utilisant la ligne de commande linux) convertir cela en un fichier GeoTIFF? (c'est-à-dire que je veux avoir un fichier, pas 2). Je suis sûr que c'est possible avec gdal ou quelque...

Je suis extrêmement nouveau dans le SIG. J'utilise gdalpour lire une carte d'utilisation des terres / couverture terrestre et j'ai besoin de choisir le lat / lng de certains types de couverture terrestre à indexer dans un ensemble de données différent qui n'est exprimé qu'en lat / lng....

J'ai un raster entier pour lequel je voudrais construire une table d'attributs raster en utilisant Python et GDAL. Je peux créer une table d'attributs raster GDAL en Python comme suit: >>> rat = gdal.RasterAttributeTable() Cela fonctionne bien, comme nous pouvons le voir: >>> rat...

J'essaie d'optimiser mon code python en utilisant pypy . Bien que je puisse importer la plupart des bibliothèques et l'exécuter avec succès, je ne suis pas en mesure d'importer des bibliothèques GDAL. Les mêmes bibliothèques lors de l'importation en Python 2.7 fonctionnent correctement. J'ai essayé...

J'écris un utilitaire simple pour recadrer des lots de fichiers raster géotiff multi-bandes dans la même zone (plus petite). En utilisant gdalwarp, je peux facilement recadrer un fichier à l'aide d'un fichier de formes de découpage à un seul polygone: gdalwarp -cutline clipper.shp -crop_to_cutline...

Je voulais fusionner environ 20 fichiers de démonstration chacun autour de 100 Mo, et j'ai utilisé gdalwarp (je connais l'option de fichier virtuel): gdalwarp $(list_of_tiffs) merged.tiff Cependant, cette commande a pris très longtemps, et après environ 15 minutes, elle traitait la 5ème image de...

Je crée régulièrement mes propres rasters GeoTIFF avec GDAL en Python, par exemple: from osgeo import gdal from numpy import random data = random.uniform(0, 10, (300, 200)) driver = gdal.GetDriverByName('GTiff') ds = driver.Create('MyRaster.tif', 200, 300) band = ds.GetRasterBand(1)...

Voir ce lien pour plus de détails. Le problème: Je veux parcourir un raster continu (un qui n'a pas de table attributaire), cellule par cellule, et obtenir la valeur de la cellule. Je veux prendre ces valeurs et exécuter des conditions dessus, en émulant les étapes d'algèbre de carte détaillées...

J'essaie de traiter certaines données raster à l'aide d'ogr / gdal et je n'arrive pas à obtenir la pleine utilisation de tous les cœurs de ma machine. Lorsque je n'exécute le processus que sur un seul cœur, j'obtiens une utilisation à 100% de ce cœur. Lorsque j'essaie de diviser en multicœurs (dans...

Dans la documentation de gdal_calc, il est indiqué Calculatrice raster en ligne de commande avec syntaxe numpy . Plus tard, il y a quelques exemples où dans l'un d'eux: gdal_calc.py -A input.tif --outfile = result.tif --calc = "A * (A> 0)" --NoDataValue = 0 - signifie que les valeurs définies de...

Je veux créer deux polygones. L'une est l'étendue rectangulaire de mon raster. Je sais que je peux obtenir l'étendue du raster en utilisant r <- raster("band5.tif") e <- extent(r) plot(e)me donne ceci: Comment puis-je créer UN polygone de la limite du raster comme indiqué ci-dessous? Voici à...

Existe-t-il un moyen de définir tous les pixels de valeur 0 d'un raster sur nodata à l'aide de GDAL? Je voudrais utiliser fillnodata pour remplir certains trous, mais les trous sont définis sur 0, je dois donc les supprimer

J'ai installé GDAL et vérifié sur la ligne de commande, cela fonctionne et lorsque j'exécute gdal2tiles.py sur la ligne de commande, cela fonctionne mais lorsque j'exécute le même fichier gdal2tiles.py avec python IDLE, j'obtiens «Aucun module nommé» osgeo '" Erreur. Pouvez-vous...

J'ai des données raster dans une table postgres que je veux entrer dans python en tant que tableau numpy. J'utilise psycopg2 pour me connecter à la base de données. Je peux télécharger les données mais elles reviennent sous forme de chaîne (probablement un binaire sérialisé). Est-ce que quelqu'un...

Je suis surtout un vecteur de vecteurs, mais pour un projet en cours, je dois travailler avec un tas de vieilles cartes papier numérisées (cartes de dégâts de bombe ww2 pour Londres, si vous êtes intéressé!) Nous avons numérisé et géoréférencé les cartes, et souhaitons maintenant créer une couche...

J'ai du mal à utiliser la fonction Polygonize en python. L'exemple de livre de cuisine pour cela peut être trouvé ici . La partie pertinente de mon code est: sourceRaster = gdal.Open('myraster.tif') band = sourceRaster.GetRasterBand(1) bandArray = band.ReadAsArray() outShapefile = "polygonized"...