2010-11-01 85 views

回答

15

使用SQLAlchemy的0.8,瓶-SQLAlchemy的和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)) 

示例查询:

计算结果转换为x和y坐标的
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() 

一种方式(转换为以GeoJSON并提取坐标):

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

如果您不仅限于使用Flask,您可能需要尝试使用MapFish,这是基于Pylons并使用GeoAlchemy。

1

您可以将它与Flask-SQLAlchemy一起使用,但您也可以将它与纯SQLAlchemy一起使用。只需将sample models from GeoAlchemy翻译为Flask-SQLAlchemy即可。类似这样的:

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

我没有测试代码,但它应该是一个公平的转录。

+2

此代码不起作用:它会给你'AttributeError:'SQLAlchemy'对象没有属性'GeometryColumn'' – jsalonen 2012-04-01 21:05:41

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