CN1302045A - 图象处理方法和记录了图象处理程序的记录媒体 - Google Patents

图象处理方法和记录了图象处理程序的记录媒体 Download PDF

Info

Publication number
CN1302045A
CN1302045A CN00136129.5A CN00136129A CN1302045A CN 1302045 A CN1302045 A CN 1302045A CN 00136129 A CN00136129 A CN 00136129A CN 1302045 A CN1302045 A CN 1302045A
Authority
CN
China
Prior art keywords
parameter
view data
image processing
pixel
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
Application number
CN00136129.5A
Other languages
English (en)
Other versions
CN1167036C (zh
Inventor
英保茂
新田库治
出立祥一
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Noritsu Koki Co Ltd
Original Assignee
Noritsu Koki Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Noritsu Koki Co Ltd filed Critical Noritsu Koki Co Ltd
Publication of CN1302045A publication Critical patent/CN1302045A/zh
Application granted granted Critical
Publication of CN1167036C publication Critical patent/CN1167036C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/73Deblurring; Sharpening
    • G06T5/75Unsharp masking
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/20Image enhancement or restoration using local operators

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)
  • Facsimile Image Signal Circuits (AREA)
  • Color Image Communication Systems (AREA)

Abstract

提供在对数字图像数据进行锐化处理时不会发生强调噪音以及发生上冲和下冲现象等问题的图象处理方法和记录了图象处理程序的记录媒体。在利用拉普拉斯算子的锐化处理中,根据原图像分别设置控制取负的拉普拉斯算子的大小的参量k、控制噪音的影响的参量θ和控制上冲和下冲现象的参量λ,根据原图像的状态变更这些参量。

Description

图象处理方法和记录了图象处理程序的记录媒体
本发明涉及对数字图像数据进行锐化处理的图象处理方法。
迄今,作为将照片等的图像印制到印相纸上的照片印相装置,已提出了通过记录了照片图像的底片将印相纸曝光的模拟打印机以及根据用扫描仪等读取记录在底片上的图像的图像数据和由数码相机等拍摄的图像数据将印相纸曝光的数字打印机等各种方案。特别是数字打印机,通过与对图像数据进行色修正及浓度修正等的图象处理的图象处理装置组合使用,可以进行模拟打印机所不能实现的色修正和浓度修正等,同时,可以容易而迅速地得到符合顾客的要求的图像,现在已广泛地应用。
另外,由于袖珍相机及一次性相机等的普及,一般人拍摄彩色照片的机会多起来了。在最近的相机中,具备了自动调焦功能及自动闪光功能等,对摄影不熟悉的人也可以以某种程度的质量进行拍摄。但是,在自动调焦功能中,特别是是否合焦了,对拍摄者是难于掌握的,有时主要的被摄体没有位于画面的中央时等,在该主要的被摄体未合焦的状态下就进行了拍摄。另外,利用自动闪光功能在逆光拍摄时等可以拍摄出某种程度的照片,但是,即使如此,有时也会成为不自然的照片。此外,由于拍摄时的轻微的手抖动原因也会成为全体模糊的照片。
如上所述,在所拍摄的图像中,由于种种原因,有时发生模糊现象。因此,在上述图象处理装置中,在输入的数字图像数据中发生模糊时,有时进行称为锐化处理的图象处理。所谓锐化处理,就是将由于模糊而引起的图像的劣化复原,或根据要求为了使图像容易看而将图像的边缘进行强调的处理。作为锐化处理的具体的方法,有利用1次微分运算的强调和2次微分运算(拉普拉斯算子)的强调等,但是,根据可以进行比较自然的强调处理,广泛采用了拉普拉斯算子的强调。
利用拉普拉斯算子的处理是,根据原图像的图像数据的各像素值,通过使该图像数据的各像素的拉普拉斯算子取负而进行边缘的强调。设原图像的图像数据的各像素值即原信号为f(i,j)(i,j表示坐标)、利用拉普拉斯算子进行处理后的各像素值即处理过的信号为F(i,j)时,则上述处理可以表为下式。F(i,j)=f(i,j)-2f(i,j)    (1)
下面,根据图7(a)~(d)说明以上的处理。
图7(a)表示在横方向取各像素的位置、在纵方向取像素值的原信号的概略情况,图7(b)表示原信号的1次微分,图7(c)表示原信号的拉普拉斯算子,图7(d)表示用上述(1)式求的(原信号)-(拉普拉斯算子)。
图7(a)所示的原信号,左侧像素值低,右侧像素值高,在中央附近存在边界区域。在该边界区域的像素值的变化略平滑缓慢,成为边缘多少模糊的状态。
对该原信号进行1次微分时,如图7(b)所示,在边界区域的中央部分成为取极大值的分布。并且,通过进一步求微分而求拉普拉斯算子时,如图7(c)所示,成为在边界区域的左侧取极大值、而在右侧取极小值的分布。将该拉普拉斯算子从原信号取负时,在原信号中从左侧的像素值低的区域变化到边界区域的部分的值,通过拉普拉斯算子的左侧的极大值取负将进一步降低。另外,在原信号中,从边界区域变化到右侧的像素值高的区域的部分的值通过使拉普拉斯算子的右侧的极小值取负(由于极小值是负的值,所以,实际上成为正的)而进一步提高。因此,如图7(d)所示,边界区域的像素值的变化与原信号相比成为陡峭变化的状态,从而可以抑制模糊,强调边缘。
下面,根据数式说明图7(a)~(d)所示的处理。上述(1)式的右边可以变形为如下形式。即
【数式1】f(i,j)-2f(i,j)=9f(i,j)-(f(i+1,j)+f(i-1,j)+f(i,j+1)+f(i,j-1)
             +f(i+1,j+1)+f(i-1,j+1)+f(i+1,j-1)+f(i-1,j-1))    (2)
为了便于说明,上述(2)式成为设想以注目像素为中心的3×3的矩阵的数学式。在该(2)式中,位于坐标(i,j)处的像素位于像素值变化的斜率的中央附近即几乎没有斜率的变化的区域或几乎没有像素值变化的区域时,则为2f(i,j)=0,从F(i,j)=f(i,j)-2f(i,j)变为F(i,j)=f(i,j),几乎没有像素值的变化。
另一方面,位于坐标(i,j)处的像素位于像素值变化的斜率的下部即位于边缘的像素值的低的一侧时,在坐标(i,j)的周围,存在很多像素值比f(i,j)大的像素。因此,(2)式右边第2项以后的坐标(i,j)周围的8个像素的像素值的总和比f(i,j)的8倍大。因此,成为f(i,j)-2f(i,j)<f(i,j),像素值小了若干。
另外,位于坐标(i,j)处的像素位于像素值变化的斜率的上部即位于边缘的像素值高的一侧时,在坐标(i,j)的周围,存在很多像素值比f(i,j)小的像素。因此,(2)式右边第2项以后的坐标(i,j)周围的8个像素的像素值的总和比f(i,j)的8倍小。因此,成为f(i,j)-2f(i,j)>f(i,j),像素值大了若干。
这样,通过利用拉普拉斯算子运算进行图象处理,可以在图像的边缘进行使像素值低的一侧更低而像素值高的一侧更高的处理,所以,边缘就进行了强调处理。
如上所述,如果利用拉普拉斯算子运算进行图象处理,通过强调图像的边缘,可以抑制图像的模糊,但是,同时也发生了噪音被放大的问题以及发生上冲和下冲的问题。
下面,先说明噪音放大的问题。图8(a)表示在发生了噪音的图像中在横方向取各像素的位置、在纵方向取像素值的原信号,图8(b)表示用上述(1)式求的(原信号)-(拉普拉斯算子)。通过比较图8(a)和(b)可知,通过利用拉普拉斯算子进行图象处理,也将噪音放大了。下面,使用具体的数值进行说明。表1
原信号A  100  100  100  100  125  100  100  100  100
拉普拉斯算子  0  0  0  25 -50  25  0  0  0
 A-B  100  100  100  75  175  75  100  100  100
在表1中,第1行的各数字表示在一维排列的原信号的各像素值,表示在像素值为100的均匀的图像数据中在1各地方包含像素值成为125的噪音的状态。第2行的各数字表示与第1行的原信号对应的拉普拉斯算子,在没有像素值的变化的像素处为0,但是,在发生了噪音的像素的两边相邻的像素处为25、而在发生了噪音的像素处则成为-50。第3行的各数字表示(原信号)-(拉普拉斯算子),在没有像素值的变化的像素处和原信号相同,为100,但是,在发生了噪音的像素的两边相邻的像素处减少为75,而在发生了噪音的像素处则增大到175。
这样,在原信号中,像素值为125的噪音通过利用拉普拉斯算子进行图象处理,像素值就增大到175,显著地放大了噪音。另外,发生了噪音的像素的两边相邻的像素的像素值减少了若干,所以,成为噪音格外显著的状态。
下面,说明发生上冲和下冲的问题。图9(a)表示在边缘比较明确的图像中在横方向取各像素的位置而在纵方向取像素值的原信号,图9(b)表示原信号的拉普拉斯算子,图9(c)表示用上述(1)式求的(原信号)-(拉普拉斯算子)。如图9(c)所示,可知对边缘比较明确的图像利用拉普拉斯算子进行图象处理时,就以边缘为中心进行了使像素值低的一侧变为更低的值、而使像素值高的一侧变为更高的值的变换,结果,就发生了上冲和下冲。下面,使用具体的数值进行说明。表2
原信号C  100  100  100  150  150  150  100  100  100
拉普拉斯算D  0  0  50  -50  0  -50  50  0  0
 C-D  100  100  50  200  150  200  50  100  100
在表2中,第1行的各数字表示一维排列的原信号的各像素值,表示像素值为100的像素区域和像素值为150的像素区域相邻的状态即在边缘处像素值急剧变化的状态。第2行的各数字表示与第1行的原信号对应的拉普拉斯算子,在没有像素值的变化的像素处为0,但是,在与像素值为150的像素相邻的像素值为100的像素处,拉普拉斯算子成为50,在与像素值为100的像素相邻的像素值为150的像素处,拉普拉斯算子成为-50。第3行的各数字表示(原信号)-(拉普拉斯算子),在没有像素值的变化的像素处,成为与原图像像素数据相同的像素值,但是,在与像素值为150的像素相邻的像素值为100的像素处,像素值减少为50,在与像素值为100的像素相邻的像素值为150的像素处,像素值上升为200。
这样,可知对边缘比较明确的图像利用拉普拉斯算子进行图象处理时,在边缘附近将发生上冲和下冲。进行这样的处理时,例如,在拍摄的暗背景中的明亮的被摄体的照片图像中,被摄体的周边部分将明亮到超过所需要的程度,而被摄体的周围的背景将发生暗的边缘,从而将成为不自然的照片图像。
本发明就是为了解决上述问题而提案的,目的旨在提供对数字图像数据进行锐化处理时不会发生噪音的强调及上冲和下冲等问题的图象处理方法和记录图象处理程序的记录媒体。
为了解决上述问题,本发明1所述的图象处理方法是根据原图像的图像数据通过将把原图像的图像数据的分布作为函数所区域的各像素的2次微分取负而进行锐化处理的图象处理方法,其特征在于:通过变更根据原图像的图像数据决定取负的2次微分的大小的第1参量,来控制锐化处理的程度。
在上述方法中,第1参量决定根据原图像的图像数据取负的2次微分的大小,该第1参量可以进行变更。即,通过根据在原图像中发生的模糊的状态变更第1参量,可以进行用以成为所希望的锐化的状态的图象处理。
本发明2所述的图象处理方法的特征在于:在本发明1所述的图象处理方法中,从输入的图像数据中抽出表示该图像数据的特征的数据,通过将这些数据输入指定的算法,来设定上述第1参量。
按照上述方法,从输入的图像数据中抽出表示该图像数据的特征的数据,通过将这些数据输入指定的算法来设定第1参量,所以,在必须处理多个图像时等,可以减少处理时间和减轻操作人员的负担。
本发明3所述的图象处理方法是根据原图像的图像数据通过将把原图像的图像数据的分布作为函数所区域的各像素的2次微分取负而进行锐化处理的图象处理方法,其特征在于:上述2次微分可以由注目像素与位于其附近的多个像素的像素值的差分的总和而得到,同时,将上述各差分乘以系数,各系数根据对应的差分的大小而变化。
2次微分基本上可以由注目像素与位于其附近的多个像素的像素值的差分的总和而得到,但是,按照上述方法,将上述各差分乘以系数,各系数根据对应的差分的大小而变化。即,例如,如果差分的大小越大,将对应的系数的大小设定得也越大,则2次微分对差分的大小是小的附近像素的影响减小。因此,例如在附近像素发生微小的噪音时等,可以抑制该噪音的影响包含在2次微分中,所以,从而可以抑制通过锐化处理而噪音得到强调的现象。
本发明4所述的图象处理方法的特征在于:在本发明3所述的图象处理方法中,上述系数在上述差分的大小大于第2参量时和上述差分的大小小于第2参量时设定为不同的值,同时,通过变更上述第2参量而控制锐化的程度。
按照上述方法,上述系数在上述差分的大小大于第2参量时和上述差分的大小小于第2参量时设定为不同的值。这里,例如,如果上述系数设定为在上述差分的大小大于第2参量时成为大的值而在上述差分的大小小于第2参量时成为小的值,则和上述一样,例如在附近像素中发生微小的噪音时等可以抑制在2次微分中包含该噪音的影响。另外,通过根据原图像的图像状态变更第2参量,可以进行用以成为所希望的锐化的状态的图象处理。
本发明5所述的图象处理方法的特征在于:在本发明4所述的图象处理方法中,从输入的图像数据中抽出表示该图像数据的特征的数据,通过将这些数据输入指定的算法而设定所述第2参量。
按照上述方法,从输入的图像数据中抽出表示该图像数据的特征的数据,通过将这些数据输入指定的算法而设定上述第2参量,所以,在必须处理多个图像时等,可以减少处理时间和减轻操作人员的负担。
本发明6所述的图象处理方法是根据原图像的图像数据通过将把原图像的图像数据的分布作为函数所区域的各像素的2次微分取负而进行锐化处理的图象处理方法,其特征在于:对所述2次微分设置了上限和/或下限,通过变更决定该上限和/或下限的第3参量来控制锐化的程度。
按照上述方法,对上述2次微分设置了上限和/或下限,所以,例如对在原图像中充分强调了边缘的区域可以防止进而进行不必要的边缘强调处理。这样,便可抑制进行锐化处理的图像发生上冲和下冲的问题。另外,所述上限和/或下限可以利用第3参量进行变更,所以,可以根据原图像的图像状态变更上限的值,从而可以进行用以成为所希望的锐化的状态的图象处理。
本发明7所述的图象处理方法的特征在于:在本发明6所述的图象处理方法中,所述上限和/或下限由用色差微分值除所述第3参量的值决定。
按照上述方法,所述上限和/或下限由用色差分值除上述第3参量的值决定。所谓色差微分值,就是边缘强的区域,所谓色差微分值小的区域,就是边缘弱的区域。即,在边缘强的区域,上限和/或下限的绝对值减小,在边缘弱的区域,上限和/或下限的绝对值增大。即,对在原图像中充分强调了边缘的区域通过减小上限和/或下限的绝对值,可以防止2次微分不必要地增大,从而可以抑制超过所需要的锐化处理。另外,对在原图像中发生某种程度的模糊的区域通过增大上限和/或下限的绝对值,可以使2次微分比较大,从而可以可靠地进行锐化处理。
本发明8所述的图象处理方法的特征在于:在本发明6或7所述的方法中,从输入的图像数据中抽出表示该图像数据的特征的数据,通过将这些数据输入指定的算法来设定上述第3参量。
按照上述方法,从输入的图像数据中抽出表示该图像数据的特征的数据,通过将这些数据输入指定的算法来设定第3参量,所以,在必须处理多个图像时等,可以减少处理时间和减轻操作人员的负担。
本发明9所述的图象处理方法的特征在于:在本发明1~8的任一项发明所述的方法中,对原图像的像素数据利用2次微分对各色成分进行锐化处理。
按照上述方法,拉普拉斯算子运算对各色成分进行,所以,例如,在原图像中,在虽然是同等程度的浓度差但由于色相发生变化而形成边缘时,也可以将该区域识别为边缘,从而可以进行边缘的强调处理。
本发明10所述的记录图象处理程序的记录媒体是记录使计算机执行根据原图像的图像数据通过将把原图像的图像数据的分布作为函数所区域的各像素的2次微分取负而进行锐化处理的图象处理的图象处理程序的记录媒体,其特征在于:通过变更根据原图像的图像数据而取负的2次微分的大小的第1参量而使计算机执行控制锐化处理的程度的处理。
在记录到上述记录媒体上的程序中,第1参量决定根据原图像的图像数据而取负的2次微分的大小,该第1参量可以变更。即,通过根据原图像中发生的模糊的状态而变更第1参量,可以进行用以成为所希望的锐化的状态的图象处理。
本发明11所述的记录图象处理程序的记录媒体的特征在于:在本发明10所述的结构中,使计算机执行从输入的图像数据中抽出表示该图像数据的特征的数据并通过将这些数据输入指定的算法而决定上述第1参量的处理。
按照记录到上述记录媒体上的程序,从输入的图像数据中抽出表示该图像数据的特征的数据,通过将这些数据输入指定的算法而设定第1参量,所以,在必须处理多个图像时等,可以减少处理时间和减轻操作人员的负担。
本发明12所述的记录图象处理程序的记录媒体是记录使计算机执行根据原图像的图像数据通过将把原图像的图像数据的分布作为函数所区域的各像素的2次微分取负而进行锐化处理的图象处理的图象处理程序的记录媒体,其特征在于:上述2次微分由注目像素与其附近的多个像素的像素值的差分的总和得到,同时使计算机执行对上述各差分乘以系数并且各系数根据对应的差分的大小而变化的处理。
2次微分基本上可以由注目像素与其附近的多个像素的像素值的差分的总和得到,但是,按照记录到上述记录媒体上的程序,对上述各差分乘以系数,而各系数根据对应的差分的大小而变化。即,例如,如果设定为差分的大小越大,使对应的系数的大小越大,则2次微分将减小差分的大小是小的附近像素的影响。因此,例如在附近像素处发生微小的噪音时等,便可抑制该噪音的影响包含在2次微分中,所以,可以抑制由于锐化处理而噪音得到强调的情况。
本发明13所述的记录图象处理程序的记录媒体的特征在于:在本发明12所述的结构中,上述系数在上述差分的大小大于第2参量时和上述差分的大小小于第2参量时设定为不同的值,同时,通过变更上述第2参量使计算机执行控制锐化的程度的处理。
按照记录到上述记录媒体上的程序,上述系数在上述差分的大小大于第2参量时和上述差分的大小小于第2参量时设定为不同的值。这里,例如,如果在上述差分的大小大于第2参量时将上述系数设定为大的值而在上述差分的大小小于第2参量时设定为小的值,就和上述一样,例如在附近像素处发生微小的噪音时等,就可以抑制该噪音的影响包含在2次微分中。另外,通过根据原图像的图像状态变更第2参量,便可进行用以成为所希望的锐化的状态的图象处理。
本发明14所述的记录图象处理程序的记录媒体的特征在于:在本发明13所述的结构中,从输入的图像数据中抽出表示该图像数据的特征的数据,通过将这些数据输入指定的算法而使计算机执行设定上述第2参量的处理。
按照记录到上述记录媒体上的程序,从输入的图像数据中抽出表示该图像数据的特征的数据,通过将这些数据输入指定的算法而设定第2参量,所以,在必须处理多个图像时等,可以减少处理时间和减轻操作人员的负担。
本发明15所述的记录图象处理程序的记录媒体是记录使计算机执行根据原图像的图像数据通过将把原图像的图像数据的分布作为函数所区域的各像素的2次微分取负而进行锐化处理的图象处理的图象处理程序的记录媒体,其特征在于:对上述2次微分设置了上限和/或下限,通过变更决定该上限和/或下限的第3参量,使计算机执行控制锐化的程度的处理。
按照记录到上述记录媒体上的程序,对上述2次微分设置了上限和/或下限,所以,可以防止例如对在原图像中充分强调了边缘的区域进而进行不必要的边缘强调处理。这样,便可抑制在进行锐化处理的图像中发生上冲和下冲的现象。另外,上述上限和/或下限可以由第3参量进行变更,所以,可以根据原图像的图像的状态来变更上限和/或下限,从而可以进行用以成为所希望的锐化的状态的图象处理。
本发明16所述的记录图象处理程序的记录媒体的特征在于:在本发明15所述的结构中,使计算机执行通过用色差微分值除上述第3参量的值决定上述上限和/或下限的处理。
按照记录到上述记录媒体上的程序,上述上限和/或下限由用色差微分值除上述第3参量的值决定。所谓色差微分值大的区域,就是边缘强的区域,所谓色差微分值小的区域,就是边缘弱的区域。即在边缘强的区域中,上限和/或下限的绝对值减小,在边缘弱的区域中,上限和/或下限的绝对值增大。即,对在原图像中充分强调了边缘的区域通过减小上限和/或下限的绝对值,防止2次微分不必要地增大,从而可以抑制超过所需要的锐化处理。另外,对在原图像中发生了某种程度的模糊的区域通过增大上限和/或下限的绝对值,2次微分则变得比较大,从而可以可靠地进行锐化处理。
本发明17所述的记录图象处理程序的记录媒体的特征在于:在本发明15或16所述的结构中,从输入的图像数据中抽出表示该图像数据的特征的数据,通过将这些数据输入指定的算法而使计算机执行设定上述第3参量的处理。
按照记录到上述记录媒体上的程序,从输入的图像数据中抽出表示该图像数据的特征的数据,通过将这些数据输入指定的算法而设定上述第3参量,所以,在必须处理多个图像时等,可以减少处理时间和减轻操作人员的负担。
本发明18所述的记录图象处理程序的记录媒体的特征在于:在本发明10~17的任一项发明所述的结构中,使计算机对原图像的像素数据执行对各色成分通过2次微分进行锐化处理的处理。
按照记录到上述记录媒体上的程序,2次微分对各色成分进行,所以,例如虽然在原图像中是同等程度的浓度差但由于色相发生变化而形成了边缘时,也可以将该区域识别为边缘,从而可以进行边缘的强调处理。
下面,简要说明附图及符号:
图1是表示本实施例利用拉普拉斯算子运算进行的锐化处理的流程的流程图。
图2是表示本发明的实施例的图像输出系统的概略结构的框图。
图3是表示色差微分值的运算的流程的流程图。
图4是表示预先将参量λ和参量θ设定为适当的值并利用与理想直线的偏离判断法设定参量k的方法的流程图。
图5是表示对横轴为原图像的色差微分平均值、纵轴为各k的值进行锐化处理的图像的色差微分平均值的曲线图。
图6是表示横轴为k的值而纵轴为各k的值的比例常数ak的曲线图。
图7(a)表示在横方向取各像素的位置而在纵方向取像素值的原图像像素数据的概略的说明图,图7(b)表示原图像的像素数据的1次微分值的说明图,图7(c)表示原图像的像素数据的拉普拉斯算子值的说明图,图7(d)是表示由(1)式求的(原图像像素数据)-(拉普拉斯算子值)的说明图。
图8(a)是表示在发生噪音的图像中在横方向取各像素的位置而在纵方向取像素值的原图像像素数据的说明图,图8(b)是表示由(1)式求的(原图像像素数据)-(拉普拉斯算子值)的说明图。
图9(a)是表示在边缘比较明确的图像中在横方向取各像素的位置而在纵方向取像素值的原图像像素数据的说明图,图9(b)是表示原图像的像素数据的拉普拉斯算子值的说明图,图9(c)是表示由(1)式求的(原图像像素数据)-(拉普拉斯算子值)的说明图。
附图中,
1-胶片扫描仪;2-图象处理装置;3-印相装置;4-图象处理运算部;5-表色系变换部;6-色差微分运算部;7-参量决定部。
实施例:
下面,根据图1~图6说明本发明的一个实施例。
图2是表示本发明的实施例的图像输出系统的概略结构的框图。该图像输出系统具有胶片扫描仪1、图象处理装置2和印相装置3。
胶片扫描仪1通过例如使光源的光照射到作为照片胶卷的底片上而由CCD(Charge Coupled Device)等接收其透过光,读取记录在底片上的图像。该胶片扫描仪1将读取的图像数据按红色成分、绿色成分和蓝色成分向图象处理装置2输出。
印相装置3根据由图象处理装置2进行了处理的图像数据通过对作为感光材料的印相纸曝光而将图像印制到印相纸上。作为将与数字图像数据相应的光照射到印相纸上的光头,使用可以根据数字图像数据而按各像素调制向印相纸照射的照射光的光调制元件。作为光调制元件,有例如PLZT曝光光头、DMD(数字微型镜元件)、LCD(液晶显示装置)、LED(LightEmitting Diode)板、激光器、FOCRT(FiberOptic Cathode Ray Tube)和CRT(Cathode Ray Tube)等。
印相装置3,也可以采用能够进行底片的扫描和印相纸的曝光这两种操作的自动打印机构成。这时,通过将照片处理装置采用将进行从图像的读取到印制的自动打印机和由PC(Personal Computer)等构成的图象处理装置2连接的结构,便可将系统简化。
图象处理装置2对从胶片扫描仪1传送来的图像数据进行锐化处理,并将进行了锐化处理的图像数据供给照片印制装置3。该图象处理装置2具有图象处理运算部4、表色系变换部5、色差微分运算部6和参量决定部7。
图象处理运算部4是对从胶片扫描仪1传送来的图像数据利用后面所述的拉普拉斯算子进行锐化处理的块。表色系变换部5是对从胶片扫描仪1传送来的图像数据和由图象处理运算部4进行了处理的图像数据进行后面所述的表色系的变换处理的块。色差微分运算部6是在由表色系变换部5变换了表色系的图像数据中对各像素进行色差微分值的运算的块。参量决定部7是决定在图象处理运算部4的运算中所使用的参量的块。
下面,详细说明上述图象处理装置2的锐化处理的方法。在本实施例中,作为锐化处理,利用以下的方法进行拉普拉斯算子的处理。
设原图像的各像素的R(红)、G(绿)、B(蓝)的像素值分别为fR(i,j)、fG(i,j)、fB(i,j)(i,j表示坐标)、进行拉普拉斯算子处理后的各像素值分别为gR(i,j)、gG(i,j)、gB(i,j)时,则本实施例的拉普拉斯算子的处理可以用下式表示。gR(i,j)=fR(i,j)-k’2fR(i,j)    (3)gG(i,j)=fG(i,j)-k’2fG(i,j)    (4)gB(i,j)=fB(i,j)-k’2fB(i,j)    (5)
在上式中,由’2f(i,j)所示的运算表示在本发明中新定义的拉普拉斯算子。关于其详细情况,后面说明。
如上述(3)~(5)式那样,本实施例的拉普拉斯算子的处理的特征在于,分别对原信号的R、G、B的像素值进行。与此相反,例如迄今即使是以彩色信息得到了原信号时也是将各像素的像素值变换灰色调,而仅根据明度进行拉普拉斯算子的处理。这时,对于明度基本上相同而由于色相发生变化而形成了边缘的部分,则不能进行边缘的强调处理。另一方面,在本实施例的拉普拉斯算子的处理中,如上述那样,是对原图像数据的R、G、B的各像素值分别进行拉普拉斯算子的处理的,所以,可以对在各色成分中像素值发生变化的部分进行强调处理。因此,对明度的变化当然可以进行边缘的强调处理,对色相的变化也可以进行边缘的强调处理。
另外,在上述(3)~(5)式中,k是用于控制边缘强调的程度的参量。取k=1时,将和先有的拉普拉斯算子的运算一样的运算分别对R、G、B的各像素值进行。如果取k<1,就减弱边缘的强调的程度,如果取k>1,就加强边缘的强调的程度。这样,通过改变k的值,便可改变边缘的强调的程度。
在上述(3)~(5)式中,对拉普拉斯算子的系数全部为相同的k,但是,也可以在各式中分别使用不同的系数。但是,在各式中分别使用不同的系数时,对各色成分的拉普拉斯算子运算的修正的强度将发生变化,所以,作为修正的对象的像素的色相将有可能发生显著的变化。
下面,说明由上述’2fK(i,j)所示的拉普拉斯算子。上述’2fK(i,j)由下式定义。
【数式2】’2fK(i,j)=σi,jmin(|"2fK(i,j)|,λ/ΔC(i,j)}    (6)
Figure A0013612900181
在上述(6)式中,由"2fK(i,j)所示的运算表示在本发明中进而新定义的拉普拉斯算子。关于其详细情况,后面说明。另外,在上式中,ΔC(i,j)表示色差微分值。关于该色差微分值,后面说明。
如前所述,在原图像中对本来边缘部分就比较明确的区域进行拉普拉斯算子的处理时,将发生上冲和下冲现象,从而将成为不自然的图像。这是因为在浓度差大的区域拉普拉斯算子的绝对值也大,从而强调边缘而引起的浓度变化量增大的缘故。
为了抑制这样的上冲和下冲现象,可以考虑设置拉普拉斯算子的绝对值的上限,但是,将该上限设定为一定值时,图像中将会发生有模糊的区域和不模糊的区域或者有模糊的程度不同的区域等不理想的情况。这是由于以下所述的理由造成的。即,拉普拉斯算子的绝对值有边缘越强而越大的倾向,所以,在边缘越强的地方,越容易发生边缘。为了抑制这种现象,必须使拉普拉斯算子的绝对值的上限与图像中边缘最强的区域一致。于是,对模糊的程度大的区域将减弱锐化的程度。
因此,如上述(6)、(7)式那样,通过设定’2fK(i,j),可以在边缘强的地方将拉普拉斯算子的绝对值的上限设定小而在边缘弱的地方设定大。
在(6)式中,色差微分值ΔC(i,j)越大的区域,即像素值的变化量越大的区域拉普拉斯算子的绝对值的上限就越小,这相当于在边缘强的地方减小拉普拉斯算子的绝对值的上限。另外,色差微分值ΔC(i,j)越小的区域即像素值的变化量越小的区域拉普拉斯算子的绝对值的上限就越大,这相当于在边缘弱的地方增大拉普拉斯算子的绝对值的上限。
另外,在(6)式中,增大λ的值时,拉普拉斯算子的绝对值的上限就增大,从而上冲和下冲现象的抑制效果就减小。另外,减小λ的值时,拉普拉斯算子的绝对值的上限就减小,从而上冲和下冲现象的抑制效果就增大。即,λ就成为控制上冲和下冲现象的抑制效果的参量。
上述(6)、(7)式是在设定拉普拉斯算子的绝对值的上限同时可以在边缘强的地方将该上限设定小而在边缘弱的地方设定大的运算的一例,只要是具有同样的功能的运算,不论使用什么样的运算式都可以。另外,在(6)、(7)式的运算中,对拉普拉斯算子设置了上限和/或下限,是具有限幅器的效果的运算,但是,也可以例如对拉普拉斯算子的绝对值设置阈值从而设定具有超过该阈值的部分就采用以某一比率压缩的值的压缩器的效果的运算式。
下面,说明由在上述(6)式中使用的"2fK(i,j)所示的拉普拉斯算子。上述"2fK(i,j)由下式定义。
【数式3】"2fK(i,j)≡δi,j(1,0){f(i+1,j)-f(i,j)}
          +δi,j (-1,0){f(i-1,j)-f(i,j)}
          +δi,j (0,1){f(i,j+1)-f(i,j)}
          +δi,j (0,-1){f(i,j-1)-f(i,j)}
          +δi,j (1,1){f(i+1,j+1)-f(i,j)}
          +δi,j (-1,1){f(i-1,j+1)-f(i,j)}
          +δi,j (1,-1){f(i+1,j-1)-f(i,j)}
          +δi,j (-1,-1){f(i-1,j-1)-f(i,j)}    (8)其中,
Figure A0013612900201
如前所述,对在原图像中发生了噪音的区域进行拉普拉斯算子的处理时,将放大该噪音,从而将使画质劣化。为了抑制这种噪音的放大,如上述(8)、(9)式那样定义拉普拉斯算子,在与相邻的像素的浓度差小于常数θ时,在对该像素的拉普拉斯算子就不包含该相邻的像素的影响。即,将判定噪音对于周边的像素,其浓度差是比较小的。
在(8)、(9)式中,如果增大θ的值,抑制噪音的放大的效果将增大,但是,将减弱拉普拉斯算子的处理的强调效果。另外,如果减小θ的值,抑制噪音的放大的效果将减小,但是,将加强拉普拉斯算子的处理的强调效果。即,θ成为控制噪音的放大的抑制效果的参量。
上述(8)、(9)式是可以抑制噪音的放大的运算的一例,只要是具有同样的功能的运算,不论使用什么样的运算式都可以。例如,在(9)式中,δi,j (1,0)取1或0占某一个值,但是,也可以是根据|f(i+1,j)-f(i,j)|的大小而变化的函数值。
下面,详细说明上述色差微分值ΔC(i,j)。首先,说明L*a*b*表色系和色差。
在2个知觉色之间所知觉的色的差异称为色差。在知觉色空间,实现色差的均等性的空间称为均等色空间。作为均等色空间而规定的有CIE(国际照明委员会:Commission International de l’Eclairage)1964U*V*W*色空间、CIE1976 L*u*v*色空间和CIE1976L*a*b*色空间等,这里,使用L*a*b*色空间。从图像数据的RGB色空间向L*a*b*色空间的变换,利用以下的式子进行。
首先,使用以下的式子将RGB表色系变换为CIEXYZ表色系。
【数式4】
Figure A0013612900211
其次,变换为(Xc,Y,Zc)数据。
【数式5】
其中,(11)式的系数是为了对于R=G=B=1而将(10)式的X,Z规格化而导入的。对于该(Xc,Y,Zc)数据,可以利用下式求明度指数L和知觉色度指数(a,b)。
【数式6】
Figure A0013612900213
其中,V(x)=11.6x1/3-1.6(x>0.01)    (13)因此,有
【数式7】
其中,成为问题的是(14)式中的第1式,如果取Y=1,则有L=92.0。由于L*是要成为Munsell值的10倍而作成的,所以,在Y=1时,最好L=100.0。因此,如果对(14)式的各式乘以100.0/92.0=1.087,则可得到下式。但是,将符号变换为L→L*、a→a*、b→b*,同时将系数取整。
【数式8】
Figure A0013612900221
利用该(15)式定义L*a*b*色空间。并且,可以利用下式求代表L*a*b*色空间的2个知觉色的2点间的直线距离即色差ΔE。
ΔE=[(ΔL*)2+(Δa*)2+(Δb*)2]1/2    (16)
并且,对各像素点(i,j),将其与所有的相邻点的色差ΔE中最大的色差作为该像素点的色差微分值ΔC(i,j)。这样定义时,在彩色图像中的边界线区域,ΔC(i,j)的值增大,而在其他的区域,ΔC(i,j)的值减小。
这里,根据上述定义试对实际的彩色图像数据进行色差微分时,可知在图像中的暗的区域中,在不是边界线区域的区域将出现大的色差微分值。由于在暗的区域即(Xc,Y,Zc)的数据值小的区域,由上述(13)式给定的V(x)在x的值小的范围内具有大的斜率,所以,可以认为这是由于微小的(Xc,Y,Zc)的数据值之差采用(L*,a*,b*)的数据之差而增大所引起的。
因此,在本实施例中,将L*a*b*表色系进行以下的修正。首先,对于V(x),为了使之在x的值小的范围内斜率不用增大,如下式所示的那样,重新定义,以使斜率的大小保持一定。V(x)=10x    (17)这样定义V(x)时,L*、a*、b*分别可以用下式表示。
【数式9】
Figure A0013612900231
如上述那样修正L*a*b*表色系后,和上述一样对实际的彩色图像数据试进行色差微分时,在图像中的暗的区域中,在不是边界线区域的区域就不会再出现大的色差微分值。但是,在图像中的明亮的区域中,可知在不是边界线区域的区域中将出现大的色差微分值。这是由于在x的值小的范围内V(x)的斜率比修正前小而在x的值大的范围内V(x)的斜率比修正前大的缘故。
因此,为了灵活地应用修正前和修正后的L*a*b*表色系的特性,进而进行以下的修正。即,在x的值大的范围内使用修正前的V(x),在x的值小的范围内使用修正后的V(x)。分界处的x的值采用在修正前的V(x)=11.6x1/3-1.6((13)式)中斜率的大小为10的点即x=0.2404。如果用公式表示,则为
【数式10】
在本实施例中,L*a*b*表色系仅用于色差的计算,色差仅由(ΔL*,Δa*,Δb*)的数据值决定,所以,(L*,a*,b*)的数据值本身并不重要。因此,仅ΔV(x)即V(x)的斜率的大小才是重要的,与V(x)本身的值无关。因此,在上述(19)式中,V(x)在x=0.2404连续,而在x≥0.2404的范围内,则给V(x)加上-3.209的偏置。
并且,通过对由(19)式求的V(x)应用下式,进一步修正L*a*b*表色系。
【数式11】
在由该(20)式定义的修正L*a*b*色空间,根据上述(16)式计算色差ΔE,对各像素点(i,j),将其与所有的相邻点的色差ΔE中最大的色差作为该像素点的色差微分值ΔC(i,j)。并且,使用这样求出的色差微分值ΔC(i,j)进行上述(6)式的运算。
下面,根据图3的流程图说明上述色差微分值的运算的流程。
首先,从输入的图像数据中设定注目像素(步骤1,以后称为S1)。其次,将由RGB表色系输入的注目像素的数据变换为XYZ表色系的数据(S2)。并且,根据由(19)式分为不同情况的式子,利用(20)式将XYZ表色系变换为修正的L*a*b*表色系(S3)。
其次,根据变换为修正的L*a*b*表色系的图像数据计算注目像素与其附近的8像素的色差(S4)。并且,将这些色差中取最大的值的色差作为色差微分值进行存储(S5)。
然后,在S6,确认对全部像素的色差微分值的计算是否结束,在未结束时就返回到S1,反复进行上述处理。
使用上述那样求出的色差微分值,也可以抽出输入的图像中的边界线。例如,在所有的像素的色差微分值中,从其值大的一方开始,预先将与15%相当的像素取为1,而将其他的像素取为0。并且,画出在图像数据中为0的像素和1的像素的分布状态,通过用线连结由1的像素所占据的区域的中心附近,便可抽出边界线。另外,以上举出了将色差微分值进行取为0和1的二进制处理而抽出边界线的例子,但是,也可以通过采用二进制以即n值化(n=3,4,…)而抽出边界线。
下面,根据图1的流程图说明本实施例利用拉普拉斯算子运算反映锐化处理的流程。
首先,在原图像的图像数据中设定注目像素(S11)。并且,对各色成分求注目像素与注目像素附近的8像素的像素值的差分(S12)。根据该差分利用上述(9)式决定上述(8)式的各项的系数,并求"2fK(i,j)(S13)。
其次,用上述(7)式判断"2fK(i,j)的正负,同时比较"2fK(i,j)的绝对值与λ/ΔC(i,j)的大小,利用上述(6)式求"2fK(i,j)(S14)。并且,根据该"2fK(i,j),利用上述(3)~(5)式生成锐化处理图像数据(S15)。
然后,在S16,确认对全部像素的锐化处理是否结束,在未结束时就返回到S11,反复进行上述处理。
如上所述,本实施例的利用拉普拉斯算子运算的锐化处理,通过分别改变在(3)~(5)式中使用的用于控制边缘强调的程度的参量k、在(6)和(7)式中使用的控制上冲和下冲现象的抑制效果的参量λ和在(8)和(9)式中使用的控制噪音的放大的抑制效果的参量θ可以改变锐化处理的程度。因此,可以根据图像的用途和目的改变锐化的作用程度。
例如,在通常作为快摄照片而拍摄的照片图像中,在应合焦的区域发生模糊时,在先有的锐化处理中,则发生噪音被放大、产生上冲和下冲现象以及对不需要合焦的区域进行了锐化处理等问题。与此相反,如果是本实施例的锐化处理,通过适当地改变上述3个参量,对已合焦的区域和不需要合焦的区域可以减弱锐化的程度,而对发生了模糊并且希望合焦的区域可以进行适度的锐化处理。
作为设定上述3个参量的数值的方法,可以考虑将进行锐化处理的图像显示在显示器等上并由操作人员通过手动改变各参量,确认锐化的程度来决定各参量的方法。该方法是由操作人员确认锐化的程度来决定参量的,所以,可以将锐化的程度设定为操作人员认为最佳的状态。但是,在必须处理多个图像时,将发生处理时间增加和对操作人员的负担增大的问题。
因此,考虑了自动地设定上述3各参量的方法。例如,可以考虑从输入的图像数据中抽出几种表示该图像数据的特征的数据并通过将这些数据输入指定的算法而设定上述3个参量的方法。作为上述算法,可以应用例如神经网络或遗传算法等以及其他最合适的探索算法等。
在实际进行本实施例的锐化处理时,在上述3参量中,参量λ和参量θ与参量k相比,对于锐化处理适当的值的范围比较大,从而伴随该变化的锐化处理的程度的变化就比较小。因此,使用预先将参量λ和参量θ设定为适当的值而根据输入的图像数据自动地设定参量k的方法,也可以充分地进行锐化处理的状态的控制。下面,说明仅自动地设定该参量k的参量设定方法。
这里,先说明参量k的值与上述色差微分值的平均值(色差微分平均值)的关系。
对共计40张采样图像,将参量k的值分别改变为k=0.200、0.400、0.600,研究色差微分平均值的变化。图5是表示横轴为原图像的色差微分平均值、纵轴为对各k的值进行了锐化处理的图像的色差微分平均值的曲线图。这样,便可知道色差微分平均值的变化对于各k的值基本上排列在通过原点的一条直线上。即,对于各k的值,处理后的图像的色差微分平均值与原图像的色差微分平均值成正比地变化。下面,试求该比例常数。
对40张采样图像分配图像号码1~40,对于图像号码i(i=1、2、…、40)的图像,设原图像的色差微分平均值为ΔCoim、对参量k进行了锐化处理的图像的色差微分平均值为ΔCk im。
用参量k进行了锐化处理的图像的色差微分平均值按比例常数ak与原图像的色差微分平均值成正比地变化时,其误差εi和误差εi的平方和S可以用下式表示。
【数式12】εi=ΔCk im-(ak·ΔCo im)    (21)S=∑εi 2=Σ(ΔCk im-(ak·ΔCo im))2 (22)
为了求比例常数ak,可以求使上述误差εi的平方和S为最小的ak。使S为最小的ak,可以如以下所示的(23)式那样,通过求解令S的1次微分为0的方程式而如(24)式那样求出。
【数式13】dS/dak=-2∑ΔCo im(ΔCk im-(ak·ΔCo im))=0    (23)ak=∑(ΔCo im·ΔCk im)/∑(ΔCo im)2           (24)
根据(24)式,如果求k=0.200、0.400、0.600时的比例常数,则得
【数式14】
Figure A0013612900271
根据上述结果可知,比例常数ak随着k的值增大而增加,但是,其增加的比率减小。即,可以认为,在增大k的值时,色差微分平均值也增加,但是,其增加的比率在减小。
图6表示横轴为k的值、纵轴为各k的值的比例常数ak的曲线图。如图6所示,在k的值比较小的范围内,比例常数ak与k的值的增大成正比地增加,但是,随着k的值的增大,比例常数ak的增加的比率减小。即,比例常数ak的增加存在饱和的k的值。
这里,将在k的值比较小的范围内根据比例关系求出的直线定义为理想直线。这时,如果k的值增大,实际的比例常数ak将偏离理想直线。下面,说明着眼于与该理想直线的偏离来决定成为饱和状态的k的值的与理想直线的偏离判断法。
可以认为,原图像的色差微分平均值越大,将饱和状态的判断所使用的与理想直线的偏离的值设定得越大的方法可以得到更自然的锐化结果。这里,为了进行规格化,利用下式定义相对色差微分平均值(ΔCknm)rel。(ΔCknm)x-1=ΔCknm/ΔCko(n=0,1,2,…)    (26)
另外,假定理想直线通过与原图像对应的点(k0=0,(ΔCkom)rel=1),并使用2点(k1,(ΔCk1m)rel)、(k2,(ΔCk2m)rel)表示其斜率a时,则得
【数式15】 a = k 1 · { ( Δ C k 1 m ) rel - 1 } + k 2 · { ( Δ C k 1 m ) rel - 1 } k 1 2 + k 2 2 ( 27 ) 并且,如下式那样定义与理想直线的偏离δkn。δkn=i(a·kn+1)-(ΔCknm)rel(n=1,2,…)  (28)
这里,在满足δkn≥d时,与理想直线的偏离就大,即判定处于了饱和的状态,设满足该关系的最小的kn为kemph。将该kemph设定为在上述锐化处理中使用的参量k。
以上所述,是利用与理想直线的偏离判断法来决定kemph的,但是,并不仅限于此,也可以应用例如斜率变化率判断法等。
所谓斜率变化率判断法,就是根据画出色差微分平均值的曲线的斜率,着眼于其变化率来判断是否为饱和状态的方法。设色差微分平均值的参量k的间隔为0.025时,则与第n(n=0、1、2、…)个对应的参量k的值kn为kn=0.025×n(n=0、1、2、…)    (29)
设用k=kn进行锐化处理的图像的色差微分平均值为ΔCknm,就用下式定义所画的曲线的斜率的大小。在下式中,取n=1、2、…。Δkn=(ΔCknm-ΔCkn-1m)/(kn-k(n-1))    (30)使用(30)式,将斜率变化率ρkn用下式定义。ρkn=(Δk1kn)/Δkn(n=1、2、…)    (31)
这里,在满足|ρkn|≥p时,就判定偏离直线并处于饱和状态,将满足该式的最小的kn作为kemph
在与理想直线的偏离判断法中,根据实验结果,确认了在作为是饱和状态的δkn≥d,δkn单调地增加。另一方面,在斜率变化率判断法中,根据实验结果,确认了在作为是饱和状态的|ρkn|≥p的范围内,|ρkn|不一定单调地增加。即,用比视为是饱和状态并且是锐化的状态大的k的值,却得到看来是不锐化的状态。例如,由于斜率变化率判断法所描绘的点的随机性高,所以,在改变所画的参量k的间隔时,有可能将完全不同的值选择为kemph了。根据以上所述,可以认为,在设定适合于自然的锐化的参量k的值kemph时,与理想直线的偏离判断法的可靠性高。
下面,说明为了利用与理想直线的偏离判断法决定kemph而按一定的步长改变k的值并通过计算其偏离量来探索成为饱和状态的k的值的最小值的方法。
这里,例如设想在k=0.500进入饱和状态的情况。这时,在上述的探索中,将改变k的值时的步长取为例如0.025时,要判定「适合于自然的锐化的参量kemph为0.500」,则必须求21次色差微分平均值,要进行很多无谓的运算。
因此,利用在与理想直线的偏离δkn在δkn≥d的范围内基本上可靠地单调增加和即使参量k的值有约0.025的差异而在处理后的图像中也不会发生肉眼可以判别的差别的情况,而将上述步长设定为0.050。但是,为了理想直线采用和步长为0.025时相同的直线,在进行k=0.025的运算后,在k≥0.050的范围内,将步长取为0.050。通过如上述那样设定步长,可以将求kemph所需要的运算次数减半。
但是,在以的方法中,求kemph所需要的运算的次数也要12次。因此,为了进一步减少求kemph所需要的运算的次数的平均值,变更按照以下的算法进行运算的顺序。但是,必须求在以下的算法中所使用的跳跃数s,下面,说明设定s=2时的算法。
首先,对k=0.000、0.025、0.050求相对色差微分平均值,按照上述(27)式求理想直线。并且,使所画的参量k的值仅跳过跳跃数s,对k=0.200、0.350、0.500、…求相对色差微分平均值,同时,求各个值的与理想直线的偏离,并将开始时满足δkn≥d的k的值作为k(2)。并且,对k(2)-0.100、k(2)-0.050求相对色差微分平均值,同时求各个值的与理想直线的偏离。在这些值中满足δkn≥d时,就将这时的k的值作为kemph,在没有满足δkn≥d的值时,就将k(2)作为kemph
按照以上的算法,将和以步长0.050顺序描绘时相同的值选择为kemph,便可减少平均的描绘次数。例如,在kemph=0.500时,求直至到该值为止的相对色差微分平均值的次数就下降到8次。
下面,作为锐化处理的一个实施例,参照图4所示的流程图说明如上述那样预先将参量λ和参量θ设定为适当的值而利用与理想直线的偏离判断法设定参量k的方法。
首先,在S21,设定参量λ和参量θ的值。这些设定值,采用根据实验结果而适当设定的固定值。
其次,在S22,根据作为进行锐化处理的图像而输入的原信号,对各像素求其色差微分值ΔC(i,j)。并且,根据对所有的像素得到的色差微分值ΔC(i,j)求该原信号的色差微分平均值ΔCom(S23)。所谓原信号的色差微分平均值,就是和将参量k取为0时的锐化处理图像的色差微分平均值同等的平均值。
然后,在S24,在k=0.025的状态,作成锐化处理图像,对各像素求其色差微分值ΔC(i,j),并求锐化处理图像的色差微分平均值ΔCom。同样,在S25,在k=0.050的状态,作成锐化处理图像,对各像素求其色差微分值ΔC(i,j),并求锐化处理图像的色差微分平均值ΔCom。
并且,根据在S23求出的原信号的、在S24求出的k=0.025的色差微分平均值和在S25求出的k=0.050的色差微分平均值计算理想直线的斜率(S26)。
然后,从k=0.200开始,以0.150的间隔改变k的值,同时对各个k的值计算与理想直线的偏离δkn(S27)。并且,将与理想直线的偏离δkn满足δkn≥d的最初的k的值设定为k(2)(S28)  。
然后,对k(2)-0.100、k(2)-0.050求相对色差微分平均值,同时求各个值的与理想直线的偏离(S29)。在这些值中满足δkn≥d时(在S210为YES(Y)),就将这时的k的值作为kemph(S211)。另一方面,在没有满足δkn≥d的值时(在S210为NO),就将k(2)作为kemph(S212)。
并且,将上述kemph设定为k,同时使用在S21设定的参量λ和参量θ的值对输入的原图像进行锐化处理(S213),然后,结束处理。
以上说明的图像输出方法,可以由图像的锐化处理功能的程序来实现。该程序存储在可以由计算机读取的记录媒体上。这时,图2所示的图象处理装置2由计算机构成。并且,作为上述记录媒体,可以是由构成图象处理装置2的计算机所具有的微处理器进行处理所需要的存储器,例如ROM本身可以是程序媒体,另外,作为外部存储装置,可以是设置程序读取装置并通过将记录媒体插入其中而可以读取的程序媒体。
不论是哪种情况,存储的程序都可以由上述计算机所具有的微处理器访问并执行,或者不论哪种情况都可以读出程序,并且读出的程序下载到微处理器的程序存储区域,从而是执行该程序的方式。假定该下载用的程序预先存储在本体装置中。
这里,上述程序媒体是构成为可以与本体分离的记录媒体,可以是磁带或盒式磁带等带系列、软盘或硬盘等磁盘或CD-ROM/MO/MD/DVD等光盘的盘系列、IC卡(包括存储卡)/光卡等卡系列,或者包含掩膜ROM、EPROM、EEPROM、闪存(ROM)等半导体存储器的固定的保持程序的媒体。
另外,在本发明中,可以采用可以与因特网的通信网络连接的系统结构,为了从通信网络中下载程序,可以是流动地保持程序的媒体。这样,在从通信网络中下载程序时,该下载用的程序可以预先存储在本体装置中,或者从别的记录媒体中获得。
作为存储在记录媒体中的内容,不限定程序,也可以是数据。
如上所述,本发明1所述的图象处理方法,是根据原图像的图像数据通过将把原图像的图像数据的分布作为函数所区域的各像素的2次微分取负而进行锐化处理的图象处理方法,是通过变更根据原图像的图像数据决定取负的2次微分的大小的第1参量来控制锐化处理的程度的方法。
这样,通过根据在原图像中发生的模糊的状态变更第1参量,便可进行用以成为所希望的锐化的状态的图象处理。
本发明2所述的图象处理方法,是从输入的图像数据中抽出表示该图像数据的特征的数据,通过将这些数据输入指定的算法而决定上述第1参量的方法。
这样,除了本发明1所述的方法的效果外,在必须处理多个图像时等,可以减少处理时间和减轻操作人员的负担。
本发明3所述的图象处理方法,是根据原图像的图像数据通过将把原图像的图像数据的分布作为函数所区域的各像素的2次微分取负而进行锐化处理的图象处理方法,是上述2次微分可以由注目像素与其附近的多个像素的像素值的差分的总和而得到,同时对上述各差分乘以系数并且各系数根据对应的差分的大小而变化的方法。
这样,例如,如果设定为差分的大小越大而对应的系数的大小也越大,则2次微分可以减小差分的大小小的附近像素的影响。因此,例如在附近像素中发生微小的噪音时等,便可抑制该噪音的影响包含在2次微分中,所以,可以抑制通过锐化处理而强调了噪音。
本发明4所述的图象处理方法,是上述系数在上述差分的大小大于第2参量时和上述差分的大小小于第2参量时设定为不同的值、同时通过变更上述第2参量来控制锐化的程度的方法。
这样,除了本发明3所述的方法的效果外,例如,如果将上述系数在上述差分的大小大于第2参量时设定为大的值而在上述差分的大小小于第2参量时设定为小的值,就和上述一样,例如在附近像素中发生微小的噪音时等,可以抑制在2次微分中包含该噪音的影响。另外,通过根据原图像的图像状态变更第2参量,可以进行用以成为所希望的锐化的状态的图象处理。
本发明5所述的图象处理方法,是从输入的图像数据中抽出表示该图像数据的特征的数据并通过将这些数据输入指定的算法而设定上述第2参量的方法。
这样,除了本发明4所述的方法的效果外,在必须处理多个图像时等,可以减少处理时间和减轻操作人员的负担。
本发明6所述的图象处理方法,是根据原图像的图像数据通过将把原图像的图像数据的分布作为函数所区域的各像素的2次微分取负而进行锐化处理的图象处理方法,是对上述2次微分设置上限和/或下限并通过变更决定该上限和/或下限的第3参量来控制锐化的程度的方法。
这样,便可抑制在进行了锐化处理的图像中发生上冲和下冲现象的问题。另外,根据原图像的图像状态变更上限和/或下限,可以进行用以成为所希望的锐化的状态的图象处理。
本发明7所述的图象处理方法,是由用色差微分值除上述第3参量的值决定上述上限和/或下限的方法。
这样,除了本发明6所述的方法的效果外,对在原图像中已充分强调了边缘的区域,通过减小上限和/或下限的绝对值,便可防止2次微分大到不必要的程度,从而可以抑制超过所需要的锐化处理。另外,对在原图像中发生了某种程度的模糊的区域,通过增大上限和/或下限的绝对值,2次微分便可比较大,从而可以可靠地进行锐化处理。
本发明8所述的图象处理方法,是从输入的图像数据中抽出表示该图像数据的特征的数据并通过将这些数据输入指定的算法而设定上述第3参量的方法。
这样,除了本发明6或7所述的方法的效果外,在必须处理多个图像时等,可以减少处理时间和减轻操作人员的负担。
本发明9所述的图象处理方法,是对原图像的像素数据按各色成分利用2次微分进行锐化处理的方法。
这样,除了本发明1~8的任一项发明所述的方法的效果外,在例如在原图像中虽然是同等程度的浓度差但由于色相发生变化而形成边缘时,也可以将该区域识别为边缘从而可以进行边缘的强调处理。
本发明10所述的记录图象处理程序的记录媒体,是记录使计算机执行根据原图像的图像数据通过将把原图像的图像数据的分布作为函数所区域的各像素的2次微分取负而进行锐化处理的图象处理的图象处理程序的记录媒体,是通过变更根据原图像的图像数据而取负的2次微分的大小的第1参量而使计算机执行控制锐化处理的程度的处理的结构。
这样,通过根据在原图像中发生的模糊的状态变更第1参量,便可进行用以成为所希望的锐化的状态的图象处理。
本发明11所述的记录图象处理程序的记录媒体,是使计算机执行从输入的图像数据中抽出表示该图像数据的特征的数据并通过将这些数据输入指定的算法而决定上述第1参量的处理的结构。
这样,除了本发明10所述的结构的效果外,在必须处理多个图像时等,可以减少处理时间和减轻操作人员的负担。
本发明12所述的记录图象处理程序的记录媒体,是记录使计算机执行根据原图像的图像数据通过将把原图像的图像数据的分布作为函数所区域的各像素的2次微分取负而进行锐化处理的图象处理的图象处理程序的记录媒体,是上述2次微分由注目像素与其附近的多个像素的像素值的差分的总和得到,同时使计算机执行对上述各差分乘以系数并且各系数根据对应的差分的大小而变化的处理的结构。
这样,例如,如果充分的大小越大将对应的系数的大小设定得越大,2次微分就可以减小充分的大小小的附近像素的影响。因此,例如在附近像素处发生微小的噪音时等,便可抑制该噪音的影响包含在2次微分中,所以,可以抑制由于锐化处理而噪音得到强调的情况。
本发明13所述的记录图象处理程序的记录媒体,是上述系数在上述差分的大小大于第2参量时和上述差分的大小小于第2参量时设定为不同的值,同时,通过变更上述第2参量使计算机执行控制锐化的程度的处理的结构。
这样,除了本发明12所述的结构的效果外,例如,如果在上述差分的大小大于第2参量时将上述系数设定为大的值而在上述差分的大小小于第2参量时设定为小的值,就和上述一样,例如在附近像素处发生微小的噪音时等,就可以抑制该噪音的影响包含在2次微分中。另外,通过根据原图像的图像状态变更第2参量,便可进行用以成为所希望的锐化的状态的图象处理。
本发明14所述的记录图象处理程序的记录媒体,是从输入的图像数据中抽出表示该图像数据的特征的数据,通过将这些数据输入指定的算法而使计算机执行设定上述第2参量的处理的结构。
这样,除了本发明13所述的结构的效果外,在必须处理多个图像时等,可以减少处理时间和减轻操作人员的负担。
本发明15所述的记录图象处理程序的记录媒体,是记录使计算机执行根据原图像的图像数据通过将把原图像的图像数据的分布作为函数所区域的各像素的2次微分取负而进行锐化处理的图象处理的图象处理程序的记录媒体,是对上述2次微分设置了上限和/或下限,通过变更决定该上限和/或下限的第3参量,使计算机执行控制锐化的程度的处理的结构。
这样,便可抑制在进行了锐化处理的图像中发生上冲和下冲现象的问题。另外,根据原图像的图像状态变更上限和/或下限的值,可以进行用以成为所希望的锐化的状态的图象处理。
本发明16所述的记录图象处理程序的记录媒体,是使计算机执行通过用色差微分值除上述第3参量的值决定上述上限和/或下限的处理的结构。
这样,除了本发明15所述的结构的效果外,对在原图像中充分强调了边缘的区域通过减小上限和/或下限的绝对值,防止2次微分不必要地增大,从而可以抑制超过所需要的锐化处理。另外,对在原图像中发生了某种程度的模糊的区域通过增大上限和/或下限的绝对值,2次微分则变得比较大,从而可以可靠地进行锐化处理。
本发明17所述的记录图象处理程序的记录媒体,是从输入的图像数据中抽出表示该图像数据的特征的数据,通过将这些数据输入指定的算法而使计算机执行设定上述第3参量的处理的结构。
这样,除了本发明15或16所述的结构的效果外,在必须处理多个图像时等,可以减少处理时间和减轻操作人员的负担。
本发明18所述的记录图象处理程序的记录媒体,是使计算机对原图像的像素数据执行对各色成分通过2次微分进行锐化处理的处理的结构。
这样,除了本发明10~17的任一项发明所述的结构的效果外,例如虽然在原图像中是同等程度的浓度差但由于色相发生变化而形成了边缘时,也可以将该区域识别为边缘,从而可以进行边缘的强调处理。

Claims (18)

1.一种图象处理方法,是根据原图像的图像数据通过将把原图像的图像数据的分布作为函数所区域的各像素的2次微分取负而进行锐化处理的图象处理方法,其特征在于:通过变更根据原图像的图像数据决定取负的2次微分的大小的第1参量,来控制锐化处理的程度。
2.根据权利要求1所述的图象处理方法,其特征在于:从输入的图像数据中抽出表示该图像数据的特征的数据,通过将这些数据输入指定的算法,来设定上述第1参量。
3.一种图象处理方法,是根据原图像的图像数据通过将把原图像的图像数据的分布作为函数所区域的各像素的2次微分取负而进行锐化处理的图象处理方法,其特征在于:所述2次微分可以由注目像素与位于其附近的多个像素的像素值的差分的总和而得到,同时,将所述各差分乘以系数,各系数根据对应的差分的大小而变化。
4.根据权利要求3所述的图象处理方法,其特征在于:所述系数在所述差分的大小大于第2参量时和上述差分的大小小于第2参量时设定为不同的值,同时,通过变更上述第2参量而控制锐化的程度。
5.根据权利要求4所述的图象处理方法,其特征在于:从输入的图像数据中抽出表示该图像数据的特征的数据,通过将这些数据输入指定的算法而设定所述第2参量。
6.一种图象处理方法,是根据原图像的图像数据通过将把原图像的图像数据的分布作为函数所区域的各像素的2次微分取负而进行锐化处理的图象处理方法,其特征在于:对所述2次微分设置了上限和/或下限,通过变更决定该上限和/或下限的第3参量来控制锐化的程度。
7.根据权利要求6所述的图象处理方法,其特征在于:所述上限和/或下限由用色差微分值除上述第3参量的值决定。
8.根据权利要求6或7所述的图象处理方法,其特征在于:从输入的图像数据中抽出表示该图像数据的特征的数据,通过将这些数据输入指定的算法来设定所述第3参量。
9.根据权利要求1~8的任一权项所述的图象处理方法,其特征在于:对原图像的像素数据利用2次微分对各色成分进行锐化处理。
10.一种记录了图象处理程序的记录媒体,是记录了使计算机执行根据原图像的图像数据通过将把原图像的图像数据的分布作为函数所区域的各像素的2次微分取负而进行锐化处理的图象处理的图象处理程序的记录媒体,其特征在于:通过变更根据原图像的图像数据而取负的2次微分的大小的第1参量而使计算机执行控制锐化处理的程度的处理。
11.根据权利要求10所述的记录了图象处理程序的记录媒体,其特征在于:使计算机执行从输入的图像数据中抽出表示该图像数据的特征的数据并通过将这些数据输入指定的算法而决定上述第1参量的处理。
12.一种记录了图象处理程序的记录媒体,是记录了使计算机执行根据原图像的图像数据通过将把原图像的图像数据的分布作为函数所区域的各像素的2次微分取负而进行锐化处理的图象处理的图象处理程序的记录媒体,其特征在于:所述2次微分由注目像素与其附近的多个像素的像素值的差分的总和得到,同时使计算机执行对所述各差分乘以系数并且各系数根据对应的差分的大小而变化的处理。
13.根据权利要求12所述的记录了图象处理程序的记录媒体,其特征在于:所述系数在所述差分的大小大于第2参量时和上述差分的大小小于第2参量时设定为不同的值,同时,通过变更所述第2参量使计算机执行控制锐化的程度的处理。
14.根据权利要求13所述的记录了图象处理程序的记录媒体,其特征在于:从输入的图像数据中抽出表示该图像数据的特征的数据,通过将这些数据输入指定的算法而使计算机执行设定所述第2参量的处理。
15.一种记录了图象处理程序的记录媒体,是记录了使计算机执行根据原图像的图像数据通过将把原图像的图像数据的分布作为函数所区域的各像素的2次微分取负而进行锐化处理的图象处理的图象处理程序的记录媒体,其特征在于:对所述2次微分设置了上限和/或下限,通过变更决定该上限和/或下限的第3参量,使计算机执行控制锐化的程度的处理。
16.根据权利要求15所述的记录了图象处理程序的记录媒体,其特征在于:使计算机执行通过用色差微分值除所述第3参量的值决定上述上限和/或下限的处理。
17.根据权利要求15或16所述的记录了图象处理程序的记录媒体,其特征在于:从输入的图像数据中抽出表示该图像数据的特征的数据,通过将这些数据输入指定的算法而使计算机执行设定所述第3参量的处理。
18.根据权利要求10~17的任一权项所述的记录了图象处理程序的记录媒体,其特征在于:使计算机对原图像的像素数据执行对各色成分通过2次微分进行锐化处理的处理。
CNB001361295A 1999-12-27 2000-12-25 图象处理方法和记录了图象处理程序的记录媒体 Expired - Fee Related CN1167036C (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP37093899A JP2001186353A (ja) 1999-12-27 1999-12-27 画像処理方法および画像処理プログラムを記録した記録媒体
JP370938/1999 1999-12-27

Publications (2)

Publication Number Publication Date
CN1302045A true CN1302045A (zh) 2001-07-04
CN1167036C CN1167036C (zh) 2004-09-15

Family

ID=18497856

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB001361295A Expired - Fee Related CN1167036C (zh) 1999-12-27 2000-12-25 图象处理方法和记录了图象处理程序的记录媒体

Country Status (4)

Country Link
US (1) US6924839B2 (zh)
EP (2) EP1113390A3 (zh)
JP (1) JP2001186353A (zh)
CN (1) CN1167036C (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100394766C (zh) * 2004-03-05 2008-06-11 佳能株式会社 图像信号处理装置和图像信号处理方法
US8982409B2 (en) 2005-12-16 2015-03-17 Thomson Licensing Method, apparatus and system for providing reproducible digital imagery products from film content

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6614474B1 (en) * 1998-08-27 2003-09-02 Polycom, Inc. Electronic pan tilt zoom video camera with adaptive edge sharpening filter
US7088392B2 (en) * 2001-08-27 2006-08-08 Ramakrishna Kakarala Digital image system and method for implementing an adaptive demosaicing method
CN1817047A (zh) * 2003-06-30 2006-08-09 株式会社尼康 用于处理具有排列的不同色彩分量的图像的图像处理装置,图像处理程序,电子照相机,以及图像处理方法
JP2005141477A (ja) * 2003-11-06 2005-06-02 Noritsu Koki Co Ltd 画像鮮鋭化処理方法とこの方法を実施する画像処理装置
US20050163401A1 (en) * 2004-01-28 2005-07-28 Techwell Inc. Display image enhancement apparatus and method using adaptive interpolation with correlation
TWI263444B (en) * 2004-12-10 2006-10-01 Hon Hai Prec Ind Co Ltd System and method for focusing images automatically
GB2423660B (en) * 2005-02-25 2010-08-18 Imagination Tech Ltd A method and apparatus for controlling overshoot in a video enhancement system
US7595911B2 (en) * 2006-02-07 2009-09-29 Sharp Laboratories Of America, Inc. Methods and systems for descreening a digital image
JP2007226520A (ja) * 2006-02-23 2007-09-06 Nippon Computer System Co Ltd 画像処理プログラム及びこれが記録されたコンピュータ読み取り可能な記録媒体、ならびに画像処理装置
JP4677376B2 (ja) * 2006-07-07 2011-04-27 キヤノン株式会社 画像処理装置、画像処理方法、画像処理プログラム及び記憶媒体
US8442344B2 (en) * 2006-08-07 2013-05-14 Qualcomm Incorporated Adaptive spatial image filter for filtering image information
US8243351B2 (en) * 2007-12-21 2012-08-14 Canon Kabushiki Kaisha Image reading apparatus, multifunction printer apparatus, and image processing method
JP4989444B2 (ja) * 2007-12-21 2012-08-01 キヤノン株式会社 画像読取装置、マルチファンクションプリンタ装置、及び画像処理方法
JP4933418B2 (ja) * 2007-12-21 2012-05-16 キヤノン株式会社 画像読取装置、マルチファンクションプリンタ装置、及び画像処理方法
TW200939154A (en) 2008-03-14 2009-09-16 Holtek Semiconductor Inc An edge extraction and enhancement
US8687912B2 (en) * 2009-03-16 2014-04-01 Nikon Corporation Adaptive overshoot control for image sharpening
WO2011021518A1 (ja) 2009-08-19 2011-02-24 日本電気株式会社 画像処理方法、画像処理装置及びプログラム
JP6953183B2 (ja) * 2017-05-24 2021-10-27 キヤノン株式会社 画像処理装置、画像処理方法およびプログラム

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4571635A (en) * 1984-02-17 1986-02-18 Minnesota Mining And Manufacturing Company Method of image enhancement by raster scanning
US4794531A (en) * 1984-11-07 1988-12-27 Hitachi, Ltd Unsharp masking for image enhancement
US5245432A (en) * 1989-07-31 1993-09-14 Imageware Research And Development Inc. Apparatus and method for transforming a digitized signal of an image to incorporate an airbrush effect
WO1991003795A1 (en) * 1989-08-28 1991-03-21 Eastman Kodak Company Digital image noise reduction of luminance and chrominance based on overlapping planar approximation
JP3437226B2 (ja) * 1993-10-20 2003-08-18 キヤノン株式会社 画像処理方法及び装置
US6714665B1 (en) * 1994-09-02 2004-03-30 Sarnoff Corporation Fully automated iris recognition system utilizing wide and narrow fields of view
JPH09326026A (ja) * 1996-06-07 1997-12-16 Advantest Corp 画像処理方法及び画像処理装置
US5880767A (en) * 1996-09-11 1999-03-09 Hewlett-Packard Company Perceptual image resolution enhancement system
US6125213A (en) * 1997-02-17 2000-09-26 Canon Kabushiki Kaisha Image processing method, an image processing apparatus, and a storage medium readable by a computer
US6252995B1 (en) * 1997-08-25 2001-06-26 Fuji Photo Film Co., Ltd. Method of and apparatus for enhancing image sharpness
JP2000032299A (ja) * 1998-07-08 2000-01-28 Sony Corp 映像信号輪郭補正回路
US6665447B1 (en) 1999-07-26 2003-12-16 Hewlett-Packard Development Company, L.P. Method for enhancing image data by sharpening
US6804408B1 (en) * 1999-12-22 2004-10-12 Eastman Kodak Company Method for enhancing a digital image with noise-dependent control of texture
JP2003526999A (ja) * 2000-03-15 2003-09-09 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ ピーキングフィルタを有するビデオ装置
US6611627B1 (en) * 2000-04-24 2003-08-26 Eastman Kodak Company Digital image processing method for edge shaping
US6650790B1 (en) * 2000-06-09 2003-11-18 Nothshore Laboratories, Inc. Digital processing apparatus for variable image-size enlargement with high-frequency bandwidth synthesis

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100394766C (zh) * 2004-03-05 2008-06-11 佳能株式会社 图像信号处理装置和图像信号处理方法
US8982409B2 (en) 2005-12-16 2015-03-17 Thomson Licensing Method, apparatus and system for providing reproducible digital imagery products from film content
US9699388B2 (en) 2005-12-16 2017-07-04 Thomson Licensing Method, apparatus and system for providing reproducible digital imagery products

Also Published As

Publication number Publication date
CN1167036C (zh) 2004-09-15
US6924839B2 (en) 2005-08-02
EP1903501A2 (en) 2008-03-26
JP2001186353A (ja) 2001-07-06
EP1113390A3 (en) 2004-04-07
EP1903501A3 (en) 2008-07-30
US20010022622A1 (en) 2001-09-20
EP1113390A2 (en) 2001-07-04

Similar Documents

Publication Publication Date Title
CN1167036C (zh) 图象处理方法和记录了图象处理程序的记录媒体
CN1207895C (zh) 图象处理方法和装置
CN1259890C (zh) 图像处理装置和方法
CN1311404C (zh) 使用了图像生成历史信息的图像的逆光处理
CN1187955C (zh) 图像文件的输出图像调整
CN1114888C (zh) 图像处理方法、设备和系统
CN1671176A (zh) 修正图像变形的图像处理装置、修正摄影图像变形的摄影装置
CN1301492C (zh) 图像处理装置
CN1691740A (zh) 放大显示装置和放大图像控制装置
CN100345160C (zh) 一种控制平均亮度的直方图均衡方法
CN1554074A (zh) 在考虑其噪声的情况下改变数字图象的方法和系统
CN1542693A (zh) 图像处理方法、图像处理装置及图像处理程序
CN101032159A (zh) 图像处理设备及其方法,图像处理程序
CN1917578A (zh) 数据处理设备、数据处理方法及程序
CN1909671A (zh) 产生用户偏爱数据的方法和装置
CN1643936A (zh) 图像处理设备和方法、以及图像处理设备的制造方法
CN1871847A (zh) 信号处理系统、信号处理方法以及信号处理程序
CN1689340A (zh) 数字静像照相机与图像校正方法
CN1744675A (zh) 图像输入装置
CN1714372A (zh) 图像信号处理
CN1649392A (zh) 具有后加工设定功能的电子相机和定制该功能的处理程序
CN1703900A (zh) 摄像系统和播放系统
CN1737852A (zh) 图像处理装置、图像处理方法、图像处理程序和记录介质
CN1878232A (zh) 图像处理装置及图像处理方法
CN1806449A (zh) 图像处理设备、图像处理方法和计算机程序

Legal Events

Date Code Title Description
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C06 Publication
PB01 Publication
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20040915

Termination date: 20111225