¿Utilizar raster2pgsql en la consola Python de QGIS?

Question

Tengo que cargar archivos GeoTiFF en una base de datos PostGIS/PostgreSQL y soy nuevo en raster2pgsql.

¿Cómo puedo ejecutar un comando raster2pgsql en la consola de Python?

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()

Answer 1

Aquí hay un script que itera sobre todos los rasters tiff en una carpeta y crea una tabla con el tamaño de mosaico automático para cada uno (basado en esto enlace ):

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)

Answer 2

+1 por la respuesta de NettaB, pero sólo quería añadir que si todos tus archivos tiff están en la misma carpeta deberías ser capaz de lograr esto usando un solo comando (ya sea desde la línea de comandos o a través de un subproceso en 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

Esto creará una nueva tabla llamada esquema.tabla_destino y empujar los datos en él. Si necesitas más información sobre los interruptores a utilizar, esto página es útil, e incluye algunos ejemplos.

Con la implementación de Python no necesitas usar Psycopg a menos que planees ejecutar consultas SQL - si sólo vas a cargar los datos directamente sólo necesitas la utilidad raster2pgsql. Así que su código puede ser adaptado a:

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)

Answer 3

Para aquellos que tienen el mismo problema en Ubuntu, Aquí está la solución que funcionó para mí, gracias a @NettaB 's solution :

import psycopg2
import subprocess
import os

input_path = "/home/projects/outputs/"
#Change to the location of pgsql bin folder
os.environ['PATH'] = r'/usr/bin'       # path to the raster2pgsql file
os.environ['PGHOST'] = 'localhost'     # localhost server 
os.environ['PGPORT'] = '5432'          # port (default is this)
os.environ['PGUSER'] = 'postgres'      # postgres username
os.environ['PGPASSWORD'] = '1234'      # my password
os.environ['PGDATABASE'] = 'p3'        # my database's name

for raster in os.listdir(input_path):
    print(raster)
    if raster.endswith(".tif"):
       name = raster.split(".tif")[0]
       # Print the found 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="p3", user="postgres", host="localhost", password="1234", port=5432)
       cursor = conn.cursor()
       cmds = 'raster2pgsql -s 3857 -t auto "' + raster + '" |psql'
       subprocess.call(cmds, shell=True)