• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    一種抗數(shù)字印刷/拍攝的數(shù)字水印算法研究及驗(yàn)證

    2022-02-22 11:48:04張明貴高靜
    包裝工程 2022年3期
    關(guān)鍵詞:通濾波數(shù)字水印魯棒性

    張明貴,高靜

    一種抗數(shù)字印刷/拍攝的數(shù)字水印算法研究及驗(yàn)證

    張明貴1,高靜2

    (1.四川寬窄印務(wù)有限責(zé)任公司,四川 德陽,618100;2.四川中煙工業(yè)有限責(zé)任公司 長城雪茄廠,四川 什邡 618400)

    數(shù)字印刷和手機(jī)拍攝是檢驗(yàn)含有數(shù)字水印的印刷產(chǎn)品版權(quán)有效性的主要途徑,解決在數(shù)字印刷和手機(jī)拍攝場景下水印信息難以提取這一問題。提出一種基于離散小波變換和奇異值分解的抗數(shù)字印刷/手機(jī)拍攝的數(shù)字水印算法。首先對宿主圖像進(jìn)行二級(jí)離散小波變換,對低頻子帶進(jìn)行奇異值分解,將Arnold置亂后的水印嵌入到奇異值分解的對角矩陣中,最后進(jìn)行小波逆變換得到嵌入水印圖像。含水印信息的峰值信噪比(PSNR)達(dá)到了39.67,在不同攻擊下的歸一化相關(guān)系數(shù)(NC)均在0.89以上,實(shí)際的印刷拍攝測試中NC值也均達(dá)到0.94,證明該水印算法在保證不可見性的同時(shí)對于印刷產(chǎn)品實(shí)際應(yīng)用中的常見攻擊具有很好的抵抗能力。對不同品牌的智能手機(jī)、印刷設(shè)備以及不同種類的紙張采集到的含水印圖像進(jìn)行測試,均能夠在印刷和拍攝場景下實(shí)現(xiàn)水印提取,證明了該算法的有效性。

    數(shù)字水印;打印拍攝;離散小波變換;奇異值分解

    隨著防偽技術(shù)的發(fā)展,信息隱藏技術(shù)在產(chǎn)品版權(quán)保護(hù)方面起到了巨大作用。它使得產(chǎn)品能具備一定的防偽功能,維護(hù)消費(fèi)者和企業(yè)的合法權(quán)益,進(jìn)一步地保證和促進(jìn)產(chǎn)品銷售,因此,數(shù)字水印技術(shù)作為信息隱藏技術(shù)研究領(lǐng)域的重要分支,是國際信息安全研究的一個(gè)重要研究熱點(diǎn)[1]。其主要的信息防偽形式是將一段文本、標(biāo)識(shí)、序列號(hào)、圖像等數(shù)字信號(hào)作為水印嵌入到載體上,同時(shí)使用一定的技術(shù)檢測方法可以將水印正確地提取出來,從而達(dá)到防偽驗(yàn)證的目的[2]。印刷包裝產(chǎn)品作為信息表達(dá)的基本載體,可以通過數(shù)字水印防偽技術(shù)的應(yīng)用,在不影響產(chǎn)品外包裝信息的條件下,起到良好的防偽效果。

    數(shù)字印刷成為了小批量印包產(chǎn)品的主流方 式[3-5]。加載水印信息的圖像后,數(shù)字印刷后圖像的低頻部分能較為完整地保存下來,但其高頻部分大部分會(huì)缺失,相當(dāng)于受到了低通濾波的攻擊。同時(shí),印刷過程中內(nèi)部各個(gè)元器件的噪聲和點(diǎn)增益都會(huì)引起印刷的圖像質(zhì)量受損,即受到了噪聲的攻擊。

    除此之外,隨著智能手機(jī)的迅速普及,通過拍照功能來進(jìn)一步識(shí)別水印信息也變成更加智能的驗(yàn)證方式,常規(guī)的圖像水印算法難以在手機(jī)拍攝后識(shí)別出完整的水印信息[6]。一些研究人員開始尋求通過移動(dòng)設(shè)備檢測水印的方法。Pramila等[7]詳細(xì)地分析了拍攝過程對于數(shù)字圖像的影響,在拍攝的過程中,相機(jī)的鏡頭會(huì)造成桶形和枕形失真。同時(shí)智能手機(jī)拍攝圖像獲取原始圖像時(shí),圖像縮放比例不同,手機(jī)拍攝圖像時(shí)存在一定的夾角,因而導(dǎo)致圖像旋轉(zhuǎn)。手機(jī)拍攝圖像不可避免地會(huì)有輕微的抖動(dòng)、壓縮,導(dǎo)致圖像模糊。

    目前數(shù)字水印算法典型的空間域水印算法是最低有效位算法(LSB)。該算法采用預(yù)先設(shè)置好的密鑰,由序列發(fā)生器產(chǎn)生隨機(jī)信號(hào),之后將產(chǎn)生的隨機(jī)信號(hào)按一定規(guī)則組合成二維形式插入到原始圖像的最低幾位,完成數(shù)字水印的嵌入[8]。最低有效位算法隱藏的信息較為脆弱,對于幾何攻擊魯棒性較差。在主流的變換域方法中,小波變換由于能較好地解決時(shí)頻分辨率的矛盾,且與新一代國際壓縮標(biāo)準(zhǔn)相兼容,從而得到了廣泛的應(yīng)用[9]。離散小波變換(DWT)是進(jìn)行信號(hào)分析非常有用的工具[10]。通過離散小波變換在原始圖像的低頻子帶中嵌入水印,可以很好地解決數(shù)字印刷高頻信息缺失這一問題。同時(shí)奇異值分解(SVD)是數(shù)值分析中基本工具之一[11],由于圖像的奇異值分解對許多圖像處理以及旋轉(zhuǎn)、縮放、平移等幾何攻擊具有一定的不變性,奇異值分解的算法可以很好地抵抗以幾何攻擊為主的拍攝攻擊。陳青等[6]提出了一種基于SIFT變換和DWT-SVD的自適應(yīng)圖像水印算法,含水印圖像在沒有受到攻擊時(shí),可以無損地提取出水印,在受到攻擊時(shí),提取出的水印品質(zhì)較好,且在面對小規(guī)模幾何攻擊時(shí),也有良好的穩(wěn)定性。

    文中考慮到數(shù)字水印在數(shù)字印刷包裝產(chǎn)品上的應(yīng)用,綜合考慮DWT和SVD的優(yōu)點(diǎn),將Arnold加密后的水印圖像嵌入原圖像中,在不改變水印圖像像素值的同時(shí)更改像素位置,實(shí)現(xiàn)一種基于離散小波變換和奇異值分解相結(jié)合的數(shù)字水印算法,從而完成抗印刷和拍攝數(shù)字水印方案,為個(gè)性化印刷包裝產(chǎn)品數(shù)字水印的實(shí)現(xiàn)提供基本的算法工具和應(yīng)用范例。

    1 水印算法基礎(chǔ)設(shè)計(jì)方法

    1.1 離散小波變換(DWT)

    離散小波變換是小波變換中的一種[12-13],兼具時(shí)、頻兩域的多分辨率分析能力,已被廣泛應(yīng)用于計(jì)算機(jī)領(lǐng)域,如數(shù)字圖像處理、數(shù)字信號(hào)處理、語音識(shí)別等,以及量子物理和石油勘探等諸多領(lǐng)域。

    研究小波變換的理論基礎(chǔ)源于傅里葉變換,具體變換公式為[14]:

    (1)

    (2)

    式(1)為傅里葉變換基本公式,其中eiwt函數(shù)可由歐拉公式直接展開,因此復(fù)雜的時(shí)域信號(hào)經(jīng)過變換之后可由一系列的正余弦波形構(gòu)成,但該信號(hào)只能夠反映頻域特性。式(2)為小波變換的基本公式,該公式存在2個(gè)變量,平移量和尺度。尺度對應(yīng)頻率,小波函數(shù)的伸縮性質(zhì),平移量對應(yīng)時(shí)間,體現(xiàn)了小波函數(shù)的平移特性,因此小波變換同時(shí)具備時(shí)頻特性。

    在圖像變換過程中,主要利用小波變換的多分辨率(及多尺度)特性分解圖像。對于數(shù)字圖像,會(huì)分別從水平方向和垂直方向進(jìn)行低通濾波和高通濾波。如圖1所示,一級(jí)離散小波變換之后圖像被分為4個(gè)子帶,其中(近似子帶)是由2個(gè)方向進(jìn)行低通濾波后產(chǎn)生的系數(shù),包含了圖像的主要信息,能夠近似地表示圖像。(水平子帶)、(垂直子帶)、(對角子帶)表示了圖像的細(xì)節(jié)或邊緣信息[15]。

    圖1 一級(jí)離散小波變換

    1.2 奇異值分解(SVD)

    奇異值分解是一種基于特征值和特征向量的數(shù)學(xué)分析工具,已被廣泛地應(yīng)用于圖像處理領(lǐng)域[16],如模式識(shí)別、圖像融合、數(shù)據(jù)壓縮、數(shù)字水印等,以及統(tǒng)計(jì)學(xué)等諸多領(lǐng)域。從矩陣論的角度分析,一幅數(shù)字圖像可以視為一個(gè)非負(fù)矩陣。設(shè)為×的數(shù)字圖像,用∈R來表示一個(gè)數(shù)字圖像矩陣,其中表示實(shí)數(shù)域。圖像經(jīng)過奇異值分解后,表示為:

    (3)

    式中:和分別為左奇異值矩陣和右奇異值矩陣;為一個(gè)非對角線上的項(xiàng)都是0的矩陣,其對角線上的元素滿足以下條件:

    (4)

    式中:為的秩,它等于非零奇異值的數(shù)量,于是有:

    (5)

    從圖像處理的角度來看,奇異值分解特性如下所述。

    1)奇異值分解對所要進(jìn)行變換的矩陣大小沒有限制,可以是正方矩陣也可以是長矩陣。

    2)一幅圖像的奇異值有很好的穩(wěn)定性,這意味著當(dāng)圖像受到輕微擾動(dòng)時(shí),它的奇異值不會(huì)發(fā)生劇烈的改變。對于裁切、縮放、旋轉(zhuǎn)等幾何攻擊干擾,圖像的奇異值依然會(huì)保持穩(wěn)定。

    3)奇異值能夠表現(xiàn)出圖像的內(nèi)在的矩陣元素關(guān)系而不是視覺特性,反映的是圖像的內(nèi)蘊(yùn)特性。所以利用奇異值分解的特性,可以提高圖像水印的抗攻擊性能。

    由于奇異值分解具備著以上特性,可以將其引入到抗數(shù)字印刷/拍攝的數(shù)字水印算法中,能夠有效地提升算法的抗打印拍攝能力。

    1.3 抗數(shù)字印刷和拍攝的印刷品數(shù)字水印算法設(shè)計(jì)

    水印模型包括了2個(gè)體系:水印嵌入和水印提取。在水印嵌入后,對含水印圖像進(jìn)行檢測,確定其含有水印,進(jìn)行水印嵌入的逆過程,即可完成水印的提取。

    如前所述,圖像經(jīng)過離散小波變換分解成1個(gè)低頻分量和3個(gè)高頻分量。低頻分量包含了圖像的大部分信息,該分量具有較高的視覺容量。相反,高頻分量能量較小,視覺容量較低,不能夠嵌入過多的水印信息,且易遭受低通濾波和噪聲的影響,如果在高頻分量中嵌入水印會(huì)使得算法不能夠有效地抵抗打印拍攝攻擊,因此,在離散小波變換分解后的低頻分量中嵌入水印,能夠更好地抵抗打印拍攝攻擊,并平衡水印的魯棒性和不可見性這一對矛盾。

    考慮到離散小波變換對旋轉(zhuǎn)、剪切等幾何攻擊魯棒性較差,而奇異值分解對于幾何攻擊的魯棒性較強(qiáng),因此將離散小波變換與奇異值分解相結(jié)合,將經(jīng)歷2次離散小波變換得到的低頻分量進(jìn)行奇異值分解,在分解的2分量中嵌入水印。

    水印嵌入算法的具體流程如下述所。

    1)讀取一幅原始圖像,對原始圖像進(jìn)行二級(jí)離散小波變換,得到2、2、2、2、1、1、17個(gè)子帶。

    2)對2低頻子帶進(jìn)行矩陣奇異值分解,得到左奇異值矩陣和右奇異值矩陣,以及對角矩陣。

    (6)

    3)讀取水印圖像,采用Arnold置亂變換算法進(jìn)行水印的加密。運(yùn)用加性水印公式(為水印的嵌入強(qiáng)度,為水印)將加密后的水印嵌入到矩陣中,之后對新產(chǎn)生的矩陣再進(jìn)行一次奇異值分解,得到新產(chǎn)生的左奇異值矩陣1和右奇異值矩陣1,以及對角矩陣1。

    (7)

    4)對、1、進(jìn)行逆奇異值分解,得到嵌有水印數(shù)據(jù)的低頻子圖2*。

    (8)

    5)對圖像進(jìn)行二級(jí)逆離散小波變換,顯示含水印圖像w,完成水印嵌入。

    水印提取算法是嵌入算法的逆過程,流程如下所述。

    1)讀取含水印圖像w,對待檢測的含水印圖像w進(jìn)行二級(jí)離散小波變換(DWT),得到2'、2'、2'、2'、1'、1'、1' 7個(gè)子帶,對低頻子帶2'進(jìn)行矩陣奇異值分解(SVD),得到左奇異值矩陣和右奇異值矩陣,以及對角矩陣

    (9)

    2)將得到的矩陣和式(10)中的1和1矩陣進(jìn)行逆奇異值分解,得到。

    (10)

    3)按照加性水印公式,根據(jù)和,得到提取的水印圖像1。

    (11)

    4)對提取的水印圖像1進(jìn)行逆Arnold置亂變換,解密得到水印圖像。

    2 實(shí)驗(yàn)結(jié)果分析與討論

    文中使用Matlab 2018b作為編程工具完成仿真實(shí)驗(yàn),選擇256像素×256像素的灰度Lena圖像(圖2a)作為原始載體圖像,選擇64像素×64像素的含有“SUST”字樣的二值圖像(圖2b)作為待嵌入的水印信息,水印的嵌入強(qiáng)度選擇為0.4。

    首先對原始圖像進(jìn)行二級(jí)離散小波變換,對低頻子帶2進(jìn)行矩陣奇異值分解,根據(jù)式(7)中的水印加性公式將置亂后的水印信息嵌入到原始圖像中,完成水印的嵌入,通過對含水印圖像進(jìn)行主觀評價(jià)和客觀評價(jià),完成水印的不可見性分析。其次對嵌入水印后的圖像進(jìn)行剪切、旋轉(zhuǎn)等幾何攻擊,以及不同程度的噪聲、低通濾波和移動(dòng)模糊等常見信號(hào)攻擊,通過對比提取出的水印圖像和原始水印圖像,采用歸一化相關(guān)系數(shù)值評價(jià)水印算法的抗攻擊能力,完成水印的魯棒性分析。最后對嵌入水印信息的圖像進(jìn)行不同打印機(jī)、不同紙張、不同手機(jī)的印刷/拍攝測試,通過觀察提取出的水印圖像,應(yīng)用歸一化相關(guān)系數(shù)值評價(jià)水印算法的魯棒性,完成算法的印刷拍攝實(shí)際測試。

    圖2 原始載體圖像和原始水印圖像

    2.1 水印的不可見性分析

    算法采用式(7)加性水印公式,將二值“SUST”水印圖像以0.4的嵌入強(qiáng)度嵌入到原始Lena圖中,得到含水印圖像。對水印不可見性而言,一般有主觀評價(jià)和客觀評價(jià)2種評價(jià)方法。主觀上多采用人的肉眼觀察的方式進(jìn)行評價(jià),客觀上衡量水印的不可見性一般采用峰值信噪比(Peak Signal To Noise Ratio, PSNR)算法,來對比含水印圖像w(,)與原始圖像(,)之間的差異,峰值信噪比的計(jì)算式為:

    (12)

    式中:M、N為原始載體圖像的大?。籌ij為原始載體圖像;I'ij為含水印圖像,一般來說,當(dāng)PSNR≥36時(shí),可以認(rèn)為水印算法的不可見性良好,滿足水印的不可見性需求。

    圖3 原始圖像和含水印圖像

    Fig.3 Original image and watermark image

    根據(jù)1.3節(jié)中的水印嵌入算法,對原始圖像(圖3a)進(jìn)行水印嵌入,得到嵌入水印后的圖像(圖3b)。通過人的視覺系統(tǒng)對比觀察圖3a和圖3b,二者并沒有明顯差異,置入水印后的圖片基本沒有變化。同時(shí),通過計(jì)算得到PSNR=39.67>36,滿足水印的不可見性客觀評價(jià)標(biāo)準(zhǔn),因此水印的不可見性良好,文中算法滿足水印的不可見性要求。

    2.2 水印的魯棒性分析

    如前所述,印刷和拍攝過程對于數(shù)字圖像主要存在2個(gè)方面的影響,一方面是印刷過程中出現(xiàn)的低通濾波攻擊、噪聲攻擊,另一方面是拍攝過程中出現(xiàn)的剪切、旋轉(zhuǎn)等幾何攻擊和移動(dòng)模糊攻擊,因此應(yīng)用各種攻擊算法對圖3b所示的含水印圖像進(jìn)行旋轉(zhuǎn)、剪切幾何攻擊以及噪聲、低通濾波和移動(dòng)模糊等攻擊,以檢測水印算法的抗印刷拍攝的能力。

    文中采用歸一化相關(guān)系數(shù)NC值來衡量水印的魯棒性,NC值的計(jì)算式為:

    (13)

    式中:w為原始水印圖像;w'為提取水印圖像;和分別為水印圖像的長和寬。通過計(jì)算得到含水印圖像遭受各種攻擊分別提取出的水印圖 像的NC值。一般情況下,當(dāng)NC值大于0.8時(shí), 就可以確認(rèn)為有效水印,否則為無效水印。攻擊圖像和提取出的水印圖像以及對應(yīng)的NC值見圖4—13。

    2.2.1 剪切攻擊下的水印提取效果分析

    文中通過剪切算法對含水印圖像的左上角、上部、中間正方形區(qū)域和中間豎直區(qū)域進(jìn)行了剪切,應(yīng)用1.3節(jié)所述的水印提取算法進(jìn)行水印提取,提取出的水印圖像以及NC值見圖4。圖4展示了不同位置剪切攻擊下的水印提取結(jié)果,可以看到,對含水印圖像的左上角、上部、中間正方形部分進(jìn)行剪切,可以完整地提取出水印圖像,提取出的水印圖像的NC值可達(dá)0.97,對含水印圖像的中間豎直部分進(jìn)行剪切,可以較為完整地提取出水印圖像,提取出的水印圖像的NC值可達(dá)0.96,可以清晰地分辨出“SUST”字樣。說明該水印算法未出現(xiàn)因剪切攻擊造成魯棒性下降的情況。相反,文中算法對各種復(fù)雜的剪切攻擊具有較強(qiáng)的抵抗能力,提取出的水印圖像清晰可見。

    2.2.2 旋轉(zhuǎn)攻擊下的水印提取效果分析

    文中通過旋轉(zhuǎn)算法對含水印圖像進(jìn)行了1°~10°的二維旋轉(zhuǎn),步長為1°。應(yīng)用1.3節(jié)中所述的水印提取算法對旋轉(zhuǎn)后的圖像進(jìn)行水印提取,提取出的水印圖像和NC值見圖5和圖6。

    圖4 不同位置裁切的Lena圖及提取的水印

    圖5 旋轉(zhuǎn)攻擊實(shí)驗(yàn)結(jié)果

    圖6 旋轉(zhuǎn)攻擊提取水印的NC值

    圖5展示了3°、5°、7°、9°旋轉(zhuǎn)攻擊下的水印提取結(jié)果,可以看到,在3°、5°旋轉(zhuǎn)攻擊下,提取出的水印清晰可見,NC值分別為0.98和0.96。由圖6可知,隨著旋轉(zhuǎn)度數(shù)逐漸提升,提取水印的NC值逐漸下降,但即使是在9°旋轉(zhuǎn)攻擊下,NC值仍然能夠達(dá)到0.89,可以清晰地分辨出“SUST”字樣。由于在人手持智能手機(jī)進(jìn)行拍攝時(shí),智能手機(jī)與印刷在紙張上的圖像的夾角一般小于5°,所以旋轉(zhuǎn)攻擊對水印圖像的提取影響甚微,該算法能夠滿足水印的抗旋轉(zhuǎn)需求。

    2.2.3 移動(dòng)模糊攻擊下的水印提取效果分析

    文中通過移動(dòng)模糊算法對含水印圖像進(jìn)行了步長為1,從1~10個(gè)像素值的移動(dòng)模糊攻擊,利用1.3節(jié)中所敘述的水印提取方法分別進(jìn)行水印提取,提取出的水印圖像以及對應(yīng)的NC值見圖7和圖8。

    圖7 移動(dòng)模糊攻擊實(shí)驗(yàn)結(jié)果

    圖8 移動(dòng)模糊攻擊提取水印的NC值

    圖7展示了部分像素值移動(dòng)模糊下的水印提取結(jié)果,可以看出,在逆時(shí)針方向移動(dòng)7個(gè)像素值的移動(dòng)模糊攻擊下,可以完整地提取出清晰的水印圖像,NC值為0.98。由圖8可知,隨著逆時(shí)針移動(dòng)像素值的提高,提取水印的NC值逐漸降低,提取出的水印圖像逐漸模糊,但即使是在10 像素值的移動(dòng)模糊攻擊下,NC值仍然可達(dá)0.94,依然可以較為清晰地分辨出“SUST”字樣,因此該算法可以有效地抵抗移動(dòng)模糊攻擊。

    2.2.4 高斯噪聲攻擊下的水印提取效果分析

    文中應(yīng)用高斯噪聲算法對含水印圖像進(jìn)行了步長為0.01,方差為0.01~0.1的高斯噪聲攻擊,采用1.3節(jié)述的提取算法對高斯噪聲攻擊后的圖像進(jìn)行水印提取,提取出的水印圖像以及NC值見圖9和圖10。

    圖9 高斯噪聲攻擊實(shí)驗(yàn)結(jié)果

    圖10 高斯噪聲攻擊提取水印的NC值

    圖9展示了方差為0.01、0.03、0.05、0.07的高斯噪聲攻擊下的水印提取結(jié)果,可以看到,隨著高斯噪聲攻擊強(qiáng)度不斷增大,圖像中的高斯噪聲逐漸增加,含水印圖像逐漸遭到破壞,不可避免地出現(xiàn)了大量的噪點(diǎn)。由圖10可知,提取水印的NC值隨著高斯噪聲方差的增加而降低。盡管受到了高斯噪聲的影響,提取水印的NC值仍然可達(dá)0.91以上,提取出的水印圖像仍然清晰可見,表明文中算法對高斯噪聲有著較強(qiáng)的魯棒性。

    2.2.5 椒鹽噪聲攻擊下的水印提取效果分析

    文中通過椒鹽噪聲算法對含水印圖像進(jìn)行了步長為1,1~15密度的椒鹽噪聲攻擊,利用1.3節(jié)所闡述的水印提取算法分別進(jìn)行水印提取,提取出的水印圖像以及NC值見圖11和圖12。

    圖11 椒鹽噪聲攻擊實(shí)驗(yàn)結(jié)果

    圖11展示了1%、5%、10%、15%密度的椒鹽噪聲攻擊下的水印提取結(jié)果。可以看出,在1%、5%、10%密度椒鹽噪聲的攻擊下,可以完整地提取出清晰的水印圖像,NC值分別為1、0.99和0.97。由圖12可知,隨著椒鹽噪聲攻擊強(qiáng)度的不斷增大,提取水印的NC值逐漸降低。當(dāng)椒鹽噪聲攻擊強(qiáng)度增大到15%時(shí),含有水印的圖像受到較大的破壞,在提取出的水印圖像上出現(xiàn)了大量的噪點(diǎn),但即使是在15%密度的椒鹽噪聲攻擊下,NC值依然高達(dá)0.95,仍然可以較為清晰地分辨出“SUST”字樣,表明該算法對椒鹽噪聲的魯棒性較強(qiáng)。

    2.2.6 低通濾波攻擊下的水印提取效果分析

    文中通過低通濾波算法對含水印圖像分別進(jìn)行了3×3、5×5、7×7、9×9的低通濾波攻擊,利用1.3節(jié)中所提出的水印提取算法分別對其進(jìn)行水印提取,提取出的水印圖像以及NC值見圖13。由圖13可以看出,在3×3和5×5的低通濾波的攻擊下,NC值為1,在7×7和9×9的低通濾波的攻擊下,NC值為0.99,可以非常清晰地分辨出“SUST”字樣,表明該算法對低通濾波的魯棒性較強(qiáng)。

    2.3 印刷拍攝實(shí)際測試

    為了驗(yàn)證水印信息在打印拍攝過程中受到像素失真和幾何失真后的穩(wěn)健不變性,使用Matlab 2018b將水印信息嵌入到原始載體圖像中,將含水印圖像保存。使用EPOSON、科美數(shù)碼等不同型號(hào)的打印機(jī)將含水印圖像打印在不同的紙張上,并使用不同型號(hào)的智能手機(jī)拍攝含水印圖像得到印刷拍攝后的圖像,修改圖像大小為256像素×256像素,使用Matlab 2018b將水印信息提取出來。打印拍攝修改大小后的含水印圖像和提取出的水印信息見圖14,提取出的水印圖像的NC值見表1。

    由圖14可以看出,經(jīng)過印刷拍攝后幾幅Lena圖像都發(fā)生了輕微的旋轉(zhuǎn),并且略微變暗。使用EPSON1390打印在噴墨打印紙上的2幅圖像提出來的水印NC值高達(dá)0.97,可以清晰地分辨出水印圖像。使用EPSON數(shù)碼打印機(jī)打印在相紙上,并用蘋果手機(jī)提取出的水印圖像NC值高達(dá)0.96,可以清晰地分辨出“SUST”字樣。使用科美數(shù)碼印刷機(jī)印刷在普通打印紙上的2幅圖像提取出的水印和使用EPSON數(shù)碼打印機(jī)打印在相紙上,并用華為mate 9提取出的水印圖像NC值均達(dá)到了0.95,可以較為清晰地分辨出水印圖像上的內(nèi)容。使用科美數(shù)碼印刷機(jī)印刷在銅版紙上的2幅圖像NC值均達(dá)到了0.94,可以較為清晰地識(shí)別出水印圖像上的“SUST”文字。

    圖13 低通濾波攻擊實(shí)驗(yàn)結(jié)果

    表1 印刷拍攝后的檢測結(jié)果

    Tab.1 Test results after printing and shooting

    由表1可知,提取出的水印的歸一化相關(guān)系數(shù)NC值均可達(dá)0.94以上。另外由圖14可以看出,提取出的水印圖像中的“SUST”字樣均較為清晰。結(jié)果表明,文中算法對于幾種不同型號(hào)的打印機(jī)打印、不同種類的紙張和不同型號(hào)的手機(jī)拍攝均具有良好的抗印刷拍攝的能力。

    3 結(jié)語

    文中通過分析數(shù)字印刷/手機(jī)拍攝過程對水印圖像的各種影響,根據(jù)印刷/拍攝特點(diǎn),提出了一種基于離散小波變換和奇異值分解的抗數(shù)字印刷/手機(jī)拍攝水印算法。通過仿真實(shí)驗(yàn)結(jié)果可以看出,該算法恢復(fù)的水印信息能夠獲得良好的圖像視覺效果,對于剪切、旋轉(zhuǎn)等幾何攻擊和噪聲、模糊、低通濾波等常規(guī)的信號(hào)處理具有很強(qiáng)的抵抗能力,符合水印算法對魯棒性和不可見性的要求。同時(shí)對于不同型號(hào)的打印機(jī)、智能手機(jī)和不同種類的紙張,均能夠在印刷和拍攝場景下實(shí)現(xiàn)水印提取,證明了該算法的有效性。

    [1] RAY A, ROY S. Recent Trends in Image Watermarking Techniques for Copyright Protection: A Survey[J]. International Journal of Multimedia Information Retrieval, 2020, 9(4): 249-270.

    WANG Jun, WAN Wen-bo. A Novel Attention-Guided JND Model for Improving Robust Image Watermarking[J]. Multimedia Tools and Applications, 2020,79(33/34): 24057-24073.

    [2] 霍佳琦. 抗打印掃描數(shù)字水印算法研究與實(shí)現(xiàn)[D]. 杭州: 浙江大學(xué), 2019: 30-34.

    HUO Jia-qi. Research and Implementation of Anti-Printing Scanning Digital Watermarking Algorithm[D]. Hangzhou: Zhejiang University, 2019: 30-34.

    [3] RIAD R, DOUZI H, El M, et al. A Print-Scan Resilient Watermarking Based on Fourier Transform and Image Restoration[J]. International Journal of Computer Applications, 2015, 128(15): 13-17.

    [4] 張雯, 孟婕. 數(shù)字水印技術(shù)在印刷防偽領(lǐng)域應(yīng)用的研究進(jìn)展[J]. 包裝工程, 2017, 38(15): 143-148.

    ZHANG Wen, MENG Jie. Progress of Applying Digital Watermarking Technology in Printing Anti-Counterfeiting[J]. Packaging Engineering, 2017, 38(15): 143-148.

    [5] 陳青, 伍東升. 基于SIFT變換和DWT-SVD的自適應(yīng)圖像水印算法[J]. 包裝工程, 2021, 42(9): 254-259.CHEN Qing, WU Dong-sheng. Adaptive Image Watermarking Algorithm Based on SIFT Transform and DWT-SVD[J]. Packaging Engineering, 2021, 42(9): 254-259.

    [6] PRAMILA A, KESKINARKAUS A, TAKALA V, et al. Extracting Watermarks from Printouts Captured with Wide Angles Using Computational Photography[J]. Multimedia Tools and Applications, 2017, 76(15): 16063-16084.

    [7] NAZARI M, MEHRABIAN M. A Novel Chaotic IWT-LSB Blind Watermarking Approach with Flexible Capacity for Secure Transmission of Authenticated Medical Images[J]. Multimedia Tools and Applications, 2021, 80(7): 10615-10655.

    [8] CHOPRA A, GUPTA S, DHALL S. Analysis of Frequency Domain Watermarking Techniques in Presence of Geometric and Simple Attacks[J]. Multimedia Tools and Applications, 2020, 79(1/2): 501-554.

    [9] 劉琰, 周理. 基于小波變換域的數(shù)字圖像嵌入和提取方法[J]. 沈陽工業(yè)大學(xué)學(xué)報(bào), 2019, 41(1): 68-72.

    LIU Yan, ZHOU Li. Digital Image Embedding and Extracting Method Based on Wavelet Transform Domain[J]. Journal of Shenyang University of Technology, 2019, 41(1): 68-72.

    [10] ZERMI N, KHALDI A, KAFI R, et al. A DWT-SVD Based Robust Digital Watermarking for Medical Image Security[J]. Forensic Science International, 2021, 320: 110691.

    [11] CHEN Bing-quan, CUI Jin-ge, XU Qing, et al. Coupling Denoising Algorithm Based on Discrete Wavelet Transform and Modified Median Filter for Medical Image[J]. Journal of Central South University, 2019, 26(1): 120-131.

    [12] 楊明鵬, 葉幗華. 基于小波變換的三維網(wǎng)格模型盲水印算法[J]. 福建師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2020, 36(3): 32-36.

    YANG Ming-peng, YE Guo-hua. A Blind Watermarking of 3D Models Based on Wavelet Transform[J]. Journal of Fujian Normal University (Natural Science Edition), 2020, 36(3): 32-36.

    [13] POONAM, ARORA S M. A DWT-SVD Based Robust Digital Watermarking for Digital Images[J]. Procedia Computer Science, 2018, 132: 1441-1448.

    [14] 李海洋. 基于離散小波變換的抗噪圖像水印處理[J]. 信息化研究, 2019, 45(5): 23-26.

    LI Hai-yang. Anti-Noise Image Watermarking Processing Based on Discrete Wavelet Transform[J]. Informatization Research, 2019, 45(5): 23-26.

    [15] WANG S, MENG X, YIN Y, et al. Optical Image Watermarking Based on Singular Value Decomposition Ghost Imaging and Lifting Wavelet Transform[J]. Optics and Lasers in Engineering, 2019, 114: 76-82.

    [16] 陳星娥, 鄭薇瑋, 吳小菁, 等. 基于奇異值分解的QR電子票數(shù)字水印防偽方法[J]. 牡丹江師范學(xué)院學(xué)報(bào)(自然科學(xué)版), 2020(1): 41-44.

    CHEN Xing-ge, ZHENG Wei-wei, WU Xiao-jing, et al. An Digital Watermarking Method for QR Electronic Ticket Based on Singular Value Decomposition[J]. Journal of Mudanjing Teachers College (Natural Sciences Edition), 2020(1): 41-44.

    Research and Verification of a Digital Watermarking Algorithm Resistant to Digital Printing/Photography

    ZHANG Ming-gui1,GAO Jing2

    (1.Sichuan Kuanzhai Printing Co., Ltd., Deyang 618100, China; 2.The Great Wall Cigar Factory, China Tobacco Sichuan Industrial Co., Ltd., Shifang 618400, China)

    Digital printing and mobile shooting are the main ways to verify the validity of copyright of printed products containing digital watermarks. Therefore, the work aims to solve the problem that it is difficult to extract watermark in digital printing and mobile shooting. A digital watermarking algorithm based on discrete wavelet transform and singular value decomposition to resist digital printing/mobile shooting was proposed. Firstly, the host image was processed by two-stage discrete wavelet transform, and singular value decomposition was performed on the low-frequency sub-bands. The watermark after Arnold was scrambled was embedded into the diagonal matrix of singular value decomposition. Finally, the embedded watermark image was obtained by inverse wavelet transform. The peak signal to noise ratio (PSNR) containing watermark information reached 39.67. The normalized correlation coefficient (NC) under different attacks was above 0.89. The NC value in actual shooting test also reached 0.94, showing that watermark algorithm not only ensured invisibility but also had good resistance to common attacks in practical application of printing products. The watermark images collected from different brands of smart phones, printing devices and different kinds of paper can be extracted in printing and shooting scenes, which proves the effectiveness of the algorithm.

    digital watermarking; printing and shooting; discrete wavelet transform; singular value decomposition

    TP391.4

    A

    1001-3563(2022)03-0274-10

    10.19554/j.cnki.1001-3563.2022.03.034

    2021-08-18

    張明貴(1973—),男,工程師,主要研究方向?yàn)閿?shù)字印刷和色彩管理。

    猜你喜歡
    通濾波數(shù)字水印魯棒性
    聲吶發(fā)射機(jī)負(fù)載阻抗變化仿真分析
    荒漠綠洲區(qū)潛在生態(tài)網(wǎng)絡(luò)增邊優(yōu)化魯棒性分析
    基于網(wǎng)屏編碼的數(shù)字水印技術(shù)
    電子制作(2019年20期)2019-12-04 03:51:30
    基于確定性指標(biāo)的弦支結(jié)構(gòu)魯棒性評價(jià)
    二階有源低通濾波電路的計(jì)算機(jī)輔助設(shè)計(jì)
    電子制作(2017年1期)2017-05-17 03:54:07
    基于非支配解集的多模式裝備項(xiàng)目群調(diào)度魯棒性優(yōu)化
    西南交通大學(xué)學(xué)報(bào)(2016年6期)2016-05-04 04:13:11
    基于復(fù)帶通濾波的智能電表量測算法及其DSP實(shí)現(xiàn)
    基于頻域分析和低通濾波的光伏并網(wǎng)逆變器諧振抑制研究
    基于數(shù)字水印的人臉與聲紋融合識(shí)別算法
    精品久久久久久久毛片微露脸| 久久精品国产亚洲av香蕉五月 | 制服诱惑二区| 这个男人来自地球电影免费观看| 女人被狂操c到高潮| 免费在线观看影片大全网站| 亚洲一区二区三区不卡视频| 一级毛片女人18水好多| 午夜久久久在线观看| 国产欧美日韩综合在线一区二区| 国产成人免费观看mmmm| 在线国产一区二区在线| 午夜视频精品福利| 99国产精品免费福利视频| 亚洲av熟女| 久久久精品国产亚洲av高清涩受| 丁香六月欧美| 久久久久久人人人人人| 精品人妻熟女毛片av久久网站| 交换朋友夫妻互换小说| av网站在线播放免费| 少妇猛男粗大的猛烈进出视频| 黄网站色视频无遮挡免费观看| 亚洲片人在线观看| 久久99一区二区三区| 老司机午夜十八禁免费视频| 看黄色毛片网站| 热re99久久精品国产66热6| 亚洲国产看品久久| 久久99一区二区三区| 色婷婷av一区二区三区视频| 亚洲专区字幕在线| 亚洲美女黄片视频| 免费女性裸体啪啪无遮挡网站| 国产熟女午夜一区二区三区| 怎么达到女性高潮| 女同久久另类99精品国产91| 亚洲,欧美精品.| 69av精品久久久久久| 十八禁人妻一区二区| 国产麻豆69| 后天国语完整版免费观看| 国产亚洲一区二区精品| 欧美人与性动交α欧美软件| 搡老熟女国产l中国老女人| 最近最新中文字幕大全电影3 | 麻豆国产av国片精品| 久久久久精品国产欧美久久久| 三上悠亚av全集在线观看| 久久午夜综合久久蜜桃| 欧美成人午夜精品| 国产欧美日韩精品亚洲av| 日韩欧美免费精品| 在线av久久热| 丝袜美足系列| 18禁观看日本| 俄罗斯特黄特色一大片| 国产成人欧美在线观看 | 成人亚洲精品一区在线观看| 老司机午夜十八禁免费视频| 免费日韩欧美在线观看| 亚洲国产毛片av蜜桃av| 午夜亚洲福利在线播放| 夜夜夜夜夜久久久久| 精品高清国产在线一区| 久久久国产欧美日韩av| 国产成人精品无人区| 18禁裸乳无遮挡动漫免费视频| 亚洲午夜理论影院| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲一码二码三码区别大吗| 免费av中文字幕在线| 欧美日韩瑟瑟在线播放| 精品人妻在线不人妻| 久久99一区二区三区| 成年女人毛片免费观看观看9 | 一级a爱片免费观看的视频| 国产一区二区激情短视频| 欧美亚洲 丝袜 人妻 在线| 夫妻午夜视频| 天堂俺去俺来也www色官网| 午夜免费观看网址| 高清黄色对白视频在线免费看| 18禁裸乳无遮挡动漫免费视频| 熟女少妇亚洲综合色aaa.| 天天躁日日躁夜夜躁夜夜| 国产精品98久久久久久宅男小说| 天天影视国产精品| 欧美精品人与动牲交sv欧美| 国产97色在线日韩免费| 高清黄色对白视频在线免费看| 自拍欧美九色日韩亚洲蝌蚪91| 国产在视频线精品| av欧美777| а√天堂www在线а√下载 | 一本一本久久a久久精品综合妖精| 国产精品久久视频播放| 国产日韩一区二区三区精品不卡| 欧美日韩亚洲高清精品| 99国产综合亚洲精品| av片东京热男人的天堂| 欧美日韩国产mv在线观看视频| av天堂在线播放| 天天添夜夜摸| 欧美日韩亚洲国产一区二区在线观看 | 欧美不卡视频在线免费观看 | 精品一区二区三卡| 18禁美女被吸乳视频| 久久久国产精品麻豆| 1024视频免费在线观看| 高清av免费在线| 一区二区三区国产精品乱码| 国产高清视频在线播放一区| 久久精品熟女亚洲av麻豆精品| 免费人成视频x8x8入口观看| 自拍欧美九色日韩亚洲蝌蚪91| 欧美乱码精品一区二区三区| 日韩欧美国产一区二区入口| 久久人妻福利社区极品人妻图片| av视频免费观看在线观看| 亚洲全国av大片| 一区二区日韩欧美中文字幕| 欧美日本中文国产一区发布| 国产精品美女特级片免费视频播放器 | av网站在线播放免费| 夜夜躁狠狠躁天天躁| 99国产精品一区二区三区| 999久久久国产精品视频| 亚洲成人国产一区在线观看| 亚洲国产欧美一区二区综合| 高清av免费在线| 亚洲成人手机| 午夜视频精品福利| 欧美老熟妇乱子伦牲交| 日本a在线网址| 欧美日韩福利视频一区二区| 精品第一国产精品| 精品福利观看| 国产真人三级小视频在线观看| 成人亚洲精品一区在线观看| 亚洲成a人片在线一区二区| 精品国产国语对白av| 国产高清视频在线播放一区| 久久国产精品大桥未久av| 99精国产麻豆久久婷婷| 在线永久观看黄色视频| 婷婷精品国产亚洲av在线 | 女人爽到高潮嗷嗷叫在线视频| 一边摸一边做爽爽视频免费| 一二三四社区在线视频社区8| 最近最新免费中文字幕在线| av天堂久久9| 不卡一级毛片| 国产欧美日韩精品亚洲av| 久久人妻熟女aⅴ| 波多野结衣一区麻豆| 亚洲中文日韩欧美视频| 老汉色∧v一级毛片| 一区福利在线观看| 一区在线观看完整版| 搡老乐熟女国产| 欧美不卡视频在线免费观看 | 欧美 亚洲 国产 日韩一| 免费在线观看影片大全网站| 他把我摸到了高潮在线观看| 99热国产这里只有精品6| 国产人伦9x9x在线观看| 国产乱人伦免费视频| 精品国产美女av久久久久小说| а√天堂www在线а√下载 | 国产精品香港三级国产av潘金莲| 亚洲一码二码三码区别大吗| 一级片'在线观看视频| 久久精品亚洲熟妇少妇任你| 国产又爽黄色视频| 久久精品熟女亚洲av麻豆精品| 一级片'在线观看视频| 久久久精品区二区三区| 免费观看精品视频网站| 久久精品国产99精品国产亚洲性色 | 热99re8久久精品国产| 亚洲人成伊人成综合网2020| 国产成人免费无遮挡视频| 亚洲专区国产一区二区| 国产91精品成人一区二区三区| 久久久久久久精品吃奶| 欧美一级毛片孕妇| 色老头精品视频在线观看| 乱人伦中国视频| 两个人免费观看高清视频| 视频区欧美日本亚洲| 国产黄色免费在线视频| av一本久久久久| 精品国产国语对白av| 婷婷精品国产亚洲av在线 | 狂野欧美激情性xxxx| 欧美日韩国产mv在线观看视频| av欧美777| 波多野结衣一区麻豆| 亚洲精品久久午夜乱码| 国产亚洲精品一区二区www | 老汉色av国产亚洲站长工具| 18禁裸乳无遮挡免费网站照片 | 99久久综合精品五月天人人| 久久精品国产综合久久久| 久久国产乱子伦精品免费另类| 极品教师在线免费播放| 国产成人欧美在线观看 | 一级a爱片免费观看的视频| 国产精品.久久久| 久9热在线精品视频| 国产精品一区二区免费欧美| 激情视频va一区二区三区| 国产欧美日韩一区二区三| 亚洲av欧美aⅴ国产| 一级片'在线观看视频| 国产免费男女视频| 天堂动漫精品| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜免费鲁丝| 自线自在国产av| 人人妻人人澡人人爽人人夜夜| 国产欧美日韩综合在线一区二区| 国产在视频线精品| 啦啦啦免费观看视频1| 色精品久久人妻99蜜桃| 夜夜躁狠狠躁天天躁| 欧美 日韩 精品 国产| 国产精品自产拍在线观看55亚洲 | 在线十欧美十亚洲十日本专区| 夜夜躁狠狠躁天天躁| 久久婷婷成人综合色麻豆| 久久久水蜜桃国产精品网| 男人操女人黄网站| 中文字幕av电影在线播放| 黄网站色视频无遮挡免费观看| 99re6热这里在线精品视频| 精品国产乱码久久久久久男人| 亚洲精品国产色婷婷电影| 一边摸一边做爽爽视频免费| 国产精品二区激情视频| 人人妻人人澡人人看| 一区二区三区国产精品乱码| 91成人精品电影| 99国产精品免费福利视频| 69精品国产乱码久久久| 脱女人内裤的视频| 国产精品影院久久| 久久99一区二区三区| 久久久久精品国产欧美久久久| 免费观看精品视频网站| 一级片免费观看大全| 国产精品永久免费网站| 亚洲精品中文字幕一二三四区| 高清av免费在线| 国产精品乱码一区二三区的特点 | 亚洲国产欧美网| 国产精品国产av在线观看| 嫁个100分男人电影在线观看| 最近最新中文字幕大全电影3 | 亚洲色图 男人天堂 中文字幕| 19禁男女啪啪无遮挡网站| 久久人妻福利社区极品人妻图片| 成人国语在线视频| 老司机深夜福利视频在线观看| 日韩成人在线观看一区二区三区| 中亚洲国语对白在线视频| 久久精品国产99精品国产亚洲性色 | 成年人免费黄色播放视频| 日本撒尿小便嘘嘘汇集6| 欧美成人午夜精品| 日韩熟女老妇一区二区性免费视频| 国产精品永久免费网站| 亚洲国产精品sss在线观看 | 久久ye,这里只有精品| 日韩大码丰满熟妇| 两个人看的免费小视频| 不卡一级毛片| 国产欧美日韩精品亚洲av| 黑人猛操日本美女一级片| 日本黄色视频三级网站网址 | 在线观看一区二区三区激情| 少妇裸体淫交视频免费看高清 | 侵犯人妻中文字幕一二三四区| 久久香蕉精品热| 韩国av一区二区三区四区| 成人黄色视频免费在线看| 18在线观看网站| 国产成人影院久久av| 极品人妻少妇av视频| 多毛熟女@视频| 黄色怎么调成土黄色| 美女扒开内裤让男人捅视频| 亚洲五月天丁香| av有码第一页| 十八禁网站免费在线| 亚洲片人在线观看| 一级毛片女人18水好多| 777米奇影视久久| 亚洲午夜理论影院| 男人舔女人的私密视频| 国产蜜桃级精品一区二区三区 | 国产成人啪精品午夜网站| 精品少妇久久久久久888优播| 777米奇影视久久| 看免费av毛片| 亚洲国产欧美网| 人妻丰满熟妇av一区二区三区 | 操美女的视频在线观看| 成人18禁在线播放| 国产亚洲一区二区精品| 国产真人三级小视频在线观看| 亚洲人成77777在线视频| 久久久久视频综合| 亚洲精品国产色婷婷电影| 一区二区三区激情视频| 在线观看免费视频日本深夜| 99国产精品免费福利视频| 热99re8久久精品国产| 很黄的视频免费| 女人爽到高潮嗷嗷叫在线视频| 国产淫语在线视频| 久久香蕉激情| 国产在视频线精品| 欧美午夜高清在线| 久久久久久久久久久久大奶| 最近最新中文字幕大全免费视频| 欧美日韩亚洲国产一区二区在线观看 | 午夜福利欧美成人| 国产午夜精品久久久久久| 亚洲成国产人片在线观看| 午夜激情av网站| 高清视频免费观看一区二区| 精品无人区乱码1区二区| 国产成人欧美| 一级片免费观看大全| 免费日韩欧美在线观看| 亚洲中文字幕日韩| 亚洲精品成人av观看孕妇| 欧美 日韩 精品 国产| 欧美色视频一区免费| 亚洲人成伊人成综合网2020| 黄片播放在线免费| 免费观看人在逋| 午夜精品久久久久久毛片777| 欧美人与性动交α欧美精品济南到| 狂野欧美激情性xxxx| 黑人操中国人逼视频| 啦啦啦在线免费观看视频4| 国产乱人伦免费视频| 悠悠久久av| 好男人电影高清在线观看| 黄片小视频在线播放| 天堂√8在线中文| 成年女人毛片免费观看观看9 | 国产极品粉嫩免费观看在线| 色老头精品视频在线观看| 乱人伦中国视频| 欧美日本中文国产一区发布| 久热这里只有精品99| 视频区欧美日本亚洲| 欧洲精品卡2卡3卡4卡5卡区| 好男人电影高清在线观看| 好男人电影高清在线观看| 国产精品一区二区在线观看99| 99久久综合精品五月天人人| 亚洲欧美激情在线| 狠狠婷婷综合久久久久久88av| 最近最新中文字幕大全电影3 | 国产免费男女视频| 亚洲国产看品久久| 亚洲精品美女久久av网站| 女人被狂操c到高潮| 国产精品亚洲一级av第二区| 精品一区二区三区视频在线观看免费 | 亚洲久久久国产精品| 99热网站在线观看| 妹子高潮喷水视频| 中文字幕人妻丝袜制服| 国产人伦9x9x在线观看| 日本五十路高清| 18禁观看日本| av有码第一页| av超薄肉色丝袜交足视频| 老司机午夜福利在线观看视频| 亚洲欧美激情在线| 丝袜美腿诱惑在线| 欧美精品av麻豆av| 久久婷婷成人综合色麻豆| 亚洲av成人不卡在线观看播放网| 他把我摸到了高潮在线观看| 制服诱惑二区| 老司机在亚洲福利影院| 极品教师在线免费播放| 一区在线观看完整版| 久久精品亚洲精品国产色婷小说| 另类亚洲欧美激情| 精品无人区乱码1区二区| 看黄色毛片网站| 精品视频人人做人人爽| 日韩有码中文字幕| 757午夜福利合集在线观看| 亚洲男人天堂网一区| 十分钟在线观看高清视频www| 亚洲精品在线美女| 亚洲精品一卡2卡三卡4卡5卡| 久久国产精品大桥未久av| 亚洲人成77777在线视频| 亚洲av日韩在线播放| 大香蕉久久成人网| 69av精品久久久久久| 婷婷丁香在线五月| 12—13女人毛片做爰片一| 亚洲av熟女| 免费在线观看黄色视频的| 午夜成年电影在线免费观看| 国产精品影院久久| 欧美日韩视频精品一区| 亚洲中文av在线| 在线观看免费午夜福利视频| 欧美人与性动交α欧美软件| 午夜精品国产一区二区电影| 午夜福利在线免费观看网站| 久久性视频一级片| 亚洲午夜精品一区,二区,三区| 国产精品自产拍在线观看55亚洲 | 啦啦啦视频在线资源免费观看| 欧美日韩黄片免| 俄罗斯特黄特色一大片| 国产精品一区二区在线不卡| 日韩成人在线观看一区二区三区| av超薄肉色丝袜交足视频| 久久国产精品人妻蜜桃| 精品第一国产精品| 建设人人有责人人尽责人人享有的| 波多野结衣一区麻豆| 欧美国产精品一级二级三级| 免费观看a级毛片全部| 久久中文字幕一级| 很黄的视频免费| 老司机深夜福利视频在线观看| av欧美777| tube8黄色片| 午夜精品国产一区二区电影| 在线观看午夜福利视频| 精品国产乱子伦一区二区三区| av网站在线播放免费| 一区二区三区国产精品乱码| 少妇 在线观看| 国产欧美日韩一区二区三区在线| 精品福利永久在线观看| 久久精品国产综合久久久| 亚洲熟妇中文字幕五十中出 | 1024香蕉在线观看| 中出人妻视频一区二区| 久久久国产成人精品二区 | 一边摸一边抽搐一进一小说 | 国产成人av激情在线播放| 香蕉久久夜色| 亚洲第一青青草原| 母亲3免费完整高清在线观看| 啦啦啦视频在线资源免费观看| 丰满的人妻完整版| 亚洲人成电影观看| 欧美乱码精品一区二区三区| 少妇猛男粗大的猛烈进出视频| 97人妻天天添夜夜摸| 午夜福利,免费看| 亚洲人成伊人成综合网2020| 亚洲精品一卡2卡三卡4卡5卡| 大陆偷拍与自拍| 亚洲专区字幕在线| 久久国产精品大桥未久av| av国产精品久久久久影院| 国产深夜福利视频在线观看| 在线观看午夜福利视频| 法律面前人人平等表现在哪些方面| 国产精品乱码一区二三区的特点 | 咕卡用的链子| 亚洲少妇的诱惑av| 满18在线观看网站| 精品欧美一区二区三区在线| 777久久人妻少妇嫩草av网站| 欧美激情极品国产一区二区三区| 亚洲七黄色美女视频| 性少妇av在线| 亚洲色图av天堂| 欧美一级毛片孕妇| 成年版毛片免费区| 日韩欧美一区视频在线观看| 亚洲一区二区三区欧美精品| 黑人猛操日本美女一级片| 亚洲va日本ⅴa欧美va伊人久久| 男人舔女人的私密视频| 露出奶头的视频| 性色av乱码一区二区三区2| 18禁美女被吸乳视频| 亚洲国产精品一区二区三区在线| 两个人看的免费小视频| 欧美日韩精品网址| 亚洲欧美激情在线| 99国产精品免费福利视频| 国产精品免费大片| 亚洲少妇的诱惑av| 捣出白浆h1v1| 成人影院久久| 窝窝影院91人妻| 亚洲精品美女久久久久99蜜臀| 丰满人妻熟妇乱又伦精品不卡| 啦啦啦 在线观看视频| 色婷婷av一区二区三区视频| 嫩草影视91久久| av福利片在线| 成年动漫av网址| 色综合欧美亚洲国产小说| 婷婷精品国产亚洲av在线 | 精品国产一区二区久久| 丁香欧美五月| 天天操日日干夜夜撸| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品美女特级片免费视频播放器 | 国产av又大| 国产99久久九九免费精品| ponron亚洲| 69精品国产乱码久久久| 免费一级毛片在线播放高清视频 | 精品一区二区三区四区五区乱码| 中文字幕av电影在线播放| 亚洲一区二区三区不卡视频| 精品一品国产午夜福利视频| 久久久久国内视频| 又紧又爽又黄一区二区| 午夜福利欧美成人| 桃红色精品国产亚洲av| 欧美精品啪啪一区二区三区| 久久精品国产亚洲av高清一级| 国产在视频线精品| 欧美亚洲日本最大视频资源| 一区在线观看完整版| 精品无人区乱码1区二区| 一边摸一边抽搐一进一出视频| 极品教师在线免费播放| 久久ye,这里只有精品| 一区二区日韩欧美中文字幕| 中国美女看黄片| 精品久久久久久久久久免费视频 | 国产色视频综合| 9热在线视频观看99| 成人亚洲精品一区在线观看| 日韩有码中文字幕| 乱人伦中国视频| 日本撒尿小便嘘嘘汇集6| 91精品国产国语对白视频| 久久精品人人爽人人爽视色| 国产成人啪精品午夜网站| 欧美精品亚洲一区二区| 大香蕉久久网| 亚洲av成人一区二区三| 日韩视频一区二区在线观看| 日日爽夜夜爽网站| 婷婷精品国产亚洲av在线 | 国内久久婷婷六月综合欲色啪| 久久 成人 亚洲| 午夜福利在线观看吧| 亚洲情色 制服丝袜| 国产主播在线观看一区二区| 欧美精品av麻豆av| 亚洲国产看品久久| 美女扒开内裤让男人捅视频| 乱人伦中国视频| 男女下面插进去视频免费观看| 国产亚洲精品一区二区www | 欧美精品人与动牲交sv欧美| 美女 人体艺术 gogo| 久久精品国产a三级三级三级| 精品久久久久久电影网| 在线国产一区二区在线| 久久精品国产亚洲av香蕉五月 | 一a级毛片在线观看| 欧美人与性动交α欧美精品济南到| 建设人人有责人人尽责人人享有的| 国产精品偷伦视频观看了| 最新美女视频免费是黄的| 18禁国产床啪视频网站| 国产高清国产精品国产三级| 久久人人97超碰香蕉20202| 在线观看午夜福利视频| 不卡一级毛片| 高清毛片免费观看视频网站 | 天天躁日日躁夜夜躁夜夜| 欧美黑人精品巨大| 欧美午夜高清在线| 巨乳人妻的诱惑在线观看| 国产成人影院久久av| 丁香六月欧美| 伦理电影免费视频| 国产亚洲av高清不卡| 两个人免费观看高清视频| 欧美黄色淫秽网站| 淫妇啪啪啪对白视频| 久热爱精品视频在线9| 精品亚洲成国产av| 91成人精品电影| 免费看十八禁软件| 亚洲精品国产区一区二| 中文字幕人妻熟女乱码| 午夜日韩欧美国产| 久久亚洲精品不卡| 深夜精品福利| 大码成人一级视频| 91国产中文字幕| 又黄又爽又免费观看的视频| 久久精品人人爽人人爽视色| 欧美黑人精品巨大| 国产av又大|