JP2006344002A - 画像処理装置、エッジ検出装置、画像処理方法、エッジ検出方法及びそのプログラム - Google Patents

画像処理装置、エッジ検出装置、画像処理方法、エッジ検出方法及びそのプログラム Download PDF

Info

Publication number
JP2006344002A
JP2006344002A JP2005169011A JP2005169011A JP2006344002A JP 2006344002 A JP2006344002 A JP 2006344002A JP 2005169011 A JP2005169011 A JP 2005169011A JP 2005169011 A JP2005169011 A JP 2005169011A JP 2006344002 A JP2006344002 A JP 2006344002A
Authority
JP
Japan
Prior art keywords
image
region
area
enlarged
edge
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
JP2005169011A
Other languages
English (en)
Other versions
JP4716255B2 (ja
Inventor
Satoshi Kubota
聡 久保田
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.)
Fujifilm Business Innovation Corp
Original Assignee
Fuji Xerox Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fuji Xerox Co Ltd filed Critical Fuji Xerox Co Ltd
Priority to JP2005169011A priority Critical patent/JP4716255B2/ja
Publication of JP2006344002A publication Critical patent/JP2006344002A/ja
Application granted granted Critical
Publication of JP4716255B2 publication Critical patent/JP4716255B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/40Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/403Edge-driven scaling; Edge-based scaling

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)
  • Editing Of Facsimile Originals (AREA)

Abstract

【課題】 軽い処理負荷で、高画質な拡大画像を高速に得ることができる画像処理装置を提供する。
【解決手段】 画像拡大プログラム4の画像ブロック特徴量算出部420は、画像ブロック中の注目領域のエッジ強度を算出し、この注目領域が特徴ブロックであるか非特徴ブロックであるかを判定し、ブロック重畳度を算出する。高画質画像ブロック生成部430は、注目領域及びその周辺領域の画像データに対して、コントラスト強調処理を施し、注目領域のエッジパターン及び注目領域の推定エッジ方向等に基づいて、注目領域に対する拡大画像ブロックを生成する。高速拡大処理部440は、入力画像データに対して拡大処理を行う。拡大画像生成部450は、算出されたブロック重畳度に基づいて、高画質画像ブロック生成部430により生成された拡大画像と、高速拡大処理部440により拡大された拡大画像を統合する。
【選択図】図19

Description

本発明は、画像を拡大する画像処理装置に関する。
画像の拡大処理は、画像の編集やファイリング、表示、印刷などを行うシステムにとって、基本的な処理の一つである。また近年、インターネットのホームページ上の画像やデジタルビデオなどのディスプレイ解像度での表示を主目的とした画像データなどの普及により、これらの低解像度画像を高解像度のプリンタなどで印刷することも頻繁に行われている。このプリンタによる印刷の際に、高画質の出力結果を得ることが望まれており、高画質の拡大処理に対する重要度が高まっている。
カラーを含む多階調で表現された画像(以下では、これを多値画像と称する)を拡大処理する既存の代表的な手法としては、最近傍法や線形補間法、キュービック・コンボリューション法などがある。最近傍法は、拡大後の各画素値として、その画素を原画像上に逆写像した際に最も距離が近い画素の画素値を使う方法である。この方法は、演算量が少ないため、高速に処理することができる。しかし、原画像の1画素がそのまま矩形形状に拡大されるため、隣り合う画素の画素値の差が小さい場合は、画質劣化の程度は小さくほとんど影響はないが、逆に大きい場合などは、斜線部やエッジ部のジャギーが目立ったり、倍率が大きい場合には、画像がモザイク状になるなど、画質劣化の程度は大きい。
線形補間法は、画素間の画素値が直線的に変化していると仮定し、拡大後の画素を逆写像した点の近傍4画素の画素値を線形に補間して画素値を求める方法である。この方法では、最近傍法よりも処理は重いが、演算量は比較的少なく、ジャギーなども発生しにくい。その一方で、直線的に変化しているという仮定に当てはまらないエッジ部分を中心に、画像全体がボケ気味になるという欠点がある。
キュービック・コンボリューション法は、標本化定理に基づいてsinc関数(sin(x)/x)を近似した補間関数を定義し、拡大後の画素を逆写像した点の近傍16画素(X、Y方向それぞれ4画素)と前記の近似補間関数との畳み込み演算により、拡大後の画素値を求める方法である。この方法では、前記2つの手法に比べて画質は比較的良いが、参照範囲が大きく、演算量が多いという欠点がある。また、高域が強調気味となる特性を持つため、エッジ部分で軽いジャギーが発生したり、ノイズ成分が強調されてしまうなどの欠点もある。
これら拡大画像の画質問題を解決する試みとして、例えば、特許文献1、特許文献2などの新規方式が提案されている。
特許文献1に記載されている技術は、原画像の変化度合いの検出法として、エッジ検出フィルタによるエッジ検出を行い、そのエッジ検出結果に基づいてエッジ画素を定義する。そして、この方法は、エッジ画素と判断された場合には、キュービック・コンボリューション法の3次関数形状を調整したM−キュービック法で拡大を行い、そうでない場合には、最近傍法で拡大する。しかし、変化度合いの大きいエッジ部分を、キュービック・コンボリューション法の3次関数形状よりもエッジを強調するように調整したM−キュービック法で行うので、エッジ部でジャギーが発生したり、エッジ部周辺のオーバーシュートあるいはアンダーシュートが増大するという欠点がある。
特許文献2に記載されている技術は、DCTの正変換、逆変換を繰り返し行いながら高周波成分を復元して拡大する手法である。さらに、この手法は、平坦部分に生じるリンギングを抑制するために、平坦部分領域においては、原画像をバイリニア法やキュービック法で補間拡大した画像を参照画像とし、参照画像の対応する領域との誤差を許容偏差に制限する処理を行う。しかし、本特許は、DCTの正変換、逆変換を繰り返し行う反復法を基本としており、画像中の大部分を占める平坦部のリンギングを抑制するための参照処理など、演算量が非常に多い。
また、画像のエッジ部分のアンダーシュート、オーバーシュートを抑制する試みとして、特許文献3の方式が提案されている。特許文献3に記載されている技術は、エッジ部を示す指標として分散値を用いて、分散値がある閾値より大きいエッジ部では、強調係数がより小さくなるように補正して高域強調処理を行うことにより、過度の強調(アンダーシュートやオーバーシュートの発生)を抑制する。しかし、この手法を画像の拡大処理に用いた場合、オーバーシュート、アンダーシュートはある程度抑制されうるが、エッジ部分がボケることは防がれない。
特開2000−188681号公報 特許第3404138号 特開平9−93439号公報
本発明は、上述した背景からなされたものであり、画像を高画質に拡大する画像処理装置を提供することを目的とする。
上記目的を達成するために、本発明にかかる画像処理装置は、入力画像の拡大処理を行う画像処理装置であって、注目画素を含む所定の大きさの画像領域の特徴量を算出する領域特徴量算出手段と、前記領域特徴量算出手段により算出された画像領域の特徴量に基づいて、この画像領域に対応する拡大画像領域を生成する拡大画像領域生成手段と、前記拡大画像領域生成手段とは異なる拡大手法を適用して、入力画像の拡大画像を生成する拡大処理手段と、前記領域特徴量算出手段により算出された特徴量が所定の条件を満たす画像領域が重畳する場合、その重畳する度合いに基づいて、前記拡大画像領域生成手段により生成された拡大画像領域の生成値を補正し、この拡大画像領域と、前記拡大処理手段により生成された拡大画像とに基づいて、入力画像に対する拡大画像を生成する拡大画像生成手段とを有する。
好適には、前記領域特徴量算出手段は、画像領域内の各画素値に基づいて、画像領域における階調変化量を算出し、この階調変化量を画像領域の特徴量とする。
好適には、前記領域特徴量算出手段は、画像領域内の各画素値に基づいて、画像領域における階調変化方向を算出し、この階調変化方向を画像領域の特徴量とする。
好適には、前記領域特徴量算出手段は、画像領域内の各画素値に基づいて、画像領域における階調変化を算出し、この階調変化が既定の画素値パターンに対応する場合に、この画素値パターンを画像領域の特徴量として選択する。
好適には、前記領域特徴量算出手段は、色空間における色成分それぞれについて特徴量を算出し、算出された特徴量に基づいて1つの色成分を選択し、この色成分について算出された特徴量を画像領域の特徴量とする。
好適には、前記拡大画像領域生成手段は、前記領域特徴量算出手段により算出された画像領域の特徴量が所定の条件を満たす場合に、画像領域の特徴量と、この画像領域の近傍領域の特徴量と、この近傍領域内の画素値とを用いて、拡大画像領域を生成する。
好適には、前記拡大画像領域生成手段は、階調変化を強調するエッジ強調処理を拡大率に応じて行って、画像領域内の画素値を補正し、補正された画素値を用いて拡大画像領域を生成する。
好適には、前記拡大画像領域生成手段は、前記領域特徴量算出手段により算出された階調変化方向に応じて、画像領域の近傍領域内の画素値を選択し、選択された画素値を用いて拡大画像領域を生成する。
好適には、前記拡大画像領域生成手段は、前記領域特徴量算出手段により算出された特徴量に対応する既定の演算式を用いて、拡大画像領域を生成する。
好適には、前記拡大処理手段は、少なくとも最近傍補間法に基づいた拡大手法を適用する。
好適には、前記拡大処理手段は、少なくとも線形補間法に基づいた拡大手法を適用する。
好適には、前記拡大処理手段は、入力画像毎あるいは入力画像中の数ライン毎に、拡大画像を生成する。
好適には、前記拡大画像生成手段は、前記拡大画像領域生成手段により生成された拡大画像領域を、前記拡大処理手段により生成された拡大画像上の対応する位置に配置する。
好適には、前記拡大画像生成手段は、前記拡大画像領域生成手段により生成された前記拡大画像領域の生成値を、前記拡大処理手段により生成された拡大画像上の対応する位置の画素値を用いて補正する。
好適には、前記拡大画像生成手段は、前記拡大画像領域生成手段により生成された拡大画像領域の生成値と、前記拡大処理手段により生成された拡大画像上の対応する位置の画素値とを、所定の割合で線形和した値を用いて、この画素値を補正する。
好適には、前記拡大画像生成手段は、複数の画像領域において前記領域特徴量算出手段により算出された画像領域の特徴量が所定の条件を満たし、これらの画像領域が互いに重畳する場合、これらの画像領域が重畳する度合いに基づいて、線形和の割合を決定する。
好適には、前記拡大画像生成手段は、前記領域特徴量算出手段により算出された画像領域の特徴量と、この画像領域と重畳する画像領域の特徴量とに基づいて、入力画像のエッジ部及びエッジ周辺部を判定する。
好適には、前記拡大画像生成手段は、複数の画像領域において前記領域特徴量算出手段により算出された画像領域の特徴量が所定の条件を満たし、これらの画像領域が互いに重畳する場合、これらの画像領域が重畳する度合いが所定値以上である画像部分を、エッジ部と判定する。
好適には、前記拡大画像生成手段は、複数の画像領域において前記領域特徴量算出手段により算出された画像領域の特徴量が所定の条件を満たし、これらの画像領域が互いに重畳する場合、これらの画像領域が重畳する度合いが所定値以下である画像部分を、エッジ周辺部と判定する。
また、好適には、前記拡大画像生成手段は、前記拡大画像領域生成手段により生成された拡大画像領域で、エッジ部に対応する画像部分の生成値を、前記拡大処理手段により生成された拡大画像の対応する位置に配置する。
また、好適には、前記拡大画像生成手段は、前記拡大画像領域生成手段により生成された拡大画像領域で、エッジ周辺部に対応する画像部分の生成値と、前記拡大処理手段により生成された拡大画像の対応する位置の画素値とを用いて、この画素値を補正する。
さらに、好適には、前記拡大画像生成手段は、前記拡大画像領域生成手段により生成された拡大画像領域で、エッジ部及びエッジ周辺部のいずれにも対応しない画像部分には、前記拡大処理手段により生成された拡大画像上の対応する位置の画素値を用いる。
また、本発明にかかるエッジ検出措置は、入力画像のエッジを検出するエッジ検出装置であって、注目画素を含む所定の大きさの画像領域の特徴量を算出する領域特徴量算出手段と、前記領域特徴量算出手段により算出された画像領域の特徴量と、この画像領域と重畳する画像領域の特徴量とに基づいて、入力画像のエッジ部及びエッジ周辺部を判定する判定手段とを有する。
好適には、前記判定手段は、複数の画像領域において前記領域特徴量算出手段により算出された画像領域の特徴量が所定の条件を満たし、これらの画像領域が互いに重畳する場合、これらの画像領域が重畳する度合いが所定値以上である画像部分を、エッジ部と判定する。
好適には、前記判定手段は、複数の画像領域において前記領域特徴量算出手段により算出された画像領域の特徴量が所定の条件を満たし、これらの画像領域が互いに重畳する場合、これらの画像領域が重畳する度合いが所定値以下である画像部分を、エッジ周辺部と判定する。
また、本発明にかかる画像処理方法は、入力画像の拡大処理を行う画像処理方法であって、注目画素を含む所定の大きさの画像領域の特徴量を算出し、前記算出された画像領域の特徴量に基づいて、この画像領域に対応する拡大画像領域を生成し、前記生成とは異なる拡大手法を適用して、入力画像の拡大画像を生成し、前記算出された特徴量が所定の条件を満たす画像領域が重畳する場合、その重畳する度合いに基づいて、前記生成された拡大画像領域の生成値を補正し、この拡大画像領域と、前記生成された拡大画像とに基づいて、入力画像に対する拡大画像を生成する。
また、本発明にかかるエッジ検出方法は、入力画像のエッジを検出するエッジ検出方法であって、注目画素を含む所定の大きさの画像領域の特徴量を算出し、前記算出された画像領域の特徴量と、この画像領域と重畳する画像領域の特徴量とに基づいて、入力画像のエッジ部及びエッジ周辺部を判定する。
また、本発明にかかる第1のプログラムは、コンピュータを含み、入力画像の拡大処理を行う画像処理装置において、注目画素を含む所定の大きさの画像領域の特徴量を算出する領域特徴量算出ステップと、前記算出された画像領域の特徴量に基づいて、この画像領域に対応する拡大画像領域を生成する拡大画像領域生成ステップと、前記生成とは異なる拡大手法を適用して、入力画像の拡大画像を生成する拡大処理ステップと、前記算出された特徴量が所定の条件を満たす画像領域が重畳する場合、その重畳する度合いに基づいて、前記生成された拡大画像領域の生成値を補正し、この拡大画像領域と、前記生成された拡大画像とに基づいて、入力画像に対する拡大画像を生成する拡大画像生成ステップとを前記画像処理装置のコンピュータに実行させる。
さらに、本発明にかかる第2のプログラムは、コンピュータを含む入力画像のエッジを検出するエッジ検出装置でおいて、注目画素を含む所定の大きさの画像領域の特徴量を算出する領域特徴量算出ステップと、前記算出された画像領域の特徴量と、この画像領域と重畳する画像領域の特徴量とに基づいて、入力画像のエッジ部及びエッジ周辺部を判定する判定ステップとを前記エッジ検出装置のコンピュータに実行させる。
本発明の画像処理装置によれば、軽い処理負荷で、高画質な拡大画像を高速に得ることができる。
[ハードウェア構成]
本実施形態における画像処理装置2のハードウェア構成を説明する。
図1は、本発明にかかる画像処理方法が適応される画像処理装置2のハードウェア構成を、制御装置20を中心に例示する図である。
図1に例示するように、画像処理装置2は、CPU202及びメモリ204などを含む制御装置20、通信装置22、HDD・CD装置などの記憶装置24、並びに、LCD表示装置あるいはCRT表示装置及びキーボード・タッチパネルなどを含むユーザインターフェース装置(UI装置)26を有する。
画像処理装置2は、例えば、画像拡大プログラム4(図2を参照して後述)がインストールされた汎用コンピュータであり、通信装置22又は記憶装置24などを介して画像データを取得し、取得された画像データを出力解像度に応じて拡大する。例えば、画像処理装置2は、プリンタ装置10に対して画像データを出力する場合には、600dpi又は2400dpiなどの解像度に変換し、UI装置26に対して画像データを出力する場合には、75dpiなどの解像度に変換する。
[画像拡大プログラム]
図2は、制御装置20(図1)により実行され、本発明にかかる画像処理方法を実現する画像拡大プログラム4の機能構成を示す図である。
図2に示すように、画像拡大プログラム4は、記憶部400、画像ブロック設定部410、画像ブロック特徴量算出部420、高画質画像ブロック生成部430、高速拡大処理部440及び拡大画像生成部450を有する。また、画像ブロック特徴量算出部420は、エッジ強度算出部422、エッジ方向推定部424及びエッジパターン選択部426を有し、拡大画像生成部450は、ブロック重畳度判定部452及び出力値補正部454を有する。
なお、画像拡大プログラム4の全部又は一部は、ASICなどのハードウェアにより実現されてもよい。
記憶部400は、画像データが拡大処理されるまで画像データを一時的に記憶し、解像度変換又は拡大処理された拡大画像データが出力装置(不図示)に出力されるまで一時的に記憶する。なお、画像データは、画像処理装置2により処理可能である画像フォーマット(例えば、BMP、TIFF、PNGなど)で記述されたデータであり、デジタルカメラ(不図示)又はスキャナ(不図示)などにより取り込まれた画像データ、又は、パーソナルコンピュータ(画像処理装置2など)等において作成又は編集等を行うアプリケーションプログラムにより作成された画像データ等である。拡大画像データ(拡大された後の画像データ)もまた、同様の画像フォーマットのデータである。
画像ブロック設定部410は、画像ブロック特徴量算出部420及び高画質画像ブロック生成部430の処理において必要とされる既定の画像ブロックサイズをそれぞれ設定する。また、画像ブロック設定部410は、記憶部400により記憶されている入力画像データから、設定されたブロックサイズの画像ブロックを順次(例えばラスタスキャン順に)切り出し、このブロックサイズの画像ブロックそれぞれを、画像ブロック特徴量算出部420及び高画質画像ブロック生成部430それぞれに対して出力する。
画像ブロック特徴量算出部420は、画像ブロック設定部410から順次入力される画像ブロックの少なくとも一部である注目領域における画像特徴量を、注目領域又はこの注目領域の周辺部を含む画像ブロック内の各画素値に基づいて算出する。画像特徴量は、例えば、注目領域のエッジ強度(階調変化量)又はエッジ角度(階調変化の方向)である。しかし、画像特徴量はこれらに限られない。画像ブロック特徴量算出部420は、例えば、注目領域の各画素値の平均値を算出し、この平均値に対する注目領域の各画素値のばらつきを表す値(例えば標準偏差又は分散)を画像特徴量として算出してもよい。
また、画像ブロック特徴量算出部420は、それぞれの注目領域について算出された画像特徴量と基準値とを比較し、注目領域の切り分けを行う。注目領域の切り分けにより、特徴のある画像ブロック(例えば、エッジを含んだ画像ブロックなど)と、特徴の少ない画像ブロック(例えば、画素値変化の少ない画像ブロックなど)とが切り分けられる。より具体的には、画像ブロック特徴量算出部420は、算出された画像特徴量が基準値以上である画像ブロックを特徴ブロックと判定し、算出された画像特徴量が基準値未満である画像ブロックを非特徴ブロックと判定する。
本例では、画像ブロック特徴量算出部420は、画像特徴量の1つとして画像ブロックのエッジ強度(数値)を算出し、算出されたエッジ強度が基準値以上である画像ブロックを特徴ブロックであると判定し、エッジ強度が基準値未満である画像ブロックを非特徴ブロックであると判定する。なお、特徴ブロックの切り分けに用いられる特徴量は、エッジ強度に限られず、例えば注目領域の平均値に対する注目領域の各画素値のばらつきを表す値である標準偏差、分散又はその他の特徴量であってもよい。
高画質画像ブロック生成部430は、画像ブロック特徴量算出部420により算出された画像特徴量を用いて、注目領域に対応する拡大画像ブロックを生成する。高画質画像ブロック生成部430における拡大処理には、注目領域に含まれる特徴を保持する拡大手法が適用されることが望ましい。また、高画質画像ブロック生成部430は、画像ブロック特徴量算出部420により算出された特徴量が基準値以上である画像ブロック(すなわち、特徴ブロック)についてのみ、拡大処理を行う。
高速拡大処理部440は、記憶部400により記憶されている入力画像データを拡大する。高速拡大処理部440は、高画質画像ブロック生成部430による拡大処理よりも処理負荷の小さい拡大アルゴリズムを適用する。具体的には、高速拡大処理部440は、例えば、最近傍補間拡大法や線形補間拡大法などを適用して、入力画像データを拡大する。また、高速拡大処理部440は、画像ブロック毎の処理ではなく、入力画像単位又は入力画像の数ライン単位で拡大処理を行う。
拡大画像生成部450は、高画質画像ブロック生成部430により拡大された拡大画像ブロックと、高速拡大処理部440により拡大された拡大画像とを用いて、原画像に対する拡大画像を生成する。より具体的には、拡大画像生成部450は、画像ブロック特徴量算出部420により算出された画像特徴量が基準値以上である画像ブロック(特徴ブロック)が重なり合う場合には、その重なり数に従い、高画質画像ブロック生成部430により拡大された拡大画像ブロックの各々の画素値を、高速拡大処理部440により拡大された拡大画像の相当する位置の画素値に基づいて補正することにより、拡大画像を生成する。
[画像ブロック特徴量算出部]
次に、画像ブロック特徴量算出部420を、より詳細に説明する。なお、注目領域が2×2画素サイズブロックであり、注目領域を含む周辺領域が4×4画素サイズブロックである場合を具体例として説明する。
エッジ強度算出部422は、画像ブロック設定部410により切り出された画像ブロックの中の注目領域のエッジ強度Gを、次の式(1)により算出する。
gx=(a+c−b−d)/2
gy=(a+b−c−d)/2
G=gx×gx+gy×gy・・・(1)
式(1)において、a、b、c及びdは、図6に例示するように、注目領域にある各画素の画素値である。gxは、主走査方向(図6の左右方向)の画素値の変化量を示し、gyは、副走査方向(図6の上下方向)の画素値の変化量を示す。
なお、エッジ強度は、上記式(1)のみにより算出される値でなく、以下の式(2)などにより算出されてもよい。
G=|gx|+|gy|・・・(2)
すなわち、エッジ強度Gは、gxの絶対値と、gyの絶対値との和として算出されてもよい。
また、エッジ強度算出部422は、算出されたエッジ強度Gと基準値(既定のしきい値Th)との比較を行い、注目領域が特徴ブロック又は非特徴ブロックのいずれであるかを判定し、判定結果を高画質画像ブロック生成部430等に対して出力する。
なお、注目領域のエッジ強度Gがしきい値Thよりも小さい場合は、エッジ方向推定部424及びエッジパターン選択部426は、この注目領域について処理を行わない。
図3は、エッジ強度算出部422において、エッジ強度Gと基準値との比較により特徴ブロックと判定された領域を例示する。
図4は、互いに重なり合う特徴ブロックを例示する。
図3に例示するように、特徴ブロックは、原画像中のエッジ付近に集中する。さらに、画像ブロック特徴量算出部420が、注目領域の画像特徴量を、1画素づつ主走査方向及び副走査方向にずらしながら算出すると、図4に例示するように、原画像中のエッジ付近においては、特徴ブロックは互いに重なり合う。
図5は、特徴ブロック同士の重なり数(ブロック重畳度)を例示する。
図5に例示するように、ブロック重畳度の範囲は、例えば1〜4であり、エッジ強度算出部422は、このブロック重畳度も算出し、算出結果を高画質画像ブロック生成部430及び拡大画像生成部450に対して出力する。
ブロック重畳度が算出されると、原画像中のエッジ付近において、エッジ領域とエッジ周辺領域とが区別される。より具体的には、ブロック重畳度が小さい領域は、エッジ周辺領域であると区別され、ブロック重畳度が大きい領域は、エッジ領域であると区別される。
エッジ方向推定部424について説明する。
図6は、注目領域及び周辺領域の具体例と、注目領域のエッジ方向の一例の説明図であり、図6(A)は、注目領域及び周辺領域を例示し、図6(B)は、この注目領域について推定されるエッジ方向を例示する図である。
図6(A)に例示するように、注目領域(画像領域)は、2×2の矩形領域(主走査方向及び副走査方向にそれぞれ2つ)を有し、周辺領域は、4×4の矩形領域(主走査方向及び副走査方向にそれぞれ4つ)を有する。各矩形は、それぞれ画素に相当し、矩形内の各数字はそれぞれの画素値を示している。すなわち、注目領域は、中心部近傍の画素{a,b,c,d}={15,104,86,203}である。
以下、この図6(A)で例示する注目領域を具体例として、エッジ方向推定部424によるエッジ方向推定処理を説明する。
図7は、エッジ方向推定部424によるエッジ方向推定処理(S10)を示すフローチャートである。
図8は、エッジ方向推定処理(S10)のS102の処理において選択される参照領域を例示する図である。なお、図8におけるハッチングされた部分は、図6に示した注目領域に相当する。
図7に示すように、ステップ100(S100)において、エッジ方向推定部424は、図6(A)に例示した注目領域のエッジ角度Θを、次の式(3)で計算する。
Θ=arctan(gy/gx)・・・(3)
図6(A)では、注目領域の画素値は、{a,b,c,d}={15,104,86,203}であり、式(1)より、
gx=−103
gy=−85
となり、これらを式(3)に代入すると、
Θ=−140.5°
となる。
このエッジ角度Θの方向は、図6(B)に示された破線方向に相当する。
さらに、エッジ方向推定部424は、算出されたエッジ角度Θが22.5°ごとに区分された方向(8方向)の角度範囲のいずれに含まれるかを判定する。本例では、エッジ角度Θが0°又は±180°を中心とした角度範囲を「方向0」とし、22.5°又は−157.5°を中心とした角度範囲を「方向1」とし、45°又は−135°を中心とした角度範囲を「方向2」とし、67.5°又は−112.5°を中心とした角度範囲を「方向3」とし、90°又は−90°を中心とした角度範囲を「方向4」とし、112.5°又は−67.5°を中心とした角度範囲を「方向5」とし、135°又は−45°を中心とした角度範囲を「方向6」とし、157.5°又は−22.5°を中心とした角度範囲を「方向7」する。これらの角度範囲は、それぞれの中心から±11.25°の範囲である。上述の具体例におけるエッジ角度Θ(=−140.5°)は、−135°±11.25°の範囲内に含まれるので、エッジ角度は「方向2」となる。
ステップ102(S102)において、エッジ方向推定部424は、算出された注目領域のエッジ角度Θに応じて、図6(A)に示した周辺領域(太線枠外の領域)の中からエッジ方向の推定に用いる参照領域を選択する。より具体的には、エッジ方向推定部424は、算出されたエッジ角度Θの方向で注目領域と隣接する可能性のある画素を含むように、参照領域を選択する。
例えば、エッジ方向推定部424は、注目領域について算出されたエッジ角度が「方向0」に含まれる場合に、図8(A)に例示する参照領域(太線で囲まれた2つの領域)を選択し、算出されたエッジ角度が「方向4」に含まれる場合に、図8(B)に例示する参照領域(太線で囲まれた2つの領域)を選択し、算出されたエッジ角度が上記以外の方向(方向1〜3、方向5〜7)に含まれる場合に、図8(C)に例示する参照領域(太線で囲まれた4つの領域)を選択する。図6に示した具体例では、エッジ角度の方向は「方向2」であるので、図8(C)に示した4つの参照領域が選択の候補となる。
なお、参照領域は、図8に例示した領域に限定されず、例えば図8(C)の場合には、参照領域数は8であったり、それぞれの方向に応じて設定されてもよい。
ステップ104(S104)において、エッジ方向推定部424は、選択された参照領域それぞれに対して、S100の処理と同様に、式(1)及び式(3)に従って、エッジ角度Θを計算する。
ステップ106(S106)において、エッジ方向推定部424は、それぞれの参照領域について算出されたエッジ角度と、注目領域について算出されたエッジ角度とを比較して、これらの差分が予め設定されている閾値Θthより小さいか否かを判断する。エッジ方向推定部424は、エッジ角度の差分が閾値Θthより小さい場合には、この参照領域のエッジ角度を適正なエッジ角度として判定してS108の処理に移行し、エッジ角度の差分が閾値Θth以上である場合には、この参照領域のエッジ角度を適正なエッジ角度ではないと判定してS110の処理に移行する。
ステップ108(S108)において、エッジ方向推定部424は、角度参照数をインクリメントする。すなわち、エッジ方向推定部424は、参照領域について算出されたエッジ角度が適正なエッジ角度であると判断された場合にのみ、角度参照数をインクリメントする。
なお、角度参照数は、エッジ角度の参照数をカウントするための変数であり、注目領域ごとに「角度参照数1」に初期化される。
ステップ110(S110)において、エッジ方向推定部424は、選択した全ての参照領域についてエッジ角度を算出したか否かを判断し、全ての参照領域についてエッジ角度が算出された場合には、S112の処理に移行し、これ以外の場合には、S104の処理に戻って次の参照領域についてエッジ角度を算出する。
ステップ112(S112)において、エッジ方向推定部424は、注目領域のエッジ角度と、適正なエッジ角度として判定された参照領域のエッジ角度との総和を計算し、算出されたエッジ角度の総和を角度参照数で割った平均エッジ角度を注目領域の推定エッジ方向とする。
なお、図6に示した具体例では、エッジ角度ΘUは、上部の参照領域{86,203,171,211}からΘU=−149.4°となり、エッジ角度ΘLは、左部の参照領域{10,15,20,86}からΘL=−131.2°となり、エッジ角度ΘDは、下部の参照領域{1,102,15,104}からΘD=−175.2°となり、エッジ角度ΘRは、右部の参照領域{104,215,203,219}からΘR=−141.0°となる。注目領域のエッジ角度Θ=−140.5°とそれぞれの参照領域のエッジ角度とが比較され、その差分が閾値Θthより小さい参照領域の数が、角度参照数としてカウントされる。
図9は、図6に示した注目領域における推定エッジ方向の一例の説明図である。例えば上述の具体例において、すべての参照領域について注目領域のエッジ角度との差分が閾値Θthより小さいとすれば、注目領域及び4つの参照領域から求められたエッジ角度の総和は−737.3°となり、角度参照数5で割ることによって、平均エッジ角度は−147.5°と求められる。この場合においても、エッジ方向推定部434は、上述の注目領域のエッジ方向と同様に、平均エッジ角度が例えば8方向のいずれかに含まれるかを判定する。本例では、平均エッジ角度が−147.5°であるため「方向1」に含まれ、「方向1」が推定エッジ方向となる。
なお、本実施形態では、1画素につき1色要素であるグレースケール画像を具体例として説明しているが、これに限定されるわけではない。例えば、1画素につき3色要素のRGB色空間のカラー画像が入力される場合には、各々の色成分のデータにおけるエッジ強度Gr、Gg、Gbの強さにより選択された色空間データにより、上記のエッジ方向推定処理が行われればよい。より具体的には、画像ブロック特徴量算出部420は、それぞれの色成分についてエッジ強度を算出し、算出されたエッジ強度Gr、Gg、Gbが最大となる色成分を選択して、選択された色成分についてのみ特徴量を算出する。このようにして、カラー画像における拡大画像データのエッジ部の色ずれなど、画質低下が抑えられることが可能となる。
また、本例では、注目領域及び参照領域について、式(1)により算出されたエッジ角度は、8方向のいずれかに分類されたが、分類は8方向に限定されず、より精度の高いエッジ方向が必要であれば、分類は、12方向(15.0°ごと)、16方向(12.25°ごと)など、さらに多数の角度領域を対象としてなされてもよい。
エッジパターン選択部426について説明する。
図10は、エッジパターン選択部426により用いられるエッジパターンテーブルを例示する図である。
図10に例示するように、エッジパターン選択部426は、推定エッジ方向とエッジパターンとを互いに対応付けたエッジパターンテーブルを有する。エッジパターンテーブルには、エッジパターンが、注目領域のパターンサイズに対応して、推定エッジ方向(例えば8方向)ごとに1つ以上登録されている。
エッジパターン選択部426は、このエッジパターンテーブルを参照して、エッジ方向推定部424によりそれぞれの注目領域について推定された推定エッジ方向に対応するエッジパターンを選択する。
本例では、図9に例示するように、注目領域に対する推定エッジ方向が、エッジ方向推定部424により「方向1」であると推定されているため、エッジパターン選択部426は、この推定エッジ方向(方向1)に従い、図10に示すエッジパターンテーブルの中から、方向1に対応する「パターン0」から「パターン3」までの4つのエッジパターンを選択し、これらをこの注目領域に対するエッジパターンの候補とする。
エッジパターン選択部426は、注目領域の画素値に基づいて、エッジパターンの候補となった1つ以上のエッジパターンの中から、1つのエッジパターンを選択する。エッジパターンの具体的な選択方法について、図11を参照しながら説明する。
図11は、図10に示した注目領域に対応するエッジパターンの選択方法を説明する図である。
本例では、推定エッジ方向が「方向1」であったため、図11(A)に例示するように、エッジパターンの候補として、「パターン0」から「パターン3」までの4つのエッジパターンが選択されている。これらのエッジパターンは、図11(B)に例示するように、ビットパターンとして表現される。具体的には、エッジパターンは、白部分を0、それ以外を1としてビットパターン化され、「ビットパターン0」から「ビットパターン3」までのビットパターンが生成される。なお、これらのビットパターンは、図10に示すエッジパターンテーブルにビットテーブルとして予め登録されていてもよい。
エッジパターン選択部426は、注目領域に相当するビットパターンを判定する。具体的には、エッジパターン選択部426は、以下に示す式(4)に従い、注目領域中の平均画素値を計算し、注目領域内の各々の画素値から平均値を引き、その符号を以て注目領域の画素値パターンとする。
Mean=(a+b+c+d)/4
a_sign=a−Mean
b_sign=b−Mean
c_sign=c−Mean
d_sign=d−Mean・・・(4)
なお、本例では、図11(C)に示すように、Mean=(15+104+86+203)/4=102であり、a_sign=−87、b_sign=2、c_sign=−16、d_sign=101となる。よって、エッジパターン選択部426は、これらの正負符号を判定して、図11(C)に示すように、注目領域のビットパターン(1010)を生成する。
エッジパターン選択部426は、図11(B)に例示するエッジパターン候補に対応するビットパターンと、図11(C)に例示する注目領域のビットパターンとのパターンマッチングを行い、最も類似するエッジパターンを選択パターンとして決定する。選択されたエッジパターンは、後述する高画質画像ブロック生成部430における拡大画像ブロック生成処理に適用される。
なお、エッジパターンは、図10に示したエッジパターンに限定されず、例えば、エッジパターン選択部426は、入力画像データの種類に応じて、エッジパターンテーブルを切り替えて、異なるエッジパターンを適用してもよい。また、エッジパターン選択部426は、各角度におけるエッジパターン候補数を増減させてもよい。
[高画質画像ブロック生成部]
次に、高画質画像ブロック生成部430をより詳細に説明する。
高画質画像ブロック生成部430は、画像ブロック特徴量算出部420において、特徴ブロックであると判定された注目領域に対して、画像ブロック特徴量算出部420により得られた注目領域に対するエッジパターン及び推定エッジ方向に基づいて、拡大画像ブロック生成処理を行う。
高画質画像ブロック生成部430は、まず、拡大処理の拡大倍率に応じたサイズ及び係数の強調カーネルを用いて、画像ブロック設定部410により切り出された画像ブロック中の注目領域およびその周辺領域の画像データのコントラストを強調する。
図12は、高画質画像ブロック生成部430により用いられる強調カーネル432(エッジ強調カーネル)を例示する図である。
図12に例示するように、第1の強調カーネル432aは、重み付け係数「1.60」及び「−0.15」を用いてコントラストを強調し、第2の強調カーネル432bは、重み付け係数「1.20」及び「−0.05」を用いてコントラストを強調する。これらの強調カーネルは、対象画像に対して既になされた拡大処理の拡大倍率に対応付けられており、互いに異なる重み付け係数を用いて、互いに異なる位置の画素値を参照する。
第1の強調カーネル432aは、図12(A)に例示するように、注目画素Pの直下画素a、直右画素b、直上画素c及び直左画素dを参照して、これらの画素の画素値にそれぞれ重み付け係数(−0.15)を掛け合わせ、重み付け係数(1.60)が掛けられた注目画素Pの画素値と合算し、合算された値を注目画素Pの画素値とする。
第2の強調カーネル432bは、第1の強調カーネル432aよりも拡大倍率の大きな画像に対して適用され、図12(B)に例示するように、注目画素Pから1画素分離間した下画素a、右画素b、上画素c及び左画素dを参照して、これらの画素の画素値にそれぞれ重み付け係数(−0.05)を掛け合わせ、重み付け係数(1.20)が掛けられた注目画素Pの画素値と合算し、合算された値を注目画素Pの画素値とする。
例えば、第1の強調カーネル432aを適用する場合に、以下の式(5)に従って、コントラスト強調後の画素値P’が算出される。
画素値P’=1.60×P−0.15×(a+b+c+d)・・・(5)
このように、強調カーネルは、既になされた拡大処理の拡大倍率に応じて異なるので、画像が、4倍、8倍と順に拡大されると、原画像の特徴を有している画素は、2画素、4画像離れた位置の画素となる。そのため、高画質画像ブロック生成部430は、図12に例示するように、拡大倍率が高いほど離れた画素を参照して強調処理を行う。例えば、高画質画像ブロック生成部430は、2倍拡大処理を2回連続して適用することにより4倍拡大が実現される場合に、最初の2倍拡大処理において第1の強調カーネル432a(図12(A))を適用し、2回目の2倍拡大処理において第2の強調カーネル432b(図12(B))を適用する。なお、8倍、16倍以降の倍率においても同様である。
また、参照する画素の位置は、図12に示すような上下左右に限られない。例えば、高画質画像ブロック生成部430は、斜め方向の画素を参照してコントラスト強調を行ったり、さらに離れた画素を参照してコントラスト強調を行ったり、又は、処理対象画像データの種類及びサイズなどにより適用する強調カーネルを切り替えたりしてもよい。
次に、高画質画像ブロック生成部430は、画像ブロック特徴量算出部420により得られた注目領域に対するエッジパターン及び推定エッジ方向と、前述したようにコントラスト強調処理が行われた注目領域とその周辺領域との画素値とを用いて、注目領域に対する拡大画像ブロックを生成する。
図13は、高画質画像ブロック生成部430における拡大画像ブロックの生成処理の具体例の説明図である。まず、高画質画像ブロック生成部430は、注目領域のエッジパターン及び推定エッジ方向基づき、コントラスト強調処理を施された画素値を用いて、3×3画像ブロックに相当する画素値を算出する。図13(A)は、図6(A)に示した注目領域及び周辺領域の一例を示す。前述したように画像ブロック特徴量算出部420において、この注目領域はエッジパターン1010に対応し、推定エッジ方向は「方向1」であると求められている。高画質画像ブロック生成部430は、(エッジパターン1010)−(推定エッジ方向「1」)の組み合わせの場合、図13(B)に示すように、3×3画像ブロックに相当するそれぞれの画素をp0〜p8とすると、図13(A)に示した注目領域の画素値{a,b,c,d}をもとに、p0〜p8の画素値を次の式によって計算する。なお、エッジパターンと推定エッジ方向との組み合わせは、「(エッジパターン)−(推定エッジ方向)」と表される。
p0=a
p1=(a+b)/2
p2=b
p3=(a+c)/2
p4=(b+c)/2
p5=(b+d)/2
p6=c
p7=(b+b)/2
p8=d
これらの計算式は、(エッジパターン)−(推定エッジ方向)の組み合わせにより一意に決定され、3×3画像ブロック相当の画素値が計算される。
図14は、他の(エッジパターン)−(推定エッジ方向)の組み合わせの場合に用いられる計算式の一例の説明図である。
図14(A)は、(エッジパターン1000)−(推定エッジ方向「1」)の場合に用いられる計算式を示し、3×3画像ブロック相当の画素値は、次の式によって計算される。
p0=a
p2=b
p3=a
p4=(b+c)/2
p5=(b+d)/2
p6=c
p7=(c+d)/2
p8=d
p1=(p4+c)/2
図14(B)は(エッジパターン1100)−(推定エッジ方向「5」)の場合に用いられる計算式を示し、3×3画像ブロック相当の画素値は、次の式によって計算される。
p0=a
p1=(a+b)/2
p2=b
p4=(a+d)/2
p6=c
p7=(c+d)/2
p8=d
p3=(p4+c)/2
p5=(p4+b)/2
図14(C)は(エッジパターン1100)−(推定エッジ方向「2」)の場合に用いられる計算式を示し、3×3画像ブロック相当の画素値は、次の式によって計算される。
p0=a
p1=a
p2=b
p3=a
p4=(b+c)/2
p5=(b+d)/2
p6=c
p7=(c+d)/2
p8=d
図14(D)は(エッジパターン0101)−(推定エッジ方向「7」)の場合に用いられる計算式を示し、3×3画像ブロック相当の画素値は、次の式によって計算される。
p0=a
p2=b
p3=(a+c)/2
p4=(a+d)/2
p5=(b+d)/2
p6=c
p8=d
p1=(p4+b)/2
p7=(p4+c)/2
なお、他のエッジパターンの場合にも、同様にそれぞれのエッジパターンに対応した計算式に従って計算を行うことによって、3×3画像ブロック相当の画素値が計算されることができる。
次に、高画質画像ブロック生成部430は、前述のように計算された3×3画像ブロック相当の画素値と、注目領域の推定エッジ方向に基づいて選択された周辺領域内の複数の参照画素とを用いて4×4画像ブロックを生成する。
図15は、注目領域における推定エッジ方向に基づく参照画素r0〜r13の選択方法の説明図である。注目領域の推定エッジ方向が方向1(22.5°)から方向3(67.5°)の場合には、図15(A)に示したように、参照画素r0〜r5は、図15(A)に太線枠で囲んだように左上から下へ3画素と右下から上へ3画素となるように選択される。また、注目領域の推定エッジ方向が方向5(112.5°)から方向7(157.5°)の場合には、図15(B)に示したように、参照画素r0〜r5は、左下から上へ3画素と右上から下へ3画素となるように選択される。参照画素r6〜r13は、推定エッジ方向に拠らず、図15(A)及び図15(B)に示すように、上下それぞれ4画素を選択される。このように、参照画素は、注目領域における推定エッジ方向に基づいて選択される。もちろん、参照画素の選択は、図15に示すように2パターンからの選択に限定されず、推定エッジ方向に従い、より多くの参照画素選択パターンが用意されてもよい。また、選択する参照画素は、推定エッジ方向に応じて変更されてもよい。
図16は、4×4画素の拡大画像ブロックの生成処理の具体例の説明図である。図16に示すように、計算された3×3画像ブロック相当の画素値p0〜p8、及び注目領域における推定エッジ方向に基づいて選択された参照画素r0〜r13を用いて、次の計算式に従って、4×4画素の拡大画像ブロックに相当する画素値(s0〜s15)が計算されて、4x4拡大画像ブロックが生成される。
s0=0.2×r6+0.16×r0+0.64×p0
s1=0.2×r7+0.32×p0+0.48×p1
s2=0.2×r8+0.48×p1+0.32×p2
s3=0.2×r9+0.64×p2+0.16×r1
s4=0.08×r0+0.32×p0+0.12×r2+0.48×p3
s5=0.16×p0+0.24×p1+0.24×p3+0.36×p4
s6=0.24×p1+0.16×p2+0.36×p4+0.24×p5
s7=0.32×p2+0.08×r1+0.48×p5+0.12×r3
s8=0.12×r2+0.48×p3+0.08×r4+0.32×p6
s9=0.24×p3+0.36×p4+0.16×p6+0.24×p7
s10=0.36×p4+0.24×p5+0.24×p7+0.16×p8
s11=0.48×p5+0.12×r3+0.32×p8+0.08×r5
s12=0.16×r4+0.64×p6+0.2×r10
s13=0.32×p6+0.48×p7+0.2×r11
s14=0.48×p7+0.32×p8+0.2×r12
s15=0.64×p8+0.16×r5+0.2×r13
このようにして、4×4画素の拡大画像ブロック(s0〜s15)は、画像ブロック特徴量算出部420により特徴ブロックと判断された注目領域に対して生成される。
以上説明したように、高画質画像ブロック生成部430が、画像ブロック特徴量算出部420により算出された画像特徴量に従って、拡大画像ブロック生成処理を行うことにより、ジャギーを抑えた高画質な拡大処理ができる。
[拡大画像生成部]
次に、拡大画像生成部450をより詳細に説明する。
拡大画像生成部450は、高画質画像ブロック生成部430により生成された注目領域に対する拡大画像ブロックと、高速拡大処理部440から出力された拡大画像とを用いて、原画像に対する拡大画像を生成する。
拡大画像生成部450において、ブロック重畳度判定部452は、エッジ強度算出部422により算出された特徴ブロックのブロック重畳度に基づいて、画像中のエッジ部分と判定された領域に対して、エッジ周辺領域(エッジ周辺部)であるのか最も確からしいエッジ領域(エッジ部)であるのかを判定する。つまり、ブロック重畳度が小さい領域は、画像中において非特徴ブロック(平坦部分)と特徴ブロック(エッジ部分)との境界領域と考えられ、エッジ周辺領域と判定される。また、ブロック重畳度が大きい領域は、周辺領域の多くが特徴ブロックであり、最も確からしいエッジ領域と判定される。例えば、ブロック重畳度1及び2の領域は、エッジ周辺領域と判定され、ブロック重畳度3及び4の領域は、最も確からしいエッジ領域と判定される。判定基準は、先述の例に限られず、例えば、ブロック重畳度の各レベル(レベル1〜レベル4)に対して、より細かくエッジを分類する判定基準が設けられてもよい。
出力値補正部454は、ブロック重畳度判定部452によるブロック重畳度に基づいた判定結果に従って、高画質画像ブロック生成部430により生成された注目領域に対する拡大画像ブロック内の画素値を補正する。
図17は、ブロック重畳度に基づいた判定結果に従う拡大画像ブロック内の画素値補正を具体的に説明する図である。
図17は、ブロック重畳度判定部452における判定の具体例として、ブロック重畳度1及び2の領域はエッジ周辺領域と判定され、ブロック重畳度3及び4の領域は最も確からしいエッジ領域と判定される例を示している。
図17に示すように、ブロック重畳度判定部452において、最も確からしいエッジ領域と判定された領域に相当する、高画質画像ブロック生成部430により生成された拡大画像ブロック内の画素値は、そのまま出力される。一方、ブロック重畳度判定部452において、エッジ周辺領域と判定された領域に相当する拡大画像ブロック内の画素値は、補正されて出力される。
図18は、画素値補正を具体的に説明した図である。
図18に示すように、エッジ周辺領域と判定された領域相当する拡大画像ブロック内の画素値は、高速拡大処理部440から出力された拡大画像上の対応する位置の画素値を用いて以下の式(6)により補正される。
補正値=α×s+(1‐α)×v ・・・(6)
ここで、sは、拡大画像ブロック内の画素値、vは、高速拡大処理部440から出力された拡大画像上の画素値、αは、変数(0≦α≦1)である。ここで、変数αは、ある所定の値であってもよいし、ブロック重畳度に基づいて可変であってもよい。
このように、画像拡大プログラム4は、ブロック重畳度に基づいて、画像中のエッジ部分を、エッジ周辺領域及び最も確からしいエッジ領域のように分類し、最も確からしいエッジ領域では、高画質画像ブロック生成部430により生成された拡大画像ブロックの画素値をそのまま出力し、エッジ周辺領域に対しては、高画質画像ブロック生成部430により生成された拡大画像ブロックの画素値を、高速拡大処理部440による拡大画像の画素値で補正することにより、エッジの鮮鋭度や滑らかさを保持しながら、エッジ周辺領域のアンダーシュートやオーバーシュートといった画質劣化の原因を抑制することが可能となる。
[全体動作]
画像処理装置2の全体動作(画像拡大処理)を説明する。
図19は、画像処理装置2の全体動作を示すフローチャート(S20)である。
図19に示すように、ステップ200(S200)において、画像ブロック設定部410は、画像ブロック特徴量算出部420及び高画質画像ブロック生成部430における処理で必要とされる既定の画像ブロックのサイズをそれぞれ設定し、記憶部400に記憶されている処理対象画像データから、設定されたブロックサイズの画像ブロックを順次(例えばラスタスキャン順に)切り出し、切り出された各画像ブロックを、画像ブロック特徴量算出部420及び高画質画像ブロック生成部430にそれぞれ出力する。
ステップ202(S202)において、画像ブロック特徴量算出部420は、入力された画像ブロック中の注目領域のエッジ強度Gを、式(1)を用いて算出する。
なお、{a,b,c,d}は、図6に例示するように注目領域内の各画素値である。入力画像データがグレースケール画像でなく、例えばRGB色空間のカラー画像である場合には、画像ブロック特徴量算出部420は、注目領域に関してR,G,Bの各色空間の色成分毎の画像ブロックそれぞれについて、式(1)を用いてエッジ強度Gr、Gg、Gbを計算し、Gr、Gg、Gbの中で最大のエッジ強度である色成分の画像ブロックを選択し、そのエッジ強度を注目領域の(すべての色成分に共通の)エッジ強度とする。
ステップ204(S204)において、画像ブロック特徴量算出部420は、注目領域のエッジ強度Gに基づいて、この注目領域が特徴ブロックであるか非特徴ブロックであるかを判定する。具体的には、画像ブロック特徴量算出部420は、算出されたエッジ強度Gが既定の閾値Th以上である場合(すなわち、注目領域内の階調変化量が大きい場合)に、この注目領域を特徴ブロックとし、算出されたエッジ強度Gが既定の閾値Th未満である場合(すなわち、注目領域内の階調変化量が小さい場合)に、この注目領域を非特徴ブロックとする。
画像処理装置2は、注目領域が特徴ブロックであると判定された場合に、この注目領域の画像特徴を保存するような拡大処理行うべくS206の処理に移行し、注目領域が非特徴ブロックであると判定された場合に、S200の処理に戻り、次の画像ブロックを処理する。
ステップ206(S206)において、画像ブロック特徴量算出部420は、特徴ブロックと判定された画像領域に対して、その特徴ブロックが既に特徴ブロックであると判定された画像領域と重なり合う場合は、ブロック重畳度を算出する。
ステップ208(S208)において、画像ブロック特徴量算出部420は、注目領域(特徴ブロック)及びその注目領域を含む1乃至複数の周辺領域中の参照領域のエッジ角度Θを、式(3)を用いて計算する。
なお、gx、gyは、式(1)において各々算出される値である。そして、画像ブロック特徴量算出部420は、算出された複数のエッジ角度Θから注目領域のエッジ方向θを推定する。例えば、画像ブロック特徴量算出部420は、得られた複数のエッジ角度Θの平均値を算出し、算出された平均値を推定エッジ方向θとする。
ステップ210(S210)において、画像ブロック特徴量算出部420は、推定されたエッジ方向θ及び注目領域(特徴ブロック)の画素分布パターンを用いて、エッジパターンを選択する。エッジパターンは、エッジ方向及び画素分布パターンに応じて、予め用意されたパターンテーブル(図10)の中から選択される。
ステップ212(S212)において、高画質画像ブロック生成部430は、拡大率に応じてサイズ及び重み付け係数が設定された強調カーネル432(図12)を用いて、画像ブロック設定部410により切り出された画像ブロック中の注目領域(特徴ブロック)及びその周辺領域の画像データに対して、コントラスト強調処理を施す。
ステップ214(S214)において、高画質画像ブロック生成部430は、画像ブロック特徴量算出部420により得られた注目領域のエッジパターン及び前記注目領域の推定エッジ方向θと、S212においてコントラスト強調が施された注目領域及び周辺領域内の画素値を用いて、エッジパターンと推定エッジ方向θとに対応する算出式(図13、図14)により、3×3画像ブロックに相当する画素値(p0〜p8)を生成する。さらに、高画質画像ブロック生成部430は、生成された画素値(p0〜p8)と推定エッジ方向θとに基づいて選択された参照画素(r0〜r13)を用いて、図16に示す算出式により、注目領域に対する拡大画像ブロックを生成する。
ステップ216(S216)において、画像処理装置2は、全ての入力画像データについて、S200からS214までの処理が完了したか否かを判定し、処理が完了していないと判定された場合に、S200の処理に戻り、次の画素ブロックに対する処理を行い、全ての入力画像データについて処理が完了していると判定された場合に、S218の処理に移行する。
ステップ218(S218)において、高速拡大処理部440は、記憶部400に記憶されている入力画像データを、画像毎または数ライン毎に入力し、最近傍補間法により拡大処理を行う。
なお、本例においては、高速拡大処理部440は、高画質画像ブロック生成部430による拡大処理(S200〜S214までの処理)が終了した後に拡大処理を行っているが、高画質画像ブロック生成部430による拡大処理と並行して拡大処理を行ってもよいし、高画質画像ブロック生成部430による拡大処理の前に拡大処理を行ってもよい。
ステップ220(S220)において、拡大画像生成部450は、S206において算出されたブロック重畳度に基づき、高画質画像ブロック生成部430により生成された拡大画像と、高速拡大処理部440により拡大された拡大画像を統合する。
拡大画像生成部450は、拡大画素の統合方法として、ブロック重畳度が大きい領域に対しては、高画質画像ブロック生成部430により生成された拡大画像の値を出力値とし、ブロック重畳度が小さい領域に対しては、高画質画像ブロック生成部430により生成された拡大画像と高速拡大処理部440により生成された拡大画像とを用いて、式(6)により算出される値を出力値とする。
以上説明したように、本実施形態における画像処理装置2は、注目画素を含む所定の大きさの画像領域の特徴量を算出し、算出された特徴量が所定の条件を満たす画像領域に対して拡大画像領域を生成し、さらに、異なる拡大手法により入力画像の拡大画像を生成し、算出された特徴量が所定の条件を満たす画像領域が重なり合う場合、その重なり数に基づいて、生成された拡大画像領域の生成値を補正し、補正された拡大画像領域と、異なる拡大手法により生成された拡大画像とに基づいて、入力画像に対する拡大画像を生成する。
このような拡大処理によって、本実施形態における画像処理装置2は、高画質拡大処理を特徴的な部分(特徴ブロック)にのみ適用することにより、全体としての処理負荷を抑え、かつ、特徴的な部分について高画質拡大処理を適用するので、特徴的な部分の特徴量を保存して、ボケやジャギーなどの画質欠陥を抑制した高画質な拡大画像を高速に得ることができる。また、本実施形態における画像処理装置2は、エッジの鮮鋭度や滑らかさを保持しながら、エッジ部周辺のアンダーシュートやオーバーシュートといった画質欠陥を抑制することができる。
本発明にかかる画像処理方法が適応される画像処理装置2のハードウェア構成を、制御装置20を中心に例示する図である。 制御装置20により実行され、本発明にかかる画像処理方法を実現する画像拡大プログラム4の機能構成を示す図である。 エッジ強度算出部422において、エッジ強度Gと基準値との比較により特徴ブロックと判定された領域を例示する。 互いに重なり合う特徴ブロックを例示する。 特徴ブロック同士の重なり数(ブロック重畳度)を例示する。 注目領域及び周辺領域の具体例と、注目領域のエッジ方向の一例の説明図であり、図6(A)は、注目領域及び周辺領域を例示し、図6(B)は、この注目領域について推定されるエッジ方向を例示する図である。 エッジ方向推定部424によるエッジ方向推定処理(S10)を示すフローチャートである。 エッジ方向推定処理(S10)のS102の処理において選択される参照領域を例示する図である。 図6に示した注目領域における推定エッジ方向の一例の説明図である。 エッジパターン選択部426により用いられるエッジパターンテーブルを例示する図である。 図10に示した注目領域に対応するエッジパターンの選択方法を説明する図である。 高画質画像ブロック生成部430により用いられる強調カーネル432(エッジ強調カーネル)を例示する図である。 高画質画像ブロック生成部430における拡大画像ブロックの生成処理の具体例の説明図である。 エッジパターンと推定エッジ方向との組み合わせの場合に用いられる計算式の一例の説明図である。 注目領域における推定エッジ方向に基づく参照画素r0〜r13の選択方法の説明図である。 4×4画素の拡大画像ブロックの生成処理の具体例の説明図である。 ブロック重畳度に基づいた判定結果に従う拡大画像ブロック内の画素値補正を具体的に説明する図である。 画素値補正を具体的に説明した図である。 画像処理装置2の全体動作を示すフローチャート(S20)である。
符号の説明
2・・・画像処理装置
10・・・プリンタ装置
20・・・制御装置
202・・・CPU
204・・・メモリ
22・・・通信装置
24・・・記憶装置
26・・・UI装置
4・・・画像拡大プログラム
400・・・記憶部
410・・・画像ブロック設定部
420・・・画像ブロック特徴量算出部
422・・・エッジ強度算出部
424・・・エッジ方向推定部
426・・・エッジパターン選択部
430・・・高画質画像ブロック生成部
432a、432b・・・強調カーネル
434・・・エッジ方向推定部
440・・・高速拡大処理部
450・・・拡大画像生成部
452・・・ブロック重畳度判定部
454・・・出力値補正部

Claims (29)

  1. 入力画像の拡大処理を行う画像処理装置であって、
    注目画素を含む所定の大きさの画像領域の特徴量を算出する領域特徴量算出手段と、
    前記領域特徴量算出手段により算出された画像領域の特徴量に基づいて、この画像領域に対応する拡大画像領域を生成する拡大画像領域生成手段と、
    前記拡大画像領域生成手段とは異なる拡大手法を適用して、入力画像の拡大画像を生成する拡大処理手段と、
    前記領域特徴量算出手段により算出された特徴量が所定の条件を満たす画像領域が重畳する場合、その重畳する度合いに基づいて、前記拡大画像領域生成手段により生成された拡大画像領域の生成値を補正し、この拡大画像領域と、前記拡大処理手段により生成された拡大画像とに基づいて、入力画像に対する拡大画像を生成する拡大画像生成手段と
    を有する画像処理装置。
  2. 前記領域特徴量算出手段は、画像領域内の各画素値に基づいて、画像領域における階調変化量を算出し、この階調変化量を画像領域の特徴量とする
    請求項1に記載の画像処理装置。
  3. 前記領域特徴量算出手段は、画像領域内の各画素値に基づいて、画像領域における階調変化方向を算出し、この階調変化方向を画像領域の特徴量とする
    請求項1に記載の画像処理装置。
  4. 前記領域特徴量算出手段は、画像領域内の各画素値に基づいて、画像領域における階調変化を算出し、この階調変化が既定の画素値パターンに対応する場合に、この画素値パターンを画像領域の特徴量として選択する
    請求項1に記載の画像処理装置。
  5. 前記領域特徴量算出手段は、色空間における色成分それぞれについて特徴量を算出し、算出された特徴量に基づいて1つの色成分を選択し、この色成分について算出された特徴量を画像領域の特徴量とする
    請求項1に記載の画像処理装置。
  6. 前記拡大画像領域生成手段は、前記領域特徴量算出手段により算出された画像領域の特徴量が所定の条件を満たす場合に、画像領域の特徴量と、この画像領域の近傍領域の特徴量と、この近傍領域内の画素値とを用いて、拡大画像領域を生成する
    請求項1に記載の画像処理装置。
  7. 前記拡大画像領域生成手段は、階調変化を強調するエッジ強調処理を拡大率に応じて行って、画像領域内の画素値を補正し、補正された画素値を用いて拡大画像領域を生成する
    請求項1に記載の画像処理装置。
  8. 前記拡大画像領域生成手段は、前記領域特徴量算出手段により算出された階調変化方向に応じて、画像領域の近傍領域内の画素値を選択し、選択された画素値を用いて拡大画像領域を生成する
    請求項3に記載の画像処理装置。
  9. 前記拡大画像領域生成手段は、前記領域特徴量算出手段により算出された特徴量に対応する既定の演算式を用いて、拡大画像領域を生成する
    請求項1に記載の画像処理装置。
  10. 前記拡大処理手段は、少なくとも最近傍補間法に基づいた拡大手法を適用する
    請求項1に記載の画像処理装置。
  11. 前記拡大処理手段は、少なくとも線形補間法に基づいた拡大手法を適用する
    請求項1に記載の画像処理装置。
  12. 前記拡大処理手段は、入力画像毎あるいは入力画像中の数ライン毎に、拡大画像を生成する
    請求項1に記載の画像処理装置。
  13. 前記拡大画像生成手段は、前記拡大画像領域生成手段により生成された拡大画像領域を、前記拡大処理手段により生成された拡大画像上の対応する位置に配置する
    請求項1に記載の画像処理装置。
  14. 前記拡大画像生成手段は、前記拡大画像領域生成手段により生成された前記拡大画像領域の生成値を、前記拡大処理手段により生成された拡大画像上の対応する位置の画素値を用いて補正する
    請求項1に記載の画像処理装置。
  15. 前記拡大画像生成手段は、前記拡大画像領域生成手段により生成された拡大画像領域の生成値と、前記拡大処理手段により生成された拡大画像上の対応する位置の画素値とを、所定の割合で線形和した値を用いて、この画素値を補正する
    請求項14に記載の画像処理装置。
  16. 前記拡大画像生成手段は、複数の画像領域において前記領域特徴量算出手段により算出された画像領域の特徴量が所定の条件を満たし、これらの画像領域が互いに重畳する場合、これらの画像領域が重畳する度合いに基づいて、線形和の割合を決定する
    請求項15に記載の画像処理装置。
  17. 前記拡大画像生成手段は、前記領域特徴量算出手段により算出された画像領域の特徴量と、この画像領域と重畳する画像領域の特徴量とに基づいて、入力画像のエッジ部及びエッジ周辺部を判定する
    請求項1に記載の画像処理装置。
  18. 前記拡大画像生成手段は、複数の画像領域において前記領域特徴量算出手段により算出された画像領域の特徴量が所定の条件を満たし、これらの画像領域が互いに重畳する場合、これらの画像領域が重畳する度合いが所定値以上である画像部分を、エッジ部と判定する
    請求項17に記載の画像処理装置。
  19. 前記拡大画像生成手段は、複数の画像領域において前記領域特徴量算出手段により算出された画像領域の特徴量が所定の条件を満たし、これらの画像領域が互いに重畳する場合、これらの画像領域が重畳する度合いが所定値以下である画像部分を、エッジ周辺部と判定する
    請求項17に記載の画像処理装置。
  20. 前記拡大画像生成手段は、前記拡大画像領域生成手段により生成された拡大画像領域で、エッジ部に対応する画像部分の生成値を、前記拡大処理手段により生成された拡大画像の対応する位置に配置する
    請求項17又は18に記載の画像処理装置。
  21. 前記拡大画像生成手段は、前記拡大画像領域生成手段により生成された拡大画像領域で、エッジ周辺部に対応する画像部分の生成値と、前記拡大処理手段により生成された拡大画像の対応する位置の画素値とを用いて、この画素値を補正する
    請求項17又は19に記載の画像処理装置。
  22. 前記拡大画像生成手段は、前記拡大画像領域生成手段により生成された拡大画像領域で、エッジ部及びエッジ周辺部のいずれにも対応しない画像部分には、前記拡大処理手段により生成された拡大画像上の対応する位置の画素値を用いる
    請求項17乃至19のいずれかに記載の画像処理装置。
  23. 入力画像のエッジを検出するエッジ検出装置であって、
    注目画素を含む所定の大きさの画像領域の特徴量を算出する領域特徴量算出手段と、
    前記領域特徴量算出手段により算出された画像領域の特徴量と、この画像領域と重畳する画像領域の特徴量とに基づいて、入力画像のエッジ部及びエッジ周辺部を判定する判定手段と
    を有するエッジ検出装置。
  24. 前記判定手段は、複数の画像領域において前記領域特徴量算出手段により算出された画像領域の特徴量が所定の条件を満たし、これらの画像領域が互いに重畳する場合、これらの画像領域が重畳する度合いが所定値以上である画像部分を、エッジ部と判定する
    請求項23に記載のエッジ検出装置。
  25. 前記判定手段は、複数の画像領域において前記領域特徴量算出手段により算出された画像領域の特徴量が所定の条件を満たし、これらの画像領域が互いに重畳する場合、これらの画像領域が重畳する度合いが所定値以下である画像部分を、エッジ周辺部と判定する
    請求項23に記載のエッジ検出装置。
  26. 入力画像の拡大処理を行う画像処理方法であって、
    注目画素を含む所定の大きさの画像領域の特徴量を算出し、
    前記算出された画像領域の特徴量に基づいて、この画像領域に対応する拡大画像領域を生成し、
    前記生成とは異なる拡大手法を適用して、入力画像の拡大画像を生成し、
    前記算出された特徴量が所定の条件を満たす画像領域が重畳する場合、その重畳する度合いに基づいて、前記生成された拡大画像領域の生成値を補正し、この拡大画像領域と、前記生成された拡大画像とに基づいて、入力画像に対する拡大画像を生成する
    画像処理方法。
  27. 入力画像のエッジを検出するエッジ検出方法であって、
    注目画素を含む所定の大きさの画像領域の特徴量を算出し、
    前記算出された画像領域の特徴量と、この画像領域と重畳する画像領域の特徴量とに基づいて、入力画像のエッジ部及びエッジ周辺部を判定する
    エッジ検出方法。
  28. コンピュータを含み、入力画像の拡大処理を行う画像処理装置において、
    注目画素を含む所定の大きさの画像領域の特徴量を算出する領域特徴量算出ステップと、
    前記算出された画像領域の特徴量に基づいて、この画像領域に対応する拡大画像領域を生成する拡大画像領域生成ステップと、
    前記生成とは異なる拡大手法を適用して、入力画像の拡大画像を生成する拡大処理ステップと、
    前記算出された特徴量が所定の条件を満たす画像領域が重畳する場合、その重畳する度合いに基づいて、前記生成された拡大画像領域の生成値を補正し、この拡大画像領域と、前記生成された拡大画像とに基づいて、入力画像に対する拡大画像を生成する拡大画像生成ステップと
    を前記画像処理装置のコンピュータに実行させるプログラム。
  29. コンピュータを含む入力画像のエッジを検出するエッジ検出装置でおいて、
    注目画素を含む所定の大きさの画像領域の特徴量を算出する領域特徴量算出ステップと、
    前記算出された画像領域の特徴量と、この画像領域と重畳する画像領域の特徴量とに基づいて、入力画像のエッジ部及びエッジ周辺部を判定する判定ステップと
    を前記エッジ検出装置のコンピュータに実行させるプログラム。
JP2005169011A 2005-06-09 2005-06-09 画像処理装置、画像処理方法、及びそのプログラム Expired - Fee Related JP4716255B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2005169011A JP4716255B2 (ja) 2005-06-09 2005-06-09 画像処理装置、画像処理方法、及びそのプログラム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2005169011A JP4716255B2 (ja) 2005-06-09 2005-06-09 画像処理装置、画像処理方法、及びそのプログラム

Publications (2)

Publication Number Publication Date
JP2006344002A true JP2006344002A (ja) 2006-12-21
JP4716255B2 JP4716255B2 (ja) 2011-07-06

Family

ID=37640929

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2005169011A Expired - Fee Related JP4716255B2 (ja) 2005-06-09 2005-06-09 画像処理装置、画像処理方法、及びそのプログラム

Country Status (1)

Country Link
JP (1) JP4716255B2 (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009140435A (ja) * 2007-12-10 2009-06-25 Sharp Corp 画像処理装置、画像表示装置、画像形成装置、画像処理方法、コンピュータプログラム及び記憶媒体
JP2010055374A (ja) * 2008-08-28 2010-03-11 Sharp Corp データ作成装置、データ作成方法、データ作成用プログラム、描画装置、描画方法、描画用プログラム、および、コンピュータ読取可能な記録媒体
JP2010055380A (ja) * 2008-08-28 2010-03-11 Sharp Corp データ作成装置、データ作成方法、データ作成用プログラム、描画装置、描画方法、描画用プログラム、および、コンピュータ読取可能な記録媒体
CN102564345A (zh) * 2010-11-20 2012-07-11 株式会社高永科技 检查方法
US8493416B2 (en) 2009-04-24 2013-07-23 Fuji Xerox Co., Ltd. Image processing apparatus, image processing method, and computer readable medium

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005123813A (ja) * 2003-10-15 2005-05-12 Fuji Xerox Co Ltd 画像処理装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005123813A (ja) * 2003-10-15 2005-05-12 Fuji Xerox Co Ltd 画像処理装置

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009140435A (ja) * 2007-12-10 2009-06-25 Sharp Corp 画像処理装置、画像表示装置、画像形成装置、画像処理方法、コンピュータプログラム及び記憶媒体
JP4510069B2 (ja) * 2007-12-10 2010-07-21 シャープ株式会社 画像処理装置、画像表示装置、画像形成装置、画像処理方法、コンピュータプログラム及び記憶媒体
JP2010055374A (ja) * 2008-08-28 2010-03-11 Sharp Corp データ作成装置、データ作成方法、データ作成用プログラム、描画装置、描画方法、描画用プログラム、および、コンピュータ読取可能な記録媒体
JP2010055380A (ja) * 2008-08-28 2010-03-11 Sharp Corp データ作成装置、データ作成方法、データ作成用プログラム、描画装置、描画方法、描画用プログラム、および、コンピュータ読取可能な記録媒体
JP4698709B2 (ja) * 2008-08-28 2011-06-08 シャープ株式会社 データ作成装置、データ作成方法、データ作成用プログラム、描画装置、描画方法、描画用プログラム、および、コンピュータ読取可能な記録媒体
JP4703695B2 (ja) * 2008-08-28 2011-06-15 シャープ株式会社 データ作成装置、データ作成方法、データ作成用プログラム、描画装置、描画方法、および、描画用プログラム
US8493416B2 (en) 2009-04-24 2013-07-23 Fuji Xerox Co., Ltd. Image processing apparatus, image processing method, and computer readable medium
CN102564345A (zh) * 2010-11-20 2012-07-11 株式会社高永科技 检查方法

Also Published As

Publication number Publication date
JP4716255B2 (ja) 2011-07-06

Similar Documents

Publication Publication Date Title
JP3915563B2 (ja) 画像処理装置および画像処理プログラム
EP2138975B1 (en) Apparatus, method and computer-readable recording medium for pixel interpolation
JP4716255B2 (ja) 画像処理装置、画像処理方法、及びそのプログラム
JP4507279B2 (ja) 画像処理装置、画像処理方法及びそのプログラム
CN109993693B (zh) 用于对图像进行插值的方法和装置
JP4126541B2 (ja) 画像処理装置及び画像処理方法、画像処理プログラム、記憶媒体
JP4868249B2 (ja) 映像信号処理装置
US7492969B2 (en) Method for magnifying image by interpolation based on edge and corner
JP4730525B2 (ja) 画像処理装置及びそのプログラム
JP4857975B2 (ja) 画像処理システムおよび画像処理プログラム
JP2006050481A (ja) 画像処理装置、画像処理方法及びそのプログラム
JP4265363B2 (ja) 画像処理装置
JP2019121061A (ja) 画像拡大装置
JP2006252108A (ja) 画像処理装置、画像処理方法及びそのプログラム
JP6444049B2 (ja) 画像処理装置、方法およびプログラム
JP5013243B2 (ja) 画像処理装置、画像処理方法、画像処理プログラム及びこれを記録した記録媒体
US9584699B2 (en) Image processing apparatus, image processing method and computer-readable storage medium with improved character display
JP4879084B2 (ja) 画像処理装置
JP4736121B2 (ja) 画像解析装置、画像処理装置、画像解析方法、画像解析プログラム及びこれを記録した記録媒体
JP2006050480A (ja) 画像処理装置、画像処理方法及びそのプログラム
JP2006252107A (ja) 画像処理装置、画像処理方法及びそのプログラム
JP4517288B2 (ja) 画像処理装置、画像処理方法及びそのプログラム
JP2015106318A (ja) 画像処理装置および画像処理方法
JP4265362B2 (ja) 画像処理装置
JP4517287B2 (ja) 画像処理装置、画像処理方法及びそのプログラム

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20080317

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20101026

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20101220

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20110215

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

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20110307

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20110320

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

Free format text: PAYMENT UNTIL: 20140408

Year of fee payment: 3

LAPS Cancellation because of no payment of annual fees