PostgreSQL , PostGIS , 球坐标 , 平面坐标 , 球面距离 , 平面距离

最新推荐文章于 2024-04-21 15:53:49 发布
最新推荐文章于 2024-04-21 15:53:49 发布
阅读量1.3k 收藏 3
点赞数 1
文章标签： postgresql 数据库
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/superhoy/article/details/127030092
版权

标签

PostgreSQL , PostGIS , 球坐标 , 平面坐标 , 球面距离 , 平面距离

背景

PostGIS中有两种常用的空间类型geometry和geography，这两种数据类型有什么差异，应该如何选择？

对于GIS来说，首先是坐标系，有两种：一种是球坐标（地理坐标），另一种是平面坐标（投影坐标）。

球坐标通常用于计算，平面坐标通常用于展示（也可以计算）。

投影坐标是从球坐标投影后展开得来(用一个圆柱将地球包起来，把地球当成会发光的光源，投影后，将圆柱展开得到)，投影的范围越大，精度就越低。范围越小，精度就越高。除了投影扇形的大小有区别，在不同的行业，也有不同的坐标系，例如用于测绘的坐标系。

目前用得最多的有SRID=4326球坐标，SRID为EPSG:3785的墨卡托投影坐标。

再来说geometry和geography两种类型，geometry支持平面对象也支持空间对象，而geography则仅支持空间对象。

geometry支持更多的函数，一些几何计算的代价更低。

geography支持的函数略少，计算代价更高。那为什么还要geography呢？因

4.2.2. When to use Geography Data type over Geometry data type    
    
The geography type allows you to store data in longitude/latitude coordinates,     
but at a cost: there are fewer functions defined on GEOGRAPHY than there are on GEOMETRY;     
those functions that are defined take more CPU time to execute.    
    
The type you choose should be conditioned on the expected working area of the application you are building.     
Will your data span the globe or a large continental area, or is it local to a state, county or municipality?    
    
If your data is contained in a small area, you might find that choosing an appropriate     
projection and using GEOMETRY is the best solution, in terms of performance and functionality available.    
    
If your data is global or covers a continental region, you may find that GEOGRAPHY     
allows you to build a system without having to worry about projection details.     
You store your data in longitude/latitude, and use the functions that have been defined on GEOGRAPHY.    
    
If you don't understand projections, and you don't want to learn about them,     
and you're prepared to accept the limitations in functionality available in GEOGRAPHY,     
then it might be easier for you to use GEOGRAPHY than GEOMETRY.     
Simply load your data up as longitude/latitude and go from there.      
    
Refer to Section 14.11, “PostGIS Function Support Matrix” for compare between     
what is supported for Geography vs. Geometry.     
For a brief listing and description of Geography functions,     
refer to Section 14.4, “PostGIS Geography Support Functions”

既然提到距离计算和投影坐标系有关，引入了本文的问题：

在不知道要计算的geometry点，在什么投影坐标系下时，往往计算距离得到的结果并不准确。

例如，下面的点，换算成2163坐标系，计算距离得到的结果并不准确。

db1=# SELECT st_distance(ST_Transform(ST_GeomFromText('POINT(120.08 30.96)', 4326), 2163 ), ST_Transform(ST_GeomFromText('POINT(120.08 30.92)', 4326), 2163 ));   
   st_distance          
------------------      
 4030.46766234184      
(1 row)

2163坐标系内容如下：

postgres=# select * from spatial_ref_sys where srid=2163;    
-[ RECORD 1 ]-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------    
srid      | 2163    
auth_name | EPSG    
auth_srid | 2163    
srtext    | PROJCS["US National Atlas Equal Area",GEOGCS["Unspecified datum based upon the Clarke 1866 Authalic Sphere",DATUM["Not_specified_based_on_Clarke_1866_Authalic_Sphere",SPHEROID["Clarke 1866 Authalic Sphere",6370997,0,AUTHORITY["EPSG","7052"]],AUTHORITY["EPSG","6052"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4052"]],PROJECTION["Lambert_Azimuthal_Equal_Area"],PARAMETER["latitude_of_center",45],PARAMETER["longitude_of_center",-100],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["X",EAST],AXIS["Y",NORTH],AUTHORITY["EPSG","2163"]]    
proj4text | +proj=laea +lat_0=45 +lon_0=-100 +x_0=0 +y_0=0 +a=6370997 +b=6370997 +units=m +no_defs     

AUTHORITY["EPSG", "9122"]指的是EPSG数据集中UNIT为degree的ID是9122；  
  
AUTHORITY["EPSG", "4326"]指的是地理坐标系WGS 84的ID是4326；  
  
AUTHORITY["EPSG", "9001"]指的是EPSG中UNIT为meter的ID是9001；  
  
AUTHORITY["EPSG", "32650"]指的是该投影坐标系WGS 84 / UTM zone 50N的ID是32650

这样会造成一个假象如下：

用st_distance函数计算出来的经纬度之间的距离，跟用java程序算出来的距离相差很大。      
      
这个例子，st_distance算出来的距离是4030，我们程序算出来的是4445，换另外两个相距很远的点，这个差距值会更大。

正确姿势

算球面距离，不要算直线距离。

1、使用球坐标，通过st_distancespheroid计算两个geometry类型的球面距离。

对于geometry类型，可以使用st_distancespheroid直接用球坐标计算，在计算时会自动设置这个椭球特性（SPHEROID["Krasovsky_1940",6378245.000000,298.299997264589] ）。

postgres=# SELECT st_distancespheroid(ST_GeomFromText('POINT(120.08 30.96)', 4326),ST_GeomFromText('POINT(120.08 30.92)', 4326), 'SPHEROID["WGS84",6378137,298.257223563]');    
 st_distancespheroid     
---------------------    
    4434.73734584354    
(1 row)

2、使用平面坐标，通过st_distance计算两个geometry类型的平面投影后的距离。

采用精准的投影坐标（小面积投影坐标系）（但是必须要覆盖到要计算的两个点）

ST_Transform(ST_GeomFromText('POINT(120.08 30.96)', 4326), 2163 )

如果geometry值的SRID不是（高精度）目标坐标系，可以使用ST_Transform函数进行转换，转换为目标投影坐标系，再计算距离。

postgres=# SELECT st_distance(ST_Transform(ST_GeomFromText('POINT(120.08 30.96)', 4326), 2390 ), ST_Transform(ST_GeomFromText('POINT(120.08 30.92)', 4326), 2390 ));   
   st_distance      
------------------  
 4547.55477647394  
(1 row)

3、使用球坐标，通过ST_Distance计算两个geography类型的球面距离。

float ST_Distance(geography gg1, geography gg2, boolean use_spheroid);   -- 适用椭球体（WGS84）    
    
    
use_spheroid设置为ture表示使用:      
  -- WGS84 椭球体参数定义    
  vspheroid := 'SPHEROID["WGS84",6378137,298.257223563]' ;

这里的XXXX就是你要选择的球坐标系SRID。在spatial_ref_sys表里可以查看各种坐标系。

postgres=# SELECT st_distance(ST_GeogFromText('SRID=xxxx;POINT(120.08 30.96)'), ST_GeogFromText('SRID=xxxx;POINT(120.08 30.92)'), true);    
  st_distance       
----------------    
 xxxxxxxxxxxxxxxxxxxx    
(1 row)    
  
例如  
  
postgres=# SELECT st_distance(ST_GeogFromText('SRID=4610;POINT(120.08 30.96)'), ST_GeogFromText('SRID=4610;POINT(120.08 30.92)'), true);    
  st_distance     
----------------  
 4434.739418211  
(1 row)

如果允许一定的偏差，可以使用全球球坐标系4326。

postgres=# SELECT st_distance(ST_GeogFromText('SRID=4326;POINT(120.08 30.96)'), ST_GeogFromText('SRID=4326;POINT(120.08 30.92)'), true);    
  st_distance       
----------------    
 4434.737345844    
(1 row)

geography只支持球坐标系，使用投影坐标会报错。

postgres=# SELECT st_distance(ST_GeogFromText('SRID=2369;POINT(120.08 30.96)'), ST_GeogFromText('SRID=2369;POINT(120.08 30.92)'), true);    
错误:  22023: Only lon/lat coordinate systems are supported in geography.  
LOCATION:  srid_is_latlong, lwgeom_transform.c:774

4、指定SPHEROID内容，通过st_distancesphereoid计算geometry的球面距离。

这种方法最为精确，但是要求了解计算距离当地的地形属性（即输入参数的spheroid的内容）。

db1=# SELECT st_distancespheroid(ST_GeomFromText('POINT(120.08 30.96)', 4326),ST_GeomFromText('POINT(120.08 30.92)', 4326), 'SPHEROID["WGS84",6378137,298.257223563]');        
 st_distancesphere       
-------------------      
    4447.803189385      
(1 row)

小结

计算距离，应该考虑到被计算的两点所在处的地球特性（spheroid）。这样计算得到的距离才是最精确的。

geometry和geography类型的选择，建议使用geometry，既能支持球坐标系，又能支持平面坐标系。主要考虑到用户是否了解位置所在处的地理特性，选择合适的坐标系。

-- 适用平面直角坐标，适用geometry类型，计算直线距离    
float ST_Distance(geometry g1, geometry g2);             
    
-- 适用椭球体坐标，适用geography类型，计算球面距离  
float ST_Distance(geography gg1, geography gg2);         
    
-- 适用椭球体坐标（WGS84），适用geography类型，计算球面距离    
float ST_Distance(geography gg1, geography gg2, boolean use_spheroid);     
    
   use_spheroid设置为ture表示使用:      
     vspheroid := 'SPHEROID["WGS84",6378137,298.257223563]' ; -- WGS84 椭球体参数定义     
    
-- 适用椭球体坐标以及投影坐标，适用geometry类型，求球面距离（需要输入spheroid）  
float ST_DistanceSpheroid(geometry geomlonlatA, geometry geomlonlatB, spheroid measurement_spheroid);

参考

1、计算球面距离

ST_DistanceSphere

2、计算球面以及平面距离(取决于输入的类型geometry还是geography)

ST_Distance

3、坐标系转换

ST_Transform

4、投影坐标和球坐标

《地理坐标系（球面坐标系）和投影坐标系（平面坐标系）》

5、PostGIS学习建议

《PostGIS 空间数据学习建议》

6、 http://blog.163.com/jey_df/blog/static/18255016120149145755781/

墨卡托 (Mercator) 投影—帮助 | ArcGIS for Desktop

键盘tops舞者
关注
  • 1
    点赞
  • 踩
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
地理空间参照系
Vonzhuxi的博客
07-10 1999
空间参照系 空间参照系（Space Reference System，SRS）是用来描述几何的定义空间的参照系。空间参照系至少定义一下内容： 基础坐标的测量单位(度，米等) 最大坐标值和最小坐标值（也称为边界） 缺省线性测量单位 数据类型，平面数据还是椭球体数据 用于将数据转换为其他SRS的投影信息 空间标识符 每个空间参照系都有一个标识符，称为空间参照系标识符（Spatial Reference Identifier，SRID）。下面是常用的几种SRID。 SRID=0 缺省空间参照系的数据默认使
postgresql14+postgis3.3+pgrouting3.4的docker镜像
08-17
3-3" --restart always -e TZ=Asia/Shanghai -p 5432:5432 -e POSTGRES_USER=postgres -e POSTGRES_PASSWORD=postgres -e PASSWORD_AUTHENTICATION="scram-sha-256" -v 自定义挂载目录:/var/lib/postgresql/data -d ...
12.（地图工具篇）PostgreSQL+PostGIS实现两坐标点之间最短路径查询算法函数.zip
09-26
下载如有问题，可私信博主。下载前建议先查看博客内容，其地址为：https://blog.csdn.net/QQ98281642/article/details/120483897
postgresql14加postgis3.2在centos离线编译安装
04-28
postgresql14加postgis3.2在centos离线编译安装
PostgreSql+PostGis创建空间数据库
05-03
PostgreSql+PostGis 创建空间数据库
postgresql自定义方法创建之坐标的自动转换
梁海江的博客
08-18 707
postgresql导入点类型的空间数据，一般会提供xy坐标(或者lon和lat)字段，一般的工具并没有提供转换为空间数据的功能。此时，我们需要提供一些方法或者触发器来实现空间字段的自动更新。而xy坐标在导入的时候也会存在两种情况，一种是以度分秒的字符串形式提供，一种是小数的方式提供。所以这里我们还要提供一个针对坐标的归一化函数。下面就是相关设计的方法方法代码。
PostGIS 距离计算规范 - 投影 与 球 坐标系, geometry 与 geography 类型
yuhui666666的博客
07-02 1844
平面坐标，球面坐标 SRID=2369 投影坐标 SRID=4326 球面坐标 ST_Transform(ST_GeomFromText('POINT(120.08 30.96)', 4326), 2390 ) 坐标系转换函数 PostGIS中有两种常用的空间类型geometry和geography，这两种数据类型有什么差异，应该如何选择？ 对于GIS来说，首先是坐标系，有两种：一种是球坐标（地理坐标），另一种是平面坐标（投影坐标）。 球坐标通常用于计算，平面坐标通常用于展示（也可以计算）...
PostGIS中geometry与geography的区别
SuperMap技术控
03-29 3481
除了geometry，PostGIS还定义了geography，用于不同的场景需求。以下通过实战对比介绍geometry与geography之间的差异、适用场景及互转方法。
PostGIS Geography Type
XiaoCai
09-26 642
PostGIS的地理几何数据类型
在pg数据库中获取俩点距离
qq_43757153的博客
09-04 411
PG数据是库是一个相对于MySQL来说更安全的开源数据库，但是在市面上对PG的使用还是比较少的，博主就结合最近的工作问题，总结一下在PG数据库中如何获取到两个坐标点之间的距离（当然，也有获取点到线，点到点的距离，下面只举例两个坐标点的距离获取方式）
超图球面与平面场景的区别和地理坐标系
bcbobo21cn的专栏
07-04 2336
前文在超图桌面版中加载了一个3dmax模型；是加载到球面场景中的；看一下能不能加载到平面场景；加载到平面场景，出现下图错误；提示数据集的坐标系是地理坐标系，不能添加到平面场景， 为什么地理坐标系不能加载到平面场景，先看一下超图帮助， 球面场景是指以模拟地球的球体三维空间形式对地球表层的场景进行展示。平面场景则使地球球面展开成平面，模拟整个大地类似一个平面的形式进行场景展示。平面场景最好在场景比较小的时候使用， 两种三维视图模式所支持的功能和操作方式大部分都相同，二者之间的差别有以下几点： 球面场
postgresql12-docker19-postgis
07-28
压缩包中包含postgresql12-docker19-postgis三个组件，每个组件都是当下比较常用的版本，并且其中还包含了现在比较流行的postgis，目前只有一个docker的安装文档，其他两个我这边还没总结，有需要的先自行百度一下...
Postgresql Geometry与wkt互转及长度计算
qq_40985985的博客
07-29 2330
这篇博客将介绍如何使用Postgresql 进行Geometry与wkt互转及长度计算。首先构建表和数据，然后sql计算。
PostGis
supernever专栏
10-22 1688
http://baike.baidu.com/link?url=VSOEQ4DFNRm9kzC5XSToxnJVKKKfrCkL5Mo-6bR_HsZ1wPAoXoQ1fVctBO-lRwvzf8Ys9ZQ1YjwzlsdFRqz9ZK postgis PostgreSQL 是一种对象-关系型数据库管理系统（ORDBMS），也是目前功能最强大、特性最丰富和最复杂的自由软件数
PostGIS空间数据库SRID背景知识 - 地理坐标系（球面坐标系）和投影坐标系（平面坐标系）...
weixin_34123613的博客
09-11 1600
背景 《PostGIS 坐标转换(SRID)的边界问题引发的背景知识 - ST_Transform》一文，介绍了使用某个坐标系时，计算国内某个点到某个点的距离出现了负数。 背景知识和坐标系有关。 什么是地理坐标系，什么是投影坐标系？ 参考此文： http://www.cnblogs.com/jetz/archive/2005/03/29/127547....
postgis基本操作
hwq317622817的博客
01-04 411
postgis
点要素生成面要素（Arcgis实操系列）
CCLasdfg的博客
10-11 5751
点转面，点要素转为面要素，Arcgis实操
【树莓派Linux内核开发】入门实操篇（虚拟机Ubuntu环境搭建+内核源码获取与配置+内核交叉编译+内核镜像挂载）
最新发布
ArsenLupin的博客
04-21 1137
本文详写记录了树莓派Linux内核开发基础入门笔记，其中包括Ubuntu虚拟机搭建，内核源码获取与配置，内核交叉编译以及内核镜像挂载。
postgresql postgis
03-16
PostgreSQL是一种开源的关系型数据库管理系统，而PostGIS是一个开源的空间数据库扩展，它允许在PostgreSQL中存储和查询地理空间数据。PostGIS提供了许多空间函数和操作符，可以用于处理和分析地理空间数据。它支持多种地理空间数据格式，包括点、线、面和多边形等。PostgreSQL和PostGIS的结合可以为地理空间数据的存储、查询和分析提供强大的支持。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值