(57)【要約】【課題】過度のオーバーシュートやアンダーシュートを生じることなく解像感のある拡大画像を得る。【解決手段】入力画像（Ｄｉｎ）を拡大して第１の拡大画像（Ｄ２Ａ）を出力するとともに、入力画像（Ｄｉｎ）の高周波数成分を拡大して第２の拡大画像（Ｄ２Ｂ）を出力し、該第２の拡大画像（Ｄ２Ｂ）の高周波数成分を取り出して第１の中間画像（Ｄ３２Ａ）を出力するとともに、第２の拡大画像（Ｄ２Ｂ）に対して非線形処理を行った第２の中間画像（Ｄ３２Ｂ）を出力する。さらに、第１の中間画像（Ｄ３２Ａ）及び第２の中間画像（Ｄ３２Ｂ）の少なくとも１つの画素値を、当該画像の画素値の符号及び第１の拡大画像（Ｄ２Ｂ）の画素値応じた増幅率で増幅する（３３Ａ、３３Ｂ）。【選択図】図１
2011048670.tif 000002
【特許請求の範囲】
【請求項１】
入力画像を拡大する画像処理装置において、
前記入力画像を拡大し、第１の拡大画像を出力する第１の画像拡大手段と、
前記入力画像の高周波数成分を取り出し、第１の高周波数成分画像を生成する第１の高周波数成分画像生成手段と、
前記第１の高周波数成分画像を拡大し、第２の拡大画像を出力する第２の画像拡大手段と、
前記第２の拡大画像を入力とし、第２の高周波数成分画像を出力する高周波数成分画像処理手段と、
前記第１の拡大画像と前記第２の高周波数成分画像を加算する第１の加算手段と
を有する画像処理装置において、
前記高周波数成分画像処理手段は、
前記第２の拡大画像の高周波数成分を取り出して第１の中間画像を出力する第２の高周波数成分画像生成手段を含む第１の補正成分生成手段と、
前記第２の拡大画像に対して非線形処理を含む処理を行った第２の中間画像を出力する非線形処理画像生成手段を含を含む第２の補正成分生成手段と、
前記第１の補正成分生成手段の出力と、前記第２の補正成分生成手段の出力を加算する第２の加算手段とを備え、
前記第２の加算手段における加算の結果が、前記高周波数成分画像処理手段の出力として用いられ、
前記第１の補正成分生成手段が、前記第１の中間画像の各画素値を前記第１の拡大画像の画素値と前記第１の中間画像の画素値の符号に応じて求められる第１の増幅率によって増幅した第３の中間画像を出力する第１の高周波数成分画像補正手段をさらに備え、
前記第３の中間画像が前記第１の補正成分生成手段の出力として用いられる
ことを特徴とする画像処理装置。
【請求項２】
入力画像を拡大する画像処理装置において、
前記入力画像を拡大し、第１の拡大画像を出力する第１の画像拡大手段と、
前記入力画像の高周波数成分を取り出し、第１の高周波数成分画像を生成する第１の高周波数成分画像生成手段と、
前記第１の高周波数成分画像を拡大し、第２の拡大画像を出力する第２の画像拡大手段と、
前記第２の拡大画像を入力とし、第２の高周波数成分画像を出力する高周波数成分画像処理手段と、
前記第１の拡大画像と前記第２の高周波数成分画像を加算する第１の加算手段と
を有する画像処理装置において、
前記高周波数成分画像処理手段は、
前記第２の拡大画像の高周波数成分を取り出して第１の中間画像を出力する第２の高周波数成分画像生成手段を含む第１の補正成分生成手段と、
前記第２の拡大画像に対して非線形処理を含む処理を行った第２の中間画像を出力する非線形処理画像生成手段を含を含む第２の補正成分生成手段と、
前記第１の補正成分生成手段の出力と、前記第２の補正成分生成手段の出力を加算する第２の加算手段とを備え、
前記第２の加算手段における加算の結果が、前記高周波数成分画像処理手段の出力として用いられ、
前記第２の補正成分生成手段が、
前記第２の中間画像の各画素値を前記第１の拡大画像の画素値と前記第２の中間画像の画素値の符号に応じて求められる第２の増幅率によって増幅した第４の中間画像を出力する第２の高周波数成分画像補正手段をさらに備え、
前記第４の中間画像が前記第２の補正成分生成手段の出力として用いられる
ことを特徴とする画像処理装置。
【請求項３】
前記第２の補正成分生成手段が、前記第２の中間画像の各画素値を前記第１の拡大画像の画素値と前記第２の中間画像の画素値の符号に応じて求められる第２の増幅率によって増幅した第４の中間画像を出力する第２の高周波数成分画像補正手段をさらに備え、
前記第４の中間画像が前記第２の補正成分生成手段の出力として用いられる
ことを特徴とする請求項１に記載の画像処理装置。
【請求項４】
入力画像を拡大する画像処理装置において、
前記入力画像を拡大し、第１の拡大画像を出力する第１の画像拡大手段と、
前記入力画像の高周波数成分を取り出し、第１の高周波数成分画像を生成する第１の高周波数成分画像生成手段と、
前記第１の高周波数成分画像を拡大し、第２の拡大画像を出力する第２の画像拡大手段と、
前記第２の拡大画像を入力とし、第２の高周波数成分画像を出力する高周波数成分画像処理手段と、
前記第１の拡大画像と前記第２の高周波数成分画像を加算する第１の加算手段と
を有する画像処理装置において、
前記高周波数成分画像処理手段は、
前記第２の拡大画像に対して非線形処理を含む処理を行った第２の中間画像を出力する非線形処理画像生成手段と、
前記第２の中間画像の各画素値を前記第１の拡大画像の画素値と前記第２の中間画像の画素値の符号に応じて求められる第２の増幅率によって増幅した第４の中間画像を出力する第２の高周波数成分画像補正手段とを有する
ことを特徴とする画像処理装置。
【請求項５】
前記非線形処理画像生成手段は、
前記第２の拡大画像に対して非線形処理を行った第１の非線形処理画像を生成する非線形処理手段と、
前記第１の非線形処理画像から高周波数成分を取り出し、前記第２の中間画像とする第３の高周波数成分画像生成手段を含む
ことを特徴とする請求項１乃至４いずれかに記載の画像処理装置。
【請求項６】
前記第１の高周波数成分画像補正手段は、
前記第１の増幅率を決定する第１の増幅率計算手段と、
前記第１の増幅率で前記第１の中間画像の各画素値を増幅する第１の画素値増幅手段を含む
ことを特徴とする請求項１又は３に記載の画像処理装置。
【請求項７】
前記第２の高周波数成分画像補正手段は、
前記第２の増幅率を決定する第１の増幅率計算手段と、
前記第２の増幅率で前記第２の中間画像の各画素値を増幅する第２の画素値増幅手段を含む
ことを特徴とする請求項２、３又は４に記載の画像処理装置。
【請求項８】
前記第１の増幅率計算手段は、前記第１の拡大画像の画素値が大きく、前記第１の中間画像の画素値の符号が正の場合、前記第１の増幅率が小さくなるよう前記第１の増幅率を求めることを特徴とする請求項６に記載の画像処理装置。
【請求項９】
前記第１の増幅率計算手段は、前記第１の拡大画像の画素値が小さく、前記第１の中間画像の画素値の符号が負の場合、前記第１の増幅率が小さくなるよう前記第１の増幅率を求めることを特徴とする請求項６に記載の画像処理装置。
【請求項１０】
前記第２の増幅率計算手段は、前記第１の拡大画像の画素値が大きく、前記第２の中間画像の画素値の符号が正の場合、前記第２の増幅率が小さくなるよう前記第２の増幅率を求めることを特徴とする請求項７に記載の画像処理装置。
【請求項１１】
前記第２の増幅率計算手段は、前記第１の拡大画像の画素値が小さく、前記第２の中間画像の画素値の符号が負の場合前記第２の増幅率が小さくなるよう前記第２の増幅率を求める
ことを特徴とする請求項７に記載の画像処理装置。
【請求項１２】
前記第１の高周波数成分画像補正手段は、
前記第１の増幅率を決定する第１の増幅率計算手段と、
前記第１の増幅率で前記第１の中間画像の各画素値を増幅する第１の画素値増幅手段を含み、
前記第２の高周波数成分画像補正手段は、
前記第２の増幅率を決定する第１の増幅率計算手段と、
前記第２の増幅率で前記第２の中間画像の各画素値を増幅する第２の画素値増幅手段を含む
ことを特徴とする請求項３に記載の画像処理装置。
【請求項１３】
前記第１の高周波数成分画像生成手段は、
前記入力画像の各画素について水平方向近傍に存在する画素の画素値を用いて水平方向高周波数成分画像を生成する第１の水平方向高周波数成分画像生成手段を有し、
前記第２の画像拡大手段は、
前記第１の水平方向高周波数成分画像を拡大し、第３の拡大画像を出力する第３の画像拡大手段を有し、
前記第２の拡大画像は、前記第３の拡大画像を含み、
前記第２の高周波数成分画像生成手段は、
前記第３の拡大画像の高周波数成分のみを取り出し、第１の水平方向中間画像を出力する第２の水平方向高周波数成分画像生成手段を有し、
前記非線形処理手段は、
前記第３の拡大画像に対して非線形処理を行った第２の非線形処理画像を出力する水平方向非線形処理手段を有し、
前記第３の高周波数成分画像生成手段は、
前記第２の非線形処理画像から高周波数成分を取り出し、第２の水平方向中間画像を出力する第３の水平方向高周波数成分画像生成手段を有し、
前記第１の中間画像は前記第１の水平方向中間画像を含み、
前記第２の中間画像は前記第２の水平方向中間画像を含み、
前記第１の増幅率計算手段は、前記第１の拡大画像の画素値と前記第１の水平方向中間画像の画素値の符号に応じて第１の水平方向増幅率を計算する第１の水平方向増幅率計算手段を有し、
前記第１の画素値増幅手段は、前記第１の水平方向増幅率で前記第１の水平方向中間画像の各画素値を増幅した第３の水平方向中間画像を出力する第１の水平方向画素値増幅手段を有し、
前記第２の増幅率計算手段は、第２の水平方向増幅率を計算する第２の水平方向増幅率計算手段を有し、
前記第２の画素値増幅手段は、前記第２の水平方向増幅率で前記第２の水平方向中間画像の各画素値を増幅した第４の水平方向中間画像を出力する第２の水平方向画素値増幅手段を有し、
前記第３の中間画像は前記第３の水平方向中間画像を含み、
前記第４の中間画像は前記第４の水平方向中間画像を含む
ことを特徴とする請求項１２に記載の画像処理装置。
【請求項１４】
前記第１の高周波数成分画像生成手段は、
前記入力画像の各画素について垂直方向近傍に存在する画素の画素値を用いて垂直方向高周波数成分画像を生成する第１の垂直方向高周波数成分画像生成手段を有し、
前記第２の画像拡大手段は、
前記第１の垂直方向高周波数成分画像を拡大し、第４の拡大画像を出力する第４の画像拡大手段を有し、
前記第２の拡大画像は、前記第４の拡大画像を含み、
前記第２の高周波数成分画像生成手段は、
前記第４の拡大画像の高周波数成分のみを取り出し、第１の垂直方向中間画像を出力する第２の垂直方向高周波数成分画像生成手段を有し、
前記非線形処理手段は、
前記第４の拡大画像に対して非線形処理を行った第３の非線形処理画像を出力する垂直方向非線形処理手段を有し、
前記第３の高周波数成分画像生成手段は、
前記第３の非線形処理画像から高周波数成分を取り出し、第２の垂直方向中間画像を出力する第３の垂直方向高周波数成分画像生成手段を有し、
前記第１の中間画像は前記第１の垂直方向中間画像を含み、
前記第２の中間画像は前記第２の垂直方向中間画像を含み、
前記第１の増幅率計算手段は、前記第１の拡大画像の画素値と前記第１の垂直方向中間画像の画素値の符号に応じて第１の垂直方向増幅率を計算する第１の垂直方向増幅率計算手段を有し、
前記第１の画素値増幅手段は、前記第１の垂直方向増幅率で前記第１の垂直方向中間画像の各画素値を増幅した第３の垂直方向中間画像を出力する第１の垂直方向画素値増幅手段を有し、
前記第２の増幅率計算手段は、第２の垂直方向増幅率を計算する第２の垂直方向増幅率計算手段を有し、
前記第２の画素値増幅手段は、前記第２の垂直方向増幅率で前記第２の垂直方向中間画像の各画素値を増幅した第４の垂直方向中間画像を出力する第２の垂直方向画素値増幅手段を有し、
前記第３の中間画像は前記第３の垂直方向中間画像を含み、
前記第４の中間画像は前記第４の垂直方向中間画像を含む
ことを特徴とする請求項１２に記載の画像処理装置。
【請求項１５】
前記第１の高周波数成分画像生成手段は、
前記入力画像の各画素について水平方向近傍に存在する画素の画素値を用いて水平方向高周波数成分画像を生成する第１の水平方向高周波数成分画像生成手段と、
前記入力画像の各画素について垂直方向近傍に存在する画素の画素値を用いて垂直方向高周波数成分画像を生成する第１の垂直方向高周波数成分画像生成手段を有し、
前記第２の画像拡大手段は、
前記第１の水平方向高周波数成分画像を拡大し、第３の拡大画像を出力する第３の画像拡大手段と、
前記第１の垂直方向高周波数成分画像を拡大し、第４の拡大画像を出力する第４の画像拡大手段を有し、
前記第２の拡大画像は、前記第３の拡大画像と前記第４の拡大画像を含み、
前記第２の高周波数成分画像生成手段は、
前記第３の拡大画像の高周波数成分のみを取り出し、第１の水平方向中間画像を出力する第２の水平方向高周波数成分画像生成手段と、
前記第４の拡大画像の高周波数成分のみを取り出し、第１の垂直方向中間画像を出力する第２の垂直方向高周波数成分画像生成手段を有し、
前記非線形処理手段は、
前記第３の拡大画像に対して非線形処理を行った第２の非線形処理画像を出力する水平方向非線形処理手段と、
前記第４の拡大画像に対して非線形処理を行った第３の非線形処理画像を出力する垂直方向非線形処理手段を有し、
前記第３の高周波数成分画像生成手段は、
前記第２の非線形処理画像から高周波数成分を取り出し、第２の水平方向中間画像を出力する第３の水平方向高周波数成分画像生成手段と、
前記第３の非線形処理画像から高周波数成分を取り出し、第２の垂直方向中間画像を出力する第３の垂直方向高周波数成分画像生成手段を有し、
前記第１の中間画像は前記第１の水平方向中間画像と前記第１の垂直方向中間画像を含み、
前記第２の中間画像は前記第２の水平方向中間画像と前記第２の垂直方向中間画像を含み、
前記第１の増幅率計算手段は、
前記第１の拡大画像の画素値と前記第１の水平方向中間画像の画素値の符号に応じて第１の水平方向増幅率を計算する第１の水平方向増幅率計算手段と、
前記第１の拡大画像の画素値と前記第１の垂直方向中間画像の画素値の符号に応じて第１の垂直方向増幅率を計算する第１の垂直方向増幅率計算手段を有し、
前記第１の画素値増幅手段は、
前記第１の水平方向増幅率で前記第１の水平方向中間画像の各画素値を増幅した第３の水平方向中間画像を出力する第１の水平方向画素値増幅手段と、
前記第１の垂直方向増幅率で前記第１の垂直方向中間画像の各画素値を増幅した第３の垂直方向中間画像を出力する第１の垂直方向画素値増幅手段を有し、
前記第２の増幅率計算手段は、
第２の水平方向増幅率を計算する第２の水平方向増幅率計算手段と、
第２の垂直方向増幅率を計算する第２の垂直方向増幅率計算手段を有し、
前記第２の画素値増幅手段は、
前記第２の水平方向増幅率で前記第２の水平方向中間画像の各画素値を増幅した第４の水平方向中間画像を出力する第２の水平方向画素値増幅手段と、
前記第２の垂直方向増幅率で前記第２の垂直方向中間画像の各画素値を増幅した第４の垂直方向中間画像を出力する第２の垂直方向画素値増幅手段を有し、
前記第３の中間画像は前記第３の水平方向中間画像と前記第３の垂直方向中間画像を含み、
前記第４の中間画像は前記第４の水平方向中間画像と前記第４の垂直方向中間画像を含む
ことを特徴とする請求項１２に記載の画像処理装置。
【請求項１６】
前記水平方向非線形形処理手段は、
前記第３の拡大画像の画素値が正から負もしくは負から正へと変化する箇所をゼロクロス点として判定する水平方向ゼロクロス判定手段と、
前記水平方向ゼロクロス判定手段の判定結果に応じて決められる増幅率で前記第３の拡大画像の画素値を増幅する水平方向信号増幅手段とを有し、
前記第３の水平方向高周波数成分画像生成手段は、
前記水平方向信号増幅手段の出力する画像から高周波数成分を取り出し、
前記垂直方向非線形形処理手段は、
前記第４の拡大画像の画素値が正から負もしくは負から正へと変化する箇所をゼロクロス点として判定する垂直方向ゼロクロス判定手段と、
前記垂直方向ゼロクロス判定手段の判定結果に応じて決められる増幅率で前記第４の拡大画像の画素値を増幅する垂直方向信号増幅手段とを有し、
前記第３の垂直方向高周波数成分画像生成手段は、
前記垂直方向信号増幅手段の出力する画像から高周波数成分を取り出す
ことを特徴とする請求項１５に記載の画像処理装置。
【請求項１７】
前記水平方向信号増幅手段は、
前記水平方向ゼロクロス判定手段で判定されたゼロクロス点を含む第１の領域内に存在する画素の画素値に対する増幅率を１より大きな値とし、それ以外の画素の画素値に対する増幅率を１とし、
前記垂直方向信号増幅手段は、
前記垂直方向ゼロクロス判定手段で判定されたゼロクロス点を含む第２の領域内に存在する画素の画素値に対する増幅率を１より大きな値とし、それ以外の画素の画素値に対する増幅率を１とする
ことを特徴とする請求項１６に記載の画像処理装置。
【請求項１８】
前記第１の領域及び前記第２の領域の少なくとも一方は前記第２の画像拡大手段における拡大率に応じて定められる
ことを特徴とする請求項１７に記載の画像処理装置。
【請求項１９】
前記第１の領域に存在する画素に対する増幅率及び前記第２の領域に存在する画素に対する増幅率の少なくとも一方は画素に応じて定められる
ことを特徴とする請求項１７又は１８に記載の画像処理装置。
【請求項２０】
請求項１乃至１９のいずれかに記載の画像処理装置を備えたことを特徴とする画像表示装置。
【請求項２１】
入力画像を拡大する画像処理方法において、
前記入力画像を拡大し、第１の拡大画像を出力する第１の画像拡大ステップと、
前記入力画像の高周波数成分を取り出し、第１の高周波数成分画像を生成する第１の高周波数成分画像生成ステップと、
前記第１の高周波数成分画像を拡大し、第２の拡大画像を出力する第２の画像拡大ステップと、
前記第２の拡大画像の高周波数成分画像に処理を加える高周波数成分画像処理ステップと、
前記第１の拡大画像と前記高周波数成分画像処理ステップによる処理の結果得られる画像を加算する第１の加算ステップと
を有する画像処理方法において、
前記高周波数成分画像処理ステップは、
前記第２の拡大画像の高周波数成分を取り出す第２の高周波数成分画像生成ステップと、
前記第２の拡大画像に対して非線形処理を含む処理を行う非線形処理画像生成ステップと、
前記第２の高周波数成分画像生成ステップにより得られる画像及び前記非線形処理画像生成ステップにより得られる画像を加算する第２の加算ステップとを含み、
前記第２の高周波数成分画像生成ステップにより得られる画像、及び前記非線形処理画像生成ステップにより得られる画像の少なくとも一つに対し、前記第２の加算ステップによる加算の前に当該画像の画素値の符号と前記第１の拡大画像の画素値に基づいて高周波数成分画像の補正を行う高周波数成分画像補正ステップをさらに有し、
前記高周波数成分画像補正ステップは、
前記補正の対象となる画像の各画素の画素値を、当該画像の画素値の符号と前記第１の拡大画像の画素値に応じて求められた増幅率によって増幅することを特徴とする
画像処理方法。
【請求項２２】
入力画像を拡大する画像処理方法において、
前記入力画像を拡大し、第１の拡大画像を出力する第１の画像拡大ステップと、
前記入力画像の高周波数成分を取り出し、第１の高周波数成分画像を生成する第１の高周波数成分画像生成ステップと、
前記第１の高周波数成分画像を拡大し、第２の拡大画像を出力する第２の画像拡大ステップと、
前記第２の拡大画像の高周波数成分画像に処理を加える高周波数成分画像処理ステップと、
前記第１の拡大画像と前記高周波数成分画像処理ステップによる処理の結果得られる画像を加算する第１の加算ステップと
を有する画像処理方法において、
前記高周波数成分画像処理ステップは、
前記第２の拡大画像に対して非線形処理を含む処理を行った第２の中間画像を出力する非線形処理画像生成ステップと、
前記第２の中間画像の各画素値を前記第１の拡大画像の画素値と前記第１の中間画像の画素値の符号に応じて求められる第２の増幅率によって増幅した第４の中間画像を出力する第２の高周波数成分画像補正ステップとを有し、
前記高周波数成分画像補正ステップは、
前記補正の対象となる画像の各画素の画素値を、当該画像の画素値の符号と前記第１の拡大画像の画素値に応じて求められた増幅率によって増幅することを特徴とする
ことを特徴とする画像処理方法。
【請求項２３】
請求項２１又は２２に記載の画像処理方法により処理された画像を表示することを特徴とする画像表示装置。
【請求項２４】
入力画像を拡大する画像処理装置において、
前記入力画像を拡大し、第１の拡大画像を出力する第１の画像拡大手段と、
前記入力画像の高周波数成分を取り出し、第１の高周波数成分画像を生成する第１の高周波数成分画像生成手段と、
前記第１の高周波数成分画像を拡大し、第２の拡大画像を出力する第２の画像拡大手段と、
前記第２の拡大画像を入力とし、第２の高周波数成分画像を出力する高周波数成分画像処理手段と、
前記第２の高周波数成分画像の各画素値を前記第１の拡大画像の画素値と前記第２の高周波数成分画像の各画素値の符号に応じて求められる第３の増幅率によって増幅した第３の高周波数成分画像を出力する第３の高周波数成分画像補正手段と、
前記第１の拡大画像と前記第３の高周波数成分画像を加算する第１の加算手段と
を有する画像処理装置において、
前記高周波数成分画像処理手段は、
前記第２の拡大画像の高周波数成分を取り出して第１の中間画像を出力する第２の高周波数成分画像生成手段と、
前記第２の拡大画像に対して非線形処理を含む処理を行った第２の中間画像を出力する非線形処理画像生成手段と、
前記第１の中間画像と前記第２の中間画像を加算する第２の加算手段とを含み、
前記第２の加算手段による加算の結果が前記第２の高周波数成分画像として前記高周波数成分画像処理手段から出力される
ことを特徴とする画像処理装置。
【請求項２５】
入力画像を拡大する画像処理装置において、
前記入力画像を拡大し、第１の拡大画像を出力する第１の画像拡大手段と、
前記入力画像の高周波数成分を取り出し、第１の高周波数成分画像を生成する第１の高周波数成分画像生成手段と、
前記第１の高周波数成分画像を拡大し、第２の拡大画像を出力する第２の画像拡大手段と、
前記第２の拡大画像に対して非線形処理を含む処理を行う非線形処理画像生成手段と、
前記非線形処理画像生成手段によって非線形処理された画像を補正対象画像とし、前記補正対象画像の画素値を、前記補正対象画像の画素値の符号と前記第１の拡大画像の画素値に応じて求められる増幅率によって増幅する高周波数成分画像補正手段と
前記高周波数成分画像補正手段による処理を受けた画像と前記第１の拡大画像を加算する加算手段を含む
ことを特徴とする画像処理装置。
【請求項２６】
入力画像を拡大する画像処理方法において、
前記入力画像を拡大し、第１の拡大画像を出力する第１の画像拡大ステップと、
前記入力画像の高周波数成分を取り出し、第１の高周波数成分画像を生成する第１の高周波数成分画像生成ステップと、
前記第１の高周波数成分画像を拡大し、第２の拡大画像を出力する第２の画像拡大ステップと、
前記第２の拡大画像を入力とし、第２の高周波数成分画像を出力する高周波数成分画像処理ステップと、
前記第２の高周波数成分画像の各画素値を前記第１の拡大画像の画素値と前記第２の高周波数成分画像の各画素値の符号に応じて求められる第３の増幅率によって増幅した第３の高周波数成分画像を出力する第３の高周波数成分画像補正ステップと、
前記第１の拡大画像と前記第３の高周波数成分画像を加算する第１の加算ステップと
を有する画像処理方法において、
前記高周波数成分画像処理ステップは、
前記第２の拡大画像の高周波数成分を取り出して第１の中間画像を出力する第２の高周波数成分画像生成ステップと、
前記第２の拡大画像に対して非線形処理を含む処理を行った第２の中間画像を出力する非線形処理画像生成ステップと、
前記第１の中間画像と前記第２の中間画像を加算する第２の加算ステップとを含み、
前記第２の加算ステップによる加算の結果を前記第２の高周波数成分画像とする
ことを特徴とする画像処理方法。
【請求項２７】
入力画像を拡大する画像処理方法において、
前記入力画像を拡大し、第１の拡大画像を出力する第１の画像拡大ステップと、
前記入力画像の高周波数成分を取り出し、第１の高周波数成分画像を生成する第１の高周波数成分画像生成ステップと、
前記第１の高周波数成分画像を拡大し、第２の拡大画像を出力する第２の画像拡大ステップと、
前記第２の拡大画像に対して非線形処理を含む処理を行う非線形処理画像生成ステップと、
前記非線形処理画像生成手段によって非線形処理された画像を補正対象画像とし、前記補正対象画像の画素値を、前記補正対象画像の画素値の符号と前記第１の拡大画像の画素値に応じて求められる増幅率によって増幅する高周波数成分画像補正ステップと
前記高周波数成分画像補正ステップによる処理を受けた画像と前記第１の拡大画像を加算する加算ステップを含む
ことを特徴とする画像処理方法。
【請求項２８】
請求項２４又は２５に記載の画像処理装置又は請求項２６又は２７に記載の画像処理方法により処理された画像を表示することを特徴とする画像表示装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、デジタル化された画像を拡大する画像処理装置及び方法、並びに画像表示装置及び方法に関するものであり、画像を拡大する際に、高周波数成分を生成することによって、解像感の高い拡大画像を得るものである。
【背景技術】
【０００２】
一般に画像処理装置は、出力画像の画素数が入力画像の画素数より多い場合、画像を拡大処理しなければならない。従来の画像処理装置では注目する画素近傍の画素がもつ画素値を重み付け加算して画像を拡大していた。
【０００３】
例えば特許文献１に記載された画像処理装置においては、各シフトレジスタからそれぞれ出力される主走査方向の隣り合った５個の画素データに所定の重み付け定数を乗算し、各画素データにおける乗算結果を加算するための演算回路を備え、画像データの拡大処理を行う場合に、演算回路での演算結果をこれらの画素データの中央の画素データとしてセレクタにて選択して出力している。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開平６−３１１３４６号公報（図１）
【発明の概要】
【発明が解決しようとする課題】
【０００５】
注目する画素近傍の画素がもつ画素値を重み付け加算することは入力画像の低周波数成分のみを通過させるローパスフィルタ処理となる。従って、上記の従来の技術では、拡大画像に対して高周波数成分を十分に与えることが出来ないため、拡大画像の解像感が失われるという問題があった。
【課題を解決するための手段】
【０００６】
本発明は上述のような課題を解消するためになされたもので、本発明の画像処理装置は、
入力画像を拡大する画像処理装置において、
前記入力画像を拡大し、第１の拡大画像を出力する第１の画像拡大手段と、
前記入力画像の高周波数成分を取り出し、第１の高周波数成分画像を生成する第１の高周波数成分画像生成手段と、
前記第１の高周波数成分画像を拡大し、第２の拡大画像を出力する第２の画像拡大手段と、
前記第２の拡大画像を入力とし、第２の高周波数成分画像を出力する高周波数成分画像処理手段と、
前記第１の拡大画像と前記第２の高周波数成分画像を加算する第１の加算手段と
を有する画像処理装置において、
前記高周波数成分画像処理手段は、
前記第２の拡大画像の高周波数成分を取り出して第１の中間画像を出力する第２の高周波数成分画像生成手段を含む第１の補正成分生成手段と、
前記第２の拡大画像に対して非線形処理を含む処理を行った第２の中間画像を出力する非線形処理画像生成手段を含を含む第２の補正成分生成手段と、
前記第１の補正成分生成手段の出力と、前記第２の補正成分生成手段の出力を加算する第２の加算手段とを備え、
前記第２の加算手段における加算の結果が、前記高周波数成分画像処理手段の出力として用いられ、
前記第１の補正成分生成手段が、前記第１の中間画像の各画素値を前記第１の拡大画像の画素値と前記第１の中間画像の画素値の符号に応じて求められる第１の増幅率によって増幅した第３の中間画像を出力する第１の高周波数成分画像補正手段をさらに備え、
前記第３の中間画像が前記第１の補正成分生成手段の出力として用いられる
ことを特徴とする。
【発明の効果】
【０００７】
本発明によれば、過度のオーバーシュートやアンダーシュートを発生することなく、高周波数成分を十分に与えることができ、解像感のある拡大画像を得ることができる。
【図面の簡単な説明】
【０００８】
【図１】本発明の実施の形態１の画像処理装置の構成を示すブロック図である。
【図２】実施の形態１の画像処理装置を用いた画像表示装置の構成例を示すブロック図である。
【図３】図１の画像処理装置の構成をより詳細に示すブロック図である。
【図４】（ａ）〜（ｄ）は、画像拡大手段２Ａの動作を示す画素配置図である。
【図５】画像拡大手段２Ａの構成例を示すブロック図である。
【図６】画像拡大手段２Ｂの構成例を示すブロック図である。
【図７】水平方向非線形処理手段３１ｈの構成例を示すブロック図である。
【図８】垂直方向非線形処理手段３１ｖの構成例を示すブロック図である。
【図９】高周波数成分処理手段３３Ａの構成例を示すブロック図である。
【図１０】（ａ）〜（ｄ）は増幅率計算手段３ＭＡの画素値と増幅率の関係を示す図である。
【図１１】高周波数成分処理手段３３Ｂの構成例を示すブロック図である。
【図１２】（ａ）〜（ｄ）は増幅率計算手段３ＨＡの画素値と増幅率の関係を示す図である。
【図１３】（ａ）〜（ｄ）は、拡大画像Ｄ２Ａを得る過程を説明するための周波数スペクトル及び周波数応答を示す図である。
【図１４】（ａ）〜（ｆ）は、中間画像Ｄ３２Ａを得る過程を説明するための周波数スペクトル及び周波数応答を示す図である。
【図１５】（ａ）〜（ｃ）は、中間画像Ｄ３２Ｂを得る過程を説明するための周波数スペクトル及び周波数応答を示す図である。
【図１６】（ａ）〜（ｅ）は、ステップエッジ信号とステップエッジ信号を異なるサンプリング周波数でサンプリングした際に得られる信号とその高周波数成分の信号強度を示す図である。
【図１７】（ａ）〜（ｆ）は、非線形処理手段３１と高周波数成分画像生成手段３２Ｂの動作を説明するための、信号強度を示す図である。
【図１８】（ａ）〜（ｄ）は、中間画像Ｄ３２Ａ、Ｄ３２Ｂを加算する効果の説明するための、信号強度を示す図である。
【図１９】（ａ）及び（ｂ）は、中間画像Ｄ３２Ａ、Ｄ３２Ｂを加算する効果及び不具合の説明するための信号強度を示す図である。
【図２０】増幅率計算手段３ＭＡ、３ＨＡの画素値と増幅率の関係を示す図である。
【図２１】拡大画像Ｄｏｕｔの周波数スペクトルの説明図である。
【図２２】（ａ）〜（ｅ）は、ステップエッジ信号とステップエッジ信号を異なるサンプリング周波数でサンプリングした際に得られる信号とその高周波数成分の信号強度を示す図である。
【図２３】（ａ）〜（ｆ）は、非線形処理手段３１と高周波数成分画像生成手段３２Ｂの動作を説明するための信号の強度を示す図である。
【図２４】本発明の実施の形態２の画像処理装置の構成例を示すブロック図である。
【図２５】実施の形態２の画像処理方法を示すフロー図である。
【図２６】高周波数成分画像生成ステップＳＴ１を示すフロー図である。
【図２７】画像拡大ステップＳＴ２Ｂを示すフロー図である。
【図２８】高周波数成分画像処理ステップＳＴ３を示すフロー図である。
【図２９】水平方向非線形処理ステップＳＴ３１ｈを示すフロー図である。
【図３０】垂直方向非線形処理ステップＳＴ３１ｖを示すフロー図である。
【図３１】高周波数成分画像補正ステップＳＴ３３Ａを示すフロー図である。
【図３２】高周波数成分画像補正ステップＳＴ３３Ｂを示すフロー図である。
【図３３】実施の形態１の画像処理装置の変形例を示すブロック図である。
【図３４】実施の形態２の高周波数成分画像処理ステップＳＴ３の変形例を示すフロー図である。
【図３５】本発明の実施の形態３の画像処理装置の構成を示すブロック図である。
【図３６】実施の形態３の画像処理装置の変形例を示すブロック図である。
【図３７】本発明の実施の形態４の画像処理方法を示すフロー図である。
【図３８】実施の形態４の画像処理方法の変形例を示すフロー図である。
【発明を実施するための形態】
【０００９】
実施の形態１．
図１は本発明の実施の形態１による画像処理装置の構成を表す図であり、例えば図２に示す画像表示装置の一部として用いることができる。ここで図２に示す画像表示装置は図１に示す画像処理装置を内部に含む画像処理装置Ｕ１及び表示部９を備えており、画像処理装置Ｕ１において画像ＤＯＲＧに対する出力として得られた画像ＤＵ１が表示部９に表示される。
【００１０】
実施の形態１による画像処理装置は、画像拡大手段２Ａと、高周波数成分画像生成手段１と、画像拡大手段２Ｂと、高周波数成分画像処理手段３と、加算手段４とを備える。
画像拡大手段２Ａは、入力画像Ｄｉｎを拡大して拡大画像Ｄ２Ａを生成する。
高周波数成分画像生成手段１は、入力画像Ｄｉｎの高周波数成分のみを取り出して高周波数成分画像Ｄ１を生成する。
【００１１】
画像拡大手段２Ｂは、高周波数成分画像生成手段１から出力される高周波数成分画像Ｄ１を拡大して拡大画像（高周波数成分拡大画像）Ｄ２Ｂを生成する。
高周波数成分画像処理手段３は、画像拡大手段２Ｂから出力される拡大画像Ｄ２Ｂに対して後述の処理を行い、高周波数成分画像（高周波数成分処理画像）Ｄ３を生成する。
【００１２】
加算手段４は、画像拡大手段２Ａから出力される拡大画像Ｄ２Ａに、高周波数成分画像処理手段３から出力される高周波数成分画像Ｄ３を加算して、その結果を最終的な拡大画像、即ち出力画像Ｄｏｕｔとして出力する。加算手段４の出力は、例えば図２に示す画像表示装置の表示部９に画像ＤＵ１として供給され、表示部９による画像表示に用いられる。
【００１３】
なお、本明細書において、拡大、高周波数成分生成、高周波数成分処理などの処理は「画像」に対して行なわれる旨記載されるが、具体的には、画像を表すデジタルデータに対して行われる。また、「画像」との記載も具体的には「画像データ」を意味する場合がある。
【００１４】
画像拡大手段２Ａ、高周波数成分画像生成手段１、画像拡大手段２Ｂ、及び高周波数成分画像処理手段３の詳細な動作については後述するが、高周波数成分画像Ｄ３のもつ周波数成分は拡大画像Ｄ２Ａがもつ周波数成分より高い周波数帯域のものとなる。従って加算手段４において拡大画像Ｄ２Ａに高周波数成分画像Ｄ３を加算することで、高周波数成分を多く含んだ拡大画像Ｄｏｕｔを得ることができる。
【００１５】
図３は図１に示した本発明の実施の形態１による画像処理装置の構成の詳細を表す図である。以下、図３を参照して、画像拡大手段２Ａ、高周波数成分画像生成手段１、画像拡大手段２Ｂ、及び高周波数成分画像処理手段３の構成をより詳細に説明する。
【００１６】
画像拡大手段２Ａは、水平方向及び垂直方向の少なくとも一方に画像を拡大するものであり、例えば、水平方向及び垂直方向に同じ倍率で拡大を行なうが、代わりに、水平方向及び垂直方向に異なる倍率で拡大を行なうものであっても良い。また、水平方向及び垂直方向の一方にのみ拡大を行なうものであっても良く、例えば入力画像に対して表示画面が横長である場合に水平方向にのみ拡大を行なうことがある。
【００１７】
高周波数成分画像生成手段１は、入力画像Ｄｉｎの高周波数成分（所定の周波数Ｆｂよりも高い成分）を取り出して、高周波数成分画像Ｄ１を生成するものであり、後述の方法でそれぞれ水平方向高周波数成分画像Ｄ１ｈ及び垂直方向高周波数成分画像Ｄ１ｖを生成する水平方向高周波数成分画像生成手段１ｈ及び垂直方向高周波数成分画像生成手段１ｖを備える。水平方向高周波数成分画像Ｄ１ｈと垂直方向高周波数成分画像Ｄ１ｖとで高周波数成分画像Ｄ１が構成されている。
【００１８】
画像拡大手段２Ｂは、水平方向高周波数成分画像Ｄ１ｈを拡大した拡大画像Ｄ２Ｂｈを生成する画像拡大手段２Ｂｈ及び垂直方向高周波数成分画像Ｄ１ｖを拡大した拡大画像Ｄ２Ｂｖを生成する画像拡大手段２Ｂｖを備える。拡大画像Ｄ２Ｂｈと拡大画像Ｄ２Ｂｖとで拡大画像Ｄ２Ｂが構成されている。
【００１９】
画像拡大手段２Ａが水平方向及び垂直方向の両方向に拡大を行なう場合、画像拡大手段２Ｂｈは、水平方向高周波数成分画像Ｄ１ｈを水平方向及び垂直方向の両方向に拡大し、画像拡大手段２Ｂｖは、垂直方向高周波数成分画像Ｄ１ｖを水平方向及び垂直方向の両方向に拡大する。画像拡大手段２Ｂｈ及び２Ｂｖによる、水平方向高周波数成分画像Ｄ１ｈ及び垂直方向高周波数成分画像Ｄ１ｖの拡大は、画像拡大手段２Ａによる拡大と水平方向及び垂直方向の各々について同じ倍率で行なわれる。
【００２０】
高周波数成分画像処理手段３は、第１の補正成分生成手段３Ａと、第２の補正成分生成手段３Ｂと、加算手段３４を備えている。第１の補正成分生成手段３Ａは、高周波数成分画像生成手段３２Ａと、高周波数成分画像補正手段３３Ａとを備え、第２の補正成分生成手段３Ｂは、非線形処理画像生成手段３０と、高周波数成分画像補正手段３３Ｂとを備える。
【００２１】
高周波数成分画像生成手段３２Ａは、拡大画像Ｄ２Ｂの高周波数成分（所定の周波数Ｆｄよりも高い成分）を取り出して中間画像（高周波数成分画像）Ｄ３２Ａを出力するものであり、拡大画像Ｄ２Ｂｈに含まれる水平方向の高周波数成分のみを取り出した水平方向中間画像Ｄ３２Ａｈを生成する水平方向高周波数成分画像生成手段３２Ａｈと、拡大画像Ｄ２Ｂｖの垂直方向の高周波数成分のみを取り出した垂直方向中間画像Ｄ３２Ａｖを生成する垂直方向高周波数成分画像生成手段３２Ａｖを備え、高周波数成分画像生成手段３２Ａからは、水平方向中間画像Ｄ３２Ａｈと垂直方向中間画像Ｄ３２Ａｖから成る中間画像Ｄ３２Ａが出力される。
【００２２】
非線形処理画像生成手段（エッジ鮮鋭化画像生成手段）３０は、拡大画像Ｄ２Ｂに対して非線形処理を含む処理を行った中間画像（エッジ鮮鋭化画像）Ｄ３２Ｂを出力するものであり、非線形処理手段３１、及び高周波数成分画像生成手段３２Ｂを備える。
非線形処理手段３１は、拡大画像Ｄ２Ｂに対して後述するエッジの鮮鋭化のための非線形処理を行った非線形処理画像Ｄ３１を生成する。
高周波数成分画像生成手段３２Ｂは非線形処理画像Ｄ３１に含まれる高周波数成分（所定の周波数Ｆｆよりも高い成分）のみを取り出した中間画像Ｄ３２Ｂを出力する。
【００２３】
非線形処理手段３１は、拡大画像Ｄ２Ｂｈに対して非線形処理した非線形処理画像Ｄ３１ｈを生成する水平方向非線形処理手段３１ｈと、拡大画像Ｄ２Ｂｖに対して非線形処理した非線形処理画像Ｄ３１ｖを生成する垂直方向非線形処理手段３１ｖを備えており、非線形処理画像Ｄ３１は非線形処理画像Ｄ３１ｈと非線形処理画像Ｄ３１ｖから成る。
【００２４】
高周波数成分画像生成手段３２Ｂは、非線形処理画像Ｄ３１ｈから高周波数成分を取り出し、水平方向中間画像Ｄ３２Ｂｈを生成する水平方向高周波数成分画像生成手段３２Ｂｈと、非線形処理画像Ｄ３１ｖから高周波数成分を取り出し、垂直方向中間画像Ｄ３２Ｂｖを生成する垂直方向高周波数成分画像生成手段３２Ｂｖを備え、中間画像Ｄ３２Ｂは水平方向中間画像Ｄ３２Ｂｈと垂直方向中間画像Ｄ３２Ｂｖから成る。
【００２５】
高周波数成分画像補正手段３３Ａは、拡大画像Ｄ２Ａと中間画像Ｄ３２Ａとを入力とし、中間画像Ｄ３２Ａを拡大画像Ｄ２Ａの画素値及び中間画像Ｄ３２Ａの画素値の符号に応じて求められる増幅率によって増幅した中間画像（補正画像）Ｄ３３Ａを出力するものであり、水平方向中間画像Ｄ３２Ａｈに対して補正処理を行って、水平方向中間画像Ｄ３３Ａｈを生成する水平方向高周波数成分画像補正手段３３Ａｈと、垂直方向中間画像Ｄ３２Ａｖに対して補正処理を行って、垂直方向中間画像Ｄ３３Ａｖを生成する垂直方向高周波数成分画像補正手段３３Ａｖを備える。水平方向中間画像Ｄ３３Ａｈと垂直方向中間画像Ｄ３３Ａｖとで中間画像Ｄ３３Ａが構成される。
【００２６】
高周波数成分画像補正手段３３Ｂは、拡大画像Ｄ２Ａと中間画像Ｄ３２Ｂとを入力とし、中間画像Ｄ３２Ｂを拡大画像Ｄ２Ａの画素値及び中間画像Ｄ３２Ｂの画素値の符号に応じて求められる増幅率によって増幅した中間画像（補正画像）Ｄ３３Ｂを出力するものであり、水平方向中間画像Ｄ３２Ｂｈに対して補正処理を行って、水平方向中間画像Ｄ３３Ｂｈを生成する水平方向高周波数成分画像補正手段３３Ｂｈと、垂直方向中間画像Ｄ３２Ｂｖに対して補正処理を行って、垂直方向中間画像Ｄ３３Ｂｖを生成する垂直方向高周波数成分画像補正手段３３Ｂｖを備える。水平方向中間画像Ｄ３３Ｂｈと垂直方向中間画像Ｄ３３Ｂｖとで中間画像Ｄ３３Ｂが構成される。
【００２７】
加算手段３４は、中間画像Ｄ３３Ａと中間画像Ｄ３３Ｂを加算し、高周波数成分画像Ｄ３として出力する。
【００２８】
以下、入力画像Ｄｉｎを水平方向、垂直方向ともに２倍に拡大した拡大画像Ｄｏｕｔを生成する場合を例にし、各構成要素の動作をさらに詳細に説明する。この説明を通じて本発明の作用、効果もより明らかなものとなるであろう。
【００２９】
まず、画像拡大手段２Ａの動作について説明する。画像拡大手段２Ａは入力画像Ｄｉｎを水平方向、垂直方向ともに２倍に拡大した拡大画像Ｄ２Ａを生成する。図４（ａ）〜（ｄ）は画像拡大手段２Ａにおける拡大画像Ｄ２Ａの生成手順の一例を模式的に示した図であり、図５は画像拡大手段２Ａの一例を示した図である。
【００３０】
画像拡大手段２Ａは、ゼロ挿入手段２１Ａと、低周波数成分数通過手段２２Ａを備える。以下、図４（ａ）〜（ｄ）を用いてゼロ挿入手段２１Ａ及び低周波数成分通過手段２２Ａの動作を説明する。図４（ａ）は入力画像Ｄｉｎ（特に画像の一部を構成する画素の配列）を、図４（ｂ）はゼロ挿入手段２１Ａで生成されるゼロ挿入画像Ｄ２１Ａを、図４（ｃ）は低周波数成分通過手段２２Ａにおいて拡大画像Ｄ２Ａを生成する際使用されるフィルタ係数を、図４（ｄ）は低周波数成分通過手段２２Ａで生成された拡大画像Ｄ２Ａを表す。図４（ｂ）、（ｄ）には画素の位置に対応させて水平座標Ｘ、垂直座標Ｙを記載している。
【００３１】
ゼロ挿入手段２１Ａでは入力画像Ｄｉｎに対して画素値０をもつ画素を水平方向には（入力画像Ｄｉｎの）１画素につき１個（隣り合う２つの画素相互間に１個）、垂直方向には（入力画像Ｄｉｎの）１ラインにつき１本（隣り合う２本のライン相互間に１本）を挿入したゼロ挿入画像Ｄ２１Ａを生成する。
「ＰＸＹ」が入力画像Ｄｉｎの座標（Ｘ，Ｙ）における画素の画素値を表し、「Ｐ’ＸＹ」がゼロ挿入画像Ｄ２１Ａの座標（Ｘ，Ｙ）の画素の画素値を表すとすると、ゼロ挿入画像Ｄ２１Ａの、Ｐ’（２Ｘ−１）（２Ｙ−１）で表わされる画素値は、入力画像ＤｉｎのＰＸＹに等しく、Ｐ’（２Ｘ−１）（２Ｙ）、Ｐ’（２Ｘ）（２Ｙ）、Ｐ’（２Ｘ）（２Ｙ−１）で表わされる画素値は、ゼロに等しい。
【００３２】
低周波数成分通過手段２２Ａではゼロ挿入画像Ｄ２１Ａに対し、図４（ｃ）に示されたフィルタ係数で表されたフィルタ演算を行うことで、図４（ｄ）に示される拡大画像Ｄ２Ａを生成する。
例えば、拡大画像Ｄ２Ａに含まれる、座標（Ｘ，Ｙ）の画素の画素値ＱＸＹは下記の式（１）のように計算される。
【００３３】
ＱＸＹ＝（４／１６）×
｛Ｐ’（Ｘ−１）（Ｙ−１）＋２Ｐ’Ｘ（Ｙ−１）＋Ｐ’（Ｘ＋１）（Ｙ−１）
＋２Ｐ’（Ｘ−１）Ｙ＋４Ｐ’ＸＹ＋２Ｐ’（Ｘ＋１）Ｙ
＋Ｐ’（Ｘ−１）（Ｙ＋１）＋２Ｐ’Ｘ（Ｙ＋１）＋Ｐ’（Ｘ＋１）（Ｙ＋１）｝
…（１）
【００３４】
なお、図４（ｃ）で表されるフィルタ係数はローパスフィルタを表すので、式（１）で表される低周波数成分通過手段２２Ａにおける処理はゼロ挿入画像Ｄ２１Ａの低周波数成分を取り出すことに対応する。
【００３５】
また、式（１）において
Ｐ’（Ｘ−１）（Ｙ−１）、２Ｐ’Ｘ（Ｙ−１）、Ｐ’（Ｘ＋１）（Ｙ−１）、２Ｐ’（Ｘ−１）Ｙ、Ｐ’ＸＹ、２Ｐ’（Ｘ＋１）Ｙ、Ｐ’（Ｘ−１）（Ｙ＋１）、Ｐ’Ｘ（Ｙ＋１）、Ｐ’（Ｘ＋１）（Ｙ＋１）
のうちいくつかはその値が０であり、それ以外は入力画像Ｄｉｎの画素値そのものになる。従って拡大処理は入力画像Ｄｉｎにおいて注目する画素の近傍の画素値を加重加算する処理と同じである。
【００３６】
次に水平方向高周波数成分画像生成手段１ｈ及び垂直方向高周波数成分画像生成手段１ｖの動作について説明する。
水平方向高周波数成分画像生成手段１ｈは入力画像Ｄｉｎに対して、入力画像Ｄｉｎの各画素及びその水平方向近傍にある、例えば所定数の画素を用いたハイパスフィルタをかけて水平方向高周波数成分画像Ｄ１ｈを生成する。
一方、垂直方向高周波数成分画像生成手段１ｖは入力画像Ｄｉｎに対して、入力画像Ｄｉｎの各画素及びその垂直方向近傍にある、例えば所定数の画素を用いたハイパスフィルタをかけて垂直方向高周波数成分画像Ｄ１ｖを生成する。
【００３７】
ハイパスフィルタをかけることは高周波数成分を取り出すことに対応し、水平方向高周波数成分画像Ｄ１ｈには、入力画像Ｄｉｎの水平方向の高周波数成分（所定の水平周波数よりも高い成分から成る）が含まれ、垂直方向高周波数成分画像Ｄ１ｖには、入力画像Ｄｉｎの垂直方向の高周波数成分（所定の垂直方向周波数よりも高い成分から成る）が含まれる。
【００３８】
水平方向高周波数成分画像生成手段１ｈで行なわれるハイパスフィルタをかける処理としては、例えば、該手段１ｈへの入力信号からその水平方向の低周波数成分（もしくは各画素に対して水平方向に整列した所定数の画素から成る局所領域における画素値の単純平均値あるいは加重平均値）を差し引くことで高周波数成分を取り出す処理を行なうことができる。
同様に、垂直方向高周波数成分画像生成手段１ｖで行なわれるハイパスフィルタをかける処理としては、例えば、該手段１ｖへの入力信号からその垂直方向の低周波数成分（もしくは各画素に対して垂直方向に整列した所定数の画素から成る局所領域における画素値の単純平均値あるいは加重平均値）を差し引くことで高周波数成分を取り出す処理を行なうことができる。
【００３９】
次に画像拡大手段２Ｂｈ及び２Ｂｖの動作について説明する。画像拡大手段２Ｂｈは水平方向高周波数成分画像Ｄ１ｈを水平方向、垂直方向とも２倍に拡大した拡大画像Ｄ２Ｂｈを生成し、画像拡大手段２Ｂｖは垂直方向高周波数成分画像Ｄ１ｖを水平方向、垂直方向とも２倍に拡大した拡大画像Ｄ２Ｂｖを生成する。
【００４０】
画像拡大手段２Ｂｈ及び画像拡大手段２Ｂｖの各々は、図５を参照して説明した画像拡大手段２Ａと同様に構成することができる。従って、画像拡大手段２Ｂｈと画像拡大手段２Ｂｖとで構成される画像拡大手段２Ｂは図６のように示すことができる。
画像拡大手段２Ｂｈの入力は水平方向高周波数成分画像Ｄ１ｈであり、出力が拡大画像Ｄ２Ｂｈとなる。画像拡大手段２Ｂｖの入力は垂直方向高周波数成分画像Ｄ１ｖであり、出力が拡大画像Ｄ２Ｂｖである。
【００４１】
画像拡大手段２Ｂｈは、ゼロ挿入手段２１Ｂｈと、低周波数成分通過手段２２Ｂｈとを備え、画像拡大手段２Ｂｖは、ゼロ挿入手段２１Ｂｖと、低周波数成分通過手段２２Ｂｖとを備える。
ゼロ挿入手段２１Ｂｈ及びゼロ挿入手段２１Ｂｖの各々は、図５のゼロ挿入手段２１Ａと同様のものであり、低周波数成分通過手段２２Ｂｈ及び低周波数成分通過手段２２Ｂｖの各々は、図５の低周波数成分通過手段２２Ａと同様のものである。
【００４２】
ゼロ挿入手段２１Ｂｈから出力されるゼロ挿入画像Ｄ２１Ｂｈとゼロ挿入手段２１Ｂｖから出力されるゼロ挿入画像Ｄ２１Ｂｖとで、ゼロ挿入手段２１Ｂの出力としてのゼロ挿入画像Ｄ２１Ｂが構成される。
低周波数成分通過手段２２Ｂｈから出力される拡大画像Ｄ２Ｂｈと低周波数成分通過手段２２Ｂｖから出力される拡大画像Ｄ２Ｂｖとで、低周波数成分通過手段２２Ｂの出力としての拡大画像Ｄ２Ｂが構成される。低周波数成分通過手段２２Ｂの出力は、ゼロ挿入画像Ｄ２１Ｂの低周波数成分（Ｆｃ以下の成分）を取り出したものである。
【００４３】
次に高周波数成分画像生成手段３２Ａの動作について説明する。
水平方向高周波数成分画像生成手段３２Ａｈは、拡大画像Ｄ２Ｂｈに水平方向のハイパスフィルタをかけて所定の水平方向周波数以上の成分から成る、高周波数成分のみを取り出し、水平方向中間画像Ｄ３２Ａｈを生成する。
一方、垂直方向高周波数成分画像生成手段３２Ａｖは、拡大画像Ｄ２Ｂｖに垂直方向のハイパスフィルタをかけて所定の垂直方向周波数以上の成分から成る、高周波数成分のみを取り出し、垂直方向中間画像Ｄ３２Ａｖを生成する。
そして水平方向中間画像Ｄ３２Ａｈと垂直方向中間画像Ｄ３２Ａｖから成る中間画像Ｄ３２Ａが高周波数成分画像生成手段３２Ａから出力される。
【００４４】
水平方向高周波数成分画像生成手段３２Ａｈで行なわれるハイパスフィルタをかける処理は、水平方向高周波数成分画像生成手段１ｈにおける処理と同様に行ない、垂直方向高周波数成分画像生成手段３２Ａｖで行なわれるハイパスフィルタをかける処理は、垂直方向高周波数成分画像生成手段１ｖにおける処理と同様に行なうことができる。
即ち、水平方向高周波数成分画像生成手段３２Ａｈで行なわれるハイパスフィルタをかける処理としては、水平方向高周波数成分画像生成手段１ｈにおける処理と同様に、例えば、該手段３２Ａｈへの入力信号からその水平方向の低周波数成分（もしくは各画素に対して水平方向に整列した所定数の画素から成る局所領域における画素値の単純平均値あるいは加重平均値）を差し引くことで高周波数成分を取り出す処理を行なうことができる。
【００４５】
同様に、垂直方向高周波数成分画像生成手段３２Ａｖで行なわれるハイパスフィルタをかける処理としては、例えば、該手段３２Ａｖへの入力信号からその垂直方向の低周波数成分（もしくは各画素に対して垂直方向に整列した所定数の画素から成る局所領域における画素値の単純平均値あるいは加重平均値）を差し引くことで高周波数成分を取り出す処理を行なうことができる。
【００４６】
次に非線形処理手段３１の動作について説明する。非線形処理手段３１は、水平方向非線形処理手段３１ｈと、垂直方向非線形処理手段３１ｖを備える。水平方向非線形処理手段３１ｈと垂直方向非線形処理手段３１ｖとは互いに同様に構成されている。但し、水平方向非線形処理手段３１ｈは水平方向の処理を行ない、垂直方向非線形処理手段３１ｖは垂直方向の処理を行なう。
【００４７】
図７は水平方向非線形処理手段３１ｈの内部構成を表す図である。図示の水平方向非線形処理手段３１ｈは、ゼロクロス判定手段３１１ｈと、信号増幅手段３１２ｈを備える。
【００４８】
ゼロクロス判定手段３１１ｈは、入力される拡大画像Ｄ２Ｂｈにおける画素値の変化を水平方向に沿って確認する。そして画素値が正の値から負の値あるいは負の値から正の値に変化する箇所をゼロクロス点として捉え、信号Ｄ３１１ｈによってゼロクロス点の前後にある画素（図示の例では、直前及び直後の各１画素）の位置を信号増幅手段３１２ｈに伝達する。
なお、水平方向非線形処理手段３１ｈではゼロクロス点の左右に位置する画素がゼロクロス点の前後にある画素として認識される。
【００４９】
水平方向信号増幅手段３１２ｈは、水平方向ゼロクロス判定手段３１１ｈの判定結果に応じて決められる増幅率で第３の拡大画像Ｄ２Ｂｈの画素値を増幅する。具体的には、信号増幅手段３１２ｈは、信号Ｄ３１１ｈをもとにゼロクロス点の前後にある画素（ゼロクロス点を含む所定の領域内に存在する画素）を特定し、ゼロクロス点の前後にある画素についてのみその画素値を増幅させた（絶対値を大きくした）非線形処理画像Ｄ３１ｈを生成する。すなわちゼロクロス点の前後にある画素の画素値に対する増幅率を１より大きな値とし、それ以外の画素の画素値に対する増幅率は１とする。
このような処理により、水平方向に並んだ画素の信号値のステップ状の変化を含むエッジの鮮鋭化が行なわれる。
【００５０】
図８は垂直方向非線形処理手段３１ｖの内部構成を表す図である。図示の垂直方向非線形処理手段３１ｖは、ゼロクロス判定手段３１１ｖと、信号増幅手段３１２ｖを備える。
【００５１】
ゼロクロス判定手段３１１ｖは、入力される拡大画像Ｄ２Ｂｖにおける画素値の変化を垂直方向に沿って確認する。そして画素値が正の値から負の値あるいは負の値から正の値に変化する箇所をゼロクロス点として捉え、信号Ｄ３１１ｖによってゼロクロス点の前後にある画素（図示の例では、直前及び直後の各１画素）の位置を信号増幅手段３１２ｖに伝達する。
なお、垂直方向非線形処理手段３１ｖではゼロクロス点の上下に位置する画素がゼロクロス点の前後にある画素として認識される。
【００５２】
垂直方向信号増幅手段３１２ｖは、垂直方向ゼロクロス判定手段３１１ｖの判定結果に応じて決められる増幅率で第４の拡大画像Ｄ２Ｂｖの画素値を増幅する。具体的には、信号増幅手段３１２ｖは、信号Ｄ３１１ｖをもとにゼロクロス点の前後にある画素（ゼロクロス点を含む所定の領域内に存在する画素）を特定し、ゼロクロス点の前後にある画素についてのみその画素値を増幅させた（絶対値を大きくした）非線形処理画像Ｄ３１ｖを生成する。すなわちゼロクロス点の前後にある画素の画素値に対する増幅率を１より大きな値とし、それ以外の画素の画素値に対する増幅率は１とする。
このような処理により、垂直方向に並んだ画素の信号値のステップ状の変化を含むエッジの鮮鋭化が行なわれる。
【００５３】
なお、図示の例では、ゼロクロス点の直前及び直後の各々１個ずつの画素についてのみ、画素値を増幅させているが、ゼロクロス点の前及び後の所定数の画素、言い換えると、ゼロクロス点を含む所定の領域内に存在する画素について、画素値を増幅させることとしても良い。また上記「所定の領域」の大きさ（所定の領域に含まれる画素の数）を画像拡大手段２Ｂにおける画像の拡大率に応じて変える（拡大率に対して適切な値に定める）こととしても良い。
【００５４】
次に高周波数成分画像生成手段３２Ｂの動作について説明する。
水平方向高周波数成分画像生成手段３２Ｂｈは、非線形処理画像Ｄ３１ｈに水平方向のハイパスフィルタをかけて所定の水平方向周波数以上の成分から成る、高周波数成分のみを取り出し、水平方向中間画像Ｄ３２Ｂｈを生成する。一方、垂直方向高周波数成分画像生成手段３２Ｂｖは、非線形処理画像Ｄ３１ｖに垂直方向のハイパスフィルタをかけて所定の垂直方向周波数以上の成分から成る、高周波数成分のみを取り出し、垂直方向中間画像Ｄ３２Ｂｖを生成する。このようにして生成された中間画像Ｄ３２Ｂが高周波数成分画像生成手段３２Ｂから出力される。
【００５５】
水平方向高周波数成分画像生成手段３２Ｂｈで行なわれるハイパスフィルタをかける処理は、水平方向高周波数成分画像生成手段１ｈにおける処理と同様に行ない、垂直方向高周波数成分画像生成手段３２Ｂｖで行なわれるハイパスフィルタをかける処理は、垂直方向高周波数成分画像生成手段１ｖにおける処理と同様に行なうことができる。
即ち、水平方向高周波数成分画像生成手段３２Ｂｈで行なわれるハイパスフィルタをかける処理としては、例えば、該手段３２Ｂｈへの入力信号からその水平方向の低周波数成分（もしくは各画素に対して水平方向に整列した所定数の画素から成る局所領域における画素値の単純平均値あるいは加重平均値）を差し引くことで高周波数成分を取り出す処理を行なうことができる。
【００５６】
同様に、垂直方向高周波数成分画像生成手段３２Ｂｖで行なわれるハイパスフィルタをかける処理としては、例えば、該手段３２Ｂｖへの入力信号からその垂直方向の低周波数成分（もしくは各画素に対して垂直方向に整列した所定数の画素から成る局所領域における画素値の単純平均値あるいは加重平均値）を差し引くことで高周波数成分を取り出す処理を行なうことができる。
【００５７】
次に、高周波数成分画像補正手段３３Ａの詳細な動作について説明する。
【００５８】
図９は高周波数成分画像補正手段３３Ａの内部構成を表す図であり、図示の高周波数成分画像補正手段３３Ａは、増幅率計算手段３ＭＡと、画素値増幅手段３ＭＢとを有する。
【００５９】
増幅率計算手段３ＭＡは、拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ａの画素値の符号をもとに増幅率Ｄ３ＭＡを決定する。
画素値増幅手段３ＭＢは、増幅率計算手段３ＭＡにより決定された増幅率Ｄ３ＭＡで中間画像Ｄ３２Ａの画素値を増幅し、その結果を中間画像Ｄ３ＭＢとして出力する。高周波数成分画像補正手段３３Ａからは中間画像Ｄ３ＭＢが中間画像Ｄ３３Ａとして出力される。
以上のように増幅率Ｄ３ＭＡによって中間画像Ｄ３２Ａの画素値が増幅されるが、この動作は画素ごとに行われる。すわなち、増幅率計算手段３ＭＡは増幅率Ｄ３ＭＡを中間画像Ｄ３２Ａの各画素値に対して計算する。すなわち、中間画像Ｄ３２Ａにおいて水平座標がｘ、垂直座標がｙで表される画素の画素値は、拡大画像Ｄ２Ａにおいて水平座標がｘ、垂直座標がｙで表される画素の画素値を用いて計算された増幅率Ｄ３ＭＡでもって増幅される。
【００６０】
先に述べた通り、増幅率Ｄ３ＭＡに基づいて中間画像Ｄ３２Ａの画素値が増幅されるが、中間画像Ｄ３２Ａは、水平方向中間画像Ｄ３２Ａｈと垂直方向中間画像Ｄ３２Ａｖから成るので、増幅率Ｄ３ＭＡとしては、水平方向中間画像Ｄ３２Ａｈに対する増幅率Ｄ３ＭＡｈと垂直方向中間画像Ｄ３２Ａｖに対する増幅率Ｄ３ＭＡｖが決定される。即ち、増幅率計算手段３ＭＡ内の水平方向増幅率計算手段３ＭＡｈでは拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ａｈの画素値の符号をもとに増幅率Ｄ３ＭＡｈが決定され、垂直方向増幅率計算手段３ＭＡｖでは拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ａｖの画素値の符号をもとに増幅率Ｄ３ＭＡｖが決定され、増幅率計算手段３ＭＡからは増幅率Ｄ３ＭＡｈ及び増幅率Ｄ３ＭＡｖが増幅率Ｄ３ＭＡとして出力される。
【００６１】
図１０（ａ）〜（ｄ）を用いて増幅率計算手段３ＭＡの動作についてさらに詳しく説明する。増幅率計算手段３ＭＡは拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ａの画素値の符号をもとに増幅率Ｄ３ＭＡを計算する。
この増幅率は拡大画像Ｄ２Ａの画素値が大きな場合、小さい値となるよう計算される。図１０（ａ）〜（ｄ）は増幅率計算手段３ＭＡにおいて計算される増幅率Ｄ３ＭＡの一例を表しており、図１０（ａ）、（ｂ）は水平方向増幅率計算手段３ＭＡｈで計算される増幅率Ｄ３ＭＡｈを、図１０（ｃ）、（ｄ）は垂直方向増幅率計算手段３ＭＡｖで計算される増幅率Ｄ３ＭＡｖを表す。
【００６２】
増幅率Ｄ３ＭＡｈは、中間画像Ｄ３２Ａｈの画素値の符号が正であれば、図１０（ａ）に示すように、拡大画像Ｄ２Ａの画素値が大きい場合、小さな値をとる。一例として、傾きｋＭｈ１、ｋＭｈ２、ｋＭｈ３がすべて負の数であるとして、拡大画像Ｄ２Ａの画素値が０の場合はある所定の値ＧＭｈｂをとり、画素値が０からある値ＡＭｈ１の間は傾きｋＭｈ１で減少していき、画素値がＡＭｈ１からある値ＡＭｈ２の間は傾きｋＭｈ２で減少していき、画素値がＡＭｈ１以上の場合は傾きｋＭｈ３で減少していくといった形が考えられる。なお、明らかに増幅率は０以上であった方がよいので、上記の決定において増幅率が負の値になった場合はその値を０とする。
【００６３】
一方、増幅率Ｄ３ＭＡｈは、中間画像Ｄ３２Ａの画素値の符号が負であれば、図１０（ｂ）に示すように、拡大画像Ｄ２Ａの画素値が小さい場合、小さな値をとる。一例として、傾きｋＭｈ１、ｋＭｈ２、ｋＭｈ３がすべて正の数であるとして、拡大画像Ｄ２Ａの画素値が０の場合はある所定の値ＧＭｈｂをとり、画素値が０からある値ＡＭｈ１の間は傾きｋＭｈ１で増加していき、画素値がＡＭｈ１からある値ＡＭｈ２の間は傾きｋＭｈ２で増加していき、画素値がＡＭｈ１以上の場合は傾きｋＭｈ３で増加していくといった形が考えられる。なお、明らかに増幅率は０以上であった方がよいので、上記の決定において増幅率が負の値になった場合はその値を０とする。
【００６４】
この関係を、増幅率Ｄ３ＭＡｈをＧ、拡大画像Ｄ２Ａの画素値をＬとして表すと
【数１】
2011048670.tif 000003
と表される。ここで、傾きｋＭｈ１、ｋＭｈ２、ｋＭｈ３は中間画像Ｄ３２Ａｈの画素値の符号が正であれば負の値を、負であれば正の値をとるものとする。
また、所定の値ＧＭｈｂは水平方向中間画像Ｄ３２Ａｈの画素値の符号が正の場合と負の場合で異なる値としてもよい。
【００６５】
なお、増幅率Ｄ３ＭＡｈは、中間画像Ｄ３２Ａｈの画素値の符号が正で拡大画像Ｄ２Ａの画素値が大きい場合、小さな値となるよう制御し、中間画像Ｄ３２Ａの画素値の符号が負で拡大画像Ｄ２Ａの画素値が小さい場合、小さな値となるよう制御することが目的であるので、傾きｋＭｈ１、ｋＭｈ２、ｋＭｈ３の正、負がすべて同じである必要はなく、上記のような制御ができるのであればどのような組み合わせでもよい。さらに言えば、増幅率Ｄ３ＭＡｈと、拡大画像Ｄ２Ａの画素値及び中間画像Ｄ３２Ａｈの画素値の符号が上記のよう制御することができるのであれば式（２）における各係数のとる値は自由である。また、増幅率Ｄ３ＭＡｈの求め方も式（２）に制限されるものではない。
【００６６】
増幅率Ｄ３ＭＡｖは、中間画像Ｄ３２Ａｖの画素値の符号が正であれば、図１０（ｃ）に示すように、拡大画像Ｄ２Ａの画素値が大きい場合、小さな値をとる。一例として、傾きｋＭｖ１、ｋＭｖ２、ｋＭｖ３がすべて負の数であるとして、拡大画像Ｄ２Ａの画素値が０の場合はある所定の値ＧＭｖｂをとり、画素値が０からある値ＡＭｖ１の間は傾きｋＭｖ１で減少していき、画素値がＡＭｖ１からある値ＡＭｖ２の間は傾きｋＭｖ２で減少していき、画素値がＡＭｖ１以上の場合は傾きｋＭｖ３で減少していくといった形が考えられる。なお、明らかに増幅率は０以上であった方がよいので、上記の決定において増幅率が負の値になった場合はその値を０とする。
【００６７】
一方、増幅率Ｄ３ＭＡｖは、中間画像Ｄ３２Ａの画素値の符号が負であれば、図１０（ｄ）に示すように、拡大画像Ｄ２Ａの画素値が小さい場合、小さな値をとる。一例として、傾きｋＭｖ１、ｋＭｖ２、ｋＭｖ３がすべて正の数であるとして、拡大画像Ｄ２Ａの画素値が０の場合はある所定の値ＧＭｖｂをとり、画素値が０からある値ＡＭｖ１の間は傾きｋＭｖ１で増加していき、画素値がＡＭｖ１からある値ＡＭｖ２の間は傾きｋＭｖ２で増加していき、画素値がＡＭｖ１以上の場合は傾きｋＭｖ３で増加していくといった形が考えられる。なお、明らかに増幅率は０以上であった方がよいので、上記の決定において増幅率が負の値になった場合はその値を０とする。
【００６８】
この関係を、増幅率Ｄ３ＭＡｖをＧ、拡大画像Ｄ２Ａの画素値をＬとして表すと
【数２】
2011048670.tif 000004
と表される。ここで、傾きｋＭｖ１、ｋＭｖ２、ｋＭｖ３は中間画像Ｄ３２Ａの画素値の符号が正であれば負の値を、負であれば正の値をとるものとする。
また、所定の値ＧＭｖｂは垂直方向中間画像Ｄ３２Ａｖの画素値の符号が正の場合と負の場合で異なる値としてもよい。
【００６９】
なお、増幅率Ｄ３ＭＡｖは、中間画像Ｄ３２Ａｖの画素値の符号が正で拡大画像Ｄ２Ａの画素値が大きい場合、小さな値となるよう制御し、中間画像Ｄ３２Ａｖの画素値の符号が負で拡大画像Ｄ２Ａの画素値が小さい場合、小さな値となるよう制御することが目的であるので、傾きｋＭｖ１、ｋＭｖ２、ｋＭｖ３の正、負がすべて同じである必要はなく、上記のような制御ができるのであればどのような組み合わせでもよい。さらに言えば、増幅率Ｄ３ＭＡｖと、拡大画像Ｄ２Ａの画素値及び中間画像Ｄ３２Ａｖの画素値の符号が上記のよう制御することができるのであれば式（３）における各係数のとる値は自由である。また、増幅率Ｄ３ＭＡｖの求め方も式（３）に制限されるものではない。
以上が増幅率計算手段３ＭＡの動作である。
【００７０】
次に画素値増幅手段３ＭＢは、増幅率Ｄ３ＭＡに基づいて中間画像Ｄ３２Ａの画素値を増幅する。中間画像Ｄ３２Ａは水平方向中間画像Ｄ３２Ａｈと垂直方向中間画像Ｄ３２Ａｖとから成るので、画素値増幅手段３ＭＢは、水平方向中間画像Ｄ３２Ａｈの画素値を増幅するための水平方向画素値増幅手段３ＭＢｈと、垂直方向中間画像Ｄ３２Ａｖの画素値を増幅するための垂直方向画素値増幅手段３ＭＢｖとを有する。
【００７１】
水平方向画素値増幅手段３ＭＢｈは、増幅率Ｄ３ＭＡｈに基づいて水平方向中間画像Ｄ３２Ａｈの画素値を増幅した画像Ｄ３ＭＢｈを出力し、垂直方向画素値増幅手段３ＭＢｖは、増幅率Ｄ３ＭＡｖに基づいて垂直方向中間画像Ｄ３２Ａｖの画素値を増幅した画像Ｄ３ＭＢｖを出力する。先にも述べたようにこの動作は画素ごとに行われる。すなわち、水平方向中間画像Ｄ３２Ａｈにおいて水平座標がｘ、垂直座標がｙで表される画素の画素値は、拡大画像Ｄ２Ａにおいて水平座標がｘ、垂直座標がｙで表される画素の画素値と水平方向中間画像Ｄ３２Ａｈにおいて水平座標がｘ、垂直座標がｙで表される画素の画素値の符号を用いて計算された増幅率Ｄ３ＭＡｈでもって増幅され、垂直方向中間画像Ｄ３２Ａｖにおいて水平座標がｘ、垂直座標がｙで表される画素の画素値は、拡大画像Ｄ２Ａにおいて水平座標がｘ、垂直座標がｙで表される画素の画素値と垂直方向中間画像Ｄ３２Ａｖにおいて水平座標がｘ、垂直座標がｙで表される画素の画素値の符号を用いて計算された増幅率Ｄ３ＭＡｖでもって増幅される。そして画素値増幅手段３ＭＢからは画像Ｄ３ＭＢｈ及び画像Ｄ３ＭＢｖが画像Ｄ３ＭＢとして出力される。
【００７２】
そして画像Ｄ３ＭＢが中間画像Ｄ３３Ａとして高周波数成分画像補正手段３３Ａから出力される。中間画像Ｄ３３Ａは画像Ｄ３ＭＢｈに相当する水平方向中間画像Ｄ３３Ａｈと画像Ｄ３ＭＢｖに相当する垂直方向中間画像Ｄ３３Ａｖから成る。
以上が画素値増幅率手段３ＭＢの動作である。
【００７３】
以上が高周波数成分補正手段３３Ａの動作であり、高周波数成分画像補正手段３３Ａでは、拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ａの画素値の符号に応じて決められる増幅率（あるいはゲイン、利得）によって、中間画像Ｄ３２Ａに補正が加えられる。
【００７４】
次に、高周波数成分画像補正手段３３Ｂの詳細な動作について説明する。
図１１は高周波数成分画像補正手段３３Ｂの内部構成を表す図であり、図示の高周波数成分画像補正手段３３Ｂは、増幅率計算手段３ＨＡと、画素値増幅手段３ＨＢとを有する。
【００７５】
増幅率計算手段３ＨＡは、拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ｂの画素値の符号をもとに増幅率Ｄ３ＨＡを決定する。
画素値増幅手段３ＨＢは、増幅率計算手段３ＨＡにより決定された増幅率Ｄ３ＨＡで中間画像Ｄ３２Ｂの画素値を増幅し、その結果を中間画像Ｄ３ＨＢとして出力する。高周波数成分画像補正手段３３Ｂからは中間画像Ｄ３ＨＢが中間画像Ｄ３３Ｂとして出力される。
以上のように増幅率Ｄ３ＨＡによって中間画像Ｄ３２Ｂの画素値が増幅されるが、この動作は画素ごとに行われる。すわなち、増幅率計算手段３ＨＡは増幅率Ｄ３ＨＡを中間画像Ｄ３２Ｂの各画素値に対して計算する。すなわち、中間画像Ｄ３２Ｂにおいて水平座標がｘ、垂直座標がｙで表される画素の画素値は、拡大画像Ｄ２Ａにおいて水平座標がｘ、垂直座標がｙで表される画素の画素値を用いて計算された増幅率Ｄ３ＨＡでもって増幅される。
【００７６】
先に述べた通り、増幅率Ｄ３ＨＡに基づいて中間画像Ｄ３２Ｂの画素値が増幅されるが、中間画像Ｄ３２Ｂは、水平方向中間画像Ｄ３２Ｂｈと垂直方向中間画像Ｄ３２Ｂｖから成るので、増幅率Ｄ３ＨＡとしては、水平方向中間画像Ｄ３２Ｂｈに対する増幅率Ｄ３ＨＡｈと垂直方向中間画像Ｄ３２Ｂｖに対する増幅率Ｄ３ＨＡｖが決定される。即ち、増幅率計算手段３ＨＡ内の水平方向増幅率計算手段３ＨＡｈでは拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ｂｈの画素値の符号をもとに増幅率Ｄ３ＨＡｈが決定され、垂直方向増幅率計算手段３ＨＡｖでは拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ｂｖの画素値の符号をもとに増幅率Ｄ３ＨＡｖが決定され、増幅率計算手段３ＨＡからは増幅率Ｄ３ＨＡｈ及び増幅率Ｄ３ＨＡｖが増幅率Ｄ３ＨＡとして出力される。
【００７７】
図１２（ａ）〜（ｄ）を用いて増幅率計算手段３ＨＡの動作についてさらに詳しく説明する。増幅率計算手段３ＨＡは拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ｂの画素値の符号をもとに増幅率Ｄ３ＨＡを計算する。
この増幅率は拡大画像Ｄ２Ａの画素値が大きな場合、小さい値となるよう計算される。図１２（ａ）〜（ｄ）は増幅率計算手段３ＨＡにおいて計算される増幅率Ｄ３ＨＡの一例を表しており、図１２（ａ）、（ｂ）は水平方向増幅率計算手段３ＨＡｈで計算される増幅率Ｄ３ＨＡｈを、図１２（ｃ）、（ｄ）は垂直方向増幅率計算手段３ＨＡｖで計算される増幅率Ｄ３ＨＡｖを表す。
【００７８】
増幅率Ｄ３ＨＡｈは、中間画像Ｄ３２Ｂｈの画素値の符号が正であれば図１２（ａ）に示すように、拡大画像Ｄ２Ａの画素値が大きい場合、小さな値をとる。一例として、傾きｋＨｈ１、ｋＨｈ２、ｋＨｈ３がすべて負の数であるとして、拡大画像Ｄ２Ａの画素値が０の場合はある所定の値ＧＨｈｂをとり、画素値が０からある値ＡＨｈ１の間は傾きｋＨｈ１で減少していき、画素値がＡＨｈ１からある値ＡＨｈ２の間は傾きｋＨｈ２で減少していき、画素値がＡＨｈ１以上の場合は傾きｋＨｈ３で減少していくといった形が考えられる。なお、明らかに増幅率は０以上であった方がよいので、上記の決定において増幅率が負の値になった場合はその値を０とする。
【００７９】
一方、増幅率Ｄ３ＨＡｈは、中間画像Ｄ３２Ｂｈの画素値の符号が負であれば、図１２（ｂ）に示すように、拡大画像Ｄ２Ａの画素値が小さい場合、小さな値をとる。一例として、傾きｋＨｈ１、ｋＨｈ２、ｋＨｈ３がすべて正の数であるとして、拡大画像Ｄ２Ａの画素値が０の場合はある所定の値ＧＨｈｂをとり、画素値が０からある値ＡＨｈ１の間は傾きｋＨｈ１で増加していき、画素値がＡＨｈ１からある値ＡＨｈ２の間は傾きｋＨｈ２で増加していき、画素値がＡＨｈ１以上の場合は傾きｋＨｈ３で増加していくといった形が考えられる。なお、明らかに増幅率は０以上であった方がよいので、上記の決定において増幅率が負の値になった場合はその値を０とする。
【００８０】
この関係を、増幅率Ｄ３ＨＡｈをＧ、拡大画像Ｄ２Ａの画素値をＬとして表すと
【数３】
2011048670.tif 000005
と表される。ここで、傾きｋＨｈ１、ｋＨｈ２、ｋＨｈ３は中間画像Ｄ３２Ｂｈの画素値の符号が正であれば負の値を、負であれば正の値をとるものとする。
また、所定の値ＧＨｈｂは水平方向中間画像Ｄ３２Ｂｈの画素値の符号が正の場合と負の場合で異なる値としてもよい。
【００８１】
なお、増幅率Ｄ３ＨＡｈは、中間画像Ｄ３２Ｂｈの画素値の符号が正で拡大画像Ｄ２Ａの画素値が大きい場合、小さな値となるよう制御し、中間画像Ｄ３２Ｂｈの画素値の符号が負で拡大画像Ｄ２Ａの画素値が小さい場合、小さな値となるよう制御することが目的であるので、傾きｋＨｈ１、ｋＨｈ２、ｋＨｈ３の正、負がすべて同じである必要はなく、上記のような制御ができるのであればどのような組み合わせでもよい。さらに言えば、増幅率Ｄ３ＨＡｈと、拡大画像Ｄ２Ａの画素値及び中間画像Ｄ３２Ｂｈの画素値の符号が上記のよう制御することができるのであれば式（４）における各係数のとる値は自由である。また、増幅率Ｄ３ＨＡｈの求め方も式（４）に制限されるものではない。
【００８２】
増幅率Ｄ３ＨＡｖは、中間画像Ｄ３２Ｂｖの画素値の符号が正であれば、図１２（ｃ）に示すように、拡大画像Ｄ２Ａの画素値が大きい場合、小さな値をとる。一例として、傾きｋＨｖ１、ｋＨｖ２、ｋＨｖ３がすべて負の数であるとして、拡大画像Ｄ２Ａの画素値が０の場合はある所定の値ＧＨｖｂをとり、画素値が０からある値ＡＨｖ１の間は傾きｋＨｖ１で減少していき、画素値がＡＨｖ１からある値ＡＨｖ２の間は傾きｋＨｖ２で減少していき、画素値がＡＨｖ１以上の場合は傾きｋＨｖ３で減少していくといった形が考えられる。なお、明らかに増幅率は０以上であった方がよいので、上記の決定において増幅率が負の値になった場合はその値を０とする。
【００８３】
一方、増幅率Ｄ３ＨＡｖは、中間画像Ｄ３２Ｂの画素値の符号が負であれば、図１２（ｄ）に示すように、拡大画像Ｄ２Ａの画素値が小さい場合、小さな値をとる。一例として、傾きｋＨｖ１、ｋＨｖ２、ｋＨｖ３がすべて正の数であるとして、拡大画像Ｄ２Ａの画素値が０の場合はある所定の値ＧＨｖｂをとり、画素値が０からある値ＡＨｖ１の間は傾きｋＨｖ１で増加していき、画素値がＡＨｖ１からある値ＡＨｖ２の間は傾きｋＨｖ２で増加していき、画素値がＡＨｖ１以上の場合は傾きｋＨｖ３で増加していくといった形が考えられる。なお、明らかに増幅率は０以上であった方がよいので、上記の決定において増幅率が負の値になった場合はその値を０とする。
【００８４】
この関係を、増幅率Ｄ３ＨＡｖをＧ、拡大画像Ｄ２Ａの画素値をＬとして表すと
【数４】
2011048670.tif 000006
と表される。ここで、傾きｋＨｖ１、ｋＨｖ２、ｋＨｖ３は中間画像Ｄ３２Ｂｖの画素値の符号が正であれば負の値を、負であれば正の値をとるものとする。
また、所定の値ＧＨｖｂは垂直方向中間画像Ｄ３２Ｂｖの画素値の符号が正の場合と負の場合で異なる値としてもよい。
【００８５】
なお、増幅率Ｄ３ＨＡｖは、中間画像Ｄ３２Ｂｖの画素値の符号が正で
拡大画像Ｄ２Ａの画素値が大きい場合、小さな値となるよう制御し、中間画像Ｄ３２Ｂの画素値の符号が負で拡大画像Ｄ２Ａの画素値が小さい場合、小さな値となるよう制御することが目的であるので、傾きｋＨｖ１、ｋＨｖ２、ｋＨｖ３の正、負がすべて同じである必要はなく、上記のような制御ができるのであればどのような組み合わせでもよい。さらに言えば、増幅率Ｄ３ＨＡｖと、拡大画像Ｄ２Ａの画素値及び中間画像Ｄ３２Ｂｖの画素値の符号が上記のよう制御することができるのであれば式（５）における各係数のとる値は自由である。また、増幅率Ｄ３ＭＡｖの求め方も式（５）に制限されるものではない。
以上が増幅率計算手段３ＨＡの動作である。
【００８６】
次に画素値増幅手段３ＨＢは、増幅率Ｄ３ＨＡに基づいて中間画像Ｄ３２Ｂの画素値を増幅する。中間画像Ｄ３２Ｂは水平方向中間画像Ｄ３２Ｂｈと垂直方向中間画像Ｄ３２Ｂｖとから成るので、画素値増幅手段３ＨＢは、水平方向中間画像Ｄ３２Ｂｈの画素値を増幅するための水平方向画素値増幅手段３ＨＢｈと、垂直方向中間画像Ｄ３２Ｂｖの画素値を増幅するための垂直方向画素値増幅手段３ＨＢｖとを有する。
【００８７】
水平方向画素値増幅手段３ＨＢｈは、増幅率Ｄ３ＨＡｈに基づいて水平方向中間画像Ｄ３２Ｂｈの画素値を増幅した画像Ｄ３ＨＢｈを出力し、垂直方向画素値増幅手段３ＨＢｖは、増幅率Ｄ３ＨＡｖに基づいて垂直方向中間画像Ｄ３２Ｂｖの画素値を増幅した画像Ｄ３ＨＢｖを出力する。先にも述べたようにこの動作は画素ごとに行われる。すなわち、水平方向中間画像Ｄ３２Ｂｈにおいて水平座標がｘ、垂直座標がｙで表される画素の画素値は、拡大画像Ｄ２Ａにおいて水平座標がｘ、垂直座標がｙで表される画素の画素値と水平方向中間画像Ｄ３２Ｂｈにおいて水平座標がｘ、垂直座標がｙで表される画素の画素値の符号を用いて計算された増幅率Ｄ３ＨＡｈでもって増幅され、垂直方向中間画像Ｄ３２Ｂｖにおいて水平座標がｘ、垂直座標がｙで表される画素の画素値は、拡大画像Ｄ２Ａにおいて水平座標がｘ、垂直座標がｙで表される画素の画素値と垂直方向中間画像Ｄ３２Ｂｖにおいて水平座標がｘ、垂直座標がｙで表される画素の画素値の符号を用いて計算された増幅率Ｄ３ＨＡｖでもって増幅される。そして画素値増幅手段３ＨＢからは画像Ｄ３ＨＢｈ及び画像Ｄ３ＨＢｖが画像Ｄ３ＨＢとして出力される。
【００８８】
そして画像Ｄ３ＨＢが中間画像Ｄ３３Ｂとして高周波数成分画像処理手段３３Ｂから出力される。中間画像Ｄ３３Ｂは画像Ｄ３ＨＢｈに相当する水平方向中間画像Ｄ３３Ｂｈと画像Ｄ３ＨＢｖに相当する垂直方向中間画像Ｄ３３Ｂｖから成る。
以上が画素値増幅率手段３ＨＢの動作である。
【００８９】
以上が高周波数成分補正手段３３Ｂの動作であり、高周波数成分画像補正手段３３Ｂでは、拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ｂの画素値の符号に応じて決められる増幅率（あるいはゲイン、利得）によって、中間画像Ｄ３２Ｂに補正が加えられる。
【００９０】
次に加算手段３４の動作について説明する。加算手段３４は、中間画像Ｄ３３Ａと中間画像Ｄ３３Ｂを加算した結果を高周波数成分画像Ｄ３として出力する。
【００９１】
ここで中間画像Ｄ３３Ａは水平方向中間画像Ｄ３３Ａｈ及び垂直方向中間画像Ｄ３３Ａｖから成り、中間画像Ｄ３３Ｂは水平方向中間画像Ｄ３３Ｂｈ及び垂直方向中間画像Ｄ３３Ｂｖから成るので、中間画像Ｄ３３Ａと中間画像Ｄ３３Ｂを加算するとは、水平方向中間画像Ｄ３３Ａｈ、垂直方向中間画像Ｄ３３Ａｖ、水平方向中間画像Ｄ３３Ｂｈ、垂直方向中間画像Ｄ３３Ｂｖを加算することを意味する。なおここでの加算処理は、単純加算に限らず、各画像に個別の重みを付けて加算する処理であっても良い。
【００９２】
最後に加算手段４の動作について説明する。加算手段４は、拡大画像Ｄ２Ａと高周波数成分画像Ｄ３を加算する。そして加算手段４において拡大画像Ｄ２Ａと高周波数成分画像Ｄ３を加算した結果得られた画像が、最終的な拡大画像Ｄｏｕｔとして画像処理装置から出力される。なおここでの加算処理は、単純加算に限らず、各画像に個別の重みを付けて加算する処理であっても良い。
【００９３】
以下、本発明における画像処理装置の作用、効果について説明する。
本発明の実施の形態では、中間画像Ｄ３２Ａ及びＤ３２Ｂをそのまま加算して高周波数成分画像Ｄ３を生成して拡大画像Ｄ２Ａに加算するわけではなく、中間画像Ｄ３２Ａ及びＤ３２Ｂを高周波数成分画像補正手段３３Ａ及び３３Ｂで補正した後に加算して高周波数成分画像Ｄ３を生成して、拡大画像Ｄ２Ａに加算しているが、以下、仮に中間画像Ｄ３２Ａ及びＤ３２Ｂをそのまま加算して高周波数成分画像Ｄ３を生成して、拡大画像Ｄ２Ａに加算した場合に得られる効果について説明し、その後で、中間画像Ｄ３２Ａ及びＤ３２Ｂの代わりに、中間画像Ｄ３３Ａ及びＤ３３Ｂを加算することによる効果について説明する。
【００９４】
拡大画像Ｄ２Ａは、入力画像Ｄｉｎのナイキスト周波数Ｆｎより低い周波数に相当する周波数成分を含み、高周波数成分画像Ｄ３は入力画像Ｄｉｎのナイキスト周波数Ｆｎ以上の周波数に相当する周波数成分を含む。従って、拡大画像Ｄ２Ａと高周波数成分画像Ｄ３を加算して生成される拡大画像Ｄｏｕｔは画像拡大後のナイキスト周波数に至る全ての周波数領域にわたって周波数成分を持つことになる。
【００９５】
まず、拡大画像Ｄ２Ａが入力画像Ｄｉｎのナイキスト周波数Ｆｎより低い周波数に相当する周波数成分を持つことについて説明する。
【００９６】
図１３（ａ）〜（ｄ）は入力画像Ｄｉｎから拡大画像Ｄ２Ａを生成する際の作用を模式的に表した図であり、図１３（ａ）は入力画像Ｄｉｎの周波数スペクトルを、図１３（ｂ）はゼロ挿入画像Ｄ２１Ａの周波数スペクトルを、図１３（ｃ）は低周波数成分通過手段２２Ａの周波数応答を、図１３（ｄ）は拡大画像Ｄ２Ａの周波数スペクトルを表している。
【００９７】
入力画像Ｄｉｎの周波数スペクトルについて説明する。入力画像Ｄｉｎからは通常、自然画などが入力されるが、これらの画像のスペクトル強度は周波数空間の原点周辺に集中している。従って入力画像Ｄｉｎの周波数スペクトルは図１３（ａ）のように表すことが出来る。ここで図１３（ａ）の縦軸はスペクトル強度を、横軸は空間周波数を、Ｆｎは入力画像Ｄｉｎのナイキスト周波数を表している。
【００９８】
なお、通常入力画像Ｄｉｎは２次元の画像のため、その周波数スペクトルも２次元の周波数空間で表されるが、その形状は図１３（ａ）に示した周波数スペクトルが原点を中心に等方的に広がったものとなる。従って周波数スペクトルについて説明するためには最低限、１次元分の形状を示せばよく、今後、特に断らない限り、周波数空間の形状は１次元分のみ示して説明を行う。
【００９９】
次にゼロ挿入画像Ｄ２１Ａの周波数スペクトルについて説明する。入力画像Ｄｉｎに対してゼロ挿入手段２１Ａで（入力画像Ｄｉｎの）１画素につき１画素、画素値０を持った画素を挿入することで周波数空間上では周波数Ｆｎを中心にした折り返しが発生する。その結果、ゼロ挿入画像Ｄ２１Ａの周波数スペクトルは図１３（ｂ）のようになる。
【０１００】
次に低周波数成分通過手段２２Ａの周波数応答について説明する。先に述べたように低周波数成分通過手段２２Ａにおける演算はローパスフィルタ処理となっているので、図１３（ｃ）に示すように低周波数成分通過手段２２Ａの周波数応答は、周波数が高くなるほど低くなる。図示の例では、低周波数成分通過手段２２Ａが第１の周波数Ｆａ（＝Ｆｎ）以下の周波数成分を主に通過させるものとしている。
【０１０１】
最後に拡大画像Ｄ２Ａの周波数スペクトルについて説明する。図１３（ｂ）に示した周波数スペクトルを持つゼロ挿入画像Ｄ２１Ａが図１３（ｃ）に示した周波数応答を持った低周波数成分通過手段２２Ａを通ることで拡大画像Ｄ２Ａが生成される。従って拡大画像Ｄ２Ａの周波数スペクトルはゼロ挿入画像Ｄ２１Ａの周波数スペクトルから、斜線で示した高周波数側の領域Ｒ２ＡＨが除かれたものとなる。
【０１０２】
従って、拡大画像Ｄ２Ａは主に入力画像Ｄｉｎのナイキスト周波数Ｆｎより低い周波数に相当する周波数成分を持つことになる。
【０１０３】
次に高周波数成分画像Ｄ３が主に入力画像Ｄｉｎのナイキスト周波数Ｆｎ以上の周波数に相当する周波数成分を持つことについて説明する。仮に高周波数成分画像補正手段３３Ａ及び３３Ｂによる補正を行わないとすれば、高周波数成分画像Ｄ３は中間画像Ｄ３２Ａと中間画像Ｄ３２Ｂを加算して得られるが、中間画像Ｄ３２Ａは特に入力画像Ｄｉｎのナイキスト周波数Ｆｎに近い周波数に相当する周波数成分を持ち、中間画像Ｄ３２Ｂは特に入力画像Ｄｉｎのナイキスト周波数Ｆｎより高い周波数に相当する周波数成分を持ち、高周波数成分画像Ｄ３では中間画像Ｄ３２Ａ、Ｄ３２Ｂがもつ周波数成分が加算されるので、入力画像Ｄｉｎのナイキスト周波数Ｆｎ以上の周波数成分を持つことになる。
【０１０４】
まず、中間画像Ｄ３２Ａの周波数スペクトルについて説明する。
図１４（ａ）〜（ｆ）は中間画像Ｄ３２Ａを生成する際の作用を模式的に表した図であり、図１４（ａ）は高周波数成分画像生成手段１の周波数応答を、図１４（ｂ）は高周波数成分画像Ｄ１（又はＤ１ｈ若しくはＤ１ｖ）の周波数スペクトルを、図１４（ｃ）は画像拡大手段２Ｂ内のゼロ挿入手段２１Ｂによって生成されるゼロ挿入画像Ｄ２１Ｂ（又はＤ２１Ｂｈ若しくはＤ２１Ｂｖ）の周波数スペクトルを、図１４（ｄ）は拡大画像Ｄ２Ｂ（又はＤ２Ｂｈ若しくはＤ２Ｂｖ）の周波数スペクトルを、図１４（ｅ）は高周波数成分画像生成手段３２Ａ（又は３２Ａｈ若しくは３２Ａｖ）の周波数応答を、図１４（ｆ）は高周波数成分画像生成手段３２Ａから出力される中間画像Ｄ３２Ａ（又はＤ３２Ａｈ若しくはＤ３２Ａｖ）の周波数スペクトルを表している。
【０１０５】
まず、高周波数成分画像生成手段１の周波数応答及び高周波数成分画像Ｄ１の周波数スペクトルについて説明する。高周波数成分画像生成手段１は入力画像Ｄｉｎのうち、主に所定の周波数Ｆｂ以上の成分を通過させるハイパスフィルタを用いて高周波数成分画像Ｄ１を生成するので、図１４（ａ）に示すように高周波数成分画像生成手段１の周波数応答は、周波数が高くなるほど高くなる。図１３（ａ）に示した周波数スペクトルを持つ入力画像Ｄｉｎが図１４（ａ）に示す周波数応答をもつハイパスフィルタを通過することで高周波数成分画像Ｄ１が得られる。図示の例では、高周波数成分画像Ｄ１の周波数スペクトルは図１４（ｂ）に示すように周波数が低い領域（周波数Ｆｂよりも低い領域）では小さくなり、周波数が高い領域（周波数Ｆｂ以上の領域）でのみある程度の強度をもつことになる。
【０１０６】
次に画像拡大手段２Ｂ内のゼロ挿入画像Ｄ２１Ｂの周波数スペクトルについて説明する。先に画像拡大手段２Ａのゼロ挿入手段２１Ａについて説明したのと同様に、ゼロ挿入手段２１Ｂによって折り返しが発生するので、画像拡大手段２Ｂ内のゼロ挿入画像Ｄ２１Ｂの周波数スペクトルは図１４（ｃ）のようになる。
【０１０７】
次に拡大画像Ｄ２Ｂの周波数スペクトルについて説明する。
拡大画像Ｄ２Ｂを生成する際、ゼロ挿入画像Ｄ２１Ｂの高周波数成分側の周波数スペクトル（例えば所定の周波数Ｆｃよりも高い領域の成分）が、低周波数成分通過手段２２Ｂによって取り除かれるので、拡大画像Ｄ２Ｂの周波数スペクトルは図１４（ｄ）に示すように高周波数側の領域（周波数Ｆｃよりも高い領域）Ｒ３２ＡＨが取り除かれたものとなる。
【０１０８】
最後に高周波数成分画像生成手段３２Ａの周波数応答及び中間画像Ｄ３２Ａの周波数スペクトルについて説明する。高周波数成分画像生成手段３２Ａは主に所定の周波数Ｆｄ以上の成分を通過させるハイパスフィルタとなっているのでその周波数応答は図１４（ｅ）に示すように周波数が高くなるほど高くなる。中間画像Ｄ３２Ａは、図１４（ｄ）に示した周波数スペクトルをもつ拡大画像Ｄ２Ｂが、図１４（ｅ）に示した周波数応答をもつハイパスフィルタを通過することで生成される。従って中間画像Ｄ３２Ａの周波数応答は図１４（ｆ）に示すように、図１４（ｄ）に示した拡大画像Ｄ２Ｂの周波数スペクトルからさらに低周波数側の領域（周波数Ｆｄよりも低い領域）Ｒ３２ＡＬが取り除かれたものとなる。
【０１０９】
従って中間画像Ｄ３２Ａは主に入力画像Ｄｉｎのナイキスト周波数Ｆｎに近い周波数に相当する周波数成分(周波数Ｆｄから周波数Ｆｃまでの周波数成分)を持つことになる。
【０１１０】
次に中間画像Ｄ３２Ｂの周波数スペクトルについて説明する。
図１５（ａ）〜（ｃ）は中間画像Ｄ３２Ｂを生成する際の作用を模式的に表した図であり、図１５（ａ）は非線形処理手段３１（又は３１ｈ若しくは３１ｖ）により高周波数成分が生成される様子を、図１５（ｂ）は高周波数成分画像生成手段３２Ｂの周波数応答を、図１５（ｃ）は中間画像Ｄ３２Ｂの周波数スペクトルを表している。
【０１１１】
後述するように、非線形処理画像Ｄ３１には入力画像Ｄｉｎのナイキスト周波数Ｆｎより高い周波数に相当する高周波数成分が生成される。図１５（ａ）はその様子を模式的に表した図である。図示の例では、周波数Ｆｅ以上の成分が生成されている。中間画像Ｄ３２Ｂは非線形処理画像Ｄ３１が高周波数成分画像生成手段３２Ｂを通過することで生成される。高周波数成分画像生成手段３２Ｂは主に所定の周波数Ｆｆ以上の成分を通過させるハイパスフィルタでありその周波数応答は図１５（ｂ）に示すように、周波数が高くなるほど高くなっている。従って中間画像Ｄ３２Ｂの周波数スペクトルは図１５（ｃ）に示すように、非線形処理画像Ｄ３１の周波数スペクトルから低周波数側の領域Ｒ３２ＢＬが取り除かれたものとなるので、入力画像Ｄｉｎのナイキスト周波数Ｆｎより高い周波数に相当するものとなる。
【０１１２】
図１６（ａ）〜（ｅ）、図１７（ａ）〜（ｆ）を用いて中間画像Ｄ３２Ｂの周波数スペクトルについてより詳しく説明を行う。なお、説明を簡単にするため各々１次元信号として記載した。
【０１１３】
図１６（ａ）〜（ｅ）は、輝度、彩度などの成分値がステップ状に変化する画像（ステップ画像）を表すステップエッジ信号と、該ステップエッジ信号を互いに異なるサンプリング周波数でサンプリングした際に得られる信号及びその高周波数成分信号の信号強度を表している。図１６（ａ）はステップエッジ信号を表す。
【０１１４】
図１６（ｂ）はステップエッジ信号をサンプリング間隔Ｓ１でサンプリングして得られる信号、図１６（ｃ）はステップエッジ信号を間隔Ｓ１でサンプリングして得られる信号の高周波数成分を表し、図１６（ｄ）はステップエッジ信号をサンプリング間隔Ｓ２でサンプリングして得られる信号、図１６（ｅ）はステップエッジ信号を間隔Ｓ１でサンプリングして得られる信号の高周波数成分を表す。
なお、サンプリング間隔Ｓ１はサンプリング間隔Ｓ２より短くなっており、サンプリング間隔を短くすることは画像を拡大することと同じである。
なお、図１６を用いて説明する、ステップエッジ信号に対するサンプリング間隔Ｓ１、Ｓ２と高周波数成分の関係は特定のサンプリング間隔の組み合わせに依存する話ではないが、以下、サンプリング間隔Ｓ２は入力画像Ｄｉｎのサンプリング間隔と同じであり、サンプリング間隔Ｓ１はサンプリング間隔Ｓ２の半分であるとする。
【０１１５】
図１６（ｂ）、（ｃ）及び図１６（ｄ）、（ｅ）に示されるようにエッジの中央は高周波数成分信号（図１６（ｃ）、（ｅ））においてゼロクロス点Ｚとして現れる。また、図１６（ｂ）、（ｃ）と図１６（ｄ）、（ｅ）を比較すると明らかなように、ゼロクロス点Ｚの前後での高周波数成分信号の傾きはサンプリング間隔を短くするにつれて（あるいは画像を拡大させるのに応じて）急になり、かつゼロクロス点Ｚの近傍で高周波数成分の局所的な最大値、最小値を与える点の位置もゼロクロス点Ｚに近づく。
【０１１６】
従って画像を拡大する際、入力画像Ｄｉｎの高周波数成分を取り出し、その変化をゼロクロス近傍で急峻にし、かつゼロクロス近傍で局所的な最大値、最小値を与える点をゼロクロス点に近づけることで入力画像Ｄｉｎの解像度には含まれない（あるいは入力画像Ｄｉｎのナイキスト周波数より高い）高周波数成分を生成し、これによりエッジの鮮鋭化が可能となる。
【０１１７】
図１７（ａ）〜（ｆ）は、高周波数成分画像生成手段１、画像拡大手段２Ｂ、非線形処理手段３１及び高周波数成分画像生成手段３２Ｂによる高周波数成分生成の手順を模式的に表した図であり、図１７（ａ）は輝度、彩度などの成分値がステップ状に変化する画像（ステップ画像）、図１７（ｂ）はステップ画像に対応した入力画像Ｄｉｎ、図１７（ｃ）は高周波数成分画像Ｄ１、図１７（ｄ）は拡大画像Ｄ２Ｂ、図１７（ｅ）は非線形処理画像Ｄ３１、図１７（ｆ）は中間画像Ｄ３２Ｂを表す。
【０１１８】
また、図１７（ａ）〜（ｅ）において座標Ｐ３はエッジ近傍において信号強度が低い値をとる領域（低レベル側）の境界に相当する画素であり、座標Ｐ４は高い値をとる領域（高レベル側）の境界に相当する画素である。
【０１１９】
ステップ画像に対応した入力画像Ｄｉｎ、高周波数成分画像Ｄ１については図１６（ａ）〜（ｅ）で説明した通りであり、その説明は省略し、まず拡大画像Ｄ２Ｂの説明を行う。
【０１２０】
なお、高周波数成分画像Ｄ１においてゼロクロス点Ｚ近傍での局所的な最大値は、高レベル側の境界に現われるので座標Ｐ４で表される画素に局所的な最大値が現われ、逆に局所的な最小値は、低レベル側の境界に現われるので座標Ｐ３で表される画素に局所的な最小値が現われる。
【０１２１】
拡大画像Ｄ２Ｂ（又はＤ２Ｂｈ若しくはＤ２Ｂｖ）はゼロ挿入手段２１Ｂで高周波数成分画像Ｄ１に対して（画像Ｄ１の）１画素につき１画素、画素値０をもった画素を挿入した後、低周波数成分通過手段２２Ｂでその低周波数成分を取り出すことで得られる。低周波数成分を取り出すことは高周波数成分画像Ｄ１（図１７（ｃ））について局所領域における平均的な画素値を求めることと同じであり、拡大画像Ｄ２Ｂ（又はＤ２Ｂｈ若しくはＤ２Ｂｖ）は、図１７（ｄ）に示したように、高周波数成分画像Ｄ１とほぼ同じ形をした、サンプリング数の増えた信号となる。
【０１２２】
なお、画像が拡大されるのでゼロクロス点Ｚと座標Ｐ３で表される画素の間に新たに座標Ｐ１で表される画素が、ゼロクロス点Ｚと座標Ｐ４で表される画素の間に新たに座標Ｐ２で表される画素が現われる。また、拡大画像Ｄ２Ｂにおいてもゼロクロス点Ｚ近傍での局所的な最大値は座標Ｐ４で表される画素に、局所的な最小値は座標Ｐ３で表される画素に現われる。
【０１２３】
次に非線形処理画像Ｄ３１の説明を行う。非線形処理画像Ｄ３１は、非線形処理手段３１が拡大画像Ｄ１中のゼロクロス点Ｚを検出し、そのゼロクロス点Ｚの前後の画素の画素値を増幅した結果として出力される。従って非線形処理画像Ｄ３１（又はＤ３１ｈ若しくはＤ３１ｖ）では、座標Ｐ１、Ｐ２で表される画素の画素値が増幅されることになり、非線形処理画像Ｄ３１は図１７（ｅ）に示したような信号となる。
【０１２４】
最後に中間画像Ｄ３２Ｂの説明を行う。中間画像Ｄ３２Ｂ（図１７（ｆ））は非線形処理画像Ｄ３１（図１７(ｅ)）のもつ高周波数成分が高周波数成分画像生成手段３２Ｂにて取り出されたものである。高周波数成分は、入力信号から入力信号の低周波数成分（もしくは局所領域における画素値の単純平均値あるいは加重平均値）を差し引くことで取り出すことができる。
【０１２５】
非線形処理画像Ｄ３１（図１７（ｅ））ではゼロクロス点Ｚの前後の画素
（座標Ｐ１、Ｐ２で表される画素）については、その画素値が信号増幅手段３１２ｈ、３１２ｖにて増幅されているため、局所領域における平均的な画素値からの差は大きくなる。一方、ゼロクロス点近傍のその他の画素については、その画素値が増幅されることはないので、局所領域における平均的な画素値からの差は小さな値となる。従って拡大画像Ｄ２Ｂ（図１７（ｄ））と比較すると中間画像Ｄ３２Ｂ（図１７（ｆ））では、ゼロクロス点Ｚの近傍での局所的な最大値、最小値を与える点はそれぞれ座標Ｐ２、Ｐ１で表される画素となり、
よりゼロクロス点Ｚへと近づく。また、局所的な最大値、最小値を与える点がゼロクロス点Ｚへと近づいた分、ゼロクロス点近傍での信号の変化も急になる。
【０１２６】
先に説明したようにこれは中間画像Ｄ３２Ｂに、入力画像Ｄｉｎの解像度には含まれない高周波数成分が含まれることを意味する。言い換えると非線形処理手段３１において、拡大画像Ｄ２Ｂのゼロクロス点前後の画素値を増幅することで、入力画像Ｄｉｎのナイキスト周波数Ｆｎより高い周波数に対応した高周波数成分を生成したことになる。
さらに言い換えると、非線形処理画像生成手段３０において拡大画像Ｄ２Ｂに対し非線形処理を行うことで、入力画像Ｄｉｎのナイキスト周波数Ｆｎより高い周波数に対応した高周波数成分を生成することが出来る。
【０１２７】
また、中間画像Ｄ３２Ｂは非線形処理手段３１において生成した高周波数成分を高周波数成分画像生成手段３２Ｂで取り出すことで生成されるので、入力画像Ｄｉｎのナイキスト周波数Ｆｎより高い周波数に対応した高周波数成分をもつ画像となる。
【０１２８】
拡大画像Ｄ２Ａ、中間画像Ｄ３２Ａ、中間画像Ｄ３２Ｂが持つ周波数成分を図示すると図２１のようになる。拡大画像Ｄ２Ａには主に入力画像Ｄｉｎのナイキスト周波数Ｆｎより低い周波数に相当する領域ＲＬに対応した周波数成分が含まれている。一方、中間画像Ｄ３２Ａには入力画像Ｄｉｎのナイキスト周波数Ｆｎに近い周波数に相当する領域ＲＭに対応した周波数成分が含まれており、中間画像Ｄ３２Ｂには入力画像Ｄｉｎのナイキスト周波数Ｆｎより高い周波数に相当する領域ＲＨに対応した周波数成分が含まれている。
【０１２９】
中間画像Ｄ３２Ａと中間画像Ｄ３２Ｂを加算して高周波数成分画像Ｄ３を生成すれば、高周波数成分画像Ｄ３には、中間画像Ｄ３２Ａが持つ周波数成分と中間画像Ｄ３２Ｂが持つ周波数成分の双方が含まれることになる。中間画像Ｄ３２Ａには入力画像Ｄｉｎのナイキスト周波数Ｆｎに近い周波数に相当する周波数成分が含まれており、中間画像Ｄ３２Ｂには入力画像Ｄｉｎのナイキスト周波数Ｆｎより高い周波数に相当する周波数成分が含まれているので、高周波数成分画像Ｄ３には入力画像Ｄｉｎのナイキスト周波数Ｆｎ以上の周波数成分が含まれることになる。そして、拡大画像Ｄ２Ａに高周波数成分画像Ｄ３を加算し出力画像Ｄｏｕｔを得ることで、出力画像Ｄｏｕｔに対して入力画像Ｄｉｎのナイキスト周波数Ｆｎ以上の周波数成分を与えることが可能になり、出力画像Ｄｏｕｔの解像感を増すことが出来る。
【０１３０】
図１８（ａ）〜（ｄ）は上記の効果を別の観点から説明するための図である。図１８（ａ）はステップエッジ信号を表している。図１８（ａ）に示すエッジではエッジ中央より左側の方が右側より輝度が低くなっている。すなわち、エッジ中央より左側が低レベル側、右側が高レベル側になる。図１８（ｂ）はステップエッジ信号をサンプリング間隔Ｓ２でサンプリングして得られる入力画像Ｄｉｎを表している。
【０１３１】
図１８（ｃ）は図１８（ｂ）に示す入力画像Ｄｉｎに対して得られる拡大画像Ｄ２Ａを表している。拡大画像Ｄ２Ａは、入力画像Ｄｉｎに対して補間演算を行って得られるので、サンプリング間隔はＳ２の半分のＳ１になるが、エッジ近傍の信号の変化はなだらかなままであり、低レベル側の境界は座標Ｐ３で表される画素であり、高レベル側の境界は座標Ｐ４で表される画素のままである。
【０１３２】
図１８（ｄ）はサンプリング間隔Ｓ１でステップエッジ信号をサンプリングした画像(図１８（ｂ）の画像と同じく符号Ｄｉｎで示す)を表している。低レベル側の境界は座標Ｐ１で表される画素であり、高レベル側の境界は座標Ｐ２で表される画素であり、拡大画像Ｄ２Ａと比較し、エッジ近傍での信号の変化が急になっている。
【０１３３】
中間画像Ｄ３２Ａと中間画像Ｄ３２Ｂを加算することで、拡大画像Ｄ２のエッジ近傍での信号の傾きが補正され、図１８（ｄ）に示すステップエッジ信号をサンプリング間隔Ｓ１でサンプリングした画像に近い画像が得られ、画像の解像感を高めることが出来る。
【０１３４】
上記の説明の様に、中間画像Ｄ３２Ａ、Ｄ３２Ｂ（あるいは高周波数成分）を拡大画像Ｄ２Ａに加算することで画像の鮮鋭感を増し、画質を向上することが可能であるが、高周波数成分の加算を過度に行うと画質低下を招くことがある。そこで、本発明の画像処理装置では、高周波数成分画像補正手段３３Ａ、３３Ｂを用いることで、その画質低下を未然に防止することとした。
【０１３５】
図１９（ａ）及び（ｂ）は高周波数成分の加算による画質低下について説明するための図である。
図１９（ａ）は高周波数成分の加算に適度に行うことによって、画像の鮮鋭感を増した場合を、図１９（ｂ）は高周波数成分の加算を過度に行った結果、画質の低下を招いた場合を表す。
【０１３６】
図１９（ａ）は、図１８（ｃ）に示された拡大画像Ｄ２Ａに対して図１７（ｄ）に示された中間画像Ｄ３２Ａ及び図１７（ｆ）に示された中間画像Ｄ３２Ｂを加算した結果を示す図であり、図１８（ｃ）において座標Ｐ３で表されたステップエッジの低輝度側の境界部分が、図１９（ａ）では座標Ｐ１で表される位置へと修正され、図１８（ｃ）において座標Ｐ４で表されたステップエッジの高輝度側の境界部分が、図１９（ａ）では座標Ｐ２で表される位置へと修正され、その結果、図１８（ｃ）と図１９（ａ）を比較すると、図１９（ａ）の方が図１６（ｂ）に示すステップエッジをサンプリング間隔Ｓ１でサンプリングした信号へと近づいていることがわかる。これは高周波数成分の加算に適度に行うことによって、画像の鮮鋭感が増したことを表す。
【０１３７】
一方、図１９（ｂ）も、図１８（ｃ）に示された拡大画像Ｄ２Ａに対して図１７（ｄ）に示された中間画像Ｄ３２Ａ及び図１７（ｆ）に示された中間画像Ｄ３２Ｂを加算した結果を示す図であるが、図１９（ａ）の場合とは異なり、高周波数成分の加算が過度に行われた場合を表している。図１９（ａ）と比較すると座標Ｐ１、Ｐ３で表される位置の輝度がその周辺と比べて不自然に低くなったり（アンダーシュート）、座標Ｐ２、Ｐ４で表される位置の輝度がその周辺と比べて不自然に高くなったり（オーバーシュート）して、画質が低下していることがわかる。
【０１３８】
中間画像Ｄ３２Ａや中間画像Ｄ３２Ｂによって加算あるいは減算される輝度の大きさ（以下、補正量）が必要以上に大きくなると、これら画質の低下要因となるアンダーシュートやオーバーシュートが発生しやすくなる。そこで中間画像Ｄ３２Ａや中間画像Ｄ３２Ｂによる補正量が必要以上に大きくならないよう調整すればよいと考えられる。
【０１３９】
その方法として例えば中間画像Ｄ３２Ａや中間画像Ｄ３２Ｂによって与えられる補正量の局所的な最大値や最小値を検出し、
検出された最大値が所定の値を超えた場合は、中間画像Ｄ３２Ａや中間画像Ｄ３２Ｂによる補正量が小さくなるよう適宜ゲインをかけることで、補正量が必要以上に大きな正の値や必要以上に小さな負の値にならないようにする方法が考えられる。
【０１４０】
ところがこのような方法をとると局所的な最大値を判定するために数画素分のデータを参照しなければならず、回路規模の増加を招く。さらに垂直方向に数画素分のデータを参照しようとするとラインメモリの追加を伴い、コストアップの要因ともなる。
【０１４１】
そこで本発明は、拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ａ若しくは中間画像Ｄ３２Ｂの画素値の符号に基づいて、中間画像Ｄ３２Ａや中間画像Ｄ３２Ｂにかける増幅率を変更することで、中間画像Ｄ３２Ａや中間画像Ｄ３２Ｂによって加算される補正量が必要以上に大きな正の値や必要以上に小さな負の値にならないようにし、オーバーシュートやアンダーシュートの発生を防止することとしたものである。実施の形態１による画像処理装置では、第１の中間画像処理手段３Ｍ及び第２の中間画像処理手段３Ｈにおいて、中間画像Ｄ３２Ａ及び中間画像Ｄ３２Ｂに対して、拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ａあるいは中間画像Ｄ３２Ｂの画素値の符号に応じて適宜異なる増幅率をかけることで、補正量を調節する。
【０１４２】
高周波数成分画像補正手段３３Ａでは図１０（ａ）〜（ｄ）に示すような増幅率Ｄ３ＭＡによって中間画像Ｄ３２Ａが補正され、高周波数成分画像補正手段３３Ｂでは図１２（ａ）〜（ｄ）に示すような増幅率Ｄ３ＨＡによって中間画像Ｄ３２Ｂが補正される。
増幅率Ｄ３ＭＡは中間画像Ｄ３２Ａの画素値の符号が正で、拡大画像Ｄ２Ａの画素値が大きい場合、小さな値をとり、中間画像Ｄ３２Ａの画素値の符号が負で、拡大画像Ｄ２Ａの画素値が小さい場合、小さな値をとる。また、増幅率Ｄ３ＨＡは中間画像Ｄ３２Ｂの画素値の符号が正で、拡大画像Ｄ２Ａの画素値が大きい場合、小さな値をとり、中間画像Ｄ３２Ｂの画素値の符号が負で、拡大画像Ｄ２Ａの画素値が小さい場合、小さな値をとる。
このような増幅率でもって補正することでオーバーシュートやアンダーシュートの発生を抑えることが出来る。理由を以下に述べる。
【０１４３】
中間画像Ｄ３２Ａは拡大画像Ｄ２Ｂに対してハイパスフィルタ処理を行うことで生成される。ここでハイパスフィルタ処理は拡大画像Ｄ２Ｂの各画素値から局所的な平均値を引くことに相当する。従って拡大画像Ｄ２Ｂにおいて画素値が大きな正の値である画素は、その画素に対して与えられるハイパスフィルタ処理後の出力値、すなわち中間画像Ｄ３２Ａの画素値も大きな正の値となる可能性が高い。
【０１４４】
一方、拡大画像Ｄ２Ｂは高周波数成分画像Ｄ１を拡大することで生成される。拡大処理は注目する画素近傍の画素値を加重加算することに相当する。
従って拡大画像Ｄ２Ｂの画素値を算出する際、高周波数成分画像Ｄ１内の画素値が大きな正の値である画素を用いる箇所では、拡大画像Ｄ２Ｂの画素値として大きな正の値が得られる可能性が高い。
【０１４５】
高周波数成分画像Ｄ１は入力画像Ｄｉｎにハイパスフィルタ処理を行うことで生成される。ここでハイパスフィルタ処理は各画素値から局所的な平均値を引くことに相当するので、入力画像Ｄｉｎにおいて画素値が大きな正の値である画素は、高周波数成分画像Ｄ１の画素値も大きな正の値になる可能性が高い。
【０１４６】
以上の話を逆から述べると、入力画像Ｄｉｎにおいて、大きな正の値を持った画素があった場合、高周波数成分画像Ｄ１中で同画素に対応する画素の画素値は大きな正の値になる可能性が高い。また、高周波数成分画像Ｄ１において、画素値が大きな正の値である画素があった場合、同画素値をもちいて算出される拡大画像Ｄ２Ｂの画素値も大きな正の値となる可能性が高いので、入力画像Ｄｉｎにおいて、大きな正の値を持った画素があった場合、拡大画像Ｄ２Ｂ上で同画素に対応する位置にある画素も大きな正の値を持つ可能性が高い。さらに、拡大画像Ｄ２Ｂにおいて、画素値が大きな正の値である画素は、その画素に対して与えられるハイパスフィルタ処理後の出力値、すなわち中間画像Ｄ３２Ａの画素値も大きな正の値となる可能性が高いので、入力画像Ｄｉｎにおいて、大きな正の値を持った画素があった場合、中間画像Ｄ３２Ａ上で同画素に対応する位置にある画素も大きな正の値を持つ可能性が高い。
【０１４７】
一方、拡大画像Ｄ２Ａは入力画像Ｄｉｎに拡大処理を行うことで生成される。拡大処理は注目する画素近傍の画素値を加重加算することに相当する。従って拡大画像Ｄ２Ａの画素値を算出する際、入力画像Ｄｉｎ内の画素値が大きな正の値である画素を用いる箇所では、拡大画像Ｄ２Ａの画素値が大きな正の値になる可能性が高い。また、先に述べた様に、入力画像Ｄｉｎにおいて、大きな正の値を持った画素があった場合、中間画像Ｄ３２Ａ上で同画素に対応する位置にある画素も大きな正の値を持つ可能性が高い。従って入力画像Ｄｉｎにおいて、大きな正の値を持った画素があった場合、中間画像Ｄ３２Ａにおいても、拡大画像Ｄ２Ａにおいても、同画素に対応する位置にある画素は大きな正の値を持つ可能性が高い。
【０１４８】
従って、拡大画像Ｄ２Ａの画素値が大きい場合、中間画像Ｄ３２Ａの画素値も大きな正の値となる可能性が高く、言い換えると中間画像Ｄ３２Ａによって補正量が過度に加算され、オーバーシュートが発生する可能性が高くなる。
【０１４９】
さらに、オーバーシュートは中間画像Ｄ３２Ａの画素値の符号が正の場合発生しうるので、結局、中間画像Ｄ３２Ａの画素値の符号が正でかつ、拡大画像Ｄ２Ａの画素値が大きい場合、中間画像Ｄ３２Ａに対する増幅率を減少してやることで、オーバーシュートが発生しにくい処理が可能となる。
【０１５０】
また、上記と逆の議論を行えば、中間画像Ｄ３２Ａの画素値の符号が負でかつ、拡大画像Ｄ２Ａの画素値が小さい場合、中間画像Ｄ３２Ａに対する増幅率を減少してやることで、アンダーシュートが発生しにくい処理が可能となることがわかる。
【０１５１】
また、中間画像Ｄ３２Ｂは拡大画像Ｄ２Ｂに対し、非線形処理手段２Ａで非線形処理を行った後、高周波数成分画像生成手段３２Ｂにおいてハイパスフィルタ処理を行うことで得られる。非線形処理手段３１ではゼロクロス点の近傍のみ中間画像Ｄ３２Ａを増幅させるので、基本的に拡大画像Ｄ２Ｂが大きな正の値を持っていると非線形処理手段３１の出力する非線形処理画像Ｄ３１も大きな正の値を持っていると考えられる。非線形処理画像Ｄ３１が大きな正の値を持っている場合、非線形処理画像Ｄ３１に対するハイパスフィルタ処理結果である中間画像Ｄ３２Ｂも大きな正の値をもつ可能性が高い。さらに、拡大画像Ｄ２Ａの画素値が大きい場合、中間画像Ｄ３２Ａの画素値も大きな正の値になる可能性が高い。
【０１５２】
さらに、オーバーシュートは中間画像Ｄ３２Ｂの画素値の符号が正の場合発生しうるので、結局、中間画像Ｄ３２Ｂの画素値の符号が正でかつ、拡大画像Ｄ２Ａの画素値が大きい場合、中間画像Ｄ３２Ｂに対する増幅率を減少してやることで、オーバーシュートが発生しにくい処理が可能となる。
【０１５３】
また、上記と逆の議論によって、中間画像Ｄ３２Ｂの画素値の符号が負でかつ、拡大画像Ｄ２Ａの画素値が小さい場合、中間画像Ｄ３２Ｂに対する増幅率を減少してやることで、アンダーシュートが発生しにくい処理が可能となることがわかる。
【０１５４】
上記のような処理は例えば、増幅率Ｄ３ＭＡを図１０（ａ）〜（ｄ）に示すように、中間画像Ｄ３２Ａの画素値の符号が正の時は拡大画像Ｄ２Ａの画素値が大きくなるほど増幅率が減少するような単調減少の関数に基づいて決定し、中間画像Ｄ３２Ａの画素値の符号が負の時は拡大画像Ｄ２Ａの画素値が小さくなるほど増幅率が減少するような単調増加の関数に基づいて決定してやる、あるいは増幅率Ｄ３ＨＡを図１２（ａ）〜（ｄ）に示すように、中間画像Ｄ３２Ｂの画素値の符号が正の時は拡大画像Ｄ２Ａの画素値が大きくなるほど増幅率が減少するような単調減少の関数に基づいて決定し、中間画像Ｄ３２Ｂの画素値の符号が負の時は拡大画像Ｄ２Ａの画素値が小さくなるほど増幅率が減少するような単調増加の関数に基づいて決定してやることで可能となる。
【０１５５】
なお、図１０（ａ）〜（ｄ）、図１２（ａ）〜（ｄ）や式（２）〜（５）に示すような関数はほんの一例であり、拡大画像Ｄ２Ａの画素値及び中間画像Ｄ３２Ａあるいは中間画像Ｄ３２Ｂの画素値の符号に応じて適宜増幅率を決定できる方法であればよい。
【０１５６】
以上のように実施の形態１による画像処理装置では、オーバーシュートやアンダーシュートの発生を抑えつつ、画像の強調処理を行うことができるという効果を奏する。画像中にオーバーシュートやアンダーシュートが過度に発生すると画像の一部のみが異様にちらついたり、暗くなったりして、視覚特性上不愉快に感じられるので、実施の形態１による画像処理装置では、視覚特性上も非常に好ましいものである。
【０１５７】
まとめると中間画像Ｄ３３Ａと画像Ｄ３３Ｂを拡大画像Ｄ２Ａに加算することで過度なオーバーシュートやアンダーシュートの発生を抑えつつ画像の解像感を高めることが可能となる。
【０１５８】
このようにして、実施の形態１による画像処理装置では、過度なオーバーシュートやアンダーシュートの発生を抑えつつ、解像感の高い出力画像Ｄｏｕｔを得ることとしている。
【０１５９】
以上に説明したように、高周波数成分画像生成手段１で生成した高周波数成分画像Ｄ１を画像拡大手段２Ｂで拡大した画像拡大手段２Ｂを、高周波数成分画像処理手段３で処理することによって、入力画像Ｄｉｎのナイキスト周波数Ｆｎ以上の周波数に相当する周波数成分を含んだ高周波数成分画像Ｄ３を得ることができる。そして、加算手段４において、入力画像Ｄｉｎのナイキスト周波数Ｆｎより低い周波数に相当する領域の周波数成分を含む拡大画像Ｄ２Ａと入力画像Ｄｉｎのナイキスト周波数Ｆｎ以上の周波数に相当する領域の周波数成分を含む高周波数成分画像Ｄ３を加算して拡大画像Ｄｏｕｔを生成することとしているので、拡大画像Ｄｏｕｔに対して高周波数成分を十分に与えることができ、解像感のある拡大画像Ｄｏｕｔを得ることができる。
【０１６０】
また、高周波数成分画像生成手段１において、水平方向の高周波数成分を取り出した水平方向高周波数成分画像Ｄ１ｈと垂直方向の高周波数成分を取り出した垂直方向高周波数成分画像Ｄ１ｖを生成することで、画像の水平方向、垂直方向のうちの任意の方向について入力画像Ｄｉｎのナイキスト周波数Ｆｎ以上の周波数に相当した周波数成分を生成することが可能となる。すなわち、水平方向高周波数成分画像Ｄ１ｈを画像拡大手段２Ｂｈで拡大した拡大画像Ｄ２Ｂｈに対し、水平方向高周波数成分画像生成手段３２Ａｈで水平方向のハイパスフィルタをかけることで、水平方向に関して入力画像Ｄｉｎのナイキスト周波数Ｆｎに近い周波数に相当した周波数成分をもった中間画像Ｄ３２Ａｈが生成され、垂直方向高周波数成分画像Ｄ１ｖを画像拡大手段２Ｂｖで拡大した拡大画像Ｄ２Ｂｖに対し、垂直方向高周波数成分画像生成手段３２Ａｖで垂直方向のハイパスフィルタをかけることで、垂直方向に関して入力画像Ｄｉｎのナイキスト周波数Ｆｎに近い周波数に相当した周波数成分をもった中間画像Ｄ３２Ａｖが生成される。
また、水平方向及び垂直方向について異なる特性のハイパスフィルタをかけることで、水平方向と垂直方向とで水平方向と垂直方向とでナイキスト周波数Ｆｎに近い周波数に相当した周波数成分を異なる程度に含むようにすることもできる。
【０１６１】
また、拡大画像Ｄ２Ｂｈに対し水平方向非線形処理手段３１ｈで非線形処理を行って生成した非線形処理画像Ｄ３１ｈに対して、水平方向高周波数成分画像生成手段３２Ｂｈでハイパスフィルタをかけることで、水平方向に関して入力画像Ｄｉｎのナイキスト周波数Ｆｎより高い周波数に相当した周波数成分をもった中間画像Ｄ３２Ｂｈが生成され、拡大画像Ｄ２Ｂｖに対し垂直方向非線形処理手段３１ｖで非線形処理を行って生成した非線形処理画像Ｄ３１ｖに対して、垂直方向高周波数成分画像生成手段３２Ｂｖでハイパスフィルタをかけることで、垂直方向に関して入力画像Ｄｉｎのナイキスト周波数Ｆｎより高い周波数に相当した周波数成分をもった中間画像Ｄ３２Ｂｖが生成される。
また、水平方向及び垂直方向について異なる特性の非線形処理及びハイパスフィルタリングを行なうことで、水平方向と垂直方向とでナイキスト周波数Ｆｎより高い周波数に相当した周波数成分を異なる程度に含むようにすることもできる。
【０１６２】
そして中間画像Ｄ３２Ａについては高周波数成分補正手段３３Ａによって、中間画像Ｄ３２Ｂについては高周波数成分補正手段３３Ｂによって補正をかけることで、過度のオーバーシュートやアンダーシュートが発生するのを防止することができる。
【０１６３】
なお、上記のように高周波数成分補正手段３３Ａ、３３Ｂを過度のオーバーシュートやアンダーシュートの発生を防止するために用いるのは一例であり、
目的に応じて拡大画像Ｄ２Ａと中間画像Ｄ３２Ａ若しくは中間画像Ｄ３２Ｂの画素値の符号と増幅率の関係は適宜設定できる。
たとえば本実施の形態１における画像処理装置のように、入力画像Ｄｉｎからハイパスフィルタを用いて高周波数成分を取り出し、それをもとに生成した成分を加算する場合、ハイパスフィルタの通過帯域を通過するようなノイズが入力画像に含まれていると、入力画像によってはノイズが強調される場合があり、このような現象は特に画像の輝度が低い場合に起こりやすい。従って、例えば図２０に示すように、中間画像Ｄ３２Ａあるいは中間画像Ｄ３２Ｂの画素値の符号が正の場合、拡大画像Ｄ２Ａの画素値ある程度が小さいかある程度大きければ、増幅率が小さくなるよう設定することも可能である。このようにすることで、画像の低輝度部分ではノイズが強調されるのを防止しつつ、過度のオーバーシュートが発生するのを防止することが可能になる。
【０１６４】
ここで図２０において記号ＧｂはＧＭｈｂ、ＧＭｖｂ、ＧＨｈｂ、ＧＨｖｂのいずれかに、記号ｋ１〜ｋ３はｋＭｈ１〜ｋＭｈ３、ｋＭｖ１〜ｋＭｖ３、ｋＨｈ１〜ｋＨｈ３、ｋＨｖ１〜ｋＨｖ３のいずれかに、記号Ａ１、Ａ２はＡＭｈ１、ＡＭｈ２、ＡＭｖ１、ＡＭｖ２、ＡＨｈ１、ＡＨｈ２、ＡＨｖ１、ＡＨｖ２のいずれかに相当する。
【０１６５】
要するに拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ａの画素値の符号に応じて中間画像Ｄ３２Ａ、Ｄ３２Ｂにかける増幅率を適宜変更してやることで、出力画像Ｄｏｕｔに発生しうる様々なアーティファクトを防止することができる。
【０１６６】
また、中間画像Ｄ３２Ａと中間画像Ｄ３２Ｂの両方を上記のような増幅率でもって補正することで最大の効果を得ることが出来るが、中間画像Ｄ３２Ａと中間画像Ｄ３２Ｂのいずれか一方のみを上記のような増幅率でもって補正するだけでも効果は得られる。
【０１６７】
また、上記の実施の形態では、高周波数成分画像処理手段３が、高周波数成分画像補正手段３３Ａ及び３３Ｂの双方を備えているが、一方のみでも良い。例えば、第２の補正成分生成手段３Ｂが、高周波数成分画像補正手段３３Ｂを備えず、非線形処理画像生成手段３０が出力する中間画像Ｄ３２Ｂを第２の補正成分生成手段３Ｂの出力として用いても良く、第１の補正成分生成手段３Ａが、高周波数成分画像補正手段３３Ａを備えず、高周波数成分画像生成手段３２Ａが出力する中間画像Ｄ３２Ａを第１の補正成分生成手段３Ａの出力として用いても良い。
【０１６８】
さらに中間画像Ｄ３３Ａと中間画像Ｄ３３Ｂを加算手段３４で加算して高周波数成分画像Ｄ３を得ることで、任意の方向について入力画像Ｄｉｎのナイキスト周波数Ｆｎ以上の周波数に相当した周波数成分を持った高周波数成分画像Ｄ３を得ることできる。
【０１６９】
なお、入力画像Ｄｉｎから拡大画像Ｄｏｕｔを生成する場合の拡大率を水平方向、垂直方向とも２倍として説明を行ったが拡大率は２倍に限定されるものではない。すなわち、画像拡大手段２Ａにおいて入力画像Ｄｉｎを水平方向、垂直方向ともに所望の倍率に拡大した拡大画像Ｄ２Ａを生成し、高周波数成分生成手段１において入力画像Ｄｉｎをもとに高周波数成分画像Ｄ１を生成し、画像拡大手段２Ｂにおいて、高周波数成分画像Ｄ１を水平方向、垂直方向ともに所望の倍率（画像拡大手段２Ａにおける拡大倍率と同じ倍率）に拡大した拡大画像Ｄ２Ｂを生成し、高周波数成分画像処理手段３において拡大画像Ｄ２Ｂをもとに高周波数成分画像Ｄ３を生成し、加算手段４において拡大画像Ｄ２Ａと高周波数成分画像Ｄ３を加算し、最終的な拡大画像Ｄｏｕｔを得ればよい。
さらに、先にも述べたように、水平方向の拡大率と垂直方向の拡大率とは同じでなくても良く、また水平方向、垂直方向の一方についてのみ拡大を行なっても良い。
【０１７０】
また、上記の説明では水平方向、垂直方向ともゼロクロス点の前後１画素についてのみ増幅率を大きくするとしたが、増幅率の制御の例はこの限りではなく、例えば拡大率に応じて適宜変化させる（拡大率に応じた値に設定する）ことも出来る。
【０１７１】
以下、拡大率が上記の例とは異なる場合について、図２２（ａ）〜（ｅ）及び図２３（ａ）〜（ｆ）を参照して説明する。
【０１７２】
図２２（ａ）にステップエッジ信号、図２２（ｂ）にステップエッジ信号をサンプリング間隔Ｓ１でサンプリングして得られる信号、図２２（ｃ）はステップエッジ信号をサンプリング間隔Ｓ１でサンプリングして得られる信号の高周波数成分を表し、図２２（ｄ）はステップエッジ信号を間隔Ｓ１の３倍の間隔Ｓ３でサンプリングして得られる信号、図２２（ｅ）はステップエッジ信号を間隔Ｓ３でサンプリングして得られる信号の高周波数成分を表す。なお、図２２（ｄ）及び（ｅ）において画素の位置ＰＬ１、ＰＲ１はステップエッジ信号の境界（輝度の明暗が変化する地点）を表す。通常、ステップエッジ信号をサンプリングした画像の高周波数成分を表す信号において、ゼロクロス点Ｚ近傍での局所的な最大値、最小値を与える画素の位置は、ステップエッジ信号の境界の位置とほぼ一致する。
【０１７３】
図２３（ａ）〜（ｆ）は、拡大率が３倍の場合の、高周波数成分画像生成手段１、画像拡大手段２Ｂ、非線形処理手段３１及び高周波数成分画像生成手段３２Ｂによる高周波数成分生成の手順を模式的に表した図であり、図２３（ａ）は輝度、彩度などの成分値がステップ状に変化する画像（ステップ画像）、図２３（ｂ）はステップ画像に対応した入力画像Ｄｉｎ、図２３（ｃ）は高周波数成分画像Ｄ１、図２３（ｄ）は拡大画像Ｄ２Ｂ、図２３（ｅ）は非線形処理画像Ｄ３１、図２３（ｆ）は中間画像Ｄ３２Ｂを表す。なお、説明を簡単にするため各々１次元信号として記載した。
【０１７４】
図２３（ｄ）に示すように、拡大画像Ｄ２Ｂにおいてゼロクロス点Ｚ近傍での局所的な最大値、最小値を与える画素の位置ＰＬ１、ＰＲ１は、拡大画像Ｄ２Ｂにおいてもステップエッジ信号の境界の位置とほぼ一致する。通常、本実施の形態の説明で用いた拡大方法ではこのＰＬ１、ＰＲ１の位置は変化せず、ＰＬ１、ＰＲ１で表す位置とゼロクロス点Ｚの間に存在する画素の数が多くなる。また、ＰＬ１、ＰＲ１で表す位置とゼロクロス点Ｚの間に存在する画素の数は拡大画像Ｄ２Ｂを生成する際の拡大率を大きくすれば（あるいはサンプリング間隔を短くすれば）多くなる。
一方、ステップエッジ信号を短いサンプリング間隔でサンプリングした画像の高周波数成分を表す信号では、ゼロクロス点Ｚ近傍での局所的な最大値、最小値を与える画素の位置はよりゼロクロス点Ｚに近づき、ゼロクロス点にＺより近い画素ほど高周波数成分を表す信号の振幅が大きくなる。
【０１７５】
従って、ゼロクロス点Ｚ前後の信号のみ増幅して非線形処理画像Ｄ３１を生成する際に、ＰＬ１及びＰＲ１よりゼロクロス点Ｚにより近い画素になるほど振幅が大きくなるよう処理してやることが好ましく、例えば位置ＰＬ１、ＰＲ１よりゼロクロス点Ｚに近い画素ではゼロクロス点Ｚにより近い画素ほど大きな増幅率で、ＰＬ１、ＰＲ１よりゼロクロス点Ｚから遠い画素については増幅率１で拡大画像Ｄ２Ｂの画素値を増幅することで、図２３（ｅ）に示すような、ゼロクロス点Ｚにより近い画素ほど大きな振幅をもった非線形処理画像Ｄ３１を生成することができる。
そしてこのようにして生成した拡大画像Ｄ２Ｂからハイパスフィルタ処理によって高周波数成分のみを取り出すことで図２３（ｆ）に示すようなサンプリング間隔Ｓ１に対応した中間画像Ｄ３２Ｂを生成できる。
【０１７６】
以上をまとめると、位置ＰＬ１、ＰＲ１とゼロクロス点Ｚの間に存在する画素の数は拡大画像Ｄ２Ｂ生成時の拡大率によって異なるので、拡大画像Ｄ２Ｂから非線形処理画像Ｄ３１を生成する際にゼロクロス点Ｚ前後において増幅率を１より大きくする画素の数を画像の拡大率に応じて変えてもよい。また、これらの画素に対する増幅率も画素に応じて、例えばゼロクロス点Ｚからの距離に応じて変えてもよい。たとえば、ゼロクロス点Ｚに近い画素ほど増幅率を大きくしてもよい。
【０１７７】
また、拡大画像Ｄ２Ａに中間画像Ｄ３３Ｂを加算するだけでもナイキスト周波数Ｆｎより高い領域の高周波数成分を与えることが出来るので、画像の解像感を増すことは可能である。すなわち、高周波数成分画像処理手段３はその内部に非線形処理を行う非線形処理画像生成手段３０と高周波数成分画像補正手段３３Ｂのみを含む構成であってもよい。
【０１７８】
また、拡大画像Ｄ２Ａに中間画像Ｄ３３Ｂを加算するだけでも、ナイキスト周波数Ｆｎより高い領域の高周波数成分を与え、画像の解像感を増すことが可能である。ナイキスト周波数Ｆｎより高い領域の高周波数成分は非線形処理画像生成手段３０による非線形処理によって生成されるので、高周波数成分画像処理手段３はその内部に非線形処理画像生成手段３０を含んでいればよく、そのような変形例として例えば図３３に示した構成も考えられる。図３３に示した構成での各構成要素のうち加算手段３５以外の動作は図３に記載されたものと同等であるのでその説明は省略する。加算手段３５は水平方向中間画像Ｄ３３Ｂｈと垂直方向中間画像Ｄ３３Ｂｖを加算する。そして加算結果が高周波数成分画像Ｄ３として高周波数成分画像処理手段３から出力される。
【０１７９】
図３３に示した画像処理装置によっても非線形処理画像生成手段３０によって生成されるナイキスト周波数Ｆｎより高い領域の高周波数成分が高周波数成分画像Ｄ３を介して拡大画像Ｄ２Ａに加算されるので、出力画像Ｄｏｕｔの解像感を高めることが出来る。また、高周波数成分画像補正手段３３Ｂによって中間画像Ｄ３２Ｂの画素値が拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ｂの画素値の符号に応じて増幅されるのでオーバーシュート、アンダーシュート等のアーティファクトを抑制できる。
【０１８０】
実施の形態２．
実施の形態１では、本発明をハードウエアにより実現するものとして説明したが、図１に示される構成の一部又は全部をソフトウエアにより、即ちプログラムされたコンピュータにより実現することも可能である。その場合の処理を図２４、並びに図２５〜図３２を参照して説明する。
【０１８１】
図２４は、実施の形態２の画像処理装置を示す。図示の画像処理装置は、ＣＰＵ１１と、プログラムメモリ１２と、データメモリ１３と、これらを接続するバス１４を有する。
ＣＰＵ１１は、プログラムメモリ１２に記憶されたプログラムに従って動作する。動作の過程で種々のデータをデータメモリ１３に記憶させる。処理の結果生成される拡大画像Ｄｏｕｔは、インターフェース１５を介して表示部９に供給され、表示部９による表示に用いられる。
以下、ＣＰＵ１１により行なわれる処理を図２５〜図３２を参照して説明する。
【０１８２】
図２５は、図２４の画像処理装置で実施される画像処理方法のフローを表す図であり、図１６に示される画像処理方法は、画像拡大ステップＳＴ２Ａ、高周波数成分画像生成ステップＳＴ１、画像拡大ステップＳＴ２Ｂ、高周波数成分画像処理ステップＳＴ３、及び加算ステップＳＴ４を有する。
【０１８３】
画像拡大ステップＳＴ２Ａは図示しない画像入力ステップにて入力された入力画像Ｄｉｎを、図１、図３の画像拡大手段２Ａと同様の処理で拡大した拡大画像Ｄ２Ａを生成する。
【０１８４】
高周波数成分画像生成ステップＳＴ１は、図２６に示すように、水平方向高周波数成分画像生成ステップＳＴ１ｈ、及び垂直方向高周波数成分画像生成ステップＳＴ１ｖを有する。水平方向高周波数成分画像生成ステップＳＴ１ｈでは、入力画像Ｄｉｎに対し、図３の水平方向高周波数成分画像生成手段１ｈと同様の処理を行い、水平方向高周波数成分画像Ｄ１ｈを生成する。一方、垂直方向高周波数成分画像生成ステップＳＴ１ｖでは、入力画像Ｄｉｎに対し、図３の垂直方向高周波数成分画像生成手段１ｖと同様の処理を行い、垂直方向高周波数成分画像Ｄ１ｖを生成する。
【０１８５】
画像拡大ステップＳＴ２Ｂは、図２７に示すように、画像拡大ステップＳＴ２Ｂｈ、及び画像拡大ステップＳＴ２Ｂｖを有する。
画像拡大ステップＳＴ２Ｂｈでは、水平方向高周波数成分画像生成ステップＳＴ１ｈで生成した水平方向高周波数成分画像Ｄ１ｈに対し、図３の画像拡大手段２Ｂｈと同様の処理を行い、拡大画像Ｄ２Ｂｈを生成する。
画像拡大ステップＳＴ２Ｂｖでは、垂直方向高周波数成分画像生成ステップＳＴ１ｖで生成した垂直方向高周波数成分画像Ｄ１ｖに対し、図３の画像拡大手段２Ｂｖと同様の処理を行い、拡大画像Ｄ２Ｂｖを生成する。
【０１８６】
次に高周波数成分画像処理ステップＳＴ３の動作を説明する。
高周波数成分画像処理ステップＳＴ３は、図２８に示すように、高周波数成分画像生成ステップＳＴ３２Ａ、非線形処理画像生成ステップＳＴ３０、高周波数成分画像補正ステップＳＴ３３Ａ、高周波数成分画像補正ステップＳＴ３３Ｂ、及び加算ステップＳＴ３４を有する。
高周波数成分通過ステップＳＴ３２Ａは、水平方向高周波数成分画像生成ステップＳＴ３２Ａｈ、垂直方向高周波数成分画像生成ステップＳＴ３２Ａｖを有する。
【０１８７】
非線形処理画像生成ステップＳＴ３０は、非線形処理ステップＳＴ３１、及び高周波数成分画像生成ステップＳＴ３２Ｂを有する。
非線形処理ステップＳＴ３１は水平方向非線形処理ステップＳＴ３１ｈ、及び垂直方向非線形処理ステップＳＴ３１ｖを有する。
高周波数成分画像生成ステップＳＴ３２Ｂは、水平方向高周波数成分画像生成ステップＳＴ３２Ｂｈ、及び垂直方向高周波数成分画像生成ステップＳＴ３２Ｂｖを有する。
【０１８８】
水平方向高周波数成分通過ステップＳＴ３２Ａｈでは、画像拡大ステップＳＴ２Ｂｈで生成した拡大画像Ｄ２Ｂｈに対し、図３の水平方向高周波数成分画像生成手段３２Ａｈと同様の処理を行い、水平方向中間画像Ｄ３２Ａｈを生成する。垂直方向高周波数成分通過ステップＳＴ３２Ａｖでは、画像拡大ステップＳＴ２Ｂｖで生成した拡大画像Ｄ２Ｂｖに対し、図３の垂直方向高周波数成分画像生成手段３２Ａｖと同様の処理を行い、垂直方向中間画像Ｄ３２Ａｖを生成する。
そして高周波数成分画像生成ステップＳＴ３２Ａでは水平方向中間画像Ｄ３２Ａｈと垂直方向中間画像Ｄ３２Ａｖから成る中間画像Ｄ３２Ａが生成される。
このように、高周波数成分通過ステップＳＴ３２Ａでは、図３の高周波数成分画像生成手段３２Ａと同様の動作が行われる。
【０１８９】
水平方向非線形処理ステップＳＴ３１ｈは、図２９に示すようにゼロクロス判定ステップＳＴ３１１ｈと信号増幅ステップＳＴ３１２ｈを有する。
水平方向非線形処理ステップＳＴ３１ｈの動作は以下のごとくである。
まず、ゼロクロス判定ステップＳＴ３１１ｈで、画像拡大ステップＳＴ２Ｂｈで生成した拡大画像Ｄ２Ｂｈにおける画素値の変化を水平方向に沿って確認する。そして画素値が正の値から負の値あるいは負の値から正の値に変化する箇所をゼロクロス点として捉え、ゼロクロス点の左右に位置する画素を特定する。信号増幅ステップＳＴ３１２ｈでは、拡大画像Ｄ２Ｂｈのうち、ゼロクロス判定ステップＳＴ３１１ｈで特定されたゼロクロス点の左右に位置する画素の画素値を増幅し、その結果得られる画像を非線形処理画像Ｄ３１ｈとして生成する。
【０１９０】
垂直方向非線形処理ステップＳＴ３１ｖは、図３０に示すようにゼロクロス判定ステップＳＴ３１１ｖと信号増幅ステップＳＴ３１２ｖを有する。
垂直方向非線形処理ステップＳＴ３１ｖの動作は以下のごとくである。
まず、ゼロクロス判定ステップＳＴ３１１ｖで、画像拡大ステップＳＴ２Ｂｖで生成した拡大画像Ｄ２Ｂｖにおける画素値の変化を垂直方向に沿って確認する。そして画素値が正の値から負の値あるいは負の値から正の値に変化する箇所をゼロクロス点として捉え、ゼロクロス点の上下に位置する画素を特定する。信号増幅ステップＳＴ３１２ｖでは、拡大画像Ｄ２Ｂｖのうち、ゼロクロス判定ステップＳＴ３１１ｖで特定されたゼロクロス点の上下に位置する画素の画素値を増幅し、その結果得られる画像を非線形処理画像Ｄ３１ｖとして生成する。
【０１９１】
そして高周波数成分画像生成ステップＳＴ３１では水平方向非線形処理画像Ｄ３１Ａｈと垂直方向非線形処理画像Ｄ３２Ａｖから成る非線形処理画像Ｄ３１が生成される。
このように、非線形処理ステップＳＴ３１では、図３の非線形処理手段３１と同様の動作が行われる。
【０１９２】
水平方向高周波数成分通過ステップＳＴ３２Ｂｈは、水平方向非線形処理ステップＳＴ３１ｈで生成した非線形処理画像Ｄ３１ｈにハイパスフィルタをかけ、水平方向中間画像Ｄ３２Ｂｈを生成する。垂直方向高周波数成分通過ステップＳＴ３２Ｂｖは、垂直方向非線形処理ステップＳＴ３１ｖで生成した非線形処理画像Ｄ３１ｖにハイパスフィルタをかけ、垂直方向中間画像Ｄ３２Ｂｖを生成する。
そして高周波数成分画像生成ステップＳＴ３２Ｂでは水平方向中間画像Ｄ３２Ｂｈと垂直方向中間画像Ｄ３２Ｂｖから成る中間画像Ｄ３２Ｂが生成される。
このように、高周波数成分通過ステップＳＴ３２Ｂでは、図３の高周波数成分画像生成手段３２Ｂと同様の動作が行われる。
【０１９３】
次に、高周波数成分画像補正ステップＳＴ３３Ａの詳細な動作について説明する。
高周波数成分画像補正ステップＳＴ３３Ａは、図３１に示すように増幅率計算ステップＳＴ３ＭＡと画素値増幅ステップＳＴ３ＭＢを含み、拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ａの画素値の符号に応じて求められる増幅率を使って、中間画像Ｄ３２Ａの各画素値を補正する。なおここで言う補正とは、拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ａの画素値の符号に応じて求められる増幅率を、中間画像Ｄ３２Ａの各画素値に掛け合わせる（あるいは拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ａの画素値の符号に応じて求められる増幅率によって画素値を増幅する）ことを意味する。
【０１９４】
ここで中間画像Ｄ３２Ａは、水平方向中間画像Ｄ３２Ａｈと垂直方向中間画像Ｄ３２Ａｖとから成るので、増幅率の計算は水平方向中間画像Ｄ３２Ａｈと垂直方向中間画像Ｄ３２Ａｖの各々について行われる。すなわち、水平方向中間画像Ｄ３２Ａｈに関しては水平方向増幅率計算ステップＳＴ３ＭＡｈにおいて各画素に対する増幅率が決定され、垂直方向中間画像Ｄ３２Ａｖに関しては垂直方向増幅率計算ステップＳＴ３ＭＡｖにおいて各画素に対する増幅率が決定される。
【０１９５】
ここで水平方向増幅率計算ステップＳＴ３ＭＡｈの動作は水平方向増幅率計算手段３ＭＡｈと、垂直方向増幅率計算ステップＳＴ３ＭＡｖの動作は垂直方向増幅率計算手段３ＭＡｖと同様であるのでその説明は省略する。
【０１９６】
次に画素値増幅ステップＳＴ３ＭＢでは、増幅率計算ステップＳＴ３ＭＡにおいて決定された増幅率に基づいて中間画像Ｄ３２Ａの各画素がもつ画素値を増幅する。ここで中間画像Ｄ３２Ａは水平方向中間画像Ｄ３２Ａｈ及び垂直方向中間画像Ｄ３２Ａｖから成るので、画素値の増幅は、水平方向中間画像Ｄ３２Ａｈと垂直方向中間画像Ｄ３２Ａｖの各々について行われる。
【０１９７】
即ち、水平方向画素値増幅ステップＳＴ３ＭＢｈでは、水平方向増幅率計算ステップＳＴ３ＭＡｈで決定された増幅率に基づいて水平方向中間画像Ｄ３２Ａｈの各画素値を増幅した画像Ｄ３ＭＢｈが生成される。また、垂直方向画素値増幅ステップＳＴ３ＭＢｖでは、垂直方向増幅率計算ステップＳＴ３ＭＡｖで決定された増幅率に基づいて垂直方向中間画像Ｄ３２Ａｖの各画素値を増幅した画像Ｄ３ＭＢｖが生成される。この動作は画素値増幅手段３ＭＢの動作と同じである。
【０１９８】
そして画像Ｄ３ＭＢｈに相当する水平方向中間画像Ｄ３３Ａｈと画像Ｄ３ＭＢｖに相当する垂直方向中間画像Ｄ３３Ａｖから成る中間画像Ｄ３３Ａが、高周波数成分画像補正ステップＳＴ３３Ａによって生成される。
【０１９９】
以上が高周波数成分画像補正ステップＳＴ３３Ａの動作であり、その動作は図９の高周波数成分補正手段３３Ａの動作と同等である。
【０２００】
次に、高周波数成分画像補正ステップＳＴ３３Ｂの詳細な動作について説明する。
高周波数成分画像補正ステップＳＴ３３Ｂは、図３２に示すように増幅率計算ステップＳＴ３ＨＡと画素値増幅ステップＳＴ３ＨＢを含み、拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ｂの画素値の符号に応じて求められる増幅率を使って、中間画像Ｄ３２Ｂの各画素値を補正する。なおここで言う補正とは、拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ｂの画素値の符号に応じて求められる増幅率を、中間画像Ｄ３２Ｂの各画素値に掛け合わせる（あるいは拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ｂの画素値の符号に応じて求められる増幅率によって画素値を増幅する）ことを意味する。
【０２０１】
ここで中間画像Ｄ３２Ｂは、水平方向中間画像Ｄ３２Ｂｈと垂直方向中間画像Ｄ３２Ｂｖとから成るので、増幅率の計算は水平方向中間画像Ｄ３２Ｂｈと垂直方向中間画像Ｄ３２Ｂｖの各々について行われる。すなわち、水平方向中間画像Ｄ３２Ｂｈに関しては水平方向増幅率計算ステップＳＴ３ＨＡｈにおいて各画素に対する増幅率が決定され、垂直方向中間画像Ｄ３２Ｂｖに関しては垂直方向増幅率計算ステップＳＴ３ＨＡｖにおいて各画素に対する増幅率が決定される。ここで水平方向増幅率計算ステップＳＴ３ＨＡｈの動作は水平方向増幅率計算手段３ＨＡｈと、垂直方向増幅率計算ステップＳＴ３ＨＡｖの動作は垂直方向増幅率計算手段３ＨＡｖと同様であるのでその説明は省略する。
【０２０２】
次に画素値増幅ステップＳＴ３ＨＢでは、増幅率計算ステップＳＴ３ＨＡにおいて決定された増幅率に基づいて中間画像Ｄ３２Ｂの各画素がもつ画素値を増幅する。ここで中間画像Ｄ３２Ｂは水平方向中間画像Ｄ３２Ｂｈ及び垂直方向中間画像Ｄ３２Ｂｖから成るので、画素値の増幅は、水平方向中間画像Ｄ３２Ｂｈと垂直方向中間画像Ｄ３２Ｂｖの各々について行われる。
【０２０３】
即ち、水平方向画素値増幅ステップＳＴ３ＨＢｈでは、水平方向増幅率計算ステップＳＴ３ＨＡｈで決定された増幅率に基づいて水平方向中間画像Ｄ３２Ｂｈの各画素値を増幅した画像Ｄ３ＨＢｈが生成される。また、垂直方向画素値増幅ステップＳＴ３ＨＢｖでは、垂直方向増幅率計算ステップＳＴ３ＭＡｖで決定された増幅率に基づいて垂直方向中間画像Ｄ３２Ａｖの各画素値を増幅した画像Ｄ３ＨＢｖが生成される。この動作は画素値増幅手段３ＨＢの動作と同じである。
【０２０４】
そして画像Ｄ３ＨＢｈに相当する水平方向中間画像Ｄ３３Ｂｈと画像Ｄ３ＨＢｖに相当する垂直方向中間画像Ｄ３３Ｂｖから成る中間画像Ｄ３３Ｂが、高周波数成分画像補正ステップＳＴ３３Ｂによって生成される。
【０２０５】
以上が高周波数成分画像補正ステップＳＴ３３Ｂの動作であり、その動作は図１１の高周波数成分補正手段３３Ｂの動作と同等である。
【０２０６】
高周波数成分画像生成ステップＳＴ３２Ａと高周波数成分画像補正ステップＳＴ３３Ａの組み合わせが、第１の補正成分生成手段３Ａの動作に対応し、非線形処理画像生成ステップＳＴ３０と高周波数成分画像補正ステップＳＴ３３Ｂの組み合わせが、第２の補正成分生成手段３Ｂの動作に対応する。
【０２０７】
加算ステップＳＴ３４は、高周波数成分画像補正ステップＳＴ３３Ａで生成した中間画像Ｄ３３Ａと高周波数成分画像補正ステップＳＴ３３Ｂで生成した中間画像Ｄ３３Ｂを加算し、高周波数成分画像Ｄ３を得る。このように、ステップＳＴ３４では、図３の加算手段３４と同様の動作が行われる。
【０２０８】
以上が、高周波数成分画像処理ステップＳＴ３の動作であり、この動作は、図１、図３の高周波数成分画像処理手段３と同じである。
【０２０９】
加算ステップＳＴ４は、画像拡大ステップＳＴ２Ａで生成した拡大画像Ｄ２Ａと高周波数成分画像処理ステップＳＴ３で生成した高周波数成分画像Ｄ３を加算した画像Ｄｏｕｔを生成する。そして生成された画像Ｄｏｕｔが図示しないステップによって、最終的な拡大画像として出力される。
この動作は、図１、図３の加算手段４と同じである。
【０２１０】
以上が実施の形態２による画像処理方法の動作である。上記の説明から明らかなように実施の形態２による画像処理方法でも実施の形態１による画像処理装置と同様の処理で画像を拡大できるため、実施の形態１による画像処理装置と同様の効果が得られる。また、実施の形態２による画像処理方法にも、実施の形態１による画像処理装置と同様の変形を行うことができ、その場合に得られる効果も実施の形態１による画像処理装置と同様である。
【０２１１】
例えば高周波数成分画像処理ステップＳＴ３を図３４に示した構成へと変形する例が考えられる。図３４中の各ステップの動作のうち加算ステップＳＴ３５以外は図２８に記載されたものと同等であるのでその説明は省略する。加算ステップＳＴ３５では水平方向中間画像Ｄ３３Ｂｈと垂直方向中間画像Ｄ３３Ｂｖを加算し、高周波数成分画像Ｄ３を生成する。このようにして得られた高周波数成分画像Ｄ３を拡大画像Ｄ２Ａに加算するだけでも、ナイキスト周波数Ｆｎより高い領域の高周波数成分を与えることが出来るので画像の解像感を増す効果が得られる。非線形処理画像生成ステップＳＴ３０による非線形処理によってナイキスト周波数Ｆｎより高い領域の高周波数成分が生成されるので、解像感を増すためには、高周波数成分画像処理ステップＳＴ３はその内部で非線形処理画像生成ステップＳＴ３０による非線形処理を実行できればよい。また、高周波数成分画像補正ステップＳＴ３３Ｂによって中間画像Ｄ３２Ｂの画素値が拡大画像Ｄ２Ａの画素値と中間画像Ｄ３２Ｂの画素値の符号に応じて増幅されるのでオーバーシュート、アンダーシュート等のアーティファクトを抑制できる。
【０２１２】
なお、変形例は上記のものに留まらず、実施の形態１による画像処理装置に加えられる変形はすべて実施の形態２による画像処理方法に適用可能である。例えば拡大率に応じて非線形処理画像生成ステップＳＴ３０の動作を適宜変化させる（設定する）ことが出来る。その際、実施の形態２による画像処理方法の各構成要素をどのように変形させればよいかは、実施の形態１との対比から明らかである。
【０２１３】
また、実施の形態２による画像処理装置は実施の形態１で説明した画像処理装置と同様画像表示装置の一部として用いることができるため、実施の形態２による画像処理装置で生成された画像Ｄｏｕｔを表示する画像表示装置も、実施の形態１で説明した画像処理装置と同様の効果が得られる。さらに、実施の形態１及び実施の形態２の画像処理装置を用いて実施される画像処理方法、及びこれを用いた画像表示方法も同様の効果が得られる。
【０２１４】
実施の形態３．
図３５は本発明の実施の形態３による画像処理装置の構成を表す図である。実施の形態３による画像処理装置は、画像拡大手段２Ａと、高周波数成分画像生成手段１と、画像拡大手段２Ｂと、高周波数成分画像処理手段３と、高周波数成分画像補正手段３３Ｃと、加算手段４とを備える。また、高周波数成分画像処理手段３は、高周波数成分画像生成手段３２Ａ、非線形処理画像生成手段３０、及び加算手段３４を備える。実施の形態３による画像処理装置も実施の形態１による画像処理装置と同様、図２に示す画像表示装置の一部として利用可能である。
【０２１５】
ここで、画像拡大手段２Ａ、高周波数成分画像生成手段１、画像拡大手段２Ｂ、高周波数成分画像生成手段３２Ａ、及び非線形処理画像生成手段３０の動作は実施の形態１による画像処理装置と同様である。
【０２１６】
加算手段３４は、中間画像Ｄ３２Ａと中間画像Ｄ３２Ｂを加算する。そして得られた画像が高周波数成分画像Ｄ３として出力される。
【０２１７】
高周波数成分画像補正手段３３Ｃは高周波数成分画像Ｄ３の各画素値を、拡大画像Ｄ２Ａの画素値及び補正対象画像（Ｄ３）の画素値の符号に応じて求められる増幅率で増幅し、高周波数成分画像Ｄ３３Ｃを出力する。
【０２１８】
高周波数成分画像補正手段３３Ｃの詳細な動作は、実施の形態１による画像処理装置において、高周波数成分画像補正手段３３Ａが中間画像Ｄ３２Ａに対して行っていた動作あるいは高周波数成分画像補正手段３３Ｂが中間画像Ｄ３２Ｂに対して行っていた動作と同様に、補正対象画像（Ｄ３）の各画素の画素値を、補正対象画像（Ｄ３）の各画素値の符号および拡大画像Ｄ２Ａ（の同じ画素の画素値）に応じて定められる増幅率で増幅するものとすることが出来るのでその説明は省略する。また、高周波数成分画像補正手段３３Ｃの詳細な構成についても同様である。
【０２１９】
加算手段４は、拡大画像Ｄ２Ａと高周波数成分画像Ｄ３３Ｃを加算し、最終的な出力画像Ｄｏｕｔを出力する。
【０２２０】
実施の形態３による画像処理装置でも実施の形態１による画像処理装置と同様、高周波数成分画像処理手段３において、入力画像Ｄｉｎのナイキスト周波数以上の高周波数成分を与えることで出力画像Ｄｏｕｔの解像感を増すことが可能である。
【０２２１】
また、高周波数成分画像補正手段３３Ｃが高周波数成分画像Ｄ３に対して行なう動作は実施の形態１による画像処理装置において、高周波数成分画像補正手段３３Ａが中間画像Ｄ３２Ａに対して行っていた動作あるいは高周波数成分画像補正手段３３Ｂが中間画像Ｄ３２Ｂに対して行っていた動作と同様である。従って、高周波数成分（高周波数成分画像Ｄ３）による補正が過度に行われ、オーバーシュートやアンダーシュート等が発生し、画質が低下するといった問題の発生を防ぐことが出来る。
【０２２２】
さらに、実施の形態１による画像処理装置と比較した場合、高周波数成分画像補正手段の数が少ないので実施の形態１より少ない回路規模で構成可能である。
【０２２３】
また実施の形態３においても、実施の形態１に対する変形と同様の変形を適宜加えることが出来る。図３６にその一例を示す。ナイキスト周波数Ｆｎより高い領域の高周波数成分は非線形処理画像生成手段３０による非線形処理によって生成されるので、高周波数成分画像処理手段３はその内部に非線形処理画像生成手段３０を含んでいればよい。図３６に示した構成での各構成要素のうち加算手段３５以外の動作は実施の形態１あるいは図３５で説明したものと同等であるのでその説明は省略する。加算手段３５は中間画像Ｄ３２Ｂに含まれる水平方向中間画像Ｄ３２Ｂｈと垂直方向中間画像Ｄ３２Ｂｖを加算する。そして加算結果が補正対象画像（Ｄ３）として高周波数成分画像処理手段３から出力される。
【０２２４】
図３６に示した画像処理装置によっても非線形処理画像生成手段３０によって生成されるナイキスト周波数Ｆｎより高い領域の高周波数成分が高周波数成分画像Ｄ３３Ｃを介して拡大画像Ｄ２Ａに加算されるので、出力画像Ｄｏｕｔの解像感を高めることが出来る。また、高周波数成分画像補正手段３３Ｃが補正対象画像（Ｄ３）に対して行なう動作によって、オーバーシュートやアンダーシュート等に起因する画質の低下を防ぐことが出来る。
【０２２５】
なお、変形例は上記のものに留まらず、例えば拡大率に応じて非線形処理画像生成手段３０の動作を適宜変化させる（設定する）ことが出来る。その際、実施の形態３による画像処理装置の各構成要素をどのように変形させればよいかは、実施の形態１との対比から明らかである。
【０２２６】
実施の形態４．
図３７は本発明の実施の形態４による画像処理方法の構成を表す図である。実施の形態４による画像処理方法は、画像拡大ステップＳＴ２Ａと、高周波数成分画像生成ステップＳＴ１と、画像拡大ステップＳＴ２Ｂと、高周波数成分画像処理ステップＳＴ３と、高周波数成分画像補正ステップＳＴ３３Ｃと、加算ステップＳＴ４とを備える。また、高周波数成分画像処理ステップＳＴ３は、高周波数成分画像生成ステップＳＴ３２Ａ、非線形処理画像生成ステップＳＴ３０、及び加算ステップＳＴ３４を含む。実施の形態４による画像処理方法も実施の形態２による画像処理方法と同様、図２４に示す画像表示装置の一部として利用可能である。
【０２２７】
ここで、画像拡大ステップＳＴ２Ａ、高周波数成分画像生成ステップＳＴ１、画像拡大ステップＳＴ２Ｂ、高周波数成分画像生成ステップＳＴ３２Ａ、及び非線形処理画像生成ステップＳＴ３０の動作は実施の形態２による画像処理方法と同様である。
【０２２８】
加算ステップＳＴ３４は、中間画像Ｄ３２Ａと中間画像Ｄ３２Ｂを加算する。そして得られた画像が高周波数成分画像Ｄ３として出力される。
【０２２９】
高周波数成分画像補正ステップＳＴ３３Ｃは高周波数成分画像Ｄ３の各画素値を、拡大画像Ｄ２Ａの画素値及び補正対象画像（Ｄ３）の画素値の符号に応じて求められる増幅率で増幅し、高周波数成分画像Ｄ３３Ｃを出力する。
【０２３０】
なお、高周波数成分画像補正ステップＳＴ３３Ｃの詳細な動作は、実施の形態２による画像処理方法において、高周波数成分画像補正ステップＳＴ３３Ａが中間画像Ｄ３２Ａに対して行っていた動作あるいは高周波数成分画像補正ステップＳＴ３３Ｂが中間画像Ｄ３２Ｂに対して行っていた動作と同様にすることが出来るのでその説明は省略する。また、高周波数成分画像補正ステップＳＴ３３Ｃの詳細な構成についても同様である。
【０２３１】
加算ステップＳＴ４は、拡大画像Ｄ２Ａと高周波数成分画像Ｄ３３Ｃを加算し、最終的な出力画像Ｄｏｕｔを出力する。
【０２３２】
実施の形態４による画像処理方法の動作は実施の形態３による画像処理装置と同様であるので、実施の形態３による画像処理装置と同様の効果を得ることが出来る。
【０２３３】
また実施の形態４においても、実施の形態１に対する変形と同様の変形を適宜加えることが出来る。図３８にその一例を示す。ナイキスト周波数Ｆｎより高い領域の高周波数成分は非線形処理画像生成ステップＳＴ３０による非線形処理によって生成されるので、高周波数成分画像処理ステップＳＴ３はその内部に非線形処理画像生成ステップＳＴ３０を含んでいればよい。図３８に示した構成での各構成要素のうち加算ステップＳＴ３５以外の動作は実施の形態２あるいは図３７で説明したものと同等であるのでその説明は省略する。加算ステップＳＴ３５は中間画像Ｄ３２Ｂに含まれる水平方向中間画像Ｄ３２Ｂｈと垂直方向中間画像Ｄ３２Ｂｖを加算する。そして加算結果が補正対象画像（Ｄ３）として高周波数成分画像補正ステップＳＴ３３Ｃで用いられる。
【０２３４】
図３８に示した画像処理方法によっても非線形処理画像生成ステップＳＴ３０によって生成されるナイキスト周波数Ｆｎより高い領域の高周波数成分が高周波数成分画像Ｄ３３Ｃを介して拡大画像Ｄ２Ａに加算されるので、出力画像Ｄｏｕｔの解像感を高めることが出来る。また、高周波数成分画像補正ステップＳＴ３３Ｃが補正対象画像（Ｄ３）に対して行なう動作によって、オーバーシュートやアンダーシュート等に起因する画質の低下を防ぐことが出来る。
【０２３５】
なお、変形例は上記のものに留まらず、例えば拡大率に応じて非線形処理画像生成ステップＳＴ３０の動作を適宜変化させる（設定する）ことが出来る。その際、実施の形態４による画像処理方法の各構成要素をどのように変形させればよいかは、実施の形態１との対比から明らかである。
【符号の説明】
【０２３６】
１  高周波数成分画像生成手段、  ２Ａ  画像拡大手段、  ２Ｂ  画像拡大手段、  ３  高周波数成分画像処理手段、  ４  加算手段、  ３０  非線形画像生成手段、  ３１  非線計処理手段、  ３２Ａ  高周波数成分画像生成手段、  ３３Ａ  高周波数成分画像補正手段、  ３３Ｂ  高周波数成分画像補正手段、  ３４  加算手段、  Ｄｉｎ  入力画像、  Ｄ１  高周波数成分画像、  Ｄ２Ａ  拡大画像、  Ｄ２Ｂ  拡大画像、  Ｄ３２Ａ  中間画像、  Ｄ３２Ｂ  中間画像、  Ｄ３３Ａ  中間画像、  Ｄ３３Ｂ  中間画像、  Ｄ３  高周波数成分画像、  Ｄｏｕｔ  出力画像。
【図１】
2011048670.tif 000007
【図２】
2011048670.tif 000008
【図３】
2011048670.tif 000009
【図４】
2011048670.tif 000010
【図５】
2011048670.tif 000011
【図６】
2011048670.tif 000012
【図７】
2011048670.tif 000013
【図８】
2011048670.tif 000014
【図９】
2011048670.tif 000015
【図１０】
2011048670.tif 000016
【図１１】
2011048670.tif 000017
【図１２】
2011048670.tif 000018
【図１３】
2011048670.tif 000019
【図１４】
2011048670.tif 000020
【図１５】
2011048670.tif 000021
【図１６】
2011048670.tif 000022
【図１７】
2011048670.tif 000023
【図１８】
2011048670.tif 000024
【図１９】
2011048670.tif 000025
【図２０】
2011048670.tif 000026
【図２１】
2011048670.tif 000027
【図２２】
2011048670.tif 000028
【図２３】
2011048670.tif 000029
【図２４】
2011048670.tif 000030
【図２５】
2011048670.tif 000031
【図２６】
2011048670.tif 000032
【図２７】
2011048670.tif 000033
【図２８】
2011048670.tif 000034
【図２９】
2011048670.tif 000035
【図３０】
2011048670.tif 000036
【図３１】
2011048670.tif 000037
【図３２】
2011048670.tif 000038
【図３３】
2011048670.tif 000039
【図３４】
2011048670.tif 000040
【図３５】
2011048670.tif 000041
【図３６】
2011048670.tif 000042
【図３７】
2011048670.tif 000043
【図３８】
2011048670.tif 000044