Le 14/05/2010 01:06, Pieren a écrit :
2010/5/13 Vincent Pottier <vpott...@gmail.com <mailto:vpott...@gmail.com>>

Je pense qu'il suffit d'ajouter une ligne dans le genre:
node,way   description     text         polygon
En effet, la colonne 'description' est ajoutée. Mais j'ai mis
node,way   description     text         linear
pour l'import de repères géodésiques.

Merci à Vincent, Émilie, Pieren. J'ai une base spacialisée qui tourne ! Et je commence à faire tourner des requêtes (il y a un an, ça m'impressionnait) :

SELECT id, lat, lon FROM geodesy_nodes
 WHERE id in (
 select DISTINCT ON (way) osm_id FROM geodesy_point
 WHERE description ILIKE '%clocher%'
  or description ILIKE '%croix%'
  or description ILIKE '%église%'
  or description ILIKE '%château%'
except
select p.osm_id
from
(select * from geodesy_point) p,
(select * from planet_osm_polygon where building is not null) c
where p.way && c.way and st_intersects(p.way,c.way) is true
 )

Il y a à peu près 40 000 repères (ou groupes de repères à la verticale) qui attendent un bâtiment dessous (église, château d'eau) !
Ce qui m'étonne, c'est la forme de lat et lon :
670596584;5333352.79653301;693719.667266716
--
FrViPofm
_______________________________________________
Talk-fr mailing list
Talk-fr@openstreetmap.org
http://lists.openstreetmap.org/listinfo/talk-fr

Répondre à