


下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
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年小自考公共事业管理保险政策解析试题及答案
- 金属切削原理与刀具知到智慧树章节测试课后答案2024年秋潍坊工程职业学院
- 行政管理语言艺术试题及答案
- 2024小自考汉语言写作提升试题及答案
- 2025年山东省事业单位初级综合类岗位公开招聘平台笔试历年典型考题及考点剖析附带答案详解
- 小自考汉语言文学常见问题与试题及答案
- (2024年秋季版)七年级历史下册 第5单元 隋唐时期 第3课《盛唐的社会景象》教学设计1 川教版
- 酒驾理论考试题及答案
- 小自考行政管理基础复习试题及答案
- 面试题及答案100例
- 成语故事杞人忧天PPT教案
- 部编版三年级上册音乐知识点汇总
- 生命体征的测量PPT幻灯片课件
- 吉林省吉林市高考报名登记表
- 质量保证体系结构图(共3页)
- 天然气长输管道的腐蚀与防护措施
- IEC60335-1(中文)
- 排沙泵检修工艺流程及验收标准
- 常见职业危害相应职业禁忌证(简表)
- Q-BQB_402-2009冷连轧碳素钢板及钢带
- 内蒙古自治区实施少数民族高层次骨干人才计划暂行办法内蒙古教育
评论
0/150
提交评论