JP2010034959A - Image processing apparatus and method, and image display device - Google Patents

Image processing apparatus and method, and image display device Download PDF

Info

Publication number
JP2010034959A
JP2010034959A JP2008196277A JP2008196277A JP2010034959A JP 2010034959 A JP2010034959 A JP 2010034959A JP 2008196277 A JP2008196277 A JP 2008196277A JP 2008196277 A JP2008196277 A JP 2008196277A JP 2010034959 A JP2010034959 A JP 2010034959A
Authority
JP
Japan
Prior art keywords
image
frequency component
intermediate image
processing
horizontal
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
JP2008196277A
Other languages
Japanese (ja)
Other versions
JP4610642B2 (en
JP2010034959A5 (en
Inventor
Shotaro Moriya
正太郎 守谷
Noritaka Okuda
悟崇 奥田
Satoshi Yamanaka
聡 山中
Koji Minami
浩次 南
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.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
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 Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP2008196277A priority Critical patent/JP4610642B2/en
Priority to EP09802651.1A priority patent/EP2320637B1/en
Priority to CN200980130164.3A priority patent/CN102113307B/en
Priority to US13/056,540 priority patent/US8498499B2/en
Priority to PCT/JP2009/003330 priority patent/WO2010013401A1/en
Publication of JP2010034959A publication Critical patent/JP2010034959A/en
Publication of JP2010034959A5 publication Critical patent/JP2010034959A5/ja
Application granted granted Critical
Publication of JP4610642B2 publication Critical patent/JP4610642B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

<P>PROBLEM TO BE SOLVED: To solve the problem that, in an image processing apparatus that properly sets an enhancement coefficient to a detailed portion image of a desired frequency band, enhancement processing may be made inappropriate or insufficient in accordance with an input image and an output image of proper image quality can not be obtained. <P>SOLUTION: An image processing apparatus includes: an intermediate image generating means (1) for generating a first intermediate image (D1) obtained by extracting a component in the vicinity of a specific frequency band from an input image (DIN); an intermediate image processing means (2) for generating a second intermediate image (D2) obtained by processing the first intermediate image (D1); and an adding means (3) for adding the input image (DIN) and the second intermediate image (D2). The intermediate image processing means (2) includes a nonlinear processing means (2A) for changing contents of processing in accordance with pixels of the intermediate image (D1). <P>COPYRIGHT: (C)2010,JPO&INPIT

Description

本発明は、入力画像に対し強調処理する画像処理装置及び方法、並びにこれらを用いた画像表示装置に関し、例えば入力画像として、元となる画像を拡大した拡大画像が入力された際に、高周波数成分の生成および加算をすることによって、解像感の高い出力画像を得るよう画像の強調処理を行うものである。   The present invention relates to an image processing apparatus and method for emphasizing an input image, and an image display apparatus using the same. For example, when an enlarged image obtained by enlarging an original image is input as an input image, a high frequency Image enhancement processing is performed so as to obtain an output image with high resolution by generating and adding components.

一般に画像を表す画像信号に対し適宜画像処理を施した後、画像を再生表示するということが行われている。   In general, an image signal representing an image is appropriately subjected to image processing and then reproduced and displayed.

例えば特許文献1に記載された画像処理装置においては、多重解像度に変換された細部画像に対して、所望の周波数帯域の細部画像に対する強調係数をその所望の周波数帯域よりも低周波数帯域の細部画像の信号に基づいて設定することにより、所望の周波数帯域を強調している。   For example, in the image processing apparatus described in Patent Document 1, with respect to a detail image converted to multi-resolution, an enhancement coefficient for a detail image in a desired frequency band is set to a detail image in a lower frequency band than the desired frequency band. The desired frequency band is emphasized by setting based on the above signal.

特開平9−44651号公報JP-A-9-44651

しかしながら、多重解像度に変換された細部画像に対して、所望の周波数帯域の細部画像に対する強調係数を適宜設定する画像処理装置では、入力画像によっては強調処理が不適切あるいは不十分となり、適正な画質の出力画像を得ることが不可能であった。   However, in an image processing apparatus that appropriately sets an enhancement coefficient for a detail image in a desired frequency band with respect to a detail image converted to multi-resolution, the enhancement processing is inappropriate or insufficient depending on the input image, and the appropriate image quality It was impossible to obtain the output image.

例えば、入力画像として拡大処理を受けた画像が入力される場合、入力画像の周波数スペクトルの高周波数成分側には、拡大処理前の画像の周波数スペクトルの一部が折り返した成分(折り返し成分)が現れる。したがって単純に高周波数成分を強調すると、この折り返し成分を強調してしまい、不適切な処理となる。また、周波数帯域を限定し、折り返し成分を含まない周波数帯域のみを強調すると、周波数スペクトルで考えた場合、高周波数成分側の強調を避けることになり、結果的に不十分な強調処理となってしまう。   For example, when an image that has undergone enlargement processing is input as an input image, a component (folding component) in which a part of the frequency spectrum of the image before enlargement processing is folded is present on the high frequency component side of the frequency spectrum of the input image. appear. Therefore, if the high frequency component is simply emphasized, this aliasing component is emphasized, which is inappropriate processing. In addition, if the frequency band is limited and only the frequency band that does not include the aliasing component is emphasized, when considering the frequency spectrum, emphasis on the high frequency component side is avoided, resulting in insufficient enhancement processing. End up.

また、入力画像としてノイズ処理を受けた画像が入力される場合、高周波数成分側の周波数スペクトルはノイズ処理によって失われている。したがって高周波数成分を取り出そうとしても、取り出すことができず、十分に画像の強調処理を行えないことがある。   When an image subjected to noise processing is input as an input image, the frequency spectrum on the high frequency component side is lost due to noise processing. Therefore, even if an attempt is made to extract a high frequency component, it cannot be extracted, and image enhancement processing may not be performed sufficiently.

本発明の画像処理装置は、
入力画像の特定の周波数帯域近傍の成分を取り出した第1の中間画像を生成する中間画像生成手段と、
前記第1の中間画像をもとに第2の中間画像を生成する中間画像処理手段と、
前記入力画像と前記第2の中間画像を加算する第1の加算手段とを有する画像処理装置において、
前記中間画像処理手段は、前記第1の中間画像の画素に応じて処理の内容を変化させる非線形処理手段を含む
ことを特徴とする。
The image processing apparatus of the present invention
Intermediate image generating means for generating a first intermediate image obtained by extracting a component in the vicinity of a specific frequency band of the input image;
Intermediate image processing means for generating a second intermediate image based on the first intermediate image;
In an image processing apparatus having a first addition means for adding the input image and the second intermediate image,
The intermediate image processing means includes non-linear processing means for changing processing contents in accordance with pixels of the first intermediate image.

本発明によれば、入力画像がその周波数スペクトルにおいて、高周波数成分側に折り返し成分を含んでいる場合や、高周波数成分を十分に含んでいない場合でも十分に画像の強調処理を行うことができる。   According to the present invention, even when the input image includes a aliasing component on the high frequency component side in the frequency spectrum or when the input image does not sufficiently include the high frequency component, the image can be sufficiently enhanced. .

実施の形態1.
図1は本発明の実施の形態1による画像処理装置の構成を表す図であり、例えば画像表示装置の一部として用いることができる。
Embodiment 1 FIG.
FIG. 1 is a diagram showing a configuration of an image processing apparatus according to Embodiment 1 of the present invention, and can be used as a part of an image display apparatus, for example.

図示の画像処理装置は、中間画像生成手段1と、中間画像処理手段2と、加算手段3とを有する。
中間画像生成手段1は、入力画像DINの特定の周波数帯域近傍の成分を取り出した中間画像(第1の中間画像)D2を生成する。
中間画像処理手段2は、第1の中間画像D2に後述する処理を行なった中間画像(第2の中間画像)D2を生成する。
加算手段3は、入力画像DINと第2の中間画像D2を加算する。加算手段3による加算結果が出力画像DOUTとして出力される。
The illustrated image processing apparatus includes an intermediate image generation unit 1, an intermediate image processing unit 2, and an addition unit 3.
The intermediate image generation means 1 generates an intermediate image (first intermediate image) D2 obtained by extracting components near a specific frequency band of the input image DIN.
The intermediate image processing means 2 generates an intermediate image (second intermediate image) D2 obtained by performing processing described later on the first intermediate image D2.
The adding means 3 adds the input image DIN and the second intermediate image D2. The addition result by the adding means 3 is output as an output image DOUT.

中間画像生成手段1は、入力画像DINから高周波数成分のみを取り出した画像(第1の高周波数成分画像)D1Aを生成する高周波数成分画像生成手段1Aと、画像D1Aの低周波数成分のみを取り出した画像D1Bを生成する低周波数成分画像生成手段1Bとを備える。高周波数成分画像生成手段1Aと低周波数成分画像生成手段1Bとで帯域通過フィルタ手段が構成されている。中間画像生成手段1からは画像D1Bが中間画像D1として出力される。   The intermediate image generation unit 1 extracts a high frequency component image generation unit 1A that generates an image (first high frequency component image) D1A obtained by extracting only a high frequency component from the input image DIN, and extracts only a low frequency component of the image D1A. And a low frequency component image generating means 1B for generating the image D1B. The high-frequency component image generating means 1A and the low-frequency component image generating means 1B constitute a band pass filter means. From the intermediate image generating means 1, the image D1B is output as the intermediate image D1.

中間画像処理手段2は、中間画像D1に対し、後述する非線形処理を行った画像D2Aを出力する非線形処理手段2Aと、画像D2Aの高周波数成分のみを取り出した画像D2Bを出力する高周波数成分画像生成手段2Bと、中間画像D1と画像D2Bを加算した画像D2Cを出力する加算手段2Cとを備える。中間画像処理手段2からは画像D2Cが中間画像D2として出力される。   The intermediate image processing means 2 outputs a non-linear processing means 2A that outputs an image D2A obtained by performing non-linear processing described later on the intermediate image D1, and a high-frequency component image that outputs an image D2B obtained by extracting only the high-frequency component of the image D2A. A generating unit 2B and an adding unit 2C that outputs an image D2C obtained by adding the intermediate image D1 and the image D2B are provided. From the intermediate image processing means 2, the image D2C is output as the intermediate image D2.

以下、実施の形態1による画像処理装置の詳細な動作について説明を行う。
まず、中間画像生成手段1の詳細な動作について説明する。
The detailed operation of the image processing apparatus according to the first embodiment will be described below.
First, the detailed operation of the intermediate image generating unit 1 will be described.

まず、中間画像生成手段1は、高周波数成分画像生成手段1Aにおいて、入力画像DINの高周波数成分のみを取り出した画像D1Aを生成する。高周波数成分の取り出しは、ハイパスフィルタ処理を行うことで可能である。なお、高周波数成分の取り出しは画像の水平方向、垂直方向それぞれについて行う。すなわち高周波数成分画像生成手段1Aは、入力画像DINに対し、水平方向のハイパスフィルタ処理を行って水平方向についてのみ高周波数成分を取り出した画像D1Ahを生成する水平方向高周波数成分画像生成手段1Ahと、垂直方向のハイパスフィルタ処理を行って垂直方向についてのみ高周波数成分を取り出した画像D1Avを生成する垂直方向高周波数成分画像生成手段1Avとを備え、画像D1Aは画像D1Ahと画像D1Avから成る。   First, the intermediate image generating unit 1 generates an image D1A obtained by extracting only the high frequency component of the input image DIN in the high frequency component image generating unit 1A. High frequency components can be extracted by performing high-pass filter processing. Note that high-frequency components are extracted in each of the horizontal and vertical directions of the image. That is, the high-frequency component image generating unit 1A performs horizontal high-pass filter processing on the input image DIN to generate an image D1Ah in which the high-frequency component is extracted only in the horizontal direction. And vertical direction high frequency component image generating means 1Av that performs high-pass filter processing in the vertical direction to generate an image D1Av obtained from high frequency components only in the vertical direction, and the image D1A is composed of an image D1Ah and an image D1Av.

次に、中間画像生成手段1は、低周波数成分画像生成手段1Bにおいて、画像D1Aの低周波数成分のみを取り出した画像D1Bを生成する。低周波数成分の取り出しは、ローパスフィルタ処理を行うことで可能である。なお、低周波数成分の取り出しは水平方向、垂直方向それぞれについて行う。すなわち低周波数成分画像生成手段1Bは、画像D1Ahに対し水平方向のローパスフィルタ処理を行った画像D1Bhを生成する水平方向低周波数成分画像生成手段1Bhと、画像D1Avに対し垂直方向のローパスフィルタ処理を行った画像D1Bvを生成する垂直方向低周波数成分画像生成手段1Bvとを備え、画像D1Bは画像D1Bhと画像D1Bvから成る。そして、中間画像生成手段1からは、画像D1Bが中間画像D1として出力される。なお、中間画像D1は、画像D1Bhに相当する画像D1h、及び画像D1Bvに相当する画像D1vから成る。   Next, the intermediate image generating unit 1 generates an image D1B obtained by extracting only the low frequency component of the image D1A in the low frequency component image generating unit 1B. The low frequency component can be extracted by performing a low pass filter process. The low frequency component is extracted in each of the horizontal direction and the vertical direction. That is, the low-frequency component image generation unit 1B performs horizontal low-frequency component image generation unit 1Bh that generates an image D1Bh obtained by performing horizontal low-pass filter processing on the image D1Ah, and performs low-pass filter processing in the vertical direction on the image D1Av. Vertical direction low frequency component image generation means 1Bv for generating the performed image D1Bv, and the image D1B is composed of an image D1Bh and an image D1Bv. Then, the intermediate image generating unit 1 outputs the image D1B as the intermediate image D1. The intermediate image D1 includes an image D1h corresponding to the image D1Bh and an image D1v corresponding to the image D1Bv.

次に、中間画像処理手段2の詳細な動作について説明する。
まず、中間画像処理手段2は、非線形処理手段2Aにおいて、中間画像D1に対して後述する非線形処理を行った画像D2Aを生成する。非線形処理は、水平方向、垂直方向それぞれについて行う。すなわち非線形処理手段2Aは、画像D1Bhに対して後述する非線形処理を行って画像D2Ahを生成する水平方向非線形処理手段2Ahと、画像D1Bvに対して後述する非線形処理を行って画像D2Avを生成する垂直方向非線形処理手段2Avとを備え、画像D2Aは画像D2Ahと画像D2Avから成る。
Next, the detailed operation of the intermediate image processing means 2 will be described.
First, the intermediate image processing unit 2 generates an image D2A obtained by performing a later-described nonlinear process on the intermediate image D1 in the nonlinear processing unit 2A. Nonlinear processing is performed for each of the horizontal direction and the vertical direction. That is, the nonlinear processing means 2A performs horizontal nonlinear processing means 2Ah that performs nonlinear processing described later on the image D1Bh to generate an image D2Ah, and vertical that performs nonlinear processing described later on image D1Bv to generate the image D2Av. The image D2A includes an image D2Ah and an image D2Av.

非線形処理手段2Aの動作についてさらに詳しく説明する。非線形処理手段2Aは同様の構成を有する水平方向非線形処理手段2Ah、及び垂直方向非線形処理手段2Avを備える。ここで、水平方向非線形処理手段2Ahは水平方向の処理を行い、垂直方向非線形処理手段2Avは垂直方向の処理を行う。   The operation of the nonlinear processing means 2A will be described in more detail. The nonlinear processing means 2A includes a horizontal nonlinear processing means 2Ah and a vertical nonlinear processing means 2Av having the same configuration. Here, the horizontal nonlinear processing means 2Ah performs horizontal processing, and the vertical nonlinear processing means 2Av performs vertical processing.

図2は水平方向非線形処理手段2Ahの構成を表す図である。水平方向非線形処理手段2Ahは、ゼロクロス判定手段311hと、信号増幅手段312hとを備える。なお、非線形処理手段2Ahには、画像D1hが入力画像DIN311hとして入力される。   FIG. 2 is a diagram showing the configuration of the horizontal nonlinear processing means 2Ah. The horizontal non-linear processing unit 2Ah includes a zero-cross determining unit 311h and a signal amplifying unit 312h. Note that the image D1h is input to the nonlinear processing means 2Ah as the input image DIN311h.

ゼロクロス判定手段311hは入力画像DIN311hにおける画素値の変化を水平方向に沿って確認する。そして画素値が正の値から負の値あるいは負の値から正の値へと変化する箇所をゼロクロス点として捉え、信号D311hによってゼロクロス点の前後にある画素(前後において隣接する画素)の位置を信号増幅手段312hに伝達する。例えばゼロクロス点の左右に位置する画素がゼロクロス点の前後に位置する画素として認識される。   The zero-cross determining unit 311h checks the change in the pixel value in the input image DIN 311h along the horizontal direction. Then, a point where the pixel value changes from a positive value to a negative value or from a negative value to a positive value is regarded as a zero cross point, and the positions of pixels (adjacent pixels in the front and rear) before and after the zero cross point by the signal D311h. The signal is transmitted to the signal amplifying means 312h. For example, pixels located on the left and right of the zero cross point are recognized as pixels located before and after the zero cross point.

信号増幅手段312hは信号D311hをもとにゼロクロス点の前後にある画素(前後において隣接する画素)を特定し、ゼロクロス点の前後にある画素についてのみその画素値を増幅させた(絶対値を大きくした)非線形処理画像D312hを生成する。すなわち、ゼロクロス点前後にある画素の画素値に対しては増幅率を1より大きな値とし、それ以外の画素の画素値に対しての増幅率は1とする。   Based on the signal D311h, the signal amplifying unit 312h identifies pixels before and after the zero cross point (adjacent pixels before and after the zero cross point), and amplifies the pixel value only for the pixels before and after the zero cross point (the absolute value is increased). A nonlinear processed image D312h is generated. That is, the amplification factor is set to a value larger than 1 for the pixel values of the pixels before and after the zero cross point, and the amplification factor is set to 1 for the pixel values of the other pixels.

そして水平方向非線形処理手段24hからは画像24hとして非線形処理画像D312hが出力される。   A non-linearly processed image D312h is output from the horizontal non-linear processing means 24h as an image 24h.

図3は垂直方向非線形処理手段2Avの構成を表す図である。垂直方向非線形処理手段2Avは、ゼロクロス判定手段311vと、信号増幅手段312vとを備える。なお、非線形処理手段2Avには、画像D1vが入力画像DIN311vとして入力される。   FIG. 3 is a diagram showing the configuration of the vertical nonlinear processing means 2Av. The vertical nonlinear processing means 2Av includes a zero cross determination means 311v and a signal amplification means 312v. Note that the image D1v is input to the nonlinear processing means 2Av as the input image DIN311v.

ゼロクロス判定手段311vは入力画像DIN311vにおける画素値の変化を垂直方向に沿って確認する。そして画素値が正の値から負の値あるいは負の値から正の値へと変化する箇所をゼロクロス点として捉え、信号D311vによってゼロクロス点の前後にある画素(前後において隣接する画素)の位置を信号増幅手段312vに伝達する。例えばゼロクロス点の上下に位置する画素がゼロクロス点の前後に位置する画素として認識される。   The zero-cross determination unit 311v confirms the change of the pixel value in the input image DIN 311v along the vertical direction. Then, a position where the pixel value changes from a positive value to a negative value or from a negative value to a positive value is regarded as a zero cross point, and the position of pixels (adjacent pixels in the front and rear) before and after the zero cross point is determined by the signal D311v. This is transmitted to the signal amplifying means 312v. For example, pixels positioned above and below the zero cross point are recognized as pixels positioned before and after the zero cross point.

信号増幅手段312vは信号D311vをもとにゼロクロス点の前後にある画素(前後において隣接する画素)を特定し、ゼロクロス点の前後にある画素についてのみその画素値を増幅させた(絶対値を大きくした)非線形処理画像D312vを生成する。すなわち、ゼロクロス点前後にある画素の画素値に対しては増幅率を1より大きな値とし、それ以外の画素の画素値に対しての増幅率は1とする。
以上が非線形処理手段2Aの動作である。
Based on the signal D311v, the signal amplifying unit 312v identifies pixels before and after the zero cross point (adjacent pixels before and after the zero cross point), and amplifies the pixel value only for the pixels before and after the zero cross point (the absolute value is increased). A nonlinear processed image D312v is generated. That is, the amplification factor is set to a value larger than 1 for the pixel values of the pixels before and after the zero cross point, and the amplification factor is set to 1 for the pixel values of the other pixels.
The above is the operation of the nonlinear processing means 2A.

次に、中間画像処理手段2は、高周波数成分画像生成手段2Bにおいて、画像D2Aの高周波数成分のみを取り出した画像D2Bを生成する。高周波数成分の取り出しは、ハイパスフィルタ処理を行うことで可能である。なお、高周波数成分の取り出しは画像の水平方向、垂直方向それぞれについて行う。すなわち高周波数成分画像生成手段2Bは、画像D2Ahに対し水平方向のハイパスフィルタ処理を行った画像D2Bhを生成する水平方向高周波数成分画像生成手段2Bhと、画像D2Avに対し垂直方向のハイパスフィルタ処理を行った画像D2Bvを生成する垂直方向高周波数成分画像生成手段2Bvとを備え、画像D2Bは画像D2Bhと画像D2Bvから成る。   Next, the intermediate image processing unit 2 generates an image D2B obtained by extracting only the high frequency component of the image D2A in the high frequency component image generation unit 2B. High frequency components can be extracted by performing high-pass filter processing. Note that high-frequency components are extracted in each of the horizontal and vertical directions of the image. That is, the high-frequency component image generation unit 2B performs horizontal high-frequency component image generation unit 2Bh that generates an image D2Bh obtained by performing horizontal high-pass filter processing on the image D2Ah, and performs vertical high-pass filter processing on the image D2Av. Vertical direction high frequency component image generation means 2Bv for generating the performed image D2Bv, and the image D2B includes the image D2Bh and the image D2Bv.

次に加算手段2Cは、中間画像D1と画像D2Bを加算して画像D2Cを生成する。なお、中間画像D1は画像D1hおよび画像D1vから成っており、画像D2Bは画像D2Bhおよび画像D2Bvから成っているので、中間画像D1と画像D2Bを加算するとは、画像D1h、D1v、D2Bh、D2Bvの全てを加算することを意味する。そして中間画像処理手段2からは、画像D2Cが中間画像D2として出力される。   Next, the adding means 2C adds the intermediate image D1 and the image D2B to generate an image D2C. The intermediate image D1 is composed of the image D1h and the image D1v, and the image D2B is composed of the image D2Bh and the image D2Bv. Therefore, when the intermediate image D1 and the image D2B are added, It means adding all. From the intermediate image processing means 2, the image D2C is output as the intermediate image D2.

ここで、加算手段2Cでの加算処理は単純加算に限らず重み付け加算を用いてもよい。すなわち、画像D1h、D1v、D2Bh、D2Bvの各々をそれぞれ異なる増幅率で増幅してから加算してもよい。   Here, the addition processing in the adding means 2C is not limited to simple addition, and weighted addition may be used. That is, the images D1h, D1v, D2Bh, and D2Bv may be added after being amplified at different amplification factors.

最後に加算手段3の動作について説明する。加算手段3は入力画像DINと中間画像D2を加算した出力画像DOUTを生成する。そして出力画像DOUTが最終的な出力画像として、画像処理装置から出力される。   Finally, the operation of the adding means 3 will be described. The adding means 3 generates an output image DOUT obtained by adding the input image DIN and the intermediate image D2. The output image DOUT is output from the image processing apparatus as a final output image.

以下、本発明における画像処理装置を画像表示装置の一部として利用する例について説明する。この説明を通じて、本発明における画像処理装置の作用、効果も明らかなものとなるであろう。なお、以下の説明では特に断らない限り、符号Fnは入力画像DINのナイキスト周波数を表す。   Hereinafter, an example in which the image processing apparatus according to the present invention is used as a part of an image display apparatus will be described. Through this description, the operation and effect of the image processing apparatus according to the present invention will become clear. In the following description, unless otherwise specified, the symbol Fn represents the Nyquist frequency of the input image DIN.

図4は本発明における画像処理装置を利用した画像表示装置であり、モニタU3上に原画DORGに対応した画像が表示される。   FIG. 4 shows an image display apparatus using the image processing apparatus according to the present invention, and an image corresponding to the original image DORG is displayed on the monitor U3.

画像拡大手段U1は、原画DORGの画像サイズがモニタU3の画像サイズより小さい場合、原画DORGを拡大した画像DU1を出力する。ここで画像を拡大する手段としては、バイキュービック法などを用いることができる。   When the image size of the original image DORG is smaller than the image size of the monitor U3, the image enlargement unit U1 outputs an image DU1 obtained by enlarging the original image DORG. Here, a bicubic method or the like can be used as means for enlarging the image.

本発明における画像処理装置U2は、画像DU1に対し、先に説明した処理を行った画像DU2を出力する。そしてモニタU8上には画像DU2が表示される。   The image processing apparatus U2 in the present invention outputs an image DU2 obtained by performing the above-described processing on the image DU1. The image DU2 is displayed on the monitor U8.

以下、原画DORGは、水平方向、垂直方向ともその画素数がモニタU3の画素数の半分であるとして、まず画像拡大手段U1の動作、作用について説明を行う。   Hereinafter, assuming that the number of pixels of the original image DORG is half that of the monitor U3 in both the horizontal direction and the vertical direction, the operation and action of the image enlarging means U1 will be described first.

図5は画像拡大手段U1の構成および動作を表す図であり、画像拡大手段U1は水平方向ゼロ挿入手段U1A、水平方向低周波数成分通過手段U1B、垂直方向ゼロ挿入手段U1C、及び垂直方向低周波数成分通過手段U1Dを備える。水平方向ゼロ挿入手段U1Aは原画DORGの水平方向に関して画素値0を持つ画素を適宜挿入した画像DU1Aを生成する。水平方向低周波数成分通過手段U1Bはローパスフィルタ処理により画像DU1Aの低周波数成分のみを取り出した画像DU1Bを生成する。垂直方向ゼロ挿入手段U1Cは画像DU1Bの垂直方向に関して画素値0を持つ画素を適宜挿入した画像DU1Cを生成する。垂直方向低周波数成分通過手段DU1Dは画像DU1Cの低周波数成分のみを取り出した画像DU1Dを生成する。そして画像DU1Dが原画DORGを水平方向、垂直方向とも2倍した画像DU1として、画像拡大手段U1から出力される。   FIG. 5 is a diagram showing the configuration and operation of the image enlarging means U1. The image enlarging means U1 includes a horizontal direction zero insertion means U1A, a horizontal direction low frequency component passing means U1B, a vertical direction zero insertion means U1C, and a vertical direction low frequency. Component passing means U1D is provided. The horizontal zero insertion means U1A generates an image DU1A in which pixels having a pixel value of 0 in the horizontal direction of the original image DORG are appropriately inserted. The horizontal direction low frequency component passing means U1B generates an image DU1B in which only the low frequency component of the image DU1A is extracted by low pass filter processing. The vertical zero insertion unit U1C generates an image DU1C in which pixels having a pixel value of 0 in the vertical direction of the image DU1B are appropriately inserted. The vertical direction low frequency component passing means DU1D generates an image DU1D obtained by extracting only the low frequency component of the image DU1C. The image DU1D is output from the image enlarging means U1 as an image DU1 obtained by doubling the original image DORG in both the horizontal and vertical directions.

図6(A)〜(E)は画像拡大手段U1の動作を詳しく説明するための図であり、図6(A)は原画DORGを、図6(B)は画像DU1Aを、図6(C)は画像DU1Bを、図6(D)は画像DU1Cを、図6(E)は画像DU1Dを表す。図6(A)〜(E)に関して、四角(各升目)は画素を表し、その中に書かれた記号あるいは数値は各画素の画素値を表す。   6A to 6E are diagrams for explaining the operation of the image enlarging means U1 in detail. FIG. 6A shows the original image DORG, FIG. 6B shows the image DU1A, and FIG. ) Represents the image DU1B, FIG. 6D represents the image DU1C, and FIG. 6E represents the image DU1D. 6A to 6E, a square (each square) represents a pixel, and a symbol or a numerical value written therein represents a pixel value of each pixel.

水平方向ゼロ挿入手段U1Aは図6(A)に示す原画DORGに対して、水平方向の1画素につき1個、画素値0をもった画素を挿入し(即ち、原画DORGの水平方向に隣接する画素列相互間に一つの、画素値0の画素から成る画素列を挿入し)、図6(B)に示す画像DU1Aを生成する。水平方向低周波数成分通過手段U1Bは図6(B)に示す画像DU1Aに対して、ローパスフィルタ処理を施し、図6(C)に示す画像DU1Bを生成する。垂直方向ゼロ挿入手段U1Cは図6(C)に示す画像DU1Bに対して、垂直方向の1画素につき1個、画素値0をもった画素を挿入し(即ち、画像DU1Bの垂直方向に隣接する画素行相互間に一つの、画素値0の画素から成る画素行を挿入し)、図6(D)に示す画像DU1Cを生成する。垂直方向低周波数成分通過手段U1Dは図6(D)に示す画像DU1Cに対して、ローパスフィルタ処理を施し、図6(E)に示す画像DU1Dを生成する。以上の処理により原画DORGを水平方向、垂直方向とも2倍に拡大した画像DU1Dが生成される。   The horizontal zero insertion means U1A inserts one pixel having a pixel value of 0 for each pixel in the horizontal direction into the original image DORG shown in FIG. 6A (that is, adjacent to the original image DORG in the horizontal direction). An image DU1A shown in FIG. 6B is generated by inserting a single pixel column composed of pixels having a pixel value of 0 between the pixel columns. The horizontal direction low frequency component passing means U1B performs a low-pass filter process on the image DU1A shown in FIG. 6B to generate an image DU1B shown in FIG. The vertical zero insertion means U1C inserts one pixel with a pixel value of 0 for each pixel in the vertical direction into the image DU1B shown in FIG. 6C (that is, adjacent to the image DU1B in the vertical direction). An image DU1C shown in FIG. 6D is generated by inserting one pixel row composed of pixels having a pixel value of 0 between the pixel rows. The vertical low frequency component passing means U1D performs low-pass filter processing on the image DU1C shown in FIG. 6D to generate an image DU1D shown in FIG. Through the above processing, an image DU1D in which the original image DORG is doubled in both the horizontal direction and the vertical direction is generated.

図7(A)〜(D)は、画像拡大手段U1による処理の作用を周波数空間上で表したものであり、図7(A)は原画DORGの周波数スペクトル、図7(B)は画像DU1Aの周波数スペクトル、図7(C)は水平方向周波数成分通過手段U1Bの周波数応答、図7(D)は画像DU1Bの周波数スペクトルを表している。なお、図7(A)〜(D)において横軸は水平方向の空間周波数を表す周波数軸であり、縦軸は周波数スペクトルもしくは周波数応答の強度を表している。なお原画DORGの画素数は入力画像DINの半分となっており、言い換えると原画DORGのサンプリング間隔は入力画像DINのサンプリング間隔の2倍になっている。したがって原画DORGのナイキスト周波数は入力画像DINのナイキスト周波数の半分すなわち、Fn/2となる。   FIGS. 7A to 7D show the effect of processing by the image enlarging means U1 on the frequency space, FIG. 7A shows the frequency spectrum of the original picture DORG, and FIG. 7B shows the image DU1A. 7C shows the frequency response of the horizontal frequency component passing means U1B, and FIG. 7D shows the frequency spectrum of the image DU1B. 7A to 7D, the horizontal axis is the frequency axis representing the spatial frequency in the horizontal direction, and the vertical axis represents the frequency spectrum or the intensity of the frequency response. Note that the number of pixels of the original image DORG is half that of the input image DIN. In other words, the sampling interval of the original image DORG is twice the sampling interval of the input image DIN. Therefore, the Nyquist frequency of the original image DORG is half of the Nyquist frequency of the input image DIN, that is, Fn / 2.

なお、図7(A)〜(D)では表記を簡素にするため、1本の周波数軸しか用いていない。しかしながら、通常、画像データは2次元平面状に並んだ画素配列上に与えられた画素値から成り、その周波数スペクトルも水平方向の周波数軸および垂直方向の周波数軸で張られる平面上に与えられるものである。したがって原画DORG等の周波数スペクトル等を正確に表すためには、水平方向の周波数軸および垂直方向の周波数軸の両方を記載する必要がある。しかしながらその周波数スペクトルの形状は通常、周波数軸上の原点を中心に等方的に広がったものであり、周波数軸1本で張られる空間上での周波数スペクトルを示しさえすれば、そこから周波数軸2本で張られる空間へ拡張して考察することは当業者にとって容易である。したがって以降の説明でも特に断らない限り、周波数空間上での説明は、1本の周波数軸で張られる空間を用いて行う。   In FIGS. 7A to 7D, only one frequency axis is used to simplify the notation. However, usually, image data consists of pixel values given on a pixel array arranged in a two-dimensional plane, and its frequency spectrum is also given on a plane stretched by a horizontal frequency axis and a vertical frequency axis. It is. Therefore, in order to accurately represent the frequency spectrum of the original picture DORG or the like, it is necessary to describe both the horizontal frequency axis and the vertical frequency axis. However, the shape of the frequency spectrum usually spreads isotropically around the origin on the frequency axis, and as long as the frequency spectrum in the space spanned by one frequency axis is shown, the frequency axis It is easy for those skilled in the art to expand and consider the space spanned by two. Therefore, unless otherwise specified in the following description, the description on the frequency space is performed using a space stretched by one frequency axis.

まず、原画DORGの周波数スペクトルについて説明する。通常、自然画像が原画DORGとして入力されるがそのスペクトル強度は周波数空間の原点周辺に集中している。したがって原画DORGの周波数スペクトルは図7(A)のように表すスペクトルSPOのようになる。   First, the frequency spectrum of the original picture DORG will be described. Normally, a natural image is input as the original image DORG, but its spectral intensity is concentrated around the origin of the frequency space. Therefore, the frequency spectrum of the original picture DORG is a spectrum SPO represented as shown in FIG.

次に、画像DU1Aのスペクトル強度について説明する。画像DU1Aは、原画DORGに対して、水平方向に1画素につき1画素、画素値0を持った画素を挿入することで生成される。このような処理を行うと周波数スペクトルには原画DORGのナイキスト周波数を中心にした折り返しが発生する。すなわち周波数±Fn/2を中心にスペクトルSPOが折り返したスペクトルSPMが発生するので、画像DU1Aの周波数スペクトルは図7(B)のように表される。   Next, the spectral intensity of the image DU1A will be described. The image DU1A is generated by inserting one pixel per pixel and a pixel value of 0 in the horizontal direction with respect to the original image DORG. When such processing is performed, aliasing around the Nyquist frequency of the original picture DORG occurs in the frequency spectrum. That is, since a spectrum SPM in which the spectrum SPO is folded around the frequency ± Fn / 2 is generated, the frequency spectrum of the image DU1A is expressed as shown in FIG.

次に、水平方向低周波数成分通過手段U1Bの周波数応答について説明する。水平方向低周波数成分通過手段U1Bはローパスフィルタによって実現されるので、その周波数応答は図7(C)に示すように周波数が高くなるほど低くなる。   Next, the frequency response of the horizontal low frequency component passing means U1B will be described. Since the horizontal low frequency component passing means U1B is realized by a low-pass filter, its frequency response becomes lower as the frequency becomes higher as shown in FIG.

最後に、画像DU1Bの周波数スペクトルについて説明する。図7(B)に示す周波数スペクトルを持った画像DU1Aに対し、図7(C)に示した周波数応答を持ったローパスフィルタ処理を行うことで、画像DU1Bが得られる。したがって画像DU1Bの周波数スペクトルは画像DU1Bに示すように、スペクトルSPMの強度がある程度落ちたスペクトルSP2と、スペクトルSPOの強度がある程度落ちたスペクトルSP1から成る。なお一般に、ローパスフィルタの周波数応答は周波数が高くなるほど低くなる。従って、スペクトルSP1の強度をスペクトルSPOと比較すると、水平方向低周波数成分通過手段U1Bによって、高周波数成分側、すなわち周波数が±Fn/2近傍でのスペクトル強度が減少したものとなる。   Finally, the frequency spectrum of the image DU1B will be described. The image DU1B is obtained by performing low-pass filter processing having the frequency response shown in FIG. 7C on the image DU1A having the frequency spectrum shown in FIG. 7B. Therefore, as shown in the image DU1B, the frequency spectrum of the image DU1B includes a spectrum SP2 in which the intensity of the spectrum SPM has dropped to some extent and a spectrum SP1 in which the intensity of the spectrum SPO has dropped to some extent. In general, the frequency response of the low-pass filter decreases as the frequency increases. Therefore, when the intensity of the spectrum SP1 is compared with the spectrum SPO, the spectrum intensity on the high frequency component side, that is, the frequency near ± Fn / 2 is reduced by the horizontal low frequency component passing means U1B.

また、画像拡大手段U1による処理のうち、垂直方向ゼロ挿入手段U1Cおよび垂直方向低周波数成分通過手段U1Dによる処理について、その周波数空間上での作用についての説明は省略するが、その処理の内容から、垂直方向の空間周波数を表す軸方向に対して、図7(A)〜(D)を用いて説明した内容と同様の作用があることは容易に理解できる。すなわち、画像DU1Dの周波数スペクトルは、図7(D)に示した周波数スペクトルが2次元上に広がったものとなる。   Of the processing by the image enlarging means U1, the description of the operation on the frequency space of the processing by the vertical zero insertion means U1C and the vertical low frequency component passing means U1D is omitted, but from the contents of the processing. It can be easily understood that there is an action similar to that described with reference to FIGS. 7A to 7D with respect to the axial direction representing the spatial frequency in the vertical direction. That is, the frequency spectrum of the image DU1D is obtained by spreading the frequency spectrum shown in FIG. 7D two-dimensionally.

また、以降の説明ではスペクトルSP2のことを折り返し成分と呼ぶ。この折り返し成分は、画像上では、比較的高い周波数成分を持ったノイズあるいは偽の信号として現れる。そのようなノイズあるいは偽の信号としてオーバーシュートやジャギーあるいはリンギング等が挙げられる。   In the following description, the spectrum SP2 is referred to as a folded component. This aliasing component appears on the image as noise or a false signal having a relatively high frequency component. Such noise or false signals include overshoot, jaggy or ringing.

以下、本発明における画像処理装置の作用、効果について説明する。
図8(A)〜(E)は入力画像DIN(もしくは画像DU1)として原画DORGを拡大して得られた画像DU1Dが入力された場合の、入力画像DINから中間画像D1を生成する際の作用、効果を模式的に表した図であり、図8(A)は入力画像DINの周波数スペクトルを、図8(B)は高周波数成分画像生成手段1の周波数応答を、図8(C)は低周波数成分画像生成手段2の周波数応答を、図8(D)は中間画像生成手段1の周波数応答を、図8(E)は中間画像D1の周波数スペクトルを表す。なお、図8(A)〜(E)においても図7(A)〜(D)と同様の理由で周波数軸は1本しか用いていない。
The operation and effect of the image processing apparatus according to the present invention will be described below.
8A to 8E show the operation when the intermediate image D1 is generated from the input image DIN when the image DU1D obtained by enlarging the original image DORG is input as the input image DIN (or image DU1). 8A schematically shows the effect, FIG. 8A shows the frequency spectrum of the input image DIN, FIG. 8B shows the frequency response of the high frequency component image generating means 1, and FIG. 8C shows the frequency response. FIG. 8D shows the frequency response of the low-frequency component image generation means 2, FIG. 8D shows the frequency response of the intermediate image generation means 1, and FIG. 8E shows the frequency spectrum of the intermediate image D1. In FIGS. 8A to 8E, only one frequency axis is used for the same reason as in FIGS. 7A to 7D.

さらに図8(A)〜(E)では、空間周波数が0以上となる範囲でのみ周波数スペクトルあるいは周波数応答の強度を表しているが、以下の説明での周波数スペクトルあるいは周波数応答は、周波数軸上の原点を中心に対称的な形状となる。したがって説明に用いる図は、空間周波数が0以上となる範囲のみを示したもので十分である。   Further, in FIGS. 8A to 8E, the intensity of the frequency spectrum or the frequency response is shown only in the range where the spatial frequency is 0 or more, but the frequency spectrum or the frequency response in the following description is on the frequency axis. It becomes a symmetric shape around the origin. Therefore, the figure used for description is sufficient to show only the range where the spatial frequency is 0 or more.

まず、入力画像DINの周波数スペクトルについて説明する。画像DU1Dが入力画像DINとして入力されるので、入力画像DINの周波数スペクトルは図8(A)に示すように、周波数スペクトルは図7(D)で説明したものと同じ形状となり、原画DORGのスペクトルSPOの強度がある程度落ちたスペクトルSP1と折り返し成分となるスペクトルSP2から成る。   First, the frequency spectrum of the input image DIN will be described. Since the image DU1D is input as the input image DIN, the frequency spectrum of the input image DIN has the same shape as that illustrated in FIG. 7D as shown in FIG. 8A, and the spectrum of the original image DORG. It consists of a spectrum SP1 in which the intensity of SPO has dropped to some extent and a spectrum SP2 that is a folded component.

次に、高周波数成分画像生成手段1Aの周波数応答について説明する。高周波数成分画像生成手段1Aはハイパスフィルタにより構成されているので、その周波数応答は図8(B)に示すように周波数が低くなるほど低くなる。   Next, the frequency response of the high frequency component image generating unit 1A will be described. Since the high frequency component image generating means 1A is constituted by a high pass filter, the frequency response thereof becomes lower as the frequency becomes lower as shown in FIG. 8B.

次に、低周波数成分画像生成手段1Bの周波数応答について説明する。低周波数成分画像生成手段1Bはローパスフィルタにより構成されているので、その周波数応答は図8(C)に示すように周波数が高くなるほど低くなる。   Next, the frequency response of the low frequency component image generation means 1B will be described. Since the low frequency component image generating means 1B is composed of a low-pass filter, the frequency response becomes lower as the frequency becomes higher as shown in FIG. 8C.

次に、中間画像生成手段1の周波数応答について説明する。入力画像DINが持つ周波数成分のうち、図8(D)に示された低周波数成分側の領域RL1の周波数成分については、中間画像生成手段1内の高周波数成分画像生成手段1Aで弱められる。一方、図8(D)に示された高周波数成分側の領域RH1の周波数成分については、中間画像生成手段1内の低周波数成分画像生成手段1Bで弱められる。したがって、中間画像生成手段1の周波数応答は、図8(D)に示すように、低周波数成分側の領域RL1と高周波数成分側の領域RH1によって帯域を制限された中間の領域RM1にピークを持ったものとなる。
この中間の領域RM1は、原画DORGに画素値0を持った画素を挿入することに伴って発生される折り返し成分を含まないものであり、原画DORGのナイキスト周波数Fn/2以下の領域の一部を占める。
Next, the frequency response of the intermediate image generating unit 1 will be described. Among the frequency components of the input image DIN, the frequency components in the low frequency component side region RL1 shown in FIG. 8D are weakened by the high frequency component image generating unit 1A in the intermediate image generating unit 1. On the other hand, the frequency component in the region RH1 on the high frequency component side shown in FIG. 8D is weakened by the low frequency component image generating unit 1B in the intermediate image generating unit 1. Therefore, as shown in FIG. 8D, the frequency response of the intermediate image generating means 1 has a peak in the intermediate region RM1 whose band is limited by the region RL1 on the low frequency component side and the region RH1 on the high frequency component side. It will have.
This intermediate region RM1 does not include the aliasing component generated when a pixel having a pixel value of 0 is inserted into the original image DORG, and is a part of the region below the Nyquist frequency Fn / 2 of the original image DORG. Occupy.

次に、中間画像D1の周波数スペクトルについて説明する。図8(A)に示す周波数スペクトルを持つ入力画像DINが、図8(D)に示した周波数応答を持つ中間画像生成手段1を通過することで、中間画像D1が得られる。そして中間画像生成手段1の周波数応答は、低周波数成分側の領域RL1と高周波数成分側の領域RH1によって帯域制限された中間の領域RM1にピークを持ったものなので、中間画像D1の周波数スペクトルは、入力画像DINの周波数スペクトルのうち、低周波数成分側の領域RL1と高周波数成分側の領域RH1に含まれる部分の強度が弱くなったものとなる。従って中間画像D1は入力画像DINの持つ高周波数成分から折り返し成分となるスペクトルSP1を取り除いたものとなる。すなわち中間画像生成手段1には、入力画像DINのもつ高周波数成分から折り返し成分となるスペクトルSP1を取り除いた中間画像D1を生成するという効果がある。   Next, the frequency spectrum of the intermediate image D1 will be described. The input image DIN having the frequency spectrum shown in FIG. 8A passes through the intermediate image generating means 1 having the frequency response shown in FIG. 8D, whereby the intermediate image D1 is obtained. The frequency response of the intermediate image generating means 1 has a peak in the intermediate region RM1 band-limited by the region RL1 on the low frequency component side and the region RH1 on the high frequency component side, so the frequency spectrum of the intermediate image D1 is In the frequency spectrum of the input image DIN, the intensity of the portion included in the low frequency component side region RL1 and the high frequency component side region RH1 is weakened. Therefore, the intermediate image D1 is obtained by removing the spectrum SP1 that is the aliasing component from the high frequency component of the input image DIN. That is, the intermediate image generating means 1 has an effect of generating the intermediate image D1 by removing the spectrum SP1 that is the aliasing component from the high frequency component of the input image DIN.

図9(A)〜(C)は中間画像処理手段2の作用、効果を表した図であり、図9(A)は非線形処理画像D2Aの周波数スペクトルを、図9(B)は高周波数成分画像2Bの周波数応答を、図9(C)は画像D2Bの周波数スペクトルを表す。なお、図9(A)〜(C)では、図8(A)〜(E)と同様の理由で、空間周波数が0以上となる範囲でのみ周波数スペクトルあるいは周波数応答の強度を表している。   FIGS. 9A to 9C are diagrams showing the operation and effect of the intermediate image processing means 2, FIG. 9A shows the frequency spectrum of the nonlinear processed image D2A, and FIG. 9B shows the high frequency component. FIG. 9C shows the frequency response of the image 2B, and FIG. 9C shows the frequency spectrum of the image D2B. In FIGS. 9A to 9C, for the same reason as in FIGS. 8A to 8E, the intensity of the frequency spectrum or frequency response is shown only in the range where the spatial frequency is 0 or more.

後述するように非線形処理画像D2Aでは、高周波数成分側の領域RH2に相当する高周波数成分が生成される。図9(A)はその様子を模式的に表した図である。画像D2Bは非線形処理画像D2Aが高周波数成分画像生成手段2Bを通過することで生成される。高周波数成分画像生成手段2Bはハイパスフィルタで構成されており、その周波数応答は図9(B)に示すように周波数が高くなるほど高いものとなる。従って画像D2Bの周波数スペクトルは図9(C)に示すように非線形処理画像D2Aの周波数スペクトルから低周波数成分側の領域RL2に相当する成分を取り除いたものとなる。言い換えると、非線形処理手段2Aには高周波数成分側の領域RH2に相当する高周波数成分を生成する効果があり、高周波数成分画像生成手段2Bには非線形処理手段2Aで生成された高周波数成分のみを取り出す効果がある。   As will be described later, in the nonlinear processed image D2A, a high frequency component corresponding to the region RH2 on the high frequency component side is generated. FIG. 9A is a diagram schematically showing the state. The image D2B is generated by passing the nonlinear processed image D2A through the high frequency component image generating means 2B. The high frequency component image generating means 2B is composed of a high-pass filter, and the frequency response becomes higher as the frequency becomes higher as shown in FIG. 9B. Therefore, the frequency spectrum of the image D2B is obtained by removing the component corresponding to the region RL2 on the low frequency component side from the frequency spectrum of the nonlinear processed image D2A as shown in FIG. 9C. In other words, the nonlinear processing means 2A has an effect of generating a high frequency component corresponding to the region RH2 on the high frequency component side, and the high frequency component image generating means 2B has only the high frequency component generated by the nonlinear processing means 2A. There is an effect to take out.

上記の作用、効果についてさらに詳しく説明する。
図10(A)〜(C)及び図11(A)〜(C)は、ステップエッジと、ステップエッジをサンプリングした際に得られる相連続する画素の信号の値を示す図である。
図10(A)はステップエッジと、サンプリング間隔S1を表しており、図10(B)はステップエッジをサンプリング間隔S1でサンプリングした際に得られる信号を表しており、図10(C)は図10(B)に表された信号の高周波数成分を表している。一方、図11(A)はステップエッジとサンプリング間隔S1より間隔の広いサンプリング間隔S2を表しており、図11(B)はステップエッジをサンプリング間隔S2でサンプリングした際に得られる信号を表しており、図11(C)は図11(B)に表された信号の高周波数成分を表している。なお、以下の説明ではサンプリング間隔S2の長さはサンプリング間隔S1の長さの半分であるとする。
The above operations and effects will be described in more detail.
FIGS. 10A to 10C and FIGS. 11A to 11C are diagrams showing step edges and signal values of pixels that are continuous when obtained by sampling the step edges.
FIG. 10A shows the step edge and the sampling interval S1, FIG. 10B shows the signal obtained when the step edge is sampled at the sampling interval S1, and FIG. 10 (B) represents the high frequency component of the signal. On the other hand, FIG. 11A shows a sampling interval S2 that is wider than the step edge and the sampling interval S1, and FIG. 11B shows a signal obtained when the step edge is sampled at the sampling interval S2. FIG. 11C shows high frequency components of the signal shown in FIG. In the following description, it is assumed that the length of the sampling interval S2 is half of the length of the sampling interval S1.

図10(C)、図11(C)に表されるようにステップエッジの中央は高周波数成分を表した信号においてゼロクロス点Zとして現れる。また、高周波数成分を表した信号のゼロクロス点Zの近傍での傾きは、サンプリング間隔が短いほど急になり、かつゼロクロス点Z近傍での局所的な最大値、最小値を与える点の位置も、サンプリング間隔が短いほどゼロクロス点Zに近づく。   As shown in FIGS. 10C and 11C, the center of the step edge appears as a zero cross point Z in the signal representing the high frequency component. In addition, the slope of the signal representing the high frequency component near the zero cross point Z becomes steeper as the sampling interval is short, and the position of the point giving the local maximum and minimum values near the zero cross point Z is also The shorter the sampling interval, the closer to the zero cross point Z.

すなわち、サンプリング間隔が変わっても、エッジ近傍において高周波数成分を表す信号のゼロクロス点の位置は変化しないが、サンプリング間隔が小さくなるほど(あるいは解像度が上がるほど)エッジ近傍での高周波数成分の傾きは急になり、局所的な最大値、最小値を与える点の位置はゼロクロス点に近づく。   That is, even if the sampling interval changes, the position of the zero crossing point of the signal representing the high frequency component does not change near the edge, but the slope of the high frequency component near the edge decreases as the sampling interval decreases (or the resolution increases). The position of the point that gives the local maximum and minimum values approaches the zero-cross point.

図12(A)〜(F)はステップエッジをサンプリング間隔S2でサンプリングした信号が2倍に拡大された後、本発明における画像処理装置に入力されるときの、中間画像生成手段1および中間画像処理手段2の動作を表している。なお、先に述べた通り、中間画像生成手段1および中間画像処理手段2の内部の処理は水平方向、垂直方向のそれぞれについて行われるのでその処理は一次元的に行われる。したがって図12(A)〜(F)では一次元信号を用いて処理の内容を表している。   12A to 12F show the intermediate image generating means 1 and the intermediate image when the signal obtained by sampling the step edge at the sampling interval S2 is doubled and then input to the image processing apparatus according to the present invention. The operation of the processing means 2 is shown. As described above, since the internal processing of the intermediate image generating unit 1 and the intermediate image processing unit 2 is performed in each of the horizontal direction and the vertical direction, the processing is performed one-dimensionally. Therefore, in FIGS. 12A to 12F, the contents of the processing are represented using a one-dimensional signal.

図12(A)は、図11(B)と同様に、ステップエッジをサンプリング間隔S2でサンプリングした信号である。図12(B)は、図12(A)に表した信号を2倍に拡大した信号である。すなわち、原画DORGに図12(A)に示すようなエッジが含まれる場合、入力画像DINとして図12(B)に示すような信号が入力される。なお、信号を2倍に拡大するとサンプリング間隔は拡大前の半分になるため、図12(B)に表した信号のサンプリング間隔は図10(A)〜(C)中のサンプリング間隔S1と同じになる。また、図12(A)において座標P3で表される位置はエッジ信号の低輝度側の境界部分であり、座標P4で表される位置はエッジ信号の高輝度側の境界部分である。   FIG. 12A shows a signal obtained by sampling the step edge at the sampling interval S2, as in FIG. 11B. FIG. 12B is a signal obtained by enlarging the signal shown in FIG. That is, when the original image DORG includes an edge as shown in FIG. 12A, a signal as shown in FIG. 12B is input as the input image DIN. Note that when the signal is doubled, the sampling interval is half that before the expansion, and therefore the sampling interval of the signal shown in FIG. 12B is the same as the sampling interval S1 in FIGS. 10A to 10C. Become. In FIG. 12A, the position represented by the coordinate P3 is the boundary portion on the low luminance side of the edge signal, and the position represented by the coordinate P4 is the boundary portion on the high luminance side of the edge signal.

図12(C)は図12(B)に表した信号の高周波数成分を表した信号、すなわち高周波数成分画像生成手段1Aから出力される画像D1Aに相当する信号である。なお、画像D1Aは、入力画像DINの高周波数成分を取り出したものなので、その中には折り返し成分も含まれている。   FIG. 12C shows a signal representing the high frequency component of the signal shown in FIG. 12B, that is, a signal corresponding to the image D1A output from the high frequency component image generating means 1A. Note that the image D1A is obtained by extracting the high-frequency component of the input image DIN, and therefore includes an aliasing component.

図12(D)は図12(C)に表した信号の低周波数成分を表した信号、すなわち低周波数成分画像生成手段1Bから出力される画像D1Bに相当する信号である。なお先に述べたとおり画像D1Bが中間画像D1として出力されるので、図12(D)は中間画像D1にも相当する。図12(D)に示すとおり、中間画像D1においてゼロクロス点Z近傍の局所的な最小値は座標P3に、局所的な最大値は座標P4に表れ、その様子は図11(C)に示した、ステップエッジをサンプリング間隔S2でサンプリングした信号から取り出した高周波数成分と一致する。また、画像D1Aに含まれていた折り返し成分は、低周波数成分画像生成手段1Bで行うローパスフィルタ処理によって取り除かれる。   FIG. 12D shows a signal representing the low frequency component of the signal shown in FIG. 12C, that is, a signal corresponding to the image D1B output from the low frequency component image generating means 1B. Since the image D1B is output as the intermediate image D1 as described above, FIG. 12D corresponds to the intermediate image D1. As shown in FIG. 12D, the local minimum value in the vicinity of the zero cross point Z in the intermediate image D1 appears at the coordinate P3, and the local maximum value appears at the coordinate P4, which is shown in FIG. 11C. The step edge coincides with the high frequency component extracted from the signal sampled at the sampling interval S2. Further, the aliasing component included in the image D1A is removed by a low-pass filter process performed by the low-frequency component image generation unit 1B.

図12(E)は、図12(D)に表した信号に対する非線形処理手段2Aに入力された際の出力信号、すなわち、中間画像D1が入力された場合に非線形処理手段2Aから出力される画像D2Aを表している。非線形処理手段2Aではゼロクロス点Zの前後の座標P1、P2の信号値が増幅される。したがって、画像D2Aは図12(E)に示すように座標P1、P2での信号値の大きさが他の値に比べ大きくなり、ゼロクロス点Z近傍で、局所的な最小値の現れる位置が座標P3からよりゼロクロス点Zに近い座標P1に、局所的な最大値の現れる位置が座標P4からよりゼロクロス点Zに近い座標P1へと変化する。これは非線形処理手段2Aにおける、ゼロクロス点Z前後の画素の値を増幅するという非線形処理によって、高周波数成分が生成されたことを意味する。このように画素ごとに適応的に増幅率を変える、あるいは画素に応じて処理の内容を適宜変えることで、高周波数成分を生成することが可能になる。すなわち非線形処理手段2Aには、中間画像D1には含まれない高周波数成分、すわなち、図9(A)に示した高周波数成分側の領域RH2に相当する高周波数成分を生成する効果がある。   FIG. 12E shows an output signal when the signal shown in FIG. 12D is input to the nonlinear processing means 2A, that is, an image output from the nonlinear processing means 2A when the intermediate image D1 is input. D2A is represented. In the nonlinear processing means 2A, the signal values of the coordinates P1 and P2 before and after the zero cross point Z are amplified. Therefore, as shown in FIG. 12E, in the image D2A, the magnitude of the signal value at the coordinates P1 and P2 is larger than the other values, and the position where the local minimum value appears near the zero cross point Z is the coordinate. The position where the local maximum value appears at the coordinate P1 closer to the zero cross point Z from P3 changes from the coordinate P4 to the coordinate P1 closer to the zero cross point Z. This means that the high frequency component is generated by the nonlinear processing in the nonlinear processing means 2A that amplifies the values of the pixels before and after the zero cross point Z. In this way, it is possible to generate a high-frequency component by adaptively changing the amplification factor for each pixel or appropriately changing the content of processing according to the pixel. In other words, the nonlinear processing means 2A has an effect of generating a high frequency component not included in the intermediate image D1, that is, a high frequency component corresponding to the region RH2 on the high frequency component side shown in FIG. 9A. is there.

図12(F)は図12(E)に表した信号の高周波数成分を表した信号、すなわち高周波数成分画像生成手段2Bから出力される画像D2Bに相当する信号である。図12(F)に示すとおり、画像D2Bにおいてゼロクロス点Z近傍の局所的な最小値は座標P1に、最大値は座標P2に表れ、その様子は図10(C)に示した、ステップエッジをサンプリング間隔S1でサンプリングした信号から取り出した高周波数成分と一致する。これは非線形処理手段2Aにおいて生成された高周波数成分が高周波数成分画像生成手段2Bによって取り出され、画像D2Bとして出力されることを意味する。また、取り出された画像D2Bはサンプリング間隔S1に対応した周波数成分を含む信号であるといえる。言い換えると、高周波数成分画像生成手段2Bには非線形処理手段2Aで生成された高周波数成分のみを取り出す効果がある。   FIG. 12F shows a signal representing the high frequency component of the signal shown in FIG. 12E, that is, a signal corresponding to the image D2B output from the high frequency component image generating means 2B. As shown in FIG. 12F, the local minimum value in the vicinity of the zero-cross point Z in the image D2B appears at the coordinate P1, and the maximum value appears at the coordinate P2. The state is shown by the step edge shown in FIG. This coincides with the high frequency component extracted from the signal sampled at the sampling interval S1. This means that the high frequency component generated in the nonlinear processing means 2A is taken out by the high frequency component image generating means 2B and output as an image D2B. Further, it can be said that the extracted image D2B is a signal including a frequency component corresponding to the sampling interval S1. In other words, the high frequency component image generation means 2B has an effect of extracting only the high frequency component generated by the nonlinear processing means 2A.

加算手段2Cでは中間画像D1と画像D2Bを加算し画像D2Cを生成する。先に述べたとおり中間画像D1は入力画像DINの持つ高周波数成分から折り返し成分を取り除いたものであり、図8(E)に示すように原画DORGのナイキスト周波数近傍の高周波数成分に対応している。図7(D)で説明したとおり、原画DORGのナイキスト周波数近傍のスペクトル強度は画像拡大手段U1での拡大処理によって弱められているので、中間画像D1を加算することで、拡大処理によって弱められたスペクトル強度を補うことができる。また、中間画像D1から、折り返し成分は取り除かれているので、オーバーシュートやジャギーあるいはリンギングといった偽の信号を強調することはない。一方、画像D2Bはサンプリング間隔S1に対応した高周波数成分である。したがって画像D2Cを加算することで原画DORGのナイキスト周波数以上の帯域の高周波数成分を与えることがでるので、画像の解像感を増すことができる。したがって中間画像D1と画像D2Bを加算した画像D2Cを入力画像DINに加算することで、折り返し成分を強調することなく高周波数成分を加算することが可能となり、画像の解像感を高めることが可能となる。   The adding means 2C adds the intermediate image D1 and the image D2B to generate an image D2C. As described above, the intermediate image D1 is obtained by removing the aliasing component from the high frequency component of the input image DIN, and corresponds to the high frequency component near the Nyquist frequency of the original image DORG as shown in FIG. Yes. As described with reference to FIG. 7D, since the spectral intensity near the Nyquist frequency of the original image DORG is weakened by the enlargement process in the image enlargement means U1, it is weakened by the enlargement process by adding the intermediate image D1. Spectral intensity can be supplemented. Further, since the aliasing component is removed from the intermediate image D1, a false signal such as overshoot, jaggy, or ringing is not emphasized. On the other hand, the image D2B is a high frequency component corresponding to the sampling interval S1. Therefore, by adding the image D2C, a high frequency component in a band equal to or higher than the Nyquist frequency of the original image DORG can be given, so that the resolution of the image can be increased. Therefore, by adding the image D2C obtained by adding the intermediate image D1 and the image D2B to the input image DIN, it is possible to add a high-frequency component without enhancing the aliasing component, and the resolution of the image can be enhanced. It becomes.

加算手段3では画像D2Cが中間画像D2として入力画像DINに加算される。したがって、折り返し成分に起因するオーバーシュートやジャギーあるいはリンギング等の増加を抑えつつ高周波数成分を加算し、画像の解像感を高めることができる。   In the adding means 3, the image D2C is added to the input image DIN as the intermediate image D2. Therefore, it is possible to add a high frequency component while suppressing an increase in overshoot, jaggy, ringing, or the like due to the aliasing component, thereby enhancing the resolution of the image.

さらに、本発明における画像処理装置では、中間画像生成手段1および中間画像処理手段2において、画像の水平方向に関する処理、垂直方向に関する処理を並列に行っているので、画像の水平方向のみ、あるいは垂直方向のみに限らず任意の方向に関して上記の効果を得ることができる。   Further, in the image processing apparatus according to the present invention, the intermediate image generation unit 1 and the intermediate image processing unit 2 perform the processing related to the horizontal direction of the image and the processing related to the vertical direction in parallel, so that only the horizontal direction of the image or the vertical direction. The above effects can be obtained not only in the direction but also in any direction.

また、本発明における画像処理装置では周波数空間で考えて原点からFnに渡る周波数帯域のうち、原画DORGのナイキスト周波数±Fn/2近傍(あるいは特定の周波数帯域)に入力画像DINが持っている成分をもとに、ナイキスト周波数±Fn近傍の高周波数成分に対応した画像D2Bを生成している。したがって、なんらかの理由で、入力画像DINにおいて、ナイキスト周波数±Fn近傍の周波数成分が失われていたとしても、画像D2Bにより、ナイキスト周波数±Fn近傍の周波数成分を与えることが可能になる。   In the image processing apparatus according to the present invention, the component of the input image DIN in the vicinity of the Nyquist frequency ± Fn / 2 (or a specific frequency band) of the original image DORG in the frequency band extending from the origin to Fn in terms of frequency space. Based on the above, an image D2B corresponding to a high frequency component near the Nyquist frequency ± Fn is generated. Therefore, even if the frequency component near the Nyquist frequency ± Fn is lost in the input image DIN for some reason, the image D2B can provide the frequency component near the Nyquist frequency ± Fn.

なお、特定の周波数帯域として用いる箇所は、±Fn/2近傍に限定されるものではない。すなわち高周波数成分画像生成手段1Aおよび低周波数成分画像生成手段1Bの周波数応答を適宜変更することで、利用する周波数帯域を変更することができる。   The location used as the specific frequency band is not limited to the vicinity of ± Fn / 2. That is, the frequency band to be used can be changed by appropriately changing the frequency response of the high-frequency component image generating unit 1A and the low-frequency component image generating unit 1B.

上記の説明ではナイキスト周波数近傍Fn近傍の周波数成分が失われる例として画像の拡大処理を挙げたが、入力画像DINに対してナイキスト周波数Fn近傍の周波数成分が失われる原因はそれに限らず、他にもノイズ除去処理等が考えられる。したがって本発明における画像処理装置の用途は画像拡大処理後に限定されるものではない。   In the above description, the image enlargement process is given as an example in which the frequency component near the Nyquist frequency Fn is lost. However, the cause of the loss of the frequency component near the Nyquist frequency Fn with respect to the input image DIN is not limited thereto. Also, noise removal processing or the like can be considered. Therefore, the application of the image processing apparatus in the present invention is not limited after the image enlargement process.

実施の形態2.
図13は本発明の実施の形態2による画像処理方法を示すフロー図であり、本発明の実施の形態2による画像処理方法は、中間画像生成ステップST1、中間画像処理ステップST2、及び加算ステップST3により実現される。
Embodiment 2. FIG.
FIG. 13 is a flowchart showing an image processing method according to the second embodiment of the present invention. The image processing method according to the second embodiment of the present invention includes an intermediate image generation step ST1, an intermediate image processing step ST2, and an addition step ST3. It is realized by.

中間画像生成ステップST1は図14に示すように、高周波数成分画像生成ステップST1A、及び低周波数成分画像生成ステップST1Bを含む。   As shown in FIG. 14, the intermediate image generation step ST1 includes a high frequency component image generation step ST1A and a low frequency component image generation step ST1B.

高周波数成分画像生成ステップST1Aは水平方向高周波数成分画像生成ステップST1Ah、及び垂直方向高周波数成分画像生成ステップST1Avを含み、低周波数成分画像生成ステップST1Bは水平方向低周波数成分画像生成ステップST1Bh、及び垂直方向高周波巣成分画像ST1Bvを含む。   The high frequency component image generation step ST1A includes a horizontal high frequency component image generation step ST1Ah and a vertical high frequency component image generation step ST1Av, and the low frequency component image generation step ST1B includes a horizontal low frequency component image generation step ST1Bh, and The vertical high-frequency nest component image ST1Bv is included.

中間画像処理ステップST2は図15に示すように、非線形処理ステップST2A、高周波数成分画像生成ステップST2B、及び加算ステップST2Cを含む。   As shown in FIG. 15, the intermediate image processing step ST2 includes a non-linear processing step ST2A, a high frequency component image generation step ST2B, and an addition step ST2C.

非線形処理ステップST2Aは水平方向非線形処理ステップST2Ah、及び垂直方向非線形処理ステップST2Avを含み、高周波数成分画像生成ステップST2Bは水平方向高周波数成分通過ステップST2Bh、及び垂直方向高周波数成分通過ステップST2Bvを含む。   The nonlinear processing step ST2A includes a horizontal nonlinear processing step ST2Ah and a vertical nonlinear processing step ST2Av, and the high frequency component image generation step ST2B includes a horizontal high frequency component passing step ST2Bh and a vertical high frequency component passing step ST2Bv. .

水平方向非線形処理ステップST2Ahは図16に示すように、ゼロクロス判定ステップST311h、及び信号増幅ステップST312hを含み、垂直方向非線形処理ステップST2Avは図17に示すように、ゼロクロス判定ステップST311v、及び信号増幅ステップST312vを含む。   As shown in FIG. 16, the horizontal non-linear processing step ST2Ah includes a zero-cross determination step ST311h and a signal amplification step ST312h. The vertical non-linear processing step ST2Av includes a zero-cross determination step ST311v and a signal amplification step as shown in FIG. Includes ST312v.

まず、図14のフローに従って中間画像生成ステップST1の動作について説明する。
高周波数成分画像生成ステップST1Aでは、図示しない画像入力ステップにて入力された入力画像DINに対し、以下のような処理が行われる。まず、水平方向高周波数成分画像生成ステップST1Ahでは、水平方向のハイパスフィルタ処理によって、入力画像DINから水平方向の高周波数成分を取り出した画像D1Ahを生成する。垂直方向高周波数成分画像ステップST1Avでは、垂直方向のハイパスフィルタ処理によって、入力画像DINから垂直方向の高周波数成分を取り出した画像D1Avを生成する。すなわち、高周波数成分画像生成ステップST1Aは、高周波数成分画像生成手段1Aと同様の処理を行い、入力画像DINから、画像D1Ahおよび画像D1Avから成る画像D1Aを生成する。この動作は高周波数成分画像生成手段1Aと同等である。
First, the operation of the intermediate image generation step ST1 will be described according to the flow of FIG.
In the high frequency component image generation step ST1A, the following processing is performed on the input image DIN input in an image input step (not shown). First, in the horizontal high-frequency component image generation step ST1Ah, an image D1Ah obtained by extracting a horizontal high-frequency component from the input image DIN is generated by horizontal high-pass filter processing. In the vertical high-frequency component image step ST1Av, an image D1Av obtained by extracting the high-frequency component in the vertical direction from the input image DIN is generated by high-pass filtering in the vertical direction. That is, the high frequency component image generation step ST1A performs the same processing as the high frequency component image generation means 1A, and generates an image D1A composed of the image D1Ah and the image D1Av from the input image DIN. This operation is equivalent to the high frequency component image generating means 1A.

低周波数成分画像生成ステップST1Bでは、画像D1Aに対し、以下のような処理が行われる。まず、水平方向低周波数成分画像生成ステップST1Bhでは、水平方向のローパスフィルタ処理によって、画像D1Ahから水平方向の低周波数成分を取り出した画像D1Bhを生成する。垂直方向低周波数成分画像生成ステップSTBvでは、垂直方向のローパスフィルタ処理によって、画像D1Avから垂直方向の低周波数成分を取り出した画像D1Bvを生成する。すなわち、低周波数成分画像生成ステップST1Bは、低周波数成分画像生成手段1Bと同様の処理を行い、画像D1Aから、画像D1Bhおよび画像D1Bvから成る画像D1Bを生成する。この動作は低周波数成分画像生成手段1Bと同等である。   In the low frequency component image generation step ST1B, the following processing is performed on the image D1A. First, in the horizontal direction low frequency component image generation step ST1Bh, an image D1Bh obtained by extracting a horizontal low frequency component from the image D1Ah is generated by a horizontal low-pass filter process. In the vertical direction low frequency component image generation step STBv, an image D1Bv obtained by extracting the low frequency component in the vertical direction from the image D1Av is generated by the low pass filter processing in the vertical direction. That is, the low frequency component image generation step ST1B performs the same processing as the low frequency component image generation means 1B, and generates an image D1B composed of the image D1Bh and the image D1Bv from the image D1A. This operation is equivalent to the low frequency component image generation means 1B.

以上が中間画像生成ステップST1の動作であり、中間画像生成ステップST1は画像D1Bhを画像D1hとし、画像D1Bvを画像D1vとし、画像D1hおよび画像D1vから成る中間画像D1を出力する。以上の動作は中間画像生成手段1と同等である。   The operation of the intermediate image generation step ST1 is as described above. In the intermediate image generation step ST1, the image D1Bh is set as the image D1h, the image D1Bv is set as the image D1v, and the intermediate image D1 composed of the image D1h and the image D1v is output. The above operation is the same as that of the intermediate image generating unit 1.

次に図15〜17に従って中間画像処理ステップST2の動作について説明する。
まず、非線形処理ステップST2Aでは中間画像D1に対し、以下のような処理を行う。
Next, the operation of the intermediate image processing step ST2 will be described with reference to FIGS.
First, in the nonlinear processing step ST2A, the following processing is performed on the intermediate image D1.

まず、水平方向非線形処理ステップST2Ahでは、図16に示すフローに従った処理で画像D1hから画像D2Ahを生成する。図16に示すフローでの処理は以下の通りである。まず、ゼロクロス判定ステップST311hでは、画像D1hにおける画素値の変化を水平方向に沿って確認する。そして画素値が正の値から負の値あるいは負の値から正の値へと変化する箇所をゼロクロス点として捉え、ゼロクロス点の左右に位置する画素を信号増幅ステップST312hに通知する。信号増幅ステップST312hでは画像D1hについて、ゼロクロス点の左右に位置すると通知された画素の画素値を増幅し、その画像を画像D2Ahとして出力する。すなわち、非線形処理ステップST2Ahは、画像D1hに対し、水平方向非線形処理手段2Ahと同様の処理を行い、画像D2Ahを生成する。   First, in the horizontal non-linear processing step ST2Ah, an image D2Ah is generated from the image D1h by processing according to the flow shown in FIG. The processing in the flow shown in FIG. 16 is as follows. First, in the zero cross determination step ST311h, a change in pixel value in the image D1h is confirmed along the horizontal direction. Then, a portion where the pixel value changes from a positive value to a negative value or from a negative value to a positive value is regarded as a zero cross point, and the pixels located on the left and right of the zero cross point are notified to the signal amplification step ST312h. In the signal amplification step ST312h, for the image D1h, the pixel value of the pixel notified to be positioned on the left and right of the zero cross point is amplified, and the image is output as the image D2Ah. That is, the non-linear processing step ST2Ah performs the same processing as the horizontal non-linear processing means 2Ah on the image D1h to generate the image D2Ah.

次に、垂直方向非線形処理ステップST2Avでは、図17に示すフローに従った処理で画像D1vから画像D2Avを生成する。図17に示すフローでの処理は以下の通りである。まず、ゼロクロス判定ステップST311vでは、画像D1vにおける画素値の変化を垂直方向に沿って確認する。そして画素値が正の値から負の値あるいは負の値から正の値へと変化する箇所をゼロクロス点として捉え、ゼロクロス点の上下に位置する画素を信号増幅ステップST312vに通知する。信号増幅ステップST312vでは画像D1vについて、ゼロクロス点の上下に位置すると通知された画素の画素値を増幅し、その画像を画像D2Avとして出力する。すなわち、非線形処理ステップST2Avは、画像D1vに対し、垂直方向非線形処理手段2Avと同様の処理を行い、画像D2Avを生成する。   Next, in the vertical non-linear processing step ST2Av, an image D2Av is generated from the image D1v by processing according to the flow shown in FIG. Processing in the flow shown in FIG. 17 is as follows. First, in the zero cross determination step ST311v, a change in pixel value in the image D1v is confirmed along the vertical direction. Then, a point where the pixel value changes from a positive value to a negative value or from a negative value to a positive value is regarded as a zero cross point, and the pixels located above and below the zero cross point are notified to the signal amplification step ST312v. In the signal amplification step ST312v, the pixel value of the pixel notified to be positioned above and below the zero cross point is amplified for the image D1v, and the image is output as the image D2Av. That is, the nonlinear processing step ST2Av performs the same processing as the vertical nonlinear processing means 2Av on the image D1v to generate the image D2Av.

以上が非線形処理ステップST2Aの動作であり、非線形処理ステップST2Aは画像D2Ahおよび画像D2Avから成る画像D2Aを生成する。その動作は非線形処理手段2Aと同等である。   The above is the operation of the non-linear processing step ST2A, and the non-linear processing step ST2A generates an image D2A composed of the image D2Ah and the image D2Av. The operation is equivalent to the nonlinear processing means 2A.

次に、高周波数成分画像生成ステップST2Bでは画像D2Aに対し、以下の様な処理を行う。   Next, in the high frequency component image generation step ST2B, the following processing is performed on the image D2A.

まず、水平方向高周波数成分画像生成ステップST2Bhでは、画像D2Ahに対し水平方向のハイパスフィルタ処理を行った画像D2Bhを生成する。すなわち、水平方向高周波数成分画像生成ステップST2Bhは、水平方向高周波数成分画像生成手段2Bhと同様の処理を行う。   First, in the horizontal direction high frequency component image generation step ST2Bh, an image D2Bh obtained by performing a high-pass filter process in the horizontal direction on the image D2Ah is generated. That is, the horizontal direction high frequency component image generation step ST2Bh performs the same processing as the horizontal direction high frequency component image generation means 2Bh.

次に、垂直方向高周波数成分画像生成ステップST2Bvでは、画像D2Avに対し垂直方向のハイパスフィルタ処理を行った画像D2Bvを生成する。すなわち、垂直方向高周波数成分画像生成ステップST2Bvは、垂直方向高周波数成分画像生成手段2Bvと同様の処理を行う。   Next, in the vertical direction high-frequency component image generation step ST2Bv, an image D2Bv obtained by performing vertical high-pass filter processing on the image D2Av is generated. That is, the vertical high frequency component image generation step ST2Bv performs the same processing as the vertical high frequency component image generation means 2Bv.

以上が高周波数成分画像生成ステップST2Bの動作であり、高周波数成分画像生成ステップST2Bは画像D2Bhおよび画像D2Bvから成る画像D2Bを生成する。その動作は高周波数成分画像生成手段2Bと同等である。   The above is the operation of the high frequency component image generation step ST2B, and the high frequency component image generation step ST2B generates an image D2B composed of the image D2Bh and the image D2Bv. The operation is the same as that of the high frequency component image generating means 2B.

加算ステップST2Cは画像D2Bと中間画像D1を加算して画像D2Cを生成する。この際、画像D2BとD1の加算は重み付け加算でもよい。この動作は加算手段2Cと同等である。   In the addition step ST2C, the image D2B and the intermediate image D1 are added to generate an image D2C. At this time, the addition of the images D2B and D1 may be weighted addition. This operation is equivalent to the adding means 2C.

以上が中間画像処理ステップST2の動作であり、中間画像処理ステップST2は画像D2Cを中間画像D2として出力する。この動作は中間画像処理手段2と同等である。   The above is the operation of the intermediate image processing step ST2, and the intermediate image processing step ST2 outputs the image D2C as the intermediate image D2. This operation is equivalent to the intermediate image processing means 2.

加算ステップST3は入力画像DINと中間画像D2を加算し、出力画像DOUTを生成する。そして出力画像DOUTが本発明における画像処理方法の最終出力画像として出力される。すなわち、加算ステップST3の動作は加算手段3の動作と同等である。
以上が本発明における画像処理方法の動作である。
In addition step ST3, the input image DIN and the intermediate image D2 are added to generate an output image DOUT. The output image DOUT is output as the final output image of the image processing method in the present invention. That is, the operation of the adding step ST3 is equivalent to the operation of the adding means 3.
The above is the operation of the image processing method according to the present invention.

本発明における画像処理方法の動作はその説明から明らかなように、本発明の実施の形態1における画像処理装置と同等である。したがって本発明における画像処理方法は、本発明の実施の形態1における画像処理装置と同様の効果を持つ。また、上記の画像処理方法で処理された画像を、図4に示す画像表示装置で表示することもできる。   As is apparent from the description, the operation of the image processing method according to the present invention is equivalent to that of the image processing apparatus according to the first embodiment of the present invention. Therefore, the image processing method according to the present invention has the same effect as the image processing apparatus according to the first embodiment of the present invention. Further, the image processed by the above-described image processing method can be displayed by the image display device shown in FIG.

本発明の実施の形態1による画像処理装置の構成を示すブロック図である。It is a block diagram which shows the structure of the image processing apparatus by Embodiment 1 of this invention. 図1の水平方向非線形処理手段2Ahの構成を示すブロック図である。It is a block diagram which shows the structure of the horizontal direction nonlinear processing means 2Ah of FIG. 図1の垂直方向非線形処理手段2Avの構成を示すブロック図である。FIG. 2 is a block diagram showing a configuration of vertical non-linear processing means 2Av in FIG. 1. 図1の画像処理装置を用いた画像表示装置の構成を示すブロック図である。It is a block diagram which shows the structure of the image display apparatus using the image processing apparatus of FIG. 図4の画像拡大手段U1の構成を示すブロック図である。It is a block diagram which shows the structure of the image expansion means U1 of FIG. (A)〜(E)は、図4の画像拡大手段U1の動作を示す画素配置図である。(A)-(E) is a pixel arrangement | positioning figure which shows operation | movement of the image expansion means U1 of FIG. (A)〜(D)は、図4の画像拡大手段U1の動作を示す周波数応答図及び周波数スペクトル図である。(A)-(D) are the frequency response figure and frequency spectrum figure which show operation | movement of the image expansion means U1 of FIG. (A)〜(E)は、図1の中間画像生成手段1の動作を示す周波数応答図及び周波数スペクトル図である。(A)-(E) are the frequency response figure and frequency spectrum figure which show the operation | movement of the intermediate | middle image generation means 1 of FIG. (A)〜(C)は、中間画像処理手段2の動作を示す周波数応答図及び周波数スペクトル図である。。(A)-(C) are the frequency response figure and frequency spectrum figure which show operation | movement of the intermediate | middle image processing means 2. FIG. . (A)〜(C)は、ステップエッジとステップエッジをサンプリング間隔S1でサンプリングしたときに得られる、相連続する画素の信号の値を示す図である。(A)-(C) is a figure which shows the value of the signal of the pixel which continues when a step edge and a step edge are sampled by sampling interval S1. (A)〜(C)は、ステップエッジとステップエッジをサンプリング間隔S2でサンプリングしたときに得られる、相連続する画素の信号の値を示す図である。(A)-(C) is a figure which shows the value of the signal of the pixel which continues when a step edge and a step edge are sampled by sampling interval S2. (A)〜(F)は、図1の中間画像生成手段1および中間画像処理手段2の動作を示す、相連続する画素の信号の値を示す図である。(A)-(F) is a figure which shows the value of the signal of the continuous pixel which shows operation | movement of the intermediate image generation means 1 and the intermediate image processing means 2 of FIG. 本発明の実施の形態2による画像処理方法を示すフロー図である。It is a flowchart which shows the image processing method by Embodiment 2 of this invention. 図13の中間画像生成ステップST1における処理を示すフロー図である。It is a flowchart which shows the process in intermediate | middle image generation step ST1 of FIG. 図13の中間画像処理ステップST2における処理を示すフロー図である。It is a flowchart which shows the process in intermediate image process step ST2 of FIG. 図15の水平方向非線形処理ステップST2Ahにおける処理を示すフロー図である。It is a flowchart which shows the process in horizontal direction nonlinear process step ST2Ah of FIG. 図15の垂直方向非線形処理ステップST2Avにおける処理を示すフロー図である。FIG. 16 is a flowchart showing processing in the vertical nonlinear processing step ST2Av in FIG. 15.

符号の説明Explanation of symbols

1 中間画像生成手段、 2 中間画像処理手段、 3 加算手段、 Din 入力画像、 Dout 出力画像。
1 intermediate image generation means, 2 intermediate image processing means, 3 addition means, Din input image, Dout output image.

Claims (10)

入力画像の特定の周波数帯域近傍の成分を取り出した第1の中間画像を生成する中間画像生成手段と、
前記第1の中間画像をもとに第2の中間画像を生成する中間画像処理手段と、
前記入力画像と前記第2の中間画像を加算する第1の加算手段とを有する画像処理装置において、
前記中間画像処理手段は、前記第1の中間画像の画素に応じて処理の内容を変化させる非線形処理手段を含む
ことを特徴とする画像処理装置。
Intermediate image generating means for generating a first intermediate image obtained by extracting a component in the vicinity of a specific frequency band of the input image;
Intermediate image processing means for generating a second intermediate image based on the first intermediate image;
In an image processing apparatus having a first addition means for adding the input image and the second intermediate image,
The image processing apparatus, wherein the intermediate image processing means includes non-linear processing means for changing processing contents in accordance with pixels of the first intermediate image.
前記中間画像生成手段は、前記入力画像の低周波数成分及び高周波数成分を除去して、所定の周波数成分のみを取り出すことを特徴とする請求項1記載の画像処理装置。   The image processing apparatus according to claim 1, wherein the intermediate image generation unit removes only a predetermined frequency component by removing a low frequency component and a high frequency component of the input image. 前記中間画像生成手段は、
前記入力画像の高周波数成分のみを取り出した第1の高周波数成分画像を生成する第1の高周波数成分画像生成手段と、
前記第1の高周波数成分画像の低周波数成分のみを取り出す低周波数成分画像生成手段とを有する
ことを特徴とする請求項1又は2に記載の画像処理装置。
The intermediate image generating means
First high frequency component image generation means for generating a first high frequency component image obtained by extracting only the high frequency component of the input image;
The image processing apparatus according to claim 1, further comprising a low-frequency component image generation unit that extracts only a low-frequency component of the first high-frequency component image.
前記非線形処理手段は、
前記第1の中間画像の各画素値を画素に応じて変化させた増幅率で増幅した非線形処理画像を生成するものであって、
前記中間画像処理手段は、
前記非線形処理画像の高周波数成分のみを取り出した第2の高周波数成分画像を生成する第2の高周波数成分画像生成手段と、
前記第1の中間画像と前記第2の高周波数成分画像を加算する第2の加算手段とをさらに有する
ことを特徴とする請求項1又は2に記載の画像処理装置。
The nonlinear processing means includes:
Generating a non-linearly processed image obtained by amplifying each pixel value of the first intermediate image with an amplification factor changed according to the pixel,
The intermediate image processing means includes
Second high frequency component image generation means for generating a second high frequency component image obtained by extracting only the high frequency component of the nonlinear processed image;
The image processing apparatus according to claim 1, further comprising second addition means for adding the first intermediate image and the second high-frequency component image.
前記第1の高周波数成分画像生成手段は、
前記入力画像の各画素の水平方向近傍に存在する画素を用いて高周波数成分を取り出した第1の水平方向高周波数成分画像を生成する第1の水平方向高周波数成分画像生成手段と、
前記入力画像の各画素の垂直方向近傍に存在する画素を用いて高周波数成分を取り出した第1の垂直方向高周波数成分画像を生成する第1の垂直方向高周波数成分画像生成手段とを有し、
前記低周波数成分画像生成手段は、
前記第1の水平方向高周波数成分画像の低周波数成分のみを取り出した第1の水平方向中間画像を生成する水平方向低周波数成分画像生成手段と、
前記第1の垂直方向高周波数成分画像の低周波数成分のみを取り出した第1の垂直方向中間画像を生成する垂直方向低周波数成分画像生成手段と
を有することを特徴とする請求項3に記載の画像処理装置。
The first high frequency component image generation means includes:
First horizontal high-frequency component image generation means for generating a first horizontal high-frequency component image obtained by extracting a high-frequency component using pixels existing in the horizontal direction of each pixel of the input image;
First vertical high-frequency component image generating means for generating a first vertical high-frequency component image obtained by extracting a high-frequency component using pixels existing in the vertical direction of each pixel of the input image; ,
The low frequency component image generation means includes
Horizontal low-frequency component image generation means for generating a first horizontal intermediate image obtained by extracting only the low-frequency component of the first horizontal high-frequency component image;
The vertical direction low frequency component image generation means which produces | generates the 1st vertical direction intermediate image which took out only the low frequency component of the said 1st vertical direction high frequency component image is characterized by the above-mentioned. Image processing device.
前記第1の中間画像が、
前記第1の水平方向中間画像および前記第1の垂直方向中間画像から成り、
前記非線形処理手段は、
前記第1の水平方向中間画像の各画素値を画素に応じて変化させた増幅率で増幅した水平方向非線形処理画像を生成する水平方向非線形処理手段と、
前記第1の垂直方向中間画像の各画素値を画素に応じて変化させた増幅率で増幅した垂直方向非線形処理画像を生成する垂直方向非線形処理手段とを有し、
前記第2の高周波数成分画像生成手段は、
前記水平方向非線形処理画像の高周波数成分のみを取り出した第2の水平方向高周波数成分画像を生成する第2の水平方向高周波数成分画像生成手段と、
前記垂直方向非線形処理画像の高周波数成分のみを取り出した第2の垂直方向高周波数成分画像を生成する第2の垂直方向高周波数成分画像生成手段とを有する
ことを特徴とする請求項5記載の画像処理装置。
The first intermediate image is
Consisting of the first horizontal intermediate image and the first vertical intermediate image,
The nonlinear processing means includes:
Horizontal non-linear processing means for generating a horizontal non-linear processed image obtained by amplifying each pixel value of the first horizontal intermediate image with an amplification factor changed according to the pixel;
Vertical direction nonlinear processing means for generating a vertical direction nonlinear processing image obtained by amplifying each pixel value of the first vertical direction intermediate image with an amplification factor changed according to the pixel,
The second high frequency component image generating means includes:
Second horizontal high frequency component image generation means for generating a second horizontal high frequency component image obtained by extracting only the high frequency component of the horizontal nonlinear processed image;
The second vertical high-frequency component image generation means for generating a second vertical high-frequency component image obtained by extracting only the high-frequency component of the vertical non-linear processing image. Image processing device.
前記水平方向非線形処理手段は、
前記第1の水平方向中間画像の画素値が正から負もしくは負から正へと変化する箇所をゼロクロス点として判定する水平方向ゼロクロス点判定手段と、
前記水平方向ゼロクロス点判定手段の判定結果に応じて前記第1の水平方向中間画像の各画素に対する増幅率を決定する水平方向信号増幅手段とを有し、
前記垂直方向非線形処理手段は、
前記第1の垂直方向中間画像の画素値が正から負もしくは負から正へと変化する箇所をゼロクロス点として判定する垂直方向ゼロクロス点判定手段と、
前記垂直方向ゼロクロス点判定手段の判定結果に応じて前記第1の垂直方向中間画像の各画素に対する増幅率を決定する垂直方向信号増幅手段とを有する
ことを特徴とする請求項6に記載の画像処理装置。
The horizontal non-linear processing means includes:
Horizontal zero-cross point determining means for determining a point where the pixel value of the first horizontal intermediate image changes from positive to negative or from negative to positive as a zero-cross point;
Horizontal signal amplification means for determining an amplification factor for each pixel of the first horizontal intermediate image according to a determination result of the horizontal zero-cross point determination means;
The vertical nonlinear processing means includes:
Vertical zero-cross point determination means for determining a point where the pixel value of the first vertical intermediate image changes from positive to negative or from negative to positive as a zero-cross point;
7. The image according to claim 6, further comprising: a vertical direction signal amplification unit that determines an amplification factor for each pixel of the first intermediate image in the first vertical direction according to a determination result of the vertical direction zero cross point determination unit. Processing equipment.
請求項1から7までのいずれかに記載の画像処理装置を備えたことを特徴とする画像表示装置。   An image display device comprising the image processing device according to claim 1. 入力画像の特定の周波数帯域近傍の成分を取り出した第1の中間画像を生成する中間画像生成ステップと、
前記第1の中間画像をもとに第2の中間画像を生成する中間画像処理ステップと、
前記入力画像と前記第2の中間画像を加算する加算ステップとを有する画像処理方法において、
前記中間画像処理ステップは、前記第1の中間画像の画素に応じて処理の内容を変化させる非線形処理ステップを含む
ことを特徴とする画像処理方法。
An intermediate image generation step of generating a first intermediate image obtained by extracting a component in the vicinity of a specific frequency band of the input image;
An intermediate image processing step of generating a second intermediate image based on the first intermediate image;
In an image processing method including an adding step of adding the input image and the second intermediate image,
The intermediate image processing step includes a non-linear processing step of changing a processing content according to a pixel of the first intermediate image.
請求項9に記載の画像処理方法で処理された画像を表示する画像表示装置。   An image display device for displaying an image processed by the image processing method according to claim 9.
JP2008196277A 2008-07-30 2008-07-30 Image processing apparatus and method, and image display apparatus Active JP4610642B2 (en)

Priority Applications (5)

Application Number Priority Date Filing Date Title
JP2008196277A JP4610642B2 (en) 2008-07-30 2008-07-30 Image processing apparatus and method, and image display apparatus
EP09802651.1A EP2320637B1 (en) 2008-07-30 2009-07-15 Image processing device and method, and image display device
CN200980130164.3A CN102113307B (en) 2008-07-30 2009-07-15 Image processing device and method, and image display device
US13/056,540 US8498499B2 (en) 2008-07-30 2009-07-15 Image processing apparatus and method and image display apparatus
PCT/JP2009/003330 WO2010013401A1 (en) 2008-07-30 2009-07-15 Image processing device and method, and image display device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2008196277A JP4610642B2 (en) 2008-07-30 2008-07-30 Image processing apparatus and method, and image display apparatus

Publications (3)

Publication Number Publication Date
JP2010034959A true JP2010034959A (en) 2010-02-12
JP2010034959A5 JP2010034959A5 (en) 2010-06-24
JP4610642B2 JP4610642B2 (en) 2011-01-12

Family

ID=41738922

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008196277A Active JP4610642B2 (en) 2008-07-30 2008-07-30 Image processing apparatus and method, and image display apparatus

Country Status (1)

Country Link
JP (1) JP4610642B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013228536A (en) * 2012-04-25 2013-11-07 Mitsubishi Electric Corp Image processing device and method, and image display device

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05344386A (en) * 1992-06-09 1993-12-24 Matsushita Electric Ind Co Ltd Image contour emphasizing device
JPH07177386A (en) * 1993-12-21 1995-07-14 Matsushita Electric Ind Co Ltd Image control emphasizing device
JPH09224186A (en) * 1995-12-12 1997-08-26 Sony Corp Video camera and control correcting device
JP2002269558A (en) * 2001-01-02 2002-09-20 Eastman Kodak Co Method of calculating noise from digital image using color cross correlation statistics
JP2003304506A (en) * 2002-04-10 2003-10-24 Matsushita Electric Ind Co Ltd Video signal processing apparatus
JP2006252508A (en) * 2005-03-08 2006-09-21 Renei Kagi Kofun Yugenkoshi Method and apparatus for noise reduction of video signal

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05344386A (en) * 1992-06-09 1993-12-24 Matsushita Electric Ind Co Ltd Image contour emphasizing device
JPH07177386A (en) * 1993-12-21 1995-07-14 Matsushita Electric Ind Co Ltd Image control emphasizing device
JPH09224186A (en) * 1995-12-12 1997-08-26 Sony Corp Video camera and control correcting device
JP2002269558A (en) * 2001-01-02 2002-09-20 Eastman Kodak Co Method of calculating noise from digital image using color cross correlation statistics
JP2003304506A (en) * 2002-04-10 2003-10-24 Matsushita Electric Ind Co Ltd Video signal processing apparatus
JP2006252508A (en) * 2005-03-08 2006-09-21 Renei Kagi Kofun Yugenkoshi Method and apparatus for noise reduction of video signal

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013228536A (en) * 2012-04-25 2013-11-07 Mitsubishi Electric Corp Image processing device and method, and image display device

Also Published As

Publication number Publication date
JP4610642B2 (en) 2011-01-12

Similar Documents

Publication Publication Date Title
JP4994400B2 (en) Image processing apparatus and method, and image display apparatus
WO2009110328A1 (en) Image processing device and method, and image display device and method
WO2010013401A1 (en) Image processing device and method, and image display device
JP4994355B2 (en) Image processing apparatus and method, and image display apparatus
JP2010028178A (en) Image quality improving device and method
JP4610642B2 (en) Image processing apparatus and method, and image display apparatus
JP5300656B2 (en) Image processing apparatus and method, and image display apparatus
JP5247628B2 (en) Image processing apparatus and method, and image display apparatus and method
JP4958876B2 (en) Image processing apparatus and method, and image display apparatus
JP5471175B2 (en) Image processing apparatus and method, and image display apparatus
JP5247633B2 (en) Image processing apparatus and method, and image display apparatus and method
JP4994409B2 (en) Image processing apparatus and method, and image display apparatus
JP5349204B2 (en) Image processing apparatus and method, and image display apparatus and method
JP5225144B2 (en) Image processing apparatus and method, and image display apparatus
JP4994354B2 (en) Image processing apparatus and method, and image display apparatus
JP2010146444A (en) Image processor and method, and image display device
JP4912418B2 (en) Image processing apparatus and method, and image display apparatus
JP5247627B2 (en) Image processing apparatus and method, and image display apparatus
JP6242431B2 (en) Image processing apparatus and method, and image display apparatus
JP5984480B2 (en) Image processing apparatus and method, and image display apparatus
JP4994407B2 (en) Image processing apparatus and method, and image display apparatus
JP2011044925A (en) Video reproducing apparatus, apparatus and method for recording and reproducing video
JP2011034396A (en) Image reproduction apparatus and image recording and reproduction apparatus
JP2011049837A (en) Device for reproducing video, and device and method for recording/reproducing video

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20100510

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20100510

A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20100527

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20100720

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20100830

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: 20100914

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20101012

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

Free format text: PAYMENT UNTIL: 20131022

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 4610642

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250