3分鐘了解Mysql空間搜GeoHash

簡單介紹：

   Mysql 内置函數方案，适合于已有業務，新增加LBS功能，增加經緯度字段方可，避免資料遷移，在5.7.5後實作更多功能實作INNODB的空間搜方法，之前版本主要是對MYISAM的支援。

    在此之前，InnoDB将幾何資料存儲為BLOB（二進制大對象）資料，在空間資料上隻能建立字首索引，當涉及空間搜尋時非常低效，尤其是在涉及複雜的幾何資料時。在大多數情況下，獲得結果的唯一方式是掃描表。

    新版本MySQL中，InnoDB支援空間索引，通過R樹來實作，使得空間搜尋變得高效，如使用内置函數(MBRWITHIN MBRCONTAINS)效率非常好。但目前空間索引隻支援兩個次元的資料。

應用場景介紹:

    公司新業務的需求，希望能夠實作目前位置快速顯示共享資訊資料，針對5.7新特性的支援度調研，利用GeoHash封裝成内置資料庫函數的簡易方案，實作前期的初期業務。

實作過程:

  1、建構表

    CREATE TABLE `tongzhou` (

      `id` INT(11) NOT NULL AUTO_INCREMENT,

      `cname` VARCHAR(255) COLLATE utf8mb4_unicode_ci DEFAULT NULL,

      `oint` POINT NOT NULL,

      PRIMARY KEY (`id`),

      SPATIAL KEY `sp_index` (`oint`)

    ) ENGINE=INNODB AUTO_INCREMENT=1 DEFAULT CHARSET=utf8mb4 COLLATE=utf8mb4_unicode_ci

  2、構造簡單資料   

    insert into `tongzhou` (`cname`, oint) values('景欣園',GEOMFROMTEXT('POINT(39.8885917679 116.6576038966)'));

    insert into `tongzhou` (`cname`, oint) values('京通羅斯福廣場',GEOMFROMTEXT('POINT(39.8890214887 116.6473661241)'));

    insert into `tongzhou` (`cname`, oint) values('中國北京市北京市通州區',GEOMFROMTEXT('POINT(39.8898670523 116.6566729546)'));

    insert into `tongzhou` (`cname`, oint) values('中國北京市北京市通州區',GEOMFROMTEXT('POINT(39.8883852752 116.655728817)'));

    insert into `tongzhou` (`cname`, oint) values('楊莊路22号院',GEOMFROMTEXT('POINT(39.8984936518 116.6339063644)'));

    insert into `tongzhou` (`cname`, oint) values('中國北京市北京市朝陽區 ',GEOMFROMTEXT('POINT(39.8975388526 116.613779068)'));

    insert into `tongzhou` (`cname`, oint) values('北京市朝陽區第三醫院 約31米',GEOMFROMTEXT('POINT(39.8873809413 116.603307724)'));

    insert into `tongzhou` (`cname`, oint) values('雙橋溫泉北裡小區',GEOMFROMTEXT('POINT(39.8918427053 116.6076636314)'));

    insert into `tongzhou` (`cname`, oint) values('京客隆配送中心',GEOMFROMTEXT('POINT(39.8916616061 116.5908622742)'));

    insert into `tongzhou` (`cname`, oint) values('世紀飯店',GEOMFROMTEXT('POINT(39.8918591688 116.6024065018)'));

    insert into `tongzhou` (`cname`, oint) values('雙橋六号井小區-北區 約48米',GEOMFROMTEXT('POINT(39.8899164443 116.6050457954)'));

    insert into `tongzhou` (`cname`, oint) values('翠屏裡小區 約190米',GEOMFROMTEXT('POINT(39.8883194176 116.6495060921)'));

    insert into `tongzhou` (`cname`, oint) values('京通羅斯福廣場',GEOMFROMTEXT('POINT(39.8894719148 116.6584646702)'));

    insert into `tongzhou` (`cname`, oint) values('life新生活廣場 約87米',GEOMFROMTEXT('POINT(39.8895542353 116.65579319)'));

    insert into `tongzhou` (`cname`, oint) values('life新生活廣場 約74米',GEOMFROMTEXT('POINT(39.8896242077 116.6568982601)'));

    insert into `tongzhou` (`cname`, oint) values('life新生活廣場',GEOMFROMTEXT('POINT(39.8901428239 116.6564154625)'));

    insert into `tongzhou` (`cname`, oint) values('漫春園',GEOMFROMTEXT('POINT(39.8913364496 116.658115983)'));

    insert into `tongzhou` (`cname`, oint) values('7天連鎖酒店(北京通州果園環島店) 約58米',GEOMFROMTEXT('POINT(39.8915998677 116.6560935974)'));

    insert into `tongzhou` (`cname`, oint) values('北京金松飯店 約77米',GEOMFROMTEXT('POINT(39.9036172391 116.6577833891)'));

    insert into `tongzhou` (`cname`, oint) values('新華南路小區-東區 約32米',GEOMFROMTEXT('POINT(39.9034444015 116.6588240862)'));

    insert into `tongzhou` (`cname`, oint) values('紅旗小區 約61米',GEOMFROMTEXT('POINT(39.9046707164 116.6598540545)'));

    insert into `tongzhou` (`cname`, oint) values('祥雲天地家園 約62米',GEOMFROMTEXT('POINT(39.9067899674 116.659129858)'));

    insert into `tongzhou` (`cname`, oint) values('梨園地鐵站 約40米',GEOMFROMTEXT('POINT(39.8832975966 116.6687965393)'));

  3、搜尋目前1公裡範圍内的位置資訊

    SET @jl=1;  //多少公裡範圍内的

    SET @jd=118.1964111328;

    SET @wd=39.5591182422;  //目前的位置資訊

    oint是實際已經存在的經緯度資訊

     SELECT *,ASTEXT(oint)FROM  tongzhou WHERE  MBRCONTAINS( LINESTRING(POINT( @wd + @jl / ( 111.12 / COS(RADIANS(@jd))), 

    @jd + @jl / 111.12  ),  

    POINT  ( @wd - @jl / ( 111.12 / COS(RADIANS(@jd))),@jd - @jl / 111.12)  ), oint) 

    ###經緯弧度（1° latitude = 111.12 kilometers）即 10/111.12，表示查找附近10公裡。

  4、測試

    通過百度地圖擷取一些資料，對

         SET @jd=118.1964111328

         SET @wd=39.5591182422

        值的更新，分析擷取範圍内的景點資訊

  5、問題：

     如果直接在WHERE采用一系列的表達式，這樣會導緻無法使用空間索引，如圖所示:

<a href="https://s2.51cto.com/oss/201711/11/2fa268136801e296920360f50ace6d02.png-wh_500x0-wm_3-wmp_4-s_2852158178.png" target="_blank"></a>

解決辦法:

     1、把WHERE後的表達式處理成一個變量進行引用，然後進行周邊搜    

    SET @aa=(SELECT LINESTRING(POINT( 39.5591182422 + 10/ ( 111.12 / COS(RADIANS(118.1964111328))), 

    118.1964111328 + 10 / 111.12  ),  

    POINT  ( 39.5591182422 - 10 / ( 111.12 / COS(RADIANS(118.1964111328))),118.1964111328 -10 / 111.12)) AS  heji)

     2、執行周邊搜

       SELECT ASTEXT(oint)FROM  tongzhou WHERE  MBRWITHIN(@aa,oint)

   順帶執行計劃效果:

其他SQL寫法補充: 

    表示3公裡範圍内的點:

      DDL：       

    CREATE TABLE `shop_id` (

      `id` int(11) NOT NULL AUTO_INCREMENT,

      `lat` double DEFAULT NULL,

      `lng` double DEFAULT NULL,

      `oint` point NOT NULL,

      SPATIAL KEY `sp_index` (`oint`)

    ) ENGINE=InnoDB AUTO_INCREMENT=275915 DEFAULT CHARSET=utf8mb4 COLLATE=utf8mb4_unicode_ci

    SELECT id, ST_Distance_Sphere(POINT(-73.951368, 40.716743), oint) AS dist, ST_AsText(oint) FROM shop_id

     WHERE ST_Contains( ST_MakeEnvelope( POINT((-73.951368+(3/111)), (40.716743+(3/111))),

      POINT((-73.951368-(3/111.12)), (40.716743-(3/111.12))) ), oint ) 

本文轉自 DBAspace 51CTO部落格，原文連結:http://blog.51cto.com/dbaspace/1980877