CN103428402A - 产生半色调影像的方法与影像处理系统及电脑产品 - Google Patents
产生半色调影像的方法与影像处理系统及电脑产品 Download PDFInfo
- Publication number
- CN103428402A CN103428402A CN2013101752754A CN201310175275A CN103428402A CN 103428402 A CN103428402 A CN 103428402A CN 2013101752754 A CN2013101752754 A CN 2013101752754A CN 201310175275 A CN201310175275 A CN 201310175275A CN 103428402 A CN103428402 A CN 103428402A
- Authority
- CN
- China
- Prior art keywords
- mentioned
- quantization error
- image
- tremble
- error value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N1/00—Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
- H04N1/40—Picture signal circuits
- H04N1/405—Halftoning, i.e. converting the picture signal of a continuous-tone original into a corresponding signal showing only two levels
- H04N1/4051—Halftoning, i.e. converting the picture signal of a continuous-tone original into a corresponding signal showing only two levels producing a dispersed dots halftone pattern, the dots having substantially the same size
- H04N1/4052—Halftoning, i.e. converting the picture signal of a continuous-tone original into a corresponding signal showing only two levels producing a dispersed dots halftone pattern, the dots having substantially the same size by error diffusion, i.e. transferring the binarising error to neighbouring dot decisions
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
Abstract
本发明公开了一种产生半色调影像的方法、影像处理系统及电脑产品,用于一影像处理系统以产生一半色调影像,上述方法包括:通过一量化器接收一具有原始像素值的原始影像,并根据一预设临界值产生一具有量化像素值的量化影像;通过一减法器将上述原始像素值减去上述量化像素值,以产生一量化误差值;通过一后抖色模块并依据上述量化误差值执行一后抖色程序,以产生一抖色值;以及通过一加法器加总上述量化影像的上述量化像素值及上述抖色值,并产生一具有半色调像素值的半色调影像。
Description
技术领域
本发明有关于一种利用产生影像的方法、影像处理系统及电脑程序产品,且特别是有关于一种用以转换一灰阶值影像至一半色调影像的方法、影像处理系统及电脑产品。
背景技术
数字半色调是一种利用数字影像处理由一输入数字影像中产生一半色调数字影像的技术。输入数字影像一般由范围由0至255的离散像素值所组成。为了在一可显示灰阶(例如,黑色)点(dot)的输出装置中重现此影像,必须转换此输入数字影像至一些使用半色调技术形式的半色调数字影像。半色调方法依据观察者的视线,通过半色调数字影像中的一些局部区域在空间中平衡视觉,使得中间的灰阶值可以通过在半色调数字影像的一些小区域“开启”一些像素值及“关闭”一些像素值而产生。部分被开启的像素值将决定显示的灰阶值。
现有的数字半色调技术可分为三大类:(1)抖色法(dithering)、(2)扩散误差法(error diffusion)以及(3)迭代最佳化(Iterative Optimization)。这些技术以及这些技术的组合皆有其自身的优点和缺点。一般来说,迭代最佳化技术要求较高的计算量,其主要用于学术方面。抖色及扩散误差此两类较实际的做法已被广泛地研究。这三种方法简要的讲述如下。
(1)抖色法
有序抖色法(Ordered dithering)可分为两类。通常分为点聚集(clustered-dot)与点扩散(dispersed-dot)有序抖色法两大类。图1a显示一点聚集有序抖色图式,图1b显示一点扩散有序抖色图式。点聚集有序抖色法在一固定间距内使用可变大小的半色调点。在点的外边缘设备像素的增加可以提高覆盖区域及点的大小。从远处看时,点的尺寸越大,覆盖区域也越大,而影像区域较暗。当显示装置可显示一独立黑色或白色像素时,点扩散有序抖色法可为较佳的选择。其使用了一固定大小、可变间距较小的点,以达到与点聚集有序抖色法相同的效果。在一给定的区域或点频率中,点间距的变化改变点的数量。在此技术中,更密集的点分布提供较暗的影像区域。在一些显示装置中,每个点包括四或五个装置像素。点扩散有序抖色法根据在原始影像中阴影的变化提供点的分布。对于特定的显示装置而言,最佳化点分布可能为最好的显示方式。
(2)错误扩散法
错误扩散法是一种依据输入数字影像值产生不同空间频率内容的适应性演算法(adaptive algorithm)。图2显示一说明基本错误扩散技术的先前技术方块图200。此技术更详细揭露于“An Adaptive Algorithm for SpatialGreyscale”Proceedings of the Society for Information Display,volume17,pp.75,1976by R.W.Floyd and L.Steinberg中。为了方便说明,将假设输入数字影像的像素范围由0至255。如图2中所示,在一临界值方块201中限制一输入数字影像Pi的像素值,以产生临界的像素值。临界值方块对于输入数字影像中低于一临界值的像素提供一信号使其像素值变更为0,对于输入数字影像中高于此临界值的像素提供一信号使其像素值变更为255。一信号差产生器202接收由临界值方块201中所传送的输入数字影像的像素值及一加法器204的输出。信号差产生器202产生表示为误差205的一信号差。信号差在一误差滤波器203中乘上一误差反馈权重,并提供给加法器204以将此权重误差加上尚未处理的邻近像素值。此可确保半色调数字影像的像素值的算术平均值被保存在一局部的影像区域中。
(3)迭代最佳化
迭代最佳化方法根据一些相关的模型,例如人类视觉系统(human visualsystem,HVS),以减少在连续色调影像及半色调影像间的误差。此错误通常通过加权最小平方误差法(Weighted Least-Square Approaches)来计算。由此种类型技术所推导的半色调影像通常具有较高的品质,但在计算复杂度较高。迭代最佳化的过程如图3所示。
当一电泳显示器(electrophorotic display,EPD)的空间解析度远低于一印表机时,在三种方法间,有序抖色法具有最低的计算量,但所产生的影像品质最差。而错误扩散法和迭代最佳化提供较佳的影像。然而,当使用手持式电子装置来执行这些方法时,其所需的计算资源将成为主要的问题。举例来说,使用8位递回运算来处理的错误扩散法,其计算复杂度在一电泳显示器系统中将是一个很大的负担。因此,一种更进一步减少计算复杂度,并提供高影像品质的系统与方法是有待开发的。
发明内容
本发明提供一种产生半色调影像的方法、影像处理系统及电脑程序产品。
本发明提出一种产生半色调影像的方法,用于一影像处理系统以产生一半色调影像,上述方法包括:通过一量化器接收一具有原始像素值的原始影像,并根据一预设临界值产生一具有量化像素值的量化影像;通过一减法器将上述原始像素值减去上述量化像素值,以产生一量化误差值;通过一后抖色模块并依据上述量化误差值执行一后抖色程序,以产生一抖色值;以及通过一加法器加总上述量化影像的上述量化像素值及上述抖色值,并产生一具有半色调像素值的半色调影像。
本发明提出一种影像处理系统,用以产生一半色调影像,上述系统包括:一量化器,接收一具有原始像素值的原始影像,并根据一预设临界值产生一具有量化像素值的量化影像;一减法器,耦接至上述量化器,用以将上述原始像素值减去上述量化像素值,以产生一量化误差值;一后抖色模块,耦接至上述减法器,用以依据上述量化误差值执行一后抖色程序,以产生一抖色值;以及一加法器,耦接至上述后抖色模块及上述量化器,用以加总上述量化影像的上述量化像素值及上述抖色值,并产生一具有半色调像素值的半色调影像。
本发明提出一种电脑程序产品,用以被一电子装置载入以于一影像处理系统中执行产生半色调影像的方法,上述影像处理系统,上述方法包括:接收一具有原始像素值的原始影像,并根据一预设临界值产生一具有量化像素值的量化影像;将上述原始像素值减去上述量化像素值,以产生一量化误差值;依据上述量化误差值执行一后抖色程序,以产生一抖色值;以及加总上述量化影像的上述量化像素值及上述抖色值,并产生一具有半色调像素值的半色调影像。
附图说明
图1a显示一4×4点聚集有序抖色图式。
图1b显示一4×4点扩散有序抖色图式。
图2显示说明执行错误扩散半色调影像的系统方块图。
图3显示一执行迭代最佳化半色调影像的系统方块图。
图4显示根据本发明一实施例用于产生一半色调影像的影像处理系统的方块图。
图5显示根据本发明一实施例用于产生抖色矩阵的虚拟程序码部分的流程图。
图6a~6c显示根据本发明一实施例的原始影像。
图6d~6f显示根据本发明一实施例的量化影像。
图6g~6i显示根据本发明一实施例的半色调影像。
其中,附图标记说明如下:
201~临界值方块;
202~信号差产生器;
203~误差滤波器;
204~加法器;
205~误差;
400~影像处理系统;
410~量化器;
420~减法器;
430~后抖色模块;
440~加法器;
S501、S502、S503、S504、S505、S506、S507、S508、S509、S510~步骤;
x~原始影像;
xq~量化影像;
e~量化误差值;
d~抖色矩阵;
y~半色调影像。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合所附图式,作详细说明如下。
I.系统架构
本发明是一种有关于具有高执行效率的后抖色(post-dithering)技术的影像处理方法,用以将原始灰阶影像转换为半色调(Halftone)影像。
图4显示根据本发明一实施例用于产生一半色调影像的影像处理系统400的方块图。如图所示,影像处理系统400处理一原始影像x以产生一半色调影像。影像处理系统400包括一量化器410、一减法器420、一后抖色模块430及一加法器440。
量化器410具有第一输入端,用以接收具有原始像素值的原始影像x,并根据一预设临界值T产生一具有量化像素值的量化影像xq。
减法器420耦接至量化器410,并具有第一输入端,用以接收原始影像x。减法器420将原始影像x的原始像素值减去量化影像xq的量化像素值,以产生一量化误差值e。
后抖色模块430耦接至减法器420,用以依据量化误差值e执行一后抖色程序,以产生具有一抖色值的一抖色矩阵d。
加法器440耦接至后抖色模块430及量化器410。加法器440用以加总量化影像的量化像素值xq及抖色矩阵d的抖色值,并产生一具有半色调像素值的半色调影像y。
在此一实施例中,本发明所提出包括一演算法的后抖色程序,强调将半色调影像的视觉误差的均方差(Mean Square Error,MSE)最小化。在数学上,视觉误差ev的均方差最小化可以下列公式表示:
1s ubject to d(i,j)={0,±T}
在此实施例中,假设原始影像尺寸为I×J,而原始影像x由二维的像素阵列(i,j)所组成,其中i及j皆为正整数,表示原始影像x的第i行及第j列。视觉误差ev为人眼v的调变转换函数(Modulation Transfer Function,MTF)及原始影像x和半色调图像y之间误差的回旋积(Convolution)总和。而抖色矩阵d的抖色值应限于相邻灰阶值之间的范围T内。
然而,由于人眼v的调变转换函数为一低通滤波器,而视觉误差ev在低频下易由眼睛察觉,但在高频下并不易察觉。因此,修改后的演算法注重于将一局部区域(Local Area)内的错误最小化,以取代的人眼v的调变转换函数。原始影像x将被分成许多非重迭、大小为M×N的局部区域,每一局部区域为一包括二维像素阵列(m,n)的矩阵,其中m及n皆为正整数且分别代表上述矩阵的第m行及第n列。公式(1)可以修改为:
subject to d(m,n)={0,±T}
其中elocal为第M×N个局部区域的误差。通过将局部误差最小化的调变转换函数,公式(2)可降低在低频中的错误数目。
II.后抖色程序(Post-dithering Process)
描述公式(2)的虚拟程序码(pseudo code)部分可由数值迭代的方法实现如下。
while k≤M×N
break;
(m,n)=max(xlocal+dlocal);
dlocal(m,n)=-T;
elocal=elocal-T;
(m,n)=max(xlocal+dlocal);
dlocal(m,n)=T;
elocal=elocal+T;
end
k=k+1;
end
请参阅虚拟程序码的部分和图5。图5显示根据本发明一实施例用于产生抖色矩阵dlocal的虚拟程序码部分的流程图。本发明所提出的后抖色程序可通过在一抖色矩阵dlocal中分配欲抖色的点(dot),以逐渐减少在原始影像xlocal的局部区域中的局部误差elocal,并减少步阶系为在相邻灰阶间的预定临界值T。
首先,给定原始影像xlocal的一局部区域,令一计数值k的初始值为1及一抖色矩阵dlocal初始值为一零矩阵(Null Matrix)。在步骤S501中,判断计数值k是否小于或等于M×N的数值(k≤M×N)。当计数值k小于或等于M×N的的数值时,执行步骤S502,反之(步骤S501中的“否”),将停止后抖色程序,完成处理量化影像xq。
在步骤S502中,取得量化误差值elocal的元素总和。当量化误差值elocal的元素总和大于负的预设临界值T的一半,或小于预设临界值T的一半(即,|elocal|<T/2)(步骤S502中的“否”)时,则执行步骤S503。当量化误差值elocal的元素总和介在负的预设临界值T的一半至预设临界值T的一半之间时(步骤S502中的“是”),则无法将量化误差值elocal的元素总和最小化,并进一步执行步骤S510。
在步骤S503中,判断量化误差值elocal的元素总和是否大于或等于预设临界值T的一半(即,elocal≥T/2)。若量化误差值elocal的元素总和大于或等于预设临界值T的一半时(在步骤S503中的“是”),在步骤S504中,找出在加总量化误差值elocal至抖色矩阵dlocal后数值最大的一矩阵元素(m,n)(即,(m,n)=max(xlocal+dlocal))。在步骤S505中,将此矩阵元素(m,n)设为在抖色矩阵dlocal中欲抖色的点。因此,进行处理此抖色矩阵dlocal的矩阵元素(m,n)以获得一负的预设临界值-T,例如,负的预设临界值-T为此抖色矩阵dlocal的矩阵元素(m,n)(即,dlocal(m,n)=-T)。在步骤S506中,令量化误差值elocal的元素总和等于量化误差值elocal的元素总和减去预设临界值T(即,elocal=elocal-T)。
在步骤S503中,若量化误差值elocal的元素总和小于或等于负的预设临界值-T的一半时(在步骤S503中的“否”),则执行步骤S507。在步骤S507中,找出在加总量化误差值elocal至抖色矩阵dlocal后数值最小的一矩阵元素(m,n)(即,(m,n)=min(xlocal+dlocal))。在步骤S508中,将此矩阵元素(m,n)设为在抖色矩阵dlocal中欲抖色的点。因此,进行处理此抖色矩阵dlocal的矩阵元素(m,n)以获得一预设临界值T,例如,预设临界值T为此抖色矩阵dlocal的矩阵元素(m,n)(即,dlocal(m,n)=T)。在步骤S509中,令量化误差值elocal的元素总和等于量化误差值elocal的元素总和加上预设临界值T(即,elocal=elocal+T)。在步骤S510中,令计数值k加1(即,k=k+1),且回到步骤S501直到计数值k大于M×N的数值为止。
III.实验结果及讨论
图6a、6b、6c分别显示每英吋300点(dots per inch,dpi)的一Lena影像、海滩影像及建筑影像。在实验中使用的这些原始影像具有256灰阶及500×500像素的解析度。在此实施例中使用16灰阶的均等量化器(UniformQuantizer)。4位的量化影像显示于图6d、6e、6f中。本发明所提出的抖色后程序是使用原始影像xlocal的一2×2局部区域。通过加上由本发明所提出抖色后程序所产生的抖色矩阵dlocal,可取得半色调影像y,如图6g、6h、6i所示。在图6a~6i中,很明显地,使用本发明所提出抖色后程序可以达到更高的影像品质。
在公式(1)中视觉误差的均方差用以估一此影像品质。在公式(1)中,人眼的调变转换函数所使用的v可表示为:
v(fr)=2.6(0.0192+0.114fr)exp{-(0.114fr)1.1},
其中fr是以角度为测量单位(Cycles Per Degree,CPD)的水平和垂直频率方均根(root-mean-square,rms)数值。针对不同大小的局部区域sx,视觉误差的实验结果如表1所示。
表1
与预期的一样,在表1中,可以发现本发明所提出的后抖色程序达到更低的视觉误差,且达到最佳的整体影像品质。值得注意的是,使用尺寸大小为局部区域的视觉误差比使用尺寸大小大于2×2局部区域的视觉误差效果较佳。因此,使用较小尺寸的局部区域可达到较佳的影像品质。
IV.结论
根据本发明的产生半色调影像的方法与影像处理系统,可执行本发明所提出的针对原始灰阶影像应用的抖色后程序以取得半色调影像。其中,为了进行影像品质的评估,在抖色后程序中,局部区域及抖色矩阵将人类的视觉特性考虑进去。因此,当人们观看所产生的半色调影像时,会感觉较为舒适。此外,本发明所提出的产生半色调影像的方法提供较高的影像品质、较低的位运算量及较低的计算复杂度。高效率的产生半色调影像的方法更可在各种电子设备中实施或运用,举例但不局限于,行动电话、无线装置、个人数字助理(Personal Digital Assistant,PDA)、手持式或可携式电脑以及电泳显示器(electrophorotic display)。
需注意的是,虽然本发明的前述模块及单元为此系统的单独组件,但此等组件可被整合至一起,因而降低系统内的组件数。同理,一或以上的组件可被分开使用,因而增加系统内的组件数。此外,本发明中的模块及单元组件可以使用任何硬体、软体方法来实现。
本发明的方法,或特定型态或其部份,可以以程序码的型态存在。程序码可以包含于实体媒体,如软盘、光盘、硬盘、或是任何其他电子设备或机器可读取(如电脑可读取)储存媒体,也或不限于外在形式的电脑程序产品,其中,当程序码被机器,如电脑载入且执行时,此机器变成用以参与本发明的装置或系统,且可执行本发明的方法步骤。程序码也可以通过一些传送媒体,如电线或电缆、光纤、或是任何传输型态进行传送,其中,当程序码被电子设备或机器,如电脑接收、载入且执行时,此机器变成用以参与本发明的系统或装置。当在一般用途处理单元实作时,程序码结合处理单元提供一操作类似于应用特定逻辑电路的独特装置。
虽然本发明已以较佳实施例揭露如上,然其并非用以限定本发明,任何熟悉此项技艺者,在不脱离本发明的精神和范围内,当可做些许更动与润饰,因此本发明的保护范围当视后附的申请专利范围所界定者为准。
Claims (18)
1.一种产生半色调影像的方法,用于一影像处理系统以产生一半色调影像,上述方法包括:
通过一量化器接收一具有原始像素值的原始影像,并根据一预设临界值产生一具有量化像素值的量化影像;
通过一减法器将上述原始像素值减去上述量化像素值,以产生一量化误差值;
通过一后抖色模块并依据上述量化误差值执行一后抖色程序,以产生一抖色值;以及
通过一加法器加总上述量化影像的上述量化像素值及上述抖色值,并产生一具有半色调像素值的半色调影像。
2.如权利要求1所述的产生半色调影像的方法,其中产生上述半色调影像的步骤还包括:
通过上述影像处理系统将上述原始影像区分为具有一局部区域尺寸的多个局部区域;以及
通过上述后抖色模块对每一局部区域执行上述后抖色程序。
3.如权利要求2所述的产生半色调影像的方法,其中上述局部区域系一包括一二维像素阵列(m,n)的矩阵,其中m及n皆为正整数且分别代表上述矩阵的第m行及第n列。
4.如权利要求3所述的产生半色调影像的方法,其中上述后抖色模块使用一具有上述局部区域尺寸且包括上述量化误差值的量化误差矩阵执行上述后抖色程序,以产生一具有上述局部区域尺寸且包括上述抖色值的抖色矩阵,其包括下列步骤:
(a)将上述抖色值初始为零;
(b)取得上述量化误差值的一总和;
(c)当上述量化误差值的上述总和大于负的上述预设临界值的一半,并小于上述预设临界值的一半时,停止上述后抖色程序;
(d)当上述量化误差值的上述总和大于或等于上述预设临界值的一半时,找出在加总上述量化误差矩阵至上述抖色矩阵的后一数值最大的矩阵元素(m,n),将上述抖色矩阵的矩阵元素(m,n)设为负的上述预设临界值,并令上述量化误差值的总和等于上述量化误差值的总和减去上述预设临界值;
(e)当上述量化误差值的上述总和小于或等于上述负的预设临界值的一半时,找出在加总上述量化误差矩阵至上述抖色矩阵的后一数值最小的矩阵元素(m,n),将上述抖色矩阵的矩阵元素(m,n)设为上述预设临界值,并令上述量化误差值的总和等于上述量化误差值的总和加上上述预设临界值;
(f)重复步骤(c)-(e),直到上述后抖色程序结束,并产生上述抖色矩阵。
5.如权利要求1所述的产生半色调影像的方法,其中上述原始影像系一灰阶影像。
6.如权利要求1所述的产生半色调影像的方法,其中上述方法用于一手持式电子装置。
7.一种影像处理系统,用以产生一半色调影像,上述系统包括:
一量化器,接收一具有原始像素值的原始影像,并根据一预设临界值产生一具有量化像素值的量化影像;
一减法器,耦接至上述量化器,用以将上述原始像素值减去上述量化像素值,以产生一量化误差值;
一后抖色模块,耦接至上述减法器,用以依据上述量化误差值执行一后抖色程序,以产生一抖色值;以及
一加法器,耦接至上述后抖色模块及上述量化器,用以加总上述量化影像的上述量化像素值及上述抖色值,并产生一具有半色调像素值的半色调影像。
8.如权利要求7所述的影像处理系统,其中上述影像处理系统将上述原始影像区分为具有一局部区域尺寸的多个局部区域(Local Area);以及上述后抖色模块对每一局部区域执行上述后抖色程序。
9.如权利要求8所述的影像处理系统,其中上述局部区域系一包括一二维像素阵列(m,n)的矩阵,其中m及n皆为正整数且分别代表上述矩阵的第m行及第n列。
10.如权利要求9所述的影像处理系统,其中上述后抖色模块使用一具有上述局部区域尺寸且包括上述量化误差值的量化误差矩阵执行上述后抖色程序,以产生一具有上述局部区域尺寸且包括上述抖色值的抖色矩阵,其包括下列步骤:
(a)将上述抖色值初始为零;
(b)取得上述量化误差值的一总和;
(c)当上述量化误差值的上述总和大于负的上述预设临界值的一半,并小于上述预设临界值的一半时,停止上述后抖色程序;
(d)当上述量化误差值的上述总和大于或等于上述预设临界值的一半时,找出在加总上述量化误差矩阵至上述抖色矩阵的后一数值最大的矩阵元素(m,n),将上述抖色矩阵的矩阵元素(m,n)设为负的上述预设临界值,并令上述量化误差值的总和等于上述量化误差值的总和减去上述预设临界值;
(e)当上述量化误差值的上述总和小于或等于上述负的预设临界值的一半时,找出在加总上述量化误差矩阵至上述抖色矩阵的后一数值最小的矩阵元素(m,n),将上述抖色矩阵的矩阵元素(m,n)设为上述预设临界值,并令上述量化误差值的总和等于上述量化误差值的总和加上上述预设临界值;
(f)重复步骤(c)-(e),直到上述后抖色程序结束,并产生上述抖色矩阵。
11.如权利要求7所述的影像处理系统,其中上述原始影像系一灰阶影像。
12.如权利要求7所述的影像处理系统,其中上述影像处理系统用于一手持式电子装置。
13.一种电脑产品,用以于一影像处理系统中执行产生半色调影像的方法,上述方法包括:
接收一具有原始像素值的原始影像,并根据一预设临界值产生一具有量化像素值的量化影像;
将上述原始像素值减去上述量化像素值,以产生一量化误差值;
依据上述量化误差值执行一后抖色程序,以产生一抖色值;以及
加总上述量化影像的上述量化像素值及上述抖色值,并产生一具有半色调像素值的半色调影像。
14.如权利要求13所述的电脑产品,其中产生上述半色调影像的步骤还包括:
将上述原始影像区分为具有一局部区域尺寸的多个区域;以及
对每一局部区域执行上述后抖色程序。
15.如权利要求14所述的电脑产品,其中上述局部区域一包括一二维像素阵列(m,n)的矩阵,其中m及n皆为正整数且分别代表上述矩阵的第m行及第n列。
16.如权利要求15所述的电脑产品,其中上述后抖色模块使用一具有上述局部区域尺寸且包括上述量化误差值的量化误差矩阵执行上述后抖色程序,以产生一具有上述局部区域尺寸且包括上述抖色值的抖色矩阵,其包括下列步骤:
(a)将上述抖色值初始为零;
(b)取得上述量化误差值的一总和;
(c)当上述量化误差值的上述总和大于负的上述预设临界值的一半,并小于上述预设临界值的一半时,停止上述后抖色程序;
(d)当上述量化误差值的上述总和大于或等于上述预设临界值的一半时,找出在加总上述量化误差矩阵至上述抖色矩阵的后一数值最大的矩阵元素(m,n),将上述抖色矩阵的矩阵元素(m,n)设为负的上述预设临界值,并令上述量化误差值的总和等于上述量化误差值的总和减去上述预设临界值;
(e)当上述量化误差值的上述总和小于或等于上述负的预设临界值的一半时,找出在加总上述量化误差矩阵至上述抖色矩阵的后一数值最小的矩阵元素(m,n),将上述抖色矩阵的矩阵元素(m,n)设为上述预设临界值,并令上述量化误差值的总和等于上述量化误差值的总和加上上述预设临界值;
(f)重复步骤(c)-(e),直到上述后抖色程序结束,并产生上述抖色矩阵。
17.如权利要求13所述的电脑产品,其中上述原始影像系一灰阶影像。
18.如权利要求13所述的电脑产品,其中上述方法用于一手持式电子装置。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/472,118 | 2012-05-15 | ||
US13/472,118 US8773456B2 (en) | 2012-05-15 | 2012-05-15 | Halftone image generation method and image processing system and computer program product thereof |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103428402A true CN103428402A (zh) | 2013-12-04 |
CN103428402B CN103428402B (zh) | 2015-09-09 |
Family
ID=49580957
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310175275.4A Expired - Fee Related CN103428402B (zh) | 2012-05-15 | 2013-05-13 | 产生半色调影像的方法与影像处理系统 |
Country Status (3)
Country | Link |
---|---|
US (1) | US8773456B2 (zh) |
CN (1) | CN103428402B (zh) |
TW (1) | TW201346838A (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6472336B2 (ja) * | 2014-06-18 | 2019-02-20 | キヤノン株式会社 | 画像処理装置、画像処理方法およびプログラム |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070133060A1 (en) * | 2005-12-14 | 2007-06-14 | Oce-Technologies B.V. | Method, apparatus and computer program for halftoning digital images |
CN101542361A (zh) * | 2007-06-15 | 2009-09-23 | 株式会社理光 | 减少电子纸显示器上的图像伪像的方法 |
-
2012
- 2012-05-15 US US13/472,118 patent/US8773456B2/en not_active Expired - Fee Related
-
2013
- 2013-05-07 TW TW102116148A patent/TW201346838A/zh unknown
- 2013-05-13 CN CN201310175275.4A patent/CN103428402B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070133060A1 (en) * | 2005-12-14 | 2007-06-14 | Oce-Technologies B.V. | Method, apparatus and computer program for halftoning digital images |
CN101542361A (zh) * | 2007-06-15 | 2009-09-23 | 株式会社理光 | 减少电子纸显示器上的图像伪像的方法 |
Also Published As
Publication number | Publication date |
---|---|
US8773456B2 (en) | 2014-07-08 |
TW201346838A (zh) | 2013-11-16 |
CN103428402B (zh) | 2015-09-09 |
US20130307867A1 (en) | 2013-11-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Orchard et al. | Color quantization of images | |
Yan et al. | Halftone visual cryptography with minimum auxiliary black pixels and uniform image quality | |
US11733648B2 (en) | Deep computational holography | |
US8897580B2 (en) | Error diffusion with color conversion and encoding | |
EP0650289A1 (en) | Method and apparatus for generating a halftone pattern for a multi-level output device | |
EP0595504A1 (en) | Model-based halftoning of color images | |
US20110142340A1 (en) | Halftone image generation method and image processing system and computer program product thereof | |
Lee et al. | A new visual cryptography with multi-level encoding | |
Guo et al. | Halftone-image security improving using overall minimal-error searching | |
JP2005286879A (ja) | 画像送信装置、画像受信装置および画像伝送システム | |
JP3154240B2 (ja) | 相関最小視覚変調パターンによるディジタル中間調調色 | |
Son | Inverse halftoning through structure-aware deep convolutional neural networks | |
Nerella et al. | Securing images using colour visual cryptography and wavelets | |
EP2187378B1 (en) | Image processing device, image processing method, and program | |
CN103428402B (zh) | 产生半色调影像的方法与影像处理系统 | |
Yu et al. | Stochastic screen halftoning for electronic imaging devices | |
JP2015104804A (ja) | 画像処理装置、画像処理方法及びプログラム | |
US7292728B1 (en) | Block quantization method for color halftoning | |
Kolpatzik et al. | Optimized error diffusion based on a human visual model | |
Su et al. | An FPGA implementation of chaotic and edge enhanced error diffusion | |
JP4293255B2 (ja) | 画像処理方法および画像処理装置 | |
Kwon et al. | Multiscale rank-based ordered dither algorithm for digital halftoning | |
Wong et al. | Feature-preserving multilevel halftoning algorithm | |
Asano et al. | Digital halftoning algorithms based on optimization criteria and their experimental evaluation | |
CN102131037A (zh) | 半色调图像产生方法及其图像处理系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150909 Termination date: 20190513 |
|
CF01 | Termination of patent right due to non-payment of annual fee |