Message posté par : HockYan
----------------------------------------
Bonjour, je relance le sujet.
J'ai un raster complet avec des valeurs différentes et un autre avec la même emprise mais découper en fonction d'une couche de masque. Sur le raster découpé, j'ai mis la valeur 10 de partout.
Ce que je cherche à réaliser c'est fusionner les deux raster en gardant la valeur 10 sur le raster découpé et les valeurs initiales sur mon raster de base.
Donc, quand mon raster découpé est no data alors il me mets la valeur de mon raster de base. Par contre, quand il y a une valeur sur mon raster découpé alors il garde cette valeur.
Pour cela, j'utilise la calculatrice raster et je mets le code suivant :
-----------------
Code :
("monraster1@1"<= -3.40282e+38)*"monraster2@1" + ("monraster1@1" > 0)*"monraster1@1"
-----------------
La valeur -3.40282e+38 est égal à la valeur no data de mon raster découpé.
Le problème quand je réalise ce code c'est qu'il me sort que mon raster découper et pas la fusion des deux.
Je ne sais pas si c'est du à mes rasters ou au code en lui-même.
Merci pour vos réponses.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=360545#p360545
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 : Yoann Lecomte
----------------------------------------
Bonjour,
La méthode SetSpatialRef n'existe pas sur le DataSource comme vous semblez être en train de faire.
Il faut que vous alliez chercher la référence à la définition du champ géométrique de votre Layer.
Voici un exemple de code qui devrait faire ce que vous souhaitez si j'ai bien compris votre question (j'ai pris un shapefile pour l'exemple):
-----------------
Code :
ds = ogr.Open('C:/some/path/to/layer.shp')
layer = ds.GetLayer()
print(layer.GetSpatialRef()) # --> None, if no spatial ref are set for this layer
srs = osr.SpatialReference()
srs.ImportFromEPSG(4326)
srs.SetAxisMappingStrategy(osr.OAMS_TRADITIONAL_GIS_ORDER)
l_defn = layer.GetLayerDefn()
geom_defn = l_defn.GetGeomFieldDefn(0) # get the definition of the geom field
geom_defn.SetSpatialRef(srs) # apply the srs
print(layer.GetSpatialRef())
-----------------
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=360499#p360499
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 : SANTANNA
----------------------------------------
Bonjour,
Je comprends que vous passez la valeur 6.78892016869880877e-05 à PIXEL_SIZE. Exact? Il me semble que la valeur min de ce champ est 0.01 (même si la bulle dans l'interface est peu logique/erronée)
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=360539#p360539
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 : nico-29
----------------------------------------
Bonjour,
J'essaie d'exécuter le process qgis:createconstantrasterlayer via PyQGIS, avec les paramètres suivants :
-----------------
Code :
extent = QgsRectangle()
extent.setXMinimum(-2.1952652257504468)
extent.setYMinimum(48.6136243196066218)
extent.setXMaximum(-1.9701446329563943)
extent.setYMaximum(48.7326340901639128)
crs = QgsCoordinateReferenceSystem('EPSG:4326')
params = {
'EXTENT': extent,
'TARGET_CRS': crs,
'PIXEL_SIZE': 0.01,
'NUMBER': 1.0,
'OUTPUT_TYPE': gdal.GDT_Float32,
'OUTPUT': QgsProcessing.TEMPORARY_OUTPUT
}
r = processing.run('qgis:createconstantrasterlayer', params)
-----------------
si je modifie avec ma valeur (6.78892016869880877e-05), j'ai un message d'erreur :
-----------------
Code :
_core.QgsProcessingException: Incapable d’exécuter l’algorithme La valeur du paramètre PIXEL_SIZE est incorrecte.
-----------------
Nico
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=360535#p360535
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 : nico-29
----------------------------------------
Bonjour,
J'ai un message d'erreur :
-----------------
Code :
AttributeError: 'NoneType' object has no attribute 'GetSpatialRef'
Traceback (most recent call last):
File "C:\Users/ngabar/AppData/Roaming/QGIS/QGIS3\profiles\default/python/plugins\s100\S100.py", line 275, in runS102
print(layer.GetSpatialRef()) # --> None, if no spatial ref are set for this layer
AttributeError: 'NoneType' object has no attribute 'GetSpatialRef'
-----------------
Dans mon cas, je ne part pas d'un shapefile, je lit un raster avec GDAL, et si je fais un print(ds), j'ai :
-----------------
Code :
<osgeo.gdal.Dataset; proxy of <Swig Object of type 'GDALDatasetShadow *' at 0x000002154EB03180> >
-----------------
Donc pourquoi ds.GetLayer() renvoie None ?
Nico
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=360534#p360534
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 : Yoann Lecomte
----------------------------------------
Bonjour,
Je veux bien faire une passe de relecture.
Vous pouvez me contacter pour me donner les détails pratiques : yoann.cartyl(a)gmail.com.
Bonne journée!
Yoann
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=360532#p360532
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 : Vottana Tep
----------------------------------------
Bonjour à tous,
Je fais suite à ce post : QGIS Plugin: Vectorisation automatique carte ancienne
https://georezo.net/forum/viewtopic.php?pid=284157
(David, merci pour ce partage. Je venais tout juste de lire l'article de P.A. Herrault auquel tu as participé avec Martin Paegelow : Vectorisation automatique des forêts dans les minutes de la carte d'Etat-Major du 19e siècle).
Ce plugin tombe du ciel, car de nombreux organismes de protection de la biodiversité travaillent actuellement sur des projets concernant les forêt anciennes.
Mais l'extension est conçue pour la version 2.0 - 2.99 de QGIS, et est incompatible avec la 3.18, avec laquelle je travaille.
J'imagine que la réponse à ma question sera "non" mais je tente quand même, est-ce qu'une mise à jour peut être envisagée ?
Sinon, j'aurais bien aimé un retour d'expérience, est-ce que des membres de la communauté Georezo a déjà utilisé le plugin ?
Je vais installer une ancienne version de QGIS et voir ce que j'arrive à faire avec l'extension.
Merci d'avance pour vos retours,
Vottana
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=360221#p360221
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 : GlaDal
----------------------------------------
Bonjour,
Je vais vous donner plusieurs pistes, à vous de choisir laquelle répondre le mieux à votre besoin :
- Extension layer menu from project (permet de créer un onglet dans la barre de menu qui requête sur des couches présentent dans un projet, consulter la doc. C'est facile d'utilisation et facile à mettre à jour pour l'ensemble des utilisateurs).
- Dans l'onglet 'Préférences', 'Options', rubrique 'Traitement', puis dans le panneau de droite, allez dans "Menus', 'Modèles', sélectionner ici le modèle que vous avez créé, vous avez la possibilité d'ajouter un bouton dans une barre d'outils avec une icone de votre choix).
- Extension 'customize toolbar'
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=360526#p360526
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 : alefebvre (alefebvre93(a)outlook.fr)
----------------------------------------
Bonjour,
A partir d'une couche de points que j’appellerais "couche S".
Et d'un point GPS que je crée pour mon projet
Je cherche à calculer sur QGIS la distance parcourue par la route comme je pourrais le faire via google maps la plus courte
Entre mon point GPS et le point le plus proche de ma couche S de ce point GPS
Je ne vois pas comment faire ... Merci par avance pour votre aide
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=360324#p360324
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 : JaiBesoinDaide (franco81996(a)hotmail.com)
----------------------------------------
Bonnjour,
Pourriez-vous m'aider si'il vous plait.
Je dispose du plan d'un camping avec plus de 600 mobil-home, j'aimerais automatiser les chemins vers chaque mobil-home afin de pouvoir proposer un plan personnalisé à chaque client. Je n'ai jamais utilisé ce type de logiciel donc je n'ai aucune expérience pourriez-vous si'il vous plait me detailler les étapes?
Je vous remercie d'avance pour votre temps, cordialement.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=360512#p360512
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