Le portail francophone de la géomatique

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

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.

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 :

intersects(  $geometry , geomFromWKT( $current_canvas_extent ))

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

Mathieu

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

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

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/

Bonsoir à tous,

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

Le filtre suivant :

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,

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

Bug corrigé ;-)

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 :

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 :

intersetcs( $geometry , $current_canvas_extent )

A tester.

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 :

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)

Diantre... l'erreur idiote...

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

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

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.

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

Bonsoir,

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

A suivre

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 ?