Add analyser analyser_osmosis_highway_zone osm-fr#326

analysers/analyser_osmosis_highway_zone.py

@@ -0,0 +1,126 @@
+#!/usr/bin/env python
+#-*- coding: utf-8 -*-
+
+###########################################################################
+##                                                                       ##
+## Copyrights Frédéric Rodrigo 2018                                      ##
+##                                                                       ##
+## This program is free software: you can redistribute it and/or modify  ##
+## it under the terms of the GNU General Public License as published by  ##
+## the Free Software Foundation, either version 3 of the License, or     ##
+## (at your option) any later version.                                   ##
+##                                                                       ##
+## This program is distributed in the hope that it will be useful,       ##
+## but WITHOUT ANY WARRANTY; without even the implied warranty of        ##
+## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         ##
+## GNU General Public License for more details.                          ##
+##                                                                       ##
+## You should have received a copy of the GNU General Public License     ##
+## along with this program.  If not, see <http://www.gnu.org/licenses/>. ##
+##                                                                       ##
+###########################################################################
+
+from .Analyser_Osmosis import Analyser_Osmosis
+
+# Cluster Zone with buffer
+sql10 = """
+DROP TABLE IF EXISTS a0;
+CREATE TEMP TABLE a0 AS
+SELECT
+  (ST_Dump(
+    ST_Union(ST_Buffer(linestring_proj, 60))
+  )).path[1] AS cid,
+  (ST_Dump(
+    ST_Union(ST_Buffer(linestring_proj, 60))
+  )).geom AS geom
+FROM
+  highways
+WHERE
+  tags?'zone:maxspeed' AND
+  highways.tags->'zone:maxspeed' LIKE '%:{0}'
+"""
+
+sql11 = """
+DROP TABLE IF EXISTS a1;
+CREATE TEMP TABLE a1 AS
+SELECT
+  cid,
+  (ST_DumpRings(geom)).path[1] AS path,
+  (ST_DumpRings(geom)).geom
+FROM
+  a0
+"""
+
+# Remove small inner hole
+sql12 = """
+DROP TABLE IF EXISTS a2;
+CREATE TEMP TABLE a2 AS
+SELECT
+  cid,
+  path,
+  ST_ExteriorRing(geom) AS geom
+FROM
+  a1
+WHERE
+  path = 0 OR
+  ST_Area(geom) > 100 * 100
+"""
+
+# Rebuild cluster and shrink the buffer back
+sql13 = """
+DROP TABLE IF EXISTS a3;
+CREATE TEMP TABLE a3 AS
+SELECT
+  ST_Buffer(CASE
+    WHEN t.inners IS NULL THEN ST_MakePolygon(a.geom)
+    ELSE ST_MakePolygon(a.geom, t.inners)
+  END, -61) AS geom
+FROM
+  (SELECT cid, geom FROM a2 WHERE path = 0) AS a
+  LEFT JOIN
+    (SELECT cid, ST_Accum(geom) AS inners FROM a2 WHERE path > 0 GROUP BY cid) AS t ON
+    a.cid = t.cid
+"""
+
+# Get candidates ways in the resulting envelope
+sql14 = """
+SELECT
+  highways.id,
+  ST_AsText(way_locate(linestring))
+FROM
+  a3
+  JOIN highways ON
+    highways.tags != ''::hstore AND
+    highways.tags?'highway' AND
+    highways.tags->'highway' IN ('residential', 'unclassified') AND
+    (NOT highways.tags?'zone:maxspeed' OR NOT highways.tags->'zone:maxspeed' LIKE '%:{0}') AND
+    (NOT highways.tags?'maxspeed' OR highways.tags->'maxspeed' = '{0}') AND
+    highways.linestring_proj && a3.geom AND
+    ST_Length(ST_Intersection(highways.linestring_proj, a3.geom)) / ST_Length(highways.linestring) > 0.8
+WHERE
+  NOT ST_IsEmpty(a3.geom)
+"""
+
+class Analyser_Osmosis_Highway_Zone(Analyser_Osmosis):
+
+    requires_tables_common = ['highways']
+
+    def __init__(self, config, logger = None):
+        Analyser_Osmosis.__init__(self, config, logger)
+        self.classs[20] = {"item":"1210", "level": 1, "tag": ["highway", "fix:survey"], "desc": T_(u"Probably missing tag zone:maxspeed=XX:%s, according to the neighborhood", 20) }
+        self.classs[30] = {"item":"1210", "level": 1, "tag": ["highway", "fix:survey"], "desc": T_(u"Probably missing tag zone:maxspeed=XX:%s, according to the neighborhood", 30) }
+        self.callback20 = lambda res: {"class":20, "data":[self.way_full, self.positionAsText]}
+        self.callback30 = lambda res: {"class":30, "data":[self.way_full, self.positionAsText]}
+
+    def analyser_osmosis_common(self):
+        self.run(sql10.format('20'))
+        self.run(sql11)
+        self.run(sql12)
+        self.run(sql13)
+        self.run(sql14.format('20'), self.callback20)
+
+        self.run(sql10.format('30'))
+        self.run(sql11)
+        self.run(sql12)
+        self.run(sql13)
+        self.run(sql14.format('30'), self.callback30)

osmose_config.py

@@ -242,6 +242,7 @@ def __init__(self, region, polygon_id=None, country_code="FR", proj=2154, analys
         self.analyser["osmosis_fantoir"] = "xxx"
         self.analyser["osmosis_highway_motorway"] = "xxx"
         self.analyser["merge_traffic_signs"] = "xxx"
+        self.analyser["osmosis_highway_zone"] = "xxx"
 
 france_region("alsace", 8636, "FR-A")
 france_region("aquitaine", 8637, "FR-B")
@@ -957,6 +958,7 @@ def __init__(self, province, polygon_id=None, country_code="DE", proj=32632, ana
         default_country.__init__(self, "europe", "germany/" + province, polygon_id, analyser_options,
                                     download_repo, download_country)
         self.analyser["merge_traffic_signs"] = "xxx"
+        self.analyser["osmosis_highway_zone"] = "xxx"
 
 de_state("baden-wuerttemberg", 62611, "DE-BW")
 #de_state("bayern", 2145268, "DE-BY")