- Qgis_fr - Listes du Forum GeoRezo

QGIS: Extraction d'images PNG, depuis Geotiff pour creation database

by qgis_fr＠georezo.net

Message posté par : GuillaumeL7 (guillaume.launey(a)gmail.com) ---------------------------------------- Bonjour, J'ai un grand nombre de cibles (environ 500 à 1500 par projet). Ces points sont référencés dans un .shp. Je dois enregistrer pour chaque contact une imagette, extraite d'un Géotiff. Pour le moment je le fait de maniére manuelle, ce qui est extrémement fastidieux. Je cherche donc le moyen d'automatiser le process sur Qgis. Les contraintes sont: - Imagettes centrées sur le point du .SHP - Imagettes nommées selon le champ "ID_Unique" de la table attributaire du .shp - Imagettes de format 500x400 px - Imagettes au format .png Au cours de mes recherches, je n'ai pas trouvé d'outil directement disponible pour faire ces extraits et me suis donc tourné vers python. Mon code actuel: ----------------- Citation : from qgis.core import QgsVectorLayer, QgsRasterLayer, QgsProject import geopandas as gpd import rasterio from rasterio.transform import from_origin from PIL import Image import numpy as np # Fonction pour extraire et enregistrer les images def extract_and_save_images(raster_path, shapefile_path, id_field_name, output_directory, width, height): raster_dataset = rasterio.open(raster_path) vector_layer = QgsVectorLayer(shapefile_path, "vector_layer", "ogr") for feature in vector_layer.getFeatures(): geom = feature.geometry() bbox = geom.boundingBox() min_x, min_y, max_x, max_y = bbox.xMinimum(), bbox.yMinimum(), bbox.xMaximum(), bbox.yMaximum() window = raster_dataset.window(min_x, min_y, max_x, max_y) # Lecture des données raster dans la fenêtre spécifiée data = raster_dataset.read(window=window) # Conversion des données en tableau NumPy data_np = np.array(data) # Création d'une nouvelle image avec Pillow image = Image.fromarray(data_np) # Redimensionnement de l'image à la taille spécifiée resized_image = image.resize((width, height)) # Enregistrement de l'image redimensionnée output_file = output_directory + f"\{feature[id_field_name]}.png" resized_image.save(output_file) raster_path = "D:\TestImages\Merge_SSS_08+11.tif" shapefile_path = "D:\TestImages\a.shp" id_field_name = "ID_Unique" output_directory = "D:\TestImages\output\" width, height = 500, 400 # Appeler la fonction pour extraire et enregistrer les images extract_and_save_images(raster_path, shapefile_path, id_field_name, output_directory, width, height) ----------------- J'enregistre ensuite ce code au format .py, puis l'appel dans la console Python avec: ----------------- Citation : exec(open('D:\Image_RasterioI.py').read()) ----------------- Cette version est la 9ème de mon code et fait appel à de plus en plus de librarie, je suis probablement entrain de faire fausse route. Pourriez vous me conseiller sur la direction à prendre? ou tout simplement m'indiquer si un outil existe et que je l'ai raté? En vous remerciant et vous souhaitant une excellente Journée, Guillaume ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=364570#p364570 Pour y répondre : qgis_fr(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

2 years, 5 months

1
0
0 0

QGIS: Enrichir un fichier geolocalisée XY par l IRIS avec un SHP

by qgis_fr＠georezo.net

Message posté par : Canarelli (pcanarelli(a)hotmail.fr) ---------------------------------------- Bonjour Je dispose d'un fichier géolocalisé XY. Je souhaiterais l'enrichir avec le code IRIS en utilisant un fichier Shapefile à l'IRIS. Etant donné que le SHP contient les infos du polygone et l'IRIS associé ca devrait être possible? J'ai cherché sur QGIS (exporter , ...) mais je ne vois pas comment faire ça? Est ce que c'est possible avec QGIS (ou un autre soft) ? Quelqu'un a une idée pour faire ça ? En vous remerciant d'avance, Cdlt ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=364462#p364462 Pour y répondre : qgis_fr(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

2 years, 5 months

1
2
0 0

QGIS/Ouverture: Lancement script python

by qgis_fr＠georezo.net

Message posté par : rep20 ---------------------------------------- Bonjour, Merci beaucoup pour votre réponse. Je vais regarder alors pour la création du plugin. ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=364561#p364561 Pour y répondre : qgis_fr(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

2 years, 5 months

1
0
0 0

QGIS/Ouverture: Lancement script python

by qgis_fr＠georezo.net

Message posté par : Yoann Lecomte ---------------------------------------- Bonjour, Via startup.py, votre code est exécuté avant que votre projet ne soit chargé. Lors de son exécution, les couches ne sont donc pas encore disponibles pour d'éventuels traitement python. Cela explique pourquoi vous n'avez aucun message d'erreur et que rien ne se passe. A mon avis, votre meilleur chance est de créer un plugin qui se connecte a https://qgis.org/pyqgis/3.2/gui/other/QgisInterface.html#qgis.gui.QgisInter… (cela peut surement être fais à l'initialisation du plugin) et qui exécute votre code lorsque le projet est chargé. Là au moins vous serez sur que vos couches sont disponibles et vous pourrez interagir avec. ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=364560#p364560 Pour y répondre : qgis_fr(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

2 years, 5 months

1
0
0 0

Lancement script python - Ouverture QGIS

by qgis_fr＠georezo.net

Message posté par : rep20 ---------------------------------------- Bonjour à tous, Je viens d'écrire un script python qui permet de supprimer les couches temporaires vides sur un projet QGIS. Certains utilisateurs oublient de les supprimer avant de fermer les différents projets QGIS. Cela permet de nettoyer automatiquement les projets. Le script fonctionne. J'aimerais faire en sorte qu'il s'éxecute automatiquement lors de l'ouverture de chaque projet. J'ai réussi à le faire fonctionner en passant par les propriétés/macros du projet. Comme il y a plusieurs projets avec le même profil je voudrais que ce script s'éxecute via les fichiers du profil. Cela permettrait de modifier une seule fois le profil plutôt que d'ouvrir tous les projets. J'ai regardé via la création d'un fichier nommé startup.py dans les dossiers sources du profil. J'ai bien pointé les traitements python sur le dossier où ce trouve ce sript et passé l'execution des scripts en automatique. Toutefois cela ne fonctionne pas. Je me demande où se trouve l'erreur car je n'ai aucun message d'erreur et quand je lance le script à la main cela fonctionne. Auriez-vous une idée ou un conseil à suivre pour executer un script python lors de l'ouverture d'un ou plusieurs QGIS ? Merci et bonne journée ! Rémi ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=364552#p364552 Pour y répondre : qgis_fr(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

2 years, 5 months

1
0
0 0

QGIS: Creation d'une BDD avec cours d'eau/affluent/effluent

by qgis_fr＠georezo.net

Message posté par : chris15000 (ref(a)gocolo.fr) ---------------------------------------- Bonjour, Je suis novice sur QGIS et j'ai quelques questions. J'ai pour projet de réaliser une base de donnée qui me permettrait de lister tous les cours d'eaux existant et pour chacun d'eux indiquer : 1/ Les villes qu'ils traversent 2/ les departements qu'ils traversent 3/ si ils sont affluent d'un autre cours d'eau lequel 4/ a l'inverse determiner les affluents de ce cours d'eau (mais ça on doit pouvoir le faire en remontant le point 3 à l'envers) POur les points 1 et 2 ça a été plutot facile en utilisant le principe d'intersection j'ai réussis aisement a sortir les données qu'il me fallait. Pour les points 3 et 4 j'ai beaucoup plus de mal. Déja sur quelle donnée utiliser. J'ai essayer avec BDTOPO de l'IGN mais je crois qu'il y a un probleme dessus du au fait que je n'ai pas le sens du cours de l'eau. En utilisant les noeuds j'ai egalement le probleme qu'un simple noeud avec un cours d'eau en arrivé et un en sortie apparait dans la BDD avec de multiples affluents qui n'apparaissent pas sur la carte. Ou alors ils sont tous sur le même tracé. Du coup je suis parti sur les conseils d'une personne du Sandre sur la BD CARTHAGE 2017 (principalement pour le sens d'ecoulement des eaux). Mais la c'est le trou noir. Je n'arrive pas à trouver comment à partir d'un cours d'eau je peux trouver le cours d'eau suivant afin de determiner si c'est un affluent. Je suis donc preneur de toutes indications me permettant d'avancer sur mon projet. ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=364521#p364521 Pour y répondre : qgis_fr(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

2 years, 5 months

1
2
0 0

QGIS: Gestionnaire BD - Erreur 6

by qgis_fr＠georezo.net

Message posté par : Max C0L1 (niloc.max(a)gmail.com) ---------------------------------------- Bonjour J'ai eu le même problème (QGIS 3.32) Ma couche d'origine est en Polygone Dans "import de couche/fichier" du gestionnaire de base de donnée de QGIS la case à cocher "Ne pas promouvoir en multi-partie" est grisée Et ma base de donnée est sur PostgreSQL 14.9 J'ai trouvé une solution avec QGIS qui n'est pas top mais elle fonctionne : - créer une nouvelle couche en Multipolygone - utiliser "refactoriser les champs" de la boite à outils de traitements pour récupérer les champs de la couche d'origine - et faire un copier/coller le la couche d'origine vers la couche refactorisée Une fois la couche en Multipolygone l'import fonctionne ! Max ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=364522#p364522 Pour y répondre : qgis_fr(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

2 years, 5 months

1
0
0 0

QGIS: Masquer des Etiquettes sous un polygone ?

by qgis_fr＠georezo.net

Message posté par : Jean Marc MICHALOWSKI ---------------------------------------- bonjour vous pouvez choisir de n'afficher que les étiquettes qui sont localisées dans votre polygone, c'est comme cela que j'ai réglé ce même problème (mais il y a bien longtemps) de mémoire, l'application du rendu "polygone inversé" ne masque pas non plus les étiquettes jean marc ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=364419#p364419 Pour y répondre : qgis_fr(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

2 years, 5 months

1
2
0 0

QGIS: Masquer des Etiquettes sous un polygone ?

by qgis_fr＠georezo.net

Message posté par : mathieu12320 ---------------------------------------- Bonjour liege67, Comment choisir de n'afficher que "les étiquettes qui sont localisées dans votre polygone" ? Je n'ai aucun champ qui me permette de le faire. Et effectivement le polygone inversé ne masque pas les étiquettes... Merci pour votre réponse. ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=364482#p364482 Pour y répondre : qgis_fr(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

2 years, 5 months

1
1
0 0

Comment calculer le pivot d'un vecteur directement ?

by qgis_fr＠georezo.net

Message posté par : piedro25 ---------------------------------------- Bonjour, Je me permet de vous solliciter pour un problème qui relève plus de la curiosité personnelle qu'autre-chose. J'ai une table avec des vecteurs (orientés, avec une longueur définie), ce sont des façades de bâtiments. J'ai extrait l'orientation de ces vecteurs via l'outil azimut de qgis, jusque là tout va bien. J'ai ensuite appliqué un pivotage du vecteur en lui appliquant -90 degrés. Il se trouve qu'à la sortie j'ai des valeurs négatives en degrés ex : -1, quand d'autres sont justes. J'ai trouvé la solution à mon problème en distinguant par un CASE les valeurs négatives pour leur ajouter 360 et retrouver la valeur normale de l'angle. Est-ce qu'il y aurait une solution qui me permettrait de passer cette dernière étape et retrouver la valeur de l'angle entre 0 et 360 degré sans avoir à effectuer ce calcul? De mon côté j'ai essayé plusieurs fois en travaillant directement en radians de résoudre mon problème, notamment en utilisant le modulo de pi ou en soustrayant directement pi/2à la conversion en degrés, je retombe toujours sur ces mêmes valeurs. De ce que j'ai compris la fonction AZIMUTH de qgis prend pour origine le Nord qui est le point Zéro, le calcul est ensuite effectué dans le sens horaire et non dans le sens trigonométrique, mais donne quand même une valeur en radians. Voilà la formule utilisée: ----------------- Code : round (degrees(azimuth(start_point($geometry), end_point($geometry))) -90) ----------------- Mes souvenirs de math commencent à dater, mais j'aimerais néanmoins vraiment comprendre comment le contourner t comprendre comment tourne cette machine. En vous remerciant pour votre aide. ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=364488#p364488 Pour y répondre : qgis_fr(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

2 years, 5 months

1
0
0 0

2026

2025

2024

2023

2022

2021

2020

2019

Qgis_fr