JP2020021168A - 不快度推定装置及び不快度推定プログラム - Google Patents

Abstract

【課題】広視野観視条件を含む様々な異なる観視距離において、様々に異なる動揺が全画面中に分布している動揺映像であっても、観視した視聴者が感知する不快度を高精度に推定できるようにする。【解決手段】不快度推定装置２は、低周波回転動揺補正量を算出する低周波回転動揺補正量算出部３０と、領域別誘目度の違い分を補正する誘目度補正量を算出する誘目度補正量算出部５０と、目立つ動き領域間の距離の違い分を観視距離を加味して補正する動き領域間距離補正量を算出する動き領域間距離補正量算出部７０と、低周波回転動揺補正量、誘目度補正量及び動き領域間距離補正量を用いて動揺認知量推定値を補正する動揺認知量補正部８０と、長時間継続する動揺に対する不快度の蓄積効果分を補正する総動揺エネルギー補正部９０と、を備える。【選択図】図９

Description

本発明は、揺れのある映像に対する動揺の不快度を推定する不快度推定装置及び不快度推定プログラムに関する。

家庭で放送、録画映像、パーソナルコンピュータ映像等を視聴する際のディスプレイが大画面化するのに伴い、画面動揺（映像の揺れ）によって誘発される不快感が増大する傾向があり、場合によっては視聴者が映像酔いを起こして健康被害に至る例もある。画面動揺の多くは、撮影時のビデオカメラの動き、すなわちグローバルモーションに起因する。かかる事態を予防するため、手持ち撮影用ビデオカメラ等には「手ぶれ補正」機能の搭載がほぼ必須となっている。また、撮影された映像に含まれるグローバルモーションを軽減する「映像スタビライズ」ソフトウェアも数多く開発されている。さらに、グローバルモーションと脈拍・血圧等との生理的な指標との関係を探ることによって、不快な状態を推定しようとする試みも行われている（非特許文献１参照）。

しかし、前記した「手ぶれ補正」では、長周期の揺れは補正ができない、又は、不十分である。また、「映像スタビライズ」では、短い期間内の映像の揺れを無くすか低減したり、滑らかな画面移動に変換したりするだけある。すなわち、「映像スタビライズ」は、どのような性質でどの程度の大きさの揺れならば残っていても良いといった客観的な基準に則っているものではない。また、この客観的な基準として脈拍・血圧等の生理的な指標を用いることができたとしても、明示できるのは健康被害に達する条件、及び、その条件に対する現在の生理的状態の程度のみである。

これに対し、放送、映画等の映像コンテンツの供給者は、映像制作段階において、映像の安全性のみならず、視聴者の大半にグローバルモーション（画面全体の動き）あるいは画面のある程度まとまった領域の動きに伴う不快感を催させないように留意する必要がある。しかしながら、前記した従来の技術は、画面の全体あるいはある程度まとまった領域の動きに対する視聴者の心理的な「不快度」を具体的に把握するものではないため、「不快度」を基準とした映像の良否を判定するには不向きであった。

かかる問題を解消するため、本願出願人は、生理的指標を用いることなく、画面の全体又はある程度まとまった領域の動きと主観評価実験で得られた不快度との関係を用いて、映像を解析して得られた物理的特徴量に基づいて画面動揺に対する不快度を推定する不快度推定装置（特許文献１参照）を開発した。

「映像酔いガイドライン検証システムの実用化に関するフィージビリティスタディ報告書」、財団法人機械システム振興協会、平成２０年３月

特開２０１６−１６２２１９号公報

さらに、本願出願人は、「同じ大きさの揺れであっても映像の内容によって不快度が異なる」と言う事象に対応するため、認知される動揺量（動揺認知量）をまず高精度に推定した上で、推定された動揺認知量と映像内容とによって不快度を推定するシステムにおける動揺認知量推定装置（特願２０１７−２３３７８、本願出願時未公開）及び不快度推定装置（特願２０１７−２３６０３、本願出願時未公開）を開発した。

ここで、前記動揺認知量推定装置及び不快度推定装置は、ハイビジョンの標準観視条件である「画面高の３倍の観視距離」で観視した際の動揺認知量及び不快度を推定するものとなっており、それ以外の観視距離で観視した場合の動揺認知量及び不快度については、「同じ動揺映像を画面高の２倍〜８倍の範囲で観視した際の動揺認知量及び不快度は、観視距離が半分になるごとにそれぞれ１ランクずつ上昇する（観視距離が２倍になるごとに１ランクずつ低下する）」という知見に基づき、画面高の３倍の観視距離（３Ｈ）で推定した動揺認知量及び不快度を補正して推定していた。

しかし、近年急激に家庭への普及が進んでいる４Ｋディスプレイや、２０１８年１２月に実用放送開始予定の８Ｋ放送における標準観視条件の観視距離は、４Ｋでは画面高の１．５倍（１．５Ｈ）、８Ｋでは画面高の０．７５倍（０．７５Ｈ）となっている。このような広視野観視条件での動揺映像に関する直近の調査では、前記の「動揺認知量及び不快度は、観視距離が半分になるごとにそれぞれ１ランクずつ上昇する」という性質が成立せず、観視距離が近くなるにつれて上昇する動揺認知量及び不快度が飽和傾向を示すことが明らかになった。さらに不快度に関しては、画面の四隅で異なる動揺を示すような場合には、観視距離が近いほど不快度が低下する割合が大きくなることも明らかになった。

したがって、前記動揺認知量推定装置及び不快度推定装置は、広視野観視条件における動揺認知量及び不快度を正確に推定することができない、という問題点を有していた。このうち動揺認知量推定装置については、動揺認知量が不快度とは異なり、観視距離が変わっても、画面の四隅で異なる動揺を示す場合と、同じ動揺が画面中央付近にまとまっている場合と、で同じ値を示すことから、観視距離が近くなるにつれて上昇する動揺認知量の飽和現象分のみを補正した、改良版の動揺認知量推定装置（特願２０１８−０１３４６０、本願出願時未公開）を既に開発している。ただし不快度推定装置については、画面の四隅で異なる動揺を示す場合、観視距離が近いほど不快度が低下する割合が大きくなることによる誤差発生に関して未解決であった。

本発明は、前記した事情に鑑みて創案されたものであり、８Ｋ映像の広視野観視条件を含む様々に異なる観視距離において、様々に異なる動揺が全画面中に分布している動揺映像を観視した際に生じる不快度を高精度に推定することが可能な不快度推定装置及び不快度推定プログラムを提供することを課題とする。

前記課題を解決するため、本発明の不快度推定装置は、映像における動揺認知量と物理的特徴量に基づいて画面動揺に対する不快度を推定する不快度推定装置であって、回転動揺成分抽出部と、デジタルフィルタ部と、低周波回転動揺補正量算出部と、領域別誘目度算出部と、誘目度補正量算出部と、動き領域間距離正規化分散値算出部と、動き領域間距離補正量算出部と、動揺認知量補正部と、総動揺エネルギー補正部と、を備える。

不快度推定装置は、回転動揺成分抽出部によって、画面を複数に分割した領域に基づき、動揺のある領域において優先的に周方向ベクトルを抽出し、デジタルフィルタ部によって、前記周方向ベクトルをデジタルフィルタリングして「低周波回転動揺に対する感度増加分」を検出する。

また、不快度推定装置は、低周波回転動揺補正量算出部によって、前記「低周波回転動揺に対する感度増加分」を画面全体での総和を取った上で動揺認知量推定用と同じ時間内加算をすることによって、不快度の上昇分に相当する低周波回転動揺補正量を算出する。

また、不快度推定装置は、領域別誘目度算出部によって、領域ごとの目立ちやすさの指標である領域別誘目度を算出し、誘目度補正量算出部によって、前記領域別誘目度と領域番号及び動き領域別短時間動揺エネルギーを用い、誘目度の偏りによる不快度の低下分を補正するための誘目度補正量を算出する。

また、不快度推定装置は、動き領域間距離正規化分散値算出部によって、前記領域別誘目度と動き領域の位置関係、及び動揺認知量を推定する際に用いられたのと同じ観視距離の値とに基づいて、動き領域間距離正規化分散値を算出する。ここで、動き領域間距離正規化分散値算出部は、前記動き領域番号、動き領域別短時間動揺エネルギー、領域ごとの誘目度及び観視距離の値を用い、想定する最も近い観視距離において、誘目度の分布状態によって当該算出値が最大となる場合に１となるように正規化した動き領域間距離正規化分散値を算出する構成であってもよい。

また、不快度推定装置は、動き領域間距離補正量算出部によって、前記動き領域間距離正規化分散値を用い、動き領域間距離の大きさによる不快度の低下分に相当する動き領域間距離補正量を算出する。すなわち、動き領域間距離補正量は、誘目度の分布状態における動き領域間距離の影響による不快度の低下分を、観視距離の影響を加味して補正するための補正量である。

また、不快度推定装置は、動揺認知量補正部によって、前記低周波回転動揺補正量と前記誘目度補正量及び前記動き領域間距離補正量を用い、入力された動揺認知量推定値を補正して、不快度との乖離分が除去された補正済み動揺認知量を得る。

また、総動揺エネルギー補正部によって、所定時間以上継続する動揺に対する不快度の蓄積効果分の補正を前記補正済み動揺認知量に施すことで補正済み総動揺エネルギーを得て、前記補正済み総動揺エネルギーを推定不快度として出力する。

かかる構成により、揺れの大きさとしては同じように感じる映像であっても、映像内で揺れが大きくかつ目立つ複数の領域間の距離が離れている場合に低下する不快度、つまり、誘目度の分布状態の違いによって動揺認知量との乖離が変化する不快度を高精度に推定することが可能になる。

また、本発明は、コンピュータを前記した不快度推定装置として機能させる不快度推定プログラムとしても具現化可能である。

本発明によると、８Ｋ映像の広視野観視条件を含む様々に異なる観視距離において、様々に異なる動揺が全画面中に分布している動揺映像を観視した際に生じる不快度を高精度に推定することが可能になる。

本発明の実施形態に係る不快度推定装置及び前処理部としての動揺認知量推定装置を備える不快度推定システムの構成例を示すブロック図である。 前処理部としての動揺認知量推定装置を示すブロック図である。 動きベクトルに対する動揺認知感度の周波数特性の一例を示すグラフである。 映像の輝度信号に対し空間周波数成分に応じて感度補正する特性の一例を示す図である。 動き領域間相関値に応じて領域別動揺エネルギーを補正する特性の一例を示す図である。 動き領域内部境界間相関値に応じて領域別動揺エネルギーを補正する特性の一例を示す図である。 動き領域の大きさと観視距離に応じて領域別動揺エネルギーを補正する特性の一例を示す図である。 動揺エネルギーを時間加算する際に乗じる時間窓の一例を示す図であり、低周波回転動揺補正量及び誘目度補正量を動揺認知量と同じ窓関数を掛けて時間積分する際の窓関数の一例を示す図である。 本発明の実施形態に係る不快度推定装置を示すブロック図である。 （ａ）は周（回転）方向の揺れを他の方向の揺れより優先して抽出する場合の周方向ベクトル算出の一例を示す図、（ｂ）は動揺のある領域の中心を回転中心とした周方向の揺れを抽出する通常の周方向ベクトル算出の一例を示す図である。 低周波回転動揺に対する感度増加分を検出するデジタルフィルタの伝達特性の一例を示す図である。 最小２次２次モーメントの算出例を説明するための図であり。（ａ）は、全画面における２次モーメントの算出例を説明するための模式図、（ｂ）は、中心から外周へ重みの大きな領域の順に同心円状に広がっていく場合の最小２次モーメントの算出例を説明するための模式図である。 動きが目立つ領域間の、観視距離を考慮した距離に応じ、動揺認知量に対して不快度との乖離分を補正する特性の一例を示す図である。

以下、本発明の実施形態について、適宜図面を参照しながら説明する。図１に示すように、本発明の実施形態に係る不快度推定システムＡは、動揺認知量推定装置１と、不快度推定装置２と、を備える。不快度推定装置２は、動揺認知量推定装置１と映像の特徴量とを用いることで、不快度推定誤差を抑制することを可能にしたものであり、機能部として、特徴量抽出部３と、不快度推定部４と、を備える。ただし、動揺認知量推定装置１の構成は、本発明の実施形態に係る不快度推定装置２の前処理を行う装置であるので、本発明の実施形態に供する上で、必ずしも下記構成例の形態を採るものではない。

＜動揺認知量推定装置＞
本発明の実施形態に供する前処理部に相当する動揺認知量推定装置１の一例は、図２に示すように、映像における動きベクトル及び設定された観視距離に基づいて画面動揺に対する動揺認知量を推定して出力するものであり、機能部として、領域別動きベクトル検出部１１０と、デジタルフィルタ部１２０と、動き領域判定部１３０と、動き領域番号付与部１４０と、水平・垂直比算出部１５０と、水平・垂直成分補正部１６０と、領域別空間周波数感度補正部１７０と、空間周波数エネルギー算出部１８０と、動き領域間相関値算出部１９０と、動き領域内部境界間相関値算出部２００と、動き領域別短時間動揺エネルギー算出部２１０と、時間窓内総動揺エネルギー算出部２２０と、を備える。ここで、動揺認知量推定装置１によって推定される動揺認知量は、動揺映像を観視した視聴者（観視者）が始めに認知する動揺量である。ここで、「始めに認知する」とは、視聴者が映像内容によって不快であるか否かというように感じ方に違いが生じる前に、視聴者が「揺れが大きいか否か」を純然と感じる、揺れの大きさの認知過程を意味している。

また、動揺認知量推定装置１は、本装置の使用者（例えば、映像の視聴者）による操作部（キーボード、マウス等）の操作結果、あるいは、センサ等による画面と視聴者との間の測距結果に基づいて、観視距離を取得する。観視距離は、映像の表示位置（スクリーン等）と視聴者の瞳孔との距離を、映像の高さ（画面高Ｈ）で割った値である。かかる観視距離を用いることによって、映像のサイズ（スクリーンサイズ）によらず、視野角による議論が可能となる。また、標準観視距離は、２Ｋ（ハイビジョン）では３Ｈ、４Ｋでは１．５Ｈ、８Ｋ（ＳＨＶ：スーパーハイビジョン）では０．７５Ｈである。観視距離は、例えば、０．７５Ｈ（ＳＨＶの標準観視距離）〜６Ｈ（ＮＴＳＣ（National Television System Committee）方式の標準観視距離）の連続的な値から適宜設定可能である。

かかる機能部のうち、領域別動きベクトル検出部１１０、デジタルフィルタ部１２０、動き領域判定部１３０、水平・垂直成分補正部１６０、領域別空間周波数感度補正部１７０及び空間周波数エネルギー算出部１８０の組み合わせに関しては、画面を分割した複数（１〜Ｎ）の領域に対応して複数（Ｎ）セットが設けられている。以下、全画面をＮに分割したそれぞれの領域を、「小領域」と呼ぶことがある。

＜領域別動きベクトル検出部＞
領域別動きベクトル検出部１１０は、当該検出部１１０に対応する領域の映像信号を取得し、取得された映像信号に基づいて、映像の時間的に隣接する画像（連続するフレーム）間における動きベクトル（フレーム間の差分）を時系列に検出する。また、領域別動きベクトル検出部１１０は、検出結果をデジタルフィルタ部１２０へ出力する。本実施形態において、領域別動きベクトル検出部１１０は、動きベクトルとして、動き量の水平方向成分、及び、垂直方向成分をそれぞれ検出する。

＜デジタルフィルタ部＞
デジタルフィルタ部１２０は、領域別動きベクトル検出部１１０によって検出された動きベクトルを取得し、取得された動きベクトルの水平方向及び垂直方向の各成分に周波数感度補正を施すことによって周波数感度補正済み動きベクトルを得た後、得られた周波数感度補正済み動きベクトルを動き領域判定部１３０、動き領域番号付与部１４０及び水平・垂直成分補正部１６０へ出力する。本実施形態において、デジタルフィルタ部１２０は、水平方向のデジタルフィルタ１２１と、垂直方向のデジタルフィルタ１２２と、を並列に備える。これは、１つのデジタルフィルタではベクトルの２成分を同時に処理できないからであり、両方とも特性は全く同じものであり、その特性は図３に示すものとなっている。

水平方向のデジタルフィルタ１２１は、動きベクトルの水平方向成分に周波数感度補正を施す。垂直方向のデジタルフィルタ１２２は、動きベクトルの垂直方向成分に周波数感度補正を施す。水平方向のデジタルフィルタ１２１及び垂直方向のデジタルフィルタ１２２は、それぞれ、水平・垂直各方向成分に対して図３に示す動揺認知量の周波数感度に相当するインパルスレスポンスを畳み込み積分するデジタルフィルタ処理を施すことによって、周波数感度補正済み動きベクトル（の水平方向成分及び垂直方向成分）を得る。ここで、各デジタルフィルタ１２１，１２２は、周波数感度に相当するインパルスレスポンスのフレーム時系列値を予め記憶している。このインパルスレスポンスのフレーム時系列値をそれぞれ、取得された動きベクトルの各成分の当該フレームの値に乗じた後、総和を取ることにより、畳み込み積分によるデジタルフィルタ処理を行うことができる。

水平方向のデジタルフィルタ１２１及び垂直方向のデジタルフィルタ１２２は、前記の畳み込み積分を用いず、離散フーリエ変換を用いて同じ処理を行うこともできる。この場合、各デジタルフィルタ１２１，１２２は、それぞれ、動揺周波数と動揺のフレーム間差分に対する感度との関係性を予め記憶している。そして、各デジタルフィルタ１２１，１２２は、取得された動きベクトルの各成分を離散フーリエ変換して全動揺周波数の各成分量を算出し、前記関係性を参照することによって、算出された全動揺周波数の各成分量に対して周波数感度補正を施す。ここでの周波数感度補正は、算出された各動揺周波数の成分量に各周波数感度の値を乗じることによって行われる。周波数感度補正後の全動揺周波数の各成分量を逆離散フーリエ変換することにより、周波数感度補正済み動きベクトル（の水平方向成分及び垂直方向成分）が得られる。ただし、離散フーリエ変換を用いる方法は、畳み込み積分を用いる方法よりも正確な周波数感度補正結果が得られるものの、計算量が畳み込み積分を用いる方法の数十倍以上となるため、１画面で多数の領域ごとに周波数補正を行う必要のある本発明の実施形態に係る動揺認知量推定装置１では、周波数感度補正誤差が数％以内に収まる畳み込み積分を用いる方法を用いる方が望ましい。

水平方向のデジタルフィルタ１２１及び垂直方向のデジタルフィルタ１２２は、それぞれ、得られた周波数感度補正済み動きベクトルの水平方向成分及び垂直方向成分を動き領域判定部１３０、動き領域番号付与部１４０及び水平・垂直成分補正部１６０へ出力する。

＜動き領域判定部＞
動き領域判定部１３０は、デジタルフィルタ部１２０から出力された周波数感度補正済み動きベクトルを取得し、周波数感度補正済み動きベクトルの大きさが予め定められた閾値以上である場合に、該当する領域において「動きあり」と判定する。本実施形態において、動き領域判定部１３０は、計算の簡略化のため、水平方向動き量（すなわち、周波数感度補正済み動きベクトルの水平方向成分の大きさ）の二乗と垂直方向動き量（すなわち、周波数感度補正済み動きベクトルの垂直方向成分の大きさ）の二乗との和が閾値（前記した閾値の二乗）以上である場合に、「動きあり」と判定することができる。動き領域判定部１３０は、「動きあり」の場合と「動きなし」の場合とで異なる値、例えば、「動きあり」の場合は「０」、「動きなし」の場合は「−１」に設定し、設定結果を動きの有無の判定結果として動き領域番号付与部１４０へ出力する。

＜動き領域番号付与部＞
動き領域番号付与部１４０は、デジタルフィルタ部１２０から出力された周波数感度補正済み動きベクトルと、動き領域判定部１３０の判定結果と、を取得し、取得された周波数感度補正済み動きベクトル及び判定結果に基づいて、動きありと判定された領域に対して同じ動き領域に属する領域ごとに動き領域番号を付与する。動き領域番号付与部１４０は、一の領域の補正済み動きベクトルと、隣接する領域の補正済み動きベクトルと、の差分ベクトルの大きさが所定値以下である場合に、一の領域と隣接する領域とが「同じ動き領域」であるとし、これらに同じ「動き領域仮番号」を付与する。かかる所定値は、画面の全体又はある程度まとまった領域の動きと主観評価実験で得られた動揺認知量との関係を用いて設定することができる。

さらに、動き領域番号付与部１４０は、隣接する「異なる動き領域」の境界において、フレーム間でいずれの領域がもう一方の領域を遮蔽しているかを調べることにより、領域の空間的な前後関係を判定する。動き領域番号付与部１４０は、領域の前後関係を全ての動き領域について調べた結果に基づいて、一番背後にある動き領域から順に「動き領域仮番号」を並べ直した後、それぞれの動き領域に並べ直した順番に沿い、「１」から順に正式な「動き領域番号」を付与する。以下、同じ動き領域番号が付与された領域（小領域）の集まりを、単に「動き領域」と呼ぶことがある。ここで、動き領域番号付与部１４０は、隣接していない複数の動き領域がいずれも同じ動き領域と隣接している場合には、複数の動き領域の前後関係を特定することができないが、全ての隣接する動き領域間で前後関係に破綻を来たさない限り、順序を入れ替えてもかまわない。ただし、隣接していない２つの動き領域ＡとＢがともに同じ動き領域Ｃと隣接するとともに同じ前後関係である場合であっても、動き領域Ｃと隣接する別の動き領域Ｄがやはり動き領域ＡとＢの両方と隣接するとともに異なる前後関係である場合には、動き領域Ａ，Ｂ，Ｃ，Ｄの前後関係の順序は一意に定まる。なお、動き領域の前後関係の順序の判定手法としては、映像における被写体のオクルージョンの検出手法等と同様の公知の手法を採用することができる。

ここで、動き領域番号付与部１４０は、一番背後にある動き領域から順に動き領域番号を付与することとしたが、動き領域番号の大小関係が前後関係を表していさえすればよいので、一番手前にある動き領域から順に動き領域番号を付与してもかまわない。

また、動き領域番号付与部１４０は、動きなしと判定された領域に関しては、動き領域判定部１３０と同様に「−１」を付与する。この値は「−１」に限定する必要はなく、動きありと判定された領域に付与された動き領域番号のいずれとも異なる特定の値でありさえすればよいが、動き領域番号が正の値であるので、最も大きい負の値である「−１」とするのが望ましい。

また、動き領域番号付与部１４０は、動きありと判定された領域であっても、隣接するいずれの領域とも異なる動きとなっている領域については、動き領域判定部１３０と同様に「０」を付与する。この値も「０」に限定する必要はなく、正の値である動き領域番号とも、動きなしと判定された領域番号である「−１」とも容易に区別できる「０」とするのが望ましい。

動き領域番号付与部１４０は、前記の手順で各領域に付与された動き領域番号を水平・垂直比算出部１５０、動き領域間相関値算出部１９０、動き領域内部境界間相関値算出部２００及び動き領域別短時間動揺エネルギー算出部２１０へ出力する。

＜水平・垂直比算出部＞
水平・垂直比算出部１５０は、動き領域番号付与部１４０から出力された動き領域番号を取得し、初めに、同じ領域番号が付与された動き領域の形状の１次モーメント（重心位置）を算出する。続いて、水平・垂直比算出部１５０は、同じ領域番号が付与された動き領域内の各領域と重心位置と各領域の位置との差の水平成分及び垂直成分のそれぞれの二乗の総和を、２次モーメントの水平成分と垂直成分として算出する。最後に、水平・垂直比算出部１５０は、２次モーメントの水平成分と垂直成分の比の平方根を計算することで、各動き領域の水平・垂直比を得て、得られた水平・垂直比を各領域の水平・垂直成分補正部１６０へ出力する。

ここで、後記する水平・垂直成分補正部１６０において動揺エネルギーの補正を行う際には、水平・垂直比の値を指数乗して用いるため、水平・垂直比算出部１５０から動き領域ごとに出力される水平・垂直比は、２次モーメントの水平成分と垂直成分の比の平方根ではなく、２次モーメントの水平成分と垂直成分の比、すなわち水平・垂直比の二乗値としてもよい。

＜水平・垂直成分補正部＞
水平・垂直成分補正部１６０は、水平・垂直比算出部１５０から出力された各領域における水平・垂直比と、デジタルフィルタ部１２０から出力された周波数感度補正済み動きベクトル、とを取得し、領域の前記周波数感度補正済み動きベクトルと前記水平・垂直比とを用いて、前記周波数感度補正済み動きベクトルの水平・垂直成分の大きさを前記水平・垂直比に応じて補正することによって水平・垂直比補正済み動きベクトルを得る。本実施形態において、水平・垂直成分補正部１６０は、水平・垂直比の値をα乗（水平・垂直成分補正部１６０からの出力が水平・垂直比ではなく水平・垂直比の二乗値である場合には０．５α乗）した後、周波数感度補正済み動きベクトルの水平成分及び垂直成分にそれぞれ乗じることにより水平・垂直比補正済み動きベクトルを得る。また、水平・垂直成分補正部１６０は、得られた水平・垂直比補正済み動きベクトルを、動き領域間相関値算出部１９０、動き領域内部境界間相関値算出部２００及び動き領域別短時間動揺エネルギー算出部２１０へ出力する。ここで、αの値については、α≒０．１８で誤差が最小となる。

＜領域別空間周波数感度補正部＞
領域別空間周波数感度補正部１７０は、領域の見えやすさの影響度に対する補正を行うのに必要となる空間周波数感度補正済みのエネルギーを算出する前処理として、領域別動きベクトル検出部１１０で画面を分割した領域ごとに、映像の輝度信号に対して図４に示す空間周波数感度補正を施す。この際、同じ映像の輝度信号であっても、観視距離によって空間周波数は異なるので、本装置の使用者が設定（入力）した、あるいはセンサによって計測された観視距離に応じて定まる空間周波数を用いる。ここで、領域別空間周波数感度補正部１７０は、観視距離と空間周波数との関係性を予め記憶している。そして、領域別空間周波数感度補正部１７０は、観視距離を取得するとともに、取得された観視距離を用いて前記関係性を参照することによって、観視距離に対応する空間周波数を読み出す。また、領域別空間周波数感度補正部１７０は、図４に示す空間周波数と相対感度との関係性を予め記憶している。そして、領域別空間周波数感度補正部１７０は、読み出された空間周波数を用いて前記関係性を参照することによって、空間周波数に対応する相対感度を読み出す。領域別空間周波数感度補正部１７０は、読み出された相対感度を用いることによって、映像の輝度信号に対して空間周波数感度補正を施す。なお、領域別空間周波数感度補正部１７０は、観視距離と相対感度との関係性を予め記憶している構成であってもよい。この場合には、領域別空間周波数感度補正部１７０は、取得された観視距離を用いて前記関係性を参照することによって、観視距離に対応する相対感度を読み出すことができる。

さらに、領域別空間周波数感度補正部１７０は、局所的な明暗コントラストの影響度を反映するために、空間周波数感度補正が施された映像の輝度信号（明暗コントラストの影響度の反映前）を、当該領域内の映像信号の直流成分の平方根で除すことにより、明暗コントラストの影響度が反映された空間周波数感度補正済み映像信号を得る。また、領域別空間周波数感度補正部１７０は、得られた空間周波数感度補正済み映像信号を空間周波数エネルギー算出部１８０へ出力する。

図４の特性は、水平・垂直・斜めのいずれ方向に対しても同等であるとしても大きな誤差を生じない（正確には斜め方向の高周波域の感度は水平・垂直方向での感度の８０％に低減するが、これによる動揺認知量の推定誤差の増加分は０．１ランク未満である）。したがって、領域別空間周波数感度補正部１７０は、グラフの原点を通る縦軸を中心軸として図４の特性を３６０°回転させた２次元空間周波数感度特性となる２次元インパルスレスポンスを、映像の輝度信号に対して畳み込み積分することによって空間周波数感度補正を施すことができる。

＜空間周波数エネルギー算出部＞
空間周波数エネルギー算出部１８０は、当該領域において領域別空間周波数感度補正部１７０から出力された空間周波数感度補正済み映像信号を取得し、取得された空間周波数感度補正済み映像信号を二乗した値の総和を求めることで、前記した領域の見えやすさの影響度に対する補正を行うのに必要となる空間周波数感度補正済みのエネルギーを算出する。また、空間周波数エネルギー算出部１８０は、算出された空間周波数感度補正済みエネルギーを動き領域別短時間動揺エネルギー算出部２１０へ出力する。

＜動き領域間相関値算出部＞
動き領域間相関値算出部１９０は、動き領域番号付与部１４０から出力された動き領域番号と、水平・垂直成分補正部１６０から出力された水平・垂直比補正済み動きベクトルと、を取得し、取得された動き領域番号及び水平・垂直比補正済み動きベクトルに基づいて動き領域間相関値を算出する。詳細には、動き領域間相関値算出部１９０は、動き領域境界における前景領域と背景領域との間の水平・垂直比補正済み動きベクトルの相関値を求めることで、前景と背景との動揺の相関の影響度に対する補正を行うのに必要となる動き領域間相関値を算出する。また、動き領域間相関値算出部１９０は、算出された動き領域間相関値を動き領域別短時間動揺エネルギー算出部２１０へ出力する。

本実施形態において、動き領域間相関値算出部１９０は、領域Ａと領域Ｂとの境界に隣接する領域Ａの画素群の動きベクトルの平均すなわち平均動きベクトルをＶ_Ａ、領域Ａと領域Ｂとの境界に隣接する領域Ｂの画素群の動きベクトルの平均すなわち平均動きベクトルをＶ_Ｂとしたとき、単位時間（例えば、２〜３秒）における動き領域間相関値Ｃ_１を、次式によって算出する。
Ｃ_１＝Σ（Ｖ_Ａ・Ｖ_Ｂ）／Σ（‖Ｖ_Ａ‖×‖Ｖ_Ｂ‖）
ここで、「・」は内積を表す。

＜動き領域内部境界間相関値算出部＞
動揺映像においては、同じ面積でかつ同じ動揺の大きさであっても、領域全体が動揺している場合の方が、領域の境界は動かず内部のみが動揺している場合よりも動揺認知量が大きくなる性質がある。動き領域内部境界間相関値算出部２００は、動き領域番号付与部１４０から出力された動き領域番号と、水平・垂直成分補正部１６０から出力された水平・垂直比補正済み動きベクトルと、を取得し、これらに基づいて動き領域内部境界間相関値を算出する。詳細には、動き領域内部境界間相関値算出部２００は、動き領域境界自体の動きベクトルと動き領域内部の水平・垂直比補正済み動きベクトルとの相関値を求めることで、前記した性質を反映する補正を施すのに必要になる動き領域内部境界間相関値を得て、動き領域別短時間動揺エネルギー算出部２１０へ出力する。

本実施形態において、動き領域内部境界間相関値算出部２００は、領域Ａの最外周（他領域との境界）の画素群の動きベクトルの平均すなわち平均動きベクトルをＶ_Ａ１、領域Ａの最外周以外（すなわち内部）の画素群の動きベクトルの平均すなわち平均動きベクトルをＶ_Ａ２としたとき、単位時間（例えば、２〜３秒）における動き領域内部境界間相関値Ｃ_２を、次式によって算出する。
Ｃ_２＝Σ（Ｖ_Ａ１・Ｖ_Ａ２）／Σ（‖Ｖ_Ａ１‖×‖Ｖ_Ａ２‖）
ここで、「・」は内積を表す。

＜動き領域別短時間動揺エネルギー算出部＞
動き領域別短時間動揺エネルギー算出部２１０は、動き領域番号付与部１４０から出力された各領域の動き領域番号と、水平・垂直成分補正部１６０から出力された水平・垂直比補正済み動きベクトルと、空間周波数エネルギー算出部１８０から出力された空間周波数感度補正済みエネルギーと、動き領域間相関値算出部１９０から出力された動き領域間相関値と、動き領域内部境界間相関値算出部２００から出力された動き領域内部境界間相関値と、本装置の使用者が設定（入力）した、あるいはセンサによって計測された観視距離と、を取得する。動き領域別短時間動揺エネルギー算出部２１０は、取得された各値に基づいて、同じ動き領域番号が付与された領域ごとに、当該領域内の水平・垂直比補正済み動きベクトルの大きさの二乗を、動きベクトルが検出された画像間の時間である短時間ごとにそれぞれ求めた個別短時間動揺エネルギーに対し、補正済み空間周波数エネルギー、動き領域間相関値、動き領域内部境界間相関値及び観視距離に応じて補正を施した後、当該動き領域内で補正が施された個別短時間動揺エネルギーの総和を求め、求められた総和を動き領域番号ごとの動き領域別短時間動揺エネルギーとして算出する。なお、補正済み空間周波数エネルギー、動き領域間相関値、及び動き領域内部境界間相関値に応じた補正は、前記した順番に行われてもよく、適宜順番を入れ替えて行われてもよい。ただし、観視距離に応じた補正は、他の補正済み個別短時間動揺エネルギーの大小関係によって補正量が異なるため、最後に補正を行う必要がある。

本実施形態において、動き領域別短時間動揺エネルギー算出部２１０は、初めに、同じ動き領域番号が付与された領域ごとに、各フレーム間において当該領域内の水平・垂直比補正済みベクトルの水平成分及び垂直成分の二乗和をそれぞれ求め、個別動揺エネルギーを得る。

次に、動き領域別短時間動揺エネルギー算出部２１０は、領域の見えやすさの影響度に対する補正として、補正済み空間周波数エネルギーをβ乗した値を個別動揺エネルギーに乗じる処理を施し、空間周波数補正済み個別動揺エネルギーを得る。ここで、β≒０．３のときに誤差が最小になる。

次に、動き領域別短時間動揺エネルギー算出部２１０は、前景と背景との動揺の相関の影響度に対する補正として、図５の特性で示される動き領域間相関値に応じた補正係数を空間周波数補正済み個別動揺エネルギーに乗じる処理を施し、動き領域間相関値補正済み個別動揺エネルギーを得る。ここで、動き領域別短時間動揺エネルギー算出部２１０は、動き領域間相関値と補正係数との関係性を予め記憶している。そして、動き領域別短時間動揺エネルギー算出部２１０は、取得された動き領域間相関値を用いて前記関係性を参照することによって、動き領域間相関値に対応する補正係数を読み出すことができる。

さらに、動き領域別短時間動揺エネルギー算出部２１０は、動き領域の境界自体がどの程度領域内部と同期して動揺しているのかによる補正として、図６の特性で示される動き領域内部境界間相関値に応じた補正係数を動き領域間相関値補正済み個別動揺エネルギーに乗じる処理を施し、動き領域内部境界間相関値補正済み個別動揺エネルギーを得る。ここで、動き領域別短時間動揺エネルギー算出部２１０は、動き領域内部境界間相関値と補正係数との関係性を予め記憶している。そして、動き領域別短時間動揺エネルギー算出部２１０は、取得された動き領域内部境界間相関値を用いて前記関係性を参照することによって、動き領域内部境界間相関値に対応する補正係数を読み出すことができる。

さらに、動き領域別短時間動揺エネルギー算出部２１０は、同じ動き領域番号の領域の大きさ（視野角）がある程度以上大きくなると、当該動き領域内での動き領域内部境界間相関値補正済み個別動揺エネルギーの加算が線形ではなく、当該エネルギーが低いものほどより加算される値が低くなる（広視野動揺に対して動揺認知量が飽和する）現象に対する補正として、まず、同じ動き領域番号の領域における動き領域内部境界間相関値補正済み個別動揺エネルギーに対して、大きいものから順に１、２、・・・と番号を付与する。次に、動き領域別短時間動揺エネルギー算出部２１０は、当該動き領域内部境界間相関値補正済み個別動揺エネルギーの番号に、設定された観視距離において画面をＮ分割した各領域の占める面積（Ｄｅｇ^２）を乗じ、その値（累積面積）に応じた補正係数を図７の特性から求め、当該動き領域内部境界間相関値補正済み個別動揺エネルギーに当該補正係数を乗じることで、観視距離補正済み個別動揺エネルギーを得る。ここで、動き領域別短時間動揺エネルギー算出部２１０は、累積面積と補正係数との関係性を予め記憶している。そして、動き領域別短時間動揺エネルギー算出部２１０は、算出された累計面積を用いて前記関係性を参照することによって、累計面積に対応する補正係数を読み出すことができる。ここで、図７に示した特性は、横軸が面積（Ｄｅｇ^２）ではなく面積の平方根（Ｄｅｇ）としている。これは、横軸が面積だと特性曲線が横軸のごく低い値から急激に減衰してしまうため視覚的に把握しにくいからであり、実際にはわざわざ平方根を計算せずに累積面積から補正係数を求めることができる横軸を面積とした特性を用いる。また、同じ映像を異なる観視距離で観視した場合には、観視距離が近い（すなわち、視野角が広い）ほど動揺振幅が大きくなる。したがって、観視距離が近いほど、観視距離補正済み個別動揺エネルギーは大きくなる。

最後に、動き領域別短時間動揺エネルギー算出部２１０は、同じ動き領域番号の動き領域内で観視距離補正済み個別動揺エネルギーの総和を算出することで、動き領域番号ごとの動き領域別短時間動揺エネルギーを得て、得られた動き領域番号ごとの動き領域別短時間動揺エネルギーを時間窓内総動揺エネルギー算出部２２０へ出力する。

＜時間窓内総動揺エネルギー算出部＞
時間窓内総動揺エネルギー算出部２２０は、動き領域別短時間動揺エネルギー算出部２１０から出力された動き領域番号ごとの動き領域別短時間動揺エネルギーを取得し、前記短時間ごとの全ての動き領域の動き領域別短時間動揺エネルギーに時間窓の伝達率を乗じたものを、時間窓分だけ非線形加算することによって動揺認知量に相当する時間窓内の総動揺エネルギーを得る。本実施形態において、時間窓内総動揺エネルギー算出部２２０は、初めに、全ての動き領域で動き領域別短時間動揺エネルギーの総和を求め、全画面短時間動揺エネルギーを得る。なお、時間窓は、動きベクトルが検出された画像間の時間である短時間よりも長い時間（例えば、後記する図８の例では、伝達率半減区間が３秒、ロールオフ区間を含めると５．５秒）に設定される。

次に、時間窓内総動揺エネルギー算出部２２０は、時刻ｔがｔ＝Ｔからｔ＝Ｔ−Ｍまで遡る時区間の全フレームにおいて、全画面短時間動揺エネルギーＰ_（ｔ）に図８の総合特性で示された時間窓の伝達率Ｗ_{（Ｔ−ｔ−ｄ）}を乗じた上でε乗して非線形加算するために、下記の計算を行い、時間窓内において動揺認知量推定値に相当する総動揺エネルギーＥ_（Ｔ）を得る。
Ｅ_（Ｔ）＝Σ（Ｗ_{（Ｔ−ｔ−ｄ）}・Ｐ_（ｔ））^ε

図８において、時間窓の伝達特性に相当する総合特性は、動きベクトルの大きさに対する不快度の減衰特性と、伝達率半減時刻０秒及び３秒のそれぞれ±ｄの時区間を遷移域とするコサインロールオフ特性、との積である合成特性の伝達率を二乗したものである。ここで、合成特性の伝達率を二乗するのは、エネルギーが動きベクトルの二乗に相当するからである。ここで、時間窓内総動揺エネルギー算出部２２０は、時間（時間窓）と総合特性の伝達率との関係性を予め記憶している。そして、時間窓内総動揺エネルギー算出部２２０は、時間窓を用いて前記関係性を参照することによって、時間窓に対応する伝達率を読み出すことができる。

ここで、−ｄは、現時点よりも未来の時刻となり、未来の時刻での動揺エネルギーは得られないため、前記の時間窓の伝達率は、Ｗ_{（Ｔ−ｔ−ｄ）}のように設定されている。これにより、動揺認知量は動揺が始まってから時間ｄだけ遅れて有効な値が得られることになるが、実際の動揺も１秒程度遅れて認知されていることから問題とはならない。また、Ｍは、Ｍ＝３＋２ｄとなる。また、εの値の誤差は、ε≒０．６７のときに最小となる。

そして、時間窓内総動揺エネルギー算出部２２０は、算出された総動揺エネルギーを動揺認知量推定値として不快度推定装置２（図１参照）へ出力する。

前記したように、動揺認知量推定装置１は、本発明の実施形態に係る不快度推定装置２の前処理を行う装置であるので、本発明の実施形態に供する上で、動揺認知量推定装置１は、前記した機能部による構成例以外の構成であってもよい。しかし、本発明の実施形態に係る不快度推定装置２では、不快度を推定するのに用いる映像の特徴量の一部として、図２に示す動揺認知量推定装置１の構成例でも用いられる、周波数感度補正済み動きベクトルの水平方向成分及び垂直方向成分、動き領域番号及び動き領域別短時間動揺エネルギーといった算出途上特徴量を流用しているので、動揺認知量推定装置１からこれらの特徴量が得られない構成である場合には、不快度推定装置２は、当該特徴量を算出するための機能部を備えることになる。

＜不快度推定装置＞
本発明の実施形態に係る不快度推定装置２は、図９に示すように、動揺認知量推定装置１から得られる動揺認知量と、動揺認知量推定装置１の算出途上特徴量であるところの、周波数感度補正済み動きベクトルの水平方向成分及び垂直方向成分、動き領域番号、及び動き領域別短時間動揺エネルギーと、映像内の部分領域において視線誘導につながる目の引きやすさの指標となる誘目度と、動揺認知量を推定する際に用いられたのと同じ観視距離の値と、に基づいて画面動揺に対する不快度を推定するものである。不快度推定装置２は、機能部として、回転動揺成分抽出部１０と、デジタルフィルタ部２０と、低周波回転動揺補正量算出部３０と、領域別誘目度算出部４０と、誘目度補正量算出部５０と、動き領域間距離正規化分散値算出部６０と、動き領域間距離補正量算出部７０と、動揺認知量補正部８０と、総動揺エネルギー補正部９０と、を備える。

かかる機能部のうち、回転動揺成分抽出部１０、デジタルフィルタ部２０、低周波回転動揺補正量算出部３０、領域別誘目度算出部４０、誘目度補正量算出部５０、動き領域間距離正規化分散値算出部６０及び動き領域間距離補正量算出部７０の組み合わせが、図１における特徴量抽出部３に相当する。また、動揺認知量補正部８０及び総動揺エネルギー補正部９０の組み合わせが、図１における不快度推定部４に相当する。

また、デジタルフィルタ部２０及び領域別誘目度算出部４０の組み合わせは、画面を分割した複数（１〜Ｎ）の領域に対応して複数（Ｎ）セットが設けられている。

＜回転動揺成分抽出部＞
回転動揺成分抽出部１０は、動揺認知量推定装置と同じように画面を複数に分割した領域中で動揺のある領域において、周方向ベクトルを抽出する。本実施形態において、回転動揺成分抽出部１０は、動揺認知量推定装置１から出力される算出途上特徴量のうち、周波数感度補正済み動きベクトルの水平方向成分及び垂直方向成分と、動き領域番号と、を取得し、全画面中で有効な動きがある領域（前記の番号付与の例では動き領域番号が０と−１ではない領域）における動きベクトルを用いて、回転角による動きベクトルの大きさの二乗の総和が最大になるとともに、回転角による動きベクトルと実際の前記補正済み動きベクトルとの二乗誤差の総和が最小となるときの回転中心と回転角とを算出する。さらに、回転動揺成分抽出部１０は、算出された回転中心及び回転角に基づいて、各領域の中心における周方向（回転中心への同心方向と直交する方向）のベクトルを算出し、当該ベクトルの水平方向成分と垂直方向成分との組を周方向ベクトルとしてデジタルフィルタ部２０へ出力する。

ここで、回転中心を定めるに際し、通常のカメラパラメータ計測時のように、並進運動（上下・左右の方向の動き）を取り除いた後の映像の中心（部分動揺映像では動きのある領域の中心）ではなく、回転エネルギーが最大かつ全ベクトルとの二乗誤差が最小となる回転中心と回転角とが用いられる。これは、手持ちカメラで傾き角が変動して撮影される場合、その回転中心は、画面中心から大きくズレている。そして、このような映像を見た人の大半は、「撮影対象の中心を移動した後に傾いた」とは感じずに、「撮影対象の中心から離れたある点を中心に傾いた」と感じることに基づいている。

ここで、図１０を参照してその違いについて説明する。領域１〜９が縦３×横３に配置された映像に関して、図１０（ａ）に示すように、アミの掛かっていない左上の４つの領域１，２，４，５のみに各動きベクトルがある場合において、回転中心は、×印のついた右下（領域５，６，８，９の境界）にあり、４つの動きベクトルは、すべて周方向ベクトルとなっている。それに対し、図１０（ｂ）に示す通常のカメラパラメータ計測では、×印のついた動きのある領域を中心にして並進運動ベクトル（×印をつけた中心位置（領域１，２，４，５の境界）に表記）と周方向ベクトルに分かれている。周方向ベクトルのエネルギーの総和である回転エネルギーは、図１０（ａ）の例では４０であるのに対し、図１０（ｂ）の例では８である。このローテートすなわち周方向（回転方向）の動揺の周期がもし６秒前後だったとしたら、これらの不快度は、０．３ランクも違ってしまうことになる。

＜デジタルフィルタ部＞
図９に示すように、デジタルフィルタ部２０は、回転動揺成分抽出部１０によって検出された周方向ベクトルを取得し、取得された周方向ベクトルをデジタルフィルタリングすることによって、低周波回転動揺に対する感度増加分ベクトルを検出する。本実施形態において、デジタルフィルタ部２０は、取得された動きベクトルすなわち周方向ベクトルの水平方向及び垂直方向の各成分に、低周波回転動揺に対する不快度の感度上昇分（周方向以外の方向の動揺に対する感度の比から１を引いた差分値）に相当する周波数感度補正を施すことによって「低周波回転動揺に対する感度増加分ベクトル」を得て、低周波回転動揺補正量算出部３０へ出力する。本実施形態において、デジタルフィルタ部２０は、水平方向のデジタルフィルタ２１と、垂直方向のデジタルフィルタ２２と、を並列に備える。これは、１つのデジタルフィルタではベクトルの２成分を同時に処理できないからであり、両方とも特性は全く同じものであり、その特性は図１１に示すものとなっている。

なお、水平方向のデジタルフィルタ２１及び垂直方向のデジタルフィルタ２２では、前記した動揺認知量推定装置１における水平方向のデジタルフィルタ１２１及び垂直方向のデジタルフィルタ１２２と同様に、畳み込み積分に代えて離散フーリエ変換を用いることも可能であるが、本発明の実施形態に係る不快度推定装置２では１画面で多数の領域ごとに周波数補正を行うため、演算量の少ない畳み込み積分を用いる方が望ましい。

＜低周波回転動揺補正量算出部＞
低周波回転動揺補正量算出部３０は、デジタルフィルタ部２０から出力された「低周波回転動揺に対する感度増加分ベクトル」を取得し、低周波回転動揺に対する感度増加分ベクトルの自乗値の画面全体での総和を取った上で時間内加算することによって低周波回転動揺補正量を算出する。ここで、低周波回転動揺補正量は、低周波回転動揺が存在することによる不快度の上昇分に相当する。本実施形態において、低周波回転動揺補正量算出部３０は、画面全体でのベクトルの自乗和の総和を取った上で、図８に示す動揺認知量推定用と同じ時間窓関数を乗じ、さらには、時間的に線形加算可能となるようなオーダー（次元）とするために、γ乗（ここでは０．６７乗）した上で時間加算をすることで、低周波回転動揺補正量を算出し、動揺認知量補正部８０へ出力する。図８の特性のうち、実際に時間窓関数として用いられるのは「総合特性」であり、これは、「減衰特性」に「ロールオフ特性」を乗じた「合成特性」を自乗したものである。低周波回転動揺補正量算出部３０は、「総合特性」を時間窓関数として用いることによって、計算量を抑えつつ、低周波回転動揺補正量を好適に算出することができる。

＜領域別誘目度算出部＞
領域別誘目度算出部４０は、映像信号を取得し、動揺認知量推定装置１と同じように分割された領域ごとに、視線誘導につながる目立ちやすさの指標である誘目度を算出して、誘目度補正量算出部５０へ出力する。ここで、誘目度はその妥当性の評価方法自体、視線の滞留率や、移動発生頻度など多数あって一意には定まっておらず、その算出法も、実際の視線誘導との相関係数が０．６以上あれば優良であるとされるものだけで数十種類にも及ぶ。本実施形態では、非営利目的で自由に利用でき、空間周波数・色分布・コントラスト・偏在率などに基づいて誘目度を算出する代表的なオープンソース・ソフトウェアである「Saliency Map Alogorithm」（Harel, Koch, Perona ： http://www.vision.caltec.edu/~harel/share/gbvs.zip）を用いて、動き領域ごとの平均誘目度を算出し、その大きさを１〜１６の１６段階に分類する。これは、動揺エネルギーが１６倍になると不快度の評定値が１ランク上がるので、同じ大きさの揺れであっても、誘目度によって不快度には最大１ランクの差が生じるように設定したからである。これは最善の設定であるわけではなく、あくまで実施形態の１つである。

＜誘目度補正量算出部＞
誘目度補正量算出部５０は、動揺認知量推定装置１から算出途上特徴量として得られる、動き領域番号及び動き領域別短時間動揺エネルギーと、領域別誘目度算出部４０から出力された領域ごとの誘目度と、を取得し、領域別誘目度に基づいて、誘目度の偏りによる不快度の低下分に相当する誘目度補正量を算出する。すなわち、画面上の各領域の誘目度に偏りがある場合には、不快度は、かかる偏りが大きくなるほど低下する。本実施形態において、誘目度補正量算出部５０は、最初に、動き領域ごとの平均誘目度を求める。次に、誘目度補正量算出部５０は、平均誘目度が最大だった動き領域の値で各動き領域の平均誘目度を除し、その値を１から差し引いた値を、各動き領域において動揺エネルギーから減じられる不快度の重み係数とする。この係数は、平均誘目度が最大だった動き領域では０、最大で１５／１６となる。続いて、誘目度補正量算出部５０は、この重み係数を動き領域別短時間動揺エネルギーに乗じた後に画面全体で総和を取った上で、かかる総和に対して、図８に示す動揺認知量推定用と同じ時間窓関数（すなわち、総合特性）を乗じる。さらに、誘目度補正量算出部５０は、時間的に線形減算可能となるようなオーダー（次元）とするために、時間窓関数が乗じられたものをγ乗（ここでは０．６７乗）した上で時間加算をすることで、誘目度補正量を算出し、算出された誘目度補正量を動揺認知量補正部８０へ出力する。

＜動き領域間距離正規化分散値算出部＞
動き領域間距離正規化分散値算出部６０は、動揺認知量推定装置１から算出途上特徴量として得られる、動き領域番号及び動き領域別短時間動揺エネルギーと、領域別誘目度算出部４０から出力された領域ごとの誘目度と、動揺認知量推定装置１で用いたものと同様に、本装置の使用者が設定（入力）した、あるいはセンサによって計測された観視距離と、を取得し、取得された動き領域番号、動き領域別短時間動揺エネルギー、誘目度及び観視距離の値に基づいて、動き領域間距離正規化分散値を算出する。ここで、動き領域間距離正規化分散値算出部６０は、前記動き領域番号、動き領域別短時間動揺エネルギー、領域ごとの誘目度及び観視距離の値を用い、想定する最も近い観視距離（後述する例では８Ｋ映像の標準観視距離である０．７５Ｈ）において、誘目度の分布状態によって当該算出値が最大となる場合に１となるように正規化した動き領域間距離正規化分散値を算出し、算出結果を動き領域間距離補正量算出部７０へ出力する。

本実施形態において、動き領域間距離正規化分散値算出部６０は、最初に、動き領域ごとの平均誘目度を求める。ここで、各小領域の中心座標を（ｘ，ｙ）（ただし、１＜ｘ＜Ｘ，１＜ｙ＜Ｙ，Ｎ＝Ｘ・Ｙ）とし、各小領域の誘目度をＳ_ｘｙとする。Ｋ個の動き領域があるとき、動き領域ごとの平均誘目度Ｓ_ｋ（１＜ｋ＜Ｋ）は、当該動き領域内にあるＭ個の小領域の各誘目度の総和をＭで除すことによって算出される。
Ｓ_ｋ＝（ΣＳ_ｘｙ）／Ｍ
（ただし、当該動き領域外の小領域では、Ｓ_ｘｙ＝０）

次に、動き領域間距離正規化分散値算出部６０は、同じ動き領域番号の動き領域において、まず、当該動き領域の平均誘目度Ｓ_ｋに当該動き領域の動揺領域別短時間動揺エネルギーＥ_ｋを乗じた上で、当該動き領域内にある小領域の数Ｍで除すことにより、当該動き領域内にある各小領域の重みＷ_ｋを算出する。ここで、Ｗ_ｋは、下記式によって算出される。
Ｗ_ｋ＝Ｓ_ｋ×Ｅ_ｋ／Ｍ

続いて、動き領域間距離正規化分散値算出部６０は、かかる重みＷ_ｋに基づいて全画面（全小領域数Ｎ）における重心位置（Ｇ_Ｘ，Ｇ_Ｙ）を算出する。距離の値（すなわち、座標値（ｘ，ｙ）の単位）は、入力された観視距離において観視した際の視角を用いる。ここで、Ｇ_ｘ，Ｇ_ｙは、下記式によって算出される。
Ｇ_Ｘ＝Σ（Ｗ_ｋ×ｘ）／Σ（Ｗ_ｋ）
Ｇ_Ｙ＝Σ（Ｗ_ｋ×ｙ）／Σ（Ｗ_ｋ）

続いて、動き領域間距離正規化分散値算出部６０は、全画面（全小領域数Ｎ）における２次モーメントＱ（重心位置からの距離の２乗に重みを乗じた値の総和）を算出する。ここで、Ｑは、下記式によって算出される。
Ｑ＝Σ（（（ｘ−Ｇ_Ｘ）^２＋（ｙ−Ｇ_Ｙ）^２）×Ｗ_ｋ）

これと並行して、動き領域間距離正規化分散値算出部６０は、同じ数と重みの小領域で全画面における２次モーメントＱの値が最小となる場合、すなわち、中心から外周へ重みの大きな領域の順に同心円状に広がっていく場合の最小２次モーメントを求める。ここで、最小２次モーメントの算出例を図１２を用いて説明する。図１２に示すように、小領域数Ｎが３６（Ｘ＝６、Ｙ＝６、Ｎ＝Ｘ×Ｙ）の画面において、動き領域数が４で、各動き領域の位置、面積、重みが図１２（ａ）となっている場合、２次モーメントＱが最小となる配置は、図１２（ｂ）のようになる。図１２（ａ）の全領域の高さが４Ｋ映像の画面高に等しく、取得された観視距離が４Ｋ映像の標準観視距離１．５Ｈであった場合、視角の単位は６Ｄｅｇ（つまり、Ｘ＝Ｙ＝３６Ｄｅｇ）となる。

ここで、面積がＡ_ｋ、重みがＷ_ｋである円形の動き領域ｋの２次モーメントＱ_ｋは、下記式によって算出される。
Ｑ_ｋ＝Ｗ_ｋ×Ａ_ｋ ^２／２π
図１２の例では、図１２（ａ）に示す２次モーメントＱが４５６であるのに対し、図１２（ｂ）に示す最小２次モーメントＱは、２９５．４となっている。

続いて、動き領域間距離正規化分散値算出部６０は、全画面における２次モーメントの値から最小２次モーメントの値を減じた値を、全画面における重みの総和で除すことによって、動き領域間距離の分散値を求める。図１２（ａ）の例では、動き領域間距離の分散値が１．６７３となっている。

最後に、動き領域間距離正規化分散値算出部６０は、動き領域間距離の分散値が最大となる、同じ重みが画面の最外周の対角位置の２点にある場合の値である「対角線距離の半分を自乗した値」で、動き領域間距離の分散値を除すことによって、最大値が１となるように正規化した動き領域間距離正規化分散値を算出する。図１２（ａ）の例において、想定する最も近い観視距離が８Ｋ映像の標準観視距離である０．７５Ｈとし、図１２（ａ）の全領域の高さが４Ｋ映像の画面高に等しく、観視距離が４Ｋ映像の標準観視距離１．５Ｈであった場合、対角線距離の半分を自乗した値が１４９．７となっている。したがって、動き領域間距離正規化分散値Ｕ’は、Ｕ’＝１．６７３／１４９．７＝０．０１１２となっている。この動き領域間距離正規化分散値の平方根が、８Ｋ映像の対角線距離で正規化した動き領域間距離に相当し、図１２（ａ）の例では０．１０６となっている。

＜動き領域間距離補正量算出部＞
動き領域間距離補正量算出部７０は、動き領域間距離正規化分散値算出部６０から出力された動き領域間距離正規化分散値を取得し、動き領域間距離の大きさによる不快度の低下を反映した動き領域間距離補正量を算出して動揺認知量補正部８０へ出力する。すなわち、画面上に複数の動き領域が存在する場合には、不快度は、動き領域間距離が大きくなるほど低下する。本実施形態において、図１３に示す動き領域間距離正規化分散値と短時間動き領域間距離補正量に相当する動揺エネルギー（振幅の二乗のオーダーの値）の相対感度との関係が求められており、動き領域間距離補正量算出部７０には、かかる関係が予め記憶されている。動き領域間距離補正量算出部７０は、取得された動き領域間距離正規化分散値に基づいて、予め記憶された前記関係から当該動き領域間距離正規化分散値に対応する相対感度を読み出し、読み出された相対感度に基づいて、動き領域間距離補正量（短時間動き領域間距離補正量）を得る。ここで、図１３の縦軸の相対感度の尺度は線形ではなく、相対感度を二乗した値が線形な間隔となる平方根スケールで表しており、横軸は値自体を動き領域間距離正規化分散値の平方根としている。これは、相対感度が動き領域間距離正規化分散値の上昇に伴って急激に減少するため、相対感度の減少具合をわかりやすく示そうとしたからである。

本実施形態において、動き領域間距離補正量算出部７０は、観視距離による影響が動き領域間距離正規化分散値算出部６０によって既に反映された動き領域間距離正規化分散値を用いて動き領域間距離補正量（短時間動き領域間距離補正量）を得る。すなわち、動き領域間距離補正量算出部７０は、動き領域間距離正規化分散値と、入力された観視距離の値に基づいて、動き領域間距離の大きさによる不快度の低下分に相当する動き領域間距離補正量を算出するものであるといえる。なお、動き領域間距離正規化分散値算出部６０が、観視距離の影響を考慮せずに動き領域間距離正規化分散値を算出し、動き領域間距離補正量算出部７０が、かかる動き領域間距離正規化分散値及び観視距離に基づいて、動き領域間距離補正量（短時間動き領域間距離補正量）を算出する構成であってもよい。

なお、観視距離が大きくなると（すなわち視角が小さくなると）、観視距離が反映された動き領域間距離正規化分散値は小さくなる。また、動き領域間距離正規化分散値が小さくなると、動き領域間距離補正量は大きくなる。また、動き領域間距離補正量が大きくなると、全視野面積が小さくなるため、不快度は小さくなる。
動き領域間距離補正量の一定値（例えば、１）からの減衰量は、観視距離が大きい程小さくなる。また、不快度のエネルギーは、例えば観視距離が２倍になると、１／１６になる。したがって、不快度は、観視距離が大きくなる程小さくなる。

続いて、動き領域間距離補正量算出部７０は、単位時間（前記動揺認知量推定装置の例では３秒プラス遷移時間）の動揺エネルギーの時間蓄積効果に合わせるために、図８に示す動揺認知量推定用の時間窓関数の各時刻の値をγ乗（ここでは０．６７乗）した値を重みとして同じくγ乗した短時間動き領域間距離補正量に乗じたものの総和を、重みの総和で除することにより、時間窓内における動き領域間距離補正量（時間的に線形加算可能なオーダーの値）を算出する。

＜動揺認知量補正部＞
動揺認知量補正部８０は、動揺認知量推定装置１から出力された動揺認知量推定値と、低周波回転動揺補正量算出部３０から出力された低周波回転動揺補正量と、誘目度補正量算出部５０から出力された誘目度補正量と、動き領域間距離補正量算出部７０から出力された動き領域間距離補正量と、を取得し、各補正量によって動揺認知量推定値を補正することにより、補正済み動揺認知量を算出し、総動揺エネルギー補正部９０へ出力する。

本実施形態において、動揺認知量補正部８０は、まず、動揺認知量推定値及び誘目度補正量のそれぞれを振幅の二乗のオーダーの値（動揺認知量が時間的に線形加算可能な値のオーダーである場合は１．５乗した値）に変換した上で、低周波回転動揺補正済み動揺認知量から誘目度補正量を減算することで、誘目度の偏りによる不快度低下の影響分を補正した誘目度補正済み動揺認知量を得る。

次に、動揺認知量補正部８０は、誘目度補正済み動揺認知量及び動き領域間距離補正量のそれぞれを振幅の二乗のオーダーの値（動揺認知量が時間的に線形加算可能な値のオーダーである場合は１．５乗した値）に変換（誘目度補正済み動揺認知量は既に変換済みなので、動き領域間距離補正量のみ変換）した上で、誘目度補正済み動揺認知量に動き領域間距離補正量を乗じることで、動き領域間距離の大きさによる不快度低下分を補正した動き領域間距離補正済み動揺認知量を得る。

続いて、動揺認知量補正部８０は、動き領域間距離補正済み動揺認知量及び低周波回転動揺補正量のそれぞれを振幅の二乗のオーダーの値（動揺認知量が時間的に線形加算可能な値のオーダーである場合は１．５乗した値）に変換（動き領域間距離補正済み動揺認知量は既に変換済みなので、低周波回転動揺補正量のみ変換）した上で線形和を求めることで、低周波回転動揺補正済み動揺認知量を得る。

最後に、動揺認知量補正部８０は、元の不快度の時間蓄積効果を線形加算で表せる動揺認知量のオーダー（振幅の二乗値の０．５乗）に変換するため、低周波回転動揺補正済み動揺認知量を０．５乗（動揺認知量が時間的に線形加算可能な値のオーダーである場合は０．７５乗）することで最終的な補正済み動揺認知量を求める。

本実施形態において動揺認知量に補正を加える処理の順番は、誘目度補正、動き領域間距離補正、低周波回転動揺補正の順になっているが、これらの処理が互いに完全に独立しているならば、順序を入れ替えても何ら問題がない。しかし、誘目度補正と低周波回転動揺補正とは完全に独立しているものの、動き領域間距離補正は誘目度補正した動揺認知量に相対感度を乗じる処理であるため、前記した順序以外の順で補正を行うと、誤った補正結果が出力されることになってしまうので、処理の順序は前記した順であることが望ましい。

＜総動揺エネルギー補正部＞
総動揺エネルギー補正部９０は、所定時間（ここでは、３秒）以上継続する動揺に対する不快度の蓄積効果分の補正を前記補正済み動揺認知量に施すことで、推定不快度に相当する補正済み総動揺エネルギーを得る。本実施形態において、総動揺エネルギー補正部９０は、所定時間（ここでは、３秒）以上の長時間動揺が継続する場合の時間蓄積効果を反映した不快度を得るためのものであって、３秒前までの総動揺エネルギー、すなわち時間蓄積効果補正後不快度が大きいほど計測時点での時間蓄積効果補正後不快度が大きくなるように、動揺認知量補正部８０から出力された計測時点の直前３秒間の不快度に相当する補正済み動揺認知量を補正して時間蓄積効果補正後不快度を算出し、算出された時間蓄積効果補正後不快度を推定不快度として出力する。本発明の実施形態において、総動揺エネルギー補正部９０は、加算部９１、記憶部９２及び乗算部９３を備えており、３秒前までの総動揺エネルギー、すなわち時間蓄積効果補正後不快度を記憶部９２で記憶しておき、計測時点において、３秒前までの総動揺エネルギーに乗算部９３で定数δ（０＜δ＜１）を乗じたものを、加算部９１で動揺認知量補正部８０から出力された計測時点の直前３秒間の不快度に相当する補正済み動揺認知量に加算することで、計測時点での時間蓄積効果補正後不快度を得る構成になっている。ここで、定数δの値は動揺が長時間継続する場合の不快度の蓄積効果と同様が止まった後の減衰効果とを同時に具現化するものであり、δ＝０．３０５程度とした場合に不快度の推定誤差が最小となる。なお、総動揺エネルギー補正部９０は、１回目の補正時には、前フレームの出力が記憶部９２に記憶されていないため、無補正の不快度を出力する。

ここで出力される推定不快度は、対数を取ると心理評価のカテゴリー値と線形な対応関係となるので、推定不快度の対数を取って心理評価のカテゴリー値を出力しても良い。また、推定不快度は、ディスプレイへの数字表示やグラフ表示をすることで、利用者への視覚的通知を行うだけではなく、推定不快度の大きさに準じて異なる音をスピーカ等の音響装置で出力することにより、聴覚的通知を行うこともできる。。

＜不快度推定装置＞
続いて、図１に戻り、不快度推定装置２の特徴量抽出部３及び不快度推定部４について説明する。

＜特徴量抽出部＞
特徴量抽出部３は、映像信号と、動揺認知量推定装置１が動揺認知量を推定する際に求めた算出途上特徴量であるところの、周波数感度補正済み動きベクトルの水平方向成分及び垂直方向成分、動き領域番号、及び動き領域別短時間動揺エネルギーと、を取得し、取得された情報に基づいて、動揺認知量以外に不快度の要因となる映像の物理的特徴量を抽出し、抽出した物理的特徴量を不快度推定部４へ出力する。動揺認知量以外の不快度の要因となる映像の物理的特徴量として前記の説明例で取り上げたのは、特定の動揺成分（低周波回転動揺）や、時間変動、そして、空間周波数・色分布・コントラスト・偏在率等によって求まる「目立ちやすさ」の指標である誘目度及び誘目度の分布状態となっている。図９の例では、回転動揺成分抽出部１０、デジタルフィルタ部２０、低周波回転動揺補正量算出部３０、領域別誘目度算出部４０、誘目度補正量算出部５０、動き領域間距離正規化分散値算出部６０及び動き領域間距離補正量算出部７０の組み合わせが特徴量抽出部３に相当する。

＜不快度推定部＞
不快度推定部４は、前記の説明例では、動揺認知量推定装置１で推定された動揺認知量と、特徴量抽出部３によって抽出された物理的特徴量（実際には補正量として出力されたもの）とを取得し、動揺認知量に補正を加えた後に、時間蓄積効果と時間減衰効果とを反映した非線形時間加算を行うことによって推定不快度を算出して、ディスプレイ、スピーカ等の外部装置（不快度を利用者へ通知する通知部）へ出力する。図９の例では、動揺認知量補正部８０及び総動揺エネルギー補正部９０の組み合わせが不快度推定部４に相当する。

＜動作例＞
続いて、本発明の実施形態に係る不快度推定装置２の動作例について、図９を参照して説明する。図１に示した動揺認知量推定装置１の構成は、本発明の実施形態に係る不快度推定装置２へ動揺認知量を出力する例を示したものなので、ここでは説明の対象外とする。

まず、回転動揺成分抽出部１０が、周波数感度補正済み動きベクトルの水平方向成分及び垂直方向成分を用いて、各領域の中心における周方向ベクトルを算出し、算出された周方向ベクトルをデジタルフィルタ部２０へ出力する。

続いて、デジタルフィルタ部２０が、周方向ベクトルの水平及び垂直の各方向成分に、周波数感度補正を施して「低周波回転動揺に対する感度増加分ベクトル」を得て、得られた「低周波回転動揺に対する感度増加分ベクトル」を低周波回転動揺補正量算出部３０へ出力する。

続いて、低周波回転動揺補正量算出部３０が、画面全体で「低周波回転動揺に対する感度増加分ベクトル」の自乗和の総和を取った値に時間窓関数を乗じ、０．６７乗した上で時間加算して低周波回転動揺補正量を算出し、算出された低周波回転動揺補正量を動揺認知量補正部８０へ出力する。

その一方、領域別誘目度算出部４０が、映像信号を用いて、領域ごとに誘目度を算出して、算出された誘目度を誘目度補正量算出部５０へ出力する。

続いて、誘目度補正量算出部５０が、領域ごとの誘目度を用いて算出した動き領域ごとの平均誘目度の大きさの順に従って、各動き領域で動揺エネルギーから減じられる不快度の重み係数を算出した上で、この重み係数を動き領域別短時間動揺エネルギーに乗じてから、さらに時間窓関数を乗じた値を０．６７乗したものを時間加算することで、誘目度補正量を算出し、算出された誘目度補正量を動揺認知量補正部８０へ出力する。

その一方、動き領域間距離正規化分散値算出部６０が、入力された観視距離における視角を距離として、当該動き領域の各小領域に等分に割り振った動揺領域別短時間動揺エネルギーに動き領域ごとの平均誘目度を乗じた値を重みとした場合の全画面における２次モーメントの値から、同じ数と重みの小領域で全画面における２次モーメントが最小となる値を減じた上で、全画面における重みの総和で除すことによって動き領域間距離の分散値を求め、さらに、動き領域間距離の分散値の最大値である「最も近い観視距離で観視した際の全画面の対角線の視角の半分を自乗した値」で動き領域間距離の分散値を除すことで、最大値が１となるように正規化した動き領域間距離正規化分散値を算出し、動き領域間距離補正量算出部７０へ出力する。

続いて、動き領域間距離補正量算出部７０が、動き領域間距離正規化分散値から求まる動揺エネルギーの相対感度を０．６７乗した値に動揺認知量推定用の時間窓関数を０．６７乗した値を重みとして乗じたものの総和を重みの総和で除すことで、動き領域間距離補正量を算出し、算出された動き領域間距離補正量を動揺認知量補正部８０へ出力する。

続いて、動揺認知量補正部８０が、動揺認知量推定値、低周波回転動揺補正量と誘目度補正量及び動き領域間距離補正量をそれぞれ振幅の二乗のオーダーの値に変換した上で、まず動揺認知量から誘目度補正量を減算して誘目度補正済み動揺認知量を得て、得られた誘目度補正済み動揺認知量に動き領域間距離補正量を乗じることで、動き領域間距離補正済み動揺認知量を得て、得られた動き領域間距離補正済み動揺認知量に低周波回転動揺補正量を加算することで、低周波回転動揺補正済み動揺認知量を得て、不快度の時間加算が線形で行えるように低周波回転動揺補正済み動揺認知量を０．５乗することで最終的な補正済み動揺認知量を得て、得られた補正済み動揺認知量を総動揺エネルギー補正部９０へ出力する。

最後に、総動揺エネルギー補正部９０が、３秒前までの時間蓄積効果補正後不快度が大きいほど計測時点での時間蓄積効果補正後不快度が大きくなるように、直前３秒間の不快度に相当する誘目度補正済み動揺認知量を補正して時間蓄積効果補正後不快度を算出し、算出された時間蓄積効果補正後不快度を推定不快度として出力する。

本発明の実施形態に係る不快度推定装置２は、動揺認知量推定装置１によって得られた動揺認知量推定値と、映像の物理的特徴量とを用いることで、広視野提示された動揺映像に対する不快度推定精度を、従来の不快度推定装置では達しえなかった水準にまで高めることができる。
また、不快度推定装置２によって得られた推定不快度は、映像コンテンツ制作者が画面動揺を低減するような映像修正を施す際に、判断基準として使用することも可能である。すなわち、不快度推定装置２は、映像コンテンツ制作者によって制作段階で用いられる場合には、映像の良否の判定、映像に含まれる画面動揺をどの程度まで低減すべきかの設定目標として、推定不快度を映像コンテンツ制作者に提示することができるので、制作に要する時間、労力及びコストの削減を図ることができるのに加え、供給される映像コンテンツの安全性を向上させることができる。
また、不快度推定装置２が視聴者側で用いられる場合には、画面動揺に関して安全であることを保証せずに制作、流通された映像に対し、視聴前又は視聴中の表示直前に不快度を推定してディスプレイ又はスピーカへ出力することで警告を発することができるので、映像酔いによる健康被害を防止することが可能になる。

以上、本発明の実施形態について説明したが、本発明は前記実施形態に限定されず、本発明の要旨を逸脱しない範囲で適宜設計変更可能である。例えば、本発明は、コンピュータを前記不快度推定装置２として機能させる不快度推定プログラムとして具現化することも可能である。あるいは、推定不快度を心理評価値と線形的に対応するようにするために、前記総動揺エネルギーの対数値を推定不快度として出力することとしてもよい。また、不快度推定装置２は、動揺認知量推定装置１と一体化された装置としても具現化可能である。また、不快度推定装置２は、観視距離を検出するセンサの検出結果を取得し、取得された検出結果に基づいて不快度を推定する構成であってもよい。この場合には、不快度推定装置２は、視聴者が個人の場合の不快度を推定するための個人観視限定ツール、又は、視聴者が複数の場合に映像に最も近い視聴者の不快度を推定するための補助的なツールとして利用することが可能である。

Ａ 不快度推定システム
１ 動揺認知量推定装置
２ 不快度推定装置
３ 特徴量抽出部
４ 不快度推定部
１０ 回転動揺成分抽出部
２０ デジタルフィルタ部
２１ 水平方向のデジタルフィルタ
２２ 垂直方向のデジタルフィルタ
３０ 低周波回転動揺補正量算出部
４０ 領域別誘目度算出部
５０ 誘目度補正量算出部
６０ 動き領域間距離正規化分散値算出部動き領域間距離補正量算出部
７０ 動き領域間距離補正量算出部
８０ 動揺認知量補正部
９０ 総動揺エネルギー補正部
９１ 加算部
９２ 記憶部
９３ 乗算部
１１０ 領域別動きベクトル検出部
１２０ デジタルフィルタ部
１２１ 水平方向のデジタルフィルタ
１２２ 垂直方向のデジタルフィルタ
１３０ 動き領域判定部
１４０ 動き領域番号付与部
１５０ 水平・垂直比演算部
１６０ 水平・垂直成分補正部
１７０ 領域別空間周波数感度補正部
１８０ 空間周波数エネルギー算出部
１９０ 動き領域間相関値算出部
２００ 動き領域内部境界間相関値算出部
２１０ 動き領域別短時間動揺エネルギー算出部
２２０ 時間窓内総動揺エネルギー算出部

Claims (2)

1. 映像における動揺認知量と物理的特徴量に基づいて画面動揺に対する不快度を推定する不快度推定装置であって、
   画面を複数に分割した領域中で動揺のある領域において、周方向ベクトルを抽出する回転動揺成分抽出部と、
   前記周方向ベクトルをデジタルフィルタリングすることによって、低周波回転動揺に対する感度増加分ベクトルを検出するデジタルフィルタ部と、
   前記低周波回転動揺に対する感度増加分ベクトルの自乗値の前記画面全体での総和を取った上で時間内加算することによって、不快度の上昇分に相当する低周波回転動揺補正量を算出する低周波回転動揺補正量算出部と、
   前記領域の目立ちやすさの指標である領域別誘目度を前記領域ごとに算出する領域別誘目度算出部と、
   前記領域別誘目度に基づいて、誘目度の偏りによる不快度の低下分に相当する誘目度補正量を算出する誘目度補正量算出部と、
   前記領域別誘目度と動き領域の位置関係とに基づいて、動き領域間距離正規化分散値を算出する動き領域間距離正規化分散値算出部と、
   前記動き領域間距離正規化分散値と、入力された観視距離の値に基づいて、動き領域間距離の大きさによる不快度の低下分に相当する動き領域間距離補正量を算出する動き領域間距離補正量算出部と、
   前記低周波回転動揺補正量、前記誘目度補正量及び前記動き領域間距離補正量を用いて、入力された動揺認知量推定値を補正し、単位時間の不快度に相当する補正済み動揺認知量を得る動揺認知量補正部と、
   所定時間以上継続する動揺に対する不快度の蓄積効果分の補正を前記補正済み動揺認知量に施すことで、推定不快度に相当する補正済み総動揺エネルギーを得る総動揺エネルギー補正部と、
   を備えることを特徴とする不快度推定装置。
2. コンピュータを請求項１に記載の不快度推定装置として機能させることを特徴とする不快度推定プログラム。

Images