Dernières discussions Sig forum

Syndiquer le contenu Forum SIG - Systèmes d'Information Géographique et Géomatique
Bienvenue sur le Forum SIG, systèmes d'information géographique et de la géomatique. Aide logiciels SIG, ESRI Gamme ArcGIS Desktop et Server, MapInfo, Géoconcept, GéoMedia, ArcGis, ArcView, QuantumGIS, QGIS, Programmation en Visual Basic, Avenue, Python, MapBasic, Bases de données Access, PostGIS, MySQL, Télédétection Erdas, Er Mapper, Cartographie en ligne avec MapServer, CartoWeb, OpenLayers, LizMap, GeoTools, Progiciels APIC, Star, DAO et CAO avec AutoCAD, MicroStation, Illustratot, Offre de stage, Offre d'emploi, Avis sur les formations universitaires et profesionnelles et de la Géomatique
Mis à jour : il y a 4 heures 6 min

[QGIS 2.x] Etiquetage multi-lignes d'une couche shapefile avec pyQGis

16/10/2017 - 22:39
Bonjour,
Je souhaite Ă©tiqueter sur deux lignes une couche shp avec pyQGis.
J'ai réussi à le faire avec un seul champ, mais avec deux champs ou plus (en utilisant concat) ça ne marche plus.
Merci de m'aider Ă  trouver la solution.
Ci-joint, le code que j'ai rédigé:

def labelLine(lyr, champNom1, champNom2):
# Cette fonction montre les etiquettes des champs "champNom1" et "champNom2" pour une couche lyr (couche Ligne)
label = QgsPalLayerSettings()
label.readFromLayer(lyr)
label.enabled = True
# on concate les deux champs
label.fieldName = "concat('BASE:',champNom1,'\\n','L=',champNom2 , 'm')"
label.placement= QgsPalLayerSettings.Line
label.bufferDraw= True
label.bufferSize= 1
label.setDataDefinedProperty(QgsPalLayerSettings.S ize,True,True,'8','')
label.writeToLayer(lyr)
Catégories: Géomatique, On en parle

[QGIS 2.x] Appliquer une projection Ă  un ensemble de couches

16/10/2017 - 14:32
Bonjour Ă  tous,

Actuellement (et assez souvent aussi), je travaille sur un projet qui contient entre 20 et 40 couches. J'utilise une projection qui me permet d'utiliser parfois d'anciennes données en même temps que les nouvelles sans trop de casse. Malheureusement, je passe mon temps à réattribuer cette projection choisie pour le projet et pour les différentes couches vectorielles et raster. Chaque fois que j'importe une couche, après redémarrage du logiciel, après création d'une nouvelle couche etc. Je pensais qu'il était possible de sélectionner l'ensemble de mes couches et de leur attribuer la bonne projection. Mais je n'ai pas encore trouvé comment faire. Est-ce que quelqu'un saurait comment faire une telle manipulation ?

Sinon, au pire, si vraiment il n'y avait pas d'autre moyen que de vérifier en permanence les propriétés d'une couche après l'autre, existe-t-il un moyen pour verrouiller une projection de sorte que la couche varie pas en cours de travail ou après rechargement dans le projet ?

Ib.
Catégories: Géomatique, On en parle

[QGIS 2.x] Déplacer des entités

15/10/2017 - 17:05
Bonjour Ă  tous,

Je ne trouve pas la commande ou l'extension pour déplacer une sélection d'entités. Il y a bien le bouton situé à gauche de l'Outil de noeud, mais il ne permet que le déplacement manuel. Existe-t-il un outil qui permettrait d'indiquer à Qgis qu'il faut déplacer ou copier la sélection de x mètres, suivant une direction particulière ?

Par exemple, ici, j'ai une polyligne que j'aimerais positionner à 15 m de sa position initiale, avec un angle de -39,38° (le zéro étant dirigé vers le nord).

J'ai fait quelques tentatives avec les outils standards mais pour l'instant je n'ai pas de résultat.
Catégories: Géomatique, On en parle

[Stage] recherche d'un stage de projet de fin d'Ă©tude (PFE)

15/10/2017 - 01:34
Actuellement étudiant en 3 éme année ingénierie en géomatique et topographe au sein de l'école supérieure de l'aéronautique et de technologie en Tunisie, je suis à la recherche d'une société qui pourrait me permettre d'effectuer mon stage de fin d'études d'une durée de 4 a 6 mois à partir de février ou mars de l'année 2018.
mail: nourallahbensalah@gmail.com Fichiers attachés
Catégories: Géomatique, On en parle

[SIG] Création d'un plan thématique avec des fenêtres

14/10/2017 - 15:31
Bonjour !

J'aimerais monter un petit projet personnel, c'est pourquoi je me tourne vers vous pour avoir des conseils de méthodes et applications.

J'aimerais à partir d'une base de données carto géoréférencée (comme OpenStreetMap), utiliser le plan d'une ville de manière à y insérer des ponctuels, qui feraient ouvrir (en cliquant dessus) une petite fenêtre dans laquelle je voudrais y insérer des éléments (photos raster et texte). A la manière de GoogleMaps (voir photo ci-dessous).



1 - Comment géoréférencer au préalable les photos raster ? sachant que je sais les situer chaque une sur un plan.
2 - Comment créer les fenêtres en lien avec les ponctuels comme j'ai expliqué ci-dessus?

Possible de faire ça sur un logiciel SIG ? (je suis débutant sur QGIS et Geoconcept).

Merci.
Catégories: Géomatique, On en parle

[ArcGIS 10.x] Des photos n apparaissent pas dans l'orthophotographie

14/10/2017 - 03:07
Bonjour,
Je cree mes ortho-images a l aide de Drone2Map. J'ai fait des prises photos aeriennes avec le phantom 4Pro le week end dernier. Les photos sont de bonnes qualites. Lorsque je les importe sur Drone2Map, il n'y a pas de problemes, l'overlap (frontal et lateral) est correct. Cependant apres traitement j'ai un trou dans mon ortho-image. Je ne comprends pas pourquoi.

Avez-vous une idee de la cause et comment la corriger ?

Merci
Catégories: Géomatique, On en parle

[QGIS 2.x] Problème avec le géo-référencement

13/10/2017 - 12:06
Bonjour,

Je vous soumet une bizarrerie concernant le géo-référencement avec la dernière version, bizarrerie que je n'ai jamais connue avec les versions précédentes :
Après avoir créé mes points, je coche la case "Charger dans QGIS quand terminé" et .... la couche raster ne se charge pas !
Pire encore, elle semble n'avoir été enregistrée nulle part donc si je souhaite ajouter cette couche raster et que je parcours jusqu'à l'endroit où je suis censé l'avoir enregistrée, elle n'y est pas présente :taistoi:
PS: Je joins la copie d'écran contenant les paramètres de transformation.
Si vous pouviez me dire ce qui se passe, je suis preneur....Sans titre QGIS.jpg D'avance, merci. Images attachées
Catégories: Géomatique, On en parle

[QGIS 2.x] Problème d'affichage des étiquettes avec le plugin CADASTRE

12/10/2017 - 13:49
J'ai toujours eu des soucis avec le plugin cadastre sous Windows (soucis qui allaitent nettement mieux sous OSX ou Ubuntu).
Sous QGis 2.18.13 et avec le plugin cadastre 1.4.1 et que ce soit suite à l'import dans un spatialite ou un postgresql j'obtiens le même résultat.

C'est à dire que je n'ai aucun affichage des étiquettes (numéro/noms).

J'ai utilisé des données toutes récentes dont voici le lien :
https://cadastre.data.gouv.fr/data/d...00AB01.tar.bz2

Je suis sous Windows 10 64bits et quelque soit le niveau de zoom, aucune Ă©tiquette ne s'affiche. Help ?! :confused:
Catégories: Géomatique, On en parle

[ArcGIS 10.x] Extraction de valeurs / plusieurs fichiers dbf

11/10/2017 - 18:42
Bonjour,

Je possede des tables au format .dbf qui ont le même champ unique. Dans ces tables qui ont la meme structure il y a 2 champs numeriques champ 1 / champ2. Je souhaiterai faire une requete pour extraire la valeurs max du champ1 de chaque identifiant unique (il faudrait donc parcourir l integralité des .dbf pour retenir que la valeur max) et récuperer la valeurs du deuxieme champ qui l'accompagne.


Merci du coup de main
Catégories: Géomatique, On en parle

[PostGIS] Tracer une bissectrice

10/10/2017 - 17:06
Bonjour,

Je cherche Ă  tracer des bissectrices via postgis.
J'ai utilisé le code ci-après, qui me permet de récupérer des informations quant à l'angle formé à la jointure de deux segments.


Code: -- 3.- Create segments from points and calculate azimuth for each line.
--    two calls of generate_series for a single function wont work (azimuth).

select
id,       
name,       
polygon_num,     
point_order as vertex,       

--       

case when point_order = 1         
  then last_value(ST_Astext(ST_Makeline(sp,ep))) over (partition by id, polygon_num)         
  else lag(ST_Astext(ST_Makeline(sp,ep)),1) over (partition by id, polygon_num order by point_order)       
end ||' - '||ST_Astext(ST_Makeline(sp,ep)) as lines,     
 
--       

abs(abs(       
  case when point_order = 1         
    then last_value(degrees(ST_Azimuth(sp,ep))) over (partition by id, polygon_num)         
    else lag(degrees(ST_Azimuth(sp,ep)),1) over (partition by id, polygon_num order by point_order)       
  end - degrees(ST_Azimuth(sp,ep))) -180 ) as ang
from (-- 2.- extract the endpoints for every 2-point line segment for each linestring     
        --    Group polygons from multipolygon     

select
id,             
name,             
coalesce(path[1],0) as polygon_num,             
generate_series(1, ST_Npoints(geom)-1) as point_order,             
ST_Pointn(geom, generate_series(1, ST_Npoints(geom)-1)) as sp,             
ST_Pointn(geom, generate_series(2, ST_Npoints(geom)  )) as ep     

from (
    -- 1.- Extract the individual linestrings and the Polygon number for later identification             

select
id,                   
name,                   
(ST_Dump(ST_Boundary(the_geom))).geom as geom,                   
(ST_Dump(ST_Boundary(the_geom))).path as path

-- To identify the polygon             

from poly_and_multipoly ) as pointlist ) as segments; Source : https://stackoverflow.com/questions/...-multipolygons

Ensuite je récupère une couche comprenant l'angle et les coordonnées des sommets à partir de la sélection précédente.

Maintenant, j'aimerais trouver le moyen de tracer une ligne qui coupe cet angle en deux parties égales (bissectrice) sur une distance définie et dans la bonne orientation. De ce que je comprends, il faut désormais que je fasse un calcul vectoriel pour récupérer les coordonnées d'un point X'Y' que je pourrai relier au point XY pour tracer la bissectrice mais je dois faire une erreur dans mon code :

Code: CREATE TABLE public.test_sommets_prime
AS SELECT
idu,
x+distance*sin(ang/2) as xp,
y+distance*cos(ang/2) as yp,
idu_ilot,
vertex,
ang
FROM public.test_sommets_principaux Mais le résultat n'est pas probant. Auriez-vous une aide à m'apporter svp ?

Catégories: Géomatique, On en parle

[QGIS 2.x] Plugin cadastre - erreur encodage recherche parcelle

10/10/2017 - 11:33
Salut Ă  tous !
Avec les outils de recherche et plus précisément la recherche de lieux du plugin Cadastre, je peux choisir la section mais à ce moment j'ai un message d'erreur et je ne peux pas choisir la parcelle. Pour info, j'ai pas eu de problème à l'import ni au chargement des données, et le plugin marche très bien pour un autre département. J'ai l'impression que c'est un problème de caractère spécial dans les données, mais je ne trouve rien, d'autant que le plugin gère apparemment les problèmes d'encodage...
Si quelqu'un a une idée...
Merci !

Le message d'erreur :

Une erreur est survenue lors de l'éxécution du code Python:

UnicodeEncodeError: 'ascii' codec can't encode character u'\xab' in position 20: ordinal not in range(128)
Traceback (most recent call last):
File "C:/Users/xx/.qgis2/python/plugins\cadastre\cadastre_dialogs.py", line 1887, in onNonSearchItemChoose
self.setupSearchCombobox(ckey, filterExpression, 'sql')
File "C:/Users/xx/.qgis2/python/plugins\cadastre\cadastre_dialogs.py", line 1570, in setupSearchCombobox
self.searchComboBoxes[combo]['orderBy']
File "C:/Users/xx/.qgis2/python/plugins\cadastre\cadastre_dialogs.py", line 1661, in getFeaturesFromSqlQuery
[header, data, rowCount, ok] = cadastre_common.fetchDataFromSqlQuery(connector, sql)
File "C:/Users/xx/.qgis2/python/plugins\cadastre\cadastre_dialogs.py", line 480, in fetchDataFromSqlQuery
print e.msg
UnicodeEncodeError: 'ascii' codec can't encode character u'\xab' in position 20: ordinal not in range(128)


Version de Python : 2.7.5 (default, May 15 2013, 22:44:16) [MSC v.1500 64 bit (AMD64)]
Version de QGIS : 2.18.9 Las Palmas, 3a16a4e
Catégories: Géomatique, On en parle

[QGIS 2.x] Création d'entités en fonction de la légende

10/10/2017 - 10:07
Bonjour Ă  tous,

Je me demandais s'il existait une manière simple (sans programmation) de créer une entité dans une couche, en fonction de la légende affectée à cette couche sous QGIS. Plus précisément, est-il possible de créer une entité avec l'un des champs rempli automatiquement, sur la base de la légende ?

Par exemple, si j'ai une shape de polylignes, légendée selon un champ "Type" dans lequel j'ai "Réseau d'eaux usées" et "Réseau d'eaux pluviales", est-il possible de créer des entités directement renseignées "Réseau d'eaux usées" ou "Réseau d'eaux pluviales" ?

Je maîtrise bien mieux ArcGIS et il est possible de le faire sous le logiciel d'ESRI, je serais étonné que QGIS ne le permette pas. Je pense être passé à côté de quelque chose.

Merci d'avance Ă  tous !

Etienne
Catégories: Géomatique, On en parle

[QGIS 2.x] Atlas : récupérer la valeur d'un champ en libellé d'un objet de la légende

09/10/2017 - 17:04
Bonjour Ă  tous,

Je fais un atlas sur des intercommunalités (anciennes et nouvelles) et des communes.

Mon atlas est paramétré pour utiliser chaque nouvelle interco (ensemble de règle dans le style avec $id=$atlasfeatureid) et montrer le regroupement des anciennes interco ainsi que les communes. Jusque là rien de très compliqué (quoique, j'ai activé l'option "Ne montrer que les entités à l'intérieur de l'entité de l'atlas" mais ma légende affiche plus d'occurrences que celles contenues dans l'interco "référence". Je pense qu'il s'agit de la discussion http://www.forumsig.org/showthread.p...-de-symbologie mais je ne comprends pas ce qui y est dit :roll:...).

Ma légende liste les anciennes interco mais je souhaiterais que le nom de la nouvelle interco soit récupérée automatiquement pour être mis en libellé "dynamique" dans la légende.

Ex : la nouvelle interco nommée AAA sert de délimitation pour la 1ère page de l'atlas et dans la légende, "AAA" apparait à côté du symbole représentant les contours de la nouvelle interco. Idem pour la page 2 et la récupération du nom de l'interco BBB en guise de libellé dans la légende...

Est-ce possible ?
Merci de vos retours.

EDIT
En me relisant, je me rends compte que je n'ai pas forcément été très clair...
La question est donc la suivante : au même titre que dans une étiquette je peux demander à Qgis de me renvoyer la valeur d'un champ comme [% "nom" %], est-il possible de faire la même chose avec le libellé des objets de la légende ?
Catégories: Géomatique, On en parle

[Dicogis 2.5] Problème d'ouverture du fichier xls

06/10/2017 - 11:39
Bonjour,

J'ai découvert l'outil Dicogis qui est très pratique.
Par contre, j'ai un souci, le fichier de résultat excel ne s'ouvre pas en général (à 2-3 exceptions près).

Le traitement se fait bien, sans erreur, l'outil affiche bien les fichiers trouvés.

Dans le log, pas d'erreur non plus.

Je ne sais pas pour quelles raisons le fichier a réussi à s'ouvrir 2 ou 3 fois.
J'ai essayé d'ouvrir excel (2010) avant ou après avoir lancé le traitement, a priori cela ne vient pas de là.

Quelqu'un aurait-il eu le problème et une solution?

Merci d'avance.

MFANGAIN Images attachées
Catégories: Géomatique, On en parle

[Dicogis 2.5] Problème d'ouverture du fichier xls

06/10/2017 - 11:39
Bonjour,

J'ai découvert l'outil Dicogis qui est très pratique.
Par contre, j'ai un souci, le fichier de résultat excel ne s'ouvre pas en général (à 2-3 exceptions près).

Le traitement se fait bien, sans erreur, l'outil affiche bien les fichiers trouvés.

Dans le log, pas d'erreur non plus.

Je ne sais pas pour quelles raisons le fichier a réussi à s'ouvrir 2 ou 3 fois.
J'ai essayé d'ouvrir excel (2010) avant ou après avoir lancé le traitement, a priori cela ne vient pas de là.

Quelqu'un aurait-il eu le problème et une solution?

Merci d'avance.

MFANGAIN Images attachées
Catégories: Géomatique, On en parle

[QGIS 2.x] Atlas et visualisation ou pas de symbologie

06/10/2017 - 10:03
Bonjour,
J'utilise la génération par atlas afin d'imprimer mes cartes de réseaux d'éclairage public par armoires de commande.
J'utilise donc une couche zone (polygone) qui dessine les contours de mon réseau par armoire, englobant donc les lampes (sur une couche), les réseaux aérien/sol (sur une autre couche).

Lorsque j'utilise la génération en atlas en utilisant ma couche zone, chaque page m'affiche la zone centrée d'une seule armoire.

MAIS ! :eek:

On voit également les lampes et réseaux des armoire voisines. (normal me direz vous)

Je cherche à ne pas afficher (masquer complètement) les symbologies de mes couches lampes et réseaux si elles ne sont pas explicitement dans le polygone (couche zone) servant à générer l'atlas.

1/ Comment faire cela ?
2/ Éventuellement je suis preneur d'une autre solution qui au lieu de masquer permettrait de les rendre en noir et blanc afin de bien les différencier de mes autres symbologies en couleurs.

Merci :cool:
Catégories: Géomatique, On en parle

[ArcGIS 9.x] Import d'un fichier .xls contenant du texte

05/10/2017 - 17:59
Bonjour, depuis peu j'ai un soucis avec l'import xls. Le fichier que je veux importer contient du texte dans une des colonnes mais n'est pas reconnu par ArcGIS (qui le considère en réel double). J'ai tout essayé: un format différent, forcer moi même à spécifier la nature de la colonne en "texte" mais rien y fait, ArcGIS ne veut pas me lire le texte qu'il y a dedans (<null>). Comment puis-je y remédier ? Merci d'avance ! (désolé pour le pavé, le retour à la ligne ne marche pas :/)
Catégories: Géomatique, On en parle

[ArcGIS 10.x] Besoin d'aide pour réaliser carte à partir d'un MNT

04/10/2017 - 16:49
Bonjour Ă  tous,

Comment déterminer une région montagneuse à partir d'un MNT (SRTM 30) à l'aide d'ArcGis ou autre logiciel SIG (Global Mapper ou Qgis) en se basant sur la définition suivante du doc 8168 "Exploitation technique des aéronefs Volume II):

Région montagneuse:Région à profil de terrain variable, où les changements d’altitude topographique dépassent
900 m (3 000 ft) à l’intérieur d’une distance de 18,5 km (10,0 NM).


Je ne demande pas forcément une solution toute prête, mais si quelqu'un à une idée je suis preneur.

Merci d'avance.
Catégories: Géomatique, On en parle

[ArcGIS 10.x] Besoin d'aide pour réaliser carte à partir d'un MNT

04/10/2017 - 16:49
Bonjour Ă  tous,

Comment déterminer une région montagneuse à partir d'un MNT (SRTM 30) à l'aide d'ArcGis ou autre logiciel SIG (Global Mapper ou Qgis) en se basant sur la définition suivante du doc 8168 "Exploitation technique des aéronefs Volume II):

Région montagneuse:Région à profil de terrain variable, où les changements d’altitude topographique dépassent
900 m (3 000 ft) à l’intérieur d’une distance de 18,5 km (10,0 NM).


Je ne demande pas forcément une solution toute prête, mais si quelqu'un à une idée je suis preneur.

Merci d'avance.
Catégories: Géomatique, On en parle

[ArcGIS 10.x] couleur de parcellles en fonction d'une caracteristique

04/10/2017 - 14:36
bonjour,

je souhaite mettre sur fond cadastrale les parcelles qui sont en assainissement non collectif d'une couleur différente de celles du cadastre principale. les données sont stockées sur une base de données à part.

cordialement
Catégories: Géomatique, On en parle