查看: 2674|回复: 0

MATLAB DCT水印源程序代码解析

字体大小: 正常放大

1198 主题	4 听众	2976 积分

该用户从未签到

电梯直达

1^#

发表于 2023-11-16 11:50 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

离散余弦变换（Discrete Cosine Transform，DCT）是一种在信号处理和数据压缩领域广泛应用的数学变换方法。其核心思想是将一系列数据点通过一组余弦函数的加权和转换，从而将信号或图像从时域/空域表示转换为频域表示。DCT在信息编码和压缩中扮演着重要的角色，尤其在图像和音频领域，它通过能量集中的方式，更高效地表示信号的频率分量。
DCT的原理包括以下几个关键点：

1.变换域： DCT将输入的序列或图像分解为一组频率成分，每个分量对应不同频率，类似于傅里叶变换，但DCT的基函数是实数。
2.余弦函数： DCT使用一组余弦函数作为基函数，这些函数在不同频率上振荡，与傅里叶变换相似但更适合实际应用。
3.变换步骤：对于一维序列，DCT将N个数据点转换为N个频率分量。对于二维图像，DCT将图像分块并对每个块应用变换，得到每个块的频率分量。
4.能量集中：与傅里叶变换不同，DCT倾向于将信号或图像的能量集中在少数低频分量上，这在图像压缩中尤为有用。
5.类型：存在不同类型的DCT变换，其中类型II是最常用的，特别是在图像压缩中得到广泛应用。

DCT的应用广泛，其中包括：

6.图像压缩： JPEG图像压缩利用DCT将图像转换为频域表示，通过去除或量化高频分量，实现了高效的图像压缩。
7.音频处理： MP3等音频压缩格式使用DCT将音频信号转换为频域表示，以减小文件大小同时保留主要信息。
8.信号处理：在通信领域，DCT用于数据传输和信号处理，有助于分析信号的频率成分。

总体而言，DCT是一种强大的数学工具，能够以更适合分析和处理的方式表示信号或图像，为多个领域提供了有力的支持。现在，我们将通过DCT的原理和应用，探讨如何将水印嵌入图像的具体过程。

现在我们使用DCT将水印提案加到图片中
这段代码执行了图像水印嵌入的过程，以下是逐行代码的解释：

clear;clc, A7 J( Q* f1 ^. ~. ]6 M
# v- d& |$ N4 l3 N: m9 a q1 Q\" Z
size=512;2 U# E\" |7 _ _# n1 W* M3 Y9 z
block=8;
5 T6 o! o+ a) P4 I' z+ x\" k+ x7 Y' w
blockno=size/block;
0 b3 M h, d3 _2 _2 h1 r/ ^
LENGTH=size*size/64;
. g* f. D1 P; S/ C) ~4 c, Y3 ^( p
Alpha1=0.02; Alpha2=0.02;1 P- @' u. P4 j6 J c! [! h5 \
T1=100;
\" R7 ]9 h% J3 L/ m( F4 W, h# g* p5 G
I=zeros(size,size);+ Q0 `: { g& H) Y, r
D=zeros(size,size);& X& @: q. t% y5 e I3 T
BW=zeros(size,size); s& J\" b z( G- X
block_dct1=zeros(block,block);
& L) d( Y* p' K: ?5 [; J; i
$ O, M, o) L1 i, x3 t( t\" g7 l4 {
i=imread('watermark64by64.png');
8 m8 t$ O7 l) w\" ]1 b; J
mark=reshape(i,1,LENGTH);
. W, o3 g: E' U
figure;imshow(i);title('水印');
9 {# }$ W& X2 `$ j6 w4 f
%subplot(3,2,1);plot(mark);title('水印');& x6 P+ d- O7 n8 N. I
%显示原图
6 B3 V% s) g$ q9 r$ N% [: j1 p
figure;I=imread('lena512512.jpg');
' q( o, I) k# y3 b* A
+ s& B* i2 Y8 ]
I0=I;
6 \: L- v Z6 w5 O V! W3 X: b
imshow(I);title('原始图像');
R& O: F# F2 y1 }, D\" o# Z) z
( J- x\" n+ W' \0 w& f; K6 O
% BW=edge(I,'canny');figure;;imshow(BW);2 i( o) t; m7 m: T& w8 D
% title('edge of origine image');! i+ j3 Q& K0 ~ s3 E8 e m
%嵌入水印
0 {\" Y* }9 h* P$ p# `3 ?9 F
k=1;9 `, T/ y8 _2 b# t' o
for m=1:blockno# l Q/ W) E7 D% k. M
for n=1:blockno
' N. X H% [8 [
x=(m-1)block+1; y=(n-1)block+1;; m& @; ]% f1 [! X6 a8 y
block_dct1=I(x:x+block-1,y:y+block-1);
b/ N9 C7 Q* S, P9 d6 \$ T, Z
block_dct1=dct2(block_dct1);( z\" X\" p E2 f) G2 h
BW_8_8=BW(x:x+block-1,y:y+block-1);! ] h3 f$ F# E/ w0 [* V
if m<=1|n<=1
* R7 V* V5 ~8 c* R* d
T=0;4 u: |# a* {4 j( x6 \! Y
else# ^# E. W: U+ T& b
T=sum(BW_8_8); T=sum(T);
7 Z; U6 z5 w. I k
end8 X\" r/ X! t7 s3 L: ^2 @
if T>T1+ b- S* b: r! L) H1 h$ A5 H
Alpha=Alpha2;3 v8 \, ]1 i }6 [
else
/ E6 l. s; N6 q5 C$ u/ }8 n
Alpha=Alpha1;4 R% }* ^% ]+ ~( x9 k
end
) u1 I' T1 `' C3 J\" Y' {
block_dct1(1,1)=block_dct1(1,1)(1+Alphamark(k));( O& g9 `3 @7 B# V4 S
block_dct1=idct2(block_dct1);4 M' J9 I% ? R$ ^6 H- c
D(x:x+block-1,y:y+block-1)=block_dct1;2 ]' u+ L% c! N7 h' p* o
k=k+1;
: i\" C( J# T. N- X. G5 u
end! }\" a1 m0 X$ g\" i$ c: G
end
/ s7 _! R0 x. M. ]* q6 A/ J
6 j8 s) D! k) c& G' A' w
figure;imshow(D,[]);title('嵌入水印后的图像');

复制代码

1.clear; clc;: 这两个命令清空了 MATLAB 的工作区和命令窗口，以确保在开始之前没有残留的变量或输出。
2.size = 512; block = 8; blockno = size / block; LENGTH = size * size / 64; Alpha1 = 0.02; Alpha2 = 0.02; T1 = 100;: 这些是设置图像处理所需的变量和参数。size 定义了图像的尺寸，block 是用于分块的大小，blockno 是图像分块后的数量，LENGTH 是图像长度，Alpha1 和 Alpha2 是用于水印嵌入的参数，T1 是一个阈值。
3.I = zeros(size, size); D = zeros(size, size); BW = zeros(size, size); block_dct1 = zeros(block, block);: 这些命令初始化了一些矩阵和数组，用于存储图像、水印、DCT 变换后的块等。
4.i = imread('watermark64by64.png'); mark = reshape(i, 1, LENGTH);: 从文件 'watermark64by64.png' 中读取水印图像，然后将其重新形状化为一个向量 mark。
5.figure; imshow(i); title('水印');: 这段代码显示了读取的水印图像。
6.figure; I = imread('lena512512.jpg'); I0 = I; imshow(I); title('原始图像');: 从文件 'lena512512.jpg' 中读取原始图像，并显示它。
7.接下来的部分是对原始图像进行边缘检测，但是注释掉了：BW = edge(I, 'canny'); figure; imshow(BW); title('edge of origine image');。因为这部分代码被注释掉了，所以不会实际执行。
8.for m = 1:blockno for n = 1:blockno ... end end: 这是一个嵌套循环，用于对图像进行分块处理和水印嵌入。
9.在嵌入水印的循环中，它首先将图像按照设定的块大小进行处理，对每个块应用 DCT 变换，然后根据水印像素的值和设定的阈值 T1 进行水印嵌入的处理。
10.最后，它显示了嵌入水印后的图像：figure; imshow(D, []); title('嵌入水印后的图像');。

%提取水印$ o0 k7 J, T$ l, P9 q
k=1;
l5 V- n1 ?8 p/ i6 q9 y4 Z q f
mark1=[];
7 P$ _8 r/ B2 B/ j: c$ {
for m=1:blockno
1 G7 t* u& _7 @( x, F
for n=1:blockno
9 y. ~5 g w& o3 D
x=(m-1)*block+1; y=(n-1)*block+1;
8 j) v9 {4 c& t1 G$ ~
block_dct1=D(x:x+block-1,y:y+block-1);1 r6 ?% e& k/ [
block_dct0=I0(x:x+block-1,y:y+block-1);
\" `$ r& R6 k2 O* L$ {2 F
block_dct1=dct2(block_dct1);
/ H( v3 s4 r# u. S/ e! ?6 _# ?9 w. E7 Q
block_dct0=dct2(block_dct0);
8 `( o; {+ k+ I4 X\" J; P0 d
BW_8_8=BW(x:x+block-1,y:y+block-1);
5 c- Y\" P! \7 f! I( D3 T3 \9 q0 l
if m<=1|n<=17 U. D+ T, f) t6 ?\" |
T=0;
2 g( z8 ?! b7 p* ]
else! p, G8 i6 b) p# }1 F& E5 [\" n& z\" [
T=sum(BW_8_8); T=sum(T);
' X2 k# J* ]# q
end+ x9 E c+ B, h' @2 y6 H
if T>T17 \% |4 h4 O; Y
Alpha=Alpha2;
3 l- y9 P/ R4 J6 ?% i
else
8 N\" }6 p l6 d8 r: y
Alpha=Alpha1;
$ B: I9 v8 l+ `! c
end
& m' @/ J( s0 g
2 l5 K, Y. J. n. n2 i- d
mark1(k)=(block_dct1(1,1)/block_dct0(1,1)-1)/Alpha;
4 n\" f5 J+ g+ Q1 w8 ?- i9 `
k=k+1;
6 M* M# _, l# K5 P; ^5 c
end
) c$ w0 u0 N0 O$ v0 r
end9 _3 V5 r8 O$ M
mark1=reshape(mark1,64,64);
! Y3 g7 i3 A5 {. O/ u
figure
+ L- t3 w& l* g* K1 u# v) q& e\" ^
imshow(mark1)
6 w4 ]. x. f( Q9 G f
title('提取的水印')

复制代码

这部分代码是关于从嵌有水印的图像中提取水印信息的过程。让我为你解释一下：

1.k=1; mark1=[];: 这两个语句初始化了一个索引变量 k 和一个空数组 mark1，用于存储提取的水印信息。
2.for m=1:blockno for n=1:blockno ... end end: 这是一个嵌套循环，用于遍历图像的分块。
3.在循环内部，它首先确定当前块的位置 x 和 y，然后分别从嵌有水印的图像 D 和原始图像 I0 中提取对应的块。
4.对提取的块应用 DCT 变换：block_dct1 = dct2(block_dct1); block_dct0 = dct2(block_dct0);。这里对嵌有水印的块和原始图像的块都进行了 DCT 变换。
5.计算了当前块的边缘信息 BW_8_8。
6.根据设定的阈值 T1，确定 Alpha 值。这与嵌入水印的过程中相似。
7.接下来的代码是用于从嵌有水印图像和原始图像中提取水印信息的关键步骤：

mark1(k) = (block_dct1(1,1) / block_dct0(1,1) - 1) / Alpha;

这个式子计算了水印信息，并将其存储到 mark1 数组中。

8.循环结束后，mark1 存储了提取出的水印信息。接下来，代码对 mark1 进行重新形状化为 64x64 的矩阵，并显示提取的水印图像。

最后的 imshow(mark1) 代码会显示提取出的水印图像。请确保之前的嵌入过程已经执行，并且使用相同的图像来提取水印，以获取正确的结果。