Matlab经纬度坐标转换xy坐标,经纬度坐标系转换为UTM坐标系(matlab)
如题所说,直接上程序。验证自己做一下,结果应该是对的。诚不我欺。
(注意:程序名与函数名保持一致!!!)
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
Matlab经纬度坐标转换xy坐标,经纬度坐标系转换为UTM坐标系(matlab)相关推荐
- 经纬度坐标系转换为UTM坐标系(matlab)
如题所说,直接上程序.验证自己做一下,结果应该是对的.诚不我欺. (注意:程序名与函数名保持一致!!!) function[shuchu]=lat_lon2utm(lat_shuru,lon_shur ...
- 经纬度坐标转换xy坐标 python_在Python中使用NewtonRaphson迭代将经纬度转换为xy Mollweide地图坐标...
我试图编写一个程序,从用户那里获取一组经度和纬度坐标,将它们转换为Mollweide投影图的x&y坐标,然后报告这些坐标处的像素值(在本例中,是噪声温度).在 我使用的地图/数据是Haslam ...
- 3DMax坐标系转换为Unity坐标系
3DMax坐标系转换为Unity坐标系 3DMax是右坐标系,并且Z轴向上 Unity 是左手坐标系,并且Y轴向上 在导出FBX文件的时候3DMax是可以指定Y,Z轴向上的 问题是unity里面X轴会 ...
- 经纬度转换XY坐标,并计算距离
import pandas as pd import numpy as np from pyproj import Transformer import matplotlib.pyplot as pl ...
- Python使用proj将GPS经纬度数据转换为utm坐标系(真实的地理位置坐标系,以米为单位),然后可以进行轨迹显示
1. 代码如下 from pyproj import Transformer import pyproj import json from matplotlib.pyplot import *if _ ...
- MATLAB极坐标与xy坐标互相转换_不改变数据形状_极坐标变量v_p(theta,r)的平面图
不改变数据形状的极坐标转换这应是原创,网上能找到的都是会改变数据形状的,比如正方形的xy图变为圆形的极坐标图. 而我的程序则是直接投影,即正方形的xy图变为正方形的数据区块投影在极坐标图上. 一共四个 ...
- 经纬度坐标转换成px_经纬度坐标转像素坐标
"```\n// 经纬度坐标转像素坐标\nfunction latlng2px(lat, lng) {\n var LatLng = (function(){\n function LatL ...
- GPS经纬度换算成XY坐标
/// <summary> /// GPS经纬度换算成x,y坐标 /// </summary> /// <para ...
- oracle经纬度换算成xy坐标,谁有全套经纬度转换成XY坐标的好的软件?
用的实施么地图?标准地形图吗? 如果是标准地形图,还是有办法转化的,因为用的是高斯-克吕格投影 该投影按照投影带中央子午线投影为直线且长度不变和赤道投影为直线的条件,确定函数的形式,从而得到高斯一克吕 ...
- 双纽线通过matlab绘图输出xy坐标,Matlab 6 (Advanced_Plotting)
Polar Chart(极坐标图): Snip20180214_15.png 例:蜗牛线(蚌线.蚶线) theta = linspace(0, 2*pi); r = 1 + 2*cos(theta); ...
最新文章
- Android studio 另一个程序正在使用此文件,进程无法访问
- php while循环控制实例讲解
- binaryreader java,C# 二进制文件的读写 | 菜鸟教程
- kubernetes问题排查
- Nginx的nginx.conf配置文件中文注释说明
- mysql shell 1.0.10_MySQL Shell(使用Shell命令管理MySQL)下载 v1.0.10 官方32位+64位Windows版 - 比克尔下载...
- Shell入门(三)之字符串
- 《四世同堂》金句摘抄(九)
- Dungeon Master——BFS
- asp.net session 如何知道是哪个浏览器客户端_小弟该如何复习 Java?
- 第三篇: UpdatePanel 控件示例
- 在ASP.NET中如何运行后台任务
- jsf服务_JSF dataTable示例
- 声卡loopback有什么用_萌新做音乐那点事 | 外置专业声卡的选择方法与推荐
- 12.Linux 高性能服务器编程 --- 高性能 IO 框架库 Libevent
- GB28181协议RTP传输
- 在计算机系统中对文件执行彻底删除,Windows系统中,教你彻底删除C盘的顽固文件,瞬间多出10个G!...
- 解决file_get_contents乱码
- python时间函数纳秒_python 时间 纳秒_Pandas时间序列(翻译)
- bert-textcnn实现多标签文本分类(基于keras+keras-bert构建)