2つのポイントが互いに1マイル以内にあるかどうかを確認する非常に効率的なチェックを実装しようとしています。
私の現在のアプローチは、Haversine距離を計算し、それが1マイル未満かどうかを確認することです。
大きなレコードセットの場合は、このyes / noフラグを計算する必要があるため、この場合は効率が重要になります。
私はそれらが1マイル以内にあるかどうかだけを気にします -距離については何も私には関係ありません。
それでは、2つの緯度/経度のポイントが互いに1マイル以内にあるかどうかを確認する最も効率的な方法は何ですか？
コメントに応えて、私はSQL Serverでこれを行っています。私のコードは以下です。
CREATE FUNCTION dbo.USR_UFN_HAVERSINE_DISTANCE
(
@LAT1 FLOAT(18)
,@LONG1 FLOAT(18)
,@LAT2 FLOAT(18)
,@LONG2 FLOAT(18)
,@UnitOfMeasure NVARCHAR(10) = 'KILOMETERS'
)
RETURNS FLOAT(18)
AS
BEGIN
DECLARE
@R FLOAT(8)
,@DLAT FLOAT(18)
,@DLON FLOAT(18)
,@A FLOAT(18)
,@C FLOAT(18)
,@D FLOAT(18)
;
SET @R =
CASE @UnitOfMeasure
WHEN 'MILES' THEN 3956.55
WHEN 'KILOMETERS' THEN 6367.45
WHEN 'FEET' THEN 20890584
WHEN 'METERS' THEN 6367450
ELSE 6367.45 --km
END
SET @DLAT = RADIANS(@LAT2 - @LAT1);
SET @DLON = RADIANS(@LONG2 - @LONG1);
SET @A = SIN(@DLAT / 2)
* SIN(@DLAT / 2)
+ COS(RADIANS(@LAT1))
* COS(RADIANS(@LAT2))
* SIN(@DLON / 2)
* SIN(@DLON / 2);
SET @C = 2 * ASIN(MIN(SQRT(@A)));
SET @D = @R * @C;
RETURN @D;
END;
あなたの研究はこれまでに可能性のある候補として何を見つけましたか？
—
PolyGeo
ソフトウェアソリューションをお探しですか、それとも独自のコードを作成していますか？これまでに何を試しましたか？
—
MaryBeth
ヘイバーシン距離を確認するだけの何が問題になっていますか 平面距離を確認するだけで、処理時間を少し節約できます。
—
トム
geomA.STDistance（geomB）<dを使用できますか？
—
Ian Turton
@Tomによって提案された「平面距離」チェックは、誤って解釈される可能性があります。正しく機能するには、注意深い解釈が必要です。1つは次のとおりです。180度子午線または極でポイントを比較する必要がないと仮定すると、ピタゴラスの公式を座標（lat1、cos（lat1）* lon1）、（lat2、cos（lat2）* lon2）に適用できます。言い換えると、（lat1-lat2）^ 2 +（cos（lat1）* lon1-cos（lat2）* lon2）^ 2を1/69 ^ 2（すべて度）と比較すると、2つの点が1マイル（パーセントの精度で）。これがHaversineより速いかどうかは不明です。
—
whuber