sql 根据经纬度查询门店地理位置

SELECT
    s.id ,
    s.name ,
    s.lng ,
    s.lat ,
    ROUND(
        6378.138 * 2 * ASIN(
            SQRT(
                POW(
                    SIN( ( @targetLat * PI() / 180 - s.lat * PI() / 180 ) / 2 ) , 2 )
                +
                COS( @targetLat * PI( ) / 180 ) * COS( s.lat * PI( ) / 180 )
                * POW( SIN( ( @targetLng * PI() / 180 - s.lng * PI() / 180 ) / 2 ) , 2 )
            )
        ) * 1000
    ) AS distance
FROM station s
ORDER BY distance ASC , s.id
LIMIT 20;