Questions marquées «numpy»

Je suis nouveau au SIG. J'ai un code qui convertit les images infrarouges de Mars en cartes d'inertie thermique, qui sont ensuite stockées sous forme de tableaux numpy 2D. J'ai enregistré ces cartes en tant que fichiers hdf5 mais j'aimerais vraiment les enregistrer en tant qu'images raster afin de...

J'ai essayé de vérifier mes filtres sur le raster DEM pour la reconnaissance des formes et il en résulte toujours des dernières lignes et colonnes manquantes (comme..20) . J'ai essayé avec la bibliothèque PIL, la charge d'image. Puis avec numpy. La sortie est la même. Je pensais que quelque chose...

J'essaie de re-projeter / rééchantillonner avec les liaisons Python GDAL, mais j'obtiens des résultats légèrement différents par rapport à ceux de l'utilitaire de ligne de commande gdalwarp. Voir la mise à jour ci-dessous pour un exemple plus court Ce script illustre l'approche Python: from osgeo...

Je souhaite apprendre à utiliser les tableaux NumPy pour optimiser le géotraitement. Une grande partie de mon travail implique des «mégadonnées», où le géotraitement prend souvent des jours pour accomplir certaines tâches. Inutile de dire que je suis très intéressé par l'optimisation de ces...

J'ai commencé à apprendre à manipuler les données LAS en python et je voulais voir comment les autres gèrent les fichiers LAS. Je voudrais lire les points (j'utilise un tableau numpy) et filtrer les classes 1 et 2 (non classifiées et au sol) dans un tableau séparé. J'ai le code suivant mais je...

Pour obtenir une valeur de pixel unique à un point d'un raster à l'aide de rasterio, il existe un exemple ici: https://github.com/mapbox/rasterio/pull/275 Cependant, existe-t-il une API directe dans rasterio (et non le cli) qui peut être utilisée pour extraire de la valeur à un seul point dans un...

J'utilise ArcGIS 10.1 et je souhaite créer un nouveau raster basé sur deux rasters préexistants. Le RasterToNumPyArray a un bon exemple que je veux adapter. import arcpy import numpy myArray = arcpy.RasterToNumPyArray('C:/data/inRaster') myArraySum = myArray.sum(1) myArraySum.shape =...

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...

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...

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'extraire des valeurs de rupture naturelles à partir d'un raster à l'aide de Python. Le processus prévu est le suivant: objet raster Arcpy vers NumPy Array (à l'aide de RasterToNumPyArray ) vers des valeurs de rupture naturelles (à l'aide de la fonction PySAL Natural Breaks ). import...

J'essaie d'apprendre les ficelles du traitement d'image de télédétection en utilisant les liaisons Python GDAL et numpy. Dans un premier temps, je lis un fichier géotiff Landsat8, je fais une manipulation simple et j'écris le résultat dans un nouveau fichier. Le code ci-dessous semble fonctionner...

Comment puis-je obtenir les coordonnées projetées ainsi que les valeurs réelles des pixels à ces coordonnées à partir d'un fichier GeoTiff, puis les enregistrer dans un tableau numpy? J'ai le fichier arsenci020l.tif, et ses coordonnées sont en mètres. Voici la sortie abrégée de gdalinfo que j'ai...

Je voudrais reclasser un fichier raster d'un raster avec 10 classes en un raster avec 8 classes utilisant pyhton, gdal et / ou numpy. Les classes sont représentées sous forme d'entiers. J'ai essayé de suivre les étapes de ce post Reclassifier les rasters en utilisant GDAL et Python , le doc...

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 ....

Quelqu'un a-t-il réussi à faire fonctionner les pandas avec l'installation Python de QGIS, sous Windows? J'ai utilisé le programme d'installation OSGeo4W pour installer QGIS2.2, fourni avec Python 2.7 et Numpy 1.5.1. Toute version raisonnablement nouvelle de pandas nécessite Numpy 1.6.1 ou une...

Aujourd'hui, j'ai mis à niveau vers ArcGIS Desktop 10.2.2 ainsi que l'outil de traitement en arrière-plan 64 bits et j'ai remarqué que numpy et matplotlib étaient installés avec Python. En fait, dans la section " Quoi de neuf? " Pour 10.2.1, il indique: ArcGIS 10.2.1 a été mis à niveau pour inclure...

J'ai un problème avec la bibliothèque GDAL lorsque j'ouvre un gros fichier GeoTiff, avec une taille d'environ 32000x32000. Je ne peux pas utiliser la ReadAsArrayfonction car la taille maximale du tableau numpy en python. Mais je me demande pourquoi QGIS peut ouvrir ce fichier facilement. Quelle est...

Le problème consiste à modéliser la propagation d'un signal (par exemple lumière ou son, etc.) à travers une série d'obstacles, comme sur la figure ci-dessous. Le signal ne peut pas traverser la surface inférieure (terrain), mais il peut traverser des obstacles. Je veux compter le nombre...