TL;DR
Le programme de conversion (travail en cours) est placé sur GitHub, veuillez donc le faire à vos propres risques.
Ceci est une continuation de Convertir le GRIB des coordonnées par polaires du radar de l'Agence météorologique GPV en netCDF (norme CF / Radial). Consultez cet article pour obtenir une description des termes CF / radial.
Données de coordonnées polaires radar de l'Agence météorologique fournies via le Centre de soutien aux entreprises météorologiques, qui a été traité la dernière fois, sont fournies en temps réel. Cependant, il n'inclut pas les informations (paramètre de polarisation) obtenues par la double polarisation du radar. Les données passées du radar à double polarisation de Tokyo, y compris ce paramètre de polarisation, se présentent sous la forme de "données météorologiques fournies pour appréhender les besoins d'utilisation" [Environnement d'utilisation des données météorologiques passées (FY2)]( Il est publié sur AWS dans le cadre de https://www.data.jma.go.jp/developer/past_data/index.html). Cette fois, je vais convertir ces données en NetCDF (CF / Radial Terms).
Données de coordonnées polaires du radar météorologique à double polarisation fournit les 8 variables suivantes:
| Le nom du paramètre | Abréviation | Numéro de paramètre(Réellement) | Numéro de paramètre(spécification) |
|---|---|---|---|
| Intensité de réflexion de polarisation horizontale | ref | 195 | 195 |
| Vitesse Doppler | vel | 228 | 2 |
| Plage de vitesse Doppler | vsw | 230 | 0 |
| Différence de facteur réflexe(Zdr) | zdr | 197 | 197 |
| Différence de phase entre les polarisations du signal reçu(ψdp) | psd | 198 | 198 |
| Taux de changement de différence de phase entre les polarisations(Kdp) | kdp | 202 | 202 |
| Coefficient de corrélation entre polarisations(ρhv) | rhv | 199 | 199 |
| Informations sur le contrôle de la qualité | qci | 215 | 206 |
Les numéros de paramètres dans les spécifications et les numéros de paramètres réellement stockés dans GRIB sont différents, ils sont donc répertoriés ensemble. La spécification indique que les numéros de paramètre 207-254 sont en attente, il semble donc qu'ils ne sont tout simplement pas répertoriés.
En ce qui concerne l'intensité de réflexion, ref (réflectivité de base) est utilisé dans le nom du fichier GRIB, mais il convient de noter que zhh (intensité de réflexion de polarisation horizontale) avec le paramètre numéro 195 est en fait stocké. Ce n'est pas un gros problème, car ils peuvent être considérés comme la même variable dans la pratique.
Je n'expliquerai pas la signification des variables car il s'agit d'une note technique, mais en japonais, vous devriez lire Fukao et Hamazu (2005) "Radio Remote Sensing of Meteorology and Atmosphere" Kyoto University Press. Vous pouvez télécharger le PDF texte intégral à partir du référentiel d'informations universitaires de l'Université de Kyoto.
La séquence de balayage du Tokyo Radar est indiquée dans le tableau ci-dessous.
Cependant, l'angle d'élévation de 90 ° (pointage vertical) du balayage numéro 20 est un balayage pour la correction du radar de polarisation, mais il n'est effectué que pour ceux dont le dernier chiffre du temps de référence est 0, et le dernier. Pas fait si le chiffre est 5. En d'autres termes, une fois tous les deux analyses de volume, toutes les 10 minutes.
Par conséquent, à l'exception du balayage numéro 20, l'espacement des distances est constant à 250 m.
| Numéro de scan | Angle d'élévation(°) | Portée maximale(km) | Intervalle de plage(m) | Espacement des angles d'aspect(°) |
|---|---|---|---|---|
| 1 | 25.0 | 64 | 250 | 0.70 |
| 2 | 17.8 | 64 | 250 | 0.70 |
| 3 | 12.8 | 64 | 250 | 0.70 |
| 4 | 9.3 | 120 | 250 | 0.70 |
| 5 | 6.7 | 120 | 250 | 0.70 |
| 6 | 4.8 | 180 | 250 | 0.70 |
| 7 | 3.5 | 180 | 250 | 0.70 |
| 8 | 2.5 | 400 | 250 | 0.35 |
| 9 | 1.7 | 400 | 250 | 0.35 |
| 10 | 1.1 | 400 | 250 | 0.35 |
| 11 | 0.7 | 400 | 250 | 0.35 |
| 12 | 0.3 | 400 | 250 | 0.35 |
| 13 | 0.0 | 400 | 250 | 0.35 |
| 14 | 0.3 | 250 | 250 | 0.70 |
| 15 | 0.3 | 250 | 250 | 0.35 |
| 16 | 0.7 | 250 | 250 | 0.35 |
| 17 | 0.7 | 250 | 250 | 0.70 |
| 18 | 1.3 | 150 | 250 | 0.70 |
| 19 | 2.2 | 150 | 250 | 0.70 |
| 20 | 90.0 | 64 | 125 | 0.70 |
Données de coordonnées polaires du radar météorologique à double polarisation a été traitée la dernière fois [données de coordonnées polaires du radar de l'Agence météorologique](https :: //www.data.jma.go.jp/add/suishin/shiyou/pdf/no13702) et le format GRIB sont identiques, mais il existe quelques différences:
Les modèles utilisés sont passés de 3.50120 à 3.50121 et 4.51022 à 4.51123. --Par exemple, l'endroit où la définition de l'angle d'élévation est décrite de la section 4 (section de définition du produit) à la section 3 (section de définition du système de treillis) a changé.
La compression des données est passée de la compression de longueur d'exécution (modèles 5.200 et 7.200) à la compression simple (faible taux de compression) (modèles 5.0 et 7.0), en supposant la compression gzip de l'ensemble du GRIB.
--En réponse au fait que l'analyse du volume prend 10 à 5 minutes, les données sont désormais séparées en fichiers séparés toutes les 5 minutes.
--La distance entre la distance et l'azimut était constante, mais il en existe maintenant trois types: (250m, 0.7 °), (250m, 0.35 °) et (125m, 0.35 °) (voir scan pour plus de détails). Voir le tableau de séquence)
Jusqu'à présent, tous les scans (angles d'élévation) étaient stockés dans un fichier pour chaque variable, mais maintenant chaque scan (angle d'élévation) est stocké dans un fichier séparé.
La vitesse Doppler n'était auparavant présente que dans certains des angles d'élévation observés, mais elle est maintenant stockée dans tous les angles d'élévation ainsi que dans l'intensité de réflexion.
La même chose que la dernière fois, Format de données de coordonnées polaires radar météorologique à double polarisation (commun à l'observation PPI / RHI) Et CfRadial Data File Format Version 1.4 et écrivez simplement le code à convertir. Les changements par rapport à la dernière fois sont les suivants:
--Depuis qu'il est devenu une simple compression, la section du document a été directement décodée par Python sans utiliser wgrib2.
Actuellement, certains défis demeurent:
――Déterminez s'il faut stocker les informations de contrôle de qualité exprimées en 8 bits telles quelles ou stocker les variables séparément pour les résultats de la qualité des données du système à polarisation simple, de la qualité des données du système à double polarisation et du MTI sélectif.
Recommended Posts