(57)【要約】【課題】入力された画像に対してボケやジャギーを生じさせること無く、鮮鋭度が高い滑らかなエッジの再現と、さらにはテクスチャなどの細部の構造を再現することができる拡大処理を、処理負荷が小さく高速で行うことが可能な画像処理システムを提供する。【解決手段】画像処理システムの周波数変換手段は、所定の大きさの画像領域を複数の周波数成分に変換し、周波数成分生成手段は、前記周波数変換手段によって変換された前記画像領域の周波数成分から所定の周波数成分を生成し、拡大画像領域生成手段は、注目領域を含む所定の大きさの画像領域と前記周波数成分生成手段によって生成された所定の周波数成分を用いて逆変換することにより前記注目領域に対する拡大画像領域を生成し、拡大画像生成手段は、前記拡大画像領域生成手段により生成された拡大画像領域を所定の方法で配置して入力画像に対する拡大画像を生成する。【選択図】図１
2008059551.tif 000002
【特許請求の範囲】
【請求項１】
所定の大きさの画像領域を複数の周波数成分に変換する周波数変換手段と、
前記周波数変換手段によって変換された前記画像領域の周波数成分から所定の周波数成
分を生成する周波数成分生成手段と、
注目領域を含む所定の大きさの画像領域と前記周波数成分生成手段によって生成された
所定の周波数成分を用いて逆変換することにより前記注目領域に対する拡大画像領域を生
成する拡大画像領域生成手段と、
前記拡大画像領域生成手段により生成された拡大画像領域を所定の方法で配置して入力
画像に対する拡大画像を生成する拡大画像生成手段
を有することを特徴とする画像処理システム。
【請求項２】
所定の大きさの画像領域を複数の周波数成分に変換する周波数変換手段と、
画像領域を分類する分類手段と、
画像領域を強調する強調手段と、
前記周波数変換手段によって変換された前記画像領域の周波数成分から所定の周波数成
分を生成する周波数成分生成手段と、
注目領域を含む所定の大きさの画像領域と前記周波数成分生成手段によって生成された
所定の周波数成分を用いて逆変換することにより前記注目領域に対する拡大画像領域を生
成する拡大画像領域生成手段と、
前記拡大画像領域生成手段とは異なる手法で前記注目領域に対する拡大画像領域を生成
する第２の拡大画像領域生成手段と、
前記拡大画像領域生成手段および前記第２の拡大画像領域生成手段により生成された拡
大画像領域を所定の方法で配置して入力画像に対する拡大画像を生成する拡大画像生成手
段
を有することを特徴とする画像処理システム。
【請求項３】
前記周波数変換手段は、ウェーブレット変換を行う
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項４】
前記周波数変換手段は、Ｈａａｒウェーブレットを用いたウェーブレット変換を行う
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項５】
前記周波数成分生成手段は、画像領域の高周波数成分を生成する
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項６】
前記周波数成分生成手段は、複数の画像領域から高周波数成分を生成する
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項７】
前記周波数成分生成手段は、注目領域の近傍の既定位置の画像領域から高周波数成分を
生成する
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項８】
前記周波数成分生成手段は、注目領域の特徴に基づいて選択された画像領域から周波数
成分を生成する
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項９】
前記周波数成分生成手段は、注目領域の特徴として注目領域のエッジ方向を算出する
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項１０】
前記周波数成分生成手段は、算出した注目領域のエッジ方向に基づいて選択された画像
領域から周波数成分を生成する
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項１１】
前記周波数成分生成手段は、算出した注目領域のエッジ方向に基づいた既定位置の画像
領域から周波数成分を生成する
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項１２】
前記周波数成分生成手段は、高周波数成分生成に用いる注目領域のエッジ方向に基づい
た複数の画像領域の各々の位置を予め位置情報として備えている
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項１３】
前記周波数成分生成手段は、高周波数成分生成に用いる複数の画像領域の既定位置を注
目領域の特徴と相関の強い複数の画像領域の位置とする
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項１４】
前記周波数成分生成手段は、画像領域の統計量に基づいて選択された画像領域から周波数成分を生成する
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項１５】
前記周波数成分生成手段は、画像領域の統計量として画像領域内の各画素と画像領域の平均値との差分である偏差を用いて、前記偏差に基づいて選択された画像領域から周波数成分を生成する
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項１６】
前記周波数成分生成手段は、画像領域の統計量として前記偏差から算出される値を用いて、前記偏差に基づいて選択された画像領域から周波数成分を生成する
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項１７】
前記周波数成分生成手段は、画像領域の統計量として画像領域の分散を用いて、前記分散に基づいて選択された画像領域から周波数成分を生成する
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項１８】
前記周波数成分生成手段は、複数の画像領域の統計量を算出し、既定の基準に基づいて、複数の画像領域の中から所定の統計量を持つ画像領域を選択し、前記選択された画像領域から周波数成分を生成する
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項１９】
前記周波数成分生成手段は、複数の画像領域の統計量を算出し、注目領域の特徴として注目領域のエッジ強度を算出し、前記エッジ強度に基づいて、複数の画像領域の中から所定の統計量を持つ画像領域を選択し、前記選択された画像領域から周波数成分を生成する
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項２０】
前記周波数成分生成手段は、複数の画像領域の統計量を算出し、注目領域中の所定画素の画素値に基づいて、複数の画像領域の中から所定の統計量を持つ画像領域を選択し、前記選択された画像領域から周波数成分を生成する
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項２１】
前記周波数成分生成手段は、注目領域の高周波数成分に近い周波数成分を持つ画像領域
から高周波数成分を生成する
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項２２】
前記周波数成分生成手段は、注目領域の高周波数成分が所定の基準値よりも小さいとき
は生成する高周波数成分値を０とする
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項２３】
前記拡大画像領域生成手段は、逆ウェーブレット変換を行う
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項２４】
前記拡大画像生成手段は、順次生成される前記拡大画像領域を重なり合うように順次配
置する
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項２５】
前記拡大画像生成手段は、順次生成される前記拡大画像領域の重なり合う画素に対して
、重なる画素値の総和を算出し、前記画素値の総和を重なった数で割ることにより前記拡
大画像領域の画素値を算出する
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項２６】
前記拡大画像生成手段は、順次生成される前記拡大画像領域を重なり合わないように順
次配置する
ことを特徴とする請求項１または２に記載の画像処理システム。
【請求項２７】
前記分類手段は、注目領域の周波数成分に基づいて算出される値により注目領域を分類
する
ことを特徴とする請求項２に記載の画像処理システム。
【請求項２８】
前記強調手段は、所定のカーネル要素およびカーネル要素間距離をもつエッジ強調カー
ネルを用いて強調処理を行う
ことを特徴とする請求項２に記載の画像処理システム。
【請求項２９】
前記強調手段は、前記分類手段による注目領域の分類結果に基づいて異なるエッジ強調
カーネルを用いて強調処理を行う
ことを特徴とする請求項２に記載の画像処理システム。
【請求項３０】
前記拡大画像領域生成手段は、前記分類手段によってテクスチャ領域と分類された画像
領域を前記注目領域とし、
前記第２の拡大画像領域生成手段は、前記分類手段によってエッジ領域と分類された画
像領域を前記注目領域とした
ことを特徴とする請求項２に記載の画像処理システム。
【請求項３１】
前記強調手段は、前記分類手段によってテクスチャ領域と分類された画像領域に対して
コントラストの強調を行い、
前記拡大画像領域生成手段は、前記強調手段によってコントラストが強調された画像領
域を前記注目領域とし、
前記第２の拡大画像領域生成手段は、前記分類手段によってエッジ領域と分類された画
像領域を前記注目領域とした
ことを特徴とする請求項２に記載の画像処理システム。
【請求項３２】
コンピュータに、請求項１ないし３１のいずれか１項に記載の画像処理システムの機能
を実現させることを特徴とする画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像処理システムおよび画像処理プログラムに関する。
【背景技術】
【０００２】
画像の拡大処理は、画像の編集やファイリング、表示、印刷などを行うシステムにとっ
て基本的な処理の一つである。また近年、インターネットのホームページ上の画像やデジ
タルビデオなどのディスプレイ解像度での表示を主目的とした画像データなどの普及によ
り、これらの低解像度画像を高解像度のプリンタなどで印刷することも頻繁に行われてい
る。このプリンタによる印刷の際に、高画質の出力結果を得ることが望まれており、高画
質の拡大処理に対する重要度が高まっている。
【０００３】
カラーを含む多階調で表現された画像（以下では、これを多値画像と称する場合もある
）を拡大処理する既存の代表的な手法としては、最近傍法や線形補間法、キュービック・
コンボリューション法などがある。最近傍法は、拡大後の各画素値として、その画素を原
画像上に逆写像した際に最も距離が近い画素の画素値を使うという方法である。この方法
は、演算量が少ないため高速に処理することができる。しかし、原画像の１画素がそのま
ま矩形形状に拡大されるため、隣り合う画素の画素値の差が小さい場合は画質劣化の程度
は小さくほとんど影響はないが、逆に大きい場合などは、斜線部やエッジ部のジャギーが
目立ったり、倍率が大きい場合には画像がモザイク状になるなど、画質劣化の程度は大き
い。
【０００４】
線形補間法は、画素間の画素値が直線的に変化していると仮定し、拡大後の画素を逆写
像した点の近傍４画素の画素値を線形に補間して画素値を求めるという方法である。この
方法では、最近傍法よりも処理は重いものの演算量は比較的少なく、ジャギーなども発生
しにくい。その一方で、直線的に変化しているという仮定に当てはまらないエッジ部分を
中心に、画像全体がボケ気味になるという欠点がある。
【０００５】
キュービック・コンボリューション法は、標本化定理に基づいてｓｉｎｃ関数（ｓｉｎ
（ｘ）／ｘ）を近似した補間関数を定義し、拡大後の画素を逆写像した点の近傍１６画素
（Ｘ、Ｙ方向それぞれ４画素）と前記の近似補間関数との畳み込み演算により、拡大後の
画素値を求める方法である。この方法は、前記２つの手法に比べて画質は比較的良いが、
拡大倍率が大きくなるとエッジ部分でジャギーが発生したり、高域が強調気味となる特性
を持つためノイズ成分が強調されてしまうなどの欠点もある。
【０００６】
近年ではこれら拡大画像の画質問題を解決する試みとして、ＤＣＴ（ｄｉｓｃｒｅｔｅ
ｃｏｓｉｎｅ  ｔｒａｎｓｆｏｒｍ、離散コサイン変換）などの直交変換を用いて、実
空間の画像信号を周波数空間の信号に変換し、何らかの高周波数成分を追加したのち逆変
換を行う事により、サンプリング時に失われた高周波数成分を復元する拡大手法が提案さ
れている。例えば特許文献１、特許文献２などの新規方式が提案されている。特許文献１
に記載されている技術では、実空間の画像信号に対してＤＣＴ、ＩＤＣＴ（ｉｎｖｅｒｓ
ｅ  ｄｉｓｃｒｅｔｅ  ｃｏｓｉｎｅ  ｔｒａｎｓｆｏｒｍ、逆離散コサイン変換）を繰
り返して行うゲルヒベルグ‐パポリスの反復法を利用して、拡大画像における高周波数成
分の復元を実現する技術が記載されている。しかし、例えば特許文献１による技術では、
ＤＣＴ、ＩＤＣＴの反復演算を行うため演算量が非常に多く、対象画像がカラー画像の場
合などではさらに処理負荷が大きくなる。
【０００７】
例えば特許文献２に記載されている技術では、画像をブロック毎にＤＣＴして周波数信
号に変換し、前記周波数信号に変換された画像ブロックの高周波部分に０を補間して、通
常ＩＤＣＴで用いる基底とは異なる基底を複数用意し、画像ブロックの特徴に基づき、用
意した基底を切替えてＩＤＣＴする手法が提案されている。しかし特許文献２による手法
では、ＩＤＣＴで用いる複数の基底を実験的に手動で作成しており、あらかじめ作成され
た基底だけでは多種多様な画像ブロックの特徴に基づくＩＤＣＴが行えず、エッジ部のジ
ャギーの発生や、リンギングなどの画質劣化が起こることが考えられる。
【０００８】
また、これらに関する技術として、例えば特許文献３には、画像平坦部におけるリンギ
ングを抑制することにより、拡大画像の高画質化を図ることを目的とし、リンギングの発
生のない、平坦部が滑らかな参照画像を用意し、処理対象画像を複数領域に分け、平坦部
領域を求め、注目画素を中心とする矩形領域に対しリンギング修正処理を行い、処理対象
画像と参照画像との誤差を画像の平坦部において少なくすることにより、リンギングの抑
制を図る技術が記載されている。
【０００９】
また、例えば特許文献４には、画像内容に逆らうことなく可視的な劣化を感じることが
ないようにしつつ、補間処理の負担を最小限に抑えることができる画像処理方法を提供す
ることを目的とし、空間周波数領域へ直交変換を用いて複数画素単位で周波数（交流）成
分を含んで圧縮された画像データを、圧縮データから非圧縮データへ復元する時に、圧縮
の単位である複数の部分画素領域についての周波数（交流）成分を含むデータを抽出する
ことを特徴とする画像処理方法において、抽出した複数の部分画素領域の周波数（交流）
成分の特徴量と、予め設定してある空間周波数特性閾値とを比較し、周波数特性として高
域周波数データブロックと低域周波数データブロックに判定し、その結果を複数の部分画
素領域の位置情報に関連付け、画像データ変倍用の補間精度の違うフィルタを少なくとも
２種類備え、画像データを変倍（拡大若しくは縮小）する時に、先に求めた結果を元に複
数の部分画素領域単位で変倍用フィルタを切り替えて使用する技術が記載されている。
【００１０】
また、例えば特許文献５には、画像の特徴抽出が高速にでき、更に、レイアウトや形状
やテクスチャといった特徴を個別的に抽出可能な柔軟性の高い特徴抽出を行うことを課題
とし、グレースケール画像又はカラー画像を入力してウエーブレット変換を施すウエーブ
レット変換部と、ウエーブレット変換後の各周波数帯域の係数に基づいて、画像のレイア
ウト、形状およびテクスチャの各々に分けて特徴量を抽出するレイアウト特徴量抽出部、
形状特徴量抽出部及びテクスチャ特徴量抽出部を備える技術が記載されている。
【特許文献１】特許第３３０２７３１号公報
【特許文献２】特許第３１９５１４２号公報
【特許文献３】特許第３４０４１３８号公報
【特許文献４】特開２００３−２７４１５４号公報
【特許文献５】特開２００１−０９２９７３号公報
【発明の開示】
【発明が解決しようとする課題】
【００１１】
本発明は、このような背景技術の状況の中でなされたもので、入力された画像に対して
ボケやジャギーを生じさせること無く、鮮鋭度が高い滑らかなエッジの再現と、さらには
テクスチャなどの細部の構造を再現することができる拡大処理を、処理負荷が小さく高速
で行うことが可能な画像処理システムおよび画像処理プログラムを提供することを目的と
するものである。
【課題を解決するための手段】
【００１２】
かかる目的を達成するための本発明の要旨とするところは、次の各項の発明に存する。
［１］  所定の大きさの画像領域を複数の周波数成分に変換する周波数変換手段と、
前記周波数変換手段によって変換された前記画像領域の周波数成分から所定の周波数成
分を生成する周波数成分生成手段と、
注目領域を含む所定の大きさの画像領域と前記周波数成分生成手段によって生成された
所定の周波数成分を用いて逆変換することにより前記注目領域に対する拡大画像領域を生
成する拡大画像領域生成手段と、
前記拡大画像領域生成手段により生成された拡大画像領域を所定の方法で配置して入力
画像に対する拡大画像を生成する拡大画像生成手段
を有することを特徴とする画像処理システム。
【００１３】
［２］  所定の大きさの画像領域を複数の周波数成分に変換する周波数変換手段と、
画像領域を分類する分類手段と、
画像領域を強調する強調手段と、
前記周波数変換手段によって変換された前記画像領域の周波数成分から所定の周波数成
分を生成する周波数成分生成手段と、
注目領域を含む所定の大きさの画像領域と前記周波数成分生成手段によって生成された
所定の周波数成分を用いて逆変換することにより前記注目領域に対する拡大画像領域を生
成する拡大画像領域生成手段と、
前記拡大画像領域生成手段とは異なる手法で前記注目領域に対する拡大画像領域を生成
する第２の拡大画像領域生成手段と、
前記拡大画像領域生成手段および前記第２の拡大画像領域生成手段により生成された拡
大画像領域を所定の方法で配置して入力画像に対する拡大画像を生成する拡大画像生成手
段
を有することを特徴とする画像処理システム。
【００１４】
［３］  前記周波数変換手段は、ウェーブレット変換を行う
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００１５】
［４］  前記周波数変換手段は、Ｈａａｒウェーブレットを用いたウェーブレット変換
を行う
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００１６】
［５］  前記周波数成分生成手段は、画像領域の高周波数成分を生成する
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００１７】
［６］  前記周波数成分生成手段は、複数の画像領域から高周波数成分を生成する
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００１８】
［７］  前記周波数成分生成手段は、注目領域の近傍の既定位置の画像領域から高周波
数成分を生成する
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００１９】
［８］  前記周波数成分生成手段は、注目領域の特徴に基づいて選択された画像領域か
ら周波数成分を生成する
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００２０】
［９］  前記周波数成分生成手段は、注目領域の特徴として注目領域のエッジ方向を算
出する
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００２１】
［１０］  前記周波数成分生成手段は、算出した注目領域のエッジ方向に基づいて選択
された画像領域から周波数成分を生成する
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００２２】
［１１］  前記周波数成分生成手段は、算出した注目領域のエッジ方向に基づいた既定
位置の画像領域から周波数成分を生成する
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００２３】
［１２］  前記周波数成分生成手段は、高周波数成分生成に用いる注目領域のエッジ方
向に基づいた複数の画像領域の各々の位置を予め位置情報として備えている
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００２４】
［１３］  前記周波数成分生成手段は、高周波数成分生成に用いる複数の画像領域の既
定位置を注目領域の特徴と相関の強い複数の画像領域の位置とする
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００２５】
［１４］  前記周波数成分生成手段は、画像領域の統計量に基づいて選択された画像領域から周波数成分を生成する
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００２６】
［１５］  前記周波数成分生成手段は、画像領域の統計量として画像領域内の各画素と画像領域の平均値との差分である偏差を用いて、前記偏差に基づいて選択された画像領域から周波数成分を生成する
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００２７】
［１６］  前記周波数成分生成手段は、画像領域の統計量として前記偏差から算出される値を用いて、前記偏差に基づいて選択された画像領域から周波数成分を生成する
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００２８】
［１７］  前記周波数成分生成手段は、画像領域の統計量として画像領域の分散を用いて、前記分散に基づいて選択された画像領域から周波数成分を生成する
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００２９】
［１８］  前記周波数成分生成手段は、複数の画像領域の統計量を算出し、既定の基準に基づいて、複数の画像領域の中から所定の統計量を持つ画像領域を選択し、前記選択された画像領域から周波数成分を生成する
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００３０】
［１９］  前記周波数成分生成手段は、複数の画像領域の統計量を算出し、注目領域の特徴として注目領域のエッジ強度を算出し、前記エッジ強度に基づいて、複数の画像領域の中から所定の統計量を持つ画像領域を選択し、前記選択された画像領域から周波数成分を生成する
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００３１】
［２０］  前記周波数成分生成手段は、複数の画像領域の統計量を算出し、注目領域中の所定画素の画素値に基づいて、複数の画像領域の中から所定の統計量を持つ画像領域を選択し、前記選択された画像領域から周波数成分を生成する
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００３２】
［２１］  前記周波数成分生成手段は、注目領域の高周波数成分に近い周波数成分を持
つ画像領域から高周波数成分を生成する
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００３３】
［２２］  前記周波数成分生成手段は、注目領域の高周波数成分が所定の基準値よりも
小さいときは生成する高周波数成分値を０とする
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００３４】
［２３］  前記拡大画像領域生成手段は、逆ウェーブレット変換を行う
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００３５】
［２４］  前記拡大画像生成手段は、順次生成される前記拡大画像領域を重なり合うよ
うに順次配置する
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００３６】
［２５］  前記拡大画像生成手段は、順次生成される前記拡大画像領域の重なり合う画
素に対して、重なる画素値の総和を算出し、前記画素値の総和を重なった数で割ることに
より前記拡大画像領域の画素値を算出する
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００３７】
［２６］  前記拡大画像生成手段は、順次生成される前記拡大画像領域を重なり合わな
いように順次配置する
ことを特徴とする［１］または［２］に記載の画像処理システム。
【００３８】
［２７］  前記分類手段は、注目領域の周波数成分に基づいて算出される値により注目
領域を分類する
ことを特徴とする［２］に記載の画像処理システム。
【００３９】
［２８］  前記強調手段は、所定のカーネル要素およびカーネル要素間距離をもつエッ
ジ強調カーネルを用いて強調処理を行う
ことを特徴とする［２］に記載の画像処理システム。
【００４０】
［２９］  前記強調手段は、前記分類手段による注目領域の分類結果に基づいて異なる
エッジ強調カーネルを用いて強調処理を行う
ことを特徴とする［２］に記載の画像処理システム。
【００４１】
［３０］  前記拡大画像領域生成手段は、前記分類手段によってテクスチャ領域と分類
された画像領域を前記注目領域とし、
前記第２の拡大画像領域生成手段は、前記分類手段によってエッジ領域と分類された画
像領域を前記注目領域とした
ことを特徴とする［２］に記載の画像処理システム。
【００４２】
［３１］  前記強調手段は、前記分類手段によってテクスチャ領域と分類された画像領
域に対してコントラストの強調を行い、
前記拡大画像領域生成手段は、前記強調手段によってコントラストが強調された画像領
域を前記注目領域とし、
前記第２の拡大画像領域生成手段は、前記分類手段によってエッジ領域と分類された画
像領域を前記注目領域とした
ことを特徴とする［２］に記載の画像処理システム。
【００４３】
［３２］  コンピュータに、［１］ないし［３１］のいずれか１項に記載の画像処
理システムの機能を実現させることを特徴とする画像処理プログラム。
【発明の効果】
【００４４】
本発明にかかる画像処理システムおよび画像処理プログラムによれば、本構成を有して
いない場合に比較して、入力された画像に対してボケやジャギーを生じさせること無く、
鮮鋭度が高い滑らかなエッジの再現と、さらにはテクスチャなどの細部の構造を再現する
ことができる拡大処理を、処理負荷が小さく高速で行うことを可能とすることができる。
【発明を実施するための最良の形態】
【００４５】
本実施の形態に関する技術は、カラーを含む多階調で表現された画像の拡大処理を行う
画像処理技術に関するものであり、特に入力された画像に対してボケやジャギーなどの画
質欠陥をなるべく生じさせること無く、高画質にしかも処理負荷が軽い拡大処理を行う画
像処理技術に関するものである。
【００４６】
（第１の実施例）
図１は、画像処理システム及び画像処理方法の第１の実施の形態を示すブロック図であ
る。図中、１は記憶部、２は画像ブロック設定部、３は拡大画像ブロック生成部、４は拡
大画像生成部である。さらに３０は画像ブロック直交変換部、３１は高周波数成分生成部
、３２は画像ブロック逆直交変換部である。
【００４７】
記憶部１は、画像データが拡大処理されるまで一時的に記憶する機能や、解像度変換ま
たは拡大処理された拡大画像データが図示しない出力装置に出力されるまで一時的に記憶
する機能などを備えたものである。なお、画像データは画像処理システムで処理可能な画
像フォーマット（例えば、ＢＭＰ、ＴＩＦＦ、ＰＮＧなど）で記述されたデータであり、
図示しないデジタルカメラやスキャナなどによって取り込まれたり、パーソナルコンピュ
ータ等において作成、編集等を処理するアプリケーションプログラムで作成された画像デ
ータ等である。拡大画像データについても、同様の画像フォーマットのデータである。
【００４８】
画像ブロック設定部２は、拡大画像ブロック生成部３における処理で必要とされる所定
の画像ブロックサイズをそれぞれ設定し、記憶部１により記憶されている入力画像データ
から、設定されたブロックサイズの画像ブロックを順次（例えばラスタスキャン順に）切
り出し、このブロックサイズの画像ブロックを出力する。
【００４９】
拡大画像ブロック生成部３は、画像ブロック設定部２から出力された画像ブロック中に
含まれる所定サイズの注目ブロックに対して、注目ブロックおよび注目ブロックの近傍に
複数ある同サイズの周辺ブロックをそれぞれ直交変換して複数の高周波数成分要素を抽出
し、前記注目ブロックと前記複数の高周波数成分要素を用いて逆直交変換を行うことによ
り拡大画像ブロックを生成する。
なお、拡大画像ブロック生成部３の詳細については後述する。
【００５０】
拡大画像生成部４は、拡大画像ブロック生成部３で生成された拡大画像ブロックを順次
配置し、解像度変換あるいは拡大された拡大画像データを記憶部１に出力する。なお、拡
大画像ブロックの配置方法については後述するが、例えば順次並べていく方法の他、拡大
画像ブロックを重なり合うように順次配置して重畳させ、重なる画素値の総和を重なった
数で割ることにより画素値を算出するように構成することもできる。
【００５１】
［拡大画像ブロック生成部の詳細説明］
次に、拡大画像ブロック生成部３についてより詳細に説明する。拡大画像ブロック生成
部３は、図１に示すように、画像ブロック直交変換部３０、高周波数成分生成部３１、画
像ブロック逆直交変換部３２で構成される。以降、本例では図２に例示するように注目ブ
ロックおよび周辺ブロックが２×２画素サイズブロックであり、注目ブロックおよび周辺
ブロックを含む画像ブロックが４×４画素サイズブロックである場合を具体例として説明
する。
【００５２】
画像ブロック直交変換部３０は、画像ブロック設定部２から出力された画像ブロック中
の注目ブロックあるいは周辺ブロックに対して、以下の式（１）に示すようにＨａａｒウ
ェーブレットを用いたウェーブレット変換をし、注目ブロックあるいは周辺ブロックをサ
ブバンド成分ＬＬ、ＬＨ、ＨＬ、ＨＨに分解する。
ＬＬ＝（ａ＋ｃ＋ｂ＋ｄ）／４
ＬＨ＝（ａ＋ｂ−ｃ−ｄ）／４
ＨＬ＝（ａ−ｂ＋ｃ−ｄ）／４
ＨＨ＝（ａ−ｂ−ｃ＋ｄ）／４            ・・・式（１）
上記式におけるａ、ｂ、ｃ及びｄは、図３に例示するように、注目ブロックあるいは周
辺ブロックにある各画素の画素値である。
なお、画像ブロック直交変換部３０における直交変換は上述したＨａａｒウェーブレッ
トを用いたウェーブレット変換に限定されるわけではなく、実空間領域の画像ブロックを
サブバンド成分に分解する変換であれば良い。
【００５３】
高周波数成分生成部３１は、画像ブロック逆直交変換部３２で必要とされる注目ブロッ
クの高周波数成分に相当する値を生成する。図４は高周波数成分生成部３１における注目
ブロックの高周波数成分生成処理のフローチャートであり、図５は高周波数成分生成部３
１における注目ブロックの高周波数成分生成処理の具体例を示した図である。
【００５４】
図４に示すように、ステップＳ４０において、高周波数成分生成部３１は、まず図５（
ａ）に示すように注目ブロック近傍の周辺ブロックを選択する。
【００５５】
ステップＳ４１において、高周波数成分生成部３１は選択された周辺ブロックを画像ブ
ロック直交変換部３０に出力し、上述の式（１）を用いて図５（ｂ）に示すようにサブバ
ンド成分ＬＬ１、ＬＨ１、ＨＬ１、ＨＨ１に分解する。
【００５６】
ステップＳ４２において、高周波数成分生成部３１は、図５（ｃ）に示すように前記周
辺ブロックのサブバンド成分のうち高周波数成分を表すＬＨ１、ＨＬ１、ＨＨ１を、図５
（ｄ）に図示するように注目ブロック（画素値ａ、ｂ、ｃ、ｄ）を多重解像度解析におけ
る低周波数成分を表すサブバンド成分ＬＬを構成する要素とみなした時のサブバンド成分
ＬＨ、ＨＬ、ＨＨを構成する所定位置の要素とみなし、ＬＨ１、ＨＬ１、ＨＨ１の値を埋
める。
【００５７】
ステップＳ４３において、高周波数成分生成部３１は、図５（ｃ）に示すように、サブ
バンド成分ＬＨ、ＨＬ、ＨＨを構成する残りすべての要素を生成し終えたかを判断し、全
ての高周波数成分の生成が終了した場合は処理を終了する。高周波数成分の生成が終了し
ていない場合はステップＳ４０からステップＳ４２までの処理を繰り返す。
【００５８】
高周波数成分生成部３１は上述したように、ステップＳ４０からステップＳ４３までの
処理を行うことで４つの周辺ブロックを選択し、最終的に図６に図示するように各注目ブ
ロックに対応した４組の高周波数成分［ＬＨ１、ＨＬ１、ＨＨ１］、［ＬＨ２、ＨＬ２、
ＨＨ２］、［ＬＨ３、ＨＬ３、ＨＨ３］、［ＬＨ４、ＨＬ４、ＨＨ４］をそれぞれ生成す
る。
なお、高周波数成分生成部３１における高周波数成分生成処理については、注目ブロッ
クのサブバンド成分ＬＨ０、ＨＬ０、ＨＨ０の絶対値である｜ＬＨ０｜、｜ＨＬ０｜、｜
ＨＨ０｜の値がすべてある所定の閾値Ｔｈ１より小さい場合のように、注目ブロックの画
素値変化が非常に小さく平坦な場合（注目ブロックに高周波数成分がほとんどない場合）
には、高周波数成分生成部３１は上述したステップＳ４０からステップＳ４３までの高周
波数成分生成処理は行わず、全ての高周波数成分［ＬＨ１、ＨＬ１、ＨＨ１］、［ＬＨ２
、ＨＬ２、ＨＨ２］、［ＬＨ３、ＨＬ３、ＨＨ３］、［ＬＨ４、ＨＬ４、ＨＨ４］を０と
するような処理を行っても良い。
【００５９】
ここで高周波数成分生成部３１における周辺ブロックの選択方法について述べる。
高周波数成分生成部３１における周辺ブロックの選択方法は、例えば図７（ａ）に示す
ように、注目ブロックの真上を第１の周辺ブロック、右を第２の周辺ブロック、左を第３
の周辺ブロック、真下を第４の周辺ブロックと選択する方法や、図７（ｂ）に示すように
、注目ブロックの左上を第１の周辺ブロック、右上を第２の周辺ブロック、左下を第３の
周辺ブロック、右下を第４の周辺ブロックと選択する方法など、注目ブロックの特徴に関
わらずあらかじめ決めておく方法が考えられる。
【００６０】
また、注目ブロックのサブバンド成分ＬＨ０、ＨＬ０を用いて次の式（２）で、注目ブ
ロックのエッジ角度Θを算出し、算出された注目ブロックのエッジ角度Θに基づいて４つ
の周辺ブロックを選択するような動的な選択方法が考えられる。
Θ＝ａｒｃｔａｎ（ＬＨ０／ＨＬ０）        ・・・式（２）
【００６１】
ここで、図１５、図１６および図１７を用いて、注目ブロックのエッジ角度Θに基づく
周辺ブロック選択の具体例について説明する。
式（２）で算出された注目ブロックのエッジ角度Θは、例えば図１５に示すように２２
．５°ごとに区分された８方向に正規化される。本例では、エッジ角度Θが０°又は±１
８０°を中心とした角度範囲を「方向（０）」とし、２２．５°又は−１５７．５°を中
心とした角度範囲を「方向（１）」とし、４５°又は−１３５°を中心とした角度範囲を
「方向（２）」とし、６７．５°又は−１１２．５°を中心とした角度範囲を「方向（３
）」とし、９０°又は−９０°を中心とした角度範囲を「方向（４）」とし、１１２．５
°又は−６７．５°を中心とした角度範囲を「方向（５）」とし、１３５°又は−４５°
を中心とした角度範囲を「方向（６）」とし、１５７．５°又は−２２．５°を中心とし
た角度範囲を「方向（７）」する。これらの角度範囲は、それぞれの中心から±１１．２
５°の範囲である。
【００６２】
ここで、先に述べた正規化されたエッジ角度Θに基づく動的な周辺ブロックの選択とい
うのは、図１６に示すように注目ブロック中の画素ａに対して、４×４画像ブロック内で
注目ブロック中の画素ａを含む周辺ブロックの選択候補である図１６（ａ）、図１６（ｂ
）、図１６（ｃ）、図１６（ｄ）の４つの位置（位置０、位置１、位置２、位置３）にあ
る候補ブロックの中から注目ブロックのエッジ角度Θに基づいて、注目ブロックに近い特
徴（ここではエッジ角度Θ）を有すると推定される１つの周辺ブロックを予め決定してお
くということである。
また他の注目ブロック中の画素ｂ、画素ｃ、画素ｄも画素ａの場合と同様に各画素毎に
注目ブロックのエッジ角度Θに基づいて注目ブロックに近い特徴を有すると推定される周
辺ブロックを予め１ブロック決定しておく。
つまり注目ブロックのエッジ角度Θに基づいて各画素（画素ａ、画素ｂ、画素ｃ、画素
ｄ）毎に選択すべき周辺ブロックの位置を「参照ブロックパターン」として予め用意して
おくことにより、高周波数成分生成部３１における高周波数成分生成で用いる周辺ブロッ
クをエッジ角度Θに基づいて動的に選択する事が可能となる。
【００６３】
具体例として、図１７に注目ブロックのエッジ角度「方向（１）」の場合に選択される
べき周辺ブロックの位置を表した「参照ブロックパターン」の１例を示す。図１７中、斜
めのハッチが掛かった２×２画像ブロックが各画素における高周波数成分生成に用いるべ
き周辺ブロックを示す。
注目ブロック中の画素ａに対して、図１７（Ａ）では選択すべき周辺ブロックを図１６
（ｂ）で示した位置１の画像ブロックとしている。同様に注目ブロック中の画素ｂに対し
ても、図１７（Ｂ）では図１６（ｂ）で示した位置１の画像ブロックとしている。また、
注目ブロック中の画素ｃに対して、図１７（Ｃ）では図１６（ｃ）で示した位置２の画像
ブロックとしている。さらに、注目ブロック中の画素ｄに対しては画素ｃと同様に、図１
７（Ｄ）では図１６（ｃ）で示した位置２の画像ブロックとしている。
つまり図１７に示した具体例においては、エッジ角度「方向（１）」を有する注目ブロ
ック（画素ａ、画素ｂ、画素ｃ、画素ｄ）の各画素が選択すべき周辺ブロックの位置を表
す「参照ブロックパターン」は［位置１、位置１、位置２、位置２］となる。
【００６４】
図１８は、注目ブロックのエッジ角度が「方向（１）」とは異なる場合の「参照ブロッ
クパターン」の具体的な例を示したものである。
図１８（ａ）は注目ブロックのエッジ角度が「方向（０）」の場合であり、「参照ブロ
ックパターン」は[位置１、位置０、位置３、位置２]である。
図１８（ｂ）は注目ブロックのエッジ角度が「方向（２）」の場合であり、「参照ブロ
ックパターン」は[位置３、位置１、位置２、位置０]である。
図１８（ｃ）は注目ブロックのエッジ角度が「方向（３）」の場合であり、「参照ブロ
ックパターン」は[位置２、位置１、位置２、位置１]である。
図１８（ｄ）は注目ブロックのエッジ角度が「方向（４）」の場合であり、「参照ブロ
ックパターン」は[位置２、位置３、位置０、位置１]である。
図１８（ｅ）は注目ブロックのエッジ角度が「方向（５）」の場合であり、「参照ブロ
ックパターン」は[位置０、位置３、位置０、位置３]である。
図１８（ｆ）は注目ブロックのエッジ角度が「方向（６）」の場合であり、「参照ブロ
ックパターン」は[位置０、位置２、位置１、位置３]である。
図１８（ｇ）は注目ブロックのエッジ角度が「方向（７）」の場合であり、「参照ブロ
ックパターン」は[位置０、位置０、位置３、位置３]である。
【００６５】
もちろん注目ブロックのエッジ角度に基づく「参照ブロックパターン」は図１７および
図１８に示した具体例に限られるわけではなく、入力画像データの種類（例えば、風景画
像、人物画像、静物画像など）、あるいは画像フォーマット（例えば、ＢＭＰ、ＪＰＥＧ
、ＴＩＦＦ）などに応じてエッジ角度毎の「参照ブロックパターン」を便宜変更するよう
にしてもよい。
【００６６】
またさらに別の周辺ブロックの選択方法として、注目ブロックの高周波数成分を表すサ
ブバンド成分ＬＨ０、ＨＬ０、ＨＨ０を用いて周辺ブロックとの類似度Δを以下の式（３
）で計算する。
ｄｘ  ＝  ｜ＨＬ０−ＨＬｘ｜
ｄｙ  ＝  ｜ＬＨ０−ＬＨｘ｜
ｄｄ  ＝  ｜ＨＨ０−ＨＨｘ｜
Δ    ＝  ｄｘ＋ｄｙ＋ｄｄ    ・・・式（３）
【００６７】
ここでＨＬｘ、ＬＨｘ、ＨＨｘ、は選択した周辺ブロックのサブバンド成分である。そしてこの類似度Δがある所定の閾値Ｔｈ０よりも小さいブロックを類似周辺ブロックとして選択するような方法をとっても良い。ただしこの場合は、画像ブロック内に類似周辺ブロックが４つ存在するとは限らないので、この場合は類似度Δの値が小さい周辺ブロックから順に選択するようにする。
【００６８】
またさらに、別の周辺ブロックの選択方法として、各周辺ブロックの統計量に相当する値に基づいて選択する方法も考えられる。
図１９（ａ）は注目ブロック（２×２サイズブロック）内のある注目画素、例えば図１６で示した注目領域（ａ、ｂ、ｃ、ｄ）中の注目画素ａが、高周波数成分生成に用いる周辺ブロック（２×２サイズブロック）の選択領域として必要な画像ブロック（３×３サイズブロック）を図示したものである。図１９に示す例では、注目画素は画素i11である。またi00、・・・、i22はそれぞれ各画素およびその画素値を表すものとする。
【００６９】
高周波数成分生成部３１は、図１９（ａ）に示す３×３画像ブロック（i00、・・・、i22）の平均値mean、および図１９（ｂ）から図１９（ｅ）に示す各周辺ブロックである周辺ブロック０（i00、i01、i10、i11）、周辺ブロック１（i01、i02、i11、i12）、周辺ブロック２（i10、i11、i20、i21）、周辺ブロック３（i11、i12、i21、i22）の分散に相当する値、ｖｅｒ０、ｖｅｒ１、ｖｅｒ２、ｖｅｒ３を以下の式（４）で算出する。
mean = (i00 + i01 + i02 + i10 + i11 + i12 + i20 + i21 + i22)/9
ｖｅｒ０= |i00 - mean| + |i01 - mean| + |i10 - mean| + |i11 - mean|
ｖｅｒ１= |i01 - mean| + |i02 - mean| + |i11 - mean| + |i12 - mean|
ｖｅｒ２= |i10 - mean| + |i11 - mean| + |i20 - mean| + |i21 - mean|
ｖｅｒ３= |i11 - mean| + |i12 - mean| + |i21 - mean| + |i22 - mean|
・・・式（４）
【００７０】
次に高周波数成分生成部３１は、式（４）で算出した分散に相当する各値ｖｅｒ０、ｖｅｒ１、ｖｅｒ２、ｖｅｒ３の中から最小値を探索し、その最小値ｖｅｒｍｉｎを持つ周辺ブロックを高周波数成分生成で用いる周辺ブロックとして選択するような方法をとってもよい。上記のような最小値に相当するブロックを選択することにより、高周波数成分生成部３１は破綻の少ない安定した高周波数成分の生成を行うことが可能となる。
ここで周辺ブロックの選択に用いる統計量は、式（４）に示した算出式でのみ算出されるわけではなく、図１９（ｂ）から図１９（ｅ）に示す各周辺ブロックのばらつき具合に相当する値を算出可能な算出式であればどのような式を用いても良い。例えば、次の式（５）のような算出式でももちろん良い。
ｖｅｒ０２=（i00 - mean）２+（i01 - mean）２+（i10 - mean）２+（i11 - mean）２
ｖｅｒ１２=（i01 - mean）２+（i02 - mean）２+（i11 - mean）２+（i12 - mean）２
ｖｅｒ２２=（i10 - mean）２+（i11 - mean）２+（i20 - mean）２+（i21 - mean）２
ｖｅｒ３２=（i11 - mean）２+（i12 - mean）２+（i21 - mean）２+（i22 - mean）２
・・・式（５）
【００７１】
また上述した例では算出した分散に相当する値の最小値を持つ周辺ブロックを高周波数成分生成で用いる周辺ブロックとして選択する例を示したが、最小値をもつ周辺ブロックではなく別の選択方法として、例えば算出した分散相当の値の大小関係と注目ブロックのエッジ強度を関連付けて、前記エッジ強度に基づいて関連付けられた分散値を持つ周辺ブロックを選択する方法も考えられる。
注目ブロックのエッジ強度Ｅは、一例として注目ブロックのサブバンド成分ＬＨ０、ＨＬ０、ＨＨ０を用いて次の式（６）で算出される。
Ｅ＝｜ＬＨ０｜＋｜ＨＬ０｜＋｜ＨＨ０｜                        ・・・式（６）
【００７２】
図２０（ａ）は式（６）で算出されたエッジ強度Ｅの強さにより、処理対象となる注目ブロックのエッジ強度領域を４領域に分割した例を示したものである。図２０（ａ）の領域０が処理対象となるエッジ強度では最もエッジ強度Ｅが弱い領域であり、図２０（ａ）の領域３が処理対象となるエッジ強度では最もエッジ強度Ｅが強い領域である事を示している。
高周波数成分生成部３１は、図２０（ａ）に示した処理対象となる注目ブロックのエッジ強度領域（領域０〜領域３）と式（４）あるいは式（５）を用いて算出された各周辺ブックの分散に相当する値とを関連付ける。具体的には図２０（ｂ）に示すように、式（４）あるいは式（５）で算出された各周辺ブロックの分散に相当する各値の大小関係と分割された注目ブロックのエッジ強度領域を関連付けて、例えば注目ブロックのエッジ強度Ｅが領域３に属する強度を持つ場合は、図２０（ｂ）に図示するように、算出された各周辺ブロックの分散に相当する値を基に最大値を持つ周辺ブロックを選択する。同様に、注目ブロックのエッジ強度Ｅが領域２に属する場合は、２番目に大きい値を持つ周辺ブロックを選択し、注目ブロックのエッジ強度Ｅが領域１に属する場合は、２番目に小さい値を持つ周辺ブロックを選択し、注目ブロックのエッジ強度Ｅが領域０に属する場合は、最小値を持つ周辺ブロックを選択する。このような選択方法をとることにより、高周波数成分生成部３１は、注目ブロックのエッジ強度に応じた高周波数成分生成が可能となり、より入力画像の特徴を高精度に再現する事ができる。
【００７３】
また上記では注目ブロックのエッジ強度Ｅと算出した分散に相当する値を関連付けた例を示したが、図２１に示すように注目ブロック中の注目画素の輝度値と算出した分散に相当する値を関連付けるようにしてもよい。図２１に図示する例では、輝度値が０〜２５５の８ビット階調で表現されているとして、輝度値［０〜６３］の場合は、最小値を持つ周辺ブロックを選択し、輝度値［６４〜１２７］の場合は、２番目に小さい値を持つ周辺ブロックを選択し、輝度値［１２８〜１９１］の場合は、２番目に大きい値を持つ周辺ブロックを選択し、輝度値［１９２〜２５５］の場合は、最大値を持つ周辺ブロックを選択するようにしている。これは入力画像中の輝度値の大きい領域は光の反射率の大きな物体の表面である場合が多く、算出された分散に相当する値が大きな周辺ブロックを選択するのが好ましいためである。
【００７４】
さらに、図２０および図２１に図示したいずれの例でも、領域を４領域に分割しているが、もちろんそれに限ったわけではなく、例えば２領域に分割して分散に相当する値の最大値と最小値を関連付けるといった手法をとってももちろん良く、要するに算出された統計量と注目ブロックあるいは注目画素の状態に基づいて、高周波数成分の生成に必要な周辺ブロックを動的に選択する手法であればどんな手法を用いても良い。
【００７５】
画像ブロック逆直交変換部３２は、図８（ａ）に示すように注目ブロックと高周波数成
分生成部３１で生成された注目ブロックの高周波数成分に相当するサブバンド成分［ＬＨ
１、ＨＬ１、ＨＨ１］、［ＬＨ２、ＨＬ２、ＨＨ２］、［ＬＨ３、ＨＬ３、ＨＨ３］、［
ＬＨ４、ＨＬ４、ＨＨ４］を用いて、以下の式（７）に示すような逆ウェーブレット変換
を行い、図８（ｂ）に示すような注目ブロックに対する拡大画像ブロックを生成する。
ａ０＝（ａ＋ＬＨ１＋ＨＬ１＋ＨＨ１）／４
ｂ０＝（ａ＋ＬＨ１−ＨＬ１−ＨＨ１）／４
ｃ０＝（ａ−ＬＨ１＋ＨＬ１−ＨＨ１）／４
ｄ０＝（ａ−ＬＨ１−ＨＬ１＋ＨＨ１）／４
ａ１＝（ｂ＋ＬＨ２＋ＨＬ２＋ＨＨ２）／４
ｂ１＝（ｂ＋ＬＨ２−ＨＬ２−ＨＨ２）／４
ｃ１＝（ｂ−ＬＨ２＋ＨＬ２−ＨＨ２）／４
ｄ１＝（ｂ−ＬＨ２−ＨＬ２＋ＨＨ２）／４
ａ２＝（ｃ＋ＬＨ３＋ＨＬ３＋ＨＨ３）／４
ｂ２＝（ｃ＋ＬＨ３−ＨＬ３−ＨＨ３）／４
ｃ２＝（ｃ−ＬＨ３＋ＨＬ３−ＨＨ３）／４
ｄ２＝（ｃ−ＬＨ３−ＨＬ３＋ＨＨ３）／４
ａ３＝（ｄ＋ＬＨ４＋ＨＬ４＋ＨＨ４）／４
ｂ３＝（ｄ＋ＬＨ４−ＨＬ４−ＨＨ４）／４
ｃ３＝（ｄ−ＬＨ４＋ＨＬ４−ＨＨ４）／４
ｄ３＝（ｄ−ＬＨ４−ＨＬ４＋ＨＨ４）／４        ・・・式（７）
【００７６】
以上説明したように、拡大画像ブロック生成部３では、注目ブロックの高周波数成分と
して、所定の複数の周辺ブロックを直交変換して生成された高周波数成分を注目ブロック
の高周波数成分と見なし、注目ブロックとこの生成された複数の高周波数成分を用いた逆
直交変換を行うことで拡大画像ブロックを生成することにより、注目ブロックが本来持つ
高周波数成分を失うことなく拡大することが可能となる。
【００７７】
［拡大画像生成部の詳細説明］
次に、拡大画像生成部４をより詳細に説明する。
拡大画像生成部４では、拡大画像ブロック生成部３で生成された注目ブロックに対する
拡大画像ブロックを所定の方法により順次配置する。図９は、拡大画像ブロック生成部３
で生成された４×４画素の拡大画像ブロックを配置する具体例の説明図である。図９に示
す例では、順次生成された拡大画像ブロック０および拡大画像ブロック１をオーバーラッ
プさせるように配置している。オーバーラップする画素は各々前画素値との平均をとるよ
うにして配置する。または、オーバーラップする画素の総和を計算し、前記画素値の総和
をオーバーラップした数で割ることにより各画素値を算出するようにしてもよい。この場
合、注目画素は例えば１画素ずつずらしながら選択し、拡大処理を行っていくことになる
。あるいは、オーバーラップさせずに並べていくこともできる。この場合、注目領域を重
ならないように選択していけばよい。
【００７８】
（第２の実施例）
図１０は、画像処理システム及び画像処理方法の第２の実施の形態を示すブロック図で
ある。
図中、１は記憶部、２は画像ブロック設定部、５は拡大画像ブロック生成部、４は拡大
画像生成部である。さらに３０は画像ブロック直交変換部、３１は高周波数成分生成部、
３２は画像ブロック逆直交変換部、５０は画像ブロック分類部、５１は画像ブロック強調
部、５２は第１拡大ブロック生成部、５３は第２拡大ブロック生成部である。
なお、記憶部１、画像ブロック設定部２、拡大画像生成部４、画像ブロック直交変換部
３０、高周波数成分生成部３１、画像ブロック逆直交変換部３２については図１に示した
第１の実施の形態におけるものと同様であるので説明を省略する。
【００７９】
図１０に例示する第２の実施の形態における画像処理システムの動作の具体例を図１２
のフローチャートを用いて説明する。
ステップＳ１２０において、画像ブロック設定部２は、注目ブロックを含む所定サイズ
の画像ブロックを切り出し、拡大画像ブロック生成部５に出力する。
【００８０】
ステップＳ１２１において、画像ブロック直交変換部３０は、画像ブロック設定部２か
ら出力された画像ブロック中に含まれる所定サイズの注目ブロックを直交変換して各サブ
バンド成分ＬＬ、ＬＨ、ＨＬ、ＨＨに分解する。
【００８１】
ステップＳ１２２において、画像ブロック分類部５０は、前記画像ブロック直交変換部
３０で得られた注目ブロックのサブバンド成分を用いて所定の基準により注目ブロックを
分類する。より具体的には、注目ブロックがエッジを含むブロックであるのか、エッジで
はなく細かな凹凸で構成されるようなテクスチャを含むブロックであるのか判定し、分類
する。
【００８２】
ステップＳ１２３において、画像ブロック強調部５１は、画像ブロック分類部５０にお
いて分類された注目ブロック毎に異なるコントラスト強調処理を行う。
【００８３】
ステップＳ１２４において、第１拡大ブロック生成部５２は、画像ブロック分類部５０
においてテクスチャを含むブロックであると判断された注目ブロックに対して拡大画像ブ
ロック生成処理を行い、また第２拡大画像ブロック生成部５３は、画像ブロック分類部５
０においてエッジを含むブロックであると判断された注目ブロックに対して拡大画像ブロ
ック生成処理を行い、前記注目ブロックの拡大画像ブロックを生成する。その後、拡大画
像生成部４は、第１拡大ブロック生成部５２および第２拡大画像ブロック生成部５３によ
り生成された拡大画像ブロックを配置する。その配置方法は、前記第１の実施の形態にお
ける拡大画像生成部４の配置方法と同様である。
【００８４】
次に、画像ブロック分類部５０、画像ブロック強調部５１、第１拡大ブロック生成部５
２、第２拡大ブロック生成部５３について説明する。
【００８５】
画像ブロック分類部５０は、注目ブロックを画像特徴量に基づいて切り分ける。ここで
画像ブロックの切り分けとは、特徴のある画像ブロック（例えば、エッジを含んだ画像ブ
ロックなど）と、特徴の少ない画像ブロック（例えば、画素値変化の少ない画像ブロック
など）とを切り分けることを意味する。より具体的には、画像ブロック直交変換部３０で
式（１）により直交変換された注目ブロックのサブバンド成分ＬＨ、ＨＬ、ＨＨの絶対値
｜ＬＨ｜、｜ＨＬ｜、｜ＨＨ｜を所定の基準値ＴＨ２とそれぞれ比較し、３つのサブバン
ド成分の絶対値が一つでも基準値ＴＨ２より大きい場合は、注目ブロックをエッジを含む
ブロック（以後、エッジブロックと呼ぶ）と判定し、３つのサブバンド成分の絶対値がす
べて基準値ＴＨ２より小さい場合は、注目ブロックをテクスチャを含むブロック（以後、
テクスチャブロックと呼ぶ）と判定する。
【００８６】
画像ブロック強調部５１は、画像ブロック設定部２から出力された注目ブロックを含む
画像ブロックに対して、先ほど述べた画像ブロック分類部５０における画像ブロックの分
類結果に基づいたコントラスト強調処理を行う。
図１１は、画像ブロック強調部５１において用いられる強調カーネル５３３（コントラ
スト強調カーネル）を例示する図である。
【００８７】
図１１に例示するように、第１の強調カーネル５３３ａは、重み付け係数「１．８０」
及び「−０．２０」を用いてコントラストを強調し、第２の強調カーネル５３３ｂは、重
み付け係数「１．２０」及び「−０．０５」を用いてコントラストを強調する。これらの
強調カーネルは、画像ブロック分類部５０の分類結果に対応付けられており、互いに異な
る重み付け係数を用いてコントラスト強調処理を行う。
具体的には図１１（Ａ）に示す第１の強調カーネル５３３ａは、画像ブロック分類部５
０でエッジブロックと分類された注目ブロックを含む画像ブロックのコントラスト強調処
理に用いられる。このようにより大きな重み付け係数を持つ強調カーネルをエッジブロッ
クに対応付ける事により、より鮮明で鮮鋭度の高いエッジを再現する事ができる。
【００８８】
また図１１（Ｂ）に示す第２の強調カーネル５３３ｂは、画像ブロック分類部５０でテ
クスチャブロックと分類された注目ブロックを含む画像ブロックのコントラスト強調処理
に用いられる。テクスチャブロックに対してエッジブロックより小さな重み付け係数を持
つ強調カーネルを対応付けることにより、テクスチャを際立たせ、テクスチャの再現性を
向上させることができる。またさらにはテクスチャの過度な強調による画像のざらつきを
抑える事も可能である。
【００８９】
第１の強調カーネル５３３ａおよび第２の強調カーネル５３３ｂは、図１１に例示する
ように、注目画素Ｐの直下画素ａ、直右画素ｂ、直上画素ｃ及び直左画素ｄを参照して、
これらの画素の画素値にそれぞれ重み付け係数（−０．２０）または（−０．０５）を掛
け合わせ、重み付け係数（１．８０）または（１．２０）が掛けられた注目画素Ｐの画素
値と合算し、合算された値を注目画素Ｐの画素値とする。
例えば、第１の強調カーネル５３３ａを適用する場合に、以下の式（８）に従って、コ
ントラスト強調後の画素値Ｐ'が算出される。
（画素値Ｐ'）＝１．８０×Ｐ−０．２０×（ａ＋ｂ＋ｃ＋ｄ）・・・式（８）
【００９０】
なお、重み付け係数は図１１（Ａ）、（Ｂ）に例示したものに限られたものではない。
また、参照する画素の位置も、図１１に示すような上下左右に限られたものではない。例
えば、斜め方向の画素を参照してコントラスト強調を行ったり、さらに離れた画素を参照
してコントラスト強調を行ったり、重み付け係数および参照する画素位置それぞれは、処
理対象画像データの種類及びサイズなどにより適宜変更してもよい。
【００９１】
第１拡大ブロック生成部５２は、高周波数成分生成部３１および画像ブロック逆直交変
換部３２で構成され、第１の実施例で述べたように、注目ブロックの高周波数成分として
、所定の複数の周辺ブロックを直交変換して生成された高周波数成分を注目ブロックの高
周波数成分と見なし、注目ブロックとこの生成された複数の高周波数成分を用いた逆直交
変換を行うことで注目ブロックが本来持つ高周波数成分を保存するような拡大画像ブロッ
ク生成処理を行う。
図１０に例示する第２の実施例では、第１拡大ブロック生成部５２は、画像ブロック分
類部５０において、テクスチャブロックと判定された注目ブロックに対して上記の拡大画
像ブロック生成処理を行う。
【００９２】
第２拡大ブロック生成部５３は、画像ブロック分類部５０において、エッジブロックと
判定された注目ブロックに対して、エッジ強度、エッジ方向、注目ブロックの画素値パタ
ーンなど特徴量を算出し、鮮鋭度が高くジャギーの発生を抑制した拡大画像ブロック生成
処理を行う。
【００９３】
上述したように、第２の実施例における画像処理システムおよび画像処理方法では、画
素値差が大きくエッジとして認識できるような入力画像データ中の領域は、第２拡大ブロ
ック生成部５３によってジャギーを抑制した鮮鋭なエッジ再現が可能な拡大処理を行うこ
とができる。さらには画素値差が大きくない領域（たとえば画素値の凹凸があまり大きく
ないような、遠景の木々や、建物の壁面、人物の髪などのテクスチャ領域）においても、
第１拡大ブロック生成部５２で入力画像データの高周波数成分を補うような拡大処理を行
うことにより、補間ボケによる細部情報の欠落などが生じず、上記テクスチャ領域の再現
性を上げることが可能となる。
【００９４】
図１３は、画像処理システムの機能または画像処理方法をコンピュータプログラムで実
現した場合におけるコンピュータプログラム及びそのコンピュータプログラムを格納した
記憶媒体の一例の説明図である。図中、３０１はプログラム、３０２はコンピュータ、３
１１は光磁気ディスク、３１２は光ディスク、３１３は磁気ディスク、３１４はメモリ、
３２１は光磁気ディスク装置、３２２は光ディスク装置、３２３は磁気ディスク装置であ
る。
【００９５】
上述の各実施の形態で説明した画像処理の機能は、コンピュータにより実行可能なプロ
グラム３０１によっても実現することが可能である。その場合、そのプログラム３０１お
よびそのプログラムが用いるデータなどは、コンピュータが読み取り可能な記憶媒体に記
憶することも可能である。記憶媒体とは、コンピュータのハードウェア資源に備えられて
いる読取装置に対して、プログラムの記述内容に応じて、磁気、光、電気等のエネルギー
の変化状態を引き起こして、それに対応する信号の形式で、読取装置にプログラムの記述
内容を伝達できるものである。例えば、光磁気ディスク３１１、光ディスク３１２（ＣＤ
やＤＶＤなどを含む）、磁気ディスク３１３、メモリ３１４（ＩＣカード、メモリカード
などを含む）等である。もちろんこれらの記憶媒体は、可搬型に限られるものではない。
【００９６】
これらの記憶媒体にプログラム３０１を格納しておき、例えばコンピュータ３０２の光
磁気ディスク装置３２１，光ディスク装置３２２，磁気ディスク装置３２３，あるいは図
示しないメモリスロットにこれらの記憶媒体を装着することによって、コンピュータから
プログラム３０１を読み出し、画像処理システムの機能または画像処理方法を実行するこ
とができる。あるいは、予め記憶媒体をコンピュータ３０２に装着しておき、例えばネッ
トワークなどを介してプログラム３０１をコンピュータ３０２に転送し、記憶媒体にプロ
グラム３０１を格納して実行させてもよい。なお、記憶部１は、コンピュータ３０２内の
メモリあるいは付属の磁気ディスク装置やその他の記憶媒体を適用することができる。も
ちろん、本実施の形態の一部の機能についてハードウェアによって構成することもできる
し、あるいは、すべてをハードウェアで構成してもよい。
【００９７】
図１４を参照して、実施の形態の画像処理システムの機能または画像処理方法を実現す
るコンピュータ３０２のハードウェア構成例について説明する。図１４に示す構成は、例
えばパーソナルコンピュータ（ＰＣ）などによって構成されるシステムであり、スキャナ
等のデータ読み取り部４１７と、プリンタなどのデータ出力部４１８を備えたハード構成
例を示している。
【００９８】
ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）４０１は、上述の実施の
形態において説明した各種の構成要素である部分、すなわち、画像ブロック設定部２、拡
大画像ブロック生成部３，５等の各構成要素の実行シーケンスを記述したコンピュータ・
プログラムに従った処理を実行する制御部である。
ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）４０２は、ＣＰＵ４０１が使用するプロ
グラムや演算パラメータ等を格納する。ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏ
ｒｙ）４０３は、ＣＰＵ４０１の実行において使用するプログラムや、その実行において
適宜変化するパラメータ等を格納する。これらはＣＰＵバスなどから構成されるホストバ
ス４０４により相互に接続されている。
【００９９】
ホストバス４０４は、ブリッジ４０５を介して、ＰＣＩ（Ｐｅｒｉｐｈｅｒａｌ  Ｃｏ
ｍｐｏｎｅｎｔ  Ｉｎｔｅｒｃｏｎｎｅｃｔ／Ｉｎｔｅｒｆａｃｅ）バスなどの外部バス
４０６に接続されている。
キーボード４０８、マウス等のポインティングデバイス４０９は、操作者により操作さ
れる入力デバイスである。ディスプレイ４１０は、液晶表示装置またはＣＲＴ（Ｃａｔｈ
ｏｄｅ  Ｒａｙ  Ｔｕｂｅ）などから成り、各種情報をテキストやイメージ情報として表
示する。
【０１００】
ＨＤＤ（Ｈａｒｄ  Ｄｉｓｋ  Ｄｒｉｖｅ）４１１は、ハードディスクを内蔵し、ハー
ドディスクを駆動し、ＣＰＵ４０１によって実行するプログラムや情報を記録または再生
させる。ハードディスクには、入力画像、拡大処理された画像などが格納される。さらに
、その他の各種のデータ処理プログラム等、各種コンピュータ・プログラムが格納される
。
ドライブ４１２は、装着されている磁気ディスク、光ディスク、光磁気ディスク、また
は半導体メモリ等のリムーバブル記録媒体４１３に記録されているデータまたはプログラ
ムを読み出して、そのデータまたはプログラムを、インタフェース４０７、外部バス４０
６、ブリッジ４０５、およびホストバス４０４を介して接続されているＲＡＭ４０３に供
給する。リムーバブル記録媒体４１３も、ハードディスクと同様のデータ記録領域として
利用可能である。
【０１０１】
接続ポート４１４は、外部接続機器４１５を接続するポートであり、ＵＳＢ、ＩＥＥＥ
１３９４等の接続部を持つ。接続ポート４１４は、インタフェース４０７、および外部バ
ス４０６、ブリッジ４０５、ホストバス４０４等を介してＣＰＵ４０１等に接続されてい
る。通信部４１６は、ネットワークに接続され、外部とのデータ通信処理を実行する。デ
ータ読み取り部４１７は、例えばスキャナであり、ドキュメントの読み取り処理を実行す
る。データ出力部４１８は、例えばプリンタであり、ドキュメントデータの出力処理を実
行する。
【０１０２】
なお、図１４に示すシステムのハードウェア構成は、１つの構成例を示すものであり、
本実施の形態は、図１４に示す構成に限らず、本実施の形態において説明した構成要素で
ある部分を実行可能な構成であればよい。例えば、一部の構成要素である部分を専用のハ
ードウェア（例えばＡＳＩＣ等）で構成してもよく、一部の構成要素である部分は外部の
システム内にあり通信回線で接続しているような形態でもよく、さらに図１４に示すシス
テムが複数互いに通信回線によって接続されていて互いに協調動作するようにしてもよい
。また、複写機、ファックス、スキャナ、プリンタ、複合機（多機能複写機とも呼ばれ、
スキャナ、プリンタ、複写機、ファックス等の機能を有している）などに組み込まれてい
てもよい。
【０１０３】
以上説明したように、第１の実施の形態は、所定の大きさの画像領域を複数の周波数成
分に変換し、周波数変換して得られた複数の画像領域の周波数成分から注目領域の高周波
数成分を生成し、注目領域と前記生成された高周波数成分とを逆変換することにより拡大
画像領域を生成してそれらを所定の方法で配置して拡大画像を生成している。このような
拡大処理によって、拡大処理によって失われる原画像に本来ある高周波数成分を補間した
拡大処理が可能となり、ボケやディテール（細部）のつぶれ等の画質欠陥を抑制すること
ができる。
また第２の実施の形態のように、注目領域を分類し、分類された注目領域毎に異なる強
調処理および画像領域拡大処理を行うことで、例えば主にテクスチャなどのボケやディテ
ールのつぶれなどが問題となる領域とエッジなどの鮮鋭性やジャギーなどが問題となる領
域それぞれに最適な画像領域拡大方法を行うことが可能となる。このような拡大処理によ
って、入力画像のボケやディテールのつぶれ、またエッジのジャギーなどの画質欠陥を抑
制した高画質な拡大画像が得ることができる。
【０１０４】
なお、説明したプログラムについては、記録媒体に格納することも可能であり、その場
合は、例えば以下の発明としても把握することができる。
コンピュータに、
所定の大きさの画像領域を複数の周波数成分に変換する周波数変換機能と、
前記周波数変換機能によって変換された前記画像領域の周波数成分から所定の周波数成
分を生成する周波数成分生成機能と、
注目領域を含む所定の大きさの画像領域と前記周波数成分生成機能によって生成された
所定の周波数成分を用いて逆変換することにより前記注目領域に対する拡大画像領域を生
成する拡大画像領域生成機能と、
前記拡大画像領域生成機能により生成された拡大画像領域を所定の方法で配置して入力
画像に対する拡大画像を生成する拡大画像生成機能
を実現させるための画像処理プログラムを記録したコンピュータ読み取り可能な記録媒
体。
【０１０５】
コンピュータに、
所定の大きさの画像領域を複数の周波数成分に変換する周波数変換機能と、
画像領域を分類する分類機能と、
画像領域を強調する強調機能と、
前記周波数変換機能によって変換された前記画像領域の周波数成分から所定の周波数成
分を生成する周波数成分生成機能と、
注目領域を含む所定の大きさの画像領域と前記周波数成分生成機能によって生成された
所定の周波数成分を用いて逆変換することにより前記注目領域に対する拡大画像領域を生
成する拡大画像領域生成機能と、
前記拡大画像領域生成機能とは異なる手法で前記注目領域に対する拡大画像領域を生成
する第２の拡大画像領域生成機能と、
前記拡大画像領域生成機能および前記第２の拡大画像領域生成機能により生成された拡
大画像領域を所定の方法で配置して入力画像に対する拡大画像を生成する拡大画像生成機
能
を実現させるための画像処理プログラムを記録したコンピュータ読み取り可能な記録媒
体。
【０１０６】
「プログラムを記録したコンピュータ読み取り可能な記録媒体」とは、プログラムのイ
ンストール、実行、プログラムの流通などのために用いられる、プログラムが記録された
コンピュータで読み取り可能な記録媒体をいう。
なお、記録媒体としては、例えば、デジタル・バーサタイル・ディスク（ＤＶＤ）であ
って、ＤＶＤフォーラムで策定された規格である「ＤＶＤ−Ｒ、ＤＶＤ−ＲＷ、ＤＶＤ−
ＲＡＭ等」、ＤＶＤ＋ＲＷで策定された規格である「ＤＶＤ＋Ｒ、ＤＶＤ＋ＲＷ等」、コ
ンパクトディスク（ＣＤ）であって、読出し専用メモリ（ＣＤ−ＲＯＭ）、ＣＤレコーダ
ブル（ＣＤ−Ｒ）、ＣＤリライタブル（ＣＤ−ＲＷ）等、光磁気ディスク（ＭＯ）、フレ
キシブルディスク（ＦＤ）、磁気テープ、ハードディスク、読出し専用メモリ（ＲＯＭ）
、電気的消去および書換可能な読出し専用メモリ（ＥＥＰＲＯＭ）、フラッシュ・メモリ
、ランダム・アクセス・メモリ（ＲＡＭ）等が含まれる。
そして、上記のプログラムまたはその一部は、上記記録媒体に記録して保存や流通等さ
せることが可能である。また、通信によって、例えば、ローカル・エリア・ネットワーク
（ＬＡＮ）、メトロポリタン・エリア・ネットワーク（ＭＡＮ）、ワイド・エリア・ネッ
トワーク（ＷＡＮ）、インターネット、イントラネット、エクストラネット等に用いられ
る有線ネットワーク、あるいは無線通信ネットワーク、さらにはこれらの組合せ等の伝送
媒体を用いて伝送することが可能であり、また、搬送波に乗せて搬送することも可能であ
る。
さらに、上記のプログラムは、他のプログラムの一部分であってもよく、あるいは別個
のプログラムと共に記録媒体に記録されていてもよい。
【図面の簡単な説明】
【０１０７】
【図１】画像処理システム及び画像処理方法の第１の実施の形態の一構成例を示すブロック図である。
【図２】注目ブロックおよび周辺ブロックを含む画像ブロックの一例についての説明図である。
【図３】注目ブロックあるいは周辺ブロックにある各画素の画素値の一例についての説明図である。
【図４】高周波数成分生成部における注目ブロックの高周波数成分生成処理例についてのフローチャートである。
【図５】高周波数成分生成部における注目ブロックの高周波数成分生成処理の具体例を示した説明図である。
【図６】各注目ブロックに対応した４組の高周波数成分の具体例を示した説明図である。
【図７】高周波数成分生成部における周辺ブロックの選択方法例を示した説明図である。
【図８】画像ブロック逆直交変換部の処理例を示した説明図である。
【図９】拡大画像ブロック生成部で生成された４×４画素の拡大画像ブロックを配置する具体例の説明図である。
【図１０】画像処理システム及び画像処理方法の第２の実施の形態の一構成例を示すブロック図である。
【図１１】画像ブロック強調部において用いられる強調カーネルを例示する図である。
【図１２】第２の実施の形態における画像処理システムの動作の具体例を示すフローチャートである。
【図１３】画像処理システムの機能または画像処理方法をコンピュータプログラムで実現した場合におけるコンピュータプログラム及びそのコンピュータプログラムを格納した記憶媒体の一例の説明図である。
【図１４】画像処理システムの機能または画像処理方法を実現するコンピュータのハードウェア構成の一例の説明図である。
【図１５】注目ブロックのエッジ角度の具体例を示す説明図である。
【図１６】注目ブロックのエッジ角度に基づく周辺ブロック選択の具体例を示す説明図である。
【図１７】参照ブロックパターンの一例を示す説明図である。
【図１８】参照ブロックパターンの具体的な例を示す説明図である。
【図１９】注目ブロック内の注目画素、周辺ブロックの一例を示す説明図である。
【図２０】エッジ強度Ｅの強さにより、処理対象となる注目ブロックのエッジ強度領域を４領域に分割した例を示す説明図である。
【図２１】注目ブロック中の注目画素の輝度値と算出した分散に相当する値を関連付けた一例を示す説明図である。
【符号の説明】
【０１０８】
１…記憶部、２…画像ブロック設定部、３…拡大画像ブロック生成部、４…拡大画像生
成部、３０…画像ブロック直交変換部、３１…高周波数成分生成部、３２…画像ブロック
逆直交変換部、３０１…プログラム、３０２…コンピュータ、３１１…光磁気ディスク、
３１２…光ディスク、３１３…磁気ディスク、３１４…メモリ、３２１…光磁気ディスク
装置、３２２…光ディスク装置、３２３…磁気ディスク装置。
【図１】
2008059551.tif 000003
【図２】
2008059551.tif 000004
【図３】
2008059551.tif 000005
【図４】
2008059551.tif 000006
【図５】
2008059551.tif 000007
【図６】
2008059551.tif 000008
【図７】
2008059551.tif 000009
【図８】
2008059551.tif 000010
【図９】
2008059551.tif 000011
【図１０】
2008059551.tif 000012
【図１１】
2008059551.tif 000013
【図１２】
2008059551.tif 000014
【図１３】
2008059551.tif 000015
【図１４】
2008059551.tif 000016
【図１５】
2008059551.tif 000017
【図１６】
2008059551.tif 000018
【図１７】
2008059551.tif 000019
【図１８】
2008059551.tif 000020
【図１９】
2008059551.tif 000021
【図２０】
2008059551.tif 000022
【図２１】
2008059551.tif 000023