JP6763241B2 - Person detection method and person detection device - Google Patents

Person detection method and person detection device Download PDF

Info

Publication number
JP6763241B2
JP6763241B2 JP2016170268A JP2016170268A JP6763241B2 JP 6763241 B2 JP6763241 B2 JP 6763241B2 JP 2016170268 A JP2016170268 A JP 2016170268A JP 2016170268 A JP2016170268 A JP 2016170268A JP 6763241 B2 JP6763241 B2 JP 6763241B2
Authority
JP
Japan
Prior art keywords
event
signal
sensor
person detection
person
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
JP2016170268A
Other languages
Japanese (ja)
Other versions
JP2018033746A (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.)
Oki Electric Industry Co Ltd
Original Assignee
Oki Electric Industry 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 Oki Electric Industry Co Ltd filed Critical Oki Electric Industry Co Ltd
Priority to JP2016170268A priority Critical patent/JP6763241B2/en
Publication of JP2018033746A publication Critical patent/JP2018033746A/en
Application granted granted Critical
Publication of JP6763241B2 publication Critical patent/JP6763241B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

この発明は、人検知方法及び人検知装置に関し、特に、検知エリア内の一人または複数の人の存在、歩行、作業などの行動を一人一人分離して検知する人検知方法及び人検知装置に関する。 The present invention relates to a person detection method and a person detection device, and more particularly to a person detection method and a person detection device that separately detects the presence of one or more people in a detection area, walking, work, and other actions.

IoT(モノのインターネット)社会への期待が増すとともに、人の動き、行動などを自動で検知できる人検知技術へのニーズが高まっている。高度な人検知技術がIT技術と結びつけば、例えば子供や要介護者の行動をセンサで自動収集して外出中スマートフォンで確認する、公共の場所で混雑状況や人の流れを把握する、あるいは不審者を発見する、店舗内で顧客がどのように移動しどこに滞留しているかの情報を得てマーケティングなどに役立てることなどが可能になる。このため、家庭、オフィス、医療、店舗など多くの場所で応用が期待されている。 With increasing expectations for the IoT (Internet of Things) society, there is an increasing need for human detection technology that can automatically detect the movements and behaviors of people. If advanced human detection technology is combined with IT technology, for example, the behavior of children and care recipients can be automatically collected by sensors and checked on a smartphone while on the go, the congestion situation and the flow of people can be grasped in public places, or suspicious. It will be possible to discover people, obtain information on how customers move and stay in the store, and use it for marketing. Therefore, it is expected to be applied in many places such as homes, offices, medical care, and stores.

人検知技術として、画像による人検知、赤外線による人検知、電磁波及び超音波照射による人検知、静電気による人検知などがある。 Human detection technology includes human detection by images, human detection by infrared rays, human detection by electromagnetic wave and ultrasonic irradiation, and human detection by static electricity.

画像による人検知では、例えば、監視カメラで得られる画像情報をソフトウエアにより画像解析する。画像による人検知には、画像情報に、必然的に、死角(ブラインドエリア)が存在するという課題がある。また、画像処理のためのシステムが高コストとなるという課題がある。さらに、画像情報の使用については、プライバシーの問題もある。 In human detection using images, for example, image information obtained by a surveillance camera is image-analyzed by software. Human detection by image has a problem that a blind spot (blind area) inevitably exists in the image information. Further, there is a problem that the system for image processing becomes expensive. Furthermore, there is also a privacy issue regarding the use of image information.

赤外線照射による人検知は、人体から放射される赤外線を赤外線センサで観測するものである。この赤外線照射による人検知は、システムが比較的低コストであり、また、プライバシーの問題も生じない。しかしながらこの方法では、検知範囲内に人が入ったかどうかはわかるものの、検知範囲内のどの位置に入ったか、何人入ったか、どう動いたかなど、より詳細な情報はわからない。 Human detection by infrared irradiation is to observe infrared rays emitted from the human body with an infrared sensor. This infrared irradiation human detection system is relatively low cost and does not cause privacy problems. However, with this method, although it is possible to know whether or not a person has entered the detection range, it is not possible to obtain more detailed information such as which position within the detection range, how many people have entered, and how the person has moved.

電磁波及び超音波照射による人検知は、近赤外線やマイクロ波、あるいは超音波など、人体にとって無害な信号を検知エリアに照射し、人体からのそれらの反射を測定するものである。この電磁波及び超音波照射による人検知は、画像による人検知に比べてデータ量が小さく解析が容易である。また、細かい情報が得られないので、プライバシーの問題が少ない。しかし、死角の問題は、監視カメラの場合より大きくなる。 Human detection by electromagnetic wave and ultrasonic irradiation irradiates the detection area with signals that are harmless to the human body, such as near infrared rays, microwaves, and ultrasonic waves, and measures their reflection from the human body. This human detection by electromagnetic wave and ultrasonic irradiation has a smaller amount of data and is easier to analyze than human detection by images. In addition, since detailed information cannot be obtained, there are few privacy issues. However, the problem of blind spots is greater than in the case of surveillance cameras.

静電気による人検知は、人体の動きがもたらす静電界の変化を測定するものである。静電界を用いた移動の検出では、画像解析を中心とした他の方法と比べて、視覚的には死角になる場所の情報も得ることができる利点がある。 Human detection by static electricity measures changes in the electrostatic field caused by the movement of the human body. The detection of movement using an electrostatic field has an advantage that it is possible to visually obtain information on a blind spot as compared with other methods centered on image analysis.

静電気による人検知について説明する。人体は導電性を持ち良好な導電体と近似される。また通常人体は地面など電気的アースと見なせるものから絶縁されている。このため、何らかの理由で帯電することによって、人体は様々な電位を持ち得る。様々な電位を持った導電体としての人体が空間中を移動するとき、空間中の静電界分布が変動する。この結果、電界の変化などを計測することにより人体の動きを知ることができる。この一例として、人の歩行に伴う電界の変化を観測することによって、個人の識別が可能となる装置が開示されている(例えば、特許文献1参照)。 A person detection by static electricity will be described. The human body has conductivity and is approximated to a good conductor. In addition, the human body is usually insulated from what can be regarded as an electrical ground such as the ground. Therefore, the human body can have various electric potentials by being charged for some reason. When the human body as a conductor with various potentials moves in space, the electrostatic field distribution in space fluctuates. As a result, the movement of the human body can be known by measuring changes in the electric field and the like. As an example of this, a device that enables individual identification by observing a change in an electric field accompanying a person's walking is disclosed (see, for example, Patent Document 1).

人体に限らず電荷を持つものにより生じる静電界の変化は、準静的すなわち変化の速度が十分小さい場合、以下の式(1)で示される、ポアソン方程式によって決定される。 The change in the electrostatic field caused by not only the human body but also an electric charge is determined by the Poisson's equation represented by the following equation (1) when it is quasi-static, that is, when the rate of change is sufficiently small.

Figure 0006763241
Figure 0006763241

ここで、Vは電位、ρは電荷密度、εは誘電率で全て空間上の位置(x,y,z)の関数である。上式(1)は空間上でρやεの空間分布が変化することによって空間上の電位の分布が変化することを示す。人体は電荷を持ち、かつ空気とは異なる誘電率を持つので、人体が移動することにより電荷や誘電率の空間分布が変化し、それによって電位分布が変化する。そこで人体が運動する空間上に複数の電極(以下センサとも称する)を置き、その電極で生じる電位を観察すれば人体の動きに関する何らかの情報が得られることがわかる。 Here, V is an electric potential, ρ is a charge density, and ε is a permittivity, which are all functions of positions (x, y, z) in space. Equation (1) above shows that the distribution of electric potential in space changes as the spatial distribution of ρ and ε changes in space. Since the human body has an electric charge and a dielectric constant different from that of air, the spatial distribution of the electric charge and the dielectric constant changes as the human body moves, and the potential distribution changes accordingly. Therefore, it can be seen that some information on the movement of the human body can be obtained by placing a plurality of electrodes (hereinafter also referred to as sensors) in the space where the human body moves and observing the potential generated by the electrodes.

なお、限られた箇所の電位を測定するだけでρやεの空間分布を決定することはできない。しかし、例えば人体を直径1mの導体球と仮定し、さらに検知範囲内には一人のみ存在し、他に移動するものはない、など単純化してモデル化すれば、少数の電極で人体の動きを推定できる。 It should be noted that the spatial distribution of ρ and ε cannot be determined only by measuring the potential at a limited location. However, if we assume that the human body is a conductor sphere with a diameter of 1 m, and if we model it simply, for example, there is only one person in the detection range and nothing else moves, the movement of the human body can be performed with a small number of electrodes. Can be estimated.

特開2004−147793号公報Japanese Unexamined Patent Publication No. 2004-147793

しかしながら、上述の従来例の静電気による人検知には、以下の課題がある。 However, the above-mentioned conventional example of detecting a person by static electricity has the following problems.

第一の問題は、複数の人体が移動している場合、それらの動きを分離してその移動を推定することが難しいことである。上式(1)は空間全ての点で成立するので、複数の人体が検知範囲内に散らばって存在する場合、センサで測定される電位分布はそれらの重ね合わせとなる。ところが測定された電位の変化が空間上のどの箇所の、人体の移動によるρとεの変化によるものかを求めるのは難しい。 The first problem is that when multiple human bodies are moving, it is difficult to separate their movements and estimate their movements. Since the above equation (1) holds at all points in the space, when a plurality of human bodies are scattered within the detection range, the potential distribution measured by the sensor is a superposition of them. However, it is difficult to determine where in space the measured potential change is due to changes in ρ and ε due to the movement of the human body.

この問題の解決策として非常に多くのセンサを設置する方法がある。しかし、センサの設置が困難となり、また多くのセンサの信号から一人一人の移動のもたらす効果を分離する複雑なアルゴリズムが必要となり、推定の誤りが増えることが予想される。 The solution to this problem is to install a large number of sensors. However, it becomes difficult to install sensors, and a complicated algorithm that separates the effects of each person's movement from the signals of many sensors is required, and it is expected that estimation errors will increase.

第二の問題は、センサがごく近傍の人の動きしか検知できないこと、および検出する運動の向きによって指向性があることである。その理由を図1を参照して説明する。図1は、静電気による人検知を説明するための模式図である。図1ではセンサとして4個の電界センサS1, S2, S3, S4を用いている。 The second problem is that the sensor can only detect the movement of a person in the immediate vicinity, and that it is directional depending on the direction of the detected movement. The reason will be described with reference to FIG. FIG. 1 is a schematic diagram for explaining human detection by static electricity. In FIG. 1, four electric field sensors S1, S2, S3, and S4 are used as sensors.

センサS1は極座標(r,θ,φ)の原点Oに、センサS2及びセンサS3はそれぞれφ=0の平面上の、O´及びO´´の位置にそれぞれ置かれているものと仮定する。測定対象の人体は点P上にあり、一定の電荷量Qを持つ空間上の一点と考える。人体(点P)は、φ=0の平面を移動するものとする。φ=0の平面は、例えば床面に対応する。点Pの位置を(L,0,0)とすると、センサS1と点Pの距離OPはLである。さらにセンサS2と点Pの距離O´PをL、センサS3と点Pの距離O´´PをLとし、OPとO´Pのなす角をΘ´、OPとO´´Pのなす角をΘ´´とする。このとき、センサS1〜S3で測定される電界の絶対値はそれぞれ以下の式(2)で表される。 It is assumed that the sensor S1 is located at the origin O in polar coordinates (r, θ, φ), and the sensor S2 and the sensor S3 are located at positions O ′ and O ″ on a plane of φ = 0, respectively. The human body to be measured is on the point P, and is considered to be a point on the space having a constant charge amount Q. It is assumed that the human body (point P) moves in a plane of φ = 0. The plane of φ = 0 corresponds to, for example, the floor surface. Assuming that the position of the point P is (L 1 , 0, 0), the distance OP between the sensor S1 and the point P is L 1 . Further, the distance O'P between the sensor S2 and the point P is L 2 , the distance O'P between the sensor S3 and the point P is L 3 , the angle formed by OP and O'P is Θ', and OP and O'P. Let the angle to be formed be Θ´´. At this time, the absolute value of the electric field measured by the sensors S1 to S3 is represented by the following equation (2), respectively.

Figure 0006763241
Figure 0006763241

ここで、εは真空の誘電率である。電荷を持った人体の移動、例えば点PがセンサS1に近づいたり遠ざかったりするときそれぞれのセンサが感じる電界強度が変化する。それぞれのセンサの感じるLに対する電界強度の変化率は上式(2)をLで微分することによって得られ、以下の式(3)で表される。 Here, ε 0 is the permittivity of the vacuum. When the charged human body moves, for example, the point P approaches or moves away from the sensor S1, the electric field strength felt by each sensor changes. The rate of change in the electric field strength with respect to L 1 felt by each sensor is obtained by differentiating the above equation (2) with L 1 , and is expressed by the following equation (3).

Figure 0006763241
Figure 0006763241

上式(3)で表される電界E〜Eの変動はそれぞれのセンサによって測定される。人体の電荷Qは一定であるという仮定の下では、Θ´、Θ´´はL〜Lが決まれば一意に決まる。従って、上式(3)の変数はL〜Lの3つである。上式(3)においては式の数と変数が等しいからL〜Lは決定可能である。つまり点Pの位置を特定できる。静電気による人検知は基本的にこのような原理に基づく。 The fluctuations of the electric fields E 1 to E 3 represented by the above equation (3) are measured by each sensor. Under the assumption that the electric charge Q of the human body is constant, Θ ′ and Θ ″ are uniquely determined if L 1 to L 3 are determined. Therefore, there are three variables in the above equation (3), L 1 to L 3 . In the above equation (3), L 1 to L 3 can be determined because the number of equations and the variable are equal. That is, the position of the point P can be specified. Human detection by static electricity is basically based on this principle.

しかしながら上式(3)を見ると、電界強度の変化はセンサからの距離の3乗に反比例している。これはセンサと検知対象の距離が大きくなるにつれ人体の移動の検知が急激に困難になることを示すものである。例えば図1でOP間の距離Lに比べて、O´P間の距離Lが10倍大きいと仮定する。するとセンサS2で測定される電界強度の変化は以下の式(4)となる。 However, looking at the above equation (3), the change in electric field strength is inversely proportional to the cube of the distance from the sensor. This indicates that the detection of the movement of the human body becomes rapidly difficult as the distance between the sensor and the detection target increases. For example, in FIG. 1, it is assumed that the distance L 2 between O'P is 10 times larger than the distance L 1 between OPs. Then, the change in the electric field strength measured by the sensor S2 is given by the following equation (4).

Figure 0006763241
Figure 0006763241

このように、センサS2で得られる値は、センサS1で得られる値の1/1000より小さくなり、その検知が困難となってしまう。上式(3)において、いずれかのセンサの変化が観測できないと、点Pの運動を決定するための方程式の数が減り、結果として人体の運動の検知ができなくなる。このように従来の電界センサによる人検知においては、センサがごく近傍の人体の動きしか検知できない問題がある。 As described above, the value obtained by the sensor S2 is smaller than 1/1000 of the value obtained by the sensor S1, and its detection becomes difficult. In the above equation (3), if the change of any of the sensors cannot be observed, the number of equations for determining the motion of the point P decreases, and as a result, the motion of the human body cannot be detected. As described above, in the conventional human detection by the electric field sensor, there is a problem that the sensor can detect only the movement of the human body in the immediate vicinity.

次に電界センサの指向性の問題について説明する。点Pの運動について、例え点PがセンサS1の近傍にいても、OPと直交するθ,φ方向のPの移動においては、センサS1 で電界の変化は生じない。このため、センサS1から点Pの動きに関する情報が得られないという問題がある。 Next, the problem of directivity of the electric field sensor will be described. Regarding the movement of the point P, even if the point P is in the vicinity of the sensor S1, the electric field does not change in the sensor S1 in the movement of P in the θ and φ directions orthogonal to the OP. Therefore, there is a problem that information on the movement of the point P cannot be obtained from the sensor S1.

この問題を解決するには、センサS2、S3と同じくOP軸から外れた位置にあるセンサS4において電界の変化を観測し、3つのセンサS2 ,S3, S4から点Pの運動を決定する必要がある。ところがセンサS4は、センサS2,S3よりも点Pからさらに離れた場所にあり、電界の変化を観測できない可能性が高い。それを避けるためには各々のセンサに非常に感度の高いものを用いるか、センサ同士の間隔を十分な観測が可能となるまで狭める必要がある。 In order to solve this problem, it is necessary to observe the change in the electric field with the sensor S4 located off the OP axis like the sensors S2 and S3, and determine the motion of the point P from the three sensors S2, S3 and S4. is there. However, the sensor S4 is located further away from the point P than the sensors S2 and S3, and there is a high possibility that the change in the electric field cannot be observed. In order to avoid this, it is necessary to use highly sensitive sensors for each sensor or to narrow the distance between the sensors until sufficient observation is possible.

以上述べたように、人体の移動により生じる静電界の変化を観測することによる人検知において、歩行など人体の大きな動きを捉えるためには、非常に感度の高いセンサを用い、かつセンサの数を高密度にしなければならない。そのため、静電気による人検知システムは、ごく近接する人体のみを検知するタッチセンサなど限られた用途でしか用いられていなかった。 As described above, in human detection by observing changes in the electrostatic field caused by the movement of the human body, in order to capture large movements of the human body such as walking, a very sensitive sensor is used and the number of sensors is increased. Must be dense. Therefore, the static electricity detection system has been used only in limited applications such as a touch sensor that detects only a human body in close proximity.

この発明は、上述の問題点に鑑みてなされたものであり、この発明の目的は、複数の人がいる場合においてもそれらの移動を個人ごとに識別して可視化できる、静電気による、人検知方法及び人検知システムを提供することにある。 The present invention has been made in view of the above-mentioned problems, and an object of the present invention is a method for detecting a person by static electricity, which can identify and visualize the movement of a plurality of people for each individual. And to provide a person detection system.

上述した目的を達成するために、この発明の人検知方法は、以下の過程を備える。先ず、複数のセンサから、信号強度を示す強度信号列を取得し、0以外の信号強度を示す強度信号列間で相関係数を計算し、相関係数を所定の閾値と比較する。比較の結果、閾値以上である場合には、信号強度比からイベント位置を特定する。閾値よりも小さい場合には、複数イベント判定を行う。複数イベント判定を行う過程は、さらに、以下の過程を備える。先ず、イベント数ξを2として、ξ個のイベント信号列Iと、ξ行M列の係数行列Aと、M行1列の強度信号列Sで与えられる、I=ASから、イベント信号列I間の相関係数r(ξ)が最小となる初期係数行列Aを決定する。次に、イベント数ξに1を加算して、ξ個のイベント信号列Iと、ξ行M列の係数行列Aと、M行1列の強度信号列Sで与えられる、I=ASから、イベント信号列I間の相関係数r(ξ)が最小となる係数行列Aを決定し、r(ξ)とr(ξ−1)を比較する。r(ξ)がr(ξ−1)以下の場合は、さらに、係数行列Aを決定する過程及びr(ξ)とr(ξ−1)を比較する過程を行い、それ以外の場合は、ξ−1以下であるときに得られた係数行列Aの各行の係数からイベント位置を決定する。 In order to achieve the above-mentioned object, the person detection method of the present invention includes the following process. First, intensity signal sequences indicating signal strength are acquired from a plurality of sensors, a correlation coefficient is calculated between intensity signal sequences indicating signal intensities other than 0, and the correlation coefficient is compared with a predetermined threshold value. As a result of comparison, if it is equal to or higher than the threshold value, the event position is specified from the signal strength ratio. If it is smaller than the threshold value, multiple event determination is performed. The process of determining a plurality of events further includes the following process. First, assuming that the number of events ξ is 2, ξ event signal sequences I, a coefficient matrix A of ξ rows and M columns, and an intensity signal sequence S of M rows and 1 column are given from I = AS, and the event signal sequence I The initial coefficient matrix A that minimizes the correlation coefficient r (ξ) between them is determined. Next, by adding 1 to the number of events ξ, from I = AS given by ξ event signal strings I, a coefficient matrix A of ξ rows and M columns, and an intensity signal sequence S of M rows and 1 column, The coefficient matrix A that minimizes the correlation coefficient r (ξ) between the event signal sequences I is determined, and r (ξ) and r (ξ-1) are compared. When r (ξ) is less than or equal to r (ξ-1), the process of determining the coefficient matrix A and the process of comparing r (ξ) and r (ξ-1) are further performed. In other cases, the process is performed. The event position is determined from the coefficients of each row of the coefficient matrix A obtained when the value is ξ-1 or less.

また、この発明の人検知装置は、複数の静電気測定器と、複数の静電気測定器の出力の高周波成分を透過させるバンドパスフィルタと、バンドパスフィルタを透過した信号から強度信号列を生成するAD変換器と、信号処理部とを備える。信号処理部は、上述の人検知方法の各手順を実行する。 Further, the human detection device of the present invention includes a plurality of electrostatic measuring instruments, a bandpass filter that transmits a high frequency component of the output of the plurality of electrostatic measuring instruments, and an AD that generates an intensity signal sequence from a signal transmitted through the bandpass filter. It includes a converter and a signal processing unit. The signal processing unit executes each procedure of the above-mentioned person detection method.

この発明の人検知方法及び人検知装置によれば、複数の電界センサより生じる信号から信号間の相関する成分を抽出することによって、人体の歩行、他物体への接触などに生じる電界の変動イベントを抽出し、イベントが生じた位置を推定して、仮想人にイベントを割り当てることによって所定の時刻の人体の位置を推定できる。また、複数のイベントをイベントごとに分離できるので、複数の人体が存在する場合においても、それらの移動を個人ごとに識別できる。 According to the human detection method and the human detection device of the present invention, an electric field fluctuation event that occurs in walking of a human body, contact with another object, etc. by extracting correlated components between signals from signals generated by a plurality of electric field sensors. Can be extracted, the position where the event occurred, and the position of the human body at a predetermined time can be estimated by assigning the event to the virtual person. Further, since a plurality of events can be separated for each event, even when a plurality of human bodies exist, their movements can be identified for each individual.

静電気による人検知を説明するための模式図である。It is a schematic diagram for demonstrating the detection of a person by static electricity. この発明の人検知装置の第1実施形態の概略構成図である。It is a schematic block diagram of the 1st Embodiment of the person detection device of this invention. 信号処理部での処理を説明するための処理フロー図である。It is a processing flow diagram for demonstrating the processing in a signal processing unit. この発明の人検知装置の第2実施形態の概略構成図である。It is a schematic block diagram of the 2nd Embodiment of the person detection device of this invention.

以下、図を参照して、この発明の実施形態について説明するが、各構成要素の形状、大きさ及び配置関係については、この発明が理解できる程度に概略的に示したものに過ぎない。また、以下、この発明の好適な構成例につき説明するが、各構成要素の材質及び数値的条件などは、単なる好適例にすぎない。従って、この発明は以下の実施の形態に限定されるものではなく、この発明の構成の範囲を逸脱せずにこの発明の効果を達成できる多くの変更又は変形を行うことができる。 Hereinafter, embodiments of the present invention will be described with reference to the drawings, but the shape, size, and arrangement of each component are only schematically shown to the extent that the present invention can be understood. Further, although a preferable configuration example of the present invention will be described below, the material and numerical conditions of each component are merely suitable examples. Therefore, the present invention is not limited to the following embodiments, and many modifications or modifications can be made that can achieve the effects of the present invention without departing from the scope of the constitution of the present invention.

(第1実施形態)
図2を参照して、人検知装置及び人検知方法の第1実施形態について説明する。図2は、この発明の人検知装置の第1実施形態の概略構成図である。人検知装置は、住居や工場のセキュリティや店舗等でのマーケティングのために設置される。
(First Embodiment)
The first embodiment of the person detection device and the person detection method will be described with reference to FIG. FIG. 2 is a schematic configuration diagram of a first embodiment of the human detection device of the present invention. The person detection device is installed for security of houses and factories and marketing in stores.

人検知装置は、センサ部1と、信号処理部2と、表示部3と、送信部4とを備えて構成される。 The human detection device includes a sensor unit 1, a signal processing unit 2, a display unit 3, and a transmission unit 4.

センサ部1は、複数の静電気測定器11を備えて構成される。静電気測定器11は、例えば、電位センサ、電界センサ、誘導電流監視センサのいずれかで構成することができる。静電気測定器11は、それぞれ、2枚の平行な電極を備えて構成される。2枚の電極の間の電位差を測定する場合は、いわゆる電界センサとして機能する。また、一方の電極を接地して他方の電極の電位を測定する場合は、いわゆる電位センサとして機能する。また、このとき接地した電極に発生する誘導電流を測定すれば、いわゆる誘導電流監視センサとして機能する。以下の説明では、静電気測定器11が電界センサである場合の例を説明する。 The sensor unit 1 is configured to include a plurality of static electricity measuring instruments 11. The static electricity measuring device 11 can be configured by, for example, any of a potential sensor, an electric field sensor, and an induced current monitoring sensor. Each of the static electricity measuring instruments 11 is configured to include two parallel electrodes. When measuring the potential difference between two electrodes, it functions as a so-called electric field sensor. Further, when one electrode is grounded and the potential of the other electrode is measured, it functions as a so-called potential sensor. Further, if the induced current generated in the grounded electrode is measured at this time, it functions as a so-called induced current monitoring sensor. In the following description, an example in which the static electricity measuring instrument 11 is an electric field sensor will be described.

電界センサである静電気測定器11では、電界に比例した大きさの電圧信号が発生する。この電圧信号は、ボルテージフォロア12、電源周波数除去フィルタ13、バンドパスフィルタ14、増幅回路15、AD変換器16を経て、信号処理部2に送られる。 In the static electricity measuring instrument 11 which is an electric field sensor, a voltage signal having a magnitude proportional to the electric field is generated. This voltage signal is sent to the signal processing unit 2 via the voltage follower 12, the power supply frequency removal filter 13, the bandpass filter 14, the amplifier circuit 15, and the AD converter 16.

電源周波数除去フィルタ13は、電源周波数(50Hz又は60Hz)、及び、その整数倍(100Hz、150Hz,200Hz・・・、又は120Hz、180Hz,240Hz・・・)の成分を取り除く。 The power supply frequency removal filter 13 removes components of the power supply frequency (50 Hz or 60 Hz) and its integral multiples (100 Hz, 150 Hz, 200 Hz ..., Or 120 Hz, 180 Hz, 240 Hz ...).

バンドパスフィルタ14は、5〜1000Hzの信号を透過させる。透過させる帯域を1000Hz以下としたのは、後段のAD変換器16でのAD変換により生じる可能性のあるエリアシングを避けるためである。一方、人の歩行時などに生じる電界変動を測定するために、透過させる帯域を5Hz以上としている。 The bandpass filter 14 transmits a signal of 5 to 1000 Hz. The transmission band is set to 1000 Hz or less in order to avoid aliasing that may occur due to AD conversion by the AD converter 16 in the subsequent stage. On the other hand, in order to measure the electric field fluctuation that occurs when a person walks, the transmission band is set to 5 Hz or more.

人体を初めとする被測定対象が空間を移動するとき、電荷や誘電率の空間分布が変化し、それによって電位分布が変化する。ところが人体が移動するときの速度は通常1m/s程度のゆっくりしたものであり、それによる電界の変化は通常1Hz以上の低周波の領域に現れる。従って、透過帯域が5〜1000Hzのバンドパスフィルタ14を透過した後の信号には、人体の単純な水平移動による電界の変化は、ほとんど現れない。 When the object to be measured, including the human body, moves in space, the spatial distribution of electric charge and permittivity changes, and the potential distribution changes accordingly. However, the speed at which the human body moves is usually slow at about 1 m / s, and the change in the electric field due to this usually appears in a low frequency region of 1 Hz or higher. Therefore, in the signal after passing through the bandpass filter 14 having a transmission band of 5 to 1000 Hz, the change in the electric field due to the simple horizontal movement of the human body hardly appears.

しかしながら、例えば人の歩行時には周波数8±2Hz程度の電界の変動が生じることが知られている。このような電界の変動は、歩行者の靴の床面を離れるとき、また着地するときに発生する。 However, it is known that, for example, when a person walks, the electric field fluctuates at a frequency of about 8 ± 2 Hz. Such fluctuations in the electric field occur when the pedestrian's shoes leave the floor and when they land.

床面は良好な導電体、電気的アースと見なせる地面の上に敷かれ、通常ある程度の電気抵抗を持つが、通常の床面素材においては床上の電荷は1ミリ秒より小さな時定数を持って地面に逃げるので、ほぼ地面と同じくアースと見なせる。 The floor is laid on a ground that can be regarded as a good conductor, electrical ground, and usually has some electrical resistance, but in normal floor materials, the charge on the floor has a time constant less than 1 millisecond. Since it escapes to the ground, it can be regarded as an earth almost like the ground.

歩行者の靴が床面を離れるとき、また着地するときに発生する電界変化は、電気的アースとしての床面と帯電した導電体としての人体の間の電気容量が、靴の離地、着地とともに増減し、結果として人体の電位が増減することが主な要因である。つまり人体の帯電量をQ、アース(床面)に対する人体の電位をV、アースと人体の間の電気容量をCとすると以下の式(5)の関係がある。 The change in electric field that occurs when a pedestrian's shoes leave the floor and land is due to the electrical capacity between the floor as an electrical ground and the human body as a charged conductor, which causes the shoes to take off and land. The main factor is that the electric potential of the human body increases or decreases as a result. That is, if the amount of charge of the human body is Q, the potential of the human body with respect to the ground (floor surface) is V, and the electric capacity between the ground and the human body is C, the following equation (5) is obtained.

Figure 0006763241
Figure 0006763241

靴面が床から離地したり着地したりする瞬間に、Cが急激に変化するので、これに従ってVが変動し、電界も変化する。 At the moment when the shoe surface takes off or lands on the floor, C changes abruptly, so V fluctuates accordingly and the electric field also changes.

他に、人体から5Hzよりも高速な電界が発生する原因としては、人体からのあるいは人体への放電がある。人体は歩行などにより数100Vから数kVに帯電しているが、別の電位を持つ物体に接触したとき両者の電位が等しくなる方向に電荷の移動が生じる。例えば、手など露出した部分が接触したとき人体の電位は急激にアース電位まで低下するが、その時定数は10ナノ秒のオーダーである。また人同士が触れ合った場合、両者の電位は等しくなるが、この現象も急激に生じる。 Another cause of the electric field generated from the human body at a speed higher than 5 Hz is a discharge from or to the human body. The human body is charged from several hundred V to several kV by walking or the like, but when it comes into contact with an object having another potential, the charge moves in the direction in which both potentials become equal. For example, when an exposed part such as a hand comes into contact, the electric potential of the human body suddenly drops to the earth potential, and the time constant is on the order of 10 nanoseconds. Also, when people come into contact with each other, their potentials become equal, but this phenomenon also occurs rapidly.

従って、静電気測定器11において発生した信号においてバンドパスフィルタ14を通過する信号は、検知範囲内における人体などの被測定対象と、測定対象とは異なる物体、例えば床面、他の人体など、の相互作用により生じるものである。この静電気測定器11で発生する信号には、以下に述べるようないくつかの特徴がある。 Therefore, in the signal generated by the static electricity measuring device 11, the signal passing through the bandpass filter 14 is the object to be measured such as a human body within the detection range and an object different from the measurement target such as a floor surface or another human body. It is caused by interaction. The signal generated by the static electricity measuring instrument 11 has some features as described below.

第一の特徴は、特に歩行時や人体の放電時において、人体の電荷や電位の変化がダイナミックであり、電界の変化量も大きいということである。歩行時の離地、着地における人体の電位は、静止時の電位の数倍となる。また人体がアースに対してゼロ電位である場合でも、靴と床面での帯電現象によって人体の電位が上昇することが知られている。人体の電位は最大で数kVにおよび、通常の電位センサ、電界センサ、電極への誘導電流の観測などで遠く離れた場所でも信号の観測が容易に可能である。 The first feature is that changes in the electric charge and potential of the human body are dynamic, and the amount of change in the electric field is also large, especially when walking or discharging the human body. The electric potential of the human body during takeoff and landing during walking is several times higher than the electric potential at rest. It is also known that even when the human body has a zero potential with respect to the earth, the potential of the human body rises due to the charging phenomenon on the shoes and the floor surface. The potential of the human body reaches a maximum of several kV, and it is possible to easily observe a signal even at a distant place by observing an induced current to an ordinary potential sensor, an electric field sensor, or an electrode.

第二の特徴は、信号強度の距離に対する減衰がゆるやかなことである。これに対して被測定対象が空間を移動するときに生じる電界の変化は、図1を用いて説明したように、センサと被測定対象の距離の3乗に反比例して小さくなる。ところが被測定対象と、それとは異なる物体、例えば床面、他の人体などとの相互作用により生じる電界の変化においては、被測定対象自体の電位が変動することによって生じるものなので、電界の変化はセンサと被測定対象の距離の2乗に反比例する。したがって、センサから離れた場所においても信号の観測が可能である。 The second feature is that the signal strength decays slowly with respect to the distance. On the other hand, the change in the electric field that occurs when the object to be measured moves in space becomes smaller in inverse proportion to the cube of the distance between the sensor and the object to be measured, as described with reference to FIG. However, the change in the electric field caused by the interaction between the object to be measured and an object different from it, such as the floor surface or another human body, is caused by the fluctuation of the potential of the object to be measured itself, so that the change in the electric field is caused. It is inversely proportional to the square of the distance between the sensor and the object to be measured. Therefore, it is possible to observe the signal even at a place away from the sensor.

第三の特徴は、あるイベントに対して得られる各々のセンサからの出力は、信号の大きさがイベントの発生位置からセンサへの距離によって異なっても、信号の波形は相関係数1で相関することである。また、十分低い周波数での電界変化は空間を瞬時に伝導すると見なせるので(準静的近似による)、それらの相関は時間差なしに生じる。第三の特徴より、各センサにより得られた信号から、時間差なしに相関する信号を抽出して解析することによって、イベントの空間的位置を推定することができる。本実施形態においてイベントの空間位置の推定は、信号処理部2によって行われる。 The third feature is that the output from each sensor obtained for a certain event has a correlation coefficient of 1 even if the signal magnitude differs depending on the distance from the event occurrence position to the sensor. It is to be. Also, since electric field changes at sufficiently low frequencies can be considered to conduct space instantaneously (by quasi-static approximation), their correlation occurs without time lag. From the third feature, the spatial position of the event can be estimated by extracting and analyzing the correlated signals without a time difference from the signals obtained by each sensor. In the present embodiment, the estimation of the spatial position of the event is performed by the signal processing unit 2.

AD変換器16は、それぞれのセンサにより得られた信号からサンプリングレート1kHzで強度信号列Sを生成する。この実施形態では強度信号列Sの時間窓を0.1ミリ秒とする。時間窓を0.1ミリ秒とすると、それぞれの強度信号列Sは以下の式(6)のように100列(0.1ms×1000Hz=100)の数列として表される。 The AD converter 16 generates an intensity signal sequence S from the signals obtained by the respective sensors at a sampling rate of 1 kHz. In this embodiment, the time window of the intensity signal sequence S is set to 0.1 milliseconds. Assuming that the time window is 0.1 ms, each intensity signal sequence S is represented as a sequence of 100 sequences (0.1 ms × 1000 Hz = 100) as shown in the following equation (6).

Figure 0006763241
Figure 0006763241

ここで、時間窓は、サンプリングレートで与えられる時間であり、時間枠は、強度信号列Sの要素の数を時間窓に乗算して与えられる。また、ここでは、センサの数Mを6とした例を説明する。 Here, the time window is the time given at the sampling rate, and the time frame is given by multiplying the time window by the number of elements of the intensity signal sequence S. Further, here, an example in which the number M of the sensors is 6 will be described.

信号処理部2は、上式(6)の強度信号列S同士の相関関係より、時間枠内で生じたイベントについて、その発生件数ξと発生時刻、発生箇所を検知する。 The signal processing unit 2 detects the number of events ξ, the time of occurrence, and the location of occurrence of the event that occurred within the time frame from the correlation between the intensity signal sequences S of the above equation (6).

図3を参照して、特に、信号処理部2での処理を説明する。図3は、信号処理部2での処理、時間枠内で生じたイベントについて、その発生件数ξと発生時刻、発生箇所を検知するための処理を説明するための処理フロー図である。 In particular, the processing in the signal processing unit 2 will be described with reference to FIG. FIG. 3 is a processing flow diagram for explaining the processing by the signal processing unit 2 and the processing for detecting the occurrence number ξ, the occurrence time, and the occurrence location of the event that occurred within the time frame.

信号処理部2が備えるイベント列取得手段21は、先ず、イベントが発生したか否かを判定する(S10)。ある時間枠について、イベント列取得手段21は、取得したどの強度信号列S1〜S6からも信号が得られない場合、すなわち(6)式において、強度信号列S1〜S6の要素が全て0である場合、時間枠にはイベントは発生していない、すなわちξ=0と判断し、この場合(No)は何も行なわず、次の時間枠での判定を行う。 The event sequence acquisition means 21 included in the signal processing unit 2 first determines whether or not an event has occurred (S10). For a certain time frame, the event sequence acquisition means 21 cannot obtain a signal from any of the acquired intensity signal sequences S1 to S6, that is, in the equation (6), all the elements of the intensity signal sequences S1 to S6 are 0. In this case, it is determined that no event has occurred in the time frame, that is, ξ = 0. In this case (No), nothing is performed and the determination in the next time frame is performed.

S10での判定の結果、強度信号列S1〜S6のいずれかの信号が得られた場合(Yes)、イベント列取得手段21は、イベントが発生したものとして強度信号列S1〜S6のうち、信号が得られていないものを除き、信号間で相関係数を計算する(S20)。Sn及びSmの相関係数rnmは以下の式(7)で計算される。上線つきSn及びSmは、それぞれSn及びSmの平均値である。 When any signal of the intensity signal sequences S1 to S6 is obtained as a result of the determination in S10 (Yes), the event sequence acquisition means 21 assumes that an event has occurred and is a signal of the intensity signal sequences S1 to S6. The correlation coefficient is calculated between the signals, except for those for which is not obtained (S20). The correlation coefficient r nm of Sn and Sm is calculated by the following equation (7). The overlined Sn and Sm are the average values of Sn and Sm, respectively.

Figure 0006763241
Figure 0006763241

次に、この相関関数rnmが、所定の閾値より大きいか否かを判定する(S30)。例えば、相関関数rnmが0.9以上など、それらの相関関数rnmが全て1と見なせる場合、時間窓内のイベントの発生件数(イベント数)ξは1件、すなわちξ=1と判断される。この場合(Yes)、以下の(8)式のように全ての強度信号列S1〜S6は一個の信号列Ivent1との比例関係で表される。 Next, it is determined whether or not this correlation function r nm is larger than a predetermined threshold value (S30). For example, when the correlation function r nm is 0.9 or more and all of the correlation functions r nm can be regarded as 1, the number of events occurring in the time window (number of events) ξ is determined to be 1, that is, ξ = 1. To. In this case (Yes), all the intensity signal sequences S1 to S6 are represented by a proportional relationship with one signal sequence Ivent1 as shown in the following equation (8).

Figure 0006763241
Figure 0006763241

上式(8)の右辺の係数行列の要素c〜cは全て0以上の実数であり、各信号間の強度比をあらわす。信号処理部2が備える位置推定手段22は、この比例定数を求め、以下の式(9)によりイベントの位置を特定する(S40)。 The elements c 1 to c 6 of the coefficient matrix on the right side of the above equation (8) are all real numbers of 0 or more, and represent the intensity ratio between each signal. The position estimation means 22 included in the signal processing unit 2 obtains this proportionality constant and specifies the position of the event by the following equation (9) (S40).

Figure 0006763241
Figure 0006763241

ここでL〜Lは、イベントの発生位置からセンサS1〜S6までの距離である。ただし、係数c〜cのうちゼロであるものに関しては、上式(9)の比例関係から除く。例えばc=0,c=0の場合は、以下の式(10)のように表して、この式(10)からイベントの位置を特定する。 Here, L 1 to L 6 are distances from the event occurrence position to the sensors S 1 to S 6 . However, the coefficients c 1 to c 6 that are zero are excluded from the proportional relationship in the above equation (9). For example, when c 5 = 0 and c 6 = 0, the position of the event is specified from this equation (10) by expressing it as the following equation (10).

Figure 0006763241
Figure 0006763241

上式(9)、上式(10)は、センサの感じる電界変化の強度が、センサとイベントの生じた距離の2乗に反比例することを利用したものである。しかしながら本実施形態においてL〜Lが、係数c〜cの2乗に反比例でなく2.1乗に反比例していると想定している。これは、床面がアースされていることの影響を考慮したためである。L〜Lに関する情報により位置推定手段22は、イベント発生位置の座標(x、y)を導き、イベントの発生時刻tとともに、第1のメモリ23に記録する。 The above equations (9) and (10) utilize the fact that the intensity of the electric field change felt by the sensor is inversely proportional to the square of the distance between the sensor and the event. However, in this embodiment, it is assumed that L 1 to L 6 are not inversely proportional to the square of the coefficients c 1 to c 6 but inversely proportional to the 2.1 power. This is because the influence of the grounding of the floor surface is taken into consideration. The position estimation means 22 derives the coordinates (x, y) of the event occurrence position based on the information regarding L 1 to L 6 , and records them in the first memory 23 together with the event occurrence time t.

さらに、信号処理部2が備える仮想人割当手段24は、イベント発生位置の座標(x、y)とイベントの発生時刻tを、既に定義された、あるいはその時点で定義した仮想人のデータとして割り当て(S200)、第2のメモリ25に記録する。既に定義された仮想人に新たなイベント(位置、発生時刻)を割り当てる場合は、既に割り当てられたイベントを上書きする形で行う。これはその仮想人の位置が更新されたことを示す。 Further, the virtual person allocating means 24 included in the signal processing unit 2 allocates the coordinates (x, y) of the event occurrence position and the event occurrence time t as virtual person data already defined or defined at that time. (S200), recording is performed in the second memory 25. When assigning a new event (position, occurrence time) to an already defined virtual person, the event that has already been assigned is overwritten. This indicates that the position of the virtual person has been updated.

仮想人の所在場所は、仕様態様に応じて表示部3により可視化して表示され、あるいは、送信部4によりインターネット上に送信される。 The location of the virtual person is visualized and displayed by the display unit 3 according to the specification mode, or is transmitted on the Internet by the transmission unit 4.

以上は、強度信号列S間の相関係数が1で、時間枠中に検知エリア内で一件のみのイベントが発生したと判断された場合である。 The above is the case where the correlation coefficient between the intensity signal sequences S is 1, and it is determined that only one event has occurred in the detection area during the time frame.

S30での判定の結果、強度信号列S間の相関係数が所定の閾値以下である、すなわち1とみなせないと判断される場合(No)、時間枠内に発生したイベントξを一件増やしてξ=2として(S50)、信号列S1〜S6を説明できるか試みる。そのため、イベント列取得手段21は、イベント信号列の探索を行い、M行ξ列(ここでは、6行2列)の行列Aのうちで、以下の式(11)の演算結果として得られた二つのイベント信号列I、Iの間の相関係数が最小となるような行列Aの係数を決定する(S60)。行列Aの係数を決定する手法については、それらの係数が複素数である場合を含み、任意好適な従来公知の方法を用いればよい。 As a result of the determination in S30, when it is determined that the correlation coefficient between the intensity signal sequences S is equal to or less than a predetermined threshold value, that is, it cannot be regarded as 1 (No), the event ξ generated within the time frame is increased by one. Let ξ = 2 (S50) and try to explain the signal sequences S1 to S6. Therefore, the event column acquisition means 21 searches for the event signal sequence, and is obtained as the calculation result of the following equation (11) in the matrix A of M rows ξ columns (here, 6 rows and 2 columns). The coefficient of the matrix A is determined so that the correlation coefficient between the two event signal sequences I 1 and I 2 is minimized (S60). As a method for determining the coefficients of the matrix A, any suitable conventionally known method may be used, including the case where those coefficients are complex numbers.

Figure 0006763241
Figure 0006763241

ここで行列Aの12個の要素は全て0か0より大きい実数であるように限定することができる。そのような限定が可能なのは、イベント信号列I、Iの変化は時間差なしに強度信号列S1〜S6に伝わり、その影響は強度信号列S1〜S6と信号源との距離によって増減するが、符号が逆転することはないからである。この限定によってより少ない計算量で行列Aの要素を決定することができる。 Here, the 12 elements of the matrix A can all be limited to 0 or a real number greater than 0. Such a limitation is possible because the changes in the event signal sequences I 1 and I 2 are transmitted to the intensity signal sequences S1 to S6 without a time difference, and the influence thereof increases or decreases depending on the distance between the intensity signal sequences S1 to S6 and the signal source. , This is because the sign is not reversed. This limitation allows the elements of matrix A to be determined with less computational complexity.

行列Aの要素及びそのときのイベント信号列I、Iが決定された後、さらにイベント列取得手段21は、そのときのイベント信号列I、Iと、その間の相関係数r(2)を記憶する(S70)。このイベント信号列I、Iは、それぞれのイベントを表すと想定される信号列である。 After the elements of the matrix A and the event signal sequences I 1 and I 2 at that time are determined, the event sequence acquisition means 21 further determines the event signal sequences I 1 and I 2 at that time and the correlation coefficient r between them ( 2) is stored (S70). The event signal sequences I 1 and I 2 are signal sequences that are supposed to represent the respective events.

なお、S30における閾値については、イベント数ξが1の場合と2の場合のサンプルを用いて、予め決定しておくことができる。 The threshold value in S30 can be determined in advance by using the samples when the number of events ξ is 1 and 2.

続けてイベント列取得手段21は、時間枠内に発生したイベントξを一件増やしてξ=3とする(S80)。そして、イベント信号列の探索を行い、M行ξ列(ここでは、6行3列)の行列Aにおいて、以下の式(12)の演算結果として得られた三つのイベント信号列I、I、Iの間の相関係数の絶対値の平均値が最小となるような行列Aの要素を決定する(S90)。 Subsequently, the event sequence acquisition means 21 increases the number of events ξ that occurred within the time frame by one so that ξ = 3 (S80). Then, the event signal sequence is searched, and in the matrix A of M rows and ξ columns (here, 6 rows and 3 columns), the three event signal sequences I 1 and I obtained as the calculation results of the following equation (12) are obtained. 2. Determine the elements of the matrix A such that the average value of the absolute values of the correlation coefficients between I 3 is minimized (S90).

Figure 0006763241
Figure 0006763241

すなわち、上式(12)で得られるイベント信号列I、I、Iから、I、I間の相関係数、I、I間の相関係数、I、I間の相関係数の3つの相関係数の絶対値の算術平均を計算し、その値が最小となる行列Aを決定する。そのような行列Aの要素及びそのときのイベント信号列I〜Iが決定された後、イベント列取得手段21は、そのときのイベント信号列I〜Iと、その間の相関係数r(3)を第1のメモリ23に記録する(S100)。 That is, from the event signal sequences I 1 , I 2 , I 3 obtained by the above equation (12), the correlation coefficient between I 1 and I 2 , the correlation coefficient between I 2 and I 3 , I 3 and I 1 The arithmetic mean of the absolute values of the three correlation coefficients between them is calculated, and the matrix A having the smallest value is determined. After the elements of such a matrix A and the event signal sequences I 1 to I 3 at that time are determined, the event sequence acquisition means 21 has the event signal sequences I 1 to I 3 at that time and the correlation coefficient between them. Record r (3) in the first memory 23 (S100).

次に、イベント列取得手段21は、r(ξ)とr(ξ−1)を比較する(S110)。ここで、r(ξ)>r(ξ−1)、すなわち、r(3)がr(2)より大きいとき(Yes)、ξ=2、すなわち、時間枠内に生じたイベントの数は2個と判断する。r(3)がr(2)より小さいとき(No)、さらに時間枠内に生じる仮定イベントの数を増やしξ=4として(S80)、任意の6行4列の行列Aにおいて、以下の式(13)の演算結果として得られた4つの信号列I〜Iの間の相関係数の絶対値の平均値が最小となるような行列Aの係数を決定する(S90)。 Next, the event sequence acquisition means 21 compares r (ξ) and r (ξ-1) (S110). Here, when r (ξ)> r (ξ-1), that is, when r (3) is greater than r (2) (Yes), ξ = 2, that is, the number of events that occurred within the time frame is 2. Judge as an individual. When r (3) is smaller than r (2) (No), the number of assumed events occurring in the time frame is further increased to ξ = 4 (S80), and in an arbitrary 6-by-4 matrix A, the following equation is used. The coefficient of the matrix A is determined so that the average value of the absolute values of the correlation coefficients between the four signal sequences I 1 to I 4 obtained as the calculation result of (13) is minimized (S90).

Figure 0006763241
Figure 0006763241

イベント列取得手段21は、最小算術平均r(4)を、第1のメモリ23に記録する(S100)。またそのときのI〜Iを記録する。r(ξ)とr(ξ−1)を比較し(S110)、r(4)がr(3)より大きいとき(Yes)、ξ=3、すなわち、時間枠内に生じるイベントの数は3個と判断する。r(4)がr(3)以下のとき(No)は、仮定するイベントを一個増やす。 The event sequence acquisition means 21 records the minimum arithmetic mean r (4) in the first memory 23 (S100). Also, record I 1 to I 4 at that time. Comparing r (ξ) and r (ξ-1) (S110), when r (4) is greater than r (3) (Yes), ξ = 3, that is, the number of events occurring in the time frame is 3. Judge as an individual. When r (4) is r (3) or less (No), the number of assumed events is increased by one.

以下同様に、S80〜S110をイベント数ξの相関関数の絶対値の平均値r(ξ)が、イベント数ξが1少ない場合の相関関数の絶対値の平均値r(ξ−1)より大きくなるまでくりかえす。r(ξ)r(ξ−1)となった時点でイベント数を決定する。つまりこの実施形態では、イベント数の増加に対して相関係数の絶対値の平均が極小値をとるとき、そのイベント数をもって時間枠内のイベント数ξと決定する。 Similarly, in S80 to S110, the average value r (ξ) of the absolute value of the correlation function of the number of events ξ is larger than the average value r (ξ-1) of the absolute values of the correlation function when the number of events ξ is 1 less. Repeat until it becomes. The number of events is determined when r (ξ) > r (ξ-1). That is, in this embodiment, when the average of the absolute values of the correlation coefficients takes the minimum value with respect to the increase in the number of events, the number of events is determined as the number of events ξ in the time frame.

以上の過程により時間内で生じたイベント数ξと、それぞれのイベントを示すイベント信号列Iが得られる。イベントの数ξがn個の場合は、以下の式(14)となる。
Through the above process , the number of events ξ that occurred in the time frame and the event signal sequence I indicating each event are obtained. When the number of events ξ is n, the following equation (14) is obtained.

Figure 0006763241
Figure 0006763241

上式(14)から、イベント信号列の相関が最小となる場合の式(15)が求まる。 From the above equation (14), the equation (15) when the correlation of the event signal sequence is minimized can be obtained.

Figure 0006763241
Figure 0006763241

さらに、もとの信号S1、S2、S3、S4、S5、S6を再構成するような6行n列の行列を求める。すなわち、以下の式(16)を満たす行列A´を求める(S120)。 Further, a 6-row, n-column matrix that reconstructs the original signals S1, S2, S3, S4, S5, and S6 is obtained. That is, the matrix A'that satisfies the following equation (16) is obtained (S120).

Figure 0006763241
Figure 0006763241

行列A´の第一列の係数の比により、第一のイベントの発生位置を求め、第二列の係数の比により、第二のイベントの発生位置を求め、以下同様にn個のイベントそれぞれの発生位置を求める(S130)。 The position where the first event occurs is obtained from the ratio of the coefficients in the first column of the matrix A', the position where the second event occurs is obtained from the ratio of the coefficients in the second column, and so on, for each of the n events. (S130).

さらに、信号処理部2が備える仮想人割当手段24は、イベント発生位置の座標(x、y)とイベントの発生時刻tを、既に定義された、あるいはその時点で定義した仮想人のデータとして割り当て(S200)、第2のメモリ25に記録する。既に定義された仮想人に新たなイベント(の位置、発生時刻)を割り当てる場合は、既に割り当てられたイベントを上書きする形で行う。これはその仮想人の位置が更新されたことを示す。 Further, the virtual person allocating means 24 included in the signal processing unit 2 allocates the coordinates (x, y) of the event occurrence position and the event occurrence time t as virtual person data already defined or defined at that time. (S200), recording is performed in the second memory 25. When assigning a new event (position, occurrence time) to an already defined virtual person, the already assigned event is overwritten. This indicates that the position of the virtual person has been updated.

以上説明したように、この第1実施形態では、複数の電界センサより生じる信号から信号間の相関する成分を抽出することによって、人体の歩行、他物体への接触などに生じる電界の変動イベントを抽出し、イベントが生じた位置を推定して、仮想人にイベントを割り当てることによって所定の時刻の人体の位置を推定できる。また、複数のイベントをイベントごとに分離できるので、複数の人体が存在する場合においても、それらの移動を個人ごとに識別できる。 As described above, in the first embodiment, by extracting the correlated components between the signals from the signals generated by the plurality of electric field sensors, the fluctuation event of the electric field generated in the walking of the human body, the contact with other objects, etc. is generated. By extracting, estimating the position where the event occurred, and assigning the event to the virtual person, the position of the human body at a predetermined time can be estimated. Further, since a plurality of events can be separated for each event, even when a plurality of human bodies exist, their movements can be identified for each individual.

また、ここで得られる信号は、人間の歩行、他物体への接触などにより生じる電界の変動イベントであり、人体の平行移動により生じる電界の変動イベントより変動量が大きい。このため、イベント発生箇所と電界センサへの距離が離れることに伴う電界変動の減少が少ないので、遠方のイベントも検知できる、従って、人検知装置において用いるセンサの数が少なくて済む、あるいは特別に高感度のセンサを用いる必要がないといった効果がある。 Further, the signal obtained here is an electric field fluctuation event caused by human walking, contact with another object, or the like, and the fluctuation amount is larger than the electric field fluctuation event caused by the parallel movement of the human body. Therefore, since the decrease in the electric field fluctuation due to the distance between the event occurrence location and the electric field sensor is small, it is possible to detect a distant event. Therefore, the number of sensors used in the human detection device can be small, or specially. There is an effect that it is not necessary to use a high-sensitivity sensor.

また、電界センサを床に垂直な方向の電界を検知するように設置すると、センサに対して床面内のどの方角で電界を変動させるイベントが生じても、センサは平等に電界を検知する、すなわち指向性がない。したがって信号処理回路によるイベントの位置推定がより精密になるという効果がある。 In addition, if the electric field sensor is installed so as to detect the electric field in the direction perpendicular to the floor, the sensor detects the electric field equally regardless of the direction in which the electric field fluctuates with respect to the sensor. That is, there is no directionality. Therefore, there is an effect that the position estimation of the event by the signal processing circuit becomes more precise.

(第2実施形態)
図4を参照して、人検知装置の第2実施形態について説明する。図4は、この発明の人検知装置の第2実施形態の概略構成図である。
(Second Embodiment)
A second embodiment of the human detection device will be described with reference to FIG. FIG. 4 is a schematic configuration diagram of a second embodiment of the human detection device of the present invention.

第2実施形態は、低周波成分を合わせて利用する点が第1実施形態と異なっている。他の構成及び動作は、第1実施形態と同様なので、重複する説明及び図示を省略することもある。 The second embodiment is different from the first embodiment in that low frequency components are used together. Since other configurations and operations are the same as those in the first embodiment, duplicate description and illustration may be omitted.

第2実施形態ではバンドパスフィルタ14により5Hzより大きな高周波数成分と5Hzより小さな低周波数成分を分離して、それぞれ異なる信号処理を行う。高周波数成分は第1実施形態と同様な処理を行う。また、低周波成分については、周波数帯域が低周波であることを除いて、増幅回路115、AD変換器116を経て低周波信号処理部122に送られるまでの構成及び動作は、高周波成分と同様である。 In the second embodiment, the bandpass filter 14 separates the high frequency component larger than 5 Hz and the low frequency component smaller than 5 Hz, and performs different signal processing for each. The high frequency component is processed in the same manner as in the first embodiment. Further, regarding the low frequency component, the configuration and operation until it is sent to the low frequency signal processing unit 122 via the amplifier circuit 115 and the AD converter 116 are the same as those of the high frequency component, except that the frequency band is low frequency. Is.

高周波信号処理部120の動作は、第1実施形態の信号処理部2と同様である。高周波成分により、人の配置、例として点P1、P2、P3にて人の存在が認知されたとき、近傍の三つのセンサの低周波成分(すなわち、DC成分)から各人の総電荷や電位を推定することができる。点P1、P2、P3での電界E、E、Eはそれぞれ以下の式(17)で表される。 The operation of the high frequency signal processing unit 120 is the same as that of the signal processing unit 2 of the first embodiment. When the presence of a person is recognized by the high frequency component, for example, at points P1, P2, and P3, the total charge and potential of each person from the low frequency components (that is, DC components) of three nearby sensors. Can be estimated. The electric fields E 1 , E 2 , and E 3 at points P1, P2, and P3 are represented by the following equations (17), respectively.

Figure 0006763241
Figure 0006763241

ここでQ、Q、QはP1、P2、P3点にいる人の電荷、LijはセンサSからPまでの距離である。Lijは、高周波成分の信号解析により既知であるから、上式(17)の連立方程式を解くことによって、Q、Q、Qを求めることができる。また人体の容量成分を与えることによって、各人の電位を推測することができる。人体の電荷や電位は、各人についての情報として、各人に結び付けて各人の識別に用いたり、また、手など露出した部分が何かに接触したなど各人の活動のバロメータとして用いたりできる。また第1実施形態において、人の歩行などを捕らえそこなったり、誤認したりした場合においても、低周波成分での電界の低下を捉えることで、高周波による情報を訂正できる。 Here Q 1, Q 2, Q 3 is the charge of the person in the P1, P2, P3 point, L ij is the distance from the sensor S i until P j. L ij is because it is known by the signal analysis of the high-frequency component, by solving the simultaneous equations of the equation (17), it can be obtained Q 1, Q 2, Q 3 . Moreover, the electric potential of each person can be estimated by giving the volume component of the human body. The electric charge and electric potential of the human body can be used as information about each person by linking them to each person to identify each person, or as a barometer of each person's activity such as when an exposed part such as a hand comes into contact with something. it can. Further, in the first embodiment, even when a person's walking or the like is missed or misidentified, the information due to the high frequency can be corrected by capturing the decrease in the electric field in the low frequency component.

第2実施形態によれば、電界センサの出力をフィルタリングし、低周波成分の信号と高周波成分のもたらす人の所在情報を併せて用いたので、各人の総電荷、電位などを得て各人の識別に用いたりするなど活用可能である。また高周波成分のもたらす人の所在情報の誤りを訂正し、より高度な人検知が可能である。 According to the second embodiment, the output of the electric field sensor is filtered, and the signal of the low frequency component and the location information of the person provided by the high frequency component are used together, so that the total charge, potential, etc. of each person are obtained and each person is obtained. It can be used for identification of. In addition, it is possible to correct the error of the location information of the person caused by the high frequency component and to detect the person at a higher level.

(その他の実施形態)
第1実施形態、第2実施形態では、電位(電界)の変化を測定するのに、電界センサを用いる例を説明したが、電位(電界)の変化を検知できるものであれば、その種類はそれらにとどまらない。例えば、電極を接地して該電極に流入する誘導電流を測ることも可能であるし、電位センサや電界センサの電位、電界に比例する出力だけでなく、それらのセンサの動作により生じる電流、電圧を測定するセンサであることも可能である。
(Other embodiments)
In the first embodiment and the second embodiment, an example in which an electric field sensor is used to measure a change in an electric potential (electric field) has been described, but if the change in the electric potential (electric field) can be detected, the type is It doesn't stop there. For example, it is possible to measure the induced current flowing into the electrode by grounding the electrode, and not only the potential of the potential sensor or the electric field sensor, the output proportional to the electric field, but also the current and voltage generated by the operation of those sensors. It can also be a sensor that measures.

例えば、電界センサを用いたのに対して、電位センサを用いると、それぞれの電位センサには電位に比例する出力と電位の時間変化に比例する出力の2つの出力がある。電位の時間変化に比例する出力は5Hzから1000Hzの間の周波数を透過するバンドパスフィルタを透過したのち、第1実施形態、第2実施形態と同様に、人間の歩行、他物体への接触などに生じるイベントを抽出する。電位に比例する出力は第2実施形態の低周波数の出力のように、前者で得た人配置に関する情報を捕捉、訂正するのに用いられる。 For example, when an electric field sensor is used, when a potential sensor is used, each potential sensor has two outputs, an output proportional to the potential and an output proportional to the time change of the potential. The output proportional to the change in potential with time passes through a bandpass filter that transmits frequencies between 5 Hz and 1000 Hz, and then, as in the first and second embodiments, human walking, contact with other objects, etc. Extract the events that occur in. The output proportional to the potential is used to capture and correct the information regarding the staffing obtained in the former, as in the low frequency output of the second embodiment.

この場合、電界ではなく電位の変動を観測するので以下の違いがある。すなわち、以下の式(18)に示すように、電界は距離Lの2乗に反比例して減少するが、電位はLに反比例して減少する。 In this case, there are the following differences because the fluctuation of the potential is observed instead of the electric field. That is, as shown in the following equation (18), the electric field decreases in inverse proportion to the square of the distance L, but the potential decreases in inverse proportion to L.

Figure 0006763241
Figure 0006763241

したがって、電位センサを用いると、センサから離れた場所と近接した場所で、センサの出力差が、第1実施形態、第2実施形態ほど大きく異ならなくなり、多くの場合位置の推定がより精密になるという利点がある。 Therefore, when the potential sensor is used, the output difference of the sensor does not differ significantly from the place away from the sensor and the place close to the sensor as in the first embodiment and the second embodiment, and the position estimation becomes more precise in many cases. There is an advantage.

イベント信号列を取得する処理については、多成分解析、独立成分分析など複数の信号から、それらをもたらした元の信号を推定する方法は多く提案されている。この発明では、その方法について限定されず、任意好適な従来公知のものを用いることができる。 Regarding the process of acquiring the event signal sequence, many methods of estimating the original signal that brought them from a plurality of signals such as multi-component analysis and independent component analysis have been proposed. In the present invention, the method is not limited, and any suitable conventionally known method can be used.

また、ある時点で生じたイベントをある仮想人に割り当てる規則、あるイベントと別のイベントが同一人によるものかどうか判断する規則は多く考えられる。例えば、歩行に関する信号が連続してある距離内で生じているとき、それは同一人に関するイベントとする、など考えられる。この具体的なアルゴリズムは任意好適な方法を用いればよい。 In addition, there are many conceivable rules for assigning an event that occurs at a certain point in time to a virtual person, and for determining whether a certain event and another event are due to the same person. For example, when a signal related to walking occurs continuously within a certain distance, it may be an event related to the same person. Any suitable method may be used for this specific algorithm.

1 センサ部
2 信号処理部
3 表示部
4 送信部
11 静電気測定器
12 ボルテージフォロア
13 電源周波数除去フィルタ
14 バンドパスフィルタ
15、115 増幅回路
16、116 AD変換器
120 高周波信号処理部
122 低周波信号処理部
1 Sensor unit 2 Signal processing unit 3 Display unit 4 Transmitter unit 11 Electrostatic measuring instrument 12 Voltage follower 13 Power frequency removal filter 14 Bandpass filter 15, 115 Amplifier circuit 16, 116 AD converter 120 High frequency signal processing unit 122 Low frequency signal processing Department

Claims (6)

複数のセンサから、信号強度を示す強度信号列を取得する過程と、
強度信号列間で相関係数を計算する過程と、
前記相関係数を所定の閾値と比較する過程と、
前記比較の結果、閾値以上である場合に行われる、信号強度比からイベント位置を特定する過程と、
閾値よりも小さい場合に行われる、複数イベント判定を行う過程と
を備え、
複数イベント判定を行う過程は、
イベント数ξを2として、
ξ個のイベント信号列Iと、ξ行M列の係数行列Aと、M行1列の強度信号列Sで与えられる、I=ASから、イベント信号列I間の相関係数r(ξ)が最小となる係数行列Aを決定する第1の過程と、
イベント数ξに1を加算して、ξ個のイベント信号列Iと、ξ行M列の係数行列Aと、M行1列の強度信号列Sで与えられる、I=ASから、イベント信号列I間の相関係数r(ξ)が最小となる係数行列Aを決定する第2の過程と、
r(ξ)とr(ξ−1)を比較する過程とを
備え、
r(ξ)がr(ξ−1)以下の場合は、さらに、前記第2の過程及びr(ξ)とr(ξ−1)を比較する過程を行い、
それ以外の場合は、ξ−1以下であるときに得られた係数行列Aの各行の係数からイベント位置を特定する
ことを特徴とする人検知方法。
The process of acquiring an intensity signal sequence indicating the signal intensity from multiple sensors,
The process of calculating the correlation coefficient between intensity signal sequences and
The process of comparing the correlation coefficient with a predetermined threshold and
As a result of the above comparison, the process of specifying the event position from the signal intensity ratio, which is performed when the value is equal to or higher than the threshold value,
It includes a process for determining multiple events, which is performed when the value is smaller than the threshold value.
The process of determining multiple events is
Let the number of events ξ be 2
Correlation coefficient r (ξ) between I = AS given by ξ event signal strings I, coefficient matrix A of ξ rows and M columns, and intensity signal sequence S of M rows and 1 column. The first process of determining the coefficient matrix A that minimizes
From I = AS, which is given by adding 1 to the number of events ξ, ξ event signal strings I, the coefficient matrix A of ξ rows and M columns, and the intensity signal sequence S of M rows and 1 column, the event signal sequence. The second process of determining the coefficient matrix A that minimizes the correlation coefficient r (ξ) between I and
With the process of comparing r (ξ) and r (ξ-1),
When r (ξ) is r (ξ-1) or less , the second process and the process of comparing r (ξ) and r (ξ-1) are further performed.
In other cases, the person detection method is characterized in that the event position is specified from the coefficient of each row of the coefficient matrix A obtained when the coefficient is ξ-1 or less.
イベント位置が特定された後、特定されたイベントを仮想人に割り当てる
ことを特徴とする請求項1に記載の人検知方法。
The person detection method according to claim 1, wherein after the event position is specified, the specified event is assigned to a virtual person.
複数の静電気測定器と、
前記複数の静電気測定器の出力の高周波成分を透過させるバンドパスフィルタと、
前記バンドパスフィルタを透過した信号から強度信号列を生成するAD変換器と、
信号処理部と
を備え、
前記信号処理部は、
請求項1又は2に記載の人検知方法の各手順を実行させる手段を備える
ことを特徴とする人検知装置。
With multiple static electricity measuring instruments
A bandpass filter that transmits high-frequency components of the outputs of the plurality of static electricity measuring instruments,
An AD converter that generates an intensity signal sequence from a signal that has passed through the bandpass filter,
Equipped with a signal processing unit
The signal processing unit
A person detection device comprising means for executing each procedure of the person detection method according to claim 1 or 2.
前記バンドパスフィルタが5〜1000Hzを透過させる
ことを特徴とする請求項3に記載の人検知装置。
The person detection device according to claim 3, wherein the bandpass filter transmits 5 to 1000 Hz.
低周波用のAD変換器及び信号処理部を備え、
前記バンドパスフィルタが5〜1000Hzを透過させるとともに、5Hz未満の低周波成分を分離して出力し、
前記低周波用のAD変換器は、低周波成分の強度信号列を生成し、
前記低周波用の信号処理部は、前記低周波成分の強度信号列から信号源の総電荷量を取得する
ことを特徴とする請求項3に記載の人検知装置。
Equipped with an AD converter for low frequencies and a signal processing unit
The bandpass filter transmits 5 to 1000 Hz and separates and outputs low frequency components of less than 5 Hz.
The AD converter for low frequencies generates an intensity signal sequence of low frequency components.
The person detection device according to claim 3, wherein the low-frequency signal processing unit acquires the total charge amount of the signal source from the intensity signal sequence of the low-frequency component.
前記静電気測定器は、電界センサ、電位センサ、誘導電流監視センサのいずれか1つ又は組み合わせである
ことを特徴とする請求項3〜5のいずれか一項に記載の人検知装置。
The person detection device according to any one of claims 3 to 5, wherein the static electricity measuring device is any one or a combination of an electric field sensor, a potential sensor, and an induced current monitoring sensor.
JP2016170268A 2016-08-31 2016-08-31 Person detection method and person detection device Active JP6763241B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016170268A JP6763241B2 (en) 2016-08-31 2016-08-31 Person detection method and person detection device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016170268A JP6763241B2 (en) 2016-08-31 2016-08-31 Person detection method and person detection device

Publications (2)

Publication Number Publication Date
JP2018033746A JP2018033746A (en) 2018-03-08
JP6763241B2 true JP6763241B2 (en) 2020-09-30

Family

ID=61564996

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016170268A Active JP6763241B2 (en) 2016-08-31 2016-08-31 Person detection method and person detection device

Country Status (1)

Country Link
JP (1) JP6763241B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109661085A (en) * 2018-12-29 2019-04-19 桂林凯歌信息科技有限公司 Dynamic and static human body sense light and control method based on induction field

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02228584A (en) * 1989-03-02 1990-09-11 Nippon Mining Co Ltd Detecting system of moving object and apparatus therefor
JP4943184B2 (en) * 2007-03-02 2012-05-30 三菱電機株式会社 Movement amount calculation device
JP5179132B2 (en) * 2007-09-27 2013-04-10 パナソニック株式会社 Human body detection device
WO2012129654A1 (en) * 2011-03-31 2012-10-04 Laurentian University Of Sudbury Multi-component electromagnetic prospecting apparatus and method of use thereof
GB201118970D0 (en) * 2011-11-02 2011-12-14 Univ Sussex Target position,movement and tracking system
JP2013186072A (en) * 2012-03-09 2013-09-19 Panasonic Corp Intruder detector

Also Published As

Publication number Publication date
JP2018033746A (en) 2018-03-08

Similar Documents

Publication Publication Date Title
Mastorakis et al. Fall detection system using Kinect’s infrared sensor
CN105426027B (en) The capacitive touch screen of adaptive touch-sensing threshold value with the acutance based on Capacitive data
JP2018526711A5 (en)
JP2015062817A5 (en) Brain activity analyzing device, brain activity analyzing method, discriminator generating device, discriminator generating method, biomarker device and program, health management device and program, and discriminator program
KR20160066927A (en) Apparatus and method for supporting computer aided diagnosis
JP2015527111A5 (en)
Prance et al. Position and movement sensing at metre standoff distances using ambient electric field
EP2764354A1 (en) Walk-through metal detection system
JP6666245B2 (en) Method, system, computer program product and apparatus for selecting operating mode of mobile platform
JP2020511285A (en) Stroke characterization by magnetic resonance imaging (MRI) fused with electroencephalogram (EEG) model
CN105303563B (en) A kind of fall detection method and device
JP6763241B2 (en) Person detection method and person detection device
Kurita et al. Noncontact detection of movements of standing up from and sitting down on a chair using electrostatic induction
JP7266599B2 (en) Devices, systems and methods for sensing patient body movement
Sowmyayani et al. Fall detection in elderly care system based on group of pictures
CN104318248B (en) Action identification method and action recognition device
US20220101709A1 (en) INDOOR OCCUPANCY ESTIMATION, TRAJECTORY TRACKING and EVENT MONITORING AND TRACKING SYSTEM
Wong et al. Enhanced classification of abnormal gait using BSN and depth
Mühlbacher-Karrer et al. Object detection based on electrical capacitance tomography
CN108981759A (en) The multi-dimensional method that system and process are imaged, are monitored or are measured using capacitance sensor
US20210038122A1 (en) Method for detecting body movements of a sleeping person
CN214180021U (en) Intelligent carpet based on electrode array
Mosca et al. Human Walking Behavior detection with a RGB-D Sensors Network for Ambient Assisted Living Applications.
Haramaki et al. A sensor fusion approach for network visualization
Król et al. Using principal component analysis and elastic net in logistic regression to identify the location of objects in EIT

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190515

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200520

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200616

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200721

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200824

R150 Certificate of patent or registration of utility model

Ref document number: 6763241

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150