计算点的方位角 (osg)

原文链接

double getAngleByLatLon(double ax, double ay, double bx, double by)
{
ax = ax * osg::PI / 180;
ay = ay * osg::PI / 180;
bx = bx * osg::PI / 180;
by = by * osg::PI / 180;

double temValue = 0;

// 经纬度相同或则纬度相同
if(ax == bx || ay == by)
{
// 经度形同
if(ax == bx)
{
if(by >= ay)
{
return 0;
}
else
{
return 180;
}
}
// 纬度相同
else
{
if(bx >= ax)
{
return 90;
}
else
{
return 270;
}
}
}

tmpValue = sin(ay) * sin(by) + cos(ay) * cos(by) * cos(bx - ax);
tmpValue = sqrt(1 - tmpValue * tmpValue);
tmpValue = cos(by) * sin(bx - ax) / tmpValue;
double angle = abs(asin(tmpValue) * 180 / osg::PI);

if(bx > ax)
{
if(by > ay)
{
// 第一象限
return angle;
}
else
{
// 第二象限
return 180 - angle;
}
}
else
{
if(by >= ay)
{
// 第四象限
return 360 - angle;
}
else
{
// 第三象限
return 180 + angle;
}
}
}

posted @   雾枫  阅读(49)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· PowerShell开发游戏 · 打蜜蜂
· 在鹅厂做java开发是什么体验
· 百万级群聊的设计实践
· WPF到Web的无缝过渡:英雄联盟客户端的OpenSilver迁移实战
· 永远不要相信用户的输入:从 SQL 注入攻防看输入验证的重要性
历史上的今天:
2022-06-22 Qt——将一个文件下的文件拷贝到指定路径文件夹下
2022-06-22 Qt——创建文件夹,指定路径下
2022-06-22 Qr——选择路径,输入文件名称,在该路径下创建该文件
2022-06-22 Qt_消息提示框
2022-06-22 Qt_获取文件路径
2022-06-22 Qt——设置窗口背景,但是不影响窗口上的控件背景
2022-06-22 Qt——模态与非模态
点击右上角即可分享
微信分享提示