数学建模暑期集训实战:经纬度转换为平面坐标

2022-06-14 11:47:59 浏览数 (1)

在一些题目中,给定目标点的经纬度,需要通过算法将其转换成平面坐标,以便更精确地计算距离。 使用墨卡托投影法将经纬度坐标投影为平面坐标。

matlab代码

代码语言:javascript复制
function [x,y]=ll_xy(lng, lat)
	earthRad = 6378137.0;
	x = ((lng .* pi) ./ 180) .* earthRad;
	a = (lat .* pi) ./ 180;
	y = (earthRad ./ 2) .* log((1.0   sin(a)) ./ (1.0 - sin(a)));
end
代码语言:javascript复制
tic
format long g
[x_p ,y_p] = ll_xy(x,y);
x_p = x_p - mean(x_p);
y_p = y_p - mean(y_p);
toc

参考文章:墨卡托投影法将经纬度转换为平面坐标

0 人点赞