Le portail francophone de la géomatique

Bonsoir,

Qgis aime bien que les couches postgis possèdent un "id" pour chaque enregistrement, c'est peut-être une piste.

Bon week-end

Bonsoir,

Une suggestion :

La ou les tables en question sont bien référencées dans la table la vue geometry_columns du schema public ? Si les données sont dans un autre schéma assurez vous que cette les tables pointent bien sur le schéma adéquate.

Bon courage et tenez nous au courant.

Dernière modification par Loic_GR (Fri 16 February 2018 22:10)

Bonjour,
Qgis a besoin de connaître la projection (SRID) de la couche et le type de géométrie.
Pour l'indiquer à posteriori (c'est mieux ce de le faire à la création ), il suffit d'utiliser
ALTER TABLE LaTable
  ALTER COLUMN geom TYPE geometry(MULTILINEPOLYGON, 2154)
    USING ST_SetSRID(geom,2154);
2154 peut être remplacé par le n° EPSG souhaité

Attention car St_collect crée des collections qui peuvent mélanger les types de géométrie,ce qui n'est pas supporté par Qgis.
Si c'est le cas, il faudra passer par ST_CollectionExtract pour ne conserver que les Polygones
http://postgis.net/docs/ST_CollectionExtract.html.

D'ailleurs, je me rends compte que :
- Si je définis le type comme MULTIPOLYGON (   ALTER COLUMN geom TYPE geometry(MULTIPOLYGON, 2154)
    USING ST_SetSRID(geom,2154) ), alors le résultat d'un ST_Collect de différentes géométries ne sera pas accepté lors d'un INSERT comme étant une géométrie valable. Par contre, le résultat d'un ST_Union le sera.

- A contrario, si je définis le type comme GeometryCollection (   ALTER COLUMN geom TYPE geometry(GeometryCollection, 2154)
    USING ST_SetSRID(geom,2154) ), alors le résultat d'un ST_Union de différentes géométries ne sera pas accepté lors d'un INSERT comme étant une géométrie valable, mais le résultat d'un ST_Collect le sera.

Correct ?

Il n'y a pas de doc pour tout ça ?
Et il n'y a pas de moyen de spécifier la colonne géométrie comme acceptant tous les sous-types de géométries ?

PS : pour la documentation, je viens de trouver ça ( http://www.postgis.fr/chrome/site/docs/ … tries.html ). Y a t'il plus complet ?

Dernière modification par Blue Monkey (Mon 19 February 2018 12:30)

Bonjour, et merci pour votre réponse.

Finalement, j'ai eu recours au code suivant pour opérer

INSERT INTO schema.matable
    ( définition, horodatage , commentaire , geom )
    VALUES
    (
        'toutes_les_zones', current_timestamp,     'Collection des zones (sous forme d''un unique multipolygone 2D)',
        (
            SELECT ST_Multi(ST_Collect(dumpedGeometries)) AS singleMultiGeometry 
            FROM
            (
                SELECT ST_Force2D((ST_Dump(geom)).geom) AS dumpedGeometries
                FROM donnees.reserve_biologique
            ) AS dumpingGeometries
        )
    ) ;

En effet, j'ai fini par comprendre les points suivants :
- que ST_Collect renvoie une géométrie de type Multi* si on lui passe en paramètres uniquement des géométries simples du même type (ex: renvoie un multipolygone si on lui passe des polygones) ; mais que il renvoie un GeometryCollection dès que on lui passe en paramètre un objet de type multi* ou bien des objets de types différents.
Dans mon cas, mes objets d'entrée étant tous des multipolygones, ST_Collect renvoie une collection de multipolygones.
- ST_Multi() permet alors de transformer alors cette collection de multipolygones en un unique multipolygone.
- Même si mes géométries d'entrée viennent d'une colonne de type geometry(multipolygone), pour alimenter ST_Collect de manière homogène avec uniquement des polygones simples, j'ai donc eu recours à ST_Dump.
- j'ai dû ajouter ST_Force2D car certains des multipolygones étaient en 3D (alors qu'ils n'avaient pas de raison de l'être) et cela me faisait des erreurs.

Si quelqu'un de plus expérimenté avait une critique (constructive) de ce code, ce serait bienvenu.

Je ne trouve toutefois toujours nulle part dans la documentation de PostGIS de liste des sous-types du type geometry.