基于DWT-DCT-SVD算法的彩色图像水印嵌入和提取matlab仿真

clc; clear; close all; warning off; addpath(genpath(pwd)); Image=imread('lena.png'); Water=imread('water.png'); %嵌入水印 [Iwm,Uw,Vw]=EmbedWatermark(Image,Water); %加入攻击方式 %JPEG2000攻击 Iwm1 = JPEG2000Attack(Iwm,90); %JPEG攻击 Iwm2 = JPEGAttack(Iwm,90); %噪声攻击 Iwm3 = noiseAttack(Iwm,3); %滤波攻击 Iwm4 = FilterAttack(Iwm,3); %整形攻击 Iwm5 = SharpenAttack(Iwm); %均衡攻击 Iwm6 = histeqAttack(Iwm); %旋转攻击 Iwm7 = RotationAttack(Iwm); figure(2) subplot(241); imshow(uint8(Iwm1)) title('JPEG2000攻击'); subplot(242); imshow(uint8(Iwm2)) title('JPEG攻击'); subplot(243); imshow(uint8(Iwm3)) title('噪声攻击'); subplot(244); imshow(uint8(Iwm4)) title('滤波攻击'); subplot(245); imshow(uint8(Iwm5)) title('整形攻击'); subplot(246); imshow(uint8(Iwm6)) title('均衡攻击'); subplot(247); imshow(uint8(Iwm7)) title('旋转攻击'); figure(1) imshow(uint8(Iwm)) title('嵌入水印图像'); figure(3) subplot(241); [WaterR,WaterG,WaterB]=ExtractWatermark(Iwm1,Image,Uw,Vw,Water); NCR=ncc(uint8(WaterR),Water); NCG=ncc(uint8(WaterG),Water); NCB=ncc(uint8(WaterB),Water); a=[NCR,NCG,NCB]; [m,p]=max(a) if p==1 imshow(uint8(WaterR)) elseif p==2 imshow(uint8(WaterG)) else imshow(uint8(WaterB)) end title('JPEG2000攻击提取水印'); subplot(242); [WaterR,WaterG,WaterB]=ExtractWatermark(Iwm2,Image,Uw,Vw,Water); NCR=ncc(uint8(WaterR),Water); NCG=ncc(uint8(WaterG),Water); NCB=ncc(uint8(WaterB),Water); a=[NCR,NCG,NCB]; [m,p]=max(a) if p==1 imshow(uint8(WaterR)) elseif p==2 imshow(uint8(WaterG)) else imshow(uint8(WaterB)) end title('JPEG攻击提取水印'); subplot(243); [WaterR,WaterG,WaterB]=ExtractWatermark(Iwm3,Image,Uw,Vw,Water); NCR=ncc(uint8(WaterR),Water); NCG=ncc(uint8(WaterG),Water); NCB=ncc(uint8(WaterB),Water); a=[NCR,NCG,NCB]; [m,p]=max(a) if p==1 imshow(uint8(WaterR)) elseif p==2 imshow(uint8(WaterG)) else imshow(uint8(WaterB)) end title('噪声攻击提取水印'); subplot(244); [WaterR,WaterG,WaterB]=ExtractWatermark(Iwm4,Image,Uw,Vw,Water); NCR=ncc(uint8(WaterR),Water); NCG=ncc(uint8(WaterG),Water); NCB=ncc(uint8(WaterB),Water); a=[NCR,NCG,NCB]; [m,p]=max(a) if p==1 imshow(uint8(WaterR)) elseif p==2 imshow(uint8(WaterG)) else imshow(uint8(WaterB)) end title('滤波攻击提取水印'); subplot(245); [WaterR,WaterG,WaterB]=ExtractWatermark(Iwm5,Image,Uw,Vw,Water); NCR=ncc(uint8(WaterR),Water); NCG=ncc(uint8(WaterG),Water); NCB=ncc(uint8(WaterB),Water); a=[NCR,NCG,NCB]; [m,p]=max(a) if p==1 imshow(uint8(WaterR)) elseif p==2 imshow(uint8(WaterG)) else imshow(uint8(WaterB)) end title('整形攻击提取水印'); subplot(246); [WaterR,WaterG,WaterB]=ExtractWatermark(Iwm6,Image,Uw,Vw,Water); NCR=ncc(uint8(WaterR),Water); NCG=ncc(uint8(WaterG),Water); NCB=ncc(uint8(WaterB),Water); a=[NCR,NCG,NCB]; [m,p]=max(a) if p==1 imshow(uint8(WaterR)) elseif p==2 imshow(uint8(WaterG)) else imshow(uint8(WaterB)) end title('均衡攻击提取水印'); subplot(247); [WaterR,WaterG,WaterB]=ExtractWatermark(Iwm7,Image,Uw,Vw,Water); NCR=ncc(uint8(WaterR),Water); NCG=ncc(uint8(WaterG),Water); NCB=ncc(uint8(WaterB),Water); a=[NCR,NCG,NCB]; [m,p]=max(a) if p==1 imshow(uint8(WaterR)) elseif p==2 imshow(uint8(WaterG)) else imshow(uint8(WaterB)) end title('旋转攻击提取水印');