


下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
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. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 预备役介绍课件
- 健康饮食产业园项目运营管理方案(范文模板)
- 城市污水管网建设工程初步设计(参考模板)
- 2025年交流调频调压牵引装置项目建议书
- 五年级上册音乐教案(全集)
- 现场施工事故应急救援预案制度
- 2025年新型节能水泵、风机和压缩机项目建议书
- 2025年防沉剂项目发展计划
- 2025年高纯工艺系统设备分条机项目建议书
- 针对2025年工业污染源全面达标排放计划的实施方案的环境风险评估与管理001
- 部编八下语文游记阅读训练题语文八年级下册能力训练(部编版)
- 保修管理控制程序
- GB/T 9117-2010带颈承插焊钢制管法兰
- GB/T 12513-2006镶玻璃构件耐火试验方法
- 人教版音乐三年级上册教材介绍-课件
- 装修改造工程施工总平面图6
- 教师的职业生涯规划与专业发展课件
- 生物安全自查表
- 广州小升初-学籍表打印版
- 天津市-解除劳动合同证明书
- 公司一年完税证明模板
评论
0/150
提交评论