1.png (43.45 KB, 下載次數: 59)
下載附件
2021-4-9 20:23 上傳
- 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('雙三次插值等高線');
- >>
復制代碼
|