Le portail francophone de la géomatique

Bonjour,

J'ai récupéré un tutoriel python que je suis entrain de l'adapter aux besoins de mon application.

Ce script crée une fenêtre et charge une image à partir d'un shapefile avec les commandes:
def addLayer(self):
"""add a (hardcoded) layer and zoom to its extent"""

layerPath = "../data/test.shp"
layerName = "test"
layerProvider = "ogr"

Existe-t-il un moyen de de faire directement la connexion à une base postgresql/postgis et de charger toutes les layes qu'elle contient ?

Dernière modification par taupinou (Fri 02 April 2010 15:54)

Bonsoir,
j'ai réussi à afficher les données de ma base  en transformant mes tables postgis en fichiers shp grace à la commande
pgsql2shp -f <nom_fichier_shp> <nom_base_localisation> <nom_table>


mais j'ai encore deux petits soucis :

1) j'ai 4 tables dans ma base donc 4 fichiers shp sont générés. Le souci c'est que je n'arrive pas à superposer les affichages des fichiers pour avoir une seule image.
Y'aurait t'il un moyen de générer un seul fichier shp avec toutes les tables ?

Pourrais tu m'expliquer d'avantage s'il te plait, car pour le moment le programme ne reconnait pas ma table geometry_column

2) Dans son état actuel, mon script génère une fenètre et affiche l'image obtenue à partir du fichier shp dans cette fenêtre.
Je voudrais qu'à la place le sript génère directement un fichier image jpg ou png et l'enregistre dans un emplacement donné.
Est ce possible de le faire ?

Merci par avance pour vos réponses.

Dernière modification par taupinou (Fri 09 April 2010 02:46)

bonjour,

La requête sert à récupérer des données d'une table spatiale. Ce n'est pas le cas de geometry_columns.

Pour récupérer automatiquement toutes les couches il faut faire un script plus complexe avec psql :

psql -d maBase -c "SELECT f_table_name,f_geometry_column FROM geometry_columns"

Traiter la réponse en bash puis faire une boucle sur ces résultats.

Pour que cela soit générique il faut que tu prennes en compte le fait qu'une table spatiale peut avoir plusieurs colonne géométrique d'où la requête SQL qui récupère le nom de la table et le nom de la colonne spatiale.

Y.

Trouvé pour la superposition !

On définit plusieurs layers à partir des shapefiles :

layerPath1 = "../chemin/fichier1.shp"
    layerName1 = "test1"
    layerProvider1 = "ogr"

layer1 = QgsVectorLayer(layerPath1, layerName1, layerProvider1)
QgsMapLayerRegistry.instance().addMapLayer(layer1);
self.canvas.setExtent(layer1.extent())
cl1 = QgsMapCanvasLayer(layer1)

enseuite, il faut prendre les cl1, cl2 ... obtenus et faire

layers= [cl1,cl3,...cln]
    self.canvas.setLayerSet(layers)

    self.canvas.saveAsImage("/home/chemin/" + "image" + ".png")


Et voilà.

Dernière modification par taupinou (Sat 17 April 2010 16:59)