JP6920942B2 - 階調補完装置及びプログラム - Google Patents

階調補完装置及びプログラム Download PDF

Info

Publication number
JP6920942B2
JP6920942B2 JP2017182496A JP2017182496A JP6920942B2 JP 6920942 B2 JP6920942 B2 JP 6920942B2 JP 2017182496 A JP2017182496 A JP 2017182496A JP 2017182496 A JP2017182496 A JP 2017182496A JP 6920942 B2 JP6920942 B2 JP 6920942B2
Authority
JP
Japan
Prior art keywords
image
frequency
gradation
frequency component
band
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2017182496A
Other languages
English (en)
Other versions
JP2019057225A (ja
Inventor
康孝 松尾
康孝 松尾
菊文 神田
菊文 神田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Japan Broadcasting Corp
Original Assignee
Japan Broadcasting Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Japan Broadcasting Corp filed Critical Japan Broadcasting Corp
Priority to JP2017182496A priority Critical patent/JP6920942B2/ja
Publication of JP2019057225A publication Critical patent/JP2019057225A/ja
Application granted granted Critical
Publication of JP6920942B2 publication Critical patent/JP6920942B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Processing (AREA)
  • Picture Signal Circuits (AREA)
  • Compression Or Coding Systems Of Tv Signals (AREA)
  • Compression Of Band Width Or Redundancy In Fax (AREA)

Description

本発明は、入力画像の階調を補完した階調補完画像を生成する階調補完装置及びプログラムに関する。
一般に、JPEGやH.265/HEVC(High Efficiency Video Coding)などの符号化方式では、直交変換処理及び量子化処理が行われている。直交変換処理では、画像信号をブロック領域毎にDCT(離散コサイン変換)して、周波数係数に変換する。そして、量子化処理では、量子化マトリクスを用いて周波数係数を量子化する。このような処理により、視覚的に目立たない高周波成分に対しては粗い量子化を施し、信号成分上で重要な低周波成分に対しては細かい量子化を施すことができ、視覚特性を考慮した圧縮処理を行うことができる(例えば、非特許文献1参照)。
また、超解像画像を生成する技術として、例えば特許文献1には、被超解像画像を含む複数のフレーム画像を周波数分解し、被超解像画像と複数の周波数分解画像の低周波成分との間でブロックマッチングを行い、周波数分解画像の低周波成分と同じ空間位相位置の高周波成分を、被超解像画像の標本化周波数を超える高周波成分に割り付けて、被超解像画像の高周波成分を推定した超解像高周波成分を生成し、低周波成分と高周波成分の周波数再構成を行うことで、超解像画像を生成する技術が開示されている。
また、特許文献2には、原画像を色成分ごとに内挿補間処理して色成分画像を生成し、色成分画像を用いて周波数分解して周波数分解画像を生成し、色成分ごとにブロックマッチングを行って色成分画像の超解像高周波成分を生成し、色成分ごとに周波数再構成を行うことにより高精細な超解像画像を生成する技術が開示されている。
さらに、特許文献3には、原画像又は原画像の変換画像を周波数分解し、ブロックマッチングを行って周波数分解画像の割り付けを行い、周波数再構成することで単一フレームの原画像から高精細な超解像画像を生成する技術が開示されている。
特開2016−134075号公報 特開2015−212905号公報 特開2015−203952号公報
大久保榮監修、「インプレス標準教科書シリーズ H.265/HEVC教科書」、株式会社インプレスジャパン、2013年10月21日
しかしながら、上述した符号化方式において、量子化は非可逆処理であるため、失われた階調値を完全に元に戻すことはできない。そのため、原画像を非可逆符号化処理して符号を生成すると、該符号を復号した際に、復号画像の画質が原画像に対して劣化してしまうという課題が残る。また、上述した超解像画像を生成する技術を用いることにより画像を高精細化できるが、劣化した画像の階調値を補完するという観点では、さらに改善の余地がある。
かかる事情に鑑みてなされた本発明の目的は、符号化処理などにより劣化した画像の階調値を高精度に補完することが可能な階調補完装置及びプログラムを提供することにある。
上記課題を解決するため、本発明に係る階調補完装置は、入力画像の階調を補完した階調補完画像を生成する階調補完装置であって、前記入力画像に対して周波数分解処理を行って複数の低周波成分画像及び複数の高周波成分画像からなる周波数分解画像を生成する周波数分解部と、前記周波数分解画像の位置合わせに用いる閾値を、前記複数の低周波成分画像における所定の低周波成分画像である参照フレームの周波数帯域が低いほど大きな値となるように設定する帯域毎閾値設定部と、前記複数の低周波成分画像間において、基準フレーム内のブロック領域と、該ブロック領域との差分値が前記閾値よりも小さい参照フレーム内のブロック領域との位置関係を示す位置合わせ情報を生成する位置合わせ部と、前記位置合わせ情報に従って前記高周波成分画像に割り付けを行い、高周波成分割付画像を生成する割付部と、前記低周波成分画像及び前記高周波成分割付画像に対して周波数再構成処理を行って階調補完画像を生成する周波数再構成部と、を備えることを特徴とする。
さらに、本発明に係る階調補完装置において、前記入力画像は、符号を復号した復号画像であり、前記帯域毎閾値設定部は、前記符号を生成する際に用いられた量子化パラメータが大きいほど、前記閾値を小さな値となるように設定すること特徴とする。
また、上記課題を解決するため、本発明に係るプログラムは、コンピュータを、上記階調補完装置として機能させることを特徴とする。
本発明によれば、符号化処理などにより劣化した画像の階調値を高精度に補完することが可能となる。
本発明の一実施形態に係る階調補完装置の構成例を示すブロック図である。 本発明の一実施形態に係る階調補完装置における位置合わせ処理の概要を示す図である。 本発明の一実施形態に係る階調補完装置における割り付け処理の概要を示す図である。 本発明の一実施形態に係る階調補完装置における周波数再構成処理の概要を示す図である。
以下、本発明の実施形態について、図面を参照して詳細に説明する。
図1に、本発明の一実施形態に係る階調補完装置の構成例を示す。図1に示す階調補完装置1は、周波数分解部11と、帯域毎閾値設定部12と、位置合わせ部13と、割付部14と、周波数再構成部15とを備える。
階調補完装置1は、復号画像などの入力画像に対して階調補完処理を行い、入力画像の階調が補完された階調補完画像を生成する。
周波数分解部11は、入力画像に対して複数階層の周波数分解(多重解像度分解)処理を行って周波数分解画像を生成する。代表的な周波数分解処理として、ウェーブレット分解処理が挙げられる。N階周波数分解処理を行うと、周波数分解階数n=1,2,…,Nの周波数分解画像が生成される。nの値が小さいほど、周波数帯域が高いことを意味する。また、同一の周波数分解階数nにおいて、さらに、低周波帯域成分の画像LL、水平方向に高周波成分となる画像LH、垂直方向に高周波成分となる画像HL、対角方向に高周波成分となる画像HHに分解される。以下、LLを低周波成分画像と称し、LH,HL,HHをまとめて高周波成分画像と称する。すなわち、周波数分解画像は、複数の低周波成分画像LL、及び複数の高周波成分画像LH,HL,HHからなる。周波数分解部11は、低周波成分画像LLを位置合わせ部13及び周波数再構成部15に出力し、高周波成分画像LH,HL,HHを割付部14に出力する。ここで、周波数分解階数Nはあらかじめ決めておいてもよいし、外部から取得してもよい。
帯域毎閾値設定部12は、周波数分解階数Nを取得する。そして、周波数分解部11により周波数分解(多重解像度分解)処理された画像の位置合わせに用いる帯域毎閾値を、後述する図2における参照フレームの周波数帯域毎(すなわち、周波数分解階数毎)に、周波数帯域が低いほど大きな値となるように設定し、設定した帯域毎閾値を位置合わせ部13に出力する。周波数帯域が低いほど大きな値となるように帯域毎閾値を設定することは、後述する位置合わせ部13において、周波数帯域が低いほど位置合わせ結果が採用されやすくなることを意味する。
本実施形態では、位置合わせ部13において、位置合わせをブロックマッチングで行い、その1画素当たりの差分値をSAD(Sam of Absolute Difference;絶対値差分和)で求めるものとする。表1に、復号画像の階調数が8bitであり、参照フレームの周波数分解階数が2〜4であるときのSADの帯域毎閾値の例を示す。なお、ブロックマッチングで1画素当たりの差分値をSSD(Sum of Squared Difference;二乗誤差和)で求める場合には、表1で示した帯域毎閾値を2乗すればよい。
Figure 0006920942
また、帯域毎閾値設定部12は、入力画像が符号を復号した復号画像である場合には、周波数分解階数Nに加えて、符号を生成する際に用いられた量子化パラメータQP(Quantization Parameter)を取得してもよい。符号化処理における量子化では、ブロック領域毎の周波数係数を、量子化マトリクスに量子化ステップを乗算した値で除算する。量子化パラメータQPの設定値が大きいほど量子化ステップは大きくなるので、圧縮率は高くなり、すなわち階調値の削減率も高くなる。
帯域毎閾値設定部12は、量子化パラメータQPを取得した場合には、量子化パラメータQPが大きいほど、帯域毎閾値を小さな値となるように設定する。特に、高周波帯域の帯域毎閾値を小さな値とする。表2に、量子化パラメータQPを取得した場合について、復号画像の階調数が8bitであり、参照フレームの周波数分解階数が2〜4であるときの、帯域毎閾値の例を示す。表2では、周波数分解階数が2,3であるとき、量子化パラメータQPが大きいほど、帯域毎閾値が小さな値となっている。なお、表2では量子化パラメータQPを40未満と40以上の2種類に分けたが、3種類以上に分けて帯域毎閾値を規定してもよい。
Figure 0006920942
位置合わせ部13は、帯域毎閾値設定部12から帯域毎閾値を入力する。そして、周波数分解部11により生成された周波数分解画像のうちの低周波成分画像間において、基準フレーム内のブロック領域と、該ブロック領域との差分値が帯域毎閾値よりも小さい参照フレーム内のブロック領域との位置関係を示す位置合わせ情報を生成し、割付部14に出力する。すなわち、位置合わせ情報は、低周波成分画像間において、相似するオブジェクトの位置関係を示すものである。ここで、差分値とは、差の絶対値の和、差の絶対値の積、差の二乗和など、種々の演算方法によって定まる差の値の総称である。
図2に、位置合わせ部13における位置合わせ処理の概要を示す。位置合わせ部13は、周波数分解画像のi階低周波成分画像LLを基準フレームとし、j階(i<j)低周波成分画像LLを参照フレームとする。図2に示す例では、i=1、j=2である。そして、基準フレームをブロック領域に分割し、任意のブロックから参照フレームへブロックマッチングを行い、探索範囲内で差分値が帯域毎閾値よりも小さいブロック、すなわち類似度(相関性)の高いブロックの位置関係を示す位置合わせ情報を生成する。差分値が帯域毎閾値より小さければよいので、類似度の高いブロックを複数(図2では3つ)検出してもよい。参照フレームの周波数分解階数jが大きいほど帯域毎閾値は大きくなるため、低周波帯域の画像ほど、基準フレームのブロックに類似するブロックとして検出されやすくなる。
なお、ブロックマッチングは、SAD,SSDなどの評価関数を用いて、既知の手法により行われる。また、ブロックマッチングは、例えば特許文献3等に開示された既知の手法である式(1)に示すパラボラフィッティング関数を用いた補間処理により、小数画素精度で行う。
Figure 0006920942
ここで、探索位置における画素位置をxとしたとき、SSD(x)は、画素位置におけるSSD値を表し、より具体的には、SSD(0)は中心位置(SSD値を最小とする位置)におけるSSD値、SSD(−1)は中心位置から−x方向又は−y方向の隣接画素の位置におけるSSD値、SSD(1)は中心位置から+x方向又は+y方向の隣接画素の位置におけるSSD値を表す。式(1)から、水平又は垂直方向の小数画素精度の画素位置(小数画素位置)をそれぞれ算出することができる。
割付部14は、位置合わせ部13により生成された位置合わせ情報に従って、周波数分解部11により生成された周波数分解画像のうちの高周波成分画像LH,HL,HHに割り付けを行い、割付後の高周波成分画像である高周波成分割付画像LH’,HL’,HH’を生成し、周波数再構成部15に出力する。
図3に、割付部14による割り付け処理の概要を示す。図3に示すように、位置合わせ情報に従って、参照フレームLLと同一階層の高周波成分画像LH,HL,HH内のブロックを、基準フレームLLと同一階層の高周波成分画像LH,HL,HHに割り付ける。図3に示す例では、i=1、j=2である。なお、図2では基準フレーム内のブロック領域に類似する参照フレーム内のブロック領域を3つ示しているが、図3では見やすくするために、基準フレーム内のブロック領域に類似する参照フレーム内のブロック領域を1つとしている。
ここで、高周波成分画像LH,HL,HHに割り付ける際には、参照フレームLL内の同じ位相位置の位置合わせ情報に従うこととする。これは、基準フレームLL内のブロックPが参照フレームLL内のブロックQに類似していれば、高周波成分画像LH,HL,HH内における、ブロックPと同じ位相位置のブロックは、高周波成分LH,HL,HH内における、ブロックQと同じ位相位置のブロックとそれぞれ類似する可能性が高いためである。
また、位置合わせ部13において、類似するブロックの位置関係を小数画素精度で求めた場合には、割付部14は、小数画素位置を通常の画素位置に合わせるために、割り付け後の高周波成分画像に対して、光学系の解像度劣化過程を模擬した点広がり関数(Point spread function;PSF)を用いた補間を行う。式(2)に、点広がり関数を示す。ここで、wはガウス関数の半値幅(分散値)である。
Figure 0006920942
割付部14は、割り付けを全て終えた後、水平、垂直、対角高周波成分として候補が複数存在する場合には、それらの値を平均するか、最大事後確率(Maximum a posteriori;MAP)再構成を行い、未知の値を推定し、最終的な高周波成分割付画像LH’,HL’,HH’として周波数再構成部15に出力する。MAP再構成の詳細については、例えば、E. Levitan and G. Herman: “A maximum a posteriori probability expectation maximization algorithm for image reconstruction in emission tomography”, IEEE Transactions on Medical Imaging, vol. 6, no. 3, pp. 185-192, Sep. 1987.を参照されたい。また、その他の方法として、ML法や、割り付けられた画素の距離に応じた重み付けにより、高周波成分画像を推定してもよい。
周波数再構成部15は、周波数分解部11により生成された低周波成分画像を低周波成分とし、割付部14により生成された高周波成分割付画像を高周波成分として、周波数再構成処理を行って階調補完画像を生成し、外部に出力する。
図4に、周波数再構成部15による周波数再構成処理の概要を示す。図4に示す例では、低周波成分画像LLを低周波成分とし、割付後の水平高周波成分画像LH’、垂直高周波成分画像HL’、対角高周波成分画像HH’を高周波成分として周波数再構成処理を行うことにより、階調補完画像が生成される。なお、周波数分解部11が周波数分解処理としてウェーブレット分解処理を行った場合には、周波数再構成部15は、同じウェーブレットを用いてウェーブレット再構成処理を行う。
なお、上述した階調補完装置1として機能させるためにコンピュータを好適に用いることができ、そのようなコンピュータは、階調補完装置1の各機能を実現する処理内容を記述したプログラムを該コンピュータの記憶部に格納しておき、該コンピュータのCPUによってこのプログラムを読み出して実行させることで実現することができる。なお、このプログラムは、コンピュータ読取り可能な記録媒体に記録可能である。
上述したように、本発明では、符号化処理を経て復号した画像では、高周波帯域の階調値の削減度合いが大きいことに着目し、位置合わせに用いる閾値を、参照フレームの周波数帯域が低いほど大きな値とし、低周波成分画像間において、差分値が該閾値よりも小さいブロック領域の位置関係を示す位置合わせ情報を生成する。さらに、量子化パラメータQPが大きいほど階調値の削減度合いが大きいことに着目し、量子化パラメータが大きいほど、閾値を小さな値としてもよい。
このため、本発明では、参照フレームうち、階調値の削減度合が小さいと考えられる低周波帯域ほど、ブロックマッチングで類似するブロックが検出されやすくなる。したがって、位置合わせの精度を向上させ、削減された階調値を高精度に補完することができるようになる。
上述の実施形態は代表的な例として説明したが、本発明の趣旨及び範囲内で、多くの変更及び置換ができることは当業者に明らかである。したがって、本発明は、上述の実施形態によって制限するものと解するべきではなく、特許請求の範囲から逸脱することなく、種々の変形や変更が可能である。例えば、実施形態の構成図に記載の複数の構成ブロックを1つに組み合わせたり、あるいは1つの構成ブロックを分割したりすることが可能である。
また、プログラムは、コンピュータ読取り可能媒体に記録されていてもよい。コンピュータ読取り可能媒体を用いれば、コンピュータにインストールすることが可能である。ここで、プログラムが記録されたコンピュータ読取り可能媒体は、非一過性の記録媒体であってもよい。非一過性の記録媒体は、特に限定されるものではないが、例えば、CD−ROMやDVD−ROMなどの記録媒体であってもよい。
1 階調補完装置
11 周波数分解部
12 帯域毎閾値設定部
13 位置合わせ部
14 割付部
15 周波数再構成部

Claims (3)

  1. 入力画像の階調を補完した階調補完画像を生成する階調補完装置であって、
    前記入力画像に対して周波数分解処理を行って複数の低周波成分画像及び複数の高周波成分画像からなる周波数分解画像を生成する周波数分解部と、
    前記周波数分解画像の位置合わせに用いる閾値を、前記複数の低周波成分画像における所定の低周波成分画像である参照フレームの周波数帯域が低いほど大きな値となるように設定する帯域毎閾値設定部と、
    前記複数の低周波成分画像間において、基準フレーム内のブロック領域と、該ブロック領域との差分値が前記閾値よりも小さい参照フレーム内のブロック領域との位置関係を示す位置合わせ情報を生成する位置合わせ部と、
    前記位置合わせ情報に従って前記高周波成分画像に割り付けを行い、高周波成分割付画像を生成する割付部と、
    前記低周波成分画像及び前記高周波成分割付画像に対して周波数再構成処理を行って階調補完画像を生成する周波数再構成部と、
    を備えることを特徴とする階調補完装置。
  2. 前記入力画像は、符号を復号した復号画像であり、
    前記帯域毎閾値設定部は、前記符号を生成する際に用いられた量子化パラメータが大きいほど、前記閾値を小さな値となるように設定すること特徴とする、請求項1に記載の階調補完装置。
  3. コンピュータを、請求項1又は2に記載の階調補完装置として機能させるためのプログラム。
JP2017182496A 2017-09-22 2017-09-22 階調補完装置及びプログラム Active JP6920942B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2017182496A JP6920942B2 (ja) 2017-09-22 2017-09-22 階調補完装置及びプログラム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2017182496A JP6920942B2 (ja) 2017-09-22 2017-09-22 階調補完装置及びプログラム

Publications (2)

Publication Number Publication Date
JP2019057225A JP2019057225A (ja) 2019-04-11
JP6920942B2 true JP6920942B2 (ja) 2021-08-18

Family

ID=66107492

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017182496A Active JP6920942B2 (ja) 2017-09-22 2017-09-22 階調補完装置及びプログラム

Country Status (1)

Country Link
JP (1) JP6920942B2 (ja)

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010218271A (ja) * 2009-03-17 2010-09-30 Tokyo Institute Of Technology パラメータ制御処理装置及び画像処理装置
JP5498970B2 (ja) * 2011-01-20 2014-05-21 日本放送協会 画像縮小装置、画像拡大装置、及びこれらのプログラム
JP5514132B2 (ja) * 2011-02-07 2014-06-04 日本放送協会 画像縮小装置、画像拡大装置、及びこれらのプログラム
JP5809920B2 (ja) * 2011-10-14 2015-11-11 東芝デジタルメディアエンジニアリング株式会社 画像処理方法及び装置
JP2015203952A (ja) * 2014-04-14 2015-11-16 日本放送協会 超解像装置及びプログラム
JP6462372B2 (ja) * 2015-01-21 2019-01-30 日本放送協会 超解像装置及びプログラム

Also Published As

Publication number Publication date
JP2019057225A (ja) 2019-04-11

Similar Documents

Publication Publication Date Title
US11272181B2 (en) Decomposition of residual data during signal encoding, decoding and reconstruction in a tiered hierarchy
US10425649B2 (en) Method and apparatus for performing graph-based prediction using optimization function
JP4920599B2 (ja) ハイブリッドビデオ圧縮の量子化雑音除去用の非線形ループ内デノイズフィルタ
US11962805B2 (en) Data processing apparatuses, methods, computer programs and computer-readable media
JP5727873B2 (ja) 動きベクトル検出装置、符号化装置、及びこれらのプログラム
JP2015203952A (ja) 超解像装置及びプログラム
Song et al. Novel near-lossless compression algorithm for medical sequence images with adaptive block-based spatial prediction
Siddeq et al. A novel 2D image compression algorithm based on two levels DWT and DCT transforms with enhanced minimize-matrix-size algorithm for high resolution structured light 3D surface reconstruction
JP7002213B2 (ja) 空間・階調超解像装置及びプログラム
JP2020028067A (ja) 画像符号化装置、画像復号装置、及びプログラム
Hou et al. Graph-based transform for data decorrelation
JP7399646B2 (ja) データ圧縮装置およびデータ圧縮方法
JP6920942B2 (ja) 階調補完装置及びプログラム
Sumalatha et al. Hierarchical lossless image compression for telemedicine applications
JP2010068084A (ja) 画像処理装置及び画像処理方法
CN107948644B (zh) 一种水下图像压缩方法及传输方法
JP6564315B2 (ja) 符号化装置、復号装置、及びプログラム
US20200311875A1 (en) Image processing device, image processing method, and recording medium
JP6864557B2 (ja) 符号化装置、復号装置、及びプログラム
JP6823492B2 (ja) 超解像補完ポストフィルタ装置及びプログラム
JP6817784B2 (ja) 超解像装置及びプログラム
JP6951163B2 (ja) 階調超解像装置及びプログラム
JP6388476B2 (ja) 符号化装置及びプログラム
Asati et al. Similarity based optimization to fractal image encoding based on multithreading parallelization
Islam et al. Optimal combining fusion on degraded compressed sensing image reconstruction

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200728

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20210629

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210727

R150 Certificate of patent or registration of utility model

Ref document number: 6920942

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250