Công trình t i [5] c a chúng tôi là qui trình n thông tin c i ti n qui trình c b n trên mi n contourlet, ý t ng đ c đ xu t b i nhóm Mohan [12]. Theo đó, quá trình nhúng đ c ti n hành qua hai giai đo n (hình 2.9), thu t toán nhúng QIM đ c s d ng trong c hai giai đo n.
- Giai đo n 1:Nhúng w trên subband t t nh t.
Theo nhóm tác gi Haifeng Li [32], các h s đ nhúng watermark t t nh t là nh ng h s subband có n ng l ng l n nh t. V i các h s trong subband này, thu t toán c n ti n hành ch n các h s nhúng nh m đ m b o b ba ràng bu c:
44
(ii) Kh n ng ph c h i trong quá trình trích thông tin. (iii) Tính trong su t c a nh nhúng.
Hình 2.9 Qui trình nhúng thông tin hai giai đo n c i ti n [5]
Ta g i các h s đ c ch n nhúng là h s kh thi.
- Giai đo n 2: Nhúng trên lowpass c a nh k t qu t giai đo n 1.
Trong quá trình nhúng d a trên hai giai đo n này, giá tr b c nh y delta c a QIM s đi u ch nh khác nhau, tùy thu c vào đ c tr ng c a h s nhúng.
Phân tích contourlet nh g c I0, ta có m t subband lowpass là IJ và các nh theo h ng t i bandpass 穴珍(,鎮賃乳), k = 0, 1, …, 2鎮乳-1; trong đó, j=1, 2, …, J và j là m c th j
nh g c I0
Bi n đ i Contourlet
Tính n ng l ng các sub-band
QIM cho subband có n ng l ng Max trên các h s kh thi
Bi n đ i Contourlet I’
nh đã nhúng I’ (G 1)
QIM trên h s lowpass c a I’
Bi n đ i Contourlet ng c Watermark w
Bi n đ i Contourlet ngu c
45
c a phân tích LP, k là nh phân tích theo h ng bandpass th k b ng dàn l c DFB lj
m c. - N ng l ng c a subband đ c tính b ng công th c 2.6: ∑∑ = = = M m N n l k j k j d m n MN E j 1 2 1 ) ( , , 1 ( , ) (2.6) v i M, N là kích th c c a dj,k
T i giá tr dmaxlà subband có n ng l ng l n nh t, ch n dmax=max(Ej,k) đ nhúng. Do tính ch t c a h s t i subband giai đo n 1 và lowpass giai đo n 2 khác nhau, nên b c nh y delta c a QIM t i hai l n nhúng s khác nhau. Khi đó, b c nh y dài (giá tr delta l n) s làm gi m giá tr đo PSNR (tính trong su t s gi m) h n so v i b c nh y delta ng n. Tuy nhiên, n u delta quá nh thì kho ng cách không đ xa đ ánh x riêng bi t các giá tr nhúng khác nhau (xem 2.1.2), nên giá tr đo NC c a watermark sau khi trích s th p đ ng ngh a v i watermark trích đ c không chính xác.
M t đ c tính n a c a h s contourlet t i các subband là tính láng gi ng. Do m c tiêu thu t toán nhúng là ch ng t n công mà không đ t yêu c u cao v dung l ng nhúng. Vì th , khi th c hi n nhúng QIM giai đo n 1, chúng tôi không nhúng m t đ dày mà th c hi n nhúng th anh hình 2.10.a, các ô đ c tô màu là nh ng h s đ c ch n đ nhúng.
Hình 2.10.a Các h s kh thi t i subband đ m b o b ba ràng bu c nhúng giai đo n 1
Trên c s d a theo qui trình t i hình 2.9, chúng tôi ti p t c c i ti n đ mô hình đ t k t qu t t h n. C th có th mô t nh sau, đ t ng đ trong su t cho mô hình, trong l n nhúng giai đo n 1, chúng tôi phân tích nh g c thành đ c t ch n và l
46
r i m i th c hi n bi n đ i contourlet. Thu t toán nhúng QIM khi đó đ c th c hi n trên contourlet c a đ c t l . Sau khi nhúng, hai t p h s contourlet c a đ c t ch n và l s đ c bi n đ i contourlet ng c, cu i cùng đem tr n hai đ c t ch n và l l i đ có nh ch a watermark Iw.
Ph ng pháp phân tích ch n l (Hình 2.10.b) đ c th c hi n nh sau:
+ T i m t pixel v trí (dòng i, c t j) s đ a vào ph n đ c t ch n n u (i+j) chia h t cho 2, ngh a là(i+j) mod 2=0.
+ Ng c l i, đ a pixel vào ph n đ c t l .
Hình 2.10.b Mô t ch n và l
tính các pixel còn thi u đ c th c hi n nh sau:
- T i đ c t l : Theo hình 2.10.b s thi u pixel (1,1), (1,3), (2,2), (3,1), (3,3). Các giá tr s tính ra b ng trung bình c a các láng gi ng l c a pixel đó. V i (1,1), láng gi ng là (1,2) và (2,1); (2,2) có 4 láng gi ng (1,2), (2,1), (2,3), (3,2).
47
Qui trình nhúng giai đo n 1 đ c th c hi n c i ti n nh trong hình 2.11 d i đây:
Hình 2.11 Qui trình nhúng thông tin giai đo n 1 theo đ c t ch n, l [26]