2010/5/14 Vincent Pottier <vpott...@gmail.com> > Le 14/05/2010 01:06, Pieren a écrit : > > 2010/5/13 Vincent Pottier <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 > ) > > Tu n'as pas besoin d'écrire P.way && c.way car ST_intersects utilise automatiquement l'index. Ca permet d'avoir une requête plus propre. Mais c'est vraiment un pinaillage de ma part.
> 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 > A moins d'utiliser le paramètre -l dans osm2pgsql, les données seront converties en projection mercator (Google). Tu n'auras donc pas de latitude, longitude. Emilie Laffray
_______________________________________________ Talk-fr mailing list Talk-fr@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-fr
