문제

가능한 복제 :
Java에서 위도/경도 값으로 작업합니다

복제하다:

두 좌표로 주어진 두 점 사이의 거리를 계산해야합니다. 내가 작업하고있는 프로젝트는 Java 프로젝트이므로 Java 코드가 훌륭 할 것이지만 의사 코드도 제공 할 수 있습니다. 그러면 직접 구현할 수 있습니다 :)

아시다시피, 좌표를 나타내는 세 가지 방법이 있습니다.

  • 도 : 분 : 초 (49 ° 30'00 "N, 123 ° 30'00"W)
  • 도 : 소수점 (49 ° 30.0 ', -123 ° 30.0'), (49d30.0m, -123d30.0 ')
  • 소수도 (49.5000 °, -123.5000 °), 일반적으로 4-6 소수점 숫자.

그것은 내 좌표가 주어진 세 번째 방법 이므로이 값에 대한 코드가 선호됩니다 :)

도움이 되었습니까?

해결책

기반 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 라이브러리, 그것은 사용합니다 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