资源描述:
《[06][ch5][大地坐标系&法截线与大地线1].ppt》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、第六讲大地坐标系法截线与大地线内容回顾球面直角三角形的纳白尔规则参考椭球定义及其作用参考椭球几何参数间的相互关系简述地球椭球、参考椭球、正常椭球、总地球椭球的区别与联系椭球面上点的法线长公式(会推)大地坐标、大地空间直角坐标的定义请思考野外测量可以获得哪些观测值,点的最终水平坐标和高程怎么得到?5.3大地坐标系与大地空间直角坐标系三、同一参考椭球下大地坐标与空间直角坐标的转换1、(B,L,H)→(X,Y,Z)参考椭球5.3大地坐标系与大地空间直角坐标系三、同一参考椭球下大地坐标与空间直角坐标的转换2、(X,Y,Z)→(B
2、,L,H)参考椭球5.3大地坐标系与大地空间直角坐标系三、同一参考椭球下大地坐标与空间直角坐标的转换2、(X,Y,Z)→(B,L,H)迭代求解法,初始值:收敛条件为:迭代收敛解为:参考椭球说明:1)ε为一小正数,如ε=5×10-10;2)J为迭代收敛时的迭代次数。5.3大地坐标系与大地空间直角坐标系三、同一参考椭球下大地坐标与空间直角坐标的转换2、(X,Y,Z)→(B,L,H)参考椭球1、(B,L,H)→(X,Y,Z)1、(B,L,H)→(X,Y,Z)已知某点1954年北京坐标系的大地坐标为:试编程求该点1954年北京坐标
3、系的大地空间直角坐标。(以米为单位输出到屏幕,保留小数点后3位)注:1954年北京坐标系使用的是克拉索夫斯基椭球。算例上机编程计算实习:第一部分要求1:从文件中读入已知数据,文件名及文件内容为BLHToXYZ.txt:6378245.0298.3334455.666771122.3335555.660要求2:每个功能模块单独编写一个子函数,如度分秒化弧度文件名aα=1:298.3BLH上机编程计算实习:第一部分1、(B,L,H)→(X,Y,Z)2、(X,Y,Z)→(B,L,H)算例已知某点1954年北京坐标系的大地空间直角
4、坐标为:试编程求该点1954年北京坐标系的大地坐标?(角度以度、分、秒形式输出到屏幕,长度以米单位输出到屏幕,均保留小数点后3位)注:1954年北京坐标系使用的是克拉索夫斯基椭球。上机编程计算实习:第一部分要求1:从文件中读入已知数据,文件名及文件内容为要求2:每个功能模块单独编写一个子函数,如弧度化度分秒XYZToBLH.txt:6378245.0298.31178143.5325181238.3903526461.538文件名aα=1:298.3XYZ上机编程计算实习:第一部分2、(X,Y,Z)→(B,L,H)3、关于
5、实习程序框架的说明构成主程序:大地测量学基础程序框架.cpp子程序:subroutine.h#include"subroutine.h″运行环境:VisualC++6.0上机编程计算实习:第一部分//C和C++的头文件#include#include#include#include//自定义常数#defineRu206264.806247096#definePI3.141592653589793//自己编的子函数,请放在头文件subroutine
6、.h中#include"subroutine.h"头部分主程序:大地测量学基础程序框架.cpp主体部分voidmain(){……//从文件中读入已知数据并检核FILE*fpin;fpin=fopen("BLHToXYZ.txt","r");doublea,alfa,B,L,H,X,Y,Z;fscanf(fpin,"%lf%lf%lf%lf%lf",&a,&alfa,&B,&L,&H);fclose(fpin);//调用子程序并显示计算结果alfa=1.0/alfa;B=RAD(B);//调用自己的子函数L=RAD(L);/
7、/调用自己的子函数BLHToXYZ(a,alfa,B,L,H,X,Y,Z);printf("==>X=%12.3fY=%12.3fZ=%12.3f",X,Y,Z);……}主程序:大地测量学基础程序框架.cppdoubleRAD(doubledfm)/*功能:度分秒化弧度子函数输入:dfm:格式为度分秒.小数秒返回值:弧度值版本:1.0作者:刘长建2007.01.30*/{……}voidBLHToXYZ(doublea,doublealfa,doubleB,doubleL,doubleH,double&X,d
8、ouble&Y,double&Z)/*……子程序:subroutine.h一、任意方向法截线曲率半径1、有关定义法截面:包含椭球面某点法线的平面。法截线:法截面与参考椭球面的交线。斜截面:不包含椭球面某点法线的平面。斜截线:斜截面与参考椭球面的交线。5.4法截线与大地线normalsectionline