Bonjour à tous
J'ai deux vecteur de et qui correspondent aux longitudes et latitudes de mes profils.
Je cherche a diviser toute ma zone en petit carreaux et chercher combien j'ai des profiles dans chaque petit carreaux. J'ai fait alors
1 2 3 4
| vec_x_14km=min(X):0.126:max(X);*%vecteur longitude de ma zone
vec_y_14km=min(Y):0.126:max(Y); %vecteur latitude de ma zone
[xx_vec,yy_vec] =meshgrid(vec_x_14km,vec_y_14km);
xx_vec=xx_vec';yy_vec=yy_vec'; |
et j'ai initialisé mon compteur comme suit
compteur_anti = 0*ones(length(vec_x_14km)-1,length(vec_y_14km)-1);
ensuite je fais une boucle sur i qui parcourt
i=1:(length(vec_x_14km)-1)
une boucle sur j qui parcourt
j=1:(length(vec_y_14km)-1)
et encore je fait une troisième boucle sur k pour parcourir mes vecteur x_lon et x_lat (1,7600)
je fais ça car je prends une première ligne et je parcourt toutes mes colonnes et pour chaque petit carreau je parcourt aussi les vecteur profils pour tester si x_lon(k) soit compris entre vec_x_14km(i) et vec_x_14km(i+1) et aussi mon y_lat(k) compris entre vec_y_14km(j) et vec_y_14km(j+1)
mais ce programme est très très long et je me demande si y a pas une méthode plus simple que la mienne
merci beaucoup
Partager