PostGIS 笔记一:简单入门
Postgis是基于PostgreSQL的空间扩展。所以你懂得PostgreSQL不代表你就懂得Postgis,你懂得Postgis也不代表你就懂得PostgreSQL。如果你看过我关于Spatialite的笔记,你具备了一定空间数据库的基础。
使用
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。因为这样更方便你去加载用经纬度数据。