2010-11-01 16 views

risposta

15

Utilizzando SQLAlchemy 0,8, Flask-SQLAlchemy e Geoalchemy 2: query

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)) 

Esempio:

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() 

Un modo di convertire il risultato a coordinate xey (convertire GeoJSON ed estrae le coordinate):

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'] 
1

Se non si è limitati all'utilizzo di Flask, si può provare a provare MapFish, che si basa su Pylon e utilizza GeoAlchemy.

1

È possibile utilizzarlo con Flask-SQLAlchemy, ma è possibile utilizzarlo anche con SQLAlchemy. Basta tradurre lo sample models from GeoAlchemy in Flask-SQLAlchemy. Qualcosa del genere:

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)) 

Non ho testato il codice, ma dovrebbe essere una buona trascrizione.

+2

Questo codice non funziona: ti dà l'attributo 'AttributeError: 'L'oggetto SQLAlchemy' non ha alcun attributo 'GeometryColumn'' – jsalonen

0
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__)