web墨卡托与墨卡托经纬度之间的转换

在相同地理位置前提下，web墨卡托投影下的经纬度和墨卡托投影下的经纬度不一样，有没有转换的公式，用MATLAB演示

经纬度与网络墨卡托(Web Mercator)投影坐标的转换_经纬度转 webmector_顽劣的石头的博客-CSDN博客本文用JavaScript实现了经纬度与Web Mercator投影坐标之间的转换..._经纬度转 webmector

https://blog.csdn.net/shaxiaozilove/article/details/54908569

我想要的是经纬度转经纬度

墨卡托投影与经纬度之间的转换公式和MATLAB演示代码如下：

转换公式：墨卡托投影经纬度转经纬度坐标系的公式如下： latitude = (y / 20037508.34) * 180; longitude = (x / 20037508.34) * 180;

经纬度坐标系经纬度转墨卡托投影的公式如下： x = longitude * 20037508.34 / 180; y = log(tan((90 + latitude) * pi / 360)) / (pi / 180); y = y * 20037508.34 / 180;

MATLAB演示代码：

% 墨卡托投影经纬度转经纬度坐标系的公式
latitude = @(y) (y / 20037508.34) * 180;
longitude = @(x) (x / 20037508.34) * 180;

% 经纬度坐标系经纬度转墨卡托投影的公式
x = @(longitude) longitude * 20037508.34 / 180;
y = @(latitude) log(tan((90 + latitude) * pi / 360)) / (pi / 180);
y = @(latitude) y(latitude) * 20037508.34 / 180;

% 示例转换
x_mercator = 5000000;
y_mercator = 10000000;

latitude_deg = latitude(y_mercator);
longitude_deg = longitude(x_mercator);

fprintf('墨卡托投影经纬度：x = %.2f, y = %.2f\n', x_mercator, y_mercator);
fprintf('经纬度坐标系经纬度：latitude = %.10f, longitude = %.10f\n', latitude_deg, longitude_deg);

以上代码可以直接运行，输入墨卡托投影的x和y坐标即可获得经纬度坐标系下的经纬度值。

请注意，这里的公式是基于中国地区的计算方法，在其他区域可能会有所不同。