名称

ST_Transform — 返回坐标已转换到不同空间参考系统的新几何。

语法

geometry ST_Transform(geometry g1, integer srid);

geometry ST_Transform(geometry geom, text to_proj);

geometry ST_Transform(geometry geom, text from_proj, text to_proj);

geometry ST_Transform(geometry geom, text from_proj, integer to_srid);

说明

返回坐标已转换到不同空间参考系统的新几何。目标空间参考 to_srid 可通过有效的 SRID 整数参数标识(即它必须存在于 spatial_ref_sys 表中)。或者,可以将定义为 PROJ.4 字符串的空间参考用于 to_proj 和/或 from_proj,但这些方法未经过优化。如果目标空间参考系统用 PROJ.4 字符串(而非 SRID)表示,则输出几何的 SRID 将设置为零。除了具有 from_proj 的函数外,输入几何必须具有已定义的 SRID。

ST_Transform 通常与 ST_SetSRID 混淆。ST_Transform 实际上会将几何的坐标从一个空间参考系统更改到另一个空间参考系统,而 ST_SetSRID() 只是更改几何的 SRID 标识符。

ST_Transform 会根据源空间参考系统和目标空间参考系统自动选择合适的转换管道。要使用特定的转换方法,请使用 ST_TransformPipeline

[Note]

需要使用 PROJ 支持编译 PostGIS。使用 PostGIS_Full_Version 确认已编译 PROJ 支持。

[Note]

如果使用多个转换,则在常用的转换上具有函数索引以利用索引使用情况非常有用。

[Note]

在 1.3.4 之前,如果与包含 CURVES 的几何一起使用,此函数会崩溃。此问题在 1.3.4+ 中已修复。

增强:2.0.0 引入了对多面体的支持。

增强:2.3.0 引入了对直接 PROJ.4 文本的支持。

此方法实现了 SQL/MM 规范。

SQL-MM 3: 5.1.6

此方法支持圆形字符串和曲线。

此函数支持多面体表面。

示例

将马萨诸塞州平面美国英尺几何图形更改为 WGS 84 长纬度

SELECT ST_AsText(ST_Transform(ST_GeomFromText('POLYGON((743238 2967416,743238 2967450,
  743265 2967450,743265.625 2967416,743238 2967416))',2249),4326)) As wgs_geom;

 wgs_geom
---------------------------
 POLYGON((-71.1776848522251 42.3902896512902,-71.1776843766326 42.3903829478009,
-71.1775844305465 42.3903826677917,-71.1775825927231 42.3902893647987,-71.177684
8522251 42.3902896512902));
(1 row)

--3D Circular String example
SELECT ST_AsEWKT(ST_Transform(ST_GeomFromEWKT('SRID=2249;CIRCULARSTRING(743238 2967416 1,743238 2967450 2,743265 2967450 3,743265.625 2967416 3,743238 2967416 4)'),4326));

         st_asewkt
--------------------------------------------------------------------------------------
 SRID=4326;CIRCULARSTRING(-71.1776848522251 42.3902896512902 1,-71.1776843766326 42.3903829478009 2,
 -71.1775844305465 42.3903826677917 3,
 -71.1775825927231 42.3902893647987 3,-71.1776848522251 42.3902896512902 4)

    

创建部分函数索引的示例。对于不确定是否会填入所有几何图形的表,最好使用部分索引,该索引会忽略空几何图形,这既可以节省空间,还可以使索引更小、更高效。

CREATE INDEX idx_geom_26986_parcels
  ON parcels
  USING gist
  (ST_Transform(geom, 26986))
  WHERE geom IS NOT NULL;
    

使用 PROJ.4 文本进行转换的示例,其中包含自定义空间参考。

-- Find intersection of two polygons near the North pole, using a custom Gnomic projection
-- See http://boundlessgeo.com/2012/02/flattening-the-peel/
 WITH data AS (
   SELECT
     ST_GeomFromText('POLYGON((170 50,170 72,-130 72,-130 50,170 50))', 4326) AS p1,
     ST_GeomFromText('POLYGON((-170 68,-170 90,-141 90,-141 68,-170 68))', 4326) AS p2,
     '+proj=gnom +ellps=WGS84 +lat_0=70 +lon_0=-160 +no_defs'::text AS gnom
 )
 SELECT ST_AsText(
   ST_Transform(
     ST_Intersection(ST_Transform(p1, gnom), ST_Transform(p2, gnom)),
   gnom, 4326))
 FROM data;
                                          st_astext
 --------------------------------------------------------------------------------
  POLYGON((-170 74.053793645338,-141 73.4268621378904,-141 68,-170 68,-170 74.053793645338))
    

配置转换行为

有时涉及网格平移的坐标转换可能会失败,例如,如果 PROJ.4 未使用网格平移文件构建,或者坐标不在定义网格平移的范围内。默认情况下,如果网格平移文件不存在,PostGIS 将引发错误,但可以通过测试不同的 PROJ.4 文本的 to_proj 值或更改 spatial_ref_sys 表中的 proj4text 值,根据每个 SRID 配置此行为。

例如,proj4text 参数 +datum=NAD87 是以下 +nadgrids 参数的简写形式

+nadgrids=@conus,@alaska,@ntv2_0.gsb,@ntv1_can.dat

@ 前缀表示如果文件不存在,则不报告错误,但如果到达列表末尾且没有找到合适的(即找到并重叠)文件,则会发出错误。

相反,如果你想确保至少存在标准文件,但如果扫描所有文件后都没有命中,则应用空转换,你可以使用

+nadgrids=@conus,@alaska,@ntv2_0.gsb,@ntv1_can.dat,null

空网格平移文件是一个有效的网格平移文件,覆盖整个世界且不应用平移。因此,作为一个完整的示例,如果你想更改 PostGIS,以便对不在正确范围内的 SRID 4267 的转换不会引发 ERROR,则可以使用以下内容

UPDATE spatial_ref_sys SET proj4text = '+proj=longlat +ellps=clrk66 +nadgrids=@conus,@alaska,@ntv2_0.gsb,@ntv1_can.dat,null +no_defs' WHERE srid = 4267;

另请参见

第 4.5 节,“空间参考系统”ST_SetSRIDST_SRIDUpdateGeometrySRIDST_TransformPipeline