(PHP)携帯電話から取得したGPS位置情報から半径xxKm内のデータを探したい

最終的にしたいことは、携帯から取得したGPS座標から任意の距離内にある施設(DBに登録)を探したいということです。

DBに置いておく施設情報には、位置情報が入っていないといけないため、緯度、経度を入れておくフィールドを用意してやります。

今回は MySQL で、"latitude" , "longitude" という DOUBLE 型のフィールドを用意。また、"id","name" 列もあるとします。

で、携帯電話で計測した位置から半径5km以内のデータを探すという処理してます。

(なお、このソースはここの続きなので、"$lat" や "$lon" は携帯から取得した座標が入っているとします。)


//$lat,$lon は時小数点形式の緯度経度が入っているものとする
//DB接続は省略
 
//クエリ生成
$query = "SELECT * FROM　store WHERE latitude IS NOT NULL AND longitude IS NOT NULL ";
//クエリ実行
$result = mysql_query($query);
//条件にあうデータを格納する配列
$findAry = array();
while ($row = mysql_fetch_assoc($result)) {
    if (isset($lat) && isset($lon)){
        //GPSからのクエリの場合距離計算(m)
        $dist = distance( $lat , $lon , $row["latitude"] ,$row["longitude"] ,true);
        //debug
        //echo "dist:" . $dist . " lat:" . $row["latitude"] . " lot:" . $row["longitude"] . "
";
        //指定距離内なら配列追加 (今回は5km)
        if ($dist < 5000){
            array_push($findAry ,array( "id" => $row["id"] , 
                                    "name" => mb_convert_encoding(htmlspecialchars($row["name"]),"sjis-win" , "UTF-8") , //DB,PHPがUFT8のため携帯向けのshift-jisに変換
                                    "lat" => $row["latitude"] , 
                                    "lon" => $row["longitude"] ,     ));
        }
    }
}
//$findAryに結果が入ってるので、これをSmartyに渡すなり、そのまま出力したりする。
 
 
/**
 * 2地点間の直線距離を求める（Hubenyの簡易式による）携帯電話GPSの測地系はwgs84なので世界測地系使用
 * @param string $a_lati A地点の緯度(時小数点形式 33.000000)
 * @param string $a_long A地点の経度(時小数点形式 132.000000)
 * @param string $b_lati B地点の緯度(時小数点形式 33.000000)
 * @param string $b_long B地点の経度(時小数点形式 132.000000)
 * @return double 直線距離（メートル）
*/
function calc_distance($a_lati , $a_long , $b_lati , $b_long ) {
    //ラジアンに変換
    $a_long = deg2rad($a_long);
    $a_lati = deg2rad($a_lati);
    $b_long = deg2rad($b_long);
    $b_lati = deg2rad($b_lati);
 
    $latave = ($a_lati + $b_lati) / 2;
    $latidiff = $a_lati - $b_lati;
    $longdiff = $a_long - $b_long;
 
    //子午線曲率半径
    //$meridian = 6334834 / sqrt(pow(1 - 0.006674 * sin($latave) * sin($latave), 3));        //日本測地系
    $meridian = 6335439 / sqrt(pow(1 - 0.006694 * sin($latave) * sin($latave), 3));      //世界測地系
 
    //卯酉線曲率半径
    //$primevertical = 6377397 / sqrt(1 - 0.006674 * sin($latave) * sin($latave));       //日本測地系
    $primevertical = 6378137 / sqrt(1 - 0.006694 * sin($latave) * sin($latave));       //世界測地系
 
    //Hubenyの簡易式
    $x = $meridian * $latidiff;
    $y = $primevertical * cos($latave) * $longdiff;
 
    return sqrt($x * $x + $y * $y);
}