postgis 计算

900913 墨卡托投影 guge

select st_area(ST_Transform(ST_GeomFromText(st_astext(data),4326),900913)) from v_shape

以下为参考：

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是墨卡托投影。