Yritän yhdistää Maanmittauslaitoksen polygon-shapefilejä ogr2ogr:n avulla, mutta lopputuloksessa on ihmeellisiä “haamurakennuksia” oikeiden rakennusten päällä. Näitä viallisia rakennuksia ei oikeasti kartalla tai yksittäisissä shapefileissä ole. Gdalin versio on 2.2.1.
Tässä on bash-skripti:
#!/bin/bash
file="rakennukset.shp"
for i in $(ls r*p.shp)
do
if [ -f "$file" ]
then
ogr2ogr -f 'ESRI Shapefile' -append -update -t_srs EPSG:4326 -nln rakennukset $file $i
else
ogr2ogr -f 'ESRI Shapefile' -t_srs EPSG:4326 $file $i
fi
done
Lopputulos näyttää Qgis:ssä tältä: https://imgur.com/a/pZZeJ
Missä vika? Kaikki rakennukset eivät ole päällekkäin, suurin osa näkyy oikein. Joissain paikoissa on rakennus sellaisessakin kohdassa, jossa ei ole lähimaillakaan muita rakennuksia.
Edit: Oikeastaan ne rakennukset, jotka eivät ole vinossa, ovat sitten tismalleen toisen päällä. Eli “oikealta” näyttävät rakennukset ovatkin kaksinkertaisina shapefilessä. Kummallista.