On en parle

[PostGIS] Requête sur un ensemble de couche présente dans un schéma

Dernières discussions Sig forum - il y a 1 heure 54 min
Salut,

J'ai une zone par table dans un schéma. Je voudrai créer une vue qui me renvoie une seule table avec l'ensemble de mes zones. Celles-ci contiennent des milliers de polygones que je fusionne via STUnion. Ensuite le but est de comparer la table théorique des zonages avec cette dernière générée à partir de x tables = au périmètre réel.

Comment les appeler automatiquement plutĂ´t que de lister Ă  la main le nom de mes 800 zones ? :confused::o:

Code: -- Dissolve
WITH dissolve_fua AS (
  SELECT
  (ROW_NUMBER() OVER())::integer AS id,
  tableoid::regclass::text AS nomtable,
  ST_Union("MA_ZONE_1".geom)::geometry(Polygon,3035) AS geom,
  SUM(ST_Area("MA_ZONE_1".geom)) AS surface
FROM
  urban_atlas_2012."MA_ZONE_1"
GROUP BY
  nomtable
)
-- Comparaison avec couche de référence des FUA
SELECT
  dissolve_fua.id,
  "MES_ZONES_DE_REF".geom,
  "MES_ZONES_DE_REF".fua_name,
  "MES_ZONES_DE_REF".country,
  "MES_ZONES_DE_REF".work_name,
  "MES_ZONES_DE_REF".commentaire,
  ST_AREA("MES_ZONES_DE_REF".geom)/ST_AREA(dissolve_fua.geom) as check_surface,
  CASE
    WHEN ST_AREA("MES_ZONES_DE_REF".geom)/ST_AREA(dissolve_fua.geom) = 1 THEN 'SURFACE OK'
    WHEN ST_AREA("MES_ZONES_DE_REF".geom)/ST_AREA(dissolve_fua.geom) != 1 THEN 'ERREUR CORRESPONDANCE'
    END
FROM
  sandbox.dissolve_fua,
  urban_atlas_2012."MES_ZONES_DE_REF"
WHERE
  SUBSTR(dissolve_fua.nomtable,19,7) = "MES_ZONES_DE_REF".fua_id;
Catégories: Géomatique, On en parle

[Autres] Problème de requêtes sous Postgresl avec les données Majic

Dernières discussions Sig forum - il y a 1 heure 58 min
Bonjour Ă  tous,

Dans le cadre de mon stage, je travail avec les données Majic d'une commune.
Ma base de données est sous Postgresql, et j'ai un soucis lors de la création d'un nouveau champ.

Exemple:

Je dois effectuer une sélection de la surface maximale de mes subdivisions fiscales (dcntsfd) par parcelle(id_par) avec le type d'occupation du sol(cgrnumd) dominant correspondant. Je cherche a crée ce champ .(cgrnumd)

Pour cela je dispose :

La table des subdivisions fiscales (suf)

Plusieurs variables:
Dcntsf, surface des subdivisions fiscales
Dcntsfd, surface dominante des subdivisions fiscales
Id_par, identifiant de parcelle
Id_suf, identifiant de la subdivision fiscale
Cgrnum, type d'occupation du sol sur la subdivision fiscale
Cgrnumd, a créer type d'occupation dominant ( sol, agriculture, terrains à bâtir...)

Important : Une parcelle peux avoir plusieurs subdivisions fiscales.


Requête testée :

SELECT distinct (id_par) id_par,MAX(dcntsf) AS dcntsfd
from suf
GROUP BY id_par
ORDER BY id_par;

Avec cette requête, j'arrive à avoir la surface de la suf la plus grande présent sur la parcelle.
Maintenant je bloque sur comment obtenir Cgrnumd, c'est à dire Cgrnum associé à la subdivision la plus grande.

Ma tentavive :

Je cherche Ă  avoir Cgrnumd, en incluant Cgrnum dans le group by, jÂ’obtiens toutes les subdivisions par parcelle et non pas seulement la subdivision la plus grande avec son type dÂ’occupation correspondant.

Je ne sais pas comment faire, merci d'avance si vous pouvez m'aiguiller.

Bonne journée à tous.

Julien

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

[Autres] Soucis de requête postgresl 9.4 sous Mac OS X (données majic)

Dernières discussions Sig forum - il y a 2 heures 5 min
Bonjour Ă  tous,

Dans le cadre de mon stage, je travail avec les données Majic d'une commune.
Ma base de données est sous Postgresql, et j'ai un soucis lors de la création d'un nouveau champ.

Exemple:

Je dois effectuer une sélection de la surface maximale de mes subdivisions fiscales (dcntsfd) par parcelle(id_par) avec le type d'occupation du sol(cgrnumd) dominant correspondant. Je cherche a crée ce champ .(cgrnumd)

Pour cela je dispose :

La table des subdivisions fiscales (suf)

Plusieurs variables:
Dcntsf, surface des subdivisions fiscales
Dcntsfd, surface dominante des subdivisions fiscales
Id_par, identifiant de parcelle
Id_suf, identifiant de la subdivision fiscale
Cgrnum, type d'occupation du sol sur la subdivision fiscale
Cgrnumd, a créer type d'occupation dominant ( sol, agriculture, terrains à bâtir...)

Important : Une parcelle peux avoir plusieurs subdivisions fiscales.


Requête testée :

SELECT distinct (id_par) id_par,MAX(dcntsf) AS dcntsfd
from suf
GROUP BY id_par
ORDER BY id_par;

Avec cette requête, j'arrive à avoir la surface de la suf la plus grande présent sur la parcelle.
Maintenant je bloque sur comment obtenir Cgrnumd, c'est à dire Cgrnum associé à la subdivision la plus grande.

Ma tentavive :

Je cherche Ă  avoir Cgrnumd, en incluant Cgrnum dans le group by, jÂ’obtiens toutes les subdivisions par parcelle et non pas seulement la subdivision la plus grande avec son type dÂ’occupation correspondant.

Je ne sais pas comment faire, merci d'avance si vous pouvez m'aiguiller.

Bonne journée à tous.

Julien

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

[ArcGIS 10.x] Quelle extraction de coordonnées sur polygone selon format légal et/ou technique, ..?

Dernières discussions Sig forum - il y a 2 heures 25 min
Bonjour Ă  tous,

Je travail en Afrique pour le classement d'une zone en aire protégée. Il me faudrait avoir accès aux coordonnées des divers polygones que constitue cette nouvelle aire protégées. Seulement il ne s'agit pas d'extraire seulement les sommets car il serait trop nombreux, d'où mes deux prochaines questions:

1) Quels sont les standards utilisés par les juristes, géomètres ou autres professionnel devant faire état d'une série de coordonnée de polygone sur un document technique, légal ou autre? Exemple complètement au hasard: tous les km a partir de l'extrémité la plus au Nord?, tous les 2 km? Les extrémités NE, N, NO, SO, S, SE et tous les 500 mètres?
J'espère que cette question est déjà assez clair. Les procédures légales dans ce genre de cas n'existe pas dans le pays où je me trouve (pour s'en inspirer), il nous faut régulièrement devoir improviser en nous basant sur d'autres standard internationaux.

2) Une fois que j'ai compris quel serait la meilleur approche d'extraction, comment pouvoir techniquement le faire sur logiciel Arcgis (ou même Qgis)? En effet, si j'ai un critère de distance, cela va compliquer les manipulations. Peut-être existe t-il un outil sur Arcgis?

Merci Ă  tous!

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

[QGIS 2.x] Enveloppe concave et couche avec trop de points

Dernières discussions Sig forum - il y a 3 heures 17 min
Bonjour !

Je cherche à réaliser une enveloppe concave à partir de pieds de vigne pour délimiter mon champ de vigne. J'ai donc trouvé l'extension Concave Hull dans QGIS qui à l'air de fonctionner plutôt bien avec des couches test de quelques points.
Souci lorsque je lance l'outil sur mes couches de travail, il a chaud et m'affiche un avertissement du genre "vous avez 21 435 entités, ca peut être long" voir très long... En fait je n'ai pas réussi à attendre jusqu'à la fin.
J'ai besoin d'une méthode beaucoup plus rapide...

Alors je me suis dit, si on délesté la couche d'un maximum de points dans l'intérieur de mon champ, et qu'on fasse tourner l'outil sur la version épurée... Problème je ne vois pas comment supprimer ces points.

Est-ce que quelqu'un aurait une idée lumineuse ?
Catégories: Géomatique, On en parle

[PostGIS] Polyligne, points et RiverGIS

Dernières discussions Sig forum - il y a 5 heures 31 min
Bonjour,

Etant en phase d'initiation à la gestion de base de données (POSTGIS), je rencontre un premier problème après la création de ma base de données.

Je n'arrive pas à importer ma couche de lignes (.shp) via pgAdmin 4 après avoir créer ma base de données. J'avais essayé sur un autre ordinateur au bureau, mais avec pgAdmin III, et j'ai eu un problème d'importation : geometry ? champs x y à créer obligatoirement ? Sur des lignes c'est pas évident..

Merci de me donner quelques pistes s'il vous plaît !
Catégories: Géomatique, On en parle

[ArcGIS 10.x] Attribut géométrique Builder

Dernières discussions Sig forum - il y a 5 heures 53 min
Bonjour, Je rencontre un problème avec le modèle Builder d'Arcgis (10.4.1) lors de la création d'attribut géométrique. Les outils déjà intégré dans des modèles builders ne marchent plus et il m'est impossible d'en créer de nouveau en utilisant l'outil (peut importe le type de données que j'indique, leur projection, les unités que je demande) J'ai fait plusieurs tests et pour le moment je sèche, du coup je me demande si d'autre personne rencontre ou non le problème Merci d'avance Mathieu Koziar
Catégories: Géomatique, On en parle

[PostGreSQL] Questions :Jointures, requĂŞtes Ă  partir d'une autre etc...

Dernières discussions Sig forum - il y a 6 heures 43 min
Bonjour,

Je souhaite aujourd'hui ajouter un champ à une table de ma BDD à partir d'une sélection sur une autre table.
Mais je me rends compte que j'ai du mal Ă  assimiler la doc.
Je me permets donc de vous solliciter afin de m'Ă©clairer.

Je ne saisis pas les différences entre les différents types de jointures.
Left join, right join, inner join etc...
A quoi servent-ils?
Dans quels cas préférer l'un plutôt que l'autre?

Qu'en est il du mot clé USING?

Une autre question concerne l'utilisation d'une requête incluse dans une autre, quelle est la syntaxe/règle générale?
Catégories: Géomatique, On en parle

[Autre] Problème d'affichage des couches avec Mapguide Maestro

Dernières discussions Sig forum - 21/05/2017 - 16:36
Bonjour,
J'ai envie de créer un interface cartographique pour afficher des des parcelles de terrain. J'utilise à cet effet MapGuide et Maestro. J'arrive bien à installer le logiciel, même si je reçois après connexion un message me disant que la "version de maestro est plus récente que celle du site sur lequel elle a été testée", Mais quand j'essaie d'afficher les couches ça maestro ne m'envoie pas dans le navigateur web comme sur le tutoriel que je suis, il m'ouvre plutôt une fenêtre qui affiche le message d'erreur suivant : "Input string was not in a correct format". Voici un aperçu du problème en image. Je travaille sur un Laptop sur lequel est installé Windows 7 professionnal en 32 bits. La version du serveur Mapguide est 2.6 et la version du client Maestro 5.1. Je remercie tous ceux qui voudront bien m'aider parce que je suis dans l'urgence d'une solution
Catégories: Géomatique, On en parle

[QGIS 2.x] Georéférencement des données raster avec QGIS

Dernières discussions Sig forum - 19/05/2017 - 19:08
Bonjour chers tous,

Je voulais solliciter de votre aide afin de connaitre"comment faire un georéférencement des données raster avec QGIS". Je veux si possible avoir la procédure à suivre pour le faire.

Cordialement, Guy.
Catégories: Géomatique, On en parle

[PostGreSQL] Update de table A avec doublons en fonction de table B - PostgreSQL

Dernières discussions Sig forum - 19/05/2017 - 09:14
Bonjour,

J'aimerais pouvoir mettre Ă  jour les lignes doublons d'une table A (ici le champ proprio) en fonction d'un champs d'une table B.
Je pense avoir balayé pas mal de forums avant de venir expliquer mon problème ici, mais je n'ai rien trouvé qui réponde à mon besoin.

Je travaille en SQL sur PG Admin 9.3.

Auriez vous des pistes de réflexion, syntaxe à me proposer?

Merci d'avance pour toutes les réponses que vous pourrez m'apporter.2017-05-18_14h57_31.png Images attachées
Catégories: Géomatique, On en parle

[QGIS 2.x] Aide pour création d'une fonction en python

Dernières discussions Sig forum - 19/05/2017 - 08:03
Bonjour, j'ai besoin d'aide pour me créer une fonction en python sous QGIS 2.18.7
J'ai 2 couches :
Armoire (CODE, Nom etc.)
Lampes (CODE,NUMERO,PUISSANCE etc.)

Nota : CODE est un STRING

Je souhaite créer un champ virtuel (sur la couche Armoire) mis à jour automatiquement afin de compter le nombre de lampes qui y sont raccordées.
Le but étant que lors de l'affichage des informations sur une armoire, j'ai un décompte des lampes qui lui sont rattachées.

La fonction dirait ceci :

Retourne moi/Compte moi le nombre de lampes (de la couche Lampes) ayant le mĂŞme CODE que l'armoire

D'avance merci
Catégories: Géomatique, On en parle

[QGIS 2.x] [Couleur du fond] Trame

Dernières discussions Sig forum - 18/05/2017 - 16:50
Bonjour,

Je souhaite modifier la couleur du fond dans les propriétés du projet.
Pour mes besoins, je souhaite que le fond soit non-pas une couleur mais une trame en hâchures (voire des symboles)...

La couleur du fond ne propose que des couleurs, est-il possible de disposer d'une trame de fond ???

Par avance merci

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

You Should Know The Benefits To Join Alibaba Clone Script

Dernières discussions Sig forum - 18/05/2017 - 14:41
Modern times bring modern techniques to buy, sale, earn and everything else. In this all Alibaba Clone Script is a unique script that provides all in one. Yes it provides all services in a single website. It is a Multi script Marketplace, where one can sale and buy any no. of products. By advertising your product here you are advertising it in whole world because Alibaba Clone Script is an international company. It will release your advertisement on international level. You can look for personal, domestic, professional, pet, baby, electronics and all other products required to a human being. It also provides B2B and B2C marketplace services. So One can start his base from here by joining. It is a platform for beginners as well as professionals. Beginners will find it a platinum platform as Alibaba Clone Script is a several years old website with higher experience and well qualified skills. It will make you learn the whole techniques of e-business. Professionals will found it a diamond opportunity as they can make their business universal popular by joining. Its B2B business supports a business with other businessmen in global. And B2C business is a business where individuals make a deal of minimum one and maximum unlimited products. E-payments are very secure and safe. In case the product you purchase from here, you found it defected, you can change it also. In order to make an order you must above 15 years old. Home services are free of cost. The products and services are supplied to all over world. You can also buy a website , business, software from here. People who want to join a membership several packages are available. You can join according to your desires. Discount coupons are announced by the category of user. Special discount coupons are announced for AlibabaÂ’s regular customers. People who want to join this business or start their own business get a better support from Alibaba Clone Script. So it is a complete package of web services. You can access it from anywhere. An order can also be taken by anytime because it is a 24 hour working website. If you are a new user and not satisfied with the services, You can also concern to admin for that or you have also a better option of feedback. By joining account, You can also concern to its members for any kind of knowledge.

For purchasing alibaba clone script one can visit the website:
http://www.b2bbusinessdirectoryscrip...script-pro.htm
Catégories: Géomatique, On en parle

[GeoServer] API REST - Accéder aux infos d'une couche

Dernières discussions Sig forum - 18/05/2017 - 14:23
Bonjour,

Je voudrais écrire un script PHP qui récupère via l'API REST les infos d'une couche Geoserver, notament le titre et l'abstract.
J'ai réussi à récupérer les info de l'entrepot, avec la requête :
Code: $url = "http://localhost:8080/geoserver/rest/workspaces/monworkspace/datastores/mondatastore/featurestype"; Ca me récupère, le nom, la description, le type de datastore, ...
Mais quelle est la requĂŞte pour descendre au niveau en dessous, c'est Ă  dire au niveau de la couche ?

Merci,
Nico

- - MISE A JOUR - -

Bon j'ai trouvé, c'est :
Code: $url = "http://localhost:8080/geoserver/rest/workspaces/monworkspace/datastores/mondatastore/featuretypes/macouche"; Mais entre les featurestype et les featuretypes, ya de quoi se mélanger les pinceaux !

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

[QGIS 2.x] Problème enregistrement simultané sur une couche PostGIS

Dernières discussions Sig forum - 18/05/2017 - 13:43
Salut,

J'ai deux utilisateurs QGIS 2.18.7 qui édite la même couche PostGIS dans des zones bien distinctes d'une agglomération.
Malheureusement, de manière aléatoire que ce soit pour 3 ou 15 enregistrements, j'ai le message d'erreur suivant :

Code: Impossible de valider les changements pour la couche MACHINBIDULE
Erreurs :
Erreur du fournisseur de données :
Erreur PostGIS lors de la modification d'attribut : ERREUR: l'instruction préparée « updatefeature » existe déjà Des idées ?
Catégories: Géomatique, On en parle

[QGIS 2.x] Convertir une image en 3 bandes de couleurs en raster classifié

Dernières discussions Sig forum - 18/05/2017 - 11:14
Bonjour,

Je dispose de données relatives aux feux de forêt sous forme d'image géo-référencées. Les variables sont représentées par classes de couleurs, l'image étant en bande de trois couleurs RVB.

Je souhaiterais classer les indices de 1 à n correspondant aux classes de couleurs et récréer un raster "analysable".

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

[ArcGIS 10.x] Python - Personnaliser paramètres Boite à outil personnalisée

Dernières discussions Sig forum - 17/05/2017 - 15:59
Bonjour,

Venant du monde logiciel libre, j'essaie de reproduire l'équivalent d'un Plugin QGIS avec une toolbox personnalisée.
Je n'ai pas trouvé réponse à ma question après avoir fouillé un peu partout, mea culpa si la question a déjà été posée...

Contexte :
Création d'un script intégré dans une tooblox perso

Problème :
Ne proposer à l'utilisateur de ne choisir QUE les couches présentes dans le document ,sans avoir à naviguer dans les dossiers pour le récupérer (avec une liste déroulante ou des cases à cocher bref, tout sauf le navigateur de dossier) :

Capture1.PNG
Capture2.PNG


Est-ce possible?
Merci d'avance pour votre aide!


-----------------------------

Pour aller plus loin :

L'objectif final est de lister les valeurs uniques présentes dans l'un des champs de la couche sélectionnée, un peu comme lorsque vous faites une requête SQL et que vous cliquez sur "Get Unique Values" après avoir sélectionné un champ. L'idée serait de lister les différentes valeurs uniques de ce champ afin que l'utilisateur puisse sélectionner celles qui l’intéresse (cases à cocher ou autre).

Capture3.PNG

En soit je pense qu'en personnalisant une requête SQL je dois pouvoir contourner le problème, mais c'est pas fou.. :roll: Images attachées
Catégories: Géomatique, On en parle

[QGIS 1.x] Problème de SCR avec OpenLayers Plugin

Dernières discussions Sig forum - 17/05/2017 - 13:39
Bonjour,

Je me permets de vous contacter n'ayant pas trouvé de solutions à mon problème. Je suis tout juste initié à Qgis et je pense que mon problème est tout bête mais je tourne en rond depuis un bon moment.

Pour résumer : Mon objectif est de superposer un Shapefile avec un fond de carte Google Satellite sur OpenLayers.
Capture dÂ’e?cran 2017-05-17 a? 12.48.22.jpgVoici le Shapefile.
Le SCR est WGS84-EPSG:4326

Lorsque je rajoute le fond de carte voici ce qu'il se passe :
Capture dÂ’e?cran 2017-05-17 a? 12.47.04.jpg
Le SCR est WGS84/Pseudo Mercator-EPSG:3857

J'imagine donc qu'il y a un problème de compatibilité dans le SCR mais lorsque je le modifie pour le fond de carte ou pour mon Shapefile j'ai des erreurs Python. Que faire?

J'espère que c'est assez clair. Je vous remercie d'avance pour votre aide !!

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