Message posté par : image95
----------------------------------------
Bonjour,
Je me permets de poster un message. Il y a quelques semaines, j'avais posté un premier message :
https://georezo.net/forum/viewtopic.php?id=136209
mais après réflexion, je me permets d'en créer un nouveau pour apporter une précision/clarification.
Pour rappel, mon but est de trouver une méthode/outil pour comparer la similarité/dissimilarité entre 2 couches vecteurs polygonales constituées de plusieurs entités. Dans mon cas d'usage, cette comparaison concerne une couche polygonale issue d'un travail de photo-interprétation (modélisation des grandes unités physionomiques/habitats). je joins une capture écran, cette couche apparait en couleur bleu).
L'autre couche vecteur polygonale est une couche vecteur polygonale issue d'un traitement de segmentation (via un outil OTB disponible via Qgis).(elle est en couleur jaune sur la capture écran jointe).
Dans mon précédent post, j'évoquais comme indicateur possible de similarité : "pour chaque entité/polygone photo-interprété, calculer le pourcentage de surface recouverte par l'entité de la couche issue de la segmentation ayant le grand recouvrement."
Pour être plus parlant, dans l'exemple de la capture écran, ce pourrait être :
"pour l'entité de la couche photo-interprétation numéro 11 (bleu), j'aimerais connaitre la part de surface recouverte par l'entité de la couche de segmentation ayant le plus grand recouvrement avec cette entité 11 (dans ce cas précis, c'est probablement l'entité 42 jaune). "
Cela pourrait être : surface de l'entité 42 en superposition avec l'entité 11 / surface totale de l'entité 11.
Mais après réflexion, je pense que c'est insuffisant. Qu'il faudrait aller un peu plus loin.
Pour ces aspects méthodes pour quantifier la précision d’une segmentation/fidélisation à la PI. Le niveau/degré de ressemblance, je pense qu'il faudrait intégrer un aspect/ une logique de "bi-similarité" : d’avoir une approche par couples de polygones les plus couvrants :
Étape 1 :
- Pour chaque entité polygone de la PI, récupérer l'id et la geom du polygone de la segmentation qui a le plus fort recouvrement. En calculer le taux de recouvrement.
(ex : surface de l'entité 42 en superposition avec l'entité 11 / surface totale de l'entité 11).
Étape 2:
-Une fois ce premier indicateur calculé, faire le calcul pour le même couple/la même séquence de polygones calculée lors de l'étape numéro 1 : effectuer le même calcul mais cette fois sur la surface totale de l'entité 42 (polygone issu de la segmentation)
(ex : surface de l'entité 42 en superposition avec l'entité 11 / surface totale de l'entité 42).
Étape 1/2 variante :
Je me disais qu'une variante/synthèse/résumé en 1 étape des étapes 1 et 2 pouvait être :
- (ex : surface de l'entité 42 en superposition avec l'entité 11 / surface totale constituée par l'entité 11 ET par l'entité 42).?
- étape 3 :
l'idée serait à partir des étapes 1/2, d'avoir une sorte de score global au niveau de la couche. Qui pourrait être une sorte moyenne/médiane pour tous les couples de polygones ? Ou une autre métrique ?
J'imagine que cet outil n'existe/,'est pas mis à disposition sous Qgis. Bien que j'ai découvert ce plugin :
https://plugins.qgis.org/plugins/similarity_plugin/#plugin-abouthttps://github.com/panickspa/SimilarityPlugin
Plugin intéressant. Mais qui à mon humble avis ne couvre qu'une partie du besoin.
Si des personnes pensent pouvoir m'aiguiller ?
Afin d'être complet, ayant fait un peu de recherche sur internet. j'ai également trouvé un indicateur potentiellement intéressant :
https://postgis.net/docs/ST_HausdorffDistance.htmlhttps://en.wikipedia.org/wiki/Hausdorff_distance
Si une personne pense pouvoir m'aiguiller. Un grand merci.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=375096#p375096
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : Alban Kraus
----------------------------------------
Bonjour,
-----------------
BHaas écrit :
Le soucis que je rencontre c'est que lors de la création des points de rencontre (via PGAdmin) mes points se retrouvent au milieu de l'océan Atlantique.
Pourtant d'après les propriétés de la couche le système de projection est bien le Lambert 93 (ESPG: 2154).
-----------------
Si vos données se retrouvent au large du golfe de Guinée, le diagnostic est clair : vous tentez d'afficher comme du Lambert-93 des données qui sont en réalité des latitudes-longitudes. Si vous associez, dans les propriétés de la couche, le SRS EPSG:4326, vos points se retrouvent-ils au bon endroit ?
-----------------
BHaas écrit :
Je me suis demandé si le problème ne venait pas de ma couche de base servant au calcul des points mais ce n'est pas le cas, les données apparaissent bien là où elles sont supposées être et avec le bon EPSG.
-----------------
Il me semble que le monde du transport aime beaucoup les données en latitude-longitude. Voudriez-vous vérifier très attentivement toutes les étapes de l'intégration en base PostGIS :
* Comment sont exprimées les coordonnées dans votre jeu de données initial ? (si c'est du GTFS, ce sont forcément des latitudes-longitudes)
* Comment les avez-vous chargées en base PostGIS ? Avez-vous bien spécifié le SRS des données sources et demandé explicitement une transformation vers Lambert-93 ?
* Dans votre base de données PostGIS, avez-vous verrouillé vos colonnes de géométrie pour n'accepter que des géométries se présentant comme étant exprimées en Lambert-93 ?
-----------------
BHaas écrit :
J'ai utilisé la requête sql: SELECT UpdateGeometrySRID('bus_lignes_pt', 'geom', 2154); pour vérifier que le SRID était bien le 2154
-----------------
Non, cette requête force les données à se présenter comme étant du Lambert-93, mais ne modifie pas les coordonnées. Un peu comme si vos données étaient capables de parler, avant : « Bonjour, je suis (0.152, 45.146) et je suis en 4326 » après « Bon, ben je suis (0.152, 45.146) et je suis en 2154 ».
Vous auriez plutôt souhaité utiliser :
-----------------
Code :
ALTER TABLE bus_lignes_pt ALTER COLUMN geom SET DATA TYPE GEOMETRY(Point, 2154);
-----------------
qui vous aurait renvoyé une erreur si les géométries n'étaient pas en 2154. Ou :
-----------------
Code :
SELECT DISTINCT ST_SRID(geom) FROM bus_lignes_pt;
-----------------
qui vous aurait affiché le (ou les différents) SRID de vos données.
Vous souhaitant un bon débogage,
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=375084#p375084
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : BHaas (b.haassig(a)gmail.com)
----------------------------------------
Bonjour à toutes et à tous,
Un peu de contexte:
Je travail à la modélisation de lignes de bus, pour éviter de me retrouver avec des lignes superposées je suis le tuto suivant: https://geotribu.fr/articles/2021/2021-04-07_carte_reseau_bus/#creer-des-po…
Le soucis que je rencontre c'est que lors de la création des points de rencontre (via PGAdmin) mes points se retrouvent au milieu de l'océan Atlantique.
Pourtant d'après les propriétés de la couche le système de projection est bien le Lambert 93 (ESPG: 2154).
Je me suis demandé si le problème ne venait pas de ma couche de base servant au calcul des points mais ce n'est pas le cas, les données apparaissent bien là où elles sont supposées être et avec le bon EPSG.
J'ai utilisé la requête sql: SELECT UpdateGeometrySRID('bus_lignes_pt', 'geom', 2154); pour vérifier que le SRID était bien le 2154 et ca n'a rien changé.
Le problème doit avoir un lien avec la première étape de calcul puisque les données de bases sont bonnes mais je ne trouve pas d'où provient l'erreur et comment la corriger.
Est-ce que quelqu'un à déjà rencontré ce genre de problème ?
PS: Je n'ai pas trouvé de section dédiée aux questions SQL désolé je poste au mauvais endroit
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=375082#p375082
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : AlineC
----------------------------------------
Ca y est, j'y suis arrivé !
Bizarrement le filtre que j'avais mis sur ma couche commune (insee=@atlas_pagename) ne marchait pas vendredi soir mais marchait ce matin !
Et après, les étiquettes étaient toujours visibles ! Je pense que c'est parce que c'était des étiquettes à partir d'une table jointe. J'ai donc rajouté mon filtre dans chaque symbologie. Et cette fois ci, le rendu est comme je le souhaitas.
Je testerais peut-être aussi la solution de conejo avec le code intersect.
Merci à toutes/tous
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=375081#p375081
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : Jean Cascalès
----------------------------------------
Bonjour,
Dans les ensembles de règles de l'onglet Symbologie, si on veut différencier la commune de l'atlas par rapport aux autres communes, il faut écrire
Code:
-----------------
Citation :
$id=@atlas_featureid
-----------------
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=375066#p375066
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : AlineC
----------------------------------------
Bonjour,
J'ai un cache commune que je voudrais activer comme cache inversée dans un atlas pour que chaque page de l'atlas ne représente que les éléments de cette commune (c'est distinct de la maille que j'ai dessiné pour caler l'emprise des données dans chaque commune).
J'y suis pas encore arrivée alors j'appelle à l'aide en ce vendredi pluvieux !
Merci
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=375053#p375053
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : AlineC
----------------------------------------
Coucou Santanna,
J'ai mis du temps à oser utiliser le bazooka, mais effectivement il n'y a que ça qui marche.
Merci
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=375052#p375052
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : Nourredine IDIR
----------------------------------------
Bonjour,
J'aurais besoin de vos retours d'expériences.
Je voudrais déployer QGIS 3.40 LTR avec la suppression (désinstallation) de l’ancienne version et mise à jour du profil, sur une trentaine d’ordinateurs de façon automatisée et silencieuse.
Nous sommes sous Windows 10 et 11, et les utilisateurs n'ont pas les droits "admin" sur leur poste.
Quelle méthode ou logiciel utilisez-vous ?
- GPO active directory, Windows server 2022.
- QGIS Deployment Toolbelt (QDT).
- WAPT ou autres logiciels payants.
- Logiciels open source.
- Autres méthodes.
Quel fichier utilisez-vous pour l’installation ?
- osgeo4w-setup.exe
- QGIS-OSGeo4W-3.40.12-1.msi
Merci par avance,
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=375045#p375045
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : ckiko13
----------------------------------------
Bonjour
J'ai une couche de points avec 2 champs:
- rue
- id_rue
Ma table ressemblerait à ceci:
RUE RUE_ID
rueA rueA_1
rueA rueA_2
rueB rueB_1
rueC rueC_1
rueC rueC_2
rueC rueC_3
rueC rueC_4
Lorsque j'ajoute un nouveau point par exemple dans la rueA, j'aimerais lorsque je remplis le champs RUE, que le champ RUE_ID prenne la suite de la 'numérotation', ici: rueA_3.
J'ai regardé du côté des fonctions ARRAY_LAST, mais je ne débouche sur rien.
Auriez-vous des pistes à me proposer?
Merci
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=375033#p375033
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : Sylvain M.
----------------------------------------
Bonsoir à tous,
J'ai une base de données au format GPKG, sur laquelle j'ai une chaine de requêtes SQL pour mettre à jour les données attributaires (des UPDATE principalement).
Dans le gestionnaire de bases de données de QGis, je suis obligé de sélectionner chaque requête une à une pour les exécuter. Si j'en sélectionne plusieurs à la suite, j'ai l'erreur :
-----------------
QGis écrit :
You can only execute one statement at a time.
-----------------
Je pense que pour enchainer les requêtes SQL en une seule fois, le mieux serait de passer par un script Python (?).
Je me débrouille un peu en Python (en dehors de QGis), mais le souci, c'est que je n'arrive pas à activer l'extension "mod_spatialite" pour que les fonctions spatiales du module sqlite3, et mes requêtes, bien qu'attributaires uniquement, concerne des couches spatiales.
J'ai par exemple ce message pour un simple UPDATE :
-----------------
Python écrit :
Erreur pour : UPDATE matable SET monchamp = mavaleur... -> no such function: ST_IsEmpty
-----------------
Comment feriez-vous pour automatiser une chaine de requêtes SQL sur un GPKG, vous ?
Je n'ai jamais utilisé la console Python de QGis (je ne sais pas faire) : peut-être serait-ce plus simple que mon installation Python Windows ?
Merci pour vos conseils !
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=375017#p375017
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net