2020-05-30

在二值图像中标识出20*20的区域

YUAN=imread('BEI.bmp');
A = rgb2gray(YUAN);
[M,N]=size(A);
t=graythresh(A); %确定灰度阈值I=imbinarize(A,t);%h=rectangle('position',[1,160,20,20],'edgecolor','b');%在二值图像中标识出20*20区域
D=uint8(I); %强行转换类型
for i = 1:M    
for j = 1:N      
  if I(i,j)           
   D(i,j) = 255;  %逻辑1即为逻辑真，对应255灰度值        else       
        D(i,j) = 0;  %逻辑0即为逻辑假，对应0灰度值     
           end 
              end
              end
              h=drawRect( D, [1,160], [20,20],  1, [255 0 0] );%输入：% src：        原始图像，可以为灰度图，可为彩色图% pt：         左上角坐标   [x1, y1]% wSize：   框的大小      [wx, wy]% lineSize： 线的宽度% color：     线的颜色      [r,  g,  b] %输出：% h：      标注好的图      画好了的图像

function [ dest ] = drawRect( src, pt, wSize,  lineSize, color )
%简介：% %将图像画上有颜色的框图，如果输入是灰度图，先转换为彩色图像，再画框图% 图像矩阵% 行向量方向  是  y% 列向量方向  是  x%---------------------------------------------------------------------
-%输入：
% src：        原始图像，可以为灰度图，可为彩色图
% pt：         左上角坐标   [x1, y1]
% wSize：   框的大小      [wx, wy]
% lineSize： 线的宽度
% color：     线的颜色      [r,  g,  b] 
%----------------------------------------------------------------------
%输出：% dest：           画好了的图像%---------------------------------------------------------------------- %flag=1: 有缺口的框
%flag=2: 无缺口的框flag = 2; 
 %判断输入参数个数
 if nargin < 5
     color = [255 255 0];
     end 
     if nargin < 4   
      lineSize = 1;
      end 
      if nargin < 3    
      disp('输入参数不够 !!!');    
      return;
      end     
      %判断框的边界问题
      [yA, xA, z] = size(src);
      x1 = pt(1);
      y1 = pt(2);
      wx = wSize(1);
      wy = wSize(2);
      if  x1>xA || ...      
        y1>yA||...     
           (x1+wx)>xA||...   
                (y1+wy)>yA     
                disp('画的框将超过图像 !!!');  
                  return;
                  end %如果是单通道的灰度图，转成3通道的图像if 1==z    
                  dest(:, : ,1) = src;    
                  dest(:, : ,2) = src;    
                  dest(:, : ,3) = src;else   
                   dest = src;end %开始画框图
                   for c = 1 : 3                 %3个通道，r，g，b分别画   
                    for dl = 1 : lineSize   %线的宽度，线条是向外面扩展的      
                      d = dl - 1;     
                         if  1==flag %有缺口的框       
                              dest(  y1-d ,            x1:(x1+wx) ,  c  ) =  color(c); %上方线条          
                                dest(  y1+wy+d ,     x1:(x1+wx) , c  ) =  color(c); %下方线条         
                                   dest(  y1:(y1+wy) ,   x1-d ,           c  ) =  color(c); %左方线条        
                                       dest(  y1:(y1+wy) ,   x1+wx+d ,    c  ) =  color(c); %左方线条      
                                         elseif 2==flag %无缺口的框            dest(  y1-d ,            (x1-d):(x1+wx+d) ,  c  ) =  color(c); %上方线条     
                                                dest(  y1+wy+d ,    (x1-d):(x1+wx+d) ,  c  ) =  color(c); %下方线条      
                                                      dest(  (y1-d):(y1+wy+d) ,   x1-d ,           c  ) =  color(c); %左方线条         
                                                         dest(  (y1-d):(y1+wy+d) ,   x1+wx+d ,    c  ) =  color(c); %左方线条    
                                                             end 
                                                                end 
                                                                   end %主循环尾  end %函数尾