


下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、学习好资料欢迎下载空间直角坐标系与大地坐标系转换程序#include<iostream>#include<cmath>#include<iomanip>usingnamespacestd;#definePI(2.0*asin(1.0)voidmain()doublea,b,c,d1,d2,f1,f2,m1,m2BL,H,X,Y,Z,W,N,e;/cout<<"请分别输入椭球的长半轴、短半轴(国际单位)"<<endl;cin>>a>>b;a=6378137;/以WGS84为例b=6356752
2、.3142;e=sqrt(a*a-b*b)/a;c=a*a/b;intx;cout<<"请输入0或1,0:大地坐标系到空间直角坐标系;1:空间直角坐标系到大地坐标系"<<endl;cin>>x;switch(x)case0:cout<<"请分别输入该点大地纬度、经度、大地高(国际单位,纬度经度请按度分秒,分别输入)"<<endl;cin>>d1>>f1>>m1>>d2>>f2>>m2>>H;B=PI*(d1+f1
3、/60+m1/3600)/180;L=PI*(d2+f2/60+m2/3600)/180;W=sqrt(1-e*e*sin(B)*sin(B);N=a/W;X=(N+H)*cos(B)*cos(L);Y=(N+H)*cos(B)*sin(L);Z=(N*(1-e*e)+H)*sin(B);cout<<"空间直角坐标系中X,YZ,坐标值(国际单位)分另1J为"<<fixed<<setprecision(6)<<X<<""<<fixed<<setprecision(6)<
4、;<Y<<""<<fixed<<setprecision(6)<<Z<<endl;break;case1:cout<<"请分别输入空间直角坐标系中X,YZ的值(国际单位)"<<endl;cin>>X>>Y>>Z;doublet,m,n,P,k,B0;m=Z/sqrt(X*X+Y*Y);t0B0=atan(m);/初值n=Z/sqrt(X*X+Y*Y);P=c*e*e/sqrt(X*X+Y*Y);k=1+(a*a-b*b)/(b*b
5、);t=m+P*n/sqrt(k+n*n);/现在为t1,之后代替t2,t3.B=atan(t);W=sqrt(1-e*e*sin(B)*sin(B);N=a/W;H=Z/sin(B)-N*(1-e*e);inti;for(i=1;fabs(B-B0)>10E-10;i+)每一次新的B与上一次计算的B比较,误差小于10E-10radBO=B;n=t;t=m+P*n/sqrt(k+n*n);/迭代B=atan(t);W=sqrt(1-e*e*sin(B)*sin(B);N=a/W;/if(X<0)&(Y>0)/L=atan(Y/X)+PI;/if(X<0)&
6、;(Y<0)/L=atan(Y/X)+PI;/if(X>0)&(Y<0)/L=2*PI-atan(Y/X);L=atan2(YX);H=sqrt(X*X+Y*Y)/cos(B)-N;intBd,Bf,Ld,Lf;doubleBm,Lm;B=180*B/PI;/B转化为度做单位Bd=B;Bf=(B-Bd)*60;Bm=(B-Bd)*60-Bf)*60;L=180*L/PI;/L转化为度做单位Ld=L;Lf=(L-Ld)*60;Lm=(L-Ld)*60-Lf)*60;cout<<"大地坐标系中纬度,经度,大地高(国际单位)分别为"«fixed«setprecision(6)«Bm«endl«Ld«""«Lf«""«fixe
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 水费征收费管理办法
- 员工月报销管理办法
- 武汉市重症管理办法
- 基金业登记管理办法
- 村墓地管理管理办法
- 外债转贷款管理办法
- 小作坊管理办法河北
- 市农业项目管理办法
- 无底薪业务管理办法
- 名师带高徒管理办法
- 甲状腺十项及临床意义
- 南京“石塘竹海”衲头庵寺重建初步思路
- 陕09N1 供暖工程标准图集
- 钢筋配筋全套表格
- 王睿-PKPD优化抗菌药物治疗
- GB/T 9254.2-2021信息技术设备、多媒体设备和接收机电磁兼容第2部分:抗扰度要求
- 中医治未病课件培训课件
- 2023年湖北农谷实业集团有限责任公司招聘笔试模拟试题及答案解析
- 急性心梗后机械并发症课件整理
- DB11-T1835-2021 给水排水管道工程施工技术规程高清最新版
- 药理学课件(全)
评论
0/150
提交评论