Verwenden Sie raster2pgsql in der Python-Konsole von QGIS?

Ich muss GeoTiFF-Dateien in eine PostGIS / PostgreSQL-Datenbank hochladen und bin neu in raster2pgsql.

Wie führe ich einen Befehl raster2pgsql in der Python-Konsole aus?

import psycopg2
import subprocess

db_name = 'enter_qgiscloud_db'
db_host = 'db.qgiscloud.com'
db_user = 'enter_qgiscloud_user'
db_password = 'enter_qgiscloud_pw'

conn = psycopg2.connect( "dbname={0} host={1} user={2} password={3}".format( db_name, db_host, db_user, db_password ) )
cursor = conn.cursor()

cmds = 'raster2pgsql -s 3857 -p -Y -I -C -M C:\qgis_cloud_data\ -F -t auto .tif | psql'
subprocess.call(cmds, shell=True)

cursor.execute(cmds)
conn.commit()

+1 für die Antwort von NettaB, wollte aber nur hinzufügen, dass Sie, wenn sich alle Ihre TIFF-Dateien im selben Ordner befinden, dies mit einem einzigen Befehl ausführen können sollten (sei es über die Befehlszeile oder über einen Unterprozess in Python):

# Set environment variables for database connection
set PGHOST=db.qgiscloud.com
set PGPORT=5432
set PGUSER=enter_qgiscloud_user
set PGPASSWORD=enter_qgiscloud_pw
set PGDATABASE=enter_qgiscloud_db

# Call the raster2pqsql utility
raster2pgsql -s 3857 -C -F -t auto C:/qgis_cloud_data/*.tif schema.target_table | psql

Dadurch wird eine neue Tabelle mit dem Namen schema.target_table erstellt und die Daten werden in diese Tabelle verschoben . Wenn Sie weitere Informationen zu den zu verwendenden Schaltern benötigen, ist diese Seite hilfreich - und enthält einige Beispiele.

Mit der Python-Implementierung müssen Sie Psycopg nur verwenden, wenn Sie SQL-Abfragen ausführen möchten. Wenn Sie nur die Daten direkt in das Programm laden, benötigen Sie nur das Dienstprogramm raster2pgsql. So kann Ihr Code angepasst werden an:

import os
import subprocess

db_name = 'enter_qgiscloud_db'
db_host = 'db.qgiscloud.com'
db_user = 'enter_qgiscloud_user'
db_password = 'enter_qgiscloud_pw'

# Set pg password environment variable - others can be included in the statement
os.environ['PGPASSWORD'] = db_password 

# Build command string
cmd = 'raster2pgsql -s 3857 -C -F -t auto C:/qgis_cloud_data/*.tif schema.target_table | psql -U {} -d {} -h {} -p 5432'.format(db_user,db_name,db_host)

# Execute
subprocess.call(cmd, shell=True)

Hier ist ein Skript, das alle TIFF-Raster in einem Ordner durchläuft und für jede Tabelle eine Tabelle mit automatischer Kachelgröße erstellt (basierend auf diesem Link ):

import psycopg2
import subprocess 
import sys, os

input_path = " C:\\qgis_cloud_data\\"
#Change to the location of pgsql bin folder
os.environ['PATH'] = r';C:\pgsql\9.6\bin'
os.environ['PGHOST'] = 'localhost'
os.environ['PGPORT'] = '9008'
os.environ['PGUSER'] = 'postgres'
os.environ['PGPASSWORD'] = 'dbpass'
os.environ['PGDATABASE'] = 'dbname'

for raster in os.listdir(input_path):
    if raster.endswith(".tif"):
       name = raster.split(".tif")[0]
       # Print the foound tiff name
       print(name)     
       raster = os.path.join(input_path, raster)                    
       # Print the full path of the tiff raster
       print(raster)
       rastername = str(name)
       rasterlayer = rastername.lower()
       conn = psycopg2.connect(database="dbname", user="postgres", host="localhost", password="dbpass", port=9008)
       cursor = conn.cursor()
       cmds = 'raster2pgsql -s 3857 -t auto "' + raster + '" |psql'
       subprocess.call(cmds, shell=True)