Ce site est optimisé pour être consulté depuis un navigateur moderne dans lequel JavaScript est activé.

opérateurs logiques GNUPLOT

BenjiBoy

Bonjour à tous,

question toute bête, je fais un graphe quotidien d'une station météo, et le numériseur s'est mis à intégrer des valeurs aberrante (-8000°C) rarement, mais régulièrement. Du coup j'avais une commande "autoscale y" qui de fait me pourri mon échelle ...
Je me suis mis en quête de faire zapper au plot les valeurs en dessous de -50 par exemple.

Donc mon code avant :

plot "meteo.csv" using 1:3 title "measured" with lines lc 4 lw 2

Mais depuis que dans les fichiers j'ai des trucs comme ça :
2020-01-29 11:40:00,59632,7.612,"(949145460, 0)",2.497,117.5,1.428,804.0,0.389,1.652,1006.0,82.2,"(949145640, 0)",482.659912109375
2020-01-29 11:50:00,59633,7.836,"(949146300, 0)",1.823,190.4,74.95,679.8,0.537,0.731,900.0,82.2,"(949146240, 0)",407.86474609375
2020-01-29 12:00:00,59634,7.682,"(949146840, 0)",2.429,162.9,56.64,549.5,0.57,1.629,735.6,79.64,"(949147200, 0)",329.7033996582031
2020-01-29 12:10:00,59635,7.816,"(949147380, 0)",1.726,61.54,55.4,558.5,0.435,0.875,699.6,81.5,"(949147260, 0)",335.1001281738281
2020-01-29 12:20:00,59636,-8190.0,"(949148160, 0)",2.497,192.2,22.94,733.4,0.557,1.264,1178.0,78.23,"(949148400, 0)",440.0456237792969
2020-01-29 12:30:00,59637,8.22,"(949149000, 0)",1.935,149.6,22.28,1132.0,0.511,1.12,1433.0,77.72,"(949148520, 0)",679.0952758789062

J'ai voulu rajouter ça

plot "meteo.csv" using 1:($3 <= -50 ? NaN : $3) title "measured" with lines lc 4 lw 2

Mais le plot ne se fait plus, l'opérateur "<=" n'est pas reconnu.

une idée ?

JBF

Et si tu utilises < au lieu de <= ?

Nasman

N'est ce pas la conséquence d'un dépassement de variable, genre donnée codée sur 16 bits et valeur supérieure à 32767 (dont la valeur qui suit 32768 est considérée comme -32768, la suivante comme -32737...) ?

BenjiBoy

à JBF :
marche pas non plus <, >, >=, ... aucun.
Ne fonctionne que "==", et "!=".

à Nasman :
tu parles du pourquoi la valeur -8190 ? Oui très certainement une erreur d'ordre de codage, d'autant que c'est assez rare ... Mais même les numériseurs les plus performants ont des enregistrements ratés ...
Et de toute façon pour que mon code soit robuste je dois intégrer cette notion de valeur max / min.

a+

JBF

D'après la doc de gnuplot 5.2 l'opérateur <= existe bien. La liste des opérateurs logiques reconnus est donnée page 33 et il y a un exemple page 34.

BenjiBoy

Suis bien d'accord, mais est-ce que ça passe dans une commande "using" ?
En tout cas, résultat de mon script :

multiplot> plot "meteo.csv" using 1:( < -50 ? NaN : ) title "measured" with lines lc 4 lw 2
                                      ^
           line 0: invalid expression

JBF

L'exemple donné page 34 de la doc gnuplot 5.2 :

For data in a file, plot the average of the data in columns 2 and 3 against the datum in column 1, but only
if the datum in column 4 is non-negative:
plot ’file’ using 1:( $4<0 ? 1/0 : ($2+$3)/2 )

La différence c'est que toi tu veux modifier la donnée à tracer, alors que dans l'exemple la condition sert juste à filtrer les données.
Est-ce que le problème n'est pas plutôt le NaN ?

BenjiBoy

Il s'agit aussi d'un filtre dans mon cas, je demande à gnuplot de ne pas interpréter toute valeur en dessous de -50.

Par ailleurs de ce que j'en ai lu sur les topics, le fait qu'il y ait "1/0" (opération interdite) à la place de "NaN" produit le même résultat, ce qui satisfait les condition est de fait ignoré par gnuplot.
Par acquis de conscience j'ai quand même essayé de remplacer NaN par 1/0, ça produit le même résultat. Il bloque sur "<". C'est troublant !

BenjiBoy

Ah, un collègue vient probablement de trouver le problème.
Et on ne peut pas le deviner au travers de mes messages.

Le script gnuplot est inséré dans un script bash au moyen des "balises" [gnuplot <<EOF] et [EOF], et il doit considérer "$3" comme une variable, ou en tout cas il ne les considère pas uniquement selon la syntaxe gnuplot, comme elles n'apparaissent pas dans l'erreur.

multiplot> plot "meteo.csv" using 1:( <-50 ? NaN :  ) 
                                      ^
           line 0: invalid expression

Et sur le web j'ai trouvé comment faire comprendre à gnuplot qu'il s'agit d'un variable linux, mais pas comment dire à gnuplot "non, il s'agit bien d'une variable gnuplot".

BenjiBoy

Et voici la bonne commande :

plot "meteo.csv" using 1:( \$3 < -50 ? NaN : \$3 )

Résolu !

kamaris

Si tu veux que le shell n'interprète rien de ce qui se trouve dans un here doc, tu mets la borne ouvrante entre guillemets :

gnuplot <<'EOF'

gnuplot <<"EOF"

BenjiBoy

Ah excellent tips aussi, merci ;-)

melixgaro

Ah oui, merci Kamaris de cette info.