rt_raster_to_gdal: не удалось загрузить выходной драйвер GDAL
Что касается первой ошибки с ST_AsTIFF , вам необходимо включить драйверы GDAL, которые по умолчанию не включены для PostGIS 2.1. Смотрите руководство о том, как это сделать. Например, у меня есть переменная окружения, настроенная на компьютере Windows с:
POSTGIS_GDAL_ENABLED_DRIVERS=GTiff PNG JPEG GIF XYZ DTED USGSDEM AAIGrid
что может быть подтверждено PostGIS с помощью:
SELECT short_name, long_name
FROM ST_GDALDrivers();
PostGIS для Numpy
Вы можете экспортировать вывод в файл GeoTIFF виртуальной памяти для чтения GDAL в массив Numpy. Для подсказок на виртуальных файлах, используемых в GDAL, см. Этот пост в блоге .
import os
import psycopg2
from osgeo import gdal
# Adjust this to connect to a PostGIS database
conn = psycopg2.connect(...)
curs = conn.cursor()
# Make a dummy table with raster data
curs.execute("""\
SELECT ST_AsRaster(ST_Buffer(ST_Point(1, 5), 10), 10, 10, '8BUI', 1) AS rast
INTO TEMP mytable;
""")
# Use a virtual memory file, which is named like this
vsipath = '/vsimem/from_postgis'
# Download raster data into Python as GeoTIFF, and make a virtual file for GDAL
curs.execute("SELECT ST_AsGDALRaster(rast, 'GTiff') FROM mytable;")
gdal.FileFromMemBuffer(vsipath, bytes(curs.fetchone()[0]))
# Read first band of raster with GDAL
ds = gdal.Open(vsipath)
band = ds.GetRasterBand(1)
arr = band.ReadAsArray()
# Close and clean up virtual memory file
ds = band = None
gdal.Unlink(vsipath)
print(arr) # this is a 2D numpy array
Показывает растеризованную буферную точку.
[[0 0 0 1 1 1 1 0 0 0]
[0 1 1 1 1 1 1 1 1 0]
[0 1 1 1 1 1 1 1 1 0]
[1 1 1 1 1 1 1 1 1 1]
[1 1 1 1 1 1 1 1 1 1]
[1 1 1 1 1 1 1 1 1 1]
[1 1 1 1 1 1 1 1 1 1]
[0 1 1 1 1 1 1 1 1 0]
[0 1 1 1 1 1 1 1 1 0]
[0 0 0 1 1 1 1 0 0 0]]
Обратите внимание, что я использовал формат 'GTiff' в примере, но другие форматы могли бы быть более подходящими. Например, если у вас большой растр, который необходимо передать по медленному интернет-соединению, попробуйте использовать «PNG» для его сжатия.