matlab军事沙盘问题

function main()
    % 提供的原始数据
    x = 0:200:1800;
    y = x;
    z = [2000,2000,2001,1992,1954,1938,1972,1995,1999,1999;
         2000,2002,2006,1908,1533,1381,1728,1959,1998,2000;
         2000,2005,2043,1921, 977, 897,1310,1930,2003,2000; 
         1997,1978,2009,2463,2374,1445,1931,2209,2050,2003;
         1992,1892,1566,1971,2768,2111,2653,2610,2121,2007;
         1991,1875,1511,1556,2221,1986,2660,2601,2119,2007;
         1996,1950,1797,2057,2849,2798,2608,2303,2052,2003;
         1999,1999,2079,2685,3390,3384,2781,2165,2016,2000;
         2000,2002,2043,2271,2668,2668,2277,2049,2003,2000;
         2000,2000,2004,2027,2067,2067,2027,2004,2000,2000];

    % 创建带色卡的等高线图和地貌图
    createContourMap(x, y, z);
    createSurfaceMap(x, y, z);

    % 图旋转以合理视角进行展示
    rotateMap(x, y, z, 2100);

    % 分析合理的区域抢占
    analyzeOptimalArea(x, y, z);
end

function createContourMap(x, y, z)
    % 最邻近插值
    x1 = 0:100:1800;
    y1 = x1';
    z1 = interp2(x, y, z, x1, y1, 'nearest');

    figure;
    colormap(mycolor);
    contourf(x1, y1, z1, 10);
    colorbar;
    title('Nearest Interpolation Contour Map');

    % 双线性插值
    x2 = 0:50:1800;
    y2 = x2';
    z2 = interp2(x, y, z, x2, y2, 'linear');

    figure;
    colormap(mycolor);
    contourf(x2, y2, z2, 10);
    colorbar;
    title('Bilinear Interpolation Contour Map');

    % 双三次插值
    x3 = 0:50:1800;
    y3 = x3';
    z3 = interp2(x, y, z, x3, y3, 'spline');

    figure;
    colormap(mycolor);
    contourf(x3, y3, z3, 10);
    colorbar;
    title('Bicubic Interpolation Contour Map');
end

function createSurfaceMap(x, y, z)
    figure;
    colormap(mycolor);
    surf(x, y, z);
    colorbar;
    title('Surface Map');
end

function rotateMap(x, y, z, elevation)
    % 找到海拔高度大于2100的山区
    [row, col] = find(z > elevation);

    % 计算山区的中心点
    center_x = mean(x(col));
    center_y = mean(y(row));

    % 计算旋转角度
    angle = atan2(center_y, center_x) * 180 / pi;

    % 旋转角度
    rotated_z = imrotate(z, angle);

    figure;
    colormap(mycolor);
    surf(rotated_z);
    colorbar;
    title('Rotated Map');
end

function analyzeOptimalArea(x, y, z)
    % 找到海拔最高的点
    [max_height, index] = max(z(:));
    [max_x, max_y] = ind2sub(size(z), index);

    % 绘制最高点及其周围区域
    figure;
    colormap(mycolor);
    surf(z);
    hold on;
    scatter3(max_x, max_y, max_height, 'r', 'filled');
    hold off;
    colorbar;
    title('Highest Point and Surrounding Area');
end