sql 根据经纬度查找附近数据

SELECTn    ROUND(n        6378.138 * 2 * ASIN(n            SQRT(n            POW( SIN(( 23.065502 * PI()/ 180-dt_lat * PI()/ 180 )/ 2 ), 2 )+ COS( 23.065502 * PI()/ 180 )* COS( dt_lat * PI()/ 180 )* POW( SIN(( 113.803962 * PI()/ 180- dt_lng * PI()/ 180 )/ 2 ), 2 )))* 1000 n    ) AS distance,n    a.id,n    a.title,n    a.price,n    a.unit,n    a.dt_lat,n    a.dt_lng nFROMn    `cms_house` AS a nHAVINGn    distance <= 1500 nORDER BYn    distance ASCnn

dt_latdt_lng 为表中 经纬度字段