Le portail francophone de la géomatique

MathieuB · Fri 27 September 2013 17:03

Bonjour à tous,

une petite colle avant le week-end :

J'utilise depuis Qgis 1.8 les possibilités de personnalisation du formulaire de renseignement des attributs.
Une nouveauté (une autre!) a fait son apparition avec la sortie de QGis 2, c'est la possibilité d'utiliser des couples clé/valeur stockés dans une table de la base de données.

J'aimerai aller un peu plus loin et appliquer un "filtre d'expression", consistant à ne proposer dans la liste déroulante que les couples clé/valeur des objets présents dans l'emprise courante de ma carte.

Concrètement, j'ai besoin de délimiter du parcellaire et je dispose pour cela du scan du cadastre et des localisants de parcelles fournis par la BD Parcellaire.

Lors du dessin d'une nouvelle parcelle, j'aimerai que les valeurs proposées pour renseigner le numéro de parcelle soient issues des valeurs de numéro de parcelle des localisants contenus dans l'emprise actuelle de ma carte, et non de l'ensemble des localisants de ma bd parcellaire.

L'éditeur de filtre ne propose rien qui soit relatif à mon objet "carte".

Pour être plus clair, je joins une capture d'écran de l'interface.

Si quelqu'un a une idée... Merci d'avance,

Fichier(s) joint(s) :
Pour accéder aux fichiers vous devez vous inscrire.

Robin · Mon 30 September 2013 09:14

Jamais testé encore, mais rien du côté de géométrie ne te conviendrait ?

MathieuB · Tue 01 October 2013 09:17

Bonjour Robin,

je pense en effet que c'est par là que ça va se passer. J'ai bien testé ceci :

Code:

intersects( $geometry, geomFromWKT( 'LINESTRING(765145.88 6274561.22, 776031.67 6284189.52)' ) )

Comme présenté par René-Luc D'Hont ici : http://www.3liz.com/blog/rldhont/index. … ns-de-QGIS

Ça fonctionne mais j'aimerai que la géométrie passée en paramètre

Code:

geomFromWKT( 'LINESTRING(765145.88 6274561.22, 776031.67 6284189.52)' )

soit une variable.

Je continue mes recherches...

MathieuB · Thu 24 October 2013 09:53

Bonjour à tous,

pour le suivi, j'ai abouti! En cherchant un peu partout une solution à mon problème, j'ai bien compris que python serait la solution.

Les 3 ressources ci-dessous m'ont bien aidé :

- http://nathanw.net/2012/11/10/user-defi … -for-qgis/
- http://www.3liz.com/blog/rldhont/index. … ns-de-QGIS
- http://gis.stackexchange.com/questions/ … n-qgis-1-8

J'ai donc créé un script python (userfunctions.py), qui contient une fonction appelée current_canvas_extent, qui retourne une géométrie de type Polygone représentant l'emprise courante de la fenêtre.

Code:

from qgis.utils import qgsfunction
from qgis.utils import iface
from qgis.core import QGis

@qgsfunction(0, "Python")
def current_canvas_extent(values, feature, parent):
    """
       retourne l etendue courante de la carte
   """
    extend = iface.mapCanvas().extent().asWktPolygon() 
    return extend

J'appelle se script au démarrage de QGis (ajouter à la commande de lancement de QGis l'option "--code chemin\vers\userfunctions.py")

J'ai donc maintenant accès dans mon éditeur d'expression à cette "variable".
Et je peux demander à QGis, grâce à l'éditeur d’expression de ne me proposer dans la liste que les parcelles présentes dans l’emprise de ma carte :

Code:

intersects(  $geometry , geomFromWKT( $current_canvas_extent ))

Je propose que "QGis Déchire" retrouve sa place dans l'aperçu d'étiquetage!

Mathieu

MathieuB · Thu 24 October 2013 10:57

Petit complément d'information de la part de Nathan Woodrow (http://plugins.qgis.org/plugins/qgsexpressionsplus/),

"If you add that into startup.py in .qgis2/python/ it will always be there when you load QGIS."

haubourg · Tue 29 October 2013 11:13

Bonjour Matthieu,
ça donne plein d'idées!

On se creuse la tête depuis un moment pour financer une fonctionnalité permettant de masquer l'affichage des objet ou des étiquettes en dehors d'une zone d'intérêt, et cette solution semble fonctionner (avec un département). Il ne reste qu'à modifier le plugin de génération de couche de masquage pour qu'il génére au passage une variable dynamique utilisateur contenant la zone d'intérêt.
Il faut aussi trouver un moyen de faire prendre ça en compte par le générateur d'Atlas.. si quelqu'un a des idées..
Régis

Cornet Jérémie · Tue 29 October 2013 11:39

haubourg a écrit:

Bonjour Matthieu,
ça donne plein d'idées!

On se creuse la tête depuis un moment pour financer une fonctionnalité permettant de masquer l'affichage des objet ou des étiquettes en dehors d'une zone d'intérêt, et cette solution semble fonctionner (avec un département). Il ne reste qu'à modifier le plugin de génération de couche de masquage pour qu'il génére au passage une variable dynamique utilisateur contenant la zone d'intérêt.
Il faut aussi trouver un moyen de faire prendre ça en compte par le générateur d'Atlas.. si quelqu'un a des idées..
Régis

Le plus simple serait de financer une amélioration du générateur d'atlas je pense.
Un fonctionnement power-user serait de pouvoir définir pour chaque couche une requête (tel que faisable à la main dans les propriétés de la couche) en pouvant utiliser les variables de l'atlas (y compris une variable contenant la géométrie en wkt qui n'existe pas pour le moment je pense).

L'autre solution est de faire ce fonctionnement à travers un plugin ou un script python.
Voire les classes suivantes
http://www.qgis.org/api/classQgsAtlasComposition.html
http://www.qgis.org/api/classQgsVectorLayer.html (voir la méthode setSubstring)

Dernière modification par Cornet Jérémie (Tue 29 October 2013 11:39)

MathieuB · Tue 29 October 2013 13:36

Bonjour à vous,

ici aussi on s'intéresse à la possibilité de masquer (avec +/- de transparence) ce qui serait hors de la zone d'impression courante de l'atlas (+/- un tampon à définir) et je suis tombé hier sur cette discussion : http://osgeo-org.1560.x6.nabble.com/Let … 21391.html !

Je suis assez favorable à l'idée de trouver de quoi financer l'ajout de cette fonctionnalité au générateur d'atlas.

Régis, tu dis que "cette solution semble fonctionner (avec un département)", qu'entends tu par là ?

Mathieu

haubourg · Tue 29 October 2013 13:58

MathieuB a écrit:

Bonjour à vous,

ici aussi on s'intéresse à la possibilité de masquer (avec +/- de transparence) ce qui serait hors de la zone d'impression courante de l'atlas (+/- un tampon à définir) et je suis tombé hier sur cette discussion : http://osgeo-org.1560.x6.nabble.com/Let … 21391.html !

Mathieu

Oui, ça fait un moment que j'essaie d'obtenir un accord dans les listes anglophone sur la meilleure façon de répondre à ce besoin.
Le problème de la couche semi transparente (déjà gérée par le plugin Mask) est qu'elle n'interfère pas avec l'étiquetage, et tous mes utilisateurs me réclame de pouvoir masquer ET cacher les étiquettes en dehors de la zone d'intérêt ayant servi à générer le masque.
- Soit on prend des solutions type mapserver - mapinfo, avec un mask image ajouté à la fin de l'étape de rendu graphique, et là on a des étiquettes coupées,
- soit on crée une vrai notion de masque dans l'API avec gestion de l'ordre de dessin des étiquettes par rapport à ce masque pour pouvoir dessiner certaines étiquettes au dessus
- soit on travaille sur la sélection des objets à afficher, et ce que permet de faire assez simplement la solution avec les variables dynamiques décrites dans ce post.

Régis, tu dis que "cette solution semble fonctionner (avec un département)", qu'entends tu par là ?

mon gros souci est la performance de cette solution. L'intersection spatiale est très vite lente et consommatrice si les objets géographiques grandissent, voir impossible si on a des géométries invalides. Un département GEOFLA croisé avec des lacs BD Carthage est déjà un assez gros objet, et je ne vois pas de ralentissement notable. L'idée du plugin "proof of concept' serait de trouver les pistes d'optimisation (nettoyage et simplification des objets à la volée, utilisation correcte des indexes spatiaux, gardes fous pour les très gros objets, utilisation d'autres opérateurs que le intersects.. ) A voir également, l'intersection prend des ressources, mais en économise sur l'affichage.. On y sera peut-être gagnants au final. A tester donc.
Régis

Cornet Jérémie · Tue 03 December 2013 23:33

A propos de l'atlas et de la possibilité de filtrer les autres couches, ça sera présent dans Qgis 2.2 :
http://nathanw.net/2013/12/02/waiting-f … s-feature/

MathieuB · Wed 04 December 2013 10:30

J'en ai rêvé, QGis le fait ;-)

MathieuB · Tue 15 July 2014 20:10

Bonsoir à tous,

Patatra, impossible de valider cette manip sur QGis 2.4...

Le filtre suivant :

Code:

intersects( $geometry, geomFromWKT( 'LINESTRING(765145.88 6274561.22, 776031.67 6284189.52)' ) )

retourne bien 1 mais la liste déroulante reste vide.

Avez vous rencontré ce souci

Merci d'avance,

MathieuB · Tue 15 July 2014 20:23

Pour précise un peu et résumer le sujet décrit ci-dessus...

J'utilise dans mes formulaire des valeurs relationnelles. Je souhaite donc alimenter la liste déroulante des communes par la couche (en lambert 93) des communes de la région. Pour limiter la liste aux seules communes présentes sur la carte, j'utilise un filtre contenant une variable python :

Code:

intersetcs(  $geometry , geomFromWKT(  $current_canvas_extent ))

La variable retourne bien l'étendue de la carte au format WKT, mais la liste reste désespérément vide. Cela fonctionne bien avec QGis 2.2.

Et j'arrive à utiliser correctement ce test par exemple pour mettre à jour une colonne...

Fichier(s) joint(s) :
Pour accéder aux fichiers vous devez vous inscrire.

MathieuB · Fri 18 July 2014 08:26

Bug remonté : http://hub.qgis.org/issues/10900

Sylvain33 · Fri 18 July 2014 11:29

MathieuB a écrit:

Bug remonté : http://hub.qgis.org/issues/10900

Merci pour le suivi j'ai le même style de requêtes grrr...

MathieuB · Thu 24 July 2014 12:20

Bug corrigé ;-)

snaileater · Wed 30 July 2014 11:02

Bonjour,
une question pour ma culture personnelle :

Quel est l'intérêt de passer par le WKT pour l'import/export de l'extent ?

Merci d'avance.

MathieuB · Thu 31 July 2014 11:53

Bonjour,

c'est trés simple : j'ai pris cette voie en explorant les objets python, ça a fonctionné, je suis passé à autre chose.

Ca meriterait peut-être d'être repris :

Code:

def current_canvas_extent(values, feature, parent):
        """ retourne l etendue courante de la carte """
        #extend = iface.mapCanvas().extent().asWktPolygon()
        extend = QgsGeometry.fromRect( iface.mapCanvas().extent())
        return extend

devrait permettre d'employer ce filtre :

Code:

intersetcs( $geometry , $current_canvas_extent )

A tester.

MathieuB · Thu 31 July 2014 13:08

Voilà qui fonctionne :

Code:

from qgis.utils import qgsfunction
from qgis.utils import iface
from qgis.core import (QGis, QgsGeometry)

@qgsfunction(0, "Python")

def current_canvas_extent(values, feature, parent):
        """ retourne l etendue courante de la carte """
        extend = QgsGeometry.fromRect( iface.mapCanvas().extent())
        return extend

http://si.cenlr.org/24-10-2013/limiter-les-propositions-d-une-liste-relationnelle-selonl-l-empise-de-la-carte

Attention, la version 2.4 comporte un bug qui empêche l'utilisation de ce filtre. Utiliser une version précédente ou une version corrigée : http://hub.qgis.org/issues/10900.

Zakh · Wed 24 August 2016 14:03

Bonjour,

Est-ce qu'il est possible de donner une géométrie différente comme élément de sélection?

Je souhaite insérer un tableau de données dans mon composeur d'impression, qui n'afficherait que les ligne pour les entités issues d'une intersections. Il s'agit en réalité des chambres du réseau d'assainissement, pour lesquels je souhaites ressortir certaines caractéristiques (profondeurs, etc) en lien à un projet, lui-même délimité par un polygone dont les sommets sont connus (fournie par le gestionnaire du projet).

Comme cette opération arrive régulièrement, je souhaitais faire une expression pour laquelle je dois juste modifier les sommets de mon polygone.

j'ai bien lu la page

http://www.3liz.com/blog/rldhont/index. … ns-de-QGIS

mais reste dans le flou. La page propose un geomFromWKT alors que qgis (2.14) dispose dans l'éditeur d'expression de geom_from_wkt. Les deux sembles reconnues par le programme. Est-ce juste?

J'ai un polygone, dont les sommets sont connus et j'essaie donc ceci :

Code:

intersects($geometry,  geom_from_wkt('POLYGONE(
536088.99 154071.52,
536096.57 154078.87,
536104.61 154098.35,
536113.97 154132.51,
536128.57 154138.56,
536149.74 154126.49,
536141.74 154225.25,
536150.51 154235.99,
536240.99 154246.94,
536251.84 154239.62,
536268.98 154176.29,
536262.51 154164.20,
536195.32 154141.66,
536204.92 154082.42,
536200.73 154072.59,
536154.54 154040.69)
'))

Mais si la requête ne retourne pas d'erreur à proprement parler, elle reste cependant vide.

Une idée?

Dernière modification par Zakh (Wed 24 August 2016 14:04)

VianneyD · Wed 24 August 2016 14:42

Bonjour,

Et en enlevant le 'E' de POLYGONE ?

Zakh · Wed 24 August 2016 14:48

Diantre... l'erreur idiote...

j'ai corrigé ce détail, mais aucun changement.

Dernière modification par Zakh (Wed 24 August 2016 14:49)

VianneyD · Wed 24 August 2016 15:12

Mmmmm... Une erreur de projection ?

Tes coordonnées ressemblent à du Lambert Zone.
C'est bien la même projection pour tes données d'entrée ?

Zakh · Wed 24 August 2016 15:58

Le système de projection est EPSG : 21781, il s'agit du système suisse. Elles sont bien dans la bonne projection.

Le polygone est un petit peu complexe (pas un simple rectangle) et est issu d'une commande de donnée réceptionnée que j'utilise en teste pour cet objectif de la table dans mon composeur.

Zakh · Wed 24 August 2016 16:06

Le diable se cache dans les détails comme on dit...

En recherchant dans ma doc Postgis, j'ai constaté que l'utilisation de POLYGON nécessite en fait deux paires de parenthèse... Là où je n'en avais mis qu'une seule.

Mea Culpa... Merci bien pour les suggestions!

amélimélo · Thu 24 November 2016 09:51

Bonjour,

Je déterre un peu ce sujet.
J'ai eu le plaisir de suivre ce tutoriel très bien fait http://si.cenlr.org/24-10-2013. Cela fonctionne très bien pour la version 2.14. Mais pour les versions supérieures (2.16 et 2.18) cela ne fonctionne plus.

Le lancement de la commande iface.mapCanvas().extent().asWktPolygon() via la console python de Qgis fonctionne bien.
Mais la requête intersects( $geometry , geomFromWKT( $current_canvas_extent )) ne retourne rien. Je pense que la fonction ne marche plus en raison de l'évolution de la version de python sur Qgis à partir de la 2.16. Simplement, je ne sais pas trop sur quoi agir. Si quelqu'un a des pistes de recherche, je le remercie d'avance.

Bonne journée,

Amélie

MathieuB · Fri 25 November 2016 13:09

Bonjour,

merci pour ce retour. J'ai regardé et c'est bizarre en effet, car les différents variables ou test ne retournent pas d’erreur, la liste reste cependant vide...

A creuser.

Mathieu

MathieuB · Wed 18 January 2017 17:40

Bonsoir,

j'ai rouvert le ticket car le "bug" est le même que pour la version 2.4...

A suivre

MathieuB · Thu 22 June 2017 13:19

Bonjour à tous,

bug résolu, si tant est que c'était un bug.
La syntaxe qui fonctionne pour le filtre est la suivante (QGIS 2.18.9) :

Code:

intersects( geometry( $currentfeature ), current_canvas_extent() )

ddtphilou · Thu 22 June 2017 14:54

bonjour à tous
je viens de tester sur une couche polygone avec comme référence la couche des communes N_COMMUNE_BDC_088.
dans la couche test j'ai créé une colonne NOM et j'ai activer "valeur relationnelle" pour cette colonne.
j'ai pris pour la sélection :
Couche = N_COMMUNE_BDC_088
Colonne clé = NOM_COMM
Colonne de valeurs = NOM_COMM

j'ai collé la requête intersects( geometry( $currentfeature ), current_canvas_extent() )

lorsque je digitalise, la liste est vide au moment de la saisie des attributs

je suis en QGIS 2.16.3

Est ce normal ?

Communauté	Ressources	L'emploi	Le marché	L'association	Fils RSS
Forum Blogs Wiki Planet Sites amis	Agenda Biblio Fiches AMO-CNIG Paris PCGIS Téléchargements	Offres d'emploi Banque de CV Formations Démarche Métiers	Geo-Entreprises Geo-Communiqué Appels d'offres	L'association Qui sommes-nous ? Nous aider Bienfaiteurs Partenaires Historique Statuts	Tout le site Job Géomatique Services GeoBlogs ...	Faire un don Mentions légales

Le portail francophone de la géomatique

Annonce

#1 Fri 27 September 2013 17:03

QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#2 Mon 30 September 2013 09:14

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#3 Tue 01 October 2013 09:17

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

Code:

Code:

#4 Thu 24 October 2013 09:53

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

Code:

Code:

#5 Thu 24 October 2013 10:57

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#6 Tue 29 October 2013 11:13

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#7 Tue 29 October 2013 11:39

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#8 Tue 29 October 2013 13:36

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#9 Tue 29 October 2013 13:58

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#10 Tue 03 December 2013 23:33

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#11 Wed 04 December 2013 10:30

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#12 Tue 15 July 2014 20:10

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

Code:

#13 Tue 15 July 2014 20:23

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

Code:

#14 Fri 18 July 2014 08:26

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#15 Fri 18 July 2014 11:29

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#16 Thu 24 July 2014 12:20

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#17 Wed 30 July 2014 11:02

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#18 Thu 31 July 2014 11:53

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

Code:

Code:

#19 Thu 31 July 2014 13:08

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

Code:

#20 Wed 24 August 2016 14:03

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

Code:

#21 Wed 24 August 2016 14:42

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#22 Wed 24 August 2016 14:48

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#23 Wed 24 August 2016 15:12

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#24 Wed 24 August 2016 15:58

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#25 Wed 24 August 2016 16:06

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#26 Thu 24 November 2016 09:51

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#27 Fri 25 November 2016 13:09

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#28 Wed 18 January 2017 17:40

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

#29 Thu 22 June 2017 13:19

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

Code:

#30 Thu 22 June 2017 14:54

Re: QGIS 2.0.1: Limiter les choix d'une liste selon l'emprise

Pied de page des forums