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" для його стиснення.