Le portail francophone de la géomatique

C'est vrai que mon explication n'est pas très claire...

Voici un schéma joint.

Je pars donc, pour simplifier, de 2 couches (rouge et verte) avec des objets qui se superposent au sein d'une même couche (rouge) et entre 2 couches.
Je souhaite découper tous les objets entre eux et obtenir le résultat dans une nouvelle couche (bleue)

De plus, les 2 couches de départ ont la même structure avec un attribut de type entier "Note". Je souhaite aussi pour chaque objet de la couche bleue obtenir la "Note" maximale des objets se superposant.

Merci d'avance.

Dernière modification par Lsam (Thu 17 April 2014 13:19)

Oui dans ce cas, polygonize peut convenir: tout dépend s'il y a des trous dans les polygones initiaux et ce qu'il faut faire de ces trous.

Nicolas

Oui pardon: j'ai oublié une étape dans ma description...

st_polygonize reconstruit les polygones élémentaires qu'il peut si le réseau de ligne est connecté ("nodé"), ce qui n'est pas le cas par défaut avec les lignes extraites des polygones.

Avec Postgis, facile, il existe st_node pour cela: extraire les frontières des polygones en lignes simples, les regrouper dans une collection, noder les lignes de la collection et passer ces lignes à st_polygonize:

pg1, pg2 et pg2 sont 3 tables de polygones de test (cf. images jointes):

with bnd as (
    select (st_dump(st_boundary(geom))).geom as geom from pg1 UNION ALL
    select (st_dump(st_boundary(geom))).geom as geom from pg2 UNION ALL
    select (st_dump(st_boundary(geom))).geom as geom from pg3
), col as (
    select st_collect(geom) as geom from bnd
) select st_polygonize(st_node(geom)) as geom 
    from col;

Ensuite identifier les trous à supprimer: le plus simple est d'avoir une couche de tous les pg des tables et de supprimer tous les polygones de la requete précédentes sils sont contenus dans un trou et n'intersectent pas d'autres polygones.

Avec Spatialite, qui ne possède pas la fonction st_node (enfin, j'ai regardé rapidement), ca doit etre possible de noder les lignes avec la fonction st_union, qui a cette propriété lorsqu'elle travaille sur des lignes.

La difficulté supplémentaire vient du fait que Spatialite ne permet pas, en SQL, "d'exploser" un multiobjet en objets simples (st_dump dans postgis), or ca serait bien pratique dans ce cas. De ce que j'ai lu sur la mailing list, il faut faire ces dump dans l'interface graphique ou en ligne de commande spatialite ( https://groups.google.com/forum/#!msg/s … CRb-mibcJ)

Je connais moins spatialite, mais ca se teste...

Nicolas

Dernière modification par Nicolas Ribot (Thu 17 April 2014 17:48)