October 2020 - Geobd - Listes du Forum GeoRezo

by geobd＠georezo.net

Message posté par : Nicolas Ribot ---------------------------------------- Bonjour, Vous auriez un exemple avec des données ? Je ne suis pas sur de comprendre le problème. Avec les fonctions de référencement linéaire, vous pouvez découper une linestring tous les x mètres sans trop de problèmes. Nicolas ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=333858#p333858 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

3 years, 7 months

1
2
0 0

Creation Trigger sur PGAdmin

by geobd＠georezo.net

Message posté par : Solène.Cottet (so.cottet(a)live.fr) ---------------------------------------- Bonjour, J'essaie de créer mon premier trigger sur l'interface PGAdmin. Ma BDD possède 2 couches de polygones, observation et site. Les observations sont incluses dans des sites. Je cherche à réaliser un trigger qui, lorsqu'une observation est ajoutée, ajoute automatiquement le nom du site dans laquelle elle se trouve, à une colonne déjà existante (nom_site, jointe précédemment). Je tâtonne un peu car je n'ai jamais réalisé ce genre de fonction. Voici le code que j'ai écris : ----------------- Code : CREATE TRIGGER update_observation BEFORE INSERT OR UPDATE ON observation FOR EACH ROW EXECUTE PROCEDURE update_nom_site_observation() CREATE OR REPLACE FUNCTION update_nom_site_observation() RETURNS TRIGGER AS BEGIN SELECT nom_site FROM site ; ON ST_Overlaps(observation.geom, site.geom) ; UPDATE observation(nom_site); END; ----------------- Je n'arrive pas à déceler où se trouve mon erreur, qui doit être une erreur de syntaxe. L'execution me renvoit l'erreur suivante : ERROR: ERREUR: erreur de syntaxe sur ou près de « CREATE » LINE 4: CREATE OR REPLACE FUNCTION update_observation() Merci d'avance pour votre aide ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=334974#p334974 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

3 years, 8 months

1
9
0 0

PG_routing : pgr_drivingdistance n'existe pas

by geobd＠georezo.net

Message posté par : Theos2000 ---------------------------------------- Bonjour, J'ai souhaité me replonger dans un travail effectué sur Pg_routing il y a 2 ans. Or aprés plusieurs sauvegardes et restaurations de cette base lorsque je lance mes requetes avec PG_routing ca me dit : "la fonction pgr_drivingdistance(unknown, integer, numeric, boolean, boolean) n'existe pas....." Est ce que quelqu'un a déja vu ca ? j'ai vérifié que le shéma topology et que les couches n'avaient pas été modifié et que les couches mobilisés pour le calcul d'isochrone existaient bien, ce qui est le cas... Du coup je comprend pas trop pourquoi un tel message.... Merci ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=334908#p334908 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

4 years

1
4
0 0

Couper une ligne tous les X mètres

by geobd＠georezo.net

Message posté par : preliator ---------------------------------------- Bonjour, Disposant d'une couche multilinestring représentant les routes d'un département, je souhaiterais couper cette couche en lignes de 500m. Voici un aperçu de ma couche de routes : https://zupimages.net/viewer.php?id=20/16/eh9e.png Après quelques recherche sur internet, je suis parti sur cette démarche : - Regrouper ma couche de route en un seul bloc (j'ai trouvé cela cohérent car la couche de base est constitués de plusieurs lignes de tailles variables) - Transformer ce bloc en linestring pour la requête suivante - Créer des points tous les 500m le long de ma linestring avec st_lineinterpolatepoints (pas possible de créer les points sur ma couche regroupée, car cette fonction n'accepte que du linestring) - Couper la couche regroupée sur cette couche de points. ----------------- Code : -- REGROUPER MA COUCHE DE ROUTE EN UN SEUL BLOC MULTILINESTRING drop table if exists ligne_regroup; create table ligne_regroup as select St_linemerge(ST_Union(geom)) as geom from public.route_echantillon; -- TRANSFORMER LE BLOC EN LINESTRING drop table if exists ligne_regroup_linestring; CREATE TABLE ligne_regroup_linestring AS SELECT (ST_Dump(ligne_regroup.geom)).geom::Geometry(LineString,2154) AS geom FROM ligne_regroup; -- CREER UN IDENTIFIANT UNIQUE SUR LA LINESTRING alter table ligne_regroup_linestring add column id serial; -- CREER LES POINT TOUS LES 500 METRES drop table if exists point_500; create table point_500 as with tmp as( select * from ligne_regroup_linestring where st_length(ligne_regroup_linestring.geom) > 500) SELECT id, st_lineinterpolatepoints((tmp.geom),500/st_length(tmp.geom)) as geom from tmp; -- ON COUPE LA ROUTE REGROUPEE SUR LES POINTS drop table if exists route_cut_500_metres; create table route_cut_500_metres as select ST_SNAP((ST_DUMP(st_difference(ligne_regroup.geom,point))).geom,all_point,0.1) as geom from (select ST_Multi(ST_Union(st_expand(point_500.geom, 0.05))) as point from point_500 ) as t1, ligne_regroup,(select ST_MULTI(ST_COLLECT(point_500.geom))as all_point from point_500) as t2; ----------------- Malheureusement, j'obtiens des résultats aberrants. Les lignes semblent bien se découper, mais il manque de nombreuses lignes. Comme on peut voir sur l'exemple suivant : https://zupimages.net/viewer.php?id=20/16/oaem.png En vert : ma couche de route regroupée (ligne_regroup) En rouge : le résultat de "route_cut_500_metres" Les points : Résultat de la requête qui me crée des points tous les 500 mètres. Merci. ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=332169#p332169 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

4 years

1
7
0 0

PostgreSQL - Extension ogr_fdw et SSL certificate problem

by geobd＠georezo.net

Message posté par : Alicemartin ---------------------------------------- Bonjour, J'ai PostgreSQL 9.5 installé en local sur un ordinateur sous Windows 10. J'y ai installé l'extension ogr_fdw afin de pouvoir importer dans ma base des couches issues de flux WFS en suivant ce tutoriel --> https://mybl0gn0tes.wordpress.com/2017/06/13/utiliser-des-flux-wfs-dans-pos…). Mon soucis est que, quand je crée un serveur distant avec des URL en https (ex : https://wxs-gpu.mongeoportail.ign.fr/externe/39wtxmgtn23okfbbs1al2lz3/wfs), j'ai un message d'erreur indiquant "SSL certificate problem: unable to get local issuer". J'ai passé la journée à essayer de régler le problème mais sans succès. Quelqu'un saurait-il quoi faire ? Je vous remercie d'avance pour votre aide et vous souhaite une belle soirée. ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=333657#p333657 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

4 years, 7 months

1
4
0 0

Trigger spatial Oracle

by geobd＠georezo.net

Message posté par : Tomapinfo (numerobisco04(a)hotmail.fr) ---------------------------------------- Bonjour, J'essaye de reproduire un trigger qui fonctionne bien dans PostgreSQL/Postgis dans Oracle. Le trigger remonte des informations lors de la saisie ou de la modification. ----------------- Code : CREATE OR REPLACE FUNCTION x.champs_auto() RETURNS trigger AS $BODY$BEGIN NEW.num_insee = x.communes.num_insee FROM x.communes WHERE st_intersects(NEW.geom, x.communes.geom); RETURN NEW; END$BODY$ LANGUAGE plpgsql VOLATILE COST 100; ----------------- Mais dans Oracle, l'ecriture n'est pas identique et au mieux j'ai réussi à faire ça. ----------------- Code : CREATE OR REPLACE TRIGGER "DBO"."Cps_auto_I" BEFORE INSERT -- ON Y -- FOR EACH ROW -- declare INSEE VARCHAR2(20) DEFAULT ''; Begin -- Select b.NUM_INSEE into INSEE from ua b WHERE SDO_RELATE(:NEW.GEOMETRIE,b.GEOMETRIE,'mask=INSIDE') = 'TRUE' and b.SOUS_TYPE = 'COMMUNES'; :NEW.NUM_INSEE := INSEE ; ; End; ----------------- Mais j'obtiens cette erreur alors que j'ai bien un index sur les 2 tables. ----------------- Citation : Erreur SQL : ORA-13226: interface non prise en charge sans index spatial ----------------- Je penses donc que l'objet n'est pas encore enregistré et donc pas dans l'Index. Car le select seul fonctionne bien. Au vue de ce message je me demande s'il y est vraiment possible de faire se genre de requête dans un trigger ? Si cela est possible avez-vous une idée pour le corriger ? Merci par avance. ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=335958#p335958 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

4 years, 9 months

1
1
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=335210#p335210 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

4 years, 9 months

1
4
0 0

PgAdminn 4 - Modification de droit

by geobd＠georezo.net

Message posté par : Bruno IRATCHET (bruno(a)georezo.net) ---------------------------------------- Sujet déplacé dans le bon forum. ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=336644#p336644 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

4 years, 11 months

1
3
0 0

Grosse bdd géométrique et QGis

by geobd＠georezo.net

Message posté par : preliator ---------------------------------------- Bonjour, Je dispose d'une bdd Postgis (nommée "STAGE") dans laquelle j'ai importé la base de données Sirene v3. (https://public.opendatasoft.com/explore/dataset/sirene_v3/information/?sort…) de 33 millions d'entrées. Etant au format CSV, je me suis servie des colonnes de Lat et Long pour créer une colonne "geom" de type geometry pour chaque ligne. Je souhaite manipuler cette bdd sur QGis, uniquement sur une commune. Après avoir crée une table "sirene_marseille" sur Postgre, filtrée sur la ville, je rencontre une très grande lenteur de QGis pour se connecter à ma BDD. Voici ce que j'obtiens durant 5 bonnes minutes : https://zupimages.net/viewer.php?id=20/42/ofd5.png Je rencontre ce problème, directement après avoir crée la géométrie de ma base Sirene. Existe t-il une alternative pour stocker une grosse bdd geometrique dans Postgis, tout en étant traitable sur QGis ? Merci. ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=336245#p336245 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

4 years, 11 months

1
3
0 0

PgAdminn 4 - Modification de droit

by geobd＠georezo.net

Message posté par : BadWolf ---------------------------------------- Je vous remercie je ne savais pas que c'était sur ce forum. ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=336647#p336647 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

4 years, 11 months

1
2
0 0

2025

2024

2023

2022

2021

2020

2019

Geobd October 2020