(57)【特許請求の範囲】
【請求項１】
直線状に配列された複数の画素のうちの欠落画素を補間する画像補間装置において、
前記直線状に配列された複数の画素のうちの前記欠落画素以外の画素は互いに所定のピッチで配列され、前記直線状に配列された複数の画素のうちの前記欠落画素を挟む画素は、相互間に前記所定のピッチの２倍の間隔を有し、
前記欠落画素以外の前記複数の画素の画素値を表すデータを格納するデータ格納部と、
前記データ格納部から読み出されたデータを受け、前記直線状に配列された複数の画素のうち、前記欠落画素の周辺に位置する複数の周辺画素の画素値の各々に、予め定められた係数を掛けて加算する積和演算を行なって、前記欠落画素のための補間値を算出するフィルタ補間部と、
前記フィルタ補間部が出力するデータを前記複数の周辺画素の画素値の範囲内に制限するリミット処理部と、
前記リミット処理部が出力するデータを、前記データ格納部が出力するデータの前記欠落画素に対応する位置に挿入するデータ挿入部と
を備え、
前記フィルタ補間部において、前記周辺画素の画素値の各々に掛けられる係数の大きさは、当該周辺画素の前記欠落画素からの距離をｄ、第１の所定の定数をｇとしたとき、{ｓｉｎ（ｇ×ｄ）｝／（ｇ×ｄ）と第２の所定の定数との積で与えられ、
前記係数は、前記欠落画素に対する係数がゼロである点を除けば、前記フィルタ補間部の前記積和演算を行う部分が低域通過フィルタである場合に用いられる係数に比例し、かつ該係数の和が１となるように定められており、
画素ピッチの逆数をＦｓ、カットオフ空間周波数をＦｃ、各周辺画素の処理対象画素からの距離を画素ピッチ数で表わした値をｓで表わすとき、
前記距離ｄがｄ＝ｓ／Ｆｓで与えられ、
前記第１の所定の定数ｇがｇ＝２π×Ｆｃで与えられ、
前記第２の所定の定数ＫＳ０がＫＳ０＝２Ｆｃ／Ｆｓで与えられる
ことを特徴とする画像補間装置。
【請求項２】
前記リミット処理部は、その入力データが、前記複数の周辺画素の画素値のうちの最大値よりも大きければ、当該入力データに対応する出力データとして、該最大値を出力し、その入力データが、前記複数の周辺画素の画素値のうちの最小値よりも小さければ、当該入力データに対応する出力データとして、該最小値を出力することを特徴とする請求項１に記載の画像補間装置。
【請求項３】
原稿に光に照射する光源と、
複数のセンサチップを直線状に配置し前記原稿からの反射光を光電変換して反射光に対応した信号を出力するセンサチップ列と、
前記センサチップ列が出力する前記反射光に対応した信号をデジタルデータに変換するＡ／Ｄ変換部と、
前記Ａ／Ｄ変換部が出力するデジタルデータを入力とする、請求項１又は２に記載の画像補間装置と
を備え、前記複数のセンサチップの各々には前記所定のピッチで画素が配列され、前記複数のセンサチップのうちの互いに隣接する前記センサチップの端部の画素が前記欠落画素を挟む画素を構成し、前記互いに隣接するセンサチップの端部の画素間に前記欠落画素が存在するものとして、前記フィルタ補間部による前記補間値を算出する処理を行うことを特徴とする画像読取装置。
【請求項４】
直線状に配列された複数の画素のうちの欠落画素を補間する画像補間方法において、
前記直線状に配列された複数の画素のうちの前記欠落画素以外の画素は互いに所定のピッチで配列され、前記直線状に配列された複数の画素のうちの前記欠落画素を挟む画素は、相互間に前記所定のピッチの２倍の間隔を有し、
前記欠落画素以外の前記複数の画素の画素値を表わすデータを格納するデータ格納ステップと、
前記データ格納ステップで格納されたデータを受け、前記直線状に配列された複数の画素のうち、前記欠落画素の周辺に位置する複数の周辺画素の画素値の各々に、予め定められた係数を掛けて加算する積和演算を行なって、前記欠落画素のための補間値を算出するフィルタ補間ステップと、
前記フィルタ補間ステップで補間により得られたデータを前記複数の周辺画素の画素値の範囲内に制限するリミット処理ステップと、
前記リミット処理ステップにより生成されたデータを、前記データ格納ステップにより格納されたデータの前記欠落画素に対応する位置に挿入するデータ挿入ステップと
を備え、
前記フィルタ補間ステップにおいて、前記周辺画素の画素値の各々に掛けられる係数の大きさは、当該周辺画素の前記欠落画素からの距離をｄ、第１の所定の定数をｇとしたとき、{ｓｉｎ（ｇ×ｄ）｝／（ｇ×ｄ）と第２の所定の定数との積で与えられ、
前記係数は、前記欠落画素に対する係数がゼロである点を除けば、前記フィルタ補間ステップの前記積和演算を行う部分が低域通過フィルタ処理を行うものである場合に用いられる係数に比例し、かつ該係数の和が１となるように定められており、
画素ピッチの逆数をＦｓ、カットオフ空間周波数をＦｃ、各周辺画素の処理対象画素からの距離を画素ピッチ数で表わした値をｓで表わすとき、
前記距離ｄがｄ＝ｓ／Ｆｓで与えられ、
前記第１の所定の定数ｇがｇ＝２π×Ｆｃで与えられ、
前記第２の所定の定数ＫＳ０がＫＳ０＝２Ｆｃ／Ｆｓで与えられる
ことを特徴とする画像補間方法。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、スキャナ、ファクシミリ、複写機等を使って原稿をデジタル的に読み取る際、何らかの原因で画像が欠落した場合に画素を補間する画像補間装置に関するものであり、特に複数の撮像素子が直線状に配列されたセンサＩＣチップが複数個直線状に配列された密着イメージセンサのセンサＩＣのつなぎ目部分に画素を補間するのに用い得る画像補間装置及び方法に関するものである。
【背景技術】
【０００２】
上記のような場合に用いられる従来の画像補間装置として、欠落画素に隣接する画素の平均値を補間データとする方法、最小自乗法を使って欠落画素に隣接する画素の回帰直線を求め、その回帰直線から補間データを計算する方法、欠落画素に隣接する４画素から４次式の曲線を求め、その４次式から補間データを計算する方法などを使って欠落画素の補間データを計算するものがあるが、この装置では、周期的なデータを補間するのが困難であった（例えば、特許文献１参照）。
【０００３】
この改善策として、欠落画素を含む複数の画素の平均値と欠落画素を含まない複数の画素の平均値が等しくなるように欠落画素の補間データを求めることにより周期的なデータを補間することも提案されている（例えば、特許文献２参照）。さらに、特性の異なる複数の補間演算回路と採点回路と管理回路と出力回路を備え、複数の補間演算回路はそれぞれ欠落画素および欠落画素近傍の画素の補間データを生成し、採点回路は複数の補間演算回路がそれぞれ出力する欠落画素近傍の画素の補間データと欠落画素近傍の画素の補間データに対応する画素の原データに基づいて採点データを生成し、管理回路は採点回路が出力する採点データに基づいて複数の補間演算回路が出力する欠落画素の補間データのうちから１つを選択するための選択信号を生成し、出力回路は前記管理回路が出力する選択信号に基づいて前記複数の補間演算回路が出力する欠落画素の補間データの１つを欠落画素のデータとして出力することにより補間方法を適応的に切り替えることも提案されている（例えば、特許文献３）。
【０００４】
【特許文献１】特開２００３−１０１７２４号公報（第５−６頁、第３〜５図）
【特許文献２】特開２００５−１４１６１１号公報（第３−７頁、第１図）
【特許文献３】特開２００５−１１７２９１号公報（第７−８頁、第１図）
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、上記の従来の技術では、欠落画素の補間精度、特に、周期的なデータを含む原稿を読み取ったデータを補間する場合の補間精度が十分でないことがあった。
【課題を解決するための手段】
【０００６】
この発明は、上述のような課題を解消するためになされたもので、本発明の画像補間装置は、
直線状に配列された複数の画素のうちの欠落画素を補間する画像補間装置において、
前記直線状に配列された複数の画素のうちの前記欠落画素以外の画素は互いに所定のピッチで配列され、前記直線状に配列された複数の画素のうちの前記欠落画素を挟む画素は、相互間に前記所定のピッチの２倍の間隔を有し、
前記欠落画素以外の前記複数の画素の画素値を表すデータを格納するデータ格納部と、
前記データ格納部から読み出されたデータを受け、前記直線状に配列された複数の画素のうち、前記欠落画素の周辺に位置する複数の周辺画素の画素値の各々に、予め定められた係数を掛けて加算する積和演算を行なって、前記欠落画素のための補間値を算出するフィルタ補間部と、
前記フィルタ補間部が出力するデータを前記複数の周辺画素の画素値の範囲内に制限するリミット処理部と、
前記リミット処理部が出力するデータを、前記データ格納部が出力するデータの前記欠落画素に対応する位置に挿入するデータ挿入部とを備え、
前記フィルタ補間部において、前記周辺画素の画素値の各々に掛けられる係数の大きさは、当該周辺画素の前記欠落画素からの距離をｄ、第１の所定の定数をｇとしたとき、{ｓｉｎ（ｇ×ｄ）｝／（ｇ×ｄ）と第２の所定の定数との積で与えられ、
前記係数は、前記欠落画素に対する係数がゼロである点を除けば、前記フィルタ補間部の前記積和演算を行う部分が低域通過フィルタである場合に用いられる係数に比例し、かつ該係数の和が１となるように定められており、
画素ピッチの逆数をＦｓ、カットオフ空間周波数をＦｃ、各周辺画素の処理対象画素からの距離を画素ピッチ数で表わした値をｓで表わすとき、
前記距離ｄがｄ＝ｓ／Ｆｓで与えられ、
前記第１の所定の定数ｇがｇ＝２π×Ｆｃで与えられ、
前記第２の所定の定数ＫＳ０がＫＳ０＝２Ｆｃ／Ｆｓで与えられることを特徴とする。
【発明の効果】
【０００７】
本発明によれば、カットオフ周波数の高い低域通過フィルタを改良した補間方法を用いたので、カットオフ周波数より低い周期的なデータを良好に補間することができる。
【発明を実施するための最良の形態】
【０００８】
実施の形態１．
図１は、本発明の実施の形態１の画像補間装置の構成を示す図である。図示の画像補間装置１は、データ格納部２、フィルタ補間部３、リミット処理部４、及び補間データ挿入部５を備える。
画像データＤＩがデータ格納部２に入力される。
データ格納部２は、縦続接続された複数のフリップフロップなどを使って欠落画素近傍の一定範囲のデータを格納し、データＤｍとして出力する。
データＤｍは、フィルタ補間部３、リミット処理部４、及び補間データ挿入部５に入力される。
【０００９】
フィルタ補間部３は、積和演算部６と係数格納部７とを有する。係数格納部７は、積和演算部６で用いられる係数を格納する。
積和演算部６は、データ格納部２から読み出されたデータＤｍと係数格納部７が出力する補間用低域通過フィルタの係数ＫＦ１〜ＫＦａを掛けて加算する積和演算を行った結果をデータＤｆとして出力する。
データＤｆは、リミット処理部４に入力される。
リミット処理部４は、データＤｆがデータＤｍの範囲内に収まるように処理を行ったデータをデータＤｃとして出力する。
データＤｃは、補間データ挿入部５に入力される。
補間データ挿入部５は、データＤｍにデータＤｃを挿入したデータをデータＤＯとして出力する。
【００１０】
図２は、本発明の画像補間装置を搭載する密着イメージセンサの構成を示す図である。
密着イメージセンサ８は、光源９、センサチップ列１０、Ａ／Ｄ変換器１１、及び画像補間装置１を備える。
なお、本発明の画像補間装置を密着イメージセンサに搭載し、以下に詳しく述べるように、センサチップ間のつなぎ目を補間する場合について説明を行うが、本発明の画像補間装置は、密着イメージセンサのつなぎ目を補間する場合だけに限定されない。例えば、１次元のラインセンサや２次元のカメラ用のイメージセンサにおいて画素欠陥により読み取れない領域を補間する場合にも有効である。
【００１１】
密着イメージセンサ８の動作について説明する。
光源９が発する光Ｌｓが原稿１２に照射される。
原稿１２からの反射光Ｌｒがセンサチップ列１０に照射される。
センサチップ列１０の各受光素子は、反射光Ｌｒを光電変換し、反射光Ｌｒに対応した信号Ｓａを出力する。
【００１２】
信号ＳａはＡ／Ｄ変換器１１に入力される。
Ａ／Ｄ変換器１１は、信号ＳａをＡ／Ｄ変換してデジタルデータＤＩを出力する。
データＤＩは画像補間装置１に入力される。
画像補間装置１は、センサチップ列１３のつなぎ目に対応する領域を補間したデータＤＯを出力する。
データＤＯは、スキャナや複写機の取り込みデータとして使用される。
【００１３】
センサチップ列と読み取りデータの関係について説明する。
センサチップ列１０は、図３（ａ）に示すように、各々直線的に整列した複数の受光素子（それぞれ画素に対応する）を備える複数のセンサチップを直線状に配置したものであり、従って、複数の受光素子が直線的に、即ち１次元的に配列されたラインセンサを構成するものである。具体的には、センサチップ列１０には、ｋ個（ｋは２以上の整数）のセンサチップ１３−１、１３−２、１３−３、…、１３−ｋが直線状に配置されたものである。センサチップ１３−１には、ｍ個（ｍは２以上の整数）のフォトダイオードなどの受光素子１４−１、１４−２、…、１４−ｍが直線状に配置され、センサチップ１３−２乃至１３−ｋも同様にｍ個の受光素子１５−１乃至１６−ｍが配置されている。
【００１４】
原稿１２は、センサチップ１３−１乃至１３−ｋが並んだ方向に対して垂直な方向に矢印ＤＤで示す方向に搬送され、原稿１２からの反射光Ｌｒが受光素子１４−１乃至１６−ｍによって読み取られる。
【００１５】
しかし、センサチップのつなぎ目部分には受光素子が存在しないため、つなぎ目部分に位置する原稿１５の反射光Ｌｒを読み取ることができない。例えば、センサチップ１３−１の右端の受光素子１４−ｍとセンサチップ１３−２の左端の受光素子１５−１の間はセンサチップの端縁部領域（額縁領域）などであって受光素子が存在しないため、原稿１２からの反射光Ｌｒを読み取ることができない。
【００１６】
そこで、画像補間装置１は、センサチップ１３−１の右端の受光素子１４−ｍとセンサチップ１３−２の左端の受光素子１５−１を含むセンサチップのつなぎ目近傍の受光素子が読み取ったデータに基づいてつなぎ目部分のデータを生成し、読み取ったデータのつなぎ目に対応する部分に挿入する。また、センサチップ１３−２と１３−３のつなぎ目部分などセンサチップ列すべてのつなぎ目部分についてセンサチップ１３−１と１３−２のつなぎ目と同様にデータを補間する。
【００１７】
より詳しく言えば、図３（ｂ）に示すように、センサーチップ１３−１乃至１３−ｋ上の受光素子１４−１乃至１４−ｍは、所定の間隔（画素ピッチ）ｐで直線上に配置されているのに対し、互いに接続される２つのセンサーチップ１３の隣接端部に位置する受光素子は互いに距離ｐ’だけ隔てられている。例えば、センサーチップ１３−１の右端の受光素子１４−ｍとセンサーチップ１３−２の左端の受光素子１５−１の距離がｐ’である。
ここで、ｐ’が、ｐのおおよそ２倍になるようにセンサーチップを配置し、センサーチップ相互間に一つの欠落画素１４−Ｘがあるものとして補間処理をすることにより、高品質な補間結果を得ることができる。そこで、以下ｐ’がｐの２倍である場合について説明するが、本発明は必ずしもこれに限定されるものではなく、距離に応じた補間処理を行うことができる。
【００１８】
次に、画像補間装置１の動作について説明する。
図４（ａ）及び（ｂ）は、データ格納部２の概略構成及び動作を説明する図である。図４（ａ）は、センサチップのつなぎ目近傍の受光素子の位置関係、図４（ｂ）は、データ格納部２の構成を示している。
【００１９】
センサチップ列１０は、１回の読み取りで原稿１２の横方向１ライン分のデータを取得できるが、Ａ／Ｄ変換器１４が１度に１画素分のデータしか処理しないため、データＤＩは１ライン分のデータの左端から順に画像補間装置１に入力される。図４（ａ）では、画素Ｌａ、…、Ｌ３、Ｌ２、Ｌ１、Ｒ１、Ｒ２、Ｒ３、…、Ｒａの順に入力される。センサチップのつなぎ目には画素がないため、位置Ｘのデータを取得できない。上記のうち、Ｌｓ（ｓ＝１〜ａ）は、位置Ｘの画素（欠落画素）に対して左側に位置し、欠落画素Ｘからｓ画素ピッチ離れた位置の画素を表し、Ｒｓ（ｓ＝１〜ａ）は、位置Ｘの画素（欠落画素）に対して右側に位置し、欠落画素Ｘからｓ画素ピッチ離れた位置の画素を表す。
【００２０】
図４（ｂ）に示すように、データ格納部２は、入力されたデータは縦続接続された複数のフリップフロップＦＦＲａ〜ＦＦＬａを使って格納し、フリップフロップに格納しているつなぎ目近傍のデータＤＬａ乃至ＤＲａを出力する。ここで、ＤＬａ乃至ＤＲａは、図１のＤｍに対応し、またそれぞれ画素Ｌａ乃至Ｒａのデータである。
なお、欠落画素Ｘに対応するデータは、データＤＩとしては得られず、データＤＬ１に対応する画素Ｌ１と、データＤＲ１に対応する画素Ｒ１との間には、２画素ピッチ分の間隔がある。
【００２１】
図５は、つなぎ目近傍の受光素子と係数格納部７が格納する係数ＫＦｓ（ｓ＝１、２、…ａ）の一例の関係を示す図である。横軸は受光素子の位置、縦軸は係数ＫＦｓの大きさである。受光素子Ｌ１とＲ１のための係数の大きさはＫＦ１、受光素子Ｌ２とＲ２のための係数の大きさはＫＦ２、受光素子Ｌ３とＲ３のための係数の大きさはＫＦ３、受光素子ＬａとＲａのための係数の大きさはＫＦａである。ここで、係数ＫＦ１乃至ＫＦａは、図１のＫＦ１〜ＫＦａに対応するものであり、これらの集合を符号ＫＦで表すこともある。
なお、欠落画素に対応する係数ＫＦ０の値がゼロとして示されているが、欠落画素のデータは得られないので、積和演算部６では、係数ＫＦ０を用いた乗算は行われない。
図５に示される例では、係数ＫＦｓ（ｓ＝１〜ａ）を表すグラフは、欠落画素Ｘの位置を中心として左右対称であり、データＤＲｓのための係数（ＫＦｓ）とデータＤＬｓのための係数（ＫＦｓ）とは互いに等しい。
【００２２】
フィルタ補間部３は、積和演算を行なう積和演算部６と、積和演算部６で用いられる係数を格納する係数格納部７とを有し、積和演算部６で用いられる係数は、欠落画素に対する係数が０であることを除き、積和演算部６が低域通過フィルタである場合に用いられる係数に比例し、かつ該係数の和が１となるように定められている。
【００２３】
フィルタ補間部（３）において、周辺画素の画素値の各々に掛けられる係数の大きさは、各周辺画素の欠落画素からの距離をｄ、第１の所定の定数をｇとしたとき、{ｓｉｎ（ｇ×ｄ）｝／（ｇ×ｄ）と第２の所定の定数との積で与えられるものとするのが望ましい。
【００２４】
積和演算部６は、データ格納部２が出力するデータＤｍ、すなわちＤＬａ〜ＤＲａと、係数格納部１０が出力する係数ＫＦ１乃至ＫＦａとに対して積和演算を行ってデータＤｆを出力する。データＤｆは、
Ｄｆ
＝ＫＦ１×ＤＬ１＋ＫＦ２×ＤＬ２＋ＫＦ３×ＤＬ３＋
…＋ＫＦａ×ＤＬａ
＋ＫＦ１×ＤＲ１＋ＫＦ２×ＤＲ２＋ＫＦ３×ＤＲ３＋
…＋ＫＦａ×ＤＲａ
で与えられる。
【００２５】
係数格納部７に格納する係数ＫＦ１乃至ＫＦａは、以下に詳しく述べるように、欠落画素乃至補間対象画素に対する係数が０である点を除けば、積和演算部６が、低域通過フィルタである場合に用いられる係数に比例し、かつ該係数の和が１となるように定められている。以下、係数格納部７に格納する係数ＫＦ１乃至ＫＦａの決定方法について説明する。
低域通過フィルタより補間用低域通過フィルタを生成し、補間用低域通過フィルタの係数を係数ＫＦ１乃至ＫＦａとする。
【００２６】
まず、低域通過フィルタについて説明する。
図６は、低域通過フィルタの係数の一例と画素の位置関係を示す図である。横軸は画素の位置（処理対象画素に対する相対位置）、縦軸は当該相対位置の画素のデータに掛けられる係数ＫＳｓの大きさである。着目画素Ｃ０の係数の大きさはＫＳ０、画素Ｌ１とＲ１の係数の大きさはＫＳ１、画素Ｌ２とＲ２の係数の大きさはＫＳ２、画素Ｌ３とＲ３の係数の大きさはＫＳ３、画素ＬａとＲａの係数の大きさはＫＳａである。
図６に示される係数ＫＳｓの値は、下記の式（１）で与えられる。
【００２７】
【数１】
0004772754.tif 000002
式（１）において、
ＫＳ０＝２Ｆｃ／Ｆｓ    …（２）
Ｆｓはサンプリング周波数（空間周波数）であり、画素ピッチｐの逆数（１／ｐ）に等しく、
Ｆｃはカットオフ周波数（空間周波数）である。
ｄは、各周辺画素（係数が掛けられる周辺画素）の処理対象画素からの距離であり、画素ピッチｐと、当該周辺画素の処理対象画素からの距離を画素ピッチ数で表した値ｓとの積、
ｄ＝ｓ×ｐ＝ｓ／Ｆｓ      …（３）
で与えられる。
定数ｇは
ｇ＝２π×Ｆｃ      …（４）
で与えられる。
式（２）、（３）及び（４）を式（１）に代入して式（１）を書き直せば、
【００２８】
【数２】
0004772754.tif 000003
或いは
【数３】
0004772754.tif 000004
となる。
【００２９】
着目画素Ｃ０のデータＤＣ０、着目画素の近傍の画素Ｌａ乃至ＲａのデータＤＬａ乃至ＤＲａと係数ＫＳ０乃至ＫＳａを積和演算すると着目画素Ｃ０の低域通過フィルタ処理後のデータＤｌｐｆ（ＤＣ０）が得られる。データＤｌｐｆ（ＤＣ０）は、
Ｄｌｐｆ（ＤＣ０）
＝ＫＳ０×ＤＣ０
＋ＫＳ１×ＤＬ１＋ＫＳ２×ＤＬ２＋ＫＳ３×ＤＬ３＋
…＋ＫＳａ×ＤＬａ
＋ＫＳ１×ＤＲ１＋ＫＳ２×ＤＲ２＋ＫＳ３×ＤＲ３＋
…＋ＫＳａ×ＤＲａ
となる。
【００３０】
積和演算に用いる係数の総和は１である。すなわち、
ＫＳ０＋ＫＳ１＋ＫＳ２＋ＫＳ３＋
…＋ＫＳａ＋ＫＳ１＋ＫＳ２＋ＫＳ３＋
…＋ＫＳａ
＝１
となる。
【００３１】
次に、補間用低域通過フィルタについて説明する。
図６で示した低域通過フィルタの係数より図５で示したような補間用低域通過フィルタの係数を算出する。図５ではつなぎ目の位置のデータを読み取ることができないため、つなぎ目の位置Ｘのデータを積和演算に使えない。そこで、図６の係数ＫＳ０を０としてつなぎ目の位置Ｘのための係数以外の係数の総和が１となるように係数ＫＦ１乃至ＫＦａを算出する。すなわち、
ＫＦ１＝ＫＳ１／ＳＫＳ
ＫＦ２＝ＫＳ２／ＳＫＳ
ＫＦ３＝ＫＳ３／ＳＫＳ
…
ＫＦａ＝ＫＳａ／ＳＫＳ
一般化すれば、
ＫＦｓ＝ＫＳｓ／ＳＫＳ
（ただし、
ＳＫＳ
＝ＫＳ１＋ＫＳ２＋ＫＳ３＋…＋ＫＳａ＋ＫＳ１＋ＫＳ２＋ＫＳ３＋…＋ＫＳａ
＝２×(ＫＳ１＋ＫＳ２＋ＫＳ３＋…＋ＫＳａ)
である）
により係数ＫＦ１〜ＫＦａ（即ち、ＫＦｓ（ｓ＝１、２、…ａ））を求める。
【００３２】
フィルタ補間部の動作について具体例を用いて説明する。まず、具体的な低域通過フィルタの周波数特性について説明し、次に、具体的な低域通過フィルタの係数より生成した補間用低域通過フィルタの周波数特性を説明する。
【００３３】
図７は、具体的な低域通過フィルタにおいて、着目画素の周辺画素の、着目画素に対する相対位置と、当該周辺画素のデータに掛けられる係数の関係を示す図である。横軸は位置、縦軸は係数の大きさである。着目画素Ｃ０の係数ＫＳ０の大きさを０．７７４とした場合であり、画素Ｌ１、Ｒ１の係数ＫＳ１の大きさは０．１９５、画素Ｌ２、Ｒ２の係数ＫＳ２の大きさは−０．１２４、画素Ｌ３、Ｒ３の係数ＫＳ３の大きさは０．０５４、画素Ｌ４、Ｒ４の係数ＫＳ４の大きさは−０．０１３である。
【００３４】
図７の低域通過フィルタを使って８画素周期、４画素周期、３画素周期、２画素周期のデータに対して行なわれる処理（フィルタリング）について、図８（ａ）乃至（ｃ）、乃至図１１（ａ）乃至（ｃ）を参照して説明する。これらのうち、図８（ａ）、図９（ａ）、図１０（ａ）、及び図１１（ａ）は、互いに同じであり、図７で示した具体的な低域通過フィルタの係数を示す。即ち、これらの図においては、位置ｉの処理対象画素に対する相対位置ｓと、低域通過フィルタにおいて、当該相対位置ｓのデータＤに掛けられる係数ＫＳｓとの関係が示されている。図８（ｂ）及び（ｃ）、図９（ｂ）及び（ｃ）、図１０（ｂ）及び（ｃ）、並びに図１１（ｂ）及び（ｃ）は、画素位置（処理対象画素の位置が符号ｉで表され、その近傍の画素の画素位置が符号ｉ＋ｓ（ｓ＝−４〜４）で表されている）と、その位置のデータＤの大きさとの関係を示す（横軸が画素位置を表し、縦軸がデータの大きさを表す）。
【００３５】
図８（ａ）乃至（ｃ）は、８画素周期のデータを処理する場合を示す。図８（ｂ）は、データの上側のピークを処理する場合、図８（ｃ）は、データの下側のピークを処理する場合を示している。
【００３６】
図８（ｂ）の位置ｉのデータを処理する場合、
Ｄ（ｉ−４）＝０、Ｄ（ｉ−３）＝０．１５、Ｄ（ｉ−２）＝０．５、Ｄ（ｉ−１）＝０．８５、Ｄ（ｉ）＝１、Ｄ（ｉ＋１）＝０．８５、Ｄ（ｉ＋２）＝０．５、Ｄ（ｉ＋３）＝０．１５、Ｄ（ｉ＋４）＝０より、
処理後のデータＤｌｐｆ（ｉ）は、
Ｄｌｐｆ（ｉ）＝−０．０１３×０＋０．０５４×０．１５−０．１２４×０．５＋０．１９５×０．８５
＋０．７７４×１
＋０．１９５×０．８５−０．１２４×０．５＋０．０５４×０．１５−０．０１３×０
＝１
となる。
【００３７】
図８（ｃ）の位置ｉのデータを処理する場合、
Ｄ（ｉ−４）＝１、Ｄ（ｉ−３）＝０．８５、Ｄ（ｉ−２）＝０．５、Ｄ（ｉ−１）＝０．１５、Ｄ（ｉ）＝０、Ｄ（ｉ＋１）＝０．１５、Ｄ（ｉ＋２）＝０．５、Ｄ（ｉ＋３）＝０．８５、Ｄ（ｉ＋４）＝１より、
処理後のデータＤｌｐｆ（ｉ）は、
Ｄｌｐｆ（ｉ）
＝−０．０１３×１＋０．０５４×０．８５−０．１２４×０．５＋０．１９５×０．１５
＋０．７７４×０
＋０．１９５×０．１５−０．１２４×０．５＋０．０５４×０．８５−０．０１３×１
＝０
となる。
【００３８】
従って、データのサンプリング周波数をＦｓとすると、８画素周期のデータのゲインＧ（Ｆｓ／８）は、
Ｇｌｐｆ（Ｆｓ／８）
＝｜処理後のデータの振幅／入力データの振幅｜
＝｜（１−０）／（１−０）｜
＝１
となる。
【００３９】
図９（ａ）乃至（ｃ）は、４画素周期のデータを処理する場合を示す。図９（ｂ）は、データの上側のピークを処理する場合、図９（ｃ）は、データの下側のピークを処理する場合を示している。
【００４０】
図９（ｂ）の位置ｉのデータを処理する場合、
Ｄ（ｉ−４）＝１、Ｄ（ｉ−３）＝０．５、Ｄ（ｉ−２）＝０、Ｄ（ｉ−１）＝０．５、Ｄ（ｉ）＝１、Ｄ（ｉ＋１）＝０．５、Ｄ（ｉ＋２）＝０、Ｄ（ｉ＋３）＝０．５、Ｄ（ｉ＋４）＝１より、
処理後のデータＤｌｐｆ（ｉ）は
Ｄｌｐｆ（ｉ）
＝−０．０１３×１＋０．０５４×０．５−０．１２４×０＋０．１９５×０．５
＋０．７７４×１
＋０．１９５×０．５−０．１２４×０＋０．０５４×０．５−０．０１３×１
＝１
となる。
【００４１】
図９（ｃ）の位置ｉのデータを処理する場合、
Ｄ（ｉ−４）＝０、Ｄ（ｉ−３）＝０．５、Ｄ（ｉ−２）＝１、Ｄ（ｉ−１）＝０．５、Ｄ（ｉ）＝０、Ｄ（ｉ＋１）＝０．５、Ｄ（ｉ＋２）＝１、Ｄ（ｉ＋３）＝０．５、Ｄ（ｉ＋４）＝０より、
処理後のデータＤｌｐｆ（ｉ）は、
Ｄｌｐｆ（ｉ）
＝−０．０１３×０＋０．０５４×０．５−０．１２４×１＋０．１９５×０．５
＋０．７７４×０
＋０．１９５×０．５−０．１２４×１＋０．０５４×０．５−０．０１３×０
＝０
となる。
【００４２】
従って、４画素周期のデータのゲインＧｌｐｆ（Ｆｓ／４）は、
Ｇｌｐｆ（Ｆｓ／４）
＝｜処理後のデータの振幅／入力データの振幅｜
＝｜（１−０）／（１−０）｜
＝１
となる。
【００４３】
図１０（ａ）乃至（ｃ）は、３画素周期のデータを処理する場合を示す。図１０（ｂ）は、データの上側のピークを処理する場合、図１０（ｃ）は、データの下側のピークを処理する場合を示している。
【００４４】
図１０（ｂ）の位置ｉのデータを処理する場合、
Ｄ（ｉ−４）＝０．２５、Ｄ（ｉ−３）＝１、Ｄ（ｉ−２）＝０．２５、Ｄ（ｉ−１）＝０．２５、Ｄ（ｉ）＝１、Ｄ（ｉ＋１）＝０．２５、Ｄ（ｉ＋２）＝０．２５、Ｄ（ｉ＋３）＝１、Ｄ（ｉ＋４）＝０．２５より、
処理後のデータＤｌｐｆ（ｉ）は、
Ｄｌｐｆ（ｉ）
＝−０．０１３×０．２５＋０．０５４×１−０．１２４×０．２５＋０．１９５×０．２５
＋０．７７４×１
＋０．１９５×０．２５−０．１２４×０．２５＋０．０５４×１−０．０１３×０．２５
＝０．９１
となる。
【００４５】
図１０（ｃ）の位置ｉのデータを処理する場合、
Ｄ（ｉ−４）＝０．７５、Ｄ（ｉ−３）＝０、Ｄ（ｉ−２）＝０．７５、Ｄ（ｉ−１）＝０．７５、Ｄ（ｉ）＝０、Ｄ（ｉ＋１）＝０．７５、Ｄ（ｉ＋２）＝０．７５、Ｄ（ｉ＋３）＝０、Ｄ（ｉ＋４）＝０．７５より、
処理後のデータＤｌｐｆ（ｉ）は、
Ｄｌｐｆ（ｉ）
＝−０．０１３×０．７５＋０．０５４×０−０．１２４×０．７５＋０．１９５×０．７５
＋０．７７４×０
＋０．１９５×０．７５−０．１２４×０．７５＋０．０５４×１−０．０１３×０．７５
＝０．０９
となる。
【００４６】
従って、３画素周期のデータのゲインＧｌｐｆ（Ｆｓ／３）は、
Ｇｌｐｆ（Ｆｓ／３）
＝｜処理後のデータの振幅／入力データの振幅｜
＝｜（０．９１−０．０９）／（１−０）｜
＝０．８２
となる。
【００４７】
図１１（ａ）乃至（ｃ）は、２画素周期のデータを処理する場合を示す。図１１（ｂ）は、データの上側のピークを処理する場合、図１１（ｃ）は、データの下側のピークを処理する場合を示している。
【００４８】
図１１（ｂ）の位置ｉのデータを処理する場合、
Ｄ（ｉ−４）＝１、Ｄ（ｉ−３）＝０、Ｄ（ｉ−２）＝１、Ｄ（ｉ−１）＝０、Ｄ（ｉ）＝１、Ｄ（ｉ＋１）＝０、Ｄ（ｉ＋２）＝１、Ｄ（ｉ＋３）＝０、Ｄ（ｉ＋４）＝１より、処理後のデータＤｌｐｆ（ｉ）は、
Ｄｌｐｆ（ｉ）
＝−０．０１３×１＋０．０５４×０−０．１２４×１＋０．１９５×０
＋０．７７４×１
＋０．１９５×０−０．１２４×１＋０．０５４×０−０．０１３×１
＝０．５
となる。
【００４９】
図１１（ｃ）の位置ｉのデータを処理する場合、
Ｄ（ｉ−４）＝０、Ｄ（ｉ−３）＝１、Ｄ（ｉ−２）＝０、Ｄ（ｉ−１）＝１、Ｄ（ｉ）＝０、Ｄ（ｉ＋１）＝１、Ｄ（ｉ＋２）＝０、Ｄ（ｉ＋３）＝１、Ｄ（ｉ＋４）＝０より、処理後のデータＤｌｐｆ（ｉ）は、
Ｄｌｐｆ（ｉ）
＝−０．０１３×０＋０．０５４×１−０．１２４×０＋０．１９５×１
＋０．７７４×０
＋０．１９５×１−０．１２４×０＋０．０５４×１−０．０１３×０
＝０．５
となる。
【００５０】
従って、２画素周期のデータのゲインＧｌｐｆ（Ｆｓ／２）は、
Ｇｌｐｆ（Ｆｓ／２）
＝｜処理後のデータの振幅／入力データの振幅｜
＝｜（０．５−０．５）／（１−０）｜
＝０
となる。
【００５１】
以上、データのピークを処理する場合について説明したが、ピーク以外の部分も同様に処理される。
【００５２】
図１２は、図７の低域通過フィルタの周波数特性を示す図である。横軸は周波数、縦軸はゲインである。図８（ａ）乃至（ｃ）、乃至図１１（ａ）乃至（ｃ）で示したように図７の低域通過フィルタは、８画素周期でゲインＧｌｐｆ（Ｆｓ／８）＝１、４画素周期でゲインＧ（Ｆｓ／４）＝１、３画素周期でゲインＧ（Ｆｓ／３）＝０．８２、２画素周期でゲインＧ（Ｆｓ／２）＝０である。
【００５３】
図１３は、図７の具体的な低域通過フィルタより算出した補間用低域通過フィルタの係数の大きさを示す図である。横軸は欠落画素（補間対象画素）に対する相対位置、縦軸は当該相対位置の画素のデータに掛けられる係数の大きさである。
【００５４】
図７の具体的な低域通過フィルタから補間用低域通過フィルタの係数を決定する方法を説明する。
図７の着目画素Ｃ０の係数を０とし、画素Ｌ１乃至Ｌ４およびＲ１乃至Ｒ４の係数の和が１となるようにする。係数ＫＦ１乃至ＫＦ４は、
ＫＦ１
＝ＫＳ１／（ＫＳ１＋ＫＳ２＋ＫＳ３＋ＫＳ４＋ＫＳ１＋ＫＳ２＋ＫＳ３＋ＫＳ４）
＝０．１９５／（０．１９５−０．１２４＋０．０５４−０．０１３
＋０．１９５−０．１２４＋０．０５４−０．０１３）
＝０．８７０
ＫＦ２
＝ＫＳ２／（ＫＳ１＋ＫＳ２＋ＫＳ３＋ＫＳ４＋ＫＳ１＋ＫＳ２＋ＫＳ３＋ＫＳ４）
＝−０．１２４／（０．１９５−０．１２４＋０．０５４−０．０１３
＋０．１９５−０．１２４＋０．０５４−０．０１３）
＝−０．５５３
ＫＦ３
＝ＫＳ３／（ＫＳ１＋ＫＳ２＋ＫＳ３＋ＫＳ４＋ＫＳ１＋ＫＳ２＋ＫＳ３＋ＫＳ４）
＝０．０５４／（０．１９５−０．１２４＋０．０５４−０．０１３
＋０．１９５−０．１２４＋０．０５４−０．０１３）
＝０．２４１
ＫＦ４
＝ＫＳ４／（ＫＳ１＋ＫＳ２＋ＫＳ３＋ＫＳ４＋ＫＳ１＋ＫＳ２＋ＫＳ３＋ＫＳ４）
＝−０．０１３／（０．１９５−０．１２４＋０．０５４−０．０１３
＋０．１９５−０．１２４＋０．０５４−０．０１３）
＝−０．０５８
となる。
【００５５】
図１３の補間用低域通過フィルタを使って８画素周期、４画素周期、３画素周期、２画素周期のデータのつなぎ目の画素ｉを補間する処理する処理について、図１４（ａ）乃至（ｃ）、乃至図１７（ａ）乃至（ｃ）を参照して説明する。これらのうち、図１４（ａ）、図１５（ａ）、図１６（ａ）、及び図１７（ａ）は、互いに同じであり、図１３で示した具体的な補間用低域通過フィルタの係数ＫＦｓを示す。即ち、これらの図においては、位置ｉの補間対象画素に対する相対位置ｓと、補間用低域通過フィルタにおいて、当該相対位置ｓのデータＤに掛けられる係数ＫＦｓとの関係が示されている。図１４（ｂ）及び（ｃ）、図１５（ｂ）及び（ｃ）、図１６（ｂ）及び（ｃ）、並びに図１７（ｂ）及び（ｃ）は、画素位置（補間対象画素の位置が符号ｉで表され、その近傍の画素の画素位置が符号ｉ＋ｓ（ｓ＝−４〜４）で表されている）と、その位置のデータＤの大きさとの関係を示す（横軸が画素位置を表し、縦軸がデータの大きさを表す）。
【００５６】
図１４（ａ）乃至（ｃ）は、８画素周期のデータを処理する場合を示す。図１４（ｂ）は、データの上側のピークを処理する場合、図１４（ｃ）は、データの下側のピークを処理する場合を示している。
【００５７】
図１４（ｂ）の位置ｉのデータを補間する場合、
Ｄ（ｉ−４）＝０、Ｄ（ｉ−３）＝０．１５、Ｄ（ｉ−２）＝０．５、Ｄ（ｉ−１）＝０．８５、Ｄ（ｉ＋１）＝０．８５、Ｄ（ｉ＋２）＝０．５、Ｄ（ｉ＋３）＝０．１５、Ｄ（ｉ＋４）＝０より、補間データＤｆ（ｉ）は、
Ｄｆ（ｉ）
＝−０．０５８×０＋０．２４１×０．１５−０．５５３×０．５＋０．８７０×０．８５
＋０．８７０×０．８５−０．５５３×０．５＋０．２４１×０．１５−０．０５８×０
＝０．９９８３
となる。
【００５８】
図１４（ｃ）の位置ｉのデータを補間する場合、
Ｄ（ｉ−４）＝１、Ｄ（ｉ−３）＝０．８５、Ｄ（ｉ−２）＝０．５、Ｄ（ｉ−１）＝０．１５、Ｄ（ｉ＋１）＝０．１５、Ｄ（ｉ＋２）＝０．５、Ｄ（ｉ＋３）＝０．８５、Ｄ（ｉ＋４）＝１より、補間データＤｆ（ｉ）は、
Ｄｆ（ｉ）
＝−０．０５８×１＋０．２４１×０．８５−０．５５３×０．５＋０．８７０×０．１５
＋０．８７０×０．１５−０．５５３×０．５＋０．２４１×０．８５−０．０５８×１
＝０．００１７
となる。
【００５９】
従って、８画素周期のデータのゲインＧｆ（Ｆｓ／８）は、
Ｇｆ（Ｆｓ／８）
＝｜処理後のデータの振幅／入力データの振幅｜
＝｜（０．９９８３−０．００１７）／（１−０）｜
＝０．９９６６
となる。
【００６０】
図１５（ａ）乃至（ｃ）は、４画素周期のデータを処理する場合を示す。図１５（ｂ）は、データの上側のピークを処理する場合、図１５（ｃ）は、データの下側のピークを処理する場合を示している。
【００６１】
図１５（ｂ）の位置ｉのデータを補間する場合、
Ｄ（ｉ−４）＝１、Ｄ（ｉ−３）＝０．５、Ｄ（ｉ−２）＝０、Ｄ（ｉ−１）＝０．５Ｄ（ｉ＋１）＝０．５、Ｄ（ｉ＋２）＝０、Ｄ（ｉ＋３）＝０．５、Ｄ（ｉ＋４）＝１より、補間データＤｆ（ｉ）は、
Ｄｆ（ｉ）
＝−０．０５８×１＋０．２４１×０．５−０．５５３×０＋０．８７０×０．５
＋０．８７０×０．５−０．５５３×０＋０．２４１×０．５−０．０５８×１
＝０．９９５
となる。
【００６２】
図１５（ｃ）の位置ｉのデータを補間する場合、
Ｄ（ｉ−４）＝０、Ｄ（ｉ−３）＝０．５、Ｄ（ｉ−２）＝１、Ｄ（ｉ−１）＝０．５、Ｄ（ｉ＋１）＝０．５、Ｄ（ｉ＋２）＝１、Ｄ（ｉ＋３）＝０．５、Ｄ（ｉ＋４）＝０より、補間データＤｆ（ｉ）は、
Ｄｆ（ｉ）
＝−０．０５８×０＋０．２４１×０．５−０．５５３×１＋０．８７０×０．５
＋０．８７０×０．５−０．５５３×１＋０．２４１×０．５−０．０５８×０
＝０．００５
となる。
【００６３】
従って、４画素周期のデータのゲインＧｆ（Ｆｓ／４）は、
Ｇｆ（Ｆｓ／４）
＝｜処理後のデータの振幅／入力データの振幅｜
＝｜（０．９９５−０．００５）／（１−０）｜
＝０．９９
となる。
【００６４】
図１６（ａ）乃至（ｃ）は、３画素周期のデータを処理する場合を示す。図１６（ｂ）は、データの上側のピークを処理する場合、図１６（ｃ）は、データの下側のピークを処理する場合を示している。
【００６５】
図１６（ｂ）の位置ｉのデータを補間する場合、
Ｄ（ｉ−４）＝０．２５、Ｄ（ｉ−３）＝１、Ｄ（ｉ−２）＝０．２５、Ｄ（ｉ−１）＝０．２５、Ｄ（ｉ＋１）＝０．２５、Ｄ（ｉ＋２）＝０．２５、Ｄ（ｉ＋３）＝１、Ｄ（ｉ＋４）＝０．２５より、補間データＤＯ（ｉ）は、
ＤＯ（ｉ）
＝−０．０５８×０．２５＋０．２４１×１−０．５５３×０．２５＋０．８７０×０．２５
＋０．８７０×０．２５−０．５５３×０．２５＋０．２４１×１−０．０５８×０．２５
＝０．６１１５
となる。
【００６６】
図１６（ｃ）の位置ｉのデータを補間する場合、
Ｄ（ｉ−４）＝０．７５、Ｄ（ｉ−３）＝０、Ｄ（ｉ−２）＝０．７５、Ｄ（ｉ−１）＝０．７５、Ｄ（ｉ＋１）＝０．７５、Ｄ（ｉ＋２）＝０．７５、Ｄ（ｉ＋３）＝０、Ｄ（ｉ＋４）＝０．７５より、補間データＤＯ（ｉ）は
ＤＯ（ｉ）
＝−０．０５８×０．７５＋０．２４１×０−０．５５３×０．７５＋０．８７０×０．７５
＋０．８７０×０．７５−０．５５３×０．７５＋０．２４１×０−０．０５８×０．７５
＝０．３８８５
となる。
【００６７】
従って、３画素周期のデータのゲインＧｆ（Ｆｓ／３）は、
Ｇｆ（Ｆｓ／３）
＝｜処理後のデータの振幅／入力データの振幅｜
＝｜（０．６１１５−０．３８８５）／（１−０）｜
＝０．２２３
となる。
【００６８】
図１７（ａ）乃至（ｃ）は、２画素周期のデータを処理する場合を示す。図１７（ｂ）は、データの上側のピークを処理する場合、図１７（ｃ）は、データの下側のピークを処理する場合を示している。
【００６９】
図１７（ｂ）の位置ｉのデータを補間する場合、補間データＤＯ（ｉ−１）を算出すると、Ｄ（ｉ−４）＝１、Ｄ（ｉ−３）＝０、Ｄ（ｉ−２）＝１、Ｄ（ｉ−１）＝０、Ｄ（ｉ＋１）＝０、Ｄ（ｉ＋２）＝１、Ｄ（ｉ＋３）＝０、Ｄ（ｉ＋４）＝１より、補間データＤｆ（ｉ）は、
Ｄｆ（ｉ）
＝−０．０５８×１＋０．２４１×０−０．５５３×１＋０．８７０×０
＋０．８７０×０−０．５５３×１＋０．２４１×０−０．０５８×１
＝−１．２２２
となる。
【００７０】
図１７（ｃ）の位置ｉのデータを補間する場合、
Ｄ（ｉ−４）＝０、Ｄ（ｉ−３）＝１、Ｄ（ｉ−２）＝０、Ｄ（ｉ−１）＝１、Ｄ（ｉ＋１）＝１、Ｄ（ｉ＋２）＝０、Ｄ（ｉ＋３）＝１、Ｄ（ｉ＋４）＝０より、補間データＤｆ（ｉ）は、
Ｄｆ（ｉ）
＝−０．０５８×０＋０．２４１×１−０．５５３×０＋０．８７０×１
＋０．８７０×１−０．５５３×０＋０．２４１×１−０．０５８×０
＝２．２２２
となる。
【００７１】
従って、２画素周期のデータのゲインＧｆ（Ｆｓ／２）は、
Ｇｆ（Ｆｓ／２）
＝｜処理後のデータの振幅／入力データの振幅｜
＝｜（−１．２２２−２．２２２）／（１−０）｜
＝３．４４４
となる。
【００７２】
図１８は、図１３の補間用低域通過フィルタの周波数特性を示す図である。横軸は周波数、縦軸はゲインである。図１４（ａ）乃至（ｃ）、乃至図１７（ａ）乃至（ｃ）で示したように図１３の補間用低域通過フィルタは、８画素周期でゲインＧ（Ｆｓ／８）＝０．９９６６、４画素周期でゲインＧ（Ｆｓ／４）＝０．９９、３画素周期でゲインＧ（Ｆｓ／３）＝０．２２３、２画素周期でゲインＧ（Ｆｓ／２）＝３．４４である。
【００７３】
補間用低域通過フィルタは、図６に示した低域通過フィルタのゲインが１である周波数範囲では、ゲインが１となるが、ゲインが１より小さくなる周波数範囲では図６に示したような低域通過フィルタとは特性が異なる。特に、周波数が高くなるとゲインが１より大きくなる。
【００７４】
フィルタ補間部３は、図６に示したような低域通過フィルタのゲインが１であるような周波数範囲では良好に補間することができる。一方、ゲインが１ではない周波数範囲では良好に補間を行なうことができない。
【００７５】
リミット処理部４の動作について説明する。リミット処理部４はデータＤｍの最大値ｍａｘ（Ｄｍ）と最小値ｍｉｎ（Ｄｍ）を算出し、フィルタ補間部３が出力するデータＤｆがｍａｘ（Ｄｍ）より大きければＤｃ＝ｍａｘ（Ｄｍ）、データＤｆがｍｉｎ（Ｄｍ）より小さければＤｃ＝ｍｉｎ（Ｄｍ）、データＤｆがｍｉｎ（Ｄｍ）以上ｍａｘ（Ｄｍ）以下であればＤｃ＝Ｄｆを出力する。
【００７６】
例えば、図１７（ａ）乃至（ｃ）で示した２画素周期のデータではｍａｘ（Ｄｍ）＝１、ｍｉｎ（Ｄｍ）＝０より、図１７（ｂ）の位置ｉを補間した場合、補間データＤｆ（ｉ）＝２．２２２なのでＤｃ（ｉ）＝１となり、図１７（ｃ）の位置ｉを補間した場合、補間データＤｆ（ｉ）＝−１．２２２よりＤｃ（ｉ）＝０となる。
【００７７】
フィルタ補間部３が出力するデータＤｆをデータＤｍの範囲内に制限することにより、周波数が高いデータを補間する場合でも補間誤差を抑えることができる。
【００７８】
以上、ハードウェアにより補間を行なう場合について説明したが、ソフトウェア又はソフトウェアとハードウェアに組合せにより補間を行なうこととしても良い。
図１９は、そのような構成における処理の手順を示すフローチャートである。
ステップＳＴ１１では、欠落画素近傍のデータＤｍを格納する。
ステップＳＴ１２では、データＤｍと、上記のようにして定めた低域通過フィルタの係数に基づいて定めた係数ＫＦを使って積和演算を行い、データＤｆを生成する。
ステップＳＴ１３では、データＤｆをデータＤｍの範囲内に制限したデータＤｃを生成する。
【００７９】
ステップＳＴ１４では、データＤｍのつなぎ目に対応する位置にデータＤｃを挿入したデータＤＯを生成する。
それぞれのステップにおける処理内容については図１で詳しく示したので、ここでは説明を省略する。
【００８０】
実施の形態２．
図２０は、本発明の実施の形態２の画像補間装置の構成を示す図である。本発明の実施の形態２の画像補間装置は、データ格納部２、フィルタ補間部３、第一の平均維持補間部１７−１、第二の平均維持補間部１７−２、第ｎの平均維持補間部１７−ｎ、選択制御部２１、出力部２０、補間データ挿入部５を備える。選択制御部２１は、採点部１８及び管理部１９を備える。
【００８１】
画像データＤＩがデータ格納部２に入力され、データ格納部２は、フリップフロップなどを使って一定範囲のデータを格納し、一定範囲のデータＤｍを出力する。データＤｍは、フィルタ補間部３、第一の平均維持補間部１７−１、第二の平均維持補間部１７−２、…第ｎの平均維持補間部１７−ｎ、採点部１８に入力される。
【００８２】
フィルタ補間部３は、データＤｍに基づいてテスト補間データＴＤ０と補間候補データＤ０を出力する。テスト補間データＴＤ０は採点部１８に入力され、補間候補データＤ０は出力部２０に入力される。第一の平均維持補間部１７−１は、データＤｍに基づいてテスト補間データＴＤ１と補間候補データＤ１を出力し、第二の平均維持補間部１７−２は、データＤｍに基づいてテスト補間データＴＤ２と補間候補データＤ２を出力し、第ｎの平均維持補間部１７−ｎは、データＤｍに基づいてテスト補間データＴＤｎと補間候補データＤｎを出力する、テスト補間データＴＤ１、ＴＤ２、ＴＤｎは採点部１８に入力され、補間候補データＤ１、Ｄ２、Ｄｎは出力部２０に入力される。
【００８３】
採点部１８はテスト補間データＴＤ０乃至ＴＤｎを、対応するデータＤｍを使って採点し、採点データＭ０乃至Ｍｎを出力する。採点データＭ０乃至Ｍｎは管理部１９に入力される。管理部１９は、採点データＭ０乃至Ｍｎに基づいて補間候補データを選択し、選択信号Ｃを出力する。選択信号Ｃは出力部２０に入力される。出力部２０は、選択信号Ｃに基づいて補間候補データＤ０乃至Ｄｎより選択し、データＤｚを出力する。データＤｚは補間データ挿入部５に入力される。補間データ挿入部５は、データＤｍにデータＤｃを挿入したデータをデータＤＯとして出力する。
【００８４】
平均維持補間部１７−１〜１７−ｎについて説明する。
図２１は平均維持補間部の動作を説明するための図である。平均維持補間部１７−１〜１７−ｎの各々は、欠落画素Ｘを含む画素列ＬＣと欠落画素Ｘを含まない画素列ＮＡおよびＮＢの平均値が等しくなるように欠落画素Ｘのデータを算出する。より詳しく言えば、欠落画素の値を未知数Ｘとして、該欠落画素の周辺の画素から成る第１群の画素（画素列ＬＣに属する画素）の値の平均値と、欠落画素の周辺に位置し、欠落画素以外の画素から成る第２群の画素（画素列ＮＡ、ＮＢに属する画素）の値の平均値とが一致するように、欠落画素の値を定める。すなわち、
（Ｌ１＋…＋Ｌｊ＋Ｘ＋Ｒ１＋…＋Ｒｊ）／ｋ
＝（（Ｌ１＋…＋Ｌｋ）／ｋ＋（Ｒ１＋…＋Ｒｋ）／ｋ）／２
を満たすようにＸを算出する（ｊは正の整数、
ここで、ｋは、画素列ＮＡの画素の数であり、このｋの値は、複数の平均維持補間部１７−１〜１７−ｎ相互間で異なる。例えば第１の平均維持補間部１７−１では、ｋ＝２であり、第２の平均維持補間部１７−２では、ｋ＝３であり、第ｎの平均維持補間部１７−ｎでは、ｋ＝ｎ＋１である。
ｊはｋとの間に
２ｊ＋１＝ｋ、即ち
ｊ＝（ｋ−１）／２
の関係を有する。
【００８５】
具体例を用いて説明する。
図２２は、ｋ＝３の平均維持補間部１７−２の動作を説明するための図である。平均維持補間部１７−２は、欠落画素Ｘを含む画素列ＬＣと欠落画素Ｘを含まない画素列ＮＡおよびＮＢの平均値が等しくなるように欠落画素Ｘのデータを算出する。すなわち、
（Ｌ１＋Ｘ＋Ｒ１）／３
＝（（Ｌ１＋Ｌ２＋Ｌ３）／３＋（Ｒ１＋Ｒ２＋Ｒ３）／３）／２
より、
Ｘ＝（（−Ｌ１＋Ｌ２＋Ｌ３）＋（−Ｒ１＋Ｒ２＋Ｒ３））／２
となる。
【００８６】
図２３は、ｋ＝３の平均維持補間部１７−２を使って３画素周期のデータを処理する場合の参照範囲を示す図である。横軸は位置、縦軸はデータの大きさである。位置ｉのデータを欠落画素Ｘとして処理する場合、位置ｉ−１、ｉ、ｉ＋１を欠落画素を含む画素列ＬＣ、位置ｉ−３、ｉ−２、ｉ−１を欠落画素を含まない画素列ＮＡ、位置ｉ＋１、ｉ＋２、ｉ＋３を欠落画素を含まない画素列ＮＢと設定する。
Ｄ（ｉ−３）＝１、Ｄ（ｉ−２）＝０．２５、Ｄ（ｉ−１）＝０．２５、Ｄ（ｉ＋１）＝０．２５、Ｄ（ｉ＋２）＝０．２５、Ｄ（ｉ＋３）＝１より、
Ｘ＝（（−Ｌ１＋Ｌ２＋Ｌ３）＋（−Ｒ１＋Ｒ２＋Ｒ３））／２
＝（（−Ｄ（ｉ−１）＋Ｄ（ｉ−２）＋Ｄ（ｉ−３））
＋（−Ｄ（ｉ＋１）＋Ｄ（ｉ＋２）＋Ｄ（ｉ＋３）））／２
＝（（−０．２５＋０．２５＋１）＋（−０．２５＋０．２５＋１））／２
＝１
となる。
【００８７】
図２４は、ｋ＝３の平均維持補間部１７−２を使って４画素周期のデータを処理する場合の参照範囲を示す図である。横軸は位置、縦軸はデータの大きさである。位置ｉのデータを欠落画素Ｘとして処理する場合、位置ｉ−１、ｉ、ｉ＋１を欠落画素を含む画素列ＬＣ、位置ｉ−３、ｉ−２、ｉ−１を欠落画素を含まない画素列ＮＡ、位置ｉ＋１、ｉ＋２、ｉ＋３を欠落画素を含まない画素列ＮＢと設定する。
Ｄ（ｉ−３）＝０．５、Ｄ（ｉ−２）＝１、Ｄ（ｉ−１）＝０．５、Ｄ（ｉ＋１）＝０．５、Ｄ（ｉ＋２）＝１、Ｄ（ｉ＋３）＝０．５より、
Ｘ＝（（−Ｌ１＋Ｌ２＋Ｌ３）＋（−Ｒ１＋Ｒ２＋Ｒ３））／２
＝（（−Ｄ（ｉ−１）＋Ｄ（ｉ−２）＋Ｄ（ｉ−３））
＋（−Ｄ（ｉ＋１）＋Ｄ（ｉ＋２）＋Ｄ（ｉ＋３）））／２
＝（（−０．５＋１＋０．５）＋（−０．５＋１＋０．５））／２
＝１
となる。
【００８８】
３画素周期のデータでは位置ｉの期待値データは１であり、ｋ＝３の平均維持補間部１７−２は正しい補間データを算出しているが、４画素周期のデータでは位置ｉの期待値データは０であり、ｋ＝３の平均維持補間部１７−２は正しい補間データを算出していない。このように平均維持補間部１７−２は、ｋの値によって良好に補間できる画素の値の変化の周期（乃至空間周波数）が変わる。逆に言えば、平均維持補間部１７−１〜１７−ｎの各々のｋの値（画素列ＬＣの画素の数が、良好に（即ち、高精度に）補間を行なうことが期待される、画素の値の変化の空間周波数に応じて定められている。
【００８９】
採点部１８の動作について説明する。
採点部１８は、上記のように、テスト補間データＴＤ０〜ＴＤｎと格納部から得られるデータＤｍに基づいて採点を行なう。
テスト補間データＴＤ０〜ＴＤｎは、欠落画素の近傍の非欠落画素について、フィルタ補間部３で算出され、テスト補間データＴＤ１〜ＴＤｎは、平均維持補間部１７−１〜１７−ｎで算出される。
【００９０】
図２５は、テスト補間データを算出する画素（テスト補間画素）と欠落画素の位置関係を示す図である。テスト補間画素を欠落画素の近傍に設定し、フィルタ補間部３及び平均維持補間部１７−１〜１７−ｎの各々は、設定されたテスト補間画素を、仮想的に欠落画素とみなし、欠落画素の補間と同じ方法で補間を行なって補間データを算出し、テスト補間データとする。採点部１８はテスト補間データＴＤ０〜ＴＤとテスト補間画素についての実際の値（データ格納部２に格納されている値）の差分を補間誤差として、補間誤差を表す採点データを出力する。
【００９１】
例えば、Ｌ５をテスト補間画素と設定し、ｋ＝３の平均維持補間部１７−２を使ってテスト補間データを算出する場合、Ｌ４乃至Ｌ６を欠落画素を含む画素列ＬＣ、Ｌ６乃至Ｌ８を欠落画素を含まない画素列ＮＡ、Ｌ２乃至Ｌ４を欠落画素を含まない画素列ＮＢと設定しＬ５の位置の補間データを算出する。
【００９２】
管理部１９の動作について説明する。
図２６は、テスト補間画素Ｔ１乃至Ｔｍと、フィルタ補間部３、第一の平均維持補間部１７−１、第二の平均維持補間部１７−２、第ｎの平均維持補間部１７−ｎの採点データＭ０乃至Ｍｎの関係を示す図である。欠落画素の近傍にテスト補間画素を複数設定した場合、採点部３によってフィルタ補間部３の採点データＭ０［Ｔ１］乃至Ｍ０［Ｔｍ］が算出される。第一の平均維持補間部１７−１乃至第ｎの平均維持補間部１７−ｎも同様に採点データＭ１乃至Ｍｎが算出される。
【００９３】
管理部１９は、フィルタ補間部の採点データを集計し、集計データＳ０を算出する。集計データＳ０は採点データＭ０［Ｔ１］乃至Ｍ０［Ｔｍ］の和である。すなわち、
Ｓ０＝Ｍ０［Ｔ１］＋Ｍ０［Ｔ２］＋…＋Ｍ０［Ｔｍ］
である。
【００９４】
第一の平均維持補間部１７−１乃至第ｎの平均維持補間部１７−ｎの集計データＳ１乃至Ｓｎも同様に
Ｓ１＝Ｍ１［Ｔ１］＋Ｍ１［Ｔ２］＋…＋Ｍ１［Ｔｍ］
Ｓ２＝Ｍ２［Ｔ１］＋Ｍ２［Ｔ２］＋…＋Ｍ２［Ｔｍ］
…
Ｓｎ＝Ｍｎ［Ｔ１］＋Ｍｎ［Ｔ２］＋…＋Ｍｎ［Ｔｍ］
となる。
【００９５】
集計データが最も小さい補間方法、すなわち、欠落画素の近傍での補間誤差（欠落画素の近傍の複数のテスト補間画素についての補間誤差の合計）が最も小さい補間方法を使って欠落画素を補間する。管理部１９は、集計データＳ０乃至Ｓｎのうち、最も小さい集計データを選択し、最も小さい集計データとなる補間方法を示す選択信号Ｃを出力する。
【００９６】
具体例を用いて本発明の実施の形態２の画像補間装置の動作を説明する。
図２７は、補間方法としてフィルタ補間部３、第一の平均維持補間部１７−１、第二の平均維持補間部１７−２を備える本発明の実施の形態２の画像補間装置の構成を示す図である。
【００９７】
第一の平均維持補間部１７−１は、ｋ＝２の平均維持補間部であり、欠落画素を含む２画素の平均値と欠落画素を含まない２画素の平均値が等しくなるように欠落画素のデータを算出する。
【００９８】
第二の平均維持補間部１７−２は、ｋ＝３の平均維持補間部であり、欠落画素を含む３画素の平均値と欠落画素を含まない３画素の平均値が等しくなるように欠落画素のデータを算出する。
【００９９】
図２８は、具体例においてテスト補間画素と欠落画素の位置関係を示す図である。欠落画素Ｘの左側に位置する画素Ｌ６をＴ１、Ｌ５をＴ２、欠落画素の右側に位置する画素Ｒ５をＴ３、Ｒ６をＴ４とする。
【０１００】
図２７の画像補間装置を用いて４画素周期のデータを補間する場合について説明する。
図２９は、位置ｉの画素を補間する場合の４画素周期のデータである。位置ｉ−６をＴ１、ｉ−５をＴ１、ｉ＋５をＴ３、ｉ＋６をＴ４とする。
【０１０１】
フィルタ補間部３は、Ｔ１乃至Ｔ４のテスト補間データＴＤ０と位置ｉの補間候補データＤ０を出力する。フィルタ補間部３の係数格納部７には、図１３の係数が格納されているものとして説明をする。テスト補間データＴＤ０［Ｔ１］乃至ＴＤ０［Ｔ４］は、
ＤＩ（ｉ−１０）＝１、ＤＩ（ｉ−９）＝０．５、ＤＩ（ｉ−８）＝０、ＤＩ（ｉ−７）＝０．５、ＤＩ（ｉ−６）＝１、ＤＩ（ｉ−５）＝０．５、ＤＩ（ｉ−４）＝０、ＤＩ（ｉ−３）＝０．５、ＤＩ（ｉ−２）＝１、ＤＩ（ｉ−１）＝０．５、ＤＩ（ｉ＋１）＝０．５、ＤＩ（ｉ＋２）＝１、ＤＩ（ｉ＋３）＝０．５、ＤＩ（ｉ＋４）＝０、ＤＩ（ｉ＋５）＝０．５、ＤＩ（ｉ＋６）＝１、ＤＩ（ｉ＋７）＝０．５、ＤＩ（ｉ＋８）＝０、ＤＩ（ｉ＋９）＝０．５、ＤＩ（ｉ＋１０）＝１より、
ＴＤ０［Ｔ１］
＝−０．０５８×１＋０．２４１×０．５−０．５５３×０＋０．８７０×０．５
＋０．８７０×０．５−０．５５３×０＋０．２４１×０．５−０．０５８×１
＝０．９９５
ＴＤ０［Ｔ２］
＝−０．０５８×０．５＋０．２４１×０−０．５５３×０．５＋０．８７０×１
＋０．８７０×０−０．５５３×０．５＋０．２４１×１−０．０５８×０．５
＝０．５
ＴＤ０［Ｔ３］
＝−０．０５８×０．５＋０．２４１×１−０．５５３×０．５＋０．８７０×０
＋０．８７０×１−０．５５３×０．５＋０．２４１×０−０．０５８×０．５
＝０．５
ＴＤ０［Ｔ４］
＝−０．０５８×１＋０．２４１×０．５−０．５５３×０＋０．８７０×０．５
＋０．８７０×０．５−０．５５３×０＋０．２４１×０．５−０．０５８×１
＝０．９９５
となる。
【０１０２】
補間候補データＤ０は、
Ｄ０
＝−０．０５８×０＋０．２４１×０．５−０．５５３×１＋０．８７０×０．５
＋０．８７０×０．５−０．５５３×１＋０．２４１×０．５−０．０５８×０
＝０．００５
となる。
【０１０３】
第一の平均維持補間部１７−１は、Ｔ１乃至Ｔ４のテスト補間データＴＤ１と位置ｉの補間候補データＤ１を出力する。テスト補間データＴＤ１［Ｔ１］乃至ＴＤ１［Ｔ４］は、
ＴＤ１［Ｔ１］＝（０＋０．５＋０−０．５）／２＝０
ＴＤ１［Ｔ２］＝（０．５＋１＋０．５−０）／２＝１
ＴＤ１［Ｔ３］＝（０．５＋０＋０．５−１）／２＝０
ＴＤ１［Ｔ４］＝（０＋０．５＋０−０．５）／２＝０
となる。
【０１０４】
補間候補データＤ１は、
Ｄ１＝（１＋０．５＋１−０．５）／２＝１
となる。
【０１０５】
第二の平均維持補間部１７−２は、Ｔ１乃至Ｔ４のテスト補間データＴＤ２と位置ｉの補間候補データＤ２を出力する。テスト補間データＴＤ２［Ｔ１］乃至ＴＤ２［Ｔ４］は、
ＴＤ２［Ｔ１］＝（０．５＋０＋０＋０．５−０．５−０．５）／２＝０
ＴＤ２［Ｔ２］＝（０＋０．５＋０．５＋１−１−０）／２＝０．５
ＴＤ２［Ｔ３］＝（１＋０．５＋０．５＋０−０−１）／２＝０．５
ＴＤ２［Ｔ４］＝（０．５＋０＋０＋０．５−０．５−０．５）／２＝０
となる。
【０１０６】
補間候補データＤ１は、
Ｄ２＝（０．５＋１＋１＋０．５−０．５−０．５）／２＝１
となる。
【０１０７】
図３０は、採点部１８が出力する採点データと管理部１９が生成する集計データをまとめた表である。フィルタ補間部３の採点データＭ０［Ｔ１］乃至Ｍ０［Ｔ４］は、
Ｍ０［Ｔ１］＝｜ＴＤ０［Ｔ１］−ＤＩ（ｉ−６）｜＝｜０．９９５−１｜＝０．００５
Ｍ０［Ｔ２］＝｜ＴＤ０［Ｔ２］−ＤＩ（ｉ−５）｜＝｜０．５−０．５｜＝０
Ｍ０［Ｔ３］＝｜ＴＤ０［Ｔ３］−ＤＩ（ｉ＋５）｜＝｜０．５−０．５｜＝０
Ｍ０［Ｔ４］＝｜ＴＤ０［Ｔ４］−ＤＩ（ｉ＋６）｜＝｜０．９９５−１｜＝０．００５
となる。
【０１０８】
従って、集計データＳ０は、
Ｓ０＝Ｍ０［Ｔ１］＋Ｍ０［Ｔ２］＋Ｍ０［Ｔ３］＋Ｍ０［Ｔ４］
＝０．００５＋０＋０＋０．００５
＝０．０１
となる。
【０１０９】
第一の平均維持補間部１７−１の採点データＭ１［Ｔ１］乃至Ｍ１［Ｔ４］は、
Ｍ１［Ｔ１］＝｜ＴＤ１［Ｔ１］−ＤＩ（ｉ−６）｜＝｜０−１｜＝１
Ｍ１［Ｔ２］＝｜ＴＤ１［Ｔ２］−ＤＩ（ｉ−５）｜＝｜１−０．５｜＝０．５
Ｍ１［Ｔ３］＝｜ＴＤ１［Ｔ３］−ＤＩ（ｉ＋５）｜＝｜０−０．５｜＝０．５
Ｍ１［Ｔ４］＝｜ＴＤ１［Ｔ４］−ＤＩ（ｉ＋６）｜＝｜０−１｜＝１
となる。
【０１１０】
従って、集計データＳ１は、
Ｓ１＝Ｍ１［Ｔ１］＋Ｍ１［Ｔ２］＋Ｍ１［Ｔ３］＋Ｍ１［Ｔ４］
＝１＋０．５＋０．５＋１
＝３
となる。
【０１１１】
第二の平均維持補間部１７−２の採点データＭ２［Ｔ１］乃至Ｍ２［Ｔ４］は、
Ｍ２［Ｔ１］＝｜ＴＤ２［Ｔ１］−ＤＩ（ｉ−６）｜＝｜０−１｜＝１
Ｍ２［Ｔ２］＝｜ＴＤ２［Ｔ２］−ＤＩ（ｉ−５）｜＝｜０．５−０．５｜＝０
Ｍ２［Ｔ３］＝｜ＴＤ２［Ｔ３］−ＤＩ（ｉ＋５）｜＝｜０．５−０．５｜＝０
Ｍ２［Ｔ４］＝｜ＴＤ２［Ｔ４］−ＤＩ（ｉ＋６）｜＝｜０−１｜＝１
となる。
【０１１２】
従って、集計データＳ１は、
Ｓ２＝Ｍ２［Ｔ１］＋Ｍ２［Ｔ２］＋Ｍ２［Ｔ３］＋Ｍ２［Ｔ４］
＝１＋０＋０＋１
＝２
となる。
【０１１３】
管理部１９は集計データが最も小さいフィルタ補間部３を選択し、フィルタ補間部３が出力する補間候補データＤ０を選択することを示す選択信号Ｃを出力する。
【０１１４】
出力部２０は、選択信号Ｃに基づいてフィルタ補間部３が出力する補間候補データＤ０＝０をデータＤｚとして出力する。
【０１１５】
補間データ挿入部５は、位置ｉにデータＤｚ＝０を補間したデータＤＯを出力する。位置ｉの期待値は０なので補間候補データＤ０＝０、Ｄ１＝１、Ｄ２＝１のうち、補間候補データＤ０が最も補間誤差が小さい。
【０１１６】
図２７の画像補間装置を用いて３画素周期のデータを補間する場合について説明する。
図３１は、位置ｉの画素を補間する場合の３画素周期のデータである。位置ｉ−６をＴ１、ｉ−５をＴ１、ｉ＋５をＴ３、ｉ＋６をＴ４とする。
【０１１７】
フィルタ補間部３は、Ｔ１乃至Ｔ４のテスト補間データＴＤ０と位置ｉの補間候補データＤ０を出力する。フィルタ補間部３の係数格納部７には、図１３の係数が格納されているものとして説明をする。テスト補間データＴＤ０［Ｔ１］乃至ＴＤ０［Ｔ４］は、
ＤＩ（ｉ−１０）＝０．７５、ＤＩ（ｉ−９）＝０、ＤＩ（ｉ−８）＝０．７５、ＤＩ（ｉ−７）＝０．７５、ＤＩ（ｉ−６）＝０、ＤＩ（ｉ−５）＝０．７５、ＤＩ（ｉ−４）＝０．７５、ＤＩ（ｉ−３）＝０、ＤＩ（ｉ−２）＝０．７５、ＤＩ（ｉ−１）＝０．７５、ＤＩ（ｉ＋１）＝０．７５、ＤＩ（ｉ＋２）＝０．７５、ＤＩ（ｉ＋３）＝０、ＤＩ（ｉ＋４）＝０．７５、ＤＩ（ｉ＋５）＝０．７５、ＤＩ（ｉ＋６）＝０、ＤＩ（ｉ＋７）＝０．７５、ＤＩ（ｉ＋８）＝０．７５、ＤＩ（ｉ＋９）＝０、ＤＩ（ｉ＋１０）＝０．７５より、
ＴＤ０［Ｔ１］
＝−０．０５８×０．７５＋０．２４１×０−０．５５３×０．７５＋０．８７０×０．７５
＋０．８７０×０．７５−０．５５３×０．７５＋０．２４１×０−０．０５８×０．７５
＝０．３８８５
ＴＤ０［Ｔ２］
＝−０．０５８×０＋０．２４１×０．７５−０．５５３×０．７５＋０．８７０×０
＋０．８７０×０．７５−０．５５３×０＋０．２４１×０．７５−０．０５８×０．７５
＝０．５５５７５
ＴＤ０［Ｔ３］
＝−０．０５８×０．７５＋０．２４１×０．７５−０．５５３×０＋０．８７０×０．７５
＋０．８７０×０−０．５５３×０．７５＋０．２４１×０．７５−０．０５８×０
＝０．５５５７５
ＴＤ０［Ｔ４］
＝−０．０５８×０．７５＋０．２４１×０−０．５５３×０．７５＋０．８７０×０．７５
＋０．８７０×０．７５−０．５５３×０．７５＋０．２４１×０−０．０５８×０．７５
＝０．３８８５
となる。
【０１１８】
補間候補データＤ０は、
Ｄ０
＝−０．０５８×０．７５＋０．２４１×０−０．５５３×０＋．７５＋０．８７０×０．７５
＋０．８７０×０．７５−０．５５３×０．７５＋０．２４１×０−０．０５８×０．７５
＝０．３８８５
となる。
【０１１９】
第一の平均維持補間部１７−１は、Ｔ１乃至Ｔ４のテスト補間データＴＤ１と位置ｉの補間候補データＤ１を出力する。テスト補間データＴＤ１［Ｔ１］乃至ＴＤ１［Ｔ４］は、
ＴＤ１［Ｔ１］＝（０．７５＋０．７５＋０．７５−０．７５）／２＝０．７５
ＴＤ１［Ｔ２］＝（０．７５＋０＋０−０．７５）／２＝０
ＴＤ１［Ｔ３］＝（０＋０．７５＋０．７５−０）／２＝０．７５
ＴＤ１［Ｔ４］＝（０．７５＋０．７５＋０．７５−０．７５）／２＝０．７５
となる。
【０１２０】
補間候補データＤ１は、
Ｄ１＝（０．７５＋０．７５＋０．７５−０．７５）／２＝０．７５
となる。
【０１２１】
第二の平均維持補間部１７−２は、Ｔ１乃至Ｔ４のテスト補間データＴＤ２と位置ｉの補間候補データＤ２を出力する。テスト補間データＴＤ２［Ｔ１］乃至ＴＤ２［Ｔ４］は、
ＴＤ２［Ｔ１］＝（０＋０．７５＋０．７５＋０−０．７５−０．７５）／２＝０
ＴＤ２［Ｔ２］＝（０．７５＋０．７５＋０＋０．７５−０−０．７５）／２＝０．７５
ＴＤ２［Ｔ３］＝（０．７５＋０＋０．７５＋０．７５−０．７５−０）／２＝０．７５
ＴＤ２［Ｔ４］＝（０＋０．７５＋０．７５＋０−０．７５−０．７５）／２＝０
となる。
【０１２２】
補間候補データＤ２は、
Ｄ２＝（０＋０．７５＋０．７５＋０−０．７５−０．７５）／２＝０
となる。
【０１２３】
図３２は、採点部１８が出力する採点データと管理部１９が生成する集計データをまとめた表である。フィルタ補間部３の採点データＭ０［Ｔ１］乃至Ｍ０［Ｔ４］は、
Ｍ０［Ｔ１］＝｜ＴＤ０［Ｔ１］−ＤＩ（ｉ−６）｜
＝｜０．３８８５−０｜＝０．３８８５
Ｍ０［Ｔ２］＝｜ＴＤ０［Ｔ２］−ＤＩ（ｉ−５）｜
＝｜０．５５５７５−０．７５｜＝０．１９４２５
Ｍ０［Ｔ３］＝｜ＴＤ０［Ｔ３］−ＤＩ（ｉ＋５）｜
＝｜０．５５５７５−０．７５｜＝０．１９４２５
Ｍ０［Ｔ４］＝｜ＴＤ０［Ｔ４］−ＤＩ（ｉ＋６）｜
＝｜０．３８８５−０｜＝０．３８８５
となる。
【０１２４】
従って、集計データＳ０は、
Ｓ０＝Ｍ０［Ｔ１］＋Ｍ０［Ｔ２］＋Ｍ０［Ｔ３］＋Ｍ０［Ｔ４］
＝０．３８８５＋０．１９４２５＋０．１９４２５＋０．３８８５
＝１．１６５５
となる。
【０１２５】
第一の平均維持補間部１７−１の採点データＭ１［Ｔ１］乃至Ｍ１［Ｔ４］は、
Ｍ１［Ｔ１］＝｜ＴＤ１［Ｔ１］−ＤＩ（ｉ−６）｜
＝｜０．７５−０｜＝０．７５
Ｍ１［Ｔ２］＝｜ＴＤ１［Ｔ２］−ＤＩ（ｉ−５）｜
＝｜０−０．７５｜＝０．７５
Ｍ１［Ｔ３］＝｜ＴＤ１［Ｔ３］−ＤＩ（ｉ＋５）｜
＝｜０．７５−０．７５｜＝０
Ｍ１［Ｔ４］＝｜ＴＤ１［Ｔ４］−ＤＩ（ｉ＋６）｜
＝｜０．７５−０｜＝０．７５
となる。
【０１２６】
従って、集計データＳ１は、
Ｓ１＝Ｍ１［Ｔ１］＋Ｍ１［Ｔ２］＋Ｍ１［Ｔ３］＋Ｍ１［Ｔ４］
＝０．７５＋０．７５＋０＋０．７５
＝２．２５
となる。
【０１２７】
第二の平均維持補間部１７−２の採点データＭ２［Ｔ１］乃至Ｍ２［Ｔ４］は、
Ｍ２［Ｔ１］＝｜ＴＤ２［Ｔ１］−ＤＩ（ｉ−６）｜＝｜０−０｜＝０
Ｍ２［Ｔ２］＝｜ＴＤ２［Ｔ２］−ＤＩ（ｉ−５）｜＝｜０．７５−０．７５｜＝０
Ｍ２［Ｔ３］＝｜ＴＤ２［Ｔ３］−ＤＩ（ｉ＋５）｜＝｜０．７５−０．７５｜＝０
Ｍ２［Ｔ４］＝｜ＴＤ２［Ｔ４］−ＤＩ（ｉ＋６）｜＝｜０−０｜＝０
となる。
【０１２８】
従って、集計データＳ１は、
Ｓ２＝Ｍ２［Ｔ１］＋Ｍ２［Ｔ２］＋Ｍ２［Ｔ３］＋Ｍ２［Ｔ４］
＝０＋０＋０＋０
＝０
となる。
【０１２９】
管理部１９は、集計データが最も小さい第二の平均維持補間部１７−２を選択し、第二の平均維持補間部１７−２が出力する補間候補データＤ２を選択することを示す選択信号Ｃを出力する。
【０１３０】
出力部２０は、選択信号Ｃに基づいて第二の平均維持補間部１７−２が出力する補間候補データＤ２＝０をデータＤｚとして出力する。
【０１３１】
補間データ挿入部５は、位置ｉにデータＤｚ＝０を補間したデータＤＯを出力する。位置ｉの期待値は０なので補間候補データＤ０＝０．３８８５、Ｄ１＝０．７５、Ｄ２＝０のうち、補間候補データＤ２が最も補間誤差が小さい。
【０１３２】
図２７の画像補間装置を用いて２画素周期のデータを補間する場合について説明する。
図３３は、位置ｉの画素を補間する場合の２画素周期のデータである。位置ｉ−６をＴ１、ｉ−５をＴ１、ｉ＋５をＴ３、ｉ＋６をＴ４とする。
【０１３３】
フィルタ補間部３は、Ｔ１乃至Ｔ４のテスト補間データＴＤ０と位置ｉの補間候補データＤ０を出力する。フィルタ補間部３の係数格納部７には、図１３の係数が格納されているものとして説明をする。テスト補間データＴＤ０［Ｔ１］乃至ＴＤ０［Ｔ４］は、
ＤＩ（ｉ−１０）＝０、ＤＩ（ｉ−９）＝１、ＤＩ（ｉ−８）＝０、ＤＩ（ｉ−７）＝１、ＤＩ（ｉ−６）＝０、ＤＩ（ｉ−５）＝１、ＤＩ（ｉ−４）＝０、ＤＩ（ｉ−３）＝１、ＤＩ（ｉ−２）＝０、ＤＩ（ｉ−１）＝１、ＤＩ（ｉ＋１）＝１、ＤＩ（ｉ＋２）＝０、ＤＩ（ｉ＋３）＝１、ＤＩ（ｉ＋４）＝０、ＤＩ（ｉ＋５）＝１、ＤＩ（ｉ＋６）＝０、ＤＩ（ｉ＋７）＝１、ＤＩ（ｉ＋８）＝０、ＤＩ（ｉ＋９）＝１、ＤＩ（ｉ＋１０）＝０より、
ＴＤ０［Ｔ１］
＝−０．０５８×０＋０．２４１×１−０．５５３×０＋０．８７０×１
＋０．８７０×１−０．５５３×０＋０．２４１×１−０．０５８×０
＝２．２２２
ＴＤ０［Ｔ２］
＝−０．０５８×１＋０．２４１×０−０．５５３×１＋０．８７０×０
＋０．８７０×０−０．５５３×１＋０．２４１×０−０．０５８×１
＝−１．２２２
ＴＤ０［Ｔ３］
＝−０．０５８×１＋０．２４１×０−０．５５３×１＋０．８７０×０
＋０．８７０×０−０．５５３×１＋０．２４１×０−０．０５８×１
＝−１．２２２
ＴＤ０［Ｔ４］
＝−０．０５８×０＋０．２４１×１−０．５５３×０＋０．８７０×１
＋０．８７０×１−０．５５３×０＋０．２４１×１−０．０５８×０
＝２．２２２
となる。
【０１３４】
補間候補データＤ０は、
Ｄ０
＝−０．０５８×０＋０．２４１×１−０．５５３×０＋０．８７０×１
＋０．８７０×１−０．５５３×０＋０．２４１×１−０．０５８×０
＝２．２２２
となる。
【０１３５】
第一の平均維持補間部１７−１は、Ｔ１乃至Ｔ４のテスト補間データＴＤ１と位置ｉの補間候補データＤ１を出力する。テスト補間データＴＤ１［Ｔ１］乃至ＴＤ１［Ｔ４］は、
ＴＤ１［Ｔ１］＝（０＋１＋０−１）／２＝０
ＴＤ１［Ｔ２］＝（１＋０＋１−０）／２＝１
ＴＤ１［Ｔ３］＝（１＋０＋１−０）／２＝１
ＴＤ１［Ｔ４］＝（０＋１＋０−１）／２＝０
となる。
【０１３６】
補間候補データＤ１は、
Ｄ１＝（０＋１＋０−１）／２＝０
となる。
【０１３７】
第二の平均維持補間部１７−２は、Ｔ１乃至Ｔ４のテスト補間データＴＤ２と位置ｉの補間候補データＤ２を出力する。テスト補間データＴＤ２［Ｔ１］乃至ＴＤ２［Ｔ４］は、
ＴＤ２［Ｔ１］＝（１＋０＋０＋１−１−１）／２＝０
ＴＤ２［Ｔ２］＝（０＋１＋１＋０−０−０）／２＝１
ＴＤ２［Ｔ３］＝（０＋１＋１＋０−０−０）／２＝１
ＴＤ２［Ｔ４］＝（１＋０＋０＋１−１−１）／２＝０
となる。
【０１３８】
補間候補データＤ１は、
Ｄ２＝（１＋０＋０＋１−１−１）／２＝０
となる。
【０１３９】
図３４は、採点部１８が出力する採点データと管理部１９が生成する集計データをまとめた表である。フィルタ補間部３の採点データＭ０［Ｔ１］乃至Ｍ０［Ｔ４］は、
Ｍ０［Ｔ１］＝｜ＴＤ０［Ｔ１］−ＤＩ（ｉ−６）｜
＝｜２．２２２−０｜＝２．２２２
Ｍ０［Ｔ２］＝｜ＴＤ０［Ｔ２］−ＤＩ（ｉ−５）｜
＝｜−１．２２２−１｜＝２．２２２
Ｍ０［Ｔ３］＝｜ＴＤ０［Ｔ３］−ＤＩ（ｉ＋５）｜
＝｜−１．２２２−１｜＝２．２２２
Ｍ０［Ｔ４］＝｜ＴＤ０［Ｔ４］−ＤＩ（ｉ＋６）｜
＝｜２．２２２−０｜＝２．２２２
となる。
【０１４０】
従って、集計データＳ０は、
Ｓ０＝Ｍ０［Ｔ１］＋Ｍ０［Ｔ２］＋Ｍ０［Ｔ３］＋Ｍ０［Ｔ４］
＝２．２２２＋２．２２２＋２．２２２＋２．２２２
＝８．８８８
となる。
【０１４１】
第一の平均維持補間部１７−１の採点データＭ１［Ｔ１］乃至Ｍ１［Ｔ４］は、
Ｍ１［Ｔ１］＝｜ＴＤ１［Ｔ１］−ＤＩ（ｉ−６）｜＝｜０−０｜＝０
Ｍ１［Ｔ２］＝｜ＴＤ１［Ｔ２］−ＤＩ（ｉ−５）｜＝｜１−１｜＝０
Ｍ１［Ｔ３］＝｜ＴＤ１［Ｔ３］−ＤＩ（ｉ＋５）｜＝｜１−１｜＝０
Ｍ１［Ｔ４］＝｜ＴＤ１［Ｔ４］−ＤＩ（ｉ＋６）｜＝｜０−０｜＝０
となる。
【０１４２】
従って、集計データＳ１は、
Ｓ１＝Ｍ１［Ｔ１］＋Ｍ１［Ｔ２］＋Ｍ１［Ｔ３］＋Ｍ１［Ｔ４］
＝０＋０＋０＋０
＝０
となる。
【０１４３】
第二の平均維持補間部１７−２の採点データＭ２［Ｔ１］乃至Ｍ２［Ｔ４］は、
Ｍ２［Ｔ１］＝｜ＴＤ２［Ｔ１］−ＤＩ（ｉ−６）｜＝｜０−０｜＝０
Ｍ２［Ｔ２］＝｜ＴＤ２［Ｔ２］−ＤＩ（ｉ−５）｜＝｜１−１｜＝０
Ｍ２［Ｔ３］＝｜ＴＤ２［Ｔ３］−ＤＩ（ｉ＋５）｜＝｜１−１｜＝０
Ｍ２［Ｔ４］＝｜ＴＤ２［Ｔ４］−ＤＩ（ｉ＋６）｜＝｜０−０｜＝０
となる。
【０１４４】
従って、集計データＳ１は、
Ｓ２＝Ｍ２［Ｔ１］＋Ｍ２［Ｔ２］＋Ｍ２［Ｔ３］＋Ｍ２［Ｔ４］
＝０＋０＋０＋０
＝０
となる。
【０１４５】
管理部１９は集計データが最も小さい第一の平均維持補間部１７−２を選択し、第一の平均維持補間部１７−２が出力する補間候補データＤ１を選択することを示す選択信号Ｃを出力する。
【０１４６】
出力部２０は、選択信号Ｃに基づいて第一の平均維持補間部１７−２が出力する補間候補データＤ１＝０をデータＤｚとして出力する。
【０１４７】
補間データ挿入部５は、位置ｉにデータＤｚ＝０を補間したデータＤＯを出力する。位置ｉの期待値は０なので補間候補データＤ０＝２．２２２、Ｄ１＝０、Ｄ２＝０のうち、補間候補データＤ１及びＤ２がＤ０よりも補間誤差が小さい。
【０１４８】
以上のようにして採点部１８と管理部１９とで構成される選択制御部２１が、フィルタ補間部３及び平均維持補間部１７−１〜１７−ｎの各々が欠落画素の近傍に位置し、値が既知である複数の画素の各々をテスト補間画素として、欠落画素の補間値を求めるのと同じ方法で算出したテスト補間値ＴＤ０〜ＴＤｎと、テスト補間画素の既知の値Ｄｍとの差を補間誤差として求め、フィルタ補間部３及び平均維持補間部１７−１〜１７−ｎのうち、複数のテスト補間画素についての補間誤差の合計が最も小さいものを選択し、出力部２０は、フィルタ補間部３及び平均維持補間部１７−１〜１７−ｎのうち、選択制御部２１により選択されたものにより、欠落画素について算出された補間値を、当該欠落画素の補間データとして出力し、データ挿入部５は、出力部２０が出力する補間データをデータ格納部２が出力するデータの欠落画素に対応する位置に挿入する。
【０１４９】
以上のように、図２０の画像補間装置は、フィルタ補間部３が不得意とする（良好な（高精度の）補間を行なうことができない）空間周波数の高いデータを補間するために平均維持補間部１７−１〜１７−ｎを備え、また、欠落画素の近傍の画素を使って補間方法が生成するテスト補間データを採点することにより欠落画素を補間するのに最適な補間方法を選択するため、周波数の低いデータはフィルタ補間部で補間し、周波数の高いデータは平均維持補間部が補間する。よって、周波数全般にわたって良好な補間結果を得ることができる。
【０１５０】
以上、ハードウェアにより補間を行なう場合について説明したが、ソフトウェア又はソフトウェアとハードウェアに組合せにより補間を行なうこととしても良い。
図３５は、そのような構成における処理の手順を示すフローチャートである。
ステップＳＴ２１では、欠落画素近傍のデータＤｍを格納する。
ステップＳＴ２２では、フィルタ補間によりテスト補間データＴＤ０を生成する。
ステップＳＴ２３では、フィルタ補間により補間候補データＤ０を生成する。
【０１５１】
ステップＳＴ２４では、平均維持補間によりテスト補間データＴＤ１乃至ＴＤｎを生成する。
ステップＳＴ２５では、平均維持補間により補間候補データＤ１乃至Ｄｎを生成する。
ステップＳＴ２６では、テスト補間データＴＤ０乃至ＴＤｎを採点し採点データＭ０乃至Ｍｎを生成する。
【０１５２】
ステップＳＴ２７では、採点データＭ０乃至Ｍｎに基づいて補間データを選択する。
ステップＳＴ２８では、データＤｍのつなぎ目に対応する位置にデータＤｃを挿入したデータＤＯを生成する。
それぞれのステップにおける処理内容については図２０で詳しく示したので、ここでは説明を省略する。
【産業上の利用可能性】
【０１５３】
本発明の活用例として、複写機やスキャナなどの密着イメージセンサに適用できる。
【図面の簡単な説明】
【０１５４】
【図１】本発明の実施の形態１の画像補間装置の構成を示す図である。
【図２】本発明の画像補間装置を搭載する密着イメージセンサの構成を示す図である。
【図３】（ａ）及び（ｂ）は、密着イメージセンサのセンサチップ列の構成を示す図である。
【図４】（ａ）及び（ｂ）は、図１のデータ格納部２の概略構成及び動作を示す図である。
【図５】つなぎ目近傍の受光素子と係数格納部７が格納する係数の関係を示す図である。
【図６】低域通過フィルタの係数と画素の位置関係を示す図である。
【図７】具体的な低域通過フィルタで用いられる、着目画素の周辺画素の、着目画素に対する相対位置と、当該周辺画素のデータに掛けられる係数との関係を示す図である。
【図８】（ａ）乃至（ｃ）は、図７の低域通過フィルタを使って８画素周期のデータの位置ｉの画素を処理（フィルタリング）する場合のデータとデータに掛けられる係数の関係を示す図である。
【図９】（ａ）乃至（ｃ）は、図７の低域通過フィルタを使って４画素周期のデータの位置ｉの画素を処理（フィルタリング）する場合のデータとデータに掛けられる係数の関係を示す図である。
【図１０】（ａ）乃至（ｃ）は、図７の低域通過フィルタを使って３画素周期のデータの位置ｉの画素を処理（フィルタリング）する場合のデータとデータに掛けられる係数の関係を示す図である。
【図１１】（ａ）乃至（ｃ）は、図７の低域通過フィルタを使って２画素周期のデータの位置ｉの画素を処理（フィルタリング）する場合のデータとデータに掛けられる係数の関係を示す図である。
【図１２】図７の低域通過フィルタの周波数特性を示す図である。
【図１３】図７の具体的な低域通過フィルタより算出した補間用低域通過フィルタの係数の大きさを示す図である。
【図１４】（ａ）乃至（ｃ）は、図１３の補間用低域通過フィルタを使って８画素周期のデータのつなぎ目の画素ｉを補間する場合のデータとデータに掛けられる係数の関係を示す図である。
【図１５】（ａ）乃至（ｃ）は、図１３の補間用低域通過フィルタを使って４画素周期のデータのつなぎ目の画素ｉを補間する場合のデータとデータに掛けられる係数の関係を示す図である。
【図１６】（ａ）乃至（ｃ）は、図１３の補間用低域通過フィルタを使って３画素周期のデータのつなぎ目の画素ｉを補間する場合のデータとデータに掛けられる係数の関係を示す図である。
【図１７】（ａ）乃至（ｃ）は、図１３の補間用低域通過フィルタを使って２画素周期のデータのつなぎ目の画素ｉを補間する場合のデータとデータに掛けられる係数の関係を示す図である。
【図１８】図１３の補間用低域通過フィルタの周波数特性を示す図である。
【図１９】ソフトウェア処理による動作を説明するフローチャートである。
【図２０】本発明の実施の形態２の画像補間装置の構成を示す図である。
【図２１】平均維持補間部の動作を説明するための図である。
【図２２】ｋ＝３の場合の平均維持補間部の動作を説明するための図である。
【図２３】ｋ＝３の平均維持補間部を使って３画素周期のデータを処理する場合の参照範囲を示す図である。
【図２４】ｋ＝３の平均維持補間部を使って４画素周期のデータを処理する場合の参照範囲を示す図である。
【図２５】テスト補間画素と欠落画素の位置関係を示す図である。
【図２６】テスト補間画素Ｔ１乃至Ｔｍと、フィルタ補間部、第一の平均維持補間部、第二の平均維持補間部、第ｎの平均維持補間部の採点データＭ０乃至Ｍｎの関係を示す図である。
【図２７】補間方法としてフィルタ補間部、第一の平均維持補間部、第二の平均維持補間部を備える本発明の実施の形態２の画像補間装置の構成を示す図である。
【図２８】具体例においてテスト補間画素と欠落画素の位置関係を示す図である。
【図２９】位置ｉの画素を補間する場合の４画素周期のデータである。
【図３０】採点部が出力する採点データと管理部が生成する集計データをまとめた表である。
【図３１】位置ｉの画素を補間する場合の３画素周期のデータである。
【図３２】採点部が出力する採点データと管理部が生成する集計データをまとめた表である。
【図３３】位置ｉの画素を補間する場合の２画素周期のデータである。
【図３４】採点部が出力する採点データと管理部が生成する集計データをまとめた表である。
【図３５】ソフトウェア処理による動作を説明するフローチャートである。
【符号の説明】
【０１５５】
１  画像補間装置、  ２  データ格納部、  ３  フィルタ補間部、  ４  リミット処理部、  ５  補間データ挿入部、  ６  積和演算部、  ７  係数格納部、  ８  密着イメージセンサ、  ９  光源、  １０  センサチップ列、  １１  Ａ／Ｄ変換器、  １２  原稿、  １３−１〜１３−ｋ  センサチップ、  １４−１〜１４−ｍ  左端のセンサチップの受光素子、  １５−１  左端から２個目のセンサチップの左端の受光素子、  １６−ｍ  右端のセンサチップの右端の受光素子、  １７−１〜１７−ｎ  平均維持補間部、  １８  採点部、  １９  管理部、  ２０  出力部。
【図１】
0004772754.tif 000005
【図２】
0004772754.tif 000006
【図３】
0004772754.tif 000007
【図４】
0004772754.tif 000008
【図５】
0004772754.tif 000009
【図６】
0004772754.tif 000010
【図７】
0004772754.tif 000011
【図８】
0004772754.tif 000012
【図９】
0004772754.tif 000013
【図１０】
0004772754.tif 000014
【図１１】
0004772754.tif 000015
【図１２】
0004772754.tif 000016
【図１３】
0004772754.tif 000017
【図１４】
0004772754.tif 000018
【図１５】
0004772754.tif 000019
【図１６】
0004772754.tif 000020
【図１７】
0004772754.tif 000021
【図１８】
0004772754.tif 000022
【図１９】
0004772754.tif 000023
【図２０】
0004772754.tif 000024
【図２１】
0004772754.tif 000025
【図２２】
0004772754.tif 000026
【図２３】
0004772754.tif 000027
【図２４】
0004772754.tif 000028
【図２５】
0004772754.tif 000029
【図２６】
0004772754.tif 000030
【図２７】
0004772754.tif 000031
【図２８】
0004772754.tif 000032
【図２９】
0004772754.tif 000033
【図３０】
0004772754.tif 000034
【図３１】
0004772754.tif 000035
【図３２】
0004772754.tif 000036
【図３３】
0004772754.tif 000037
【図３４】
0004772754.tif 000038
【図３５】
0004772754.tif 000039