2

我正在尝试使用 SQLAlchemy(0.8)/GeoAlchemy2(0.2.1) 复制嵌套的原始 PostGreSQL/PostGIS 栅格查询，但不知道如何访问geomval数据类型的组件。它是一种返回“geom”和“val”的复合数据类型。

这是有效的原始查询：

SELECT (dap).val, (dap).geom
FROM (SELECT ST_DumpAsPolygons(rast) as dap FROM my_raster_table) thing

我目前正在使用的 SQLAlchemy 查询：

import geoalchemy2 as ga2
from sqlalchemy import *
from sqlalchemy.orm import sessionmaker

metadata = MetaData()
my_raster_table = Table('my_raster_table', metadata,
                  Column('rid', Integer),
                  Column('rast', ga2.Raster))

engine = create_engine(my_conn_str)
session = sessionmaker(engine)()

q = session.query(ga2.func.ST_DumpAsPolygons(my_raster_table.c.rast).label('dap'))

然后我想在子查询中访问它，如下所示：

q2 = session.query(ga2.func.ST_Area(q.subquery().c.dap.geom))

但是这种语法不起作用，或者我不会发布这个问题；）。有人有想法吗？谢谢！

4

1 回答 1

2

解决方案最终相当简单：

首先，定义一个自定义GeomvalType，继承 geoalchemy2 的CompositeType并指定一个特定于geomval的类型图：

class GeomvalType(ga2.types.CompositeType):
    typemap = {'geom':ga2.Geometry('MULTIPOLYGON'),'val':Float}

接下来，使用type_coerce将ST_DumpAsPolygons结果转换为初始查询中的GeomvalType ：

q = session.query(type_coerce(ga2.func.ST_DumpAsPolygons(my_raster_table.c.rast), GeomvalType()).label('dap'))

最后，像我之前尝试的那样从子查询中访问它（成功！）：

q2 = session.query(ga2.func.ST_Area(q.subquery().c.dap.geom))
于 2013-07-18T12:48:51.763 回答