PostgresqlAndPostGis

Postgresql&PostGis

问题及解决方法

关于postGIS没有template_postgis模版的问题解决

// 1、建立普通数据库
  create database xx;
  
// 2、然后输入官网给的这几条添加扩展语句(官网可找到)：
-- Enable PostGIS (includes raster)
CREATE EXTENSION postgis;
-- Enable Topology
CREATE EXTENSION postgis_topology;
-- Enable PostGIS Advanced 3D 
-- and other geoprocessing algorithms
-- sfcgal not available with all distributions
CREATE EXTENSION postgis_sfcgal;
-- fuzzy matching needed for Tiger
CREATE EXTENSION fuzzystrmatch;
-- rule based standardizer
CREATE EXTENSION address_standardizer;
-- example rule data set
CREATE EXTENSION address_standardizer_data_us;
-- Enable US Tiger Geocoder
CREATE EXTENSION postgis_tiger_geocoder;

关于PostGresql中没有postgis函数问题

create extension postgis	// postgis 相关函数

相关函数的使用

二进制geometry转geojson

select st_asgeojson(the_geom) from cities

geojson转二进制geometry

select st_geomfromgeojson('{"type":"Point","coordinates":[27.91162491,-33.01529]}')

将二进制的geometry格式的转换成 4326(经纬度)的格式

update test3 set geom = st_setsrid(geom,4326)

将wkt转成geometry

SELECT   ST_GeomFromText(wkt)  from cities

创建空间索引

create index idx_test_shenz on test_shenz(geom) using gist

将geometry转换成wkt格式：

select st_astext(geometry) from cities limit 3;

将数据库中的经纬度转为geometry格式

// ALTER TABLE ajxx ADD geom geometry;   
update ajxx set geom=st_geomfromtext('Point('||dqjd||' '||dqwd||')',4326) // 经度/纬度

插入数据返回Id

insert into point(pointtype,pointx,pointy,pointval)values(1,2,3,4) RETURNING id;

将查询的结果导成新表

CREATE TABLE films_recent AS
 SELECT * FROM films WHERE date_prod >= '2002-01-01';

将shp文件转成sql 直接入pg库

D:\Program Files\PostgreSQL9.6\bin>shp2pgsql -s 4326 -W "GBK" D:\FZ_PCS_PG_ALL.shp fz_pcs >  D:\fz_pcs.sql

异常

ERROR:  relation "public.kafka_swan_vehicle_id_seq" does not exist

解决

新建查询，执行 CREATE SEQUENCE IF NOT EXISTS XXXX_id_seq;

表中含有自增主键导出sql文件 再导入另外一个数据库时自增主键报错问题解决

//有自增字段的表导出成sql文件之后，执行的时候，要在建表之前加上 digitalhub_country_gid_seq  代表主键
DROP SEQUENCE if EXISTS "public"."digitalhub_country_gid_seq"; 
CREATE SEQUENCE "public"."digitalhub_country_gid_seq"  
 INCREMENT 1  
MINVALUE 1  
MAXVALUE 9223372036854775807  
START 1  
CACHE 1;  
就是把之前的删掉了再建

判断两点之间函数

select ST_Length(Geography(ST_GeomFromText('LINESTRING(116.47 39.89,116.47 39.95)')));

查看sql正在运行的进程、查看表是否被锁、解锁表

SELECT
    procpid,
    START,
    now() - START AS lap,
    current_query
FROM
    (
        SELECT
            backendid,
            pg_stat_get_backend_pid (S.backendid) AS procpid,
            pg_stat_get_backend_activity_start (S.backendid) AS START,
            pg_stat_get_backend_activity (S.backendid) AS current_query
        FROM
            (
                SELECT
                    pg_stat_get_backend_idset () AS backendid
            ) AS S
    ) AS S
WHERE
    current_query <> '<IDLE>'
ORDER BY
    lap DESC;

• procpid：进程id
• start：进程开始时间
• lap：经过时间
• current_query：执行中的sql
• 怎样停止正在执行的sql
• SELECT pg_cancel_backend(进程id);
• 或者用系统函数
• kill -9 进程id;

查看表是否被锁，解锁表

select oid from pg_class where relname='m_ss_kjcx_tech_inovate_talent' 
select pid from pg_locks where relation='4384913' 
--如果查询到了结果，表示该表被锁 则需要释放锁定 
select pg_cancel_backend(140650678843136)
通过命令： 
=# select pg_cancel_backend(线程id); 
来kill掉指定的SQL语句。(这个函数只能 kill Select 查询，而updae,delete DML不生效)
使用可以kill 各种DML(SELECT,UPDATE,DELETE,DROP)操作= 
=# select pg_terminate_backend(pid int)
-- PGSQL
select oid from pg_class where relname='m_ss_kjcx_tech_inovate_talent' 
select pid from pg_locks where relation='4384913' 
--如果查询到了结果，表示该表被锁 则需要释放锁定 
select pg_cancel_backend(140650678843136)
--- 来查看有哪些SQL正在执行。 waiting='t'
select pid,backend_start,application_name,query_start,waiting,state ,query from pg_stat_activity where query <>'' order by query_start asc
通过命令： 
=# select pg_cancel_backend(线程id); 
来kill掉指定的SQL语句。(这个函数只能 kill Select 查询，而updae,delete DML不生效)
使用可以kill 各种DML(SELECT,UPDATE,DELETE,DROP)操作= 
=# select pg_terminate_backend(pid int)

常用函数

PostGIS中的常用函数
以下内容包括比较多的尖括号，发布到blogger的时候会显示不正常，内容太多我也无暇一个个手动改代码，因此如有问题就去参考PostGIS官方文档。
首先需要说明一下，这里许多函数是以ST_[X]yyy形式命名的，事实上很多函数也可以通过xyyy的形式访问，在PostGIS的函数库中我们可以看到这两种函数定义完全一样。
1. OGC标准函数
管理函数：
添加几何字段 AddGeometryColumn(, , , , , )
删除几何字段 DropGeometryColumn(, , )
检查数据库几何字段并在geometry_columns中归档 Probe_Geometry_Columns()
给几何对象设置空间参考(在通过一个范围做空间查询时常用)ST_SetSRID(geometry, integer)
几何对象关系函数：
获取两个几何对象间的距离                                                                                        ST_Distance(geometry, geometry)
如果两个几何对象间距离在给定值范围内，则返回                                                    TRUEST_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)
几何对象处理函数：
获取几何对象的中心                                                                                                      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, geometryB)
返回两个几何对象的合并结果                                                                                       ST_Union(geometry, geometry)
返回一系列几何对象的合并结果                                                                                    ST_Union(geometry set)
用较少的内存和较长的时间完成合并操作，结果和ST_Union相同 ST_MemUnion(geometry set)
几何对象存取函数：
获取几何对象的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)
几何对象构造函数：
参考语义：
Text：WKT
WKB：WKB
Geom:Geometry
M:Multi
Bd:BuildArea
Coll:Collection ST_GeomFromText(text,[])    将wkt转为geometry
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)