postgis基础应用

江城老金

已于 2022-11-09 09:07:28 修改

阅读量597

点赞数 1

分类专栏：数据库文章标签： postgresql

于 2021-04-26 10:40:13 首次发布

本文链接：https://blog.csdn.net/joshho/article/details/116145297

版权

数据库专栏收录该内容

8 篇文章 1 订阅

订阅专栏

部分函数操作示例

添加空间字段

SELECT AddGeometryColumn ('someschema','sometable','somecolumn', 4490,'POINT', 2);

设置geometry字段的值

update sometable set shape=ST_GeomFromText('POINT('||longitude||' '||latitude||')', 4490) where xxx=xxx;

更新geometry字段SRID

select UpdateGeometrySRID('some_table_name', 'shape', 4490)

SRID简介

SRID WKID 空间参考简介： https://blog.csdn.net/bigbigtreewhu/article/details/52162277

postgis一些管理函数

1、AddGeometryColumn

添加一个几何类型字段到已经存在的可用表中，默认情况下使用类型修饰符定义约束。该i函数实现版本较早，支持3D但不能删除z，也支持添加圆和曲面。
函数：

text AddGeometryColumn(varchar table_name, varchar column_name, integer srid, varchar type, integer dimension, booleanuse_typmod=true);
text AddGeometryColumn(varchar schema_name, varchar table_name, varchar column_name, integer srid, varchar type, integer dimension, boolean use_typmod=true);
text AddGeometryColumn(varchar catalog_name, varchar schema_name, varchar table_name, varchar column_name, integer srid, varchar type, integer dimension, boolean use_typmod=true);

注意事项：
 （1）srid必须在SPATIAL_REF_SYS表中
 （2）添加到的表必须存在
 （3）添加的类型存在

例如：

CREATE SCHEMA my_schema;
CREATE TABLE my_schema.my_spatial_table (id serial);
SELECT AddGeometryColumn ('my_schema','my_spatial_table','geom',4326,'POINT',2);
SELECT AddGeometryColumn ('my_schema','my_spatial_table','geom_c',4326,'POINT',2, false);
SELECT AddGeometryColumn('my_schema','my_spatial_table','geomcp_c',4326,'CURVEPOLYGON',2,false);

2、DropGeometryColumn

从空间表中删除某一几何类型字段。（与AddGeometryColumn类似）
函数：

text DropGeometryColumn(varchar table_name, varchar column_name);
text DropGeometryColumn(varchar schema_name, varchar table_name, varchar column_name);
text DropGeometryColumn(varchar catalog_name, varchar schema_name, varchar table_name, varchar column_name);

例如：

SELECT DropGeometryColumn ('my_schema','my_spatial_table','geom');
--2.0版本里也可用该方法删除字段。
ALTER TABLE my_schema.my_spatial_table DROP column geom;

3、DropGeometryTable

删除一个表和所有的空间字段。

函数：

boolean DropGeometryTable(varchar table_name);
boolean DropGeometryTable(varchar schema_name, varchar table_name);
boolean DropGeometryTable(varchar catalog_name, varchar schema_name, varchar table_name);

例如：

SELECT DropGeometryTable ('my_schema','my_spatial_table');
DROP TABLE my_schema.my_spatial_table;

4、获取PostGIS版本相关信息函数（以下函数均返回text）

PostGIS_Full_Version 返回所有版本信息
PostGIS_GEOS_Version 返回GEOS类库的版本信息，及是否可用。
PostGIS_LibXML_Version 返回libxml2库的版本信息
PostGIS_Lib_Build_Date 返回PostGIS类库构建时间
PostGIS_Lib_Version 返回PostGIS类库版本
PostGIS_PROJ_Version 返回PROJ4类库的版本信息，及是否可用。
PostGIS_Scripts_Build_Date 返回Scripts构建时间
PostGIS_Scripts_Installed 返回安装脚本版本
PostGIS_Scripts_Released 返回卸载脚本版本
PostGIS_Version 返回PostGIS版本信息

例如：

SELECT PostGIS_Full_Version();

SELECT PostGIS_GEOS_Version();

SELECT PostGIS_LibXML_Version();

SELECT PostGIS_Lib_Build_Date();

SELECT PostGIS_Lib_Version();

SELECT PostGIS_PROJ_Version();

SELECT PostGIS_Scripts_Build_Date();

SELECT PostGIS_Scripts_Installed();

SELECT PostGIS_Scripts_Released();

SELECT PostGIS_Version();定义的

5、Populate_Geometry_Columns

使用时，须确保定义的字段可修改或者空间参考已经注册。use_typmod=false时，获得过去行为。

函数：

text Populate_Geometry_Columns(boolean use_typmod=true);
int Populate_Geometry_Columns(oid relation_oid, boolean use_typmod=true);

注意：
2.0.0版本改变：现在使用类型修饰符替换检查约束约束几何类型。你仍可使用检查修饰替换新use_typmod，设置为 false。

2.0.0版本加强：介绍了use_typmod可选参数允许被控制，如果列被创建typmodifiers或检查约束。

例如：

CREATE TABLE public.myspatial_table_cs(gid serial, geom geometry);
INSERT INTO myspatial_table_cs(geom) VALUES(ST_GeomFromText('LINESTRING(1 2, 3 4)',4326) );
SELECT Populate_Geometry_Columns('public.myspatial_table_cs'::regclass, false);

6、UpdateGeometrySRID

修改空间参考SRID，也可用ST_SetSRID，注意修改的SRID必须是模式中存在的。

函数：

text UpdateGeometrySRID(varchar table_name, varchar column_name, integer srid);
text UpdateGeometrySRID(varchar schema_name, varchar table_name, varchar column_name, integer srid);
text UpdateGeometrySRID(varchar catalog_name, varchar schema_name, varchar table_name, varchar column_name, integer srid);

postgis操作geometry知识二

1、WKT定义几何对象格式：

POINT(0 0) ——点
LINESTRING(0 0,1 1,1 2) ——线
POLYGON((0 0,4 0,4 4,0 4,0 0),(1 1, 2 1, 2 2, 1 2,1 1)) ——面
MULTIPOINT(0 0,1 2) ——多点
MULTILINESTRING((0 0,1 1,1 2),(2 3,3 2,5 4)) ——多线
MULTIPOLYGON(((0 0,4 0,4 4,0 4,0 0),(1 1,2 1,2 2,1 2,1 1)), ((-1 -1,-1 -2,-2 -2,-2 -1,-1 -1))) ——多面
GEOMETRYCOLLECTION(POINT(2 3),LINESTRING((2 3,3 4))) ——几何集合

2、常用函数：

wkt转geometry st_geomfromtext(wkt,wkid)

geometry转wkt st_astext(geom)

获取点对象x、y坐标值 st_x(geom)、st_y(geom)

获取线/面对象四至 st_xmin(geom)、st_ymin(geom)、st_xmax(geom)、st_ymax(geom)

计算两点之间距离 st_distance(geom,geom) / st_distance(wkt,wkt)

计算线的长度 st_length(geom) / st_length(wkt)

计算面积 st_area(geom) / st_area(wkt)

缓冲区计算 st_buffer(geom,distance) / st_buffer(wkt,distance)

3、管理函数：

添加几何字段 AddGeometryColumn(, , , , , )
删除几何字段 DropGeometryColumn(, , )
检查数据库几何字段并在geometry_columns中归档 Probe_Geometry_Columns()
给几何对象设置空间参考（在通过一个范围做空间查询时常用） ST_SetSRID(geometry, integer)

4、几何对象关系函数：

获取两个几何对象间的距离 ST_Distance(geometry, geometry)
如果两个几何对象间距离在给定值范围内，则返回TRUE ST_DWithin(geometry, geometry, float)
判断两个几何对象是否相等
（比如LINESTRING(0 0, 2 2)和LINESTRING(0 0, 1 1, 2 2)是相同的几何对象） ST_Equals(geometry, geometry)
判断两个几何对象是否分离 ST_Disjoint(geometry, geometry)
判断两个几何对象是否相交 ST_Intersects(geometry, geometry)
判断两个几何对象的边缘是否接触 ST_Touches(geometry, geometry)
判断两个几何对象是否互相穿过 ST_Crosses(geometry, geometry)
判断A是否被B包含 ST_Within(geometry A, geometry B)
判断两个几何对象是否是重叠 ST_Overlaps(geometry, geometry)
判断A是否包含B ST_Contains(geometry A, geometry B)
判断A是否覆盖 B ST_Covers(geometry A, geometry B)
判断A是否被B所覆盖 ST_CoveredBy(geometry A, geometry B)
通过DE-9IM 矩阵判断两个几何对象的关系是否成立 ST_Relate(geometry, geometry, intersectionPatternMatrix)
获得两个几何对象的关系（DE-9IM矩阵） ST_Relate(geometry, geometry)

5、几何对象处理函数：

获取几何对象的中心 ST_Centroid(geometry)
面积量测 ST_Area(geometry)
长度量测 ST_Length(geometry)
返回曲面上的一个点 ST_PointOnSurface(geometry)
获取边界 ST_Boundary(geometry)
获取缓冲后的几何对象 ST_Buffer(geometry, double, [integer])
获取多几何对象的外接对象 ST_ConvexHull(geometry)
获取两个几何对象相交的部分 ST_Intersection(geometry, geometry)
将经度小于0的值加360使所有经度值在0-360间 ST_Shift_Longitude(geometry)
获取两个几何对象不相交的部分（A、B可互换） ST_SymDifference(geometry A, geometry B)
从A去除和B相交的部分后返回 ST_Difference(geometry A, geometry B)
返回两个几何对象的合并结果 ST_Union(geometry, geometry)
返回一系列几何对象的合并结果 ST_Union(geometry set)
用较少的内存和较长的时间完成合并操作，结果和ST_Union相同 ST_MemUnion(geometry set)

6、几何对象存取函数：

获取几何对象的WKT描述 ST_AsText(geometry)
获取几何对象的WKB描述 ST_AsBinary(geometry)
获取几何对象的空间参考ID ST_SRID(geometry)
获取几何对象的维数 ST_Dimension(geometry)
获取几何对象的边界范围 ST_Envelope(geometry)
判断几何对象是否为空 ST_IsEmpty(geometry)
判断几何对象是否不包含特殊点（比如自相交） ST_IsSimple(geometry)
判断几何对象是否闭合 ST_IsClosed(geometry)
判断曲线是否闭合并且不包含特殊点 ST_IsRing(geometry)
获取多几何对象中的对象个数 ST_NumGeometries(geometry)
获取多几何对象中第N个对象 ST_GeometryN(geometry,int)
获取几何对象中的点个数 ST_NumPoints(geometry)
获取几何对象的第N个点 ST_PointN(geometry,integer)
获取多边形的外边缘 ST_ExteriorRing(geometry)
获取多边形内边界个数 ST_NumInteriorRings(geometry)
同上 ST_NumInteriorRing(geometry)
获取多边形的第N个内边界 ST_InteriorRingN(geometry,integer)
获取线的终点 ST_EndPoint(geometry)
获取线的起始点 ST_StartPoint(geometry)
获取几何对象的类型 GeometryType(geometry)
类似上，但是不检查M值，即POINTM对象会被判断为point ST_GeometryType(geometry)
获取点的X坐标 ST_X(geometry)
获取点的Y坐标 ST_Y(geometry)
获取点的Z坐标 ST_Z(geometry)
获取点的M值 ST_M(geometry)

7、几何对象构造函数：

参考语义：
Text：WKT
WKB：WKB
Geom:Geometry
M:Multi
Bd:BuildArea
Coll:Collection ST_GeomFromText(text,[])

ST_PointFromText(text,[])
ST_LineFromText(text,[])
ST_LinestringFromText(text,[])
ST_PolyFromText(text,[])
ST_PolygonFromText(text,[])
ST_MPointFromText(text,[])
ST_MLineFromText(text,[])
ST_MPolyFromText(text,[])
ST_GeomCollFromText(text,[])
ST_GeomFromWKB(bytea,[])
ST_GeometryFromWKB(bytea,[])
ST_PointFromWKB(bytea,[])
ST_LineFromWKB(bytea,[])
ST_LinestringFromWKB(bytea,[])
ST_PolyFromWKB(bytea,[])
ST_PolygonFromWKB(bytea,[])
ST_MPointFromWKB(bytea,[])
ST_MLineFromWKB(bytea,[])
ST_MPolyFromWKB(bytea,[])
ST_GeomCollFromWKB(bytea,[])
ST_BdPolyFromText(text WKT, integer SRID)
ST_BdMPolyFromText(text WKT, integer SRID)