calais-data

Des données sur le passage à Calais - retour accueil

git clone git://bebou.netlib.re/calais-data
Log | Files | Refs | 

commit 772327b95216067e9da34fb6264b34c12fb0bd03
parent d45b00634e6b568e792811210af7658bcd2243bb
Auterice: Arthur Pons <arthur.pons@unistra.fr>
Date:   Thu, 10 Oct 2024 15:58:31 +0200

Plus de progression

Diffstat:

Mfigures/fig1.gp | 20+++++++++++---------


Mfigures/fig1.sh | 5++---


Mfigures/fig10.gp | 46+++++++++++++++++++++++++++++++++++++++++++---


Mfigures/fig10.sh | 14+++++++++-----


Mfigures/fig11.gp | 3++-


Mfigures/fig13.gp | 3++-


Mfigures/fig2.gp | 3++-


Mfigures/fig3.gp | 16++++++++++------


Mfigures/fig3.sh | 4+++-


Mfigures/fig4.gp | 14++++++--------


Mfigures/fig4.sh | 3+--


Mfigures/fig5.gp | 14++++++++------


Mfigures/fig5.sh | 5++---


Mfigures/fig6.gp | 31++++++++++---------------------


Mfigures/fig6.sh | 20++------------------


Mfigures/fig7.gp | 32++++++++++++++++----------------


Mfigures/fig7.sh | 26++++++++++++++++++--------


Mfigures/fig8.gp | 19+++++++++++--------


Mfigures/fig8.sh | 10+++++-----


Mfigures/fig9.gp | 27+++++++++++++++++++--------


Mfigures/fig9.sh | 9++++++++-


Mplan | 1+


22 files changed, 191 insertions(+), 134 deletions(-)

diff --git a/figures/fig1.gp b/figures/fig1.gp
@@ -1,12 +1,14 @@
 #!/usr/bin/gnuplot
 
-set term png size 1000,600
+set term pngcairo size 1000,600
+set linetype 1 lc rgb "#0d41a1"
 
-set style fill solid border -1
-set boxwidth 0.9
-set style histogram gap 0
-
-set title "Nombre de personnes passées par vitesse de vent"
-set ylabel "Nombre de personnes"
-set xlabel "Vitesse de vent (miles/heure)"
-plot 'figures/fig1.tsv' u 2:xtic(1) w histo notitle
+set title "La vitesse du vent en fonction des jours"
+set ylabel "Vitesse de vent moyenne (miles/heure)"
+set xlabel "Nombre de jours depuis le 1er janvier 2018"
+set yrange [0:40]
+set xr [0:2500]
+set xtics autofreq 364
+plot '+' u (2225):(35) w p ps 1 pt 7 lc 1 title "Passage"        ,\
+     '+' u (2225):(35) w p ps 1 pt 7 lc 4 title "Pas de passage" ,\
+'figures/fig1.tsv' u :3:::($2>0 ? 1:4) w p pt 7 lc variable notitle
diff --git a/figures/fig1.sh b/figures/fig1.sh
@@ -1,6 +1,5 @@
 #! /bin/sh
 
-< full.tsv cut -f 2,5 |
-	grep -v nan | cut -d'.' -f1 | sort -t'	' -nk2 |
-	tsv-summarize -H -g 2 --sum 1
+< full.tsv cut -f1,2,5 |
+	tail -n+2
 
diff --git a/figures/fig10.gp b/figures/fig10.gp
@@ -1,5 +1,45 @@
 #!/usr/bin/gnuplot
 
-set term png size 1000,600
-set key autotitle columnhead
-plot 'figures/fig10.tsv' u 2
+# https://stackoverflow.com/questions/63771600/is-there-a-way-to-have-3-different-y-axes-on-one-graph-using-gnuplot
+
+set term pngcairo size 1000,600
+set linetype 1 lc rgb "#0d41a1"
+
+set margins screen 0.1, screen 0.77, screen 0.1, screen 0.90
+
+set multiplot
+
+set boxwidth 1
+set style fill solid 1 border -1
+set style data histograms
+set style histogram gap 1
+
+set ylabel "Nombre de jours"
+set title "Personnes, jours de passage et personnes par bateau par année"
+
+set yrange [0:365]
+set ytics autofreq 26 nomirror
+set y2tics autofreq 3900
+set y2range [0:54600]
+set y2label "Nombre de personnes"
+set xtics scale 0
+set grid
+
+plot 'figures/fig10.tsv' using 2:xtic(1) t "Nombre de personnes" axes x1y2 lc 1,\
+'' using 3:xtic(1) t "Nombre de jours avec passage" lc 4,\
+'' u ($4*600):xtic(1) axes x1y2 t "Nombre moyen de personnes par bateau" lc 2
+
+unset xlabel
+unset ylabel
+unset y2label
+unset tics
+
+set rmargin at screen 0.91
+set title " "
+set border 8
+set y2label "Nombre moyen de personnes par bateau" offset 0,0
+set y2tics nomirror offset 0,0
+set y2range [0:100]
+set y2tics autofreq 10
+plot NaN notitle
+
diff --git a/figures/fig10.sh b/figures/fig10.sh
@@ -1,8 +1,12 @@
 #! /bin/sh
 
-< full.tsv cut -f1,4 |
-	grep -v '0$' |
-	sed -E 's,/[0-9]+	,	,' |
-	tsv-summarize -H -g 1 --mean 2 |
-	sed '1 s/_/-/'
+tmpd=$(mktemp -d)
+< full.tsv cut -f 1,2,4 |
+	sed -E 's,/[0-9]+/[0-9]+	,	,' |
+	sed -E 's,	0	,		,;s,0$,,' > $tmpd/data
+< /$tmpd/data tsv-summarize -H -g 1 --not-missing-count 2 > $tmpd/jpassage
+< /$tmpd/data tsv-summarize -H -g 1 --sum 2 --exclude-missing > $tmpd/nbpassage
+< /$tmpd/data tsv-summarize -H -g 1 --mean 3 --exclude-missing > $tmpd/densite
+< $tmpd/jpassage tsv-join -f - -k 1 -d 1 -a 2 $tmpd/nbpassage | tail -n+2 > $tmpd/a
+< $tmpd/densite tsv-join -f - -k 1 -d 1 -a 2 $tmpd/a
 
diff --git a/figures/fig11.gp b/figures/fig11.gp
@@ -1,6 +1,7 @@
 #!/usr/bin/gnuplot
 
-set term png size 1000,600
+set term pngcairo size 1000,600
+set linetype 1 lc rgb "#0d41a1"
 set key autotitle columnhead
 set style fill solid
 set style histo gap 2
diff --git a/figures/fig13.gp b/figures/fig13.gp
@@ -1,6 +1,7 @@
 #!/usr/bin/gnuplot
 
-set term png size 1000,600
+set term pngcairo size 1000,600
+set linetype 1 lc rgb "#0d41a1"
 set key autotitle columnhead
 set style fill solid
 set style histo gap 1
diff --git a/figures/fig2.gp b/figures/fig2.gp
@@ -1,7 +1,8 @@
 #!/usr/bin/gnuplot
 
-set term png size 1000,600
+set term pngcairo size 1000,600
 set style fill solid
+set linetype 1 lc rgb "#0d41a1"
 
 set boxwidth 0.7
 set style fill solid 1.0 border -1
diff --git a/figures/fig3.gp b/figures/fig3.gp
@@ -1,9 +1,13 @@
 #!/usr/bin/gnuplot
 
-set term png size 1000,600
-set key autotitle columnhead
-set title "Nombre de personnes exilées arrivées en Angleterre par jour"
-set xtics autofreq 365
+set term pngcairo size 1000,600
+set linetype 1 lc rgb "#0d41a1"
+
+set style fill solid border -1
+set boxwidth 0.9
+set style histogram gap 0
+
+set title "Nombre de personnes passées par vitesse de vent"
 set ylabel "Nombre de personnes"
-set xlabel "Nombre de jours depuis le 1er janvier 2018"
-plot 'figures/fig3.tsv' u 2 w p pt 7
+set xlabel "Vitesse de vent (miles/heure)"
+plot 'figures/fig3.tsv' u 2:xtic(1) w histo notitle
diff --git a/figures/fig3.sh b/figures/fig3.sh
@@ -1,4 +1,6 @@
 #! /bin/sh
 
-< full.tsv cut -f1,2
+< full.tsv cut -f 2,5 |
+	grep -v nan | cut -d'.' -f1 | sort -t'	' -nk2 |
+	tsv-summarize -H -g 2 --sum 1
 
diff --git a/figures/fig4.gp b/figures/fig4.gp
@@ -1,12 +1,10 @@
 #!/usr/bin/gnuplot
 
-set term png size 1000,600
+set term pngcairo size 1000,600
+set linetype 1 lc rgb "#0d41a1"
 set key autotitle columnhead
-set title "Nombre de personnes exilées arrivées en Angleterre par semaines depuis janvier 2018"
-set xtics autofreq 51
+set title "Nombre de personnes exilées arrivées en Angleterre par jour"
+set xtics autofreq 365
 set ylabel "Nombre de personnes"
-set xlabel "Nombre de semaines depuis janver 2018"
-set style fill solid
-set style histo gap 0
-set xrange [0:370]
-plot 'figures/fig4.tsv' u 2 t "" w histo
+set xlabel "Nombre de jours depuis le 1er janvier 2018"
+plot 'figures/fig4.tsv' u 2 w p pt 7
diff --git a/figures/fig4.sh b/figures/fig4.sh
@@ -1,5 +1,4 @@
 #! /bin/sh
 
-< weekdays.tsv cut -f1,9
-
+< full.tsv cut -f1,2
 
diff --git a/figures/fig5.gp b/figures/fig5.gp
@@ -1,12 +1,14 @@
 #!/usr/bin/gnuplot
 
-set term png size 1000,600
+set term pngcairo size 1000,600
+set linetype 1 lc rgb "#0d41a1"
 set key autotitle columnhead
-set title "Nombre de personnes exilées arrivées en Angleterre par mois depuis janvier 2018"
-set xtics autofreq 12
+set title "Nombre de personnes exilées arrivées en Angleterre par semaines depuis janvier 2018"
+set xtics autofreq 51
 set ylabel "Nombre de personnes"
-set xlabel "Nombre de mois depuis janver 2018"
-set style fill solid border -1
+set xlabel "Nombre de semaines depuis janver 2018"
+set style fill solid
 set style histo gap 0
-set xrange [0:85]
+set xrange [0:355]
+do for [i=4:355:4]{ set arrow from i,0 to i,4000 lt 0 nohead }
 plot 'figures/fig5.tsv' u 2 t "" w histo
diff --git a/figures/fig5.sh b/figures/fig5.sh
@@ -1,6 +1,5 @@
 #! /bin/sh
 
-< full.tsv cut -f1,2 |
-	sed -E 's,/[0-9]+	,	,' |
-	tsv-summarize -H -g 1 --sum 2
+< weekdays.tsv cut -f1,9
+
 
diff --git a/figures/fig6.gp b/figures/fig6.gp
@@ -1,24 +1,13 @@
 #!/usr/bin/gnuplot
 
-set term png size 1000,1800
-
-set multiplot layout 3, 1
-set tmargin 4
-
+set term pngcairo size 1000,600
+set linetype 1 lc rgb "#0d41a1"
+set key autotitle columnhead
+set title "Nombre de personnes exilées arrivées en Angleterre par mois depuis janvier 2018"
+set xtics autofreq 12
+set ylabel "Nombre de personnes"
+set xlabel "Nombre de mois depuis janver 2018"
 set style fill solid border -1
-set boxwidth 0.9
-set style histogram gap 0
-
-set ylabel "Nombre de personnes passées"
-set xlabel
-set xtics scale 0
-
-set title "Nombre de personnes exilées passées par années"
-plot 'figures/fig6.tsv' u 2:xtic(1) index 0 w histo notitle
-set title "Nombre de personnes exilées passées par mois"
-plot 'figures/fig6.tsv' u 2:xtic(1) index 1 w histo notitle
-set title "Nombre de personnes exilées passées par jour de la semaine"
-set yrange [0:]
-plot 'figures/fig6.tsv' u 2:xtic(1) index 2 w histo notitle
-
-unset multiplot
+set style histo gap 0
+set xrange [0:85]
+plot 'figures/fig6.tsv' u 2 t "" w histo
diff --git a/figures/fig6.sh b/figures/fig6.sh
@@ -1,22 +1,6 @@
 #! /bin/sh
 
 < full.tsv cut -f1,2 |
-	sed -E 's,/[0-9]+/[0-9]+	,	,' |
-	tail -n+2 |
-	tsv-summarize -g 1 --sum 2
-
-echo;echo
-
-< full.tsv cut -f1,2 |
-	sed -E 's,[0-9]+/([0-9]+)/[0-9]+	,\1	,' |
-	tail -n+2 |
-	tsv-summarize -g 1 --sum 2
-
-echo;echo
-
-a=$(mktemp);b=$(mktemp)
-< weekdays.tsv head -n1 | cut -f2-8 | tr '\t' '\n' > $a
-< weekdays.tsv tail -n1 | tr '\t' '\n' | head -n-1 | grep -v '^$' > $b
-paste $a $b
-rm $a $b
+	sed -E 's,/[0-9]+	,	,' |
+	tsv-summarize -H -g 1 --sum 2
 
diff --git a/figures/fig7.gp b/figures/fig7.gp
@@ -1,25 +1,25 @@
 #!/usr/bin/gnuplot
 
-set term png size 1000,1200
+set term pngcairo size 1000,1000
+set linetype 1 lc rgb "#0d41a1"
 
-set multiplot layout 2, 1
+set multiplot layout 3, 1
 set tmargin 4
 
 set style fill solid border -1
+set boxwidth 0.9
+set style histogram gap 0
 
-set ylabel "Nombre moyen de personnes dans un bateau"
-set xlabel "Nombre de jours avec passage depuis le 1er janvier 2018"
-set title "Evolution du nombre moyen de personnes dans un bateau par jour de passage"
-set xtics autofreq 364
-set xrange [0:850]
-plot 'figures/fig7.tsv' index 0 pt 7 notitle
-
-unset xrange
-set style fill solid 0.50 border lt -1;
-set pointsize 0.5;
-set title "Répartition du nombre de personne par bateau les jours de passage par année"
-set ylabel "Nombre moyen de personnes dans un bateau"
-unset xlabel
+set ylabel "Nombre de personnes passées"
+set xlabel
 set xtics scale 0
-plot 'figures/fig7.tsv' using (2):2:(0):1 index 1 with boxplot lc variable notitle
 
+set title "Nombre de personnes exilées passées par années"
+plot 'figures/fig7.tsv' u 2:xtic(1) index 0 w histo notitle
+set title "Nombre de personnes exilées passées par mois"
+plot 'figures/fig7.tsv' u 2:xtic(1) index 1 w histo notitle
+set title "Nombre de personnes exilées passées par jour de la semaine"
+set yrange [0:]
+plot 'figures/fig7.tsv' u 2:xtic(1) index 2 w histo notitle
+
+unset multiplot
diff --git a/figures/fig7.sh b/figures/fig7.sh
@@ -1,12 +1,22 @@
 #! /bin/sh
 
-< full.tsv cut -f4 |
-	grep -v '^0$' |
-	tail -n+2
-echo;echo
-< full.tsv cut -f1,4 |
-	grep -v '0$' |
+< full.tsv cut -f1,2 |
 	sed -E 's,/[0-9]+/[0-9]+	,	,' |
-	#sed -E 's,/[0-9]+	,	,' |
-	tail -n+2
+	tail -n+2 |
+	tsv-summarize -g 1 --sum 2
+
+echo;echo
+
+< full.tsv cut -f1,2 |
+	sed -E 's,[0-9]+/([0-9]+)/[0-9]+	,\1	,' |
+	tail -n+2 |
+	tsv-summarize -g 1 --sum 2
+
+echo;echo
+
+a=$(mktemp);b=$(mktemp)
+< weekdays.tsv head -n1 | cut -f2-8 | tr '\t' '\n' > $a
+< weekdays.tsv tail -n1 | tr '\t' '\n' | head -n-1 | grep -v '^$' > $b
+paste $a $b
+rm $a $b
 
diff --git a/figures/fig8.gp b/figures/fig8.gp
@@ -1,22 +1,25 @@
 #!/usr/bin/gnuplot
 
-set term png size 1000,1200
+set term pngcairo size 1000,1200
+set linetype 1 lc rgb "#0d41a1"
 
 set multiplot layout 2, 1
 set tmargin 4
 
-set title "Evolution de la vitesse de vent moyenne par jour de passage"
-set xlabel "Nombre de jours avec passage depuis le 1er janvier 2018"
-set ylabel "Vitesse du vent (miles/h)"
-set xrange [0:850]
+set style fill solid border -1
 
-plot 'figures/fig8.tsv' index 0 notitle
+set ylabel "Nombre moyen de personnes dans un bateau"
+set xlabel "Nombre de jours avec passage depuis le 1er janvier 2018"
+set title "Evolution du nombre moyen de personnes dans un bateau par jour de passage"
+set xtics autofreq 364
+set xrange [0:2550]
+plot 'figures/fig8.tsv' u ($1>0 ? $1:NaN) index 0 pt 7 notitle
 
 unset xrange
-
 set style fill solid 0.50 border lt -1;
 set pointsize 0.5;
-set title "Répartition des vitesses de vent moyenne par jour de passage par année"
+set title "Répartition du nombre de personne par bateau les jours de passage par année"
+set ylabel "Nombre moyen de personnes dans un bateau"
 unset xlabel
 set xtics scale 0
 plot 'figures/fig8.tsv' using (2):2:(0):1 index 1 with boxplot lc variable notitle
diff --git a/figures/fig8.sh b/figures/fig8.sh
@@ -1,11 +1,11 @@
-#!/bin/sh
+#! /bin/sh
 
-< full.tsv cut -f2,5 |
-	grep -v '^0	' | cut -f2 |
+< full.tsv cut -f4 |
+#	grep -v '^0$' |
 	tail -n+2
 echo;echo
-< full.tsv cut -f1,2,5 |
-	grep -v '	0	' | cut -f1,3 |
+< full.tsv cut -f1,4 |
+	grep -v '0$' |
 	sed -E 's,/[0-9]+/[0-9]+	,	,' |
 	#sed -E 's,/[0-9]+	,	,' |
 	tail -n+2
diff --git a/figures/fig9.gp b/figures/fig9.gp
@@ -1,13 +1,24 @@
 #!/usr/bin/gnuplot
 
-set term png size 1000,600
+set term pngcairo size 1000,1000
+set linetype 1 lc rgb "#0d41a1"
 
-set title "La vitesse du vent en fonction des jours"
-set ylabel "Vitesse de vent moyenne (miles/heure)"
+set multiplot layout 2, 1
+set tmargin 4
+
+set title "Evolution de la vitesse de vent moyenne les jours de passage"
 set xlabel "Nombre de jours depuis le 1er janvier 2018"
-set yrange [0:40]
-set xr [0:2500]
+set ylabel "Vitesse du vent (miles/h)"
+set xrange [0:2550]
 set xtics autofreq 364
-plot '+' u (2225):(35) w p ps 1 pt 7 lc 1 title "Passage"        ,\
-     '+' u (2225):(35) w p ps 1 pt 7 lc 4 title "Pas de passage" ,\
-'figures/fig9.tsv' u :3:::($2>0 ? 1:4) w p pt 7 lc variable notitle
+plot 'figures/fig9.tsv' u ($1>0?$2:NaN) index 0 pt 7 notitle
+
+unset xrange
+
+set style fill solid 0.50 border lt -1;
+set pointsize 0.5;
+set title "Répartition des vitesses de vent moyenne par jour de passage par année"
+unset xlabel
+set xtics scale 0
+plot 'figures/fig9.tsv' using (2):2:(0):1 index 1 with boxplot lc variable notitle
+
diff --git a/figures/fig9.sh b/figures/fig9.sh
@@ -1,5 +1,12 @@
-#! /bin/sh
+#!/bin/sh
 
+< full.tsv cut -f2,5 |
+	#grep -v '^0	' | cut -f2 |
+	tail -n+2
+echo;echo
 < full.tsv cut -f1,2,5 |
+	grep -v '	0	' | cut -f1,3 |
+	sed -E 's,/[0-9]+/[0-9]+	,	,' |
+	#sed -E 's,/[0-9]+	,	,' |
 	tail -n+2
 
diff --git a/plan b/plan
@@ -7,3 +7,4 @@ nombre de passages par années mois et jour de la semaine
 évolution de la densité des bateaux avec le temps par jour et par tranche de un an
 vitesse de vent moyenne lors des jours de passage par jour et par année
 vitesse de vent moyenne par jour, en couleur 1 les passages, en couleur 2 les non passage
+Personnes passées, jours de passage et personnes par bateau par année