Le portail francophone de la géomatique

Bonjour et bienvenue,

si vos données sont stockées dans une base spatialite, vous pouvez faire une sélection SQL correspondant à votre requête.

https://sites.google.com/site/sgbdspati … couvrement

Si vos données sont dans un fichier shp, je ne vois pas de solution.

Pour créer une base de donnée spatialite, vous pouvez simplement faire comme si vous alliez créer une nouvelle couche de données spatialite, . Créer cette base où bon vous semble et quitter l'assistant de création de couche.

Une fois la base de donnée créée, vous ouvrez le gestionnaire de bases de données (votre base doit apparaitre).

Vous glissez déposez vos shp depuis l'explorateur windows sur votre base de données et l'assistant d'import va s'ouvrir.

merci j'ai réussi à travailler efficacement avec QSpatialite, une extension de QGIS.

Bien sûr
voilà j’ai installé l'extension QSpatialite,
j'ai ensuite crée ma base de données en enregistrant le fichier qui la contient à l'emplacement de mon projet QGIS
j'ai importé les couches dont j'avais besoin dans la base Spatialite et
ensuite j’ai simplement rédigé ma requête SQL dans l'interface, puis quand elle m'a convenu j'ai réimporté la nouvelle table comme couche dans QGIS
A tout hasard je joins la dite requête, inspirée du lien ci dessus indiqué par MathieuB que je remercie un fois de plus:

select st_area( st_intersection(c.geometry,t.geometry)) 'aire de recouvrement',
st_area( st_intersection(c.geometry,t.geometry))/st_area(c.geometry) 'pourcentage de recouvrement',
c.'PKUID', c.'Geometry',  c.'ID_GEOFLA',  c.'CODE_COMM' , c.'INSEE_COM', c.'NOM_COMM',  c.'STATUT',  c.'SUPERFICIE',  c.'POPULATION' , c.'CODE_CANT', c.'CODE_ARR',  c.'CODE_DEPT', c.'NOM_DEPT',  c.'CODE_REG',  c.'NOM_REGION'
from  Communes3950 c,  PointsRalliement_tampon t
where
st_intersects(c.geometry,t.geometry) and
st_area( st_intersection(c.geometry,t.geometry))/st_area(c.geometry) > 0.4

Dernière modification par Chicon_power (Thu 21 August 2014 18:12)