前言
前段时间做实验的时候用到经纬度坐标转笛卡尔坐标,在CSDN找了一下并整理成转换函数,现在找不到文章出处了,如有知道的麻烦告诉我,我好注明出处,谢谢。
一、转换公式笔记
为方便编公式,把整个转换公式拆分成几个部分:A、BB、C、D。
公式如下:
二、度数转弧度函数
例:180°转弧度制
rad=deg2rad(180)
三、坐标系转换函数
新建脚本,复制代码,保存文件命名为“LL3XY”至matlab搜索路径下。
% X:水平直角坐标,单位为千米(km);
% Y:纵向直角坐标,单位为千米(km);
% B:纬度,单位为弧度(rad);
% L:经度,单位为弧度(rad);
% Bo:投影基准纬度,Bo =0,单位为弧度((rad);
% Lo:坐标原点的经度,Lo =0,单位为弧度(rad);
% a:地球椭球体长半轴,a=6378137.0000,单位为米(m);
% b:地球椭球体短半轴,b=6356752.3142,单位为米(m);
% e:第一偏心率;
% e’:第二偏心率。
function [X,Y]=LL3XY(B,L)
a=6378137.0000;
b=6356752.3142;
e=0.00669438444204258280946884516695;
e1=0.00673950125438707813648942896063;
Lo =0;
Bo =0;
% 输入B、L
A=tan(pi/4+B/2);
BB=(1-e*sin(B))/(1+e*sin(B));
C=a^2/b;
D=sqrt(1+e1*e1*cos(Bo)*cos(Bo));
K=C*cos(Bo)/D;
Y=K*log(A*BB^(e/2))/1000;
X=K*(L-Lo)/1000;
end
四、总结
整体转换步骤:
(若经纬度的度分秒→度,可用dms2degrees函数;matlab自带,可自行查询用法)
1、把经纬度转化为以°为单位的小数;
2、用deg2rad函数转化为弧度制;
3、以弧度制的经纬度输入LL3XY转换函数中,输出以赤道和本初子午线交点为原点的XY坐标。
分享代码若有错漏之处,请大家留言批评指正!!
如果本篇文章对您有用的话,欢迎点赞收藏噢,谢谢谢谢,哈哈哈哈哈哈!!
主页还有更加丰富的内容噢 O(∩_∩)O :
分享代码若有错漏之处,请大家留言批评指正!!
如果本篇文章对您有用的话,欢迎点赞收藏噢,谢谢谢谢,哈哈哈哈哈哈!!
主页还有更加丰富的内容噢 O(∩_∩)O :
Matlab 土法求航海DCPA和TCPA,并根据DCPA正负判断目标船过本船船首or船尾
基于Matlab雷达视窗的来船运动矢量绘制
Matlab 沿着曲线的动态图制作
Matlab改进埃尔米特(Hermite)分段三次插值——(可在pchip函数中自定义导数值)
Matlab 四元素船舶领域代码复现