如何从IPython Notebook中查询存储的PostGIS光栅数据

Question

测试数据

• 1个TIF文件（159KB）

三重目标:

1. 使用raster2pgsql将光栅加载到PostGIS并在QGIS中可视化
2. 在我的IPython Notebook中连接到PostGIS并将栅格加载到NumPy数组中
3. 在我的IPyhton Notebook中，使用Pandas加载存储在PostGIS中的具有不同时间步长的一个像素光栅的时间序列

到目前为止，我已经使用< code>raster2pgsql命令将一幅光栅图像导入到PostGIS光栅中，并使用数据库管理器在QGIS中将其可视化:

raster2pgsql -s 4326 -d -I -C -M -R -l 4 D:\Downloads\raster//modis.tif -F -t 100x100 public.ndvi | psql -U postgres -d rastertest -h localhost -p 5432

但是如何从IPython笔记本中访问/查询此光栅？

我找到了这个演示文稿，它是关于SQLALchemy和GeoAlchemy2的。并且提到它也支持PostGIS栅格。好像很有意思！但是使用文档，我看不出如何将其应用于栅格数据

我想我可以使用以下代码连接到我的PostGIS数据库，其中postgres=user、密码=admin和数据库=rastertest：

from sqlalchemy import create_engine
engine = create_engine('postgresql://postgres:admin@localhost/rastertest', echo=True)

但是…任何建议都非常感谢。

Answer 1

您应该使用psycopg模块从python连接到postgres db。一些代码示例:

import psycopg2


def connect_db():

try:
    conn = psycopg2.connect("dbname='rastertest' user='admin' host='localhost' password='password'")
    conn.set_session(autocommit=True) #if you want your updates to take effect without being in a transaction and requiring a commit, for a beginner, I would set this to True
    return conn
except:
    print "I am unable to connect to the database"
    return None

def get_raster(raster_id,conn):

    query= "SELECT ST_AsText(geom) from raster_table where id={}".format(raster_id)
    conn.cursor().execute(query)
    res = cur.fetchall()

    return res[0][0]

也许您可以使用栅格的文本表示。或者，看看这里的http://postgis.net/docs/RT_reference.html，看看是否有任何函数返回您想要的numpy数组，并相应地替换get_raster中的查询。(可能是这个http://postgis.net/docs/RT_ST_DumpValues.html)