电脑桌面
添加小米粒文库到电脑桌面
安装后可以在桌面快捷访问

基于matlab的大地坐标与直角坐标间的转换VIP免费

基于matlab的大地坐标与直角坐标间的转换_第1页
1/5
基于matlab的大地坐标与直角坐标间的转换_第2页
2/5
基于matlab的大地坐标与直角坐标间的转换_第3页
3/5
测量程序设计实验报告实验名称:大地坐标与空间直角坐标的换算实验四大地坐标与空间直角坐标的换算一、实验目的编写大地坐标与空间直角坐标相互转换的程序,并对格式化文件数据进行计算,验证程序。二、实验内容:1、大地坐标向空间直角坐标换算转换公式:x=(N+h)cosBcosLy=(N+h)cosBsinLz=[N(1−e2)+h]sinB(1)其中:L为经度,B为纬度,h为大地高,N=a√1−e2sin2B为卯酉圈曲率半径,e=√a2−b2a为第一偏心率,a为旋转椭球长半轴,b为短半轴。WGS84椭球参数:长半轴a=6378137扁率f=1/298.257223563根据上式创建以geo2xyz命名的函数,函数输入输出格式为[x,y,z]=geo2xyz(L,B,h)2、空间直角坐标向大地坐标换算根据式(1)推导大地坐标向空间直角坐标转换公式:L=arctan(y/x)B=arctan(z+Ne2sinB√x2+y2)h=√x2+y2cosB−N注意计算纬度时需要用到迭代,可用B=arctan(azb√x2+y2)作为初始值。创建以xyz2geo命名的函数,函数输入输出格式为[L,B,h]=xyz2geo(x,y,z)三、实验步骤1、大地坐标向空间直角坐标换算主程序:%%大地坐标向空间直角坐标换算%函数的输入输出格式为[x,y,z]=geo2xyz(L,B,h)[filename,pathname]=uigetfile('*.txt','请选择打开的数据文件');file=[pathname,filename];data=importdata(file);L=data.data(:,1);B=data.data(:,2);h=data.data(:,3);[x,y,z]=geo2xyz(L,B,h);A=[x,y,z];A=A';[filename_out,pathname_out]=uiputfile('*.txt','请选择要输出数据文件');fileout=[pathname_out,filename_out];fid=fopen(fileout,'wt');fprintf(fid,'xyz\n');fprintf(fid,'%15.7f%15.7f%15.7f\n',A);close('all');函数:function[x,y,z]=geo2xyz(L,B,h)%大地坐标经纬度转换成空间直角坐标B=dms2rad(B);L=dms2rad(L);a=6378137;%a是长半轴f=1/298.257223563;%f是扁率b=a-a*f;e=sqrt(a^2-b^2)/a;N=a./(sqrt(1-e^2.*(sin(B)).^2));%N为卯酉圈半径率,e为第一偏心率x=(N+h).*cos(B).*cos(L);y=(N+h).*cos(B).*sin(L);z=(N*(1-e^2)+h).*sin(B);endfunctionrad=dms2rad(jiaodu)%度分秒->弧度(rad)degree=fix(jiaodu);mimute=fix((jiaodu-degree)*100);second=(jiaodu-degree-mimute/100)*10000;degree=degree+mimute/60+second/3600;rad=degree/180*pi;end2、空间直角坐标向大地坐标换算主程序:%%将文件data.2.txt中的空间直角坐标系转换为大地坐标,并将计算结果按照格式存储在文件data3.txt中%data3.txt格式为:经度(ddmmss)纬度(ddmmss)大地高[filename,pathname]=uigetfile('*.txt','请选择打开的数据文件');file=[pathname,filename];data=importdata(file);x=data.data(:,1);y=data.data(:,2);z=data.data(:,3);[L,B,H]=xyz2geo(x,y,z);[filename_out,pathname_out]=uiputfile('*.txt','请选择要输出数据文件');fileout=[pathname_out,filename_out];fid=fopen(fileout,'wt');fprintf(fid,'经度(ddmmss)纬度(ddmmss)大地高\n');fprintf(fid,'%10.7f%10.7f%7.3f\n',[L,B,H]');fclose('all');函数:function[L,B,H]=xyz2geo(x,y,z)%将直角坐标转换为大地坐标%%已知:WGS-84椭球参数f=1/298.257223563;%扁率f=(a-b)/aa=6378137;%长半轴b=a*(1-f);%短半轴e=(sqrt(a^2-b^2))/a;%第一偏心率%%经度L的计算L=atan2(y,x);L=rad2dms(L);%%纬度B的计算B0=atan2((a*z),(b.*sqrt(x.^2+y.^2)));%B初始值while1N=a./(sqrt(1-(e^2).*(sin(B0).^2)));%卯酉圈曲率半径NB=atan2(z+N.*(e^2).*sin(B0),(sqrt(x.^2+y.^2)));ifabs(B0-B)<10^-6break;endabs(B0-B)B0=B;end%%大地高H的计算H=(sqrt(x.^2+y.^2))./cos(B)-N;B=degree2dms(B.*180/pi);%纬度BEndfunctiondms=degree2dms(jiaodu)%度->度分秒(dd.mmss)degree=fix(jiaodu);mimute=fix((jiaodu-degree)*60);second=((jiaodu-degree)*60-mimute)*60;dms=degree+mimute/100+second/10000;

1、当您付费下载文档后,您只拥有了使用权限,并不意味着购买了版权,文档只能用于自身使用,不得用于其他商业用途(如 [转卖]进行直接盈利或[编辑后售卖]进行间接盈利)。
2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。
3、如文档内容存在违规,或者侵犯商业秘密、侵犯著作权等,请点击“违规举报”。

碎片内容

基于matlab的大地坐标与直角坐标间的转换

您可能关注的文档

确认删除?
VIP
微信客服
  • 扫码咨询
会员Q群
  • 会员专属群点击这里加入QQ群
客服邮箱
回到顶部