Discussion :

[GlmGil]

Bonsoir, Je travaille avec des matrices 2D contenant, pour un certain nombre de particule: les latitudes et longitudes de chaque particule en fonction des pas de temps :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
 
    Tabtime = np.zeros(0)
    Tabpart = np.zeros(0)
    Tabtime = np.arange(0, 87,1)
    Tabpart = np.arange(0, 15,1)
    for i in Tabtime:
    for j in Tabpart:     
        if l[i][j] <= 120 and l[i][j] >= 26:
             if L[i][j] >= 15 and L[i][j] <= 45:   
                lt[i][j] = np.append(lt, l[i][j])
                Lt[i][j] = np.append(Lt, L[i][j])

Tabtime contient les 87 pas de temps, et Tabpart contient les particules, ce que je cherche à faire est de garder les coordonnées des particules qui se situent uniquement entre 26° et 120° de longitude et 15° à 45° de Latitude.
J'essaye de ranger ces coordonnées selectionnés dans deux nouvelles matrices mais avec mon code je n'arrive pas à les remplir :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
 
 
print(lt)                
 
>> 
[[0. 0. 0. ... 0. 0. 0.]
 [0. 0. 0. ... 0. 0. 0.]
 [0. 0. 0. ... 0. 0. 0.]
 ...
 [0. 0. 0. ... 0. 0. 0.]
 [0. 0. 0. ... 0. 0. 0.]
 [0. 0. 0. ... 0. 0. 0.]]

Si vous auriez une idée pour cela je veux bien

Merci!

[charliemtx]

Bonjour,

Je ne suis pas sur de comprendre votre code ; les variables l,lt,L,Lt ne sont jamais initialisées et pourtant vous utilisez leur valeur ? D'où viennent ces valeurs ? Un fichier excel ?

Si je comprends bien, à chaque pas de temps, vous avez une matrice qui vous donne la latitude et la longitude de chaque particule, est-ce que ce genre de chose vous convient ?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
import pandas as pd
 
# Create fake dataframe
df = pd.DataFrame(columns=['particle','latitude','longitude'])
df.set_index('particle', inplace=True) # Set particle as index
 
# Create fake data
import random
for i in range(100) : df.loc[str(i)]=[random.random(),random.random()]
 
# Extract data, filter over latitude AND longitude
filtered_data = df[(df['latitude']>0.5) & (df['longitude']<0.5)]
print(filtered_data)

[lg_53]

Et sans utiliser pandas, la fonction np.where devrait suffire à répondre à votre besoin !

[GlmGil]

Bonjour,
merci pour la réponse,

Toutes ces variables ont été initialisées plus tôt dans le code, voici le code depuis le début :

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
32
33
34
35
36
37
38
 
 
 
import netCDF4 as ndf
import matplotlib.pyplot as plt
import numpy as np
import glob
import math
 
 
 
ext="*.nc4"
for files in sorted(glob.glob(ext),reverse=True):
    print(files)  
    part_ds = ndf.Dataset(files, "r", format="NETCDF4")  
    part = part_ds.variables['part']    
    t = part_ds.variables['time']      
    l = part_ds.variables['lon']  
    L = part_ds.variables['lat']
    O = part_ds.variables['O3']
    T = part_ds.variables['T']
    PV = part_ds.variables['PV']
    p = part_ds.variables['p']   
    Tabtime = np.zeros(0)
    Tabpart = np.zeros(0)
    Tabtime = np.arange(0, len(t),1)
    Tabpart = np.arange(0, 15,1)
    lt=np.zeros((len(t),15))
    Lt=np.zeros((len(t),15))
 
    for i in Tabtime:
        for j in Tabpart:     
            if l[i][j] <= 120 and l[i][j] >= 26:
                 if L[i][j] >= 15 and L[i][j] <= 45:   
                    lt[i][j] = np.append(lt, l[i][j])
                    Lt[i][j] = np.append(Lt, L[i][j])
 
print(lt)

Les données sont importées depuis un ensemble de fichier NetCDF, dans lesquels elles sont sont la forme time(particule, longitude, latitude).

Donc vous avez bien compris les données que je traite. Par contre moi je ne sais pas exactement comment utiliser votre code dans mon cas.

D'ailleurs je cherche plus globalement à récupérer toutes les valeurs (O , T, PV, p, L, l) uniquement dans le cas des particules qui sont à : 15°<L<45° et 26°<l<120°.

[wiztricks]

Salut,

Déjà votre code n'affiche que les données du dernier fichier puisque les variables sont réinitialisées à chaque itération.
Ensuite, vérifiez que les conditions

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
            if l[i][j] <= 120 and l[i][j] >= 26:
                 if L[i][j] >= 15 and L[i][j] <= 45:

sont vraies de temps en temps car arrivé à:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

                    lt[i][j] = np.append(lt, l[i][j])

çà devrait se vautrer.

Pour le reste vous avez déjà un problème de cohérence "à la base".

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
    lt=np.zeros((len(t),15))
    Lt=np.zeros((len(t),15))

crée des tableaux avec de tas de zéros.... et la condition ne sera jamais vraie.

Après on peut utiliser des fonctionnalités un peu plus avancées mais pour l'instant, çà ne marche pas parce que vous n'êtes pas assez attentif à ce que vous avez écrit.

- W