MATLAB点云处理--三维曲面构建

Halcom

1. function ysw
   
2. clc;clear all;close all
   
3. x=[0.4,0.4,0.4,0.4,0.4,...
   
4.     0.5,0.5,0.5,0.5,0.5,...
   
5.     0.6,0.6,0.6,0.6,0.6];
   
6. y=[0.15,0.2,0.28,0.35,0.4,...
   
7.     0.15,0.2,0.28,0.35,0.4,...
   
8.     0.15,0.2,0.28,0.35,0.4,];
   
9. z1=[566,566,566,907,1351,...
   
10.     472,472,472,907,1351,...
    
11.     377,377,377,907,1351];
    
12. z2=[810,810,810,2000,3000,...
    
13.     675,675,675,2000,3000,...
    
14.     540,540,540,2000,3000];
    
15. sf1=fit([x',y'],z1','poly23');
    
16. sf2=fit([x',y'],z2','poly23');
    
17. 
    
18. xx = min(x):0.01:max(x);
    
19. yy = min(y):0.01:max(y);
    
20. [xx,yy] = meshgrid(xx,yy);
    
21. zz1 = sf1(xx,yy);
    
22. zz2 = sf2(xx,yy);
    
23. % figure(1),hold on;
    
24. % plot3(xx,yy,zz1,'ro');
    
25. % plot3(xx,yy,zz2,'go');
    
26. % hold off
    
27. 
    
28. p(:,1) = [xx(:);xx(:)]/2;
    
29. p(:,2) = [yy(:);yy(:)]/2;
    
30. p(:,3) = [zz1(:);zz2(:)]/5000000;
    
31. [t]=MyCrust(p);
    
32. % plot3(p(:,1),p(:,2),p(:,3),'g.')
    
33. %% plot of the oyput triangulation
    
34. figure(2)
    
35. hold on
    
36. % axis equal
    
37. % trisurf(t,p(:,1),p(:,2),p(:,3),'facecolor','c','edgecolor','b')%plot della superficie trattata
    
38. trisurf(t,p(:,1)*2,p(:,2)*2,p(:,3)*5000000,'edgecolor','b')
    
39. view(-37.5,30);
    
40. % axis tight;

复制代码

具体的代码下载链接：

游客，如果您要查看本帖隐藏内容请回复

z863338313

谢谢楼主，学习学习，看看哪个对我有用

舞动青春911

谢谢楼主分享，正在学习中