Interpolation de valeurs dans un fichier

Bonjour,
J'ai un fichier contenant deux colonnes :

12.050 .0000001617893
12.149 .0000001560011
12.248 .0000001521674
12.348 0
12.449 0
...

Et j'ai des valeurs de références qui commencent à 12.020, et avancent de 0.150 :

12.020
12.170
12.320
12.470
...


Tout en gardant les valeurs de la colonne 2, j'aimerais juste que les valeurs de la colonne 1 du fichier soient remplacées par les valeurs de référence les plus proche :

12.050 sera remplacé par 12.020
12.149 par 12.170
12.248 par 12.320
12.348 par 12.320
12.449 par 12.470
...


Sauriez-vous comment faire ?
J'ai essayé avec cette boucle mais ce n'est pas précis parce que les valeurs qui ne sont pas dans le pas de 0.150 sont oubliées :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16

for ((i=1202; i <= 2207; i += 15))   do alt=$(echo "scale=3; ${i}/100" | bc) # création valeur de référence echo $alt >> tab_$i # ecriture dans un fichier fic2=tab_$i awk -vp=10 'NR==FNR{arr[int($1*p)/p]++;next} (int($1*p)/p) in arr' {$fic2,$fichier} >> AvAMA_$i # comparaison valeur de réf/valeur du fichier col2=$(awk '{ print $2 }' AvAMA_$i) # selection colonne 2 à garder   echo $alt $col2 >> fichier_final   rm tab_$i AvAMA_$i     done

---

Et créer des fichiers pour des opérations intermédiaires n'est pas très malin j'imagine...

Moi, je le ferais bêtement avec awk, un truc du genre:

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17

awk_script=' BEGIN { ref0 = 12.020 pas = 0.150 pas2 = pas/2 } { col1 = $1 col2 = $2   k = int((col1 + pas2 - ref0) / pas) ncol1 = ref0 + k * pas printf "%.3f ", ncol1   if ( col2 == 0 ) { printf "0\n"} else { printf "%.13f\n", col2} } '

---

Code:

---

1 2 3 4 5 6

awk "${awk_script}" "plop" 12.020 0.0000001617893 12.170 0.0000001560011 12.320 0.0000001521674 12.320 0 12.470 0

---

Rq: je connais int qui tronque, mais pas de fonction qui arrondisse, d'où l'ajout de pas/2.
Rq: je ne maîtrise pas très bien les formats de printf...

Cette méthode marche parfaitement!

Je n'avais jamais pensé aux scripts awk, ça m'évite un tas d'opérations intermédiaires

Merci à toi!

on notera bien toutefois, que, dans awk "${awk_script}" "plop", les accolades sont totalement superflues.

Citation:

---

Envoyé par N_BaH

on notera bien toutefois, que, dans awk "${awk_script}" "plop", les accolades sont totalement superflues.

---

Tout à fait exact !
C'est une vieille habitude (sur laquelle nous ne reviendrons pas) et qui m'aide pour la lecture et me libère, pour les rares (je l'accorde) cas où elles sont nécessaires, de la charge mentale de vérifier leur nécessité sans risquer de tomber sur un bug pas toujours évident à cerner au premier coup d'oeil:

Code:

---

1 2 3 4 5 6 7

a=A; b=B; a_b=ab   echo "ex1: $a_$b $a_b ${a}_${b}" ex1: B ab A_B   echo "ex2: $a-$b $a-b ${a}-${b}" ex2: A-B A-b A-B

---

Avec les accolades, j'obtiens ce que je veux dans tous les cas, même si elles ne sont nécessaires que dans l'exemple 1.
C'est pour ça que j'ai pris l'habitude de les utiliser comme séparateur syntaxique pour ne pas me poser la question de leur nécessité.

C'est aussi pour limiter les erreurs de frappe que j'utilise souvent l'option set -u en début de script.

Citation:

---

Envoyé par GlmGil

Je n'avais jamais pensé aux scripts awk, ça m'évite un tas d'opérations intermédiaires

---

Lorsque j'ai à traiter un fichier avec des colonnes et avec un peu de traitement algorithmique, j'aime bien utiliser un "vrai" langage de programmation, comme "awk", qui est plus puissant que certains ne l'imaginent quand ils voient de simples awk '{print $2}' par exemple.

PS: Je ne nie pas que "bash" est aussi un "vrai" langage de programmation...

Citation:

---

[...]Avec les accolades, j'obtiens ce que je veux dans tous les cas, même si elles ne sont nécessaires que dans l'exemple 1.
C'est pour ça que j'ai pris l'habitude de les utiliser comme séparateur syntaxique pour ne pas me poser la question de leur nécessité.

C'est aussi pour limiter les erreurs de frappe que j'utilise souvent l'option set -u en début de script.

---

tout ça n'inspire pas confiance. Ça ne met pas en valeur tes compétences. :)

Citation:

---

Envoyé par jack-ft

PS: Je ne nie pas que "bash" est aussi un "vrai" langage de programmation...

---

Sauf qu'en bash pure, on ne peut pas tout faire :aie:

Citation:

---

Envoyé par jack-ft

Moi, je le ferais bêtement avec awk, un truc du genre:

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17

awk_script=' BEGIN { ref0 = 12.020 pas = 0.150 pas2 = pas/2 } { col1 = $1 col2 = $2   k = int((col1 + pas2 - ref0) / pas) ncol1 = ref0 + k * pas printf "%.3f ", ncol1   if ( col2 == 0 ) { printf "0\n"} else { printf "%.13f\n", col2} } '

---

---

Bonjour,

je me permets de rebondir dans cette discussion pour une petite question sur le même sujet.
Je vois que ce script marche très bien mais j'aurais besoin de savoir combien de valeur de mon fichier initial sont 'interpolable' par chaque valeur de référence.

Par exemple :

j'ai ce fichier ci :

Code:

---

1 2 3 4 5 6 7 8 9

12.03100 12.05900 12.49500 12.49800 12.50000 12.50500 12.78700 12.81200

---

et j'aimerais savoir combien de valeur exactement sont 'interpolable' par la valeur ref 12.020 ?


En fait mon but est de faire ces comparaisons sur plusieurs fichiers et d'interpoler aux valeurs ref le fichier où il y a le plus grand nombre de valeurs interpolables

C'est peut-être pas très clair dit comme ça mais voilà un exemple :

premier fichier :

Code:

---

1 2 3 4 5

12.03100 12.05900 12.49500 14.00599

---

deuxième fichier :

Code:

---

1 2 3 4 5 6 7 8

12.02100 12.08900 12.23564 12.45500 12.55555 13.56600 14.00245

---

Donc ce sera le fichier 2 auquel je vais interpoler ma valeur ref 12.020.

Vous sauriez comment faire ça ?

Merci pour votre aide!

Bonjour :coucou:

Compte-les, au lieu de les afficher.

Ensuite, tu peux utiliser la variable intégrée FILENAME comme clé de tableau associatif, pour écrire les fichiers les uns à la suite des autres dans la ligne de commande.

Code:

---

awk '... nb[FILENAME]++; ...' fichier1.txt fichier.txt fichier2.txt fichier3.txt fichier4.txt fichiern.txt

---

Il ne reste plus qu'à trouver le max du tableau nb[], et son indice associé.

Bonjour !

Je ne suis pas sûr d'avoir compris, en gros il faudrait ranger les valeurs interpolées dans un fichier temporaire pour chaque valeur 'ref' , et de compter combien il y en a dans chaque fichier c'est ça ?

Pas sûr de bien savoir m'y prendre, je ne comprends pas vraiment comment fonctionne le script awk :lol: