[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

JP6644877B2 - Image processing device - Google Patents

Image processing device Download PDF

Info

Publication number
JP6644877B2
JP6644877B2 JP2018512795A JP2018512795A JP6644877B2 JP 6644877 B2 JP6644877 B2 JP 6644877B2 JP 2018512795 A JP2018512795 A JP 2018512795A JP 2018512795 A JP2018512795 A JP 2018512795A JP 6644877 B2 JP6644877 B2 JP 6644877B2
Authority
JP
Japan
Prior art keywords
pixel
pixel value
value
image
color
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.)
Active
Application number
JP2018512795A
Other languages
Japanese (ja)
Other versions
JPWO2017183273A1 (en
Inventor
忠和 角戸
忠和 角戸
恵太 中山
恵太 中山
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Kokusai Electric Inc
Original Assignee
Hitachi Kokusai Electric Inc
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 Hitachi Kokusai Electric Inc filed Critical Hitachi Kokusai Electric Inc
Publication of JPWO2017183273A1 publication Critical patent/JPWO2017183273A1/en
Application granted granted Critical
Publication of JP6644877B2 publication Critical patent/JP6644877B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N23/00Cameras or camera modules comprising electronic image sensors; Control thereof
    • H04N23/10Cameras or camera modules comprising electronic image sensors; Control thereof for generating image signals from different wavelengths
    • H04N23/12Cameras or camera modules comprising electronic image sensors; Control thereof for generating image signals from different wavelengths with one sensor only

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Color Television Image Signal Generators (AREA)

Description

本発明は、画像処理装置及び画像処理方法に係り、例えば単板式センサなどから得られるRAW画像からRGB色情報を補間する画像処理装置に関する。   The present invention relates to an image processing apparatus and an image processing method, and more particularly to an image processing apparatus that interpolates RGB color information from a RAW image obtained from a single-plate sensor or the like.

近年、ビデオ監視の分野において一般的に使用される画像の解像度は、これまでのアナログ解像度(720pix×480pix)からFull HD解像度(1920pix×1080pix)へ移行し、画質が飛躍的に向上した。今後は、さらに高解像度化する流れがある。また、安全・安心が求められる社会情勢の中で、監視分野ではカメラで撮影した映像を忠実に再現することに加えて、カメラで撮影した映像から、人・動物・物(例えば自動車や二輪車)など様々な情報を検出する画像処理技術も向上している。そうした画像処理の中でも画像の視認性や忠実性を高める技術が近年重要視されている。   In recent years, the resolution of images generally used in the field of video surveillance has shifted from the conventional analog resolution (720 pix × 480 pix) to the Full HD resolution (1920 pix × 1080 pix), and the image quality has been dramatically improved. In the future, there is a trend to further increase the resolution. In addition, in a social situation where safety and security are required, in the field of surveillance, in addition to faithfully reproducing images captured by cameras, images captured by cameras can be used to detect people, animals, objects (for example, automobiles and motorcycles). Image processing techniques for detecting various types of information have also been improved. Among such image processing, a technique for improving the visibility and fidelity of an image has recently been regarded as important.

一般に、CCDセンサやCMOSセンサには、1つのセンサでRGBの各色を取り込む単板方式と、光学的にRGBを分離した多板方式がある。カラー画像は、画素毎にRGBの色情報(R画素値、G画素値、B画素値)を使って生成されるが、単板方式のセンサでは、原色フィルタを通して値が得られることから、各画素位置で1色の色情報しか出力しない。このセンサから読み出されたままの、画素情報が欠落している画像はRAW画像と呼ばれる。このRAW画像を生成するセンサには、原色フィルタとしてBayerパターンが用いられることが多い。   In general, a CCD sensor or a CMOS sensor includes a single-plate system in which each color of RGB is captured by one sensor, and a multi-plate system in which RGB is optically separated. A color image is generated using RGB color information (R pixel value, G pixel value, B pixel value) for each pixel. However, in a single-plate sensor, a value is obtained through a primary color filter. Only one color information is output at the pixel position. An image that is read from this sensor and lacks pixel information is called a RAW image. A Bayer pattern is often used as a primary color filter in a sensor that generates this RAW image.

図1に示す様に、RAW画像はそのままでは、3つの色情報のうち2つの色情報が欠落している。そこで、欠落した2つの色情報を補うために、デモザイキング処理と呼ばれるRGBの色情報を補間する処理を行う。デモザイキング処理として、最も単純な手法は線形補間であり、補間対象の欠落している色情報を周辺の同色の色情報から補間する。すなわち、Rの色情報が欠落している場合には、周囲のR画素の平均値にて補間することになる。画像の周波数成分の少ない単一的な領域であれば、線形補間を用いて良好な補間を行うことができる。しかしながら、周波数成分の高い画像のエッジ部分が存在する領域では、偽色と呼ばれる、被写体が本来持っていない色情報を補間してしまう。   As shown in FIG. 1, two color information out of the three color information is missing when the RAW image is left as it is. Therefore, in order to compensate for the two missing pieces of color information, a process of interpolating RGB color information called demosaicing is performed. As the demosaicing process, the simplest method is linear interpolation, in which missing color information to be interpolated is interpolated from surrounding color information of the same color. That is, when the color information of R is missing, interpolation is performed using the average value of the surrounding R pixels. In a single region having a small frequency component of an image, good interpolation can be performed using linear interpolation. However, in an area where an edge portion of an image having a high frequency component exists, color information called false color, which the subject does not originally have, is interpolated.

図2を用いて、偽色について説明する。画素位置2の色情報が欠落し、それを線形補間によって補間する場合を想定する。線形補間によって画素位置1の色情報P1と画素位置3の色情報P3を用いて補間する。結果として補間された色情報P2aは、2つの色情報P1、P3の平均値となる。このとき、図示のように、例えば本来被写体が持つ色情報P2bが、補間された色情報P2aよりも大きな値である場合、補間された色情報P2aは本来被写体が持たない色情報となり、偽色が発生していることになる。この様に、偽色はエッジ付近で生じやすい。   The false color will be described with reference to FIG. Assume that color information at pixel position 2 is missing and is interpolated by linear interpolation. Interpolation is performed using the color information P1 at the pixel position 1 and the color information P3 at the pixel position 3 by linear interpolation. As a result, the interpolated color information P2a becomes an average value of the two pieces of color information P1 and P3. At this time, as shown in the drawing, for example, if the color information P2b originally possessed by the subject is a value larger than the interpolated color information P2a, the interpolated color information P2a becomes color information originally possessed by the subject and a false color Has occurred. Thus, false colors are likely to occur near the edges.

偽色を防ぐために、画像のエッジ方向に応じて補間処理の方法を変える手法があり、例えばACPI法(適応型カラーブレーン補間法)が知られる。ACPI法では、RAW画像からエッジ方向を評価し、色情報の変化の少ない方向に並んでいる画素を用いて補間を行う。これにより偽色を抑制することはできるが、限られた数個の補間方向のうち、エッジの勾配方向にできるだけ直交する一つを選んでいるに過ぎず、偽色を完全に防ぐことはできない。   In order to prevent false colors, there is a method of changing an interpolation processing method according to an edge direction of an image. For example, an ACPI method (adaptive color brain interpolation method) is known. In the ACPI method, an edge direction is evaluated from a RAW image, and interpolation is performed using pixels arranged in a direction in which color information changes little. Although false colors can be suppressed by this, one of only a limited number of interpolation directions that is as orthogonal as possible to the edge gradient direction is selected, and false colors cannot be completely prevented. .

更に単板式センサでは、偽色以外にも色モアレという現象が発生する。色モアレは、高周波成分がある被写体映像で発生する現象である。単板式のセンサの場合、RGGBの4画素を1組としてデモザイキング処理を行うが、このとき、R/Bの画素数は画素数に対して1/4の解像度しかない。そのため高周波成分に対してエイリアシングが発生し、色情報を正しく再現できない、或いは、周期的な濃淡模様が生じてしまう。   Furthermore, in the single-plate sensor, a color moiré phenomenon occurs in addition to the false color. Color moire is a phenomenon that occurs in a subject image having high frequency components. In the case of a single-plate sensor, the demosaicing process is performed with four R, G, G and B pixels as one set. At this time, the number of R / B pixels is only 1/4 the resolution of the number of pixels. For this reason, aliasing occurs with respect to the high-frequency component, and the color information cannot be correctly reproduced, or a periodic shading pattern occurs.

なお、デモザイキング処理には、これまで様々な手法が提案されており、更に、それらで生じるノイズ等を低減し画像品位を向上させる技術も知られる(例えば特許文献1及び2参照)。   Various techniques have been proposed for the demosaicing process, and techniques for reducing noise and the like generated thereby and improving image quality are known (for example, see Patent Documents 1 and 2).

特開2005−354610号公報JP 2005-354610 A 特許第4520886号公報Japanese Patent No. 4520886

Jim S Jimmy Li and Sharmil Randhawa, "Color filter array demosaicking using high-order interpolation techniques with a weighted median filter for sharp color edge preservation",IEEE Transactions on Image Processing, 2009年9月, vol.18, No.9, pp.1946-1957Jim S Jimmy Li and Sharmil Randhawa, "Color filter array demosaicking using high-order interpolation techniques with a weighted median filter for sharp color edge preservation", IEEE Transactions on Image Processing, September 2009, vol.18, No.9, pp.1946-1957 Xin Li, "Demosaicing by Successive Approximation", IEEE Transactions on Image Processing, 2005年3月, vol.14, No.3, pp.370-379Xin Li, "Demosaicing by Successive Approximation", IEEE Transactions on Image Processing, March 2005, vol.14, No.3, pp.370-379

従来のデモザイキング処理では、視覚的に目立たない程度に偽色を確実に抑圧することができなかった。なお、確実な抑圧という点では、全ての色チャネルに対して共通の空間LPF(低域通過フィルタ)を施す方法があるが、解像感が低下してしまう。最近では、超解像やディープラーニングを用いた手法も報告されているが、ごく一部の画像の領域で生じている不具合への対処としては、処理負荷が大きすぎるという問題がある。   In the conventional demosaicing process, it is not possible to reliably suppress false colors to such an extent that they are not visually noticeable. In terms of reliable suppression, there is a method of applying a common spatial LPF (low-pass filter) to all the color channels, but the resolution is reduced. Recently, techniques using super-resolution or deep learning have been reported, but there is a problem that a processing load is too large as a measure against a problem occurring in a very small part of an image area.

本発明は、このような状況に鑑みなされたもので、比較的少ない処理量で、視覚的に目立つ偽色を十分に抑圧することを目的とする。   The present invention has been made in view of such a situation, and has as its object to sufficiently suppress visually conspicuous false colors with a relatively small processing amount.

本発明は、RAW画像を取得し、デモザイキング処理を施すデモザイキング処理器と、前記デモザイキング処理器による処理結果に対して、R画素値及びB画素値に関して所定の画素領域を用いてメディアンフィルタを複数回実施するフィルタリング部とを備える。
また、前記RAW画像は、ある単位パターンが周期的に配置された(2次元上で並進対称性を有する)RGB配列を有しており、前記フィルタリング部は、前記メディアンフィルタを施す前段の処理として、各画素の色情報としてデモザイキング処理で得られたR画素値、G画素値及びB画素値に対して、それぞれの画素において、(R画素値−G画素値)、及び(B画素値−G画素値)を算出し、それら算出された値に対して前記メディアンフィルタを施し、フィリタリング後の値に対して差分に用いた前記G画素値を加算してもよい。
また、前記デモザイキング処理器は、前記デモザイキング処理として適応型カラーブレーン補間法を用いて処理してもよい。
また、デモザイキング処理器とフィルタリング部の間に一つもしくは複数種類の画質調整処理を実行する画質調整処理器を備えていても良い。ここで、画質調整処理と呼んでいるものは、γ補正のような一般的にカメラに採用される画像処理だけでなく、夜間監視映像の暗部補正など画像の特徴を大きく強調するようなものも対象とする。
The present invention provides a demosaicing processor that acquires a RAW image and performs a demosaicing process, and a median filter using a predetermined pixel area for an R pixel value and a B pixel value with respect to a processing result of the demosaicing processor. And a filtering unit for performing the filtering a plurality of times.
Further, the RAW image has an RGB array in which certain unit patterns are periodically arranged (having translational symmetry in two dimensions), and the filtering unit performs processing before the median filter is applied. For the R pixel value, the G pixel value, and the B pixel value obtained by the demosaicing process as the color information of each pixel, (R pixel value−G pixel value) and (B pixel value− G pixel value), the median filter may be applied to the calculated values, and the G pixel value used for the difference may be added to the filtered value.
In addition, the demosaicing processor may perform the demosaicing process using an adaptive color brain interpolation method.
Further, an image quality adjustment processor that executes one or more types of image quality adjustment processing may be provided between the demosaicing processor and the filtering unit. Here, what is called image quality adjustment processing includes not only image processing generally used in cameras such as γ correction, but also processing that greatly emphasizes image features such as dark area correction of night monitoring video. set to target.

本発明によると、RAW画像のデモザイキング処理により生じる色モアレ等を確実に低減することができる。   ADVANTAGE OF THE INVENTION According to this invention, the color moire etc. which arise by the demosaicing process of a RAW image can be reduced reliably.

背景技術に係る、単板方式センサ出力画素の補間プロセスを説明する図。The figure explaining the interpolation process of the single board type sensor output pixel concerning background art. 背景技術に係る、線形補間による偽色を説明する図。The figure explaining the false color by linear interpolation concerning background art. 実施形態に係る、画像処理システム1のブロック図。FIG. 1 is a block diagram of an image processing system 1 according to an embodiment. 実施形態に係る、偽色除去器14の処理を説明するフローチャート。6 is a flowchart illustrating processing of the false color remover 14 according to the embodiment. 偽色除去器14によるメディアンフィルタ処理のハードウェア実装例。7 is a hardware implementation example of median filter processing by the false color remover 14. 偽色除去器14によるメディアンフィルタ処理の別のハードウェア実装例。7 shows another hardware implementation example of the median filter processing by the false color remover 14.

次に、本発明を実施するための形態(以下、単に「実施形態」という)を、図面を参照して具体的に説明する。本実施形態では、デモザイキング処理の後段でメディアンフィルタを用いた非線形反復法によって色モアレを低減する。メディアンフィルタは任意のm×n領域の画素を用いる。   Next, embodiments for implementing the present invention (hereinafter, simply referred to as “embodiments”) will be specifically described with reference to the drawings. In the present embodiment, color moiré is reduced by a non-linear iterative method using a median filter at a stage subsequent to the demosaicing process. The median filter uses pixels in an arbitrary m × n area.

図3は、本実施形態に係る画像処理システム 1の概略構成を示すブロック図である。ここでは、主に本実施形態に特徴的な処理に着目して示している。   FIG. 3 is a block diagram illustrating a schematic configuration of the image processing system 1 according to the present embodiment. Here, the description focuses mainly on processing characteristic of the present embodiment.

画像処理システム 1は、カメラ20と画像処理装置10とを備える。カメラ20は、単板方式のセンサを有し、センサ面には上述したBayerパターンのオンチップカラーフィルタを被せている。このカメラ20は、RAW画像出力機能を有し、画像処理装置10へRAW画像を出力することができる。すなわち、RAW画像は、BayerパターンのRGB配列の画素データを有する。   The image processing system 1 includes a camera 20 and an image processing device 10. The camera 20 has a single-plate sensor, and the sensor surface is covered with the above-described Bayer pattern on-chip color filter. The camera 20 has a RAW image output function and can output a RAW image to the image processing device 10. That is, the RAW image has the pixel data of the RGB array of the Bayer pattern.

画像処理装置10は、RAW画像入力部11と、デモザイキング処理器12と、画質調整処理器13と、偽色除去器14と、画像出力部15とを備える。   The image processing device 10 includes a RAW image input unit 11, a demosaicing processor 12, an image quality adjustment processor 13, a false color remover 14, and an image output unit 15.

RAW画像入力部11は、カメラ20からRAW画像を取得する。RAW画像は、デモザイキング処理器12に渡す前に、黒レベル補正や固定パターン雑音補正等がされていることを前提とする。   The RAW image input unit 11 acquires a RAW image from the camera 20. It is assumed that the RAW image has been subjected to black level correction, fixed pattern noise correction, and the like before being passed to the demosaicing processor 12.

デモザイキング処理器12は、RAW画像入力部11からのRAW画像に対してデモザイキング処理を施し、色分離されたRGB画像を出力する。本例ではデモザイキング処理としてACPI法を用いるが、VNG(Variable Number of Gradients)法のような他の補間法を任意に適用することができる。   The demosaicing processor 12 performs a demosaicing process on the RAW image from the RAW image input unit 11 and outputs a color-separated RGB image. In this example, the ACPI method is used as the demosaicing processing, but another interpolation method such as a VNG (Variable Number of Gradients) method can be arbitrarily applied.

画質調整処理器13は、γ補正、暗部(ブラック)補正、ハイライト(ニー)補正や、トーンマッピング等の輝度領域の操作や、独立色マスキング、かすみ除去などの色相や彩度の領域の操作や、空間領域の操作や、時間領域の操作を含む、画質や視認性の向上に有用な各種の画像処理を実行する。画質調整処理器13の処理は、偽色を強調し或いは復活させる場合があるため、メディアンフィルタによる偽色除去処理は画質調整処理後に適用する方が、効果が確かであると考えられる。   The image quality adjustment processor 13 performs operations on luminance areas such as gamma correction, dark area (black) correction, highlight (knee) correction, tone mapping, and operation on hue and saturation areas such as independent color masking and haze removal. And various kinds of image processing useful for improving image quality and visibility, including an operation in a space area and an operation in a time area. Since the processing of the image quality adjustment processor 13 may enhance or restore a false color, it is considered that applying the false color removal processing by the median filter after the image quality adjustment processing is more effective.

偽色除去器14(反復フィルタリング部)は、反復法処理として、デモザイキング処理後の画像データのRGBの色情報において、R画素とB画素に対して、G画素の差分値を算出し、算出後の差分画素データに対してm×n画素のメディアンフィルタをユーザ指定回数だけ繰り返し施す。これによって、デモザイキング処理によって生じる色モアレを低減する。   The false color remover 14 (iterative filtering unit) calculates a difference value between a G pixel and an R pixel and a B pixel in the RGB color information of the image data after the demosaicing process as an iterative process. A median filter of m × n pixels is repeatedly applied to the subsequent difference pixel data by the number of times designated by the user. As a result, color moire caused by the demosaicing process is reduced.

画像出力部15は、偽色除去器14における処理後の画像データを、所定の記憶装置または表示装置に出力する。なお、画像出力部15はRGB出力(4:4:4フォーマット)の場合には、特に変換処理を行わず出力し、Y/C出力をする場合には、RGB→YC変換処理を実施して出力する。変換処理に使用する演算式はBT.709など規格化された演算式を用いる。   The image output unit 15 outputs the processed image data in the false color remover 14 to a predetermined storage device or display device. In the case of RGB output (4: 4: 4 format), the image output unit 15 performs output without performing any particular conversion process, and when performing Y / C output, performs an RGB → YC conversion process. Output. The arithmetic expression used for the conversion process uses an arithmetic expression standardized such as BT.709.

図4に偽色除去器14による反復法処理の詳細が示される。偽色除去器14は、ユーザの指示を取得し、反復を行う回数を設定する(S10)。既に設定されている場合は、この設定処理はスキップされる。   FIG. 4 shows the details of the iterative processing by the false color remover 14. The false color remover 14 acquires the user's instruction and sets the number of times of repetition (S10). If it has already been set, this setting process is skipped.

次に、偽色除去器14は、デモザイキング処理後のRGB画像の各画素位置において、R画素値とB画素値に対して、G画素値との差分値を算出する(S11)。すなわち、例えば、ある画素の画素値X(R画素値、G画素値、B画素値)がXi(Ri, Gi, Bi)とした場合、RG差分画素データΔXi(Ri-Gi, Bi-Gi)を算出する。Next, the false color remover 14 calculates a difference value between the R pixel value and the B pixel value and the G pixel value at each pixel position of the RGB image after the demosaicing process (S11). That is, for example, when the pixel value X (R pixel value, G pixel value, B pixel value) of a certain pixel is X i (R i , G i , B i ), RG difference pixel data ΔX i (R i − G i , B i -G i ) are calculated.

次に、偽色除去器14は、算出後のRG差分画素データΔXi(Ri-Gi, Bi-Gi)で表される画像データに対して、m×n画素のメディアンフィルタ処理を実施する(S12)。なお、メディアンフィルタ処理とは処理対象画素の周辺m×n画素(ウィンドウ)のデータの中央値を算出し、算出した中央値を処理結果として採用する処理である。このとき、画像端では周辺画素がm×n存在しないケースが発生する。その場合には、存在しない画素は中央値算出の対象から除外して処理を行う。もしくは、パディング処理として画像端の画素値をm×nになるようにコピーして処理を行う。パディング処理の有無によって、中央値算出の結果は変わらない。Next, the false color remover 14 performs a median filter process of m × n pixels on the image data represented by the calculated RG difference pixel data ΔX i (R i -G i , B i -G i ). Is performed (S12). The median filter process is a process of calculating a median of data of m × n pixels (window) around a pixel to be processed, and employing the calculated median as a processing result. At this time, there is a case where m × n peripheral pixels do not exist at the image end. In this case, non-existing pixels are excluded from the median value calculation and the processing is performed. Alternatively, as a padding process, a process is performed by copying the pixel value at the end of the image so as to be m × n. The result of the median calculation does not change depending on whether or not padding is performed.

その後、偽色除去器14は、メディアンフィルタ処理後の各画素データΔXf(Rf, Bf)に、S11で減算したG画素値(Gi)を加算し(S13)、加算処理後のR、Bの画素値を処理結果として保持する(S14)。すなわち、(Ri, Gi, Bi)を、(Rf+Gi, Gi, Bf+Gi)で上書きする。なお、G画素値については、処理を行わず、入力される画素値をそのまま保持する。Thereafter, the false color remover 14 adds the G pixel value (Gi) subtracted in S11 to each pixel data ΔXf (Rf, Bf) after the median filter processing (S13), The pixel value is held as a processing result (S14). That is, (R i , G i , B i ) is overwritten with (R f + G i , G i , B f + G i ). The G pixel value is not processed and the input pixel value is held as it is.

ここで、反復回数が設定回数に達していない場合(S15のN)、S11の処理に戻る。このとき、S11以降の処理では、前回処理で得られたR画素値、G画素値及びB画素値が用いられる。反復回数が設定回数に達した場合(S15のY)、反復法処理は終了する。なお、S11からS14の各ステップの処理は、画像内の全画素について行われる。RAW画像のときから存在した画素値については、処理の対象から除外することもできるが、効果に大きな違いは無い。   If the number of repetitions has not reached the set number (N of S15), the process returns to S11. At this time, in the processing after S11, the R pixel value, the G pixel value, and the B pixel value obtained in the previous processing are used. When the number of repetitions has reached the set number (Y in S15), the iterative method processing ends. Note that the processing of each step from S11 to S14 is performed for all pixels in the image. Pixel values that have existed since the time of the RAW image can be excluded from processing, but there is no significant difference in the effect.

図5に、偽色除去器14によるメディアンフィルタ処理のハードウェア実装例が示される。図5の構成は、図4で示した反復処理の内のS11からS13についての、R画素値への操作に相当し、B画素値についても同様である。なおメディアンフィルタのウィンドウサイズは3×3としてある。また複数回の反復処理を行うには、図5の構成を単にカスケード接続すればよい。   FIG. 5 shows a hardware implementation example of the median filter processing by the false color remover 14. The configuration in FIG. 5 corresponds to the operation for the R pixel values from S11 to S13 in the repetitive processing shown in FIG. 4, and the same applies to the B pixel values. The window size of the median filter is 3 × 3. Further, in order to perform the repetitive processing a plurality of times, the configuration in FIG. 5 may be simply cascaded.

最初に、j番目の反復におけるR画素値(Ri j)、G画素値(Gi j)が減算器40に入力され、減算器40はそれらの差分を算出し、フリップフロップ(FF)31と1H遅延線41へ出力される。フリップフロップ(FF)31〜39は、それぞれ1つの画素値を保持するレジスタであり、メディアンの算出元となる3×3画素分の値を保持する。1H遅延線41、42は、1水平ライン分の画素値を保持するデジタル遅延線であり、先入れ先出し(FIFO)メモリ等で構成され得る。減算器40からの差分は負になり得るので、減算前よりもビット深度が1つ増える。First, R pixel values in the j-th iteration (R i j), G pixel value (G i j) is input to the subtracter 40, the subtracter 40 calculates the difference between them, the flip-flop (FF) 31 Is output to the 1H delay line 41. Each of the flip-flops (FF) 31 to 39 is a register that holds one pixel value, and holds a value of 3 × 3 pixels from which a median is calculated. The 1H delay lines 41 and 42 are digital delay lines that hold pixel values for one horizontal line, and can be configured by a first-in first-out (FIFO) memory or the like. Since the difference from the subtractor 40 can be negative, the bit depth is increased by one compared to before the subtraction.

メディアンフィルタ43は、入力された9個の画素値の中央値を取り出して、フィルタされた画素値(Rf j)として出力する。メディアンフィルタ43は、例えば、バブルソートのアルゴリズムをハードコーディングしたもので、隣同士で値を比較できる9個並列のレジスタを9段設けて、パイプライン動作させることで実現されうる。The median filter 43 extracts the median value of the nine input pixel values and outputs the median value as a filtered pixel value (R f j ). The median filter 43 is, for example, a hard-coded algorithm of a bubble sort, and can be realized by providing nine stages of parallel registers that can compare values between adjacent units and performing a pipeline operation.

2H+2タップ遅延線44は、画素値(Ri j)が減算器40に入力されてから、対応するフィルタされた画素値(Rf j)がメディアンフィルタ43から出力されるまでの遅延と同じ遅延を、G画素値(Gi j)に施す。加算器45は、フィルタされた画素値(Rf j)と遅延された画素値(Ri j)を加算して、出力する。2H + 2 tapped delay line 44, from the pixel values (R i j) is input to a subtracter 40, a delay until the corresponding filtered pixel value (R f j) is output from the median filter 43 The same delay is applied to the G pixel value (G i j ). The adder 45 adds the filtered pixel value (R f j) and delayed pixel values (R i j), and outputs.

図5の構成をFPGAで実現する場合、1H遅延線41、42や2H+2タップ遅延線44は、FPGA内部の組み込みブロックRAM(EBR)を用いて構成されうる。1つのEBRは18Kbit、36Kbit程度の容量を有し、外部RAMに比較して非常に高い転送帯域幅が利用可能である。それ以外のメディアンフィルタ43等は、FPGAの基本的な単位であるProgrammable Function Unit(PFU)或いはスライスを組合せて構成されうる。この様に、本例の偽色除去器14は乗算器(DSPスライス)を必要としないため、FPGA内部の余ったリソースを利用して、容易に組み込むことができる場合がある。アスペクト比が高い、つまり横長の画像では、その両側付近で横方向の倍率色収差が目立つようになる。単一の種類のウィンドウを用いてそのような収差による偽色をも補正しようとする場合、アスペクト比に合せてウィンドウも横長にした方がよい場合がある。   When the configuration of FIG. 5 is implemented by an FPGA, the 1H delay lines 41 and 42 and the 2H + 2 tap delay line 44 can be configured using an embedded block RAM (EBR) inside the FPGA. One EBR has a capacity of about 18 Kbit or 36 Kbit, and can use a very high transfer bandwidth as compared with an external RAM. Other median filters 43 and the like can be configured by combining Programmable Function Units (PFU) or slices, which are basic units of the FPGA. As described above, since the false color remover 14 of the present example does not require a multiplier (DSP slice), it may be possible to easily incorporate the false color remover 14 using extra resources inside the FPGA. In an image having a high aspect ratio, that is, a horizontally long image, lateral chromatic aberration becomes conspicuous near both sides thereof. When trying to correct even the false color due to such aberration using a single type of window, it may be better to make the window horizontally long in accordance with the aspect ratio.

図4及び図5で説明した偽色除去器14による、色モアレ低減効果の実測値を表1に示す。

Figure 0006644877
ここでは、客観的評価値を導入している。評価結果には、CIELAB表色系の色差式を用いて評価しており、差分値(ΔEab)が小さいほど本来の被写体に近い画素値が得られていることを意味する。反復処理を複数回実施することで効果があることが客観的評価からも確認できる。具体的には、デモザイキング処理(ACPI)のみでは、差分値(ΔEab)は1.94である。デモザイキング処理を行った結果に対して反復処理を1回のみ実施した結果、差分値(ΔEab)は1.3095となり、色モアレの低減効果が確認された。そしてさらに反復処理を3回実施した結果、差分値(ΔEab)は1.2663となり、複数回実施することで、さらに色モアレ低減の効果を確認することができた。Table 1 shows the measured values of the color moiré reduction effect by the false color remover 14 described with reference to FIGS.
Figure 0006644877
Here, an objective evaluation value is introduced. The evaluation result is evaluated using a color difference formula of the CIELAB color system, and a smaller difference value (ΔEab) means that a pixel value closer to the original subject is obtained. It can be confirmed from an objective evaluation that the effect is obtained by performing the repetitive processing a plurality of times. Specifically, only in the demosaicing process (ACPI), the difference value (ΔEab) is 1.94. As a result of performing the iterative process only once on the result of the demosaicing process, the difference value (ΔEab) was 1.3095, and the effect of reducing color moiré was confirmed. Further, as a result of performing the repetition processing three more times, the difference value (ΔEab) was 1.2663. By performing the repetition processing a plurality of times, the effect of reducing the color moiré could be further confirmed.

図6に、偽色除去器14によるメディアンフィルタ処理の別のハードウェア実装例が示される。この例は、メディアンフィルタ43に代えて、スイッチングメディアンフィルタ48を備えた点などにおいて、図5と異なる。なお図5と同一の構成には同一の符号を付して、説明を省略する。この例では、フィルタウィンドウの中心に対応する画素の値を出力するFF35の出力が取り出され、比較器46と、減算器49に与えられる。   FIG. 6 shows another example of hardware implementation of the median filter processing by the false color remover 14. This example differs from FIG. 5 in that a switching median filter 48 is provided instead of the median filter 43. The same components as those in FIG. 5 are denoted by the same reference numerals, and description thereof will be omitted. In this example, the output of the FF 35 that outputs the value of the pixel corresponding to the center of the filter window is extracted and provided to the comparator 46 and the subtractor 49.

比較器46は、中心画素値を、所定の閾値と比較し、その結果のブール論理値を出力する。画素アドレス生成器47は、図6のハードウェア全体の動作タイミングを支配しており、FF35から出力されている中心画素値が、画像中のどこに位置するかを示すアドレスを出力することができる。   Comparator 46 compares the center pixel value with a predetermined threshold and outputs the resulting Boolean logic value. The pixel address generator 47 controls the operation timing of the entire hardware in FIG. 6 and can output an address indicating where the center pixel value output from the FF 35 is located in the image.

本例のスイッチングメディアンフィルタ(SMF)48は、比較器46や画素アドレス生成器47から与えられた条件などによって、動作を切替えることができる順序統計フィルタ或いはスタックフィルタである。一例として、中心画素値が閾値より小さい場合、SMF48はα-トリムド平均値フィルタとして動作し、そこでは、ソートされたN個(N=m×n。本例では9個)の画素値の列における両端のαN個が破棄され、中心寄りの(1-α)N個の画素値の算術平均が出力される。また、中心画素が、Bayer配列中に存在した画素である(つまり補間によって生成された画素ではない)場合、メディアンフィルタなしで中心画素値がそのまま出力される。その他の場合は、通常のメディアンフィルタとして動作する。トリムド平均値フィルタを用いた場合、相対的に感度が低いRやBチャンネルにおいて、低輝度領域でのS/Nを改善し、色の階調を緻密に再現することができる。また、Bayer配列からの原信号を維持することで、画像のエッジを保存できる。   The switching median filter (SMF) 48 of the present example is an order statistical filter or a stack filter whose operation can be switched according to conditions given from the comparator 46 and the pixel address generator 47 and the like. As an example, if the central pixel value is less than the threshold, SMF48 operates as an α-trimmed average filter, where a sorted N (N = m × n, in this example, 9) column of pixel values Are discarded, and the arithmetic average of (1-α) N pixel values near the center is output. When the center pixel is a pixel existing in the Bayer array (that is, not a pixel generated by interpolation), the center pixel value is output without a median filter. In other cases, it operates as a normal median filter. When the trimmed average filter is used, in the R and B channels having relatively low sensitivity, the S / N in the low luminance area can be improved, and the color gradation can be precisely reproduced. Also, by maintaining the original signal from the Bayer array, the edges of the image can be preserved.

この例では、Bayer配列からの原信号はそのまま出力するとしたが、フィルタの出力と比較して原信号の方が大きい場合は、原信号を破棄してフィルタの出力を用いるようにしてもよい。なお、コンディショナルなフィルタ動作は上述の組合せに限られず、条件として、雑音や高周波成分の評価結果を用いてもよく、フィルタとしてDouble Window Modified Trimmed Meanフィルタ、スタックフィルタ、εフィルタなどを用いてもよい。   In this example, the original signal from the Bayer array is output as it is. However, if the original signal is larger than the output of the filter, the original signal may be discarded and the output of the filter may be used. Note that the conditional filter operation is not limited to the above-described combinations, and as a condition, an evaluation result of noise or a high-frequency component may be used, and a Double Window Modified Trimmed Mean filter, a stack filter, an ε filter, etc. Good.

減算器49は、SMF48からのフィルタされた画素値から、FF35からの中心画素値を減算し、修正量として乗算器50へ出力する。乗算器50は、修正量に所定の固定定数を乗算し、G修正量として減算器51へ出力する。減算器51は、2H+2タップ遅延線44からのG画素値から、乗算器50からのG修正量を減算して、出力する。これらの構成は、R画素値の修正による輝度の変化を補償するために設けてあり、必須ではない。補償の強さは、所定の固定定数を適切に選ぶことによって決定される。固定定数を乗算する乗算器50は、PFUやスライスによって構成されうる。   The subtracter 49 subtracts the center pixel value from the FF 35 from the filtered pixel value from the SMF 48, and outputs the result to the multiplier 50 as a correction amount. The multiplier 50 multiplies the correction amount by a predetermined fixed constant and outputs the result to the subtractor 51 as a G correction amount. The subtractor 51 subtracts the G correction amount from the multiplier 50 from the G pixel value from the 2H + 2 tap delay line 44 and outputs the result. These configurations are provided to compensate for a change in luminance due to the correction of the R pixel value, and are not essential. The strength of the compensation is determined by appropriately choosing a predetermined fixed constant. The multiplier 50 that multiplies the fixed constant can be configured by a PFU or a slice.

本例では、平均化を伴うフィルタを用いたことにより、小数を持つ画素値が得られ、その固定小数点の値を持つ画素値を偽色除去器14内部の反復計算の間、使用することができ、更に偽色除去器14の出力とすることができる。例えば色空間の変換の際に、固定小数点の画素値を、BT.709 の変換係数と乗算することができる。   In this example, by using a filter with averaging, a pixel value having a decimal value is obtained, and the pixel value having the fixed-point value can be used during the iterative calculation inside the false color remover 14. And the output of the false color remover 14. For example, during color space conversion, a fixed-point pixel value can be multiplied by a BT.709 conversion coefficient.

以上説明した様に、本実施形態の画像処理システム 1によると、メディアンフィルタのような非線形フィルタを反復処理することで色モアレ低減効果を得ることができる。反復処理は、必ずしも同一処理を連続して行うものに限らず、複数回の反復処理の一部を、画質調整処理器13の処理の前に行ってもよく、或いは画像出力部14からの出力された後にパソコン等によるソフトウェア処理で行ってもよい。また、オンチップカラーフィルタは、ベイヤーフィルタに限らず、4色を1組とする補色フィルタや、ベイヤーフィルタのG画素の1つを全透過或いは赤外透過フィルタに置換した場合であっても、補間によってG画素値(もしくは擬似輝度値)を全画素位置で生成すれば適用できる。   As described above, according to the image processing system 1 of the present embodiment, a color moiré reduction effect can be obtained by repeatedly processing a non-linear filter such as a median filter. The repetitive processing is not limited to performing the same processing continuously, and a part of the plurality of repetitive processing may be performed before the processing of the image quality adjustment processor 13 or the output from the image output unit 14 may be performed. After that, the processing may be performed by software processing using a personal computer or the like. In addition, the on-chip color filter is not limited to the Bayer filter, even when a complementary color filter having four colors as one set, or when one of the G pixels of the Bayer filter is replaced with a total transmission or infrared transmission filter, This is applicable if G pixel values (or pseudo luminance values) are generated at all pixel positions by interpolation.

1 画像処理システム、 10 画像処理装置、 11 RAW画像入力部、 12 デモザイキング処理器、 13 画質調整処理器、 14 偽色除去器(反復フィルタリング部)、 15 画像出力部、 20 カメラ 1 Image processing system, 10 Image processing device, 11 RAW image input unit, 12 Demosaicing processor, 13 Image quality adjustment processor, 14 False color remover (iterative filtering unit), 15 Image output unit, 20 camera

本発明は、ディジタルカメラ、ビデオカメラ等の機器に組みこまれて利用されるほか、RAW現像もしくはポストプロダクション等で用いられる画像編集ソフトウェアや、画像を保管・編集するクラウドサービスにおいても利用することができる。   The present invention can be used by being incorporated into devices such as digital cameras and video cameras, as well as image editing software used in RAW development or post-production, and cloud services for storing and editing images. it can.

Claims (5)

2次元上で並進対称性を有するカラーフィルタを設けた撮像素子から取得したRAW画像に、デモザイキング処理を施すデモザイキング処理器と、
前記デモザイキング処理器による処理結果に対して、各R画素値及びB画素値に関して所定のウィンドウを用いてメディアンフィルタを複数回実施する偽色除去器と、
前記偽色除去器は、前記メディアンフィルタを施す前段の処理として、各画素の色情報としてデモザイキング処理で得られたR画素値、G画素値及びB画素値に対して、それぞれの画素において、(R画素値−G画素値)、及び(B画素値−G画素値)を算出し、それら算出された値に対して前記メディアンフィルタを施し、フィリタリング後の値に対して差分に用いた前記G画素値を加算するものであって、
さらに、前記偽色除去器によって複数回実施されるメディアンフィルタは、スイッチングメディアンフィルタであって、前記ウィンドウの中心画素の値が所定の閾値より小さい場合、α−トリムド平均値フィルタとして動作し、或いは、前記ウィンドウの中心画素が、Bayer配列中に存在した画素である場合、前記中心画素の値をそのまま出力する動作をし、
前記メディアンフィルタが複数回実施される間、画素値のビット深度は前記偽色除去器に入力される前に比べて増加することを特徴とする画像処理装置。
A demosaicing processor for performing a demosaicing process on a RAW image obtained from an image sensor provided with a color filter having translational symmetry in two dimensions;
A false color remover that performs a median filter a plurality of times using a predetermined window for each R pixel value and B pixel value with respect to the processing result by the demosaicing processor;
The false color eliminator performs, as a process prior to applying the median filter, an R pixel value, a G pixel value, and a B pixel value obtained in the demosaicing process as color information of each pixel. (R pixel value−G pixel value) and (B pixel value−G pixel value) were calculated, the calculated value was subjected to the median filter, and the value after filtering was used as a difference. G pixel values are added,
Further, the median filter performed by the false color remover a plurality of times is a switching median filter, and operates as an α-trimmed average filter when a value of a center pixel of the window is smaller than a predetermined threshold value, or When the center pixel of the window is a pixel existing in the Bayer array, the operation of outputting the value of the center pixel as it is is performed,
The image processing apparatus according to claim 1, wherein the bit depth of a pixel value is increased as compared to before being input to the false color remover, while the median filter is performed a plurality of times.
前記RAW画像は、BayerパターンのRGB配列を有しており、
前記デモザイキング処理器は、前記デモザイキング処理として適応型カラーブレーン補間法を用いて処理するものであり、
前記デモザイキング処理器による処理結果に対し、前記偽色除去器に渡される前に、画質や視認性を向上させる、少なくとも輝度領域もしくは色領域の画像処理を行う画質調整処理器を備えることを特徴とする請求項1に記載の画像処理装置。
The RAW image has an RGB array of a Bayer pattern,
The demosaicing processor performs processing using an adaptive color brain interpolation method as the demosaicing processing,
Before passing the processing result by the demosaicing processing unit to the false color elimination unit , the image processing unit includes an image quality adjustment processing unit that performs image processing of at least a luminance region or a color region to improve image quality and visibility. The image processing apparatus according to claim 1.
前記ウィンドウは、前記RAW画像のアスペクト比に対応させて横長の矩形状に設定されることを特徴とする請求項1または2に記載の画像処理装置。   The image processing apparatus according to claim 1, wherein the window is set in a horizontally long rectangular shape corresponding to an aspect ratio of the RAW image. 2次元上で並進対称性を有するカラーフィルタを設けた撮像素子から取得したRAW画像に、デモザイキング処理を施すデモザイキング処理手段と、
前記デモザイキング処理手段による処理結果に対して、各R画素値及びB画素値に関して所定のウィンドウを用いてメディアンフィルタを複数回実施する偽色除去手段と、
前記偽色除去手段は、前記メディアンフィルタを施す前段の処理として、各画素の色情報としてデモザイキング処理で得られたR画素値、G画素値及びB画素値に対して、それぞれの画素において、(R画素値−G画素値)、及び(B画素値−G画素値)を算出し、それら算出された値に対して前記メディアンフィルタを施し、フィリタリング後の値に対して差分に用いた前記G画素値を加算するものであって、
さらに、前記偽色除去手段によって複数回実施されるメディアンフィルタは、前記ウィンドウの中心画素の値が所定の閾値より小さい場合、α−トリムド平均値フィルタとして動作し、或いは、前記ウィンドウの中心画素が、Bayer配列中に存在した画素である場合、前記中心画素の値をそのまま出力する動作をし、
前記メディアンフィルタが複数回実施される間、画素値のビット深度は前記偽色除去手段に入力される前に比べて増加することを特徴とする画像処理装置における画像処理方法。
Demosaicing processing means for performing demosaicing processing on a RAW image obtained from an image sensor provided with a color filter having translational symmetry in two dimensions;
A false color removing unit that performs a median filter a plurality of times using a predetermined window for each of the R pixel value and the B pixel value with respect to a processing result obtained by the demosaicing processing unit;
The false color removing unit performs, as a process prior to applying the median filter, an R pixel value, a G pixel value, and a B pixel value obtained by demosaicing as color information of each pixel. (R pixel value−G pixel value) and (B pixel value−G pixel value) were calculated, the calculated value was subjected to the median filter, and the value after filtering was used as a difference. G pixel values are added,
Further, the median filter executed multiple times by the false color removing means operates as an α-trimmed average filter when the value of the center pixel of the window is smaller than a predetermined threshold, or the center pixel of the window is , When the pixel exists in the Bayer array, the operation of outputting the value of the central pixel as it is is performed,
An image processing method according to an image processing apparatus, wherein the bit depth of a pixel value increases before the median filter is executed a plurality of times as compared to before being input to the false color removing unit.
前記ウィンドウは、前記RAW画像のアスペクト比に対応させて横長の矩形状に設定されることを特徴とする請求項4に記載の画像処理方法。
The image processing method according to claim 4, wherein the window is set in a horizontally long rectangular shape corresponding to an aspect ratio of the RAW image.
JP2018512795A 2016-04-20 2017-02-16 Image processing device Active JP6644877B2 (en)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2016084343 2016-04-20
JP2016084343 2016-04-20
PCT/JP2017/005696 WO2017183273A1 (en) 2016-04-20 2017-02-16 Image processing device

Publications (2)

Publication Number Publication Date
JPWO2017183273A1 JPWO2017183273A1 (en) 2019-01-17
JP6644877B2 true JP6644877B2 (en) 2020-02-12

Family

ID=60115985

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018512795A Active JP6644877B2 (en) 2016-04-20 2017-02-16 Image processing device

Country Status (2)

Country Link
JP (1) JP6644877B2 (en)
WO (1) WO2017183273A1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US12136183B2 (en) 2020-10-23 2024-11-05 Samsung Electronics Co., Ltd. Demosaicing method and demosaicing device

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113840124B (en) * 2021-10-12 2023-08-18 锐芯微电子股份有限公司 Image processing method and system
JP7299398B1 (en) 2022-10-14 2023-06-27 株式会社 日立産業制御ソリューションズ Image processing device and image processing method

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000217123A (en) * 1999-01-20 2000-08-04 Canon Inc Image pickup device and image processing method for image pickup device
JP2001022923A (en) * 1999-07-08 2001-01-26 Sanyo Electric Co Ltd Image data processor
WO2006043414A1 (en) * 2004-10-19 2006-04-27 Megachips Lsi Solutions Inc. Image processing method and image processing device
JP2009077051A (en) * 2007-09-19 2009-04-09 Toshiba Corp Imaging apparatus and imaging method thereof

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US12136183B2 (en) 2020-10-23 2024-11-05 Samsung Electronics Co., Ltd. Demosaicing method and demosaicing device

Also Published As

Publication number Publication date
JPWO2017183273A1 (en) 2019-01-17
WO2017183273A1 (en) 2017-10-26

Similar Documents

Publication Publication Date Title
US9013611B1 (en) Method and device for generating a digital image based upon a selected set of chrominance groups
US8295595B2 (en) Generating full color images by demosaicing noise removed pixels from images
KR100786931B1 (en) Image signal processing apparatus and image signal processing method
KR101327789B1 (en) Method and apparatus for reducing various noises of image simultaneously
JP4290193B2 (en) Image processing device
US8773593B2 (en) Noise reduction filter circuitry and method
KR20110029001A (en) Image porcessing method and apparatus in the bayer image
JP2012256202A (en) Image processing apparatus and method, and program
US8253829B2 (en) Image processing apparatus, imaging apparatus, and image processing method
JP4952251B2 (en) Image processing method, image processing program, recording medium, and image processing apparatus
JP2009303139A (en) Solid-state imaging apparatus
WO2008119228A1 (en) Low complexity color de-noising filter
JP6644877B2 (en) Image processing device
JP3309941B2 (en) Noise detection circuit, noise elimination circuit and contour enhancement circuit
JP2011151739A (en) Image processing apparatus
US8385671B1 (en) Digital camera and method
US20140037207A1 (en) System and a method of adaptively suppressing false-color artifacts
US8049788B2 (en) Color difference correction and imaging device
JP2007323635A (en) Recursive filtering of video image
CN113129389B (en) Method for judging moire, method for inhibiting moire and circuit system
JP3633561B2 (en) Image processing apparatus, image processing method, and program
KR20080035440A (en) Image processing apparatus and method
JP4380399B2 (en) Imaging apparatus, noise reduction apparatus, noise reduction method, and program
US9007494B2 (en) Image processing apparatus, method for controlling the same and storage medium
JP5103580B2 (en) Image processing apparatus and digital camera

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180911

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20180911

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20191016

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20191210

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200108

R150 Certificate of patent or registration of utility model

Ref document number: 6644877

Country of ref document: JP

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