JP6412032B2 - Space recognition device - Google Patents

Space recognition device Download PDF

Info

Publication number
JP6412032B2
JP6412032B2 JP2016016075A JP2016016075A JP6412032B2 JP 6412032 B2 JP6412032 B2 JP 6412032B2 JP 2016016075 A JP2016016075 A JP 2016016075A JP 2016016075 A JP2016016075 A JP 2016016075A JP 6412032 B2 JP6412032 B2 JP 6412032B2
Authority
JP
Japan
Prior art keywords
time
distribution
change
index value
value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2016016075A
Other languages
Japanese (ja)
Other versions
JP2017135656A (en
Inventor
徳見 修
修 徳見
黒川 高晴
高晴 黒川
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Secom Co Ltd
Original Assignee
Secom Co Ltd
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 Secom Co Ltd filed Critical Secom Co Ltd
Priority to JP2016016075A priority Critical patent/JP6412032B2/en
Publication of JP2017135656A publication Critical patent/JP2017135656A/en
Application granted granted Critical
Publication of JP6412032B2 publication Critical patent/JP6412032B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本発明は、所定空間における物理量(例えば輝度)を連続計測した計測信号(例えば輝度信号)から当該空間における変化を認識する空間認識装置に関する。   The present invention relates to a space recognition device that recognizes a change in a space from a measurement signal (for example, luminance signal) obtained by continuously measuring a physical quantity (for example, luminance) in a predetermined space.

防犯等を目的として、監視カメラが撮像した画像を背景差分処理するなどして侵入者等による変化領域を検出することが従来から行われている。このような変化領域の検出に際し、草木の揺れや照明変動による誤検出が問題とされてきた。   For the purpose of crime prevention or the like, it has been conventionally performed to detect a change area by an intruder or the like by performing background difference processing on an image captured by a surveillance camera. In detecting such a change area, erroneous detection due to shaking of a plant or illumination variation has been a problem.

従来、草木の揺れや照明変動による誤検出を低減するために背景モデルと呼ばれる手法を採る場合がある。   Conventionally, there is a case where a technique called a background model is adopted in order to reduce false detection due to shaking of a plant or illumination variation.

例えば、特許文献1に記載の侵入物体検知装置においては、複数の監視画像の各画素における画素値の出現確率分布を表す背景モデルを生成し、新たに撮像された監視画像と背景モデルを比較して各画素における相関値を求めることによって侵入者による変化領域を検出することが記載されている。   For example, in the intruding object detection device described in Patent Literature 1, a background model that represents an appearance probability distribution of pixel values in each pixel of a plurality of monitoring images is generated, and the newly captured monitoring image is compared with the background model. In addition, it is described that a change area by an intruder is detected by obtaining a correlation value in each pixel.

特開2011−123742号公報JP 2011-123742 A

しかしながら、従来技術においては、背景の画像情報側における揺らぎは軽減されるものの、背景の画像情報と比較される側(上記説明において侵入物体が写っている監視画像の側)の揺らぎの影響が依然として残り、これらの揺らぎにより変化領域の検出精度が低下していた。   However, in the prior art, fluctuations on the background image information side are reduced, but the influence of fluctuations on the side compared with the background image information (the monitoring image side in which the intruding object is reflected in the above description) still remains. In addition, the detection accuracy of the change region is reduced due to these fluctuations.

すなわち、通常用いられる監視カメラでは、各撮像素子が一定時間(例えば1/5秒間)の輝度信号の平均値を、その撮像素子に対応する画素の輝度値として出力する。そして、このように出力される輝度値は、上記一定時間内に発生した外乱やカメラ内部のノイズの影響で常に揺らいでおり、不安定である。そのため、この輝度値を用いて変化領域を検出する限り、変化領域の一部が欠ける、または変化領域の一部に背景が混入するなど揺らぎによる誤差の発生を避けることが困難であった。   In other words, in a normally used monitoring camera, each image sensor outputs an average value of luminance signals for a fixed time (for example, 1/5 second) as the luminance value of the pixel corresponding to the image sensor. The luminance value output in this way is unstable because it always fluctuates due to the disturbance generated within the predetermined time and the noise inside the camera. Therefore, as long as the change area is detected using this luminance value, it has been difficult to avoid the occurrence of an error due to fluctuation such as a part of the change area missing or a background mixed in a part of the change area.

また、従来技術において変化領域を検出するために算出する相関値は、揺らぎを含む複数の輝度値どうしを演算して得られた結果であり、スカラである。揺らぎを含む複数の輝度値どうしの演算結果である相関値においては揺らぎの影響が倍増し得る上、スカラとなってしまっているため演算結果から非揺らぎ成分を分離して揺らぎ成分による誤差を排除することが困難となっていた。相関値を算出する代わりに差分値を算出する例もあるが、差分値も相関値と同様に揺らぎの影響が倍増し得るため、非揺らぎ成分は分離困難となる。   In addition, the correlation value calculated for detecting the change region in the prior art is a result obtained by calculating a plurality of luminance values including fluctuation, and is a scalar. In the correlation value, which is the calculation result of multiple luminance values including fluctuations, the influence of fluctuation can be doubled, and since it is a scalar, non-fluctuation components are separated from the calculation results and errors due to fluctuation components are eliminated. It was difficult to do. Although there is an example in which a difference value is calculated instead of calculating a correlation value, since the influence of fluctuation can be doubled similarly to the correlation value, the non-fluctuation component is difficult to separate.

そこで、本発明は上記問題点に鑑みてなされたものであり、所定空間の物理量を計測した計測信号に揺らぎが生じていても、当該空間における変化を高精度に認識できる空間認識装置を提供することを目的とするものである。   Therefore, the present invention has been made in view of the above problems, and provides a space recognition device that can recognize a change in the space with high accuracy even if fluctuation occurs in a measurement signal obtained by measuring a physical quantity in a predetermined space. It is for the purpose.

上記した目的を達成するために、本発明に係る空間認識装置は、所定の計測面上に配置され、所定の空間の物理量を計測して計測信号を出力する計測部と、
前記計測面上に設定された小領域において、互いに異なる時間区間の前記計測信号のそれぞれから特徴量ごとの時間頻度を算出して前記時間区間ごとの特徴量分布を生成する特徴量分布生成部と、
前記時間区間が異なる前記特徴量分布間で、前記小領域における前記特徴量どうしの相違性または類似性を表す変化指標値を求め、前記変化指標値ごとの時間頻度を要素とする変化指標値分布を前記小領域ごとに生成する変化指標値分布生成部と、
前記小領域における前記変化指標値分布を用いて前記空間における変化を認識する空間変化認識部と、
を備えたことを特徴とする。
In order to achieve the above object, a space recognition device according to the present invention is arranged on a predetermined measurement surface, measures a physical quantity of a predetermined space, and outputs a measurement signal;
A feature amount distribution generation unit that calculates a time frequency for each feature amount from each of the measurement signals in different time intervals in a small region set on the measurement surface and generates a feature amount distribution for each time interval ; ,
A change index value distribution that obtains a change index value representing a difference or similarity between the feature quantities in the small region between the feature quantity distributions having different time intervals, and uses a time frequency for each change index value as an element A change index value distribution generation unit for generating each small region ,
A spatial change recognition unit that recognizes a change in the space using the change index value distribution in the small region;
It is provided with.

また、本発明に係る空間認識装置は、前記変化指標値分布生成部が、前記変化指標値ごとに、前記時間区間が異なる前記特徴量分布間における前記特徴量の前記時間頻度の組み合わせのうち当該変化指標値が算出された組み合わせごとの時間頻度の積を総和して、当該変化指標値の時間頻度を算出してもよい。 Further, in the space recognition device according to the present invention, the change index value distribution generation unit includes the time frequency combination of the feature quantities between the feature quantity distributions having different time intervals for each change index value. It is also possible to calculate the time frequency of the change index value by summing up the products of the time frequencies for each combination for which the change index value is calculated.

さらに、本発明に係る空間認識装置は、前記空間変化認識部が、前記変化指標値を変化らしさの度合いを表す評価値に変換する変換関数を用いて、前記小領域の前記変化指標値分布を前記各変化指標値に対応する前記評価値ごとの時間頻度を要素とする評価値分布に変換し、前記小領域における前記評価値分布を用いて前記空間における変化を認識してもよい。   Furthermore, in the space recognition device according to the present invention, the space change recognition unit converts the change index value distribution of the small region using a conversion function that converts the change index value into an evaluation value representing a degree of change. A change in the space may be recognized using the evaluation value distribution in the small region by converting the time frequency for each evaluation value corresponding to each change index value into an evaluation value distribution.

また、本発明に係る空間認識装置は、前記特徴量分布生成部が、
第一時間区間の前記計測信号から前記特徴量ごとの前記時間頻度を算出して短時間特徴量分布を生成する短時間特徴量分布生成部と、
前記第一時間区間よりも長い第二時間区間の前記計測信号から前記特徴量ごとの前記時間頻度を算出して長時間特徴量分布を生成する長時間特徴量分布生成部と、
を備え、
前記変化指標値分布生成部が、前記短時間特徴量分布と前記長時間特徴量分布の間で前記小領域における前記変化指標値分布を生成してもよい。
Further, in the space recognition device according to the present invention, the feature amount distribution generation unit includes:
A short-time feature value distribution generation unit that calculates the time frequency for each feature value from the measurement signal in a first time interval and generates a short-time feature value distribution;
A long-time feature amount distribution generation unit that calculates the time frequency for each feature amount from the measurement signal in a second time interval longer than the first time interval and generates a long-time feature amount distribution;
With
The change index value distribution generation unit may generate the change index value distribution in the small region between the short time feature value distribution and the long time feature value distribution.

さらに、本発明に係る空間認識装置は、前記計測部が、前記計測面を撮像面として該撮像面上に配置され、前記空間を撮像して前記計測信号としての輝度信号を出力する光電変換部からなる構成としてもよい。   Furthermore, in the space recognition apparatus according to the present invention, the measurement unit is arranged on the imaging surface with the measurement surface as an imaging surface, and the photoelectric conversion unit that images the space and outputs a luminance signal as the measurement signal It is good also as composition which consists of.

本発明の空間認識装置によれば、特徴量分布および変化指標値分布においては、時間頻度の低い揺らぎ成分と、時間頻度の高い非揺らぎ成分の間に有意な違いが生じる。そのため、非揺らぎ成分による変化を正しく認識でき、認識精度が格段に向上する。   According to the space recognition apparatus of the present invention, a significant difference occurs between a fluctuation component with a low time frequency and a non-fluctuation component with a high time frequency in the feature amount distribution and the change index value distribution. Therefore, it is possible to correctly recognize changes due to non-fluctuating components, and the recognition accuracy is greatly improved.

本発明に係る空間認識装置の概略構成を示す図である。It is a figure which shows schematic structure of the space recognition apparatus which concerns on this invention. 本発明に係る空間認識装置によって得られる分布画像の模式図である。It is a schematic diagram of the distribution image obtained by the space recognition apparatus according to the present invention. 本発明に係る空間認識装置の機能ブロック図である。It is a functional block diagram of the space recognition device concerning the present invention. (a)時間と輝度値との関係の一例を示す図である。 (b)輝度値と時間頻度との関係の一例を示す図である。(A) It is a figure which shows an example of the relationship between time and a luminance value. (B) It is a figure which shows an example of the relationship between a luminance value and time frequency. 短時間特徴量分布と長時間特徴量分布の説明図である。It is explanatory drawing of a short time feature-value distribution and a long-time feature-value distribution. (a)〜(e)監視空間内に変化があった場合の短時間特徴量分布、長時間特徴量分布、変化指標値分布、指標値/評価値変換関数、評価値分布を示す図である。(A)-(e) It is a figure which shows the short time feature-value distribution, long-time feature-value distribution, change index value distribution, index value / evaluation value conversion function, and evaluation value distribution when there is a change in the monitoring space. . (a)〜(e)監視空間内に変化がない場合の短時間特徴量分布、長時間特徴量分布、変化指標値分布、指標値/評価値変換関数、評価値分布を示す図である。(A)-(e) It is a figure which shows the short time feature-value distribution, long-time feature-value distribution, change index value distribution, index value / evaluation value conversion function, and evaluation value distribution when there is no change in the monitoring space. (a),(b)人が撮像されている撮像画像とその撮像画像における変化領域らしさの最頻値Rの一例を示す図である。(A), (b) It is a figure which shows an example of the mode value R of the change area | region likelihood in the captured image and the captured image in which the person was imaged. 本発明に係る空間認識装置における撮像部と処理部の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the imaging part and process part in the space recognition apparatus which concerns on this invention.

以下、本発明を実施するための形態について、添付した図面の図1〜9を参照しながら詳細に説明する。   Hereinafter, embodiments for carrying out the present invention will be described in detail with reference to FIGS. 1 to 9 of the accompanying drawings.

[空間認識装置の構成について]
図1に示すように、本実施の形態の空間認識装置1は、撮像部2、記憶部3、処理部4、出力部5を含んで構成され、所定空間における物理量(例えば輝度)を連続計測した計測信号(例えば輝度信号)から当該空間における変化を認識する。
[Configuration of space recognition device]
As shown in FIG. 1, the space recognition device 1 according to the present embodiment includes an imaging unit 2, a storage unit 3, a processing unit 4, and an output unit 5, and continuously measures a physical quantity (for example, luminance) in a predetermined space. The change in the space is recognized from the measured signal (for example, luminance signal).

撮像部2は、所定の空間(監視空間)を高速撮像して所定周期で図2に示すような、特徴量分布の集まりである分布画像を生成し、処理部4に出力する。なお、図2におけるx方向(水平方向)およびy方向(垂直方向)の一つ一つの矩形部分が画素に対応している。   The imaging unit 2 captures a predetermined space (monitoring space) at high speed, generates a distribution image that is a collection of feature amount distributions as shown in FIG. 2 at a predetermined cycle, and outputs the distribution image to the processing unit 4. Note that each rectangular portion in the x direction (horizontal direction) and the y direction (vertical direction) in FIG. 2 corresponds to a pixel.

撮像部2は、図3に示す複数の光電変換部20(20−1,20−2,…,20−N)を備える。複数の光電変換部20−1,20−2,…,20−Nのそれぞれは計測部であり、それぞれ監視空間に存在する物体からの反射光を光電変換して、反射光の輝度に応じた輝度信号を出力する。   The imaging unit 2 includes a plurality of photoelectric conversion units 20 (20-1, 20-2,..., 20-N) illustrated in FIG. Each of the plurality of photoelectric conversion units 20-1, 20-2,..., 20 -N is a measurement unit, and photoelectrically converts reflected light from an object existing in the monitoring space, according to the brightness of the reflected light. Outputs a luminance signal.

また、撮像部2は、図3に示す特徴量分布生成部21を備える。特徴量分布生成部21は、複数の光電変換部20−1,20−2,…,20−Nからの輝度信号を一定時間保持し、保持した所定時間幅の輝度信号から、1または複数の光電変換部20を単位とする特徴量ごとの時間頻度を特徴量分布として算出する。   The imaging unit 2 includes a feature amount distribution generation unit 21 illustrated in FIG. The feature quantity distribution generation unit 21 holds the luminance signals from the plurality of photoelectric conversion units 20-1, 20-2,... The time frequency for each feature quantity with the photoelectric conversion unit 20 as a unit is calculated as a feature quantity distribution.

本例では、1または複数の光電変換部20を単位とする特徴量分布を複数算出し、それらの集まりによる画像を分布画像としている。   In this example, a plurality of feature amount distributions with one or a plurality of photoelectric conversion units 20 as a unit are calculated, and an image obtained by the collection is used as a distribution image.

記憶部3は、ROM(Read Only Memory)、RAM(Random Access Memory)等のメモリ装置で構成され、各種プログラムや各種データを記憶する。また、記憶部3は、処理部4と接続されて処理部4との間でこれらの情報を入出力する。   The storage unit 3 includes a memory device such as a ROM (Read Only Memory) and a RAM (Random Access Memory), and stores various programs and various data. The storage unit 3 is connected to the processing unit 4 and inputs / outputs such information to / from the processing unit 4.

処理部4は、CPU(Central Processing Unit )、DSP(Digital Signal Processor)、MCU(Micro Control Unit)等の演算装置で構成される。   The processing unit 4 is configured by an arithmetic device such as a CPU (Central Processing Unit), a DSP (Digital Signal Processor), or an MCU (Micro Control Unit).

処理部4は、記憶部3と接続され、記憶部3からプログラムを読み出して実行することにより各種処理・制御手段として動作する。また、処理部4は、各種データを記憶部3に記憶させ、各種データの読み出しを行う。   The processing unit 4 is connected to the storage unit 3 and operates as various processing / control units by reading out and executing a program from the storage unit 3. The processing unit 4 stores various data in the storage unit 3 and reads out the various data.

さらに、処理部4は、撮像部2および出力部5とも接続され、撮像部2からの分布画像を処理して監視空間内の変化領域を認識するとともに、認識した変化領域が侵入者によるものか否かを判定し、侵入者による変化領域と判定した場合にアラーム信号を出力部5に出力する。   Further, the processing unit 4 is also connected to the imaging unit 2 and the output unit 5, processes the distribution image from the imaging unit 2 to recognize the change area in the monitoring space, and whether the recognized change area is caused by an intruder. If it is determined that the area has been changed by an intruder, an alarm signal is output to the output unit 5.

出力部5は、処理部4と接続され、処理部4からのアラーム信号を外部出力する。具体的には、出力部5は、通信装置で構成され、所定のネットワークを介して遠隔の警備センターにアラーム信号を伝送する。   The output unit 5 is connected to the processing unit 4 and outputs an alarm signal from the processing unit 4 to the outside. Specifically, the output unit 5 includes a communication device and transmits an alarm signal to a remote security center via a predetermined network.

さらに、図3の機能ブロック図を参照しながら各部の構成について説明する。   Further, the configuration of each unit will be described with reference to the functional block diagram of FIG.

複数の光電変換部20−1,20−2,…,20−Nは、それぞれが撮像素子、電子シャッター、増幅器、A/D変換器を含む。電子シャッター、増幅器およびA/D変換器は2以上の光電変換部20により共有させてもよい。   Each of the plurality of photoelectric conversion units 20-1, 20-2,..., 20-N includes an image sensor, an electronic shutter, an amplifier, and an A / D converter. The electronic shutter, amplifier, and A / D converter may be shared by two or more photoelectric conversion units 20.

撮像素子は、CCD(Charge Coupled Devices)またはC−MOS(Complementary Metal Oxide Semiconductor )などで構成され、監視空間の可視光線の強度を計測して、計測した強度を電圧値とする輝度信号を出力する。   The imaging device is composed of a CCD (Charge Coupled Devices) or a C-MOS (Complementary Metal Oxide Semiconductor), and measures the intensity of visible light in the monitoring space and outputs a luminance signal having the measured intensity as a voltage value. .

例えば、各光電変換部20は、1/500秒ごとに電子シャッターを作動させて、作動時の撮像素子の電圧を増幅器で増幅し、A/D変換器でサンプリングすることにより、1/500秒ごとに輝度値を0〜255の256階調に離散化したディジタルの輝度信号を出力する。   For example, each photoelectric conversion unit 20 operates an electronic shutter every 1/500 seconds, amplifies the voltage of the imaging device at the time of operation with an amplifier, and samples with an A / D converter, thereby 1/500 seconds. A digital luminance signal in which the luminance value is discretized into 256 gradations from 0 to 255 is output every time.

複数の撮像素子は、所定の撮像面上に配置されて面をなし、当該面に監視空間を投影する。例えば、撮像素子を平面の基板上に、基板のX方向に等間隔で640個、Y方向に等間隔で480個というように、アレー状に配置し、1/500秒ごとに640×480画素の濃淡画像を撮像する。   The plurality of image sensors are arranged on a predetermined imaging surface to form a surface, and project the monitoring space onto the surface. For example, the image sensors are arranged in an array on a flat substrate such that 640 elements are equally spaced in the X direction of the substrate and 480 elements are equally spaced in the Y direction, and 640 × 480 pixels every 1/500 second. The gray image is taken.

また、複数の撮像素子は、それぞれがR、G、Bの各帯域に適合したカラーフィルタを備えた3つの撮像素子を単位とする撮像素子の組を、平面の基板上に、X方向に等間隔で640組、Y方向に等間隔で480組というように、アレー状に配置し、640×480画素のカラー画像を撮像することもできる。   In addition, a plurality of image pickup devices each include a set of image pickup devices each having three image pickup devices each having a color filter suitable for each of the R, G, and B bands on a flat substrate in the X direction. 640 sets at intervals and 480 sets at equal intervals in the Y direction can be arranged in an array to pick up a color image of 640 × 480 pixels.

さらに、RGBの輝度信号をHSVなど他の表色系に変換して出力してもよい。   Further, RGB luminance signals may be converted into another color system such as HSV and output.

光電変換部20の撮像周期(計測周期)は、検出対象とする変化よりも十分短ければよく、例えば人による変化を主たる検出対象とする場合は、1/100秒、または1/200秒などでもよい。   The imaging cycle (measurement cycle) of the photoelectric conversion unit 20 only needs to be sufficiently shorter than a change to be detected. For example, when a change by a person is a main detection target, it may be 1/100 second or 1/200 second. Good.

また、階調数、光電変換部20の個数も上記例に限らず、検出対象や空間の大きさに応じた設定とすることができる。   The number of gradations and the number of photoelectric conversion units 20 are not limited to the above example, and can be set according to the detection target and the size of the space.

以上のように、複数の光電変換部20−1,20−2,…,20−Nは、所定の撮像面上に配置され、所定の空間を撮像してそれぞれが輝度信号を出力する。すなわち、複数の計測部は、所定の計測面上に配置され、所定の空間の物理量を計測してそれぞれが計測信号を出力する。   As described above, the plurality of photoelectric conversion units 20-1, 20-2,..., 20-N are arranged on a predetermined imaging surface, image a predetermined space, and each output a luminance signal. That is, the plurality of measurement units are arranged on a predetermined measurement surface, measure physical quantities in a predetermined space, and each output a measurement signal.

なお、本実施の形態では、撮像部2が特徴量分布生成部21を含んで機能し、処理部4が変化指標値分布生成部22、空間変化認識部24、侵入検知部25を含んで機能し、記憶部3が指標値/評価値変換関数記憶部23を含んで機能する。   In the present embodiment, the imaging unit 2 functions including the feature amount distribution generation unit 21, and the processing unit 4 functions including the change index value distribution generation unit 22, the space change recognition unit 24, and the intrusion detection unit 25. The storage unit 3 functions including the index value / evaluation value conversion function storage unit 23.

特徴量分布生成部21は、バッファ、演算回路を含む。バッファは、例えばリングバッファで構成され、複数の光電変換部20−1,20−2,…,20−Nからの輝度信号を一定時間保持し、撮像周期ごとの輝度値を循環記憶する。   The feature amount distribution generation unit 21 includes a buffer and an arithmetic circuit. The buffer is composed of, for example, a ring buffer, holds the luminance signals from the plurality of photoelectric conversion units 20-1, 20-2,..., 20-N for a certain period, and circulates and stores the luminance values for each imaging cycle.

本例では、予め定めた頻度算出時間幅以上のバッファを光電変換部20の数だけ備える。   In this example, the number of buffers equal to or greater than a predetermined frequency calculation time width is provided as many as the number of photoelectric conversion units 20.

例えば、撮像周期が1/500秒、後述する第二時間幅が1分間であれば、バッファ長が30000以上のバッファを光電変換部20の数だけ用意する。   For example, if the imaging cycle is 1/500 seconds and the second time width described later is 1 minute, buffers having a buffer length of 30000 or more are prepared for the number of photoelectric conversion units 20.

演算回路は、バッファに保持した頻度算出時間幅の輝度信号から特徴量ごとの時間頻度を頻度算出周期で算出して、特徴量分布を生成し、変化指標値分布生成部22に出力する。   The arithmetic circuit calculates the time frequency for each feature amount from the luminance signal of the frequency calculation time width held in the buffer in the frequency calculation cycle, generates a feature amount distribution, and outputs the feature amount distribution to the change index value distribution generation unit 22.

例えば演算回路は、撮像部2の撮像面上に予め設定された複数の小領域のそれぞれにおいて特徴量の正規化ヒストグラムを算出する。   For example, the arithmetic circuit calculates a normalized histogram of feature amounts in each of a plurality of small regions set in advance on the imaging surface of the imaging unit 2.

例えば、小領域を上記濃淡画像の画素、特徴量を上記濃淡画像の輝度値とすることができる。つまり、光電変換部20と小領域が1対1対応し、特徴量はスカラである。   For example, the small area can be the pixel of the gray image, and the feature amount can be the luminance value of the gray image. That is, the photoelectric conversion unit 20 and the small area have a one-to-one correspondence, and the feature amount is a scalar.

また、例えば、小領域を上記カラー画像の画素、特徴量をR、G、Bの輝度値としてもよい。この場合、光電変換部20と小領域が3対1対応し、特徴量は3次元ベクトルである。   Further, for example, the small area may be the pixel of the color image, and the feature amount may be the R, G, and B luminance values. In this case, the photoelectric conversion unit 20 and the small region have a one-to-one correspondence, and the feature amount is a three-dimensional vector.

さらに、例えば、上記濃淡画像またはカラー画像における近傍の複数画素(2×2画素、3×3画素など)を小領域としてもよい。この場合、解像度変換を行いつつ、各小領域における時間頻度算出の元データを増やせる。   Furthermore, for example, a plurality of neighboring pixels (2 × 2 pixels, 3 × 3 pixels, etc.) in the grayscale image or color image may be set as a small region. In this case, the original data for calculating the time frequency in each small area can be increased while performing resolution conversion.

頻度算出周期は、後述する第一時間幅1/5秒に対して、1/5秒であり、つまりオーバーラップなしの状態である。   The frequency calculation cycle is 1/5 second with respect to a first time width 1/5 second described later, that is, there is no overlap.

なお、頻度算出周期を第一時間幅よりも短く設定してオーバーラップを設けてもよいし、周期を第一時間幅よりも長く設定して間引きしてもよい。   Note that the frequency calculation cycle may be set shorter than the first time width and an overlap may be provided, or the cycle may be set longer than the first time width and thinned.

以上のように、特徴量分布生成部21は、計測面上(撮像部2の撮像面上)に設定された複数の小領域のそれぞれにおいて、順次、所定時間幅の計測信号(輝度信号)から特徴量ごとの時間頻度を算出して特徴量分布を生成する。   As described above, the feature amount distribution generation unit 21 sequentially starts from a measurement signal (luminance signal) having a predetermined time width in each of a plurality of small regions set on the measurement surface (on the image pickup surface of the image pickup unit 2). A time frequency for each feature amount is calculated to generate a feature amount distribution.

ここで、図4(a)はある画素における輝度信号の時間と輝度値の関係の一例を示す。図4(a)に示すように、輝度値が0〜255の範囲で時間の経過とともに変化したものとする。頻度算出時間幅Tの輝度値の数が例えば100個である場合、各輝度値の度数を計数し、計数結果のそれぞれを100で除して時間頻度を算出すると、当該画素における輝度値と時間頻度は図4(b)に示す正規化ヒストグラムのような関係になる。   Here, FIG. 4A shows an example of the relationship between the luminance signal time and the luminance value in a certain pixel. As shown in FIG. 4A, it is assumed that the luminance value changes with time in the range of 0-255. When the number of luminance values of the frequency calculation time width T is 100, for example, the frequency value of each luminance value is counted, and when the time frequency is calculated by dividing each of the counting results by 100, the luminance value and time in the pixel are calculated. The frequency has a relationship like a normalized histogram shown in FIG.

特徴量分布生成部21は、図3に示すように、短時間特徴量分布生成部210および長時間特徴量分布生成部211を備え、時間幅の異なる2種類の特徴量分布を生成する。   As shown in FIG. 3, the feature amount distribution generation unit 21 includes a short-time feature amount distribution generation unit 210 and a long-time feature amount distribution generation unit 211, and generates two types of feature amount distributions having different time widths.

短時間特徴量分布生成部210は、複数の小領域のそれぞれにおいて、図5に示す第一時間幅Taの輝度信号から特徴量ごとの時間頻度を算出して短時間特徴量分布を生成し、変化指標値分布生成部22に出力する。   The short time feature amount distribution generation unit 210 calculates a time frequency for each feature amount from the luminance signal having the first time width Ta shown in FIG. 5 in each of the plurality of small regions, and generates a short time feature amount distribution. It outputs to the change index value distribution generation unit 22.

例えば、1/5秒ごとに、各画素について、バッファから最新1/5秒間の輝度信号(時系列の輝度値)を読み出し、読み出した輝度信号から輝度値ごとの時間頻度を算出することにより、短時間特徴量分布を生成する。   For example, by reading out the latest 1/5 second luminance signal (time series luminance value) from the buffer for each pixel every 1/5 second, and calculating the time frequency for each luminance value from the read luminance signal, A short-time feature distribution is generated.

長時間特徴量分布生成部211は、複数の小領域のそれぞれにおいて、図5に示す第一時間幅Taよりも長い第二時間幅Tbの輝度信号から特徴量ごとの時間頻度を算出して長時間特徴量分布を生成し、変化指標値分布生成部22に出力する。   The long-time feature value distribution generation unit 211 calculates the time frequency for each feature value from the luminance signal having the second time width Tb longer than the first time width Ta shown in FIG. A temporal feature quantity distribution is generated and output to the change index value distribution generation unit 22.

例えば、1/5秒ごとに、各画素について、バッファから最新1分間の輝度信号(時系列の輝度値)を読み出し、読み出した輝度信号から輝度値ごとの時間頻度を算出することにより、長時間特徴量分布を生成する。   For example, by reading out the latest 1 minute luminance signal (time series luminance value) from the buffer every 1/5 second and calculating the time frequency for each luminance value from the read luminance signal, Generate a feature distribution.

なお、短時間特徴量分布を生成する時間区間が長時間特徴量分布を生成する時間区間に含まれていても構わない。   The time interval for generating the short-time feature value distribution may be included in the time interval for generating the long-time feature value distribution.

各特徴量分布においては、非揺らぎ成分である真の要素(真の特徴量)の時間頻度が有意に高い値を示し、これと比較して、揺らぎ成分である偽の要素(偽の特徴量)の時間頻度は十分に低いものとなる。   In each feature quantity distribution, the time frequency of the true element (true feature quantity) that is a non-fluctuating component shows a significantly high value, and compared with this, a false element (fake feature quantity) that is a fluctuation component ) Time frequency is sufficiently low.

例えば、突発的なノイズが発生し、その振幅が大きかったとしても(つまり真の輝度値から大きく外れた値だったとしても)、特徴量分布においてはその時間頻度が十分に低いものとなり、真の輝度値とノイズを区別可能となる。   For example, even if sudden noise occurs and its amplitude is large (that is, even if it is a value greatly deviating from the true luminance value), the time frequency in the feature quantity distribution is sufficiently low, and true It is possible to distinguish the luminance value and noise.

そして、短時間特徴量分布では、建物等の静止物体や草木等の振動/回転物体といった背景物体、および人や車両のような移動物体について、非揺らぎ成分である真の特徴量の時間頻度が有意に高い値を示し、揺らぎ成分である偽の特徴量の時間頻度は十分に低いものとなる。   In the short-time feature quantity distribution, the time frequency of the true feature quantity as a non-fluctuating component is obtained for background objects such as stationary objects such as buildings, vibration / rotation objects such as plants, and moving objects such as people and vehicles. It shows a significantly high value, and the time frequency of the fake feature amount that is a fluctuation component is sufficiently low.

これに対し、長時間特徴量分布では、背景物体について、非揺らぎ成分である真の特徴量の時間頻度が有意に高い値を示し、移動物体、および揺らぎ成分である偽の特徴量の時間頻度は十分に低いものとなる。特に、長時間特徴量分布によって、静止物体のみならず振動/回転物体についても背景物体の非揺らぎ成分として扱うことができる。   On the other hand, in the long-time feature amount distribution, the time frequency of the true feature amount that is a non-fluctuation component is significantly high for the background object, and the time frequency of the false feature amount that is a moving object and the fluctuation component Is sufficiently low. In particular, due to the long-time feature distribution, not only a stationary object but also a vibration / rotation object can be treated as a non-fluctuating component of a background object.

変化指標値分布生成部22は、特徴量分布生成部21が順次生成した特徴量分布間で小領域における、特徴量どうしの相違性または類似性を表す変化指標値を求め、変化指標値ごとの時間頻度を要素とする変化指標値分布を算出する。   The change index value distribution generation unit 22 obtains a change index value representing a difference or similarity between feature quantities in a small area between the feature quantity distributions sequentially generated by the feature quantity distribution generation unit 21, and determines the change index value for each change index value. A change index value distribution having time frequency as an element is calculated.

変化指標値は、変化領域を認識するための指標となる値である。例えば、特徴量を濃淡画像の輝度値、変化指標値を輝度差とすることができる。その場合、変化指標値である輝度差は、特徴量どうしの相違性を表す値である。したがって、輝度差が大きな小領域どうしほど相違性が高く、変化領域である確度が高いこととなり、輝度差が0または0近傍である小領域どうしは相違性が低く、変化領域である確度が低いこととなる。   The change index value is a value serving as an index for recognizing the change area. For example, the feature value can be a luminance value of a grayscale image, and the change index value can be a luminance difference. In this case, the luminance difference that is the change index value is a value that represents the difference between the feature amounts. Accordingly, the difference between the small areas having a large luminance difference is high, and the probability of being a change area is high. The difference between the small areas having a luminance difference of 0 or near 0 is low, and the accuracy of being a change area is low. It will be.

変化指標値分布は、時間頻度を加味して変化領域を認識するための基礎となる情報である。具体的には、変化指標値ごとに、特徴量分布間における特徴量の時間頻度の組み合わせのうち当該変化指標値が算出された組み合わせごとの時間頻度の積を総和して、当該変化指標値の時間頻度を算出する。   The change index value distribution is information that serves as a basis for recognizing a change area in consideration of time frequency. Specifically, for each change index value, the product of the time frequency for each combination for which the change index value is calculated among the combinations of the time frequency of the feature quantities between the feature quantity distributions is summed, and the change index value Calculate the time frequency.

例えば、特徴量を濃淡画像の輝度値、変化指標値を輝度差とした場合、輝度差ごとに、輝度値分布間における輝度値の時間頻度の組み合わせのうち当該輝度差が算出された組み合わせごとの時間頻度の積を総和して、当該輝度差の時間頻度を算出し、輝度差分布を求める。   For example, when the feature value is the luminance value of the grayscale image and the change index value is the luminance difference, for each luminance difference, for each combination for which the luminance difference is calculated among the luminance frequency temporal combinations between the luminance value distributions. The product of the time frequency is summed to calculate the time frequency of the brightness difference, and the brightness difference distribution is obtained.

また、振動/回転物体を変化領域と誤認識しないように、短時間特徴量分布生成部210が生成した短時間特徴量分布と、長時間特徴量分布生成部211が生成した長時間特徴量分布の間で小領域ごとに変化指標値分布を算出するのが好ましい。   Further, the short-time feature value distribution generated by the short-time feature value distribution generation unit 210 and the long-time feature value distribution generated by the long-time feature value distribution generation unit 211 so that the vibration / rotation object is not mistakenly recognized as a change region. It is preferable to calculate the change index value distribution for each small region between the two.

例えば、短時間輝度値分布を構成する任意の画素(x,y)における輝度値iの時間頻度Pa(i)と、長時間輝度値分布の対応画素(x,y)における輝度値jの時間頻度Pb(j)との輝度差分布Ps(|i−j|)は下記式(1)のようになる。   For example, the time frequency Pa (i) of the luminance value i at an arbitrary pixel (x, y) constituting the short-time luminance value distribution and the time of the luminance value j at the corresponding pixel (x, y) of the long-time luminance value distribution The luminance difference distribution Ps (| i−j |) with the frequency Pb (j) is expressed by the following formula (1).

Ps(|i−j|)= ∫j・∫i・Pb(j)・Pa(i)・di・dj…式(1)   Ps (| i−j |) = ∫j · ∫i · Pb (j) · Pa (i) · di · dj Expression (1)

なお、輝度値i,jおよびiとjの輝度差|i−j|の値域はいずれも0〜255である。また、∫j∫ididjは、輝度差が|i−j|である輝度値iとjの組み合わせについて時間頻度の積を総和することを意味する。   The value ranges of the luminance values i, j and the luminance difference | i−j | between i and j are all 0 to 255. Further, ∫j∫ididj means summing up the products of time frequencies for the combination of luminance values i and j whose luminance difference is | i−j |.

ここで、輝度差分布Psのうち輝度差が0の要素Ps(0)について説明する。
長時間輝度値分布側の輝度値j=0,1,…,255に対して、輝度差|i−j|が0となる短時間輝度値分布側の輝度値iはそれぞれ0,1,…,255であるから、長時間輝度値分布側の時間頻度Pb(0)と短時間輝度値分布側の時間頻度Pa(0)の積、長時間輝度値分布側の時間頻度Pb(1)と短時間輝度値分布側の時間頻度Pa(1)の積、…、長時間輝度値分布側の時間頻度Pb(255)と短時間輝度値分布側の時間頻度Pa(255)の積を総和して要素Ps(0)が求まる。
Here, the element Ps (0) having a luminance difference of 0 in the luminance difference distribution Ps will be described.
The luminance value i on the short-time luminance value distribution side where the luminance difference | i−j | is 0 for the luminance value j = 0, 1,..., 255 on the long-time luminance value distribution side is 0, 1,. , 255, the product of the time frequency Pb (0) on the long-time luminance value distribution side and the time frequency Pa (0) on the short-time luminance value distribution side, and the time frequency Pb (1) on the long-time luminance value distribution side, The product of the time frequency Pa (1) on the short-time luminance value distribution side, ..., the product of the time frequency Pb (255) on the long-time luminance value distribution side and the time frequency Pa (255) on the short-time luminance value distribution side Thus, the element Ps (0) is obtained.

次に、輝度差分布Psのうち輝度差が1の要素Ps(1)について説明する。
長時間輝度値分布側の輝度値j=0,1,2,…,255に対して、輝度差|i−j|が1となる短時間輝度値分布側の輝度値iはそれぞれ1,0と2,1と3,…,254であるから、Pb(0)とPa(1)の積、Pb(1)とPa(0)の積、Pb(1)とPa(2)の積、Pb(2)とPa(1)の積、Pb(2)とPa(3)の積、…、Pb(255)とPa(254)の積を総和して要素Ps(1)が求まる。Ps(2)からPs(255)についても同様に求める。
Next, the element Ps (1) having a luminance difference of 1 in the luminance difference distribution Ps will be described.
The luminance value i on the short-time luminance value distribution side where the luminance difference | i−j | is 1 for the luminance value j = 0, 1, 2,. , 2, 1 and 3, ..., 254, the product of Pb (0) and Pa (1), the product of Pb (1) and Pa (0), the product of Pb (1) and Pa (2), The product of Pb (2) and Pa (1), the product of Pb (2) and Pa (3),..., The product of Pb (255) and Pa (254) is summed to obtain the element Ps (1). The same applies to Ps (2) to Ps (255).

ここで、図6は人が現れた場合など、監視空間内に変化があった場合における短時間特徴量分布、長時間特徴量分布、変化指標値分布、指標値/評価値変換関数、評価値分布の一例を示す。また、図7は監視空間内に変化が無い場合における短時間特徴量分布、長時間特徴量分布、変化指標値分布、指標値/評価値変換関数、評価値分布の一例を示す。図7を見ると、監視空間内に変化が無くても輝度信号の揺らぎは発生していることが判る。また、図6と図7を比較すると、監視空間内に変化があった場合には、短時間特徴量分布と長時間特徴量分布においてピークの時間頻度を示す輝度値が異なっており、監視空間内に変化が無い場合には、短時間特徴量分布と長時間特徴量分布においてピークの時間頻度を示す輝度値が略一致している。   Here, FIG. 6 shows a short-time feature value distribution, a long-time feature value distribution, a change index value distribution, an index value / evaluation value conversion function, an evaluation value when there is a change in the monitoring space, such as when a person appears. An example of distribution is shown. FIG. 7 shows an example of the short-time feature value distribution, long-time feature value distribution, change index value distribution, index value / evaluation value conversion function, and evaluation value distribution when there is no change in the monitoring space. FIG. 7 shows that the luminance signal fluctuates even if there is no change in the monitoring space. Further, comparing FIG. 6 and FIG. 7, when there is a change in the monitoring space, the brightness value indicating the peak time frequency is different in the short-time feature distribution and the long-time feature distribution. When there is no change, the luminance values indicating the peak time frequency in the short-time feature amount distribution and the long-time feature amount distribution substantially coincide with each other.

変化指標値分布においては、小領域間に真の変化があれば真に相違性がある要素または真に類似性がある要素の時間頻度が有意に高い値を示し、これと比較して、揺らぎ成分の混入した偽の相違性または偽の類似性が現れている要素の時間頻度は十分に低いものとなる。   In the change index value distribution, if there is a true change between small areas, the time frequency of elements that are truly different or elements that are truly similar shows a significantly high value. The time frequency of an element showing a false difference or false similarity mixed with components is sufficiently low.

指標値/評価値変換関数記憶部23は、変化指標値を、変化らしさを表す評価値に変換する指標値/評価値変換関数を予め記憶する。   The index value / evaluation value conversion function storage unit 23 stores in advance an index value / evaluation value conversion function for converting the change index value into an evaluation value representing the likelihood of change.

例えば、評価値Rは小領域の変化領域らしさを表し、その値域は0〜1とする。また、変化領域らしさが高い変化指標値ほど高い評価値Rに対応させる。さらに、変化領域らしさが全く無ければ評価値R=0、確実に変化領域らしければ評価値R=1、あいまいさがあれば0<R<1とする。   For example, the evaluation value R represents the change area like small area, and the value range is 0 to 1. Further, a change index value having a higher change area likelihood is associated with a higher evaluation value R. Further, the evaluation value R = 0 if there is no change area, the evaluation value R = 1 if the change area is sure, and 0 <R <1 if there is ambiguity.

これを実現する連続関数の一例として、0近傍の輝度差kには揺らぎによるあいまいさがあることを考慮し、k=0であれば評価値R=0、0<k<TkであればRをkについて0<R<1の単純増加関数とし、k≧TkであればR=1とする。   As an example of a continuous function for realizing this, considering that there is ambiguity due to fluctuation in the luminance difference k near 0, if k = 0, the evaluation value R = 0, and if 0 <k <Tk, R Is a simple increase function of 0 <R <1 for k, and R = 1 if k ≧ Tk.

なお、Tkは想定する揺らぎの大きさに応じて設定する。例えば、事前の実験に基づいて、値域が0〜255である輝度差kに対してTk=40などと定めておく。またはこれを実現するステップ関数の一例として、輝度差kが0であれば評価値R=0、輝度差kが0より大きければR=1とする。   Note that Tk is set according to the assumed fluctuation magnitude. For example, based on a prior experiment, Tk = 40 or the like is determined for the luminance difference k whose value range is 0 to 255. Alternatively, as an example of a step function for realizing this, if the luminance difference k is 0, the evaluation value R = 0, and if the luminance difference k is greater than 0, R = 1.

つまり、上記例はいずれも、大きな輝度差ほど変化領域らしさを表す評価値を高く割り当てる変換関数である。   That is, all of the above examples are conversion functions that assign a higher evaluation value that represents the change area as the luminance difference increases.

空間変化認識部24は、各小領域における変化指標値分布に基づいて監視空間における変化を認識する。例えば、空間変化認識部24は監視空間における変化領域を認識する。   The space change recognition unit 24 recognizes a change in the monitoring space based on the change index value distribution in each small area. For example, the space change recognition unit 24 recognizes a change region in the monitoring space.

具体的には、空間変化認識部24は、まず、指標値/評価値変換関数記憶部23に記憶されている変換関数を読み出し、下記式(2)に従い変化指標値の時間頻度を当該変化指標値に対応する評価値の時間頻度に変換することによって、各小領域の変化指標値分布を評価値分布に変換する。   Specifically, the space change recognition unit 24 first reads the conversion function stored in the index value / evaluation value conversion function storage unit 23, and calculates the time frequency of the change index value according to the following equation (2). The change index value distribution of each small area is converted into an evaluation value distribution by converting the time frequency of the evaluation value corresponding to the value.

P(R)=∫k・Ps(R(k))・dk …式(2)   P (R) = ∫k · Ps (R (k)) · dk (2)

但し、kは変化指標値、Psは変化指標値の時間頻度である。また、∫k dkは変化指標値kについての総和を意味する。   Here, k is the change index value, and Ps is the time frequency of the change index value. Further, ∫k dk means the sum total regarding the change index value k.

変換関数R(k)により変化指標値は評価値に変換され、割り当てられた評価値が同一である変化指標値の時間頻度が積算されて当該評価値の時間頻度に変換される。   The change index value is converted into an evaluation value by the conversion function R (k), and the time frequencies of the change index values having the same assigned evaluation value are integrated and converted into the time frequency of the evaluation value.

例えばkは輝度差であり、k=|i−j|となる。
つまり、変換関数を上記連続関数とした例では、P(0)=Ps(R(0))、…、P(0.05)=Ps(R(1))、…、P(1)=Ps(R(Tk))+Ps(R(Tk+1))+…+Ps(R(255))となる。
For example, k is a luminance difference, and k = | i−j |.
That is, in the example in which the conversion function is the continuous function, P (0) = Ps (R (0)),..., P (0.05) = Ps (R (1)),. Ps (R (Tk)) + Ps (R (Tk + 1)) +... + Ps (R (255)).

空間変化認識部24は、次に、最頻値(最も高い時間頻度が算出された評価値)が0より大きな評価値分布が算出された小領域を変化領域と認識する。   Next, the space change recognition unit 24 recognizes, as a change area, a small area in which an evaluation value distribution with a mode value (an evaluation value for which the highest time frequency is calculated) is greater than 0 is calculated.

また、変換関数を上記連続関数とした場合は、輝度差kが0近傍であり揺らぎによるあいまいさを含んだ輝度差が評価値0<R<1の領域に分散するため、不当にR=1に積算されることを防げる。そのため、非揺らぎ成分による変化指標値を揺らぎ成分による変化指標値から精度良く分離して変化領域を高精度に認識できる。   Further, when the conversion function is the continuous function, the luminance difference k is in the vicinity of 0 and the luminance difference including the ambiguity due to the fluctuation is dispersed in the region where the evaluation value 0 <R <1, and therefore, R = 1 is inappropriately set. Can be prevented from being accumulated. Therefore, it is possible to separate the change index value based on the non-fluctuation component from the change index value based on the fluctuation component with high accuracy and recognize the change region with high accuracy.

または、変換関数を上記ステップ関数とした場合は、輝度差が0である時間頻度が、輝度差が0でない時間頻度よりも高い変化指標値分布が算出された小領域が変化領域であると認識することとなる。このようにしても、非揺らぎ成分による変化指標値を揺らぎ成分による変化指標値から精度良く分離して変化領域を高精度に認識できる。   Alternatively, when the conversion function is the above step function, a small region in which a change index value distribution having a higher time frequency when the luminance difference is 0 is higher than a time frequency when the luminance difference is not 0 is recognized as the changed region. Will be. Even in this case, the change index value based on the non-fluctuation component can be accurately separated from the change index value based on the fluctuation component, and the change region can be recognized with high accuracy.

ここで、図8(a),(b)は人が撮像されている分布画像とその分布画像における評価値Rの最頻値Rfの一例を示す。図8(a)に示す撮像画像のA−A線においては、図8(b)に示すように、人の胴体部分と両腕部分に該当する位置の最頻値Rfが1を示し、他の位置の最頻値Rfが0を示している。   Here, FIGS. 8A and 8B show an example of a distribution image in which a person is imaged and the mode value Rf of the evaluation value R in the distribution image. In the AA line of the captured image shown in FIG. 8A, as shown in FIG. 8B, the mode value Rf at the position corresponding to the human torso and both arms shows 1 and others. The mode value Rf at the position of 0 indicates 0.

なお、空間変化認識部24は、変換関数を介さずに変化指標値分布を直接解析して変化を認識してもよい。この場合、空間変化認識部24は、例えば小領域における変化指標値分布における最頻値が非0であれば当該小領域は変化領域と認識し、0であれば変化領域ではないと認識することができる。このようにしても、非揺らぎ成分による変化指標値を揺らぎ成分による変化指標値から精度良く分離して変化領域を高精度に認識できる。   The space change recognition unit 24 may recognize the change by directly analyzing the change index value distribution without using the conversion function. In this case, for example, the spatial change recognition unit 24 recognizes the small region as a change region if the mode value in the change index value distribution in the small region is non-zero, and recognizes that the small region is not a change region if it is 0. Can do. Even in this case, the change index value based on the non-fluctuation component can be accurately separated from the change index value based on the fluctuation component, and the change region can be recognized with high accuracy.

なお、変化領域と認識した小領域のうち近傍の小領域同士を連結することによって、これらをまとめた変化領域としてもよい。例えば、変化領域と認識した画素のうち8近傍の画素どうしに順次同一ラベルを付与し、ラベルごとの画素群を変化領域の情報とする。   In addition, it is good also as a change area | region which put these together by connecting adjacent small areas among the small areas recognized as a change area. For example, among the pixels recognized as the change region, the same label is sequentially given to the pixels in the vicinity of eight, and the pixel group for each label is used as the change region information.

空間変化認識部24は、上記のように認識した各変化領域の座標情報を侵入検知部25に出力する。   The space change recognition unit 24 outputs the coordinate information of each change region recognized as described above to the intrusion detection unit 25.

侵入検知部25は、変化領域の大きさが予め定めた範囲である場合に、監視空間内に侵入者の存在を検知したとして所定のアラーム信号を生成し、生成したアラーム信号を出力部5に出力する。   When the size of the change area is in a predetermined range, the intrusion detection unit 25 generates a predetermined alarm signal by detecting the presence of an intruder in the monitoring space, and sends the generated alarm signal to the output unit 5. Output.

なお、空間変化認識部24は、各変化領域の短時間特徴量分布、長時間特徴量分布、評価値を加えて出力してもよい。   The space change recognition unit 24 may add and output the short-time feature value distribution, long-time feature value distribution, and evaluation value of each change region.

その場合、侵入検知部25は、変化領域の大きさに加え、変化領域における短時間特徴量分布の特徴や変化領域における短時間特徴量分布の特徴と長時間特徴量分布の特徴の比較結果、変化領域の形状や、変化領域の動きを検知条件としてもよい。   In that case, in addition to the size of the change area, the intrusion detection unit 25 compares the characteristics of the short-time feature quantity distribution in the change area and the comparison results of the features of the short-time feature quantity distribution and the long-time feature quantity distribution in the change area, The shape of the change area or the movement of the change area may be used as the detection condition.

また、その場合、侵入検知部25は、大きさ、特徴の比較結果、形状、動きについて基準値との差を求め、評価値で重み付けた差がしきい値を超えた場合に侵入者が存在すると判定することもできる。   In that case, the intrusion detection unit 25 obtains a difference from the reference value for the size, feature comparison result, shape, and movement, and an intruder exists when the difference weighted by the evaluation value exceeds the threshold value. Then it can be determined.

次に、上記のように構成される空間認識装置1における撮像部2と処理部4の動作について図9のフローチャートを参照しながら説明する。   Next, operations of the imaging unit 2 and the processing unit 4 in the space recognition device 1 configured as described above will be described with reference to the flowchart of FIG.

撮像部2の各光電変換部20(20−1,20−2,…,20−N)は定期的に電子シャッターが作動するたびに、撮像素子の電圧をA/D変換して輝度値をサンプリングし、サンプリングされた各輝度値を撮像部2の特徴量分布生成部21がバッファに書き込むことで輝度信号をバッファリングする(S1)。   Each photoelectric conversion unit 20 (20-1, 20-2,..., 20-N) of the image pickup unit 2 performs A / D conversion on the voltage of the image pickup device every time the electronic shutter is operated to obtain a luminance value. Sampling is performed, and the luminance value sampled is written into the buffer by the feature amount distribution generation unit 21 of the imaging unit 2 to buffer the luminance signal (S1).

特徴量分布生成部21は、バッファへの書き込み回数あるいはバッファにおける最新の書き込み位置などを基に頻度算出周期が到来したか否かを判定する(S2)。   The feature amount distribution generation unit 21 determines whether or not the frequency calculation cycle has arrived based on the number of times of writing to the buffer or the latest writing position in the buffer (S2).

例えば、撮像周期を1/500秒、頻度算出周期を1/5秒とすれば、100個の輝度値を書き込むたびに頻度算出周期が到来する。   For example, if the imaging cycle is 1/500 seconds and the frequency calculation cycle is 1/5 second, the frequency calculation cycle comes every time 100 luminance values are written.

頻度算出周期が到来していなければ(S2にてNO)、特徴量分布生成部21は処理をステップS1へ戻し、バッファリングを続ける。   If the frequency calculation cycle has not arrived (NO in S2), the feature amount distribution generation unit 21 returns the process to step S1 and continues buffering.

他方、頻度算出周期が到来すると(S2にてYES)、特徴量分布生成部21は処理をステップS3へ進める。   On the other hand, when the frequency calculation cycle arrives (YES in S2), the feature amount distribution generation unit 21 advances the process to step S3.

特徴量分布生成部21は、短時間特徴量分布生成部210として動作し、バッファから最新の書き込み位置を終端とする第一時間幅の輝度信号を小領域ごとに読み出し、小領域ごとに輝度値の時間頻度を算出して短時間特徴量分布を生成する(S3)。   The feature amount distribution generation unit 21 operates as the short time feature amount distribution generation unit 210, reads a luminance signal having a first time width that ends at the latest writing position from the buffer for each small region, and sets a luminance value for each small region. Is calculated to generate a short-time feature amount distribution (S3).

例えば、撮像周期を1/500秒、第一時間幅を1/5秒間、輝度値を0〜255の整数、小領域を画素とすれば、画素ごとに、100個の輝度値から輝度値0,1,…,255の度数を計数して、計数結果のそれぞれを100で除した正規化ヒストグラムを短時間輝度値分布として算出する。   For example, if the imaging cycle is 1/500 seconds, the first time width is 1/5 seconds, the luminance value is an integer from 0 to 255, and the small area is a pixel, the luminance value is 0 from 100 luminance values for each pixel. , 1,..., 255 are counted, and a normalized histogram obtained by dividing each count result by 100 is calculated as a short-time luminance value distribution.

次に、特徴量分布生成部21は、長時間特徴量分布生成部211として動作し、バッファから最新の書き込み位置を終端とする第二時間幅の輝度信号を小領域ごとに読み出し、小領域ごとに輝度値の時間頻度を算出して長時間特徴量分布を生成する(S4)。   Next, the feature amount distribution generation unit 21 operates as the long-time feature amount distribution generation unit 211, reads out a luminance signal having a second time width that ends at the latest writing position from the buffer for each small region, The luminance frequency time frequency is calculated to generate a long-time feature amount distribution (S4).

例えば、撮像周期を1/500秒、第二時間幅を1分間、輝度値を0〜255の整数、小領域を画素とすれば、画素ごとに、30000個の輝度値から輝度値0,1,…,255の度数を計数して、計数結果のそれぞれを30000で除した正規化ヒストグラムを長時間輝度値分布として算出する。   For example, if the imaging cycle is 1/500 seconds, the second time width is 1 minute, the luminance value is an integer of 0 to 255, and the small area is a pixel, the luminance value is 0,1 from 30000 luminance values for each pixel. ,..., 255 are counted, and a normalized histogram obtained by dividing each count result by 30000 is calculated as a long-time luminance value distribution.

あるいは、長時間の輝度信号をバッファリングする代わりに短時間特徴量分布をバッファリングしておき、短時間特徴量分布を小領域ごとに積分および正規化して長時間特徴量分布を生成することもできる。   Alternatively, instead of buffering a long-time luminance signal, the short-time feature distribution can be buffered, and the short-time feature distribution can be integrated and normalized for each small region to generate a long-time feature distribution. it can.

撮像部2は、生成した短時間特徴量分布および長時間特徴量分布を処理部4に出力する。   The imaging unit 2 outputs the generated short-time feature value distribution and long-time feature value distribution to the processing unit 4.

なお、図9のフローチャートでは撮像部2と処理部4の処理を一連の処理として表現しているが、特徴量分布を出力した撮像部2は処理部4によるステップS3〜S9の処理終了を待たずに、ステップS1へ処理を戻してもよい。   In the flowchart of FIG. 9, the processes of the imaging unit 2 and the processing unit 4 are expressed as a series of processes. However, the imaging unit 2 that has output the feature amount distribution waits for the processing unit 4 to finish the processes in steps S3 to S9. Instead, the process may be returned to step S1.

続いて、処理部4の変化指標値分布生成部22は、入力された短時間特徴量分布および長時間特徴量分布の対応する小領域ごとに変化指標値分布を算出して指標値分布画像を生成する(S5)。   Subsequently, the change index value distribution generation unit 22 of the processing unit 4 calculates the change index value distribution for each of the input short-time feature quantity distribution and corresponding small regions of the long-time feature quantity distribution, and generates an index value distribution image. Generate (S5).

例えば、変化指標値分布生成部22は、短時間輝度値分布および長時間輝度値分布の対応する画素ごとに式(1)を適用して輝度差分布を算出する。   For example, the change index value distribution generation unit 22 calculates the luminance difference distribution by applying Equation (1) for each corresponding pixel of the short-time luminance value distribution and the long-time luminance value distribution.

続いて、処理部4の空間変化認識部24は、指標値/評価値変換関数記憶部23から変換関数を読み出し、各小領域の変化指標値分布に変換関数を適用して当該小領域の評価値分布を算出する(S6)。例えば、空間変化認識部24は、各画素の輝度差分布に式(2)を適用して当該画素の評価値分布を算出する。   Subsequently, the spatial change recognition unit 24 of the processing unit 4 reads the conversion function from the index value / evaluation value conversion function storage unit 23, applies the conversion function to the change index value distribution of each small region, and evaluates the small region. A value distribution is calculated (S6). For example, the spatial change recognition unit 24 calculates the evaluation value distribution of the pixel by applying Equation (2) to the luminance difference distribution of each pixel.

全ての小領域において評価値分布の最頻値が0(最も時間頻度の高い評価値が0)の場合、空間変化認識部24は、変化領域が検出されなかったとして(ステップS7にてNO)、処理をステップS1へ戻す。処理部4は特徴量分布の入力待ちとなる。   When the mode value of the evaluation value distribution is 0 in all the small regions (the evaluation value with the highest time frequency is 0), the spatial change recognition unit 24 assumes that no change region has been detected (NO in step S7). The process returns to step S1. The processing unit 4 waits for input of the feature amount distribution.

他方、評価値分布の最頻値が0より大きな小領域が存在している場合、空間変化認識部24は、変化領域が検出されたとして(ステップS7にてYES)、変化領域の情報を生成し、処理をステップS8に進める。   On the other hand, if there is a small region where the mode value of the evaluation value distribution is greater than 0, the space change recognition unit 24 generates information on the changed region, assuming that the changed region is detected (YES in step S7). Then, the process proceeds to step S8.

例えば、空間変化認識部24は、評価値分布の最頻値が0より大きな画素のうち座標が隣接している画素を一つの変化領域としてまとめ、各変化領域の短時間特徴量分布、長時間特徴量分布および変化領域らしさを侵入検知部25に出力する。   For example, the spatial change recognizing unit 24 puts together pixels whose coordinates are adjacent to each other among the pixels having the mode value of the evaluation value distribution greater than 0 as one change area. The feature amount distribution and the change area likelihood are output to the intrusion detection unit 25.

続いて、処理部4の侵入検知部25は、空間変化認識部24が生成した変化領域の情報を侵入検知条件と比較する(S8)。   Subsequently, the intrusion detection unit 25 of the processing unit 4 compares the information on the change area generated by the space change recognition unit 24 with the intrusion detection condition (S8).

侵入検知条件を満たす変化領域が生成された場合は(S8にてYES)、所定のアラーム信号を生成して出力部5に入力する(S9)。   If a change area that satisfies the intrusion detection condition is generated (YES in S8), a predetermined alarm signal is generated and input to output unit 5 (S9).

そして、出力部5は、侵入検知部25からアラーム信号が入力されると、そのアラーム信号を警備センターへ送信する。   When the alarm signal is input from the intrusion detection unit 25, the output unit 5 transmits the alarm signal to the security center.

なお、侵入検知条件を満たす変化領域が生成されなかった場合は(S8にてNO)、ステップS9の処理はスキップされる。   If a change area that satisfies the intrusion detection condition is not generated (NO in S8), the process of step S9 is skipped.

以上の処理を終えると処理はステップS1へ戻され、撮像部2は次の分布画像の生成を行い、処理部4は撮像部2からの分布画像の入力を待つ。   When the above processing is completed, the processing returns to step S1, the imaging unit 2 generates the next distribution image, and the processing unit 4 waits for the input of the distribution image from the imaging unit 2.

[変形例]
上述した実施の形態においては、特徴量どうしの相違性を表す変化指標値として輝度差、すなわち特徴量どうしの差を用いる例を示した。別の実施の形態においては、特徴量どうしの相違性を表す変化指標値として特徴量間の距離値を用いることもできる。この場合も、変化指標値分布生成部22は、距離値ごとに、特徴量分布間における特徴量の時間頻度の組み合わせのうち当該距離値が算出された組み合わせごとの時間頻度の積を総和して、当該距離値の時間頻度を算出することで距離値ごとの時間頻度を要素とする距離値分布を求める。また、指標値/評価値変換関数記憶部23は、大きな変化指標値ほど変化らしさを表す評価値を高く割り当てる変換関数を記憶する。そして、空間変化認識部24は、その変換関数を用いて距離値の時間頻度を当該距離値に対応する評価値の時間頻度に変換することによって各小領域の距離値分布を評価値の時間頻度を要素とする評価値分布に変換し、各小領域における評価値分布の最頻値に基づいて空間における変化を認識する。
[Modification]
In the above-described embodiment, an example in which a luminance difference, that is, a difference between feature amounts is used as a change index value representing a difference between feature amounts has been described. In another embodiment, a distance value between feature quantities can be used as a change index value representing a difference between feature quantities. Also in this case, the change index value distribution generation unit 22 sums, for each distance value, the product of the time frequencies for each combination for which the distance value is calculated among the combinations of the time frequencies of the feature amounts between the feature amount distributions. Then, by calculating the time frequency of the distance value, a distance value distribution having the time frequency for each distance value as an element is obtained. Further, the index value / evaluation value conversion function storage unit 23 stores a conversion function that assigns a higher evaluation value representing the degree of change to a larger change index value. Then, the space change recognition unit 24 uses the conversion function to convert the time frequency of the distance value into the time frequency of the evaluation value corresponding to the distance value, thereby converting the distance value distribution of each small region into the time frequency of the evaluation value. Is converted into an evaluation value distribution having elements as elements, and changes in the space are recognized based on the mode value of the evaluation value distribution in each small region.

また別の実施の形態においては、特徴量どうしの類似性を表す変化指標値を用いることができる。この場合、変化指標値は類似度または相関値などとすることができる。例えば、上述した輝度差を255から減じた値は類似度の例である。この場合も、変化指標値分布生成部22は、類似度(または相関値)ごとに、特徴量分布間における特徴量の時間頻度の組み合わせのうち当該類似度(または相関値)が算出された組み合わせごとの時間頻度の積を総和して、当該類似度(または相関値)の時間頻度を算出することで類似度(または相関値)ごとの時間頻度を要素とする類似度分布(または相関値分布)を求める。また、この場合、指標値/評価値変換関数記憶部23は小さな変化指標値ほど変化らしさを表す評価値を高く割り当てる変換関数を記憶する。そして、空間変化認識部24は、その変換関数を用いて類似度(または相関値)の時間頻度を当該類似度(または相関値)に対応する評価値の時間頻度に変換することによって各小領域の類似度分布(または相関値分布)を評価値の時間頻度を要素とする評価値分布に変換し、各小領域における評価値分布の最頻値に基づいて空間における変化を認識する。   In another embodiment, a change index value representing similarity between feature quantities can be used. In this case, the change index value can be a similarity or a correlation value. For example, a value obtained by subtracting the above-described luminance difference from 255 is an example of similarity. Also in this case, the change index value distribution generation unit 22 calculates, for each similarity (or correlation value), the similarity (or correlation value) in which the similarity (or correlation value) is calculated among the time frequency combinations of the feature amounts between the feature amount distributions. Similarity distribution (or correlation value distribution) with the time frequency of each similarity (or correlation value) as an element by calculating the time frequency of the similarity (or correlation value) by summing up the products of time frequency for each ) In this case, the index value / evaluation value conversion function storage unit 23 stores a conversion function that assigns a higher evaluation value that represents the degree of change as the change index value becomes smaller. Then, the space change recognition unit 24 uses the conversion function to convert the time frequency of the similarity (or correlation value) into the time frequency of the evaluation value corresponding to the similarity (or correlation value), thereby converting each small region. The similarity distribution (or correlation value distribution) is converted into an evaluation value distribution using the time frequency of the evaluation value as an element, and a change in the space is recognized based on the mode value of the evaluation value distribution in each small region.

また、上述した実施の形態および各変形例において、特徴量は輝度値、輝度値のヒストグラムまたは動きベクトルなどとすることができる。なお、特徴量を、動きベクトルとする場合、変化指標値分布生成部22は、分布画像間で近傍の小領域を含めた範囲で変化指標値である動きベクトルを算出することとなる。   In the above-described embodiment and each modification, the feature amount may be a luminance value, a luminance value histogram, a motion vector, or the like. When the feature quantity is a motion vector, the change index value distribution generation unit 22 calculates a motion vector that is a change index value in a range including a small area nearby between distribution images.

また、上述した実施の形態および各変形例において、図3の変化指標値分布生成部22は、長時間特徴量分布生成部211を設けずに、短時間特徴量分布生成部210にて前後して生成された短時間特徴量分布間で小領域ごとに変化指標値分布を算出してもよい。   Further, in the above-described embodiment and each modification, the change index value distribution generation unit 22 in FIG. 3 does not include the long-time feature amount distribution generation unit 211, and is changed by the short-time feature amount distribution generation unit 210. The change index value distribution may be calculated for each small region between the short-time feature amount distributions generated in this way.

また、上述した実施の形態および各変形例において、特徴量分布生成部21は、処理部4内で実現してもよい。その場合、撮像部2は、撮像周期ごとに輝度値を並べた撮像画像を処理部4に出力する。処理部4は、特徴量分布生成部21として動作し、撮像部2からの撮像画像を記憶部3にバッファリングして頻度算出周期ごとに特徴量分布を生成する。   In the above-described embodiment and each modification, the feature amount distribution generation unit 21 may be realized in the processing unit 4. In that case, the imaging unit 2 outputs a captured image in which luminance values are arranged for each imaging cycle to the processing unit 4. The processing unit 4 operates as the feature amount distribution generation unit 21, buffers the captured image from the imaging unit 2 in the storage unit 3, and generates a feature amount distribution for each frequency calculation period.

上述した実施の形態および各変形例では、空間変化認識部24は、変化領域を認識したが、これに限らず、空間における明るさの変化、物体の動き(速さまたは/および移動方向)、物体数の増減など種々の変化の認識に用いることができる。   In the above-described embodiment and each modified example, the space change recognition unit 24 recognizes the change region, but is not limited to this, the brightness change in the space, the movement of the object (speed or / and moving direction), It can be used for recognition of various changes such as increase / decrease in the number of objects.

上述した実施の形態および各変形例では、計測部が監視空間の可視光線の強度に応じた振幅を有する輝度信号を出力する構成としているが、これに限定されるものではない。   In the above-described embodiment and each modification, the measurement unit outputs a luminance signal having an amplitude corresponding to the intensity of visible light in the monitoring space, but is not limited thereto.

例えば、計測部が監視空間の近赤外線の強度に応じた振幅を有する輝度信号を出力する構成、計測部が監視空間の赤外線の強度に応じた振幅を有する輝度信号を出力する構成、計測部が監視空間の紫外線の強度に応じた振幅を有する輝度信号を出力する構成、計測部が監視空間のX線の強度に応じた振幅を有する輝度信号を出力する構成、計測部としての距離センサが当該計測部から監視空間に存在する物体までの距離に応じた振幅を有する距離信号を出力する構成など、計測部が空間の所定物理量を計測して計測信号を出力する構成であってもよい。なお、上記距離センサは、例えば、レーザー光を投受光する方式、マイクロ波またはミリ波などの電磁波を送受信する方式、超音波を送受信する方式などが考えられる。   For example, a configuration in which the measurement unit outputs a luminance signal having an amplitude corresponding to the intensity of near infrared light in the monitoring space, a configuration in which the measurement unit outputs a luminance signal having an amplitude according to the intensity of infrared light in the monitoring space, and the measurement unit A configuration for outputting a luminance signal having an amplitude corresponding to the intensity of ultraviolet light in the monitoring space, a configuration for the measurement unit to output a luminance signal having an amplitude corresponding to the intensity of the X-ray in the monitoring space, and a distance sensor as the measurement unit A configuration in which the measurement unit measures a predetermined physical quantity in the space and outputs a measurement signal, such as a configuration that outputs a distance signal having an amplitude corresponding to the distance from the measurement unit to an object existing in the monitoring space, may be used. As the distance sensor, for example, a method of projecting and receiving laser light, a method of transmitting and receiving electromagnetic waves such as microwaves and millimeter waves, and a method of transmitting and receiving ultrasonic waves are conceivable.

上述した実施の形態および各変形例では、複数の小領域を計測面の全体に設定する例を示したが、計測面の一部に複数の小領域を設定してもよいし、計測面の一部に1つの小領域を設定してもよい。   In the above-described embodiment and each modification, an example in which a plurality of small areas are set on the entire measurement surface has been described. However, a plurality of small areas may be set on a part of the measurement surface. One small area may be set for a part.

以上、本発明に係る空間認識装置の最良の形態について説明したが、この形態による記述及び図面により本発明が限定されることはない。すなわち、この形態に基づいて当業者等によりなされる他の形態、実施例及び運用技術などはすべて本発明の範疇に含まれることは勿論である。   The best mode of the space recognition device according to the present invention has been described above, but the present invention is not limited by the description and drawings according to this mode. That is, it is a matter of course that all other forms, examples, operation techniques, and the like made by those skilled in the art based on this form are included in the scope of the present invention.

1 空間認識装置
2 撮像部
3 記憶部
4 出力部
20(20−1,20−2,…,20−N) 光電変換部
21 特徴量分布生成部
22 変化指標値分布生成部
23 指標値/評価値変換関数記憶部
24 空間変化認識部
25 侵入検知部
210 短時間特徴量分布生成部
211 長時間特徴量分布生成部
DESCRIPTION OF SYMBOLS 1 Spatial recognition apparatus 2 Imaging part 3 Memory | storage part 4 Output part 20 (20-1, 20-2, ..., 20-N) Photoelectric conversion part 21 Feature-value distribution generation part 22 Change index value distribution generation part 23 Index value / evaluation Value conversion function storage unit 24 Spatial change recognition unit 25 Intrusion detection unit 210 Short time feature amount distribution generation unit 211 Long time feature amount distribution generation unit

Claims (5)

所定の計測面上に配置され、所定の空間の物理量を計測して計測信号を出力する計測部と、
前記計測面上に設定された小領域において、互いに異なる時間区間の前記計測信号のそれぞれから特徴量ごとの時間頻度を算出して前記時間区間ごとの特徴量分布を生成する特徴量分布生成部と、
前記時間区間が異なる前記特徴量分布間で、前記小領域における前記特徴量どうしの相違性または類似性を表す変化指標値を求め、前記変化指標値ごとの時間頻度を要素とする変化指標値分布を前記小領域ごとに生成する変化指標値分布生成部と、
前記小領域における前記変化指標値分布を用いて前記空間における変化を認識する空間変化認識部と、
を備えたことを特徴とする空間認識装置。
A measurement unit arranged on a predetermined measurement surface and measuring a physical quantity of a predetermined space and outputting a measurement signal;
A feature amount distribution generation unit that calculates a time frequency for each feature amount from each of the measurement signals in different time intervals in a small region set on the measurement surface and generates a feature amount distribution for each time interval ; ,
A change index value distribution that obtains a change index value representing a difference or similarity between the feature quantities in the small region between the feature quantity distributions having different time intervals, and uses a time frequency for each change index value as an element A change index value distribution generation unit for generating each small region ,
A spatial change recognition unit that recognizes a change in the space using the change index value distribution in the small region;
A space recognition device comprising:
前記変化指標値分布生成部は、前記変化指標値ごとに、前記時間区間が異なる前記特徴量分布間における前記特徴量の前記時間頻度の組み合わせのうち当該変化指標値が算出された組み合わせごとの時間頻度の積を総和して、当該変化指標値の時間頻度を算出する請求項1に記載の空間認識装置。 The change index value distribution generation unit is configured to calculate, for each change index value, a time for each combination in which the change index value is calculated among the time frequency combinations of the feature quantities between the feature quantity distributions having different time intervals. The space recognition apparatus according to claim 1, wherein the frequency products are summed to calculate the time frequency of the change index value. 前記空間変化認識部は、前記変化指標値を変化らしさの度合いを表す評価値に変換する変換関数を用いて、前記小領域の前記変化指標値分布を前記各変化指標値に対応する前記評価値ごとの時間頻度を要素とする評価値分布に変換し、前記小領域における前記評価値分布を用いて前記空間における変化を認識する請求項1または2に記載の空間認識装置。 The spatial change recognition unit uses the conversion function for converting the change index value into an evaluation value representing a degree of change, and the evaluation value corresponding to each change index value in the change index value distribution of the small region The space recognition device according to claim 1, wherein a change in the space is recognized using the evaluation value distribution in the small area by converting the evaluation value distribution having each time frequency as an element. 前記特徴量分布生成部は、
第一時間区間の前記計測信号から前記特徴量ごとの前記時間頻度を算出して短時間特徴量分布を生成する短時間特徴量分布生成部と、
前記第一時間区間よりも長い第二時間区間の前記計測信号から前記特徴量ごとの前記時間頻度を算出して長時間特徴量分布を生成する長時間特徴量分布生成部と、
を備え、
前記変化指標値分布生成部は、前記短時間特徴量分布と前記長時間特徴量分布の間で前記小領域における前記変化指標値分布を生成する請求項1〜3の何れか一つに記載の空間認識装置。
The feature quantity distribution generation unit
A short-time feature value distribution generation unit that calculates the time frequency for each feature value from the measurement signal in a first time interval and generates a short-time feature value distribution;
A long-time feature amount distribution generation unit that calculates the time frequency for each feature amount from the measurement signal in a second time interval longer than the first time interval and generates a long-time feature amount distribution;
With
The said change index value distribution production | generation part produces | generates the said change index value distribution in the said small area | region between the said short time feature-value distribution and the said long-time feature-value distribution. Spatial recognition device.
前記計測部は、前記計測面を撮像面として該撮像面上に配置され、前記空間を撮像して前記計測信号としての輝度信号を出力する光電変換部からなることを特徴とする請求項1〜4の何れか一つに記載の空間認識装置。 The said measurement part is arrange | positioned on this imaging surface by making the said measurement surface into an imaging surface, and consists of a photoelectric conversion part which images the said space and outputs the luminance signal as said measurement signal. The space recognition device according to any one of 4.
JP2016016075A 2016-01-29 2016-01-29 Space recognition device Active JP6412032B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016016075A JP6412032B2 (en) 2016-01-29 2016-01-29 Space recognition device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016016075A JP6412032B2 (en) 2016-01-29 2016-01-29 Space recognition device

Publications (2)

Publication Number Publication Date
JP2017135656A JP2017135656A (en) 2017-08-03
JP6412032B2 true JP6412032B2 (en) 2018-10-24

Family

ID=59503790

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016016075A Active JP6412032B2 (en) 2016-01-29 2016-01-29 Space recognition device

Country Status (1)

Country Link
JP (1) JP6412032B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11195025B2 (en) 2017-11-27 2021-12-07 Nec Corporation Information processing device, information processing method, and storage medium for temporally dividing time-series data for analysis
JP2020149158A (en) * 2019-03-11 2020-09-17 Necエンベデッドプロダクツ株式会社 Crime prevention device, crime prevention method, and program

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3123587B2 (en) * 1994-03-09 2001-01-15 日本電信電話株式会社 Moving object region extraction method using background subtraction
JP5218906B2 (en) * 2008-10-16 2013-06-26 国立大学法人東京農工大学 Smoke detection device and smoke detection method
JP5388829B2 (en) * 2009-12-11 2014-01-15 セコム株式会社 Intruder detection device
JP5715863B2 (en) * 2011-03-25 2015-05-13 セコム株式会社 Image processing device
JP2013070129A (en) * 2011-09-21 2013-04-18 Panasonic Corp Image information extraction apparatus, image transmission apparatus using the same, image receiving apparatus, and image transmission system
JP6020199B2 (en) * 2013-01-24 2016-11-02 株式会社ソシオネクスト Image processing apparatus, method, program, and imaging apparatus

Also Published As

Publication number Publication date
JP2017135656A (en) 2017-08-03

Similar Documents

Publication Publication Date Title
JP4429298B2 (en) Object number detection device and object number detection method
JP4389602B2 (en) Object detection apparatus, object detection method, and program
US20100046802A1 (en) Distance estimation apparatus, distance estimation method, storage medium storing program, integrated circuit, and camera
GB2475104A (en) Detecting movement of 3D objects using a TOF camera
US8923552B2 (en) Object detection apparatus and object detection method
JP6412032B2 (en) Space recognition device
JP4079136B2 (en) Motion detection device and motion detection method
JP4887540B2 (en) Vehicle periphery monitoring device, vehicle, vehicle periphery monitoring program, and vehicle periphery monitoring method
JP3409538B2 (en) Human body detection device
JP4918615B2 (en) Object number detection device and object number detection method
JP6214425B2 (en) Specific motion detection device
JP5679760B2 (en) Intruder detection device
KR101733668B1 (en) Object detection apparatus, surveillance camera, and object detection method
JP4580189B2 (en) Sensing device
JPH03242592A (en) Monitoring device for image
JP6332703B2 (en) Space recognition device
JP2693586B2 (en) Image identification / tracking device
JP6332704B2 (en) Object detection device
JP4674920B2 (en) Object number detection device and object number detection method
JPWO2018168182A1 (en) MOBILE BODY DETECTING DEVICE, MOBILE BODY DETECTING METHOD, AND MOBILE BODY DETECTING PROGRAM
JP6332702B2 (en) Space recognition device
JP2018166308A (en) Camera calibration device
JP4645321B2 (en) Moving object detection device using images
JP4509654B2 (en) Sensing device
KR100744158B1 (en) Probabilistic visual surveillance system and method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20170516

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20180213

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180320

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180502

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20180927

R150 Certificate of patent or registration of utility model

Ref document number: 6412032

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250