资源描述：
《基于matlab的大地坐标与直角坐标间的转换.docx》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯最新料推荐⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯测量程序设计实验报告实验名称：大地坐标与空间直角坐标的换算1⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯最新料推荐⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯实验四大地坐标与空间直角坐标的换算一、实验目的编写大地坐标与空间直角坐标相互转换的程序，并对格式化文件数据进行计算，验证程序。二、实验内容：1、大地坐标向空间直角坐标换算转换公式：x(Nh)cosBcosLy(Nh)cosBsinL（1）z[N(1e2)h]sinB其中：L为经度，B
2、为纬度，h为大地高，Na为卯酉圈曲率半径，1e2sin2Bea2b2为第一偏心率，a为旋转椭球长半轴，b为短半轴。aWGS84椭球参数：长半轴a=6378137扁率f=1/298.257223563根据上式创建以geo2xyz命名的函数，函数输入输出格式为[x,y,z]=geo2xyz(L,B,h)2、空间直角坐标向大地坐标换算根据式（1）推导大地坐标向空间直角坐标转换公式：Larctan(y/x)BzNe2sinB)arctan(x2y2hx2y2NcosB注意计算纬度时需要用到迭代，可用Baz)作为初始值。arctan(b
3、x2y2创建以xyz2geo命名的函数，函数输入输出格式为[L,B,h]=xyz2geo(x,y,z)2⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯最新料推荐⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯三、实验步骤1、大地坐标向空间直角坐标换算主程序：%%大地坐标向空间直角坐标换算%函数的输入输出格式为[x,y,z]=geo2xyz(L,B,h)[filename,pathname]=uigetfile('*.txt','请选择打开的数据文件');file=[pathname,filename];data=importdata(f
4、ile);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');fprintf(fid,'%15.7f%15.7f%15.7f
5、',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-
6、e^2)+h).*sin(B);endfunctionrad=dms2rad(jiaodu)%度分秒－>弧度（rad）degree=fix(jiaodu);mimute=fix((jiaodu-degree)*100);3⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯最新料推荐⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯second=(jiaodu-degree-mimute/100)*10000;degree=degree+mimute/60+second/3600;rad=degree/180*pi;end2、空间直角坐标向大地
7、坐标换算主程序：%%将文件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]=xyz2ge
8、o(x,y,z);[filename_out,pathname_out]=uiputfile('*.txt','请选择要输出数据文件');fileout=[pathname_out,filename_out];fid=fopen(fileout,'wt');fprintf(f
