Apparition de valeurs aberrantes lors de manipulation de masked array numpy

Bonjour,

Dans le cadre de mon stage, je travaille sur des données climatiques. Je travaille donc sur des fichiers netCDF, et manipule des masked arrays.
En particulier, j'ai récupéré les données de température E-OBS de 1950 à 2020 (mesures journalières de température dans toute l'Europe).

Par la suite, j'ai réalisé la climatologie de cette période (moyenne des 1er janviers, moyenne des 2 janvier, etc. afin de pouvoir tracer des anomalies de température), et ai enregistré ces moyennes dans un fichier netCDF dont voici l'en-tête :

Pièce jointe 596717

Des tracés de cycles saisonniers montrent que ces données sont valides.
Cependant, j'ai par la suite cherché à lisser ces moyennes poiint par point. Mais la manipulation de ce tableau 3D fait apparaître des valeurs aberrantes dans un nombre non négligeable de points. Pour lisser mes valeurs, j'utilise une moyenne glissante sur 21 jours, et j'ai concaténé mon tableau de données trois fois en ne travaillant que sur celui du milieu pour ne pas avoiir de problèmes aux bords (fin décembre et début janvier).

Voici le bout de code qui suffit à faire apparaître ces valeurs (lors de la création de la variable extended_temp). Les "print" m'ont permis de vérifier que des valeurs de 9*10³⁶ apparaissent lors de cette manipulation. J'ai aussi vérifié que le mask était bien conservé, ce qui est le cas, et j'ai cherché à le forcer (for i in range(1098)...) mais cela n'a rien changé.

Pièce jointe 596718

Par ailleurs, dans un autre script dans lequel je cherche à récupérer seulement les données d'hiver, et à réordonner sous la forme Décembre-Janvier-Février, le même problème intervient (lors de la création de la variable T_mean_ano).

Pièce jointe 596719

Auriez-vous une idée de ce qui peut causer ceci ?
Je m'arrache les cheveux depuis plusieurs jours...

J'espère n'avoir pas été trop confus !

Merci d'avance,

Bonjour

En l'état on ne peut rien vous dire :
1) Il faut poster votre code entre balise code (le bouton #). Sans quoi on ne peut absolument rien tester sans se taper la recopie de votre code. On veut bien aider mais de là à recopier du code .... Donc bannissez les captures d'écran
2) Pour tester il faudrait qu'on ait les données sur lesquelles vous vous appuyez. Donc soit le fichier d'input, mais le mieux serait que vous définissiez un petit jeu de test à la main (en rentrant direct le tableau dans le code), pour illustrer votre souci sur un problème réduit. Car pareil, votre problème porte sur un point précis de vos données, et personne n'ira se farcir la lecture/compréhension de la manière dont sont structurées vos données et ou il faut regarder précisément.

Que vous ayez à poster un message sur un forum ou non, tout résolution de bug commence par là : définir un encadrement à l'intérieur duquel on sait qu'il y a le problème, et le resserer jusqu'à parvenir à qqch de plus simple possible. C'est comme ca que tout débogage commence.