File:BAG import graph.png

Summary

Progress of BAG import in the Netherlands. The number of buildings with source tags '3dShapes' and 'BAG' is shown, as well as the number of buildings with other source tags. Buildings without source tag are included in the Other/Overig category.

Data

Computing the data

Current method

• OSM data for the Netherlands (netherlands-latest.osm.pbf) downloaded from Geofabrik.
• Convert to o5m format with Osmconvert:

osmconvert netherlands-latest.osm.pbf -o=nl.o5m

• Filter with Osmfilter, keeping only buildings:

osmfilter nl.o5m --keep="building=" -o=buildings.o5m

• The Geofabrik file contains a small strip of land across the border. Buildings that are not inside the Netherlands need to be discarded. Download the border of the Netherlands from Openstreetmap and convert to a .poly file "Netherlands.poly". Then clip to this polygon:

osmconvert buildings.o5m -B=Netherlands.poly -o=buildings-clipped.o5m

• Filter buildings and save statistics:

osmfilter buildings-clipped.o5m --keep="building=" --ignore-dependencies --drop-nodes --out-count=source > buildingsource.txt
osmfilter buildings-clipped.o5m --keep="building=" --ignore-dependencies --drop-nodes --out-count > buildingkeys.txt

Old method, using a PostgreSQL database

• Buildings filtered with Osmosis:

osmosis --read-pbf file="netherlands-latest.osm.pbf" --tf reject-nodes 
    --tf accept-ways building=* --tf accept-relations building=* --write-pgsql-dump

• Loaded in database using 'pgsnapshot_load_0.6.sql' from Osmosis
• SQL query:

select source, count(*) as "num" from 
 (select lower(tags->'source') as "source" from ways union all 
  select lower(tags->'source') as "source" from relations) as "combine"
group by source

Gnuplot script

set term pngcairo size 728, 410 font "Helvetica,12"
set output "BAG_import_graph.png"
set xdata time
set timefmt "%d-%b-%y"
set xrange ["1-Nov-13":"30-Sep-14"]
set format x "%b-%y"
unset mxtics
set yrange [0:10]
set tics out nomirror scale 0.75
set xtics 2592000 textcolor rgb "black"
set ytics 2 textcolor rgb "black"
set ylabel "Aantal (in miljoen)"
set style data linespoints
set border lt rgb "#757575"
set grid ls 1 lt rgb "#D9D9D9"
set key inside c r samplen 2 spacing 1.65  reverse Left
plot "data.txt" using 1:($3/1e6) lw 2.5 lt rgb "#C0504D" pt 18 ps 0.8 title "BAG",\
     "data.txt" using 1:($2/1e6) lw 2.5 lt rgb "#4F81BD" pt 26 ps 1.1 title "3dShapes",\
     "data.txt" using 1:($4/1e6) lw 2.5 lt rgb "#9BBB59" pt 22 ps 0.8 title "Overig"

Licensing