JP2002245464A - Device and method for image evaluation, and program making computer evaluate image - Google Patents

Device and method for image evaluation, and program making computer evaluate image

Info

Publication number
JP2002245464A
JP2002245464A JP2001040855A JP2001040855A JP2002245464A JP 2002245464 A JP2002245464 A JP 2002245464A JP 2001040855 A JP2001040855 A JP 2001040855A JP 2001040855 A JP2001040855 A JP 2001040855A JP 2002245464 A JP2002245464 A JP 2002245464A
Authority
JP
Japan
Prior art keywords
image
evaluation
spatial
filter
image evaluation
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
JP2001040855A
Other languages
Japanese (ja)
Other versions
JP4040259B2 (en
Inventor
Kenji Kagitani
賢治 鎰谷
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.)
Ricoh Co Ltd
Original Assignee
Ricoh 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 Ricoh Co Ltd filed Critical Ricoh Co Ltd
Priority to JP2001040855A priority Critical patent/JP4040259B2/en
Publication of JP2002245464A publication Critical patent/JP2002245464A/en
Application granted granted Critical
Publication of JP4040259B2 publication Critical patent/JP4040259B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Image Analysis (AREA)

Abstract

PROBLEM TO BE SOLVED: To output an evaluated value which is correlated to the subjective evaluation of an image. SOLUTION: An image evaluating device 100 is equipped with an image data input part 101 which inputs spatial optical information of the image, a convolutional arithmetic part 103 which performs convolutional arithmetic between a spatial arithmetic filter stored in an arithmetic filter storage part 106 and a spatial arithmetic filter for the spatial optical information of the image inputted by the image data input part 101, and an image quality evaluation part 105 which predicts the image quality that a person feels according to the output of the convolutional arithmetic part 103.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は、ハードコピー画像
の品質を評価する画像評価装置、画像評価方法およびコ
ンピュータに画像を評価させるプログラムに関し、特
に、マルチチャンネルモデルに基づいたロジスティック
回帰分析により重み(係数)の決定された評価式を用い
て画像の品質を評価する画像評価装置、画像評価方法お
よびコンピュータに画像を評価させるプログラムに関す
る。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an image evaluation device for evaluating the quality of a hard copy image, an image evaluation method, and a program for causing a computer to evaluate an image. More particularly, the present invention relates to a logistic regression analysis based on a multi-channel model. The present invention relates to an image evaluation device, an image evaluation method, and a program for causing a computer to evaluate an image by using the evaluation formula for which the coefficient has been determined.

【0002】[0002]

【従来の技術】従来、人間の視覚機能においては、空間
周波数選択性を有する複数のチャンネル(以降において
適宜マルチチャンネルと称する)が存在することが知ら
れている。神経生理学的および解剖学的研究により、網
膜神経節細胞(ganglion cell)には、ほ
ぼ円対称の受容野(receptive field)
をもち、受容野中心部に興奮性領域を有しつつ周辺部に
抑制性領域を有するオン中心型細胞と、その逆の特性を
有するオフ中心型細胞が存在することが明らかにされて
いる。このオン中心型細胞とオフ中心型細胞が空間周波
数選択性のチャンネルの視覚経路をなしていると考えら
れている。図17は、チャンネルの受容野の概略形状を
示した説明図であり、受容野平面上に各受光体からの出
力をz軸方向にプロットしたものを水平方向から表示し
た図である。図はオン中心型細胞を示している。
2. Description of the Related Art It has been known that a plurality of channels having spatial frequency selectivity (hereinafter, appropriately referred to as multi-channels) exist in human visual functions. According to neurophysiological and anatomical studies, retinal ganglion cells have a substantially circularly symmetric receptive field.
It has been clarified that there are on-centre cells having an excitatory region in the central part of the receptive field and an inhibitory region in the peripheral part, and off-centre cells having the opposite characteristics. It is thought that the on-centered cell and the off-centered cell form the visual pathway of a spatial frequency selective channel. FIG. 17 is an explanatory view showing the schematic shape of the receptive field of the channel, and is a diagram in which the output from each photoreceptor is plotted on the plane of the receptive field in the z-axis direction and displayed from the horizontal direction. The figure shows on-centre cells.

【0003】一方、ハードコピー画像の画像品質を評価
する方法として、特開平7−325922号に示される
技術が挙げられる。すなわち、人間の視覚特性を考慮す
るために、物理量である空間周波数領域(パワースペク
トラム)に変換された画像情報に対して、人間の視覚の
MTF特性(VTF)を乗じた量を用いて画像品質を評
価するという方法が用いられている。
On the other hand, as a method of evaluating the image quality of a hard copy image, there is a technique disclosed in Japanese Patent Application Laid-Open No. 7-325922. That is, in order to take into account human visual characteristics, image quality is converted by using the amount obtained by multiplying the MTF characteristic (VTF) of human vision with image information converted into a spatial frequency domain (power spectrum), which is a physical quantity. Is used.

【0004】[0004]

【発明が解決しようとする課題】しかしながら、従来の
技術では以下の問題点があった。すなわち、従来の画像
品質の評価方法では、デジタル画像のテクスチャの知覚
と粒状性を統合することが困難であるという問題点があ
った。また、テクスチャのみによるエッジの検出が困難
であるという問題点があった。
However, the prior art has the following problems. That is, the conventional image quality evaluation method has a problem that it is difficult to integrate the perception of texture and the granularity of a digital image. In addition, there is a problem that it is difficult to detect an edge using only a texture.

【0005】これは、画質を知覚する人間の視覚系が非
常に非線形なシステムであるため、空間周波数特性は近
似的にしか有効でないことに起因すると考えられる。ま
た、視覚のMTF特性は人間の正弦波状の刺激に対する
コントラスト感度の閾値のデータを基に得られたもので
あるため、マルチチャンネルによる視覚現象の1側面を
記述したものにすぎないことに起因すると考えられる。
[0005] This is considered to be due to the fact that the spatial frequency characteristic is only approximately effective because the human visual system that perceives image quality is a very nonlinear system. In addition, since the visual MTF characteristic is obtained based on the data of the threshold value of the contrast sensitivity to a human sine wave-like stimulus, the visual MTF characteristic is merely one aspect of the visual phenomenon by multi-channel. Conceivable.

【0006】たとえば、濃度の同じ画像を二つ並べて見
た場合に、人間はそのテクスチャの差によって両画像の
見え方が異なると感じる。一方、従来の空間周波数特性
を用いた方法では同じ見え方をすると評価してしまうと
いう問題点があった。
For example, when two images having the same density are viewed side by side, a human perceives that the appearance of both images is different due to the difference in texture. On the other hand, the conventional method using the spatial frequency characteristic has a problem that the same appearance is evaluated as being evaluated.

【0007】換言すると、画像に対する主観評価と相関
のとれた評価値を出力することができない場合があると
いう問題点があった。
In other words, there is a problem that an evaluation value correlated with the subjective evaluation of an image may not be output.

【0008】本発明は上記に鑑みてなされたものであっ
て、画像に対する主観評価と相関のとれた評価値を出力
することを目的とする。
The present invention has been made in view of the above, and has as its object to output an evaluation value correlated with the subjective evaluation of an image.

【0009】[0009]

【課題を解決するための手段】上記の目的を達成するた
めに、請求項1に記載の画像評価装置は、画像の空間的
光学的情報を入力する画像入力手段と、前記画像入力手
段により入力された空間的光学的情報に対して使用する
空間演算フィルタ空間演算フィルタと、前記画像入力手
段により入力された画像の空間的光学的情報に対する前
記空間演算フィルタの畳み込み演算をおこなう畳込演算
手段と、前記畳込演算手段の出力を基に、人間が知覚す
る画像品質を予測する画質予測手段と、を具備したこと
を特徴とする。すなわち、請求項1にかかる発明は、人
間の視覚系の特性に基づく画像品質評価をおこなう。
According to another aspect of the present invention, there is provided an image evaluation apparatus, comprising: an image input unit for inputting spatial and optical information of an image; A spatial operation filter used for the obtained spatial optical information, a spatial operation filter, and a convolution operation unit for performing a convolution operation of the spatial operation filter on the spatial optical information of the image input by the image input unit. And image quality prediction means for predicting the image quality perceived by humans based on the output of the convolution operation means. That is, the invention according to claim 1 performs image quality evaluation based on characteristics of a human visual system.

【0010】また、請求項2に記載の画像評価装置は、
請求項1記載の画像評価装置において、前記空間演算フ
ィルタが等方性微分演算フィルタであることを特徴とす
る。すなわち、請求項2にかかる発明は、2次元の画像
構造に起因する画質要因を人間の視覚系の特性に基づい
て定量化する。
[0010] The image evaluation apparatus according to claim 2 is
2. The image evaluation device according to claim 1, wherein the spatial operation filter is an isotropic differential operation filter. That is, the invention according to claim 2 quantifies an image quality factor caused by a two-dimensional image structure based on characteristics of a human visual system.

【0011】また、請求項3に記載の画像評価装置は、
請求項1または2に記載の画像評価装置において、前記
空間演算フィルタが各々特性の相異なる複数のフィルタ
からなり、前記画質予測手段が前記複数のフィルタにそ
れぞれに対応した前記畳み込み演算手段からの出力すべ
てを用いて、人間が知覚する画像品質を予測することを
特徴とする。すなわち、請求項3にかかる発明は、様々
な尺度の画像構造に起因する画質要因を人間の視覚系の
特性に基づいて定量化し、画質を総合的に定量化する。
Further, the image evaluation device according to the third aspect is characterized in that:
3. The image evaluation device according to claim 1, wherein the spatial operation filter includes a plurality of filters having different characteristics, and the image quality prediction unit outputs from the convolution operation unit corresponding to each of the plurality of filters. It is characterized by predicting the image quality perceived by humans using everything. That is, the invention according to claim 3 quantifies the image quality factors caused by the image structures of various scales based on the characteristics of the human visual system, and quantifies the image quality comprehensively.

【0012】また、請求項4に記載の画像評価装置は、
請求項3に記載の画像評価装置において、前記複数のフ
ィルタが人間の視覚系の空間周波数特性に対して最適化
されていることを特徴とする。すなわち、請求項4にか
かる発明は、人間の視覚系の空間周波数特性に基づく、
より精緻な画質の定量化をおこなう。
[0012] The image evaluation device according to claim 4 is
4. The image evaluation device according to claim 3, wherein the plurality of filters are optimized for spatial frequency characteristics of a human visual system. That is, the invention according to claim 4 is based on the spatial frequency characteristics of the human visual system.
Perform more precise image quality quantification.

【0013】また、請求項5に記載の画像評価装置は、
請求項1〜4の何れか一つに記載の画像評価装置におい
て、前記空間演算フィルタとして各々特性の相異なる5
個のフィルタを具備したことを特徴とする。すなわち、
請求項5にかかる発明は、人間の視覚系のマルチチャン
ネルモデルに基づく、より精緻な画質の定量化をおこな
う。
[0013] The image evaluation apparatus according to claim 5 is
5. The image evaluation device according to claim 1, wherein the spatial calculation filters have different characteristics. 5.
It is characterized by having a number of filters. That is,
The invention according to claim 5 performs more precise quantification of image quality based on a multi-channel model of the human visual system.

【0014】また、請求項6に記載の画像評価装置は、
請求項3、4または5に記載の画像評価装置において、
前記複数のフィルタのうち必要な数個を選択するフィル
タ選択手段を具備したことを特徴とする。すなわち、請
求項6にかかる発明は、画像構造が各々相異なる視覚的
尺度に分布するような複数種の画像欠陥に対応した知覚
を人間の視覚系のマルチチャンネルモデルに基づいて、
より高速かつ精緻に定量化する。
Further, the image evaluation device according to claim 6 is
The image evaluation device according to claim 3, 4 or 5,
The image processing apparatus further includes a filter selection unit that selects a required number of the plurality of filters. That is, the invention according to claim 6 provides a perception corresponding to a plurality of types of image defects in which image structures are distributed on different visual scales based on a multi-channel model of a human visual system.
Quantify faster and more precisely.

【0015】また、請求項7に記載の画像評価装置は、
請求項1〜6の何れか一つに記載の画像評価装置におい
て、前記空間演算フィルタのサイズを各々可変とするフ
ィルタサイズ変更手段を具備したことを特徴とする。す
なわち、請求項7にかかる発明は、被評価画像の画像密
度に依存しない画質評価をおこなう。これは、マルチチ
ャンネルモデルに基づく画像品質評価をおこなうために
は被評価画像の画像密度が網膜上の結像密度にほぼ等し
いことが望ましいからである。
Further, the image evaluation device according to claim 7 is
7. The image evaluation apparatus according to claim 1, further comprising: a filter size changing unit configured to change a size of the spatial operation filter. That is, the invention according to claim 7 performs image quality evaluation independent of the image density of the evaluated image. This is because, in order to perform image quality evaluation based on the multi-channel model, it is desirable that the image density of the evaluated image be approximately equal to the image density on the retina.

【0016】また、請求項8に記載の画像評価装置は、
評価すべき画像の画像データを入力する画像データ入力
手段と、前記画像データ入力手段により入力された画像
データに対して使用する空間演算フィルタを選択する空
間演算フィルタ選択手段と、前記画像データ入力手段に
より入力された画像データに対して前記空間演算フィル
タ選択手段により選択された空間演算フィルタを用いて
畳み込み演算をおこなう畳込演算手段と、前記畳込演算
手段によりおこなわれた畳み込み演算の結果に基づいて
前記画像の品質を評価する画像品質評価手段と、を具備
したことを特徴とする。すなわち、請求項8にかかる発
明は、人間の視覚系の特性に基づく画像品質評価をおこ
なう。
Further, the image evaluation device according to the present invention is characterized in that:
Image data input means for inputting image data of an image to be evaluated, spatial calculation filter selection means for selecting a spatial calculation filter to be used for the image data input by the image data input means, and the image data input means A convolution operation unit that performs a convolution operation on the image data input by using the spatial operation filter selected by the spatial operation filter selection unit, based on a result of the convolution operation performed by the convolution operation unit. Image quality evaluation means for evaluating the quality of the image. That is, the invention according to claim 8 performs image quality evaluation based on characteristics of a human visual system.

【0017】また、請求項9に記載の画像評価装置は、
請求項8記載の画像評価装置において、前記空間演算フ
ィルタを等方性微分演算フィルタとしたことを特徴とす
る。すなわち、請求項9にかかる発明は、2次元の画像
構造に起因する画質要因を人間の視覚系の特性に基づい
て定量化する。
Further, the image evaluation apparatus according to claim 9 is
9. The image evaluation device according to claim 8, wherein the spatial operation filter is an isotropic differential operation filter. That is, the invention according to claim 9 quantifies an image quality factor caused by a two-dimensional image structure based on characteristics of a human visual system.

【0018】また、請求項10に記載の画像評価装置
は、請求項9に記載の画像評価装置において、前記等方
性微分演算フィルタを、関数LG(r)=(−1/(π
σ4))・(1−r2/(2σ2))・exp(−r2
(2σ2))(但し、rは、r=sqrt(x2+y2
で与えられる距離であり、x,yは前記評価すべき画像
の画像データの或る位置を中心とした座標系における座
標であり、σは視覚系のマルチチャンネルモデルでいう
5種類のチャンネルに対応した定数である)としたこと
を特徴とする。すなわち、請求項10にかかる発明は、
人間の視覚系のマルチチャンネルモデルに基づく、より
精細な画質の定量化をおこなう。
According to a tenth aspect of the present invention, in the image evaluation apparatus of the ninth aspect, the isotropic differential operation filter is provided with a function LG (r) = (− 1 / (π
σ 4 )) · (1−r 2 / (2σ 2 )) · exp (−r 2 /
(2σ 2 )) (where r is r = sqrt (x 2 + y 2 )
Where x and y are coordinates in a coordinate system centered on a certain position of the image data of the image to be evaluated, and σ corresponds to five types of channels in a multi-channel model of a visual system. This is a characteristic constant). That is, the invention according to claim 10 is
Quantify finer image quality based on multi-channel model of human visual system.

【0019】また、請求項11に記載の画像評価装置
は、二つの画像を比較して同じに見えるか差があるよう
に見えるかの主観的な評価に沿った評価を出力する画像
評価装置であって、評価すべき二つの画像の画像データ
を入力する画像データ入力手段と、前記画像データ入力
手段により入力されたデータに基づいて、テクスチャ量
の差の絶対値ΔTとテクスチャ量の平均値Tmeanと
を算出するテクスチャ量算出手段と、前記評価すべき二
つの画像の網点面積率の差の絶対値ΔMと網点面積率の
平均値Mmeanとを取得する網点面積率取得手段と、
前記テクスチャ量算出手段により算出されたテクスチャ
量の差の絶対値ΔTとテクスチャ量の平均値Tmea
n、および、前記網点面積率取得手段により取得された
網点面積率の差の絶対値ΔMと網点面積率の平均値Mm
eanとを用いて、式(1)により画像評価値Dを算出
する画像評価値算出手段と、を具備したことを特徴とす
る(但しp1〜p9は、予めおこなわれた主観評価実験に
より決定された定数とする)。すなわち、請求項11に
かかる発明は、巨視的には同じ濃度である画像のテクス
チャによる見え方の差を算出する。
An image evaluation apparatus according to an eleventh aspect is an image evaluation apparatus for comparing two images and outputting an evaluation according to a subjective evaluation of whether the two images look the same or appear to have a difference. An image data input unit for inputting image data of two images to be evaluated; and an absolute value ΔT of a difference between texture amounts and an average value Tmean of texture amounts based on the data input by the image data input unit. Texture amount calculating means for calculating the halftone dot area ratio of the two images to be evaluated, and a halftone dot area ratio obtaining means for obtaining an average value Mmean of the halftone dot area ratio,
The absolute value ΔT of the difference between the texture amounts calculated by the texture amount calculating means and the average value Tmea of the texture amounts
n, the absolute value ΔM of the difference between the halftone dot area ratios obtained by the halftone dot area ratio obtaining means and the average value Mm of the halftone dot area ratios
and an image evaluation value calculating means for calculating an image evaluation value D according to equation (1) using E.ean (where p 1 to p 9 are based on a subjective evaluation experiment conducted in advance. The determined constant). That is, the invention according to claim 11 calculates a difference in appearance due to the texture of an image having the same density macroscopically.

【0020】また、請求項12に記載の画像評価装置
は、請求項11に記載の画像評価装置において、さら
に、前記画像評価値算出手段により算出された画像評価
値Dを、式(2)に基づいて画像評価値PDに補正する
画像評価値補正手段を具備したことを特徴とする(但し
q1〜q4は、回帰分析により決定された定数とす
る)。すなわち、請求項12にかかる発明は、画像に対
する主観評価と高精度に相関のとれた評価値を出力す
る。
According to a twelfth aspect of the present invention, in the image evaluation apparatus according to the eleventh aspect, the image evaluation value D calculated by the image evaluation value calculation means is further expressed by the following equation (2). An image evaluation value correction unit for correcting the image evaluation value PD based on the image evaluation value PD is provided (however, q1 to q4 are constants determined by regression analysis). That is, the invention according to claim 12 outputs an evaluation value highly correlated with the subjective evaluation of the image.

【0021】また、請求項13に記載の画像評価装置
は、請求項1〜12に記載の画像評価装置において、前
記画像データ入力手段により入力される画像データの入
力密度を、人間が当該画像を観察する際の網膜上の結像
密度と略等しくなるように調整する入力密度調節手段を
備えたことを特徴とする。すなわち、請求項13にかか
る発明は、被評価画像の画像密度と網膜上の画像密度を
略等しくする。
According to a thirteenth aspect of the present invention, in the image evaluation apparatus according to the first aspect, the input density of the image data input by the image data input means is determined by a person. Input density adjusting means for adjusting the image density on the retina at the time of observation so as to be substantially equal to the density is provided. That is, in the invention according to claim 13, the image density of the evaluated image is made substantially equal to the image density on the retina.

【0022】また、請求項14に記載の画像評価方法
は、画像の空間的光学的情報に対して空間演算フィルタ
を畳み込んで得られる出力を基に、人間が知覚する画像
品質を予測することを特徴とする。すなわち、請求項1
4にかかる発明は、人間の視覚系の特性に基づく画像品
質評価をおこなう。
According to a fourteenth aspect of the present invention, there is provided an image evaluation method for predicting an image quality perceived by a human based on an output obtained by convolving a spatial and optical information of an image with a spatial operation filter. It is characterized by. That is, claim 1
The invention according to 4 performs image quality evaluation based on characteristics of a human visual system.

【0023】また、請求項15に記載の画像評価方法
は、請求項14記載の画像評価方法において、前記空間
演算フィルタが等方性微分演算フィルタであることを特
徴とする。すなわち、請求項15にかかる発明は、2次
元の画像構造に起因する画質要因を人間の視覚系の特性
に基づいて定量化する。
According to a fifteenth aspect of the present invention, in the image evaluation method of the fourteenth aspect, the spatial operation filter is an isotropic differential operation filter. That is, the invention according to claim 15 quantifies the image quality factor caused by the two-dimensional image structure based on the characteristics of the human visual system.

【0024】また、請求項16に記載の画像評価方法
は、請求項14または15に記載の画像評価方法におい
て、空間演算フィルタとして各々特性の相異なる複数の
空間演算フィルタ群を用い、当該空間演算フィルタ群に
より畳み込んで得られる各出力のすべてを用いて、人間
が知覚する画像品質を予測することを特徴とする。すな
わち、請求項16にかかる発明は、様々な尺度の画像構
造に起因する画質要因を人間の視覚系の特性に基づいて
定量化し、画質を総合的に定量化する。
The image evaluation method according to claim 16 is the image evaluation method according to claim 14 or 15, wherein a plurality of spatial operation filter groups having different characteristics are used as the spatial operation filters. It is characterized in that the image quality perceived by a human is predicted by using all of the outputs obtained by convolution by the filter group. That is, the invention according to claim 16 quantifies the image quality factors caused by the image structures of various scales based on the characteristics of the human visual system, and comprehensively quantifies the image quality.

【0025】また、請求項17に記載の画像評価方法
は、請求項16に記載の画像評価方法において、前記空
間演算フィルタ群が人間の視覚系の空間周波数特性に対
して最適化されていることを特徴とする。すなわち、請
求項17にかかる発明は、人間の視覚系の空間周波数特
性に基づく、より精緻な画質の定量化をおこなう。
According to a seventeenth aspect of the present invention, in the image evaluation method of the sixteenth aspect, the spatial operation filter group is optimized for a spatial frequency characteristic of a human visual system. It is characterized by. That is, the invention according to claim 17 performs more precise quantification of image quality based on the spatial frequency characteristics of the human visual system.

【0026】また、請求項18に記載の画像評価方法
は、請求項14〜17の何れか一つに記載の画像評価方
法において、前記空間演算フィルタとして各々特性の相
異なる5個の空間演算フィルタを用いることを特徴とす
る。すなわち、請求項18にかかる発明は、人間の視覚
系のマルチチャンネルモデルに基づく、より精緻な画質
の定量化をおこなう。
The image evaluation method according to claim 18 is the image evaluation method according to any one of claims 14 to 17, wherein five spatial operation filters having different characteristics are used as the spatial operation filters. Is used. That is, the invention according to claim 18 performs more precise quantification of image quality based on a multi-channel model of the human visual system.

【0027】また、請求項19に記載の画像評価方法
は、請求項16、17または18に記載の画像評価方法
において、前記空間演算フィルタのうち必要な数個を選
択して使用することを特徴とする。すなわち、請求項1
9にかかる発明は、画像構造が各々相異なる視覚的尺度
に分布するような複数種の画像欠陥に対応した知覚を人
間の視覚系のマルチチャンネルモデルに基づいて、より
高速かつ精緻に定量化する。
According to a nineteenth aspect of the present invention, in the image evaluation method of the sixteenth, seventeenth, or eighteenth aspect, a necessary number of the spatial operation filters are selected and used. And That is, claim 1
The invention according to 9 quantifies perception corresponding to a plurality of types of image defects in which image structures are distributed on different visual scales faster and more precisely based on a multi-channel model of the human visual system. .

【0028】また、請求項20に記載の画像評価方法
は、請求項14〜19の何れか一つに記載の画像評価方
法において、前記空間演算フィルタのサイズが各々可変
であることを特徴とする。すなわち、請求項7にかかる
発明は、被評価画像の画像密度に依存しない画質評価を
おこなう画像評価装置を提供する。これは、マルチチャ
ンネルモデルに基づく画像品質評価をおこなうためには
被評価画像の画像密度が網膜上の結像密度にほぼ等しい
ことが望ましいからである。
According to a twentieth aspect of the present invention, in the image evaluation method according to any one of the fourteenth to nineteenth aspects, the size of the spatial operation filter is variable. . That is, the invention according to claim 7 provides an image evaluation device that performs image quality evaluation independent of the image density of the evaluated image. This is because, in order to perform image quality evaluation based on the multi-channel model, it is desirable that the image density of the evaluated image be approximately equal to the image density on the retina.

【0029】また、請求項21に記載のプログラムは、
コンピュータに画像評価を実行させるプログラムであっ
て、コンピュータを、前記請求項1〜請求項13の何れ
か一つに記載の各手段として機能させることを特徴とす
る。すなわち、請求項21にかかる発明は、請求項1〜
請求項13の何れか一つに記載の各手段をコンピュータ
に実行させるプログラムであるので、コンピュータを用
いて、画像に対する主観評価と相関のとれた評価値を出
力することができる。
Further, the program according to claim 21 is:
A program for causing a computer to execute image evaluation, wherein the program causes a computer to function as each unit according to any one of claims 1 to 13. That is, the invention according to claim 21 is based on claims 1 to
Since the computer is a program that causes a computer to execute the means described in any one of the above aspects, it is possible to output an evaluation value that is correlated with the subjective evaluation of the image using the computer.

【0030】[0030]

【発明の実施の形態】以下、本発明の実施の形態を図面
を参照しながら詳細に説明する。 実施の形態1.実施の形態1では、まず、視覚系のモデ
ルについて概説し、つぎに、画像評価装置をパーソナル
コンピュータ(以降適宜PCと称する)に適用した例を
説明し、つづいて実験例を、最後に処理流れについて説
明する。
Embodiments of the present invention will be described below in detail with reference to the drawings. Embodiment 1 FIG. In the first embodiment, first, an outline of a model of a visual system will be described. Next, an example in which the image evaluation apparatus is applied to a personal computer (hereinafter, appropriately referred to as a PC) will be described. Will be described.

【0031】(視覚系のモデルについて)ここでは、視
覚のマルチチャンネル理論と受容野について概説する。
1968年にCampbellとRobsonが心理物
理学的実験の結果をもとに、人間の視覚機能における空
間周波数選択性を有するチャンネルの存在を結論づけて
以来、この空間周波数選択性を有するチャンネルについ
ては様々な心理物理学的研究がおこなわれてきた。
(About the Model of the Visual System) Here, the multichannel theory of vision and the receptive field will be outlined.
Since Campbell and Robson concluded the existence of a channel with spatial frequency selectivity in human visual function based on the results of psychophysical experiments in 1968, various channels with this spatial frequency selectivity have been Psychophysical studies have been performed.

【0032】一方、網膜神経節細胞には、略円対称の受
容野をもち、受容野中心部に興奮性領域を有し周辺部に
抑制領域を有するオン中心型細胞と、その逆のオフ中心
型細胞が存在することが、神経生理学的、および解剖学
的研究の結果明らかにされており、これらが空間周波数
選択性のチャンネルの視覚系路をなしていると考えられ
ている。
On the other hand, retinal ganglion cells have a substantially circularly symmetric receptive field, an on-center type cell having an excitatory region in the central portion of the receptive field and an inhibitory region in the peripheral portion, and the opposite off-center type cell. The presence of type cells has been demonstrated by neurophysiological and anatomical studies and is thought to form the visual pathway of spatially frequency-selective channels.

【0033】さらに、Wilsonらは視覚の検出域の
データを基礎にして、人間におけるマルチチャンネルメ
カニズムの定量的モデルを構築した。Wilsonらは
これらのチャンネルに対応する4種に分類される受容野
をそれぞれ二つのガウス関数の差(DOG:Diffe
rence of Gaussians)として記述
し、各チャンネルが刺激を検出する確率が空間的確率加
算の法則(Low ofprobability su
mmation)に従うならば、0.25〜16.0
[cycle/deg]の範囲の1次元の空間的パター
ンの閾検出に関する心理物理学的データを基本的に説明
できることを証明した。
In addition, Wilson et al. Built a quantitative model of the multi-channel mechanism in humans based on data in the visual detection zone. Wilson et al. Described the receptive field classified into four types corresponding to these channels as the difference between two Gaussian functions (DOG: Diffe
The probability of detecting a stimulus in each channel is determined by the law of spatial probability addition (Low of probability su).
MMation), 0.25 to 16.0
It has been proved that the psychophysical data concerning the threshold detection of a one-dimensional spatial pattern in the range of [cycle / deg] can be basically explained.

【0034】すなわち、Wilsonらは、視野の各点
において存在する1次元の4種類のチャンネル(マルチ
チャンネル)を提案した。これらは受容野の中心的領域
の大きさの小さい順にN,S,T,Uと名付けられた。
これらのチャンネルの大きさは偏心度、すなわち、視角
で表された中心窩からの距離と共に線形に増大する。ま
た、1次元的にはこの様な受容野を有するチャンネルは
網膜に写った外界の光の強度変化の検出器の役割を担っ
ていると考えられている。
That is, Wilson et al. Proposed four types of one-dimensional channels (multi-channels) existing at each point of the visual field. These were named N, S, T, U in ascending order of the size of the central area of the receptive field.
The size of these channels increases linearly with eccentricity, the distance from the fovea expressed in viewing angles. One-dimensionally, a channel having such a receptive field is considered to play a role of a detector for a change in the intensity of external light reflected on the retina.

【0035】さらに、この定量的モデルはMarらによ
り2次元的に拡張、補強されている。ここで、Mar
は、Wilsonらのマルチチャンネル理論を2次元に
展開するに当たり、受容野中心部が1個の錐体によって
構成されるミジェット神経節細胞に対応する最小チャン
ネルの存在を予測し、計5つのチャンネルでモデルを構
築した。
Further, this quantitative model is extended and reinforced two-dimensionally by Mar et al. Where Mar
Predicts the existence of the smallest channel corresponding to midget ganglion cells in which the central part of the receptive field is formed by a single cone in developing Wilson et al.'S multi-channel theory in two dimensions. The model was built.

【0036】Marらは、刺激の強度分布を検出する最
良のフィルタはΔGフィルタであるとする。ここでΔは
ラプラス演算子であり、Gは標準偏差σ(但しσは各チ
ャンネルによって固有の5つの値をとる)をもつ二次元
ガウス分布である。GおよびΔGは下式(3)、(4)
により表記される。
Assume that the best filter for detecting the intensity distribution of a stimulus is a ΔG filter. Here, Δ is a Laplace operator, and G is a two-dimensional Gaussian distribution having a standard deviation σ (where σ takes five unique values for each channel). G and ΔG are expressed by the following equations (3) and (4).
Is represented by

【数3】 (Equation 3)

【数4】 (Equation 4)

【0037】以降においては、MarらによるGもしく
はΔGの関数プロファイルを用いた当方性微分演算フィ
ルタにより、マルチチャンネルモデルに裏打ちされた画
像評価装置、画像評価方法およびコンピュータに画像を
評価させるプログラムについて説明する。
In the following, an image evaluation apparatus, an image evaluation method, and a program for causing a computer to evaluate an image, which are backed by a multi-channel model, by an isotropic differential operation filter using a G or ΔG function profile by Mar et al. Will be described. I do.

【0038】(画像評価装置について)図1は、実施の
形態1の画像評価装置をPCに適用した場合の機能ブロ
ック図であり、図2は、実施の形態1の画像評価装置の
ハードウェア構成の一例を示した説明図である。
(Regarding Image Evaluation Apparatus) FIG. 1 is a functional block diagram when the image evaluation apparatus of the first embodiment is applied to a PC, and FIG. 2 is a hardware configuration of the image evaluation apparatus of the first embodiment. FIG. 4 is an explanatory diagram showing an example of the embodiment.

【0039】PC100は、評価対象である画像の画像
データを入力する画像データ入力部101と、画像デー
タ入力部101により入力された画像データに対して使
用する空間演算フィルタを選択する演算フィルタ選択部
102と、画像データ入力部101により入力された画
像データに対して演算フィルタ選択部102により選択
された空間演算フィルタを用いて畳み込み演算をおこな
う畳込演算部103と、畳込演算部103によりおこな
われた畳み込み演算の結果に基づいて評価対象の画像の
品質の評価値を出力する画像品質評価部104と、を有
する。なお、以降においてはPC100を画像評価装置
100と称する場合がある。
The PC 100 includes an image data input unit 101 for inputting image data of an image to be evaluated, and an operation filter selection unit for selecting a spatial operation filter to be used for the image data input by the image data input unit 101. 102, a convolution operation unit 103 that performs a convolution operation on the image data input by the image data input unit 101 using the spatial operation filter selected by the operation filter selection unit 102, and a convolution operation unit 103. And an image quality evaluation unit 104 that outputs an evaluation value of the quality of the image to be evaluated based on the result of the convolution operation. In the following, the PC 100 may be referred to as the image evaluation device 100.

【0040】なお、画像品質評価部104では、評価す
べき画像の網点面積率を取得する網点面積率取得部10
5からの出力値を用いて画像の品質の評価値を出力す
る。なお、使用の態様によっては網点面積率以外の画像
の物理的属性、たとえば反射率などを用いてもよい。ま
た、演算フィルタ選択部102が選択する演算フィルタ
は演算フィルタ格納部106に格納されているものとす
る。演算フィルタ格納部106には、空間演算フィル
タ、たとえば、等方性微分演算フィルタやラプラシアン
などを格納する。以降においては、前述のΔGを用いた
例を説明するが、使用の態様によっては、DOGを用い
てもよい。
The image quality evaluation unit 104 obtains a dot area ratio of an image to be evaluated.
The evaluation value of the quality of the image is output using the output value from 5. Depending on the mode of use, a physical attribute of the image other than the dot area ratio, for example, a reflectance may be used. It is assumed that the operation filter selected by the operation filter selection unit 102 is stored in the operation filter storage unit 106. The operation filter storage unit 106 stores a spatial operation filter, for example, an isotropic differential operation filter, Laplacian, or the like. Hereinafter, an example using ΔG described above will be described, but DOG may be used depending on the mode of use.

【0041】PC100は、ハードウェア構成として、
畳み込み演算を実行することを初めとし画像の品質の評
価値を算出する演算をおこなうCPU201(図2参
照)と、CPU201のワークエリアであるRAM20
2と、OSを含み様々なプログラムないしソフトウェア
を格納し、また、評価対象の画像の画像データや、畳み
込み演算をおこなう際の空間演算フィルタや演算の際の
係数を格納するハードディスク203を有する。
The PC 100 has the following hardware configuration.
A CPU 201 (see FIG. 2) for performing operations for calculating an evaluation value of image quality including execution of a convolution operation, and a RAM 20 which is a work area of the CPU 201
2 and a hard disk 203 that stores various programs or software including an OS, and stores image data of an image to be evaluated, a spatial operation filter for performing a convolution operation, and a coefficient for an operation.

【0042】また、PC100は、評価すべき画像を表
示するCRT204と、CRT204の出力制御をおこ
なうビデオカード205と、評価すべき画像を印刷する
プリンタ206と、各種の指示をおこなうキーボード2
07およびマウス208と、プリンタ206、キーボー
ド207およびマウス208の入出力を制御するI/F
209と、評価すべき画像を読み込むスキャナ210
と、以上の各部を接続するバス211を有する。
The PC 100 includes a CRT 204 for displaying an image to be evaluated, a video card 205 for controlling the output of the CRT 204, a printer 206 for printing the image to be evaluated, and a keyboard 2 for giving various instructions.
I / F that controls input and output of a printer 07, a mouse 208, and a printer 206, a keyboard 207, and a mouse 208.
209 and a scanner 210 for reading an image to be evaluated
And a bus 211 for connecting the above components.

【0043】なお、CRT204、プリンタ206は、
評価すべき所定の画像を表示もしくは出力するが、これ
は、各種の係数決めをおこなう際もしくは回帰分析をお
こなうためのサンプリングのために用いてもよい。同様
に、スキャナ210は、サンプリングの際に用いるもの
であってもよい。換言すれば、PC100により、これ
らの入力機器および出力機器の評価をおこなうことがで
きる。
Note that the CRT 204 and the printer 206
A predetermined image to be evaluated is displayed or output, which may be used when determining various coefficients or for sampling for performing regression analysis. Similarly, the scanner 210 may be used for sampling. In other words, the input device and the output device can be evaluated by the PC 100.

【0044】ここで、ハードディスク203の格納対象
について説明する。ハードディスク203は、PC10
0の基本動作を制御するOS231と、畳み込み演算や
各種のフィルタや係数の読み出しをおこない画像の品質
の評価値を出力するプログラムである画像品質評価アプ
リケーション232と、を有する。
Here, the storage target of the hard disk 203 will be described. The hard disk 203 is a PC 10
It has an OS 231 for controlling the basic operation of 0, and an image quality evaluation application 232 which is a program for performing a convolution operation, reading various filters and coefficients, and outputting an evaluation value of image quality.

【0045】ハードディスク203は、評価すべき画像
データを格納する画像データ格納部233と、畳み込み
演算をおこなう際に参照するフィルタ格納部234と、
畳み込み演算をおこなう際に使用する係数や、画像の品
質を評価する際に使用する係数を格納する係数格納部2
35と、を有する。図に示したように、画像データ格納
部233には評価対象画像データ1、評価対象画像デー
タ2、・・・が格納され、フィルタ格納部234には、
微分演算フィルタ1、微分演算フィルタ2、・・・が格
納され、係数格納部235には、係数1、係数2、・・
・が格納されている。なお、係数には、後述するt1
5、A1〜A5、p1〜p9、q1〜q4等を格納する。
The hard disk 203 includes an image data storage unit 233 for storing image data to be evaluated, a filter storage unit 234 to be referred to when performing a convolution operation,
Coefficient storage unit 2 for storing coefficients used when performing a convolution operation and coefficients used when evaluating the quality of an image
35. As shown in the figure, the image data storage unit 233 stores evaluation target image data 1, evaluation target image data 2,..., And the filter storage unit 234 stores
The differential operation filter 1, the differential operation filter 2,... Are stored, and the coefficient storage unit 235 stores the coefficient 1, the coefficient 2,.
・ Is stored. It is to be noted that the coefficient, which will be described later t 1 ~
t 5, and stores the A 1 ~A 5, p 1 ~p 9, q 1 ~q 4 or the like.

【0046】つぎに、各部の内容を説明する。 (画像評価装置100:画像データ入力部101の内
容)画像データ入力部101は、評価すべき画像の画像
データを入力する。画像データ入力部101は、たとえ
ば300mm等所定の観察距離にて人間が被評価画像を
観察する際の人間の網膜上の結像密度と近似的に等しい
画像入力密度を有する画像を入力する。なお、予め電子
的に作成された画像を入力してもよい。画像データ入力
部101は、たとえば、スキャナ210やハードディス
ク203(特に、画像データ格納部233)によりその
機能を実現することができる。なお、使用の態様によっ
ては、マイクロ濃度計(図示せず)により実現すること
もできる。
Next, the contents of each section will be described. (Image Evaluation Device 100: Contents of Image Data Input Unit 101) The image data input unit 101 inputs image data of an image to be evaluated. The image data input unit 101 inputs an image having an image input density approximately equal to the image formation density on the human retina when a human observes the evaluated image at a predetermined observation distance such as 300 mm. Note that an image created electronically in advance may be input. The function of the image data input unit 101 can be realized by, for example, the scanner 210 or the hard disk 203 (in particular, the image data storage unit 233). In addition, depending on the mode of use, it can also be realized by a microdensitometer (not shown).

【0047】(画像評価装置100:演算フィルタ選択
部102の内容)演算フィルタ選択部102は、畳み込
み演算に使用する演算フィルタを選択し、畳込演算部1
03に引き渡す。演算フィルタは前述したGもしくはΔ
G(式(3)もしくは式(4)を参照)を用いる。な
お、式GもしくはΔG中のσは、各チャンネルに従って
定められた定数であり、5つのチャンネルがあるとする
Marのモデルの場合はσは5つあるということにな
る。以降においては、適宜σ1〜σ5(σ1最小、σ5
大)と表記することとする。
(Image Evaluation Apparatus 100: Contents of Arithmetic Filter Selection Unit 102) The arithmetic filter selection unit 102 selects an arithmetic filter to be used for the convolution operation.
Hand over to 03. The operation filter is G or Δ
G (see Equation (3) or Equation (4)) is used. In the expression G or ΔG, σ is a constant determined according to each channel. In the case of a Mar model in which there are five channels, there are five σ. Hereinafter, they will be referred to as σ 1 to σ 51 minimum, σ 5 maximum) as appropriate.

【0048】従って、演算フィルタ選択部102は、評
価すべき画像に応じて、すべてのσに対応したフィルタ
を選択する場合もあれば、一つのσ、たとえばσ2に対
応したフィルタを選択する場合もある。演算フィルタ選
択部102は、たとえば、ハードディスク203(特に
画像品質評価アプリケーション232とフィルタ格納部
234)、CPU201によりその機能を実現すること
ができる。なお、選択される演算フィルタは演算フィル
タ格納部106に格納されており、演算フィルタ格納部
106は、ハードディスク203(特にフィルタ格納部
234)によりその機能を実現することができる。
Therefore, the calculation filter selection unit 102 may select all the filters corresponding to σ depending on the image to be evaluated, or select one σ, for example, a filter corresponding to σ 2 , depending on the image to be evaluated. There is also. The function of the calculation filter selection unit 102 can be realized by, for example, the hard disk 203 (in particular, the image quality evaluation application 232 and the filter storage unit 234) and the CPU 201. The selected operation filter is stored in the operation filter storage unit 106, and the function of the operation filter storage unit 106 can be realized by the hard disk 203 (in particular, the filter storage unit 234).

【0049】(画像評価装置100:畳込演算部103
の内容)畳込演算部103は、入力された画像データに
選択された演算フィルタを施し画像の品質の評価値を算
出する際に使用される基礎データを出力する。すなわ
ち、畳込演算部103は、画像データに畳み込み演算を
施す。後述する例では、基礎データとして画像のテクス
チャを表現する物理量を求めている。畳込演算部103
は、たとえば、ハードディスク203(特に画像品質評
価アプリケーション232)、CPU201、RAM2
02によりその機能を実現することができる。
(Image evaluation device 100: convolution operation unit 103)
The convolution operation unit 103 outputs the basic data used when applying the selected operation filter to the input image data and calculating the evaluation value of the image quality. That is, the convolution operation unit 103 performs a convolution operation on the image data. In an example described later, a physical quantity expressing a texture of an image is obtained as basic data. Convolution operation unit 103
Are, for example, the hard disk 203 (especially the image quality evaluation application 232), the CPU 201, the RAM 2
02 can realize the function.

【0050】(画像評価装置100:画像品質評価部1
04)画像品質評価部104は、畳み込み演算の結果に
基づいて画像の品質の評価値を出力する。後述する例で
は、テクスチャ量と網点面積率を用いて、画像に対する
主観評価と相関のとれた評価値を出力する。具体的には
式(1)に表される評価値Dを出力する。この際使用す
る係数p1〜p9は、予め求めておいた係数であり、係数
格納部235に格納されている。画像品質評価部104
は、たとえば、ハードディスク203(特に画像品質評
価アプリケーション232と係数格納部235)と、C
PU201とRAM202によりその機能を実現するこ
とができる。
(Image Evaluation Apparatus 100: Image Quality Evaluation Unit 1)
04) The image quality evaluation unit 104 outputs an image quality evaluation value based on the result of the convolution operation. In an example described later, an evaluation value correlated with the subjective evaluation of the image is output using the texture amount and the dot area ratio. Specifically, it outputs an evaluation value D represented by Expression (1). The coefficients p 1 to p 9 used at this time are coefficients determined in advance, and are stored in the coefficient storage unit 235. Image quality evaluation unit 104
Is, for example, the hard disk 203 (especially the image quality evaluation application 232 and the coefficient storage unit 235) and the C
The function can be realized by the PU 201 and the RAM 202.

【0051】(画像評価装置100:網点面積率取得部
105の内容)網点面積率取得部105は、評価すべき
画像の網点面積率を入手し、適宜差の絶対値や平均値を
算出し、画像品質評価部104に出力する。ここで、取
得とは、画像データ入力部101から入力した画像デー
タに基づいて網点面積率を算出する場合と、キーボード
207やテンキーなどを介して数値を入力する場合の何
れであってもよいことを意味する。
(Image Evaluation Apparatus 100: Contents of Dot Area Ratio Acquisition Unit 105) The dot area ratio acquisition unit 105 obtains the dot area ratio of the image to be evaluated and appropriately calculates the absolute value or average value of the difference. Calculated and output to the image quality evaluation unit 104. Here, the acquisition may be either a case where the halftone dot area ratio is calculated based on the image data input from the image data input unit 101 or a case where a numerical value is input via the keyboard 207 or numeric keys. Means that.

【0052】また、差の値や平均値とは、たとえば、A
社のプリンタとB社のプリンタを用いて同一画像を出力
し、比較する場合にその網点面積率の差の値や平均値を
いう。なお、網点面積率取得部105が、差の絶対値や
平均値を算出してもよいし、別途画像品質評価部104
が算出する態様であってもよい。網点面積率取得部10
5は、たとえば、ハードディスク203(特に画像品質
評価アプリケーション232)や、スキャナ210、キ
ーボード207およびCPU201によりその機能を実
現することができる。
The difference value and the average value are, for example, A
When the same image is output using a printer of the company B and a printer of the company B and compared, the difference value or the average value of the halftone dot area ratios is used. Note that the dot area ratio acquisition unit 105 may calculate the absolute value or the average value of the difference, or the image quality evaluation unit 104
May be calculated. Halftone dot area acquisition unit 10
5 can be realized by the hard disk 203 (especially the image quality evaluation application 232), the scanner 210, the keyboard 207, and the CPU 201, for example.

【0053】画像評価装置であるPC100は、以上の
ような構成を有することにより、たとえば特定の空間周
波数帯域に生じるバンディングのような画像欠陥の知覚
される大きさを定量化でき、人間の視覚系の特性に基づ
く画像品質評価をおこなうこと、すなわち、画像に対す
る主観評価と相関のとれた評価値を出力することが可能
となる。
The PC 100, which is an image evaluation apparatus, can quantify the perceived size of an image defect such as banding that occurs in a specific spatial frequency band by having the above-described configuration, and can realize a human visual system. It is possible to perform an image quality evaluation based on the above characteristics, that is, to output an evaluation value correlated with the subjective evaluation of the image.

【0054】(実験例)つぎに、画像評価装置100を
適用した実験例について説明する。ここでは画像の作
成、物理量の算出、主観評価実験、評価値と主観評価実
験との相関の順について説明する。ここでは、簡単のた
めに、中間調処理された画像における階調変化に伴うテ
クスチャ変化による画質劣化を定量化する(評価する)
実験について説明する。
(Experimental Example) Next, an experimental example to which the image evaluation apparatus 100 is applied will be described. Here, the order of creation of an image, calculation of a physical quantity, a subjective evaluation experiment, and a correlation between an evaluation value and a subjective evaluation experiment will be described. Here, for the sake of simplicity, image quality deterioration due to a texture change due to a gradation change in a halftone-processed image is quantified (evaluated).
The experiment will be described.

【0055】(実験例:画像の作成)実験は、テクスチ
ャ構造を有する画像として、300dpiおよび600
dpiでそれぞれ150線でスクリーン角30度の2値
中間調画像の2種類を用いた。150線とは、網点が1
インチの中に150個入ることを示す画像の一属性値で
ある。図3は、実験に用いたスクリーン角30度で15
0線の2値中間調画像(300dpiおよび600dp
i)の拡大図である。
(Experimental Example: Creation of Image) The experiment was conducted at 300 dpi and 600 dpi as an image having a texture structure.
Two types of binary halftone images with a screen angle of 30 degrees and 150 lines each at dpi were used. 150 lines means that the halftone dot is 1
This is one attribute value of an image indicating that 150 pieces are included in an inch. FIG. 3 shows that the screen angle used in the experiment was 30 degrees and the screen angle was 15 degrees.
0 line binary halftone image (300 dpi and 600 dpi)
It is an enlarged view of i).

【0056】物理量算出および主観評価実験に用いる中
間調画像データはコンピュータを用いて作成した。ここ
で、物理量を算出するための画像データとして均一入力
値の矩形画像データを用い、主観評価に用いる画像デー
タとして各々異なる入力値をもち、隣接する2個の矩形
部分からなる画像データを用いた。図4は、物理量を算
出するための画像データとして、1画素8ビットの矩形
の画像で画像全域にわたって値が等しい画像の例を示し
た図である。
The halftone image data used for the physical quantity calculation and the subjective evaluation experiment was created using a computer. Here, rectangular image data having uniform input values were used as image data for calculating physical quantities, and image data having two different input values and image data consisting of two adjacent rectangular portions were used as image data used for subjective evaluation. . FIG. 4 is a diagram showing an example of an image of a rectangular image of 8 bits per pixel and having the same value over the entire image area as image data for calculating a physical quantity.

【0057】図5は、主観評価に用いる画像データとし
て各々異なる入力値をもち、隣接する2個の矩形部分か
らなる画像の例を示した図である。図6は、図5に示し
た画像の拡大図である。なお、図5は見やすさの観点か
ら、矩形部分の左側は160を、右側を128にしてあ
る。
FIG. 5 is a diagram showing an example of an image composed of two adjacent rectangular portions having different input values as image data used for subjective evaluation. FIG. 6 is an enlarged view of the image shown in FIG. In FIG. 5, from the viewpoint of legibility, the left side of the rectangular portion is 160 and the right side is 128.

【0058】これらコンピュータ上で作成された主観評
価用の画像をKodak社のAPPROVALで出力し
た。なお、APPROVALからは図6に示したよう
に、二値化した画像を出力した。
Images for subjective evaluation created on these computers were output by Kodak's APPROVAL. In addition, a binarized image was output from APPROVAL as shown in FIG.

【0059】(実験例:物理量の算出)Marが提唱し
た視覚系のマルチチャンネル理論に基づいて、5種類の
チャンネルを、画質劣化を定量化する基となる物理量の
算出に用いた。Marによれば、受容野(すなわち、あ
る1個の神経節細胞に連結している網膜上の錐体(桿
体)の集合で表される網膜上の領域)の中心部の直径
は、各々、0.022°、0.073°、0.146
°、0.276°、0.495°である(なお、網膜上
の単位面積に結像する見込みの視野角は視距離に関係な
く一定であるので、上記のごとく°で表すことができ
る)。
(Experimental Example: Calculation of Physical Amount) Based on the multi-channel theory of the visual system proposed by Mar, five types of channels were used to calculate a physical amount as a basis for quantifying image quality deterioration. According to Mar, the diameter of the central part of the receptive field (ie, the area on the retina represented by the collection of cones (rods) on the retina connecting to one ganglion cell) is , 0.022 °, 0.073 °, 0.146
°, 0.276 °, and 0.495 ° (note that the viewing angle expected to form an image on the unit area on the retina is constant regardless of the viewing distance, and can be represented by ° as described above). .

【0060】従って、眼球光学系の接点が網膜から17
mmの位置にあり、錐体の直径が2[μm]であるとす
ると、この見込み角度は、錐体直径のそれぞれ1.65
倍、5.42倍、10.84倍、20.46倍、36.
72倍に相当する。また、Marの定義によるチャンネ
ルの数学的記述は、受容野中心からの距離をr(r=s
qrt(x2+y2))とすると、上述した式(4)で与
えられる。
Therefore, the point of contact of the eyeball optical system is 17
mm, and the diameter of the cone is 2 [μm], the expected angles are 1.65 each of the cone diameter.
Times, 5.42 times, 10.84 times, 20.46 times, 36.
Equivalent to 72 times. The mathematical description of the channel according to Mar defines that the distance from the center of the receptive field is r (r = s
qrt (x 2 + y 2 )), which is given by the above equation (4).

【0061】受容野中心部の境界においては、式(4)
の左辺が0となるので(図17参照)、ここで用いる5
種類のチャンネルのσの値は、錐体の直径を1cone
と表すこととすると、1.166cone、3.832
cone、7.665cone、14.464con
e、25.962coneとなる。
At the boundary of the central part of the receptive field, equation (4)
Is 0 on the left side (see FIG. 17).
The value of σ for each type of channel is 1 cone
In this case, 1.166 cone, 3.832
cone, 7.665 cone, 14.464 con
e, 25.962 cone.

【0062】画像の観察距離を350mmとして、解像
度600dpiすなわち、600分の1インチの物体の
像は網膜上では0.0021mmであり、ほぼ錐体の直
径と一致する。従って、600dpiの画像は近似的に
そのまま網膜像として用いることが可能となる。
Assuming that the observation distance of the image is 350 mm, the image of an object having a resolution of 600 dpi, that is, a 1/600 inch is 0.0021 mm on the retina, which is almost equal to the diameter of the cone. Therefore, an image of 600 dpi can be used as a retinal image approximately as it is.

【0063】このようなチャンネルで画像を畳み込む
と、各々のチャンネルのスケールにおける強度変化が検
出される。均一な強度(濃度)の画像に対してはその出
力はすべての画素において0であるが、画像内に強度の
変化があると変化に応じた正または負の値が出力され
る。畳み込みの一例として下式(5)、(6)を用い
て、テクスチャを表現する物理量Tを採用し、後述する
評価値Dを算出することとした。
When an image is convolved with such a channel, a change in intensity on the scale of each channel is detected. For an image with uniform intensity (density), the output is 0 for all pixels, but if there is a change in intensity in the image, a positive or negative value corresponding to the change is output. As an example of convolution, using the following equations (5) and (6), a physical quantity T representing a texture is adopted, and an evaluation value D described later is calculated.

【0064】[0064]

【数5】 (Equation 5)

【0065】但し、式(5)および式(6)の添字iは
各チャンネルに対応しており、Gi(r)は、チャンネ
ルiにおけるG(r)を表すものとする。Nは被演算画
素数を示す。線形和Tを求める際の係数Aiは、1次元
の空間周波数λの正弦波の強度分布をもつ多値画像ti
を求め、様々なλに対するtiの線形和を人間の視覚系
のコントラスト感度空間周波数特性(VTF:Visu
al TransferFunction)に対して回
帰をおこなうことにより決定した。VTFに対して回帰
をおこなったのは、VTFは人間の視覚特性のうち最も
基本的な特性を表現する関数と考えられるからである。
すなわち、Aiは、人間の視覚系の空間周波数特性に最
適化されたされた係数ということができる。
Here, the suffix i in the equations (5) and (6) corresponds to each channel, and G i (r) represents G (r) in the channel i. N indicates the number of pixels to be operated. The coefficient A i for obtaining the linear sum T is a multi-valued image t i having a one-dimensional sine wave intensity distribution with a spatial frequency λ.
Is calculated, and the linear sum of t i for various λ is converted into the contrast sensitivity spatial frequency characteristic (VTF: Visu) of the human visual system.
al Transfer Function). The reason why the regression is performed on the VTF is that the VTF is considered to be a function expressing the most basic characteristic among human visual characteristics.
That, A i can be said to be optimized coefficients in the spatial frequency characteristics of the human visual system.

【0066】図7は、VTFの関数形を示した説明図で
あり、図8は、回帰をおこなうことによって決定された
係数Aiを示した図表である。図8から明らかなよう
に、tiおよびTは画像の強度の変動量に対する増加関
数となっている。また、VTFに対して回帰させたの
で、Aiはσiに一対一に対応する定数であり、評価する
画像の種類によらず使用することができる。
FIG. 7 is an explanatory diagram showing a function form of VTF, and FIG. 8 is a table showing coefficients A i determined by performing regression. As is clear from FIG. 8, t i and T are increasing functions with respect to the variation amount of the intensity of the image. In addition, since regression is performed on the VTF, A i is a constant corresponding to σ i on a one-to-one basis, and can be used regardless of the type of image to be evaluated.

【0067】(実験例:主観評価実験)主観的評価実験
は、入力値の相異なる2個の矩形図形を隣接させたサン
プル群を被験者に提示しおこなった。被験者は、知覚さ
れるテクスチャの差異に応じて、サンプル群を5レベル
に分類した。すなわち、レベル1は「はっきり差が分か
る」、レベル2は「結構差が分かる」、レベル3は、
「少し差があるように見える」、レベル4は「何となく
差があるように見える」、レベル5は「同一に見える」
として、二つの隣接した画像を評価した。図9は、実験
に用いたAPPROVALからの出力例を示した図であ
る。前述したように、主観評価実験では図に示した30
0dpi、150線、スクリーン角30°の二値中間調
画像の何れか二つのサンプルを隣接させておこなった。
(Experimental Example: Subjective Evaluation Experiment) In the subjective evaluation experiment, a sample group in which two rectangular figures having different input values were adjacent was presented to the subject. Subjects classified the sample group into five levels according to the difference in perceived texture. That is, level 1 is "clearly understand the difference", level 2 is "clearly understand the difference", and level 3 is
"Looks a bit different", Level 4 "Looks somewhat different", Level 5 "Looks the same"
As two adjacent images were evaluated. FIG. 9 is a diagram showing an example of output from APPROVAL used in the experiment. As described above, in the subjective evaluation experiment, 30
Any two samples of a binary halftone image at 0 dpi, 150 lines, and a screen angle of 30 ° were placed adjacent to each other.

【0068】(実験例:評価値と主観評価実験との相
関)各サンプルに対して知覚されるテクスチャの差異の
大きさを式(1)で表される関数により与えられる値D
により予測した。上述したように、ΔTは、式(6)で
与えられるテクスチャの差、すなわち、隣接させた片方
のサンプルについてのテクスチャ量T(テクスチャ量T
Lとする)と他方のサンプルについてのテクスチャ量T
(テクスチャ量TRとする)との差(TL−TR)であ
る。また、Tmeanは、Tmean=(TL+TR)/
2である。なお、Tを求める際のR(x,y)として
は、600dpiでスキャンした画像をそのまま用い
た。
(Experimental Example: Correlation between Evaluation Value and Subjective Evaluation Experiment) The magnitude of the difference in texture perceived for each sample is represented by a value D given by a function represented by equation (1).
Predicted by As described above, ΔT is the difference between the textures given by Expression (6), that is, the texture amount T (texture amount T
L ) and the texture amount T for the other sample
Which is the difference between (a texture amount T R) (T L -T R ). Also, Tmean is Tmean = (T L + T R ) /
2. Note that an image scanned at 600 dpi was used as it was as R (x, y) when obtaining T.

【0069】同様に、ΔMは予め分かっているサンプル
の網点面積率の差、すなわち、隣接させた片方のサンプ
ルについての網点面積率M(面積率MLとする)と他方
のサンプルについての網点面積率M(面積率MRとす
る)との差(ML−MR)である。また、Mmeanは、
Mmean=(ML+MR)/2である。
Similarly, ΔM is the difference between the dot area ratios of the samples known in advance, that is, the dot area ratio M (referred to as the area ratio M L ) of one adjacent sample and the dot ratio of the other sample. a halftone dot area rate M (the area ratio M R) the difference between the (M L -M R). Also, M mean
Mmean = a (M L + M R) / 2.

【0070】本願発明者は、テクスチャの量の差(Δ
T)を平均テクスチャ量(Tmean)で補正した量
と、画像の強度(明度)の差(ΔM)を平均強度(明
度)(Mmean)で補正した量の線形和となる形で式
(1)を案出した。各係数p1〜p9については、主観評
価実験の結果に対して回帰をおこない決定した。図10
は、最適化により求められた係数p1〜p9の値を示した
図表である。
The inventor of the present application has found that the difference in the amount of texture (Δ
Equation (1) in the form of a linear sum of an amount obtained by correcting T) with the average texture amount (Tmean) and an amount obtained by correcting the difference (ΔM) between the intensity (brightness) of the image with the average intensity (brightness) (Mmean). Devised. Each of the coefficients p 1 to p 9 was determined by performing a regression on the result of the subjective evaluation experiment. FIG.
Is a table showing the values of the coefficients p 1 ~p 9 obtained by optimization.

【0071】係数p3が負であるため式(1)の第1項
は平均テクスチャ量(Tmean)に対する増加関数で
あり、画像の強度変動が大きいほどテクスチャを知覚し
やすいという現象に対応しているものと考えられる。ま
た、明度は網点面積率に対して減少関数であるが、係数
p7が負であるため第2項は平均強度(明度)に対する
増加関数であり、画像が明るいほど明るさの差を知覚し
やすいという現象に対応しているものと考えられる。従
って、式(1)により算出される評価値Dは、実際の心
理物理的現象と矛盾しない予測式ないし評価式であると
いえる。
Since the coefficient p3 is negative, the first term of the equation (1) is an increasing function with respect to the average texture amount (Tmean), and corresponds to the phenomenon that the greater the fluctuation of the image intensity, the easier the texture is perceived. It is considered something. The brightness is a decreasing function with respect to the halftone dot area ratio, but since the coefficient p7 is negative, the second term is an increasing function with respect to the average intensity (brightness). It is thought that this corresponds to the phenomenon of being easy. Therefore, it can be said that the evaluation value D calculated by Expression (1) is a prediction expression or evaluation expression that does not contradict the actual psychophysical phenomenon.

【0072】(画像評価装置100の処理流れ)最後
に、画像評価装置100の処理流れについて説明する。
図11は画像評価装置100の処理流れの一例を示した
フローチャートである。画像の品質の評価をおこなうた
めに、まず、評価すべき画像の画像データを入力する
(ステップS1101)。この入力は、スキャナ210
(図2参照)から入力してもよいが、電子的に作成した
理想的なデータを入力してもよい。
(Processing Flow of Image Evaluation Apparatus 100) Finally, the processing flow of the image evaluation apparatus 100 will be described.
FIG. 11 is a flowchart illustrating an example of the processing flow of the image evaluation device 100. In order to evaluate the quality of an image, first, image data of an image to be evaluated is input (step S1101). This input is input to the scanner 210
(See FIG. 2), or ideal electronically created data may be input.

【0073】つぎに、ステップS1101で入力された
画像データに対して使用する空間演算フィルタを選択す
る(ステップS1102)。空間演算フィルタは、式
(4)に示したMarのモデルのような等方性微分演算
フィルタが挙げられ、使用するσi(i=channe
l1〜channel5)の全部もしくはいくつかを選
択する態様であってもよいし、また、式(4)とは異な
る複数の空間演算フィルタから使用すべき空間演算フィ
ルタを選択してもよい。
Next, a spatial operation filter to be used for the image data input in step S1101 is selected (step S1102). Examples of the spatial operation filter include an isotropic differential operation filter such as the Mar model shown in Expression (4), and σ i (i = channel
11 to 5) may be selected, or a spatial operation filter to be used may be selected from a plurality of spatial operation filters different from Expression (4).

【0074】つぎに、ステップS1101で入力された
画像データに対してステップS1102で選択された空
間演算フィルタを用いて畳み込み演算をおこなう(ステ
ップS1103)。具体的な一例として、式(5)およ
び式(6)により表される演算をおこなうことが挙げら
れる。つづいて、ステップS1103でおこなわれた畳
み込み演算の結果に基づいてステップS1101で入力
された画像データを有する画像の品質を評価する(ステ
ップS1104)。具体的な一例として、式(1)に表
される評価式に基づいて評価値Dを求めることが挙げら
れる。すなわち、ステップS1103で求めたテクスチ
ャ量と、別途取得した網点面積率とに基づいて評価値D
を求める。
Next, a convolution operation is performed on the image data input in step S1101 using the spatial operation filter selected in step S1102 (step S1103). As a specific example, performing an operation represented by Expressions (5) and (6) can be cited. Subsequently, the quality of the image having the image data input in step S1101 is evaluated based on the result of the convolution operation performed in step S1103 (step S1104). As a specific example, obtaining the evaluation value D based on the evaluation expression represented by Expression (1) is given. That is, the evaluation value D is calculated based on the texture amount obtained in step S1103 and the dot area ratio separately acquired.
Ask for.

【0075】以上説明した様に実施の形態1の画像評価
装置は、人間の視覚系の特性に基づく画像品質評価をお
こなうことができ、これにより、画像に対する主観評価
と相関のとれた評価値を出力することが可能な画像評価
装置を提供することができる。特に、空間演算フィルタ
としてマルチチャンネルモデルに基づいた等方性微分演
算フィルタΔGを用い、各チャンネルからの出力値ti
に対する係数Aiも、VTFに回帰させて得た値なの
で、人間の視覚特性を反映させた評価値を得ることがで
きる。
As described above, the image evaluation apparatus according to the first embodiment can perform image quality evaluation based on the characteristics of the human visual system, and can thereby obtain an evaluation value correlated with the subjective evaluation of an image. An image evaluation device capable of outputting can be provided. In particular, an isotropic differential operation filter ΔG based on a multi-channel model is used as a spatial operation filter, and output values t i from each channel are used.
Factor A i for even, so values obtained by regression to VTF, it is possible to obtain an evaluation value that reflects the human visual characteristics.

【0076】実施の形態2.実施の形態2では、評価値
Dを補正してより評価の精度を向上させる、すなわち、
人間主観評価と高精度に相関のとれた評価値を出力する
画像評価装置について説明する。なお、実施の形態2で
は、実施の形態1と同様の構成部分については同一の符
号を付し、その説明を省略するものとする。
Embodiment 2 In the second embodiment, the evaluation value D is corrected to further improve the accuracy of the evaluation.
An image evaluation device that outputs an evaluation value highly accurately correlated with human subjective evaluation will be described. In the second embodiment, the same components as those in the first embodiment are denoted by the same reference numerals, and description thereof will be omitted.

【0077】図12は、実施の形態2の画像評価装置を
PCに適用した場合の機能ブロック図である。画像評価
装置1200は、実施の形態1の画像評価装置100
に、さらに、画像評価値補正部1201を含ませたもの
である。つぎに、画像評価値補正部1201の内容につ
いて説明する。
FIG. 12 is a functional block diagram when the image evaluation device according to the second embodiment is applied to a PC. The image evaluation device 1200 is an image evaluation device 100 according to the first embodiment.
Further, an image evaluation value correction unit 1201 is included. Next, the contents of the image evaluation value correction unit 1201 will be described.

【0078】(画像評価装置1200:画像評価値補正
部1201の内容)画像評価値補正部1201は、画像
品質評価部104から出力された評価値Dを式(2)に
基づいて画像評価値PDに補正する。本願発明者は、ロ
ジスティック回帰をおこなうため式(2)の関数形を採
用し、ロジスティック回帰をおこなって係数q1〜q4
決定した。図13は、ロジスティック回帰により求めら
れた係数q1〜q4の値を示した図表である。
(Image Evaluation Device 1200: Contents of Image Evaluation Value Correction Unit 1201) The image evaluation value correction unit 1201 converts the evaluation value D output from the image quality evaluation unit 104 into an image evaluation value PD based on the equation (2). To be corrected. The inventor of the present application adopted the function form of Expression (2) to perform logistic regression, and determined coefficients q 1 to q 4 by performing logistic regression. Figure 13 is a table showing the value of the coefficient q 1 to q 4 obtained by logistic regression.

【0079】図14は、式(2)による主観評価結果の
予測精度、すなわち、補正された評価値の予測精度を示
した説明図である。図に示したように、予測精度は相関
係数0.901(寄与率0.811)であり、画像評価
装置1200は、画像に対する主観評価と非常に高い相
関のとれた評価値を出力することが可能であることが分
かった。画像評価値補正部1201は、たとえば、ハー
ドディスク203(特に画像品質評価アプリケーション
232、係数格納部235)、CPU201およびRA
M202によりその機能を実現することができる。
FIG. 14 is an explanatory diagram showing the prediction accuracy of the subjective evaluation result by equation (2), that is, the prediction accuracy of the corrected evaluation value. As shown in the figure, the prediction accuracy is a correlation coefficient of 0.901 (contribution ratio of 0.811), and the image evaluation device 1200 outputs an evaluation value that is highly correlated with the subjective evaluation of the image. Turned out to be possible. The image evaluation value correction unit 1201 includes, for example, the hard disk 203 (in particular, the image quality evaluation application 232, the coefficient storage unit 235), the CPU 201, and the RA.
The function can be realized by M202.

【0080】(画像評価装置1200の処理流れ)最後
に、画像評価装置100の処理流れについて説明する。
図15は画像評価装置1200の処理流れの一例を示し
たフローチャートである。画像の品質の評価をおこなう
ために、まず、評価すべき画像の画像データを入力する
(ステップS1501)。この入力は、スキャナ210
(図2参照)から入力してもよいが、電子的に作成した
理想的なデータを入力してもよい。
(Processing Flow of Image Evaluation Apparatus 1200) Finally, the processing flow of the image evaluation apparatus 100 will be described.
FIG. 15 is a flowchart illustrating an example of the processing flow of the image evaluation device 1200. In order to evaluate the quality of an image, first, image data of an image to be evaluated is input (step S1501). This input is input to the scanner 210
(See FIG. 2), or ideal electronically created data may be input.

【0081】つぎに、ステップS1501で入力された
画像データに対して使用する空間演算フィルタを選択す
る(ステップS1502)。空間演算フィルタは、式
(4)に示したMarのモデルのように等方性微分演算
フィルタが挙げられ、使用するσi(i=channe
l1〜channel5)の全部もしくはいくつかを選
択する態様であってもよいし、また、式(4)とは異な
る複数の空間演算フィルタから使用すべき空間演算フィ
ルタを選択してもよい。
Next, a spatial operation filter to be used for the image data input in step S1501 is selected (step S1502). An example of the spatial operation filter is an isotropic differential operation filter as in the Mar model shown in Expression (4), and σ i (i = channel
11 to 5) may be selected, or a spatial operation filter to be used may be selected from a plurality of spatial operation filters different from Expression (4).

【0082】つぎに、ステップS1501で入力された
画像データに対してステップS1502で選択された空
間演算フィルタを用いて畳み込み演算をおこなう(ステ
ップS1503)。具体的な一例として、式(5)およ
び式(6)により表される演算をおこなうことが挙げら
れる。つづいて、ステップS1503でおこなわれた畳
み込み演算の結果に基づいてステップS1501で入力
された画像データを有する画像の品質を評価する(ステ
ップS1504)。具体的な一例として、式(1)に表
される評価式に基づいて評価値Dを求めることが挙げら
れる。すなわち、ステップS1503で求めたテクスチ
ャ量と、別途取得した網点面積率とに基づいて評価値D
を求める。
Next, a convolution operation is performed on the image data input in step S1501 using the spatial operation filter selected in step S1502 (step S1503). As a specific example, performing an operation represented by Expressions (5) and (6) can be cited. Subsequently, the quality of the image having the image data input in step S1501 is evaluated based on the result of the convolution operation performed in step S1503 (step S1504). As a specific example, obtaining the evaluation value D based on the evaluation expression represented by Expression (1) is given. That is, the evaluation value D is calculated based on the texture amount obtained in step S1503 and the halftone dot area ratio acquired separately.
Ask for.

【0083】画像評価装置1200は、ステップS15
04で算出された評価値Dを式(2)に基づいて評価値
PDに補正する(ステップS1505)。補正する際の
係数q1〜q4は係数格納部235に格納しておく。
The image evaluation apparatus 1200 proceeds to step S15
The evaluation value D calculated in step 04 is corrected to an evaluation value PD based on equation (2) (step S1505). The coefficients q 1 to q 4 for correction are stored in the coefficient storage unit 235.

【0084】以上説明したように、実施の形態2の画像
評価装置は、画像に対する主観評価と高精度に相関のと
れた評価値を出力することが可能となる。たとえば、従
来では人間には差が知覚されるが評価値ではうまく差が
表れないような、テクスチャのみによるエッジの検出を
相関よく評価することが可能となる。
As described above, the image evaluation device according to the second embodiment can output an evaluation value highly correlated with the subjective evaluation of an image. For example, in the related art, it is possible to evaluate the detection of an edge using only a texture with good correlation so that a difference is perceived by a human but the difference is not well represented by an evaluation value.

【0085】実施の形態3.実施の形態3では、入力さ
れる画像データの密度を調節する画像評価装置について
説明する。なお、実施の形態3では、実施の形態1と同
様の構成部分については同一の符合を付し、その説明を
省略するものとする。
Embodiment 3 Embodiment 3 describes an image evaluation device that adjusts the density of input image data. In the third embodiment, the same components as those in the first embodiment are denoted by the same reference numerals, and description thereof will be omitted.

【0086】図16は、実施の形態3の画像評価装置を
PCに適用した場合の機能ブロック図である。画像評価
装置1600は、実施の形態1の画像評価装置100
に、さらに、入力密度調節部1601を含ませたもので
ある。つぎに、入力密度調節部1601の内容について
説明する。
FIG. 16 is a functional block diagram when the image evaluation apparatus according to the third embodiment is applied to a PC. The image evaluation device 1600 is an image evaluation device 100 according to the first embodiment.
In addition, an input density adjusting unit 1601 is further included. Next, the contents of the input density adjustment unit 1601 will be described.

【0087】(画像評価装置1600:入力密度調節部
1601の内容)入力密度調節部1601は、画像デー
タ入力部101により入力される画像データの入力密度
を、人間が当該画像を観察する際の網膜上の結像密度と
略等しくなるように調節する。調節は以下の方法により
おこなう。まず、人間が被評価画像を観察する際の観察
距離を入力する。この入力された観察距離から、人間が
被評価画像を観察する際の網膜上の結像密度を演算す
る。つぎに、演算された網膜上の結像密度と近似的に等
しくなるように、画像入力密度を調節する。
(Image Evaluation Apparatus 1600: Contents of Input Density Adjusting Unit 1601) The input density adjusting unit 1601 determines the input density of the image data input by the image data input unit 101 by using the retina when a human observes the image. Adjust so as to be approximately equal to the above imaging density. The adjustment is performed by the following method. First, an observation distance when a human observes the evaluated image is input. From the input observation distance, the imaging density on the retina when a human observes the evaluated image is calculated. Next, the image input density is adjusted so as to be approximately equal to the calculated image density on the retina.

【0088】たとえば、人間が被評価画像を観察する際
の観察距離をD[mm]とすると、人間の眼球光学系の
接点は通常網膜から17[mm]離れた位置にあり、且
つ、網膜上錐体の直径は略0.002[mm]であると
いわれているので、錐体1個に結像する物体の直径I
[mm]は、I=D・0.002/17で表される。従
って、入力密度調節部1601は、画像入力密度が1画
素当たりI[mm]となるように調整する。
For example, assuming that the observation distance when a human observes an image to be evaluated is D [mm], the contact point of the human eye optical system is usually located at a distance of 17 [mm] from the retina and is located on the retina. Since the diameter of the cone is said to be approximately 0.002 [mm], the diameter I of the object imaged on one cone is
[Mm] is represented by I = D · 0.002 / 17. Therefore, the input density adjusting unit 1601 adjusts the image input density to be I [mm] per pixel.

【0089】この様に、入力密度を調整することによ
り、実際の人間の視覚系の特性に基づく画像品質評価を
おこなうことが可能となる。すなわち、被評価画像を人
間が観察する際の観察距離が異なると、網膜上の結像密
度も異なってしまうが、入力密度調節部1601は、様
々な相異なる観察距離にて人間が被評価画像を観察する
場合に対応でき、相関の高い評価値を出力することが可
能となる。なお、仕様の態様によっては、各σを同一比
をもって拡大もしくは縮小調節してもよい。
As described above, by adjusting the input density, it is possible to evaluate the image quality based on the characteristics of the actual human visual system. That is, if the observation distance at which the evaluated image is observed by a human is different, the image formation density on the retina is also different. However, the input density adjusting unit 1601 determines that the human can be evaluated at various different observation distances. Can be observed, and an evaluation value having a high correlation can be output. Note that, depending on the mode of the specification, each σ may be adjusted to be enlarged or reduced with the same ratio.

【0090】以上説明したように、本発明は、心理物理
的側面、神経生理学的側面、解剖学的側面からもその存
在が確実視されている視覚系のチャンネルについてのモ
デルを用いている。これにより、実在する視覚メカニズ
ムに基づく演繹法的手法を画質評価に導入することがで
き、画像に対する主観評価と相関のとれた評価値を出力
することが可能となる。すなわち、本発明により、本質
的に網膜上の強度変化を検出する機構であるチャネルを
用いて、粒状性、鮮鋭性、バンディングといった画像の
濃度変動に起因する画質劣化の要因の定量化をおこなう
ことが可能となり、画像の品質の評価値として出力する
ことが可能となった。
As described above, the present invention uses a model of a visual system channel whose existence is surely recognized from the psychophysical, neurophysiological, and anatomical aspects. As a result, a deductive method based on a real visual mechanism can be introduced into image quality evaluation, and an evaluation value correlated with the subjective evaluation of an image can be output. In other words, according to the present invention, using a channel, which is a mechanism for detecting a change in intensity on the retina, quantifying factors such as graininess, sharpness, and banding that cause image quality deterioration due to image density fluctuation. Can be output as an evaluation value of image quality.

【0091】また、この評価装置、評価手法を用いて、
各種の画像出力機器(プリンタ、スキャナ、CRTモニ
タ)を評価することができる。例えば、A社とB社のプ
リンタが同じ2400dpiとしていても、同じ画像デ
ータからプリントアウトされた画像は微妙に感じが異な
っていることがある。この場合、本発明の画像評価装
置、画像評価方法を用いることにより、どちらが「よ
い」プリンタであるかという、人間の感覚的な尺度に沿
った評価を定量的におこうことが可能となる。このよう
な、いわば量的変数(画像データ)から質的変数(画像
がよい悪いの人間の判断)を予測可能とするのは、本装
置もしくは本方法が、ロジスティック回帰分析を利用し
ているからであるといえる。
Also, using this evaluation device and evaluation method,
Various image output devices (printer, scanner, CRT monitor) can be evaluated. For example, even if the printers of Company A and Company B use the same 2400 dpi, images printed out from the same image data may have slightly different feelings. In this case, by using the image evaluation device and the image evaluation method of the present invention, it is possible to quantitatively evaluate which is a “good” printer in accordance with a human sensory scale. The reason that such a qualitative variable (human judgment of good or bad image) can be predicted from such a quantitative variable (image data) is that the present apparatus or the present method uses logistic regression analysis. You can say that.

【0092】なお、本実施の形態で説明した画像処理方
法は、あらかじめ用意されたプログラムをパーソナル・
コンピューターやワークステーション等のコンピュータ
で実行することにより実現することができる。このプロ
グラムは、ハードディスク、フロッピー(登録商標)デ
ィスク、CD−ROM、MO、DVD等のコンピュータ
で読み取り可能な記録媒体に記録され、コンピュータに
よって記録媒体から読み出されることによって実行され
る。またこのプログラムは、上記記録媒体を介して、イ
ンターネット等のネットワークを介して配布することが
できる。
The image processing method described in the present embodiment uses a program prepared in advance for a personal computer.
It can be realized by executing on a computer such as a computer or a workstation. This program is recorded on a computer-readable recording medium such as a hard disk, a floppy (registered trademark) disk, a CD-ROM, an MO, and a DVD, and is executed by being read from the recording medium by the computer. This program can be distributed via the recording medium and a network such as the Internet.

【0093】[0093]

【発明の効果】以上説明したように、本発明の画像評価
装置(請求項1)は、画像入力手段が画像の空間的光学
的情報を入力し、畳込演算手段が、画像入力手段により
入力された画像の空間的光学的情報に対する空間演算フ
ィルタの畳み込み演算をおこない、画質予測手段が、畳
込演算手段の出力を基に、人間が知覚する画像品質を予
測するので、人間の視覚系の特性に基づく画像品質評価
をおこなうことができ、これにより、画像に対する主観
評価と相関のとれた評価値を出力することが可能な画像
評価装置を提供することができる。
As described above, according to the image evaluation apparatus of the present invention (claim 1), the image input means inputs the spatial and optical information of the image, and the convolution operation means inputs the image by the image input means. The spatial image processing unit performs a convolution operation on the spatial and optical information of the obtained image, and the image quality prediction unit predicts the image quality perceived by a human based on the output of the convolution operation unit. It is possible to perform image quality evaluation based on characteristics, and thereby to provide an image evaluation device capable of outputting an evaluation value correlated with a subjective evaluation of an image.

【0094】また、本発明の画像評価装置(請求項2)
は、請求項1記載の画像評価装置において、空間演算フ
ィルタが等方性微分演算フィルタであるので、2次元の
画像構造に起因する画質要因を人間の視覚系の特性に基
づいて定量化することができ、これにより、画像に対す
る主観評価と相関のとれた評価値を出力することが可能
な画像評価装置を提供することができる。
Further, the image evaluation device of the present invention (Claim 2)
In the image evaluation device according to claim 1, since the spatial operation filter is an isotropic differential operation filter, an image quality factor caused by a two-dimensional image structure is quantified based on characteristics of a human visual system. Accordingly, it is possible to provide an image evaluation device capable of outputting an evaluation value correlated with the subjective evaluation of an image.

【0095】また、本発明の画像評価装置(請求項3)
は、請求項1または2に記載の画像評価装置において、
空間演算フィルタが各々特性の相異なる複数のフィルタ
からなり、画質予測手段が複数のフィルタにそれぞれに
対応した畳み込み演算手段からの出力すべてを用いて、
人間が知覚する画像品質を予測するので、様々な尺度の
画像構造に起因する画質要因を人間の視覚系の特性に基
づいて定量化し、画質を総合的に定量化することがで
き、これにより、画像に対する主観評価と相関のとれた
評価値を出力することが可能な画像評価装置を提供する
ことができる。
Further, the image evaluation apparatus of the present invention (Claim 3)
Is the image evaluation device according to claim 1 or 2,
The spatial calculation filter includes a plurality of filters having different characteristics, and the image quality prediction unit uses all outputs from the convolution calculation units corresponding to the plurality of filters, respectively.
Since it predicts the image quality perceived by humans, it is possible to quantify the image quality factors caused by the image structure of various scales based on the characteristics of the human visual system and comprehensively quantify the image quality. An image evaluation device capable of outputting an evaluation value correlated with a subjective evaluation of an image can be provided.

【0096】また、本発明の画像評価装置(請求項4)
は、請求項3に記載の画像評価装置において、複数のフ
ィルタが人間の視覚系の空間周波数特性に対して最適化
されているので、人間の視覚系の空間周波数特性に基づ
く、より精緻な画質の定量化をおこなうことができ、こ
れにより、画像に対する主観評価と相関のとれた評価値
を出力することが可能な画像評価装置を提供することが
できる。
Further, the image evaluation device of the present invention (Claim 4)
In the image evaluation device according to claim 3, since the plurality of filters are optimized for the spatial frequency characteristics of the human visual system, a more precise image quality based on the spatial frequency characteristics of the human visual system Therefore, it is possible to provide an image evaluation apparatus capable of outputting an evaluation value correlated with the subjective evaluation of an image.

【0097】また、本発明の画像評価装置(請求項5)
は、請求項1〜4の何れか一つに記載の画像評価装置に
おいて、空間演算フィルタとして各々特性の相異なる5
個のフィルタを具備したので、人間の視覚系のマルチチ
ャンネルモデルに基づく、より精緻な画質の定量化をお
こなうことができ、これにより、画像に対する主観評価
と相関のとれた評価値を出力することが可能な画像評価
装置を提供することができる。
Further, the image evaluation device of the present invention (Claim 5)
Is the image evaluation device according to any one of claims 1 to 4, wherein the spatial operation filters have different characteristics.
With the number of filters, it is possible to perform more precise quantification of image quality based on a multi-channel model of the human visual system, thereby outputting an evaluation value correlated with the subjective evaluation of the image. It is possible to provide an image evaluation device capable of performing the following.

【0098】また、本発明の画像評価装置(請求項6)
は、請求項3、4または5に記載の画像評価装置におい
て、フィルタ選択手段が、複数のフィルタのうち必要な
数個を選択するので、画像構造が各々相異なる視覚的尺
度に分布するような複数種の画像欠陥に対応した知覚を
人間の視覚系のマルチチャンネルモデルに基づいて、よ
り高速かつ精緻に定量化することができ、これにより、
画像に対する主観評価と相関のとれた評価値を出力する
ことが可能な画像評価装置を提供することができる。
Further, the image evaluation apparatus of the present invention (claim 6)
In the image evaluation device according to the third, fourth or fifth aspect, since the filter selecting means selects a required number of the plurality of filters, the image structure is distributed on different visual scales. Perceptions corresponding to multiple types of image defects can be quantified faster and more precisely based on a multi-channel model of the human visual system.
An image evaluation device capable of outputting an evaluation value correlated with a subjective evaluation of an image can be provided.

【0099】また、本発明の画像評価装置(請求項7)
は、請求項1〜6の何れか一つに記載の画像評価装置に
おいて、フィルタサイズ変更手段が、空間演算フィルタ
のサイズを各々可変とするので、被評価画像の画像密度
に依存しない画質評価をおこなうことができ、これによ
り、画像に対する主観評価と相関のとれた評価値を出力
することが可能な画像評価装置を提供することができ
る。
Further, the image evaluation device of the present invention (claim 7)
In the image evaluation apparatus according to any one of claims 1 to 6, the filter size changing unit changes the size of the spatial operation filter, so that the image quality evaluation independent of the image density of the image to be evaluated can be performed. This makes it possible to provide an image evaluation device capable of outputting an evaluation value correlated with the subjective evaluation of an image.

【0100】また、本発明の画像評価装置(請求項8)
は、画像データ入力手段が、評価すべき画像の画像デー
タを入力し、空間演算フィルタ選択手段が、画像データ
入力手段により入力された画像データに対して使用する
空間演算フィルタを選択し、畳込演算手段が、画像デー
タ入力手段により入力された画像データに対して空間演
算フィルタ選択手段により選択された空間演算フィルタ
を用いて畳み込み演算をおこない、画像品質評価手段
が、畳込演算手段によりおこなわれた畳み込み演算の結
果に基づいて画像の品質を評価するので、人間の視覚系
の特性に基づく画像品質評価をおこなうことができ、こ
れにより、画像に対する主観評価と相関のとれた評価値
を出力することが可能な画像評価装置を提供することが
できる。
An image evaluation apparatus according to the present invention (claim 8)
The image data input means inputs image data of an image to be evaluated, and the spatial calculation filter selection means selects a spatial calculation filter to be used for the image data input by the image data input means, and performs convolution. The arithmetic means performs a convolution operation on the image data input by the image data input means using the spatial operation filter selected by the spatial operation filter selection means, and the image quality evaluation means is performed by the convolution operation means. Since the quality of the image is evaluated based on the result of the convolution operation, the image quality can be evaluated based on the characteristics of the human visual system, thereby outputting an evaluation value correlated with the subjective evaluation of the image. It is possible to provide an image evaluation device capable of performing the above.

【0101】また、本発明の画像評価装置(請求項9)
は、請求項8に記載の画像評価装置において、空間演算
フィルタを等方性微分演算フィルタとしたので、2次元
の画像構造に起因する画質要因を人間の視覚系の特性に
基づいて定量化することができ、これにより、画像に対
する主観評価と相関のとれた評価値を出力することが可
能な画像評価装置を提供することができる。
An image evaluation apparatus according to the present invention (claim 9)
In the image evaluation device according to the eighth aspect, since the spatial operation filter is an isotropic differential operation filter, the image quality factor due to the two-dimensional image structure is quantified based on the characteristics of the human visual system. Accordingly, it is possible to provide an image evaluation device capable of outputting an evaluation value correlated with the subjective evaluation of an image.

【0102】また、本発明の画像評価装置(請求項1
0)は、請求項9に記載の画像評価装置において、等方
性微分演算フィルタを、関数LG(r)=(−1/(π
σ4))・(1−r2/(2σ2))・exp(−r2
(2σ2))(但し、rは、r=sqrt(x2+y2
で与えられる距離であり、x,yは評価すべき画像の画
像データの或る位置を中心とした座標系における座標で
あり、σは視覚系のマルチチャンネルモデルでいう5種
類のチャンネルに対応した定数である)としたので、人
間の視覚系のマルチチャンネルモデルに基づく、より精
細な画質の定量化をおこなうことができ、これにより、
画像に対する主観評価と相関のとれた評価値を出力する
ことが可能な画像評価装置を提供することができる。
The image evaluation device of the present invention (claim 1)
0) is the image evaluation device according to claim 9, wherein the isotropic differential operation filter is set to a function LG (r) = (− 1 / (π)
σ 4 )) · (1−r 2 / (2σ 2 )) · exp (−r 2 /
(2σ 2 )) (where r is r = sqrt (x 2 + y 2 )
X and y are coordinates in a coordinate system centered on a certain position of the image data of the image to be evaluated, and σ corresponds to five types of channels in a multi-channel model of a visual system. Is a constant), which allows for more detailed image quality quantification based on a multi-channel model of the human visual system,
An image evaluation device capable of outputting an evaluation value correlated with the subjective evaluation of an image can be provided.

【0103】また、本発明の画像評価装置(請求項1
1)は、二つの画像を比較して同じに見えるか差がある
ように見えるかの主観的な評価に沿った評価を出力する
画像評価装置であって、画像データ入力手段が、評価す
べき二つの画像の画像データを入力し、テクスチャ量算
出手段が、画像データ入力手段により入力されたデータ
に基づいて、テクスチャ量の差の絶対値ΔTとテクスチ
ャ量の平均値Tmeanとを算出し、網点面積率取得手
段が、評価すべき二つの画像の網点面積率の差の絶対値
ΔMと網点面積率の平均値Mmeanとを取得し、画像
評価値算出手段が、テクスチャ量算出手段により算出さ
れたテクスチャ量の差の絶対値ΔTとテクスチャ量の平
均値Tmean、および、網点面積率取得手段により取
得された網点面積率の差の絶対値ΔMと網点面積率の平
均値Mmeanとを用いて、式(1)により画像評価値
Dを算出するので、巨視的には同じ濃度である画像のテ
クスチャによる見え方の差を算出することができ、これ
により、画像に対する主観評価と相関のとれた評価値を
出力することが可能な画像評価装置を提供することがで
きる。
The image evaluation apparatus of the present invention (claim 1)
1) is an image evaluation device that outputs an evaluation according to a subjective evaluation of whether two images look the same or there is a difference by comparing the two images, wherein the image data input means should evaluate The image data of the two images is input, and the texture amount calculation means calculates the absolute value ΔT of the difference between the texture amounts and the average value Tmean of the texture amounts based on the data input by the image data input means. The dot area ratio obtaining means obtains the absolute value ΔM of the difference between the halftone dot area ratios of the two images to be evaluated and the average value Mmean of the halftone dot area ratios. The calculated absolute value ΔT of the difference in the texture amount and the average value Tmean of the texture amount, and the absolute value ΔM of the difference in the halftone dot area ratio and the average value Mmean of the halftone dot area ratio obtained by the halftone dot area ratio obtaining means. And Then, since the image evaluation value D is calculated by the equation (1), it is possible to calculate a difference in appearance due to the texture of an image having the same density macroscopically, thereby obtaining a correlation with the subjective evaluation of the image. And an image evaluation device capable of outputting the evaluated value.

【0104】また、本発明の画像評価装置(請求項1
2)は、請求項11に記載の画像評価装置において、さ
らに、画像評価値補正手段が、画像評価値算出手段によ
り算出された画像評価値Dを、式(2)に基づいて画像
評価値PDに補正するので、画像に対する主観評価と高
精度に相関のとれた評価値を出力することが可能な画像
評価装置を提供することができる。
The image evaluation device of the present invention (claim 1)
2) The image evaluation device according to claim 11, wherein the image evaluation value correction unit further converts the image evaluation value D calculated by the image evaluation value calculation unit into an image evaluation value PD based on Expression (2). Therefore, it is possible to provide an image evaluation apparatus capable of outputting an evaluation value highly accurately correlated with the subjective evaluation of an image.

【0105】また、本発明の画像評価装置(請求項1
3)は、請求項1〜12に記載の画像評価装置におい
て、入力密度調節手段が、画像データ入力手段により入
力される画像データの入力密度を、人間が当該画像を観
察する際の網膜上の結像密度と略等しくなるように調整
するので、被評価画像の画像密度と網膜上の画像密度を
略等しくすることができ、これにより、画像に対する主
観評価と相関のとれた評価値を出力することが可能な画
像評価装置を提供することができる。
The image evaluation device of the present invention (claim 1)
3) The image evaluation device according to claim 1, wherein the input density adjusting unit determines the input density of the image data input by the image data input unit on the retina when a human observes the image. Since the image density is adjusted so as to be substantially equal to the image density, the image density of the image to be evaluated and the image density on the retina can be made approximately equal, thereby outputting an evaluation value correlated with the subjective evaluation of the image. It is possible to provide an image evaluation device capable of performing the above.

【0106】また、本発明の画像評価方法(請求項1
4)は、画像の空間的光学的情報に対して空間演算フィ
ルタを畳み込んで得られる出力を基に、人間が知覚する
画像品質を予測するので、人間の視覚系の特性に基づく
画像品質評価をおこなうことができ、これにより、画像
に対する主観評価と相関のとれた評価値を出力すること
が可能な画像評価方法を提供することができる。
The image evaluation method of the present invention (claim 1)
4) predicts the image quality perceived by a human based on the output obtained by convolving the spatial and optical information of the image with the spatial operation filter, and thus evaluates the image quality based on the characteristics of the human visual system. Accordingly, it is possible to provide an image evaluation method capable of outputting an evaluation value correlated with the subjective evaluation of an image.

【0107】また、本発明の画像評価方法(請求項1
5)は、請求項14に記載の画像評価方法において、空
間演算フィルタが等方性微分演算フィルタであるので、
2次元の画像構造に起因する画質要因を人間の視覚系の
特性に基づいて定量化でき、これにより、画像に対する
主観評価と相関のとれた評価値を出力することが可能な
画像評価方法を提供することができる。
The image evaluation method of the present invention (claim 1)
5) In the image evaluation method according to claim 14, the spatial operation filter is an isotropic differential operation filter.
Provided is an image evaluation method capable of quantifying an image quality factor caused by a two-dimensional image structure based on characteristics of a human visual system, thereby outputting an evaluation value correlated with a subjective evaluation of an image. can do.

【0108】また、本発明の画像評価方法(請求項1
6)は、請求項14または15に記載の画像評価方法に
おいて、空間演算フィルタとして各々特性の相異なる複
数の空間演算フィルタ群を用い、当該空間演算フィルタ
群により畳み込んで得られる各出力のすべてを用いて、
人間が知覚する画像品質を予測するので、様々な尺度の
画像構造に起因する画質要因を人間の視覚系の特性に基
づいて定量化し、画質を総合的に定量化することがで
き、これにより、画像に対する主観評価と相関のとれた
評価値を出力することが可能な画像評価方法を提供する
ことができる。
The image evaluation method of the present invention (claim 1)
6) In the image evaluation method according to claim 14 or 15, all of the outputs obtained by convolving with the spatial operation filter group using a plurality of spatial operation filter groups each having a different characteristic as the spatial operation filter. Using,
Since it predicts the image quality perceived by humans, it is possible to quantify the image quality factors caused by the image structure of various scales based on the characteristics of the human visual system and comprehensively quantify the image quality. An image evaluation method capable of outputting an evaluation value correlated with a subjective evaluation of an image can be provided.

【0109】また、本発明の画像評価方法(請求項1
7)は、請求項16に記載の画像評価方法において、空
間演算フィルタ群が人間の視覚系の空間周波数特性に対
して最適化されているので、人間の視覚系の空間周波数
特性に基づくより精緻な画質の定量化をおこなうことが
でき、これにより、画像に対する主観評価と相関のとれ
た評価値を出力することが可能な画像評価方法を提供す
ることができる。
The image evaluation method of the present invention (claim 1)
7) In the image evaluation method according to claim 16, since the spatial operation filter group is optimized for the spatial frequency characteristics of the human visual system, more refinement based on the spatial frequency characteristics of the human visual system. Thus, it is possible to provide an image evaluation method capable of outputting an evaluation value correlated with a subjective evaluation of an image.

【0110】また、本発明の画像評価方法(請求項1
8)は、請求項14〜17の何れか一つに記載の画像評
価方法において、空間演算フィルタとして各々特性の相
異なる5個の空間演算フィルタを用いるので、人間の視
覚系のマルチチャンネルモデルに基づくより精緻な画質
の定量化をおこなうことができ、これにより、画像に対
する主観評価と相関のとれた評価値を出力することが可
能な画像評価方法を提供することができる。
The image evaluation method of the present invention (claim 1)
8) In the image evaluation method according to any one of claims 14 to 17, since five spatial operation filters having different characteristics are used as the spatial operation filters, a multi-channel model of a human visual system is used. It is possible to perform a more precise quantification of the image quality based on the image evaluation, thereby providing an image evaluation method capable of outputting an evaluation value correlated with the subjective evaluation of the image.

【0111】また、本発明の画像評価方法(請求項1
9)は、請求項16、17または18に記載の画像評価
方法において、空間演算フィルタのうち必要な数個を選
択して使用するので、画像構造が各々相異なる視覚的尺
度に分布するような複数種の画像欠陥に対応した知覚を
人間の視覚系のマルチチャンネルモデルに基づいて、よ
り高速かつ精緻に定量化することができ、これにより、
画像に対する主観評価と相関のとれた評価値を出力する
ことが可能な画像評価方法を提供することができる。
The image evaluation method of the present invention (claim 1)
In 9), in the image evaluation method according to claim 16, 17 or 18, a necessary number of spatial operation filters are selected and used, so that the image structure is distributed on different visual scales. Perceptions corresponding to multiple types of image defects can be quantified faster and more precisely based on a multi-channel model of the human visual system.
An image evaluation method capable of outputting an evaluation value correlated with a subjective evaluation of an image can be provided.

【0112】また、本発明の画像評価方法(請求項2
0)は、請求項14〜19の何れか一つに記載の画像評
価方法において、空間演算フィルタのサイズが各々可変
であるので、被評価画像の画像密度に依存しない画質評
価をおこなうことができ、これにより、画像に対する主
観評価と相関のとれた評価値を出力することが可能な画
像評価方法を提供することができる。
The image evaluation method of the present invention (claim 2)
0) In the image evaluation method according to any one of claims 14 to 19, since the size of the spatial operation filter is variable, image quality evaluation independent of the image density of the evaluated image can be performed. Accordingly, it is possible to provide an image evaluation method capable of outputting an evaluation value correlated with the subjective evaluation of an image.

【0113】また、本発明のプログラム(請求項21)
は、コンピュータに画像評価を実行させるプログラムで
あって、コンピュータを、前記請求項1〜請求項13の
何れか一つに記載の各手段として機能させるので、請求
項1〜請求項13の何れか一つに記載の各手段をコンピ
ュータに実行させることができる。
Further, the program of the present invention (claim 21)
Is a program that causes a computer to execute image evaluation, and causes the computer to function as each unit according to any one of claims 1 to 13. Each means described in one can be executed by a computer.

【図面の簡単な説明】[Brief description of the drawings]

【図1】実施の形態1の画像評価装置をPCに適用した
場合の機能ブロック図である。
FIG. 1 is a functional block diagram when an image evaluation device according to a first embodiment is applied to a PC.

【図2】実施の形態1の画像評価装置のハードウェア構
成の一例を示した説明図である。
FIG. 2 is an explanatory diagram illustrating an example of a hardware configuration of the image evaluation device according to the first embodiment;

【図3】実験に用いたスクリーン角30度で150線の
2値中間調画像(300dpiおよび600dpi)の
拡大図である。
FIG. 3 is an enlarged view of a binary halftone image (300 dpi and 600 dpi) with a screen angle of 30 degrees and 150 lines used in the experiment.

【図4】物理量を算出するための画像データとして、1
画素8ビットの矩形の画像で画像全域にわたって値が等
しい画像の例を示した図である。
FIG. 4 shows 1 as image data for calculating a physical quantity.
FIG. 7 is a diagram illustrating an example of an image of a rectangular image having 8 bits of pixels and having the same value over the entire image.

【図5】主観評価に用いる画像データとして各々異なる
入力値をもち、隣接する2個の矩形部分からなる画像の
例を示した図である。
FIG. 5 is a diagram illustrating an example of an image including two adjacent rectangular portions having different input values as image data used for subjective evaluation.

【図6】図5に示した画像の拡大図である。FIG. 6 is an enlarged view of the image shown in FIG. 5;

【図7】VTFの関数形を示した説明図である。FIG. 7 is an explanatory diagram showing a function form of a VTF.

【図8】実施の形態1の画像評価装置に用いるテクスチ
ャ量Tを算出する際に使用する係数であって、回帰をお
こなうことによって決定された係数Aiを示した図表で
ある。
FIG. 8 is a table showing coefficients A i used in calculating the texture amount T used in the image evaluation device according to the first embodiment, which coefficients are determined by performing regression.

【図9】実施の形態1で説明した主観評価実験に用いた
APPROVALからの出力例を示した図である。
FIG. 9 is a diagram showing an output example from APPROVAL used in the subjective evaluation experiment described in the first embodiment.

【図10】実施の形態1の画像評価装置に用いる評価値
Dを算出する際に使用する係数であって、最適化をする
ことにより求められた係数p1〜p9の値を示した図表で
ある。
FIG. 10 is a chart showing coefficients p 1 to p 9 which are coefficients used when calculating an evaluation value D used in the image evaluation apparatus according to the first embodiment and which are obtained by optimization; It is.

【図11】実施の形態1の画像評価装置の処理流れの一
例を示したフローチャートである。
FIG. 11 is a flowchart illustrating an example of a processing flow of the image evaluation device according to the first embodiment;

【図12】実施の形態2の画像評価装置をPCに適用し
た場合の機能ブロック図である。
FIG. 12 is a functional block diagram when the image evaluation apparatus according to the second embodiment is applied to a PC.

【図13】実施の形態2の画像評価装置に用いる評価値
PDを算出する際に使用する係数であって、ロジスティ
ック回帰により求められた係数q1〜q4の値を示した図
表である。
FIG. 13 is a table showing coefficients q 1 to q 4 obtained by logistic regression, which are coefficients used when calculating an evaluation value PD used in the image evaluation apparatus according to the second embodiment.

【図14】実施の形態2の画像評価装置で補正された評
価値の予測精度を示した説明図である。
FIG. 14 is an explanatory diagram showing prediction accuracy of an evaluation value corrected by the image evaluation device according to the second embodiment.

【図15】実施の形態2の画像評価装置の処理流れの一
例を示したフローチャートである。
FIG. 15 is a flowchart illustrating an example of a processing flow of the image evaluation device according to the second embodiment;

【図16】実施の形態3の画像評価装置をPCに適用し
た場合の機能ブロック図である。
FIG. 16 is a functional block diagram when the image evaluation device according to the third embodiment is applied to a PC.

【図17】チャンネルの受容野の概略形状を示した説明
図であり、受容野平面上に各受光体からの出力をz軸方
向にプロットしたものを水平方向から表示した図であ
る。
FIG. 17 is an explanatory diagram showing a schematic shape of a receptive field of a channel, and is a diagram in which outputs from respective photoreceptors are plotted in a z-axis direction on a receptive field plane and displayed in a horizontal direction.

【符号の説明】[Explanation of symbols]

100 画像評価装置 101 画像データ入力部 102 演算フィルタ選択部 103 畳込演算部 104 画像品質評価部 105 網点面積率取得部 106 演算フィルタ格納部 203 ハードディスク 206 プリンタ 210 スキャナ 232 画像品質評価アプリケーション 233 画像データ格納部 234 フィルタ格納部 235 係数格納部 1200 画像評価装置 1201 画像評価値補正部 1600 画像評価装置 1601 入力密度調節部 D 画像評価値 M 網点面積率 PD 補正後の評価値 T テクスチャ量 ΔG 等方性微分演算フィルタ REFERENCE SIGNS LIST 100 image evaluation device 101 image data input unit 102 operation filter selection unit 103 convolution operation unit 104 image quality evaluation unit 105 halftone dot area ratio acquisition unit 106 operation filter storage unit 203 hard disk 206 printer 210 scanner 232 image quality evaluation application 233 image data Storage unit 234 Filter storage unit 235 Coefficient storage unit 1200 Image evaluation device 1201 Image evaluation value correction unit 1600 Image evaluation device 1601 Input density adjustment unit D Image evaluation value M Halftone dot area ratio PD Evaluation value after correction T Texture amount ΔG Isotropic Sex differential operation filter

Claims (21)

【特許請求の範囲】[Claims] 【請求項1】 画像の空間的光学的情報を入力する画像
入力手段と、 前記画像入力手段により入力された空間的光学的情報に
対して使用する空間演算フィルタと、 前記画像入力手段により入力された画像の空間的光学的
情報に対する前記空間演算フィルタの畳み込み演算をお
こなう畳込演算手段と、 前記畳込演算手段の出力を基に、人間が知覚する画像品
質を予測する画質予測手段と、 を具備したことを特徴とする画像評価装置。
1. An image input unit for inputting spatial and optical information of an image, a spatial operation filter used for the spatial and optical information input by the image input unit, and an input by the image input unit Convolution operation means for performing a convolution operation of the spatial operation filter on the spatial and optical information of the image, and image quality prediction means for predicting image quality perceived by humans based on the output of the convolution operation means. An image evaluation device, comprising:
【請求項2】 前記空間演算フィルタが等方性微分演算
フィルタであることを特徴とする請求項1記載の画像評
価装置。
2. The image evaluation device according to claim 1, wherein the spatial operation filter is an isotropic differential operation filter.
【請求項3】 前記空間演算フィルタは各々特性の相異
なる複数のフィルタからなり、前記画質予測手段は前記
複数のフィルタにそれぞれに対応した前記畳み込み演算
手段からの出力すべてを用いて、人間が知覚する画像品
質を予測することを特徴とする請求項1または2に記載
の画像評価装置。
3. The spatial operation filter comprises a plurality of filters having different characteristics, and the image quality prediction means uses all the outputs from the convolution operation means respectively corresponding to the plurality of filters, so that a human can perceive. The image evaluation device according to claim 1, wherein an image quality to be predicted is predicted.
【請求項4】 前記複数のフィルタが人間の視覚系の空
間周波数特性に対して最適化されていることを特徴とす
る請求項3に記載の画像評価装置。
4. The apparatus according to claim 3, wherein the plurality of filters are optimized for spatial frequency characteristics of a human visual system.
【請求項5】 前記空間演算フィルタとして各々特性の
相異なる5個のフィルタを具備したことを特徴とする請
求項1〜4の何れか一つに記載の画像評価装置。
5. The image evaluation device according to claim 1, further comprising five filters having different characteristics, respectively, as said spatial operation filter.
【請求項6】 前記複数のフィルタのうち必要な数個を
選択するフィルタ選択手段を具備したことを特徴とする
請求項3、4または5に記載の画像評価装置。
6. The image evaluation apparatus according to claim 3, further comprising a filter selection unit that selects a required number of the plurality of filters.
【請求項7】 前記空間演算フィルタのサイズを各々可
変とするフィルタサイズ変更手段を具備したことを特徴
とする請求項1〜6の何れか一つに記載の画像評価装
置。
7. The image evaluation apparatus according to claim 1, further comprising: a filter size changing unit that changes a size of the spatial operation filter.
【請求項8】 評価すべき画像の画像データを入力する
画像データ入力手段と、 前記画像データ入力手段により入力された画像データに
対して使用する空間演算フィルタを選択する空間演算フ
ィルタ選択手段と、 前記画像データ入力手段により入力された画像データに
対して前記空間演算フィルタ選択手段により選択された
空間演算フィルタを用いて畳み込み演算をおこなう畳込
演算手段と、 前記畳込演算手段によりおこなわれた畳み込み演算の結
果に基づいて前記画像の品質を評価する画像品質評価手
段と、 を具備したことを特徴とする画像評価装置。
8. An image data input unit for inputting image data of an image to be evaluated, a spatial operation filter selecting unit for selecting a spatial operation filter to be used for the image data input by the image data input unit, Convolution operation means for performing a convolution operation on the image data input by the image data input means using the spatial operation filter selected by the spatial operation filter selection means; convolution operation performed by the convolution operation means An image quality evaluation device, comprising: image quality evaluation means for evaluating the quality of the image based on a result of the calculation.
【請求項9】 前記空間演算フィルタを等方性微分演算
フィルタとしたことを特徴とする請求項8記載の画像評
価装置。
9. The image evaluation device according to claim 8, wherein the spatial operation filter is an isotropic differential operation filter.
【請求項10】 前記等方性微分演算フィルタを、関数
LG(r)=(−1/(πσ4))・(1−r2/(2σ
2))・exp(−r2/(2σ2))(但し、rは、r
=sqrt(x2+y2)で与えられる距離であり、x,
yは前記評価すべき画像の画像データの或る位置を中心
とした座標系における座標であり、σは視覚系のマルチ
チャンネルモデルでいう5種類のチャンネルに対応した
定数である)としたことを特徴とする請求項9に記載の
画像評価装置。
10. The isotropic differential operation filter is defined by a function LG (r) = (− 1 / (πσ 4 )) · (1−r 2 / (2σ)
2 )). Exp (−r 2 / (2σ 2 )) (where r is r
= Sqrt (x 2 + y 2 ), where x,
y is a coordinate in a coordinate system centered on a certain position of the image data of the image to be evaluated, and σ is a constant corresponding to five types of channels in a multi-channel model of a visual system.) The image evaluation device according to claim 9, wherein:
【請求項11】 二つの画像を比較して同じに見えるか
差があるように見えるかの主観的な評価に沿った評価を
出力する画像評価装置であって、 評価すべき二つの画像の画像データを入力する画像デー
タ入力手段と、 前記画像データ入力手段により入力されたデータに基づ
いて、テクスチャ量の差の絶対値ΔTとテクスチャ量の
平均値Tmeanとを算出するテクスチャ量算出手段
と、 前記評価すべき二つの画像の網点面積率の差の絶対値Δ
Mと網点面積率の平均値Mmeanとを取得する網点面
積率取得手段と、 前記テクスチャ量算出手段により算出されたテクスチャ
量の差の絶対値ΔTとテクスチャ量の平均値Tmea
n、および、前記網点面積率取得手段により取得された
網点面積率の差の絶対値ΔMと網点面積率の平均値Mm
eanとを用いて、下式(1)により画像評価値Dを算
出する画像評価値算出手段と、 を具備したことを特徴とする画像評価装置。 【数1】 (但しp1〜p9は、予めおこなわれた主観評価実験によ
り決定された定数とする。)
11. An image evaluation device for comparing two images and outputting an evaluation in accordance with a subjective evaluation of whether the two images look the same or different from each other, comprising: Image data input means for inputting data; texture amount calculation means for calculating an absolute value ΔT of a difference between texture amounts and an average value Tmean of texture amounts based on the data input by the image data input means; Absolute value Δ of the difference between the dot area ratios of two images to be evaluated
A dot area ratio acquiring unit for acquiring M and an average value Mmean of the dot area ratio; an absolute value ΔT of a difference between the texture amounts calculated by the texture amount calculating unit; and an average value Tmea of the texture amounts
n, the absolute value ΔM of the difference between the halftone dot area ratios obtained by the halftone dot area ratio obtaining means and the average value Mm of the halftone dot area ratios
An image evaluation device comprising: an image evaluation value calculating unit that calculates an image evaluation value D according to the following equation (1) using E.an. (Equation 1) (However, p 1 to p 9 are constants determined by a subjective evaluation experiment performed in advance.)
【請求項12】 さらに、前記画像評価値算出手段によ
り算出された画像評価値Dを、下式(2)に基づいて画
像評価値PDに補正する画像評価値補正手段を具備した
ことを特徴とする請求項11に記載の画像評価装置。 【数2】 (但しq1〜q4は、回帰分析により決定された定数と
する。)
12. An image evaluation value correction means for correcting the image evaluation value D calculated by the image evaluation value calculation means to an image evaluation value PD based on the following equation (2). The image evaluation device according to claim 11, wherein: (Equation 2) (However, q1 to q4 are constants determined by regression analysis.)
【請求項13】 前記画像データ入力手段により入力さ
れる画像データの入力密度を、人間が当該画像を観察す
る際の網膜上の結像密度と略等しくなるように調整する
入力密度調節手段を備えたことを特徴とする請求項1〜
12に記載の画像評価装置。
13. An input density adjusting means for adjusting an input density of image data input by said image data input means so as to be substantially equal to an image formation density on a retina when a human observes the image. Claim 1 characterized by the above-mentioned.
13. The image evaluation device according to 12.
【請求項14】 画像の空間的光学的情報に対して空間
演算フィルタを畳み込んで得られる出力を基に、人間が
知覚する画像品質を予測することを特徴とする画像評価
方法。
14. An image evaluation method, wherein an image quality perceived by a human is predicted based on an output obtained by convolving a spatial operation filter with spatial and optical information of an image.
【請求項15】 前記空間演算フィルタが等方性微分演
算フィルタであることを特徴とする請求項14に記載の
画像評価方法。
15. The image evaluation method according to claim 14, wherein the spatial operation filter is an isotropic differential operation filter.
【請求項16】 空間演算フィルタとして各々特性の相
異なる複数の空間演算フィルタ群を用い、当該空間演算
フィルタ群により畳み込んで得られる各出力のすべてを
用いて、人間が知覚する画像品質を予測することを特徴
とする請求項14または15に記載の画像評価方法。
16. A plurality of spatial operation filters having different characteristics are used as spatial operation filters, and image quality perceived by humans is predicted using all outputs obtained by convolution by the spatial operation filters. The image evaluation method according to claim 14, wherein the evaluation is performed.
【請求項17】 前記空間演算フィルタ群が人間の視覚
系の空間周波数特性に対して最適化されていることを特
徴とする請求項16に記載の画像評価方法。
17. The image evaluation method according to claim 16, wherein the spatial operation filter group is optimized for a spatial frequency characteristic of a human visual system.
【請求項18】 前記空間演算フィルタとして各々特性
の相異なる5個の空間演算フィルタを用いることを特徴
とする請求項14〜17の何れか一つに記載の画像評価
方法。
18. The image evaluation method according to claim 14, wherein five spatial operation filters having different characteristics are used as the spatial operation filters.
【請求項19】 前記空間演算フィルタのうち必要な数
個を選択して使用することを特徴とする請求項16、1
7または18に記載の画像評価方法。
19. The apparatus according to claim 16, wherein a necessary number of the spatial operation filters are selected and used.
19. The image evaluation method according to 7 or 18.
【請求項20】 前記空間演算フィルタのサイズが各々
可変であることを特徴とする請求項14〜19の何れか
一つに記載の画像評価方法。
20. The image evaluation method according to claim 14, wherein the size of each of the spatial operation filters is variable.
【請求項21】 コンピュータに画像評価を実行させる
プログラムであって、 コンピュータを、 前記請求項1〜請求項13の何れか一つに記載の各手段
として機能させることを特徴とするプログラム。
21. A program for causing a computer to execute image evaluation, wherein the program causes a computer to function as each unit according to any one of claims 1 to 13.
JP2001040855A 2001-02-16 2001-02-16 Image evaluation device Expired - Fee Related JP4040259B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2001040855A JP4040259B2 (en) 2001-02-16 2001-02-16 Image evaluation device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2001040855A JP4040259B2 (en) 2001-02-16 2001-02-16 Image evaluation device

Publications (2)

Publication Number Publication Date
JP2002245464A true JP2002245464A (en) 2002-08-30
JP4040259B2 JP4040259B2 (en) 2008-01-30

Family

ID=18903377

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2001040855A Expired - Fee Related JP4040259B2 (en) 2001-02-16 2001-02-16 Image evaluation device

Country Status (1)

Country Link
JP (1) JP4040259B2 (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007043655A (en) * 2005-06-29 2007-02-15 Ntt Docomo Inc Video evaluation apparatus, method, and program, and spatio-temporal frequency analyzing apparatus, method, and program
JP2009037634A (en) * 2001-11-16 2009-02-19 Mitsutoyo Corp Operation method of machine vision system and case-specific boundary positioning system
JP2009199611A (en) * 2006-04-26 2009-09-03 Aware Inc Fingerprint preview quality and segmentation
JP2010055404A (en) * 2008-08-28 2010-03-11 Canon Inc Image processing method and image processing apparatus
JP2010134704A (en) * 2008-12-04 2010-06-17 Canon Inc Image processing apparatus, image processing method, and computer program
WO2016175923A1 (en) * 2015-04-28 2016-11-03 Qualcomm Incorporated Filter specificity as training criterion for neural networks
JP2018160863A (en) * 2017-03-23 2018-10-11 京セラドキュメントソリューションズ株式会社 Dither matrix creation method, dither matrix creation device, image processing apparatus, and dither matrix creation program
CN112990156A (en) * 2021-05-12 2021-06-18 深圳市安软科技股份有限公司 Optimal target capturing method and device based on video and related equipment
JP7428249B2 (en) 2020-06-02 2024-02-06 コニカミノルタ株式会社 Image quality evaluation device, image forming device, and image quality evaluation method

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102637301B (en) * 2012-03-20 2014-08-13 武汉大学 Method for automatically evaluating color quality of image during aerial photography in real time

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0415882A (en) * 1990-05-09 1992-01-21 Matsushita Electric Ind Co Ltd Picture quality evaluation device
JPH0696218A (en) * 1992-06-30 1994-04-08 Eastman Kodak Co Method and device for determining difference between images, which can be visually recognized
JPH06168315A (en) * 1992-06-15 1994-06-14 Dainippon Printing Co Ltd Image inspection device
JPH0720042A (en) * 1993-07-07 1995-01-24 Kubota Corp Device for spectroscopic analysis

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0415882A (en) * 1990-05-09 1992-01-21 Matsushita Electric Ind Co Ltd Picture quality evaluation device
JPH06168315A (en) * 1992-06-15 1994-06-14 Dainippon Printing Co Ltd Image inspection device
JPH0696218A (en) * 1992-06-30 1994-04-08 Eastman Kodak Co Method and device for determining difference between images, which can be visually recognized
JPH0720042A (en) * 1993-07-07 1995-01-24 Kubota Corp Device for spectroscopic analysis

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4727703B2 (en) * 2001-11-16 2011-07-20 株式会社ミツトヨ Method of operating machine vision system and boundary positioning system
JP2009037634A (en) * 2001-11-16 2009-02-19 Mitsutoyo Corp Operation method of machine vision system and case-specific boundary positioning system
JP2007043655A (en) * 2005-06-29 2007-02-15 Ntt Docomo Inc Video evaluation apparatus, method, and program, and spatio-temporal frequency analyzing apparatus, method, and program
US8452060B2 (en) 2006-04-26 2013-05-28 Aware, Inc. Fingerprint preview quality and segmentation
US9152843B2 (en) 2006-04-26 2015-10-06 Aware, Inc. Fingerprint preview quality and segmentation
US7936907B2 (en) 2006-04-26 2011-05-03 Aware, Inc. Fingerprint preview quality and segmentation
US10325137B2 (en) 2006-04-26 2019-06-18 Aware, Inc. Fingerprint preview quality and segmentation
US8238621B2 (en) 2006-04-26 2012-08-07 Aware, Inc. Fingerprint preview quality and segmentation
US10776604B2 (en) 2006-04-26 2020-09-15 Aware, Inc. Fingerprint preview quality and segmentation
US9031291B2 (en) 2006-04-26 2015-05-12 Aware, Inc. Fingerprint preview quality and segmentation
US11250239B2 (en) 2006-04-26 2022-02-15 Aware, Inc. Fingerprint preview quality and segmentation
US9405957B2 (en) 2006-04-26 2016-08-02 Aware, Inc. Fingerprint preview quality and segmentation
JP2009199611A (en) * 2006-04-26 2009-09-03 Aware Inc Fingerprint preview quality and segmentation
US9626548B2 (en) 2006-04-26 2017-04-18 Aware, Inc. Fingerprint preview quality and segmentation
US9792483B2 (en) 2006-04-26 2017-10-17 Aware, Inc. Fingerprint preview quality and segmentation
US10083339B2 (en) 2006-04-26 2018-09-25 Aware, Inc. Fingerprint preview quality and segmentation
JP2010055404A (en) * 2008-08-28 2010-03-11 Canon Inc Image processing method and image processing apparatus
JP2010134704A (en) * 2008-12-04 2010-06-17 Canon Inc Image processing apparatus, image processing method, and computer program
WO2016175923A1 (en) * 2015-04-28 2016-11-03 Qualcomm Incorporated Filter specificity as training criterion for neural networks
US10515304B2 (en) 2015-04-28 2019-12-24 Qualcomm Incorporated Filter specificity as training criterion for neural networks
JP2018160863A (en) * 2017-03-23 2018-10-11 京セラドキュメントソリューションズ株式会社 Dither matrix creation method, dither matrix creation device, image processing apparatus, and dither matrix creation program
JP7428249B2 (en) 2020-06-02 2024-02-06 コニカミノルタ株式会社 Image quality evaluation device, image forming device, and image quality evaluation method
CN112990156A (en) * 2021-05-12 2021-06-18 深圳市安软科技股份有限公司 Optimal target capturing method and device based on video and related equipment
CN112990156B (en) * 2021-05-12 2021-08-13 深圳市安软科技股份有限公司 Optimal target capturing method and device based on video and related equipment

Also Published As

Publication number Publication date
JP4040259B2 (en) 2008-01-30

Similar Documents

Publication Publication Date Title
Zheng et al. A new metric based on extended spatial frequency and its application to DWT based fusion algorithms
Mantiuk et al. HDR-VDP-2: A calibrated visual metric for visibility and quality predictions in all luminance conditions
Chen et al. A human perception inspired quality metric for image fusion based on regional information
US7747099B2 (en) System and method for comparing images with different contrast levels
JP4959760B2 (en) Method for evaluating noise appearance in images
Peli Contrast sensitivity function and image discrimination
Zhu et al. A no-reference objective image quality metric based on perceptually weighted local noise
TW200831887A (en) Image processing apparatus, data processing apparatus, parameter adjusting method and record medium
EP2884457B1 (en) Image evaluation device and image evaluation program
JP2002245464A (en) Device and method for image evaluation, and program making computer evaluate image
Liu et al. Free energy adjusted peak signal to noise ratio (FEA-PSNR) for image quality assessment
JP2017060737A (en) Processing for creating transmission image without artificial noises
JP4383352B2 (en) Histological evaluation of nuclear polymorphism
Zhang et al. A perceptually relevant channelized joint observer (PCJO) for the detection-localization of parametric signals
CN110889820A (en) Image quality assessment
Renuka et al. An objective measure for assessing the quality of contrast enhancement on magnetic resonance images
Alaguselvi et al. Quantitative analysis of fundus image enhancement in the detection of diabetic retinopathy using deep convolutional neural network
Qureshi et al. A comprehensive performance evaluation of objective quality metrics for contrast enhancement techniques
US20240104724A1 (en) Radiomics standardization
Nagpal et al. Performance analysis of diabetic retinopathy using diverse image enhancement techniques
CN114862762A (en) Quality evaluation method and device for human body scanning image
de Arruda et al. Multiple sclerosis segmentation method in magnetic resonance imaging using fuzzy connectedness, binarization, mathematical morphology, and 3D reconstruction
Farrell et al. Optimizing subpixel rendering using a perceptual metric
Huberman et al. Reducing lateral visual biases in displays
Jigo et al. How voluntary and involuntary attention differentially shape spatial resolution

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20050128

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20070801

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20070807

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20071009

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20071106

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20071107

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101116

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111116

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111116

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121116

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20131116

Year of fee payment: 6

LAPS Cancellation because of no payment of annual fees