Frasco con GeoAlchemy código de ejemplo [cerrado]
Pregunta
Puede cualquier persona proporcionar cualquier código de ejemplo para Frasco con GeoAlchemy?
Solución
Uso SQLAlchemy 0,8, matraz de SQLAlchemy y Geoalchemy 2:
from app import db
from geoalchemy2.types import Geometry
class Point(db.Model):
"""represents an x/y coordinate location."""
__tablename__ = 'point'
id = db.Column(db.Integer, primary_key=True)
geom = db.Column(Geometry(geometry_type='POINT', srid=4326))
Consulta de ejemplo:
from geoalchemy2.elements import WKTElement
from app import models
def get_nearest(lat, lon):
# find the nearest point to the input coordinates
# convert the input coordinates to a WKT point and query for nearest point
pt = WKTElement('POINT({0} {1})'.format(lon, lat), srid=4326)
return models.Point.query.order_by(models.Point.geom.distance_box(pt)).first()
Una forma de convertir el resultado a coordenadas x e y (convertido a GeoJSON y coordenadas extracto):
import geoalchemy2.functions as func
import json
from app import db
def point_geom_to_xy(pt):
# extract x and y coordinates from a point geometry
geom_json = json.loads(db.session.scalar(func.ST_AsGeoJSON(pt.geom)))
return geom_json['coordinates']
Otros consejos
Si no está limitado a utilizar el frasco, es posible que desee dar MapFish un intento, que se basa en pilones y utiliza GeoAlchemy.
Se puede utilizar con matraz de SQLAlchemy , pero se puede utilizar con SQLAlchemy llanura también. Sólo traducir los modelos de ejemplo href="http://geoalchemy.org/tutorial.html#model-definition" de GeoAlchemy para matraz de SQLAlchemy. Algo como esto:
class Spot(db.Model):
__tablename__ = 'spots'
id = db.Column(Integer, primary_key=True)
name = db.Column(Unicode, nullable=False)
height = db.Column(Integer)
created = db.Column(DateTime, default=datetime.now())
geom = db.GeometryColumn(Point(2))
No he probado el código, pero debe ser una transcripción justo.
from myapp import db
from geoalchemy import GeometryColumn, Point
class FixXX(db.Model):
__tablename__ = 'fixXX'
fix_pk = db.Column(db.Integer, primary_key=True)
fix = db.Column(db.String)
geometry = GeometryColumn(Point(2, srid=4326))
GeometryDDL(FixXX.__table__)