説明

geometry型については、二つのジオメトリの、2次元の最小デカルト距離 (空間参照系に基づきます)を、投影法の単位で返します。geography型については、デフォルトでは、二つのジオメトリ間の測地距離をメートル単位で返します。use_spheroidをFALSEにすると、回転楕円体の替わりに球面で計算します。

This method implements the OpenGIS Simple Features Implementation Specification for SQL 1.1.

This method implements the SQL/MM specification. SQL-MM 3: 5.1.23

This method supports Circular Strings and Curves

This method is also provided by SFCGAL backend.

Availability: 1.5.0 1.5でジオグラフィ対応が導入されました。大きいジオメトリや頂点の多いジオメトリについての速度が改善しました。

Enhanced: 2.1.0 ジオグラフィでの速度が改善されました。詳細はMaking Geography fasterをご覧ください。

Enhanced: 2.1.0 - 曲線ジオメトリ対応が導入されました。

Enhanced: 2.2.0 - 精度とロバスト性の向上のためにGeographicLibを使って回転楕円体面上での計測を行うようにしています。この新機能を使うには、Proj 4.9.0以上が必要です。

基本的なジオメトリの例

-- ジオメトリの例 - 度単位 4326はWGS84経度緯度で単位は度
SELECT ST_Distance(
                'SRID=4326;POINT(-72.1235 42.3521)'::geometry,
                'SRID=4326;LINESTRING(-72.1260 42.45, -72.123 42.1546)'::geometry
        );
st_distance
-----------------
0.00150567726382282

-- ジオメトリの例 - メートル単位 (SRID: 3857 一般的なWeb地図のピクセルに比例)
-- 値が外れていても近傍は正しく比較でき、
-- KNNやKMeansのようなアルゴリズムに良い選択となります
SELECT ST_Distance(
                        ST_Transform('SRID=4326;POINT(-72.1235 42.3521)'::geometry, 3857),
                        ST_Transform('SRID=4326;LINESTRING(-72.1260 42.45, -72.123 42.1546)'::geometry, 3857)
                );
st_distance
-----------------
167.441410065196

-- ジオメトリの例 - メートル単位 (SRID: 3857 上と同じですが、歪みを考慮してcos(lat)で修正)
SELECT ST_Distance(
                        ST_Transform('SRID=4326;POINT(-72.1235 42.3521)'::geometry, 3857),
                        ST_Transform('SRID=4326;LINESTRING(-72.1260 42.45, -72.123 42.1546)'::geometry, 3857)
                ) * cosd(42.3521);
st_distance
-----------------
123.742351254151

-- ジオメトリの例 - メートル単位 (SRID: 26986 メートル単位マサチューセッツ州平面) (マサチューセッツで最も精度が高い)
SELECT ST_Distance(
                        ST_Transform('SRID=4326;POINT(-72.1235 42.3521)'::geometry, 26986),
                        ST_Transform('SRID=4326;LINESTRING(-72.1260 42.45, -72.123 42.1546)'::geometry, 26986)
                );
st_distance
-----------------
123.797937878454

-- ジオメトリの例 - メートル単位 (SRID: 2163 米国ナショナルアトラス正積図法)
(least accurate)
SELECT ST_Distance(
                        ST_Transform('SRID=4326;POINT(-72.1235 42.3521)'::geometry, 2163),
                        ST_Transform('SRID=4326;LINESTRING(-72.1260 42.45, -72.123 42.1546)'::geometry, 2163)
                );

st_distance
------------------
126.664256056812

ジオグラフィの例

-- ジオメトリの例と同じですが、単位がメートルである点に注意して下さい。
-- 若干の速度向上と精度低下となる球面を使った計測です。
SELECT ST_Distance(gg1, gg2) As spheroid_dist, ST_Distance(gg1, gg2, false) As sphere_dist
FROM (SELECT
        'SRID=4326;POINT(-72.1235 42.3521)'::geography as gg1,
        'SRID=4326;LINESTRING(-72.1260 42.45, -72.123 42.1546)'::geography as gg2
        ) As foo  ;

  spheroid_dist   |   sphere_dist
------------------+------------------
 123.802076746848 | 123.475736916397

関連情報

ST_3DDistance, ST_DWithin, ST_DistanceSphere, ST_DistanceSpheroid, ST_MaxDistance, ST_HausdorffDistance, ST_FrechetDistance, ST_Transform