J'essaie de fusionner 60 fichiers géotiff en un seul fichier énorme afin de pouvoir produire des tuiles à partir de ce fichier géotiff énorme résultant en utilisant gdal2tiles. Le processus de fusion a fonctionné, mais pour une raison quelconque, il produit un fichier d'une taille de 78 Go au lieu...
J'essaie de configurer un virtualenv python isolé pour travailler sur des projets SIG sur mon Mac OS X. On dirait que je devrai installer le GEOS, PROJ, GDAL / OGR à partir du site kyngchaos ici http://www.kyngchaos.com/ logiciels / frameworks et pas easy_install dans mon virtualenv. Est-ce exact?...
Je mosaïque certaines images avec gdal et j'aimerais améliorer le résultat final en utilisant une couche d'altération progressive / graduelle vers le bord de chaque image pour supprimer les bords nets au milieu de la mosaïque. Le problème que j'ai, c'est que la partie de chaque image individuelle...
J'ai essayé d'utiliser le code de QGis Save Raster as Rendered Image . iest une couche raster: pipelayer = i pipeextent = pipelayer.extent() pipewidth, pipeheight = (pipelayer.width(), pipelayer.height()) piperenderer = pipelayer.renderer() pipeprovider = pipelayer.dataProvider() crs =...
Je travaille sur une projection de sélection de site où je ne pourrai effectuer aucune observation in situ. Un élément clé du processus de sélection est la robustesse du terrain sous-jacent. Je prévois d'utiliser GDALdem pour créer les cartes de robustesse dérivées. J'ai lu à la fois le Wilson et...
J'écris parfois de courts scripts / instructions de travail pour des collègues qui ne sont pas très compétents en GDAL / programmation en général. Pour éviter tout problème avec les scripts, lors de la fusion / conversion de géotiffes avec GDAL, j'essaie simplement d'utiliser l'option BIGTIFF = YES...
J'utilise gSimplify(package rgeos) pour simplifier les géométries d'un fichier de formes. La fonction fonctionne bien, mais maintenant je ne peux pas écrire la sortie dans un nouveau fichier de formes. J'ai essayé plusieurs façons: writeOGR(simplyshape, file, driver="ESRI Shapefile", layer='test')...
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"...
J'essaie de convertir des données raster en données polygonales. Auparavant, j'utilisais le fichier shapefile.py pour créer des données ponctuelles sur des lat / lon individuels sur le raster. Mais cela créera une énorme donnée et l'interroger sera un gros problème. Donc, ce serait mieux si je suis...
Je travaille avec GDAL en Python et j'ai voulu que GDAL lève des exceptions pour que je puisse faire quelque chose d'utile avec elles. Par exemple, si un fichier ne peut pas être ouvert avec gdal.Open (), je voudrais lever une exception au lieu de gdal renvoyant None. Je voudrais juste un moyen de...
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...
Lorsque j'essaie d'installer Shapely sur mon ordinateur Windows 64 bits, je ne parviens pas à faire fonctionner la bibliothèque GEOS. Jusqu'à présent, j'ai exécuté le programme d'installation OSGeo4W à partir duquel j'ai installé GDAL (je pense que la bibliothèque geos est incluse dans ce package)....
J'ai une grille binaire Arc / Info --- spécifiquement, un raster d'accumulation de flux ArcGIS --- et je voudrais identifier toutes les cellules ayant une valeur spécifique (ou dans une plage de valeurs). En fin de compte, je voudrais un fichier de formes de points représentant ces cellules. Je...
Existe-t-il un moyen d'effectuer la même tâche que l'utilitaire gdalbuildvrt à l'aide des liaisons GDAL Python? Jusqu'à présent, je n'ai trouvé aucun autre moyen de le faire que de créer un vrt d'un seul ensemble de données et de modifier manuellement le xml. Je voudrais créer un vrt à partir de...
J'utilise Python 2.6.5 (32 bits) avec Numpy 1.3 et Gdal 1.9.1 installés sur Windows 7 64 bits. J'essaie de lire un ensemble de données raster de 800 Mo Imagine (.img) dans un tableau Numpy pour faire de l'algèbre raster, mais dès que j'exécute le code suivant, Python.exe se bloque. from osgeo...
J'essaie d'utiliser ogrinfo pour obtenir des détails sur un fichier de formes que j'ai téléchargé. Actuellement, la seule façon pour moi de le faire est de le charger dans QGIS et de cliquer manuellement pour trouver des informations dessus, comme ouvrir la table attributaire. Je veux juste pouvoir...
J'ai un ensemble de rasters (.vrt) avec des données quotidiennes sur l'humidité du sol. Je veux additionner les pixels de tous les rasters afin d'avoir une mesure par mois. Cependant, le monde entier n'est pas couvert chaque jour, ce qui se traduit par des valeurs de nodata aux endroits où des...
Il y a eu des tonnes de questions à ce sujet, mais aucune n'a résolu mon problème. J'ai une machine avec: Windows 7 x64 Python 3.4.3 x64 gdal installé: gdalinfo --version --> GDAL 1.11.4, released 2016/01/25; gdal a été installé à partir d'une roue GDAL-1.11.4-cp34-none-win_amd64.whld' ici Le...
ArcGIS Desktop a un ensemble de données de mosaïque raster , QGIS et ses amis ont des tables raster virtuelles . Un survol de la documentation révèle beaucoup de similitudes. Ils créent ou permettent tous deux une représentation alternative des données qui se trouvent ailleurs sans avoir à créer...
Passer de la liste de diffusion gdal-dev: Le lundi 2 septembre 2013 à 19 h 09, David Shean a écrit: Hi list, j'essaie de regrouper une série temporelle de rasters GTiff avec une projection / étendue / résolution identique à un seul fichier NetCDF pour la distribution. J'ai passé la dernière heure à...