一、山区地貌:在某山区测得一些地点的高程如下表:(平面区域1200≤x ≤4000,1200≤y ≤3600),试作出该山区的地貌图和等高线图。
clc;clear;clf;
x = 1200:400:4000;
y = 1200:400:3600;
highth = [1480 1500 1550 1510 1430 1300 1200 980
1500 1550 1600 1550 1600 1600 1600 1550
1500 1200 1100 1550 1600 1550 1380 1070
1500 1200 1100 1350 1450 1200 1150 1010
1390 1500 1500 1400 900 1100 1060 950
1320 1450 1420 1400 1300 700 900 850
1130 1250 1280 1230 1040 900 500 700];
xi = 1200:10:4000;
yi = 1200:10:3600;
% 进行三次插值
zi = interp2(x,y,highth,xi',yi,'cubic');
% 绘制三维图和等高线图
mc = meshc(xi,yi,zi);
% 调整等高线图得属性
mc(2).Fill = 'on';
mc(2).LineColor = 'b';
mc(2).ShowText = 'on';
view([-155 17]);
title("山区地貌图");
% 单独绘制等高线图
[M,c] = contour(xi,yi,zi,8);
c.Fill = 'on';
c.LineColor = 'b';
c.ShowText = 'on';colorbar;
title("山区等高线图");