Vergelijking tussen Wadden en AIS data op de Waddenzee

Introductie

Onderstaande kaarten geven een impressie en eerste aanzet tot discussie over het nut van het gebruik van radar data en ais data voor het in kaart brengen van scheepsbewegingen op het wad.

Een paar methodologische punten:

De data die te zien is in onderstaande kaarten betreft de data van juni 2017.

De punten zijn geaggregeerd tot vlakken van ~ 100 bij 100 meter.

De data is weergegeven als het natuurlijke logarithme van de puntendichtheid. Reden hiervoor is dat dit het beste past bij de distributie van de puntendichtheid (~ Poisson-distributie). Hieronder staat de distributie over de gehele breedte van de Wadden, voor zowel radar als AIS op 4 juni.

## Warning: Removed 1876994 rows containing non-finite values (stat_bin).
## Warning: Removed 1903131 rows containing non-finite values (stat_bin).

Deze diagram laat vrijwel niets zien, aangezien er een aantal hoge outliers zijn met hele hoge puntendichtheden (maar die weinig voorkomen). Door in te zoomen op de distributie van alle vlakken met minder dan 100 punten zien we beter wat er gebeurt qua schependichtheid op het wad:

## Warning: Removed 1876994 rows containing non-finite values (stat_bin).
## Warning: Removed 1903131 rows containing non-finite values (stat_bin).

Deze histogram laat zien dat de meeste meetvlakken op het wad slechts 1 keer een boot registreren (en de Radar registreert iets meer vlakken met 1 punt). De AIS laat een mooie verloop zien die hoort bij een poisson verdeling, voor de radar loopt dit iets anders. Deze heeft, relatief aan wat we zouden verwachten, veel vlakken met rondom de 5 tot 10 markers. Dit zijn naar alle waarschijnlijkheid punten die constant aanwezig zijn op het wad, te denken valt aan boeien en de zaadinvanginstallaties.

De reden dat ik 4 juni laat zien is omdat dit een dag was met relatief veel pleziervaart (pinksterzondag), en een dag met mooi weer (geen regen, windkracht 3 tot 4, KNMI: Daggegevens weer Lauwersoog). Hierdoor levert de radar data nog een redelijk natuurlijk beeld op.

Als we de data van 4 juni naast bijvoorbeeld een regenachtige dag, met aardig wat wind leggen (7 juni, woensdag, 9.5 uur regen, windkracht 5 tot 7), dan zien we een meer piekend verloop van de radar data. De ruis tussen de 5 tot 10 punten per vlak is behoorlijk toegenomen, en ook verder in de distributie komen er meer onwaarschijnlijkheden voor.

Vergelijkende kaarten

In onderstaande kaart vergelijken we twee dagen, 4 juni 2017 (Pinkerzondag, zonnig), en 7 juni (woensdag, regenachtig). Het gaat hierbij om de Radar data over deze twee dagen genomen.

Wat bovenstaande kaart duidelijk laat zien is dat er behoorlijk wat valide data in de radar gegevens zit. De wadden zijn aanzienlijk drukker op Pinksterzondag,

## Warning in colors(.): Some values were outside the color scale and will be
## treated as NA
mapaisradar