標(biāo)題: matla山區(qū)地貌繪制問題 [打印本頁]
作者: 嘿嘿嘿223    時(shí)間: 2021-4-9 19:36
標(biāo)題: matla山區(qū)地貌繪制問題

  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('原始數(shù)據(jù)地貌圖');
  17. >> subplot(4,2,2);
  18. >> contour(x1,y1,z1,20)
  19. >> title('原始數(shù)據(jù)等高線');
  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. >>
復(fù)制代碼






歡迎光臨 (http://www.torrancerestoration.com/bbs/) Powered by Discuz! X3.1