[求助]：各位可否幫幫忙.會不會把花瓶變成3D.還是運行不了..這個星期要交了.

kaiyim

                               
登录/注册后可看大图

各位可否幫幫忙........會不會把上面那個花瓶變成3D的MODE....我已經做了很久.....還是運行不了,可以可幫幫忙, 看看我的代码有沒有錯....謝謝謝謝謝謝

2. [FONT =courier new]
   
3. function vase;
   
4. %/*this function generate a vector model of the profile of the bottle*/
   
5. close all;
   
6. 
   
7. %read in the image and turn to gery scale image
   
8. I = rgb2gray(imread('vase.jpeg'));
   
9. imshow(I);
   
10. 
    
11. %adjust the contrast to avoid edge the showdaw of the bottle
    
12. I = imadjust(I,[],[], 0.15);
    
13. figure, imshow(I)
    
14. 
    
15. %set up the best threshold
    
16. level = 0.1
    
17. BW1 = edge(I,'canny', level);
    
18. figure, imshow(BW1);
    
19. 
    
20. % remove spur pixels
    
21. BW2 = bwmorph(BW1, 'spur', 2);
    
22. figure, imshow(BW2), title('bwmorph--spur')
    
23. 
    
24. %using the closing technique to fill in the bottle
    
25. BW3 = bwmorph(BW2, 'dilate', 40);
    
26. figure, imshow(BW3), title('bwmorph--dilate')
    
27. BW4 = bwmorph(BW3, 'erode', 40);
    
28. figure, imshow(BW4), title('bwmorph--erode')
    
29. 
    
30. %remove all the piexls expect the boundary of the bottle
    
31. BW = bwmorph(BW4, 'remove', inf);
    
32. figure, imshow(BW), title('bwmorph--remove')
    
33. 
    
34. 
    
35. 
    
36. %rotate the image by 180 to let the sweep generator working properly
    
37. BW = rot90(rot90(BW));
    
38. figure, imshow(BW);
    
39. 
    
40. %find out the centroid of the bottle
    
41. theCentroid = regionprops(BW, 'Centroid');
    
42. theCentroid = round([theCentroid.Centroid]);
    
43. 
    
44. %find all the nodes of the whole bottle
    
45. [Y X] = find(BW);
    
46. 
    
47. % %select the peixels of the left half bottle
    
48. bufferIndex = find(X >=theCentroid(1));
    
49. X = X(bufferIndex);
    
50. Y = Y(bufferIndex);
    
51. 
    
52. %adjust the coordinates
    
53. %subtracting the min.vaule of X element to set the orgin start from 0
    
54. X = abs(X - min(X));
    
55. %adjust the position of the bottle to touch the ground
    
56. Y = Y - min(Y);
    
57. 
    
58. %Sort elements of X Y in ascending order of Y
    
59. [bufferY, bufferIndex] = sort(Y);
    
60. X = X(bufferIndex);
    
61. Y = Y(bufferIndex);
    
62. 
    
63. %reduce the unncessary piexls
    
64. numOfCoordinate = size(X,1);
    
65. lastIndex = numOfCoordinate;
    
66. currIndex = 1;
    
67. while (currIndex < lastIndex)
    
68.     currIndex2 = 1;
    
69.     while (currIndex2 < lastIndex)
    
70.         %remove the pixel if the X or Y is same as the neighborhood
    
71.         %as it means it is not the essential pixel to describe the bottle
    
72.         if (X(currIndex2)==X(currIndex2+1)) | (Y(currIndex2)==Y(currIndex2+1))
    
73.             %remove the pixel
    
74.             X(currIndex2) = [];   
    
75.             Y(currIndex2) = [];
    
76.             lastIndex = lastIndex-1;
    
77.         else
    
78.             currIndex2 = currIndex2 +1;
    
79.         end
    
80.     end
    
81.     lastIndex = lastIndex-currIndex;
    
82. end   
    
83. %repair the hole of the bottom of the bottle
    
84. X = [X; 0];
    
85. Y = [Y; min(Y)+1];
    
86. 
    
87. numOfCoordinate = size(X,1);
    
88. 
    
89. %call part3 function to sweep generate a 3D bottle
    
90. sweep(X,Y,numOfCoordinate);
    
91. 
    
92. function sweep(X,Y,numOfCoordinate);   
    
93. %/*this function sweeps generate a 3D bottle*/
    
94. 
    
95. %set the num of nodes for the sweep generator
    
96. n=50;
    
97. %calculate the number of nodes and patches
    
98. numnodes = numOfCoordinate*n;
    
99. numpatch = (numOfCoordinate-1)*n;
    
100. 
     
101. %calculate the node values
     
102. for i = 0:n-1
     
103.     for j = 1:numOfCoordinate
     
104.         theta = 2*i*pi/n;
     
105.         Node(j+numOfCoordinate*i,1) = X(j)*cos(theta);
     
106.         Node(j+numOfCoordinate*i,2) = X(j)*sin(theta);
     
107.         Node(j+numOfCoordinate*i,3) = Y(j);
     
108.     end
     
109. end
     
110. 
     
111. %assign the nodes to the patches
     
112. for i = 1:n
     
113.     term = numOfCoordinate*(i-1);
     
114.     for j = 1:(numOfCoordinate-1)
     
115.         FaceNode(((numOfCoordinate-1)*(i-1))+j,:) = [(term+j) (term+j+1) (term+j+1+numOfCoordinate) (term+j+numOfCoordinate)];
     
116.     end
     
117. end
     
118. 
     
119. %connects the last profile to the first profile
     
120. for i = numpatch-(numOfCoordinate-2):numpatch
     
121.     FaceNode(i,3) = FaceNode(i,3)-numnodes;
     
122.     FaceNode(i,4) = FaceNode(i,4)-numnodes;
     
123. end
     
124. 
     
125. %displays the 3D bottle
     
126. figure, patch('Vertices',Node,'Faces',FaceNode,'facevertexCdata',[0.2 0.2 1],'facecolor','black', 'EdgeColor','none')
     
127. light('Position',[100 -100 0],'Style','infinite','color','white');
     
128. 
     
129. %the lighting effects
     
130. lighting phong
     
131. 
     
132. axis square;
     
133. axis([-170, 170, -170, 170, 0, 350]);
     
134. rotate3d on;
     
135. [/FONT]

[FONT =courier new]<br /> function vase;<br /> %/*this function generate a vector model of the profile of the bottle*/<br /> close all;<br /> <br /> %read in the image and turn to gery scale image<br /> I = rgb2gray(imread('vase.jpeg'));<br /> imshow(I);<br /> <br /> %adjust the contrast to avoid edge the showdaw of the bottle<br /> I = imadjust(I,[],[], 0.15);<br /> figure, imshow(I)<br /> <br /> %set up the best threshold<br /> level = 0.1<br /> BW1 = edge(I,'canny', level);<br /> figure, imshow(BW1);<br /> <br /> % remove spur pixels<br /> BW2 = bwmorph(BW1, 'spur', 2);<br /> figure, imshow(BW2), title('bwmorph--spur')<br /> <br /> %using the closing technique to fill in the bottle<br /> BW3 = bwmorph(BW2, 'dilate', 40);<br /> figure, imshow(BW3), title('bwmorph--dilate')<br /> BW4 = bwmorph(BW3, 'erode', 40);<br /> figure, imshow(BW4), title('bwmorph--erode')<br /> <br /> %remove all the piexls expect the boundary of the bottle<br /> BW = bwmorph(BW4, 'remove', inf);<br /> figure, imshow(BW), title('bwmorph--remove')<br /> <br /> <br /> <br /> %rotate the image by 180 to let the sweep generator working properly<br /> BW = rot90(rot90(BW));<br /> figure, imshow(BW);<br /> <br /> %find out the centroid of the bottle<br /> theCentroid = regionprops(BW, 'Centroid');<br /> theCentroid = round([theCentroid.Centroid]);<br /> <br /> %find all the nodes of the whole bottle<br /> [Y X] = find(BW);<br /> <br /> % %select the peixels of the left half bottle<br /> bufferIndex = find(X >=theCentroid(1));<br /> X = X(bufferIndex);<br /> Y = Y(bufferIndex);<br /> <br /> %adjust the coordinates<br /> %subtracting the min.vaule of X element to set the orgin start from 0<br /> X = abs(X - min(X));<br /> %adjust the position of the bottle to touch the ground<br /> Y = Y - min(Y);<br /> <br /> %Sort elements of X Y in ascending order of Y<br /> [bufferY, bufferIndex] = sort(Y);<br /> X = X(bufferIndex);<br /> Y = Y(bufferIndex);<br /> <br /> %reduce the unncessary piexls<br /> numOfCoordinate = size(X,1);<br /> lastIndex = numOfCoordinate;<br /> currIndex = 1;<br /> while (currIndex < lastIndex) <br />     currIndex2 = 1;<br />     while (currIndex2 < lastIndex)<br />         %remove the pixel if the X or Y is same as the neighborhood<br />         %as it means it is not the essential pixel to describe the bottle<br />         if (X(currIndex2)==X(currIndex2+1)) | (Y(currIndex2)==Y(currIndex2+1))<br />             %remove the pixel<br />             X(currIndex2) = [];    <br />             Y(currIndex2) = [];<br />             lastIndex = lastIndex-1;<br />         else<br />             currIndex2 = currIndex2 +1;<br />         end<br />     end<br />     lastIndex = lastIndex-currIndex;<br /> end    <br /> %repair the hole of the bottom of the bottle<br /> X = [X; 0];<br /> Y = [Y; min(Y)+1];<br /> <br /> numOfCoordinate = size(X,1);<br /> <br /> %call part3 function to sweep generate a 3D bottle<br /> sweep(X,Y,numOfCoordinate);<br /> <br /> function sweep(X,Y,numOfCoordinate);   <br /> %/*this function sweeps generate a 3D bottle*/<br /> <br /> %set the num of nodes for the sweep generator<br /> n=50;<br /> %calculate the number of nodes and patches<br /> numnodes = numOfCoordinate*n;<br /> numpatch = (numOfCoordinate-1)*n;<br /> <br /> %calculate the node values<br /> for i = 0:n-1<br />     for j = 1:numOfCoordinate<br />         theta = 2*i*pi/n;<br />         Node(j+numOfCoordinate*i,1) = X(j)*cos(theta); <br />         Node(j+numOfCoordinate*i,2) = X(j)*sin(theta);<br />         Node(j+numOfCoordinate*i,3) = Y(j);<br />     end<br /> end<br /> <br /> %assign the nodes to the patches<br /> for i = 1:n<br />     term = numOfCoordinate*(i-1);<br />     for j = 1:(numOfCoordinate-1)<br />         FaceNode(((numOfCoordinate-1)*(i-1))+j,:) = [(term+j) (term+j+1) (term+j+1+numOfCoordinate) (term+j+numOfCoordinate)];<br />     end<br /> end<br /> <br /> %connects the last profile to the first profile<br /> for i = numpatch-(numOfCoordinate-2):numpatch<br />     FaceNode(i,3) = FaceNode(i,3)-numnodes;<br />     FaceNode(i,4) = FaceNode(i,4)-numnodes;<br /> end<br /> <br /> %displays the 3D bottle<br /> figure, patch('Vertices',Node,'Faces',FaceNode,'facevertexCdata',[0.2 0.2 1],'facecolor','black', 'EdgeColor','none')<br /> light('Position',[100 -100 0],'Style','infinite','color','white');<br /> <br /> %the lighting effects<br /> lighting phong<br /> <br /> axis square;<br /> axis([-170, 170, -170, 170, 0, 350]);<br /> rotate3d on;<br /> [/FONT]