GoPro & Linux

https://github.com/progweb/gpx2video/commits/main/

Pour la télémétrie, j'ai corrigé hier soir.

Ffprobe sur ta vidéo en input ?

Ok pour la liste des évolutions/corrections

ffprobe de la vidéo :

Code:: ffprobe Jonzac09.mp4
ffprobe version 4.4.2-0ubuntu0.22.04.1 Copyright (c) 2007-2021 the FFmpeg developers
built with gcc 11 (Ubuntu 11.2.0-19ubuntu1)
configuration: --prefix=/usr --extra-version=0ubuntu0.22.04.1 --toolchain=hardened --libdir=/usr/lib/x86_64-linux-gnu --incdir=/usr/include/x86_64-linux-gnu --arch=amd64 --enable-gpl --disable-stripping --enable-gnutls --enable-ladspa --enable-libaom --enable-libass --enable-libbluray --enable-libbs2b --enable-libcaca --enable-libcdio --enable-libcodec2 --enable-libdav1d --enable-libflite --enable-libfontconfig --enable-libfreetype --enable-libfribidi --enable-libgme --enable-libgsm --enable-libjack --enable-libmp3lame --enable-libmysofa --enable-libopenjpeg --enable-libopenmpt --enable-libopus --enable-libpulse --enable-librabbitmq --enable-librubberband --enable-libshine --enable-libsnappy --enable-libsoxr --enable-libspeex --enable-libsrt --enable-libssh --enable-libtheora --enable-libtwolame --enable-libvidstab --enable-libvorbis --enable-libvpx --enable-libwebp --enable-libx265 --enable-libxml2 --enable-libxvid --enable-libzimg --enable-libzmq --enable-libzvbi --enable-lv2 --enable-omx --enable-openal --enable-opencl --enable-opengl --enable-sdl2 --enable-pocketsphinx --enable-librsvg --enable-libmfx --enable-libdc1394 --enable-libdrm --enable-libiec61883 --enable-chromaprint --enable-frei0r --enable-libx264 --enable-shared
WARNING: library configuration mismatch
avcodec configuration: --prefix=/usr --extra-version=0ubuntu0.22.04.1 --toolchain=hardened --libdir=/usr/lib/x86_64-linux-gnu --incdir=/usr/include/x86_64-linux-gnu --arch=amd64 --enable-gpl --disable-stripping --enable-gnutls --enable-ladspa --enable-libaom --enable-libass --enable-libbluray --enable-libbs2b --enable-libcaca --enable-libcdio --enable-libcodec2 --enable-libdav1d --enable-libflite --enable-libfontconfig --enable-libfreetype --enable-libfribidi --enable-libgme --enable-libgsm --enable-libjack --enable-libmp3lame --enable-libmysofa --enable-libopenjpeg --enable-libopenmpt --enable-libopus --enable-libpulse --enable-librabbitmq --enable-librubberband --enable-libshine --enable-libsnappy --enable-libsoxr --enable-libspeex --enable-libsrt --enable-libssh --enable-libtheora --enable-libtwolame --enable-libvidstab --enable-libvorbis --enable-libvpx --enable-libwebp --enable-libx265 --enable-libxml2 --enable-libxvid --enable-libzimg --enable-libzmq --enable-libzvbi --enable-lv2 --enable-omx --enable-openal --enable-opencl --enable-opengl --enable-sdl2 --enable-pocketsphinx --enable-librsvg --enable-libmfx --enable-libdc1394 --enable-libdrm --enable-libiec61883 --enable-chromaprint --enable-frei0r --enable-libx264 --enable-shared --enable-version3 --disable-doc --disable-programs --enable-libaribb24 --enable-libopencore_amrnb --enable-libopencore_amrwb --enable-libtesseract --enable-libvo_amrwbenc --enable-libsmbclient
libavutil 56. 70.100 / 56. 70.100
libavcodec 58.134.100 / 58.134.100
libavformat 58. 76.100 / 58. 76.100
libavdevice 58. 13.100 / 58. 13.100
libavfilter 7.110.100 / 7.110.100
libswscale 5. 9.100 / 5. 9.100
libswresample 3. 9.100 / 3. 9.100
libpostproc 55. 9.100 / 55. 9.100
Input #0, mov,mp4,m4a,3gp,3g2,mj2, from 'Jonzac09.mp4':
Metadata:
major_brand : isom
minor_version : 512
compatible_brands: isomiso2avc1mp41
encoder : Lavf58.76.100
Duration: 00:46:48.94, start: 0.000000, bitrate: 6178 kb/s
Stream #0:0(und): Video: h264 (High 4:2:2) (avc1 / 0x31637661), yuvj422p(pc, bt470bg/unknown/unknown), 1280x720 [SAR 1:1 DAR 16:9], 6044 kb/s, 25 fps, 25 tbr, 12800 tbn, 50 tbc (default)
Metadata:
handler_name : VideoHandler
vendor_id : [0][0][0][0]
Stream #0:1(und): Audio: aac (LC) (mp4a / 0x6134706D), 48000 Hz, stereo, fltp, 127 kb/s (default)
Metadata:
handler_name : SoundHandler
vendor_id : [0][0][0][0]

J'essaye la version d'hier soir et je reviens
Pulsar33

Re,

Je n'ai encodé que 3 minutes mais le problème est toujours là. On dirait que les valeurs calculées par delta semblent 3 fois plus élevées qu'elles ne devraient.

J'ai regardé plus en détail le gpx que j'utilise et qui est généré par le site VisuGPX pour injecter l'élévation dans ma trace qui en est dépourvue, il y a 2 points par seconde

Code:: <time>2024-05-05T10:37:23+02:00</time>
<time>2024-05-05T10:37:23.665000+02:00</time>

Même si ça peut perturber, je ne suis pas certain que ce soit l'explication ?

@+ Puls

Bon, j'ai recodé une vidéo de l'an dernier et il n'y a pas le problème avec ta version d'hier soir.

Il y a parfois seulement plusieurs poits par seconde mais souvent un seul. Voici un extrait du gpx.

Code:: <trkpt lat="45.4449" lon="-0.4469">
   <ele>31.3</ele>
   <time>2023-04-19T16:28:40.053000+0200</time>
   </trkpt>
   <trkpt lat="45.44485" lon="-0.44692">
   <ele>31.3</ele>
   <time>2023-04-19T16:28:40.993000+0200</time>
   </trkpt>
   <trkpt lat="45.44482" lon="-0.44694">
   <ele>31.3</ele>
   <time>2023-04-19T16:28:42.010000+0200</time>
   </trkpt>
   <trkpt lat="45.44477" lon="-0.44697">
   <ele>31.1</ele>
   <time>2023-04-19T16:28:43.010000+0200</time>
   </trkpt>
   <trkpt lat="45.44472" lon="-0.447">
   <ele>30.8</ele>
   <time>2023-04-19T16:28:44.002000+0200</time>
   </trkpt>
   <trkpt lat="45.44468" lon="-0.44702">
   <ele>30.7</ele>
   <time>2023-04-19T16:28:44.997000+0200</time>
   </trkpt>
   <trkpt lat="45.44463" lon="-0.44705">
   <ele>30.7</ele>
   <time>2023-04-19T16:28:45.988000+0200</time>
   </trkpt>
   <trkpt lat="45.44459" lon="-0.44708">
   <ele>30.7</ele>
   <time>2023-04-19T16:28:46.991000+0200</time>
   </trkpt>
   <trkpt lat="45.44454" lon="-0.44711">
   <ele>30.9</ele>
   <time>2023-04-19T16:28:47.976000+0200</time>

A l'inverse, le gpx du cas qui ne marche pas contient systématiquement deux points par seconde :

Code:: <trkpt lat="45.43648" lon="-0.433">
   <ele>37.1</ele>
   <time>2024-05-05T10:37:14.664000+02:00</time>
   </trkpt>
   <trkpt lat="45.43648" lon="-0.433">
   <ele>37.1</ele>
   <time>2024-05-05T10:37:15+02:00</time>
   </trkpt>
   <trkpt lat="45.43648" lon="-0.433">
   <ele>37.1</ele>
   <time>2024-05-05T10:37:15.666000+02:00</time>
   </trkpt>
   <trkpt lat="45.43648" lon="-0.433">
   <ele>37.1</ele>
   <time>2024-05-05T10:37:16+02:00</time>
   </trkpt>
   <trkpt lat="45.43648" lon="-0.433">
   <ele>37.1</ele>
   <time>2024-05-05T10:37:16.665000+02:00</time>
   </trkpt>
   <trkpt lat="45.43648" lon="-0.433">
   <ele>37.1</ele>
   <time>2024-05-05T10:37:17+02:00</time>
   </trkpt>
   <trkpt lat="45.43648" lon="-0.433">
   <ele>37.1</ele>
   <time>2024-05-05T10:37:17.666000+02:00</time>
   </trkpt>
   <trkpt lat="45.43648" lon="-0.433">
   <ele>37.1</ele>
   <time>2024-05-05T10:37:18+02:00</time>
   </trkpt>
   <trkpt lat="45.43648" lon="-0.433">
   <ele>37.1</ele>
   <time>2024-05-05T10:37:18.666000+02:00</time>
   </trkpt>
   <trkpt lat="45.43648" lon="-0.433">
   <ele>37.1</ele>
   <time>2024-05-05T10:37:19+02:00</time>
   </trkpt>

J'espère que ça peut aider
@+ Puls

Pulsar33 a écrit:: A l'inverse, le gpx du cas qui ne marche pas contient systématiquement deux points par seconde :

J'espère que ça peut aider
@+ Puls

Je suis rassuré... je n’ai pas tout cassé Smile

Bon tu as mis le doigt sur le vrai problème.

De mémoire, j’ai tout codé à la seconde... du coup entre 2 points, tu bouges et la durée est nulle... donc ça fait une vitesse proche de l’infinie Wink

Je vais reprendre tout cela.

Merci pour ce temps passer à tester.

Avec plaisir !
Bon courage pour la correction
Bonne soirée
Pulsar33

En attendant de faire mieux (parce que ça va me prendre un peu de temps), j'ai fait un fix pour gérer ton cas.

Ca devrait être mieux.

Bonjour et merci
Cependant, c'est une journée de voyage pour moi aujourd'hui. Je te tiens au courant dès que j'aurai pu essayer
Bonne journée
Pulsar33

Bonjour,

J'ai enfin pu essayer, merci !
Le Fix fonctionne en effet, même si bien sûr mes remarques antérieures sur la nécessité d'un filtrage optimisé demeurent.

Bon courage pour la suite.
Je devrais pouvoir être un peu plus réactif maintenant.
Pulsar33

Bonjour,

Je ne me souviens plus : est-ce qu'on peut injecter un csv (à ton format) en lieu et place du gpx maintenant ?
J'aimerais bien faire des essais de filtrage en attendant que tu aies codé ce qui va bien ...

Bon week-end
Pulsar33

Pas encore Smile

Bonjour,
Des progrès depuis "GPX skip point if same time" ?
Bonne journée
Pulsar33

Bonjour,

Que puis-je faire pour te motiver ? C'est dommage de rester dans cet entre-deux
Si régler le problème du filtrage te semble trop lourd, ajoute au moins l'injection de csv en entrée
Je pourrai ainsi faire des essais de filtrage

Allez, courage !
Bonne journée
Pulsar33

Pulsar33 a écrit:: Bonjour,

Je ne me souviens plus : est-ce qu'on peut injecter un csv (à ton format) en lieu et place du gpx maintenant ?
J'aimerais bien faire des essais de filtrage en attendant que tu aies codé ce qui va bien ...

Bon week-end
Pulsar33

Bonjour,

Pour me servir d'un gps depuis quelques temps (sécurisé par une feuille de papier Smile

), il m'arrive d'utiliser gpsbabel, qui visiblement sait convertir un fichier "csv" (à formater correctement j'imagine) en "gpx".

Quel genre de filtre imagines-tu essayer ? Kalman ?

__________________________________
"Dieu ne réclame aucun droit d’auteur sur la réalité." - M. Anglade

Bonjour O'Leave,

Ton idée est intéressante.

Le soucis en l'absence de poursuite du développement de progweb est que le filtrage actuel des coordonnées d'où découlent la vitesse et la pente ne semble pas adapté. Je sais que progweb s'est cassé les dents sur différents filtres, dont Kalman. Je ne suis pas certain que ce soit judicieux dans notre contexte à partir du moment où la trace est figée et non dynamique. On peut tout à fait utiliser un filtre du deuxième ordre bien calibré et déphaser les coordonnées pour rattraper le retard. C'est ce que nous avions fait en restitution de mission pour un avion d'entrainement et ça donnait de très bons résultats.

J'essaierais bien cette solution et c'est pourquoi j'avais proposé à progweb de pouvoir injecter en lieu et place du gpx un csv modifié au même format que celui que crée son programme afin de ne pas perturber son modèle de données. Il n'a toujours pas trouvé "le temps" de mettre cette option en place. Néanmoins je ne sais pas quel serait l'impact d'un nouveau ré-échantillonnage de ces données lissées par gpx2video.

Du coup je peux utiliser la chaîne suivante :
- uploader mon gpx sur visugps (j'utilise geovelo qui affiche le profil vertical mais n'exporte pas l'altitude dans le fichier gpx, malgré plusieurs rapports de bug à ce sujet de ma part)
- downloader le gpx avec altitude renseignée par visugps puis le convertir en csv compatible de gpsbabel
- modifier le csv pour intégrer les données filtrées en lieu et place des données initiales
- utiliser gpsbabel pour régénérer un gpx tout beau tout propre et le donner à manger à gpx2video
Bon, c'est un peu lourd ... mais pourquoi pas à court terme.

Bonne soirée
Pulsar33

PS : j'avais aussi en tête de faire un bout de code gpxFilter que j'ai d'ailleurs commencé mais pas terminé. On verra bien ...
Le mieux serait quand même que le filtrage de gpx2video soit amélioré

Bonjour

@ O'Leave
Es-tu parti des sources pour gpsbabel ou quelle version toute prête as-tu installée ?

Bonne journée
Pulsar33

Bonjour,

Finalement, j'ai codé le petit outil de filtrage du gpx et le résultat me semble tout à fait satisfaisant
J'ai mis un filtre d'ordre deux sur chaque coordonnée avec un coefficient d'amortissement de 0.7 sur 6 cycles
Il me semble inutile de filtrer en vectoriel ni avec des prédictions pour le résultat attendu et compte tenue de la trajectoire intégralement connue.
Les 6 cycles de décalage du filtre sont éliminés dans le fichier de sortie au prix de la perte des 6 derniers points du fichier ce qui n'est pas un problème.

A toutes fins utiles, voici la partie du code qui s'occupe du filtrage.
Dans mon implémentation, à l'issue du parsing du gpx entrant, les points (trkpt) sont stockés dans les noeuds (node) d'une liste chaînée (trklist)

Code:: //-Traiter un fichier Data
int ProcessData(MyConf *conf){
if (conf->trklist->nbn > 2){
// Initialiser les deux premiers points
listnode *node = conf->trklist->head;
trkpt *pt = (trkpt *)node->ptr;
pt->flat = pt->lat;
pt->flon = pt->lon;
pt->fele = pt->ele;
node = node->succ; pt = (trkpt *)node->ptr;
pt->flat = pt->lat;
pt->flon = pt->lon;
pt->fele = pt->ele;
// Filtrer les coordonnées
double a = conf->cycles;
double a2 = a*a;
double z = conf->damping;
double z2 = z*z;
trkpt *pt1, *pt2;
node = node->succ;
while (node->succ){
pt = (trkpt *)node->ptr;
pt1 = (trkpt *)node->pred->ptr;
pt2 = (trkpt *)node->pred->pred->ptr;
pt->flat = (pt->lat+(pt1->flat*(a+a2/(2*z2)))-(pt2->flat*a2/(4*z2)))/(1+a+a2/(4*z2));
pt->flon = (pt->lon+(pt1->flon*(a+a2/(2*z2)))-(pt2->flon*a2/(4*z2)))/(1+a+a2/(4*z2));
pt->fele = (pt->ele+(pt1->fele*(a+a2/(2*z2)))-(pt2->fele*a2/(4*z2)))/(1+a+a2/(4*z2));
node = node->succ;
}
return conf->trklist->nbn;
}

J'essayerai de publier une vidéo montrant le résultat ...
Bon week-end
Pulsar33

Bonjour,

Je vois que tu ne perds pas la main. Pour répondre à ta question, j'utilise "GPSBabel Version 1.8.0" simplement issue de la distribution linux de l'ordi, LMDE 6.

__________________________________
"Dieu ne réclame aucun droit d’auteur sur la réalité." - M. Anglade

Merci O'Leave

Ça fait un moment que je me demande si je ne vais pas passer de MINT pur à LMDE ...
Ceci dit GPSBabel est dispo pour les deux en effet.
Il y avait plusieurs mois que je n'avais pas touché à du code mais c'est revenu en quelques heures ! C'est cool, ça fait du bien à la tête jocolor

J'ai intégré l'offset de démarrage qui vire les premiers points inutiles aisément (ma vidéo et mon gpx ne sont pas synchronisés nativement)
Il me reste à pouvoir offseter l'heure car gpx2video ne semble pas tenir compte correctement de ce que lui dit le gpx.
Je croyais l'avoir déjà signalé à progweb, mais je ne retrouve pas mon message
Il faudra que j'explique clairement le problème.

Bonne soirée
Pulsar33

Bonjour,

Voici une vidéo dont la partie haute est l'incrustation générée par gpx2video avec les données brutes et la partie basse celle générée avec les données filtrées.
La synchronisation entre les deux vidéos est à quelques dixèmes de seconde près, d'où l'écho léger qu'on entend

Première moitié de la vidéo : la trajectoire dans les ronds-points
On voit clairement que le filtre réduit l'amplitude des coordonnées. Il en résulte un aplatissement de la trajectoire dans les ronds-points assez visible
Il faudrait donc conserver les positions non filtrées car à l'échelle de l'affichage, leurs variations non désirées ne sont pas visibles

Deuxième moitié de la vidéo : le filtrage de la pente (et de la vitesse)
Sur un faux-plat assez long, on voit nettement que la pente non filtrée dit n'importe quoi
C'est aussi visible dans la première séquence mais moins flagrant
La pente non filtrée varie de -7% à +12%, avant et après les ponts, et même à -12% fugitivement
De telles variations se font parfois en une à deux secondes
A l'inverse, la pente filtrée reste sagement autour de 1% ce qui correspond très bien à la réalité
Même si c'est moins facile à voir, le rendu de la vitesse filtrée est aussi nettement meilleur
Il faudrait donc utiliser les positions filtrées pour calculer la pente, la vitesse et ce qui en découle

Hélas, de l'extérieur du programme, je ne fournis que les positions, donc filtrées ou non filtrées.
Seul gpx2video pourrait en interne utiliser les non filtrées pour les positions et la trajectoire, et les filtrées pour les valeurs dérivées
Je peux jouer sur le nombre de cycles et le coefficient d'amortissement mais je ne pense pas trouver de meilleur compromis

progweb, la balle est dans ton camp !
Bon dimanche
Pulsar33

Ça va arriver, juste un peu de patience.

Bonjour,

OK ...

Tant que tu y es, pourrais-tu afficher cette erreur même si -q est actif stp ?
ERROR: Parsing of 'Jonzac09VGO_F.gpx' failed due to mismatched tag on line 12 and column 4
WARNING: Can't read GPS data, skip map widget

Quand on bidouille les gpx, ce serait bien d'être averti même si tous les autres affichages console ne sont pas souhaitables

Bon courage
Pulsar33

Et un détail pour les puristes dont je suis :
Les pentes à "-0%" c'est un peu choquant

Bonne journée
Pulsar33

Une petite remarque sur le filtre, pour utiliser gpx2video, ce n'est pas les positions lat et l'on que tu dois corriger, mais le timestamp de chaque point et surtout calculer un timestamp à la seconde (tant que gpx2video ne tiendra pas compte des milisecondes).
Ensuite une fois le timestamp corrigé, tu oeux corriger l'altitude.

Mais cela va arriver....

» GNU/Linux!!
» CAO sous Linux
» GPS compatible linux Ubuntu
» Mac ou Pc ?
» navigateur links sous Linux