JP2011238120A - Discomfort degree estimation device and discomfort degree estimation program - Google Patents

Discomfort degree estimation device and discomfort degree estimation program Download PDF

Info

Publication number
JP2011238120A
JP2011238120A JP2010110401A JP2010110401A JP2011238120A JP 2011238120 A JP2011238120 A JP 2011238120A JP 2010110401 A JP2010110401 A JP 2010110401A JP 2010110401 A JP2010110401 A JP 2010110401A JP 2011238120 A JP2011238120 A JP 2011238120A
Authority
JP
Japan
Prior art keywords
discomfort
global motion
unit
types
degree
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.)
Granted
Application number
JP2010110401A
Other languages
Japanese (ja)
Other versions
JP5391145B2 (en
Inventor
Makoto Tadenuma
眞 蓼沼
Hisaya Morita
寿哉 森田
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
Nippon Hoso Kyokai NHK
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 Nippon Hoso Kyokai NHK filed Critical Nippon Hoso Kyokai NHK
Priority to JP2010110401A priority Critical patent/JP5391145B2/en
Publication of JP2011238120A publication Critical patent/JP2011238120A/en
Application granted granted Critical
Publication of JP5391145B2 publication Critical patent/JP5391145B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Studio Devices (AREA)
  • Image Analysis (AREA)

Abstract

PROBLEM TO BE SOLVED: To provide a discomfort degree estimation device for estimating a discomfort degree to be sensed by a viewer by global motions, based on the global motions.SOLUTION: The discomfort degree estimation device 1 includes: a global motion detection part 10 for detecting a plurality of kinds of global motions of an image; a Fourier transformation part 20 for calculating frequency components of the plurality of kinds of global motions by performing discrete Fourier transformation for a longer time as a frequency is lower; a weighting part 30 for performing weighting in the frequency components of the plurality of kinds of global motions; a discomfort degree calculation part 40 for calculating a discomfort degree by calculating a nonlinear sum of the weighted frequency components of the plurality of kinds of global motions and transforming the calculated nonlinear sum into a logarithm; and a discomfort degree correction part 50 for correcting the discomfort degree by multiplying a previous time discomfort degree by a constant which is larger than zero and smaller than one, and adding the previous time discomfort degree after the multiplication of the constant to a this time discomfort degree.

Description

本発明は、映像の不快度を推定する不快度推定装置及び不快度推定プログラムに関する。   The present invention relates to a discomfort level estimation apparatus and a discomfort level estimation program for estimating a discomfort level of an image.

家庭で放送、録画映像、パーソナルコンピュータ映像等を視聴する際のディスプレイが大画面化するのに伴い、画面動揺(映像の揺れ)、いわゆるグローバルモーションによって誘発される不快感が増大する傾向があり、場合によっては視聴者が映像酔いを起こして健康被害に至る例もある。かかる事態を予防するため、手持ち撮影用ビデオカメラ等には「手ぶれ補正」機能の搭載がほぼ必須となっている。また、撮影された映像に含まれるグローバルモーションを軽減する「映像スタビライズ」ソフトウェアも数多く開発されている。さらに、グローバルモーションと脈拍・血圧等との生理的な指標との関係を探ることによって、不快な状態を推定しようとする試みも行われている(非特許文献1参照)。   As the display when watching broadcasts, recorded videos, personal computer videos, etc. at home increases in screen size, there is a tendency for discomfort induced by screen motion (video motion), so-called global motion, to increase. In some cases, viewers may experience video sickness resulting in health damage. In order to prevent such a situation, it is almost essential to install a “camera shake correction” function in a hand-held video camera or the like. In addition, a lot of “Video Stabilization” software has been developed to reduce the global motion contained in the captured video. Furthermore, an attempt has been made to estimate an unpleasant state by searching for a relationship between global motion and physiological indices such as pulse and blood pressure (see Non-Patent Document 1).

「映像酔いガイドライン検証システムの実用化に関するフィージビリティスタディ報告書」、財団法人機械システム振興協会、平成20年3月“Feasibility Study Report on Practical Use of Visual Sickness Guidelines Verification System”, Mechanical Systems Promotion Association, March 2008

しかし、前記した「手ぶれ補正」では長周期の揺れは補正ができない、あるいは不十分であり、「映像スタビライズ」では、短い期間内の映像の揺れを無くすか低減したり、滑らかな画面移動に変換したりするだけで、どのような性質でどの程度の大きさの揺れならば残っていても良いといった客観的な基準に則っているものではない。また、この客観的な基準として脈拍・血圧等の生理的な指標を用いることができたとしても、明示できるのは健康被害に達する条件および、その条件に対する現在の生理的状態の程度のみである。   However, the above-mentioned “camera shake correction” cannot correct or inadequate long-period shaking, and “image stabilization” eliminates or reduces image shaking within a short period or converts it to smooth screen movement. It is not based on an objective standard that it can be left as long as it is shaken by what kind of property and how much. In addition, even if physiological indices such as pulse and blood pressure can be used as this objective standard, it is possible to specify only the conditions for reaching a health hazard and the current physiological state for those conditions. .

これに対し、放送、映画等の映像コンテンツの供給者は、映像制作段階において、映像の安全性のみならず、視聴者の大半にグローバルモーションに伴う不快感を催させないように留意する必要がある。しかしながら、前記した従来の技術は、グローバルモーションに対する視聴者の心理的な「不快度」を具体的に把握するものではなく、グローバルモーションによる「不快度」を基準として映像の良否を判定するには不向きであった。   On the other hand, suppliers of video content such as broadcasts and movies need to be careful not to cause discomfort associated with global motion to the majority of viewers in addition to video safety during the video production stage. . However, the above-described conventional technology does not specifically grasp the viewer's psychological “discomfort” with respect to global motion, but to determine the quality of a video based on the “discomfort” due to global motion. It was unsuitable.

本発明は、前記した事情に鑑みて創案されたものであり、グローバルモーションによって視聴者が感じる不快度をグローバルモーションに基づいて推定することが可能な不快度推定装置及び不快度推定プログラムを提供することを課題とする。   The present invention has been developed in view of the above-described circumstances, and provides a discomfort level estimation device and a discomfort level estimation program capable of estimating a discomfort level felt by a viewer based on global motion based on the global motion. This is the issue.

前記課題を解決するため、本発明の不快度推定装置は、映像における画面動揺を示すグローバルモーションに基づいて不快度を推定する不快度推定装置であって、グローバルモーション検出部と、離散フーリエ変換部と、重み付け部と、不快度算出部と、不快度補正部と、を備える。   In order to solve the above-mentioned problem, a discomfort degree estimation apparatus according to the present invention is a discomfort degree estimation apparatus that estimates a discomfort degree based on a global motion indicating screen motion in a video, and includes a global motion detection unit and a discrete Fourier transform unit. And a weighting unit, a discomfort calculation unit, and a discomfort correction unit.

かかる構成により、不快度推定装置は、グローバルモーション検出部によって、前記映像の時間的に隣接する画像間における複数種類のグローバルモーションを検出し、離散フーリエ変換部によって、前記複数種類のグローバルモーションに対して、それぞれ、周波数が低いほど長い期間で離散フーリエ変換を施すことによって、前記複数種類のグローバルモーションの周波数成分を算出し、重み付け部によって、前記複数種類のグローバルモーションの周波数成分に対して、それぞれ重み付けを施し、不快度算出部によって、前記複数種類のグローバルモーションの周波数成分の非線形和を算出することによって、不快度を算出し、不快度補正部によって、前回の前記不快度が大きいほど今回の前記不快度が大きくなるように今回の前記不快度を補正する。   With this configuration, the discomfort level estimation device detects a plurality of types of global motion between temporally adjacent images of the video by a global motion detection unit, and performs a discrete Fourier transform unit on the plurality of types of global motion. The frequency components of the plurality of types of global motion are calculated by performing discrete Fourier transform in a longer period as the frequency is lower, and the frequency components of the plurality of types of global motion are respectively calculated by the weighting unit. Weighting is performed, and the discomfort degree calculation unit calculates a non-linear sum of the frequency components of the plurality of types of global motions, thereby calculating the discomfort degree. In order to increase the discomfort level, Degrees to correct for.

この不快度推定装置では、視聴者の生理的指標によらず、グローバルモーションに基づいて不快度を推定することができる。また、重み付け部が複数種類のグローバルモーションに対してそれぞれ重み付けを施すので、複数種類のグローバルモーションの周波数成分と、それに起因する不快度と、の関係を調整し、複数種類のグローバルモーションに基づく不快度を同列に扱うことができる。また、不快度補正部が前回の不快度に基づいて今回の不快度を補正するので、グローバルモーションが長時間継続する場合の蓄積効果を反映した不快度を出力することができる。   With this discomfort level estimation device, the discomfort level can be estimated based on global motion regardless of the viewer's physiological index. In addition, since the weighting unit weights multiple types of global motion, the relationship between the frequency components of multiple types of global motion and the level of discomfort resulting from them is adjusted, and unpleasantness based on multiple types of global motion. Degrees can be handled in the same row. Further, since the discomfort level correcting unit corrects the current discomfort level based on the previous discomfort level, it is possible to output the discomfort level reflecting the accumulation effect when the global motion continues for a long time.

前記グローバルモーション検出部は、前記複数種類のグローバルモーションとして、上下移動、左右移動、回転、及び、ズームに関するグローバルモーションを検出する構成であってもよい。   The global motion detection unit may be configured to detect global motion related to vertical movement, horizontal movement, rotation, and zoom as the plurality of types of global motion.

前記重み付け部は、前記フーリエ変換部によって算出された前記複数種類のグローバルモーションの周波数成分に対して、それぞれ、前記複数種類のグローバルモーションの周波数ごとに予め求められた不快閾値の逆数を、前記複数種類のグローバルモーションの周波数成分に対して乗算することによって重み付けを施す構成であってもよい。   The weighting unit, for the frequency components of the plurality of types of global motion calculated by the Fourier transform unit, respectively, calculates the reciprocal number of the discomfort threshold obtained in advance for each frequency of the plurality of types of global motion. A configuration may be used in which weighting is performed by multiplying frequency components of various types of global motion.

前記不快度算出部は、前記非線形和を対数に変換したものを前記不快度として前記不快度補正部へ出力する構成であってもよい。   The discomfort degree calculation unit may output the non-linear sum converted to a logarithm as the discomfort degree to the discomfort degree correction unit.

この不快度推定装置では、等間隔に増加するような値として補正前の不快度を算出することができる。   In this discomfort degree estimation apparatus, the uncomfortable degree before correction can be calculated as a value that increases at equal intervals.

前記不快度補正部は、前回の前記不快度に0よりも大きく1よりも小さい定数を乗じ、今回の前記不快度に、定数が乗じられた前回の前記不快度を加算することによって、不快度を算出する構成であってもよい。   The discomfort level correcting unit multiplies the previous discomfort level by a constant larger than 0 and smaller than 1, and adds the previous discomfort level multiplied by the constant to the current discomfort level, thereby obtaining the discomfort level. May be configured to calculate.

また、本発明は、コンピュータを前記した不快度推定装置として機能させる不快度推定プログラムとしても具現化可能である。   The present invention can also be embodied as a discomfort level estimation program that causes a computer to function as the above-described discomfort level estimation device.

本発明によれば、グローバルモーションによって視聴者が感じる不快度をグローバルモーションに基づいて推定することができる。   According to the present invention, it is possible to estimate the degree of discomfort felt by the viewer by global motion based on the global motion.

本発明の実施形態に係る不快度推定装置を示すブロック図である。It is a block diagram which shows the discomfort degree estimation apparatus which concerns on embodiment of this invention. 離散フーリエ変換部によって算出される周波数成分の周波数と、当該周波数成分を算出するための時間窓の時間窓長(フレーム数)と、の関係を示す表である。It is a table | surface which shows the relationship between the frequency of the frequency component calculated by a discrete Fourier-transform part, and the time window length (frame number) of the time window for calculating the said frequency component. 不快閾値の逆数の周波数特性の一例を示すグラフである。It is a graph which shows an example of the frequency characteristic of the reciprocal number of a discomfort threshold. 図1の不快度算出部を示す詳細ブロック図である。It is a detailed block diagram which shows the discomfort degree calculation part of FIG.

以下、本発明の実施形態について、適宜図面を参照しながら説明する。同様の部分には同一符号を付し、重複する説明を省略する。本発明において、「グローバルモーション」は、映像の動きの中における画面全体の動き(画面動揺)を指し、映像の中で被写体が局所的に動く「ローカルモーション」とは区別される。   Hereinafter, embodiments of the present invention will be described with reference to the drawings as appropriate. Similar parts are denoted by the same reference numerals, and redundant description is omitted. In the present invention, “global motion” refers to movement of the entire screen (motion of the screen) in the motion of the video, and is distinguished from “local motion” in which the subject moves locally in the video.

図1に示すように、本発明の実施形態に係る不快度推定装置1は、機能部として、グローバルモーション検出部10と、離散フーリエ変換部20と、重み付け部30と、不快度算出部40と、不快度補正部50と、を備える。   As shown in FIG. 1, the discomfort level estimation apparatus 1 according to the embodiment of the present invention includes a global motion detection unit 10, a discrete Fourier transform unit 20, a weighting unit 30, and a discomfort level calculation unit 40 as functional units. The discomfort degree correction unit 50 is provided.

<グローバルモーション検出部>
グローバルモーション検出部10は、映像信号を取得し、取得された映像信号に基づいて、映像の時間的に隣接する画像間における複数種類のグローバルモーション(フレーム間の差分)を時系列に検出し、検出結果を離散フーリエ変換部20へ出力する。本実施形態において、グローバルモーション検出部10は、複数種類のグローバルモーションとして、上下移動に関するグローバルモーション(撮像装置の上下方向(チルト)移動量)、左右移動量に関するグローバルモーション(撮像装置の左右方向(パン)移動量)、回転に関するグローバルモーション(撮像装置の光学軸まわり(ロール)の回転量)、及び、ズームに関するグローバルモーション(撮像装置のズーム倍率の変化量)を検出する。
<Global motion detector>
The global motion detection unit 10 acquires a video signal, and based on the acquired video signal, detects a plurality of types of global motion (difference between frames) between temporally adjacent images of the video in time series, The detection result is output to the discrete Fourier transform unit 20. In this embodiment, the global motion detection unit 10 includes, as a plurality of types of global motion, a global motion related to vertical movement (vertical direction (tilt) movement amount of the imaging device) and a global motion related to horizontal movement amount (horizontal direction of the imaging device ( Pan) movement amount), global motion related to rotation (rotation amount around the optical axis (roll) of the image pickup apparatus), and global motion related to zoom (change amount of zoom magnification of the image pickup apparatus).

なお、上下移動に関するグローバルモーションに関しては、画像中の被写体が上向きに動いた場合をプラスとし、下向きに動いた場合をマイナスとする。また、左右移動に関するグローバルモーションに関しては、画像中の被写体が右向きに動いた場合をプラスとし、左向きに動いた場合をマイナスとする。また、回転に関するグローバルモーションに関しては、画像中の被写体が反時計回りに回転した場合をプラスとし、時計回りに回転した場合をマイナスとする。また、ズームに関するグローバルモーションに関しては、画像中の被写体のサイズが拡大した場合をプラスとし、縮小した場合をマイナスとする。   In addition, regarding the global motion related to the vertical movement, the case where the subject in the image moves upward is positive, and the case where the subject moves downward is negative. As for the global motion related to the left / right movement, the positive is when the subject in the image moves to the right, and the negative when the subject moves to the left. In addition, regarding global motion related to rotation, a positive value is obtained when the subject in the image rotates counterclockwise, and a negative value when the object rotates clockwise. As for the global motion related to zoom, a positive value is given when the size of the subject in the image is enlarged, and a negative value is given when the size is reduced.

また、グローバルモーション検出部10は、前記した複数種類のグローバルモーションを、視聴者が標準観視条件で映像を視聴した際の画面全体における平均移動視角量として正規化する構成であってもよい。ここで、標準観視条件とは、例えば、ハイビジョン映像の場合において、画面の高さの3倍の距離から視聴するような視聴条件のことを言う。かかる構成により、画素数(上下移動及び左右移動に関するグローバルモーション)、角度(回転に関するグローバルモーション)、サイズ比(ズームに関するグローバルモーション)と単位が異なっていた各グローバルモーションを同等に扱うことができる。   Further, the global motion detection unit 10 may be configured to normalize the plurality of types of global motion as the average moving viewing angle amount in the entire screen when the viewer views the video under the standard viewing condition. Here, the standard viewing condition refers to a viewing condition for viewing from a distance three times the height of the screen in the case of a high-definition video, for example. With such a configuration, it is possible to treat each global motion having different units with the same number of pixels (global motion related to vertical movement and horizontal movement), angle (global motion related to rotation), and size ratio (global motion related to zoom).

グローバルモーション検出部10は、一般的なグローバルモーションの検出手法として、画面を複数のブロックに分割し、各ブロックの画像のフレーム間での移動量をブロックマッチングによって算出し、ブロック位置を補正した後の対応自乗誤差が最小となる場合のグローバルモーションの値を算出する。かかる算出手法において、グローバルモーション検出部10は、RANSAC(RANdom SAmple Consensus)法、MPDE(MultiPartial Differential Equation)法等といった公知の手法を適宜採用することができる。   As a general global motion detection method, the global motion detection unit 10 divides the screen into a plurality of blocks, calculates the amount of movement of each block image between frames by block matching, and corrects the block position. The global motion value is calculated when the corresponding square error of is minimum. In such a calculation method, the global motion detection unit 10 can appropriately adopt a known method such as a RANSAC (RANdom SAmple Consensus) method, an MPDE (MultiPartial Differential Equation) method, or the like.

<離散フーリエ変換部>
離散フーリエ変換部20は、グローバルモーション検出部10によって検出された複数種類のグローバルモーションを取得し、取得された複数種類のグローバルモーションに対して、それぞれ、周波数が低いほど長い期間で離散フーリエ変換を施すことによって、複数種類のグローバルモーションの周波数成分(振幅の絶対値に比例する値)を算出し、算出結果を重み付け部30へ出力する。
<Discrete Fourier transform unit>
The discrete Fourier transform unit 20 acquires a plurality of types of global motion detected by the global motion detection unit 10, and performs discrete Fourier transform on the acquired plurality of types of global motion in a longer period as the frequency is lower. As a result, frequency components (values proportional to the absolute value of the amplitude) of a plurality of types of global motion are calculated, and the calculation result is output to the weighting unit 30.

本実施形態において、離散フーリエ変換部20は、1[秒]ごとに離散フーリエ変換を行う。これは、大半の視聴者(評価実験における平均的な評定者)において不快度が閾値に達するまでには早くても2[秒]を要することから、1[秒]間隔で不快度を推定すれば十分な精度が得られるためである。   In the present embodiment, the discrete Fourier transform unit 20 performs a discrete Fourier transform every 1 [second]. This is because most viewers (average raters in the evaluation experiment) take 2 [seconds] at the earliest to reach the threshold, and therefore, the discomfort can be estimated at intervals of 1 [second]. This is because sufficient accuracy can be obtained.

本実施形態において、離散フーリエ変換部20は、上下移動に関するグローバルモーションに対して離散フーリエ変換を施す上下移動用離散フーリエ変換部21と、左右移動に関するグローバルモーションに対して離散フーリエ変換を施す左右移動用離散フーリエ変換部22と、回転に関するグローバルモーションに対して離散フーリエ変換を施す回転用離散フーリエ変換部23と、ズームに関するグローバルモーションに対して離散フーリエ変換を施すズーム用離散フーリエ変換部24と、を備える。   In the present embodiment, the discrete Fourier transform unit 20 includes a discrete Fourier transform unit 21 that performs a discrete Fourier transform on a global motion related to vertical movement, and a left-right move that performs a discrete Fourier transform on a global motion related to a horizontal movement. A discrete Fourier transform unit 22 for rotation, a discrete discrete Fourier transform unit 23 for performing discrete Fourier transform on a global motion related to rotation, a discrete Fourier transform unit 24 for zooming performing discrete Fourier transform on a global motion related to zoom, Is provided.

図2は、離散フーリエ変換部によって算出される周波数成分の周波数と、当該周波数成分を算出するための時間窓の時間窓長(フレーム数)と、の関係を示す表である。図2に示すように、離散フーリエ変換部20は、グローバルモーションが揺れとしてではなく、ほぼ定速運動として認識される0.06[Hz]以下の周波数成分(時間窓長34.13[秒]の0〜2周期)を、周波数成分の算出の対象外としている。   FIG. 2 is a table showing the relationship between the frequency component frequency calculated by the discrete Fourier transform unit and the time window length (number of frames) of the time window for calculating the frequency component. As shown in FIG. 2, the discrete Fourier transform unit 20 has a frequency component of 0.06 [Hz] or less (time window length 34.13 [seconds]) that is recognized as a substantially constant speed motion, not as a global motion. 0 to 2 periods) is excluded from the calculation of frequency components.

また、グローバルモーションが不快なものとして認識されるには、3〜5周期を要することが評価実験によって明らかになっているので、離散フーリエ変換部20は、各時間窓において、最低周波数(1周期)の約3〜5倍の周波数におけるグローバルモーションの周波数成分を算出する。   Moreover, since it is clear from the evaluation experiment that it takes 3 to 5 cycles for the global motion to be recognized as unpleasant, the discrete Fourier transform unit 20 has the lowest frequency (1 cycle) in each time window. ) To calculate the frequency component of global motion at a frequency about 3 to 5 times that of).

すなわち、最低周波数の2倍の周波数(2周期)が0.06[Hz]よりも若干小さくなる時間窓を最長の時間窓とすればよいので、本実施形態では、離散フーリエ変換部20は、標準的なフレーム周波数である30[Hz]の映像において、最長の時間窓を1024フレーム(34.13[秒])に設定して離散フーリエ変換を行う。   That is, since the time window in which the frequency (two cycles) twice the lowest frequency is slightly smaller than 0.06 [Hz] may be used as the longest time window, in this embodiment, the discrete Fourier transform unit 20 is In an image of 30 [Hz] that is a standard frame frequency, the longest time window is set to 1024 frames (34.13 [seconds]), and discrete Fourier transform is performed.

そして、離散フーリエ変換部20は、時間窓長を1/2にしつつ、それぞれの時間窓長において最低周波数の約3〜5倍の周波数におけるグローバルモーションの周波数成分を算出する。なお、10[Hz]以上のグローバルモーションに関しては、画面の動揺ではなく二重像として認識されるため、離散フーリエ変換部20は、時間窓を16フレーム(0.53[秒])に設定して離散フーリエ変換を行う。   Then, the discrete Fourier transform unit 20 calculates the frequency component of global motion at a frequency of about 3 to 5 times the lowest frequency in each time window length while halving the time window length. In addition, since the global motion of 10 [Hz] or more is recognized as a double image rather than a screen shake, the discrete Fourier transform unit 20 sets the time window to 16 frames (0.53 [seconds]). To perform a discrete Fourier transform.

上下移動用離散フーリエ変換部21は、フレーム数を1024に設定して、0.088[Hz]、0.117[Hz]、0.146[Hz]の上下移動に関するグローバルモーションの周波数成分を算出し、フレーム数を512に設定して、0.176[Hz]、0.234[Hz]、0.293[Hz]の上下移動に関するグローバルモーションの周波数成分を算出し、フレーム数を256に設定して、0.352[Hz]、0.469[Hz]、0.586[Hz]の上下移動に関するグローバルモーションの周波数成分を算出し、フレーム数を128に設定して、0.703[Hz]、0.938[Hz]、1.172[Hz]の上下移動に関するグローバルモーションの周波数成分を算出し、フレーム数を64に設定して、1.406[Hz]、1.875[Hz]、2.344[Hz]の上下移動に関するグローバルモーションの周波数成分を算出し、フレーム数を32に設定して、2.813[Hz]、3.750[Hz]、4.688[Hz]の上下移動に関するグローバルモーションの周波数成分を算出し、フレーム数を16に設定して、5.625[Hz]、7.500[Hz]、9.375[Hz]、11.250[Hz]、13.125[Hz]、15.000[Hz]の上下移動に関するグローバルモーションの周波数成分を算出する。なお、左右移動用離散フーリエ変換部22、回転用離散フーリエ変換部23及びズーム用離散フーリエ変換部24も、上下移動用離散フーリエ変換部21と同様に、各グローバルモーションの周波数成分を算出する。   The discrete Fourier transform unit 21 for vertical movement sets the number of frames to 1024 and calculates the frequency component of global motion related to vertical movement of 0.088 [Hz], 0.117 [Hz], and 0.146 [Hz]. Then, the number of frames is set to 512, the frequency component of global motion related to vertical movement of 0.176 [Hz], 0.234 [Hz], and 0.293 [Hz] is calculated, and the number of frames is set to 256. Then, the frequency component of the global motion related to the vertical movement of 0.352 [Hz], 0.469 [Hz], 0.586 [Hz] is calculated, the number of frames is set to 128, and 0.703 [Hz ], 0.938 [Hz], 1.172 [Hz], the global motion frequency component for vertical movement is calculated, the number of frames is set to 64, and 1.40 The frequency component of the global motion related to the vertical movement of [Hz], 1.875 [Hz], 2.344 [Hz] is calculated, the number of frames is set to 32, and 2.813 [Hz], 3.750 [ The frequency component of the global motion related to the vertical movement of [Hz] and 4.688 [Hz] is calculated, the number of frames is set to 16, 5.625 [Hz], 7.500 [Hz], and 9.375 [Hz]. ] The frequency component of the global motion regarding the vertical movement of 11.250 [Hz], 13.125 [Hz], and 15.000 [Hz] is calculated. Note that the left-right moving discrete Fourier transform unit 22, the rotating discrete Fourier transform unit 23, and the zooming discrete Fourier transform unit 24 also calculate the frequency components of each global motion in the same manner as the up-and-down moving discrete Fourier transform unit 21.

本実施形態では、各時間窓において離散フーリエ変換が行われる周波数に重なりも抜けも無いので、離散フーリエ変換部20は、グローバルモーションの総エネルギーが保存された状態で周波数成分を算出することができる。したがって、本実施形態に係る不快度推定装置1は、ウェーブレット変換等を行うフィルタバンクを備える必要がない。   In the present embodiment, since there is no overlap or omission in the frequency at which the discrete Fourier transform is performed in each time window, the discrete Fourier transform unit 20 can calculate the frequency component in a state where the total energy of the global motion is stored. . Therefore, the discomfort degree estimation apparatus 1 according to the present embodiment does not need to include a filter bank that performs wavelet transform or the like.

また、離散フーリエ変換部20は、フレーム数(時間窓長)を繰り返し1/2にするのに都合の良い2のべき乗に設定しているが、離散フーリエ変換を施して算出される周波数成分は、各時間窓において3個(最短の時間窓においてのみ6つ)であるので、高速フーリエ変換(FFT:Fast Fourier Transform)を行う必要がなく、フレーム数を2のべき乗に設定しなくてもよい。   The discrete Fourier transform unit 20 sets the number of frames (time window length) to a power of 2 that is convenient for repeatedly halving, but the frequency component calculated by performing the discrete Fourier transform is Since there are three in each time window (six only in the shortest time window), there is no need to perform Fast Fourier Transform (FFT), and the number of frames need not be set to a power of two. .

<重み付け部>
重み付け部30は、離散フーリエ変換部20から出力された複数種類のグローバルモーションの周波数成分を取得し、取得された複数種類のグローバルモーションの周波数成分に対して、それぞれ、重み付けを施すものであり、本実施形態では、重み付けとして、周波数ごとに予め求められた不快閾値の逆数を、複数種類のグローバルモーションの周波数成分に対して乗算し、乗算結果を不快度算出部40へ出力する。本実施形態において、重み付け部30は、上下移動用重み付け部31と、左右移動用重み付け部32と、回転用重み付け部33と、ズーム用重み付け部34と、を備える。
<Weighting section>
The weighting unit 30 acquires a plurality of types of global motion frequency components output from the discrete Fourier transform unit 20, and applies weights to the acquired plurality of types of global motion frequency components, In the present embodiment, as the weighting, the reciprocal number of the discomfort threshold obtained in advance for each frequency is multiplied by the frequency components of a plurality of types of global motion, and the multiplication result is output to the discomfort degree calculation unit 40. In the present embodiment, the weighting unit 30 includes a vertical movement weighting unit 31, a horizontal movement weighting unit 32, a rotation weighting unit 33, and a zoom weighting unit 34.

図3は、不快閾値の逆数の周波数特性の一例を示すグラフである。ここで、不快閾値とは、各グローバルオーションに関して、周波数ごとに、視聴者が不快であると感じるグローバルモーションの周波数成分の閾値(不快閾値)を多数の評定者の幾何平均によって求めた値である。図3における縦軸(不快閾値の逆数)は、リニアスケールであり、単位としては、例えば、平均視覚速度の逆数である[sec/deg]等を採用することができる。また、不快閾値の逆数は、周波数間の値の比(相対値)が意味を有するものであって、その値自体に関しては、当業者によって適宜設定可能である。   FIG. 3 is a graph showing an example of the frequency characteristic of the reciprocal of the unpleasant threshold. Here, the unpleasant threshold value is a value obtained by calculating the frequency component threshold value (unpleasant threshold value) of the global motion that the viewer feels unpleasant for each frequency by the geometric average of many evaluators. . The vertical axis in FIG. 3 (the reciprocal of the unpleasant threshold) is a linear scale, and for example, [sec / deg] that is the reciprocal of the average visual speed can be adopted. Further, the reciprocal of the unpleasant threshold is meaningful by the ratio (relative value) of values between frequencies, and the value itself can be appropriately set by those skilled in the art.

上下移動用重み付け部31は、上下移動用離散フーリエ変換部21から出力された上下移動に関するグローバルモーションの周波数成分を取得し、取得された複数の周波数成分に重み付けを施し、重み付けの結果を不快度算出部40へ出力する。本実施形態では、0.088[Hz]の上下移動に関する周波数成分をa、0.117[Hz]の上下移動に関する周波数成分をa、・・・、15.000[Hz]の上下移動に関する周波数成分をa24とし、0.088[Hz]の上下移動に関する不快閾値の逆数をk、0.117[Hz]の上下移動に関する不快閾値の逆数をk、・・・、15.000[Hz]の上下移動に関する不快閾値の逆数をk24としたとき、上下移動用重み付け部31は、重み付けの結果として、k・a,k・a,・・・,k24・a24を不快度算出部40へ出力する。 The weighting unit for vertical movement 31 acquires the frequency component of the global motion related to the vertical movement output from the discrete Fourier transform unit 21 for vertical movement, weights the obtained plurality of frequency components, and sets the weighted result as the discomfort level. Output to the calculation unit 40. In the present embodiment, the frequency component related to the vertical movement of 0.088 [Hz] is a 1 , the frequency component related to the vertical movement of 0.117 [Hz] is the vertical movement of a 2 ,..., 15.000 [Hz]. 15. The frequency component for A is 24 , the reciprocal of the unpleasant threshold for up / down movement of 0.088 [Hz] is k 1 , the reciprocal of the unpleasant threshold for up / down movement of 0.117 [Hz] is k 2 ,. 000 when the reciprocal of the discomfort threshold for vertical movement of the [Hz] set to k 24, elevating weighting unit 31, as a result of the weighting, k 1 · a 1, k 2 · a 2, ···, k 24 Output a 24 to the discomfort level calculation unit 40.

左右移動用重み付け部32は、左右移動用離散フーリエ変換部22から出力された左右移動に関するグローバルモーションの周波数成分を取得し、取得された複数の周波数成分に重み付けを施し、重み付けの結果を不快度算出部40へ出力する。本実施形態では、0.088[Hz]の左右移動に関する周波数成分をa25、0.117[Hz]の左右移動に関する周波数成分をa26、・・・、15.000[Hz]の左右移動に関する周波数成分をa48とし、0.088[Hz]の左右移動に関する不快閾値の逆数をk25、0.117[Hz]の左右移動に関する不快閾値の逆数をk26、・・・、15.000[Hz]の左右移動に関する不快閾値の逆数をk48としたとき、左右移動用重み付け部32は、予め記憶された不快閾値の逆数(重み値)と周波数成分とを乗算することによってk25・a25,k26・a26,・・・,k48・a48を算出し、重み付けの結果として、かかる算出結果を不快度算出部40へ出力する。 The left / right moving weighting unit 32 acquires the frequency component of the global motion related to the left / right movement output from the left / right moving discrete Fourier transform unit 22, weights the acquired plurality of frequency components, and determines the weighted result as the discomfort level. Output to the calculation unit 40. In the present embodiment, the frequency component related to the horizontal movement of 0.088 [Hz] is a 25 , the frequency component related to the horizontal movement of 0.117 [Hz] is the horizontal movement of a 26 ,..., 15.000 [Hz]. 15. The frequency component related to A is 48 , the reciprocal of the discomfort threshold relating to 0.088 [Hz] left-right movement is k 25 , the reciprocal of the discomfort threshold relating to 0.117 [Hz] left-right movement is k 26 ,. 000 when the reciprocal of the discomfort threshold for lateral movement of the [Hz] set to k 48, lateral movement weighting unit 32, k 25 by multiplying the reciprocal of the pre-stored discomfort threshold value (weight value) and a frequency component A 25 , k 26 · a 26 ,..., K 48 · a 48 are calculated, and the calculation result is output to the discomfort degree calculation unit 40 as a weighting result.

回転用重み付け部33は、回転用離散フーリエ変換部23から出力された回転に関するグローバルモーションの周波数成分を取得し、取得された複数の周波数成分に重み付けを施し、重み付けの結果を不快度算出部40へ出力する。本実施形態では、0.088[Hz]の回転に関する周波数成分をa49、0.117[Hz]の回転に関する周波数成分をa50、・・・、15.000[Hz]の回転に関する周波数成分をa72とし、0.088[Hz]の回転に関する不快閾値の逆数をk49、0.117[Hz]の回転に関する不快閾値の逆数をk50、・・・、15.000[Hz]の回転に関する不快閾値の逆数をk72としたとき、回転用重み付け部33は、予め記憶された不快閾値の逆数(重み値)と周波数成分とを乗算することによってk49・a49,k50・a50,・・・,k72・a72を算出し、重み付けの結果として、かかる算出結果を不快度算出部40へ出力する。 The rotation weighting unit 33 acquires the frequency component of the global motion related to the rotation output from the rotation discrete Fourier transform unit 23, weights the plurality of acquired frequency components, and sets the weighted result as the discomfort degree calculation unit 40. Output to. In the present embodiment, a frequency component related to rotation of 0.088 [Hz] is a 49 , a frequency component related to rotation of 0.117 [Hz] is a frequency component related to rotation of a 50 ,..., 15.000 [Hz]. was a a 72, 0.088 k 49 the inverse of the discomfort threshold for rotation of [Hz], 0.117 [Hz] k 50 the inverse of the discomfort threshold for rotation, ..., of 15.000 [Hz] when the reciprocal of the discomfort threshold for rotating and k 72, rotation weighting unit 33 in advance reciprocal of the stored discomfort threshold (weight value) and k 49 · a 49 by multiplying the frequency components, k 50 · a 50 ,..., k 72 · a 72 are calculated, and the calculation result is output to the discomfort degree calculation unit 40 as a result of weighting.

ズーム用重み付け部34は、ズーム用離散フーリエ変換部24から出力されたズームに関するグローバルモーションの周波数成分を取得し、取得された複数の周波数成分に重み付けを施し、重み付けの結果を不快度算出部40へ出力する。本実施形態では、0.088[Hz]のズームに関する周波数成分をa73、0.117[Hz]のズームに関する周波数成分をa74、・・・、15.000[Hz]のズームに関する周波数成分をa96とし、0.088[Hz]のズームに関する不快閾値の逆数をk73、0.117[Hz]のズームに関する不快閾値の逆数をk74、・・・、15.000[Hz]のズームに関する不快閾値の逆数をk96としたとき、ズーム用重み付け部34は、予め記憶された不快閾値の逆数(重み値)と周波数成分とを乗算することによってk73・a73,k74・a74,・・・,k96・a96を算出し、重み付けの結果として、かかる算出結果を不快度算出部40へ出力する。 The zoom weighting unit 34 acquires the frequency component of the global motion related to the zoom output from the zoom discrete Fourier transform unit 24, weights the acquired plurality of frequency components, and sets the weighted result as the discomfort degree calculation unit 40. Output to. In the present embodiment, a frequency component related to a zoom of 0.088 [Hz] is a 73 , a frequency component related to a zoom of 0.117 [Hz] is a frequency component related to a zoom of a 74 ,. Is a 96 , the reciprocal of the discomfort threshold for zoom at 0.088 [Hz] is k 73 , and the reciprocal of the discomfort threshold for zoom at 0.117 [Hz] is k 74 ,..., 15.000 [Hz]. when the reciprocal of the discomfort threshold for zooming and k 96, zooming weighting unit 34, k 73 · a 73, k 74 · by multiplying the reciprocal of the pre-stored discomfort threshold value (weight value) and a frequency component a 74 ,..., k 96 · a 96 are calculated, and the calculation result is output to the discomfort degree calculation unit 40 as a weighted result.

<不快度算出部>
不快度算出部40は、重み付け部30によって重み付けされた周波数成分k・a〜k96・a96を取得し、かかる重み付けされた周波数成分k・a〜k96・a96に基づいて、不快度cを算出する。本実施形態において、不快度算出部40は、不快度cを1秒ごとに繰り返し算出して不快度補正部50へ出力する。図4は、図1の不快度算出部を示す詳細ブロック図である。図4に示すように、不快度算出部40は、成分加算部41と、対数変換部42と、所定時間内最大値検出部43と、を備える。
<Discomfort calculation unit>
Discomfort calculating unit 40 obtains the frequency component k 1 · a 1 ~k 96 · a 96 weighted by the weighting unit 30, based on such a weighted frequency components k 1 · a 1 ~k 96 · a 96 Then, the discomfort degree c is calculated. In the present embodiment, the discomfort degree calculation unit 40 repeatedly calculates the discomfort degree c every second and outputs it to the discomfort degree correction unit 50. FIG. 4 is a detailed block diagram illustrating the discomfort level calculation unit of FIG. As shown in FIG. 4, the discomfort degree calculation unit 40 includes a component addition unit 41, a logarithmic conversion unit 42, and a maximum value detection unit 43 within a predetermined time.

成分加算部41は、重み付け部30から出力された、重み付けされた周波数成分k1・a1〜k96・a96を取得し、これらを非線形加算することによって、非線形和bを算出し、算出結果を成分加算部42へ出力する。
b={Σ(k・aγ1/γ
The component addition unit 41 obtains the weighted frequency components k1 · a1 to k96 · a96 output from the weighting unit 30, calculates the nonlinear sum b by performing nonlinear addition, and adds the calculated result to the component To the unit 42.
b = {Σ (k n · a n) γ} 1 / γ

すなわち、成分加算部41は、重み付けされた周波数成分k・a(n=1,2,・・・96)をそれぞれγ乗したものを全て加算し、その1/γ乗を算出することよって、非線形和bを得る。ここで、非線形和bは、γ=1の場合には、重み付けされた周波数成分の和(振幅和)となり、γ=2の場合には、グローバルモーションのエネルギーの和となる。主観評価実験の各評定者について調べると、γが約2である場合に後記する不快値dが主観評価実験に適合する。しかし、評定による不快値が同じ2つの種類のグローバルモーションを合成した場合に、評定者ごとに見ると平均値のまわりにそれぞれ30%程度の標準偏差が存在するため、電力和の場合よりも不快値が4%ほど大きくなる。この場合には、γは約1.8となる。したがって、γの値は、
1.8≦γ≦2.0
を満たすように設定されることが望ましい。
That is, the component adding unit 41 adds all the weighted frequency components k n · a n (n = 1, 2,... 96) to the γ power and calculates the 1 / γ power. Therefore, the nonlinear sum b is obtained. Here, the non-linear sum b is a sum of weighted frequency components (amplitude sum) when γ = 1, and is a sum of global motion energy when γ = 2. When examining each evaluator of the subjective evaluation experiment, when γ is about 2, an unpleasant value d described later is suitable for the subjective evaluation experiment. However, when two types of global motion with the same unpleasant value by rating are synthesized, there is a standard deviation of about 30% around the average value when viewed by each evaluator. The value increases by about 4%. In this case, γ is about 1.8. Therefore, the value of γ is
1.8 ≦ γ ≦ 2.0
It is desirable to set so as to satisfy.

対数変換部42は、成分加算部41から出力された非線形和bを対数に変換することによって、不快度cを算出し、算出結果を所定時間内最大値検出部43へ出力する。
c=log
The logarithmic conversion unit 42 calculates the discomfort degree c by converting the non-linear sum b output from the component addition unit 41 into a logarithm, and outputs the calculation result to the maximum value detection unit 43 within a predetermined time.
c = log X b

このように、非線形和bを対数変換することによって不快度cを算出するのは、主観評価実験における評定に用いた「不快ではない」、「やや不快」、「不快」、「非常に不快」、「耐え難く不快」のカテゴリーの境界における、非線形和bを対数変換した値が、ほぼ等間隔になることに基づく措置である。ここで、非線形和bが約6.5倍になると、不快度が1カテゴリー上昇するので、対数の底Xは、
X≒6.5
に設定されることが望ましい。
なお、前記した不快閾値は、主観評価実験のカテゴリーにおける、「やや不快」と「不快」との境界を示す閾値である。
In this way, the degree of discomfort c is calculated by logarithmically transforming the non-linear sum b. This is a measure based on the fact that the values obtained by logarithmically transforming the nonlinear sum b at the boundary of the category of “unbearable and uncomfortable” are almost equally spaced. Here, when the nonlinear sum b is about 6.5 times, the discomfort level increases by one category.
X≈6.5
It is desirable to be set to.
The discomfort threshold described above is a threshold indicating a boundary between “somewhat discomfort” and “discomfort” in the category of the subjective evaluation experiment.

また、非線形和bには、ある値以下では長時間蓄積しても不快閾値に達しない臨界値があり、非線形和bが臨界値に等しい場合に不快度cが0となるように対数の底Xを設定し、対数変換部42は、不快度cがマイナスの場合には、当該不快度cを出力しない(又は、不快度cとして0を出力する)ように構成される。かかる臨界値は、主観評価実験によると、対数変換前の値として、不快閾値の約1/50の値に設定されることが望ましい。   Further, the nonlinear sum b has a critical value that does not reach the discomfort threshold even if it accumulates for a long time below a certain value. When the nonlinear sum b is equal to the critical value, the base of the logarithm is such that the discomfort degree c becomes zero. X is set, and the logarithmic conversion unit 42 is configured not to output the discomfort level c (or output 0 as the discomfort level c) when the discomfort level c is negative. Such a critical value is desirably set to a value about 1/50 of the discomfort threshold as a value before logarithmic conversion according to a subjective evaluation experiment.

所定時間内最大値検出部43は、離散フーリエ変換部20による離散フーリエ変換処理時点の直前の所定期間(例えば、5〜7[秒])の間における不快度cの最大値を検出し、当該最大値を不快度cとして不快度補正部50の加算部51へ出力する。これは、間欠的にグローバルモーションが無い期間が3秒程度である映像に対する評定不快度が、連続して同じ振幅のグローバルモーションが発生している映像に対する評定不快度とほぼ同様となること(すなわち、グローバルモーションが3秒程度無かった場合でも、その前後にグローバルモーションが発生していれば、ヒトは先のグローバルモーションによる影響を受けた不快度を感じてしまう)に基づく措置である。   The maximum value detection unit 43 within a predetermined time detects the maximum value of the discomfort degree c during a predetermined period (for example, 5 to 7 [seconds]) immediately before the discrete Fourier transform processing time point by the discrete Fourier transform unit 20, The maximum value is output to the adding unit 51 of the discomfort level correcting unit 50 as the discomfort level c. This means that the rating discomfort for a video in which the period of no global motion is intermittent is about 3 seconds is almost the same as the rating discomfort for a video in which global motion of the same amplitude is continuously generated (that is, Even if there is no global motion for about 3 seconds, if global motion has occurred before and after that, humans will feel discomfort influenced by the previous global motion).

本実施形態において、所定時間内最大値検出部43は、記憶部43aと、最大値検出部43bと、を備える。   In the present embodiment, the maximum value detection unit 43 within the predetermined time includes a storage unit 43a and a maximum value detection unit 43b.

記憶部43aは、対数変換部42から出力された不快度cを、直近の所定個数(例えば、5〜7個)だけ記憶し、記憶された所定個数の不快度cを最大値検出部43bへ出力するシフトレジスタである。   The storage unit 43a stores the most recent predetermined number (for example, 5 to 7) of the discomfort degree c output from the logarithmic conversion unit 42, and stores the predetermined number of discomfort levels c to the maximum value detection unit 43b. This is a shift register to output.

最大値検出部43bは、記憶部43aから出力された所定個数の不快度cを比較し、最大の不快度cを検出して加算部51へ出力する。   The maximum value detection unit 43b compares a predetermined number of discomfort levels c output from the storage unit 43a, detects the maximum discomfort level c, and outputs it to the addition unit 51.

<不快度補正部>
図1に戻り、不快度補正部50は、前回の(1秒前の)不快度c(以下、前回の不快度cとする)が大きいほど今回の不快度c(以下、今回の不快度cm+1とする)が大きくなるように今回の不快度cm+1を補正することによって、補正後の不快度dを算出する。本実施形態において、不快度算出部50は、加算部51と、記憶部52と、乗算部53と、を備える。
<Discomfort correction unit>
Returning to Figure 1, discomfort degree correcting unit 50, the previous (one second ago) unpleasantness c (hereinafter, the previous and unpleasantness c m) is too large this unpleasantness c (hereinafter, this unpleasantness The uncomfortable degree d after correction is calculated by correcting the current uncomfortable degree cm + 1 so that (c m + 1 ) becomes large. In the present embodiment, the discomfort degree calculation unit 50 includes an addition unit 51, a storage unit 52, and a multiplication unit 53.

加算部51は、不快度算出部40から出力された前回の不快度cを取得し、取得された前回の不快度cを記憶部52へ出力する。 The adding unit 51 acquires the previous discomfort level cm output from the discomfort level calculation unit 40, and outputs the acquired previous discomfort level cm to the storage unit 52.

記憶部52は、加算部51から出力された前回の不快度cを、1個だけ記憶し、記憶された前回の不快度cを乗算部53へ出力するシフトレジスタである。 The storage unit 52 is a shift register that stores only one previous discomfort degree cm output from the addition unit 51 and outputs the stored previous discomfort degree cm to the multiplication unit 53.

乗算部53は、記憶部52から出力された前回の不快度cを取得し、取得された前回の不快度cに定数p(0<p<1)を乗じ、算出結果のp・cを加算部51へ出力する。この定数pは、グローバルモーションが長時間継続する場合における不快度の蓄積効果を表すものであり、不快度の推定誤差を最小とするには、p=0.005程度とすることが望ましい。 The multiplication unit 53 acquires the previous discomfort degree cm output from the storage unit 52, multiplies the acquired previous discomfort degree cm by a constant p (0 <p <1), and calculates the calculation result p · c m is output to the adder 51. The constant p represents an effect of accumulating discomfort when the global motion continues for a long time. In order to minimize an estimation error of discomfort, it is desirable that p is about 0.005.

そして、加算部51は、不快度算出部40から出力された今回の不快度cm+1を取得し、今回の不快度cm+1に、定数が乗じられた前回の不快度p・cを加算することによって、補正後の不快度dを算出し、算出された不快度dをディスプレイ、スピーカ等の外部装置(不快度dを利用者へ通知する通知部)へ出力するとともに、新たな前回の不快度として記憶部52へ出力する。
d=cm+1+p・c
かかる加算部51は、離散フーリエ変換部20による離散フーリエ変換のタイミングである1[秒]ごとに補正後の不快度dを算出して出力する。
Then, the adding unit 51 obtains the current discomfort level cm + 1 that is output from the discomfort level calculation unit 40, and adds the previous discomfort level p · cm multiplied by a constant to the current discomfort level cm + 1. As a result, the corrected discomfort level d is calculated, and the calculated discomfort level d is output to an external device (notification unit for notifying the user of the discomfort level d) such as a display or a speaker, and a new previous discomfort level is output. It outputs to the storage unit 52 as the degree.
d = c m + 1 + p · c m
The adding unit 51 calculates and outputs the corrected discomfort degree d every 1 [second] which is the timing of the discrete Fourier transform by the discrete Fourier transform unit 20.

本発明の実施形態に係る不快度推定装置1は、複数種類のグローバルモーションの周波数成分が時間的に変化するような画面動揺を伴う映像に対しても、主観評価実験における平均値に近い不快度dを推定することができる。また、推定された不快度dは、映像コンテンツ製作者が画面動揺を低減するような映像修正を施す際に、映像修正の度合いの指標として使用可能である。   The discomfort degree estimation apparatus 1 according to the embodiment of the present invention has a discomfort degree that is close to the average value in the subjective evaluation experiment, even for images with screen fluctuations in which frequency components of multiple types of global motion change over time. d can be estimated. Further, the estimated discomfort degree d can be used as an index of the degree of video correction when the video content producer performs video correction to reduce screen shaking.

すなわち、不快度推定装置1は、映像コンテンツ制作者によって制作段階で用いられる場合には、映像の良否の判定、映像に含まれる画面動揺をどの程度まで低減すべきかの目標設定等に好適な不快度dを映像コンテンツ製作者に提示することができるので、制作に要する時間、労力及びコストの削減が図られるだけでなく、供給される映像コンテンツの安全性及び快適性も高められる。
また、不快度推定装置1は、映像の視聴者側で用いられる場合には、画面動揺に関して安全、快適であることを保証せずに制作、流通された映像に対して、視聴前又は視聴中の表示直前に不快度dを推定してディスプレイ又はスピーカへ出力することによって、視聴時に警告を発することができるので、映像酔いによる健康被害及び不快感の誘発を防止することが可能になる。
That is, the discomfort level estimation apparatus 1 is suitable for determining whether a video is good or not, and for setting a target for how much the screen shaking included in the video should be reduced, when used by the video content creator at the production stage. Since the degree d can be presented to the video content producer, not only can the time, labor, and cost required for production be reduced, but also the safety and comfort of the supplied video content can be improved.
Further, when used on the viewer side of a video, the discomfort level estimation device 1 is pre-viewing or watching a video that has been produced and distributed without guaranteeing that the screen is safe and comfortable with respect to shaking. Since the degree of discomfort d is estimated and output to the display or the speaker immediately before the display of the image, a warning can be issued at the time of viewing, so that it is possible to prevent health damage and inducing discomfort due to video sickness.

以上、本発明の実施形態について実施形態を参照して説明したが、本発明は前記実施形態に限定されず、本発明の要旨を逸脱しない範囲で適宜設計変更可能である。例えば、本発明は、コンピュータを前記不快度推定装置1として機能させる不快度推定プログラムとして具現化することも可能である。また、不快度決定部40によって決定された不快度dが閾値を超えた場合に、映像が不快なグローバルモーションを含んでいると判定し、判定結果をディスプレイ等の外部装置へ出力する構成であってもよい。また、対数変換部42及び所定時間内最大値検出部43は省略可能である。   As mentioned above, although embodiment of this invention was described with reference to embodiment, this invention is not limited to the said embodiment, A design change is possible suitably in the range which does not deviate from the summary of this invention. For example, the present invention can also be embodied as a discomfort level estimation program that causes a computer to function as the discomfort level estimation device 1. In addition, when the discomfort degree d determined by the discomfort degree determination unit 40 exceeds a threshold, it is determined that the video includes an unpleasant global motion, and the determination result is output to an external device such as a display. May be. Further, the logarithmic conversion unit 42 and the maximum value detection unit 43 within a predetermined time can be omitted.

1 不快度推定装置
10 グローバルモーション検出部
20 離散フーリエ変換部
30 重み付け部
40 不快度算出部
50 不快度補正部
DESCRIPTION OF SYMBOLS 1 Discomfort degree estimation apparatus 10 Global motion detection part 20 Discrete Fourier transform part 30 Weighting part 40 Discomfort degree calculation part 50 Discomfort degree correction part

Claims (6)

映像における画面動揺を示すグローバルモーションに基づいて不快度を推定する不快度推定装置であって、
前記映像の時間的に隣接する画像間における複数種類のグローバルモーションを検出するグローバルモーション検出部と、
検出された前記複数種類のグローバルモーションに対して、それぞれ、周波数が低いほど長い期間で離散フーリエ変換を施すことによって、前記複数種類のグローバルモーションの周波数成分を算出する離散フーリエ変換部と、
算出された前記複数種類のグローバルモーションの周波数成分に対して、それぞれ重み付けを施す重み付け部と、
重み付けされた前記複数種類のグローバルモーションの周波数成分の非線形和を算出することによって、前記不快度を算出する不快度算出部と、
前回の前記不快度が大きいほど今回の前記不快度が大きくなるように今回の前記不快度を補正する不快度補正部と、
を備えることを特徴とする不快度推定装置。
A discomfort level estimation device that estimates a discomfort level based on global motion indicating a screen shake in a video,
A global motion detector that detects a plurality of types of global motion between temporally adjacent images of the video;
A discrete Fourier transform unit that calculates a frequency component of the plurality of types of global motion by performing discrete Fourier transform in a longer period as the frequency is lower for each of the detected types of global motion,
A weighting unit that weights the calculated frequency components of the plurality of types of global motion,
A discomfort degree calculating unit that calculates the discomfort degree by calculating a nonlinear sum of weighted frequency components of the plurality of types of global motion;
A discomfort level correction unit that corrects the current discomfort level so that the current discomfort level increases as the previous discomfort level increases;
A discomfort degree estimation device comprising:
前記グローバルモーション検出部は、前記複数種類のグローバルモーションとして、上下移動、左右移動、回転、及び、ズームに関するグローバルモーションを検出する
ことを特徴とする請求項1に記載の不快度推定装置。
The said global motion detection part detects the global motion regarding a vertical movement, a left-right movement, a rotation, and a zoom as said multiple types of global motion. The discomfort degree estimation apparatus of Claim 1 characterized by the above-mentioned.
前記重み付け部は、前記フーリエ変換部によって算出された前記複数種類のグローバルモーションの周波数成分に対して、それぞれ、前記複数種類のグローバルモーションの周波数ごとに予め求められた不快閾値の逆数を、前記複数種類のグローバルモーションの周波数成分に対して乗算することによって重み付けを施す
ことを特徴とする請求項1又は請求項2に記載の不快度推定装置。
The weighting unit, for the frequency components of the plurality of types of global motion calculated by the Fourier transform unit, respectively, calculates the reciprocal number of the discomfort threshold obtained in advance for each frequency of the plurality of types of global motion. 3. The discomfort degree estimation apparatus according to claim 1, wherein weighting is performed by multiplying frequency components of different types of global motion.
前記不快度算出部は、前記非線形和を対数に変換したものを前記不快度として前記不快度補正部へ出力する
ことを特徴とする請求項1から請求項3のいずれか一項に記載の不快度推定装置。
The unpleasantness degree calculation unit outputs the conversion of the non-linear sum into a logarithm as the unpleasantness level to the unpleasantness degree correction unit. The unpleasantness according to any one of claims 1 to 3 Degree estimation device.
前記不快度補正部は、前回の前記不快度に0よりも大きく1よりも小さい定数を乗じ、今回の前記不快度に、定数が乗じられた前回の前記不快度を加算することによって、補正後の今回の前記不快度を得る
ことを特徴とする請求項1から請求項4のいずれか一項に記載の不快度推定装置。
The discomfort level correcting unit multiplies the previous discomfort level by a constant larger than 0 and smaller than 1, and adds the previous discomfort level multiplied by the constant to the current discomfort level after correction. The discomfort degree estimation apparatus according to any one of claims 1 to 4, wherein the discomfort degree of the current time is obtained.
映像における画面動揺を示すグローバルモーションに基づいて不快度を推定する不快度推定プログラムであって、
コンピュータを、
前記映像の時間的に隣接する画像間における複数種類のグローバルモーションを検出するグローバルモーション検出部、
検出された前記複数種類のグローバルモーションに対して、それぞれ、周波数が低いほど長い期間で離散フーリエ変換を施すことによって、前記複数種類のグローバルモーションの周波数成分を算出する離散フーリエ変換部、
算出された前記複数種類のグローバルモーションの周波数成分に対して、それぞれ重み付けを施す重み付け部、
重み付けされた前記複数種類のグローバルモーションの周波数成分の非線形和を算出することによって、前記不快度を算出する不快度算出部、及び、
前回の前記不快度が大きいほど今回の前記不快度が大きくなるように今回の前記不快度を補正する不快度補正部、
として機能させることを特徴とする不快度推定プログラム。
A discomfort degree estimation program that estimates discomfort based on global motion indicating screen sway in a video,
Computer
A global motion detector for detecting a plurality of types of global motion between temporally adjacent images of the video,
A discrete Fourier transform unit for calculating the frequency components of the plurality of types of global motion by performing discrete Fourier transform on the detected types of global motions in a longer period as the frequency is lower,
A weighting unit for weighting the calculated frequency components of the plurality of types of global motion,
A discomfort degree calculating unit for calculating the discomfort degree by calculating a nonlinear sum of weighted frequency components of the plurality of types of global motion; and
A discomfort level correction unit that corrects the current discomfort level so that the current discomfort level increases as the previous discomfort level increases;
An unpleasantness degree estimation program characterized by functioning as
JP2010110401A 2010-05-12 2010-05-12 Discomfort degree estimation apparatus and discomfort degree estimation program Expired - Fee Related JP5391145B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2010110401A JP5391145B2 (en) 2010-05-12 2010-05-12 Discomfort degree estimation apparatus and discomfort degree estimation program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2010110401A JP5391145B2 (en) 2010-05-12 2010-05-12 Discomfort degree estimation apparatus and discomfort degree estimation program

Publications (2)

Publication Number Publication Date
JP2011238120A true JP2011238120A (en) 2011-11-24
JP5391145B2 JP5391145B2 (en) 2014-01-15

Family

ID=45326012

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2010110401A Expired - Fee Related JP5391145B2 (en) 2010-05-12 2010-05-12 Discomfort degree estimation apparatus and discomfort degree estimation program

Country Status (1)

Country Link
JP (1) JP5391145B2 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013021455A (en) * 2011-07-08 2013-01-31 Nippon Hoso Kyokai <Nhk> Discomfort degree estimation apparatus and discomfort degree estimation program
JP2016162219A (en) * 2015-03-02 2016-09-05 日本放送協会 Discomfort degree estimation apparatus and discomfort degree estimation program
JP2021101324A (en) * 2019-12-24 2021-07-08 本田技研工業株式会社 Gesture detection device, gesture detection method, and program
JP2022533282A (en) * 2019-10-18 2022-07-22 グーグル エルエルシー Multi-tier scalable media analytics

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007108458A1 (en) * 2006-03-23 2007-09-27 Matsushita Electric Industrial Co., Ltd. Content imaging apparatus
WO2008032739A1 (en) * 2006-09-12 2008-03-20 Panasonic Corporation Content imaging device
JP2008299241A (en) * 2007-06-04 2008-12-11 Sharp Corp Image processing apparatus and display
JP2010004113A (en) * 2008-06-18 2010-01-07 Olympus Imaging Corp Imaging apparatus, imaging method, image processing device and image processing method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007108458A1 (en) * 2006-03-23 2007-09-27 Matsushita Electric Industrial Co., Ltd. Content imaging apparatus
WO2008032739A1 (en) * 2006-09-12 2008-03-20 Panasonic Corporation Content imaging device
JP2008299241A (en) * 2007-06-04 2008-12-11 Sharp Corp Image processing apparatus and display
JP2010004113A (en) * 2008-06-18 2010-01-07 Olympus Imaging Corp Imaging apparatus, imaging method, image processing device and image processing method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JPN7013003438; 財団法人 機械システム振興協会: 映像酔いガイドライン検証システムの実用化に関するフィージビリティスタディ報告書(要旨) , 200803, pp.17-33 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013021455A (en) * 2011-07-08 2013-01-31 Nippon Hoso Kyokai <Nhk> Discomfort degree estimation apparatus and discomfort degree estimation program
JP2016162219A (en) * 2015-03-02 2016-09-05 日本放送協会 Discomfort degree estimation apparatus and discomfort degree estimation program
JP2022533282A (en) * 2019-10-18 2022-07-22 グーグル エルエルシー Multi-tier scalable media analytics
JP7234356B2 (en) 2019-10-18 2023-03-07 グーグル エルエルシー Multi-tier scalable media analytics
JP2021101324A (en) * 2019-12-24 2021-07-08 本田技研工業株式会社 Gesture detection device, gesture detection method, and program
JP7175257B2 (en) 2019-12-24 2022-11-18 本田技研工業株式会社 GESTURE DETECTION DEVICE, GESTURE DETECTION METHOD, AND PROGRAM

Also Published As

Publication number Publication date
JP5391145B2 (en) 2014-01-15

Similar Documents

Publication Publication Date Title
JP5680524B2 (en) Image processing device
JP5391145B2 (en) Discomfort degree estimation apparatus and discomfort degree estimation program
US20100328532A1 (en) Image generating device, static text detecting device and method thereof
US20150206289A1 (en) Joint Video Deblurring and Stabilization
JP2012513719A (en) Generating image scaling curves
JP5530399B2 (en) Image processing method, image processing apparatus, and image processing program
US20110293172A1 (en) Image processing apparatus, image processing method, and image display apparatus
JP5410366B2 (en) Discomfort degree estimation apparatus and discomfort degree estimation program
JP2010268525A5 (en)
US9602762B2 (en) Fallback in frame rate conversion system
JP5395737B2 (en) Discomfort degree estimation apparatus and discomfort degree estimation program
JP2012014627A5 (en)
JP5727886B2 (en) Discomfort degree estimation apparatus and discomfort degree estimation program
US20200410645A1 (en) Image processing apparatus, image processing method, and storage medium
JP5997575B2 (en) Discomfort degree estimation apparatus and discomfort degree estimation program
JP6076110B2 (en) Discomfort degree estimation apparatus and discomfort degree estimation program
JP6487718B2 (en) Discomfort degree estimation apparatus and discomfort degree estimation program
JP2014121079A5 (en)
JP2015106192A (en) Discomfort degree estimation device and discomfort degree estimation program
JP2016111473A (en) Objective image quality evaluation device and program
JP2017029401A (en) Determination device for change in pupil diameter by light stimulus, determination method and program
JP2020021168A (en) Discomfort degree estimation device and discomfort degree estimation program
JP5364927B2 (en) Motion compensation apparatus and method
JP6196541B2 (en) Motion perception amount estimation device, discomfort level estimation device, and motion perception amount estimation program
JP6867822B2 (en) Discomfort estimation device and discomfort estimation program

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20130104

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20130906

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: 20130917

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20131011

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees