SQL 2座標距離計算後のソート

2つの座標の列が(x 1,y 1)、(x 2,y 2)である場合、それらの間の距離:
SQRT((X1-X2)*(X1-X2)+(Y1-Y2)*(Y1-Y2))

SELECT * FROM m_store 
ORDER BY 
SQRT((121.517759-`longitude`)*(121.517759-`longitude`)+(31.178469-`latitude`)*(31.178469-`latitude`))  

longitude経度
义齿

/**
*             ,    
*@param lng1,lng2   
*@param lat1,lat2   
*@return float   ，   
*@author www.Alixixi.com
**/
function getdistance($lng1,$lat1,$lng2,$lat2){
	//       
	$radLat1=deg2rad($lat1);//deg2rad()          
	$radLat2=deg2rad($lat2);
	$radLng1=deg2rad($lng1);
	$radLng2=deg2rad($lng2);
	$a=$radLat1-$radLat2;
	$b=$radLng1-$radLng2;
	$s=2*asin(sqrt(pow(sin($a/2),2)+cos($radLat1)*cos($radLat2)*pow(sin($b/2),2)))*6378.137*1000;
	return $s;
}


echo getdistance(121.477551,31.270338, 118.782347,32.072796).' ';