Discussion :

[soma12]

Bonjour,

j'ai un fichier netcdf4 que j'aimerai afficher tous les valeurs des pixel ou des valeurs du variable air_dep sous forme d'un tableau que chaque ligne-colonne correspond à une valeurs du variable :



cela juste un exemple , je pense que je devrai utiliser quelque chose comme Var_air[x-5 : x+5,y-5:y+5]

mais je sais pas trop comment le faire les x y sont les coordonnées ligne colonne

merci d'avance

[lg_53]

Bonjour

Comment ferait tu déjà pour afficher la valeur d'un seul pixel ? Il faut toujours raisonner sur des problèmes simples d'abord.

Montre du code.

[soma12]

Bonjour ,

voilà le code soi par cette méthode sois par l'affichage avec imshow on peut voir la valeur du variable dans un pixel de coordonnées (x, y) choisi :

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
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
import numpy as np
from netCDF4 import Dataset
 
 
x_coord, y_coord = 12.6109, 23.9767
 
path = '/data/air2012.nc'
 
 
def ExtractVarsFromNetcdf(x_coord, y_coord, path, varnames):
    """
 
    @params:
       x_coord - Required : the x coordinate of the point
       x_coord - Required : the y coordinate of the point
       path+fic - Required : The directory of the netcdf file.
       varnames - Required : The netcdf variables
    """
    with Dataset(path, "r") as nc:
        # Get the nc lat and lon from the point's x, y
        #lon = np.ma.getdata(nc.variables["lon"][int(round(x_coord))])
        #lat = np.ma.getdata(nc.variables["lat"][int(round(y_coord))])
        # Return a np.array with the netcdf data
        air_dep = np.ma.getdata(nc.variables['air_dep'][0,x_coord , y_coord])
 
        return air_dep
 
 
air_dep = ExtractVarsFromNetcdf(x_coord, y_coord, path, varnames = ['air_dep'])
 
print((int(air_dep))

ca donne :

-38


voilà mais ce que je cherche c'est comme une matrice 5x5 mais sous forme d'un tableau !

merci

[wiztricks]

Salut,

voilà mais ce que je cherche c'est comme une matrice 5x5 mais sous forme d'un tableau !

Il doit y avoir des façons plus intelligente d'utiliser l'API de netCDF4 que de récupérer une valeur après l'autre mais il y a certainement des forums d'utilisateurs ou des exemples qui traînent qui pourraient être instructifs.

A défaut, côté Python, votre question se résume à ranger 25 valeurs dans une liste des listes 5x5. Ce qui est une difficulté que tout débutant apprend à résoudre...
Dit autrement, si vous avez pris le temps de débuter, quelle difficulté rencontrez vous pour faire une chose qui paraît si simple?

- W

[soma12]

les difficultés c'est au niveau du fichier netcdf4 c'est pas simple a codé avec ce type de fichier surtout que pour mes connaissances j'ai pas trop des connaissances de ce type des fichiers

[wiztricks]

Salut,

j'ai pas trop des connaissances de ce type des fichiers

C'est une base de données un peu spéciale.
Et comme toute base de données, il vaut mieux savoir l'utiliser avant de vouloir programmer. Ce serait une base de données "populaire", je vous suggèrerais d'aller faire un tour dans la rubrique SGDB...
Mais dans ce cas, je ne peux que modifier le titre de la discussion en y mettant NetCDF4 histoire que si un programmeur Python qui connaît passe dans le coin...
Ceci dit, il y a certainement des tutos. et des articles à consulter sur Internet.

- W