Message posté par : Lsam ---------------------------------------- Bonjour, Pour ma culture personnelle, je souhaite acquérir des connaissances en Python pour manipuler des données spatiales. Pas de souci pour l'instant pour me connecter à mon serveur PostgreSQL, et pour charger des couches vecteur ou non spatiales. Par contre, je ne trouve pas de documentation claire pour charger et manipuler des couches raster issus de serveurs PostGIS. Pour l'instant, j'essaie avec le module rasterio : ----------------- Code : import psycopg2 import rasterio connection = psycopg2.connect(host="hote", database="base", user="utilisateur") curseur = connection.cursor() curseur.execute("""SELECT * FROM monschema.matable_raster ;""") # il n'y a pas moyen de charger la couche directement ? relief = curseur.fetchall() # est-ce correct ? relief2 = rasterio.open(relief, drive='PostGISRaster') # cela me renvoie le binaire du raster visiblement, donc je suis contraint d'interrompre la commande ----------------- Auriez-vous une meilleure piste pour pouvoir manipuler ma couche PostGIS comme tout objet raster ? Je veux en fait faire la même chose qu'un script sur R, où le chargement de la couche raster est très simple : ----------------- Code : library(RPostgreSQL) library(rpostgis) pilote <- dbDriver("PostgreSQL") connection <- dbConnect(pilote, host="hote", dbname="mabase", user="utilisateur") relief <- pgGetRast(connection, c("monschema","matable_raster")) ----------------- Merci ! ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=343200#p343200 Pour y répondre : geolibre(a)ml.georezo.net ou reply de votre messagerie Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement -- Association GeoRezo - le portail géomatique https://georezo.net