用matlab写一段墨卡托投影的程序

在MATLAB中，可以使用基础地图工具箱（Mapping Toolbox）来创建墨卡托投影（Mercator Projection）。墨卡托投影是一种等角圆柱投影，广泛用于海洋和航空导航。下面是一个简单的示例，展示如何将经纬度坐标转换为墨卡托平面上的点：

% 导入Mapping Toolbox
if ~isToolboxInstalled(' Mapping')
    error('需要安装Mappping Toolbox才能运行此代码');
end

% 定义纬度和经度范围
lat = -90:10:90; % 纬度范围
lon = -180:10:180; % 经度范围

% 创建地球模型
R = earth; % 地球半径，近似值6371公里

% 转换到墨卡托投影
[merc_lat, merc_lon] = pol2cart(lon, lat, R); % 使用pol2cart函数进行转换

% 打印转换后的结果
disp(['纬度: ', num2str(merc_lat), ' 经度: ', num2str(merc_lon)])
% 或者你可以直接画出墨卡托投影的地图
figure;
geoplot(lat, lon, 'b.'); % 用蓝色线绘制原始经纬度
geoplot(merc_lat, merc_lon, 'r.'); % 用红色线绘制墨卡托投影的点
title('Mercator Projection');
xlabel('Mercator Longitude');
ylabel('Mercator Latitude');