大地坐标系法截线与大地线.pptxVIP

第1页/共30页大地坐标系法截线与大地线第2页/共30页5.3 大地坐标系与大地空间直角坐标系三、同一参考椭球下大地坐标与空间直角坐标的转换1、（B，L，H）→（X，Y，Z）参考椭球第3页/共30页5.3 大地坐标系与大地空间直角坐标系三、同一参考椭球下大地坐标与空间直角坐标的转换2、（X，Y，Z） → （B，L，H）参考椭球第4页/共30页迭代求解法，初始值：收敛条件为：迭代收敛解为：5.3 大地坐标系与大地空间直角坐标系三、同一参考椭球下大地坐标与空间直角坐标的转换2、（X，Y，Z） → （B，L，H）说明：1）ε为一小正数，如ε=5×10-10 ；2）J 为迭代收敛时的迭代次数。参考椭球第5页/共30页5.3 大地坐标系与大地空间直角坐标系三、同一参考椭球下大地坐标与空间直角坐标的转换1、（B，L，H）→（X，Y，Z）2、（X，Y，Z） → （B，L，H）参考椭球第6页/共30页上机编程计算实习：第一部分1、 （B，L，H） →（X，Y，Z） 算例已知某点1954年北京坐标系的大地坐标为：试编程求该点1954年北京坐标系的大地空间直角坐标。（以米为单位输出到屏幕，保留小数点后3位）注：1954年北京坐标系使用的是克拉索夫斯基椭球。第7页/共30页 要求1：从文件中读入已知数据，文件名及文件内容为BLHToXYZ.txt:6378245.0298.3334455.666771122.3335555.660文件名aα=1：298.3BLH 要求2：每个功能模块单独编写一个子函数，如度分秒 化弧度上机编程计算实习：第一部分1、 （B，L，H） →（X，Y，Z）第8页/共30页 算例已知某点1954年北京坐标系的大地空间直角坐标为：试编程求该点1954年北京坐标系的大地坐标？（角度以度、分、秒形式输出到屏幕，长度以米单位输出到屏幕，均保留小数点后3位）注：1954年北京坐标系使用的是克拉索夫斯基椭球。上机编程计算实习：第一部分2、 （X，Y，Z） →（B，L，H）第9页/共30页 要求1：从文件中读入已知数据，文件名及文件内容为XYZToBLH.txt:6378245.02985325181238.3903526461.538文件名aα=1：298.3XYZ 要求2：每个功能模块单独编写一个子函数，如弧度化度分秒上机编程计算实习：第一部分2、 （X，Y，Z） →（B，L，H）第10页/共30页 运行环境：Visual C++6.0 构成主程序：大地测量学基础程序框架.cpp#includesubroutine.h″子程序：subroutine.h上机编程计算实习：第一部分3、 关于实习程序框架的说明第11页/共30页主程序：大地测量学基础程序框架.cpp//C和C++的头文件#include math.h#include stdlib.h#include stdio.h#include iostream.h//自定义常数#define Ru 206264.806247096#define PI 3.141592653589793//自己编的子函数，请放在头文件subroutine.h中#include subroutine.h头部分第12页/共30页void main(){ …… //从文件中读入已知数据并检核 FILE *fpin; fpin=fopen(BLHToXYZ.txt,r); double a,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); //调用自己的子函数 BLHToXYZ(a,alfa,B,L,H,X,Y,Z); printf(==\nX=%12.3f\nY=%12.3f\nZ=%12.3f\n,X,Y,Z); …… }主体部分主程序：大地测量学基础程序框架.cpp第13页/共30页子程序：subroutine.hdouble RAD(double dfm)/* 功能： 度分秒化弧度子函数 输入：dfm： 格式为 度分秒.小数秒 返回值：弧度值 版本：1.0 作者： 刘长建2007.01.30*/{ ……}void BLHToXYZ(double a,double alfa,double B,double L,double H,double X,double Y,double Z)/*……第14页/共30页5.4 法截线与大地线normal section line and geodetic line子午圈：赤 道：其