25 Aug
2023
25 Aug
'23
9:09 a.m.
Message posté par : image95
----------------------------------------
Bonjour,
Avec postgis, je souhaite générer une nouvelle table étant le résultat de
l'intersection d'une table de geometries mixtes et de ma couche cadastre
polygonales. Voici ma tentative. En toute modestie, je pense que cette requête est bonne.
Mais j'ai une erreur en lien à la topologie.
-----------------
Code :
create TABLE schema.resultat as
select
gid, suivi,ouvrage_type, source,code,dr,convention,
"t2.parcelle, "Cigeo_t2.section,
st_intersection(t1.geom, t2.geom ) as geom
from schema.t1
join schema.t2
on st_intersects (t1.geom,t2.geom)
-----------------
-----------------
Code :
ERROR: ERREUR: GEOSIntersects: TopologyException: side location conflict at
872675.09999999998 6823686.2999999998
SQL state: XX000
-----------------
J'ai fait des tests avec ST_MakeValid(geom); Mais je ne trouve pas de solution.
Une personne saurait elle m'éclairer ? Un grand merci.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=362238#p362238
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
25 Aug
25 Aug
10:58 a.m.
Message posté par : Nicolas Ribot
----------------------------------------
Bonjour,
Vous travaillez avec quelle version de Postgis (select postgis_full_version();) ? Les
dernieres versions de geos (3.10+) règlent pas mal de problèmes de topologie.
Pour corriger ces problèmes de topologie, on peut tenter de modifier légèrement les coord
des géométries avec st_snapToGrid, ou st_reducePrecision (nouvelle fonction, plus
puissante je crois) et d'utiliser les nouvelles variantes de st_intersection,
st_union, st_difference, st_symDifference qui prennent un paramètre de précision (le meme
qu'on utilisera pour st_reducePrecision)
La premiere étape est de corriger les géométries invalides dans les tables sources:
-----------------
Code :
update matable set geom = st_makeValid(geom)
where not st_isvalid(geom);
-----------------
Puis de tester si l'utilsation de st_reducePrecision sur les geom corrige le pb,
directement dans la requete:
-----------------
Code :
select gid, suivi,ouvrage_type, source,code,dr,convention,
t2.parcelle, t2.section, st_intersection(st_reducePrecision(t1.geom, 0.01),
st_reducePrecision(t2.geom, 0.01) ) as geom
from schema.t1 join schema.t2 on st_intersects (t1.geom,t2.geom);
-----------------
Si c'est le cas, on peut mettre à jour les geom des tables si nécessaire (ou juste
appliquer le reducePrecision dans la requete de calcul d'intersection)
Nicolas
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=362249#p362249
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
11:22 a.m.
Message posté par : Nicolas Ribot
----------------------------------------
Ok merci,
Ca peut valoir le coup de tester sur une version récente de postgis/geos (postgis 3.3 ou
3.4, geos 3.11 ou plus)
(vous pouvez facilement installer une ou plusieurs versions de pg/postgis sur votre
machine (https://www.postgresql.org/download/) pour tester
Nicolas
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=362254#p362254
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
5:24 p.m.
Message posté par : image95
----------------------------------------
Je viens de faire un test avec une nouvelle version de postgres et postgis.
-----------------
Code :
select PostGIS_Full_Version();
-----------------
-----------------
Code :
POSTGIS="3.4.0 3.4.0" [EXTENSION] PGSQL="150"
GEOS="3.12.0-CAPI-1.18.0" PROJ="8.2.1 NETWORK_ENABLED=OFF
URL_ENDPOINT=https://cdn.proj.org
USER_WRITABLE_DIRECTORY=C:\Windows\ServiceProfiles\NetworkService\AppData\Local/proj"
LIBXML="2.9.14" LIBJSON="0.12" LIBPROTOBUF="1.2.1"
WAGYU="0.5.0 (Internal)"
-----------------
Puis
c
-----------------
Code :
reate TABLE schema.resultat as
select
gid, suivi,ouvrage_type, source,code,dr,convention,
t2.parcelle, t2.section,
st_intersection(t1.geom, t2.geom ) as geom
from schema.t1
join schema.t2
on st_intersects (t1.geom,t2.geom)
-----------------
-----------------
Code :
ERROR: lwgeom_intersection_prec: GEOS Error: TopologyException: side location conflict at
875253.84683544014 6826777.4367088638. This can occur if the input geometry is invalid.
ERREUR: lwgeom_intersection_prec: GEOS Error: TopologyException: side location conflict at
875253.84683544014 6826777.4367088638. This can occur if the input geometry is invalid.
SQL state: XX000
---
update andra_humain."Cigeo_CadastreUF_Po" set geom = st_makeValid(geom)
where not st_isvalid(geom);
-----------------
Puis,
-----------------
Code :
update matable set geom = st_makeValid(geom)
where not st_isvalid(geom);
-----------------
NOTICE: Self-intersection at or near point 845991.14978354971 6839069.5350649348 NOTICE:
Self-intersection at or near point 846524.08181817865 6837960.9654545467 NOTICE:
Self-intersection at or near point 882943.5153532261 6830835.5559041062 NOTICE: Hole lies
outside shell at or near point 884955.69999999995 6831733.5999999996 NOTICE: Too few
points in geometry component at or near point 883230.09999999998 6829270.4000000004 ERROR:
Geometry type (MultiLineString) does not match column type (MultiPolygon) ERREUR: Geometry
type (MultiLineString) does not match column type (MultiPolygon) SQL state: 22023
Et ,
-----------------
Code :
select gid, suivi,ouvrage_type, source,code,dr,convention,
t2.parcelle, t2.section, st_intersection(st_reducePrecision(t1.geom, 0.01),
st_reducePrecision(t2.geom, 0.01) ) as geom
from schema.t1 join schema.t2 on st_intersects (t1.geom,t2.geom);
-----------------
-----------------
Code :
ERROR: lwgeom_reduceprecision: GEOS Error: TopologyException: unable to assign free hole
to a shell at 877338.90000000002 6835526.5 ERREUR: lwgeom_reduceprecision: GEOS Error:
TopologyException: unable to assign free hole to a shell at 877338.90000000002 6835526.5
SQL state: XX000
-----------------
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=362265#p362265
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
11:04 a.m.
Message posté par : image95
----------------------------------------
Bonjour Nicolas,
Merci. Avant de tester votre proposition, je réponds à votre question concernant la
version de postgis :
-----------------
Code :
POSTGIS="3.1.0 3.1.0" [EXTENSION] PGSQL="130"
GEOS="3.9.0-CAPI-1.14.1" PROJ="7.1.1" GDAL="GDAL 3.2.0, released
2020/10/26" LIBXML="2.9.9" LIBJSON="0.12"
LIBPROTOBUF="1.2.1" WAGYU="0.5.0 (Internal)" TOPOLOGY RASTER
-----------------
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=362251#p362251
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
1:43 p.m.
Message posté par : image95
----------------------------------------
Merci. POur info, lorsque je tente votre proposition avec st_makevalid et
l'utilisation de st_reducePrecision :
-----------------
Code :
ERROR: ERREUR: lwgeom_reduceprecision: GEOS Error: TopologyException: unable to assign
free hole to a shell at 877338.90000000002 6835526.5
SQL state: XX000
-----------------
Je vais en attendant essayer de faire un test d'upgrade de postgis.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=362257#p362257
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
7:43 p.m.
Message posté par : Nicolas Ribot
----------------------------------------
Bonsoir,
Peut etre utiliser st_snapToGrid ?
Sinon regarder l'endroit ou ca plante (le message d'erreur donne les coord) et
corriger les data.
Nicolas
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=362267#p362267
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
777
days inactive
777
days old
6 comments
1 participants
participants (1)
-
geobd＠georezo.net