- x1=1200:10:4000;
- >> y1=1200:10:3600;
- >> z=[1130 1250 1280 1230 1040 900 500 700
- ;1320 1450 1420 1400 1300 700 900 850
- ;1390 1500 1500 1400 900 1100 1060 950
- ;1500 1200 1100 1350 1450 1200 1150 1010
- ;1500 1200 1100 1550 1600 1550 1380 1070
- ;1500 1550 1600 1550 1600 1600 1600 1550
- ;1480 1500 1550 1510 1430 1300 1200 980
- ];
- >> x1=1200:10:4000;
- >> y1=1200:10:3600;
- >> [x1,y1]=meshgrid(x1,y1);
- >> subplot(4,2,1);
- >> mesh(x,y,z)
- >> title('原始数据地貌图');
- >> subplot(4,2,2);
- >> contour(x1,y1,z1,20)
- >> title('原始数据等高线');
- >> subplot(4,2,3);
- >> z1= interp2(x,y,z,x1,y1,'nearest');
- >> mesh(x1,y1,z1)
- >> title('最邻近值地貌图');
- >> subplot(4,2,4);
- >> contour(x1,y1,z1,20)
- >> title('最邻近值等高线')
- >> subplot(4,2,5);
- >> z1= interp2(x,y,z,x1,y1,'linear');
- >> mesh(x1,y1,z1)
- >> title('双线性插值地貌图');
- >> subplot(4,2,6);
- >> contour(x1,y1,z1,20)
- >> title('双线性插值等高线');
- >> subplot(4,2,7);
- >> z1= interp2(x,y,z,x1,y1,'cubic');
- >> mesh(x1,y1,z1)
- >> title('双三次插值地貌图');
- >> subplot(4,2,8);
- >> contour(x1,y1,z1,20)
- >> title('双三次插值等高线');
- >>
复制代码
|