我有一个多对多关系,OsmAdminUnit(多边形几何)的实例被分组到OsmAdminAgg实例中。
模型定义基本上是:
class OsmAdminUnit(db.Model):
__tablename__ = 'osm_admin'
id = db.Column(db.Integer, primary_key=True)
geometry = db.Column(Geometry(
geometry_type='GEOMETRY',
srid=3857), nullable=False)
agg_units = db.relationship('OsmAdminAgg',
secondary=aggregations,
backref=db.backref('osm_admin', lazy='dynamic'))
class OsmAdminAgg(db.Model):
__tablename__ = 'admin_agg'
id = db.Column(db.Integer, primary_key=True)
name = db.Column(db.String(100), unique=True, nullable=False)
现在我正在努力做的是选择属于某个OsmAdminAgg的OsmAdminUnits,并通过应用GeoAlchemy的ST_Union来合并polgyons。
选择属于 id = 1 的管理员 agg 的所有管理单元有效:
units = OsmAdminUnit.query.filter(OsmAdminUnit.agg_units.any(id=1)).all()
但我不明白如何ST_Union该结果。到目前为止,我的方法是:
union = db.session.query(
OsmAdminUnit.geometry.ST_Union().ST_AsGeoJSON().label('agg_union')
).filter(OsmAdminUnit.agg_units.any(id=1)).subquery()
那么,如何获得这些几何图形的联合,并将其作为GeoJSON获得呢?
顺便说一句,我正在使用SQLAlchemy,Flask-SQLAlchemy,Geoalchemy2在Flask之上构建它。
试试这个:
from sqlalchemy.sql.functions import func
union = db.session.query(func.ST_AsGeoJSON(func.ST_Union(
OsmAdminUnit.geometry)).label('agg_union')
).filter(OsmAdminUnit.agg_units.any(id=1)).subquery()
您可以在GeoAlchemy 2文档中看到此的基本模板。实质上,您需要将func
传递给查询(而不是模型)以选择联合本身。
在您的情况下,如下所示:
import sqlalchemy
union = db.session.query(
sqlalchemy.func.ST_AsGeoJSON(
sqlalchemy.func.ST_Union(OsmAdminUnit.geometry)
).label('agg_union')
).filter(
OsmAdminUnit.agg_units.any(id=1)
).all()
这将获取与过滤器匹配的 OsmAdminUnit 记录的geometry
值的并集,并将其作为字符串化的 GeoJSON 返回。
接受的答案对我不起作用,我认为导入在我使用的 sqlalchemy 版本中可能有所不同。