Le portail francophone de la géomatique

bonjour,

Un anneau est bien un trou pour vous?

Je doute qu'une enveloppe convexe soit dans tous les cas la solution, si le polygone est concave le résultat ne sera surement pas celui attendu.
vos pouvez utiliser un buffer de 0 éventuellement, sinon plutot une combinaison de st_exteriorring et st_makepolygon.

Dernière modification par vincentp (Wed 20 February 2013 01:44)

En effet un anneau est bien un trou, mais c'est le terme de QGIS. J'ai déjà essayé le buffer à 0 mais cela n'efface pas les "véritables anneaux" (exemple : mer intérieure, lac...). En faite j'ai une couche mondiale des bassins hydrographiques et je voudrais la simplifier au maximum (enlever les petites îles et les trous à l'intérieur des polygones)

Ok, merci beaucoup je testerai cette méthode et je vous tiendrai au courant

[EDIT] Encore merci pour l'aide, ça fonctionne nickel.

Dernière modification par gaetanpru (Thu 21 February 2013 10:40)

Bonjour,

Vu le message d'erreur renvoyé, je vérifierai que ton champ geom est bien de type POLYGON.
Que renvoie la requete suivante ?

SELECT ST_GeometryType(geom) FROM public.ma_table

Attention, ST_ExteriorRing ne fonctionne pas avec les MULTIPOLYGON : http://postgis.net/docs/manual-2.0/ST_ExteriorRing.html

Bonjour

Essayez sous Postgis la requête suivante

ST_GeometryType((st_dump(multigeom)).geom))

Pour vérifier que votre multi polygone est bien "explosé" en polygones.

J'ai testé ta méthode Vianney,
La première requête me renvoi bien 0 : multi_geom bigint 0.
Cela signifie donc bien que je n'ai pas d'entité avec plusieurs parties.
Par contre la deuxième me retourne : La requête a été exécutée avec succès en 515 ms, mais ne renvoie aucun résultat.

En fait j'ai ré-importé ma table dans PostGIS en cochant l'option Créer des géométries simples au lieu de multiparties et ça marche.

Néanmoins, quand je retourne à mon problème initial (remplir les trous de mes polygones), cette requête ne change rien :

SELECT st_makepolygon(st_exteriorring(geom))
FROM public.ma_table
;

Dernière modification par Xavier Julian (Fri 29 January 2016 14:32)

Voilà le truc évident à côté du quel je passais, merci !
Excusez moi, je suis autodidacte sur PostGIS, donc il me manque certains éléments basiques.

Maintenant, comment ajouter une condition pour que la création de polygone ne se fasse que sur les trous inférieurs à une certaine taille ?
(j'essaye de corriger des erreurs de topologie, mais j'ai aussi des polygones avec des anneaux que je veux garder).

Bon début de semaine à tous !

Bonjour,

st_area(geom)<20

renvoie un boolean (vrai ou faux).



il faut écrire

INSERT INTO monschema.matable SELECT st_makepolygon(stexteriorring(geom)) FROM public.matabledepart WHERE st_area(geom)<20;

Bonjour,

Je reviens sur ce post un an plus tard car les astuces qui y sont données m'ont beaucoup aidée mais il n'y a pas eu de suite concernant le remplissage des anneaux inférieurs à tant de m².

En effet, ChristopheV, vous nous proposez la requête suivante :

INSERT INTO monschema.matable SELECT st_makepolygon(stexteriorring(geom)) FROM public.matabledepart WHERE st_area(geom)<20;

éventuellement via un WITH.

Mais le st_area(geom)<20 ne va bien chercher que les géométries non nulles? Il faudrait donc que les anneaux soient "remplis" pour qu'on puisse connaitre leur surface?

En effet, si j'ai bien compris ce que fait le st_makepolygon(st_exteriorring((st_dump(geom)).geom)), on cherche à sélectionner le périmètre extérieur de nos polygones, périmètre qu'on remplit ensuite entièrement sans se soucier de savoir si on souhaite conserver certaines anneaux.

Pour ma part, j'ai cherché dans ma condition non pas à sélectionner les géométries de ma table de départ qui font moins de x m², mais à remplir les trous de ma couche de départ via un st_interiorringN(geom,n), à sélectionner ceux qui font moins de x m² et à redécouper ma couche remplie via st_exteriorring avec ces nouveaux polygones obtenus via st_interiorring... Le problème avec st_interiorringN(geom,n), c'est ce fameux n qui consiste à désigner quels trous on souhaite remplir. Or moi, je souhaite tous les remplir sans cette condition n!

Est-ce que je me complique la vie et, surtout, est-ce que ma méthode est applicable?

En vous remerciant par avance pour votre aide!

Bien cordialement,

Lucie D.

Dernière modification par Lucie D. (Wed 03 May 2017 15:19)

Bonjour,

Comme le dit Alain, la réponse a été apportée sur notre liste interne, mais il ne s'agit pas tout à fait de la requête ci-dessus (en tout cas pas pour ce que je souhaitais faire).

Voici comment je procède (merci aux personnes qui m'ont débloquée!) :

1°/ requête qui me permet d'obtenir un objet sans trous :

drop table if exists tache_sans_trous;
create temp table tache_sans_trous (id bigint);

SELECT AddGeometryColumn ('tache_sans_trous','geom',2154,'POLYGON',2);

INSERT INTO tache_sans_trous
SELECT id, st_makepolygon(st_exteriorring(st_geometryn(geom,1))) as geom
FROM tache_avec_trous;

2°/ requête qui me permet de créer une couche des anneaux supérieurs à 1000 m² (correspond aux anneaux que je souhaite conserver)

DROP TABLE IF EXISTS anneaux_sup_1000;
with tache_anneaux_poly as(
select row_number() over () as id, st_geometryN(geom,1) AS geom
from tache_avec_trous),

decompte_trou as (
select id, geom, ST_NumInteriorRings(geom) as nb_trou
from tache_anneaux_poly),

creation_des_anneaux as (
select id, n as num_trou, st_makepolygon(st_interiorringn(geom,n)) as geom
from decompte_trou
CROSS JOIN generate_series(1,nb_trou) as n
where nb_trou>0)

select *
into temp anneaux_sup_1000
from creation_des_anneaux
where st_area(geom)>1000;

CREATE INDEX tache_sans_trous_geom_gist ON tache_sans_trous USING gist(geom);
CREATE INDEX anneaux_sup_1000_geom_gist ON anneaux_sup_1000 USING gist(geom);

3°/ découpage de la tache_sans_trous avec anneaux_sup_1000

DROP TABLE IF EXISTS tache_anneaux_sup1000;
CREATE TABLE tache_anneaux_sup1000
AS
SELECT st_multi(st_collectionextract(st_forcecollection(st_makevalid(st_difference((st_multi(st_collectionextract(st_forcecollection(st_makevalid(st_union(t1.geom))),3))),(st_multi(st_collectionextract(st_forcecollection(st_makevalid(st_union(t2.geom))),3)))))),3))::geometry (MultiPolygon,2154) as geom
FROM tache_sans_trous as t1, anneaux_sup_1000 as t2
WHERE st_intersects(t1.geom,t2.geom);

4°/ insertion des derniers polygones qui ne sont pas concernés par le st_intersects pour obtenir l'ensemble de la tache

INSERT INTO tache_anneaux_sup1000
select st_multi(t1.geom)
from tache_sans_trous as t1, tache_anneaux_sup1000 as t2
where st_intersects(t1.geom,t2.geom) IS FALSE;

Dernière modification par Lucie D. (Wed 10 May 2017 13:50)