2010-11-01 18 views

Répondre

15

Utilisation SQLAlchemy 0,8, Flask-SQLAlchemy et 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)) 

Exemple de requête:

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

Une façon de convertir le résultat de coordonnées x et y (convertir en GeoJSON et extraire les coordonnées):

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

Si vous n'êtes pas limité à l'utilisation de Flask, vous pouvez essayer MapFish, basé sur Pylons et utilisant GeoAlchemy.

1

Vous pouvez l'utiliser avec Flask-SQLAlchemy, mais vous pouvez également l'utiliser avec SQLAlchemy. Il suffit de traduire le sample models from GeoAlchemy en Flask-SQLAlchemy. Quelque chose comme ceci:

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

Je n'ai pas testé le code, mais il devrait être une transcription équitable.

+2

Ce code ne fonctionne pas: il vous donne 'AttributeError: l'objet 'SQLAlchemy' n'a pas d'attribut '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__)