public static JWD GetJWDB(JWD A, double distance, double angle)
{
double dx = distance*1000 * Math.Sin(angle * Math.PI /180.0);
double dy = distance*1000 * Math.Cos(angle * Math.PI /180.0);
double BJD = (dx/A.Ed + A.m_RadLo) * 180.0/Math.PI;
double BWD = (dy/A.Ec + A.m_RadLa) * 180.0/Math.PI;
JWD B = new JWD(BJD, BWD);
return B;
}
//! 计算点A 和 点B的经纬度，求他们的距离和点B相对于点A的方位
/*!
* \param A A点经纬度
* \param B B点经纬度
* \param angle B相对于A的方位, 不需要返回该值，则将其设为空
* \return A点B点的角度
*/
public static double angle(JWD A, JWD B)
{
double dx = (B.m_RadLo - A.m_RadLo) * A.Ed;
double dy = (B.m_RadLa - A.m_RadLa) * A.Ec;
double angle = 0.0;
angle = Math.Atan(Math.Abs(dx/dy))*180.0/Math.PI;
// 判断象限
double dLo = B.m_Longitude - A.m_Longitude;
double dLa = B.m_Latitude - A.m_Latitude;
if(dLo > 0 && dLa <= 0) {
angle = (90.0 - angle) + 90.0;
}
else if(dLo <= 0 && dLa < 0) {
angle = angle + 180.0;
}
else if(dLo < 0 && dLa >= 0) {
angle = (90.0 - angle) + 270;
}
return angle;
}
//! 已知点A经纬度，根据B点据A点的距离，和方位，求B点的经纬度
/*!
* /param longitude 已知点A经度
* /param latitude 已知点A纬度
* /param distance B点到A点的距离
* /param angle B点相对于A点的方位
* /return B点的经纬度坐标
*/
public static JWD GetJWDB(double longitude, double latitude, double distance, double angle)
{
JWD A = new JWD(longitude, latitude);
return GetJWDB(A, distance, angle);
}