Как конвертировать SDO_GEOMTRY в GeoJSON

Question

Я работаю с sqlalchemy и geoalchemy и конвертирую свои результаты в geojson. С нормальным способом, как это:

print json.dumps([dict(r) for r in connection.execute(query)])

это невозможно, потому что cx_Oracle.Objets не сериализуем! Я могу иметь доступ через отдельные атрибуты, как этот:

result = connection.execute(query)
result2 = result.fetchone()[0]
print result2.SDO_ORDINATES

Вот моя программа:

#!/usr/bin/env python
# coding: utf8
#from __future__ import absolute_import, division, print_function
from sqlalchemy import create_engine
from sqlalchemy import Table, MetaData
from sqlalchemy.sql import and_, select
from geoalchemy import Geometry, GeometryExtensionColumn
from geoalchemy import *
from geoalchemy.oracle import oracle_functions
from geoalchemy.oracle import OracleComparator
import cx_Oracle
import json
import sdo

#def main():
engine = create_engine('oracle+cx_oracle://TEST_3D:limo1013@10.40.33.160:1521/sdetest')
metadata = MetaData(engine)

# Loading tables 
building = Table(
    'building',
    metadata,
    GeometryExtensionColumn('centroid_geom', Geometry(2, srid= 431467)),
    autoload=True,
    autoload_with=engine
)
GeometryDDL(building)

thematic_surface = Table('thematic_surface', metadata, autoload=True)
surface_geometry = Table('surface_geometry', metadata, autoload=True)
objectclass = Table('objectclass', metadata, autoload=True)

connection = engine.connect()

# define the query
query = select([(surface_geometry.c.geometry)]  #building.c.id, surface_geometry.c.geometry, objectclass.c.classname
).where(
    and_(
        building.c.grid_id_400 == 4158,
        building.c.id == thematic_surface.c.building_id,
        thematic_surface.c.lod2_multi_surface_id == surface_geometry.c.root_id,
        surface_geometry.c.geometry != None,
        thematic_surface.c.objectclass_id == objectclass.c.id,
    )
)
# Execute and print the result of the query
#print json.dumps([dict(r) for r in connection.execute(query)])
result = connection.execute(query)

Я преобразую все мои объекты cx_Oracle.Object в GeoJSON, но как? В интернете есть функция sdo2geojson, которая прекрасно работает в sql developer, но, конечно, эта функция не известна для python...

Я надеюсь, что кто-то может помочь мне???

Answer 1

Это использует (пока еще не выпущенную) версию cx_Oracle, которая поддерживает связывание объектов и другие более сложные способы использования объектов. Используя пример, предоставленный cx_Oracle для демонстрации вставки геометрии, следующий код преобразует созданный таким образом объект в JSON. Функция ObjectRepr(), включенная ниже, должна работать для любого объекта, возвращаемого из Oracle. Он просто читает метаданные об объекте и превращает объект в словарь атрибутов или список значений.

import cx_Oracle
import json

connection = cx_Oracle.Connection("user/pw@tns")
typeObj = connection.gettype("SDO_GEOMETRY")
cursor = connection.cursor()
cursor.execute("""
        select Geometry
        from TestGeometry
        where IntCol = 1""")
obj, = cursor.fetchone()

def ObjectRepr(obj):
    if obj.type.iscollection:
        returnValue = []
        for value in obj.aslist():
            if isinstance(value, cx_Oracle.Object):
                value = ObjectRepr(value)
            returnValue.append(value)
    else:
        returnValue = {}
        for attr in obj.type.attributes:
            value = getattr(obj, attr.name)
            if value is None:
                continue
            elif isinstance(value, cx_Oracle.Object):
                value = ObjectRepr(value)
            returnValue[attr.name] = value
    return returnValue

print("JSON:", json.dumps(ObjectRepr(obj)))