python - 无法将几何图形转换为 geojson

Question

我的问题是这样的。我正在为一些数据创建模型。

class Cables(Base):
    __tablename__ = 'cables'

    id = Column(Integer, nullable=False)
    route = Column(Geometry(geometry_type='LINESTRING', srid=4326), nullable=False)

现在，我想将这样的路线转换为 GeoJSON。

我尝试过的事情

@app.route("/api/cable/<int:id>", methods=['GET'])
def get_cable(id):
    cable = session.query(Cables).filter(Cables.id == id).first()
    return str(geoalchemy2.functions.ST_AsGeoJSON(cable.route))

返回ST_AsGeoJSON(ST_GeomFromEWKB(:ST_GeomFromEWKB_1))

如果我更改返回:

return geoalchemy2.functions.ST_AsGeoJSON(cable.route)

返回TypeError:'ST_AsGeoJSON'对象不可调用

return str(cable.route)

返回0102000020e610000002000000b34fd4d9bca351c032e14d5134c240c0d24f8055e0a351c0dedea9f4dcbf40c0 这表明我确实有一个几何对象。

return cable.route

返回类型错误:“WKBElement”对象不可调用

如果我打印路线类型，

print(type(cable.route))

返回

<class 'geoalchemy2.elements.WKBElement'>

我认为它应该返回此类的对象，而不是类本身。我现在很困惑，我不知道我现在应该做什么。

有什么建议吗？

Answer 1

调用 ST_AsGeoJSON 的正确方法似乎是在查询内部。 例如，

ruta = session.query(Cables.route.ST_AsGeoJSON()).filter(Cables.id == id).first()

我最终做的是安装一个新的库(shapely)，以便读取十六进制字节串，然后将其转换为字典，因为字典可以以简单的方式转换为 json。

def ewkb_route_to_dict(self):
    """
    returns the cable's route as a dictionary.
    It uses shapely.wkb.loads. The first argument is the route as a bytestring,
    and the second one (True) is to tell loads that
    the input wkb will be represented as a hex string
    """
    return mapping(shapely.wkb.loads(str(self.route), True))

然后在我的 to string 方法中:

def __str__(self):
    d = dict()
    ...
    d["route"] = self.ewkb_route_to_dict()
    ...
    return dumps(d)

这将正确地将几何图形转换为 GeoJSON。