图像匹配与识别

甘肃政法学院本科生实验报告(五）实验课程名称:数字图像处理技实验日期: 2016年 月日开课时间：2015-2016学年第二 学期甘肃政法学院实验管理中心印制实验题目图像匹配与识别图像Matlab数字水印设计小组合作姓名班级学号一、实验目的1、 图像匹配与识别2、 图像Matlab数字水印设计二．实验环境Matlab7.0三、实验内容与步骤matlab分为可感知和不易感知两种。1、可感知水印。这是一种可以看见的水印，就像插入或覆盖在图像上的标Internet2、不易感知水印。这是一种应用更加广泛的水印，主要用来鉴别产品的真者可以从中提取出标记，从而证明物品为某人所有。二、图像匹配与识别图像识别技术是人工智能的一个重要领域，人们提出了不同的图像识别模是一种程。四、实验过程与分析水印处理：clearall;closeall;clc;M=256;%原图像长度N=64;%水印长度[filename1,pathname]=uigetfile('*.*','selecttheimage');image1=imread(num2str(filename1));subplot(2,2,1);imshow(image1);title('originalimage'); %imageforwatermarkingimage1=double(image1);imagew=imread('dmg2.tif');subplot(2,2,2);imshow(imagew);title('originalwatermark'); %originalwatermark%嵌入水印[ca,ch,cv,cd]=dwt2(image1,'db1');[cas,chs,cvs,cds]=dwt2(ca,'db1');fori=1:Nforj=1:Nifimagew(i,j)==0a=-1;elsea=1;endCa(i,j)=cas(i,j)*(1+a*0.03);endendIM=idwt2(Ca,chs,cvs,cds,'db1');markedimage=double(idwt2(IM,ch,cv,cd,'db1'));%显示嵌入后水印图像subplot(2,2,3);colormap(gray(256));image(markedimage);title('markedimage');imwrite(markedimage,gray(256),'watermarked.bmp','bmp');%提取水印image1=imread(num2str(filename1));image1=double(image1);imaged=imread('watermarked.bmp');[ca,ch,cv,cd]=dwt2(image1,'db1');imaged=imread('watermarked.bmp');[ca,ch,cv,cd]=dwt2(image1,'db1');[cas,chs,cvs,cds]=dwt2(ca,'db1');[caa,chh,cvv,cdd]=dwt2(imaged,'db1');[caas,chhs,cvvs,cdds]=dwt2(caa,'db1');forp=1:Nforq=1:Na=caas(p,q)/cas(p,q)-1;ifa<0W(p,q)=0;elseW(p,q)=255;endendend%显示提取的水印subplot(2,2,4);colormap(gray(256));image(W);title('从含水印图像中提取的水印');imwrite(W,gray(256),'watermark.bmp','bmp');显示图像的傅立叶变换FFT及逆变换IFFT。I=imread('lena.bmp');subplot(1,3,1),imshow(uint8(I)); %显示原始图J=fft2(I); %图像FFT变换B=fftshift(J);%显示FFT变换结果subplot(1,3,2),imshow(log(abs(B)),[]);H=ifft2(J);%使用逆变换复原图像并显示结果subplot(1,3,3),imshow(uint8(H));计算并显示图像的离散余弦变换DCT及逆变换IDCT。程序代码如下：I=imread('lena.bmp');subplot(1,4,1),imshow(uint8(I);)%显示原始图像J=dct2(I);DCTJ(abs(J)<100)=0;%将较小的系数设置为0N=idct2(J);subplot(1,4,4),imshow(uint8(N));二值图像匹配m=imread(filepath); %读取原图figure,imshow(im,[]);title('Raw'); %显示原图im=im2bw(im); %转二值图像figure,imshow(im,[]),title('BW'); %显示二值图im2=imfill(im,'holes'); %填充im3=bwperim(im2); %轮廓提取figure,imshow(im2,[]);title('') %显figure,imshow(im3,[]);图像识别：clear;clc;reference_img=imread('optical1.tif');%提取原始图片作为基准target_img=imread('optical2.tif');%要求进行匹配的图像subplot(2,3,1);%显示原始图像imshow(reference_img);title('原始图像');subplot(2,3,2);imshow(target_img);title('目标图像');[Ix,Iy]=size(target_img);%x=[160,103,102,209];%自己选定的特征点y=[116,246,160,299];u=[87,35,34,141];v=[21,151,66,204];Px=polyfit(x,u,1);%Py=polyfit(y,v,1);fori=1:Ixforj=1:Iym=i*Px(1,1)+Px(1,2);%n=j*Py(1,1)+Py(1,2);m_integer=floor(m);%n_integer=floor(n);m_decimal=m-m_integer;%对映射值取小数n_decimal=n-n_integer;if(ge(m_integer,Ix)||ge(n_integer,Iy)||lt(m_integer,0)||lt(n_integer,0)||m_integer==0||n_integer==0)%判断映射点是否满足条件result_img(i,j)=0;result_img3(i,j)=0;elseresult_img(i,j)=target_img(m_integer,n_integer);%没有进行插值result_img3(i,j)=(1-m_decimal)*(n_decimal*target_img(m_integer,n_integer+1)+(1-n_decimal)*target_img(m_integer,n_integer))+m_decimal*(n_decimal*target_img(m_integer+1,n_integer+1)+(1-n_decimal)*target_img(m_integer+1,n_integer));%双线性插值endifresult_img(i,j)==0%把没有进行插值图像和原始图像进行组合result_img2(i,j)=reference_img(i,j);elseresult_img2(i,j)=result_img(i,j);endifresult_img3(i,j)==0%把插值后图像和原始图像进行组合result_img4(i,j)=reference_img(i,j);ifresult_img3(i,j)==0%把插值后图像和原始图像进行组合result_img4(i,j)=reference_img(i,j);elseresult_img4(i,j)=result_img3(i,j);endendendtitle('未插值图像匹配');titl