求助：如何求地球上2直线间夹角？

已知地球上3点Point3D pt1, Point3D pt2, Point3D pt3,(Point3D.X=经度，Point3D.Y=纬度，Point3D.Z=高度)如何求这3点形成的角度呢（pt1-pt2-pt3, pt2是角度的顶点）？

有什么算法大家谁能提供下吗？谁能帮助填下下面的函数吗？

/**

*求地球上3点形成的角度

*

*param[in]  pt1 

*param[in]  pt2  角度顶点

*param[in]  pt3 

*

*reture 直线角度

*/

pulic double GetEarth3DAngle(Point3D pt1, Point3D pt2, Point3D pt3)

{

 

 

}