Message posté par : YasmineB
----------------------------------------
Update : Pour le Lambert 2 étendu j'ai exécuté cette commande dans la calculatrice de champs :
x(transform($geometry, 'EPSG:2154', 'EPSG:27572')).
Par contre pour les coordonnées en WGS84 en degré minute seconde je sèche... Je suis preneuse si vous avez des astuces !
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332228#p332228
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
Message posté par : lylambda
----------------------------------------
Bonjour,
Ayant un besoin similaire à Pascal153, je me permets de déterrer ce fil.
J'utilise une action python qui importe une dalle raster d'un Modèle Numérique de terrain (MNT, à savoir le RGE-ALTI) à partir d'un quadrillage Shapefile.
Par défaut, le style d'une dalle de MNT se base sur un dégradé de gris en fonction de l’élévation minimum et maximum au sein de la dalle.
Ainsi, en important plusieurs dalles, chacune aura son propre dégradé de gris et l'affichage globale ne sera pas cohérent.
Je souhaite donc que mon action python importe la dalle ET applique un style à la volée à partir d'un qml que j'ai préparé.
C'est sur cette seconde partie que je sollicite une aide du forum.
Deux alternatives existent à la solution que je recherche. Je les ai écarté pour l'instant :
- Copier et renommer mon qml autant de fois que j'ai de dalles. Vu leur nombre, cela ne me semble pas une bonne solution.
- Réaliser un VRT. Cette option est trop lente à grande échelle pour un usage efficace.
Voici le script en question :
-----------------
Code :
import os
def getVectorLayerByName(NomCouche):
layermap=QgsProject.instance().mapLayers()
for name, layer in layermap.items():
if layer.name()==NomCouche:
if layer.isValid():
return layer
else:
return None
mypath="[% "DOSSIER" %]".replace("\\","/")+ "[% "NOM2" %]"
instRegistry = QgsProject.instance()
macouche = instRegistry.mapLayersByName( "[% "NOM2" %]")
courant=qgis.utils.iface.mapCanvas().currentLayer()
raster_ouvert=getVectorLayerByName("[% "NOM2" %]")
if raster_ouvert is not None:
QgsProject.instance().removeMapLayer(raster_ouvert.id())
qgis.utils.iface.mapCanvas().refresh()
else:
qgis.utils.iface.addRasterLayer(mypath,"[% "NOM2" %]")
qgis.utils.iface.setActiveLayer(courant)
-----------------
Cette action permet donc d'ouvrir ou fermer une dalle raster en fonction d'un carroyage vecteur et des Champs virtuels "DOSSIER" et "NOM2". Elle est idéale pour l'Orthophoto HD, qui est trop lourde pour en faire un VRT (https://georezo.net/forum/viewtopic.php?pid=321243#p321243 (cf ce message) ).
Pour l'utiliser sur les dalles du RGE-ALTI, je souhaite ajouter une fonction supplémentaire : l'application d'un style à la volée (stocké en qml) à la dalle ouverte .
En fouillant le web, je trouve des indications sur l'utilisation de fonction tel que "layer.loadNamedStyle(CHEMIN/STYLE.qml)". Mais je n'ai pas su les mettre en application, mes essais se concluant par de belles erreurs.
Un habitué de python saurait-il m'aider ?
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332256#p332256
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
Message posté par : Augustin_P
----------------------------------------
Bonjour,
Je souhaiterais extraire, sans modifier la forme, l'ensemble des polygones d'une couche qui sont compris, partiellement ou intégralement, dans les polygones d'une autre couche.
Pour expliquer le contexte. J'ai le cadastre sur QGIS, parcelles.shp, et j'ai défini une zone d'étude, site.shp. Je souhaiterais extraire toutes les parcelles (polygones) qui sont comprises partiellement ou intégralement dans ma zone d'étude (polygone) et cela sans couper les parcelles.
Ci-joint un exemple ce que je souhaiterais faire.
Merci.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332249#p332249
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
Message posté par : JP Beaudoin
----------------------------------------
bonjour
je suis sur qgis 3.6.0
je tests présentement et j'essaie de mettre une infobulle d'une photo
regarder les infos dans la pièces jointe
si mon code n'est pas bon .. retranscrire svp
je suis peu familier avec ce type de code
merci!
----------------------------------------
Ce message est accompagné de fichiers, pour les télécharger, suivre le lien ci-dessous.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332253#p332253
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
Message posté par : SALES P
----------------------------------------
Bonjour
Je cherche à générer un raster virtuel pour un SCAN 25 sur un département. Je dispose de fichiers .jp2
Lorsque je lance la construction en pointant le répertoire source et en indiquant un nom de fichier en sortie, j'obtiens une couche toute noire !!
Faut-il cocher certains options de construction dans le panneau de commande ?
Merci de votre aide
P SALES
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332247#p332247
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
Message posté par : Sylvain - GeoSigNature
----------------------------------------
Cela doit être possible en SQL spatial.
Voir par exemple ce sujet (via PostGis) : [PostGIS] : Fusionner les parcelles qui se touchent (https://georezo.net/forum/viewtopic.php?pid=227926)
Il "suffit" d'ajouter une condition "WHERE" pour les conditions attributaires.
Si vous êtres partants pour la solution "SQL", vous pouvez peut-être déplacer dans la section "Géo'BD" où vous pourrez croisez quelques pointures du SQL Spatial ;)
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332238#p332238
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
Message posté par : Sylvain - GeoSigNature
----------------------------------------
Bonjour,
Il est nécessaire que l'ensemble des fichiers décompressés soient dans le même dossier.
Vous devez, a minima, avoir au moins un fichier d'extension SHP, DBF, PRJ (+ SHX, CPG, ...), et ceux-ci doivent avoir le même nom.
Quelles sont les extensions des fichiers décompressés ?
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332222#p332222
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
Message posté par : lylambda
----------------------------------------
Bonjour Yugath,
As-tu tenté de décompresser les sous dossiers/fichiers qui présentent l'extension shp.7z ?
Il est possible que toute l'arborescence soit compressée, comme des poupée gigogne.
Sinon, tu peux tenter de supprimer manuellement l'extension .7z si c'est seulement un souci de nommage. Mais j'y crois moins.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332219#p332219
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
Message posté par : Vincent
----------------------------------------
Ça devrait être la même opération, mais le résultat ira dans un champ textuel.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332232#p332232
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
Message posté par : Vincent
----------------------------------------
-----------------
Citation :
Il faudrait que j'ordonne de compter et trouver le meilleur moyen de faire un groupe de 28. A condition bien sûr que ces regions soient côte à côte ou séparées par un océan.
-----------------
Ce n'est vraiment pas clair. Si vous pouvez rephraser, vous aurez plus de chances d'avoir d'autres réponses. Merci.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332230#p332230
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