PostGIS 笔记一：简单入门

Postgis是基于PostgreSQL的空间扩展。所以你懂得PostgreSQL不代表你就懂得Postgis，你懂得Postgis也不代表你就懂得PostgreSQL。如果你看过我关于Spatialite的笔记，你具备了一定空间数据库的基础。

• 1简单的入门
• 2CLI Cheat Sheet
• 3python for Spatialite
• 4NET平台使用spatilite扩展
• 5Spatiliate2GeoJson数据的转换

使用

Postgis分为两个大版本，它们之间有较大幅度的更新，网上很多教程都是基于旧版本的postgis。在安装的阶段当中，新版本的postgis会出现“template_postgis模板不存在”的问题。解决方案有两个：

• 一是参考postgis不具备template_postgis模板自建模板。

• 二是创建PostGIS扩展

CREATE EXTENSION postgis

同样，很多第三方库也会区分postgis的版本，专门针对2.0版本后进行定制。

运用

空间字段可以分为geometry和geography两种，大部分情况使用满足opengis标准的geometry。

创建空间字段

方法一

打开空间扩展后，前提是你首先创建一张表，运用opengis标准推荐的方法创建空间字段。

SELECT AddGeometryColumn ('scheme'，'cities', 'the_geom', 4326, 'POINT', 2);

缺点在于这是一个扩展的方法，表名、字段名等参数是区分大小写的。

方法二

也可以使用geometry type的方法，同时创建表。

CREATE TABLE ROADS ( ID int4, ROAD_NAME varchar(25), geom geometry(LINESTRING,4326) );

或者修改表

ALTER TABLE roads ADD COLUMN geom2 geometry(LINESTRINGZ,4326);

插入一条数据

由于空间数据的描述相对复杂，空间数据的插入需要依靠别的数据格式进行转换。ST_GeomFromText是从文本进行转换。

INSERT INTO roads (road_id, roads_geom, road_name) VALUES (3,ST_GeomFromText('LINESTRING(192783 228138,192612 229814)',-1),'Paul St');

或者用wkt，可以参考WKT&WKB 笔记一：格式介绍。

INSERT INTO table ( SHAPE, NAME )VALUES ( GeomFromEWKT('SRID=4326;POINTM(116.39 39.9 10)'), '北京' )

读取一条数据

显示数据可以将空间字段保存成不同的格式。

 SELECT id, ST_AsText(the_geom), ST_AsGeoJson(the_geom), ST_AsEwkt(the_geom), ST_X(the_geom), ST_Y(the_geom) FROM cities;

用Geography数据格式不用Geometry数据格式

新型的GEOGRAPHY格式允许你用longitude/latitude coordinates存数据，但是有一个代价：针对GEOGRAPHY很少方法被定义，而针对GEOMETRY相对较多。而且这些方法都相对来说更消耗CPU。

你选择这个数据类型必须是因为你应用的情况。你的数据如果是全球或者大陆范围或者是一个洲，一个国家或者一个城市。

• 如果数据在一个很小的范围内，你可能需要选择一个恰当的投影，用Geometry是最好的解决方案，因为性能优先和功能可用。

• 你的数据是全球或者大范围的，你会发现GEOGRAPHY让你建立一个系统不用担心投影信息。你可以用经纬度存起你的数据。

• 你不懂得投影，你也不想去学投影，你可以接受GEOGRAPHY的方法缺失，你应该选择GEOGRAPHY而不是GEOMETRY。因为这样更方便你去加载用经纬度数据。

来自PostGIS_Geography VS Geometry

空间分析方法

PostGIS_GeographyFunctions

参考

PostgreSQL的postgres_fdw跨库使用