Picard

bonjour @meteo-melin, j'ai réussi à décoder correctement les bufr grâce à la librairie @theperk (un grand merci). Maintenant, comme toi, je me confronte à un pb de localisation des données bufr (IPRN21 ou IMFR27). - la loc dans les IPRN20 en format h5 est correct, dans le sens où j'arrive à générer une grille lon/lat qui matche météociel et que je valide par ailleurs (colloc avec mon pluvio) - la loc des bufr est décalée (*) , quelque soit la solution que j'emploie (et j'ai testé un grand nb de solutions ces derniers jours) j'en suis là (voir le code ci-dessous) mais cette solution est décalée (*) par rapport à IPRN20 (je cherche une précision au km) (*) les cellules de pluie sont décalées si quelqu'un.e à une idée, je suis preneur. Je continue de chercher de mon coté. PS: Est-ce qu'il pourrait y avoir un pb de parallaxe entre IPRN20 et IPRN21 ? datas_messages = lib_bufr.deco_bufr(fname) dm = datas_messages[0] # ~ print('\n'.join(sorted(list(datas_messages[0].keys())))) rows = int(dm["Number of pixels per column"][0]) cols = int(dm["Number of pixels per row"][0]) dx = float(dm["Pixel size on horizontal - 1"][0]) dy = float(dm["Pixel size on horizontal - 2"][0]) lat_ref = float(dm["Latitude de reference"][0]) lon0 = float(dm["Longitude du meridien parallele a l'axe des Y"][0]) ptype = int(dm["Projection type"][0]) # Projection OPERA : LCC (ptype=1) ou Stéréographique (ptype=3) if ptype == 1: projstr = f"+proj=lcc +lat_1=44 +lat_2=49 +lat_0={lat_ref} +lon_0={lon0} +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs" elif ptype == 3: projstr = f"+proj=stere +lat_0=90 +lon_0={lon0} +lat_ts={lat_ref} +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs" else: raise ValueError(f"Projection type non géré: {ptype}") crs_geo = CRS.from_epsg(4326) crs_proj = CRS.from_proj4(projstr) fwd = Transformer.from_crs(crs_geo, crs_proj, always_xy=True) inv = Transformer.from_crs(crs_proj, crs_geo, always_xy=True) lats = np.asarray(dm["Latitude (high accuracy)"], float) lons = np.asarray(dm["Longitude (high accuracy)"], float) lons = np.where(lons > 180, lons - 360, lons) # Coin UL = (lat max, lon min) parmi les tie-points i_ul = np.argmin((lats - lats.max())**2 + (lons - lons.min())**2) x_ul, y_ul = fwd.transform(lons[i_ul], lats[i_ul]) # ~ DECAL_X = +85 # ~ DECAL_Y = +35 DECAL_X = 0 DECAL_Y = 0 # Grille aux centres de pixels xs = x_ul + (np.arange(cols) + DECAL_X*0.5) * dx ys = y_ul - (np.arange(rows) + DECAL_Y*0.5) * dy # vers le sud XX, YY = np.meshgrid(xs, ys) lon, lat = inv.transform(XX, YY) lon = np.where(lon > 180, lon - 360, lon)

est-ce qu'il possible de modifier le titre de ce topic en [RESOLU] ?

nooooooon mais quel boulet je mettais l'adresse de mon routeur, erreur de newbie ... désolé... maintenant, ça marche parfaitement, merci pour le support @libertykite d'avoir répondu, ça m'a obligé à regarder tout ça de plus prés.

Bonjour, je n'arrive pas à faire fonctionner le script d'exemple de récupération des données via l'api open data. J'ai bien généré une ip via la page dédiée. J'ai même généré 3 clés différentes: 2 pour les 2 ip différentes selon la machine depuis laquelle j'éxecute le code, et 1 clé redondante pour une des 2 machines j’exécute le code suivant mais quelque soit la machine et la clé, j'ai l'erreur suivante et effectivement, l'objet response contient Question subsidiaire: je n'arrive pas à supprimer les clés via la page dédiée... je clique sur la poubelle, un pop-up d'avertissement apparaît (pour confirmer) et pourtant les 3 clés sont toujours là.... J'ai vu que certains d'entre vous ont réussi, je dois me planter qq part ^ ^. une idée: @Marc09, @Fred59_, @libertykite ?

Merci beaucoup, @Fred59_

merci @Fred59_! je suis bien conscient du contexte associatif, étant moi-même bien impliqué dans un tiers-lieu où j'ai mon bureau 😁 https://leprevert81.fr/ et j'ai bien l'intention de faire en sorte que ma boîte adhère à l'asso (on est 1 dans ma boîte, ça devrait se faire facilement) @+ B.

@Fred59_, un petit up sur ce sujet ? merci, Bruno

Bonjour à tou.te.s, Je m'appelle Bruno Picard et je viens tout juste de connaître ce site et vos activités: tout est très enthousiasmant. En deux mots de présentation, je travaille dans l'observation de l'atmosphère par radiomètre micro-ondes, en collaboration avec des scientifiques, des expert.e.s du climat et les agences spatiales (CNES, ESA, Eumetsat). Par ailleurs, j'ai vu que vous aviez une activité autour de l'aspect climatique des séries temporelles, j'espère pouvoir y contribuer. bonne journée à tou.te.s Bruno PS: quelques exemples de mon activité ici https://www.linkedin.com/in/bruno-picard-874533155/ ou là https://www.linkedin.com/company/fluctus-sas/

Merci pour ta réponse rapide @Fred59_! Est-ce qu'il y aurait une liste offline par hasard ? Ou un fichier kmz ?

Bonjour à tou.te.s, Je m'appelle Bruno Picard et je viens tout juste de connaître ce site et vos activités: tout est très enthousiasmant. En deux mots de présentation, je travaille dans l'observation de l'atmosphère par radiomètre micro-ondes, en collaboration avec des scientifiques, des experts du climat et les agences spatiales (CNES, ESA, Eumetsat). Je souhaiterais utiliser votre réseau pour la validation des observations de la vapeur d'eau par satellite et leurs comparaisons aux mesures des stations GPS: les événements de pluie perturbent ces derniers. J'ai commencé à coder qq interfaces avec votre base open data. Ma question est la suivante: existe-t-il une méthode pour récupérer automatique la liste des codes des stations ? cela me permettrait de faciliter la récupération automatique des données puisque la requête demande un stations[]=CODE. Par ailleurs, j'ai vu que vous aviez une activité autour de l'aspect climatique des séries temporelles, j'espère pouvoir y contribuer. bonne journée à tou.te.s Bruno PS: quelques exemples de mon activité ici https://www.linkedin.com/in/bruno-picard-874533155/ ou là https://www.linkedin.com/company/fluctus-sas/