Javaで経度と緯度がわかっている場合、メートル単位で距離を計算する[複製]

StackOverflow https://stackoverflow.com/questions/837872

質問

  

可能な重複:
   Javaで緯度/経度の値を操作する

複製:

2つの座標で与えられる2点間の距離を計算する必要があります。私が取り組んでいるプロジェクトはJavaプロジェクトですので、Javaコードは素晴らしいでしょうが、擬似コードを与えることもできますし、自分で実装することもできます:)

おそらくご存知のように、座標を表すには3つの方法があります:

  • 度:分:秒(49° 30'00" N、123° 30'00" W)
  • 度:分(49° 30.0 '、-123° 30.0')、(49d30.0m、-123d30.0 ')
  • 10進数(49.5000°、-123.5000°)、通常は4〜6個の10進数。

これは私の座標が与えられる3番目の方法なので、この値のコードが優先されます:)

役に立ちましたか?

解決

stackoverflowに関する別の質問に基づいて、このコードを取得しました。これは、マイルではなくメートルで結果を計算します。 )

 public static float distFrom(float lat1, float lng1, float lat2, float lng2) {
    double earthRadius = 6371000; //meters
    double dLat = Math.toRadians(lat2-lat1);
    double dLng = Math.toRadians(lng2-lng1);
    double a = Math.sin(dLat/2) * Math.sin(dLat/2) +
               Math.cos(Math.toRadians(lat1)) * Math.cos(Math.toRadians(lat2)) *
               Math.sin(dLng/2) * Math.sin(dLng/2);
    double c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a));
    float dist = (float) (earthRadius * c);

    return dist;
    }

他のヒント

GPS用のJava Geodesy Library 、地球表面の曲率を考慮した Vincentyの公式を使用します。

実装は次のようになります。

import org.gavaghan.geodesy.*;

...

GeodeticCalculator geoCalc = new GeodeticCalculator();

Ellipsoid reference = Ellipsoid.WGS84;  

GlobalPosition pointA = new GlobalPosition(latitude, longitude, 0.0); // Point A

GlobalPosition userPos = new GlobalPosition(userLat, userLon, 0.0); // Point B

double distance = geoCalc.calculateGeodeticCurve(reference, userPos, pointA).getEllipsoidalDistance(); // Distance between Point A and Point B

結果の距離はメートル単位です。

C ++では、次のように実行されます。

#define LOCAL_PI 3.1415926535897932385 

double ToRadians(double degrees) 
{
  double radians = degrees * LOCAL_PI / 180;
  return radians;
}

double DirectDistance(double lat1, double lng1, double lat2, double lng2) 
{
  double earthRadius = 3958.75;
  double dLat = ToRadians(lat2-lat1);
  double dLng = ToRadians(lng2-lng1);
  double a = sin(dLat/2) * sin(dLat/2) + 
             cos(ToRadians(lat1)) * cos(ToRadians(lat2)) * 
             sin(dLng/2) * sin(dLng/2);
  double c = 2 * atan2(sqrt(a), sqrt(1-a));
  double dist = earthRadius * c;
  double meterConversion = 1609.00;
  return dist * meterConversion;
}
ライセンス: CC-BY-SA帰属
所属していません StackOverflow
scroll top