已知點的經緯度 與另一個點的距離 以及方向角 求另一點的經緯度
阿新 • • 發佈：2019-01-24
public class AngleUtil {
/**
* 求B點經緯度
* @param A 已知點的經緯度，
* @param distance AB兩地的距離 單位km
* @param angle AB連線與正北方向的夾角（0~360）
* @return B點的經緯度
*/
final static double Rc=6378137;
final static double Rj=6356725;
double m_LoDeg,m_LoMin,m_LoSec;
double m_LaDeg,m_LaMin,m_LaSec;
double m_Longitude,m_Latitude;
double m_RadLo,m_RadLa;
double Ec;
double Ed;
public AngleUtil(double longitude,double latitude){
m_LoDeg=(int)longitude;
m_LoMin=(int)((longitude-m_LoDeg)*60);
m_LoSec=(longitude-m_LoDeg-m_LoMin/60.)*3600;
m_LaDeg=(int)latitude;
m_LaMin=(int)((latitude-m_LaDeg)*60);
m_LaSec=(latitude-m_LaDeg-m_LaMin/60.)*3600;
m_Longitude=longitude;
m_Latitude=latitude;
m_RadLo=longitude*Math.PI/180.;
m_RadLa=latitude*Math.PI/180.;
Ec=Rj+(Rc-Rj)*(90.-m_Latitude)/90.;
Ed=Ec*Math.cos(m_RadLa);
}
public static String getMyLatLng(AngleUtil A,double distance,double angle){//方法
double dx = distance*1000*Math.sin(Math.toRadians(angle));
double dy= distance*1000*Math.cos(Math.toRadians(angle));
double bjd=(dx/A.Ed+A.m_RadLo)*180./Math.PI;
double bwd=(dy/A.Ec+A.m_RadLa)*180./Math.PI;
String lnglat=bjd+","+bwd;
return lnglat;
}
}