Oracle根据经纬度查询一定范围内的数据

指定一个经纬度，给定一个范围值(单位:千米)，查出在经纬度周围这个范围内的数据。 
经度:23.33 
纬度:34.66 
范围:2000km 
SA_LONGITUDE为数据表经度字段 
SA_LATITUDE为数据表纬度字段 

 

select * from SY_STORE_ADDRESS where sqrt(  

    (  
     ((23.33-SA_LONGITUDE)*ACOS(-1)*12656*cos(((23.33+SA_LONGITUDE)/2)*ACOS(-1)/180)/180)  
     *  
     ((23.33-SA_LONGITUDE)*ACOS(-1)*12656*cos (((23.33+SA_LONGITUDE)/2)*ACOS(-1)/180)/180)  
    )  
    +  
    (  
     ((34.66-SA_LATITUDE)*ACOS(-1)*12656/180)  
     *  
     ((34.66-SA_LATITUDE)*ACOS(-1)*12656/180)  
    )  
)<2000