本發(fā)明涉及一種圖像水印方法,尤其是涉及一種對(duì)抗色調(diào)映射的高動(dòng)態(tài)范圍圖像水印方法。
背景技術(shù):
:近年來(lái),高動(dòng)態(tài)范圍(highdynamicrange,hdr)圖像作為數(shù)字圖像領(lǐng)域的重要發(fā)展和突破,已在數(shù)碼攝影、超高清電影與電視、視頻游戲、遙感探測(cè)和醫(yī)學(xué)成像等各方面日益受到關(guān)注。不同于傳統(tǒng)的低動(dòng)態(tài)范圍(lowdynamicrange,ldr)圖像,高動(dòng)態(tài)范圍圖像采用浮點(diǎn)數(shù)據(jù)來(lái)記錄像素信息,能夠更加準(zhǔn)確記錄真實(shí)場(chǎng)景的全部色彩范圍值,能夠表現(xiàn)出豐富的色彩細(xì)節(jié)和明暗層次。低動(dòng)態(tài)范圍圖像能夠直接在常規(guī)顯示設(shè)備上再現(xiàn),而高動(dòng)態(tài)范圍圖像只能在專用顯示設(shè)備上再現(xiàn),但是由于專用顯示設(shè)備的研制過(guò)程復(fù)雜、成本較高、技術(shù)還不夠成熟而普及面有限,因此為了能夠在常規(guī)顯示設(shè)備上再現(xiàn)高動(dòng)態(tài)范圍圖像豐富的顏色和階調(diào)細(xì)節(jié),研究人員設(shè)計(jì)了各種能將高動(dòng)態(tài)范圍圖像轉(zhuǎn)化成低動(dòng)態(tài)范圍圖像的色調(diào)映射算子(tonemappingoperator,tmo)。由于數(shù)字圖像本身具有易復(fù)制、易傳播、易篡改的特點(diǎn),因此高動(dòng)態(tài)范圍圖像的知識(shí)產(chǎn)權(quán)保護(hù)問(wèn)題亟待解決。雖然低動(dòng)態(tài)范圍圖像數(shù)字水印技術(shù)已經(jīng)相當(dāng)成熟,但由于高動(dòng)態(tài)范圍圖像獨(dú)特的存儲(chǔ)方式和色調(diào)映射算子的存在,現(xiàn)有的低動(dòng)態(tài)范圍圖像數(shù)字水印技術(shù)不能直接移植到高動(dòng)態(tài)范圍圖像領(lǐng)域,因此設(shè)計(jì)一種對(duì)抗色調(diào)映射算子的魯棒性高動(dòng)態(tài)范圍圖像水印算法顯得尤為重要。近年來(lái),高動(dòng)態(tài)范圍圖像水印技術(shù)的研究日益得到重視。c.m.yu,k.c.wu,c.m.wang,"adistortion-freedatahidingschemeforhighdynamicrangeimages",displays32(5),pp.225-236,2011.(余、吳等人,一種無(wú)失真的高動(dòng)態(tài)范圍圖像信息隱藏方案,顯示器32(5),225-236,2011)、z.h.wang,c.c.chang,t.y.lin,c.c.lin,"anoveldistortion-freedatahidingschemeforhighdynamicrangeimages",fourthinternationalconferenceondigitalhome,guangzhou,china,pp.33-38,2012.(王、常等人,一種新型無(wú)失真的高動(dòng)態(tài)范圍圖像信息隱藏方案,第四屆數(shù)字之家國(guó)際會(huì)議,中國(guó)廣州,33-38,2012)、c.c.chang,t.s.nguyen,c.c.lin,"distortion-freedataembeddingschemeforhighdynamicrangeimages",journalofelectronicscienceandtechnology11(1),pp.20-26,2013.(常、nguyen等人,一種無(wú)失真的高動(dòng)態(tài)范圍圖像數(shù)據(jù)嵌入方案,電子科學(xué)技術(shù)11(1),20-26,2013)、c.c.chang,t.s.nguyen,c.c.lin,"anewdistortion-freedataembeddingschemeforhigh-dynamicrangeimages",multimediatoolsandapplications.75(1),pp.145-163,2016.(常、nguyen等人,一種新型的無(wú)失真高動(dòng)態(tài)范圍圖像數(shù)據(jù)嵌入方案,多媒體工具及應(yīng)用,75(1),145-163,2016),這些方法利用了rgbe格式的高動(dòng)態(tài)范圍圖像的存儲(chǔ)特點(diǎn),即rgb三通道和e通道存在多種同樣的映射關(guān)系,實(shí)現(xiàn)了水印的無(wú)損嵌入。y.m.cheng,c.m.wang,"anovelapproachtosteganographyinhighdynamicrangeimages".ieeemultimedia16(3),pp.70-80,2009.(陳、王等人,一種新型的高動(dòng)態(tài)范圍圖像數(shù)據(jù)隱寫(xiě)方法,ieee多媒體16(3),70-80,2009)、m.t.li,n.c.huang,c.m.wang,"adatahidingschemeforhdrimages".internationaljournalofinnovativecomputinginformationandcontrol7(5a),pp.2021-2035,2011.(李、黃等人,高動(dòng)態(tài)范圍圖像信息隱藏方案,國(guó)際創(chuàng)新計(jì)算信息與控制雜志7(5a),2021-2035,2011),這些方法采用了最低有效位(leastsignificantbits,lsb)算法在rgbe格式和logluv(tiff)格式的高動(dòng)態(tài)范圍圖像上進(jìn)行水印嵌入。y.t.lin,c.m.wang,w.s.chen,f.p.lin,w.lin,"anoveldatahidingalgorithmforhighdynamicrangeimages",ieeetransactionsonmultimedia19(1),pp.196-211,2017.(林、王等人,一種新型的高動(dòng)態(tài)范圍圖像數(shù)據(jù)隱藏算法,ieee多媒體匯刊19(1),196-211,2017),該方法利用了openexr格式的高動(dòng)態(tài)范圍圖像中尾數(shù)對(duì)圖像的質(zhì)量影響較小這一存儲(chǔ)特點(diǎn),進(jìn)行水印嵌入。上述這些方法主要依靠高動(dòng)態(tài)范圍圖像的存儲(chǔ)格式特點(diǎn),在不可見(jiàn)性和嵌入容量方面取得了一定的效果,但是在魯棒性,尤其是對(duì)抗色調(diào)映射算子方面無(wú)能為力;此外,上述各種方法均不能適用于所有格式的高動(dòng)態(tài)范圍圖像處理,存在使用局限性。現(xiàn)有的對(duì)抗色調(diào)映射算子的高動(dòng)態(tài)范圍圖像水印算法主要有以下幾種。f.guerrini,m.okuda,n.adami,r.leonardi,"highdynamicrangeimagewatermarkingrobustagainsttone-mappingoperators".ieeetrans.informationforensicsandsecurity6(2),pp.283-295,2011.(guerrini、okuda等人,對(duì)抗色調(diào)映射的魯棒性高動(dòng)態(tài)范圍圖像水印,ieee信息取證和安全匯刊6(2),283-295,2011),該方法在小波變換的低頻區(qū)域上,對(duì)低感知區(qū)域分塊并以峰度作為特征值,利用量化索引調(diào)制(quantizationindexmodulation,qim)方法進(jìn)行水印嵌入,該方法表現(xiàn)出了良好的色調(diào)映射算子(7種)魯棒性和不可見(jiàn)性,但是水印誤碼率高于20%。x.xue,t.jinno,x.jin,m.okuda,s.goto,"watermarkingforhdrimagerobusttotonemapping".ieicetransactions94-a(11),pp.2334-2341,2011.(薛、jinno等人,針對(duì)色調(diào)映射的魯棒性高動(dòng)態(tài)范圍圖像水印,ieice匯刊94-a(11),2334-2341,2011),該方法首先對(duì)高動(dòng)態(tài)范圍圖像進(jìn)行對(duì)數(shù)預(yù)處理,然后利用雙邊濾波器生成基本層和細(xì)節(jié)層,最后在細(xì)節(jié)層利用小波變換,對(duì)橫向、縱向的高頻區(qū)域進(jìn)行水印嵌入,該方法可以有效地對(duì)抗色調(diào)映射算子(4種)攻擊,但是測(cè)試圖像數(shù)量相對(duì)較少,水印誤碼率仍然較高。e.maiorana,v.solachidis,p.campisi,y.lou,"robustmulti-bitwatermarkingforhdrimagesintheradon-dctdomain".int.symposiumonimageandsignalprocessingandanalysis,trieste,italy,pp.284-289,2013.(maiorana、solachidis等人,radon-dct域的魯棒性多位高動(dòng)態(tài)范圍圖像水印,圖像信號(hào)處理和分析國(guó)際專題會(huì)議,意大利的里雅斯特,284-289,2013),該方法對(duì)高動(dòng)態(tài)范圍圖像進(jìn)行對(duì)數(shù)預(yù)處理后,進(jìn)行離散小波變換和radon-dct變換,采用量化索引調(diào)制方法在最大能量區(qū)域嵌入水印,該方法的不可見(jiàn)性較好,但水印誤碼率達(dá)到22%。v.solachidis,e.maiorana,p.campisi,"hdrimagemulti-bitwatermarkingusingbilateral-filtering-basedmasking".imageprocessing:algorithmsandsystems,vol.865505,2013.(solachidis、maiorana等人,基于雙邊濾波器掩膜的多位高動(dòng)態(tài)范圍圖像水印,圖像處理:算法和系統(tǒng),卷號(hào)865505,2013),該方法利用恰可感知失真模型和對(duì)比敏感度函數(shù)設(shè)計(jì)了感知掩模,小波變換后嵌入水印,該方法在7種色調(diào)映射算子上表現(xiàn)出了較好的魯棒性,不足之處在于嵌入量只有128bit。v.solachidis,e.maiorana,p.campisi,f.banterle,"hdrimagewatermarkingbasedonbracketingdecomposition".int.conf.ondigitalsignalprocessing,santorini,greece,pp.1-6,2013.(solachidis、maiorana等人,基于曝光分解的高動(dòng)態(tài)范圍圖像水印,數(shù)字信號(hào)處理國(guó)際會(huì)議,希臘圣托里尼,1-6,2013.),該方法將高動(dòng)態(tài)范圍圖像分解為一系列不同亮度范圍的低動(dòng)態(tài)范圍圖像,并分別進(jìn)行小波變換嵌入水印,這種多次嵌入造成了該方法的不可見(jiàn)性較差。j.l.wu,"robustwatermarkingframeworkforhighdynamicrangeimagesagainsttone-mappingattacks".watermarking2,229-242,2012.(吳,對(duì)抗色調(diào)映射攻擊的魯棒性高動(dòng)態(tài)范圍圖像水印結(jié)構(gòu),水印2,229-242,2012),該方法首先選用一種特定的色調(diào)映射算子對(duì)高動(dòng)態(tài)范圍圖像進(jìn)行預(yù)處理,對(duì)得到的低動(dòng)態(tài)范圍圖像在dct域進(jìn)行水印嵌入,然后重建為高動(dòng)態(tài)范圍圖像,該方法對(duì)色調(diào)映射算子的魯棒性較好,但需要預(yù)先確定色調(diào)映射算子,顯然這并不現(xiàn)實(shí)。以上對(duì)抗色調(diào)映射算子的高動(dòng)態(tài)范圍圖像水印算法從不同角度出發(fā),利用高動(dòng)態(tài)范圍圖像的特性,在對(duì)抗少數(shù)色調(diào)映射算子方面取得了一定的效果,然而,以上算法并未從根源上找出高動(dòng)態(tài)范圍圖像水印和低動(dòng)態(tài)范圍圖像水印的不同之處,特別是色調(diào)映射算子對(duì)高動(dòng)態(tài)范圍圖像水印的影響方面,分析不夠深入和準(zhǔn)確,因而在對(duì)抗色調(diào)映射算子的數(shù)量和水印誤碼率方面仍有待進(jìn)一步提高。技術(shù)實(shí)現(xiàn)要素:本發(fā)明所要解決的技術(shù)問(wèn)題是提供一種對(duì)抗色調(diào)映射的高動(dòng)態(tài)范圍圖像水印方法,其能夠在高動(dòng)態(tài)范圍圖像中嵌入水印信息,可以在不可見(jiàn)性和嵌入容量之間進(jìn)行權(quán)衡,靈活性和實(shí)用性高,且能夠在受常規(guī)攻擊后得到的高動(dòng)態(tài)范圍圖像或在受色調(diào)映射攻擊后得到的低動(dòng)態(tài)范圍圖像中提取出水印信息,能夠?qū)顾猩{(diào)映射攻擊,魯棒性好、水印誤碼率低。本發(fā)明解決上述技術(shù)問(wèn)題所采用的技術(shù)方案為:一種對(duì)抗色調(diào)映射的高動(dòng)態(tài)范圍圖像水印方法,其特征在于包括水印嵌入和水印提取兩部分;所述的水印嵌入部分的具體步驟為:①_1、將待嵌入的原始二值水印序列記為w;將待嵌入水印的rgbe、openexr或logluv格式的原始高動(dòng)態(tài)范圍圖像記為iorg,iorg的高度為rorg,iorg的寬度為corg;并人為設(shè)定iorg嵌入水印后得到的高動(dòng)態(tài)范圍圖像相對(duì)于iorg的圖像質(zhì)量值,記為q;其中,q∈[70,100];①_2、將iorg從rgb空間轉(zhuǎn)化到hsv空間,得到iorg的h層圖像、s層圖像和v層圖像,對(duì)應(yīng)記為iorg-h、iorg-s和iorg-v;然后判斷rorg和corg是否能夠被4整除,如果rorg和corg均能夠被4整除,則對(duì)iorg-v不進(jìn)行像素點(diǎn)填充,并將iorg-v重新記為iorg-v-new,iorg-v-new的高度為rnew,iorg-v-new的寬度為cnew,rnew=rorg,cnew=corg;如果rorg能夠被4整除,而corg不能夠被4整除,則在iorg-v的最右側(cè)填充4-mod(corg,4)列像素點(diǎn),且令填充的所有像素點(diǎn)的像素值為0,并將填充后得到的圖像記為iorg-v-new,iorg-v-new的高度為rnew,iorg-v-new的寬度為cnew,rnew=rorg,cnew=4-mod(corg,4)+corg;如果rorg不能夠被4整除,而corg能夠被4整除,則在iorg-v的最下側(cè)填充4-mod(rorg,4)行像素點(diǎn),且令填充的所有像素點(diǎn)的像素值為0,并將填充后得到的圖像記為iorg-v-new,iorg-v-new的高度為rnew,iorg-v-new的寬度為cnew,rnew=4-mod(rorg,4)+rorg,cnew=corg;如果rorg和corg均不能夠被4整除,則在iorg-v的最右側(cè)填充4-mod(corg,4)列像素點(diǎn),在iorg-v的最下側(cè)填充4-mod(rorg,4)行像素點(diǎn),且令填充的所有像素點(diǎn)的像素值為0,并將填充后得到的圖像記為iorg-v-new,iorg-v-new的高度為rnew,iorg-v-new的寬度為cnew,rnew=4-mod(rorg,4)+rorg,cnew=4-mod(corg,4)+corg;其中,rnew=rorg和cnew=corg中的“=”為賦值符號(hào),mod()為取余函數(shù);①_3、對(duì)iorg-v-new進(jìn)行一級(jí)冗余小波變換,得到iorg-v-new的低頻近似分量、水平方向高頻分量、垂直方向高頻分量和對(duì)角線方向高頻分量,對(duì)應(yīng)記為llorg-v-new、hlorg-v-new、lhorg-v-new和hhorg-v-new;然后將llorg-v-new劃分成個(gè)互不重疊的尺寸大小為4×4的分塊,將llorg-v-new中坐標(biāo)位置為(x,y)的分塊記為bll,org(x,y);其中,①_4、對(duì)iorg-v-new中的每個(gè)像素點(diǎn)的像素值進(jìn)行對(duì)數(shù)歸一化處理,得到對(duì)數(shù)歸一化處理后的圖像,記為i'org-v-new;然后將i'org-v-new劃分成個(gè)互不重疊的尺寸大小為4×4的分塊,將i'org-v-new中坐標(biāo)位置為(x,y)的分塊記為borg-v-new(x,y);再對(duì)i'org-v-new進(jìn)行下采樣處理,得到下采樣處理后的圖像,記為valueorg-v-new,valueorg-v-new的高度為valueorg-v-new的寬度為將valueorg-v-new中坐標(biāo)位置為(x,y)的像素點(diǎn)的像素值記為valueorg-v-new(x,y),valueorg-v-new(x,y)的值等于borg-v-new(x,y)中的所有像素點(diǎn)的像素值中的最小值;①_5、對(duì)valueorg-v-new進(jìn)行二值化處理,得到iorg-v-new的感知掩膜,記為markorg-v-new,將markorg-v-new中坐標(biāo)位置為(x,y)的像素點(diǎn)的像素值記為markorg-v-new(x,y),然后根據(jù)markorg-v-new中的每個(gè)像素點(diǎn)的像素值,判定llorg-v-new中對(duì)應(yīng)坐標(biāo)位置的分塊是否嵌入水印,對(duì)于bll,org(x,y),如果markorg-v-new(x,y)=0,則判定bll,org(x,y)不需嵌入水?。蝗绻鹠arkorg-v-new(x,y)=1,則判定bll,org(x,y)需嵌入水?。黄渲?,β表示掩膜亮度閾值,β∈[0,1],β的初始值為1;①_6、對(duì)llorg-v-new中判定為需嵌入水印的每個(gè)分塊進(jìn)行奇異值分解變換,得到llorg-v-new中判定為需嵌入水印的每個(gè)分塊的維數(shù)為4×4的u矩陣分量、維數(shù)為4×4的s矩陣分量、維數(shù)為4×4的v矩陣分量;然后在llorg-v-new中判定為需嵌入水印的每個(gè)分塊的u矩陣分量中的第2行第1列元素和第4行第1列元素中嵌入水印,假設(shè)bll,org(x,y)判定為需嵌入水印,則在bll,org(x,y)的u矩陣分量中的第2行第1列元素u21和第4行第1列元素u41中嵌入水印,如果從w中讀取的當(dāng)前待嵌入的二值水印為1,則當(dāng)u21-u41<0且|u21-u41|>t時(shí),對(duì)u21和u41不進(jìn)行修改,至此完成bll,org(x,y)的水印嵌入;當(dāng)|u21-u41|<t時(shí),對(duì)u21和u41進(jìn)行修改,令并令至此完成bll,org(x,y)的水印嵌入;如果從w中讀取的當(dāng)前待嵌入的二值水印為0,則當(dāng)u21-u41>0且|u21-u41|>t時(shí),對(duì)u21和u41不進(jìn)行修改,至此完成bll,org(x,y)的水印嵌入;當(dāng)|u21-u41|<t時(shí),對(duì)u21和u41進(jìn)行修改,令并令至此完成bll,org(x,y)的水印嵌入;其中,符號(hào)“||”為取絕對(duì)值符號(hào),t表示水印嵌入閾值,sign()為取符號(hào)函數(shù),和中的“=”為賦值符號(hào);①_7、對(duì)嵌有水印的每個(gè)分塊進(jìn)行奇異值分解逆變換,重建低頻近似分量,將重建得到的低頻近似分量記為ll*org-v-new;然后將ll*org-v-new、hlorg-v-new、lhorg-v-new和hhorg-v-new合成完整的嵌有水印的圖像;接著對(duì)嵌有水印的圖像進(jìn)行一級(jí)冗余小波變換的逆變換,得到逆變換后的嵌有水印的圖像,記為i*org-v-new;之后丟棄i*org-v-new的最右側(cè)的4-mod(corg,4)列像素點(diǎn)和最下側(cè)的4-mod(rorg,4)行像素點(diǎn),得到嵌有水印的v層圖像,記為i*org-v;再將iorg-h、iorg-s和i*org-v合成完整的嵌有水印的hsv空間圖像,記為i*hsv;最后將i*hsv從hsv空間轉(zhuǎn)化到rgb空間,得到嵌有水印的rgb空間圖像,記為i*org;①_8、利用hdr-vdp-2方法,計(jì)算i*org相對(duì)于iorg的圖像質(zhì)量值,記為q*;然后判斷q*>q或|q*-q|<6是否成立,如果成立,則將markorg-v-new確定為水印秘鑰k,將i*org確定為最終的嵌有水印的高動(dòng)態(tài)范圍圖像,水印嵌入過(guò)程結(jié)束;否則,令β=β-δstep,然后返回步驟①_5繼續(xù)執(zhí)行;其中,β=β-δstep中的“=”為賦值符號(hào),δstep表示調(diào)整步長(zhǎng);所述的水印提取部分的具體步驟為:②_1、讀取嵌有水印的高動(dòng)態(tài)范圍圖像受常規(guī)攻擊后得到的高動(dòng)態(tài)范圍圖像或受色調(diào)映射后得到的低動(dòng)態(tài)范圍圖像,記為istego,istego的高度為rorg,istego的寬度為corg;②_2、將istego從rgb空間轉(zhuǎn)化到hsv空間,得到istego的h層圖像、s層圖像和v層圖像,對(duì)應(yīng)記為istego-h、istego-s和istego-v;然后判斷rorg和corg是否能夠被4整除,如果rorg和corg均能夠被4整除,則對(duì)istego-v不進(jìn)行像素點(diǎn)填充,并將istego-v重新記為istego-v-new,istego-v-new的高度為rnew,istego-v-new的寬度為cnew,rnew=rorg,cnew=corg;如果rorg能夠被4整除,而corg不能夠被4整除,則在istego-v的最右側(cè)填充4-mod(corg,4)列像素點(diǎn),且令填充的所有像素點(diǎn)的像素值為0,并將填充后得到的圖像記為istego-v-new,istego-v-new的高度為rnew,istego-v-new的寬度為cnew,rnew=rorg,cnew=4-mod(corg,4)+corg;如果rorg不能夠被4整除,而corg能夠被4整除,則在istego-v的最下側(cè)填充4-mod(rorg,4)行像素點(diǎn),且令填充的所有像素點(diǎn)的像素值為0,并將填充后得到的圖像記為istego-v-new,istego-v-new的高度為rnew,istego-v-new的寬度為cnew,rnew=4-mod(rorg,4)+rorg,cnew=corg;如果rorg和corg均不能夠被4整除,則在istego-v的最右側(cè)填充4-mod(corg,4)列像素點(diǎn),在istego-v的最下側(cè)填充4-mod(rorg,4)行像素點(diǎn),且令填充的所有像素點(diǎn)的像素值為0,并將填充后得到的圖像記為istego-v-new,istego-v-new的高度為rnew,istego-v-new的寬度為cnew,rnew=4-mod(rorg,4)+rorg,cnew=4-mod(corg,4)+corg;②_3、對(duì)istego-v-new進(jìn)行一級(jí)冗余小波變換,得到istego-v-new的低頻近似分量、水平方向高頻分量、垂直方向高頻分量和對(duì)角線方向高頻分量,對(duì)應(yīng)記為llstego-v-new、hlstego-v-new、lhstego-v-new和hhstego-v-new;然后將llstego-v-new劃分成個(gè)互不重疊的尺寸大小為4×4的分塊,將llstego-v-new中坐標(biāo)位置為(x,y)的分塊記為bll,stego(x,y);其中,②_4、讀取水印秘鑰k;然后根據(jù)水印秘鑰k中的每個(gè)像素點(diǎn)的像素值,判定llstego-v-new中對(duì)應(yīng)坐標(biāo)位置的分塊是否提取水印,對(duì)于bll,stego(x,y),如果markorg-v-new(x,y)=0,則判定bll,stego(x,y)不需提取水??;如果markorg-v-new(x,y)=1,則判定bll,stego(x,y)需提取水??;②_5、對(duì)llstego-v-new中判定為需提取水印的每個(gè)分塊進(jìn)行奇異值分解變換,得到llstego-v-new中判定為需提取水印的每個(gè)分塊的維數(shù)為4×4的u矩陣分量、維數(shù)為4×4的s矩陣分量、維數(shù)為4×4的v矩陣分量;然后在llstego-v-new中判定為需提取水印的每個(gè)分塊的u矩陣分量中的第2行第1列元素和第4行第1列元素中提取水印,假設(shè)bll,stego(x,y)判定為需提取水印,則在bll,stego(x,y)的u矩陣分量中的第2行第1列元素u*21和第4行第1列元素u*41中提取水印,如果u*21>u*41,則提取出的二值水印為0;如果u*21≤u*41,則提取出的二值水印為1;再根據(jù)提取出的所有二值水印,得到提取出的二值水印序列,記為w*。所述的步驟①_6中取t=0.08。所述的步驟①_8中取δstep=0.01。與現(xiàn)有技術(shù)相比,本發(fā)明的優(yōu)點(diǎn)在于:1)本發(fā)明方法不依賴于高動(dòng)態(tài)范圍圖像的存儲(chǔ)格式,將現(xiàn)有的rgbe、openexr或logluv(tiff)三種公知格式的高動(dòng)態(tài)范圍圖像從rgb空間轉(zhuǎn)化為hsv空間,直接在原始的浮點(diǎn)數(shù)值上進(jìn)行水印嵌入,因而,現(xiàn)有的rgbe、openexr和logluv(tiff)三種公知存儲(chǔ)格式之間的讀寫(xiě)和相互轉(zhuǎn)化,對(duì)水印的嵌入和提取沒(méi)有任何影響,因而本發(fā)明方法具有極高的通用性。2)本發(fā)明方法為對(duì)抗色調(diào)映射的高動(dòng)態(tài)范圍圖像水印方法,對(duì)所有色調(diào)映射算子都具有較高的魯棒性,即從嵌有水印的高動(dòng)態(tài)范圍圖像受常規(guī)攻擊后得到的高動(dòng)態(tài)范圍圖像或受色調(diào)映射攻擊后得到的低動(dòng)態(tài)范圍圖像中,都能對(duì)水印信息進(jìn)行有效的提取,滿足高動(dòng)態(tài)范圍圖像數(shù)字水印的魯棒性要求,因而具有極高的實(shí)用性。3)本發(fā)明方法在水印嵌入過(guò)程中,人工設(shè)定嵌入水印的高動(dòng)態(tài)范圍圖像相對(duì)于原始高動(dòng)態(tài)范圍圖像的圖像質(zhì)量值,這樣用戶可以在不可見(jiàn)性和嵌入容量之間進(jìn)行權(quán)衡,具有較高的靈活性和實(shí)用性。4)本發(fā)明方法在水印嵌入過(guò)程中,獲取感知掩膜的算法復(fù)雜度極低,可以根據(jù)人工設(shè)定的圖像質(zhì)量值實(shí)時(shí)進(jìn)行更新;同時(shí),感知掩膜的尺寸近似為原始高動(dòng)態(tài)范圍圖像的尺寸的1/16,且為二值圖像,所需存儲(chǔ)空間小,便于作為水印秘鑰使用,也增強(qiáng)了高動(dòng)態(tài)范圍圖像數(shù)字水印的安全性。5)本發(fā)明方法在水印嵌入過(guò)程中,利用了冗余小波變換的擴(kuò)頻特性和冗余性增強(qiáng)了本發(fā)明方法的魯棒性和嵌入容量;同時(shí),利用奇異值分解修改圖像的結(jié)構(gòu)信息增強(qiáng)了本發(fā)明方法的嵌入強(qiáng)度;最后,設(shè)計(jì)了低復(fù)雜度的亮度感知掩模自適應(yīng)調(diào)節(jié)水印的不可見(jiàn)性;在水印提取過(guò)程中,可以從嵌有水印的高動(dòng)態(tài)范圍圖像受常規(guī)攻擊后得到的高動(dòng)態(tài)范圍圖像或受色調(diào)映射攻擊后得到的低動(dòng)態(tài)范圍圖像中,對(duì)水印信息進(jìn)行有效的提取,能夠?qū)顾猩{(diào)映射攻擊,魯棒性好、水印誤碼率低。附圖說(shuō)明圖1為本發(fā)明方法的水印嵌入部分的流程框圖;圖2為本發(fā)明方法的水印提取部分的流程框圖。具體實(shí)施方式以下結(jié)合附圖實(shí)施例對(duì)本發(fā)明作進(jìn)一步詳細(xì)描述。本發(fā)明提出的一種對(duì)抗色調(diào)映射的高動(dòng)態(tài)范圍圖像水印方法,其包括水印嵌入和水印提取兩部分。所述的水印嵌入部分的流程框圖如圖1所示,其具體步驟為:①_1、將待嵌入的原始二值水印序列記為w;將待嵌入水印的rgbe、openexr或logluv格式的原始高動(dòng)態(tài)范圍圖像記為iorg,iorg的高度為rorg,iorg的寬度為corg;并人為設(shè)定iorg嵌入水印后得到的高動(dòng)態(tài)范圍圖像相對(duì)于iorg的圖像質(zhì)量值,記為q;其中,q∈[70,100],在本實(shí)施例中將q值設(shè)為90,w的長(zhǎng)度可隨q值進(jìn)行調(diào)節(jié),q值越大則w的長(zhǎng)度越小。①_2、高動(dòng)態(tài)范圍圖像有三種公知格式,分別為rgbe、openexr或logluv(tiff)格式,三種格式的高動(dòng)態(tài)范圍圖像均為rgb空間圖像。將iorg從rgb空間轉(zhuǎn)化到hsv空間,得到iorg的h層圖像、s層圖像和v層圖像,對(duì)應(yīng)記為iorg-h、iorg-s和iorg-v,選取iorg-v作為嵌入水印的宿主圖像;然后判斷rorg和corg是否能夠被4整除,如果rorg和corg均能夠被4整除,則對(duì)iorg-v不進(jìn)行像素點(diǎn)填充,并將iorg-v重新記為iorg-v-new,iorg-v-new的高度為rnew,iorg-v-new的寬度為cnew,rnew=rorg,cnew=corg;如果rorg能夠被4整除,而corg不能夠被4整除,則在iorg-v的最右側(cè)填充4-mod(corg,4)列像素點(diǎn),且令填充的所有像素點(diǎn)的像素值為0,并將填充后得到的圖像記為iorg-v-new,iorg-v-new的高度為rnew,iorg-v-new的寬度為cnew,rnew=rorg,cnew=4-mod(corg,4)+corg;如果rorg不能夠被4整除,而corg能夠被4整除,則在iorg-v的最下側(cè)填充4-mod(rorg,4)行像素點(diǎn),且令填充的所有像素點(diǎn)的像素值為0,并將填充后得到的圖像記為iorg-v-new,iorg-v-new的高度為rnew,iorg-v-new的寬度為cnew,rnew=4-mod(rorg,4)+rorg,cnew=corg;如果rorg和corg均不能夠被4整除,則在iorg-v的最右側(cè)填充4-mod(corg,4)列像素點(diǎn),在iorg-v的最下側(cè)填充4-mod(rorg,4)行像素點(diǎn),且令填充的所有像素點(diǎn)的像素值為0,并將填充后得到的圖像記為iorg-v-new,iorg-v-new的高度為rnew,iorg-v-new的寬度為cnew,rnew=4-mod(rorg,4)+rorg,cnew=4-mod(corg,4)+corg;其中,rnew=rorg和cnew=corg中的“=”為賦值符號(hào),mod()為取余函數(shù)。①_3、對(duì)iorg-v-new進(jìn)行一級(jí)冗余小波變換(redundantdiscretewavelettransform,rdwt),得到iorg-v-new的低頻近似分量、水平方向高頻分量、垂直方向高頻分量和對(duì)角線方向高頻分量,對(duì)應(yīng)記為llorg-v-new、hlorg-v-new、lhorg-v-new和hhorg-v-new;然后將llorg-v-new劃分成個(gè)互不重疊的尺寸大小為4×4的分塊,將llorg-v-new中坐標(biāo)位置為(x,y)的分塊記為bll,org(x,y);其中,①_4、對(duì)iorg-v-new中的每個(gè)像素點(diǎn)的像素值進(jìn)行對(duì)數(shù)歸一化處理,得到對(duì)數(shù)歸一化處理后的圖像,記為i'org-v-new;然后將i'org-v-new劃分成個(gè)互不重疊的尺寸大小為4×4的分塊,將i'org-v-new中坐標(biāo)位置為(x,y)的分塊記為borg-v-new(x,y);再對(duì)i'org-v-new進(jìn)行下采樣處理,得到下采樣處理后的圖像,記為valueorg-v-new,valueorg-v-new的高度為valueorg-v-new的寬度為將valueorg-v-new中坐標(biāo)位置為(x,y)的像素點(diǎn)的像素值記為valueorg-v-new(x,y),valueorg-v-new(x,y)的值等于borg-v-new(x,y)中的所有像素點(diǎn)的像素值中的最小值。①_5、對(duì)valueorg-v-new進(jìn)行二值化處理,得到iorg-v-new的感知掩膜,記為markorg-v-new,將markorg-v-new中坐標(biāo)位置為(x,y)的像素點(diǎn)的像素值記為markorg-v-new(x,y),然后根據(jù)markorg-v-new中的每個(gè)像素點(diǎn)的像素值,判定llorg-v-new中對(duì)應(yīng)坐標(biāo)位置的分塊是否嵌入水印,對(duì)于bll,org(x,y),如果markorg-v-new(x,y)=0,則判定bll,org(x,y)不需嵌入水??;如果markorg-v-new(x,y)=1,則判定bll,org(x,y)需嵌入水??;其中,β表示掩膜亮度閾值,β∈[0,1],β的初始值為1。①_6、對(duì)llorg-v-new中判定為需嵌入水印的每個(gè)分塊進(jìn)行奇異值分解(singularvaluedecomposition,svd)變換,得到llorg-v-new中判定為需嵌入水印的每個(gè)分塊的維數(shù)為4×4的u矩陣分量、維數(shù)為4×4的s矩陣分量、維數(shù)為4×4的v矩陣分量。然后在llorg-v-new中判定為需嵌入水印的每個(gè)分塊的u矩陣分量中的第2行第1列元素和第4行第1列元素中嵌入水印,假設(shè)bll,org(x,y)判定為需嵌入水印,則在bll,org(x,y)的u矩陣分量中的第2行第1列元素u21和第4行第1列元素u41中嵌入水印,如果從w中讀取的當(dāng)前待嵌入的二值水印為1,則當(dāng)u21-u41<0且|u21-u41|>t時(shí),對(duì)u21和u41不進(jìn)行修改,至此完成bll,org(x,y)的水印嵌入;當(dāng)|u21-u41|<t時(shí),對(duì)u21和u41進(jìn)行修改,令并令至此完成bll,org(x,y)的水印嵌入。如果從w中讀取的當(dāng)前待嵌入的二值水印為0,則當(dāng)u21-u41>0且|u21-u41|>t時(shí),對(duì)u21和u41不進(jìn)行修改,至此完成bll,org(x,y)的水印嵌入;當(dāng)|u21-u41|<t時(shí),對(duì)u21和u41進(jìn)行修改,令并令至此完成bll,org(x,y)的水印嵌入。其中,符號(hào)“||”為取絕對(duì)值符號(hào),t表示水印嵌入閾值,t的值越大表示嵌入強(qiáng)度越高,在本實(shí)施例中取t=0.08,sign()為取符號(hào)函數(shù),和中的“=”為賦值符號(hào)。①_7、對(duì)嵌有水印的每個(gè)分塊進(jìn)行奇異值分解(svd)逆變換,重建低頻近似分量,將重建得到的低頻近似分量記為ll*org-v-new;然后將ll*org-v-new、hlorg-v-new、lhorg-v-new和hhorg-v-new合成完整的嵌有水印的圖像;接著對(duì)嵌有水印的圖像進(jìn)行一級(jí)冗余小波變換(rdwt)的逆變換,得到逆變換后的嵌有水印的圖像,記為i*org-v-new;之后丟棄i*org-v-new的最右側(cè)的4-mod(corg,4)列像素點(diǎn)和最下側(cè)的4-mod(rorg,4)行像素點(diǎn),得到嵌有水印的v層圖像,記為i*org-v;再將iorg-h、iorg-s和i*org-v合成完整的嵌有水印的hsv空間圖像,記為i*hsv;最后將i*hsv從hsv空間轉(zhuǎn)化到rgb空間,得到嵌有水印的rgb空間圖像,記為i*org。①_8、利用現(xiàn)有的hdr-vdp-2方法(參考文獻(xiàn)為:m.narwaria,r.k.mantiuk,m.p.d.silva,p.l.callet."hdr-vdp-2.2:acalibratedmethodforobjectivequalitypredictionofhighdynamicrangeandstandardimages".injournalofelectronicimaging,24(1),2015.(narwaria、mantiuk等人,hdr-vdp-2.2:一種高動(dòng)態(tài)范圍圖像和標(biāo)準(zhǔn)圖像的客觀質(zhì)量預(yù)測(cè)基準(zhǔn)方法,電子成像期刊,24(1),2015)),計(jì)算i*org相對(duì)于iorg的圖像質(zhì)量值,記為q*;然后判斷q*>q或|q*-q|<6是否成立,如果成立,則將markorg-v-new確定為水印秘鑰k,將i*org確定為最終的嵌有水印的高動(dòng)態(tài)范圍圖像,水印嵌入過(guò)程結(jié)束;否則,令β=β-δstep,然后返回步驟①_5繼續(xù)執(zhí)行;其中,β=β-δstep中的“=”為賦值符號(hào),δstep表示調(diào)整步長(zhǎng),在本實(shí)施例中取δstep=0.01。所述的水印提取部分的流程框圖如圖2所示,其具體步驟為:②_1、讀取嵌有水印的高動(dòng)態(tài)范圍圖像受常規(guī)攻擊后得到的高動(dòng)態(tài)范圍圖像或受色調(diào)映射后得到的低動(dòng)態(tài)范圍圖像,記為istego,istego的高度為rorg,istego的寬度為corg。②_2、將istego從rgb空間轉(zhuǎn)化到hsv空間,得到istego的h層圖像、s層圖像和v層圖像,對(duì)應(yīng)記為istego-h、istego-s和istego-v,選取istego-v作為提取水印的圖像;然后判斷rorg和corg是否能夠被4整除,如果rorg和corg均能夠被4整除,則對(duì)istego-v不進(jìn)行像素點(diǎn)填充,并將istego-v重新記為istego-v-new,istego-v-new的高度為rnew,istego-v-new的寬度為cnew,rnew=rorg,cnew=corg;如果rorg能夠被4整除,而corg不能夠被4整除,則在istego-v的最右側(cè)填充4-mod(corg,4)列像素點(diǎn),且令填充的所有像素點(diǎn)的像素值為0,并將填充后得到的圖像記為istego-v-new,istego-v-new的高度為rnew,istego-v-new的寬度為cnew,rnew=rorg,cnew=4-mod(corg,4)+corg;如果rorg不能夠被4整除,而corg能夠被4整除,則在istego-v的最下側(cè)填充4-mod(rorg,4)行像素點(diǎn),且令填充的所有像素點(diǎn)的像素值為0,并將填充后得到的圖像記為istego-v-new,istego-v-new的高度為rnew,istego-v-new的寬度為cnew,rnew=4-mod(rorg,4)+rorg,cnew=corg;如果rorg和corg均不能夠被4整除,則在istego-v的最右側(cè)填充4-mod(corg,4)列像素點(diǎn),在istego-v的最下側(cè)填充4-mod(rorg,4)行像素點(diǎn),且令填充的所有像素點(diǎn)的像素值為0,并將填充后得到的圖像記為istego-v-new,istego-v-new的高度為rnew,istego-v-new的寬度為cnew,rnew=4-mod(rorg,4)+rorg,cnew=4-mod(corg,4)+corg。②_3、對(duì)istego-v-new進(jìn)行一級(jí)冗余小波變換(rdwt),得到istego-v-new的低頻近似分量、水平方向高頻分量、垂直方向高頻分量和對(duì)角線方向高頻分量,對(duì)應(yīng)記為llstego-v-new、hlstego-v-new、lhstego-v-new和hhstego-v-new;然后將llstego-v-new劃分成個(gè)互不重疊的尺寸大小為4×4的分塊,將llstego-v-new中坐標(biāo)位置為(x,y)的分塊記為bll,stego(x,y);其中,②_4、讀取水印秘鑰k;然后根據(jù)水印秘鑰k中的每個(gè)像素點(diǎn)的像素值,判定llstego-v-new中對(duì)應(yīng)坐標(biāo)位置的分塊是否提取水印,對(duì)于bll,stego(x,y),如果markorg-v-new(x,y)=0,則判定bll,stego(x,y)不需提取水?。蝗绻鹠arkorg-v-new(x,y)=1,則判定bll,stego(x,y)需提取水印。②_5、對(duì)llstego-v-new中判定為需提取水印的每個(gè)分塊進(jìn)行奇異值分解(svd)變換,得到llstego-v-new中判定為需提取水印的每個(gè)分塊的維數(shù)為4×4的u矩陣分量、維數(shù)為4×4的s矩陣分量、維數(shù)為4×4的v矩陣分量;然后在llstego-v-new中判定為需提取水印的每個(gè)分塊的u矩陣分量中的第2行第1列元素和第4行第1列元素中提取水印,假設(shè)bll,stego(x,y)判定為需提取水印,則在bll,stego(x,y)的u矩陣分量中的第2行第1列元素u*21和第4行第1列元素u*41中提取水印,如果u*21>u*41,則提取出的二值水印為0;如果u*21≤u*41,則提取出的二值水印為1;再根據(jù)提取出的所有二值水印,得到提取出的二值水印序列,記為w*。為進(jìn)一步說(shuō)明本發(fā)明方法的可行性和有效性,對(duì)本發(fā)明方法進(jìn)行測(cè)試。在本實(shí)施例中,選取公開(kāi)的高動(dòng)態(tài)范圍圖像庫(kù)和現(xiàn)有的色調(diào)映射算子進(jìn)行測(cè)試,其來(lái)源分別為:1)高動(dòng)態(tài)范圍圖像庫(kù)來(lái)源:①gregward的網(wǎng)站:http://www.anyhere.com/gward/hdrenc/pages/originals.html,共計(jì)33幅高動(dòng)態(tài)范圍圖像,動(dòng)態(tài)范圍的數(shù)量級(jí)為2.0~8.9,尺寸范圍為512×346~6144×6144。②imagegallery的網(wǎng)站:http://resources.mpi-inf.mpg.de/hdr/gallery.html,共計(jì)7幅高動(dòng)態(tài)范圍圖像,動(dòng)態(tài)范圍的數(shù)量級(jí)為3.0~8.6,尺寸范圍為760×1016~3270×1396。2)色調(diào)映射算子來(lái)源:hdrtoolboxformatlab,版本為1.1.0,共計(jì)27個(gè)色調(diào)映射算子。在本實(shí)施例中選用測(cè)試數(shù)據(jù)庫(kù)中的15幅高動(dòng)態(tài)范圍圖像,包含gregward的網(wǎng)站中的11幅高動(dòng)態(tài)范圍圖像,imagegallery的網(wǎng)站中的4幅高動(dòng)態(tài)范圍圖像,具體信息如表1所列;并選用hdrtoolbox中的26個(gè)色調(diào)映射算子(其中banterle色調(diào)映射算子要求高動(dòng)態(tài)范圍圖像的數(shù)值范圍為0.015~3000cd/m2,不具有通用性,因此本發(fā)明不再討論),具體信息如表2所列。表1高動(dòng)態(tài)范圍圖像的具體信息表2色調(diào)映射算子的具體信息序號(hào)名稱日期序號(hào)名稱日期1ashikhmin200214logarithmic2bruceexpoblend201315mertens20073chiu199316normalize4drago200317pattanaik20005durand200218ramantmo20096exponential19reinharddevlin20057fattal200220reinhard20028ferwerda199621schlick19959gamma22tumblinrushmeier199310kimkautzconsistent200823vanhateren200611krawczyk200524wardglobal199412kuang200725wardhistadj199713lischinski200626yee2003對(duì)上述15幅高動(dòng)態(tài)范圍圖像,按本發(fā)明方法的水印嵌入過(guò)程嵌入二值水印信息,并按本發(fā)明方法的水印提取過(guò)程提取出二值水印信息;然后采用信噪比(snr)、結(jié)構(gòu)相似度(ssim)和hdr-vdp-2方法來(lái)評(píng)價(jià)嵌有水印的高動(dòng)態(tài)范圍圖像相對(duì)于原始高動(dòng)態(tài)范圍圖像的圖像質(zhì)量值,即對(duì)水印不可見(jiàn)性的評(píng)價(jià);同時(shí),采用誤碼率(ber)來(lái)衡量提取的二值水印信息和原始二值水印信息的相似度,即對(duì)水印魯棒性的評(píng)價(jià)。在嵌有水印的高動(dòng)態(tài)范圍圖像相對(duì)于原始高動(dòng)態(tài)范圍圖像的圖像質(zhì)量值達(dá)到人工設(shè)定的圖像質(zhì)量值時(shí),水印信息的嵌入容量(bit)越大,誤碼率(ber)越小,則數(shù)字水印性能越好;反之,則越差。表3列出了人工設(shè)定的圖像質(zhì)量值為90時(shí),獲取嵌有水印的高動(dòng)態(tài)范圍圖像的相關(guān)測(cè)試數(shù)據(jù);表4列出了嵌有水印的高動(dòng)態(tài)范圍圖像受色調(diào)映射攻擊后得到的低動(dòng)態(tài)范圍圖像中提取的二值水印信息的誤碼率。從表3和表4所列的數(shù)據(jù),可知人工設(shè)定的圖像質(zhì)量值為90時(shí),本發(fā)明方法的平均嵌入容量為65035bit,嵌有水印的高動(dòng)態(tài)范圍圖像受色調(diào)映射攻擊后平均誤碼率為12.12%,表現(xiàn)出了優(yōu)異的嵌入容量和色調(diào)映射算子魯棒性。表3人工設(shè)定的圖像質(zhì)量值為90時(shí),獲取嵌有水印的高動(dòng)態(tài)范圍圖像的相關(guān)測(cè)試數(shù)據(jù)表4嵌有水印的高動(dòng)態(tài)范圍圖像受色調(diào)映射攻擊后得到的低動(dòng)態(tài)范圍圖像中提取的二值水印信息的平均誤碼率f.guerrini,m.okuda,n.adami,r.leonardi,"highdynamicrangeimagewatermarkingrobustagainsttone-mappingoperators".ieeetrans.inf.forensicssecur.6(2),pp.283-295,2011.(guerrini、okuda等人,對(duì)抗色調(diào)映射的魯棒性高動(dòng)態(tài)范圍圖像水印,ieee信息取證和安全匯刊6(2),283-295,2011)、e.maiorana,v.solachidis,p.campisi,y.lou,"robustmulti-bitwatermarkingforhdrimagesintheradon-dctdomain".in:proc.int.symposiumonimageandsignalprocessingandanalysis,trieste,italy,pp.284-289,2013.(maiorana、solachidis等人,radon-dct域的魯棒性多位高動(dòng)態(tài)范圍圖像水印,圖像信號(hào)處理和分析專題國(guó)際會(huì)議,意大利的里雅斯特,284-289,2013)、j.l.wu,"robustwatermarkingframeworkforhighdynamicrangeimagesagainsttone-mappingattacks".watermarking2,229-242,2012.(吳,對(duì)抗色調(diào)映射攻擊的魯棒性高動(dòng)態(tài)范圍圖像水印結(jié)構(gòu),水印2,229-242,2012)是現(xiàn)有的對(duì)抗色調(diào)映射算子的頂尖算法,在此對(duì)現(xiàn)有的三種對(duì)抗色調(diào)映射算子的頂尖算法和本發(fā)明方法分別在色調(diào)映射算子魯棒性、平均誤碼率、嵌入容量以及不可感知性進(jìn)行對(duì)比。為了公平起見(jiàn),本發(fā)明方法分別選取上述第一種和第二種對(duì)抗色調(diào)映射算子的頂尖算法中的測(cè)試圖像和色調(diào)映射算子(tmo)進(jìn)行實(shí)驗(yàn),對(duì)比數(shù)據(jù)如表5、表6、表7所列。鑒于第三種對(duì)抗色調(diào)映射算子的頂尖算法需要預(yù)先確定色調(diào)映射算子,缺乏實(shí)用性,故不進(jìn)行重點(diǎn)對(duì)比。從表5和表6所列的數(shù)據(jù)可以看出,本發(fā)明方法在設(shè)定q=90的情況下,嵌入容量在第一種和第二種對(duì)抗色調(diào)映射算子的頂尖算法的20倍以上;同時(shí),在受色調(diào)映射(tm)攻擊前后進(jìn)行水印提取,本發(fā)明方法的誤碼率也明顯優(yōu)于第一種和第二種對(duì)抗色調(diào)映射算子的頂尖算法。通過(guò)表7的綜合對(duì)比可知,本發(fā)明方法在對(duì)抗色調(diào)映射算子的魯棒性上優(yōu)于上述三種對(duì)抗色調(diào)映射算子的頂尖算法,表現(xiàn)出了優(yōu)異的性能和普遍的實(shí)用性。表5本發(fā)明方法和第一種對(duì)抗色調(diào)映射算子的頂尖算法在嵌入容量(bit)和受色調(diào)映射攻擊前后水印提取的誤碼率(%)對(duì)比表6本發(fā)明方法和第二種對(duì)抗色調(diào)映射算子的頂尖算法在嵌入容量(bit)和受色調(diào)映射攻擊前后水印提取的誤碼率(%)對(duì)比表7本發(fā)明方法與上述三種對(duì)抗色調(diào)映射算子的頂尖算法在綜合性能上的對(duì)比當(dāng)前第1頁(yè)12