Le portail francophone de la géomatique

Bonjour,
j'ai une solution mais elle n'est pas des plus évidentes.
Il est nécessaires d'installer Postgres/Postgis, d'importer tous tes fichiers vecteurs dans cette bdd et de faire une requête spatiale dans postgis telle qu'on peut en voir dans MapInfo à savoir une intersection entre ton rectangle et tes données vecteur.
Autre possibilité, utiliser Grass pour cela.
@+

Bonjour,

Il faut créer un Mapset tout d'abord (Grass n'utilise pas directement les couches de Quantum GIS) puis utiliser la commande d'import vectoriel dans les outils de Grass.

@+
Philippe

Bonjour noxia,

on séloigne du sujet initial mais en double-cliquant sur la table que tu veux sélectionner dans QGIS, tu ouvres l'éditeur de requête qui te permet d'appliquer un ensemble de critères à ta sélection.

Si tu as besoin de manière récurrente de ta sélection, tu peux créer une vue dans postgresql et l'afficher dans QGIS

par exemple je veux créer la vue des communes de l'Aude :

CREATE VIEW commune_aude AS SELECT * FROM commune WHERE num_insee LIKE '11%';

Buzuc, as tu essayé de simplifier tes chemins avec Inkscape (onglet chemin => simplifier ou Ctrl + L) ?

Par ailleurs, je crois que GVSIG te permet de faire des intersections entre couches (geoprocess -> intersection).


Mathieu Bossaert

Dernière modification par Mathieu BOSSAERT (Wed 10 September 2008 14:35)

Désolé ça ne répond pas à ma question. Je vais la poser ailleurs car c'est vrai que ça sort du contexte apparement.

Merci !
j'étudie tout ça, la simplification de polygones semble interessante maintenant il faut que je m'y mette...

La simplification des polygones par requete (SIMPLIFY) n'a pas l'air de faire grand chose sur la couche!

Mais c'est peut-être parce qu'il faut sauvegarder les valeurs que me donne la requête ?
En effet j'ai fait un select et cela ne fait qu'afficher le résultat de la requête non?
Mon code :

SELECT simplify(the_geom,500) as "53simpl" 
    FROM "53"

Auriez-vous une idée ?

Ca avance à la vitesse grand V !
J'ai laissé tombé la simplification de polygone : trop de déformations même avec un paramètre à 100.
J'ai en fait décomposé le fichier CLC(Corine Land Cover) en plusieurs couches postGIS. une pour les forets, une pour l'habitat urbain et une pour les cultures.
Aucun soucis pour les deux premiers : j'ai deux fichiers d'environ 5Mo parfaitement exploitables dans InkScape.
La couche des cultures demande par contre à être simplifiée (svg de 18Mo !).

Je vais déja voir ce que ça donne en faisant une intersection entre la couche et les limites de ma carte comme suggéré par Xavier78.
Le problème est que je n'arrive pas à écrire la requête qui va bien.

Pourriez-vous me donner un exemple à partir duquel je peux extrapoler ?
J'ai aussi cherché la doc des fonctions PostGis mais je ne trove que des descriptions sommaires, y'a-t-il ça quelque part ?

Encore merci à tous

J'ai testé ta requête et j'ai un message d'erreur :

La requête va renvoyer 0 ligne(s).

ERREUR:  Un SELECT * sans table spécifiée n'est pas valide

********** Erreur **********

ERREUR: Un SELECT * sans table spécifiée n'est pas valide
État SQL :42601

La même requête en ajoutant un from :

SELECT * FROM "clc53" WHERE contains("clc53".the_geom, "limites".the_geom)

me renvoi

La requête va renvoyer 0 ligne(s).

ERREUR:  entrée manquante de la clause FROM pour la table « limites »
LIGNE 1 : ... SELECT * FROM "clc53" WHERE contains("clc53".the_geom, "limites"....
                                                               ^

********** Erreur **********

ERREUR: entrée manquante de la clause FROM pour la table « limites »
État SQL :42P01
Caractère : 58

Je maîtrise assez  bien SQL mais là j'avoue que les commandes "contains" ou l'opérateur d'intersection "&&" me donnent du fil à retordre.

De plus, mes deux couches n'ont pas le même SRID : clc est en Lambert II et les limites en WGS84.
Ca devient un peu le b..... mais je pense qu'il va falloir que je transforme la couche clc en WGS84 (SRID final que j'exploiterai)

Dernière modification par buzuc (Fri 12 September 2008 15:33)

Merci Yves, je viens de comprendre ça aussi.

Par contre pourquoi est-ce que les shp de CLC n'ont pas de SRID défini ? Normalement ils devraient être en "Lambert Zone II" mais quand je les ouvre dans QGIS c'est WGS84 qui est défini !
J'ai tenté de changer le SRID en enregistrant la couche sous un autre nom mais ça n'a pas l'air de marcher...

Et je n'ai pas trouvé non plus comment le faire avec PostGIS.

Bien, bien...

J'ai donc transformer mes deux couches pour qu'elles aient le même SRID (WGS84).
J'ai lancé la requête

SELECT * FROM "clc53",limites WHERE contains("clc53".the_geom, "limites".the_geom) AND "clc53".the_geom && "limites".the_geom;

Et le résultat c'est que je n'ai aucun enregistrement renvoyé
Cela veut dire que je n'ai pas d'intersection entre mes deux couches ? C'est pas glop !