Landuse overlap / Erstes experiment / Frage

Florian Lohoff f at zz.de
Di Apr 10 14:38:11 CEST 2018


Hi,

On Tue, Apr 10, 2018 at 10:29:15AM +0200, Hartmut Holzgraefe wrote:
> On 10.04.2018 08:51, Florian Lohoff wrote:
> 
> > Hier ist mal eine Karte mit (Multi)Polygonen wo das nicht der
> > fall ist d.h. wo landuses übereinanderlappen ohne entsprechende
> > multipolygon relations.
> > 
> > https://osm.zz.de/dbview/?db=landuseoverlap-owl&layer=overlap#51.87342,8.54548,12z
> 
> Ich werde da gerade nicht ganz schlau draus.
> 
> * die Links in den Popups funktionieren nicht, enthalten weder OSM-IDs
>   noch Koordinaten
> 
> * bei manchen der Polygone verstehe ich überhaupt nicht wo die herkommen
> 
> zB: https://osm.zz.de/dbview/?db=landuseoverlap-owl&layer=overlap#52.03128,8.52797,17z
> 
> Das kleinere Polygon da ist klar, da liegt ein "landuse=construction" oben
> auf "landuse=residential". Dort ist kürzlich ein Haus abgerissen
> worden und wird nun neu bebaut.
> 
> Das größere Polygon ist mir aber ein absolutes Rätsel, das finde ich
> in den OSM-Daten nicht ...?
> 
> Oder hat da schon jemand kürzlich aufgeräumt?

Jein - ich schraube noch - Die Popups waren eben kaputt. Beim nächsten
Rechnen stehen sinnvolle dinge drin.

Das die Geometrieen manchmal nicht mit den OSM dingen übereinstimmen
liegt daran das ich ein OGRIntersection bilde aus den beiden
überlappenden Flächen. Da kommt aber zum teil so ein schrott bei raus:

GEOMETRYCOLLECTION (POINT (6.3797753 50.5918261),LINESTRING (6.3803743 50.5915808,6.3809734 50.591472),LINESTRING (6.3809734 50.591472,6.3815604 50.591392),LINESTRING (6.3815604 50.591392,6.3818183 50.5912774),LINESTRING (6.3818183 50.5912774,6.3820918 50.5910856),LINESTRING (6.3820918 50.5910856,6.3825175 50.5905126),LINESTRING (6.3825175 50.5905126,6.3826134 50.5903178),LINESTRING (6.3826134 50.5903178,6.38262 50.5901483),LINESTRING (6.38262 50.5901483,6.3826064 50.5899329),LINESTRING (6.3826064 50.5899329,6.3825387 50.5898106),POLYGON ((6.38219007319966 50.5893783153194,6.3825387 50.5898106,6.3821985 50.5893667,6.38219007319966 50.5893783153194)))

Sowas kann die SQLite scheinbar nicht. Ich kann das als Geometrycollection da
reinschreiben, aber raus kommt mit AsGeoJSON was syntaktisch kaputtes.

Dazu kommt das ja auch Leaflet mit einer Geometrycollection nicht viel
anfangen wird können (Wie render ich einen Node wenn ich eine Stroke Breite habe?)

Deshalb gehe ich bei solchen konstrukten hin und mache ein OGRConvexHull über
die Geometrycollection. Damit liegen alle Fehlstellen innerhalb des Polygons.

Ich hoffe es wird klarer warum die teilweise so aussehen.

In den Popups stehen jetzt auch brauchbare Informationen.

Niedersachsen:
https://osm.zz.de/dbview/?db=landuseoverlap-nds&layer=overlap#52.13077,8.11516,14z

NRW:
https://osm.zz.de/dbview/?db=landuseoverlap-nrw&layer=overlap#52.08309,7.86587,11z

OWL:
https://osm.zz.de/dbview/?db=landuseoverlap-owl&layer=overlap#51.99133,8.41707,12z

Im moment berechne ich die alle 20 Minuten. Das klappt gerade nicht mehr immer
deshalb evtl die großen NRW/NDS demnächst nicht mehr so häufig.

Flo
-- 
Florian Lohoff                                                 f at zz.de
             UTF-8 Test: The 🐈 ran after a 🐁, but the 🐁 ran away


Mehr Informationen über die Mailingliste OSM