下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
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. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2024学年泰安市高三语文(上)12月考试卷附答案解析
- 《电工电子技术 》课件-第10章 逻辑代数基础与组合逻辑电路
- 2025年白城货运从业资格证考试题库
- 分析风险的报告范文
- 2025年吉林市货运从业资格证模拟考试题库下载
- 2025年滨州货运资格证模拟考试题库
- 2025年衡水货物从业资格证考试题
- 《摩擦力受力分析》课件
- 《设计暴雨》课件
- 2025企业融资典当借款合同
- 麓湖营销体系及逻辑
- 九年级历史上册 第19课《巴黎公社》导学案 中华书局版-中华书局版初中九年级上册历史学案
- 《9加几》评课稿
- CTCS列控系统及车载设备介绍
- 某某单位关于开展谈心谈话活动的情况报告情况统计五篇范文
- 无线铁塔及天馈线安装专项施工方案
- 气动夯管技术在管道施工中的应用
- ARAMCO阿美认证检验员考试题及答案(共56页)
- 仪器自检自效校验记录 2
- 浅谈脱口秀的语言幽默感——以《脱口秀大会第二季》为例
- 聚合物改性教案(1-2)(课堂PPT)
评论
0/150
提交评论