图像融合 MATLAB.rar_matlab_matlab 图像_图像融合

clear clc bkg_img =imread('1.jpg'); edge_margin=30; insert_flag=0; bkg_img=rgb2gray(bkg_img); subplot(231),imshow(bkg_img); title('原图像'); [X0,Y0]=size(bkg_img); %得到多次膨胀后的最大连通区域坐标 B=[1 1 1 1 1 1 1 1 1]; img2=im2bw(bkg_img); iimg2=~img2; A2=imdilate(iimg2,B);%图像A1被结构元素B膨胀 A3=imdilate(A2,B); A4=imdilate(A3,B); A5=imdilate(A4,B); A6=imdilate(A5,B); A7=imdilate(A6,B); subplot(232),imshow(A7); title('膨胀后图像'); img_reg = regionprops(A7, 'area', 'boundingbox'); areas = [img_reg.Area]; rects = cat(1, img_reg.BoundingBox); %显示最大连通区域 [~, max_id] = max(areas); max_rect = rects(max_id, :); %读取threat img 旋转T1 img=imread('2.jpg'); subplot(233),imshow(img); title('原图像'); if ndims(img)==3 img=rgb2gray(img); end [X1,Y1]=size(img); for i_img=1:X1 for j_img=1:Y1 img(i_img,j_img)=255-img(i_img,j_img); end end subplot(234),imshow(img); title('危险品图像'); angle=360*rand(1,1); img=imrotate(img, angle,'bilinear','loose'); [X2,Y2]=size(img); for i_img=1:X2 for j_img=1:Y2 img(i_img,j_img)=255-img(i_img,j_img); end end subplot(235),imshow(img); title('图像'); scale_down=1.05-0.1*rand; imgscale=imresize(img,scale_down); [Xscale,Yscale]=size(imgscale); [x1,y1,x2,y2]=tip_in_maxrect(max_rect(2),max_rect(1),max_rect(4),max_rect(3),Xscale,Yscale,edge_margin); bkg_img_double=im2double(bkg_img); imgscale_double=im2double(imgscale); if x1==0&&x2==0 insert_flag=0; else for i_img=x1+1:x2 for j_img=y1+1:y2 bkg_img_double(i_img,j_img)= bkg_img_double(i_img,j_img)*imgscale_double(i_img-x1,j_img-y1); end img_syn=cat(3,bkg_img_double,bkg_img_double,bkg_img_double); end subplot(236),imshow(img_syn); title('融合图像'); insert_flag=1; end