September 2020 - Geomatique - Listes du Forum GeoRezo

Mise a jour des donnees SHP sous postrges

by geomatique＠georezo.net

Message posté par : ktl (r.njaraherizo(a)gmail.com) ---------------------------------------- Bonjour à toutes et tous, J'aimerais faire une mise à jour (ajout d'autre géomètrie) des shapefile qui sont déjà stocker sous postegresql via une formulaire html. Mais je n'arrive pas de trouver de requête pour effectuer cette tache. Quelqu'un peut m'aider à ce genre de situation? Merci! ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=335425#p335425 Pour y répondre : geomatique(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
1
0 0

Cites mayas: trouver le systeme de projection d'un SHP sans PRJ

by geomatique＠georezo.net

Message posté par : Nicolas89 (occis.nicolas(a)gmail.com) ---------------------------------------- Bonjour, Étant complètement novice en la matière, j'espère que vous serez bienveillant avec ce message. :) Si vous avez quelques instants, je sollicite votre aide sur le sujet suivant. Récemment, la technologie LIDAR a permis de découvrir des milliers de vestiges mayas inconnus au Nord du Guatemala : https://decryptageo.fr/la-complexite-des-anciennes-cites-mayas-desormais-de… Je souhaiterais les représenter sur une carte via GoogleMaps ou QGis. Le jeu de données est public, et peut être téléchargé à droite de https://datadryad.org/stash/dataset/doi:10.5061/dryad.k51j708 Mais les fichiers SHP n'ont pas de PRJ associés. Comment peut-on connaître le système de projection utilisé ? Dans QGis, j'ai essayé d'appliquer les projections - SR-ORG:6731: GTM - SR-ORG:6866: Guatemala GTM (42500) - EPSG:32061: NAD27 / Guatemala Norte - EPSG:32062: NAD27 / Guatemala Sur proposées par https://spatialreference.org/ref/?search=guatemala Mais la zone archéologique apparaît systématiquement dans l'océan trop au Sud-Ouest. Merci d'avance pour votre aide :) Cordialement, Nicolas ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=335328#p335328 Pour y répondre : geomatique(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
5
0 0

transformer polylignes en polygones

by geomatique＠georezo.net

Message posté par : régisj ---------------------------------------- Bonjour, Quelqu'un connait-il un moyen simple des transformer des milliers de polylignes (triangles) en polygones depuis Postgis ? J'ai essayé : DROP TABLE IF EXISTS icpe.test; CREATE TABLE icpe.test AS (SELECT ST_MakePolygon (secteur_polyligne_r44.geom)as geom FROM icpe.secteur_polyligne_r44); et j'ai le message : ERREUR: Shell is not a line ********** Erreur ********** Est-ce à dire que mes polylignes (secteur_polyligne_r44) ne sont pas fermées? Merci pour votre aide. Cordialement ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=335336#p335336 Pour y répondre : geomatique(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
1
0 0

QGIS Extension: differents points de vues ?

by geomatique＠georezo.net

Message posté par : tony13160 (tony.guyot(a)resodetection.com) ---------------------------------------- Bonjour, Je ne sais pas si ici je peux poser ma question ? mais je la pose quand meme : J'aimerai connaitre le nom de l'extension Qgis afin de pouvoir observer un plan selon differents axes : X et Y : vue normale du dessus et par ex X et Z pour une vue de cote. Merci de votre aide, Cordialement. ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=335292#p335292 Pour y répondre : geomatique(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
0
0 0

Flux WMS Géoportail de l'urbanisme

by geomatique＠georezo.net

Message posté par : jlavaure ---------------------------------------- @IGN_GPU Je vous remercie pour votre réponse. Mais si on prend la couche "Zonages", peut importe le niveau de zoom il n'y a rien dans mon QGIS et peut-importe la version. On a essayé sur plusieurs postes dans notre entreprise et c'est le même problème. Avez-vous une solution à nous apporter ? ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=330171#p330171 Pour y répondre : geomatique(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

QGIS 3.10.1 Acorona: Remplissage automatisé des champs

by geomatique＠georezo.net

Message posté par : KIENDREBEOGO Yembi ---------------------------------------- Bonsoir. Besoin d'aide. Je travaille avec QGIS 3.10.1 Acorona. J'aimerais nettoyer de façon automatisée du "texte" ou du "nombre" dans les enregistrements d'un champ. Quelles sont les formules adéquates? Merci de me venir en aide SVP! ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=335112#p335112 Pour y répondre : geomatique(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

[Academique] [Enquete] Representations SIG 3D pour la participation

by geomatique＠georezo.net

Message posté par : tibo.c (thibaud.chassin(a)gmail.com) ---------------------------------------- Bonsoir à tous, Je suis actuellement doctorant en géovisualisation au Laboratoire LaSIG de l'EPFL. Ma thèse porte sur l'utilisation de SIGs 3D dans les démarches participatives de projets d'aménagement du territoire. Je lance aujourd'hui une phase d'expérimentation visant à évaluer l'efficacité de différentes représentations 3D. Je vous propose donc une enquête d'environ 20 min sous la forme d'un petit jeu. Les informations que vous me communiquerez resteront anonymes et confidentielles et ne seront analysées que globalement avec l’ensemble des réponses des autres personnes qui auront répondu. C'est maintenant à vous de jouer ! Merci d'avance pour votre aide et n'hésitez pas à partager :) Lien de l'enquête : http://dev-mockup3d.heig-vd.ch/ (Un pc vous sera nécessaire et la langue se change en haut à droite) ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=334622#p334622 Pour y répondre : geomatique(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

Optimisation du ST_DIFFERENCE

by geomatique＠georezo.net

Message posté par : hadri45 (hadrien.siney(a)me.com) ---------------------------------------- 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=335100#p335100 Pour y répondre : geomatique(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

2025

2024

2023

2022

2021

2020

2019

Geomatique September 2020