[OSM-dev-fr] [SQL] utilisation de ST_Difference

[Nicolas Dumoulin]

Le mercredi 29 décembre 2010 16:37:17 Vincent de Chateau-Thierry, vous avez 
écrit :
> Une piste pour voir : si dans ta concaténation de points (ta variable
> 'points') tu mets les valeurs de bbox que tu cites avant, il s'agit de
> longitudes & latitudes, alors que dans ST_GeomFromText tu indiques que ce
> sont non pas des lon & lat, mais des mètres dans la projection "Google"
> (900913). Si c'est le cas, essaies plutôt avec dans ton ST_GeomFromText
> "4326" plutôt que "900913", pour signifier que ce sont des longitudes et
> latitudes.

Merci de ta réponse, mais j'obtiens alors l'erreur suivante :
RuntimeError: PSQL error:
ERREUR:  Operation on mixed SRID geometries
Full sql was: 'SELECT AsBinary("way") AS geom from (select 
ST_DIFFERENCE(ST_GeomFromText('POINT(3.098 45.739 3.1543 45.764)', 4326),way) 
as way from planet_osm_polygon where name='Aubière') as roads WHERE "way" && 
SetSRID('BOX3D(344867.7824775615 5738534.792254916,351135.0698092228 
5742697.566128345)'::box3d, 900913)'

Pour info, j'ajoute mes couches postgis avec le srs "+proj=latlong 
+datum=WGS84" (copié-collé des exemples).
-- 
Nicolas Dumoulin
http://wiki.openstreetmap.org/wiki/User:NicolasDumoulin