标题: matla山区地貌绘制问题 [打印本页]

---

作者: 嘿嘿嘿223    时间: 2021-4-9 19:36
标题: matla山区地貌绘制问题

1. x1=1200:10:4000;
   
2. >> y1=1200:10:3600;
   
3. >> z=[1130 1250 1280 1230 1040 900 500 700
   
4. ;1320 1450 1420 1400 1300 700 900 850
   
5. ;1390 1500 1500 1400 900 1100 1060 950
   
6. ;1500 1200 1100 1350 1450 1200 1150 1010
   
7. ;1500 1200 1100 1550 1600 1550 1380 1070
   
8. ;1500 1550 1600 1550 1600 1600 1600 1550
   
9. ;1480 1500 1550 1510 1430 1300 1200 980
   
10. ];
    
11. >> x1=1200:10:4000;
    
12. >> y1=1200:10:3600;
    
13. >> [x1,y1]=meshgrid(x1,y1);
    
14. >> subplot(4,2,1);
    
15. >> mesh(x,y,z)
    
16. >> title('原始数据地貌图');
    
17. >> subplot(4,2,2);
    
18. >> contour(x1,y1,z1,20)
    
19. >> title('原始数据等高线');
    
20. >> subplot(4,2,3);
    
21. >> z1= interp2(x,y,z,x1,y1,'nearest');
    
22. >> mesh(x1,y1,z1)
    
23. >> title('最邻近值地貌图');
    
24. >> subplot(4,2,4);
    
25. >> contour(x1,y1,z1,20)
    
26. >> title('最邻近值等高线')
    
27. >> subplot(4,2,5);
    
28. >> z1= interp2(x,y,z,x1,y1,'linear');
    
29. >>  mesh(x1,y1,z1)
    
30. >> title('双线性插值地貌图');
    
31. >> subplot(4,2,6);
    
32. >> contour(x1,y1,z1,20)
    
33. >> title('双线性插值等高线');
    
34. >> subplot(4,2,7);
    
35. >> z1= interp2(x,y,z,x1,y1,'cubic');
    
36. >> mesh(x1,y1,z1)
    
37. >> title('双三次插值地貌图');
    
38. >> subplot(4,2,8);
    
39. >> contour(x1,y1,z1,20)
    
40. >> title('双三次插值等高线');
    
41. >>
    

复制代码

---

欢迎光临 (http://www.51hei.com/bbs/)

Powered by Discuz! X3.1