- 积分
- 2853
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-12-16
- 最后登录
- 1970-1-1
|
发表于 2022-7-11 15:07:10
|
显示全部楼层
回帖奖励 +3 金钱
unction [lat,lon] =lambert72toWGS84(x,y)
%UNTITLED 此处显示有关此函数的摘要
% 此处显示详细说明
n = 0.77164219;
F = 1.81329763;
thetaFudge = 0.00014204;
e = 0.08199189;
a = 6378388;
xDiff = 149910;
yDiff = 5400150;
theta0 = 0.07604294;
xReal = xDiff - x;
yReal = yDiff - y;
rho = sqrt(xReal * xReal + yReal * yReal);
theta =atan(xReal /(-yReal));
lon = (theta0 + (theta + thetaFudge)/n) * 180/pi;
lat=0;
for i=1:5
lat= 2*atan((F*a/rho).^(1/n)*((1+e*sin(lat))/(1-sin(lat))).^(e/2))-pi/2;
end
lat=lat*180/pi;
end
参考:https://www.bram.us/2012/03/14/convert-lambert72-to-wgs84/ |
|