I=imread('brand02.jpg');
HS_I=double(I);
tic
G=gradient(HS_I);  % 計(jì)算梯度值
[y,x]=size(I);
T=60;
%%%%%%% 商標(biāo)文字分割 %%%%%%%%%%
n=30;
GY=zeros(y,1);                 
for j=1:y   
    for i=1:x-n
         Max=max(G(j,i:i+n));
         Min=min(G(j,i:i+n));
         MGD(j,i)=Max-Min;       % 計(jì)算最大梯度差
         if (MGD(j,i)>T)         
             GY(j,1)=GY(j,1)+1;  % MGD大于閾值的次數(shù)統(tǒng)計(jì)
         end    
    end
end    
[temp MaxY]=max(GY);
PY1=MaxY;
while ((GY(PY1,1)>=n)&&(PY1>1))  % 次數(shù)閾值檢測(cè)
       PY1=PY1-1;
end    
PY2=MaxY;
while ((GY(PY2,1)>=n)&&(PY2<y))
       PY2=PY2+1;
end
%%%%%%%%%% X 方向的分割 %%%%%%%%%%%
n=15;
GX=zeros(1,x);
for i=1:x
    for j=1:y-n
         Max=max(G(j:j+n,i));
         Min=min(G(j:j+n,i));
         MGD(j,i)=Max-Min;
         if (MGD(j,i)>T)
             GX(1,i)=GX(1,i)+1;
         end
     end
 end
 PX1=1;
while (((GX(1,PX1)<=n)||(GX(1,PX1+1)<=n))&&(PX1<x)) % 列方向保證同時(shí)相鄰的兩列滿足次數(shù)要求
       PX1=PX1+1;
end    
PX2=x;
while (((GX(1,PX2)<=n)||(GX(1,PX2-1)<=n))&&(PX2>PX1))
       PX2=PX2-1;
end    

HS_I=uint8(HS_I);
IY=HS_I(PY1:PY2,PX1:PX2);
t=toc
%%%%%%%%%%%%%%%%%%%%%%%%%%
figure,imshow(HS_I);
figure,plot(MGD(40,:));grid
figure,plot(MGD(140,:));grid
figure,imshow(IY);