标签:phi lamda0 lon UTM shuru 180 matlab lat 坐标系
如题所说,直接上程序。验证自己做一下,结果应该是对的。诚不我欺。
(注意:程序名与函数名保持一致!!!)
function[shuchu]=lat_lon2utm(lat_shuru,lon_shuru)
%地理经纬度坐标转换为UTM坐标
size_shuzu=size(lat_shuru);
for i=1:size_shuzu(2)
%输入经纬度
% lat=29.819206;
% lon=116.133243;
lat=lat_shuru(i);
lon=lon_shuru(i);
%% %%unit:km
a=6378.137;
e=0.0818192;
k0=0.9996;
E0=500;
N0=0;
Zonenum=fix(lon/6)+31;
lamda0=(Zonenum-1)*6-180+3;%degree
lamda0=lamda0*pi/180;%radian
phi=lat*pi/180;
lamda=lon*pi/180;%radian
v=1/sqrt(1-e^2*(sin(phi)^2));
A=(lamda-lamda0)*cos(phi);
T=tan(phi)*tan(phi);
C=e^2*cos(phi)*cos(phi)/(1-e^2);
s=(1-e^2/4-3*e^4/64-5*e^6/256)*phi-(3*e^2/8+3*e^4/32+45*e^6/1024)*sin(2*phi)+(15*e^4/256+45*e^6/1024)*sin(4*phi)-35*e^6/3072*sin(6*phi);
UTME=E0+k0*a*v*(A+(1-T+C)*A^3/6+(5-18*T+T^2)*A^5/120);
UTMN=N0+k0*a*(s+v*tan(phi)*(A^2/2+(5-T+9*C+4*C^2)*A^4/24+(61-58*T+T^2)*A^6/720));
shuchu(1,i)=UTME;
shuchu(2,i)=UTMN;
end
标签:phi,lamda0,lon,UTM,shuru,180,matlab,lat,坐标系 来源: https://blog.csdn.net/weixin_43955546/article/details/113408099
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。