MySQL空間函數(shù)ST_Distance_Sphere()的使用方式

更新時(shí)間：2022年11月17日 15:37:24 作者：站在墻頭上

這篇文章主要介紹了MySQL空間函數(shù)ST_Distance_Sphere()的使用方式，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。如有錯(cuò)誤或未考慮完全的地方，望不吝賜教

空間函數(shù)ST_Distance_Sphere()的使用

返回球體上兩個(gè)點(diǎn)和/或多點(diǎn)之間的最小球面距離（以米為單位）

官網(wǎng)（https://dev.mysql.com/doc/refman/5.7/en/spatial-convenience-functions.html）

測(cè)試

表結(jié)構(gòu)如下：

create table t_geo_test
(
    ID    int auto_increment
        primary key,
    NAME  varchar(64) not null,
    SHAPE geometry    not null
)
    engine = InnoDB
    charset = utf8;

測(cè)試數(shù)據(jù)：(插入坐標(biāo)時(shí)，以下兩種寫(xiě)法都可以)

INSERT INTO t_geo_test VALUES (1, '測(cè)試', POINT(121.590347, 31.388094));
INSERT INTO t_geo_test VALUES (2, 'somewhere', ST_GeomFromText('POINT(121.366961 31.190049)'));

計(jì)算距離：

SELECT st_distance_sphere(POINT(121.590347, 31.388094),SHAPE) AS distant FROM t_geo_test;

查詢到某點(diǎn)的距離小于一個(gè)值的所有數(shù)據(jù)：

SELECT *, ST_Distance_Sphere(POINT(121.590347, 31.388094),SHAPE) AS distant FROM t_geo_test WHERE ST_Distance_Sphere(POINT(121.590347, 31.388094),SHAPE) < 1000 ORDER BY distant;

注意：

坐標(biāo)不能隨意寫(xiě)測(cè)試數(shù)據(jù)，注意官網(wǎng)中的說(shuō)的條件：

幾何參數(shù)應(yīng)由指定（經(jīng)度，緯度）坐標(biāo)值的點(diǎn)組成：
經(jīng)度和緯度分別是該點(diǎn)的第一和第二坐標(biāo)。
兩個(gè)坐標(biāo)均以度為單位。
經(jīng)度值必須在（-180，180]范圍內(nèi)。正值位于本初子午線以東。
緯度值必須在[-90，90]范圍內(nèi)。正值位于赤道以北。

否則則會(huì)報(bào)錯(cuò)：

[HY000][1210] Incorrect arguments to st_distance_sphere

st_distance_sphere計(jì)算兩坐標(biāo)點(diǎn)距離

最近項(xiàng)目中需要計(jì)算一個(gè)坐標(biāo)點(diǎn)與多個(gè)點(diǎn)的距離，發(fā)現(xiàn)用Python來(lái)實(shí)現(xiàn)效率很低。經(jīng)同事推薦，將這些坐標(biāo)點(diǎn)存入了Mysql數(shù)據(jù)庫(kù)，然后用數(shù)據(jù)庫(kù)自帶的方法st_distance_sphere計(jì)算距離。經(jīng)過(guò)比較發(fā)現(xiàn)確實(shí)效率提高了很多，特此記錄一下。

注意：使用st_distance_sphere需要Mysql數(shù)據(jù)庫(kù)版本為5.7及以上。

例子

現(xiàn)數(shù)據(jù)庫(kù)表中有如下坐標(biāo)點(diǎn)數(shù)據(jù)：

需要找出與第一個(gè)坐標(biāo)點(diǎn)（113.8064049, 22.7300434）相距小于3500米的坐標(biāo)點(diǎn)（表中共有3500個(gè)點(diǎn)）。

sql實(shí)現(xiàn)：

select t.num,t.city,t.wgs84_lng,t.wgs84_lat,
  TRUNCATE(st_distance_sphere(point (113.8064049, 22.7300434),point(t.wgs84_lng,t.wgs84_lat)),2) as distance
	from moran_point t
	where t.city = '深圳' 
	HAVING distance > 0 and distance < 3500
	ORDER BY distance;

查詢結(jié)果：