本發(fā)明涉及圖像處理領域，具體而言，涉及一種數字圖像的半色調處理方法和系統。

背景技術：

數字圖像半色調處理是利用有限的色調對原始連續(xù)色調圖像進行表示的技術，利用了人類視覺系統低通濾波的屬性，使得半色調處理的結果與原始圖像在感官上比較類似，但是，目前如何提升對數字圖像進行處理的效率和處理效果，成為了亟待解決的問題。

相關技術中，主要是基于優(yōu)化迭代的圖像半色調算法，以與輸入圖像整體灰度值相同的二值噪聲圖像為初始半色調結果，然后隨機交換當前結果中的黑白像素點，判斷新的結果在色調、結構、對比度上是否與輸入圖像更接近，當滿足判斷條件時，則交換兩個像素點；否則，保持現有像素點不變繼續(xù)測試交換其他黑白像素點。

但是，現有技術中的方法，黑白像素點交換過程是通過模擬退火算法進行優(yōu)化直至收斂的，該方法對圖像處理時需要時間比較長，導致對圖像的處理效率比較低，并且圖像的處理效果不能很好的保持原圖像的特征。

技術實現要素：

有鑒于此，本發(fā)明實施例的目的在于提供一種數字圖像的半色調處理方法和系統，以實現有效的保持輸入圖像的亮度和對比度信息，同時節(jié)約了對圖像處理的時間。

第一方面，本發(fā)明實施例提供了一種數字圖像的半色調處理方法，所述方法包括：

獲取待處理圖像，并對所述待處理圖像進行分析得到所述待處理圖像的特征信息，其中，所述特征信息包括：像素點的局部亮度信息、對比度信息和所述待處理圖像的結構信息；

根據所述像素點的局部亮度信息和對比度信息，計算所述待處理圖像的量化閾值；

計算各所述像素點的量化誤差，并將所述量化誤差擴散到各所述像素點的鄰域內，計算量化誤差擴散后的各所述像素點的量化誤差灰度值；

根據所述量化閾值和各所述像素點的量化誤差灰度值，對所述待處理圖像進行二值化，得到二值化后所述待處理圖像中各像素點的二值化灰度值，并對所述二值化后的像素點的所述二值化灰度值進行誤差傳遞，得到半色調處理圖像。

結合第一方面，本發(fā)明實施例提供了第一方面的第一種可能的實施方式，其中，獲取待處理圖像，并對所述待處理圖像進行分析得到所述待處理圖像的特征信息，包括：

計算所述待處理圖像的梯度，并根據所述待處理圖像的梯度得到所述待處理圖像的結構信息；

根據所述待處理圖像的結構信息，通過以下公式計算所述像素點的局部亮度信息和對比度信息：

其中，w_i,j表示鄰域內像素點的權重，σ＝1.5，L[m,n]表示所述像素點的局部亮度信息，C[m,n]表示對比度信息，m和n分別表示像素點的橫縱坐標，f[m+i,n+j]表示所述待處理圖像的初始灰度值。

結合第一方面，本發(fā)明實施例提供了第一方面的第二種可能的實施方式，其中，根據所述像素點的局部亮度信息和對比度信息，計算所述待處理圖像的量化閾值，包括：

通過以下公式計算所述量化閾值：

t[m,n]＝127.5-f[m,n]×C[m,n]×(f[m,n]-L[m,n])，

其中，t[m,n]表示量化閾值，127.5為量化閾值的默認閾值，L[m,n]表示所述像素點的局部亮度信息，C[m,n]表示對比度信息，m和n分別表示像素點的橫縱坐標，f[m,n]表示所述待處理圖像的初始灰度值。

結合第一方面，本發(fā)明實施例提供了第一方面的第三種可能的實施方式，其中，在所述獲取待處理圖像，并對所述待處理圖像進行分析得到所述待處理圖像的特征信息步驟之后，所述方法還包括：根據所述待處理圖像的結構信息計算誤差擴散系數；

其中，通過以下公式計算誤差擴散系數：

其中，s和t分別表示目標像素點的橫縱坐標，w[s,t]表示誤差擴散系數，d_st為當前像素點與目標像素點之間的距離，I_st表示目標像素點的灰度值，表示結構信息對d_st調整后的指數。

結合第一方面，本發(fā)明實施例提供了第一方面的第四種可能的實施方式，其中，計算各所述像素點的量化誤差，并將所述量化誤差擴散到各所述像素點的鄰域內，計算量化誤差擴散后的各所述像素點的量化誤差灰度值，包括：

通過以下公式計算各像素點的量化誤差灰度值：

u[m,n]＝f[m,n]+∑_s,tw[s,t]×e[m-s,n-t]，

其中，m和n分別表示像素點的橫縱坐標，s和t分別表示目標像素點的橫縱坐標，u[m,n]表示量化誤差灰度值，f[m,n]表示所述待處理圖像的初始灰度值，w[s,t]表示誤差擴散系數，e[m-s,n-t]表示量化誤差。

結合第一方面，本發(fā)明實施例提供了第一方面的第五種可能的實施方式，其中，根據所述量化閾值和各所述像素點的量化誤差灰度值，對所述待處理圖像進行二值化，得到二值化后所述待處理圖像中各像素點的二值化灰度值，并對所述二值化后的像素點的所述二值化灰度值進行誤差傳遞，得到半色調處理圖像，包括：

判斷各所述像素點的量化誤差灰度值是否大于所述量化閾值；

如果是，則確認所述像素點為白色；

如果否，則確認所述像素點為黑色；

根據所述誤差擴散系數，將二值化后產生的新的量化誤差擴散到鄰域內未處理的像素點，得到半色調處理后的圖像。

第二方面，本發(fā)明實施例提供了一種數字圖像的半色調處理系統，所述系統包括：

特征信息獲取模塊，用于獲取待處理圖像，并對所述待處理圖像進行分析得到所述待處理圖像的特征信息，其中，所述特征信息包括：像素點的局部亮度信息、對比度信息和所述待處理圖像的結構信息；

量化閾值計算模塊，根據所述像素點的局部亮度信息和對比度信息，計算所述待處理圖像的量化閾值；

量化誤差灰度值計算模塊，用于計算各所述像素點的量化誤差，并將所述量化誤差擴散到各所述像素點的鄰域內，計算量化誤差擴散后的各所述像素點的量化誤差灰度值；

半色調處理模塊，用于根據所述量化閾值和各所述像素點的量化誤差灰度值，對所述待處理圖像進行二值化，得到二值化后所述待處理圖像中各像素點的二值化灰度值，并對所述二值化后的像素點的所述二值化灰度值進行誤差傳遞，得到半色調處理圖像。

結合第二方面，本發(fā)明實施例提供了第二方面的第一種可能的實施方式，其中，所述系統還包括：

誤差擴散系數計算模塊，用于通過以下公式計算誤差擴散系數：

其中，s和t分別表示目標像素點的橫縱坐標，w[s,t]表示誤差擴散系數，d_st為當前像素點與目標像素點之間的距離，I_st表示目標像素點的灰度值，表示結構信息對d_st調整后的指數。

結合第二方面，本發(fā)明實施例提供了第二方面的第二種可能的實施方式，其中，所述量化誤差灰度值計算模塊包括：

量化誤差灰度值計算單元，用于通過以下公式計算各像素點的量化誤差灰度值：

u[m,n]＝f[m,n]+∑_s,tw[s,t]×e[m-s,n-t]，

其中，m和n分別表示像素點的橫縱坐標，s和t分別表示目標像素點的橫縱坐標，u[m,n]表示量化誤差灰度值，f[m,n]表示所述待處理圖像的初始灰度值，w[s,t]表示誤差擴散系數，e[m-s,n-t]表示量化誤差。

結合第二方面，本發(fā)明實施例提供了第二方面的第三種可能的實施方式，其中，所述半色調處理模塊包括：

判斷模塊，用于判斷各所述像素點的量化誤差灰度值是否大于所述量化閾值；

第一確認模塊，用于當所述判斷模塊的判斷結果為是時，則確認所述像素點為白色；

第二確認模塊，用于當所述判斷模塊的判斷結果為否時，則確認所述像素點為黑色；

半色調處理單元，用于根據所述誤差擴散系數，將二值化后產生的新的量化誤差擴散到鄰域內未處理的像素點，得到半色調處理后的圖像。

本發(fā)明實施例提供的一種數字圖像的半色調處理方法和系統，通過對待處理圖像的像素點的局部亮度信息、對比度信息和結構信息進行提取，并根據局部亮度信息和對比度信息計算出量化閾值，根據量化閾值和各像素點的量化誤差灰度值對待處理圖像進行二值化，并對二值化后的像素點進行誤差傳遞，從而得到半色調處理圖像，該方法有效的保持了輸入圖像的亮度和對比度信息，同時節(jié)約了對圖像處理的時間。

附圖說明

為了更清楚地說明本發(fā)明具體實施方式或現有技術中的技術方案，下面將對具體實施方式或現有技術描述中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖是本發(fā)明的一些實施方式，對于本領域普通技術人員來講，在不付出創(chuàng)造性勞動的前提下，還可以根據這些附圖獲得其他的附圖。

圖1示出了本發(fā)明實施例所提供的一種數字圖像的半色調處理方法的流程圖；

圖2示出了本發(fā)明實施例所提供的一種數字圖像的半色調處理方法的另一流程圖；

圖3示出了本發(fā)明實施例所提供的一種數字圖像的半色調處理系統的結構示意圖；

圖4示出了本發(fā)明實施例所提供的半色調處理模塊的結構示意圖。

圖標：

附圖3中，各標號所代表的部件列表如下：

30-特征信息獲取模塊； 31-量化閾值計算模塊；

32-量化誤差灰度值計算模塊； 33-半色調處理模塊。

附圖4中，各標號所代表的部件列表如下：

40-判斷模塊； 41-第一確認模塊；

42-第二確認模塊； 43-半色調處理單元。

具體實施方式

下面將結合本發(fā)明實施例中附圖，對本發(fā)明實施例中的技術方案進行清楚、完整地描述，顯然，所描述的實施例僅僅是本發(fā)明一部分實施例，而不是全部的實施例。通常在此處附圖中描述和示出的本發(fā)明實施例的組件可以以各種不同的配置來布置和設計。因此，以下對在附圖中提供的本發(fā)明的實施例的詳細描述并非旨在限制要求保護的本發(fā)明的范圍，而是僅僅表示本發(fā)明的選定實施例。基于本發(fā)明的實施例，本領域技術人員在沒有做出創(chuàng)造性勞動的前提下所獲得的所有其他實施例，都屬于本發(fā)明保護的范圍。

考慮到相關技術中，主要是基于優(yōu)化迭代的圖像半色調算法，以與輸入圖像整體灰度值相同的二值噪聲圖像為初始半色調結果，然后隨機交換當前結果中的黑白像素點，判斷新的結果在色調、結構、對比度上是否與輸入圖像更接近，當滿足判斷條件時，則交換兩個像素點；否則，保持現有像素點不變繼續(xù)測試交換其他黑白像素點。但是，現有技術中的方法，黑白像素點交換過程是通過模擬退火算法進行優(yōu)化直至收斂的，該方法對圖像處理時需要時間比較長，導致對圖像的處理效率比較低，并且圖像的處理效果不能很好的保持原圖像的特征?；诖耍景l(fā)明實施例提供了一種數字圖像的半色調處理方法和系統，下面通過實施例進行描述。

實施例1

本實施例提供的數字圖像的半色調處理方法，有效的保持了輸入圖像的亮度和對比度信息，同時節(jié)約了對圖像處理的時間，提高了圖像處理的效率。

參見圖1，本實施例提供了一種數字圖像的半色調處理方法，該方法包括：

步驟S101，獲取待處理圖像，并對待處理圖像進行分析得到待處理圖像的特征信息，其中，特征信息包括：像素點的局部亮度信息、對比度信息和待處理圖像的結構信息。

參見圖2，步驟S101進一步地包括步驟S201至步驟S202：

步驟S201,計算待處理圖像的梯度，并根據待處理圖像的梯度得到待處理圖像的結構信息；

步驟S202,根據待處理圖像的結構信息，通過以下公式計算像素點的局部亮度信息和對比度信息：

其中，w_i,j表示鄰域內像素點的權重，σ＝1.5，L[m,n]表示像素點的局部亮度信息，C[m,n]表示對比度信息，m和n分別表示像素點的橫縱坐標，f[m+i,n+j]表示待處理圖像的初始灰度值。

其中，計算輸入圖像每個像素的局部亮度信息：以每個像素為中心，取其周圍7x7鄰域內像素灰度的平均值作為其局部亮度信息，鄰域內像素的權重w由歸一化后的二維高斯函數確定，即

σ＝1.5為高斯函數的標準差。

計算輸入圖像每個像素的局部對比度信息：取像素周圍7x7的鄰域，計算鄰域內像素的灰度值與當前像素的局部亮度之間的平方差，利用平方差的平均值作為當前像素的局部對比度信息，即

其中w_i,j為歸一化的二維高斯函數。

步驟S102，根據像素點的局部亮度信息和對比度信息，計算待處理圖像的量化閾值。

步驟S102具體包括下列步驟(1)：

(1)通過以下公式計算量化閾值：

t[m,n]＝127.5-f[m,n]×C[m,n]×(f[m,n]-L[m,n])，

其中，t[m,n]表示量化閾值，127.5為量化閾值的默認閾值，L[m,n]表示所述像素點的局部亮度信息，C[m,n]表示對比度信息，m和n分別表示像素點的橫縱坐標，f[m,n]表示所述待處理圖像的初始灰度值。

利用計算得到的像素局部亮度信息和對比度信息對誤差擴散的量化閾值進行調節(jié)；假設輸入灰度圖像的每個像素用8個bit進行表示，其灰度值取值范圍為0-255，因此默認二值量化的閾值為127.5，因此，引入像素的局部亮度和對比度信息后，調節(jié)以后的量化閾值為：

t[m,n]＝127.5-f[m,n]×C[m,n]×(f[m,n]-L[m,n])。

步驟S103，計算各像素點的量化誤差，并將量化誤差擴散到各像素點的鄰域內，計算量化誤差擴散后的各像素點的量化誤差灰度值。

步驟S103具體包括：通過以下公式計算各像素點的量化誤差灰度值：

u[m,n]＝f[m,n]+∑_s,tw[s,t]×e[m-s,n-t]，

其中，m和n分別表示像素點的橫縱坐標，s和t分別表示目標像素點的橫縱坐標，u[m,n]表示量化誤差灰度值，f[m,n]表示所述待處理圖像的初始灰度值，w[s,t]表示誤差擴散系數，e[m-s,n-t]表示量化誤差。

量化誤差灰度值u[m,n]包括兩部分內容，一是初始輸入圖像的灰度值f[m,n]，另一就是擴散到當前像素的量化誤差e，從而得到經過量化誤差擴散的當前像素量化誤差灰度值u[m,n]。

步驟S104，根據量化閾值和各像素點的量化誤差灰度值，對待處理圖像進行二值化，得到二值化后待處理圖像中各像素點的二值化灰度值，并對二值化后的像素點的二值化灰度值進行誤差傳遞，得到半色調處理圖像。

步驟S104進一步地包括：步驟(1)至步驟(4)：

(1)判斷各像素點的量化誤差灰度值是否大于量化閾值，當判斷結果為是時，執(zhí)行步驟(1)，當判斷結果為否時，執(zhí)行步驟(3)；

(2)確認像素點為白色；

(3)確認像素點為黑色；

(4)根據誤差擴散系數，將二值化后產生的新的量化誤差擴散到鄰域內未處理的像素點，得到半色調處理后的圖像。

根據量化閾值和各經過量化誤差擴散后的像素點的量化誤差灰度值，對該待處理圖像進行二值化，具體過程如下：

按照掃描線順序對待處理圖像像素點進行二值量化：在處理某個像素點時，比較當前像素點的量化誤差灰度值u[m,n]與閾值t[m,n]之間的大小，并根據大小關系判斷在輸出結果中當前像素應該為黑色還是白色，即當u[m,n]≥t[m,n]二值化灰度值為g[m,n]＝255，否則則為g[m,n]＝0。

為了保持圖像的整體色調，當前像素量化操作完成以后，即二值化完成之后，需要將產生的新的量化誤差e[m,n]＝u[m,n]-g[m,n]擴散到鄰域內未處理像素。

除此之外，該方法還包括：根據待處理圖像的結構信息計算誤差擴散系數；

其中，通過以下公式計算誤差擴散系數：

其中，s和t分別表示目標像素點的橫縱坐標，w[s,t]表示誤差擴散系數，d_st為當前像素點與目標像素點之間的距離，I_st表示目標像素點的灰度值，表示結構信息對d_st調整后的指數。

為了使鄰域內的像素接受不同程度的量化誤差，引入了擴散系數w[s,t]，擴散系數的設置遵從原則：與距離較遠的像素比臨近的像素應該吸收較多的誤差。為了保證圖像的色彩關系：量化誤差為正時，灰度值較大的像素應該吸收較多的誤差；量化誤差為負時，灰度值較小的像素應該吸收較多的誤差。因此，擴散系數為

其中d_st為誤差擴散源像素與目標像素之間的距離；量化誤差為正時，I_st＝f[s,t]，為負時，I_st＝255-f[s,t]。此外，為了保持低對比度區(qū)域的結構信息，利用圖像結構信息對擴散系數中距離d_st的指數進行調整。調整后的指數k為：

綜上所述，本實施例提供了一種數字圖像的半色調處理方法，通過對待處理圖像的像素點的局部亮度信息、對比度信息和結構信息進行提取，并根據局部亮度信息和對比度信息計算出量化閾值，根據量化閾值和各像素點的量化誤差灰度值對待處理圖像進行二值化，并對二值化后的像素點進行誤差傳遞，從而得到半色調處理圖像，該方法有效的保持了輸入圖像的亮度和對比度信息，同時節(jié)約了對圖像處理的時間。

實施例2

參見圖3，本實施例提供了一種數字圖像的半色調處理系統，該系統包括：特征信息獲取模塊30，用于獲取待處理圖像，并對待處理圖像進行分析得到待處理圖像的特征信息，其中，特征信息包括：像素點的局部亮度信息、對比度信息和待處理圖像的結構信息；

量化閾值計算模塊31，根據像素點的局部亮度信息和對比度信息，計算待處理圖像的量化閾值；

量化誤差灰度值計算模塊32，用于計算各像素點的量化誤差，并將量化誤差擴散到各像素點的鄰域內，計算量化誤差擴散后的各像素點的量化誤差灰度值；

半色調處理模塊33，用于根據量化閾值和各像素點的量化誤差灰度值，對待處理圖像進行二值化，得到二值化后待處理圖像中各像素點的二值化灰度值，并對二值化后的像素點的二值化灰度值進行誤差傳遞，得到半色調處理圖像。

該系統還包括：誤差擴散系數計算模塊，用于通過以下公式計算誤差擴散系數：

其中，s和t分別表示目標像素點的橫縱坐標，w[s,t]表示誤差擴散系數，d_st為當前像素點與目標像素點之間的距離，I_st表示目標像素點的灰度值，表示結構信息對d_st調整后的指數。

量化誤差灰度值計算模塊32包括：

量化誤差灰度值計算單元，用于通過以下公式計算各像素點的量化誤差灰度值：

u[m,n]＝f[m,n]+∑_s,tw[s,t]×e[m-s,n-t]，

其中，其中，m和n分別表示像素點的橫縱坐標，s和t分別表示目標像素點的橫縱坐標，u[m,n]表示量化誤差灰度值，f[m,n]表示所述待處理圖像的初始灰度值，w[s,t]表示誤差擴散系數，e[m-s,n-t]表示量化誤差。

參見圖4，半色調處理模塊33包括：

判斷模塊40，用于判斷各像素點的量化誤差灰度值是否大于量化閾值；

第一確認模塊41，用于當判斷模塊的判斷結果為是時，則確認像素點為白色；

第二確認模塊42，用于當判斷模塊的判斷結果為否時，則確認像素點為黑色；

半色調處理單元43，用于根據誤差擴散系數，將二值化后產生的新的量化誤差擴散到鄰域內未處理的像素點，得到半色調處理后的圖像。

綜上所述，本實施例提供了一種數字圖像的半色調處理系統，通過對待處理圖像的像素點的局部亮度信息、對比度信息和結構信息進行提取，并根據局部亮度信息和對比度信息計算出量化閾值，根據量化閾值和各像素點的量化誤差灰度值對待處理圖像進行二值化，并對二值化后的像素點進行誤差傳遞，從而得到半色調處理圖像，該方法有效的保持了輸入圖像的亮度和對比度信息，同時節(jié)約了對圖像處理的時間。

本發(fā)明實施例所提供的數字圖像的半色調處理方法的計算機程序產品，包括存儲了程序代碼的計算機可讀存儲介質，程序代碼包括的指令可用于執(zhí)行前面方法實施例中的方法，具體實現可參見方法實施例，在此不再贅述。

所屬領域的技術人員可以清楚地了解到，為描述的方便和簡潔，上述描述的系統、裝置和單元的具體工作過程，可以參考前述方法實施例中的對應過程，在此不再贅述。

以上所述，僅為本發(fā)明的具體實施方式，但本發(fā)明的保護范圍并不局限于此，任何熟悉本技術領域的技術人員在本發(fā)明揭露的技術范圍內，可輕易想到變化或替換，都應涵蓋在本發(fā)明的保護范圍之內。因此，本發(fā)明的保護范圍應所述以權利要求的保護范圍為準。