mysql计算距离

SELECT 
  d.*,
  round(
    1000 * 2 * 6378.137 * ASIN(
      SQRT(
        POW(
          SIN(
            PI() * (113.92474300 - longitude) / 360
          ),
          2
        ) + COS(PI() * 22.67503500 / 180) * COS(latitude * PI() / 180) * POW(
          SIN(
            PI() * (22.67503500 - latitude) / 360
          ),
          2
        )
      )
    )
  , 0) AS juli 
FROM
  `d_device`  d
  where longitude is not null and latitude is not null
ORDER BY juli ASC 
LIMIT 0, 20 ;