Changeset 663

Timestamp:

02/10/10 22:02:57 (2 years ago)

Author:

taborts

Message:

fixed problem with all figures not appearing

Location:

branches/2010-image-rec

Files:

• hist_test.asv (modified) (2 diffs)
• hist_test.m (modified) (2 diffs)

branches/2010-image-rec/hist_test.asv

@@ -106 +106 @@
 mask = imclose(mask, se);
 
-textureactual=texture;
-for threshold = .3:.02:.5
-    texture(textureactual > threshold & mask == 0) = 1;
-    texture(texture ~= 1) = 0;
-    imtool(texture);
+line_mask = texture;
+for threshold = .42
+    line_mask(texture > threshold & mask == 0) = 1;
+    line_mask(texture <= threshold | mask ~= 0) = 0;
+    imtool(line_mask);
 end
 
@@ -124 +124 @@
 line([left_x right_x right_x left_x left_x], ...
      [top_y top_y bottom_y bottom_y top_y]);
+ 
+ %%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
+ 
+ figure(3);
+ 
+ [H,T,R] = hough(line_mask);
+       imshow(H,[],'XData',T,'YData',R,'InitialMagnification','fit');
+       xlabel('\theta'), ylabel('\rho');
+       axis on, axis normal, hold on;
+       P  = houghpeaks(H,5,'NHoodSize',[161,161])%ceil(0.3*max(H(:))));
+       x = T(P(:,2)); 
+       y = R(P(:,1));
+       plot(x,y,'s','color','white');
+ 
+  lines = houghlines(line_mask,T,R,P,'FillGap',5,'MinLength',7);
+       figure, imshow(im_rgb), hold on
+       max_len = 0;
+       for k = 1:length(lines)
+         xy = [lines(k).point1; lines(k).point2];
+         plot(xy(:,1),xy(:,2),'LineWidth',2,'Color','blue');
+ 
+         % plot beginnings and ends of lines
+         plot(xy(1,1),xy(1,2),'x','LineWidth',2,'Color','yellow');
+         plot(xy(2,1),xy(2,2),'x','LineWidth',2,'Color','red');
+ 
+         % determine the endpoints of the longest line segment 
+         len = norm(lines(k).point1 - lines(k).point2);
+         if ( len > max_len)
+           max_len = len;
+           xy_long = xy;
+         end
+       end
+ 
+       % highlight the longest line segment
+       plot(xy_long(:,1),xy_long(:,2),'LineWidth',2,'Color','cyan');
+ 
+ 
+ %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
 fprintf('Done\n');

branches/2010-image-rec/hist_test.m

@@ -127 +127 @@
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
  
+ figure(3);
+ 
  [H,T,R] = hough(line_mask);
        imshow(H,[],'XData',T,'YData',R,'InitialMagnification','fit');
@@ -141 +143 @@
        for k = 1:length(lines)
          xy = [lines(k).point1; lines(k).point2];
-         plot(xy(:,1),xy(:,2),'LineWidth',2,'Color','green');
+         plot(xy(:,1),xy(:,2),'LineWidth',2,'Color','blue');
  
          % plot beginnings and ends of lines