Le portail francophone de la géomatique

Bonjour,

je bloque sur une requête:
j'ai 2 couches, toutes 2 de type multipolygone (1 EPCI et 1 couche occupation du sol corine land cover). Je souhaite découper la couche CLC au niveau des limites de la couche EPCI. Je parviens bien à découper la couche, mais il y un bug sur une partie des objets à cheval sur la limite de l'EPCI. Certains sont ignorés par la requête et ne sont pas découpés.
Je ne sais pas si je suis très clair... Un capture est souvent plus efficace que de longs discours...

DROP TABLE IF EXISTS "10_P024_theix".ocs_clc_test;
CREATE TABLE "10_P024_theix".ocs_clc_test AS
   SELECT s.* 
FROM 
 "10_P024_theix".epci r,
 "10_P024_theix".ocs_clc s
   WHERE s.the_geom && r.the_geom 
   AND ST_Intersects(s.the_geom, r.the_geom);

Si quelqu'un a une idée...

Merci

########
EDIT

J'ai testé une autre requête, mais le résultat est identique à la première :

 DROP TABLE IF EXISTS "10_P024_theix".ocs_clc_3;
CREATE TABLE "10_P024_theix".ocs_clc_3 AS
SELECT s.gid, (ST_Intersection(r.the_geom, s.the_geom)) As the_geom
FROM "10_P024_theix".epci r
    INNER JOIN "10_P024_theix".ocs_clc s
    ON ST_intersects(r.the_geom, s.the_geom)
    WHERE ST_intersects(r.the_geom, s.the_geom)
;

Après quelques tests, il apparaît que ce sont les polygones devant être découpées en plusieurs entités qui posent problème (cf 2e pièce jointe)

Dernière modification par matieu_dumo (Wed 16 March 2011 09:55)