提取感兴趣区域

%%对于一幅图像的感兴趣区域获取
% Input_Image1=imread('G:\红外图片\tank300.bmp');
%Input_Image=imread('E:\魏晗论文\m文件(加上程序应用图片及说明)\m文件(加上程序应用图片及说明)\原车辆图像\a19.bmp');
Input_Image=imread('E:\09.3.1图像融合论文与资料\小论文\程序与图片1\shiyan\mt1r07081608viscrop.jpg');
% Input_Image=imread('E:\picture\lena.bmp');
% Input_Image=Input_Image1(3:end-3,3:end-3);
figure(1);
imshow(Input_Image);
%title('原始图像');
[h g]=imhist(Input_Image); % h 为像素个数，g 为灰度级
ff=double(Input_Image);
%/////////选择感兴趣区域具体程序段///////////////////
[M N]=size(Input_Image);

%%%%%%%%%%%%%%%%%%%%%%%%%%//输入300,得到初始阈值tg
tg=0;
for a=1:255
if (h(a)>300)&(h(a+1)<300)
if tg=tg
X(i,j)=1;
else
X(i,j)=0;
end
end
end
X=X~=0;

%标记连接成分（8邻域）
[LX XNum]=bwlabeln(X,8);
size(LX);

%以伪彩色的形式显示标记图像
RGBX=label2rgb(LX,@jet,'k');

figure(2);
subplot(221);imshow(Input_Image);title('原始图像');
subplot(222);imhist(Input_Image);title('原始图像直方图');
subplot(223);imshow(X);title('原始图像二值图');
subplot(224);imshow(RGBX);title('标记图像');
%/////////////////////////////////////////////////////////////////
%寻找最大区域作为参考区域
[r1 c1]=find(LX==1);
[a b]=size(r1);
MaxRegion=a;
Max_Region_flag=1;
Max_Region_r=r1;
Max_Region_c=c1;
for i=2:XNum
[r,c] = find(LX == i);
[a b]=size(r);
if a>MaxRegion
MaxRegion=a;
Max_Region_flag=i;
Max_Region_r=r;
Max_Region_c=c;
end
end
%确定参考区域的四个顶点坐标

rmin=min(min(Max_Region_r));
rmax=max(max(Max_Region_r));
cmin=min(min(Max_Region_c));
cmax=max(max(Max_Region_c));

%确定参考区域的行数和列数
Ref_row=rmax-rmin+1;
Ref_column=cmax-cmin+1;

%将参考区域图像赋给一个新的图像矩阵Ref_Image
Ref_Image=zeros(M,N);
for i=1:M
for j=1:N
if(i>rmin & icmin & jM
ROImaxr=M;
end
if ROIminc<0
ROIminc=1;
end
if ROImaxc>N
ROImaxc=N;
end

%将感兴趣区域图像赋给一个新的图像矩阵ROI_Image
ROI_Image=zeros(M,N);
for i=1:M
for j=1:N
if(i>ROIminr & iROIminc & jROIminr & iROIminc & j=tg
TWO_Image(i,j)=255;
else
TWO_Image(i,j)=0;
end
else
TWO_Image(i,j)=0;
end
end
end
figure(5);
imshow(uint8(TWO_Image));
title('二值化感兴趣区域');