答えて
2点(GPXファイルの各ウェイポイントのペア)間の距離を計算する伝統的な方法は、Haversine式です。
私は、アルゴリズムを実装するSQL Server機能を持っています。これは他の言語にも簡単に翻訳できるはずです:
create function dbo.udf_Haversine(@lat1 float, @long1 float,
@lat2 float, @long2 float) returns float begin
declare @dlon float, @dlat float, @rlat1 float,
@rlat2 float, @rlong1 float, @rlong2 float,
@a float, @c float, @R float, @d float, @DtoR float
select @DtoR = 0.017453293
select @R = 3959 -- Earth radius
select
@rlat1 = @lat1 * @DtoR,
@rlong1 = @long1 * @DtoR,
@rlat2 = @lat2 * @DtoR,
@rlong2 = @long2 * @DtoR
select
@dlon = @rlong1 - @rlong2,
@dlat = @rlat1 - @rlat2
select @a = power(sin(@dlat/2), 2) + cos(@rlat1) *
cos(@rlat2) * power(sin(@dlon/2), 2)
select @c = 2 * atn2(sqrt(@a), sqrt([email protected]))
select @d = @R * @c
return @d
end
これは距離をマイルで返します。キロメートルについては、地球半径をkm相当量に置き換えてください。
Hereは、より詳細な説明です。
編集:この機能は、郵便番号データベースで半径検索を行うのに十分速く、正確です。 this siteについては素晴らしい仕事をしてきました(しかし、今はリンクが壊れているので、これ以上はありません)。
Mike Gavaghan has an algorithm彼のサイトの距離計算に使用します。 C#とJAVA版のコードがあります。
Vincenty formulaeのDelphi実装はhereです。
更新されたリンクがありますか? –
@Adam Carter:URLを確認しても問題ありません。 – menjaraz
ここにScalaの実装があります。
3958.761は、マイル単位でmean radius of the Earthです。 km(またはその他の単位)で結果を取得するには、この番号を変更します。
// The Haversine formula
def haversineDistance(pointA: (Double, Double), pointB: (Double, Double)): Double = {
val deltaLat = math.toRadians(pointB._1 - pointA._1)
val deltaLong = math.toRadians(pointB._2 - pointA._2)
val a = math.pow(math.sin(deltaLat/2), 2) + math.cos(math.toRadians(pointA._1)) * math.cos(math.toRadians(pointB._1)) * math.pow(math.sin(deltaLong/2), 2)
val greatCircleDistance = 2 * math.atan2(math.sqrt(a), math.sqrt(1 - a))
3958.761 * greatCircleDistance
}
// A sequence of gpx trackpoint lat,long pairs parsed from the track GPX data
val trkpts: Seq[(Double, Double)] = {
val x = scala.xml.XML.loadString(track)
(x \\ "trkpt").map(trkpt => ((trkpt \ "@lat").text.toDouble, (trkpt \ "@lon").text.toDouble))
}
// Distance of track in miles using Haversine formula
val trackDistance: Double = {
trkpts match {
case head :: tail => tail.foldLeft(head, 0.0)((accum, elem) => (elem, accum._2 + haversineDistance(accum._1, elem)))._2
case Nil => 0.0
}
}
この質問はかなり古いですが、私は完全性のためにpythonオプションを追加したいと思います。 GeoPyはgreat-circle distance
とVincenty distance
の両方を持っています。
- 1. 加速度計からの距離の計算方法
- 2. パスからの最短距離の計算方法は?
- 3. 対象距離の計算方法
- 4. MySQLでの距離計算の距離計算
- 5. 重心からの距離の計算
- 6. 個別のdfからの距離の計算方法R
- 7. gpxファイルからの正確なカウント距離 - ph
- 8. ポリゴンから地理距離を計算する方法は?
- 9. 距離計算のバグandroid
- 10. クラウドでの距離計算
- 11. NS-3の距離計算
- 12. (正方形)チェックイン距離の計算
- 13. ビデオの動きからの距離を計算する方法は?
- 14. OpenLayers:2点間の距離の計算方法は?
- 15. OpenCVの焦点距離の計算方法は?
- 16. 5つのマーカー間の距離を計算する方法は?
- 17. 2点以上の距離の計算方法は? C#
- 18. ドキュメントと各重心(k-means)の距離の計算方法は?
- 19. ポストギスによる都市間の距離の計算方法は?
- 20. 距離計算建物
- 21. Python - 疎ベクトル/距離計算
- 22. Android GPS距離計算
- 23. ショートパス。距離計算。 java。
- 24. 道路距離計算mysql
- 25. 計算ユークリッド距離が
- 26. Estimoteビーコン距離計算
- 27. Geokit/Geocoderでの距離と距離の計算
- 28. 旅行総距離の計算方法は?
- 29. 距離とは別の場所を計算する方法
- 30. 超球面の中心からの距離の計算方法(白色化PCA)
ありがとうございます。私はそれをjavaに移植し、ここに投稿します。 @DtoRとはどういう意味ですか?半径までの距離? – guerda
これは、度をラジアン、π/ 180に変換する要因です。 – cdonner
リンクをありがとう。すぐにJavaバージョンを投稿します – guerda