Bonjour,
J'ai fait une requête sur une base de type osmisis pour obtenir les vides
entre les limites communales. Ils sont issues pour la plus part de problèmes
du cadastre lui même ou de ratés à de l'outil d'aide à l'import.
Voilà le résultat :
http://f.rodrigo.free.fr/tmp/osm-map/trous.html
La mise à jour n'est pas automatique.
Pour ceux qui veulent regarder sous le capot, la requête est en PJ.
Fred
SELECT
ST_AsKML(ST_Envelope(geom)),
ST_NPoints(geom)
FROM
(SELECT
(ST_Dump(ST_Polygonize(geom))).geom AS geom
FROM
(SELECT
id
FROM
ways
JOIN relation_members ON
id=member_id AND
member_type='W'
JOIN relation_tags ON
relation_tags.relation_id=relation_members.relation_id AND
relation_tags.k='admin_level' AND relation_tags.v='8'
GROUP BY
id
HAVING
COUNT(id)=1
) AS foo
JOIN way_geometry ON
id=way_id
AND NOT ST_IsClosed(geom) -- retire les polygones (îles et communes isolés)
) AS bar
WHERE
ST_NPoints(geom) < 100 -- Valeur exp. determiné sur l'Aquitaine pour ne pas avoir de faux positifs
;
_______________________________________________
Talk-fr mailing list
[email protected]
http://lists.openstreetmap.org/listinfo/talk-fr
