- Geobd - Listes du Forum GeoRezo

[Postgres 12] regex_replace: enlever espaces vides

by geobd＠georezo.net

Message posté par : ckiko13 ---------------------------------------- Bonjour Dans une table, je souhaiterais enlever des espaces vides entre des virgules ou à la fin de la ligne. Pour ce faire, j'utilise cette expression: ----------------- Code : select REGEXP_REPLACE (nom_prop,'/\s+/',' ') from parc_group_proprio; ----------------- Où nom_prop = champ concerné & parc_group_proprio = ma table Par exemple, de cet enregistrement: ----------------- Citation : BONNA/AUGUSTE , ,019 RUE COTE ,26000 VALENCE ----------------- Je souhaite obtenir cela: ----------------- Citation : BONNA/AUGUSTE,019 RUE COTE,26000 VALENCE ----------------- J'essaye plusieurs variantes de mon expression régulière ----------------- Citation : '\\s+' '\\s+' '/ss+/' ----------------- , mais rien ne semble fonctionner. Auriez-vous une idée pour reformuler correctement mon expression? Merci ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=335845#p335845 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 years

1
2
0 0

Postgis union et decoupe de polygones dans une même couche

by geobd＠georezo.net

Message posté par : Tonio456 (antoine.prevors(a)gmail.com) ---------------------------------------- Bonsoir à tous, Je cherche à lancer une requête via postgis pour découper des polygones d'une même couche qui se superposent et donc créer de nouvelles géométries tout en récupérant la valeur de mon champ "nombre" pour en faire la moyenne. Ma table ressemble à cela, j'ai a peu 500 cercles https://ibb.co/kc19b4j Si je prends cet exemple https://ibb.co/MhTvTtr J'ai 3 cercles avec un champ nombre de valeurs 1,2 et3 et je cherche à obtenir comme résultat https://ibb.co/Vpm7H9W c'est à dire 7 polygones qui reprennent la moyenne des attributs des polygones avec lesquels ils se superposent J'espère etre assez clair dans mon explication J'etais parti sur une requete du genre ----------------- Code : SELECT a.label,b.label,st_intersection(a.geom,b.geom), avg (matable.nombre) as moyenne FROM matable as a, matable as b WHERE st_intersects(a.geom,b.geom) AND a.idb.id ----------------- mais je pense m'égarer complètement... si quelqu'un a une piste je suis preneur !:) Je vous remercie d'avance ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=335620#p335620 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 years

1
3
0 0

Utilisation de PostgreSQL comme BD pour projet Lizmap

by geobd＠georezo.net

Message posté par : Francis YAMEOGO ---------------------------------------- Bonjour à tous. J'ai élaboré un projet Webmapping a partie de Lizmap qui est une solution qui permet de générer des cartes interactive identique a votre projet QGIS. A présent je compte utiliser PostgreSQL comme Base de donnée de mon projet. A ce niveau je suis bloqué. Dans la documentation Lizmap , il est dit de créer un lien dans le fichier localConfig.uni.PHP en enregistrant les username, password, hôte de la base de donnée PostgreSQL. Chose faite mais j'arrive toujours pas a créer ce lien. Quelqu'un aurait une expérience en la matière. Bien vouloir me faire savoir. Merci ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=335433#p335433 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 years

1
6
0 0

Combinaison de requetes

by geobd＠georezo.net

Message posté par : staubli ---------------------------------------- Salut tout le monde J'ai trois tables differentes que je voudrais combiner au travers d'un VIEW. Mais j'ai beau tout essayer je n'y arrive pas. J'ai une table Avec mes objets (Points) qui contient leur ID, leur geom et une Denomination Sous forme de Code. Une table Avec toutes les denominations de tous mes objets (tel Code correspond a telle Denomination pour tel objet, exemple SY12=poteau electrique). Une 3e table qui indique des poins kilometriques sur une route. Je veux donc pour chaque objet lui associer sa denomination et son kilometrage (son plus proche voisin). Je suis arrivé d'un coté à associer un objet à son nom et d'un autre coté à associer un objet à son kilometrage mais je n'arrive donc pas à combiner les deux. J'imagine que ca doit etre tout bete mais je n'y arrive pas. Je commence juste à migrer les données vers Postgre et je n'ai pas fait beaucoup de SQL, ce qui explique mes difficultés. Requete associer une Denomination à un Code ----------------- Code : SELECT p."ID", p.geom, p.name, p.angle, p.linie, i.bezeichnung FROM infra_visu._prellbockdist p JOIN infra_visu._symbolindex_bahnplan i ON p.name::text = i.symbolnr::text; ----------------- Requete associer un Point kilometrique ----------------- Code : SELECT DISTINCT ON (_prellbockdist."ID") _prellbockdist."ID", _prellbockdist.name, _prellbockdist.angle, _prellbockdist.linie, _kilometrierung.linie, _kilometrierung.km FROM infra_visu._prellbockdist, infra_visu._kilometrierung WHERE ST_DWithin(_kilometrierung.geom , _prellbockdist.geom, 100) Order by _prellbockdist."ID", ST_Distance (_kilometrierung.geom , _prellbockdist.geom) ----------------- ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=335393#p335393 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 years

1
8
0 0

Pgrouting - Le plus court chemin d'un point A à un point B

by geobd＠georezo.net

Message posté par : preliator ---------------------------------------- Bonjour, Je m'initie depuis peu sur Pgrouting, notamment grâce à la document en ligne. Disposant d'une couche de route (source : OSM), je souhaite calculer le chemin le plus court entre deux points. J'utilise Windows 10, ainsi que postgreSQL et pgAdmin (dernière version). La requête finale fonctionne, mais est malheureusement vide. Voici ma démarche : ----------------- Code : -- Sur QGis, conversion de la couche route de Multilinestring à linestring avec l'outil "De morceaux multiples à morceaux uniques" -- Importation de la couche route (nommée "route_049") dans PostgreSQL -- Création des colonnes Source et Target ALTER TABLE route_049 ADD COLUMN "source" integer; ALTER TABLE route_049 ADD COLUMN "target" integer; --Création que de la topologie (fonctionne bien, et me renvoie "OK") SELECT pgr_createtopology('route_049', 0.001, 'geom', 'id', 'source', 'target'); --Tentative de création d'un trajet le plus court entre deux points (https://zupimages.net/viewer.php?id=20/37/2piv.png --> Les points proviennent de "route_049_vertices_pgr") SELECT * FROM pgr_dijkstra('select id, source, target, st_length(geom) as cost, st_length(geom) as cost FROM route_049', 30, 4, false) --0 ligne retournée ----------------- J'ai du adapter mon code, car les formules que j'entrais n'était pas compatibles avec la doc. Aurais-je loupé quelque chose ? Merci. ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=335270#p335270 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 years

1
3
0 0

Optimisation du ST_DIFFERENCE

by geobd＠georezo.net

Message posté par : hadri45 ---------------------------------------- Bonjour, Je dois réaliser des croisements de couches forestières afin d'enrichir la base de données des Forêts officielle. Pour cette opération, le mieux est d'utiliser la fonction Union d'ArcGIS mais j'ai l'obligation de réaliser cela sous Postgres/PostGIS. Pour se faire, j'ai donc vu qu'il fallait réaliser cela en 2 étapes avec d'abord l'intersection pour les deux couches qui s'intersectent puis avec INSERT INTO, récupérer celle qui ne s'intersecte pas avec le ST_DIFFERENCE. Voici le code où je dois réaliser deux croisements tout en gardant les attributs (code_psg,code_onf et tfv) pour des analyses complémentaires. Ce code me donne le bon résultat mais vu que je dois réaliser ce traitement pour la France entière il me faut l'optimiser le plus possible sachant que sur un département cela prend plus de 2h30 alors que sur ArcGIS le traitement prend 4 minutes. ----------------- Code : -- TABLE DU CROISEMENT PRO / PSG DROP TABLE IF EXISTS surface_unitaire_01.propsg_01; -- 1. Intersection CREATE TABLE surface_unitaire_01.propsg_01 AS SELECT o.code_onf, p.code AS code_psg, ST_SetSRID((ST_Dump(ST_Intersection(o.geom, p.geom))).geom, 932006) AS geom FROM surface_unitaire_01.pro_2018_simple_01 o INNER JOIN surface_unitaire_01.psg_2018_simple_01 p ON (o.geom && p.geom AND ST_Intersects(o.geom, p.geom)); -- Nettoyage des géométries invalides DELETE FROM surface_unitaire_01.propsg_01 WHERE ST_GeometryType(geom) != 'ST_Polygon'; CREATE INDEX sp_idx_propsg_01_geom ON surface_unitaire_01.propsg_01 USING GIST (geom); ANALYZE surface_unitaire_01.propsg_01; -- 2. Différences PRO - PSG et PSG - PRO CREATE TEMP TABLE pro_union AS SELECT code_onf, ST_Union(geom) AS geom FROM surface_unitaire_01.propsg_01 GROUP BY code_onf; CREATE INDEX sp_idx_pro_union_geom ON pro_union USING GIST (geom); ANALYZE pro_union; CREATE TEMP TABLE psg_union AS SELECT code_psg, ST_Union(geom) AS geom FROM surface_unitaire_01.propsg_01 GROUP BY code_psg; CREATE INDEX sp_idx_psg_union_geom ON psg_union USING GIST (geom); ANALYZE psg_union; DROP INDEX surface_unitaire_01.sp_idx_propsg_01_geom; INSERT INTO surface_unitaire_01.propsg_01 SELECT o.code_onf, NULL::VARCHAR(3) AS code_psg, (ST_Dump(ST_Difference(o.geom, ST_SetSRID(COALESCE(t1.geom, 'GEOMETRYCOLLECTION EMPTY'::GEOMETRY), 932006)))).geom AS geom FROM surface_unitaire_01.pro_2018_simple_01 o LEFT JOIN pro_union t1 ON o.code_onf = t1.code_onf UNION SELECT NULL::VARCHAR(3) AS code_onf, p.code AS code_psg, (ST_Dump(ST_Difference(p.geom, ST_SetSRID(COALESCE(t2.geom, 'GEOMETRYCOLLECTION EMPTY'::GEOMETRY), 932006)))).geom AS geom FROM surface_unitaire_01.psg_2018_simple_01 p LEFT JOIN psg_union t2 ON p.code = t2.code_psg; -- Nettoyage des géométries invalides SELECT DISTINCT ST_GeometryType(geom) FROM surface_unitaire_01.propsg_01; SELECT DISTINCT ST_SRID(geom) FROM surface_unitaire_01.propsg_01; SELECT COUNT(*) FROM surface_unitaire_01.propsg_01 WHERE NOT ST_IsValid(geom); DELETE FROM surface_unitaire_01.propsg_01 WHERE ST_GeometryType(geom) != 'ST_Polygon'; CREATE INDEX sp_idx_propsg_01_geom ON surface_unitaire_01.propsg_01 USING GIST (geom); ANALYZE surface_unitaire_01.propsg_01; DROP TABLE pro_union; DROP TABLE psg_union; -- contrôles de surfaces SELECT 'couche onf' AS origine, code_onf, SUM(ST_Area(geom)) AS surface_total FROM surface_unitaire_01.pro_2018_simple_01 GROUP BY origine, code_onf UNION SELECT 'croisement' AS origine, code_onf, SUM(ST_Area(geom)) AS surface_total FROM surface_unitaire_01.propsg_01 GROUP BY origine, code_onf ORDER BY code_onf, origine; SELECT 'couche psg' AS origine, code AS code_psg, SUM(ST_Area(geom)) AS surface_total FROM surface_unitaire_01.psg_2018_simple_01 GROUP BY origine, code_psg UNION SELECT 'croisement' AS origine, code_psg, SUM(ST_Area(geom)) AS surface_total FROM surface_unitaire_01.propsg_01 GROUP BY origine, code_psg ORDER BY code_psg, origine; SELECT DISTINCT st_geometrytype(geom) FROM surface_unitaire_01.bd_foret_01_simple SELECT id_foret, ST_NPoints(geom), ST_MemSize(geom) FROM surface_unitaire_01.bd_foret_01_simple ORDER BY ST_NPoints DESC; SELECT COUNT(*), COUNT(*) FILTER (WHERE ST_MemSize(geom) > 8192) FROM surface_unitaire_01.bd_foret_01_simple; CREATE TABLE surface_unitaire_01.bd_foret_01_simple_sub AS SELECT id_foret, tfv, ST_SubDivide(geom) AS geom FROM surface_unitaire_01.bd_foret_01_simple; CREATE INDEX idx_foret_geom_sub_01 ON surface_unitaire_01.bd_foret_01_simple_sub USING GIST (geom); CREATE TABLE surface_unitaire_01.propsg_01_sub AS SELECT ST_SUBDIVIDE(geom) as geom, code_onf,code_psg FROM surface_unitaire_01.propsg_01 -- TABLE DU CROISEMENT BD_FORET / PRO_PSG DROP TABLE IF EXISTS surface_unitaire_01.final_cedric_01; -- 1. Intersection CREATE TABLE surface_unitaire_01.final_cedric_01 AS SELECT p.code_onf, p.code_psg, b.tfv, ST_SetSRID((ST_Dump(ST_Intersection(ST_Buffer(p.geom, 0.0000001), b.geom))).geom, 932006) AS geom FROM surface_unitaire_01.propsg_01 p INNER JOIN surface_unitaire_01.bd_foret_01_simple_sub b ON (p.geom && b.geom AND ST_Intersects(p.geom, b.geom)); -- Nettoyage des géométries invalides SELECT DISTINCT ST_GeometryType(geom) FROM surface_unitaire_01.final_cedric_01; DELETE FROM surface_unitaire_01.final_cedric_01 WHERE ST_GeometryType(geom) != 'ST_Polygon'; CREATE INDEX sp_idx_final_cedric_01_geom ON surface_unitaire_01.final_cedric_01 USING GIST (geom); ANALYZE surface_unitaire_01.final_cedric_01; -- 2. Différences BD_FORET - PRO_PSG et PRO_PSG - BD_FORET DROP TABLE IF EXISTS pro_psg; CREATE TEMP TABLE pro_psg AS SELECT code_onf, code_psg,geom FROM surface_unitaire_01.final_cedric_01; CREATE INDEX sp_idx_pro_psg_geom ON pro_psg USING GIST (geom); ANALYZE pro_psg; DROP TABLE IF EXISTS tfv; CREATE TEMP TABLE tfv AS SELECT tfv, geom FROM surface_unitaire_01.final_cedric_01; CREATE INDEX sp_idx_tfv_geom ON tfv USING GIST (geom); ANALYZE tfv; DROP INDEX surface_unitaire_01.sp_idx_final_cedric_01_geom; INSERT INTO surface_unitaire_01.final_cedric_01 SELECT p.code_onf, p.code_psg, NULL::VARCHAR(10) AS tfv, (ST_Dump(ST_Difference(p.geom, ST_SetSRID(COALESCE(t1.geom, 'GEOMETRYCOLLECTION EMPTY'::GEOMETRY), 932006)))).geom AS geom FROM surface_unitaire_01.propsg_01_sub p LEFT JOIN pro_psg as t1 ON p.code_onf = t1.code_onf AND p.code_psg = t1.code_psg UNION SELECT NULL::VARCHAR(3) AS code_onf, NULL::VARCHAR(3) AS code_psg, t.tfv, (ST_Dump(ST_Difference(t.geom, ST_SetSRID(COALESCE(t2.geom, 'GEOMETRYCOLLECTION EMPTY'::GEOMETRY), 932006)))).geom AS geom FROM surface_unitaire_01.bd_foret_01_simple_sub t LEFT JOIN tfv t2 ON t.tfv = t2.tfv; CREATE INDEX sp_idx_final_cedric_01_geom ON surface_unitaire_01.final_cedric_01 USING GIST (geom); ANALYZE surface_unitaire_01.propsg_01; DROP TABLE pro_psg_union; DROP TABLE IF EXISTS tfv_union; DROP TABLE tfv; ----------------- J'ai donc essayé de subdiviser ce qui accélère les temps de traitements mais le plus gros problème reste le dernier ST_DIFFERENCE entre la bd_foret et pro_psg qui est toujours aussi long malgré cette subdivision et ne permet pas de valider le code. J'ai également essayé de transformer les géométries en tuile (MapVector Tile) avec ST_ASMVTGeom mais cela est encore plus long. J'ai également essayé de faire sans le ST_UNION ou ST_COLLECT mais cela ne change pas beaucoup le temps de traitement. Ma couche principale est de l'ordre de plus de 25000 entités et les autres couches à croiser sont autour de 2000. Pour la configuration, je suis sous PostgreSQL 12 et PostGIS 3.0. Auriez-vous donc une idée afin d'optimiser le traitement ? Merci, ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=335212#p335212 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 years, 1 month

1
7
0 0

architecture bdd

by geobd＠georezo.net

Message posté par : Ricola62 ---------------------------------------- Bonjour, Je m'interroge concernant l'architecture de ma base de donnée à adopter afin d'avoir le maximum de sécurité. Il me faut réaliser d'un coté une webapp pour la consultation et permettant la génération de document pdf et de l'autre une partie saisie qui enrichira la base de donnée. Je me posais la question d'avoir une base de donnée temporaire (pour la saisie) et basculer les données dans bdd consultation lorsque celles ci sont validées (par l'admninistrateur) via une automatisation chaque début de semaine par exemple, ce qui permettrait d'éviter les coquilles lors de la consultation sur la webapp. J'aimerais avoir vos avis sur la question et quelques suggestions. Merci d'avance ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=335300#p335300 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 years, 1 month

1
0
0 0

Optimisation du ST_DIFFERENCE

by geobd＠georezo.net

Message posté par : hadri45 ---------------------------------------- Bonjour, Je dois réaliser des croisements de couches forestières afin d'enrichir la base de données des Forêts officielle. Pour cette opération, le mieux est d'utiliser la fonction Union d'ArcGIS mais j'ai l'obligation de réaliser cela sous Postgres/PostGIS. Pour se faire, j'ai donc vu qu'il fallait réaliser cela en 2 étapes avec d'abord l'intersection pour les deux couches qui s'intersectent puis avec INSERT INTO, récupérer celle qui ne s'intersecte pas avec le ST_DIFFERENCE. Voici le code où je dois réaliser deux croisements tout en gardant les attributs (code_psg,code_onf et tfv) pour des analyses complémentaires. Ce code me donne le bon résultat mais vu que je dois réaliser ce traitement pour la France entière il me faut l'optimiser le plus possible sachant que sur un département cela prend plus de 2h30 alors que sur ArcGIS le traitement prend 4 minutes. ----------------- Code : -- TABLE DU CROISEMENT PRO / PSG DROP TABLE IF EXISTS surface_unitaire_01.propsg_01; -- 1. Intersection CREATE TABLE surface_unitaire_01.propsg_01 AS SELECT o.code_onf, p.code AS code_psg, ST_SetSRID((ST_Dump(ST_Intersection(o.geom, p.geom))).geom, 932006) AS geom FROM surface_unitaire_01.pro_2018_simple_01 o INNER JOIN surface_unitaire_01.psg_2018_simple_01 p ON (o.geom && p.geom AND ST_Intersects(o.geom, p.geom)); -- Nettoyage des géométries invalides DELETE FROM surface_unitaire_01.propsg_01 WHERE ST_GeometryType(geom) != 'ST_Polygon'; CREATE INDEX sp_idx_propsg_01_geom ON surface_unitaire_01.propsg_01 USING GIST (geom); ANALYZE surface_unitaire_01.propsg_01; -- 2. Différences PRO - PSG et PSG - PRO CREATE TEMP TABLE pro_union AS SELECT code_onf, ST_Union(geom) AS geom FROM surface_unitaire_01.propsg_01 GROUP BY code_onf; CREATE INDEX sp_idx_pro_union_geom ON pro_union USING GIST (geom); ANALYZE pro_union; CREATE TEMP TABLE psg_union AS SELECT code_psg, ST_Union(geom) AS geom FROM surface_unitaire_01.propsg_01 GROUP BY code_psg; CREATE INDEX sp_idx_psg_union_geom ON psg_union USING GIST (geom); ANALYZE psg_union; DROP INDEX surface_unitaire_01.sp_idx_propsg_01_geom; INSERT INTO surface_unitaire_01.propsg_01 SELECT o.code_onf, NULL::VARCHAR(3) AS code_psg, (ST_Dump(ST_Difference(o.geom, ST_SetSRID(COALESCE(t1.geom, 'GEOMETRYCOLLECTION EMPTY'::GEOMETRY), 932006)))).geom AS geom FROM surface_unitaire_01.pro_2018_simple_01 o LEFT JOIN pro_union t1 ON o.code_onf = t1.code_onf UNION SELECT NULL::VARCHAR(3) AS code_onf, p.code AS code_psg, (ST_Dump(ST_Difference(p.geom, ST_SetSRID(COALESCE(t2.geom, 'GEOMETRYCOLLECTION EMPTY'::GEOMETRY), 932006)))).geom AS geom FROM surface_unitaire_01.psg_2018_simple_01 p LEFT JOIN psg_union t2 ON p.code = t2.code_psg; -- Nettoyage des géométries invalides SELECT DISTINCT ST_GeometryType(geom) FROM surface_unitaire_01.propsg_01; SELECT DISTINCT ST_SRID(geom) FROM surface_unitaire_01.propsg_01; SELECT COUNT(*) FROM surface_unitaire_01.propsg_01 WHERE NOT ST_IsValid(geom); DELETE FROM surface_unitaire_01.propsg_01 WHERE ST_GeometryType(geom) != 'ST_Polygon'; CREATE INDEX sp_idx_propsg_01_geom ON surface_unitaire_01.propsg_01 USING GIST (geom); ANALYZE surface_unitaire_01.propsg_01; DROP TABLE pro_union; DROP TABLE psg_union; -- contrôles de surfaces SELECT 'couche onf' AS origine, code_onf, SUM(ST_Area(geom)) AS surface_total FROM surface_unitaire_01.pro_2018_simple_01 GROUP BY origine, code_onf UNION SELECT 'croisement' AS origine, code_onf, SUM(ST_Area(geom)) AS surface_total FROM surface_unitaire_01.propsg_01 GROUP BY origine, code_onf ORDER BY code_onf, origine; SELECT 'couche psg' AS origine, code AS code_psg, SUM(ST_Area(geom)) AS surface_total FROM surface_unitaire_01.psg_2018_simple_01 GROUP BY origine, code_psg UNION SELECT 'croisement' AS origine, code_psg, SUM(ST_Area(geom)) AS surface_total FROM surface_unitaire_01.propsg_01 GROUP BY origine, code_psg ORDER BY code_psg, origine; SELECT DISTINCT st_geometrytype(geom) FROM surface_unitaire_01.bd_foret_01_simple SELECT id_foret, ST_NPoints(geom), ST_MemSize(geom) FROM surface_unitaire_01.bd_foret_01_simple ORDER BY ST_NPoints DESC; SELECT COUNT(*), COUNT(*) FILTER (WHERE ST_MemSize(geom) > 8192) FROM surface_unitaire_01.bd_foret_01_simple; CREATE TABLE surface_unitaire_01.bd_foret_01_simple_sub AS SELECT id_foret, tfv, ST_SubDivide(geom) AS geom FROM surface_unitaire_01.bd_foret_01_simple; CREATE INDEX idx_foret_geom_sub_01 ON surface_unitaire_01.bd_foret_01_simple_sub USING GIST (geom); CREATE TABLE surface_unitaire_01.propsg_01_sub AS SELECT ST_SUBDIVIDE(geom) as geom, code_onf,code_psg FROM surface_unitaire_01.propsg_01 -- TABLE DU CROISEMENT BD_FORET / PRO_PSG DROP TABLE IF EXISTS surface_unitaire_01.final_cedric_01; -- 1. Intersection CREATE TABLE surface_unitaire_01.final_cedric_01 AS SELECT p.code_onf, p.code_psg, b.tfv, ST_SetSRID((ST_Dump(ST_Intersection(ST_Buffer(p.geom, 0.0000001), b.geom))).geom, 932006) AS geom FROM surface_unitaire_01.propsg_01 p INNER JOIN surface_unitaire_01.bd_foret_01_simple_sub b ON (p.geom && b.geom AND ST_Intersects(p.geom, b.geom)); -- Nettoyage des géométries invalides SELECT DISTINCT ST_GeometryType(geom) FROM surface_unitaire_01.final_cedric_01; DELETE FROM surface_unitaire_01.final_cedric_01 WHERE ST_GeometryType(geom) != 'ST_Polygon'; CREATE INDEX sp_idx_final_cedric_01_geom ON surface_unitaire_01.final_cedric_01 USING GIST (geom); ANALYZE surface_unitaire_01.final_cedric_01; -- 2. Différences BD_FORET - PRO_PSG et PRO_PSG - BD_FORET DROP TABLE IF EXISTS pro_psg; CREATE TEMP TABLE pro_psg AS SELECT code_onf, code_psg,geom FROM surface_unitaire_01.final_cedric_01; CREATE INDEX sp_idx_pro_psg_geom ON pro_psg USING GIST (geom); ANALYZE pro_psg; DROP TABLE IF EXISTS tfv; CREATE TEMP TABLE tfv AS SELECT tfv, geom FROM surface_unitaire_01.final_cedric_01; CREATE INDEX sp_idx_tfv_geom ON tfv USING GIST (geom); ANALYZE tfv; DROP INDEX surface_unitaire_01.sp_idx_final_cedric_01_geom; INSERT INTO surface_unitaire_01.final_cedric_01 SELECT p.code_onf, p.code_psg, NULL::VARCHAR(10) AS tfv, (ST_Dump(ST_Difference(p.geom, ST_SetSRID(COALESCE(t1.geom, 'GEOMETRYCOLLECTION EMPTY'::GEOMETRY), 932006)))).geom AS geom FROM surface_unitaire_01.propsg_01_sub p LEFT JOIN pro_psg as t1 ON p.code_onf = t1.code_onf AND p.code_psg = t1.code_psg UNION SELECT NULL::VARCHAR(3) AS code_onf, NULL::VARCHAR(3) AS code_psg, t.tfv, (ST_Dump(ST_Difference(t.geom, ST_SetSRID(COALESCE(t2.geom, 'GEOMETRYCOLLECTION EMPTY'::GEOMETRY), 932006)))).geom AS geom FROM surface_unitaire_01.bd_foret_01_simple_sub t LEFT JOIN tfv t2 ON t.tfv = t2.tfv; CREATE INDEX sp_idx_final_cedric_01_geom ON surface_unitaire_01.final_cedric_01 USING GIST (geom); ANALYZE surface_unitaire_01.propsg_01; DROP TABLE pro_psg_union; DROP TABLE IF EXISTS tfv_union; DROP TABLE tfv; ----------------- J'ai donc essayé de subdiviser ce qui accélère les temps de traitements mais le plus gros problème reste le dernier ST_DIFFERENCE entre la bd_foret et pro_psg qui est toujours aussi long malgré cette subdivision et ne permet pas de valider le code. J'ai également essayé de transformer les géométries en tuile (MapVector Tile) avec ST_ASMVTGeom mais cela est encore plus long. J'ai également essayé de faire sans le ST_UNION ou ST_COLLECT mais cela ne change pas beaucoup le temps de traitement. Ma couche principale est de l'ordre de plus de 25000 entités et les autres couches à croiser sont autour de 2000. Pour la configuration, je suis sous PostgreSQL 12 et PostGIS 3.0. Auriez-vous donc une idée afin d'optimiser le traitement ? Merci, ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=335211#p335211 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 years, 1 month

1
0
0 0

trigger update apres Editeur d'attributs portal

by geobd＠georezo.net

Message posté par : benoit111 ---------------------------------------- Bonjour, sgbd : postgressql 9.2 Je mets à jour une table 'ESSAI' via portal et avec l'outil Editeur d'attributs (https://enterprise.arcgis.com/fr/portal/latest/use/widget-batch-attribute-e…) J'ai créer un trigger sur cette table ----------------- Code : CREATE TRIGGER essai_update AFTER INSERT OR UPDATE OR DELETE ON ESSAI FOR EACH ROW EXECUTE PROCEDURE fonction_update(); ----------------- qui appelle une fonction pour remplir une table report à chaque update. Si je réalise un update sur ESSAI à partir de pg_admin, le mécanisme fonctionne. Par contre sur portal avec le widget Editeur d'attributs impossible de remplir la table essai. Je ne trouve pas d'info la dessus alors si quelqu'un à eut le même problème et peut me donner quelques pistes ? merci par avance ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=335063#p335063 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 years, 1 month

1
1
0 0

Problème de mot de passe pour l'import d'un raster sur PostgreSQL

by geobd＠georezo.net

Message posté par : Killian Jaumain ---------------------------------------- Bonjour, J'ai PostgreSQL12.3 avec PostGIS 3.0.2 sur Windows. J'essaie d'importer des raster et des géo package sur PostgreSQL mais quand j’écris le code : ----------------- Code : C:\Program Files\PostgreSQL\12\bin>raster2pgsql -s 2154 -C -M C:\Users\jaumaik24d\Desktop\scan_ign\scan.gpkg | psql -d Generale -U postgres -h localhost -p 5432 ----------------- J'obtiens : ----------------- Code : Processing 1/1: C:\Users\jaumaik24d\Desktop\scan_ign\scan.gpkg Mot de passe pour l'utilisateur postgres : ----------------- Et il est impossible d'écrire le mot de passe même par copier coller. Pouvez vous m'aider à trouver ou se situe le problème ? Cordialement jau19 ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=334995#p334995 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 years, 1 month

1
3
0 0

2025

2024

2023

2022

2021

2020

2019

Geobd