上一篇博文PostGIS匯入匯出ESRI Shapefile資料介紹瞭如何匯入空間向量資料到PostgreSQL中，緊接上一篇，本文將介紹如何使用PostGIS匯入匯出空間柵格資料。
實驗環境和上一篇的博文一樣： Ubuntu 16.04 LTS 。
資料為全球影像（ESRI ArcGIS提供的示例資料），下載連線：全球影像百度網盤下載
可以使用GDAL的gdalinfo命令檢視其詳細資訊：

在QGIS中檢視如下：

PostGIS提供了raster2pgsql工具用於柵格資料的匯入（可以使用man raster2pgsql命令檢視幫助文件）。
使用如下命令進行影像資料的插入，具體引數的含義這裡不在累贅，上篇文章有說明。關於raster2pgsql的引數可以使用man命令進行檢視。

raster2pgsql -s 4326 -C ~/Downloads/gisdata/wsiearth.tif staging.wsiearth | psql -h localhost -p 5432 -U postgres -d postgis_in_action -W

可以使用psql進入資料庫檢視：

還可以使用\d staging.wsiearth檢視關係表的結構：

可以看到有一個rid和rast的列，rid是以一個整形對插入的柵格資料進行的標示，rast列的型別是raster儲存了具體資料。

如果想要匯出資料，可以使用PostGIS提供的內建函式。其中ST_AsGDALRaster 是一個通用性的函式，可以匯出為GDAL支援的任意格式。此外，還有ST_AsPNG，ST_AsJPEG和ST_AsTIFF等函式，提供了直接匯出到指定格式，並且接受更少的引數。
注意：PostGIS的預設匯出GDAL支援的格式資料是關閉的，啟用的話需要修改PostgreSQL的配置檔案。修改的方法比較多，而且不同的PostgreSQL版本的支援情況不一樣。我使用如下命令進行修改：

sudo vim /etc/postgresql/9.5/main/postgresql.conf

在該配置檔案後面新增：

postgis.gdal_enabled_drivers = 'ENABLE_ALL'
postgis.enable_outdb_rasters = True

重啟PostgreSQL：

sudo service postgresql restart

第二個問題是：這些內建函式提供的匯出結果是PostgreSQL的bytea (byte array)資料型別，我們需要自己寫程式進行轉換成實際的影象檔案。

# -*- coding: utf-8 -*-
 


import psycopg2

# Connect to an existing database
conn = psycopg2.connect('host=localhost port=5432 user=postgres password=password dbname=postgis_in_action')

# Open a cursor to perform database operations
cur = conn.cursor()

# Execute SQL query
cur.execute("SELECT ST_AsTIFF(rast, 'LZW') AS rasttiff FROM staging.wsiearth WHERE rid=1;")

# Fetch data as Python objects
rasttiff = cur.fetchone()

# Write data to file
if rasttiff is not None:
    open('/home/theone/Desktop/wsiearth.tif', 'wb').write(str(rasttiff[0]))

# Close communication with the database
cur.close()
conn.close()

在QGIS中檢視匯出結果如下：

可以發現匯出的結果，丟失了影象的colormap，但是每個畫素值都是正確的。