專利名稱：數(shù)字水印的嵌入方法、提取方法及其裝置的制作方法
技術(shù)領(lǐng)域：
本發(fā)明涉及數(shù)字水印信息技術(shù)領(lǐng)域，尤其涉及數(shù)字水印的嵌入方法、提取方法及
其裝置。
背景技術(shù)：
在全球信息數(shù)字化進程的推動下，二值圖像作為信息流通的一種重要方式，廣泛 的分布于各種應用環(huán)節(jié)。我們平時接觸最多的圖像、文字等一般都需要通過印刷到紙質(zhì)媒 體上輸出，如個人檔案、學歷證書、專利證件、商標、產(chǎn)品包裝、護照等等。而隨著信息技術(shù) 的發(fā)展和高質(zhì)量圖像輸入輸出設(shè)備的進步，數(shù)碼相機、錄像機、掃描儀和打印機等設(shè)備和功 能強大的軟件，特別是高精度彩色噴墨、激光打印機和高精度掃描儀的普及，為人們的工作 和生活帶來了許多便利，但同時也使各種二值圖像的偽造及篡改變得非常容易。出于版權(quán) 保護、數(shù)據(jù)完整和數(shù)據(jù)隱寫的考慮，對二值圖像數(shù)字水印技術(shù)的研究提出了更為嚴峻的要 求。 由于二值圖像本身的特殊性，g卩能夠利用的圖像信息僅由黑白兩種像素點構(gòu)成， 隨意改變其中任何一種像素點都會造成比較明顯的視覺影響，因此在二值圖像中實施數(shù)字 水印技術(shù)存在一定的難度。目前，針對二值圖像提出的數(shù)字水印技術(shù)中，主要有如下幾種
現(xiàn)有技術(shù)一提出了一種基于像素游程的水印嵌入方法，主要用于傳真文件的信息 隱藏，該方法的具體技術(shù)原理是通過修改游程邊緣的像素點，即修改游程的長度，將游程 加1或減l，最終用游程的奇偶性來攜帶信息。該方法在游程較多時會出現(xiàn)較大的失真，很 容易引起視覺上的感知，產(chǎn)生明顯的凹凸不平。并且。該方法嵌入的水印信息對打印、掃描 等抵抗能力很差，很容易由于打印、掃描之類的攻擊而改變游程的奇偶性，對水印信息的提 取造成一定的困難，無法準確地檢測出攜帶的水印信息。 現(xiàn)有技術(shù)二提出了一種調(diào)制行間距的數(shù)字水印嵌入方法，類似地也可以通過調(diào)制 字間距實現(xiàn)數(shù)字水印的嵌入，但是這種方法在提取數(shù)字水印信息時需要原始文檔的參與， 即通過與原始文檔的行間距或字間距進行對比，確定文檔中行間距或字間距的變化，進而 根據(jù)確定結(jié)果提取水印信息，這種方法最大的缺點是在提取水印信息時需要原始文檔參 與，無法實現(xiàn)盲檢測。并且，在水印信息嵌入過程中，若對行間距或字間距改變較多，則容易 引起視覺上的感知，若對行間距或字間距改變較少，則不易檢測出來。 現(xiàn)有技術(shù)三提出了一種基于文字特征的信息隱藏方法，該方法通過改變文字的局
部特征以嵌入水印信息，例如，對文字中個別筆畫進行上移、下移或伸縮等形成不同的變形
字，并通過這些變形字攜帶水印信息。由于需要對文字之間的筆畫連接關(guān)系做變形，因此容
易引起視覺上的感知，并且該方法嵌入的水印信息對打印、掃描等抵抗能力較差，很容易由
于打印、掃描之類的攻擊而改變筆畫的連接關(guān)系，因此無法準確地檢測出變形字。 綜上所述，現(xiàn)有二值圖像的數(shù)字水印嵌入技術(shù)，容易引起視覺上的感知，對視覺效
果影響明顯，使通過對原始文檔的修改或變形而攜帶的數(shù)字水印容易被發(fā)現(xiàn)，導致數(shù)字水
印的隱藏性差；并且，基于上述現(xiàn)有技術(shù)嵌入的數(shù)字水印抵抗打印、掃描等攻擊的能力差，
5使得嵌入的數(shù)字水印的魯棒性差，進而給提取數(shù)字水印信息帶來困難。

發(fā)明內(nèi)容
本發(fā)明提供一種數(shù)字水印的嵌入方法及其裝置，以提高圖像中數(shù)字水印信息的隱 藏性以及魯棒性。 本發(fā)明提供一種數(shù)字水印的提取方法及其裝置，以提高圖像中提取數(shù)字水印信息 的準確度。 本發(fā)明實施例通過如下技術(shù)方案實現(xiàn) 本發(fā)明實施例提供了一種數(shù)字水印的嵌入方法，包括如下步驟 根據(jù)分塊規(guī)則將圖像劃分為設(shè)定數(shù)目的圖像塊，統(tǒng)計每個圖像塊的特征值，并根
據(jù)排列規(guī)則將統(tǒng)計出的特征值組成第一矩陣； 將所述第一矩陣進行離散余弦變換得到第二矩陣，并根據(jù)預先定義的與數(shù)字水印 對應的閾值替換所述第二矩陣中的指定特征值得到第三矩陣；
將所述第三矩陣進行離散余弦逆變換得到第四矩陣； 根據(jù)所述第四矩陣中的特征值相對于所述第一矩陣中的相應特征值的差值，確定
待修改的圖像塊以及該圖像塊的修改方式，并根據(jù)確定結(jié)果修改所述圖像。 本發(fā)明實施例還提供了一種數(shù)據(jù)水印的提取方法，包括如下步驟 根據(jù)與對圖像嵌入數(shù)字水印時相同的分塊規(guī)則將所述圖像劃分為設(shè)定數(shù)目的圖
像塊，統(tǒng)計每個圖像塊的特征值，并根據(jù)與對所述圖像嵌入數(shù)字水印時相同的排列規(guī)則將
統(tǒng)計出的特征值組成第一矩陣； 將所述第一矩陣進行離散余弦變換得到第二矩陣，并根據(jù)所述第二矩陣中的特征 值確定對所述圖像嵌入數(shù)字水印時用于替換指定特征值的閾值； 根據(jù)確定的所述閾值以及設(shè)定的待嵌入數(shù)字水印與閾值的對應關(guān)系，提取嵌入的 數(shù)字水印。
本發(fā)明實施例還提供了一種數(shù)字水印的嵌入裝置，包括 圖像劃分單元，用于根據(jù)分塊規(guī)則將圖像劃分為設(shè)定數(shù)目的圖像塊； 特征值統(tǒng)計單元，用于統(tǒng)計圖像劃分單元劃分得到的每個圖像塊的特征值； 矩陣處理單元，用于根據(jù)排列規(guī)則將特征值統(tǒng)計單元統(tǒng)計出的特征值組成第一矩
陣，將所述第一矩陣進行離散余弦變換得到第二矩陣，根據(jù)預先定義的與數(shù)字水印對應的
閾值替換所述第二矩陣中的指定特征值得到第三矩陣，并將所述第三矩陣進行離散余弦逆
變換得到第四矩陣； 水印嵌入單元，用于根據(jù)矩陣處理單元得到的第四矩陣中的特征值相對于第一矩 陣中的相應特征值的差值，確定待修改的圖像塊以及該圖像塊的修改方式，并根據(jù)確定結(jié) 果修改所述圖像。 本發(fā)明實施例還提供了一種數(shù)字水印的提取裝置，包括 圖像劃分單元，用于根據(jù)與對圖像嵌入數(shù)字水印時相同的分塊規(guī)則將所述圖像劃 分為設(shè)定數(shù)目的圖像塊； 特征值統(tǒng)計單元，用于統(tǒng)計圖像劃分單元劃分得到的每個圖像塊的特征值；
矩陣處理單元，用于根據(jù)與對所述圖像嵌入數(shù)字水印時相同的排列規(guī)則將特征值統(tǒng)計單元統(tǒng)計出的特征值組成第一矩陣，將所述第一矩陣進行離散余弦變換得到第二矩 陣； 水印提取單元，用于根據(jù)矩陣處理單元得到的第二矩陣中的特征值確定對所述圖
像嵌入數(shù)字水印時用于替換指定特征值的閾值，并根據(jù)確定的所述閾值以及設(shè)定的待嵌入
數(shù)字水印與閾值的對應關(guān)系，提取嵌入的數(shù)字水印。
通過上述技術(shù)方案，本發(fā)明實施例能夠產(chǎn)生如下有益效果 在對圖像嵌入數(shù)字水印時，首先在空域上對該圖像劃分并組織成空域特征值矩陣 (即第一矩陣)，進一步進行離散余弦變換得到頻域特征值矩陣(即第二矩陣)，修改該頻域 特征值矩陣中的指定特征值后，進行離散余弦逆變換得到空域特征值矩陣(即第四矩陣)， 并根據(jù)經(jīng)過上述變化的空域特征值矩陣中的特征值相對于變化前的空域特征值矩陣中的 相應特征值的差值，確定待修改的圖像塊以及該圖像塊的修改方式，并基于該確定結(jié)果修 改圖像。該方案充分利用圖像在空域和頻域上的轉(zhuǎn)換關(guān)系，將圖像的空域特征和頻域特征 相結(jié)合，使得嵌入的數(shù)字水印的視覺隱藏效果好，并且抵抗打印、掃描等攻擊的魯棒性高。
在對圖像提取數(shù)字水印時，根據(jù)與對圖像嵌入數(shù)字水印時相同的分塊規(guī)則以及排 列規(guī)則，對圖像劃分并組織成空域特征值矩陣(即第一矩陣)，進一步進行離散余弦變換得 到頻域特征值矩陣(即第二矩陣)，并根據(jù)該頻域特征值矩陣中的特征值確定對圖像嵌入 數(shù)字水印時用于替換指定特征值的閾值，從而根據(jù)確定出的閾值以及設(shè)定的待嵌入數(shù)字水 印與閾值的對應關(guān)系，提取嵌入的數(shù)字水印。該方案充分利用頻域特征值矩陣相對穩(wěn)定的 特征，根據(jù)頻域特征值矩陣中的特征值提取嵌入的數(shù)字水印，準確度高。


圖1為本發(fā)明實施例中對圖像嵌入數(shù)字水印的方法流程圖； 圖2為本發(fā)明實施例中對圖像進行分塊的示意圖一 ； 圖3為本發(fā)明實施例中各個矩陣的關(guān)系示意圖； 圖4為本發(fā)明實施例中對圖像提取數(shù)字水印的方法流程圖； 圖5為本發(fā)明實施例中對圖像進行分塊的示意圖二 ； 圖6為本發(fā)明實施例中對圖5所示圖像嵌入數(shù)字水印后的圖像示意圖； 圖7為本發(fā)明實施例中數(shù)字水印的嵌入裝置的結(jié)構(gòu)示意圖； 圖8為本發(fā)明實施例中數(shù)字水印的提取裝置的結(jié)構(gòu)示意圖。
具體實施例方式
本發(fā)明基于大量圖像在打印過程中呈現(xiàn)的基本規(guī)律，結(jié)合頻域DCT (Discrete Cosine Transformation,離散余弦變換)的系數(shù)特征，提出數(shù)字水印的嵌入方法、提取方法 及其裝置，以提高圖像中嵌入的數(shù)字水印的視覺隱藏效果和魯棒性，以及提高數(shù)字水印提 取的準確度。 通過大量的試驗，本發(fā)明實施例總結(jié)出圖像在打印過程中呈現(xiàn)的基本規(guī)律，以及
頻域DCT的系數(shù)特征，具體如下 1、圖像在打印過程中呈現(xiàn)的基本規(guī)律 在相同環(huán)境下采用高質(zhì)量的打印機輸出，分別對中文楷體、中文宋體、中文黑體以
7及英文、日文等做了大量的測試，發(fā)現(xiàn)打印前后，圖像中每個區(qū)域的黑點數(shù)變化量是隨機 的，是不可預見的，但是在每個區(qū)域范圍內(nèi)，黑點數(shù)的變化方向是一定的，要么增加要么減 少，也就是說，打印對整個圖像中的黑像素點數(shù)目的影響趨勢是一致的。
2、頻域DCT的系數(shù)特征 針對 一 個實數(shù)矩陣，變化該矩陣中的 一 個固定值后進行D C T ，經(jīng)過多次試驗表 明在DCT系數(shù)上反映的僅僅是直流分量的變換，而交流系數(shù)基本不變。進一步地，對該 實數(shù)矩陣進行DCT后，對DCT系數(shù)的低頻、中頻和高頻部分分別做相應的修改，然后進行 IDCT(Inverse Discrete Consine Transform,離散余弦逆變換)，經(jīng)過多次試驗表明修改 高頻系數(shù)對原實數(shù)矩陣(即進行DCT前的實數(shù)矩陣)產(chǎn)生的影響要遠遠小于修改低頻系數(shù) 對原實數(shù)矩陣產(chǎn)生的影響，影響程度分別按照低、中、高頻的順序依次降低。并且，影響會因 為被修改系數(shù)的個數(shù)增加而增加，被修改的系數(shù)越多，原實數(shù)矩陣變化越大。最后，固定待 修改系數(shù)位置，并用不同的閾值替換待修改系數(shù)，發(fā)現(xiàn)該閾值的取值直接關(guān)系著原實數(shù)矩 陣的變化，閾值的取值越大，對原實數(shù)矩陣產(chǎn)生的影響越明顯，反之越小。
基于上述試驗結(jié)論，本發(fā)明實施例提出了在圖像中嵌入數(shù)字水印和提取數(shù)字水印 的方法，其主要思想包括在對圖像嵌入數(shù)字水印的過程中，首先在空域上對圖像劃分并組 織成空域特征值矩陣，通過DCT將該空域特征值矩陣轉(zhuǎn)換為頻域特征值矩陣，修改該頻域 特征值矩陣中的指定特征值后，進行IDCT返回空域特征值矩陣，并根據(jù)經(jīng)過上述變化的空 域特征值矩陣中的特征值相對于變化前的空域特征值矩陣中的相應特征值的差值，確定待 修改的圖像塊以及該圖像塊的修改方式，并基于該確定結(jié)果修改圖像，從而完成了數(shù)字水 印的嵌入。在提取過程中，統(tǒng)計相對應嵌入數(shù)字水印時劃分的各個圖像塊的特征值，組成特 征值矩陣并進行DCT得到頻域特征值矩陣，根據(jù)該頻域特征值矩陣中的特征值，即可提取 出數(shù)字水印。 下面結(jié)合說明書附圖對本發(fā)明實施例的主要實現(xiàn)原理、具體實施過程及其對應能 夠達到的有益效果進行詳細的闡述 如圖1所示，本發(fā)明一個實施例中，對圖像嵌入數(shù)字水印時，執(zhí)行如下步驟
步驟101、根據(jù)設(shè)定分塊規(guī)則，將待嵌入數(shù)字水印的圖像劃分為設(shè)定數(shù)目的圖像 塊。 該步驟中，分塊規(guī)則可以根據(jù)實際需要靈活設(shè)定，例如分塊規(guī)則中可以包括要將
圖像劃分的圖像塊數(shù)目以及每個圖像塊的大小。對于要將圖像劃分的圖像塊數(shù)目，一般以
能組成一個矩陣并且該矩陣為方陣為最佳，即將圖像劃分為NXS塊，其中N等于S為最佳，
在N不等于S的情況下，N和S越接近越好。對于每個圖像塊的大小，一般以在提取數(shù)字水
印時劃分的圖像塊的邊界區(qū)域能夠和嵌入數(shù)字水印時劃分的圖像塊的邊界區(qū)域吻合為最
佳，如圖2所示，為一種以文字邊界劃分圖像塊的方法示例。 步驟102、統(tǒng)計待嵌入數(shù)字水印的圖像中每個圖像塊的特征值Ei.。 該步驟中，統(tǒng)計圖像塊的特征值即統(tǒng)計該圖像塊中的黑像素點數(shù)在該待嵌入數(shù)字
水印的圖像包含的總像素點數(shù)中所占的比例，具體過程為 統(tǒng)計每個圖像塊包含的黑像素點數(shù)和待嵌入數(shù)字水印的圖像中包含的總像素點 數(shù)； 將統(tǒng)計得到的各個圖像塊中的黑像素點數(shù)與總像素點數(shù)的比值確定為各個圖像
8塊的特征值。 步驟103、根據(jù)設(shè)定排列規(guī)則將統(tǒng)計出的特征值組成特征值矩陣M。 該步驟中，排列規(guī)則表示圖像塊在待嵌入數(shù)字水印的圖像中的位置與該圖像塊的
特征值在特征值矩陣M中的位置之間的對應關(guān)系，具體可以根據(jù)需要靈活設(shè)定，例如，可以
根據(jù)圖像塊在圖像中從左到右以及從上到下的順序，將各個圖像塊對應的特征值從矩陣的
第一行第一列開始按照從左到右以及從上到下的順序依次排列，形成特征值矩陣。 步驟104、對特征值矩陣M進行DCT得到頻域特征值矩陣C。 步驟105、根據(jù)預先定義的與數(shù)字水印對應的閾值替換頻域特征值矩陣C中的指 定特征值得到頻域特征值矩陣C'。 該步驟中，預先定義與數(shù)字水印對應的閾值的過程為 首先，確定一個理想的閾值K，該閾值K是根據(jù)經(jīng)驗而設(shè)定的，具體地，要保證視覺 效果和魯棒性之間的平衡，即為了保證視覺效果，根據(jù)閾值K影響的待翻轉(zhuǎn)像素點要盡量 少，但也不能太少，否則會給數(shù)字水印的提取帶來困難。 其次，設(shè)定待嵌入數(shù)字水印與閾值K的對應關(guān)系，例如，需要嵌入的二進制序列由
{0， 1}組成，則用K代表0，用-K代表1 ;或者，用K代表1，用-K代表0。 步驟106、將頻域特征值矩陣C'進行IDCT得到新的特征值矩陣M'。 為了使上述矩陣間的轉(zhuǎn)換關(guān)系更加明確，如圖3所述，列出了特征值矩陣M、M'和
經(jīng)過DCT的頻域特征值矩陣C、C'的關(guān)系示意圖。 步驟107、根據(jù)特征值矩陣M'中的特征值相對于特征值矩陣M中的相應特征值的 差值X(j,k)，確定待修改的圖像塊以及該圖像塊的修改方式。
該步驟中，差值X(j,k)的計算公式如下
X(j,k)=M' (j,k)-M(j,k) 其中，j表示該特征值在特征值矩陣M'以及在特征值矩陣M中所在的行數(shù)，k表 示該特征值在特征值矩陣M'以及在特征值矩陣M中所在的列數(shù)。 確定待修改的圖像塊以及該圖像塊的修改方式的具體過程視X(j, k)的正負而有所 不同，具體地 當差值X(j,k)為正數(shù)時，確定存在該差值的特征值對應的圖像塊為待修改的圖像 塊，并根據(jù)該差值確定該圖像塊的修改方式為將與該差值相應數(shù)目的白像素點翻轉(zhuǎn)為黑 像素點，即將X(j,k)個白像素點翻轉(zhuǎn)為黑像素點； 當差值X(j,k)為負數(shù)時，確定存在該差值的特征值對應的圖像塊為待修改的圖像 塊，并根據(jù)該差值確定該圖像塊的修改方式為將與該差值相應數(shù)目的黑像素點翻轉(zhuǎn)為白 像素點，即將lx(j,k)l個黑像素點翻轉(zhuǎn)為白像素點。 步驟108、根據(jù)確定出的待修改圖像塊以及每個待修改圖像塊的修改方式，修改該 待嵌入數(shù)字水印的圖像中的像素點。 該步驟中，根據(jù)待修改的圖像塊以及該圖像塊的修改方式，修改該圖像塊中的像 素點時，像素點的翻轉(zhuǎn)優(yōu)先級可以參考Min Wu[2004]提出的翻轉(zhuǎn)圖像中黑白像素點的方 法，該方法首先通過對局部光滑行和局部連通性的考察，確定出圖像中的可修改點，并為每 一個可修改點給出了一個評分標準，修改分數(shù)越高的點對圖像的視覺影響就越??；通過置 亂、分塊等方法均勻分布每個塊中的可修改點，采用奇偶量化的方法來嵌入數(shù)字水印信息。
進一步地，上述流程中的步驟105中，頻域特征值矩陣C中的指定特征值是根據(jù)實 際效果預先確定的，具體確定過程包括 首先，確定特征值矩陣C中要被替換的特征值的數(shù)目T，具體如下
特征值矩陣C中要修改的特征值的個數(shù)T，決定了信息量嵌入的多少，同時也影響 著視覺效果。修改的個數(shù)越多，嵌入的信息量(即數(shù)據(jù)的bit容量)越大，但需翻轉(zhuǎn)的像素 點越多，對圖像的破壞也越大。T的取值主要依據(jù)源圖像和修改后圖像的相關(guān)度(即圖像的 相似度)以及變化的總像素點數(shù)目決定，假設(shè)修改T個，即意味著嵌入了 T個bit數(shù)據(jù)。
其中，源圖像和修改后圖像的相關(guān)度可以通過如下公式得到
<formula>formula see original document page 10</formula> 源圖像和修改后圖像中變化的總像素點數(shù)目，可以通過如下公式得到
<formula>formula see original document page 10</formula> x = 根據(jù)得到的變化前后兩個矩陣的相關(guān)度以及變化的總點數(shù)，選取合適的T值，具 體地，兩個矩陣的相關(guān)度越大，變化的總點數(shù)越少，則K的取值越佳。 確定出特征值矩陣C中要被替換的特征值的數(shù)目T后，在具體替換時，需要對要替 換的特征值定位，根據(jù)上述總結(jié)出的頻域DCT系數(shù)的特征，即修改高頻系數(shù)對原實數(shù)矩陣 產(chǎn)生的影響相對最小，本實施例采用逆Zig-Zag的順序，從矩陣的最右下角的特征值開始 依次修改。 本發(fā)明又一實施例中，對應于上述流程，在步驟102中，統(tǒng)計圖像塊的特征值的方 法還可以包括 統(tǒng)計每個圖像塊包含的白像素點數(shù)和待嵌入數(shù)字水印的圖像包含的總像素點數(shù)， 并將統(tǒng)計得到的各個圖像塊中的白像素點數(shù)與總像素點數(shù)的比值確定為各個圖像塊的特 征值。 相應地，對應于上述流程，在步驟107中，根據(jù)特征值矩陣M'中的特征值相對于 特征值矩陣M中的相應特征值的差值，確定待修改的圖像塊以及該圖像塊的修改方式的過 程包括 當差值X(j,k)為正數(shù)時，確定存在該差值的特征值對應的圖像塊為待修改的圖像 塊，并根據(jù)該差值確定該圖像塊的修改方式為將與該差值相應數(shù)目的黑像素點翻轉(zhuǎn)為白 像素點，即將X(j,k)個黑像素點翻轉(zhuǎn)為白像素點； 當差值X(j,k)為負數(shù)時，確定存在該差值的特征值對應的圖像塊為待修改的圖像 塊，并根據(jù)該差值確定該圖像塊的修改方式為將與該差值相應數(shù)目的白像素點翻轉(zhuǎn)為黑 像素點，即將lx(j,k)l個白像素點翻轉(zhuǎn)為黑像素點。 如圖4所示，本發(fā)明一個實施例中，對圖像提取數(shù)字水印時，執(zhí)行如下步驟
步驟401、根據(jù)與對圖像嵌入數(shù)字水印時相同的分塊規(guī)則，將該圖像劃分為設(shè)定數(shù) 目的圖像塊。
該步驟中，圖像可以為掃描得到的圖像，所采用的分塊規(guī)則與對該圖像嵌入數(shù)字
水印時采用的分塊規(guī)則相同，即要將圖像劃分的圖像塊數(shù)目與對該圖像嵌入數(shù)字水印時劃
分的圖像塊數(shù)目相同，并且每個圖像塊的大小與對該圖像嵌入數(shù)字水印時劃分的相應圖像
塊的大小統(tǒng)一，這樣可以提高提取數(shù)字水印的精度。 步驟402、統(tǒng)計劃分后得到的每個圖像塊的特征值Ei'。 該步驟統(tǒng)計特征值的方法與對圖像嵌入數(shù)字水印時統(tǒng)計特征值的方法一致，即與 上述步驟102中統(tǒng)計每個圖像塊所采用的方法一致，此處不再詳細描述。
步驟403、根據(jù)與對圖像嵌入數(shù)字水印時相同的排列規(guī)則，將統(tǒng)計出的特征值組成 特征值矩陣M〃 。 該步驟與上述步驟103所采用的排列規(guī)則一致，此處不再詳細描述。 步驟404、將征值矩陣M〃進行DCT得到頻域特征值矩陣C〃 。 步驟405、根據(jù)得到的頻域特征值矩陣C〃中的特征值，確定對該圖像嵌入數(shù)字水
印時用于替換指定特征值的閾值。 該步驟中，得到頻域特征值矩陣C〃后，根據(jù)對圖像嵌入數(shù)字水印時修改的特征
值的數(shù)目和修改順序，依次確定所選用的閾值。例如，嵌入數(shù)字水印時，采用逆Zig-Zag的
順序，從矩陣最右下角的特征值開始依次修改了 T個特征值，則提取數(shù)字水印時，也采用逆
Zig-Zag的順序，從頻域特征值矩陣C〃最右下角的特征值開始確定閾值。 步驟406、根據(jù)確定出的閾值以及設(shè)定的待嵌入數(shù)字水印與閾值的對應關(guān)系，確定
嵌入的數(shù)字水印。 該步驟中，若設(shè)定的待嵌入數(shù)字水印與閾值的對應關(guān)系為用K代表O，用-K代表 l，并且確定出的閾值依次為K、K、-K、K、-K、K，則據(jù)此確定出的數(shù)字水印即為{0，0，1，0，1，0}。 通過上述過程可以準確地提取出數(shù)字水印，進一步地，本發(fā)明又一實施例中，在得 到頻域特征值矩陣C〃后，以及在根據(jù)該矩陣中的特征值提取數(shù)字水印之前，可以對該矩陣 的可靠性進行驗證，具體地 對照得到的頻域特征值矩陣C〃與對該圖像嵌入數(shù)字水印時得到頻域特征值矩陣 C'(即上述步驟105中得到的得到頻域特征值矩陣C')，按照約定好的被替換特征值的 位置，檢查C〃中被替換的特征值與C'相應特征值的符號是否一致。通常具體的值會有不 同程度的改變，但是正負號會保持相對的穩(wěn)定。即當C'中某特征值被修改為K，那么C〃中 相應位置的符號應該保持不變，仍為"+ ";相反的，若C'中某系數(shù)被修改為-K，那么對應的
C〃中相應位置的符號應該保持"-"，允許值在一定程度上變化，但不會發(fā)生正負號翻轉(zhuǎn)。 根據(jù)上述過程，若被修改的所有特征值或大于設(shè)定閾值個特征值在C〃和C'中 的符號一致，則確定得到頻域特征值矩陣C〃可靠，可以根據(jù)該矩陣中的特征值提取數(shù)字水 印；若被修改的所有特征值或大于設(shè)定閾值個特征值在C〃和C'中的符號不一致，則確定 得到頻域特征值矩陣C〃不可靠，需要重新提取。其中，閾值可以根據(jù)被修改的特征值的總 數(shù)以及要達到的精度靈活設(shè)定。 根據(jù)上述過程提供的方法， 一個具體實施例如下 A、參見圖5，要對該圖所示的二值文本圖像嵌入數(shù)字水印時，將該圖像劃分為64 塊，其中圖像塊的邊界如虛線框所示(圖中僅示出了兩個圖像塊的邊界)。
11
B、統(tǒng)計每個圖像塊的特征值Ei(i = 1、2......64)，具體公式為Ei = (NBi/NF) XT ; 其中，Ei為第i個圖像塊的特征值；NBi為第i個圖像塊包含的黑像素點數(shù)；NF為 圖5所示的圖像包含的所有像素點數(shù)；T是比例精度，根據(jù)希望達到的精度取值，例如可以 取100或100000 (本實施例中取T = 100000)。C、按照圖像塊從左到右從上到下的順序，將統(tǒng)計得到的64個特征值依次排列成矩陣M。該步驟得到的矩陣M如下1. 0e+00#0. 96531. 00331. 11650. 99751. 04401. 06370.93401. 08971. 04131. 06501. 07221. 12631. 15181. 12010.84231. 16031. 00380. 98451. 07311. 19831.06150. 86110.97291.01631. 10660. 77971. 24750. 94831. 07530.73760. 71300.99950. 86330. 98090. 91521. 01851. 11471. 12900.78141. 20520. 98771. 29671. 03421. 06061. 13711. 03061. 08871. 03081. 00291. 03281. 07131. 00471. 10221.01721.21350.99951. 16061. 02711. 24680. 96861. 00540. 91090.98360.9582D、對得到的矩陣M進行DCT，得到矩陣c。該步驟得到的矩陣C如下1. 0e+00#8.24400.1136"0.1333"0.19130.1574"0.16630.12600.1230"0. 0507"0. 0377"0. 0895"0. 07380.股5"0. 07080.1438"0.10870.15290. 02530. 06770. 0560"0.1868a 1009"0. 01930. 02870.1135"0. 21040. 0081a 1115"0.1460"0. 2091"0.1091"0.1731"0. 22130. 05680. 0333"0. 05820. 0802"0. 03880. 27810. 0993"0.10410.04600. 0165"0. 0432■0.04040.1902"0. 03320. 22550. 02780.10750. 0173"0. 0558"0. 07750. 0560"0.1504"0. 0132"0. 0182"0. 2193"0. 03660. 0351"0. 0107a 0131"0. 0657"0.1369E、修改矩陣C中的指定特征值，取閾值K=500，修改個數(shù)T=10。修改后的矩陣c'如下1. 0e+00#8.24400.1136"0.1333"0.19130.1574"0.16630.12600.1230"0. 0507"0. 0377"0. 0895"0. 07380.股5"0. 07080.1438"0.10870.15290. 02530. 06770. 0560"0.1868a誦"0. 01930. 02870.1135"0. 21040. 0081a 1115"0.1460"0. 2091"0.1091"0.1731"0. 22130. 05680. 0333"0. 05820. 0802"0. 03880. 2781"0. 0500"0.10410.04600. 0165"0. 0432■0.04040.1902a 05000. 05000. 02780.10750. 0173"0. 0558"0. 0775a 0500"0. 0500a 0500"0. 0182"0. 2193"0. 03660. 0351"0. 0500a 0500"0. 0500a 0500
12
F、對矩陣C'進行IDCT，得到矩陣M'。
修改后的矩陣M'如下
1.0e+004*0.96611.00301.10871.0237a99831.1147a89651.10301.03121.09111.04931.11631.20461.0506a89311.14311.0270a91641.16751.11621.1087a8345a99631.00481.0878a84131.14721.0591a9752a8314a63061.0349a8627a9719a9441a9796i.1576i.0712a84671.1745a9966i.27921.05021.0423i.1548i.0323i.06541.0456i.0013i.03381.07151.0116i.0879i.0236i.2214a9929i.1586i.03341.2383a9740i.0047a9118a9796a9608 G、根據(jù)矩陣M'中特征值與矩陣M中相應特征值的差值，確定每個圖像塊中需要翻
轉(zhuǎn)的像素點。 確定結(jié)果如下2_1-1759-102114-8429-2358-51-22118-155114-3852-152211-184106-5952-25-42138-224248-224210-18478_1-2065-8796-129146-6820-3936-41404-5233_32016-321418-15_414-1912_12_96 H、根據(jù)上述確定結(jié)果，翻轉(zhuǎn)每個圖像塊中的像素點，其中 若差值為正數(shù)，則將該圖像塊中相應數(shù)目的白像素點翻轉(zhuǎn)為黑像素點，若差值為
負數(shù)，則將該圖像塊中相應數(shù)目的黑像素點翻轉(zhuǎn)為白像素點。 根據(jù)上述過程嵌入數(shù)字水印后得到的二值文本圖像如圖6所示。 基于上述流程，本發(fā)明實施例還提供了一種數(shù)字水印的嵌入裝置，其結(jié)構(gòu)如圖7
所示，包括 圖像劃分單元701、特征值統(tǒng)計單元702、矩陣處理單元703以及水印嵌入單元 704。具體地: 圖像劃分單元701，用于根據(jù)設(shè)定分塊規(guī)則將圖像劃分為設(shè)定數(shù)目的圖像塊.
特征值統(tǒng)計單元702，用于統(tǒng)計圖像劃分單元701劃分得到的每個圖像塊的特征值。 矩陣處理單元703，用于根據(jù)設(shè)定排列規(guī)則將特征值統(tǒng)計單元702統(tǒng)計出的特征 值組成第一矩陣，將該第一矩陣進行離散余弦變換得到第二矩陣，根據(jù)預先定義的與數(shù)字 水印對應的閾值替換得到的第二矩陣中的指定特征值得到第三矩陣，并將得到的第三矩陣 進行離散余弦逆變換得到第四矩陣。 水印嵌入單元704，用于根據(jù)矩陣處理單元703得到的第四矩陣中的特征值相對 于第一矩陣中的相應特征值的差值，確定待修改的圖像塊以及該圖像塊的修改方式，并根據(jù)確定結(jié)果修改圖像。
進一步地，一個實施例中，上述特征值統(tǒng)計單元702可以包括 第一統(tǒng)計模塊，用于統(tǒng)計圖像塊包含的黑像素點數(shù)和圖像包含的總像素點數(shù)，并
將統(tǒng)計出的圖像塊的黑像素點數(shù)與總像素點數(shù)的比值確定為該圖像塊的特征值； 或者，第二統(tǒng)計模塊，用于統(tǒng)計圖像塊包含的白像素點數(shù)和圖像包含的總像素點
數(shù)，并將統(tǒng)計得到的白像素點數(shù)與總像素點數(shù)的比值確定為該圖像塊的特征值。 相應地，若圖像塊的特征值由上述第一統(tǒng)計模塊統(tǒng)計得到，則水印嵌入單元704
進一步用于 若第四矩陣中的特征值相對于第一矩陣中的相應特征值的差值為正數(shù)，則確定存 在該差值的特征值對應的圖像塊為待修改的圖像塊，并根據(jù)該差值確定圖像塊的修改方式 為將與該差值相應數(shù)目的白像素點翻轉(zhuǎn)為黑像素點； 若第四矩陣中的特征值相對于第一矩陣中的相應特征值的差值為負數(shù)，確定存在 該差值的特征值對應的圖像塊為待修改的圖像塊，并根據(jù)該差值確定圖像塊的修改方式 為將與該差值相應數(shù)目的黑像素點翻轉(zhuǎn)為白像素點。 若圖像塊的特征值由上述第二統(tǒng)計模塊統(tǒng)計得到，則水印嵌入單元704進一步用 于 若第四矩陣中的特征值相對于第一矩陣中的相應特征值的差值為正數(shù)，確定存在 該差值的特征值對應的圖像塊為待修改的圖像塊，并根據(jù)該差值確定圖像塊的修改方式 為將與該差值相應數(shù)目的黑像素點翻轉(zhuǎn)為白像素點； 若第四矩陣中的特征值相對于第一矩陣中的相應特征值的差值為負數(shù)，確定存在 該差值的特征值對應的圖像塊為待修改的圖像塊，并根據(jù)該差值確定圖像塊的修改方式 為將與該差值相應數(shù)目的白像素點翻轉(zhuǎn)為黑像素點。 基于上述流程，本發(fā)明實施例提供了一種數(shù)字水印的提取裝置，如圖8所示，該裝 置包括 圖像劃分單元801、特征值統(tǒng)計單元802、矩陣處理單元803以及水印提取單元 804。具體地 圖像劃分單元801，用于根據(jù)與對圖像嵌入數(shù)字水印時相同的分塊規(guī)則將圖像劃 分為設(shè)定數(shù)目的圖像塊。 特征值統(tǒng)計單元802，用于統(tǒng)計圖像劃分單元801劃分得到的每個圖像塊的特征值。 矩陣處理單元803，用于根據(jù)與對圖像嵌入數(shù)字水印時相同的排列規(guī)則將特征值 統(tǒng)計單元802統(tǒng)計出的特征值組成第一矩陣，將組成的第一矩陣進行離散余弦變換得到第 二矩陣。 水印提取單元804，用于根據(jù)矩陣處理單元803得到的第二矩陣中的特征值確定
對圖像嵌入數(shù)字水印時用于替換指定特征值的閾值，并根據(jù)確定出的閾值以及設(shè)定的待嵌
入數(shù)字水印與閾值的對應關(guān)系，確定嵌入的數(shù)字水印。 進一步地，一個實施例中，上述特征值統(tǒng)計單元802可以包括 第一統(tǒng)計模塊，用于統(tǒng)計圖像塊包含的黑像素點數(shù)和圖像包含的總像素點數(shù)，并
將統(tǒng)計出的圖像塊的黑像素點數(shù)與總像素點數(shù)的比值確定為該圖像塊的特征值；
14
或者，第二統(tǒng)計模塊，用于統(tǒng)計圖像塊包含的白像素點數(shù)和圖像包含的總像素點
數(shù)，并將統(tǒng)計得到的白像素點數(shù)與總像素點數(shù)的比值確定為該圖像塊的特征值。 在對圖像嵌入數(shù)字水印時，首先在空域上對該圖像劃分并組織成空域特征值矩陣
(即第一矩陣)，進一步進行DCT得到頻域特征值矩陣(即第二矩陣)，修改該頻域特征值矩
陣中的指定特征值后，進行IDCT得到空域特征值矩陣(即第四矩陣)，并根據(jù)經(jīng)過上述變化
的空域特征值矩陣中的特征值相對于變化前的空域特征值矩陣中的相應特征值的差值，確
定待修改的圖像塊以及該圖像塊的修改方式，并基于該確定結(jié)果修改圖像。該方案充分利
用圖像在空域和頻域上的轉(zhuǎn)換關(guān)系，將圖像的空域特征和頻域特征相結(jié)合，使得嵌入的數(shù)
字水印的視覺隱藏效果好，并且抵抗打印、掃描等攻擊的魯棒性高。 在對圖像提取數(shù)字水印時，根據(jù)與對圖像嵌入數(shù)字水印時相同的分塊規(guī)則以及排 列規(guī)則，對圖像劃分并組織成空域特征值矩陣(即第一矩陣)，進一步進行DCT得到頻域特 征值矩陣(即第二矩陣)，并根據(jù)該頻域特征值矩陣中的特征值確定對圖像嵌入數(shù)字水印 時用于替換指定特征值的閾值，從而根據(jù)確定出的閾值以及設(shè)定的待嵌入數(shù)字水印與閾值 的對應關(guān)系，確定嵌入的數(shù)字水印。該方案充分考慮到頻域特征值矩陣相對穩(wěn)定的特征，根 據(jù)頻域特征值矩陣中的特征值提取嵌入的數(shù)字水印，準確度高。 顯然，本領(lǐng)域的技術(shù)人員可以對本發(fā)明進行各種改動和變型而不脫離本發(fā)明的精 神和范圍。這樣，倘若本發(fā)明的這些修改和變型屬于本發(fā)明權(quán)利要求及其等同技術(shù)的范圍 之內(nèi)，則本發(fā)明也意圖包含這些改動和變型在內(nèi)。
權(quán)利要求
一種數(shù)字水印的嵌入方法，其特征在于，包括根據(jù)分塊規(guī)則將圖像劃分為設(shè)定數(shù)目的圖像塊，統(tǒng)計每個圖像塊的特征值，并根據(jù)排列規(guī)則將統(tǒng)計出的特征值組成第一矩陣；將所述第一矩陣進行離散余弦變換得到第二矩陣，并根據(jù)預先定義的與數(shù)字水印對應的閾值替換所述第二矩陣中的指定特征值得到第三矩陣；將所述第三矩陣進行離散余弦逆變換得到第四矩陣；根據(jù)所述第四矩陣中的特征值相對于所述第一矩陣中的相應特征值的差值，確定待修改的圖像塊以及該圖像塊的修改方式，并根據(jù)確定結(jié)果修改所述圖像。
2. 如權(quán)利要求1所述的方法，其特征在于，統(tǒng)計圖像塊的特征值，包括 統(tǒng)計所述圖像塊包含的黑像素點數(shù)和所述圖像包含的總像素點數(shù)，并將所述黑像素點數(shù)與所述總像素點數(shù)的比值確定為所述圖像塊的特征值。
3. 如權(quán)利要求2所述的方法，其特征在于，當所述差值為正數(shù)時，所述確定待修改的圖 像塊以及該圖像塊的修改方式，包括確定存在所述差值的特征值對應的圖像塊為待修改的圖像塊，并根據(jù)所述差值確定所 述圖像塊的修改方式為將與所述差值相應數(shù)目的白像素點翻轉(zhuǎn)為黑像素點；當所述差值為負數(shù)時，所述確定待修改的圖像塊以及該圖像塊的修改方式，包括 確定存在所述差值的特征值對應的圖像塊為待修改的圖像塊，并根據(jù)所述差值確定所 述圖像塊的修改方式為將與所述差值相應數(shù)目的黑像素點翻轉(zhuǎn)為白像素點。
4. 如權(quán)利要求1所述的方法，其特征在于，統(tǒng)計圖像塊的特征值，包括 統(tǒng)計所述圖像塊包含的白像素點數(shù)和所述圖像包含的總像素點數(shù)，并將所述白像素點數(shù)與所述總像素點數(shù)的比值確定為所述圖像塊的特征值。
5. 如權(quán)利要求4所述的方法，其特征在于，當所述差值為正數(shù)時，所述確定待修改的圖 像塊以及該圖像塊的修改方式，包括確定存在所述差值的特征值對應的圖像塊為待修改的圖像塊，并根據(jù)所述差值確定所 述圖像塊的修改方式為將與所述差值相應數(shù)目的黑像素點翻轉(zhuǎn)為白像素點；當所述差值為負數(shù)時，所述確定待修改的圖像塊以及該圖像塊的修改方式，包括 確定存在所述差值的特征值對應的圖像塊為待修改的圖像塊，并根據(jù)所述差值確定所 述圖像塊的修改方式為將與所述差值相應數(shù)目的白像素點翻轉(zhuǎn)為黑像素點。
6. 如權(quán)利要求1所述的方法，其特征在于，所述排列規(guī)則包括所述圖像塊在所述圖像中的位置與所述圖像塊的特征值在所述第一矩陣中的位置之 間的對應關(guān)系。
7. —種數(shù)字水印的提取方法，其特征在于，包括根據(jù)與對圖像嵌入數(shù)字水印時相同的分塊規(guī)則將所述圖像劃分為設(shè)定數(shù)目的圖像塊， 統(tǒng)計每個圖像塊的特征值，并根據(jù)與對所述圖像嵌入數(shù)字水印時相同的排列規(guī)則將統(tǒng)計出 的特征值組成第一矩陣；將所述第一矩陣進行離散余弦變換得到第二矩陣，并根據(jù)所述第二矩陣中的特征值確 定對所述圖像嵌入數(shù)字水印時用于替換指定特征值的閾值；根據(jù)確定的所述閾值以及設(shè)定的待嵌入數(shù)字水印與閾值的對應關(guān)系，提取嵌入的數(shù)字 水印。
8. 如權(quán)利要求7所述的方法，其特征在于，統(tǒng)計圖像塊的特征值，包括 統(tǒng)計所述圖像塊包含的黑像素點數(shù)和所述圖像包含的總像素點數(shù)，并將所述黑像素點數(shù)與所述總像素點數(shù)的比值確定為所述圖像塊的特征值；或者，統(tǒng)計所述圖像塊包含的白像素點數(shù)和所述圖像包含的總像素點數(shù)，并將所述白 像素點數(shù)與所述總像素點數(shù)的比值確定為所述圖像塊的特征值。
9. 一種數(shù)字水印的嵌入裝置，其特征在于，包括圖像劃分單元，用于根據(jù)分塊規(guī)則將圖像劃分為設(shè)定數(shù)目的圖像塊； 特征值統(tǒng)計單元，用于統(tǒng)計圖像劃分單元劃分得到的每個圖像塊的特征值； 矩陣處理單元，用于根據(jù)排列規(guī)則將特征值統(tǒng)計單元統(tǒng)計出的特征值組成第一矩陣， 將所述第一矩陣進行離散余弦變換得到第二矩陣，根據(jù)預先定義的與數(shù)字水印對應的閾值 替換所述第二矩陣中的指定特征值得到第三矩陣，并將所述第三矩陣進行離散余弦逆變換 得到第四矩陣；水印嵌入單元，用于根據(jù)矩陣處理單元得到的第四矩陣中的特征值相對于第一矩陣中 的相應特征值的差值，確定待修改的圖像塊以及該圖像塊的修改方式，并根據(jù)確定結(jié)果修 改所述圖像。
10. 如權(quán)利要求9所述的裝置，其特征在于，所述特征值統(tǒng)計單元，包括 第一統(tǒng)計模塊，用于統(tǒng)計所述圖像塊包含的黑像素點數(shù)和所述圖像包含的總像素點數(shù)，并將所述黑像素點數(shù)與所述總像素點數(shù)的比值確定為所述圖像塊的特征值；或者，第二統(tǒng)計模塊，用于統(tǒng)計所述圖像塊包含的白像素點數(shù)和所述圖像包含的總像 素點數(shù)，并將所述白像素點數(shù)與所述總像素點數(shù)的比值確定為所述圖像塊的特征值。
11. 如權(quán)利要求io所述的裝置，其特征在于，若所述圖像塊的特征值由所述第一統(tǒng)計模塊統(tǒng)計得到，所述水印嵌入單元進一步用于若所述差值為正數(shù)，確定存在所述差值的特征值對應的圖像塊為待修改的圖像塊，并 根據(jù)所述差值確定所述圖像塊的修改方式為將與所述差值相應數(shù)目的白像素點翻轉(zhuǎn)為黑 像素點；若所述差值為負數(shù)，確定存在所述差值的特征值對應的圖像塊為待修改的圖像塊，并 根據(jù)所述差值確定所述圖像塊的修改方式為將與所述差值相應數(shù)目的黑像素點翻轉(zhuǎn)為白 像素點。
12. 如權(quán)利要求IO所述的裝置，其特征在于，若所述圖像塊的特征值由所述第二統(tǒng)計 模塊統(tǒng)計得到，所述水印嵌入單元進一步用于若所述差值為正數(shù)，確定存在所述差值的特征值對應的圖像塊為待修改的圖像塊，并 根據(jù)所述差值確定所述圖像塊的修改方式為將與所述差值相應數(shù)目的黑像素點翻轉(zhuǎn)為白 像素點；若所述差值為負數(shù)，確定存在所述差值的特征值對應的圖像塊為待修改的圖像塊，并 根據(jù)所述差值確定所述圖像塊的修改方式為將與所述差值相應數(shù)目的白像素點翻轉(zhuǎn)為黑 像素點。
13. —種數(shù)字水印的提取裝置，其特征在于，包括圖像劃分單元，用于根據(jù)與對圖像嵌入數(shù)字水印時相同的分塊規(guī)則將所述圖像劃分為 設(shè)定數(shù)目的圖像塊；特征值統(tǒng)計單元，用于統(tǒng)計圖像劃分單元劃分得到的每個圖像塊的特征值； 矩陣處理單元，用于根據(jù)與對所述圖像嵌入數(shù)字水印時相同的排列規(guī)則將特征值統(tǒng)計單元統(tǒng)計出的特征值組成第一矩陣，將所述第一矩陣進行離散余弦變換得到第二矩陣；水印提取單元，用于根據(jù)矩陣處理單元得到的第二矩陣中的特征值確定對所述圖像嵌入數(shù)字水印時用于替換指定特征值的閾值，并根據(jù)確定的所述閾值以及設(shè)定的待嵌入數(shù)字水印與閾值的對應關(guān)系，提取嵌入的數(shù)字水印。
14.如權(quán)利要求13所述的裝置，其特征在于，所述特征值統(tǒng)計單元，包括 第一統(tǒng)計模塊，用于統(tǒng)計所述圖像塊包含的黑像素點數(shù)和所述圖像包含的總像素點數(shù)，并將所述黑像素點數(shù)與所述總像素點數(shù)的比值確定為所述圖像塊的特征值；或者，第二統(tǒng)計模塊，用于統(tǒng)計所述圖像塊包含的白像素點數(shù)和所述圖像包含的總像素點數(shù)，并將所述白像素點數(shù)與所述總像素點數(shù)的比值確定為所述圖像塊的特征值。
全文摘要
本發(fā)明公開了數(shù)字水印的嵌入方法、提取方法及裝置，包括嵌入水印時，根據(jù)分塊規(guī)則將圖像劃分并組成第一矩陣，進行離散余弦變換得到第二矩陣，替換第二矩陣中的指定特征值得到第三矩陣，進行離散余弦逆變換得到第四矩陣，根據(jù)第四矩陣中的特征值相對于第一矩陣中的相應特征值的差值，確定待修改的圖像塊以及該圖像塊的修改方式，并根據(jù)確定結(jié)果修改圖像。提取水印時，根據(jù)與嵌入水印時相同的分塊規(guī)則將圖像劃分并組成第一矩陣，進行離散余弦變換得到第二矩陣，根據(jù)第二矩陣確定對圖像嵌入水印時用于替換指定特征值的閾值，進而提取數(shù)字水印。根據(jù)該方案，提高了圖像中數(shù)字水印信息的隱藏性和魯棒性，以及提高了圖像中提取數(shù)字水印信息的準確度。
文檔編號G06T1/00GK101751655SQ20081024031
公開日2010年6月23日 申請日期2008年12月17日 優(yōu)先權(quán)日2008年12月17日
發(fā)明者于權(quán), 吳於茜, 和斌, 康凱 申請人:北京大學;北大方正集團有限公司;北京方正電子政務信息科技有限公司