首页 > 数据库 > 其他综合 >

PostGIS根据已有列创建Geometry,并计算距离(教程)

2018-02-01

PostGIS根据已有列创建Geometry,并计算距离(教程)。1 创建包含坐标系信息的Geometry列;2 从已有字段生成POINT列;3 更新已有的Geometry字段并设置srid;4 计算距离

1.创建包含坐标系信息的Geometry列

AddGeometryColumn的用法如下:

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

通过如下语句生成一个坐标系为srid=4326的POINT类型的Geometry列:

SELECT AddGeometryColumn ('public','road_srid','geom_srid',4326,'POINT',2);

2.从已有字段生成POINT列

通过ST_MakePoint函数生成POINT,x、y为表中原有字段,通过ST_SetSRID函数设置srid,使其匹配第一步中的srid

update public.road_srid set xy_srid=ST_SetSRID(ST_MakePoint(x,y),4326)

3.更新已有的Geometry字段并设置srid

geom为未设置srid的Geometry字段,将geom转化为WKB,使用ST_GeomFromWKB函数生成Geometry:

UPDATE public.road_srid SET geom = ST_GeomFromWKB(ST_AsBinary(geom),4326);

4.计算距离

使用ST_Distance计算两个Geometry的距离;如果要计算的两个Geometry的srid不同则无法进行计算,因此需要保证计算的两个Geometry的srid相同;

SELECT ST_Distance(
        ST_GeomFromText('POINT(114.1235 34.3521)',4326),
        ST_GeomFromText('LINESTRING(114.1260 34.45, 114.123 34.1546)', 4326)
    );

以上计算出的结果为0.00150567726382282,单位为度,由于此srid代表wgs84大地坐标系,单位为度。

SELECT ST_Distance(
            ST_Transform(ST_GeomFromText('POINT(-72.1235 42.3521)',4326),900913),
            ST_Transform(ST_GeomFromText('LINESTRING(-72.1260 42.45, -72.123 42.1546)', 4326),900913)
        );

以上计算出的结果为123.797937878428,通过ST_Transform将大地坐标系进行投影,得到的是投影坐标系下的坐标,此时计算出的结果单位为米;其中srid 为900913是墨卡托投影。

相关文章
最新文章
热点推荐