在一些题目中,给定目标点的经纬度,需要通过算法将其转换成平面坐标,以便更精确地计算距离。 使用墨卡托投影法将经纬度坐标投影为平面坐标。
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
参考文章:墨卡托投影法将经纬度转换为平面坐标