一、定义空间字段：
因为该版本支持创建Geography字段类型，因此创建空间字段不同于之前版本，采用了如下形式：
geom =Column(Geometry('POLYGON’))
二、空间函数的使用：
需要注意的是所有函数都带有ST_前缀，这也是为了与postgis写法一致.
函数的使用方法：
第一种是使用func，例如在orm中你可以：
Session.query(func.ST_Area(Lake.geom))或Session.query(Lake.geom.ST_Area()) 进行面积查询
不使用orm你可以：
select([func.ST_Area(lake_table.c.geom)])
三、支持wkt和wkb格式
例如WKTElement('POINT(5 45)’)和WKTElement('POINT(5 45)',srid=4326)对应着Point类型的Geometry，和postgis中的ST_GeomFromText效果相同。
同时支持与shapely库的空间类型互相转换，详情请查阅geoalchemy2.shape.from_shape和geoalchemy2.shape.to_shape方法
空间函数整理：（谷歌机翻加上了自己的一点调整）
函数名 返回类型
2.1 geometry访问相关
（'ST_Envelope'，types.Geometry，
'返回提供的geometry的外边框，双精度（float8）’）
（'ST_GeometryN'，无，
'如果几何是``GEOMETRYCOLLECTION``，``（MULTI）POINT``，``（MULTI）LINESTRING``，``MULTICURVE``或``（MULTI）POLYGON``，``POLYHEDRALSURFACE``，则返回基于1的第N个几何， 否则返回``None``。），
（'ST_GeometryType'，无，
'返回``ST_Geometry``返回值的几何类型。'），
（'ST_IsValid'，无，
如果``ST_Geometry``返回的空间数据有效，则返回``True``。'），
（'ST_NPoints'，无，
'返回Geometry数据中的节点数。'），
（'ST_SRID'，无，
‘返回Geometry数据的SRID数值，该数值参考spatial_ref_sys表中的定义）
（'ST_X'，无，
返回点的X坐标，如果不可用，则返回“无”。输入必须是一个Point类型数据），
（'ST_Y'，无，返回点的Y坐标，如果不可用，则返回“无”。输入必须是一个Point类型数据）
（'ST_Z'，无，返回点的Z坐标，如果不可用，则返回“无”。输入必须是一个Point类型数据）
2.2 geometry编辑相关
（'ST_Transform'，types.Geometry，
'返回一个新的几何体，其坐标转换为输入的SRID'），
2.3 geometry输出
（'ST_AsBinary'，无，转化成WKB），
（'ST_AsEWKB'，无，转化成EWKB），
（'ST_AsTWKB'，无，转化成TWKB），
（'ST_AsGeoJSON'，无，'将geometry转化成GeoJSON元素返回。'），
（'ST_AsGML'，无，'将geometry转化成GML版本2或3。'），
（'ST_AsKML'，无，
'将geometry转化为KML元素返回。多个类型。默认参数：'version = 2，default precision= 15'），
（'ST_AsSVG'，无，
'返回输入geometry或geography的SVG格式数据'），
（'ST_AsText'，无，返回输入geometry的WKT格式，没有SRID信息）
（'ST_AsEWKT'，无，
'返回输入geometry的EWKT格式数据，具有SRID信息。'），
2.4 空间关系和测量
（'ST_Area'，无，如果是polygon或multi-polygon，则返回平面面积；对于geometry类型的面积单位取决于SRID。对于geography类型面积单位为平方米），
（'ST_Azimuth'，无，返回pointb和pointb之间的角度的弧度值，顺时针计算，12 = 0; 3 = PI / 2; 6 = PI; 9 = 3PI / 2），
（'ST_Centroid'，types.Geometry，返回geometry的几何中心），
（'ST_Contains'，无，A是否包含B'），
（'ST_ContainsProperly'，A是否完全包含B，不内切），
（'ST_Covers'，无，A是否包含B，可以内切，
（'ST_CoveredBy'，没有，B是否包含A），
（'ST_Crosses'，无，
如果提供的几何有一些但不是全部，则返回``True``
'内部点共同点。'），
（'ST_Disjoint'，无，
如果几何不“空间相交”，则返回``True`` - 如果'
“他们不共享任何空间。”），
（'ST_Distance'，无，
'对于几何类型返回二维笛卡尔最小距离'
'（基于空间参考）在投影单位中的两个几何之间。为''
'地理类型默认返回'之间的球形最小距离'
'米的两个地理位置。'），
（'ST_Distance_Sphere'，无，
'返回两个lon / lat几何之间的最小距离（以米为单位）。使用'
'球形地球，半径6370986米。比...快 '
'``ST_Distance_Spheroid``，但不太准确。 PostGIS版本'
'在1.5之前只实现积分。'），
（'ST_DFullyWithin'，没有，
'如果所有几何都在指定的范围内，则返回``True``
'彼此的距离'）
'ST_DWithin'，没有，
如果空间图形在另一个空间图形指定的距离内，则返回True。单位视空间坐标系而定，geography的单位为米。默认为球面，栅格数据默认为非球面
（'ST_Intersects'，无，
如果geometry/geography相交，则返回true。
（'ST_Length'，无，
'如果它是一个LINESTRING或MULTILINESTRING，则返回几何的2d长度。几何以空间参考为单位，geography以米为单位（默认球体）'），
（'ST_LineLocatePoint'，无，
'返回一个介于0和1之间的浮点数，表示LineString上最接近给定Point的点的距离占LineString总长的百分比），
（'ST_OrderingEquals'，无，
如果给定的几何表示相同的几何，则返回``True``
'和点的方向顺序相同。'），
（'ST_Overlaps'，无，
如果Geometries共享空间，则返回``True``，是相同的'
'维度，但彼此不完全包含。'），
（'ST_Perimeter'，无，
'返回ST_Surface边界的长度测量值或'
'ST_MultiSurface几何或地理。 （多边形，多边形）。 “
'几何测量以空间参考为单位，地理是'
'米。'），
（'ST_Project'，types.Geography，
'返回一个从起始点投射的``POINT``使用距离'
'米和轴承（方位角）以弧度表示。'），
（'ST_Relate'，无，
'如果此几何与空间相关'，则返回``True``
'另一个几何，通过测试内部之间的交叉点，'
'由值指定的两个几何的边界和外部'
'在intersectionMatrixPattern中。如果没有intersectionMatrixPattern'
'传入，然后返回最大的intersectionMatrixPattern'
'涉及两种几何形状。'），
（'ST_Touches'，无，
如果几何有至少一个共同点，'返回``True``
'但他们的内部不相交。'），
（'ST_Within'，没有，
'如果几何A完全在几何B'内，则返回``True``，
＃
＃几何处理
＃
（'ST_Buffer'，types.Geometry，
'对于几何体：返回表示其所有点的几何体'
'与此几何的距离小于或等于距离。 “
'计算在此几何的空间参照系中。\ n \ n'
'对于地理：使用平面变换包装。介绍1.5'
'支持不同的端盖和斜接设置来控制形状。'），
（'ST_Difference'，types.Geometry，
'返回表示几何A的那部分的几何'
'与几何B不相交''），
（'ST_Dump'，types.GeometryDump，
'返回一组geometry_dump（geom，path）行，组成''
'geometry g1。'），
（'ST_DumpPoints'，types.GeometryDump，
'返回一组所有点的geometry_dump（geom，path）行'
'构成一个几何。'），
（'ST_Intersection'，types.Geometry，
'返回表示geomA和'的共享部分的几何
“geomB。地理位置实现了对几何体的转换'
'交叉点然后转换回WGS84。'），
（'ST_LineMerge'，types.Geometry，
'返回一个（一组）LineString（s）形成的缝合在一起'
'MULTILINESTRING的组成线工作。'），
（'ST_LineSubstring'，types.Geometry，
'返回一个线串作为输入的一个子串开始和'
'以总共2d长度的给定分数结束。第二和第三'
'参数是介于0和1之间的float8值。这仅适用于'
“线串。要使用连续的MULTILINESTRING，请使用'
'与ST_LineMerge结合。'
“”
'如果\'开始\'和\'结尾\'具有相同的值，这相当于'
'到ST_LineInterpolatePoint。'），
（'ST_Union'，types.Geometry，
'返回表示'的点集联合的几何'
'的几何形状。'），
（'ST_Simplify'，types.Geometry，
'使用'返回给定几何的“简化”版本
'Douglas-Peucker算法。'），
＃
#Raster构造函数
＃
（'ST_GeomFromText'，types.Geometry，
'从OGC着名文本构造一个PostGIS ST_Geometry对象'
'表示。’），
（'ST_GeomFromEWKT'，types.Geometry，
'从OGC Extended Well-Known构造一个PostGIS ST_Geometry对象'
'text（EWKT）表示。'），
（'ST_GeomFromEWKB'，types.Geometry，
'从OGC Extended Well-Known构造一个PostGIS ST_Geometry对象'
'二进制（EWKB）表示。'），
（'ST_GeogFromText'，types.Geography，
'从众所周知的文本中返回地理对象或扩展众所周知的'
'表示。'），
（'ST_AsRaster'，types.Raster，
（'将PostGIS几何体转换为PostGIS栅格。'，'RT_ST_AsRaster'）），
＃
#Raster访问者
＃
（'ST_Height'，无，
（'返回栅格的高度，以像素为单位。'，'RT_ST_Height'）），
（'ST_Width'，无，
（'返回栅格的宽度，以像素为单位。'，'RT_ST_Width'）），
＃
#Raster Pixel Accessors和Setters
＃
（'ST_Value'，无，
（'返回给定columnx，rowy pixel或at'中给定band的值
'一个特定的几何点。乐队编号从1开始并假设为'
'如果没有指定则为1。如果``exclude_nodata_value``设置为'
'``false``，然后所有像素都包含nodata像素被认为是'
'相交并返回值。如果没有传递``exclude_nodata_value``'
'然后从栅格的元数据中读取它。'，'RT_ST_Value’））
翻译不动了，贴上原文：
('ST_Envelope', types.Geometry,
'Returns a geometry representing the double precision (float8) bounding'
'box of the supplied geometry.'),
('ST_GeometryN',None,
'Return the 1-based Nth geometry if the geometry is a '
'``GEOMETRYCOLLECTION``, ``(MULTI)POINT``, ``(MULTI)LINESTRING``, '
'``MULTICURVE`` or ``(MULTI)POLYGON``, ``POLYHEDRALSURFACE`` Otherwise, '
'return ``None``.'),
('ST_GeometryType',None,
'Return the geometry type of the ``ST_Geometry`` value.'),
('ST_IsValid',None,
'Returns ``True`` if the ``ST_Geometry`` is well formed.'),
('ST_NPoints',None,
'Return the number of points (vertices) in a geometry.'),
('ST_SRID',None,
'Returns the spatial reference identifier for the ``ST_Geometry`` as '
'defined in ``spatial_ref_sys`` table.'),
('ST_X',None,
'Return the X coordinate of the point, or ``None`` if not available. '
'Input must be a point.'),
('ST_Y',None,
'Return the Y coordinate of the point, or ``None`` if not available. '
'Input must be a point.'),
('ST_Z',None,
'Return the Z coordinate of the point, or ``None`` if not available. '
'Input must be a point.'),
#
# Geometry Editors
#
('ST_Transform', types.Geometry,
'Return a new geometry with its coordinates transformed to the SRID '
'referenced by the integer parameter.'),
#
# Geometry Outputs
#
('ST_AsBinary',None,
'Return the Well-Known Binary (WKB) representation of the geometry/'
'geography without SRID meta data.'),
('ST_AsEWKB',None,
'Return the Well-Known Binary (WKB) representation of the geometry/'
'geography with SRID meta data.'),
('ST_AsTWKB',None,
'Returns the geometry as TWKB, aka "Tiny Well-Known Binary"'),
('ST_AsGeoJSON',None,'Return the geometry as a GeoJSON element.'),
('ST_AsGML',None,'Return the geometry as a GML version 2 or 3 element.'),
('ST_AsKML',None,
'Return the geometry as a KML element. Several variants. Default '
'version=2, default precision=15'),
('ST_AsSVG',None,
'Returns a Geometry in SVG path data given a geometry or geography '
'object.'),
('ST_AsText',None,
'Return the Well-Known Text (WKT) representation of the geometry/'
'geography without SRID metadata.'),
('ST_AsEWKT',None,
'Return the Well-Known Text (WKT) representation of the geometry/'
'geography with SRID metadata.'),
#
# Spatial Relationships and Measurements
#
('ST_Area',None,
'Returns the area of the surface if it is a polygon or multi-polygon. '
'For ``geometry`` type area is in SRID units. For ``geography`` area is '
'in square meters.'),
('ST_Azimuth',None,
'Returns the angle in radians from the horizontal of the '
'vector defined by pointA and pointB. Angle is computed clockwise from '
'down-to-up: on the clock: 12=0; 3=PI/2; 6=PI; 9=3PI/2.'),
('ST_Centroid', types.Geometry,
'Returns the geometric center of a geometry.'),
('ST_Contains',None,
'Returns ``True`` if and only if no points of B lie in the exterior of '
'A, and at least one point of the interior of B lies in the interior '
'of A.'),
('ST_ContainsProperly',None,
'Returns ``True`` if B intersects the interior of A but not the boundary '
'(or exterior). A does not contain properly itself, but does contain '
'itself.'),
('ST_Covers',None,
'Returns ``True`` if no point in Geometry B is outside Geometry A'),
('ST_CoveredBy',None,
'Returns ``True`` if no point in Geometry/Geography A is outside Geometry'
'/Geography B'),
('ST_Crosses',None,
'Returns ``True`` if the supplied geometries have some, but not all, '
'interior points in common.'),
('ST_Disjoint',None,
'Returns ``True`` if the Geometries do not "spatially intersect" - if '
'they do not share any space together.'),
('ST_Distance',None,
'For geometry type Returns the 2-dimensional cartesian minimum distance '
'(based on spatial ref) between two geometries in projected units. For '
'geography type defaults to return spheroidal minimum distance between '
'two geographies in meters.'),
('ST_Distance_Sphere',None,
'Returns minimum distance in meters between two lon/lat geometries. Uses '
'a spherical earth and radius of 6370986 meters. Faster than '
'``ST_Distance_Spheroid``, but less accurate. PostGIS versions '
'prior to 1.5 only implemented for points.'),
('ST_DFullyWithin',None,
'Returns ``True`` if all of the geometries are within the specified '
'distance of one another'),
('ST_DWithin',None,
'Returns ``True`` if the geometries are within the specified distance of '
'one another. For geometry units are in those of spatial reference and '
'for geography units are in meters and measurement is defaulted to '
'``use_spheroid=True`` (measure around spheroid), for faster check, '
'``use_spheroid=False`` to measure along sphere.'),
('ST_Equals',None,
'Returns ``True`` if the given geometries represent the same geometry. '
'Directionality is ignored.'),
('ST_Intersects',None,
'Returns ``True`` if the Geometries/Geography "spatially intersect in '
'2D" - (share any portion of space) and ``False`` if they don\'t (they '
'are Disjoint). For geography -- tolerance is 0.00001 meters (so any '
'points that close are considered to intersect)'),
('ST_Length',None,
'Returns the 2d length of the geometry if it is a linestring or '
'multilinestring. geometry are in units of spatial reference and '
'geography are in meters (default spheroid)'),
('ST_LineLocatePoint',None,
'Returns a float between 0 and 1 representing the location of the '
'closest point on LineString to the given Point, as a fraction of '
'total 2d line length.'),
('ST_OrderingEquals',None,
'Returns ``True`` if the given geometries represent the same geometry '
'and points are in the same directional order.'),
('ST_Overlaps',None,
'Returns ``True`` if the Geometries share space, are of the same '
'dimension, but are not completely contained by each other.'),
('ST_Perimeter',None,
'Return the length measurement of the boundary of an ST_Surface or '
'ST_MultiSurface geometry or geography. (Polygon, Multipolygon). '
'geometry measurement is in units of spatial reference and geography is '
'in meters.'),
('ST_Project', types.Geography,
'Returns a ``POINT`` projected from a start point using a distance in '
'meters and bearing (azimuth) in radians.'),
('ST_Relate',None,
'Returns ``True`` if this Geometry is spatially related to '
'anotherGeometry, by testing for intersections between the Interior, '
'Boundary and Exterior of the two geometries as specified by the values '
'in the intersectionMatrixPattern. If no intersectionMatrixPattern is '
'passed in, then returns the maximum intersectionMatrixPattern that '
'relates the 2 geometries.'),
('ST_Touches',None,
'Returns ``True`` if the geometries have at least one point in common, '
'but their interiors do not intersect.'),
('ST_Within',None,
'Returns ``True`` if the geometry A is completely inside geometry B'),
#
# Geometry Processing
#
('ST_Buffer', types.Geometry,
'For geometry: Returns a geometry that represents all points whose '
'distance from this Geometry is less than or equal to distance. '
'Calculations are in the Spatial Reference System of this Geometry.\n\n'
'For geography: Uses a planar transform wrapper. Introduced in 1.5 '
'support for different end cap and mitre settings to control shape.'),
('ST_Difference', types.Geometry,
'Returns a geometry that represents that part of geometry A that does '
'not intersect with geometry B.'),
('ST_Dump', types.GeometryDump,
'Returns a set of geometry_dump (geom,path) rows, that make up a '
'geometry g1.'),
('ST_DumpPoints', types.GeometryDump,
'Returns a set of geometry_dump (geom,path) rows of all points that '
'make up a geometry.'),
('ST_Intersection', types.Geometry,
'Returns a geometry that represents the shared portion of geomA and '
'geomB. The geography implementation does a transform to geometry to do '
'the intersection and then transform back to WGS84.'),
('ST_LineMerge', types.Geometry,
'Returns a (set of) LineString(s) formed by sewing together the '
'constituent line work of a MULTILINESTRING.'),
('ST_LineSubstring', types.Geometry,
'Return a linestring being a substring of the input one starting and '
'ending at the given fractions of total 2d length. Second and third '
'arguments are float8 values between 0 and 1. This only works with '
'LINESTRINGs. To use with contiguous MULTILINESTRINGs use in '
'conjunction with ST_LineMerge.'
''
'If\'start\'and\'end\'have the same value this is equivalent '
'to ST_LineInterpolatePoint.'),
('ST_Union', types.Geometry,
'Returns a geometry that represents the point set union of the '
'Geometries.'),
('ST_Simplify', types.Geometry,
'Returns a "simplified" version of the given geometry using the '
'Douglas-Peucker algorithm.'),
#
# Raster Constructors
#
('ST_GeomFromText', types.Geometry,
'Constructs a PostGIS ST_Geometry object from the OGC Well-Known text '
'representation.'),
('ST_GeomFromEWKT', types.Geometry,
'Constructs a PostGIS ST_Geometry object from the OGC Extended Well-Known '
'text (EWKT) representation.'),
('ST_GeomFromEWKB', types.Geometry,
'Constructs a PostGIS ST_Geometry object from the OGC Extended Well-Known '
'binary (EWKB) representation.'),
('ST_GeogFromText', types.Geography,
'Returns a geography object from the well-known text or extended well-known '
'representation.'),
('ST_AsRaster', types.Raster,
('Converts a PostGIS geometry to a PostGIS raster.','RT_ST_AsRaster')),
#
# Raster Accessors
#
('ST_Height',None,
('Returns the height of the raster in pixels.','RT_ST_Height')),
('ST_Width',None,
('Returns the width of the raster in pixels.','RT_ST_Width')),
#
# Raster Pixel Accessors and Setters
#
('ST_Value',None,
('Returns the value of a given band in a given columnx, rowy pixel or at '
'a particular geometric point. Band numbers start at 1 and assumed to '
'be 1 if not specified. If ``exclude_nodata_value`` is set to '
'``false``, then all pixels include nodata pixels are considered to '
'intersect and return value. If ``exclude_nodata_value`` is not passed '
'in then reads it from metadata of raster.','RT_ST_Value'))
