% 溢油地点的经纬度坐标(示例数据,替换为实际数据)
oilSpillLocations = [
-80.1234, 25.6789; % 经度, 纬度
-75.4321, 30.9876;
% 添加更多溢油地点的经纬度坐标...
];
% 使用米勒投影函数转换经纬度坐标
[x, y] = millerProjection(oilSpillLocations(:, 1), oilSpillLocations(:, 2));
% 绘制地图并标记溢油地点
figure;
plot(x, y, 'ro'); % 以红色圆圈标记溢油地点
title('Oil Spill Locations in Miller Projection');
% 定义米勒投影函数
function [x, y] = millerProjection(lon, lat)
R = 6371000; % 地球半径,单位:米
% 将经纬度转换为弧度
lonRad = deg2rad(lon);
latRad = deg2rad(lat);
% 米勒投影的转换公式
x = lonRad;
y = log(tan(pi/4 + latRad/2)) * R;
end
这段代码将在米勒投影下绘制地图,并以红色圆圈标记溢油地点的位置。