erreur lors d'ectraction des valeurs du pixel , lat, long d'un point dans un fichier netCDF4

Version imprimable

bonjour;

je suis nouvelle sur les fichier netCDF4 j'ai une application avec wxpython ou je peux lire les fichiers netCDF4 et les visualiser et de sélectionner des points x y des lignes et colonnes : x = col & y = row : x_coord = 89.9877 , y_coord = 35.8141 ( cordonnées sous curseur )

et les affichers sur une autre fenetres

ce que j'aimerai faire ces d'avoir des informations sur la long , lat & valeur numérique du pixel dans ce point

j'ai un code que je comprend pas comment fonctionne

voilà le code:
Code:

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 32 33 34 35 36 import numpy as np import netCDF4 from netCDF4 import Dataset fic='air.departure.sig995.2012.nc' ncdir='/home/data/' nc_data = netCDF4.Dataset(path+fic,'r') varnames = nc_data.variables.keys() x_coord = 89.9877 y_coord = 35.8141 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 ncdir - Required : The directory of the netcdf file. varnames - Required : The netcdf variables """ with Dataset(ncdir, "r") as nc: # Get the nc lat and lon from the point's x, y lon = nc.variables["lon"][int(round(x_coord))] lat = nc.variables["lat"][int(round(y_coord))] # Return a np.array with the netcdf data air_dep = np.ma.getdata( [nc.variables[air_dep][:, x_coord, y_coord] for air_dep in varnames] ) return air_dep, lon, lat
quand je l’exécute y'a pas d'erreur mais quand d exécuter cette fonction:

Code:

nc_data, lon, lat = ExtractVarsFromNetcdf(89.9877, 35.8141, path+fic, [ air_dep])

je reçois un erreur que air_dep , lat, mon sont pas définis

j'aimerai comprendre comment marche ce code svp et ou se trouve l'erreur

merci d'avance :)

merci

08/07/2019, 11h08
nekcorp

Bonjour,

Citation:

j'ai un code que je comprend pas comment fonctionne

D'où vient ce code ? Peut être voir avec la personne qui a écrit ce code pour lui poser des questions.

Citation:

je reçois un erreur que air_dep , lat, mon sont pas définis

Pouvez vous mettre le message d'erreur en entier, car dans votre code air_dep intervient à différents endroits, ça permettra de mieux comprendre pourquoi vous avez une erreur sur ce paramètre.

Code:

nc_data, lon, lat = ExtractVarsFromNetcdf(89.9877, 35.8141, path+fic, [ air_dep])

Pourquoi mettez vous des crochets autour de air_dep ?

Bonjour,

bon la personne je pense qu'elle a copie de quelque part puisque il arrive pas a m'expliquer comment ca fonctionne

voilà l'erreur :

Code:

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
32
33
In [19]: 
    ...: 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
    ...:         ncdir      - Required : The directory of the netcdf file.
    ...:         varnames   - Required : The netcdf variables
    ...:     """
    ...: 
    ...:     with Dataset(ncdir, "r") as nc:
    ...: 
    ...:         # Get the nc lat and lon from the point's x, y
    ...:         lon = nc.variables["lon"][int(round(x_coord))]
    ...:         lat = nc.variables["lat"][int(round(y_coord))]
    ...: 
    ...:         # Return a np.array with the netcdf data
    ...:         air_dep = np.ma.getdata(
    ...:             [nc.variables[air_dep][:, x_coord, y_coord] for air_dep in 
    ...: varnames]
    ...:         )
    ...: 
    ...:         return air_dep, lon, lat
    ...:     
 
In [20]: nc, lon, lat = ExtractVarsFromNetcdf(89.9877, 35.8141, ncdir, [ air_dep
    ...: ])
---------------------------------------------------------------------------
NameError                                 Traceback (most recent call last)
<ipython-input-20-e3b0979e7612> in <module>()
----> 1 nc, lon, lat = ExtractVarsFromNetcdf(89.9877, 35.8141, ncdir, [ air_dep])
 
NameError: name 'air_dep' is not defined

voilà merci

08/07/2019, 11h54
nekcorp
Code:

1 2 nc, lon, lat = ExtractVarsFromNetcdf(89.9877, 35.8141, ncdir, [ air_dep ...: ])
Lorsque vous appelez air_dep dans la fonction il doit être déclaré préalablement, comme vous avez déclaré x_coord, y_coord, path.

Bonjour,

je viens de comprendre comment ça fonctionne merci :

voilà :

Code:

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
32
33
 
import numpy as np
 
from netCDF4 import Dataset
 
 
x_coord, y_coord = 12.6109, 23.9767
 
path = 'D:/data/netCDF4/air.departure.sig995.2012.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, lon, lat
 
 
air_dep, lon, lat = ExtractVarsFromNetcdf(x_coord, y_coord, path, varnames = ['air_dep'])
 
print(air_dep, lon, lat)

Merci à vous