JP2010258914A - Prominent area image generating method, prominent area image generating device, program, and recording medium - Google Patents

Prominent area image generating method, prominent area image generating device, program, and recording medium Download PDF

Info

Publication number
JP2010258914A
JP2010258914A JP2009108474A JP2009108474A JP2010258914A JP 2010258914 A JP2010258914 A JP 2010258914A JP 2009108474 A JP2009108474 A JP 2009108474A JP 2009108474 A JP2009108474 A JP 2009108474A JP 2010258914 A JP2010258914 A JP 2010258914A
Authority
JP
Japan
Prior art keywords
image
saliency area
feature amount
saliency
likelihood
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
JP2009108474A
Other languages
Japanese (ja)
Other versions
JP5235770B2 (en
Inventor
Shogo Kimura
昭悟 木村
Junji Yamato
淳司 大和
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.)
Nippon Telegraph and Telephone Corp
Original Assignee
Nippon Telegraph and Telephone 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 Nippon Telegraph and Telephone Corp filed Critical Nippon Telegraph and Telephone Corp
Priority to JP2009108474A priority Critical patent/JP5235770B2/en
Publication of JP2010258914A publication Critical patent/JP2010258914A/en
Application granted granted Critical
Publication of JP5235770B2 publication Critical patent/JP5235770B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/25Determination of region of interest [ROI] or a volume of interest [VOI]

Abstract

<P>PROBLEM TO BE SOLVED: To provide a technology for enabling an area to be divided even when no information on an object area and a background area is given in advance. <P>SOLUTION: An area image generating device 1000 includes: a degree-of-interest image extraction unit 1 which extracts a degree-of-interest image from input image, a prominent area prior probability image extraction unit 2 which extracts a prominent area prior probability image showing a probability that each position of an input image as each of frames constituting the input image is in a prominent area; a feature quantity likelihood calculation unit 3 which calculates a feature quantity likelihood showing a likelihood of the image feature quantity that the input image is included in the prominent area and an area other than the prominent area; a prominent area image extraction unit 4 which extracts a prominent area image showing a prominent area of the input image from the input image, the prominent area prior probability image, and the feature quantity likelihood; and a prominent area image generation unit 5 which generates the prominent area image from the each prominent area image obtained by executing the above processes on each input image. <P>COPYRIGHT: (C)2011,JPO&INPIT

Description

本発明は、顕著領域映像生成方法、顕著領域映像生成装置、プログラムおよび記録媒体に関する。特に、本発明は、人間の視覚特性を利用した信号処理によって顕著領域(物体領域)、背景領域(非顕著領域)に関する事前情報を算出し、算出した事前情報を利用して、入力映像から顕著領域と非顕著領域とを高精度に分割する、顕著領域映像生成方法、顕著領域映像生成装置、プログラムおよび記録媒体に関する。   The present invention relates to a saliency image generation method, a saliency image generation device, a program, and a recording medium. In particular, the present invention calculates prior information about a saliency area (object area) and a background area (non-salience area) by signal processing using human visual characteristics, and uses the calculated prior information to make a significant difference from the input video. The present invention relates to a saliency area image generation method, a saliency area image generation apparatus, a program, and a recording medium that divide an area and a non-salience area with high accuracy.

画像や映像から、人物・動物・物体など興味の対象である領域(以下、「物体領域」という)を、背景などの物体領域以外の領域(以下、「背景領域」という)と区別して抽出する領域分割技術は、クロマキーなどを利用しない自由な画像映像合成、背景領域の変動に頑健な物体認識・画像映像検索、領域の重要性に応じてビットレートを調整可能な画像映像符号化など、応用範囲の広い重要な技術である。領域分割技術として、画像領域分割をある種の統計モデルに対する事後確率最大化問題として定式化し、事後確率最大化問題を統計モデルと等価なグラフの最小カットを求めることによって解決する方法が知られている(例えば、非特許文献1参照)。また、非特許文献1に記載の方法を映像信号に拡張し、映像信号の時間的な連続性を利用してグラフの最小カットを効率的に求めることで高速に映像の領域分割を実現する方法も知られている(例えば、非特許文献2参照)。   Extracting regions of interest (hereinafter referred to as “object regions”) such as people, animals, and objects from images and video separately from regions other than object regions such as backgrounds (hereinafter referred to as “background regions”) Area segmentation technology can be applied to free image synthesis without using chroma keys, object recognition and image retrieval robust to changes in the background area, and image video coding that can adjust the bit rate according to the importance of the area. It is an important technology with a wide range. As a region segmentation technique, a method is known in which image region segmentation is formulated as a posterior probability maximization problem for a certain statistical model, and the posterior probability maximization problem is solved by finding the minimum cut of the graph equivalent to the statistical model. (For example, refer nonpatent literature 1). In addition, the method described in Non-Patent Document 1 is extended to a video signal, and a video segmentation is realized at high speed by efficiently obtaining a minimum cut of a graph using temporal continuity of the video signal. Is also known (see, for example, Non-Patent Document 2).

Y.Boykov and G.F.Lea,“Graph cuts and efficient N-D image segmentation,” International Journal of Computer Vision,Vol.70,No.2,pp.109-131,2006.Y. Boykov and G.F.Lea, “Graph cuts and efficient N-D image segmentation,” International Journal of Computer Vision, Vol. 70, No. 2, pp. 109-131, 2006. P.Kohli and P. Torr,“Dynamic graph cuts for efficient inference in Markov random fields,”IEEE Transactions on Pattern Analysis and Machine Intelligence,Vol.29,No.12,pp.2079-2088,2007.P. Kohli and P. Torr, “Dynamic graph cuts for efficient inference in Markov random fields,” IEEE Transactions on Pattern Analysis and Machine Intelligence, Vol. 29, No. 12, pp. 2079-2088, 2007.

領域分割技術の応用では、物体領域や背景領域に関する事前情報(手がかり)が与えられない状況においても正しく物体領域と背景領域とを分割できる手法、換言すれば、人手を介せず領域分割を含めて全ての手順を自動で実行できる手法が望まれている。領域分割技術の応用範囲が拡大するからである。しかしながら、非特許文献1、2に記載の方法は、物体領域・背景領域に関する事前情報が部分的に手動で与えられることを想定しているため、物体領域・背景領域に関する事前情報が全く与えられない場合には利用することができないという問題がある。即ち、領域分割技術の応用範囲を著しく制限しているという問題がある。   In the application of region segmentation technology, a method that can correctly segment the object region and the background region even in the situation where prior information (cue) about the object region and the background region is not given, in other words, including region segmentation without human intervention. Therefore, a method that can automatically execute all the procedures is desired. This is because the application range of the area division technique is expanded. However, since the methods described in Non-Patent Documents 1 and 2 assume that prior information on the object region / background region is partially manually given, no prior information on the object region / background region is given at all. If not, there is a problem that it cannot be used. That is, there is a problem that the application range of the region division technique is remarkably limited.

本発明は、上述した課題に鑑みてなされたものであって、物体領域、背景領域に関する事前情報が全く与えられない場合においても領域分割を可能にするための技術を提供することを目的とする。   The present invention has been made in view of the above-described problems, and an object of the present invention is to provide a technique for enabling region division even when prior information regarding an object region and a background region is not given at all. .

上記問題を解決するために、本発明の一態様である顕著領域映像生成方法は、入力映像から、人間が注意を向けやすい度合いである注目度を示す注目度映像を抽出する注目度映像抽出過程と、入力映像を構成する各フレームである入力画像の各位置が顕著領域である確率を示す顕著領域事前確率画像を抽出する顕著領域事前確率画像抽出過程と、前記入力画像の顕著領域および顕著領域外の領域にそれぞれ含まれる画像特徴量の尤度を示す特徴量尤度を算出する特徴量尤度算出過程と、前記入力画像、前記顕著領域事前確率画像および前記特徴量尤度から、前記入力画像の顕著領域を示す顕著領域画像を抽出する顕著領域画像抽出部と、前記各入力画像に対し、前記注目度映像抽出過程、前記顕著領域事前確率画像抽出過程、前記特徴量尤度算出過程および前記顕著領域画像抽出過程を実行して得られる前記各顕著領域画像から前記顕著領域映像を生成する顕著領域映像生成過程とを有し、前記顕著領域事前確率画像抽出過程は、一の前記入力画像の各位置が顕著領域である確率を示す顕著領域事前確率画像を、前記注目度映像抽出過程によって抽出された前記注目度映像内の前記入力画像に対応する画像である注目度画像および前記顕著領域画像に基づいて抽出し、前記特徴量尤度算出過程は、前記特徴量尤度を、前記入力画像、前記注目度画像、前記顕著領域事前確率画像、前記顕著領域画像および前回迄に算出した前記特徴量尤度の少なくとも1つに基づいて算出することを特徴とする。   In order to solve the above problem, a remarkable area image generation method according to an aspect of the present invention includes an attention level image extraction process for extracting an attention level image indicating a degree of attention that is a degree to which a person is likely to pay attention from an input image. A saliency area prior probability image extraction process for extracting a saliency area prior probability image indicating the probability that each position of the input image, which is each frame constituting the input video, is a saliency area, and the saliency area and saliency area of the input image From the feature amount likelihood calculation process for calculating the feature amount likelihood indicating the likelihood of the image feature amount included in each of the outer regions, the input image, the saliency region prior probability image, and the feature amount likelihood, the input A saliency area image extracting unit for extracting a saliency area image indicating a saliency area of the image; and for each of the input images, the attention degree video extraction process, the saliency area prior probability image extraction process, and the feature amount likelihood. And a saliency area image generation process for generating the saliency area image from the saliency area images obtained by executing the saliency area image extraction process and the saliency area prior probability image extraction process. A remarkable area prior probability image indicating the probability that each position of the input image is a remarkable area, an attention degree image that is an image corresponding to the input image in the attention degree video extracted by the attention degree video extraction process, and The feature amount likelihood calculation process extracts the feature amount likelihood from the input image, the attention degree image, the saliency area prior probability image, the saliency area image, and the previous time. It is calculated based on at least one of the calculated feature amount likelihoods.

上記顕著領域映像生成方法において、前記注目度映像抽出過程は、前記入力画像に基づいて、該入力画像の中で顕著な特性を持つ空間領域を表示した画像である基礎注目度画像を算出する基礎注目度画像抽出過程と、現在の前記入力画像の各位置における顕著性を確率的な表現を用いて表示した画像である確率的基礎注目度画像を、前記基礎注目度画像抽出過程によって算出された基礎注目度画像と、前回の前記入力画像から該確率的基礎注目度画像抽出過程によって算出された前記確率的基礎注目度画像と、逐次更新され、視線位置推定に用いる第1の母数である確率的基礎注目度母数とに基づいて算出する確率的基礎注目度画像抽出過程と、現在の前記入力画像における前記視線位置確率密度映像のフレームである視線位置確率密度画像を、前記確率的基礎注目度画像抽出過程によって算出された確率的基礎注目度画像と、前回の前記入力画像から該視線位置確率密度画像抽出過程によって算出された前記視線位置確率密度画像と、逐次更新され、視線位置推定に用いる第2の母数である視線位置確率密度母数とに基づいて算出する視線位置確率密度画像抽出過程と、前記基礎注目度画像抽出過程と、前記確率的基礎注目度画像抽出過程と、前記視線位置確率密度画像抽出過程とを、それぞれの前記入力画像に対して順に繰り返し行うことにより算出される時系列の前記視線位置確率密度画像を前記視線位置確率密度映像として出力する視線位置確率密度映像出力過程とを有し、前記視線位置確率密度画像抽出過程は、視線移動の大きさを制御する確率変数である視線移動状態変数を、前回の前記入力画像から該視線位置確率密度画像抽出過程によって算出された前記視線位置確率密度画像と、前回の前記入力画像から該視線移動状態変数更新過程によって算出された視線移動状態変数と、前記視線位置確率密度母数とに基づいて更新し、該視線移動状態変数の集合である視線移動状態変数集合を出力する視線移動状態変数更新過程と、視線移動を考慮した代表的な視線位置を示す代表視線位置の集合である代表視線位置集合を、前記確率的基礎注目度画像抽出過程によって算出された確率的基礎注目度画像と、前回の前記入力画像から該代表視線位置更新過程によって更新された代表視線位置集合と、前記視線移動状態変数集合と、前記視線位置確率密度母数とに基づいて更新する代表視線位置更新過程と、それぞれの前記代表視線位置に関連付けられた重みからなる代表視線位置重み係数の集合である代表視線位置重み係数集合を、前記確率的基礎注目度画像抽出過程によって算出された確率的基礎注目度画像と、前記代表視線位置更新過程によって更新された代表視線位置集合と、前記視線移動状態変数更新過程から出力された視線移動状態変数集合と、前記視線位置確率密度母数とに基づいて算出する代表視線位置重み係数算出過程と、前記代表視線位置更新過程によって更新された代表視線位置集合と、代表視線位置重み係数算出過程によって算出された代表視線位置重み係数集合とに基づいて、前記視線位置確率密度画像を算出する視線位置確率密度画像出力過程とを有し、前記代表視線位置集合と、前記代表視線位置重み係数集合とを含む前記視線位置確率密度画像を算出するようにしてもよい。   In the saliency area video generation method, the attention level image extraction step is a basis for calculating a basic attention level image that is an image displaying a spatial area having a remarkable characteristic in the input image based on the input image. Attention level image extraction process and a stochastic basic attention level image that is an image displaying prominence at each position of the current input image using a probabilistic expression was calculated by the basic attention level image extraction process The basic attention degree image, the probabilistic basic attention degree image calculated by the probabilistic basic attention degree image extraction process from the previous input image, and a first parameter that is sequentially updated and used for gaze position estimation. A stochastic basic attention degree image extraction process that is calculated based on a stochastic basic attention degree parameter, and a gaze position probability density image that is a frame of the gaze position probability density image in the current input image The probabilistic basic attention degree image calculated by the probabilistic basic attention level image extraction process, and the gaze position probability density image calculated by the gaze position probability density image extraction process from the previous input image, and sequentially updated A gaze position probability density image extraction process calculated based on a gaze position probability density parameter that is a second parameter used for gaze position estimation, the basic attention degree image extraction process, and the stochastic basic attention degree The time-series gaze position probability density image calculated by sequentially repeating the image extraction process and the gaze position probability density image extraction process for each of the input images is output as the gaze position probability density image. A line-of-sight position probability density image output process, and the line-of-sight position probability density image extraction process is a line-of-sight movement state that is a random variable that controls the magnitude of line-of-sight movement The line-of-sight position probability density image calculated by the line-of-sight position probability density image extraction process from the previous input image, and the line-of-sight movement state variable calculated by the line-of-sight movement state variable update process from the previous input image And a line-of-sight movement state variable update process for outputting a line-of-sight movement state variable set that is a set of the line-of-sight movement state variables, and a representative line of sight in consideration of line-of-sight movement. A representative line-of-sight position set, which is a set of representative line-of-sight positions indicating a position, is obtained by a probabilistic basic attention level image calculated by the probabilistic basic attention level image extraction process and a representative gaze position update process from the previous input image. Representative visual line position update process for updating based on the updated representative visual line position set, the visual line movement state variable set, and the visual line position probability density parameter, respectively. A representative gaze position weight coefficient set, which is a set of representative gaze position weight coefficients composed of weights associated with the representative gaze position, and a probabilistic basic attention image calculated by the probabilistic basic attention image extraction process; The representative gaze position calculated based on the representative gaze position set updated by the representative gaze position update process, the gaze movement state variable set output from the gaze movement state variable update process, and the gaze position probability density parameter. Based on the weight coefficient calculation process, the representative eye position position set updated by the representative eye position update process, and the representative eye position weight coefficient set calculated by the representative eye position weight coefficient calculation process, the eye position probability density image A visual line position probability density image output process for calculating the representative visual line position set and the representative visual line position weighting coefficient set It may be calculated the linear position probability density image.

上記顕著領域映像生成方法において、前記顕著領域事前確率画像抽出過程は、前記注目度画像のみを用いて前記顕著領域事前確率画像を生成する顕著領域事前確率画像生成過程と、前記顕著領域画像を用いて前記顕著領域事前確率画像生成過程によって生成された前記顕著領域事前確率画像を更新する顕著領域事前確率画像更新過程とから構成されるようにしてもよい。   In the saliency area image generation method, the saliency area prior probability image extraction process uses the saliency area prior probability image generation process for generating the saliency area prior probability image using only the attention level image, and the saliency area image. And a saliency area prior probability image update process for updating the saliency area prior probability image generated by the saliency area prior probability image generation process.

上記顕著領域映像生成方法において、前記特徴量尤度算出過程は、顕著領域に含まれる画像特徴量の尤度を示す顕著領域特徴量尤度を、前記入力画像、前記顕著領域事前確率画像、前記顕著領域画像および前回迄に算出した前記顕著領域特徴量尤度のうち少なくとも1つに基づいて算出する顕著領域特徴量尤度算出過程と、顕著領域外の領域に含まれる画像特徴量の尤度を示す非顕著領域特徴量尤度を、前記入力画像、前記顕著領域事前確率画像、前記顕著領域画像および前回迄に算出した前記非顕著領域特徴量尤度のうち少なくとも1つに基づいて算出する非顕著領域特徴量尤度算出過程と、前記顕著領域特徴量尤度および前記非顕著領域特徴量尤度を加算して特徴量尤度として出力する特徴量尤度出力過程とから構成されるようにしてもよい。   In the saliency area video generation method, the feature amount likelihood calculation step includes the saliency area feature amount likelihood indicating the likelihood of the image feature amount included in the saliency area, the input image, the saliency area prior probability image, A saliency area feature amount likelihood calculation process based on at least one of the saliency area image and the saliency area feature amount likelihood calculated up to the previous time, and the likelihood of the image feature amount included in the area outside the saliency area Is calculated based on at least one of the input image, the saliency area prior probability image, the saliency area image, and the non-salience area feature amount likelihood calculated until the previous time. A non-significant area feature quantity likelihood calculation process and a feature quantity likelihood output process of adding the saliency area feature quantity likelihood and the non-salience area feature quantity likelihood to output as a feature quantity likelihood. even if There.

上記顕著領域映像生成方法において、前記顕著領域特徴量尤度算出過程は、前記入力画像、前記顕著領域事前確率画像および前記顕著領域画像に基づいて、前記顕著領域特徴量尤度を生成する顕著領域特徴量尤度生成過程と、前記顕著領域特徴量尤度生成過程によって生成された前記顕著領域特徴量尤度を更新する顕著領域特徴量尤度更新過程とから構成され、非顕著領域特徴量尤度算出過程は、前記入力画像、前記顕著領域事前確率画像および前記顕著領域画像に基づいて、前記非顕著領域特徴量尤度を生成する非顕著領域特徴量尤度生成過程と、前記非顕著領域特徴量尤度生成過程によって生成された前記非顕著領域特徴量尤度を更新する非顕著領域特徴量尤度更新過程とから構成され、前記顕著領域特徴量尤度更新過程は、前記入力画像、前記顕著領域画像および前回迄に更新した更新後の前記顕著領域特徴量尤度のうち少なくとも1つに基づいて前記顕著領域特徴量尤度を更新し、前記非顕著領域特徴量尤度更新過程は、前記入力画像、非顕著領域画像および前回迄に更新した更新後の前記非顕著領域特徴量尤度のうち少なくとも1つに基づいて前記非顕著領域特徴量尤度を更新するようにしてもよい。   In the saliency area video generation method, the saliency area feature amount likelihood calculating step includes generating a saliency area feature amount likelihood based on the input image, the saliency area prior probability image, and the saliency area image. A feature amount likelihood generation process and a saliency area feature amount likelihood update process for updating the saliency area feature amount likelihood generated by the saliency area feature amount likelihood generation process. The degree calculation process includes: a non-salience area feature amount likelihood generation process for generating the non-salience area feature amount likelihood based on the input image, the saliency area prior probability image, and the saliency area image; A non-salience area feature amount likelihood update process for updating the non-salience area feature amount likelihood generated by the feature amount likelihood generation process, and the saliency area feature amount likelihood update process includes the input image Updating the saliency area feature quantity likelihood based on at least one of the saliency area image and the updated saliency area feature quantity likelihood updated up to the previous time, and updating the non-salience area feature quantity likelihood May update the non-significant region feature value likelihood based on at least one of the input image, the non-significant region image, and the updated non-significant region feature value likelihood updated so far. Good.

上記顕著領域映像生成方法は、前記入力画像を異なる解像度によってそれぞれ平滑化した複数の平滑化画像からなる平滑化画像群を生成する平滑化画像群生成過程と、前記平準化画像群に対し、前記顕著領域事前確率画像抽出過程、前記特徴量尤度算出過程、前記顕著領域画像抽出過程を実行し、前記入力画像の前記顕著領域画像を確定する顕著領域画像確定過程とを更に有し、前記特徴量尤度算出過程および顕著領域画像抽出過程は、前記入力画像に代えて前記平滑化画像を用い、前記顕著領域映像生成過程は、前記各入力画像に対し、前記注目度映像抽出過程、前記顕著領域事前確率画像抽出過程、前記特徴量尤度算出過程、前記顕著領域画像抽出過程、前記平滑化画像群生成過程および顕著領域画像確定過程を実行して得られる前記顕著領域画像から前記顕著領域映像を生成するようにしてもよい。   The saliency area image generation method includes a smoothed image group generation process for generating a smoothed image group including a plurality of smoothed images obtained by smoothing the input image at different resolutions, and the leveled image group, A saliency area prior probability image extraction process, the feature amount likelihood calculation process, and the saliency area image extraction process, and further including a saliency area image determination process for determining the saliency area image of the input image, The quantity likelihood calculation process and the saliency area image extraction process use the smoothed image instead of the input image, and the saliency area image generation process performs the attention level image extraction process and the saliency for each input image. The saliency obtained by executing an area prior probability image extraction process, the feature likelihood calculation process, the saliency area image extraction process, the smoothed image group generation process, and the saliency area image determination process. It may generate the salient region image from the area image.

上記問題を解決するために、本発明の他の態様である顕著領域映像生成装置は、入力映像から、人間が注意を向けやすい度合いである注目度を示す注目度映像を抽出する注目度映像抽出部と、入力映像を構成する各フレームである前記入力画像の各位置が顕著領域である確率を示す顕著領域事前確率画像を抽出する顕著領域事前確率画像抽出部と、前記入力画像の顕著領域および顕著領域外の領域にそれぞれ含まれる画像特徴量の尤度を示す特徴量尤度を算出する特徴量尤度算出部と、前記入力画像、前記顕著領域事前確率画像および前記特徴量尤度から、前記入力画像の顕著領域を示す顕著領域画像を抽出する顕著領域画像抽出部と、前記各入力画像に対し、前記注目度映像抽出部、前記顕著領域事前確率画像抽出部、前記特徴量尤度算出部および前記顕著領域画像抽出部を実行して得られる前記各顕著領域画像から前記顕著領域映像を生成する顕著領域映像生成部とを備え、前記顕著領域事前確率画像抽出部は、一の前記入力画像の各位置が顕著領域である確率を示す顕著領域事前確率画像を、前記注目度映像抽出部によって抽出された前記注目度映像内の前記入力画像に対応する画像である注目度画像および前記顕著領域画像に基づいて抽出し、前記特徴量尤度算出部は、前記特徴量尤度を、前記入力画像、前記注目度画像、前記顕著領域事前確率画像、前記顕著領域画像および前回迄に算出した前記特徴量尤度の少なくとも1つに基づいて算出することを特徴とする。   In order to solve the above-described problem, a remarkable area image generation device according to another aspect of the present invention extracts an attention level image that extracts an attention level image indicating a degree of attention that is a degree to which a person is likely to pay attention from an input image. A saliency area prior probability image extraction unit that extracts a saliency area prior probability image indicating a probability that each position of the input image that is each frame constituting the input video is a saliency area, and a saliency area of the input image and From the feature amount likelihood calculating unit that calculates the likelihood of the image feature amount included in each of the regions outside the saliency region, the input image, the saliency region prior probability image, and the feature amount likelihood, A saliency area image extraction unit that extracts a saliency area image indicating a saliency area of the input image; and, for each of the input images, the attention level video extraction unit, the saliency area prior probability image extraction unit, and the feature amount likelihood calculation. And a saliency area image generation unit that generates the saliency area image from each saliency area image obtained by executing the saliency area image extraction unit, and the saliency area prior probability image extraction unit includes the one input image A noticeable area prior probability image indicating a probability that each position of the noticeable area is a noticeable area, an attention degree image that is an image corresponding to the input image in the attention degree video extracted by the attention degree video extraction unit, and the salient area The feature amount likelihood calculating unit extracts the feature amount likelihood from the input image, the attention degree image, the saliency area prior probability image, the saliency area image, and the previous time. The calculation is based on at least one of the feature amount likelihoods.

上記顕著領域映像生成装置において、前記注目度映像抽出部は、前記入力画像に基づいて、該入力画像の中で顕著な特性を持つ空間領域を表示した画像である基礎注目度画像を算出する基礎注目度画像抽出部と、現在の前記入力画像の各位置における顕著性を確率的な表現を用いて表示した画像である確率的基礎注目度画像を、前記基礎注目度画像抽出部によって算出された基礎注目度画像と、前回の前記入力画像から該確率的基礎注目度画像抽出部によって算出された前記確率的基礎注目度画像と、逐次更新され、視線位置推定に用いる第1の母数である確率的基礎注目度母数とに基づいて算出する確率的基礎注目度画像抽出部と、現在の前記入力画像における前記視線位置確率密度映像のフレームである視線位置確率密度画像を、前記確率的基礎注目度画像抽出部によって算出された確率的基礎注目度画像と、前回の前記入力画像から該視線位置確率密度画像抽出部によって算出された前記視線位置確率密度画像と、逐次更新され、視線位置推定に用いる第2の母数である視線位置確率密度母数とに基づいて算出する視線位置確率密度画像抽出部と、前記基礎注目度画像抽出部と、前記確率的基礎注目度画像抽出部と、前記視線位置確率密度画像抽出部とを、それぞれの前記入力画像に対して順に繰り返し行うことにより算出される時系列の前記視線位置確率密度画像を前記視線位置確率密度映像として出力する視線位置確率密度映像出力部とを備え、前記視線位置確率密度画像抽出部は、視線移動の大きさを制御する確率変数である視線移動状態変数を、前回の前記入力画像から該視線位置確率密度画像抽出部によって算出された前記視線位置確率密度画像と、前回の前記入力画像から該視線移動状態変数更新部によって算出された視線移動状態変数と、前記視線位置確率密度母数とに基づいて更新し、該視線移動状態変数の集合である視線移動状態変数集合を出力する視線移動状態変数更新部と、視線移動を考慮した代表的な視線位置を示す代表視線位置の集合である代表視線位置集合を、前記確率的基礎注目度画像抽出部によって算出された確率的基礎注目度画像と、前回の前記入力画像から該代表視線位置更新部によって更新された代表視線位置集合と、前記視線移動状態変数集合と、前記視線位置確率密度母数とに基づいて更新する代表視線位置更新部と、それぞれの前記代表視線位置に関連付けられた重みからなる代表視線位置重み係数の集合である代表視線位置重み係数集合を、前記確率的基礎注目度画像抽出部によって算出された確率的基礎注目度画像と、前記代表視線位置更新部によって更新された代表視線位置集合と、前記視線移動状態変数更新部から出力された視線移動状態変数集合と、前記視線位置確率密度母数とに基づいて算出する代表視線位置重み係数算出部と、前記代表視線位置更新部によって更新された代表視線位置集合と、代表視線位置重み係数算出部によって算出された代表視線位置重み係数集合とに基づいて、前記視線位置確率密度画像を算出する視線位置確率密度画像出力部とを有し、前記代表視線位置集合と、前記代表視線位置重み係数集合とを含む前記視線位置確率密度画像を算出するようにしてもよい。   In the saliency area video generation device, the attention level video extraction unit calculates a basic attention level image that is an image displaying a spatial area having a remarkable characteristic in the input image based on the input image. The basic attention degree image extraction unit calculates a probabilistic basic attention degree image that is an image displaying the saliency at each position of the current input image using a probabilistic expression. The basic attention image, the probabilistic basic attention image calculated by the probabilistic basic attention image extraction unit from the previous input image, and a first parameter that is sequentially updated and used for gaze position estimation. A probabilistic basic attention degree image extraction unit that is calculated based on a probabilistic basic attention degree parameter, and a gaze position probability density image that is a frame of the gaze position probability density image in the current input image; The probabilistic basic attention level image calculated by the static basic attention level image extraction unit and the gaze position probability density image calculated by the gaze position probability density image extraction unit from the previous input image are sequentially updated, and the line of sight A gaze position probability density image extraction unit that is calculated based on a gaze position probability density parameter that is a second parameter used for position estimation, the basic attention level image extraction unit, and the probabilistic basic attention level image extraction unit And a line-of-sight position probability density image extracting unit that outputs the line-of-sight position probability density image in time series calculated as the line-of-sight position probability density image as the line-of-sight position probability density image. A probability density image output unit, and the line-of-sight position probability density image extraction unit converts a line-of-sight movement state variable, which is a random variable that controls the amount of line-of-sight movement, into the previous input image. The line-of-sight position probability density image calculated by the line-of-sight position probability density image extraction unit, the line-of-sight movement state variable calculated by the line-of-sight movement state variable update unit from the previous input image, and the line-of-sight position probability density mother A line-of-sight movement state variable update unit that outputs a line-of-sight movement state variable set that is a set of the line-of-sight movement state variables, and a set of representative line-of-sight positions indicating representative line-of-sight positions that take line-of-sight movement into account A representative gaze position set calculated by the probabilistic basic attention level image extraction unit, and a representative gaze position set updated by the representative gaze position update unit from the previous input image, , A representative line-of-sight position update unit for updating based on the line-of-sight movement state variable set and the line-of-sight position probability density parameter, and weights associated with the representative line-of-sight positions The representative gaze position weighting coefficient set, which is a set of representative gaze position weighting coefficients, is updated by the probabilistic basic attention degree image calculated by the probabilistic basic attention degree image extracting unit and the representative gaze position updating unit. A representative gaze position weight coefficient calculating unit that calculates a representative gaze position set, a gaze movement state variable set output from the gaze movement state variable update unit, and the gaze position probability density parameter; and the representative gaze position A gaze position probability density image output unit that calculates the gaze position probability density image based on the representative gaze position weight set updated by the update unit and the representative gaze position weight coefficient set calculated by the representative gaze position weight coefficient calculation unit The visual line position probability density image including the representative visual line position set and the representative visual line position weighting coefficient set may be calculated.

上記顕著領域映像生成装置において、前記顕著領域事前確率画像抽出部は、前記注目度画像のみを用いて前記顕著領域事前確率画像を生成する顕著領域事前確率画像生成部と、 前記顕著領域画像を用いて前記顕著領域事前確率画像生成部によって生成された前記顕著領域事前確率画像を更新する顕著領域事前確率画像更新部とから構成されるようにしてもよい。   In the saliency area video generation device, the saliency area prior probability image extraction unit uses the saliency area prior probability image generation part that generates the saliency area prior probability image using only the attention level image, and the saliency area image. And a saliency area prior probability image update unit that updates the saliency area prior probability image generated by the saliency area prior probability image generation unit.

上記顕著領域映像生成装置において、前記特徴量尤度算出部は、顕著領域に含まれる画像特徴量の尤度を示す顕著領域特徴量尤度を、前記入力画像、前記顕著領域事前確率画像、前記顕著領域画像および前回迄に算出した前記顕著領域特徴量尤度のうち少なくとも1つに基づいて算出する顕著領域特徴量尤度算出部と、顕著領域外の領域に含まれる画像特徴量の尤度を示す非顕著領域特徴量尤度を、前記入力画像、前記顕著領域事前確率画像、前記顕著領域画像および前回迄に算出した前記非顕著領域特徴量尤度のうち少なくとも1つに基づいて算出する非顕著領域特徴量尤度算出部と、前記顕著領域特徴量尤度および前記非顕著領域特徴量尤度を加算して特徴量尤度として出力する特徴量尤度出力部とから構成されるようにしてもよい。   In the saliency area video generation device, the feature amount likelihood calculating unit calculates the saliency area feature amount likelihood indicating the likelihood of the image feature amount included in the saliency area, the input image, the saliency area prior probability image, the The saliency area feature amount likelihood calculation unit that calculates based on at least one of the saliency area image and the saliency area feature amount likelihood calculated up to the previous time, and the likelihood of the image feature amount included in the area outside the saliency area Is calculated based on at least one of the input image, the saliency area prior probability image, the saliency area image, and the non-salience area feature amount likelihood calculated until the previous time. A non-significant region feature amount likelihood calculating unit and a feature amount likelihood output unit that adds the saliency region feature amount likelihood and the non-significant region feature amount likelihood and outputs the result as a feature amount likelihood. It may be.

上記顕著領域映像生成装置において、前記顕著領域特徴量尤度算出部は、前記入力画像、前記顕著領域事前確率画像および前記顕著領域画像に基づいて、前記顕著領域特徴量尤度を生成する顕著領域特徴量尤度生成部と、前記顕著領域特徴量尤度生成部によって生成された前記顕著領域特徴量尤度を更新する顕著領域特徴量尤度更新部とから構成され、前記非顕著領域特徴量尤度算出部は、前記入力画像、前記顕著領域事前確率画像および前記顕著領域画像に基づいて、前記非顕著領域特徴量尤度を生成する非顕著領域特徴量尤度生成部と、前記非顕著領域特徴量尤度生成部によって生成された前記非顕著領域特徴量尤度を更新する非顕著領域特徴量尤度更新部とから構成され、前記顕著領域特徴量尤度更新部は、前記入力画像、前記顕著領域画像および前回迄に更新した更新後の前記顕著領域特徴量尤度のうち少なくとも1つに基づいて前記顕著領域特徴量尤度を更新し、前記非顕著領域特徴量尤度更新部は、前記入力画像、非顕著領域画像および前回迄に更新した更新後の前記非顕著領域特徴量尤度のうち少なくとも1つに基づいて前記非顕著領域特徴量尤度を更新するようにしてもよい。   In the saliency area video generation device, the saliency area feature amount likelihood calculating unit generates the saliency area feature amount likelihood based on the input image, the saliency area prior probability image, and the saliency area image. The non-significant area feature quantity includes a feature quantity likelihood generation section and a saliency area feature quantity likelihood update section that updates the saliency area feature quantity likelihood generated by the saliency area feature quantity likelihood generation section. The likelihood calculation unit includes a non-salience region feature amount likelihood generation unit that generates the non-salience region feature amount likelihood based on the input image, the saliency region prior probability image, and the saliency region image; A non-significant region feature amount likelihood update unit that updates the non-salience region feature amount likelihood generated by the region feature amount likelihood generation unit, and the saliency region feature amount likelihood update unit includes the input image The prominent territory The saliency area feature quantity likelihood is updated based on at least one of the image and the updated saliency area feature quantity likelihood updated up to the previous time, and the non-salience area feature quantity likelihood update unit is configured to input the input The non-salience area feature amount likelihood may be updated based on at least one of the image, the non-salience area image, and the updated non-salience area feature amount likelihood updated so far.

上記顕著領域映像生成装置は、前記入力画像を異なる解像度によってそれぞれ平滑化した複数の平滑化画像からなる平滑化画像群を生成する平滑化画像群生成部と、前記平準化画像群に対し、前記顕著領域事前確率画像抽出部、前記特徴量尤度算出部、前記顕著領域画像抽出部の処理を実行し、前記入力画像の前記顕著領域画像を確定する顕著領域画像確定部とを更に備え、前記特徴量尤度算出部および顕著領域画像抽出部は、前記入力画像に代えて前記平滑化画像を用い、前記顕著領域映像生成部は、前記各入力画像に対し、前記注目度映像抽出部、前記顕著領域事前確率画像抽出部、前記特徴量尤度算出部、前記顕著領域画像抽出部、前記平滑化画像群生成部および前記顕著領域画像確定部の各処理を実行して得られる前記顕著領域画像から前記顕著領域映像を生成するようにしてもよい。   The saliency area video generation device includes a smoothed image group generation unit that generates a smoothed image group including a plurality of smoothed images obtained by smoothing the input image with different resolutions, and the leveled image group, A saliency area prior probability image extraction unit, the feature amount likelihood calculation unit, and a saliency area image extraction unit, and a saliency area image determination unit that determines the saliency area image of the input image. The feature amount likelihood calculation unit and the saliency area image extraction unit use the smoothed image instead of the input image, and the saliency area image generation unit performs the attention level image extraction unit, The saliency area image obtained by executing each process of the saliency area prior probability image extraction unit, the feature amount likelihood calculation unit, the saliency area image extraction unit, the smoothed image group generation unit, and the saliency area image determination unit It may generate an al the salient region image.

上記問題を解決するために、本発明の他の態様であるプログラムは、入力映像から、人間が注意を向けやすい度合いである注目度を示す注目度映像を抽出する注目度映像抽出ステップと、入力映像を構成する各フレームである前記入力画像の各位置が顕著領域である確率を示す顕著領域事前確率画像を抽出する顕著領域事前確率画像抽出ステップと、前記入力画像の顕著領域および顕著領域外の領域にそれぞれ含まれる画像特徴量の尤度を示す特徴量尤度を算出する特徴量尤度算出ステップと、前記入力画像、前記顕著領域事前確率画像および前記特徴量尤度から、前記入力画像の顕著領域を示す顕著領域画像を抽出する顕著領域画像抽出ステップと、前記各入力画像に対し、前記注目度映像抽出ステップ、前記顕著領域事前確率画像抽出ステップ、前記特徴量尤度算出ステップおよび前記顕著領域画像抽出ステップを実行して得られる前記各顕著領域画像から前記顕著領域映像を生成する顕著領域映像生成ステップとをコンピュータに実行させるプログラムであって、前記顕著領域事前確率画像抽出ステップは、一の前記入力画像の各位置が顕著領域である確率を示す顕著領域事前確率画像を、前記注目度映像抽出ステップによって抽出された前記注目度映像内の前記入力画像に対応する画像である注目度画像および前記顕著領域画像に基づいて抽出し、前記特徴量尤度算出ステップは、前記特徴量尤度を、前記入力画像、前記注目度画像、前記顕著領域事前確率画像、前記顕著領域画像および前回迄に算出した前記特徴量尤度の少なくとも1つに基づいて算出することを特徴とする。   In order to solve the above problem, a program according to another aspect of the present invention includes an attention level video extraction step of extracting an attention level video indicating a level of attention that is a degree to which a person is likely to pay attention from an input video, and an input A saliency area prior probability image extracting step for extracting a saliency area prior probability image indicating a probability that each position of the input image which is each frame constituting the video is a saliency area; The feature amount likelihood calculating step for calculating the feature amount likelihood indicating the likelihood of the image feature amount included in each region, and the input image, the saliency region prior probability image, and the feature amount likelihood A saliency area image extraction step for extracting a saliency area image indicating a saliency area; and, for each of the input images, the attention degree video extraction step and the saliency area prior probability image extraction stage. And a saliency area image generation step for generating the saliency area image from the saliency area images obtained by executing the feature amount likelihood calculation step and the saliency area image extraction step. The saliency area prior probability image extraction step includes a saliency area prior probability image indicating the probability that each position of the one input image is a saliency area in the attention level video extracted by the attention level video extraction step. Based on the attention level image that is an image corresponding to the input image and the saliency area image, and the feature amount likelihood calculating step includes calculating the feature amount likelihood as the input image, the attention level image, It is calculated based on at least one of the saliency area prior probability image, the saliency area image, and the feature amount likelihood calculated until the previous time. That.

上記問題を解決するために、本発明の他の態様である記録媒体は、入力映像から、人間が注意を向けやすい度合いである注目度を示す注目度映像を抽出する注目度映像抽出ステップと、入力映像を構成する各フレームである前記入力画像の各位置が顕著領域である確率を示す顕著領域事前確率画像を抽出する顕著領域事前確率画像抽出ステップと、前記入力画像の顕著領域および顕著領域外の領域にそれぞれ含まれる画像特徴量の尤度を示す特徴量尤度を算出する特徴量尤度算出ステップと、前記入力画像、前記顕著領域事前確率画像および前記特徴量尤度から、前記入力画像の顕著領域を示す顕著領域画像を抽出する顕著領域画像抽出ステップと、前記各入力画像に対し、前記注目度映像抽出ステップ、前記顕著領域事前確率画像抽出ステップ、前記特徴量尤度算出ステップおよび前記顕著領域画像抽出ステップを実行して得られる前記各顕著領域画像から前記顕著領域映像を生成する顕著領域映像生成ステップとをコンピュータに実行させるためのプログラムを記録したコンピュータ読み取り可能な記録憶媒体であって、前記顕著領域事前確率画像抽出ステップは、一の前記入力画像の各位置が顕著領域である確率を示す顕著領域事前確率画像を、前記注目度映像抽出ステップによって抽出された前記注目度映像内の前記入力画像に対応する画像である注目度画像および前記顕著領域画像に基づいて抽出し、前記特徴量尤度算出ステップは、前記特徴量尤度を、前記入力画像、前記注目度画像、前記顕著領域事前確率画像、前記顕著領域画像および前回迄に算出した前記特徴量尤度の少なくとも1つに基づいて算出することを特徴とする。   In order to solve the above problem, a recording medium according to another aspect of the present invention includes an attention level video extraction step for extracting an attention level video indicating a degree of attention that is a degree to which a human is easily directed from an input video; A saliency area prior probability image extracting step for extracting a saliency area prior probability image indicating a probability that each position of the input image which is each frame constituting the input video is a saliency area; From the feature amount likelihood calculating step for calculating the feature amount likelihood indicating the likelihood of the image feature amount included in each of the regions, the input image, the saliency region prior probability image, and the feature amount likelihood, the input image A saliency area image extracting step for extracting a saliency area image indicating a saliency area of the image, a focus level image extraction step, and a saliency area prior probability image extraction step for each input image. A program for causing a computer to execute a saliency area image generation step for generating the saliency area image from each saliency area image obtained by executing the feature amount likelihood calculation step and the saliency area image extraction step. In the recorded computer-readable recording medium, the remarkable area prior probability image extraction step includes: a remarkable area prior probability image indicating a probability that each position of the one input image is a remarkable area; Extracting based on the attention level image and the saliency area image, which are images corresponding to the input image in the attention level video extracted in the extraction step, and the feature amount likelihood calculating step includes calculating the feature amount likelihood. , The input image, the attention degree image, the saliency area prior probability image, the saliency area image, and the feature amount calculated up to the previous time And calculating, based on at least one degree.

本発明によれば、物体領域・背景領域に関する事前情報が全く与えられない場合においても領域分割が可能になる。従って、物体領域・背景領域に関する事前知識がない場合でも、精度良く物体領域と背景領域を分割して、注目している領域(物体領域)を抽出することができるようになる。   According to the present invention, it is possible to divide a region even when no prior information regarding the object region / background region is given. Therefore, even if there is no prior knowledge about the object region / background region, the object region and the background region can be divided with high accuracy and the region of interest (object region) can be extracted.

本発明の第1の実施形態に係る顕著領域映像生成装置1000による顕著領域事前確率画像の算出過程の模式図である。It is a schematic diagram of the calculation process of the saliency area prior probability image by the saliency area image generation device 1000 according to the first embodiment of the present invention. 顕著領域映像生成装置1000の機能ブロック図の一例である。3 is an example of a functional block diagram of a saliency area image generation device 1000. FIG. 注目度映像抽出部1の機能ブロック図である。It is a functional block diagram of the attention level video extraction unit 1. 顕著領域事前確率画像抽出部2、特徴量尤度算出部3および顕著領域画像抽出部4の機能ブロック図である。4 is a functional block diagram of a saliency area prior probability image extraction unit 2, a feature amount likelihood calculation unit 3, and a saliency area image extraction unit 4. FIG. 特徴量尤度の算出過程の模式図である。It is a schematic diagram of the feature amount likelihood calculation process. 顕著領域抽出グラフの例である。It is an example of a remarkable area | region extraction graph. 顕著領域映像生成装置1100の機能ブロック図の一例である。3 is an example of a functional block diagram of a saliency area image generation device 1100. FIG. 顕著領域抽出の例である。It is an example of remarkable area | region extraction. 顕著領域抽出の比較であるComparison of salient region extraction

(第1の実施形態)
以下、本発明の第1の実施形態に係る顕著領域映像生成装置1000について図面を参照して説明する。なお、本実施形態(後述する第2の実施形態も同様)においては、領域分割を映像顕著性に基づいて実現するため、以下、「顕著領域抽出」と「領域分割」とを同義で用いる。同様に、「顕著領域」と「物体領域」とを同義で用い、「非顕著領域」と「背景領域」とを同義で用いる。また、以下の説明において、式中の文字上部に ̄が付いた文字は、文中において文字の前に ̄を記載して示す。例えば、式中の文字(下記式1)は文中において ̄xと記載して示す。
(First embodiment)
Hereinafter, a saliency image generating apparatus 1000 according to a first embodiment of the present invention will be described with reference to the drawings. In the present embodiment (the same applies to the second embodiment described later), in order to realize area division based on video saliency, hereinafter, “significant area extraction” and “area division” are used synonymously. Similarly, “saliency area” and “object area” are used synonymously, and “non-salience area” and “background area” are used synonymously. In the following description, a character with a  ̄ in the upper part of the character in the formula is indicated with a  ̄ in front of the character in the sentence. For example, a character in the formula (the following formula 1) is indicated as  ̄x in the sentence.

また、式中の文字上部に〜が付いた文字は、文中において文字の前に〜を記載して示す。例えば、式中の文字(下記式2)は文中において〜ηと記載して示し、式中の文字(下記式3)は文中において〜Σと記載して示す。   In addition, a character with “˜” attached to the upper part of the character in the formula is indicated by “˜” before the character in the sentence. For example, a character in the formula (Formula 2 below) is indicated by ˜η in the sentence, and a character in the formula (Formula 3 below) is indicated by ˜Σ in the sentence.

なお、式中の文字(下記式4)と文中のgは同じである。   In addition, the character in a formula (following formula 4) and g in a sentence are the same.

顕著領域映像生成装置1000は、図1に示すように、外部から入力映像を取得し、当該入力映像を構成する各入力フレーム(各入力画像)からそれぞれの顕著領域を抽出した各顕著領域フレーム(各顕著領域画像)から構成される顕著領域映像を生成し、外部に出力する。   As shown in FIG. 1, the saliency area video generation apparatus 1000 obtains an input video from the outside, and extracts each saliency area frame (extracting each saliency area from each input frame (each input image) constituting the input video ( A saliency area image composed of each saliency area image) is generated and output to the outside.

顕著領域映像生成装置1000は、図2に示すように、注目度映像抽出部1、顕著領域事前確率画像抽出部2、特徴量尤度算出部3、顕著領域画像抽出部4および顕著領域映像生成部5を備える。注目度映像抽出部1は、図3に示すように、基礎注目度画像抽出部11、確率的基礎注目度画像抽出部12、確率的基礎注目度母数逐次推定部13、視線位置確率密度画像抽出部14、視線位置確率密度映像出力部15を備える。視線位置確率密度画像抽出部14は、視線移動状態変数更新部141、代表視線位置更新部142、代表視線位置重み係数算出部143、視線位置確率密度画像出力部144、代表視線位置集合再構成部145を備える。顕著領域事前確率画像抽出部2は、図4(a)に示すように、顕著領域事前確率画像生成部21および顕著領域事前確率画像更新部22を備える。特徴量尤度算出部3は、顕著領域特徴量尤度算出部31、図4(b)に示すように、非顕著領域特徴量尤度算出部32および特徴量尤度出力部33を備える。顕著領域特徴量尤度算出部31は、顕著領域特徴量尤度生成部311および顕著領域特徴量尤度更新部312を備える。非顕著領域特徴量尤度算出部32は、非顕著領域特徴量尤度生成部321および非顕著領域特徴量尤度更新部322を備える。顕著領域画像抽出部4は、図4(c)に示すように、顕著領域抽出グラフ生成部41および顕著領域抽出グラフ分割部42を備える。   As shown in FIG. 2, the saliency area image generation apparatus 1000 includes an attention degree image extraction unit 1, a saliency area prior probability image extraction unit 2, a feature amount likelihood calculation unit 3, a saliency area image extraction unit 4, and a saliency area image generation. Part 5 is provided. As shown in FIG. 3, the attention level video extraction unit 1 includes a basic attention level image extraction unit 11, a stochastic basic attention level image extraction unit 12, a stochastic basic attention level parameter sequential estimation unit 13, and a gaze position probability density image. An extraction unit 14 and a line-of-sight position probability density video output unit 15 are provided. The gaze position probability density image extraction unit 14 includes a gaze movement state variable update unit 141, a representative gaze position update unit 142, a representative gaze position weight coefficient calculation unit 143, a gaze position probability density image output unit 144, and a representative gaze position set reconstruction unit. 145. As shown in FIG. 4A, the saliency area prior probability image extraction unit 2 includes a saliency area prior probability image generation unit 21 and a saliency area prior probability image update unit 22. The feature amount likelihood calculation unit 3 includes a saliency region feature amount likelihood calculation unit 31, a non-salience region feature amount likelihood calculation unit 32, and a feature amount likelihood output unit 33, as shown in FIG. The saliency area feature quantity likelihood calculation unit 31 includes a saliency area feature quantity likelihood generation unit 311 and a saliency area feature quantity likelihood update unit 312. The non-significant region feature quantity likelihood calculating unit 32 includes a non-significant region feature amount likelihood generating unit 321 and a non-significant region feature amount likelihood updating unit 322. As shown in FIG. 4C, the saliency area image extraction unit 4 includes a saliency area extraction graph generation unit 41 and a saliency area extraction graph division unit 42.

注目度映像抽出部1は、入力映像を取得する。注目度映像抽出部1は、入力映像の各フレームの中で人間が注意を向けやすい度合いである注目度を示す映像である注目度映像を抽出する。注目度映像抽出部1は、抽出した注目度映像を顕著領域事前確率画像抽出部2に出力(供給)する。   The attention level video extraction unit 1 acquires an input video. The attention level video extracting unit 1 extracts a attention level video that is a video showing a level of attention, which is the degree to which a person is likely to pay attention in each frame of the input video. The attention level video extraction unit 1 outputs (supplies) the extracted attention level video to the saliency prior probability image extraction unit 2.

具体的には、注目度映像抽出部1は、視線位置推定の対象となる入力映像、視線位置推定に必要となる第1の母数である確率的基礎注目度母数Θ(t)、及び視線位置推定に必要となる第2の母数である視線位置確率密度母数Θ(t)が入力され、入力映像に含まれる時系列の各入力画像(各フレーム)内の各位置において、人間が視線を向ける確率を示した視線位置確率密度画像X(t)を算出する。更に、視線位置推定装置100は、算出した視線位置確率密度画像X(t)の時系列の映像である視線位置確率密度映像を出力する。 Specifically, the attention level video extraction unit 1 includes an input video that is a target of eye gaze position estimation, a stochastic basic attention degree parameter Θ s (t) that is a first parameter necessary for gaze position estimation, And a gaze position probability density parameter Θ x (t), which is a second parameter necessary for gaze position estimation, is input at each position in each time-series input image (each frame) included in the input video. Then, a line-of-sight position probability density image X (t) showing the probability that the person turns the line of sight is calculated. Furthermore, the gaze position estimation apparatus 100 outputs a gaze position probability density image that is a time-series image of the calculated gaze position probability density image X (t).

基礎注目度画像抽出部11は、入力された入力映像から視線位置推定を行う入力画像(フレーム)を取り出す。また、基礎注目度画像抽出部11は、取り出した入力画像の中で顕著な特性を持つ空間領域を示した画像である基礎注目度画像を抽出する。そして、抽出した基礎注目度画像を確率的基礎注目度画像抽出部12及び確率的基礎注目度母数逐次推定部13に出力する。   The basic attention level image extraction unit 11 extracts an input image (frame) for estimating the line-of-sight position from the input video input. Further, the basic attention level image extraction unit 11 extracts a basic attention level image that is an image showing a spatial region having a remarkable characteristic in the extracted input image. Then, the extracted basic attention level image is output to the probabilistic basic attention level image extraction unit 12 and the stochastic basic attention level parameter sequential estimation unit 13.

基礎注目度画像抽出部11における基礎注目度画像の抽出処理は、特許文献1に記載されている基礎注目度画像抽出部11の処理と同様であるため、処理内容の詳細な説明は省略する。ただし、本実施形態では、時刻iの入力画像から算出する基礎注目度画像を下記式(5)(以下、「基礎注目度画像 ̄S(i)」と表す)とする。   Since the basic attention level image extraction processing in the basic attention level image extraction unit 11 is the same as the processing of the basic attention level image extraction unit 11 described in Patent Document 1, detailed description of the processing content is omitted. However, in the present embodiment, the basic attention degree image calculated from the input image at time i is represented by the following equation (5) (hereinafter referred to as “basic attention degree image  ̄S (i)”).

確率的基礎注目度画像抽出部12は、現在の入力画像の各位置における顕著性を確率的な表現を用いて表示した画像である確率的基礎注目度画像S(t)を抽出する。なお、確率的基礎注目度画像抽出部12による確率的基礎注目度画像S(t)の抽出は、基礎注目度画像抽出部11から入力された基礎注目度画像 ̄S(i)、該確率的基礎注目度画像抽出部12が算出したこれまでの確率的基礎注目度画像S(t)、及び確率的基礎注目度母数Θ(t)に基づいて行われる。
また、確率的基礎注目度画像抽出部12によって抽出される確率的基礎注目度画像S(t)は、各位置yにおける確率的基礎注目度s(t,y)の期待値である下記式(6)(以下、「期待値^s(t,y|t)」と表す)及び標準偏差σs(t,y|t)を保持する画像である。
The probabilistic basic attention level image extraction unit 12 extracts a probabilistic basic attention level image S (t) that is an image in which saliency at each position of the current input image is displayed using a probabilistic expression. Note that the stochastic basic attention level image S (t) is extracted by the probabilistic basic attention level image extraction unit 12 based on the basic attention level image  ̄S (i) input from the basic attention level image extraction unit 11 and the probabilistic basic attention level image extraction unit 12. This is performed based on the previous stochastic basic attention level image S (t) calculated by the basic attention level image extraction unit 12 and the stochastic basic attention level parameter Θ s (t).
The probabilistic basic attention image S (t) extracted by the probabilistic basic attention image extracting unit 12 is an expected value of the probabilistic basic attention s (t, y) at each position y (the following formula ( 6) An image holding a standard deviation σs (t, y | t) (hereinafter, expressed as “expected value ^ s (t, y | t)”).

また、確率的基礎注目度画像抽出部12は、抽出した確率的基礎注目度画像S(t)を視線位置確率密度画像抽出部14及び確率的基礎注目度母数逐次推定部13に出力する。
なお、確率的基礎注目度画像抽出部12には、確率的基礎注目度母数逐次推定部13によって更新された確率的基礎注目度母数Θ(t+1)が入力される。
The stochastic basic attention level image extraction unit 12 outputs the extracted probabilistic basic attention level image S (t) to the gaze position probability density image extraction unit 14 and the stochastic basic attention level parameter sequential estimation unit 13.
The stochastic basic attention degree image extraction unit 12 receives the stochastic basic attention degree parameter Θ s (t + 1) updated by the stochastic basic attention degree parameter sequential estimation unit 13.

確率的基礎注目度画像抽出部12における確率的基礎注目度画像S(t)の抽出は、非特許文献1、2に記載されている手法によって算出することができる。また、確率的基礎注目度画像抽出部12における確率的基礎注目度画像S(t)の抽出方法は、特に限定されるものではないが、一例として、カルマンフィルタを用いた推定方法について述べる。   The extraction of the stochastic basic attention level image S (t) in the stochastic basic attention level image extraction unit 12 can be calculated by the methods described in Non-Patent Documents 1 and 2. The method of extracting the stochastic basic attention level image S (t) in the stochastic basic attention level image extracting unit 12 is not particularly limited, but an estimation method using a Kalman filter will be described as an example.

まず、現在(時刻t)の確率的基礎注目度画像S(t)(確率変数)の位置yにおける画素値s(t,y)(確率変数)が、現在の基礎注目度画像である下記式(7)(以下、「基礎注目度画像 ̄S(t)」と表す)の位置yにおける画素値である下記式(8)(以下、「画素値 ̄s(t,y)」と表す)、及び1時点前(時刻t−1)の確率的基礎注目度画像S(t−1)の位置yにおける画素値s(t−1,y)について、下記式(9)、(10)のような関係式を満たしているものとする。   First, the pixel value s (t, y) (probability variable) at the position y of the current (time t) stochastic basic attention image S (t) (probability variable) is the current basic attention image. (7) (hereinafter referred to as “pixel value  ̄ s (t, y)”) (7) (hereinafter referred to as “pixel value  ̄ s (t, y)”) , And the pixel value s (t−1, y) at the position y of the stochastic basic attention level image S (t−1) before one time point (time t−1), the following formulas (9) and (10) It is assumed that the relational expression is satisfied.

ここで、確率的基礎注目度母数Θ(t)は、時刻t及び位置yに依存する形で下記式(11)のように与えられているものとする。 Here, it is assumed that the stochastic basic attention degree parameter Θ s (t) is given by the following equation (11) in a form depending on the time t and the position y.

また、上記式(9)、(10)において、p(a|b)は、bが与えられたときのaの確率密度を示す。また、下記式(12)は期待値が下記式(13)で標準偏差がσである正規分布に従うsの確率密度を示し、下記式(14)のように表される。   In the above formulas (9) and (10), p (a | b) indicates the probability density of a when b is given. Also, the following formula (12) shows the probability density of s according to the normal distribution whose expected value is the following formula (13) and the standard deviation is σ, and is expressed as the following formula (14).

以降の説明においては、画素値 ̄s(t,y)を位置yにおける基礎注目度と呼ぶこととする。また、同様に、確率的基礎注目度s(t,y)を位置yにおける確率的基礎注目度と呼ぶこととする。また、特に必要な場合を除いて、位置yを省略するものとする。例えば、s(t,y)をs(t)と表す。   In the following description, the pixel value  ̄s (t, y) is referred to as the basic attention level at the position y. Similarly, the stochastic basic attention level s (t, y) is referred to as the stochastic basic attention level at the position y. Further, the position y is omitted unless particularly necessary. For example, s (t, y) is represented as s (t).

続いて、1時点前の確率的基礎注目度s(t−1)が、これまでの確率的基礎注目度画像抽出部12の処理により、下記式(15)のような確率密度を用いた表現にて抽出されているものとする。   Subsequently, the stochastic basic attention level s (t−1) one point before is expressed by using the probability density as in the following formula (15) by the processing of the probabilistic basic attention level image extraction unit 12 so far. It is assumed that it has been extracted at

上記式(15)において、下記式(16)は、時刻tから時刻tまでの基礎注目度の系列、下記式(17)は時刻1から時刻tまでの基礎注目度である下記式(18)が与えられているときの時刻tの確率的基礎注目度s(t)の期待値、σ(t|t)は、この時の標準偏差を示す。 In the above formula (15), the following formula (16) is a series of basic attention degrees from time t 1 to time t 2, and the following formula (17) is the following formula that is basic attention degrees from time 1 to time t 2. The expected value σ s (t 1 | t 2 ) of the stochastic basic attention degree s (t 1 ) at time t 1 when (18) is given indicates the standard deviation at this time.

このとき、確率的基礎注目度画像抽出部12は、現在の確率的基礎注目度s(t)の下記式(19)に示す確率密度における期待値である下記式(20)(以下、「期待値^s(t|t)」と表す)及び標準偏差σ(t|t)を、下記式(21)、(22)のようにして更新する。 At this time, the stochastic basic attention level image extraction unit 12 uses the following formula (20) (hereinafter, “expectation”) which is an expected value in the probability density shown in the following formula (19) of the current stochastic basic attention level s (t). Value ^ s (t | t) ") and standard deviation σ s (t | t) are updated as in the following equations (21) and (22).

なお、上述の確率的基礎注目度画像抽出部12における期待値^s(t|t)及び標準偏差σ(t|t)の更新は、画像中の各位置で独立して実行することができる。 Note that the update of the expected value ^ s (t | t) and the standard deviation σ s (t | t) in the above-described stochastic basic attention level image extraction unit 12 may be executed independently at each position in the image. it can.

また、上記式(10)に替えて、下記式(23)のような関係式を用いることによって、入力画像の各位置における動き成分を考慮した実施形態とすることも可能である。   Further, by using a relational expression such as the following expression (23) instead of the above expression (10), an embodiment in which a motion component at each position of the input image is taken into consideration can be obtained.

上記式(23)において、Δy(t)は時刻t、位置yにおけるオプティカルフローであり、例えば、特許文献1に記載されている運動特徴画像抽出部115と同様の方法を用いる。   In the above equation (23), Δy (t) is an optical flow at time t and position y. For example, the same method as the motion feature image extraction unit 115 described in Patent Document 1 is used.

上記に述べたカルマンフィルタを用いた推定方法では、各位置の確率的基礎注目度s(t,y)が空間方向で独立に抽出されていたが、確率的基礎注目度の空間的な連続性を導入することもできる。以下、動的マルコフ確率場と呼ばれる統計モデルに基づく確率的基礎注目度の記述を行い、平均場近似と呼ばれる統計解析手法により確率的基礎注目度s(t,y)を解析的に導出する場合について述べる。   In the estimation method using the Kalman filter described above, the stochastic basic attention degree s (t, y) at each position is extracted independently in the spatial direction. It can also be introduced. Hereinafter, a description of a stochastic basic attention degree based on a statistical model called a dynamic Markov random field will be described, and a stochastic basic attention degree s (t, y) will be derived analytically by a statistical analysis technique called mean field approximation. Is described.

まず、現在(時刻t)の確率的基礎注目度画像S(t)(確率変数)の位置yにおける画素値s(t,y)(確率変数)が、現在の基礎注目度画像 ̄S(t)の位置yにおける画素値 ̄s(t,y)、1時点前(時刻t−1)の確率的基礎注目度画像S(t−1)の位置yにおける画素値s(t−1,y)、及び位置yの近傍D(y)に含まれる各位置である下記式(24)における現在の確率的基礎注目度画像S(t)の画素値である下記式(25)について、下記式(26)〜(30)のような関係式を満たしているものとする。   First, the pixel value s (t, y) (probability variable) at the position y of the current (time t) probabilistic basic attention image S (t) (probability variable) is the current basic attention image  ̄S (t ) At the position y of the probabilistic basic attention level image S (t−1) one time before (time t−1). ), And the following formula (25) that is the pixel value of the current probabilistic basic attention level image S (t) in the following formula (24) that is each position included in the vicinity D (y) of the position y: It is assumed that the relational expressions (26) to (30) are satisfied.

ここで、確率的基礎注目度母数Θ(t)は、時刻t及び位置yに依存する形で下記式(31)のように再定義されているものとする。 Here, it is assumed that the stochastic basic attention degree parameter Θ s (t) is redefined as shown in the following formula (31) in a form depending on the time t and the position y.

近傍D(y)の決定方法としては、例えば、位置yの上下左右の4点、もしくはさらに斜め位置4点を加えた8点、などが考えられる。   As a method for determining the neighborhood D (y), for example, four points on the top, bottom, left, and right of the position y, or eight points including four additional oblique positions may be considered.

続いて、上述したカルマンフィルタを用いた推定方法と同様に、1時点前の確率的基礎注目度s(t−1,y)が、これまでの確率的基礎注目度画像抽出部12の処理により、下記式(32)のような確率密度を用いた表現にて抽出されているものとする。   Subsequently, as in the estimation method using the Kalman filter described above, the stochastic basic attention level s (t−1, y) one point before is obtained by the processing of the probabilistic basic attention level image extraction unit 12 so far. It is assumed that it is extracted by an expression using a probability density such as the following formula (32).

上記式(32)において、下記式(33)は、時刻tまでの基礎注目度画像である下記式(34)が与えられているときの時刻t・位置yの確率的基礎注目度s(t,y)の期待値、σ(t,y|t)は、この時の標準偏差を示す。 In the above equation (32), the following equation (33) is a stochastic basic attention degree s at time t 1 and position y when the following expression (34) which is a basic attention degree image up to time t 2 is given. (t 1, y) of the expected value, σ s (t 1, y | t 2) represents the standard deviation at this time.

このとき、確率的基礎注目度画像抽出部12は、位置yにおける現在の確率的基礎注目度s(t,y)の下記式(35)に示す確率密度における期待値^s(t,y|t)及び標準偏差σ(t,y|t)を更新することが目的となる。 At this time, the stochastic basic attention level image extraction unit 12 expects the current value of the probabilistic basic attention level s (t, y) at the position y in the probability density represented by the following equation (35) ^ s (t, y | The objective is to update t) and the standard deviation σ s (t, y | t).

また、確率的基礎注目度画像抽出部12による更新は、繰り返し計算を用いた下記式(36)〜(41)の方法によって行われる。   The update by the probabilistic basic attention level image extraction unit 12 is performed by the following formulas (36) to (41) using repetitive calculation.

上記式(36)〜(41)において、|D(y)|は集合D(y)の要素数を示す。また、上記式(36)〜(41)を用いた計算において、上記式(39)に示すような無限ステップの繰り返しは不可能であるため、実際には、第l+1ステップの出力である下記式(42)と第lステップの出力である下記式(43)との差が十分小さくなった時点で繰り返し計算を打ち切ることとする。   In the above formulas (36) to (41), | D (y) | indicates the number of elements of the set D (y). In addition, in the calculation using the above formulas (36) to (41), it is impossible to repeat the infinite steps as shown in the above formula (39). When the difference between (42) and the following equation (43) which is the output of the l-th step becomes sufficiently small, the calculation is repeatedly iterated.

また、ステップに関するインデックスlを固定したとき、上記式(38)に示す更新は、画像の各位置において独立に算出することができる。また、それ以外の更新式については、時刻tを固定することによって、同様に画像の各位置において独立に算出することができる。
このことによって、上述の確率的基礎注目度画像抽出部12における期待値^s(t,y|t)及び標準偏差σ(t,y|t)の更新は、上述したカルマンフィルタを用いた推定方法と同様に、画像中の各位置で独立して実行することができ、これらの更新処理を容易に並列化できる。
Further, when the index 1 relating to the step is fixed, the update shown in the above equation (38) can be calculated independently at each position of the image. Further, other update formulas can be calculated independently at each position of the image by fixing the time t.
As a result, the update of the expected value ^ s (t, y | t) and the standard deviation σ s (t, y | t) in the above-described stochastic basic attention level image extraction unit 12 is estimated using the above-described Kalman filter. Similar to the method, it can be executed independently at each position in the image, and these update processes can be easily parallelized.

確率的基礎注目度母数逐次推定部13は、基礎注目度画像抽出部11から入力された基礎注目度画像 ̄S(i)、確率的基礎注目度画像抽出部12から入力された確率的基礎注目度画像S(t)、及び事前に与えられた母数である確率的基礎注目度母数Θ(t)に基づいて、確率的基礎注目度母数Θ(t)を逐次的に更新する。
また、確率的基礎注目度母数逐次推定部13は、更新した確率的基礎注目度母数Θ(t+1)を確率的基礎注目度画像抽出部12に出力する。
The stochastic basic attention degree parameter sequential estimation unit 13 includes a basic attention degree image  ̄S (i) input from the basic attention degree image extraction unit 11 and a stochastic basis input from the stochastic basic attention degree image extraction unit 12. Based on the attention degree image S (t) and the probabilistic basic attention degree parameter Θ s (t) which is a parameter given in advance, the probabilistic basic attention degree parameter Θ s (t) is sequentially obtained. Update.
The stochastic basic attention degree parameter sequential estimation unit 13 outputs the updated stochastic basic attention degree parameter Θ s (t + 1) to the stochastic basic attention degree image extraction unit 12.

なお、確率的基礎注目度母数逐次推定部13は、確率的基礎注目度母数Θ(t)を更新していない場合は、事前に与えられた母数である確率的基礎注目度母数Θ(t)を確率的基礎注目度母数Θ(t+1)として確率的基礎注目度画像抽出部12に出力する。即ち、確率的基礎注目度画像抽出部12から確率的基礎注目度画像S(t)が入力されていない初期段階では、確率的基礎注目度母数Θ(t)の更新ができないため、入力された確率的基礎注目度母数Θ(t)をそのまま確率的基礎注目度画像抽出部12に出力する。 Note that the stochastic basic attention degree parameter sequential estimation unit 13 does not update the stochastic basic attention degree parameter Θ s (t), and the probabilistic basic attention degree mother, which is a parameter given in advance. The number Θ s (t) is output to the stochastic basic attention level image extraction unit 12 as a stochastic basic attention level parameter Θ s (t + 1). That is, since the stochastic basic attention degree parameter Θ s (t) cannot be updated in the initial stage where the stochastic basic attention degree image S (t) is not input from the stochastic basic attention degree image extraction unit 12, probabilistic basis prominence population parameter theta s a (t) is output as a probability basis saliency image extracting section 12.

確率的基礎注目度母数逐次推定部13における確率的基礎注目度母数Θ(t+1)の推定方法は特に限定されるものではないが、本実施形態においては、適応カルマンフィルタを用いた推定方法について述べる。 The estimation method of the stochastic basic attention degree parameter Θ s (t + 1) in the stochastic basic attention degree parameter sequential estimation unit 13 is not particularly limited, but in this embodiment, the estimation method using an adaptive Kalman filter Is described.

確率的基礎注目度母数逐次推定部13において、次の時刻t+1で用いられる確率的基礎注目度母数Θ(t+1)を、下記式(44)に示す。 In probabilistic basis attention parametric sequential estimation unit 13, the next time t + 1 with probability basis prominence population parameter theta s used the (t + 1), represented by the following formula (44).

確率的基礎注目度母数逐次推定部13は、既に基礎注目度画像抽出部11によって算出されている基礎注目度画像 ̄S(i)、及び確率的基礎注目度画像抽出部12によって算出されている確率的基礎注目度画像S(t)を構成する確率的基礎注目度の期待値及び標準偏差を用いて、下記式(45)〜(52)のように計算する。   The stochastic basic attention level parameter sequential estimation unit 13 is calculated by the basic attention level image  ̄S (i) that has already been calculated by the basic attention level image extraction unit 11 and the probabilistic basic attention level image extraction unit 12. Using the expected value and standard deviation of the probabilistic basic attention degree that constitutes the probabilistic basic attention degree image S (t), calculation is performed as in the following formulas (45) to (52).

上記式(45)〜(52)において、下記式(53)(以下、「 ̄σs1」と表す)及び下記式(54)(以下、「 ̄σs2」と表す)は基本確率的基礎注目度母数であり、予め定めておく、もしくは事前に学習によって算出しておくものである。 In the above formulas (45) to (52), the following formula (53) (hereinafter referred to as “ ̄σ s1 ”) and the following formula (54) (hereinafter referred to as “ ̄σ s2 ”) are basic stochastic basic attentions. It is a degree parameter and is determined in advance or calculated by learning in advance.

また、λs1及びλs2は予め定められた母数の混合比であり、これらの数値を適切に定めることで、逐次更新で獲得した母数である下記式(55),下記式(56)と、予め定めておいた母数 ̄σs1, ̄σs2とのバランスを制御することができる。 Also, λ s1 and λ s2 are predetermined mixing ratios of the parameters, and by appropriately determining these numerical values, the parameters obtained by successive updating are expressed by the following equations (55) and (56). And the predetermined parameters  ̄σ s1 and  ̄σ s2 can be controlled.

なお、λs1=λs2=0とすることにより、確率的基礎注目度母数逐次推定部13による確率的基礎注目度母数Θ(t+1)の推定処理を行わないことと等価になる。また、Nは過去の情報を保持するバッファの時間長である。 Note that setting λ s1 = λ s2 = 0 is equivalent to not performing the process of estimating the stochastic basic attention degree parameter Θ s (t + 1) by the stochastic basic attention degree parameter sequential estimation unit 13. N s is the time length of a buffer that holds past information.

視線位置確率密度画像抽出部14は、視線移動状態変数更新部141、代表視線位置更新部142、代表視線位置重み係数算出部143、視線位置確率密度画像出力部144、代表視線位置集合再構成部145から構成される。
視線位置確率密度画像抽出部14は、視線位置確率密度映像を構成するフレームである視線位置確率密度画像X(t)を抽出する。なお、視線位置確率密度画像抽出部14による視線位置確率密度画像X(t)の抽出は、確率的基礎注目度画像抽出部12から入力された確率的基礎注目度画像S(t)、該視線位置確率密度画像抽出部14が抽出したこれまでの視線位置確率密度画像X(t)、及び事前に与えられた母数である視線位置確率密度母数Θ(t)に基づいて行われる。
また、視線位置確率密度画像抽出部14は、視線位置確率密度画像X(t)を視線位置確率密度映像出力部5に出力する。
The gaze position probability density image extraction unit 14 includes a gaze movement state variable update unit 141, a representative gaze position update unit 142, a representative gaze position weight coefficient calculation unit 143, a gaze position probability density image output unit 144, and a representative gaze position set reconstruction unit. 145.
The line-of-sight position probability density image extraction unit 14 extracts a line-of-sight position probability density image X (t) that is a frame constituting the line-of-sight position probability density image. Note that the eye-gaze position probability density image X (t) is extracted by the eye-gaze position probability density image extraction unit 14 from the probabilistic basic attention level image S (t) input from the probabilistic basic attention level image extraction unit 12 and the line of sight. This is performed based on the line-of-sight position probability density image X (t) extracted by the position probability density image extraction unit 14 and the line-of-sight position probability density parameter Θ x (t) given in advance.
Further, the line-of-sight position probability density image extraction unit 14 outputs the line-of-sight position probability density image X (t) to the line-of-sight position probability density video output unit 5.

視線移動状態変数更新部141は、代表視線位置集合再構成部145から出力されたこれまでの出力である視線位置確率密度画像X(t)、及び事前に与えられた母数である視線位置確率密度母数Θ(t)に基づいて、これまでの視線位置確率密度画像X(t)に含まれている視線移動の大きさを制御する確率変数である視線移動状態変数u(t)を更新する。
また、視線移動状態変数更新部141は、更新した視線移動状態変数u(t)の集合である視線移動状態変数集合U(t)を代表視線位置更新部142及び代表視線位置集合再構成部145に出力する。
The line-of-sight movement state variable updating unit 141 includes a line-of-sight position probability density image X (t) that has been output from the representative line-of-sight position set reconstruction unit 145 and a line-of-sight position probability that is a parameter given in advance. Based on the density parameter Θ x (t), the gaze movement state variable u (t), which is a random variable for controlling the magnitude of the gaze movement included in the gaze position probability density image X (t) so far, is obtained. Update.
Also, the line-of-sight movement state variable update unit 141 converts the line-of-sight movement state variable set U (t), which is a set of the updated line-of-sight movement state variables u (t), to the representative line-of-sight position update unit 142 and the representative line-of-sight position set reconstruction unit 145. Output to.

視線移動状態変数更新部141における視線移動状態確率変数集合U(t)の更新方法は特に限定されるものではないが、本実施形態による方法について述べる。   The method of updating the line-of-sight movement state random variable set U (t) in the line-of-sight movement state variable updating unit 141 is not particularly limited, but the method according to this embodiment will be described.

まず、代表視線位置集合再構成部145の出力の一部として、1時点前(時刻t−1)の視線移動状態変数集合U(t−1)が、下記式(57)のように与えられているものとする。   First, as part of the output of the representative line-of-sight position set reconstruction unit 145, the line-of-sight movement state variable set U (t-1) one point before (time t-1) is given as in the following equation (57). It shall be.

上記式(57)において、Nは視線移動状態変数集合の要素数、即ち、視線移動状態変数のサンプル数を示す。また、各視線移動状態変数はm通りの値(1,2,・・・,m)のいずれかを取るものとする。 In the above equation (57), Nu represents the number of elements of the line-of-sight movement state variable set, that is, the number of samples of the line-of-sight movement state variable. Each eye movement state variable shall take one of the values of the street m u (1,2, ···, m u).

このとき、各視線移動状態変数のサンプルu(t−1)から、視線位置確率密度母数Θ(t)の1つである視線移動遷移確率行列Φ={φ(i,j)(i,j)に基づいて、現在の視線移動状態変数のサンプルu(t)をランダムに生成する。視線移動遷移確率行列はm行m列の行列によって表現され、i行j列の要素φ(i,j)によって、状態jから状態iへ遷移する確率を表現している。よって、Φは下記式(58)に示すような性質を満たす。 At this time, a line-of-sight movement transition probability matrix Φ = {φ (i, j) } , which is one of the line-of-sight position probability density parameters Θ x (t), from the samples u n (t−1) of the line-of-sight movement state variables. Based on (i, j) , a sample u n (t) of the current line-of-sight movement state variable is randomly generated. The line-of-sight movement transition probability matrix is expressed by a matrix of m u rows and m u columns, and expresses a probability of transition from the state j to the state i by an element φ (i, j) of i rows and j columns. Therefore, Φ satisfies the properties shown in the following formula (58).

即ち、上述のようにして生成した視線移動状態変数のサンプルu(t)の集合である下記式(59)が、現在の視線移動状態変数集合U(t)となる。
なお、現在の視線移動状態変数集合U(t)に含まれる視線移動状態変数のサンプルu(t)の経験確率分布が、視線移動状態変数の生起確率の近似となっている。
In other words, the following equation (59), which is a set of samples u n (t) of the line-of-sight movement state variables generated as described above, is the current line-of-sight movement state variable set U (t).
It should be noted that the empirical probability distribution of the line-of-sight movement state variable sample u n (t) included in the current line-of-sight movement state variable set U (t) is an approximation of the occurrence probability of the line-of-sight movement state variable.

別の実施形態として、視線移動状態変数更新部141では何の処理もしないことも可能である。しかし、これは、前述の視線移動状態変数更新部141における視線移動状態確率変数集合U(t)の更新処理において、m=1、即ち、視線移動状態が1通りしかないことと等価である。 As another embodiment, the line-of-sight movement state variable updating unit 141 can perform no processing. However, this is equivalent to m u = 1, that is, there is only one line-of-sight movement state in the update process of the line-of-sight movement state random variable set U (t) in the line-of-sight movement state variable update unit 141 described above. .

代表視線位置更新部142は、代表視線位置集合再構成部145から出力されたこれまでの出力である視線位置確率密度画像X(t)、視線移動状態変数更新部141から入力された視線移動状態変数集合U(t)、及び事前に与えられた母数である視線位置確率密度母数Θ(t)に基づいて、視線移動状態変数u(t)によって制御される視線移動を考慮し、これまでの視線位置確率密度画像X(t)に含まれている代表的な視線位置を表現する代表視線位置の集合である代表視線位置集合V(t)を更新する。
また、代表視線位置更新部142は、更新した代表視線位置集合V(t)を代表視線位置重み係数算出部143、視線位置確率密度画像出力部144、及び代表視線位置集合再構成部145に出力する。
The representative line-of-sight position update unit 142 includes the line-of-sight position probability density image X (t) output from the representative line-of-sight position set reconstruction unit 145 and the line-of-sight movement state input from the line-of-sight movement state variable update unit 141. Based on the variable set U (t) and the gaze position probability density parameter Θ x (t), which is a parameter given in advance, the gaze movement controlled by the gaze movement state variable u (t) is considered, A representative line-of-sight position set V (t), which is a set of representative line-of-sight positions representing representative line-of-sight positions included in the line-of-sight position probability density image X (t), is updated.
Further, the representative visual line position update unit 142 outputs the updated representative visual line position set V (t) to the representative visual line position weight coefficient calculation unit 143, the visual line position probability density image output unit 144, and the representative visual line position set reconstruction unit 145. To do.

代表視線位置更新部142における代表視線位置集合V(t)の更新方法は特に限定されるものではないが、本実施形態による方法について述べる。   The method of updating the representative line-of-sight position set V (t) in the representative line-of-sight position update unit 142 is not particularly limited, but the method according to the present embodiment will be described.

まず、視線移動状態変数u(t)によって視線位置x(t)が制御されるモデルとして、以下の2通りのモデルについて述べる。   First, the following two models will be described as models in which the line-of-sight position x (t) is controlled by the line-of-sight movement state variable u (t).

(モデル1):1時点前(時刻t−1)における視線位置が、視線位置x(t−1)として与えられているときの、現在(時刻t)における視線位置x(t)の生起確率を、現在の視線移動状態変数u(t)に依存する形で、下記式(60)のように与える。 (Model 1): Occurrence probability of the line-of-sight position x (t) at the current time (time t) when the line-of-sight position before the time point (time t-1) is given as the line-of-sight position x (t-1) Is given by the following equation (60) in a form depending on the current line-of-sight movement state variable u (t).

上記式(60)において、γxi及びσxi(i=0,1,・・・,m−1)はそれぞれ視線位置確率密度母数Θ(t)を構成する定数、下記式(61)(以下、「確率密度Q(x; ̄x,γ,σ)」と表す)は中心が下記式(62)、最頻距離がγ、最頻距離からの標準偏差に相当する母数をσとする下記式(63)に示す確率密度関数を表す。 In the above equation (60), γ xi and σ xi (i = 0, 1,..., M u −1) are constants constituting the line-of-sight position probability density parameter Θ x (t), respectively, ) (Hereinafter referred to as “probability density Q (x;  ̄x, γ, σ)”), the center is the following formula (62), the most frequent distance is γ, and the parameter corresponding to the standard deviation from the most frequent distance is The probability density function shown in the following formula (63) as σ is represented.

上記式(63)において、‖x‖はベクトルxのノルム、Zは確率密度Q(x; ̄x,γ,σ)の全定義域における積分値を1にするための下記式(64)で表される正規化定数を示す。 In the above equation (63), ‖x‖ is the norm of the vector x, Z L is the following equation (64) for setting the integral value in the entire domain of probability density Q (x;  ̄x, γ, σ) to 1. The normalization constant represented by

(モデル2):1時点前(時刻t−1)における視線位置が、視線位置x(t−1)として与えられているときの、現在(時刻t)における視線位置x(t)の生起確率として、ベータ分布を用いる。なお、1次元変量xに対するベータ分布は、定義域を下記式(65)とする下記式(66)によって定義される。 (Model 2): Occurrence probability of the gaze position x (t) at the present time (time t) when the gaze position before the time point (time t-1) is given as the gaze position x (t-1) A beta distribution is used. The beta distribution for the one-dimensional variable x is defined by the following formula (66) in which the domain of definition is the following formula (65).

上記式(66)において、a,bはそれぞれベータ分布を特徴付ける母数を示す。また、B(a,b)はベータ関数と呼ばれ、ベータ分布の全定義域における積分値を1にするための下記式(67)に示す正規化定数を表す。   In the above equation (66), a and b each indicate a parameter that characterizes the beta distribution. B (a, b) is called a beta function, and represents a normalization constant represented by the following formula (67) for setting the integral value in the entire domain of the beta distribution to 1.

なお、本実施形態においては、正規化定数の変量として位置xと、予め定められた原点xとの距離を用い、定義域を下記式(68)とするベータ分布である下記式(69)を用いる。   In the present embodiment, the following equation (69), which is a beta distribution using the distance between the position x and a predetermined origin x as a variable of the normalization constant, and having a defined region as the following equation (68), Use.

即ち、上述のベータ分布である下記式(70)は、下記式(71)によって与えられる。   That is, the following formula (70) which is the above-mentioned beta distribution is given by the following formula (71).

上記式(71)のように正規化されたベータ分布を用いることにより、1時点前(時刻t−1)の視線位置が視線位置x(t−1)として与えられているときの、現在(時刻t)の視線位置x(t)の生起確率は、現在の視線移動状態変数u(t)に依存する形で、下記式(72)のように与えられる。   By using the normalized beta distribution as in the above formula (71), the current position when the line-of-sight position one time before (time t−1) is given as the line-of-sight position x (t−1) ( The occurrence probability of the line-of-sight position x (t) at time t) is given by the following equation (72), depending on the current line-of-sight movement state variable u (t).

上記式(72)において、axi及びbxi(i=0,1,m−1)は、それぞれ視線位置確率密度母数Θ(t)を構成する定数を表す。 In the above formula (72), a xi and b xi (i = 0, 1, m u −1) represent constants constituting the line-of-sight position probability density parameter Θ x (t), respectively.

代表視線位置更新部142は、上述したモデルのいずれかを用いる方法によって、次に述べるように、代表視線位置集合V(t)を更新する。   The representative line-of-sight position update unit 142 updates the representative line-of-sight position set V (t) by a method using any of the above-described models as described below.

まず、代表視線位置集合再構成部145の出力の一部として、1時点前(時刻t−1)の代表視線位置集合V(t−1)が、下記式(73)のように与えられているものとする。   First, as a part of the output of the representative line-of-sight position set reconstruction unit 145, a representative line-of-sight position set V (t−1) one point before (time t−1) is given as in the following equation (73). It shall be.

上記式(73)において、Nは代表視線位置集合V(t)の要素数、即ち、代表視線位置のサンプル数を示す。なお、一般的な実施形態においては、代表視線位置集合V(t)の要素数Nを、視線移動状態変数集合U(t)の要素数Nと同一にしておく。 In the above formula (73), N x is the number of elements of the representative line of sight position set V (t), that is, the number of samples of the representative line-of-sight position. In a general embodiment, the number of elements N x of the representative line-of-sight position set V (t) is made the same as the number of elements N u of the line-of-sight movement state variable set U (t).

また、代表視線位置更新部142は、1時点前(時刻t−1)における代表視線位置の各サンプルx(t−1)から、上述したいずれかのモデルによって下記式(74)のように示す確率密度関数を用いて、現時点(時刻t)における代表視線位置のサンプルx(t)をランダムに生成する。 Also, the representative line-of-sight position update unit 142 calculates, from one sample x n (t−1) of the representative line-of-sight position one time before (time t−1), as shown in the following formula (74) using one of the models described above. A sample x n (t) of the representative line-of-sight position at the present time (time t) is randomly generated using the probability density function shown.

なお、上記式(74)のように示す確率密度関数を用いてランダムにサンプルを生成する方法は、サンプル生成に用いる確率密度関数が複雑である。このため、直接的な方法によってランダムなサンプルを生成することは困難である。しかし、上述のようなランダムのサンプル生成は、例えば、マルコフ連鎖モンテカルロ法に基づいたサンプル生成の方法を用いることができる。   Note that, in the method of randomly generating a sample using the probability density function represented by the above formula (74), the probability density function used for sample generation is complicated. For this reason, it is difficult to generate a random sample by a direct method. However, for the random sample generation as described above, for example, a sample generation method based on the Markov chain Monte Carlo method can be used.

次に、一般にMetropolis−Hastingsアルゴリズムと呼ばれるマルコフ連鎖モンテカルロ法に基づいたサンプル生成の詳細な方法について述べる。   Next, a detailed method of sample generation based on the Markov chain Monte Carlo method generally called the Metropolis-Hastings algorithm will be described.

まず、代表視線位置更新部142は、代表視線位置の仮のサンプルの初期値である下記式(75)として、1時点前(時刻t−1)における代表視線位置のサンプルを下記式(76)のように与える。   First, the representative line-of-sight position update unit 142 sets the sample of the representative line-of-sight position one time before (time t−1) as the following expression (75) that is the initial value of the temporary sample of the representative line-of-sight position: Give like.

次に、原点対称な確率密度関数を用いて、2次元ベクトルである下記式(77)を生成し、この2次元ベクトルである下記式(78)を、第k−1ステップの代表視線位置の仮のサンプルである下記式(79)に加えることによって、第kステップの代表視線位置の仮のサンプルである下記式(80)を下記式(81)のように生成する。   Next, the following equation (77), which is a two-dimensional vector, is generated using a probability density function that is symmetric with respect to the origin, and the following equation (78), which is a two-dimensional vector, is expressed as the representative line-of-sight position in the (k-1) th step. By adding to the following formula (79) that is a temporary sample, the following formula (80) that is a temporary sample of the representative line-of-sight position of the k-th step is generated as the following formula (81).

この原点対称な確率密度関数は、原点に対する対称性を満足しているのみで良く、例えば、原点を中心とする2次元正規分布、原点を中心とする各要素±δの範囲内の一様分布、などが考えられる。 This origin-symmetric probability density function only needs to satisfy the symmetry with respect to the origin, for example, a two-dimensional normal distribution centered on the origin, and uniform within the range of each element ± δ x centered on the origin. Distribution, etc. can be considered.

そして、第kステップの仮の代表視線位置のサンプルである上記式(80)の生起確率と第k−1ステップの代表視線位置の仮のサンプルである上記式(79)の生起確率との比である下記式(82)を下記式(83)に基づいて計算する。   Then, the ratio of the occurrence probability of the above formula (80), which is a sample of the temporary representative line-of-sight position in the k-th step, and the occurrence probability of the above-described expression (79), which is a temporary sample of the representative line-of-sight position in the k-1 step The following formula (82) is calculated based on the following formula (83).

最後に、下記式(84)の一様乱数である下記式(85)を発生させ、下記式(86)の場合のみ、第kステップの代表視線位置の仮のサンプルである上述の数式(80)を棄却して第k−1ステップの代表視線位置の仮のサンプルである上述の数式(79)に置き換える。   Finally, the following formula (85), which is a uniform random number of the following formula (84), is generated. Only in the case of the following formula (86), the above formula (80), which is a temporary sample of the representative line-of-sight position in the k-th step. ) Is replaced with the above-described mathematical expression (79), which is a temporary sample of the representative line-of-sight position in the (k−1) -th step.

その後、上述の仮のサンプルの生成ステップをあらかじめ定められた回数(K回)繰り返し、第Kステップの仮のサンプルである下記式(87)を下記式(88)に示すような時刻tの代表視線位置のサンプルとする。 Thereafter, the above-described tentative sample predetermined number (K x times) a generation step of repeating, the K x formula is a sample of tentative step (87) the following equation (88) shows such a time t The representative line-of-sight position is a sample.

上記に述べたとおり、マルコフ連鎖モンテカルロ法に基づいて、サンプルを生成する。この生成したサンプルの集合である下記式(89)が、即ち、現時点の代表視線位置集合V(t)となる。また、現時点の代表視線位置集合V(t)に含まれる代表視線位置サンプルの経験確率分布が、視線位置の生起確率の近似となっている。   As described above, a sample is generated based on the Markov chain Monte Carlo method. The following expression (89), which is a set of the generated samples, is the current representative gaze position set V (t). Also, the experience probability distribution of representative visual line position samples included in the current representative visual line position set V (t) is an approximation of the occurrence probability of the visual line position.

代表視線位置重み係数算出部143は、確率的基礎注目度画像抽出部12から入力された確率的基礎注目度画像S(t)、代表視線位置更新部142から入力された代表視線位置集合V(t)、及び事前に与えられた母数である視線位置確率密度母数Θ(t)に基づいて、各代表視線位置に関連付けられた重みである代表視線位置重み係数を算出する。
また、代表視線位置重み係数算出部143は、算出した代表視線位置重み係数の集合である代表視線位置重み係数集合である下記式(90)を視線位置確率密度画像出力部144及び代表視線位置集合再構成部145に出力する。
The representative gaze position weighting coefficient calculation unit 143 includes the probabilistic basic attention level image S (t) input from the probabilistic basic attention level image extraction unit 12 and the representative gaze position set V ( Based on t) and a gaze position probability density parameter Θ x (t) that is a parameter given in advance, a representative gaze position weight coefficient that is a weight associated with each representative gaze position is calculated.
In addition, the representative gaze position weight coefficient calculating unit 143 converts the following gaze position probability density image output unit 144 and the representative gaze position set as the representative gaze position weight coefficient set which is a set of the calculated representative gaze position weight coefficients. The data is output to the reconstruction unit 145.

代表視線位置重み係数算出部143における代表視線位置重み係数集合W(t)の抽出方法は特に限定されるものではないが、本実施形態においては、信号検出理論に基づいた代表視線位置重み係数集合W(t)の抽出方法について述べる。   The method of extracting the representative line-of-sight position weighting coefficient set W (t) in the representative line-of-sight position weighting coefficient calculation unit 143 is not particularly limited, but in this embodiment, the representative line-of-sight position weighting coefficient set based on the signal detection theory is used. A method for extracting W (t) will be described.

代表視線位置のサンプルx(t)(n=1,2,・・・,N)に関連付けられる代表視線位置重み係数w(t)は、下記式(91)及び下記式(92)によって算出される。なお、下記式(91)及び下記式(92)は、位置x(t)における確率的基礎注目度s(t,y)の実現値が、ある位置集合D(x(t))以外の位置yにおける確率的基礎注目度s(t,y)の実現値以上となる確率を算出するものである。 The representative line-of-sight position weighting coefficient w n (t) associated with the representative line-of-sight position sample x n (t) (n = 1, 2,..., N x ) is expressed by the following expressions (91) and (92). Is calculated by In addition, the following formula (91) and the following formula (92) indicate that the actual value of the probabilistic basic attention level s (t, y) at the position x n (t) is a certain position set D x (x n (t)). The probability of being equal to or higher than the actual value of the probabilistic basic attention degree s (t, y) at the position y other than is calculated.

なお、上記式(91)および上記式(92)においてのみs=s(t,x(t))の表記を用いている。また、上記式(91)および上記式(92)において、下記式(93)は、位置yにおける現在の確率的基礎注目度s(t,y)の確率分布関数を示し、位置yにおける現在の確率密度p(s(t,x))に対応して、下記式(94)のように定義される。 Note that the notation of s = s (t, x n (t)) is used only in the above formula (91) and the above formula (92). In the above formula (91) and the above formula (92), the following formula (93) represents the probability distribution function of the current probabilistic basic attention level s (t, y) at the position y, and the current distribution at the position y. Corresponding to the probability density p (s (t, x)), it is defined as the following formula (94).

ある位置集合D(x)の与え方は種々の方法が考えられるが、例えば、位置x以外の任意位置の集合、位置x以外で基礎注目度である下記式(95)が局所的に最大となる位置yの集合、位置x以外で確率的基礎注目度s(t,y)の期待値である下記式(96)が局所的に最大となる位置yの集合、などが考えられる。 Various methods can be considered for giving a certain position set D x (x). For example, a set of arbitrary positions other than the position x, and the following formula (95) which is a basic attention degree other than the position x is locally maximum. And a set of positions y where the following expression (96), which is an expected value of the probabilistic basic attention degree s (t, y), other than the position x is locally maximized.

上記に述べた信号検出理論に基づいた代表視線位置重み係数集合W(t)の抽出方法では、サンプリングにより代表視線位置集合V(t)及び代表視線位置重み係数集合W(t)を抽出していたが、サンプリングを用いずに代表視線位置重み係数集合W(t)を抽出することもできる。以下、サンプリングを用いずに代表視線位置重み係数集合W(t)を抽出する方法について述べる。   In the method of extracting the representative gaze position weight coefficient set W (t) based on the signal detection theory described above, the representative gaze position set V (t) and the representative gaze position weight coefficient set W (t) are extracted by sampling. However, the representative line-of-sight position weighting coefficient set W (t) can also be extracted without using sampling. Hereinafter, a method for extracting the representative gaze position weighting coefficient set W (t) without using sampling will be described.

サンプリングを用いずに代表視線位置重み係数集合W(t)を抽出する方法では、代表視線位置更新部142による代表視線位置集合V(t)の更新と、代表視線位置重み係数算出部143による代表視線位置重み係数集合W(t)の抽出とが同時に行われる。   In the method of extracting the representative gaze position weight coefficient set W (t) without using sampling, the representative gaze position update unit 142 updates the representative gaze position set V (t) and the representative gaze position weight coefficient calculation unit 143 represents the representative. The line-of-sight position weighting coefficient set W (t) is extracted at the same time.

まず、上述した信号検出理論に基づいた代表視線位置重み係数集合W(t)の抽出方法と同様に下記式(97)によって、位置x(t)における確率的基礎注目度s(t,y)の実現値が、ある位置集合D(x(t))以外の位置yにおける確率的基礎注目度s(t,y)の実現値以上になる確率を、入力画像中の各位置において算出する。 First, similarly to the method of extracting the representative gaze position weighting coefficient set W (t) based on the signal detection theory described above, the probabilistic basic attention s (t, y) at the position x (t) is obtained by the following equation (97). Is calculated at each position in the input image at a position in the input image that is equal to or greater than the actual value of the probabilistic basic attention level s (t, y) at a position y other than a position set D x (x (t)). .

続いて、この上記式(97)によって算出した確率分布を、下記式(98)〜(102)に示すようなEMアルゴリズムを用いて混合ガウス分布でモデル化する。即ち、混合ガウス分布の各パラメータであるガウス分布の混合比π(t)(n=1,2,・・・,M)、各ガウス分布の平均ベクトルである下記式(103)、及び共分散行列S(t)を、下記式(98)〜(102)によるモデル化のステップをk=1,2,・・・について各パラメータが収束するまで繰り返すことによって混合ガウス分布のモデルを導出する。なお、混合ガウス分布のモデル導出の際、位置xがどのガウス分布に所属するものであるかを表現する確率変数zを導入する。 Subsequently, the probability distribution calculated by the above equation (97) is modeled by a mixed Gaussian distribution using an EM algorithm as shown in the following equations (98) to (102). That is, the mixing ratio π n (t) (n = 1, 2,..., M y ) of the Gaussian distribution, which is each parameter of the mixed Gaussian distribution, the following equation (103) that is an average vector of each Gaussian distribution, and A model of the mixed Gaussian distribution is obtained by repeating the modeling step according to the following equations (98) to (102) for the covariance matrix S n (t) until each parameter converges for k = 1, 2,. To derive. When deriving a model of the mixed Gaussian distribution, a random variable z expressing which Gaussian distribution the position x belongs to is introduced.

上記式(98)〜(102)において、α(n=1,2,・・・,M)は、下記式(104)を満たすようにあらかじめ定められた定数を示す。 In the above formula (98) ~ (102), α n (n = 1,2, ···, M y) denotes a predetermined constant so as to satisfy the following equation (104).

なお、ガウス分布の混合比π(t)が予め定められた定数よりも小さいガウス分布については、ガウス分布の混合比π(t)に対する寄与が小さいものとして除去し、最終的に残ったN個のガウス分布によって混合ガウス分布を構成する。そして、この混合ガウス分布の各平均位置である上記式(103)(n=1,2,・・・,N)を、現時点(時刻t)の代表視線位置v(t)(n=1,2,・・・,N)として決定する。
このことから、サンプリングを用いずに代表視線位置重み係数集合W(t)を抽出する方法では、代表視線位置集合V(t)の要素数Nは、予め与えられるものではなく、入力画像によって異なることがわかる。
Note that the smaller Gaussian distribution than constant mixing ratio of the Gaussian distribution [pi n (t) is predetermined, and removed as contribution to the mixing ratio of the Gaussian distribution [pi n (t) is small, remaining finally A mixed Gaussian distribution is formed by N x Gaussian distributions. Then, the above equation (103) (n = 1, 2,..., N x ), which is each average position of the mixed Gaussian distribution, is converted into the representative line-of-sight position v n (t) (n = 1, 2,..., N x ).
From this, in the method of extracting the representative gaze position weight coefficient set W (t) without using sampling, the number of elements N x of the representative gaze position set V (t) is not given in advance, and depends on the input image. I can see that they are different.

一方、代表視線位置重み係数w(t)(n=1,2,・・・,N)については、1時点前(時刻t−1)の代表視線位置集合である下記式(105)、1時点前(時刻t−1)の代表視線位置重み係数集合である下記式(106)及び上記の混合ガウス分布の混合比π(t)に基づいて、下記式(107)のようにして算出する。 On the other hand, the representative line of sight position weight coefficient w n (t) (n = 1,2, ···, N x) for the 1 point before (time t-1) of the representative line of sight position set at a following formula (105) Based on the following formula (106) which is a representative gaze position weighting coefficient set one time before (time t−1) and the mixture ratio π n (t) of the above mixed Gaussian distribution, the following formula (107) is obtained. To calculate.

これは、即ち、1時点前の代表視線位置集合V(t)及び代表視線位置重み係数w(t)から構成される混合ガウス分布を、上記式(98)〜(102)によってモデル化した混合ガウス分布に、視線移動に関する確率密度p(s(t,x))を考慮して遷移させたものである。 In other words, a mixed Gaussian distribution composed of the representative line-of-sight position set V (t) and the representative line-of-sight position weighting coefficient w n (t) before one time point is modeled by the above equations (98) to (102). The transition is made to the mixed Gaussian distribution in consideration of the probability density p (s (t, x)) regarding the line-of-sight movement.

上記に述べたとおり、サンプリングを用いずに代表視線位置重み係数集合W(t)を抽出する方法では、代表視線位置更新部142及び代表視線位置重み係数算出部143は、代表視線位置集合である下記式(108)及び代表視線位置重み係数集合である上記式(90)を抽出し、視線位置確率密度画像出力部144に出力する。   As described above, in the method of extracting the representative gaze position weight coefficient set W (t) without using sampling, the representative gaze position update unit 142 and the representative gaze position weight coefficient calculation unit 143 are a representative gaze position set. The following formula (108) and the above formula (90) which is a representative gaze position weighting coefficient set are extracted and output to the gaze position probability density image output unit 144.

視線位置確率密度画像出力部144は、代表視線位置更新部142から入力された代表視線位置集合V(t)、及び代表視線位置重み係数算出部143から入力された代表視線位置重み係数集合W(t)に基づいて、代表視線位置確率密度画像H(t)を抽出する。
また、視線位置確率密度画像出力部144は、抽出した代表視線位置確率密度画像H(t)を代表視線位置集合再構成部145に出力する。
The line-of-sight position probability density image output unit 144 includes a representative line-of-sight position set V (t) input from the representative line-of-sight position update unit 142 and a representative line-of-sight position weight coefficient set W ( Based on t), a representative gaze position probability density image H (t) is extracted.
The line-of-sight position probability density image output unit 144 outputs the extracted representative line-of-sight position probability density image H (t) to the representative line-of-sight position set reconstruction unit 145.

視線位置確率密度画像出力部144による代表視線位置確率密度画像H(t)の算出方法は特に限定されるものではないが、本実施形態による方法について述べる。   The method for calculating the representative eye-gaze position probability density image H (t) by the eye-gaze position probability density image output unit 144 is not particularly limited, but the method according to the present embodiment will be described.

視線位置確率密度画像出力部144は、現時点(時刻t)の代表視線位置確率密度画像H(t)の位置x(t)における画素値を、代表視線位置集合V(t)及び代表視線位置重み集合W(t)に基づいて、下記式(109)のように算出する。   The gaze position probability density image output unit 144 uses the representative gaze position set V (t) and the representative gaze position weight as the pixel value at the position x (t) of the representative gaze position probability density image H (t) at the current time (time t). Based on the set W (t), calculation is performed as in the following formula (109).

上記式(109)において、f(・)は予め定められた関数であり、例えば、下記式(110)に示すデルタ関数や、下記式(111)に示す2次元正規分布などが考えられる。   In the above equation (109), f (•) is a predetermined function, and for example, a delta function represented by the following equation (110), a two-dimensional normal distribution represented by the following equation (111), and the like are conceivable.

代表視線位置集合再構成部145は、代表視線位置更新部142から入力された代表視線位置集合V(t)、視線移動状態変数更新部141から入力された視線移動状態変数集合U(t)、及び代表視線位置重み係数算出部143から入力された代表視線位置重み係数集合W(t)に基づいて、代表視線位置集合V(t)及び視線移動状態変数集合U(t)を、代表視線位置重み係数集合W(t)の示す重み配分に従って再構成する。
また、代表視線位置集合再構成部145は、代表視線位置重み係数集合W(t)を再構成する。
また、代表視線位置集合再構成部145は、再構成された代表視線位置集合V(t)、視線移動状態変数集合U(t)、及び代表視線位置重み係数集合W(t)に基づいた、視線位置確率密度画像X(t)を視線位置確率密度映像出力部5に出力する。
また、代表視線位置集合再構成部145は、視線位置確率密度画像X(t)を視線移動状態変数更新部141及び代表視線位置更新部142に出力する。
The representative line-of-sight position set reconstruction unit 145 includes a representative line-of-sight position set V (t) input from the representative line-of-sight position update unit 142, a line-of-sight movement state variable set U (t) input from the line-of-sight movement state variable update unit 141, And the representative gaze position weight coefficient set W (t) input from the representative gaze position weight coefficient calculation unit 143, the representative gaze position set V (t) and the gaze movement state variable set U (t) Reconfiguration is performed according to the weight distribution indicated by the weight coefficient set W (t).
Further, the representative gaze position set reconstruction unit 145 reconstructs the representative gaze position weight coefficient set W (t).
The representative gaze position set reconstruction unit 145, the reconstructed representative gaze position set V * (t), eye movement state variable set U * (t), and the representative line of sight position weight coefficient set W * (t) The line-of-sight position probability density image X (t) is output to the line-of-sight position probability density image output unit 5.
In addition, the representative gaze position set reconstruction unit 145 outputs the gaze position probability density image X (t) to the gaze movement state variable update unit 141 and the representative gaze position update unit 142.

代表視線位置集合再構成部145における代表視線位置集合V(t)及び視線移動状態変数集合U(t)の再構成方法は特に限定されるものではないが、本実施形態による方法について述べる。   The method for reconstructing the representative line-of-sight position set V (t) and the line-of-sight movement state variable set U (t) in the representative line-of-sight position set reconstruction unit 145 is not particularly limited, but the method according to this embodiment will be described.

まず、代表視線位置重み係数w(t)(n=1,2,・・・,N)の累積和c(t)を下記式(112)によって算出する。なお、累積和c(t)を算出する際に必要に応じて、代表視線位置重み係数w(t)の大きい順に代表視線位置v(t)、視線移動状態変数u(t)及び代表視線位置重み係数w(t)の並べ替えを行う。 First, a cumulative sum c n (t) of representative line-of-sight position weighting factors w n (t) (n = 1, 2,..., N x ) is calculated by the following equation (112). If necessary in calculating the cumulative sum c n (t), representative gaze position weight coefficient w n descending order representative line-of-sight position of the (t) v n (t) , eye movement state variable u (t) and The representative line-of-sight position weighting coefficient w n (t) is rearranged.

以降の処理のため、c(t)=0と定める。 For the subsequent processing, it is determined that c 0 (t) = 0.

次に、ある数κを下記式(113)の範囲でランダムに定め、以降、n=2,3,・・・,Nについて、κを下記式(114)のように定める。 Then, randomly determined a certain number kappa 1 in the range of the following formula (113), and later, n = 2,3, · · ·, for N x, defining a kappa n as the following equation (114).

そして、n=1,2,・・・,Nのそれぞれについて、下記式(115)の条件を満たす整数nを求める。 Then, for each of n = 1, 2,..., N x , an integer n * that satisfies the condition of the following formula (115) is obtained.

そして、新しい代表視線位置である下記式(116)を下記式(117)のように定める。   Then, the following expression (116), which is a new representative line-of-sight position, is defined as the following expression (117).

また、新しい視線移動状態変数である下記式(118)を下記式(119)のように定める。   Further, the following formula (118), which is a new line-of-sight movement state variable, is defined as the following formula (119).

なお、新しい代表視線位置重み係数である下記式(120)は、全て1/Nとする。 Incidentally, the new representative sight position weight coefficient formula (120) are all 1 / N x.

なお、上述した代表視線位置集合V(t)及び視線移動状態変数集合U(t)の再構成は、必ずしも全ての時刻において実施する必要はなく、例えば、ある一定時間間隔を置いて実施することや、全く実施しないこともできる。
また、例えば、下記式(121)に示す代表視線位置重み係数の偏りに関する条件を満たさないときのみ実施することもできる。
Note that the above-described reconstruction of the representative line-of-sight position set V (t) and the line-of-sight movement state variable set U (t) does not necessarily have to be performed at all times, for example, at certain time intervals. Or you can do nothing at all.
Further, for example, it can be performed only when the condition regarding the bias of the representative gaze position weighting coefficient represented by the following formula (121) is not satisfied.

上記式(121)において、Neffは下記式(122)を満たすように予め定められた定数である。 In the above formula (121), N eff is a constant determined in advance so as to satisfy the following formula (122).

また、代表視線位置集合再構成部145は、上記に述べた代表視線位置集合V(t)、視線移動状態変数集合U(t)、及び代表視線位置重み係数集合W(t)の再構成方法によって再構成した新しい代表視線位置集合である下記式(123)、新しい視線移動状態変数集合である下記式(124)、及び新しい代表視線位置重み係数集合である下記式(125)に基づいて、視線位置確率密度画像出力部144から入力された代表視線位置確率密度画像H(t)を再構成した視線位置確率密度画像X(t)を、視線位置確率密度画像抽出部14の出力として視線位置確率密度映像出力部5に出力する。   Further, the representative gaze position set reconstruction unit 145 reconstructs the representative gaze position set V (t), the gaze movement state variable set U (t), and the representative gaze position weight coefficient set W (t) described above. Based on the following formula (123) that is a new representative gaze position set reconstructed by the following formula, the following formula (124) that is a new gaze movement state variable set, and the following formula (125) that is a new representative gaze position weight coefficient set, The line-of-sight position probability density image X (t) reconstructed from the line-of-sight position probability density image H (t) input from the line-of-sight position probability density image output unit 144 is used as the output of the line-of-sight position probability density image extraction unit 14 Output to the probability density video output unit 5.

なお、代表視線位置集合再構成部145による代表視線位置集合V(t)、視線移動状態変数集合U(t)、及び代表視線位置重み係数集合W(t)の再構成を全く実施しない場合は、視線位置確率密度画像出力部144から入力された代表視線位置確率密度画像H(t)を視線位置確率密度画像抽出部14の出力である視線位置確率密度画像X(t)として視線位置確率密度映像出力部5に出力する。   In the case where the representative visual line position set reconstruction unit 145 does not reconstruct the representative visual line position set V (t), the visual line movement state variable set U (t), and the representative visual line position weight coefficient set W (t) at all. The visual line position probability density image H (t) input from the visual line position probability density image output unit 144 is used as the visual line position probability density image X (t) that is the output of the visual line position probability density image extraction unit 14. Output to the video output unit 5.

視線位置確率密度映像出力部5は、入力映像に含まれる時系列の各入力画像から、基礎注目度画像抽出部11、確率的基礎注目度画像抽出部12、確率的基礎注目度母数逐次推定部13、及び視線位置確率密度画像抽出部14の処理によって抽出される視線位置確率密度画像X(t)の時系列である視線位置確率密度映像を抽出して、出力する。   The line-of-sight position probability density video output unit 5 sequentially estimates a basic attention level image extraction unit 11, a probabilistic basic attention level image extraction unit 12, and a stochastic basic attention level parameter from each time-series input image included in the input video. The line-of-sight position probability density image X (t) extracted by the processing of the unit 13 and the line-of-sight position probability density image extracting unit 14 is extracted and output.

上記に述べたとおり、第1の実施形態によれば、視線位置推定の対象となる入力映像、確率的基礎注目度母数Θ(t)、及び視線位置確率密度母数Θ(t)に基づいて、視線位置確率密度映像を出力する場合に、確率的基礎注目度母数Θ(t)を逐次更新することができる。 As described above, according to the first embodiment, the input image, the probabilistic basic attention degree parameter Θ s (t), and the sight position probability density parameter Θ x (t), which are targets for eye gaze position estimation. On the basis of the above, the stochastic basic attention degree parameter Θ s (t) can be sequentially updated when the line-of-sight position probability density image is output.

また、第1の実施形態によれば、確率的基礎注目度画像抽出部12における期待値及び標準偏差の更新を、入力画像中の各位置で独立して実行することができる。その結果、確率的基礎注目度画像抽出部12による期待値及び標準偏差の更新処理を、複数コアを持つ計算機やGraphic Processor Unit(GPU)などの並列処理が可能な計算機上で容易に並列化することができ、処理を高速化することができる。   Further, according to the first embodiment, the expected value and the standard deviation in the probabilistic basic attention level image extraction unit 12 can be independently updated at each position in the input image. As a result, the expected value and standard deviation update processing by the probabilistic basic attention level image extraction unit 12 is easily parallelized on a computer capable of parallel processing, such as a computer having multiple cores or a graphic processor unit (GPU). And the processing speed can be increased.

なお、注目度映像抽出部1は、確率的基礎注目度画像抽出部12、確率的基礎注目度母数逐次推定部13、視線位置確率密度画像抽出部14、視線位置確率密度映像出力部15による上述の手法に代えて、下記非特許文献3、4に記載の手法でも、注目度映像を抽出することもできる。但し、確率的基礎注目度画像抽出部12、確率的基礎注目度母数逐次推定部13、視線位置確率密度画像抽出部14、視線位置確率密度映像出力部15による手法は、上述の如く並列処理に適しているため、非特許文献3、4に記載の手法に比べ、高速に計算(高速に注目度映像を抽出)することができる。   Note that the attention level video extraction unit 1 includes a probabilistic basic attention level image extraction unit 12, a stochastic basic attention level parameter sequential estimation unit 13, a gaze position probability density image extraction unit 14, and a gaze position probability density video output unit 15. Instead of the above-described method, the attention degree video can also be extracted by the methods described in Non-Patent Documents 3 and 4 below. However, the probabilistic basic attention level image extraction unit 12, the probabilistic basic attention level parameter sequential estimation unit 13, the gaze position probability density image extraction unit 14, and the gaze position probability density video output unit 15 employ parallel processing as described above. Therefore, compared to the methods described in Non-Patent Documents 3 and 4, it is possible to calculate at high speed (extract attention level video at high speed).

(非特許文献3)Derek Pang,Akisato Kimura,Tatsuto Takeuchi,Junji Yamato and Kunio Kashino,”A stochastic model of selective visual attention with a dynamic Bayesian network,”Proc.International Conference on Multimedia and Expo (ICME2008),pp.1073−1076, Hannover,Germany,June 2008.
(非特許文献4)Akisato Kimura,Derek Pang,Tatsuto Takeuchi, Junji Yamato and Kunio Kashino,”Dynamic Markov random fields for stochastic modeling of visual attention,”Proc.International Conference on Pattern Recognition (ICPR2008), Mo.BT8.35,Tampa, Florida,USA,December 2008.
(Non-Patent Document 3) Derek Pang, Akisato Kimura, Tatsuto Takeuchi, Junji Yamato and Kunio Kashino, “A stochastic model of selective visual attention with a dynamic Bayesian network,” Proc. International Conference on Multimedia and Expo (ICME2008), pp. 1073-1076, Hannover, Germany, June 2008.
(Non-Patent Document 4) Akisato Kimura, Derek Pang, Tatsuto Takeuchi, Junji Yamato and Kunio Kashino, “Dynamic Markov random fields for stochastic modeling of visual attention,” Proc. International Conference on Pattern Recognition (ICPR2008), Mo.BT8.35 , Tampa, Florida, USA, December 2008.

顕著領域事前確率画像抽出部2は、入力映像を構成する各フレームである入力画像の各位置が顕著領域である確率を示す顕著領域事前確率画像を抽出する。具体的には、顕著領域事前確率画像抽出部2は、注目度映像の一のフレームである注目度画像および顕著領域画像抽出部4によって抽出された顕著領域画像から、入力映像中の対応するフレームである入力画像の各位置が顕著領域である確率を表示する顕著領域事前確率画像を抽出する。換言すれば、顕著領域事前確率画像抽出部2は、一の入力画像の各位置が顕著領域である確率を示す顕著領域事前確率画像を、注目度映像抽出部1によって抽出された注目度映像内の当該入力画像に対応する画像である注目度画像および顕著領域画像抽出部4によって抽出された当該入力画像に対応する顕著領域画像に基づいて抽出する。顕著領域事前確率画像抽出部2は、抽出した顕著領域事前確率画像を特徴量尤度算出部3および顕著領域画像抽出部4に出力する。顕著領域事前確率画像抽出部2が顕著領域事前確率画像を抽出する方法は特に限定しないが、本実施形態においては、顕著領域事前確率画像生成部21と顕著領域事前確率画像更新部22とによって抽出する方法について説明する。   The saliency area prior probability image extraction unit 2 extracts a saliency area prior probability image indicating the probability that each position of the input image, which is each frame constituting the input video, is a saliency area. Specifically, the saliency area prior probability image extraction unit 2 corresponds to a corresponding frame in the input video from the attention level image which is one frame of the attention level video and the saliency area image extracted by the saliency area image extraction unit 4. A saliency area prior probability image that displays the probability that each position of the input image is a saliency area is extracted. In other words, the saliency area prior probability image extraction unit 2 extracts the saliency area prior probability image indicating the probability that each position of one input image is a saliency area in the attention level video extracted by the attention level video extraction unit 1. Are extracted based on the attention level image which is an image corresponding to the input image and the saliency area image corresponding to the input image extracted by the saliency area image extraction unit 4. The saliency area prior probability image extraction unit 2 outputs the extracted saliency area prior probability image to the feature amount likelihood calculation unit 3 and the saliency area image extraction unit 4. The method by which the saliency area prior probability image extraction unit 2 extracts the saliency area prior probability image is not particularly limited. In this embodiment, the saliency area prior probability image generation unit 21 and the saliency area prior probability image update unit 22 perform extraction. How to do will be described.

顕著領域事前確率画像生成部21は、注目度画像を入力し、注目度画像のみから顕著領域事前確率画像を生成する。顕著領域事前確率画像生成部21が注目度画像から顕著領域事前確率画像を生成する方法は特に限定しないが、本実施形態においては、ガウス混合分布モデルを利用した方法について説明する。   The saliency area prior probability image generation unit 21 receives the attention level image and generates a saliency area prior probability image only from the attention level image. The method for generating the saliency area prior probability image from the attention level image by the saliency area prior probability image generation unit 21 is not particularly limited. In the present embodiment, a method using a Gaussian mixture distribution model will be described.

顕著領域事前確率画像生成部21は、まず、時刻tの注目度画像(即ち、基礎注目度画像 ̄S(t)若しくは視線位置確率密度画像X(t))が、それぞれ中心位置〜x(t)・共分散行列〜Σs,j(t)(j=1,2,・・・,M)を持ち、混合比が〜ηs,j(t)であるM個のガウス分布の混合によって構成されていると仮定し、そのモデルパラメータ(即ち、M個の中心位置、共分散行列、混合比)を注目度画像(b)から推定する。推定方法の具体例は、以下の2つである。 The saliency area prior probability image generation unit 21 first calculates the attention level image at time t (that is, the basic attention level image  ̄S (t) or the gaze position probability density image X (t)) from the center position to x j ( t) M s Gaussian distributions with covariance matrix ~ Σ s, j (t) (j = 1,2, ..., M s ) and mixing ratio ~ η s, j (t) Model parameters (ie, M s center positions, covariance matrix, mixture ratio) are estimated from the attention degree image (b). Specific examples of the estimation method are the following two.

(推定方法1)
EMアルゴリズムを用いて導出する。このとき、EMアルゴリズムに与える各サンプルは注目度画像(b)のある特定の位置xに対応し、位置xにおける画素値と等しい値の重みを持つことに注意する。EMアルゴリズムによるガウス混合分布パラメータの推定は、下記式(126)〜下記式(129)をk=1,2,・・・で繰り返すことによって行われ、各パラメータが収束した時点で手順を打ち切り、パラメータを固定する。
(Estimation method 1)
Derived using EM algorithm. At this time, it should be noted that each sample given to the EM algorithm corresponds to a specific position x of the attention degree image (b) and has a weight equal to the pixel value at the position x. The estimation of the Gaussian mixture distribution parameter by the EM algorithm is performed by repeating the following formula (126) to the following formula (129) with k = 1, 2,..., And the procedure is terminated when each parameter converges. Fix the parameters.

ここで、g(x;〜x,Σ)は多次元正規分布であり、次元数がDとするときには下記式(130)で定義される。   Here, g (x; ˜x, Σ) is a multidimensional normal distribution, and is defined by the following formula (130) when the number of dimensions is D.

また、視線位置確率密度画像X(t)の位置xにおける画素値を、位置xをEMアルゴリズムのサンプルとみなしたときの重みとして利用していることから、ここではw(t)と表現している。 Further, since the pixel value at the position x of the line-of-sight position probability density image X (t) is used as a weight when the position x is regarded as a sample of the EM algorithm, it is expressed here as w x (t). ing.

(推定方法2)
注目度画像の画素値の極大値をM個検出し、極大値となる位置を中心位置〜x(t)(j=1,2,・・・,M)として定め、その位置の注目度画像の画素値を混合比〜ηs,j(t)とする。共分散行列〜Σs,j(t)については、第1の推定方法と同様にして求めるか、予め定めておいた値を利用する。
(Estimation method 2)
M s maximum values of the pixel values of the attention level image are detected, and a position where the maximum value is obtained is determined as a center position to x j (t) (j = 1, 2,..., M s ). Let the pixel value of the attention level image be a mixture ratio ~ η s, j (t). The covariance matrix ~ Σ s, j (t) is obtained in the same manner as in the first estimation method or a predetermined value is used.

以上のようにして、顕著領域事前確率画像生成部21は、注目度画像からガウス混合分布のモデルパラメータを推定し、顕著領域事前確率画像を生成する。具体的には、顕著領域事前確率画像生成部21は、モデルパラメータの1つである混合比を、その最大値が1と等しくなるように正規化し、その後に各位置におけるガウス混合分布の確率を計算して、当該位置の顕著領域事前確率画像 ̄Ξ(t)の画素値ξ(x,t)とする(下記式(131))。 As described above, the saliency area prior probability image generation unit 21 estimates the model parameter of the Gaussian mixture distribution from the attention degree image, and generates the saliency area prior probability image. Specifically, the saliency area prior probability image generation unit 21 normalizes the mixture ratio, which is one of the model parameters, so that its maximum value is equal to 1, and then calculates the probability of the Gaussian mixture distribution at each position. The pixel value ξ 1 (x, t) of the saliency area prior probability image 1 1 (t) at the position is calculated (the following formula (131)).

上記の実施形態では、全ての位置においてガウス混合分布を用いた方法によって顕著領域事前確率画像を生成しているが、顕著領域が画像の中心位置に存在しやすいことを考慮した後述する第2の実施形態も可能である。即ち、この場合には、顕著領域事前確率画像を第1の実施形態と同様に生成した後、画像の左右両端もしくは上下左右の端の一定領域のピクセル値を強制的に0とする。当該方法は、画像の端に顕著領域が存在する可能性を排除することを意味している。若しくは、顕著領域事前確率画像を第1の実施形態と同様に生成した後、画像の中心位置からの距離に比例する重みを顕著領域事前確率画像に掛け合わせ、その出力を新たに顕著領域事前確率とする実施形態も考えられる。上記の通り、顕著領域事前確率画像生成部21は、顕著領域事前確率画像 ̄Ξ(t)を生成し出力する。 In the above-described embodiment, the saliency area prior probability image is generated by the method using the Gaussian mixture distribution at all positions. However, the second area described later in consideration of the fact that the saliency area tends to exist at the center position of the image. Embodiments are possible. That is, in this case, after generating the saliency area prior probability image in the same manner as in the first embodiment, the pixel values of the constant areas at the left and right ends or the upper and lower left and right ends of the image are forcibly set to zero. This method means eliminating the possibility of a significant area at the edge of the image. Alternatively, after the saliency area prior probability image is generated in the same manner as in the first embodiment, the weight proportional to the distance from the center position of the image is multiplied by the saliency area prior probability image, and the output is newly generated as the saliency area prior probability. An embodiment is also conceivable. As described above, the saliency area prior probability image generation unit 21 generates and outputs the saliency area prior probability image  ̄Ξ 1 (t).

顕著領域事前確率画像更新部22は、顕著領域画像を用いて顕著領域事前確率画像生成部21によって生成された顕著領域事前確率画像を更新する。即ち、顕著領域事前確率画像更新部22は、顕著領域事前確率画像生成部21によって生成された顕著領域事前確率画像、および、顕著領域画像抽出部4によって抽出された顕著領域画像を入力し、顕著領域画像を用いて顕著領域事前確率画像を更新する。顕著領域事前確率画像更新部22が顕著領域事前確率画像を更新する方法は特に限定しないが、本実施形態においては、カルマンフィルタの原理を利用する方法について説明する。   The saliency area prior probability image update unit 22 updates the saliency area prior probability image generated by the saliency area prior probability image generation unit 21 using the saliency area image. That is, the saliency area prior probability image update unit 22 inputs the saliency area prior probability image generated by the saliency area prior probability image generation unit 21 and the saliency area image extracted by the saliency area image extraction unit 4, and The saliency prior probability image is updated using the region image. The method by which the saliency area prior probability image update unit 22 updates the saliency area prior probability image is not particularly limited, but in the present embodiment, a method using the principle of the Kalman filter will be described.

現時点(時刻t)の顕著領域事前確率画像Ξ(t)(確率変数)の位置xにおける画素値ξ(x,t)(確率変数)が、現時点の更新前顕著領域事前確率画像 ̄Ξ1(t)の位置xにおける画素値 ̄ξ1(x,t)、および、1時点前(時刻t−1)の顕著領域画像A(t−1)の位置xにおける画素値a(x,t−1)について、下記式(132)(133)を満たしているものとする。 The pixel value ξ 1 (x, t) (probability variable) at the position x of the saliency area prior probability image Ξ 1 (t) (probability variable) at the current time (time t) is the current saliency area prior probability image  ̄Ξ1 before update. The pixel value  ̄ξ1 (x, t) at the position x in (t) and the pixel value a (x, t− at the position x in the saliency area image A (t−1) one time before (time t−1). It is assumed that the following formulas (132) and (133) are satisfied for 1).

ここで、θ=(σ,σ)はあらかじめ与えられるパラメータである。またf(・) は、顕著領域画像の画素値を実数値に変換する関数であり、例えば、下記式(134)(135のように設定する。 Here, θ = (σ 1 , σ 2 ) is a parameter given in advance. Further, f (•) is a function for converting the pixel value of the saliency area image into a real value, and is set as in the following formulas (134) and (135), for example.

このとき、顕著領域事前確率画像更新部22は、現時点の顕著領域事前確率画像Ξ(t)の位置xにおける画素値ξ(x,t)を、カルマンフィルタの原理を利用することにより、下記式(136)(137)によって更新する。 At this time, the saliency area prior probability image update unit 22 uses the Kalman filter principle to calculate the pixel value ξ 1 (x, t) at the position x of the current saliency area prior probability image 1 1 (t) as follows. It is updated by the equations (136) and (137).

上記の実施形態では、各時刻の顕著領域事前確率を保持しているが、この分散を次の時刻での更新の際に利用しなくてもよい。即ち、上記式(136)(137)に、下記式(138)を追加してもよい。   In the above embodiment, the saliency prior probability at each time is held, but this distribution may not be used when updating at the next time. That is, the following formula (138) may be added to the above formulas (136) and (137).

上記の通り、顕著領域事前確率画像更新部22は、顕著領域事前確率画像 ̄Ξ(t) をΞ(t)に更新し、更新後の顕著領域事前確率画像Ξ1(t)を出力する。 As described above, salient region prior probability image update section 22 updates the marked area prior probability image ¯Ξ 1 (t) to Ξ 1 (t), and outputs the salient region prior probability image .xi.1 (t) the updated .

上記の通り、顕著領域事前確率画像抽出部2は、図5に示すように、顕著領域事前確率画像Ξ1(t)を抽出(生成、更新)し、出力する。   As described above, the saliency area prior probability image extraction unit 2 extracts (generates and updates) and outputs the saliency area prior probability image Ξ1 (t) as shown in FIG.

特徴量尤度算出部3は、入力画像の顕著領域および顕著領域外の領域にそれぞれ含まれる画像特徴量の尤度を示す特徴量尤度を算出する。具体的には、特徴量尤度算出部3は、特徴量尤度を、入力画像、注目度画像、顕著領域事前確率画像抽出部2によって抽出された顕著領域事前確率画像、顕著領域画像抽出部4によって抽出された顕著領域画像および前回迄に算出した特徴量尤度の少なくとも1つに基づいて算出する。例えば、特徴量尤度算出部3は、入力画像、顕著領域事前確率画像、顕著領域画像および前回迄に算出した特徴量尤度から、特徴量尤度を算出する。特徴量尤度算出部3は、算出した特徴量尤度を顕著領域画像抽出部4に出力する。特徴量尤度算出部3が徴量尤度を算出する方法は特に限定しないが、本実施形態においては、顕著領域特徴量尤度算出部31と非顕著領域特徴量尤度算出部32と特徴量尤度出力部33とによって算出する方法について説明する。   The feature amount likelihood calculating unit 3 calculates a feature amount likelihood indicating the likelihood of the image feature amount included in each of the saliency area and the area outside the saliency area of the input image. Specifically, the feature amount likelihood calculating unit 3 determines the feature amount likelihood as the input image, the attention degree image, the saliency area prior probability image extracted by the saliency area prior probability image extracting unit 2, and the saliency area image extracting unit. 4 is calculated based on at least one of the salient region image extracted by 4 and the feature amount likelihood calculated until the previous time. For example, the feature amount likelihood calculation unit 3 calculates the feature amount likelihood from the input image, the saliency area prior probability image, the saliency area image, and the feature amount likelihood calculated until the previous time. The feature amount likelihood calculating unit 3 outputs the calculated feature amount likelihood to the saliency area image extracting unit 4. The method by which the feature amount likelihood calculation unit 3 calculates the collection likelihood is not particularly limited, but in the present embodiment, the saliency region feature amount likelihood calculation unit 31, the non-salience region feature amount likelihood calculation unit 32, and the feature A method of calculation by the quantity likelihood output unit 33 will be described.

顕著領域特徴量尤度算出部31は、顕著領域に含まれる画像特徴量の尤度を示す顕著領域特徴量尤度を、入力画像、顕著領域事前確率画像、顕著領域画像および前回迄に算出した顕著領域特徴量尤度のうち少なくとも1つに基づいて算出する。顕著領域特徴量尤度算出部31が顕著領域特徴量尤度を算出する方法は特に限定しないが、本実施形態においては、顕著領域特徴量尤度生成部311と顕著領域特徴量尤度更新部312とによって算出する方法について説明する。   The saliency area feature amount likelihood calculation unit 31 calculates the saliency area feature amount likelihood indicating the likelihood of the image feature amount included in the saliency area by the input image, the saliency area prior probability image, the saliency area image, and the previous time. Calculation is made based on at least one of the saliency feature amount likelihoods. The method by which the saliency area feature amount likelihood calculating unit 31 calculates the saliency area feature amount likelihood is not particularly limited, but in the present embodiment, the saliency area feature amount likelihood generating unit 311 and the saliency area feature amount likelihood updating unit are included. 312 will be described.

顕著領域特徴量尤度生成部311は、入力画像、顕著領域事前確率画像および顕著領域画像に基づいて顕著領域特徴量尤度を新たに生成(算出)し、出力する。顕著領域特徴量尤度生成部311が顕著領域特徴量尤度を生成する方法は、特に限定しないが、本実施形態においては、ガウス混合分布モデルを利用した方法について説明する。   The saliency area feature amount likelihood generation unit 311 newly generates (calculates) and outputs a saliency area feature amount likelihood based on the input image, the saliency area prior probability image, and the saliency area image. The method for generating the saliency feature amount likelihood by the saliency feature amount likelihood generation unit 311 is not particularly limited, but in the present embodiment, a method using a Gaussian mixture distribution model will be described.

顕著領域特徴量尤度生成部311は、まず、時刻tにおいて、顕著領域に特有の特徴量の確率分布である顕著領域特徴量確率が、それぞれ平均〜c(t)・共分散行列〜Σf,j(t)(j= 1,2,・・・,M)を持ち、混合比が〜ηf,j(t)であるM個のガウス分布の混合によって構成されていると仮定し、これらのモデルパラメータを、顕著領域事前確率画像の画素値で重み付けした入力画像の画素値から推定する。モデルパラメータの推定には、例えばEMアルゴリズムを用いる。具体的には、下記式(139)〜下記式(142)をk=1,2,・・・ で繰り返すことによって行われ、各パラメータが収束した時点で手順を打ち切り、パラメータを固定する。 The saliency area feature amount likelihood generation unit 311 first calculates, at time t, the saliency area feature quantity probabilities, which are probability distributions of the characteristic quantities peculiar to the saliency area, from the mean to c j (t) and the covariance matrix to Σ, respectively. f, j (t) (j = 1, 2,..., M f ), and a mixture ratio of ˜η f, j (t) is composed of a mixture of M f Gaussian distributions. Assuming these model parameters are estimated from the pixel values of the input image weighted with the pixel values of the saliency prior probability image. For example, an EM algorithm is used for estimating the model parameters. Specifically, the following formula (139) to the following formula (142) are repeated by k = 1, 2,..., And when each parameter converges, the procedure is terminated and the parameters are fixed.

ここで、入力画像の位置xにおける画素値は、RGBの3次元ベクトルとしてc(x, t)で表現される。上記のようにして、顕著領域特徴量尤度生成部311は、推定したガウス混合分布のモデルパラメータから、顕著領域特徴量尤度を算出する。具体的には、推定したモデルパラメータで特徴付けられるガウス混合分布を尤度とする下記式(143)によって算出する。   Here, the pixel value at the position x of the input image is expressed as c (x, t) as a three-dimensional RGB vector. As described above, the saliency area feature amount likelihood generation unit 311 calculates the saliency area feature amount likelihood from the model parameters of the estimated Gaussian mixture distribution. Specifically, it is calculated by the following equation (143) with a Gaussian mixture distribution characterized by the estimated model parameters as the likelihood.

上記の通り、顕著領域特徴量尤度生成部311は、顕著領域特徴量尤度 ̄ψ(c, t)を生成(算出)し、出力する。 As described above, the saliency area feature amount likelihood generation unit 311 generates (calculates) and outputs the saliency area feature amount likelihood  ̄ψ 1 (c, t).

顕著領域特徴量尤度更新部312は、顕著領域特徴量尤度生成部311によって生成された顕著領域特徴量尤度を更新する。具体的には、顕著領域特徴量尤度更新部312は、入力画像、顕著領域画像および前回迄に更新した更新後の顕著領域特徴量尤度のうち少なくとも1つに基づいて、顕著領域特徴量尤度生成部311によって生成された顕著領域特徴量尤度を更新する。顕著領域特徴量尤度更新部312が顕著領域特徴量尤度(e1)を更新する方法は特に限定しないが、本実施形態においては、以下の2通りの方法を説明する。   The saliency area feature amount likelihood update unit 312 updates the saliency area feature amount likelihood generated by the saliency area feature amount likelihood generation unit 311. Specifically, the saliency area feature amount likelihood updating unit 312 is based on at least one of the input image, the saliency area image, and the updated saliency area feature amount likelihood updated until the previous time, and the saliency area feature amount. The saliency area feature amount likelihood generated by the likelihood generation unit 311 is updated. The method by which the saliency feature amount likelihood update unit 312 updates the saliency feature amount likelihood (e1) is not particularly limited, but in the present embodiment, the following two methods will be described.

(更新方法1)
2種類の顕著領域特徴量尤度を混合することによって更新する。具体的には、求めるべき現時点(時刻t)の顕著領域特徴量尤度ψ(c,t)を、顕著領域特徴量尤度生成部311から出力された更新前の顕著領域特徴量尤度ψ1(c,t)、および、1時点前(時刻t−1)の顕著領域特徴量尤度 ̄ψ(c,t−1)を、あらかじめ定められた混合比λにて混合する下記式(144)によって計算する。
(Update method 1)
It is updated by mixing two types of saliency area feature likelihoods. Specifically, the saliency area feature amount likelihood ψ 1 (c, t) at the current time (time t) to be obtained is determined as the saliency area feature amount likelihood before update output from the saliency area feature amount likelihood generation unit 311. ψ1 (c, t) and the saliency area feature likelihood  ̄ψ 1 (c, t-1) one point before (time t-1) are mixed at a predetermined mixing ratio λ c Calculated according to equation (144).

(更新方法2)
1時点前(時刻t−1)の顕著領域画像A(t−1)に基づいて、1時点前の顕著領域特徴量尤度ψ(c,t−1)を更新した上で第1の実施形態と同様の方法を実施する。具体的には、1時点前の顕著領域画像A(t−1)において顕著領域であるとされた領域(下記式(145)によって示される顕著領域Aobj(t))を取り出し、顕著領域Aobj(t)にある入力画像中の画素値から、顕著領域特徴量尤度生成部311に示した方法と同様の方法で顕著領域特徴量尤度ψ(c,t−1)を再学習する。但し、本実施形態では、重みとして、顕著領域事前確率画像に代えて顕著領域画像を用いるものとする。1時点前の顕著領域特徴量尤度ψ1(c,t−1)を再学習した後、第1の実施形態と同様の方法により、現在の顕著領域特徴量尤度ψ(c,t)を生成する。
(Update method 2)
Based on the saliency area image A (t−1) before one time point (time t−1), the saliency area feature amount likelihood ψ 1 (c, t−1) before one time point is updated and the first A method similar to that of the embodiment is performed. Specifically, a salient region and areas at one time before the salient region image A (t-1) to (formula (salient region A obj indicated by 145) (t)) is taken out, salient region A The saliency feature amount likelihood ψ 1 (c, t−1) is relearned from the pixel value in the input image at obj (t) by the same method as the method shown in the saliency feature amount likelihood generation unit 311. To do. However, in the present embodiment, the saliency area image is used as the weight instead of the saliency area prior probability image. After re-learning the saliency area feature amount likelihood ψ1 (c, t-1) one point before, the current saliency region feature amount likelihood ψ 1 (c, t) is obtained by the same method as in the first embodiment. Is generated.

上記の通り、顕著領域特徴量尤度更新部312は、顕著領域特徴量尤度 ̄ψ(c,t)をψ(c,t)に更新し、出力する。上記の通り、顕著領域特徴量尤度算出部31は、顕著領域特徴量尤度ψ(c,t) を算出し、出力する。 As described above, the saliency area feature quantity likelihood update unit 312 updates the saliency area feature quantity likelihood  ̄ψ 1 (c, t) to ψ 1 (c, t) and outputs it. As described above, the saliency area feature amount likelihood calculating unit 31 calculates and outputs the saliency area feature amount likelihood ψ 1 (c, t).

非顕著領域特徴量尤度算出部32は、顕著領域外の領域に含まれる画像特徴量の尤度を示す非顕著領域特徴量尤度を、入力画像、顕著領域事前確率画像、顕著領域画像および前回迄に算出した非顕著領域特徴量尤度のうち少なくとも1つに基づいて算出する。非顕著領域特徴量尤度算出部32が非顕著領域特徴量尤度を算出する方法は特に限定しないが、本実施形態においては、非顕著領域特徴量尤度生成部321と非顕著領域特徴量尤度更新部322とによって算出する方法について説明する。   The non-significant region feature amount likelihood calculating unit 32 converts the non-significant region feature amount likelihood indicating the likelihood of the image feature amount included in the region outside the saliency region into an input image, a saliency region prior probability image, a saliency region image, and The calculation is performed based on at least one of the non-significant region feature amount likelihoods calculated so far. The method by which the non-significant region feature amount likelihood calculating unit 32 calculates the non-significant region feature amount likelihood is not particularly limited, but in the present embodiment, the non-significant region feature amount likelihood generating unit 321 and the non-significant region feature amount likelihood are calculated. A method of calculation by the likelihood update unit 322 will be described.

非顕著領域特徴量尤度生成部321は、入力画像、顕著領域事前確率画像および顕著領域画像に基づいて非顕著領域特徴量尤度を新たに生成(算出)し、出力する。非顕著領域特徴量尤度生成部321が非顕著領域特徴量尤度を生成する方法は特に限定しないが、本実施形態においては、ガウス混合分布モデルを利用した方法について説明する。当該方法は、前述の顕著領域特徴量尤度生成部311の方法とほぼ同様であるが、顕著領域事前確率画像に代えて、顕著領域事前確率画像Ξ(t)の各画素値ξ(x,t)をある規則に従って変換して生成した画像である非顕著領域事前確率画像Ξ(t)を用いるものとする。当該変換規則として、例えば以下の2つの方法が考えられる。 The non-salience area feature amount likelihood generation unit 321 newly generates (calculates) and outputs a non-salience area feature amount likelihood based on the input image, the saliency area prior probability image, and the saliency area image. The method by which the non-salience feature amount likelihood generation unit 321 generates the non-salience feature amount likelihood is not particularly limited, but in the present embodiment, a method using a Gaussian mixture distribution model will be described. The method is almost the same as the method of salient region feature amount likelihood generator 311 described above, instead of the salient region prior probability image, each pixel value of the notable region prior probability image Ξ 1 (t) ξ 1 ( Assume that a non-salience area prior probability image Ξ 2 (t), which is an image generated by converting x, t) according to a certain rule, is used. For example, the following two methods can be considered as the conversion rule.

(方法1)
非顕著領域事前確率画像(f)の位置xにおける画素値ξ(x,t)を、1−ξ(x,t)に変換する。
(方法2)
ξ(x,t)=0である位置xのみ、その位置の非顕著領域事前確率画像(f)の画素値を1とする。それ以外の位置は、画素値を0とする。
(Method 1)
The pixel value ξ 2 (x, t) at the position x of the non-salience area prior probability image (f) is converted into 1−ξ 1 (x, t).
(Method 2)
Only in the position x where ξ 1 (x, t) = 0, the pixel value of the non-salience area prior probability image (f) at that position is set to 1. For other positions, the pixel value is 0.

上記の通り、非顕著領域特徴量尤度生成部321は、非顕著領域特徴量尤度 ̄ψ(c,t)を生成(算出)し、出力する。 As described above, the non-salient region feature quantity likelihood generating unit 321, a non-salient region feature quantity likelihood ¯ψ 2 (c, t) generates (calculates), outputs.

非顕著領域特徴量尤度更新部322は、非顕著領域特徴量尤度生成部321によって生成された非顕著領域特徴量尤度を更新する。具体的には、非顕著領域特徴量尤度更新部322は、入力画像、非顕著領域画像および前回迄に更新した更新後の非顕著領域特徴量尤度のうち少なくとも1つに基づいて、非顕著領域特徴量尤度生成部321によって生成された非顕著領域特徴量尤度を更新する。なお、非顕著領域画像は、顕著領域事前確率画像抽出部2によって抽出される顕著領域外の領域に係る画像である。非顕著領域特徴量尤度更新部322が非顕著領域特徴量尤度を更新する方法は、顕著領域特徴量尤度更新部312の方法と同様である。但し、顕著領域事前確率画像に代えて非顕著領域事前確率画像を、顕著領域特徴量尤度に代えて非顕著領域特徴量尤度を、顕著領域に代えて非顕著領域(下記式(146)によって示される顕著領域Abkg(t))を用いる。 The non-significant region feature quantity likelihood update unit 322 updates the non-significant region feature amount likelihood generated by the non-significant region feature amount likelihood generation unit 321. Specifically, the non-salience area feature amount likelihood update unit 322 performs non-remarkable based on at least one of the input image, the non-salience area image, and the updated non-salience area feature amount likelihood updated up to the previous time. The non-salience area feature amount likelihood generated by the saliency area feature amount likelihood generation unit 321 is updated. The non-significant area image is an image relating to an area outside the saliency area extracted by the saliency area prior probability image extraction unit 2. The method by which the non-salience area feature amount likelihood update unit 322 updates the non-salience area feature amount likelihood is the same as the method of the saliency area feature amount likelihood update unit 312. However, instead of the saliency area prior probability image, the non-salience area prior probability image is replaced with the saliency area feature amount likelihood, and the non-salience area feature amount likelihood is replaced with the saliency area. The saliency area A bkg (t)) indicated by is used.

上記の通り、非顕著領域特徴量尤度更新部322は、非顕著領域特徴量尤度ψ(c,t)を更新し、更新後の非顕著領域特徴量尤度を出力する。上記の通り、非顕著領域特徴量尤度算出部32は、非顕著領域特徴量尤度ψ(c,t)を抽出し、出力する。 As described above, the non-significant region feature value likelihood update unit 322 updates the non-significant region feature value likelihood ψ 2 (c, t), and outputs the updated non-salience region feature value likelihood. As described above, the non-salience area feature amount likelihood calculating unit 32 extracts and outputs the non-salience area feature amount likelihood ψ 2 (c, t).

特徴量尤度出力部33は、顕著領域特徴量尤度および非顕著領域特徴量尤度を加算して特徴量尤度として出力する。   The feature quantity likelihood output unit 33 adds the saliency area feature quantity likelihood and the non-salience area feature quantity likelihood and outputs the result as a feature quantity likelihood.

顕著領域画像抽出部4は、入力画像、顕著領域事前確率画像および特徴量尤度から、入力画像の顕著領域を示す顕著領域画像を抽出する。顕著領域画像抽出部4は、抽出した顕著領域画像を顕著領域事前確率画像抽出部2、特徴量尤度算出部3および顕著領域映像生成部5に出力する。顕著領域画像抽出部4が顕著領域画像を抽出する方法は特に限定しないが、本実施形態においては、非特許文献1に記載の方法を基礎としたグラフカットを用いる方法について説明する。当該方法は、顕著領域抽出グラフ生成部41と、顕著領域抽出グラフ分割部42とによって顕著領域画像を抽出する。   The saliency area image extraction unit 4 extracts a saliency area image indicating the saliency area of the input image from the input image, the saliency area prior probability image, and the feature amount likelihood. The saliency area image extraction unit 4 outputs the extracted saliency area image to the saliency area prior probability image extraction part 2, the feature amount likelihood calculation part 3, and the saliency area video generation part 5. The method of extracting the saliency area image by the saliency area image extracting unit 4 is not particularly limited, but in this embodiment, a method using a graph cut based on the method described in Non-Patent Document 1 will be described. In this method, the saliency area extraction graph generation unit 41 and the saliency area extraction graph division unit 42 extract the saliency area image.

顕著領域抽出グラフ生成部41は、入力画像、顕著領域事前確率画像および特徴量尤度を入力し、顕著領域画像を抽出するためのグラフである顕著領域抽出グラフを生成し、出力する。   The saliency area extraction graph generation unit 41 inputs the input image, the saliency area prior probability image, and the feature amount likelihood, and generates and outputs a saliency area extraction graph that is a graph for extracting the saliency area image.

具体的には、顕著領域抽出グラフ生成部41は、まず、時刻tの顕著領域抽出グラフG(t)の頂点として、入力画像の各位置x∈Iに対応する頂点と、顕著領域・非顕著領域のラベルにそれぞれ対応する頂点の2種類の頂点を用意する。即ち、頂点は総計画素数+2個となる。以降、簡単のため、各位置xに対応する頂点をv、顕著領域のラベルに対応する頂点をSOURCES、非顕著領域のラベルに対応する頂点をSINKTとして、それぞれ表現する。また、顕著領域抽出グラフの辺として、近傍位置に対応する頂点の間相互に配置される有向辺であるn−linkと、SOURCEから各頂点・各頂点からSINKにそれぞれ配置される有向辺であるt−linkの2種類の辺を用意する。近傍としては、例えば上下左右の4近傍、若しくはさらに斜め方向を加えた8近傍を考える。このようにして、顕著領域抽出グラフは有向グラフとして、例えば図6に示すような形で構成される。 Specifically, the saliency area extraction graph generation unit 41 first sets vertices corresponding to the respective positions x∈I of the input image as vertices of the saliency area extraction graph G (t) at time t, and saliency areas / non-salience. Two types of vertices corresponding to the label of the area are prepared. That is, the total number of vertices is +2. Hereinafter, for simplicity, the vertex corresponding to each position x is represented as v x , the vertex corresponding to the label of the saliency area is SOURCES, and the vertex corresponding to the label of the non-salience area is SINKT, respectively. Also, as the sides of the saliency area extraction graph, n-links, which are directed edges between vertices corresponding to neighboring positions, and directed edges, which are arranged from SOURCE to each vertex and from each vertex to SINK, respectively. Two types of t-link sides are prepared. As the neighborhood, for example, 4 neighborhoods in the upper, lower, left, and right directions, or 8 neighborhoods including an oblique direction are considered. In this way, the saliency area extraction graph is configured as a directed graph, for example, in the form shown in FIG.

次に、顕著領域抽出グラフ生成部41は、顕著領域抽出グラフの各辺にコストを与える。t−linkのコストは、顕著領域事前確率画像および特徴量尤度から算出される。具体的には、SOURCESから頂点vへのt−linkのコストC(S,v;t)は対応する非顕著領域事前確率と非顕著領域特徴量尤度の和で、頂点vからSINKTへのt−linkのコストC(T,v;t)は、対応する顕著領域事前確率と顕著領域特徴量尤度を用いて、下記式(147)(148)のように与えられる。 Next, the saliency area extraction graph generation unit 41 gives a cost to each side of the saliency area extraction graph. The cost of t-link is calculated from the significant area prior probability image and the feature amount likelihood. Specifically, the cost C (S, v x ; t) of t-link from SOURCES to vertex v x is the sum of the corresponding non-salience area prior probability and non-salience area feature likelihood, and from vertex v x The cost C (T, v x ; t) of t-link to SINKT is given by the following equations (147) and (148) using the corresponding saliency area prior probability and saliency area feature likelihood.

一方、n−linkのコストは、近接画素間の輝度値の類似性に基づいて算出される。具体的には、ある2点vとvとの間のn−linkのコストC(v,v)は、下記式(149)で与えられる。 On the other hand, the cost of n-link is calculated based on the similarity of luminance values between adjacent pixels. Specifically, the cost C (v x, v y) of n-link between one 2-point v x and v y is given by the following formula (149).

顕著領域抽出グラフ分割部42は、顕著領域抽出グラフを入力し、顕著領域抽出グラフを分割することで顕著領域画像を生成し、出力する。   The saliency area extraction graph dividing unit 42 receives the saliency area extraction graph, divides the saliency area extraction graph, and generates and outputs a saliency area image.

具体的には、顕著領域抽出グラフ分割部42は、まず、顕著領域抽出グラフに含まれる頂点を、SOURCEを含む部分集合とSINKを含む部分集合に分割することを考える。このとき、SOURCE側の頂点の部分集合からSINK側の頂点の部分集合へまたがる辺のコストが最も小さくなるように分割する。逆向き、即ち、SINK側の頂点の部分集合からSOURCE側の頂点の部分集合へまたがる辺のコストは考慮しないことに注意する。このような問題は、グラフの最小カット問題と呼ばれ、グラフの最大フロー問題と等価であることが知られている。この最大フロー問題を解く方法として、非特許文献1の他、下記非特許文献5に記載の「Ford-Fulkerson algorithm」、下記非特許文献6に記載の「Goldberg-Tarjan algorithm」などが一般に広く知られている。   Specifically, the saliency area extraction graph dividing unit 42 first considers dividing the vertices included in the saliency area extraction graph into a subset including SOURCE and a subset including SINK. At this time, the division is performed so that the cost of the side extending from the subset of the vertexes on the SOURCE side to the subset of the vertexes on the SINK side is minimized. Note that it does not take into account the cost in the opposite direction, ie, the edge spanning from a subset of vertices on the SINK side to a subset of vertices on the SOURCE side. Such a problem is called a graph minimum cut problem, and is known to be equivalent to a graph maximum flow problem. As a method for solving this maximum flow problem, in addition to Non-Patent Document 1, “Ford-Fulkerson algorithm” described in Non-Patent Document 5 below, “Goldberg-Tarjan algorithm” described in Non-Patent Document 6 below, etc. are generally widely known. It has been.

(非特許文献5)L.R.Ford,D.R.Fulkerson:“Maximal flow through a network,”Canadial Journal of Mathematics,Vol.8,pp.399-404,1956.
(非特許文献6)A.V.Goldberg,R.E.Tarjan:“A new approach to the maximum-flow problem,”Journal of the ACM,Vol.35,pp.921-940,1988.
(Non-Patent Document 5) LRFord, DRFulkerson: “Maximal flow through a network,” Canadial Journal of Mathematics, Vol. 8, pp. 399-404, 1956.
(Non-Patent Document 6) AVGoldberg, REtarjan: “A new approach to the maximum-flow problem,” Journal of the ACM, Vol. 35, pp. 921-940, 1988.

上記の方法で顕著領域抽出グラフを分割した結果、SOURCEを含む部分グラフに属した頂点に対応する画素位置を顕著領域Aobj(t)に、SINKを含む部分グラフに属した頂点に対応する画素位置を非顕著領域Abkg(t)に所属させる。顕著領域画像は、下記式(150)に示すように、顕著領域に属する位置の画素値を1、非顕著領域に属する位置の画素値を0とする画像である。 As a result of dividing the salient region extraction graph by the above method, the pixel position corresponding to the vertex belonging to the subgraph including SOURCE is set to the salient region A obj (t), and the pixel corresponding to the vertex belonging to the subgraph including SINK The position is made to belong to the non-salience area A bkg (t). As shown in the following formula (150), the saliency area image is an image in which the pixel value of the position belonging to the saliency area is 1 and the pixel value of the position belonging to the non-salience area is 0.

上記の通り、顕著領域抽出グラフ分割部42は、顕著領域画像A(x,t) を抽出し、出力する。つまり、顕著領域画像抽出部4は、顕著領域画像A(x,t)を抽出し、この顕著領域画像を出力する。   As described above, the saliency area extraction graph dividing unit 42 extracts and outputs the saliency area image A (x, t). That is, the saliency area image extraction unit 4 extracts the saliency area image A (x, t) and outputs the saliency area image.

顕著領域映像生成部5は、各入力画像に対し、注目度映像抽出部1、顕著領域事前確率画像抽出部2、特徴量尤度算出部3および顕著領域画像抽出部4を実行して得られる各顕著領域画像から顕著領域映像を生成する。換言すれば、顕著領域映像生成部5は、注目度映像抽出部1、顕著領域事前確率画像抽出部2、特徴量尤度算出部3および顕著領域画像抽出部4を順に各入力画像に対して繰り返し実行することにより抽出された顕著領域画像の集合を入力し、集合内の顕著領域画像から構成される顕著領域映像を生成する(集合内の顕著領域画像を時系列に並べて顕著領域映像を生成する)。顕著領域映像生成部5は、生成した顕著領域映像を外部に出力する。   The saliency area video generation unit 5 is obtained by executing the attention level video extraction unit 1, the saliency area prior probability image extraction unit 2, the feature amount likelihood calculation unit 3, and the saliency area image extraction unit 4 for each input image. A saliency area image is generated from each saliency area image. In other words, the saliency area video generation unit 5 applies the attention level video extraction unit 1, the saliency area prior probability image extraction unit 2, the feature amount likelihood calculation unit 3, and the saliency area image extraction unit 4 in order to each input image. A set of salient area images extracted by repeated execution is input and a salient area video composed of salient area images in the set is generated (a salient area video is generated by arranging the salient area images in the set in time series) To do). The saliency area image generation unit 5 outputs the generated saliency area image to the outside.

(第2の実施形態)
以下、本発明の第2の実施形態に係る顕著領域映像生成装置1100について図面を参照して説明する。顕著領域映像生成装置1100は、外部から入力映像を取得し、当該入力映像を構成する各入力フレーム(各入力画像)からそれぞれの顕著領域を抽出した各顕著領域フレーム(各顕著領域画像)から構成される顕著領域映像を生成し、外部に出力する。
(Second Embodiment)
Hereinafter, a saliency image generating apparatus 1100 according to a second embodiment of the present invention will be described with reference to the drawings. The saliency area image generation device 1100 is configured from each saliency area frame (each saliency area image) obtained by acquiring an input image from outside and extracting each saliency area from each input frame (each input image) constituting the input image. Generated saliency area video and output it to the outside.

顕著領域映像生成装置1100は、図7に示すように、注目度映像抽出部1、顕著領域事前確率画像抽出部2、特徴量尤度算出部3、顕著領域画像抽出部4、顕著領域映像生成部5、平滑化画像群生成部6、顕著領域画像確定部7を備える。注目度映像抽出部1、顕著領域画像抽出部4および顕著領域映像生成部5は、第1の実施形態と同様であるため説明を省略する。   As shown in FIG. 7, the saliency area image generation device 1100 includes an attention degree image extraction unit 1, a saliency area prior probability image extraction unit 2, a feature amount likelihood calculation unit 3, a saliency area image extraction unit 4, and a saliency area image generation unit. Unit 5, smoothed image group generation unit 6, and saliency area image determination unit 7. The attention level video extraction unit 1, the saliency area image extraction unit 4 and the saliency area image generation unit 5 are the same as those in the first embodiment, and thus the description thereof is omitted.

平滑化画像群生成部6は、入力画像を異なる解像度によってそれぞれ平滑化した複数の平滑化画像からなる平滑化画像群を生成する。つまり、平滑化画像群生成部6は、入力画像を入力し、当該入力画像を異なる解像度によってそれぞれ平滑化した平滑化画像群を生成する。平滑化画像群生成部6は、生成した平滑化画像群を特徴量尤度算出部3、顕著領域画像抽出部4および顕著領域画像確定部7に出力する。平滑化画像群生成部6が平滑化画像群を生成する方法は特に限定しないが、本実施形態では、入力画像に対し平滑化と縮小とを繰り返す方法について説明する。   The smoothed image group generation unit 6 generates a smoothed image group composed of a plurality of smoothed images obtained by smoothing the input image with different resolutions. That is, the smoothed image group generation unit 6 inputs an input image and generates a smoothed image group obtained by smoothing the input image with different resolutions. The smoothed image group generation unit 6 outputs the generated smoothed image group to the feature amount likelihood calculation unit 3, the saliency area image extraction unit 4, and the saliency area image determination unit 7. A method for generating the smoothed image group by the smoothed image group generation unit 6 is not particularly limited. In the present embodiment, a method for repeating smoothing and reduction on the input image will be described.

平滑化画像群生成部6は、時刻tの平滑化画像の初期値H(t)とし、ある整数kについて平滑化画像Hk−1(t)が与えられているとき、入力画像を所定の標準偏差パラメータσを有するガウス平滑化フィルタを用いて平滑化する。平滑化画像群生成部6は、ガウス平滑化フィルタを用いて平滑化した画像を、下記式(151)を満たす所定の倍率aを用いて縮小し、平滑化画像H(t)を生成する。 Smoothed image group generating section 6, the initial value H 0 of the time t of the smoothed image (t), when the smoothed image H k-1 (t) is given for some integer k, given an input image Is smoothed using a Gaussian smoothing filter having a standard deviation parameter σ g of. Smoothed image group generating section 6, an image obtained by smoothing with a Gaussian smoothing filter, reduced with a predetermined magnification a g which satisfies the following expression (151), generating a smoothed image H k (t) To do.

整数kが平滑化画像における平滑化の度合いと対応していることから、以下、整数kを平滑化係数と呼ぶ。上記の過程をk=1,2,・・・,n−2で繰り返すことにより、平滑化画像群を形成する(下記式(152))。 Since the integer k corresponds to the degree of smoothing in the smoothed image, the integer k is hereinafter referred to as a smoothing coefficient. A smoothed image group is formed by repeating the above process at k = 1, 2,..., N g −2 (the following formula (152)).

このとき、特にσ=0,a=1とすると、各平滑化画像が全て入力画像と同一となる。以上のように、平滑化画像群生成部6は、平滑化画像群を抽出し、出力する。 At this time, if σ g = 0 and a g = 1 in particular, all the smoothed images are the same as the input image. As described above, the smoothed image group generation unit 6 extracts and outputs a smoothed image group.

顕著領域事前確率画像抽出部2は、注目度映像の一のフレームである注目度画像および顕著領域画から、入力映像中の対応するフレームである入力画像の各位置が顕著領域である確率を表示する顕著領域事前確率画像を抽出する。顕著領域事前確率画像抽出部2が顕著領域事前確率画像を抽出する方法は特に限定しないが、本実施形態においては、顕著領域事前確率画像生成部21と、顕著領域事前確率画像更新部22とによって抽出する方法について説明する。   The saliency area prior probability image extraction unit 2 displays the probability that each position of the input image that is a corresponding frame in the input video is a saliency area from the attention level image and the saliency area image that are one frame of the attention level video. The saliency prior probability image to be extracted is extracted. The method by which the saliency area prior probability image extraction unit 2 extracts the saliency area prior probability image is not particularly limited, but in this embodiment, the saliency area prior probability image generation unit 21 and the saliency area prior probability image update unit 22 A method of extracting will be described.

顕著領域事前確率画像生成部21は、第1の実施形態と同様であるため説明を省略する。顕著領域事前確率画像更新部22も、第1の実施形態とほぼ同様である。但し、以下の点が第1の実施形態と異なる。   Since the saliency area prior probability image generation unit 21 is the same as that of the first embodiment, the description thereof is omitted. The saliency area prior probability image update unit 22 is substantially the same as that in the first embodiment. However, the following points are different from the first embodiment.

1.ある時刻tにおいて本処理を初めて実行する際、即ち、以降の特徴量尤度算出部3および顕著領域画像抽出部4において、平滑化係数最大の平滑化画像(下記式(153))が用いられる場合には、第1の実施形態と同様の方法で顕著領域事前確率画像(d)を更新する。 1. When this process is executed for the first time at a certain time t, that is, in the subsequent feature amount likelihood calculation unit 3 and saliency area image extraction unit 4, a smoothed image with the maximum smoothing coefficient (the following equation (153)) is obtained. When used, the saliency prior probability image (d) is updated by the same method as in the first embodiment.

2.ある時刻tにおいて本処理を再度実行する際、即ち、以降の特徴量尤度算出部3および顕著領域画像抽出部4において平滑化画像H(t)(k=n−2,n−3,・・・・,0)が用いる場合には、以下の変更を行った上で第1の実施形態と同様の処理を行う。
(1)第1の実施形態に記載の更新式のパラメータの一つであるσ、および、1時点前(時刻t−1)顕著領域事前確率の分散(下記式(154))を強制的に0に置き換える。
2. When this process is executed again at a certain time t, that is, in the subsequent feature amount likelihood calculation unit 3 and the saliency region image extraction unit 4, the smoothed image H k (t) (k = ng- 2, ng- 3,..., 0), the same processing as in the first embodiment is performed after making the following changes.
(1) σ 1 , which is one of the parameters of the update formula described in the first embodiment, and the variance of the saliency prior probability (timed t−1) before one time point (time (t−1)) (the following formula (154)) are forced Replace with 0.

(2)1時点前(時刻t−1)の顕著領域画像A(t−1)に代えて、平滑化係数が1つ大きい平滑化画像Hk+1(t)を用いて生成された顕著領域画像A(t;k+1)を用いる。
(3)顕著領域事前確率の分散(下記式(155))を更新せずに、平均ξ(x,t)のみを第1の実施形態と同様の方法で更新する。
(2) A saliency area image generated using a smoothed image H k + 1 (t) having one larger smoothing coefficient instead of the saliency area image A (t−1) one point before (time t−1). A (t; k + 1) is used.
(3) Only the average ξ 1 (x, t) is updated by the same method as in the first embodiment without updating the variance of the saliency prior probability (the following formula (155)).

3.入力として平滑化係数kの平滑化画像H(t)を用いたことを明確にするため、出力である顕著領域事前確率画像をΞ(t;k)と表記する。 3. In order to clarify that the smoothed image H k (t) having the smoothing coefficient k is used as an input, the saliency area prior probability image as an output is denoted as Ξ 1 (t; k).

特徴量尤度算出部3も、第1の実施形態とほぼ同様である。但し、以下の点が第1の実施形態と異なる。
1.入力画像に代えて平滑化画像H(t)(k=n−1,n−2,・・・,0) のうち1つを用いてもよい。このとき、時刻tにおいて本処理をj(j=1,2,・・・,n)回目に実行するときには、平滑化係数k=n−jの平滑化画像(下記式156)が用いられる。このことは、平滑化係数が大きい平滑化画像から順に用いられることを意味する。
The feature amount likelihood calculating unit 3 is substantially the same as that in the first embodiment. However, the following points are different from the first embodiment.
1. Instead of the input image, one of the smoothed images H k (t) (k = ng− 1, ng− 2,..., 0) may be used. At this time, when this process is executed j (j = 1, 2,..., N g ) times at time t, a smoothed image (equation 156 below) with a smoothing coefficient k = ng− j is used. It is done. This means that the smoothed image having the largest smoothing coefficient is used in order.

2.ある時刻tにおいて本処理を初めて実行する際、即ち、平滑化係数n−1の平滑化画像Hng−1(t)が入力として用いられる場合には、第1項目以外は第1の実施形態と同様である。
3.ある時刻tにおいて本処理を再度実行する際、即ち、平滑化係数k(k=n−2,n−1,・・・,0)の平滑化画像H(t)が入力として用いられる場合には、以下を用いる。
(1)1時点前(時刻t−1)の顕著領域画像A(t−1)に代えて、現時点(時刻t)で平滑化係数が1つ大きい平滑化画像Hk+1(t)を用いて生成された顕著領域画像A(t;k+1)を用いる。
(2)1時点前の顕著領域特徴量尤度ψ(c,t−1)に代えて、現時点で平滑化係数が1つ大きい平滑化画像Hk+1(t)を用いて生成された顕著領域特徴量尤度ψ(c,t;k+1)を用いる。
(3)1時点前の非顕著領域特徴量尤度ψ(c,t−1)に代えて、現時点で平滑化係数が1つ大きい平滑化画像Hk+1(t)を用いて生成された非顕著領域特徴量尤度ψ(c,t;k+1)を用いる。
4.入力として平滑化係数kの平滑化画像H(t)を用いたことを明確にするため、出力である顕著領域特徴量尤度をψ(t;k)、非顕著領域特徴量尤度をψ(t;k)と表記する。
2. When first run the process at a certain time t, that is, when the smoothing factor n g -1 smoothed image H ng-1 (t) is used as an input, except the first item first embodiment It is the same as the form.
3. When performing the process at a certain time t again, i.e., used as a smoothing coefficient k (k = n g -2, n g -1, ···, 0) of the smoothed image H k (t) is input If so, use:
(1) Instead of the saliency area image A (t−1) one time before (time t−1), a smoothed image H k + 1 (t) having a larger smoothing coefficient at the current time (time t) is used. The generated saliency area image A (t; k + 1) is used.
(2) The saliency generated by using the smoothed image H k + 1 (t) having one larger smoothing coefficient at the current time instead of the saliency area feature likelihood ψ 1 (c, t−1) one time before The region feature amount likelihood ψ 1 (c, t; k + 1) is used.
(3) Instead of the non-significant region feature amount likelihood ψ 2 (c, t−1) one point before, it is generated using a smoothed image H k + 1 (t) having one larger smoothing coefficient at the present time. The non-salience area feature likelihood ψ 2 (c, t; k + 1) is used.
4). In order to clarify that the smoothed image H k (t) having the smoothing coefficient k is used as an input, the saliency area feature quantity likelihood as an output is ψ 1 (t; k), and the non-salience area feature quantity likelihood. Is denoted as ψ 2 (t; k).

顕著領域画像確定部7は、平準化画像群に対し、顕著領域事前確率画像抽出部2、特徴量尤度算出部3、顕著領域画像抽出部4の処理を実行し、入力画像の顕著領域画像を確定する。即ち、顕著領域画像確定部7は、平滑化係数kの平滑化画像H(t)に対し、顕著領域事前確率画像抽出部2、特徴量尤度算出部3および顕著領域画像抽出部4を順に実行し、抽出された顕著領域画像A(t;k)を入力し、1つ前のステップで抽出された顕著領域画像A(t;k+1)からの変化がない場合には、現時点(時刻t)の入力画像に対しての最終的な顕著領域画像を確定し、この顕著領域画像A(t)= A(t;k)を出力し、変化がある場合には、kを1つ小さくして、再度顕著領域事前確率画像抽出部2に戻る。 The saliency area image determination unit 7 executes the processes of the saliency area prior probability image extraction unit 2, the feature amount likelihood calculation unit 3, and the saliency area image extraction unit 4 for the leveled image group, and the saliency area image of the input image Confirm. That is, the saliency area image determination unit 7 applies the saliency area prior probability image extraction unit 2, the feature amount likelihood calculation unit 3, and the saliency area image extraction unit 4 to the smoothed image H k (t) with the smoothing coefficient k. When the extracted saliency area image A (t; k) is input and there is no change from the saliency area image A (t; k + 1) extracted in the previous step, the current time (time The final saliency image for the input image of t) is determined, and this saliency image A (t) = A (t; k) is output. If there is a change, k is decreased by one. Then, the process returns to the saliency area prior probability image extraction unit 2 again.

図8は、顕著領域映像生成装置1000による顕著領域映像の生成結果(顕著領域の抽出結果)である。なお、入力映像は、640×480ピクセル、30〜90秒、各パラメータの数値は、σ=0.0.31、σ=0.037、M=5、λ=0.25、λ=100、σ=0、1である。図8において、第1行および第3行は入力映像、第2行は第1行に対応する顕著領域の抽出結果、第4行は第3行に対応する顕著領域の抽出結果である。図8に示すように、顕著領域映像生成装置1000によれば、適切に顕著領域が抽出され、その境界もほぼ完全である。 FIG. 8 shows a saliency area image generation result (salience area extraction result) by the saliency area image generation apparatus 1000. The input video is 640 × 480 pixels, 30 to 90 seconds, and the numerical values of the parameters are σ 1 = 0.0.31, σ 2 = 0.037, M f = 5, λ c = 0.25, λ i = 100, σ = 0, 1. In FIG. 8, the first row and the third row are the input video, the second row is the extraction result of the salient region corresponding to the first row, and the fourth row is the extraction result of the salient region corresponding to the third row. As shown in FIG. 8, according to the saliency area image generation device 1000, the saliency area is appropriately extracted, and the boundary is almost perfect.

図9は、顕著領域映像生成装置1000による方法と、顕著領域映像生成装置1000から顕著領域事前確率画像更新部22、顕著領域特徴量尤度更新部312および非顕著領域特徴量尤度更新部322を取り除いた方法、即ち、顕著領域事前確率画像および特徴量尤度を逐次的に更新しない方法との比較である。図9に示すように、顕著領域事前確率画像および特徴量尤度を逐次的に更新しない方法では分割結果がフレームによって大きく異なるが、顕著領域映像生成装置1000による方法では安定して画像分割ができる。   FIG. 9 shows a method using the saliency area video generation apparatus 1000 and a saliency area prior probability image update unit 22, a saliency area feature amount likelihood update part 312, and a non-salience area feature amount likelihood update part 322 from the saliency area image generation apparatus 1000. This is a comparison with a method in which the saliency area prior probability image and the feature amount likelihood are not sequentially updated. As shown in FIG. 9, in the method in which the saliency area prior probability image and the feature amount likelihood are not sequentially updated, the division result varies greatly depending on the frame, but the method using the saliency area image generation device 1000 can stably divide the image. .

以上の説明したように、本発明では、主に以下の2点により、上記の映像顕著領域抽出方法を実現している。   As described above, in the present invention, the above-described video salient region extraction method is realized mainly by the following two points.

(1)注目度映像抽出部1による、人間の視覚機構を模擬したモデルに基づく映像顕著性の算出、並びに、顕著領域事前確率画像生成部21および顕著領域特徴量尤度算出部31による、映像顕著性に基づく顕著領域、非顕著領域に関する事前情報の生成
(2)顕著領域事前確率画像更新部22および非顕著領域特徴量尤度算出部32による、顕著領域、非顕著領域に関する事前情報の逐次更新
物体領域・背景領域に関する事前情報が全く与えられない場合においても領域分割が可能になる。従って、物体領域、背景領域に関する事前知識がない場合でも、精度良く物体領域と背景領域を分割して、注目している領域(物体領域)を抽出することができるようになる。
(1) Image saliency calculation based on a model simulating a human visual mechanism by the attention level video extraction unit 1, and video by the saliency area prior probability image generation unit 21 and the saliency feature amount likelihood calculation unit 31 Generation of prior information on saliency areas and non-salience areas based on saliency (2) Sequential generation of prior information on saliency areas and non-salience areas by the saliency area prior probability image update unit 22 and non-salience area feature amount likelihood calculation unit 32 Even when no prior information regarding the updated object area / background area is given, the area can be divided. Therefore, even when there is no prior knowledge about the object region and the background region, the object region and the background region can be divided with high accuracy and the region of interest (object region) can be extracted.

これにより、物体領域・背景領域に関する事前情報が全く与えられない場合においても領域分割が可能になる。従って、物体領域・背景領域に関する事前知識がない場合でも、精度良く物体領域と背景領域を分割して、注目している領域(物体領域)を抽出することができるようになる。   As a result, it is possible to divide the region even when no prior information regarding the object region / background region is given. Therefore, even if there is no prior knowledge about the object region / background region, the object region and the background region can be divided with high accuracy and the region of interest (object region) can be extracted.

なお、顕著領域映像生成装置1000、1100の各処理を実行するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、当該記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより、顕著領域映像生成装置1000、1100に係る上述した種々の処理を行ってもよい。なお、ここでいう「コンピュータシステム」とは、OSや周辺機器等のハードウェアを含むものであってもよい。また、「コンピュータシステム」は、WWWシステムを利用している場合であれば、ホームページ提供環境(あるいは表示環境)も含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ROM、フラッシュメモリ等の書き込み可能な不揮発性メモリ、CD−ROM等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。   A program for executing each process of the saliency image generating apparatuses 1000 and 1100 is recorded on a computer-readable recording medium, and the program recorded on the recording medium is read into a computer system and executed. The above-described various processes related to the saliency image generation apparatuses 1000 and 1100 may be performed. Here, the “computer system” may include an OS and hardware such as peripheral devices. Further, the “computer system” includes a homepage providing environment (or display environment) if a WWW system is used. The “computer-readable recording medium” means a flexible disk, a magneto-optical disk, a ROM, a writable nonvolatile memory such as a flash memory, a portable medium such as a CD-ROM, a hard disk built in a computer system, etc. This is a storage device.

さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムが送信された場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリ(例えばDRAM(Dynamic Random Access Memory))のように、一定時間プログラムを保持しているものも含むものとする。また、上記プログラムは、このプログラムを記憶装置等に格納したコンピュータシステムから、伝送媒体を介して、あるいは、伝送媒体中の伝送波により他のコンピュータシステムに伝送されてもよい。ここで、プログラムを伝送する「伝送媒体」は、インターネット等のネットワーク(通信網)や電話回線等の通信回線(通信線)のように情報を伝送する機能を有する媒体のことをいう。また、上記プログラムは、前述した機能の一部を実現するためのものであっても良い。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル(差分プログラム)であっても良い。   Further, the “computer-readable recording medium” means a volatile memory (for example, DRAM (Dynamic DRAM) in a computer system that becomes a server or a client when a program is transmitted through a network such as the Internet or a communication line such as a telephone line. Random Access Memory)), etc., which hold programs for a certain period of time. The program may be transmitted from a computer system storing the program in a storage device or the like to another computer system via a transmission medium or by a transmission wave in the transmission medium. Here, the “transmission medium” for transmitting the program refers to a medium having a function of transmitting information, such as a network (communication network) such as the Internet or a communication line (communication line) such as a telephone line. The program may be for realizing a part of the functions described above. Furthermore, what can implement | achieve the function mentioned above in combination with the program already recorded on the computer system, and what is called a difference file (difference program) may be sufficient.

以上、この発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。   The embodiment of the present invention has been described in detail with reference to the drawings. However, the specific configuration is not limited to this embodiment, and includes designs and the like that do not depart from the gist of the present invention.

1 注目度映像抽出部 2 顕著領域事前確率画像抽出部 3 特徴量尤度算出部 4顕著領域画像抽出部 5 顕著領域映像生成部 6 平滑化画像群生成部 7 顕著領域画像確定部 11 基礎注目度画像抽出部 12 確率的基礎注目度画像抽出部 13 確率的基礎注目度母数逐次推定部 14 視線位置確率密度画像抽出部 15 視線位置確率密度映像出力部 21 顕著領域事前確率画像生成部 22 顕著領域事前確率画像更新部 31 顕著領域特徴量尤度算出部 32 非顕著領域特徴量尤度算出部 33 特徴量尤度出力部 41 顕著領域抽出グラフ生成部 42 顕著領域抽出グラフ分割部 141 視線移動状態変数更新部 142 代表視線位置更新部 143 代表視線位置重み係数算出部 144 視線位置確率密度画像出力部 145代表視線位置集合再構成部 311 顕著領域特徴量尤度生成部 312 顕著領域特徴量尤度更新部 321 非顕著領域特徴量尤度生成部 322 非顕著領域特徴量尤度更新部 1000、1100 顕著領域映像生成装置   DESCRIPTION OF SYMBOLS 1 Attention degree video extraction part 2 Saliency area prior probability image extraction part 3 Feature quantity likelihood calculation part 4 Saliency area image extraction part 5 Saliency area image generation part 6 Smoothed image group generation part 7 Saliency area image determination part 11 Basic attention degree Image extraction unit 12 Probabilistic basic attention level image extraction unit 13 Probabilistic basic attention level parameter sequential estimation unit 14 Gaze position probability density image extraction unit 15 Gaze position probability density video output unit 21 Saliency area prior probability image generation unit 22 Saliency area Prior probability image update unit 31 Saliency region feature amount likelihood calculation unit 32 Non-salience region feature amount likelihood calculation unit 33 Feature amount likelihood output unit 41 Saliency region extraction graph generation unit 42 Salience region extraction graph division unit 141 Gaze movement state variable Update unit 142 Representative eye gaze position update unit 143 Representative eye gaze position weight coefficient calculation unit 144 Eye gaze position probability density image output unit 14 Representative line-of-sight position set reconstruction unit 311 Striking region feature amount likelihood generating unit 312 Striking region feature amount likelihood updating unit 321 Non-salience region feature amount likelihood generating unit 322 Non-salience region feature amount likelihood updating unit 1000, 1100 Striking region Video generation device

Claims (14)

入力映像から、人間が注意を向けやすい度合いである注目度を示す注目度映像を抽出する注目度映像抽出過程と、
入力映像を構成する各フレームである入力画像の各位置が顕著領域である確率を示す顕著領域事前確率画像を抽出する顕著領域事前確率画像抽出過程と、
前記入力画像の顕著領域および顕著領域外の領域にそれぞれ含まれる画像特徴量の尤度を示す特徴量尤度を算出する特徴量尤度算出過程と、
前記入力画像、前記顕著領域事前確率画像および前記特徴量尤度から、前記入力画像の顕著領域を示す顕著領域画像を抽出する顕著領域画像抽出過程と、
前記各入力画像に対し、前記注目度映像抽出過程、前記顕著領域事前確率画像抽出過程、前記特徴量尤度算出過程および前記顕著領域画像抽出過程を実行して得られる前記各顕著領域画像から前記顕著領域映像を生成する顕著領域映像生成過程と
を有し、
前記顕著領域事前確率画像抽出過程は、
一の前記入力画像の各位置が顕著領域である確率を示す顕著領域事前確率画像を、前記注目度映像抽出過程によって抽出された前記注目度映像内の前記入力画像に対応する画像である注目度画像および前記顕著領域画像に基づいて抽出し、
前記特徴量尤度算出過程は、
前記特徴量尤度を、前記入力画像、前記注目度画像、前記顕著領域事前確率画像、前記顕著領域画像および前回迄に算出した前記特徴量尤度の少なくとも1つに基づいて算出する
ことを特徴とする顕著領域映像生成方法。
Attention level video extraction process for extracting attention level video indicating the level of attention that is the degree to which human attention is likely to be drawn from the input video,
A saliency area prior probability image extraction process for extracting a saliency area prior probability image indicating the probability that each position of the input image that is each frame constituting the input video is a saliency area; and
A feature amount likelihood calculating process for calculating a feature amount likelihood indicating the likelihood of the image feature amount included in each of the saliency area and the non-salience area of the input image;
A saliency area image extraction process for extracting a saliency area image indicating a saliency area of the input image from the input image, the saliency area prior probability image, and the feature amount likelihood;
From each saliency area image obtained by executing the attention level video extraction process, the saliency area prior probability image extraction process, the feature amount likelihood calculation process, and the saliency area image extraction process for each input image. A saliency image generation process for generating a saliency image,
The saliency area prior probability image extraction process includes:
The degree of attention is an image corresponding to the input image in the attention level video extracted by the attention level video extraction process, wherein the remarkable area prior probability image indicating the probability that each position of the one input image is a remarkable area Extracting based on the image and the salient region image;
The feature amount likelihood calculation process includes:
The feature quantity likelihood is calculated based on at least one of the input image, the attention degree image, the saliency area prior probability image, the saliency area image, and the feature quantity likelihood calculated up to the previous time. A saliency area image generation method.
前記注目度映像抽出過程は、
前記入力画像に基づいて、該入力画像の中で顕著な特性を持つ空間領域を表示した画像である基礎注目度画像を算出する基礎注目度画像抽出過程と、
現在の前記入力画像の各位置における顕著性を確率的な表現を用いて表示した画像である確率的基礎注目度画像を、前記基礎注目度画像抽出過程によって算出された基礎注目度画像と、前回の前記入力画像から該確率的基礎注目度画像抽出過程によって算出された前記確率的基礎注目度画像と、逐次更新され、視線位置推定に用いる第1の母数である確率的基礎注目度母数とに基づいて算出する確率的基礎注目度画像抽出過程と、
現在の前記入力画像における前記視線位置確率密度映像のフレームである視線位置確率密度画像を、前記確率的基礎注目度画像抽出過程によって算出された確率的基礎注目度画像と、前回の前記入力画像から該視線位置確率密度画像抽出過程によって算出された前記視線位置確率密度画像と、逐次更新され、視線位置推定に用いる第2の母数である視線位置確率密度母数とに基づいて算出する視線位置確率密度画像抽出過程と、
前記基礎注目度画像抽出過程と、前記確率的基礎注目度画像抽出過程と、前記視線位置確率密度画像抽出過程とを、それぞれの前記入力画像に対して順に繰り返し行うことにより算出される時系列の前記視線位置確率密度画像を前記視線位置確率密度映像として出力する視線位置確率密度映像出力過程と
を有し、
前記視線位置確率密度画像抽出過程は、
視線移動の大きさを制御する確率変数である視線移動状態変数を、前回の前記入力画像から該視線位置確率密度画像抽出過程によって算出された前記視線位置確率密度画像と、前回の前記入力画像から該視線移動状態変数更新過程によって算出された視線移動状態変数と、前記視線位置確率密度母数とに基づいて更新し、該視線移動状態変数の集合である視線移動状態変数集合を出力する視線移動状態変数更新過程と、
視線移動を考慮した代表的な視線位置を示す代表視線位置の集合である代表視線位置集合を、前記確率的基礎注目度画像抽出過程によって算出された確率的基礎注目度画像と、前回の前記入力画像から該代表視線位置更新過程によって更新された代表視線位置集合と、前記視線移動状態変数集合と、前記視線位置確率密度母数とに基づいて更新する代表視線位置更新過程と、
それぞれの前記代表視線位置に関連付けられた重みからなる代表視線位置重み係数の集合である代表視線位置重み係数集合を、前記確率的基礎注目度画像抽出過程によって算出された確率的基礎注目度画像と、前記代表視線位置更新過程によって更新された代表視線位置集合と、前記視線移動状態変数更新過程から出力された視線移動状態変数集合と、前記視線位置確率密度母数とに基づいて算出する代表視線位置重み係数算出過程と、
前記代表視線位置更新過程によって更新された代表視線位置集合と、代表視線位置重み係数算出過程によって算出された代表視線位置重み係数集合とに基づいて、前記視線位置確率密度画像を算出する視線位置確率密度画像出力過程と
を有し、前記代表視線位置集合と、前記代表視線位置重み係数集合とを含む前記視線位置確率密度画像を算出する
ことを特徴とする請求項1に記載の顕著領域映像生成方法。
The attention level video extraction process includes:
A basic attention degree image extraction process for calculating a basic attention degree image, which is an image displaying a spatial region having a remarkable characteristic in the input image based on the input image;
A stochastic basic attention image that is an image displaying prominence at each position of the current input image using a probabilistic expression, a basic attention image calculated by the basic attention image extraction process, and a previous The probabilistic basic attention degree image calculated by the probabilistic basic attention degree image extraction process from the input image and the first basic parameter that is sequentially updated and used for eye-gaze position estimation is the probabilistic basic attention degree parameter. A stochastic basic attention level image extraction process calculated based on
A gaze position probability density image that is a frame of the gaze position probability density image in the current input image is obtained from the probabilistic basic attention image calculated by the probabilistic basic attention image extraction process and the previous input image. A line-of-sight position calculated based on the line-of-sight position probability density image calculated by the line-of-sight position probability density image extraction process and a line-of-sight position probability density parameter that is sequentially updated and used for line-of-sight position estimation Probability density image extraction process;
A time series calculated by sequentially repeating the basic attention level image extraction process, the probabilistic basic attention level image extraction process, and the gaze position probability density image extraction process for each of the input images. A line-of-sight position probability density image output process for outputting the line-of-sight position probability density image as the line-of-sight position probability density image;
The line-of-sight position probability density image extraction process includes:
A line-of-sight movement state variable, which is a random variable that controls the amount of line-of-sight movement, is calculated from the line-of-sight position probability density image calculated by the line-of-sight position probability density image extraction process from the previous input image and the previous input image. Gaze movement that is updated based on the gaze movement state variable calculated by the gaze movement state variable update process and the gaze position probability density parameter and outputs a gaze movement state variable set that is a set of the gaze movement state variables State variable update process;
A representative gaze position set, which is a set of representative gaze positions indicating representative gaze positions considering gaze movement, a stochastic basic attention image calculated by the stochastic basic attention image extraction process, and the previous input A representative gaze position update process that is updated based on the representative gaze position set updated from the image by the representative gaze position update process, the gaze movement state variable set, and the gaze position probability density parameter;
A representative gaze position weight coefficient set, which is a set of representative gaze position weight coefficients composed of weights associated with each of the representative gaze positions, and a probabilistic basic attention image calculated by the probabilistic basic attention image extraction process, A representative line-of-sight calculated based on the representative line-of-sight position set updated by the representative line-of-sight position update process, the line-of-sight movement state variable set output from the line-of-sight movement state variable update process, and the line-of-sight position probability density parameter Position weighting factor calculation process,
Gaze position probability for calculating the gaze position probability density image based on the representative gaze position set updated by the representative gaze position update process and the representative gaze position weight coefficient set calculated by the representative gaze position weight coefficient calculation process The saliency area image generation according to claim 1, further comprising: calculating a gaze position probability density image including the representative gaze position set and the representative gaze position weighting coefficient set. Method.
前記顕著領域事前確率画像抽出過程は、
前記注目度画像のみを用いて前記顕著領域事前確率画像を生成する顕著領域事前確率画像生成過程と、
前記顕著領域画像を用いて前記顕著領域事前確率画像生成過程によって生成された前記顕著領域事前確率画像を更新する顕著領域事前確率画像更新過程と
から構成されることを特徴とする請求項1または請求項2の何れか1項に記載の顕著領域映像生成方法。
The saliency area prior probability image extraction process includes:
A saliency area prior probability image generation process for generating the saliency area prior probability image using only the attention level image; and
2. The saliency area prior probability image update process in which the saliency area prior probability image generated by the saliency area prior probability image generation process is updated using the saliency area image. Item 3. The saliency image generation method according to any one of Items 2 to 3.
前記特徴量尤度算出過程は、
顕著領域に含まれる画像特徴量の尤度を示す顕著領域特徴量尤度を、前記入力画像、前記顕著領域事前確率画像、前記顕著領域画像および前回迄に算出した前記顕著領域特徴量尤度のうち少なくとも1つに基づいて算出する顕著領域特徴量尤度算出過程と、
顕著領域外の領域に含まれる画像特徴量の尤度を示す非顕著領域特徴量尤度を、前記入力画像、前記顕著領域事前確率画像、前記顕著領域画像および前回迄に算出した前記非顕著領域特徴量尤度のうち少なくとも1つに基づいて算出する非顕著領域特徴量尤度算出過程と、
前記顕著領域特徴量尤度および前記非顕著領域特徴量尤度を加算して特徴量尤度として出力する特徴量尤度出力過程と
から構成されることを特徴とする請求項1から請求項3の何れか1項に記載の顕著領域映像生成方法。
The feature amount likelihood calculation process includes:
The saliency area feature amount likelihood indicating the likelihood of the image feature amount included in the saliency area is calculated based on the input image, the saliency area prior probability image, the saliency area image, and the saliency area feature amount likelihood calculated up to the previous time. A saliency area feature likelihood calculation process to be calculated based on at least one of them,
The non-salience area feature amount likelihood indicating the likelihood of the image feature amount included in the area outside the saliency area is calculated as the input image, the saliency area prior probability image, the saliency area image, and the non-salience area calculated until the previous time. A non-significant region feature amount likelihood calculation process based on at least one of the feature amount likelihoods;
4. The feature amount likelihood output process of adding the saliency region feature amount likelihood and the non-salience region feature amount likelihood and outputting the feature amount likelihood as a feature amount likelihood. The saliency area | region video generation method of any one of these.
前記顕著領域特徴量尤度算出過程は、
前記入力画像、前記顕著領域事前確率画像および前記顕著領域画像に基づいて、前記顕著領域特徴量尤度を生成する顕著領域特徴量尤度生成過程と、
前記顕著領域特徴量尤度生成過程によって生成された前記顕著領域特徴量尤度を更新する顕著領域特徴量尤度更新過程と
から構成され、
前記非顕著領域特徴量尤度算出過程は、
前記入力画像、前記顕著領域事前確率画像および前記顕著領域画像に基づいて、前記非顕著領域特徴量尤度を生成する非顕著領域特徴量尤度生成過程と、
前記非顕著領域特徴量尤度生成過程によって生成された前記非顕著領域特徴量尤度を更新する非顕著領域特徴量尤度更新過程と
から構成され、
前記顕著領域特徴量尤度更新過程は、
前記入力画像、前記顕著領域画像および前回迄に更新した更新後の前記顕著領域特徴量尤度のうち少なくとも1つに基づいて前記顕著領域特徴量尤度を更新し、
前記非顕著領域特徴量尤度更新過程は、
前記入力画像、非顕著領域画像および前回迄に更新した更新後の前記非顕著領域特徴量尤度のうち少なくとも1つに基づいて前記非顕著領域特徴量尤度を更新する
ことを特徴とする請求項4に記載の顕著領域映像生成方法。
The saliency area feature likelihood calculation process includes:
Based on the input image, the saliency area prior probability image, and the saliency area image, a saliency area feature amount likelihood generation process for generating the saliency area feature amount likelihood;
A saliency area feature amount likelihood update process for updating the saliency area feature amount likelihood generated by the saliency area feature amount likelihood generation process,
The non-significant region feature amount likelihood calculation process includes:
Based on the input image, the saliency area prior probability image, and the saliency area image, a non-salience area feature amount likelihood generation process for generating the non-salience area feature amount likelihood;
A non-salience area feature amount likelihood update process for updating the non-salience area feature amount likelihood generated by the non-salience area feature amount likelihood generation process,
The saliency area feature likelihood update process includes:
Updating the saliency area feature amount likelihood based on at least one of the input image, the saliency area image and the updated saliency area feature amount likelihood updated until the previous time,
The non-significant region feature amount likelihood update process includes:
The non-significant region feature value likelihood is updated based on at least one of the input image, the non-significant region image, and the updated non-significant region feature value likelihood updated up to the previous time. Item 5. The saliency image generation method according to Item 4.
前記入力画像を異なる解像度によってそれぞれ平滑化した複数の平滑化画像からなる平滑化画像群を生成する平滑化画像群生成過程と、
前記平準化画像群に対し、前記顕著領域事前確率画像抽出過程、前記特徴量尤度算出過程、前記顕著領域画像抽出過程を実行し、前記入力画像の前記顕著領域画像を確定する顕著領域画像確定過程と
を更に有し、
前記特徴量尤度算出過程および顕著領域画像抽出過程は、
前記入力画像に代えて前記平滑化画像を用い、
前記顕著領域映像生成過程は、
前記各入力画像に対し、前記注目度映像抽出過程、前記顕著領域事前確率画像抽出過程、前記特徴量尤度算出過程、前記顕著領域画像抽出過程、前記平滑化画像群生成過程および前記顕著領域画像確定過程を実行して得られる前記顕著領域画像から前記顕著領域映像を生成する
ことを特徴とする請求項1から請求項5の何れか1項に記載の顕著領域映像生成方法。
A smoothed image group generation process for generating a smoothed image group composed of a plurality of smoothed images obtained by smoothing the input image with different resolutions, respectively;
Perform the saliency area prior probability image extraction process, the feature amount likelihood calculation process, and the saliency area image extraction process on the leveled image group, and determine the saliency area image of the input image. And further comprising a process,
The feature amount likelihood calculation process and the saliency area image extraction process include:
Using the smoothed image instead of the input image,
The saliency image generation process includes:
For each input image, the attention level video extraction process, the saliency area prior probability image extraction process, the feature amount likelihood calculation process, the saliency area image extraction process, the smoothed image group generation process, and the saliency area image 6. The saliency area image generation method according to claim 1, wherein the saliency area image is generated from the saliency area image obtained by executing a determination process.
入力映像から、人間が注意を向けやすい度合いである注目度を示す注目度映像を抽出する注目度映像抽出部と、
入力映像を構成する各フレームである前記入力画像の各位置が顕著領域である確率を示す顕著領域事前確率画像を抽出する顕著領域事前確率画像抽出部と、
前記入力画像の顕著領域および顕著領域外の領域にそれぞれ含まれる画像特徴量の尤度を示す特徴量尤度を算出する特徴量尤度算出部と、
前記入力画像、前記顕著領域事前確率画像および前記特徴量尤度から、前記入力画像の顕著領域を示す顕著領域画像を抽出する顕著領域画像抽出部と、
前記各入力画像に対し、前記注目度映像抽出部、前記顕著領域事前確率画像抽出部、前記特徴量尤度算出部および前記顕著領域画像抽出部を実行して得られる前記各顕著領域画像から前記顕著領域映像を生成する顕著領域映像生成部と
を備え、
前記顕著領域事前確率画像抽出部は、
一の前記入力画像の各位置が顕著領域である確率を示す顕著領域事前確率画像を、前記注目度映像抽出部によって抽出された前記注目度映像内の前記入力画像に対応する画像である注目度画像および前記顕著領域画像に基づいて抽出し、
前記特徴量尤度算出部は、
前記特徴量尤度を、前記入力画像、前記注目度画像、前記顕著領域事前確率画像、前記顕著領域画像および前回迄に算出した前記特徴量尤度の少なくとも1つに基づいて算出する
ことを特徴とする顕著領域映像生成装置。
An attention level video extraction unit that extracts a degree of attention video indicating a degree of attention that is a degree that a person is likely to pay attention from an input video;
A saliency area prior probability image extraction unit that extracts a saliency area prior probability image indicating the probability that each position of the input image that is each frame constituting the input video is a saliency area; and
A feature amount likelihood calculating unit that calculates a feature amount likelihood indicating the likelihood of the image feature amount included in each of the saliency area and the area outside the saliency area of the input image;
A saliency area image extracting unit that extracts a saliency area image indicating a saliency area of the input image from the input image, the saliency area prior probability image, and the feature amount likelihood;
From each saliency area image obtained by executing the attention level video extraction unit, the saliency area prior probability image extraction unit, the feature amount likelihood calculation unit, and the saliency area image extraction unit for each input image A saliency image generating unit for generating a saliency image,
The saliency area prior probability image extraction unit includes:
Attention degree which is an image corresponding to the input image in the attention degree video extracted by the attention degree video extracting unit, a saliency area prior probability image indicating the probability that each position of the one input image is a saliency area Extracting based on the image and the salient region image;
The feature amount likelihood calculating unit includes:
The feature quantity likelihood is calculated based on at least one of the input image, the attention degree image, the saliency area prior probability image, the saliency area image, and the feature quantity likelihood calculated up to the previous time. A saliency area image generating device.
前記注目度映像抽出部は、
前記入力画像に基づいて、該入力画像の中で顕著な特性を持つ空間領域を表示した画像である基礎注目度画像を算出する基礎注目度画像抽出部と、
現在の前記入力画像の各位置における顕著性を確率的な表現を用いて表示した画像である確率的基礎注目度画像を、前記基礎注目度画像抽出部によって算出された基礎注目度画像と、前回の前記入力画像から該確率的基礎注目度画像抽出部によって算出された前記確率的基礎注目度画像と、逐次更新され、視線位置推定に用いる第1の母数である確率的基礎注目度母数とに基づいて算出する確率的基礎注目度画像抽出部と、
現在の前記入力画像における前記視線位置確率密度映像のフレームである視線位置確率密度画像を、前記確率的基礎注目度画像抽出部によって算出された確率的基礎注目度画像と、前回の前記入力画像から該視線位置確率密度画像抽出部によって算出された前記視線位置確率密度画像と、逐次更新され、視線位置推定に用いる第2の母数である視線位置確率密度母数とに基づいて算出する視線位置確率密度画像抽出部と、
前記基礎注目度画像抽出部と、前記確率的基礎注目度画像抽出部と、前記視線位置確率密度画像抽出部とを、それぞれの前記入力画像に対して順に繰り返し行うことにより算出される時系列の前記視線位置確率密度画像を前記視線位置確率密度映像として出力する視線位置確率密度映像出力部と
を備え、
前記視線位置確率密度画像抽出部は、
視線移動の大きさを制御する確率変数である視線移動状態変数を、前回の前記入力画像から該視線位置確率密度画像抽出部によって算出された前記視線位置確率密度画像と、前回の前記入力画像から該視線移動状態変数更新部によって算出された視線移動状態変数と、前記視線位置確率密度母数とに基づいて更新し、該視線移動状態変数の集合である視線移動状態変数集合を出力する視線移動状態変数更新部と、
視線移動を考慮した代表的な視線位置を示す代表視線位置の集合である代表視線位置集合を、前記確率的基礎注目度画像抽出部によって算出された確率的基礎注目度画像と、前回の前記入力画像から該代表視線位置更新部によって更新された代表視線位置集合と、前記視線移動状態変数集合と、前記視線位置確率密度母数とに基づいて更新する代表視線位置更新部と、
それぞれの前記代表視線位置に関連付けられた重みからなる代表視線位置重み係数の集合である代表視線位置重み係数集合を、前記確率的基礎注目度画像抽出部によって算出された確率的基礎注目度画像と、前記代表視線位置更新部によって更新された代表視線位置集合と、前記視線移動状態変数更新部から出力された視線移動状態変数集合と、前記視線位置確率密度母数とに基づいて算出する代表視線位置重み係数算出部と、
前記代表視線位置更新部によって更新された代表視線位置集合と、代表視線位置重み係数算出部によって算出された代表視線位置重み係数集合とに基づいて、前記視線位置確率密度画像を算出する視線位置確率密度画像出力部と
を有し、前記代表視線位置集合と、前記代表視線位置重み係数集合とを含む前記視線位置確率密度画像を算出する
ことを特徴とする請求項7に記載の顕著領域映像生成装置。
The attention level video extraction unit
A basic attention level image extraction unit that calculates a basic attention level image, which is an image displaying a spatial region having a remarkable characteristic in the input image, based on the input image;
A stochastic basic attention image, which is an image displaying prominence at each position of the current input image using a probabilistic expression, a basic attention image calculated by the basic attention image extraction unit, and a previous time The probabilistic basic attention level image calculated by the probabilistic basic attention level image extraction unit from the input image and a probabilistic basic attention level parameter that is sequentially updated and used as a first parameter for eye-gaze position estimation A stochastic basic attention level image extraction unit that is calculated based on
A gaze position probability density image, which is a frame of the gaze position probability density image in the current input image, is obtained from the probabilistic basic attention image calculated by the probabilistic basic attention image extraction unit and the previous input image. A line-of-sight position calculated based on the line-of-sight position probability density image calculated by the line-of-sight position probability density image extraction unit and a line-of-sight position probability density parameter that is sequentially updated and used for line-of-sight position estimation A probability density image extraction unit;
A time series calculated by sequentially repeating the basic attention level image extraction unit, the probabilistic basic attention level image extraction unit, and the gaze position probability density image extraction unit for each of the input images. A line-of-sight position probability density image output unit that outputs the line-of-sight position probability density image as the line-of-sight position probability density image;
The line-of-sight position probability density image extraction unit,
A line-of-sight movement state variable, which is a random variable that controls the magnitude of line-of-sight movement, is calculated from the line-of-sight position probability density image calculated by the line-of-sight position probability density image extraction unit from the previous input image and the previous input image. A line-of-sight movement that is updated based on the line-of-sight movement state variable calculated by the line-of-sight movement state variable update unit and the line-of-sight position probability density parameter and outputs a line-of-sight movement state variable set that is a set of the line-of-sight movement state variables A state variable updater;
A representative gaze position set, which is a set of representative gaze positions indicating representative gaze positions considering gaze movement, a stochastic basic attention level image calculated by the probabilistic basic attention level image extraction unit, and the previous input A representative gaze position update unit that updates the representative gaze position set updated by the representative gaze position update unit based on the gaze movement state variable set and the gaze position probability density parameter;
A representative gaze position weighting coefficient set, which is a set of representative gaze position weighting coefficients composed of weights associated with each of the representative gaze positions, and a probabilistic basic attention degree image calculated by the probabilistic basic attention degree image extraction unit, The representative line-of-sight calculated based on the representative line-of-sight position set updated by the representative line-of-sight position update unit, the line-of-sight movement state variable set output from the line-of-sight movement state variable update unit, and the line-of-sight position probability density parameter A position weight coefficient calculation unit;
The line-of-sight position probability of calculating the line-of-sight position probability density image based on the representative line-of-sight position set updated by the representative line-of-sight position update unit and the representative line-of-sight position weight coefficient set calculated by the representative line-of-sight position weight coefficient calculating unit The saliency area image generation according to claim 7, further comprising: a density image output unit, wherein the sight line position probability density image including the representative sight line position set and the representative sight position weighting coefficient set is calculated. apparatus.
前記顕著領域事前確率画像抽出部は、
前記注目度画像のみを用いて前記顕著領域事前確率画像を生成する顕著領域事前確率画像生成部と、
前記顕著領域画像を用いて前記顕著領域事前確率画像生成部によって生成された前記顕著領域事前確率画像を更新する顕著領域事前確率画像更新部と
から構成されることを特徴とする請求項7または請求項8の何れか1項に記載の顕著領域映像生成装置。
The saliency area prior probability image extraction unit includes:
A saliency area prior probability image generation unit that generates the saliency area prior probability image using only the attention level image; and
8. The saliency area prior probability image update unit configured to update the saliency area prior probability image generated by the saliency area prior probability image generation unit using the saliency area image. Item 9. The saliency area image generation device according to any one of Items 8 to 9.
前記特徴量尤度算出部は、
顕著領域に含まれる画像特徴量の尤度を示す顕著領域特徴量尤度を、前記入力画像、前記顕著領域事前確率画像、前記顕著領域画像および前回迄に算出した前記顕著領域特徴量尤度のうち少なくとも1つに基づいて算出する顕著領域特徴量尤度算出部と、
顕著領域外の領域に含まれる画像特徴量の尤度を示す非顕著領域特徴量尤度を、前記入力画像、前記顕著領域事前確率画像、前記顕著領域画像および前回迄に算出した前記非顕著領域特徴量尤度のうち少なくとも1つに基づいて算出する非顕著領域特徴量尤度算出部と、
前記顕著領域特徴量尤度および前記非顕著領域特徴量尤度を加算して特徴量尤度として出力する特徴量尤度出力部と
から構成されることを特徴とする請求項7から請求項9の何れか1項に記載の顕著領域映像生成装置。
The feature amount likelihood calculating unit includes:
The saliency area feature amount likelihood indicating the likelihood of the image feature amount included in the saliency area is calculated based on the input image, the saliency area prior probability image, the saliency area image, and the saliency area feature amount likelihood calculated up to the previous time. A saliency area feature amount likelihood calculating unit that calculates based on at least one of them,
The non-salience area feature amount likelihood indicating the likelihood of the image feature amount included in the area outside the saliency area is calculated as the input image, the saliency area prior probability image, the saliency area image, and the non-salience area calculated until the previous time. A non-significant region feature amount likelihood calculating unit that calculates based on at least one of the feature amount likelihoods;
10. The feature amount likelihood output unit configured to add the saliency region feature amount likelihood and the non-salience region feature amount likelihood and output as a feature amount likelihood. The saliency area | region video generation apparatus of any one of these.
前記顕著領域特徴量尤度算出部は、
前記入力画像、前記顕著領域事前確率画像および前記顕著領域画像に基づいて、前記顕著領域特徴量尤度を生成する顕著領域特徴量尤度生成部と、
前記顕著領域特徴量尤度生成部によって生成された前記顕著領域特徴量尤度を更新する顕著領域特徴量尤度更新部と
から構成され、
前記非顕著領域特徴量尤度算出部は、
前記入力画像、前記顕著領域事前確率画像および前記顕著領域画像に基づいて、前記非顕著領域特徴量尤度を生成する非顕著領域特徴量尤度生成部と、
前記非顕著領域特徴量尤度生成部によって生成された前記非顕著領域特徴量尤度を更新する非顕著領域特徴量尤度更新部と
から構成され、
前記顕著領域特徴量尤度更新部は、
前記入力画像、前記顕著領域画像および前回迄に更新した更新後の前記顕著領域特徴量尤度のうち少なくとも1つに基づいて前記顕著領域特徴量尤度を更新し、
前記非顕著領域特徴量尤度更新部は、
前記入力画像、非顕著領域画像および前回迄に更新した更新後の前記非顕著領域特徴量尤度のうち少なくとも1つに基づいて前記非顕著領域特徴量尤度を更新する
ことを特徴とする請求項10に記載の顕著領域映像生成装置。
The saliency area feature likelihood calculating unit
A saliency area feature amount likelihood generating unit that generates the saliency area feature amount likelihood based on the input image, the saliency area prior probability image, and the saliency area image;
A saliency area feature amount likelihood update unit configured to update the saliency area feature amount likelihood generated by the saliency area feature amount likelihood generation unit;
The non-significant region feature amount likelihood calculating unit
A non-salience area feature amount likelihood generating unit that generates the non-salience area feature amount likelihood based on the input image, the saliency area prior probability image, and the saliency area image;
A non-salience area feature amount likelihood update unit configured to update the non-salience area feature amount likelihood generated by the non-salience area feature amount likelihood generation unit;
The saliency area feature likelihood update unit includes:
Updating the saliency area feature amount likelihood based on at least one of the input image, the saliency area image and the updated saliency area feature amount likelihood updated until the previous time,
The non-significant region feature amount likelihood update unit
The non-significant region feature value likelihood is updated based on at least one of the input image, the non-significant region image, and the updated non-significant region feature value likelihood updated up to the previous time. Item 13. The saliency image generating device according to Item 10.
前記入力画像を異なる解像度によってそれぞれ平滑化した複数の平滑化画像からなる平滑化画像群を生成する平滑化画像群生成部と、
前記平準化画像群に対し、前記顕著領域事前確率画像抽出部、前記特徴量尤度算出部、前記顕著領域画像抽出部の処理を実行し、前記入力画像の前記顕著領域画像を確定する顕著領域画像確定部と
を更に備え、
前記特徴量尤度算出部および顕著領域画像抽出部は、
前記入力画像に代えて前記平滑化画像を用い、
前記顕著領域映像生成部は、
前記各入力画像に対し、前記注目度映像抽出部、前記顕著領域事前確率画像抽出部、前記特徴量尤度算出部、前記顕著領域画像抽出部、前記平滑化画像群生成部および前記顕著領域画像確定部の各処理を実行して得られる前記顕著領域画像から前記顕著領域映像を生成する
ことを特徴とする請求項7から請求項11の何れか1項に記載の顕著領域映像生成装置。
A smoothed image group generation unit that generates a smoothed image group composed of a plurality of smoothed images obtained by smoothing the input image at different resolutions;
A saliency area that executes processing of the saliency area prior probability image extraction unit, the feature amount likelihood calculation unit, and the saliency area image extraction unit for the leveled image group to determine the saliency area image of the input image An image determination unit;
The feature amount likelihood calculating unit and the salient region image extracting unit are:
Using the smoothed image instead of the input image,
The saliency area image generation unit includes:
For each input image, the attention level video extraction unit, the saliency area prior probability image extraction unit, the feature amount likelihood calculation unit, the saliency area image extraction unit, the smoothed image group generation unit, and the saliency area image 12. The saliency area image generation device according to claim 7, wherein the saliency area image is generated from the saliency area image obtained by executing each process of the determination unit.
入力映像から、人間が注意を向けやすい度合いである注目度を示す注目度映像を抽出する注目度映像抽出ステップと、
入力映像を構成する各フレームである前記入力画像の各位置が顕著領域である確率を示す顕著領域事前確率画像を抽出する顕著領域事前確率画像抽出ステップと、
前記入力画像の顕著領域および顕著領域外の領域にそれぞれ含まれる画像特徴量の尤度を示す特徴量尤度を算出する特徴量尤度算出ステップと、
前記入力画像、前記顕著領域事前確率画像および前記特徴量尤度から、前記入力画像の顕著領域を示す顕著領域画像を抽出する顕著領域画像抽出ステップと、
前記各入力画像に対し、前記注目度映像抽出ステップ、前記顕著領域事前確率画像抽出ステップ、前記特徴量尤度算出ステップおよび前記顕著領域画像抽出ステップを実行して得られる前記各顕著領域画像から前記顕著領域映像を生成する顕著領域映像生成ステップと
をコンピュータに実行させるプログラムであって、
前記顕著領域事前確率画像抽出ステップは、
一の前記入力画像の各位置が顕著領域である確率を示す顕著領域事前確率画像を、前記注目度映像抽出ステップによって抽出された前記注目度映像内の前記入力画像に対応する画像である注目度画像および前記顕著領域画像に基づいて抽出し、
前記特徴量尤度算出ステップは、
前記特徴量尤度を、前記入力画像、前記注目度画像、前記顕著領域事前確率画像、前記顕著領域画像および前回迄に算出した前記特徴量尤度の少なくとも1つに基づいて算出する
ことを特徴とするプログラム。
Attention level video extraction step for extracting an attention level video indicating a degree of attention that is the degree to which a human is likely to pay attention from an input video;
A saliency area prior probability image extraction step for extracting a saliency area prior probability image indicating a probability that each position of the input image that is each frame constituting the input video is a saliency area; and
A feature amount likelihood calculating step of calculating a feature amount likelihood indicating the likelihood of the image feature amount included in each of the saliency area and the area outside the saliency area of the input image;
A saliency area image extracting step of extracting a saliency area image indicating a saliency area of the input image from the input image, the saliency area prior probability image, and the feature amount likelihood; and
From each saliency area image obtained by executing the attention level video extraction step, the saliency area prior probability image extraction step, the feature amount likelihood calculation step and the saliency area image extraction step for each input image. A program for causing a computer to execute a saliency area image generation step for generating a saliency area image,
The saliency area prior probability image extraction step includes:
Attention degree which is an image corresponding to the input image in the attention degree video extracted by the attention degree video extraction step, a saliency area prior probability image indicating the probability that each position of the one input image is a saliency area Extracting based on the image and the salient region image;
The feature amount likelihood calculating step includes:
The feature quantity likelihood is calculated based on at least one of the input image, the attention degree image, the saliency area prior probability image, the saliency area image, and the feature quantity likelihood calculated up to the previous time. Program.
入力映像から、人間が注意を向けやすい度合いである注目度を示す注目度映像を抽出する注目度映像抽出ステップと、
入力映像を構成する各フレームである前記入力画像の各位置が顕著領域である確率を示す顕著領域事前確率画像を抽出する顕著領域事前確率画像抽出ステップと、
前記入力画像の顕著領域および顕著領域外の領域にそれぞれ含まれる画像特徴量の尤度を示す特徴量尤度を算出する特徴量尤度算出ステップと、
前記入力画像、前記顕著領域事前確率画像および前記特徴量尤度から、前記入力画像の顕著領域を示す顕著領域画像を抽出する顕著領域画像抽出ステップと、
前記各入力画像に対し、前記注目度映像抽出ステップ、前記顕著領域事前確率画像抽出ステップ、前記特徴量尤度算出ステップおよび前記顕著領域画像抽出ステップを実行して得られる前記各顕著領域画像から前記顕著領域映像を生成する顕著領域映像生成ステップと
をコンピュータに実行させるためのプログラムを記録したコンピュータ読み取り可能な記録憶媒体であって、
前記顕著領域事前確率画像抽出ステップは、
一の前記入力画像の各位置が顕著領域である確率を示す顕著領域事前確率画像を、前記注目度映像抽出ステップによって抽出された前記注目度映像内の前記入力画像に対応する画像である注目度画像および前記顕著領域画像に基づいて抽出し、
前記特徴量尤度算出ステップは、
前記特徴量尤度を、前記入力画像、前記注目度画像、前記顕著領域事前確率画像、前記顕著領域画像および前回迄に算出した前記特徴量尤度の少なくとも1つに基づいて算出する
ことを特徴とする記録媒体。
Attention level video extraction step for extracting an attention level video indicating a degree of attention that is the degree to which a human is likely to pay attention from an input video;
A saliency area prior probability image extraction step for extracting a saliency area prior probability image indicating a probability that each position of the input image that is each frame constituting the input video is a saliency area; and
A feature amount likelihood calculating step of calculating a feature amount likelihood indicating the likelihood of the image feature amount included in each of the saliency area and the area outside the saliency area of the input image;
A saliency area image extracting step of extracting a saliency area image indicating a saliency area of the input image from the input image, the saliency area prior probability image, and the feature amount likelihood; and
From each saliency area image obtained by executing the attention level video extraction step, the saliency area prior probability image extraction step, the feature amount likelihood calculation step and the saliency area image extraction step for each input image. A computer-readable storage medium storing a program for causing a computer to execute a saliency area image generation step for generating a saliency area image,
The saliency area prior probability image extraction step includes:
Attention degree which is an image corresponding to the input image in the attention degree video extracted by the attention degree video extraction step, a saliency area prior probability image indicating the probability that each position of the one input image is a saliency area Extracting based on the image and the salient region image;
The feature amount likelihood calculating step includes:
The feature quantity likelihood is calculated based on at least one of the input image, the attention degree image, the saliency area prior probability image, the saliency area image, and the feature quantity likelihood calculated up to the previous time. A recording medium.
JP2009108474A 2009-04-27 2009-04-27 Striking area image generation method, saliency area image generation apparatus, program, and recording medium Active JP5235770B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2009108474A JP5235770B2 (en) 2009-04-27 2009-04-27 Striking area image generation method, saliency area image generation apparatus, program, and recording medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2009108474A JP5235770B2 (en) 2009-04-27 2009-04-27 Striking area image generation method, saliency area image generation apparatus, program, and recording medium

Publications (2)

Publication Number Publication Date
JP2010258914A true JP2010258914A (en) 2010-11-11
JP5235770B2 JP5235770B2 (en) 2013-07-10

Family

ID=43319300

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009108474A Active JP5235770B2 (en) 2009-04-27 2009-04-27 Striking area image generation method, saliency area image generation apparatus, program, and recording medium

Country Status (1)

Country Link
JP (1) JP5235770B2 (en)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPWO2011061905A1 (en) * 2009-11-20 2013-04-04 日本電気株式会社 Object region extraction device, object region extraction method, and program
JP2013097369A (en) * 2011-11-03 2013-05-20 Kotatsu Kokusai Denshi Kofun Yugenkoshi Method for displaying background wallpaper and one or more user interface elements on display unit of electrical apparatus at the same time, computer program product for implementing method thereof and electrical apparatus implementing method thereof
JP2013142991A (en) * 2012-01-10 2013-07-22 Nippon Telegr & Teleph Corp <Ntt> Object area detection device, method and program
EP2797052A2 (en) 2013-04-26 2014-10-29 Omron Corporation Detecting a saliency region in an image
EP2797051A2 (en) 2013-04-26 2014-10-29 Omron Corporation Image saliency map determination device, method, program, and recording medium
JP2015011526A (en) * 2013-06-28 2015-01-19 Kddi株式会社 Action recognition system, method, and program, and recognizer construction system
EP2879080A2 (en) 2013-12-02 2015-06-03 Omron Corporation Image processing device and method, and computer readable medium
EP2919162A1 (en) 2014-03-13 2015-09-16 Omron Corporation Image processing apparatus and image processing method
KR20160009477A (en) * 2014-05-05 2016-01-26 시아오미 아이엔씨. Method and device for image segmentation
EP3151160A1 (en) 2015-09-29 2017-04-05 Omron Corporation Visual attention detector and visual attention detection method
US9633444B2 (en) 2014-05-05 2017-04-25 Xiaomi Inc. Method and device for image segmentation
CN109934825A (en) * 2019-03-01 2019-06-25 辽宁工程技术大学 A kind of Remote Sensing Image Segmentation of binding hierarchy gauss hybrid models and M-H

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009003615A (en) * 2007-06-20 2009-01-08 Nippon Telegr & Teleph Corp <Ntt> Attention region extraction method, attention region extraction device, computer program, and recording medium

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009003615A (en) * 2007-06-20 2009-01-08 Nippon Telegr & Teleph Corp <Ntt> Attention region extraction method, attention region extraction device, computer program, and recording medium

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPWO2011061905A1 (en) * 2009-11-20 2013-04-04 日本電気株式会社 Object region extraction device, object region extraction method, and program
JP2013097369A (en) * 2011-11-03 2013-05-20 Kotatsu Kokusai Denshi Kofun Yugenkoshi Method for displaying background wallpaper and one or more user interface elements on display unit of electrical apparatus at the same time, computer program product for implementing method thereof and electrical apparatus implementing method thereof
US8943426B2 (en) 2011-11-03 2015-01-27 Htc Corporation Method for displaying background wallpaper and one or more user interface elements on display unit of electrical apparatus at the same time, computer program product for the method and electrical apparatus implementing the method
JP2013142991A (en) * 2012-01-10 2013-07-22 Nippon Telegr & Teleph Corp <Ntt> Object area detection device, method and program
US9189861B2 (en) 2013-04-26 2015-11-17 Omron Corporation Image processing device, image processing method, program, and recording medium
EP2797052A2 (en) 2013-04-26 2014-10-29 Omron Corporation Detecting a saliency region in an image
EP2797051A2 (en) 2013-04-26 2014-10-29 Omron Corporation Image saliency map determination device, method, program, and recording medium
US9020270B2 (en) 2013-04-26 2015-04-28 Omron Corporation Image processing device, image processing method, image processing control program, and recording medium
JP2015011526A (en) * 2013-06-28 2015-01-19 Kddi株式会社 Action recognition system, method, and program, and recognizer construction system
EP2879080A2 (en) 2013-12-02 2015-06-03 Omron Corporation Image processing device and method, and computer readable medium
KR20150063924A (en) 2013-12-02 2015-06-10 오므론 가부시키가이샤 Image processing apparatus, image processing method and computer readable medium
US9367762B2 (en) 2013-12-02 2016-06-14 Omron Corporation Image processing device and method, and computer readable medium
EP2919162A1 (en) 2014-03-13 2015-09-16 Omron Corporation Image processing apparatus and image processing method
KR20150107579A (en) 2014-03-13 2015-09-23 오므론 가부시키가이샤 Image processing apparatus and image processing method
US9600746B2 (en) 2014-03-13 2017-03-21 Omron Corporation Image processing apparatus and image processing method
KR20160009477A (en) * 2014-05-05 2016-01-26 시아오미 아이엔씨. Method and device for image segmentation
KR101670004B1 (en) * 2014-05-05 2016-10-27 시아오미 아이엔씨. Method, device, program and recording medium for image segmentation
US9633444B2 (en) 2014-05-05 2017-04-25 Xiaomi Inc. Method and device for image segmentation
EP3151160A1 (en) 2015-09-29 2017-04-05 Omron Corporation Visual attention detector and visual attention detection method
KR20170038144A (en) 2015-09-29 2017-04-06 오므론 가부시키가이샤 Attention detection apparatus and attention detection method
US9904868B2 (en) 2015-09-29 2018-02-27 Omron Corporation Visual attention detector and visual attention detection method
CN109934825A (en) * 2019-03-01 2019-06-25 辽宁工程技术大学 A kind of Remote Sensing Image Segmentation of binding hierarchy gauss hybrid models and M-H
CN109934825B (en) * 2019-03-01 2023-02-28 辽宁工程技术大学 Remote sensing image segmentation method combining hierarchical Gaussian mixture model and M-H

Also Published As

Publication number Publication date
JP5235770B2 (en) 2013-07-10

Similar Documents

Publication Publication Date Title
JP5235770B2 (en) Striking area image generation method, saliency area image generation apparatus, program, and recording medium
US10600185B2 (en) Automatic liver segmentation using adversarial image-to-image network
Ostrovski et al. Autoregressive quantile networks for generative modeling
US20200401916A1 (en) Systems and methods for training generative machine learning models
Svensén et al. Robust Bayesian mixture modelling
Doucet et al. Marginal maximum a posteriori estimation using Markov chain Monte Carlo
JP5766620B2 (en) Object region detection apparatus, method, and program
Li et al. Exploring compositional high order pattern potentials for structured output learning
US8229221B2 (en) Image processing using masked restricted boltzmann machines
Papastamoulis Handling the label switching problem in latent class models via the ECR algorithm
Yu et al. Modeling spatial extremes via ensemble-of-trees of pairwise copulas
Nguyen Gaussian mixture model based spatial information concept for image segmentation
CN115526223A (en) Score-based generative modeling in a potential space
Wenzel Generative adversarial networks and other generative models
CN113284151A (en) Pancreas segmentation method and system based on deep convolutional neural network
JP6623184B2 (en) Global structure extraction apparatus, method, and program for multilayer neural network
Adams et al. Dynamic trees for image modelling
Segal Synthesising Clinically Realistic Chest X-Rays using Generative Adversarial Networks
US20190156182A1 (en) Data inference apparatus, data inference method and non-transitory computer readable medium
JP2010257344A (en) Sight line position estimating method, sight line position estimating device, program, and recording medium
Lu et al. Generalized competitive learning of Gaussian mixture models
Chin et al. Efficient data augmentation for multivariate probit models with panel data: an application to general practitioner decision making about contraceptives
JP2018195200A (en) Whole situation structure extraction device of multilayer neutral net, method, and program
JP2014149788A (en) Object area boundary estimation device, object area boundary estimation method, and object area boundary estimation program
Baumgartner et al. A new approach to image segmentation with two-dimensional hidden Markov models

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20111111

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20130326

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

Ref document number: 5235770

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20160405

Year of fee payment: 3

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350