IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

Calcul scientifique Python Discussion :

Extract all pixels values from GeoTIFF with Python


Sujet :

Calcul scientifique Python

  1. #1
    Futur Membre du Club
    Homme Profil pro
    Étudiant
    Inscrit en
    Juillet 2019
    Messages
    5
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France

    Informations professionnelles :
    Activité : Étudiant
    Secteur : High Tech - Multimédia et Internet

    Informations forums :
    Inscription : Juillet 2019
    Messages : 5
    Points : 6
    Points
    6
    Par défaut Extract all pixels values from GeoTIFF with Python
    bonjour,
    j'ai un petit soucis pour récupérer toutes les valeurs des pixels d'un Geotiff d'une seule bande. D'abord j'ai récupéré les coordonnées x et y et je voudrais les exploiter pour accéder à la valeur du pixel. voici mon code:
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    raster = rasterio.open('.../g.tif')
       br_pixel = (raster.width * raster.height) 
     
       liste3 = []
       i = 0   
       while i <= nbr_pixel:
          (x,y) = raster.xy(i,i)
          listX.append(int(x))
          listY.append(int(y))
       i = i + 1 
       l = list(zip(listX,listY))
       liste3 = raster.sample(l)
       print(liste3)
     
       <generator object sample_gen at 0x7f7ef5821c50>

  2. #2
    Expert éminent

    Homme Profil pro
    Inscrit en
    Octobre 2008
    Messages
    4 298
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : Belgique

    Informations forums :
    Inscription : Octobre 2008
    Messages : 4 298
    Points : 6 778
    Points
    6 778
    Par défaut
    Salut,

    Tu as créé une boucle infinie là. La ligne i = i + 1doit avoir une indentation de plus.

    D'autre part, tu incrémentes x et y en même temps donc, en fait, tu parcours ton tableau en diagonale. (0, 0), (1, 1), (2, 2), (3, 3), etc.

    Ta boucle devrait être comme ceci:
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
     
    for y in range(raster.height):
        for x in range(raster.width):
            pixel = raster.xy(x, y)
            ...
    après, je ne suis pas trop sûre de ce que tu veux faire.

    Tu as vu la doc ?
    https://rasterio.readthedocs.io/en/latest/

  3. #3
    Futur Membre du Club
    Homme Profil pro
    Étudiant
    Inscrit en
    Juillet 2019
    Messages
    5
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France

    Informations professionnelles :
    Activité : Étudiant
    Secteur : High Tech - Multimédia et Internet

    Informations forums :
    Inscription : Juillet 2019
    Messages : 5
    Points : 6
    Points
    6
    Par défaut
    Merci pour l'astuce. la boucle parcourt ligne par ligne en récupérant x et y. Maintenant il me faut la valeur du pixel afin d'obtenir x,y,valpixel.
    Merci encore une fois
    zouhaib

Discussions similaires

  1. Réponses: 5
    Dernier message: 08/12/2009, 13h24
  2. Réponses: 1
    Dernier message: 03/05/2009, 22h43
  3. Réponses: 1
    Dernier message: 19/01/2009, 16h36
  4. Réponses: 1
    Dernier message: 12/11/2008, 20h57
  5. Start with python
    Par elghadi_mohamed dans le forum Général Python
    Réponses: 4
    Dernier message: 21/04/2007, 00h33

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo