数学建模暑期集训18：经纬度转换为平面坐标

Posted 2021-09-07 Z|Star

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

matlab代码

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

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

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