On 18/05/2010 22:56, Vincent Pottier wrote:
> Bonsoir,
> Je commence à bricoler sur postgres, mais j'ai quelques soucis.
> * Je suis très vexé. Le Doubs n'est pas importé dans la base
> J'ai fais plusieurs imports de l'extract France de Geofabric via 
> osm2pgsql. Je ne trouve pas le Doubs avec la requête :
> SELECT ref, name from france_polygon WHERE admin_level='6' ORDER BY ref;
> alors qu'il semble bien formé dans osmose
> http://www.openstreetmap.org/browse/relation/7462
>
> osm2pgsql v. 0.66
> postgres v. 8.4.3
> Bon, au dela de mon sentiment personnel finalement peu affecté, ce qui 
> m'inquiète c'est que du coup l'import que j'ai n'est pas fiable...
>
>   
L'import Geofabrik ne complete pas les ways ou les relations qui ont des
points manquants. Il se peut que s'il manque un point le polygone ne
soit pas cree. Mais la je suis sure que Sylvain saura repondre mieux que
moi.

> * par ailleurs, une requête pour obtenir les communes sans voirie me 
> retourne des faux négatifs :
> SELECT      c."ref:INSEE", c.name, ST_AsText(ST_Centroid(c.way))
> FROM        france_polygon AS c
>                  LEFT OUTER JOIN france_roads AS r
>                      ON ST_Crosses(c.way, r.way)
>                          AND ( r.highway='residential'
>                             OR r.highway='unclassified'
>                             OR r.highway='tertiary'
>                             OR r.highway='secondary'
>                             )
> WHERE c.admin_level='8'
>      AND r.osm_id IS NULL
> ORDER BY c."ref:INSEE"
>
> dont Beure : http://osm.org/go/0CUkJeSY--
> ou Deluz : http://osm.org/go/0CU4hPy4-
> D'après le schéma sur 
> http://postgis.refractions.net/documentation/manual-svn/ST_Crosses.html
> et notamment le troisième schéma, c'est bien cette fonction qui 
> permettrait de trouver les communes ayant de la voirie qui pénètre.
> Mais déjà deux faux négatifs, c'est trop.
>
>   
Pour debugger, j'aurais tendance a dire que ca serait bien que tu
recuperes tous les id pour un de tes faux negatifs. Ensuite, tu affiches
l'id des routes que Postgis pense faire une intersection avec. Utilises
ensuite un openjump ou QGis pour visualiser ces ids et voir pourquoi il
te retourne ces faux positifs.
Je pense que c'est le plus simple.
Peut etre que ST_Intersects est plus interessant.

> * Enfin, j'ai créé des index gist(machin.way) sur les tables. Le gain de 
> vitesse n'est pas fulgurant et j'ai toujours le fameux message : 
> NOTICE:  LWGEOM_gist_joinsel called with incorrect join type
>   
Le notice n'est pas un probleme, tu peux l'ignorer.

Emilie Laffray

Attachment: signature.asc
Description: OpenPGP digital signature

_______________________________________________
Talk-fr mailing list
Talk-fr@openstreetmap.org
http://lists.openstreetmap.org/listinfo/talk-fr

Répondre à