Android 2.1 SQLite3で緯度経度から距離検索

Android 2.1のSQLiteにはsin, cos関数がないため例示されているMySQLでの経度緯度からの距離検索はそのままの形では利用できない。解決のヒントをもとに実装したので要点まとめ。

▼元の公式からクエリに利用できる形へ変換
距離 = C * acos ( sin(lat)*sin(qlat)+cos(lat)*cos(qlat)*cos(lng-qlng) )

[lat, lng] 座標1の緯度経度(DB内を想定) …青字はDB格納時に予め算出できる部分
[qlat, qlng] 座標2の緯度経度(基準位置を想定) …緑字はクエリ投入時に算出できる部分
C=6371 (距離の単位がkmの場合) または 3959 (単位がmiの場合)

cos(a-b) = cos(a)*cos(b)+sin(a)*sin(b) なので、

距離 = C * acos ( sin(lat)*sin(qlat)+cos(lat)*cos(qlat)*(cos(lng)*cos(qlng)+sin(lng)*sin(qlng)) )

acos(x) = y は cos(y)=x なので、

cos(距離/C) = sin(lat)*sin(qlat)+cos(lat)*cos(qlat)*(cos(lng)*cos(qlng)+sin(lng)*sin(qlng))

SQLite内のエントリには緯度経度それぞれのsin, cos値を追加で格納しておけば、検索条件とする緯度経度のsin, cos、cos(距離/C)を使った四則演算で距離を評価できる形になる。

このときのcos(距離/C)値は距離として-1が最も遠く+1が最も近い状態。

遠い(地球の裏側) -1 … 0 (地球の1/4) … +1 (同一の緯度経度)

▼SQLiteへのエントリ登録
1つの座標について sin(lat), cos(lat), sin(lng), cos(lng) をそれぞれ real型で登録しておく。元の lat, lng 値はそれぞれ atan(sin/cos) で求まるが、誤差が問題になる場合はあわせて記録。

// テーブル登録例
create table location ( _id integer primary key autoincrement, dat0 integer, ... ,
 sinlat real not null, coslat real not null, sinlng real not null, coslng real not null);

※JavaのMath関数がラジアン基準なので記録の際には注意。

Math.sin(Math.toRadians(loc.getLatitude()));

▼SQLite用 Location距離検索クエリの作成
上述のテーブルを距離条件で検索する場合は、基準とするLocation (latitude, longitude)を元に検索用のクエリを動的に作成する。

// SQLiteDatabase.rawQuery用クエリ文
public String searchNearQuery(Location loc, double range_km){
	double km_cos=Math.cos(range_km/6371);	// 距離基準cos値
	double radlat=Math.toRadians(loc.getLatitude()), radlong=Math.toRadians(loc.getLongitude());
	double qsinlat=Math.sin(radlat), qcoslat=Math.cos(radlat);
	double qsinlng=Math.sin(radlong), qcoslng=Math.cos(radlong);

	StringBuilder sb=new StringBuilder();
	sb.append("SELECT _id, dat0, ..., ");
	sb.append("(sinlat*"+qsinlat+" + coslat*"+qcoslat+"*(coslng*"+qcoslng+"+sinlng*"+qsinlng+")) AS distcos ");
	sb.append(" FROM location ");
	sb.append(" WHERE distcos > "+km_cos);	// 値が大きい方が近い
	sb.append(" ORDER BY distcos DESC ");	// 近い順に出力
//	sb.append(" LIMIT 10");	// ←↑↓必要な場合追加
//	Log.d("searchNearQuery", sb.toString());
	return sb.toString();
}

▼検索結果からの距離の取得
distcos値をacos角度x単位で距離へ変換。

if(distcos>=1.0) distance_km = 0.0;	// 誤差の都合、同一座標で1.0を超える場合があるため(※)
else distance_km = (6371*Math.acos(distcos));	// miの場合は6371の代わりに3959
// ※厳密には地球の裏側で-1を超える場合もある。