JP7440131B2 - 生体信号測定用のレーダにおけるレーダ収集信号のランダムノイズの除去方法及びその装置 - Google Patents
生体信号測定用のレーダにおけるレーダ収集信号のランダムノイズの除去方法及びその装置 Download PDFInfo
- Publication number
- JP7440131B2 JP7440131B2 JP2022554166A JP2022554166A JP7440131B2 JP 7440131 B2 JP7440131 B2 JP 7440131B2 JP 2022554166 A JP2022554166 A JP 2022554166A JP 2022554166 A JP2022554166 A JP 2022554166A JP 7440131 B2 JP7440131 B2 JP 7440131B2
- Authority
- JP
- Japan
- Prior art keywords
- signal
- signals
- radar
- index
- random noise
- 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
Links
- 238000000034 method Methods 0.000 title claims description 16
- 230000000737 periodic effect Effects 0.000 claims description 24
- 238000001914 filtration Methods 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 14
- 230000003595 spectral effect Effects 0.000 claims description 14
- 230000010355 oscillation Effects 0.000 claims description 12
- 238000001228 spectrum Methods 0.000 claims description 11
- 238000005070 sampling Methods 0.000 claims description 8
- 238000005259 measurement Methods 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims 1
- 230000000694 effects Effects 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 239000000969 carrier Substances 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000001747 exhibiting effect Effects 0.000 description 1
- 230000001121 heart beat frequency Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000002250 progressing effect Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/288—Coherent receivers
- G01S7/2883—Coherent receivers using FFT processing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
- G01S7/2921—Extracting wanted echo-signals based on data belonging to one radar period
- G01S7/2922—Extracting wanted echo-signals based on data belonging to one radar period by using a controlled threshold
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/0507—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves using microwaves or terahertz waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/52—Discriminating between fixed and moving objects or between objects moving at different speeds
- G01S13/56—Discriminating between fixed and moving objects or between objects moving at different speeds for presence detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
- G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/295—Means for transforming co-ordinates or for evaluating data, e.g. using computers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/415—Identification of targets based on measurements of movement associated with the target
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
- A61B5/0816—Measuring devices for examining respiratory frequency
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Signal Processing (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Physiology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Power Engineering (AREA)
- Radar Systems Or Details Thereof (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Description
本発明は、マイクロ波ドップラレーダを用いて人体の生体信号を取得する間に外部から取り込まれる各種の振動ノイズを効果的に減衰しかつ除去するための方法とその装置に関する。
一般に、センシング技術の一つとしてのレーダは、物体までの正確な距離と観測個所に対する物体の相対速度とを正確に測定することができる。レーダ装置は、概ねマイクロ波の電磁気波を物体に発射して、その物体から反射される電磁気波を受信して作動する。処理された信号は、オペレータやレーダによって制御される周辺装置が使用可能な形態に変換される。レーダを用いて動きや生体信号などを感知する技術の商用化が進んでいく一方である。実際に、レーダが用いられる環境下で生じ得る外部の影響によるレーダ受信信号に含まれている各種のノイズ問題を解決するための方案が提案されている。
また、ドップラレーダ（Ｄｏｐｐｌｅｒ Ｒａｄａｒ）は、電波のドップラ効果を用いたものであり、目標に向かって送信するレーダ電波の周波数と反射する周波数との差分にて移動する目標を検出する。気象レーダや航空機の自律航法装置、軍事用のレーダに用いられる。気象用のものにおいては、雲の内部において起きている風の速度の変化を測定する。自律航法装置においては、地面に達する電波の速度を測定して現在の位置を算出する。軍事用のレーダにおいては、通常、単一パルス信号を用いて地表及び海水面上の反射波中において移動する目標のみを捕捉し、かつ追跡するパルスドップラレーダが主流をなしている。
図１中、従来のマイクロ波レーダを用いて生体信号を測定する方法においては、心拍動、呼吸によって生じる人体の振動や動きに伴い生じるドップラ効果をデータとして取り集めることによりセンシングするのが普通である。
しかしながら、実際に人体の生体信号を測定する過程においては、人体の動きや筋肉の動きまたは周りの環境要因によるさらなるドップラ効果が生じる。したがって、図２中、実際にレーダがセンシングする収集信号には、生体信号データの他に、生体信号が存在する時間的、周波数的な同一の領域にランダムなノイズが一緒に載ってしまう。このようなランダムノイズは、周期性に乏しい故に、経時的に種々に変化しながら現れる特性があるのに対し、生体信号は、所定の周期的なパターンを保持するという特性がある。そのため、人体の心拍スペクトル領域に混ざったランダムノイズとの差分を演算して区別できる取り組みが行われなければならない。
本発明と関わる先行技術であって、特許文献１は、心拍を測定しようとするターゲットから反射されたレーダパルスが重なり合った単一のフレームを予め定められた受信時間に従って累積させてフレームセット（Ｆｒａｍｅ Ｓｅｔ）を生成するステップと、前記生成されたフレームセットに含まれている単一のフレームのそれぞれに対して最大のピークを示す第１のサンプラインデックス及び前記第１のサンプラインデックスのうち最も多い最大のピークに対応するサンプラインデックスである第２のサンプラインデックスを用いて、前記ターゲットの動きによって前記最大のピークの遷移区間が生じた単一のフレームを除去するステップと、前記ターゲットの動きによって前記最大のピークの遷移区間が生じた単一のフレームが除去されて空き空間が含まれているフレームセットに対して、不完全データに対する周期的なパターンを探知してピーク周波数を識別するアルゴリズムを適用することにより、前記ターゲットの心拍周波数を決めるステップと、を含み、前記ターゲットから反射されたレーダパルスは、前記ターゲットの動きが反映されたレーダ信号であり、前記遷移区間が生じた部分を除去するステップは、前記フレームセットに含まれている単一のフレームのそれぞれに対して最大のピークを示す第１のサンプラインデックスを取り出すステップと、前記最大のピークを示す第１のサンプラインデックスのうち、最も多い最大のピークに対応するサンプラインデックスを動きプロファイルを生成するための基準を示す第２のサンプラインデックスとして決めるステップと、前記第１のサンプラインデックスと第２のサンプラインデックスとの差分を用いて、前記ターゲットの動きに基づく動きプロファイルを生成するステップと、前記動きプロファイルを用いて、前記フレームセットに含まれている単一のフレームを整列するステップと、を含むターゲットの生体情報の決定方法が開示されている。
また、特許文献２には、基底帯域信号を送信した後、送信した基底帯域信号が物体に反射されて戻ってくる基底帯域信号に基づいてドップラを測定するドップラレーダにおいて、前記ドップラレーダは、前記基底帯域信号として、多重キャリアを含む基底帯域信号を用い、前記ドップラレーダは、受信した基底帯域信号をフーリエ変換を用いてＮ個の周波数信号に変換し、変換されたＮ個の周波数のうち、任意のＭ個の周波数信号を選び、選ばれたＭ個の周波数信号のそれぞれにキャリブレーション定数を乗算して出力されるＭ個の周波数信号のそれぞれの位相値を取り出し、取り出されたＭ個の位相値のそれぞれに距離変換定数を乗算して出力されるＭ個の距離情報のそれぞれにフィルタを適用して出力される値を加算してドップラ測定値を出力する信号処理部を備える多重キャリアドップラレーダが開示されている。
本発明は、上記の問題を解決するためのものであって、ランダムノイズが比較的に周期性を有する生体信号領域に混ざっているとき、ランダムノイズを効果的に除去することを目的としている。
また、本発明は、マイクロ波レーダを用いて非接触式にて生体信号を測定する正確度を改善することを他の目的としている。
本発明は、上記の目的を達成するために、生体信号測定用のレーダにおけるレーダ収集信号のランダムノイズの除去方法において、（ａ）ＩＦ信号生成器の電圧制御発振部（ＶＣＯ）において外部から印加された電圧Ｖ（ｔ）にて所定の周期の発振周波数を生成して、単一の送信アンテナを介して動的ターゲットの主な部位に送出し、前記動的ターゲットから反射された信号と動的ターゲットの周りにおいて生じた各種の信号をｎ個の受信アンテナにおいて受信するステップと、（ｂ）前記ＩＦ信号生成器は、受信されたｎ個の信号からそれぞれドップラＩＦ信号を生成するステップと、（ｃ）アナログ／デジタル変換器は、前記ＩＦ信号生成器から入力された複数のドップラＩＦ信号をそれぞれのデジタルデータに変換するステップと、（ｄ）前記アナログ／デジタル変換器において単位時間の間に収集されたデジタル信号を高速フーリエ変換器においてサンプリング時間間隔にてシンボルを有するデータ列として構成し、このデータ列を複数のインデックス（Ｉｎｄｅｘ）の周波数成分シンボル列（ｓｙｍｂｏｌ ｓｅｔ）に切り換えるステップと、（ｅ）前記信号演算器の演算部において周波数成分シンボル列に切り換えられたインデックスシンボルを各インデックス列に加算して、ｎ個の受信アンテナの数で割った値を演算するステップと、（ｆ）前記演算部において演算された値に対して共通して生じる周期的な信号と共通的ではない非周期的な信号とのスペクトル成分の偏差を所定の基準のしきい（Ｔｈｒｅｓｈｏｌｄ）値に応じて分類して周期的な信号のみをフィルタリングするステップと、を含んでなる、生体信号測定用のレーダにおけるレーダ収集信号のランダムノイズの除去方法を提供するところに特徴がある。
また、本発明は、生体信号測定用のレーダにおけるレーダ収集信号のランダムノイズの除去装置において、電圧制御発振部（Ｖｏｌｔａｇｅ Ｃｏｎｔｒｏｌｌｅｄ Ｏｓｃｉｌｌａｔｏｒ；ＶＣＯ）において外部から印加された電圧Ｖ（ｔ）にて所定の周期の発振周波数を生成して単一の送信アンテナＴｘから送出し、所定の距離の動的ターゲットを有するレーダ信号収集領域から反射された反射波を複数の受信アンテナＲｘ－１～Ｒｘ－ｎにおいてそれぞれ受信し、受信されたｎ個の信号からドップラＩＦ信号を生成するＩＦ信号生成器と、前記ＩＦ信号生成器において生成されたｎ個のアナログ信号をデジタル信号に変換するアナログ／デジタル変換器と、前記アナログ／デジタル変換器において単位時間の間に収集されたデジタル信号を高速フーリエ変換器においてサンプリング時間間隔にてシンボルを有するデータ列として構成し、構成されたデータ列をインデックス（Ｉｎｄｅｘ）の周波数成分シンボル列（ｓｙｍｂｏｌ ｓｅｔ）に切り換えた後、演算部においてインデックスシンボルを各インデックス列に加算して、ｎ個の受信アンテナＲｘの数で割った値を求める信号演算器と、前記信号演算器において求めた値に対して共通して生じる周期的な信号と共通的ではない非周期的な信号とのスペクトル成分の偏差を所定の基準のしきい（Ｔｈｒｅｓｈｏｌｄ）値に応じて分類してフィルタリングするデジタルフィルタリング信号処理器と、を備える、生体信号測定用のレーダにおけるレーダ収集信号のランダムノイズの除去装置を提供するところに特徴がある。
さらに、本発明において、前記受信アンテナＲｘ－１～Ｒｘ－ｎは、動的ターゲットとそれ以外の領域から同時に信号を受信できる位置及び方向に並べられてもよい。
さらにまた、本発明において、前記デジタルフィルタリング信号処理器は、基準しきい値にてフィルタリングするために所定の帯域のバンドパスフィルタ（ＢＰＦ）や高域パスフィルタ（ＨＰＦ）または低域パスフィルタ（ＬＰＦ）のいずれか一つが適用されてもよい。
本発明によれば、単一の送信部と複数の受信部とからなる非接触式生体信号測定用のドップラレーダを用いて、生体信号領域に含まれているランダムノイズが比較的に周期性を有する生体信号領域に混ざっているとき、ランダムノイズを効果的に減衰しかつ除去することにより、生体信号をセンシングする正確度を向上させることができるというメリットがある。
以下、添付図面に基づいて、本発明に係る生体信号測定用のレーダにおけるレーダ収集信号のランダムノイズの除去装置に関する実施形態について詳述する。
図３及び図４中、マイクロ波（Ｍｉｃｒｏｗａｖｅ）ドップラレーダ１の中間周波数（ＩＦ：Ｉｎｔｅｒｍｅｄｉａｔｅ Ｆｒｅｑｕｅｎｃｙ）信号生成器１０は、電圧制御発振部（Ｖｏｌｔａｇｅ Ｃｏｎｔｒｏｌｌｅｄ Ｏｓｃｉｌｌａｔｏｒ；ＶＣＯ）１１において生じた発振周波数と受信された周波数との差分から生じるＩＦ信号を生成するものである。ＩＦ信号生成器１０には、電圧制御発振部１１から出力された所定の周期の発振周波数を外部に送出する送信アンテナが組み込まれる。そして、ＩＦ信号生成器１０には、送信アンテナから送出された信号が所定の距離の動的ターゲット、例えば、人体から反射された信号と人体の周りのレーダ信号収集領域から信号を受信する受信アンテナが組み込まれる。ここで、送信アンテナＴｘは、単一個から構成され、受信アンテナＲｘ－１～Ｒｘ－ｎは、複数から構成される。受信アンテナＲｘ－１～Ｒｘ－ｎは、動的ターゲットとそれ以外の領域から同時に信号を受信できるように適切な位置及び方向に並べられることが好ましい。
そして、ＩＦ信号生成器１０には、各受信アンテナＲｘ－１～Ｒｘ－ｎのそれぞれにミキサＭ１～Ｍｎが結合される。ミキサＭ１～Ｍｎは、電圧制御発振部１１から出力された発振周波数と各受信アンテナＲｘ－１～Ｒｘ－ｎにおいて受信されたｎ個の信号から生じる差分にてドップラＩＦ信号を生成する。
ＩＦ信号生成器１０は、信号が送受信可能なドップラレーダトランシーバ（Ｔｒａｎｃｅｉｖｅｒ）と同じ構造であり、かつ、単一の送信アンテナＴｘと複数の受信アンテナＲｘ－１～Ｒｘ－ｎが配備される。したがって、ＩＦ信号生成器１０は、電圧制御発振部１１において外部から印加された電圧Ｖ（ｔ）にて所定の周期の発振周波数を生成して単一の送信アンテナＴｘから送出し、所定の距離の動的ターゲットを有するレーダ信号収集領域から反射された反射波を複数の受信アンテナＲｘ－１～Ｒｘ－ｎにおいてそれぞれ受信し、受信されたｎ個の信号からドップラＩＦ信号を生成する。
アナログ／デジタル変換器（Ａｎａｌｏｇ／Ｄｉｇｉｔａｌ Ｃｏｎｖｅｒｔｅｒ）２０は、入力されたアナログ信号をデジタル信号に変換して出力するものであり、ＩＦ信号生成器１０において生成されたｎ個のアナログ信号をデジタル信号にそれぞれ変換して出力する。このとき、デジタル信号を周波数シンボルに切り換えるために、高速フーリエ変換を行わなければならないため、単位時間の間に収集されたデジタル信号のデータを用いてサンプリング時間間隔にてデータ列を構成する。
信号演算器３０は、高速フーリエ変換器（Ｆａｓｔ Ｆｏｕｒｉｅｒ Ｔｒａｎｓｆｏｒｍｅｒ；ＦＦＴ）と演算部３２を備える。高速フーリエ変換器３１は、アナログ／デジタル変換器２０において単位時間の間に収集されたデジタル信号のデータ列を周波数成分のシンボルに変換し、これは、ｎ個のスペクトルインデックスデータ列として構成する。そして、演算部３２は、構成されたスペクトルインデックスシンボルを各インデックス列に加算して、ｎ個の受信アンテナＲｘ－１～Ｒｘ－ｎの数で割った値を求める。
デジタルフィルタリング信号処理器４０は、信号演算器３０において単位時間の間に求められた値に対して共通して生じる周期的な信号と共通的ではない非周期的な信号とのスペクトル成分の偏差を所定の基準のしきい（Ｔｈｒｅｓｈｏｌｄ）値に応じて分類してフィルタリングするものである。すなわち、デジタルフィルタリング信号処理器４０は、単位時間の間に共通して生じる周期性を有する信号の演算結果は所定の大きさの値に保持されるものの、共通的ではない非周期性を有する信号は格段に低い値に下がってしまうため、周期的な信号と非周期的な信号とのスペクトル成分の偏差を大きくした後、適切な基準しきい値に応じて分類してフィルタリングして信号を処理する。基準しきい値にてフィルタリングするためのフィルタは、所定の帯域のバンドパスフィルタ（ＢＰＦ）や高域パスフィルタ（ＨＰＦ）または低域パスフィルタ（ＬＰＦ）などが適宜に選択できる筈である。
このような構成を有する本発明に係る生体信号測定用のレーダにおけるレーダ収集信号のランダムノイズの除去方法について説明する。
まず、図５中、ＩＦ信号生成器１０の電圧制御発振部１１において外部から印加された電圧Ｖ（ｔ）にて所定の周期の発振周波数を生成して、送信アンテナＴｘを介して動的ターゲットの主な部位に送出される（Ｓ１）。さらに、送信アンテナＴｘから送出される信号は、動的ターゲットである生体信号を最も上手く収集できる人体の部位に向かって送出することが好ましい。
動的ターゲットから反射された信号と動的ターゲットの周りの各種の信号は、ｎ個の受信アンテナＲｘ－１～Ｒｘ－ｎにおいて受信される（Ｓ２）。このとき、ｎ個の受信アンテナＲｘ－１～Ｒｘ－ｎは、生体信号を最も多く含む信号を受け取る受信アンテナＲｘ－１と、生体信号とともに、人体の他の動きや周りにおいて生じた効果音などが含まれている信号を受け取る複数の受信アンテナＲｘ－２～Ｒｘ－ｎと、を備える。
複数の受信アンテナＲｘ－１～Ｒｘ－ｎにおいて受信されたｎ個の信号は、ＩＦ信号生成器１０においてそれぞれドップラＩＦ信号として生成される（Ｓ３）。そして、アナログ／デジタル変換器２０は、ＩＦ信号生成器１０において生成された後に入力された複数のドップラＩＦ信号をそれぞれのデジタルデータに変換する（Ｓ４）。デジタルデータは、複数の受信アンテナＲｘ－１～Ｒｘ－ｎにおいて信号を受信する間に時間的に持続的に繰り返しが行われるため、リアルタイムにて単位時間別にデータを連続して収集することができる筈である。
アナログ／デジタル変換器２０において単位時間の間に収集されたデジタル信号を信号演算器３０の高速フーリエ変換器３１においてサンプリング時間間隔にてシンボルを有するデータ列として生成し、このデータ列は、複数のインデックス（Ｉｎｄｅｘ）の周波数成分シンボル列（ｓｙｍｂｏｌ ｓｅｔ）に切り換えられる（Ｓ５）。
一方、ドップラＩＦ信号のスペクトル変換は、次の過程により行われる。すなわち、周波数成分のスペクトル変換においては、所定の時間の間に収集された時間領域の信号を、周期性を見出して周波数成分に変換する。したがって、高速フーリエ変換のために所定の時間の間に収集された時間領域のデジタルサンプリングデータを必要とする。さらに、サンプリング時間間隔のシンボルは、時間領域においてデジタル信号にサンプリングされたオリジナルの信号成分のシンボルである。例えば、３０秒間の時間領域の信号データを取り集めてフーリエ変換を行うと、３０秒の間に周期性が現れたすべての信号の周波数成分のスペクトルが表わされる。このとき、スペクトルグラフ上の横軸が周波数軸となり、これがスペクトルインデックスである。
複数のインデックスは、図６の表において、複数の受信信号は高速フーリエ変換器３１を経由しながらＩｍの周波数成分シンボル列に切り換えられる。そして、演算部３２は、複数の受信アンテナＲｘ－１～Ｒｘ－ｎからそれぞれ取得されたインデックスシンボルを各インデックスＩ（１），Ｉ（２），Ｉ（３），．．．Ｉ（ｍ）列に加算して、受信アンテナＲｘ－１～Ｒｘ－ｎの数であるｎで割った値を演算する（Ｓ６）。
ここで、ｎ個の受信アンテナＲｘ－１～Ｒｘ－ｎ信号は、ｎ個のスペクトルインデックスデータとして表わされ、演算は、各アンテナの同じスペクトル成分であるインデックス列へと進む。なお、Ｉ（１），Ｉ（２），…Ｉ（ｍ）と表わされたものは、実際の周波数（Ｈｚ）成分であり、このシンボルデータは、スペクトルの強さを示す。
さらに、図７のスペクトルグラフにおいて、共通して生じる周期的な信号は、所定の値を保持するものの、共通的ではない非周期的な信号は格段に低い値に下がってしまうため、周期的な信号と非周期的な信号とのスペクトル偏差を作り出すことが可能になる。すなわち、周波数ドメインインデックスＩ（ｍ）列のＩ（１）＝ａ１＋ｂ１＋ｃ１…（ｎ）１、Ｉ（２）＝ａ２＋ｂ２＋ｃ２…（ｎ）２、Ｉ（３）＝ａ３＋ｂ３＋ｃ３…（ｎ）３、Ｉ（ｍ）＝ａ（ｍ）＋ｂ（ｍ）＋ｃ（ｍ）…（ｎ）（ｍ）などといったように計算するが、受信アンテナＲｘ－１～Ｒｘ－ｎの数がｎ個であれば、各インデックスの演算結果は、Ｉ（１）／ｎ，Ｉ（２）／ｎ，Ｉ（３）／ｎ，…Ｉ（ｍ）／ｎといったように計算される。
そして、デジタルフィルタリング信号処理器４０は、演算部３２において演算された値に対して共通して生じる周期的な信号と共通的ではない非周期的な信号とのスペクトル成分の偏差を所定の基準のしきい（Ｔｈｒｅｓｈｏｌｄ）値に応じて分類して周期的な信号のみをフィルタリングする（Ｓ７）。このとき、しきい値を基準として分類すれば、基準点以上の高い値のみを残すことができ、これは、通常、時間的に周期的な成分を有している生体信号である可能性が高くなる。
このような方法を用いて、各受信アンテナＲｘ－１～Ｒｘ－ｎにおいて受信されたデータのうち、共通して受信された部分を取り、共通していない部分をランダムノイズとみなしてフィルタリングして除去することができる。
このように、単一の送受信部からなる生体信号測定用のドップラレーダを用いて生体信号を測定する従来の方法に比べて、本発明は、単一の送信部と複数の受信部とからなる生体信号測定用のドップラレーダを用いて、生体信号領域に含まれているランダムノイズを効果的に除去することにより、非接触方式において生体信号をセンシングする正確度を改善することができるというメリットがある。
以上の説明において、本発明については、特定の実施形態と関連付けて図示及び説明したが、特許請求の範囲によって示された発明の思想及び領域から逸脱しない範囲内において種々の改造及び変化を行うことが可能であるということは、この技術分野における通常の知識を有する者であれば誰でも理解できる筈である。
Claims (4)
- 生体信号測定用のレーダにおけるレーダ収集信号のランダムノイズの除去方法であって、
（ａ）ＩＦ信号生成器の電圧制御発振部（ＶＣＯ）において外部から印加された電圧（Ｖ（ｔ））にて所定の周期の発振周波数を有する信号を連続的に生成して、単一の送信アンテナを介して動的ターゲットの主な部位に送出し、前記動的ターゲットから反射された信号と前記動的ターゲットの周りにおいて生じた各種の信号をｎ個（ｎは２以上の整数）の受信アンテナにおいて受信するステップと、
（ｂ）前記ＩＦ信号生成器は、前記発信周波数を有する信号と受信されたｎ個の信号との差分からそれぞれドップラＩＦ信号を生成するステップと、
（ｃ）アナログ／デジタル変換器が、前記ＩＦ信号生成器から入力された複数のドップラＩＦ信号をそれぞれのデジタルデータに変換するステップと、
（ｄ）前記アナログ／デジタル変換器（２０）において単位時間の間に収集されたデジタル信号を高速フーリエ変換器においてサンプリング時間間隔にてシンボルを有するデータ列として構成し、前記データ列を複数のインデックス（Ｉｎｄｅｘ）の周波数成分シンボル列（ｓｙｍｂｏｌ ｓｅｔ）に切り換えるステップと、
（ｅ）信号演算器の演算部において前記周波数成分シンボル列に切り換えられたインデックスシンボルを各インデックス列に加算して、ｎ個の受信アンテナの数で割った値を演算するステップと、
（ｆ）前記演算部において演算された値に対して共通して生じる周期的な信号と共通的ではない非周期的な信号とのスペクトル成分の偏差を所定の基準のしきい（Ｔｈｒｅｓｈｏｌｄ）値に応じて分類して周期的な信号のみをフィルタリングするステップと、
を含み、
前記周波数成分シンボル列のシンボルデータはスペクトルの強さを示す、生体信号測定用のレーダにおけるレーダ収集信号のランダムノイズの除去方法。 - 生体信号測定用のレーダにおけるレーダ収集信号のランダムノイズの除去装置であって、
電圧制御発振部（Ｖｏｌｔａｇｅ Ｃｏｎｔｒｏｌｌｅｄ Ｏｓｃｉｌｌａｔｏｒ；ＶＣＯ）において外部から印加された電圧（Ｖ（ｔ））にて所定の周期の発振周波数を有する信号を連続的に生成して単一の送信アンテナ（Ｔｘ）から送出し、所定の距離の動的ターゲットを有するレーダ信号収集領域から反射された反射波を複数の受信アンテナ（Ｒｘ－１～Ｒｘ－ｎ）（ｎは２以上の整数）においてそれぞれ受信し、前記発信周波数を有する信号と受信されたｎ個の信号との差分からそれぞれドップラＩＦ信号を生成するＩＦ信号生成器と、
前記ＩＦ信号生成器において生成されたｎ個のアナログ信号をデジタル信号に変換するアナログ／デジタル変換器と、
前記アナログ／デジタル変換器において単位時間の間に収集されたデジタル信号を高速フーリエ変換器においてサンプリング時間間隔にてシンボルを有するデータ列として構成し、構成された前記データ列をインデックス（Ｉｎｄｅｘ）の周波数成分シンボル列（ｓｙｍｂｏｌ ｓｅｔ）に切り換えた後、演算部においてインデックスシンボルを各インデックス列に加算して、ｎ個の受信アンテナ（Ｒｘ）の数で割った値を求める信号演算器と、
前記信号演算器において求めた値に対して共通して生じる周期的な信号と共通的ではない非周期的な信号とのスペクトル成分の偏差を所定の基準のしきい（Ｔｈｒｅｓｈｏｌｄ）値に応じて分類してフィルタリングするデジタルフィルタリング信号処理器と、
を備え、
前記周波数成分シンボル列のシンボルデータはスペクトルの強さを示す、生体信号測定用のレーダにおけるレーダ収集信号のランダムノイズの除去装置。 - 前記受信アンテナ（Ｒｘ－１～Ｒｘ－ｎ）は、動的ターゲットとそれ以外の領域から同時に信号を受信できる位置及び方向に並べられた、請求項２に記載の生体信号測定用のレーダにおけるレーダ収集信号のランダムノイズの除去装置。
- 前記デジタルフィルタリング信号処理器は、基準しきい値にてフィルタリングするために所定の帯域のバンドパスフィルタ（ＢＰＦ）や高域パスフィルタ（ＨＰＦ）または低域パスフィルタ（ＬＰＦ）のいずれか一つが適用される、請求項２に記載の生体信号測定用のレーダにおけるレーダ収集信号のランダムノイズの除去装置。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| KR1020190143162A KR102122758B1 (ko) | 2019-11-11 | 2019-11-11 | 생체신호 측정용 레이더에서 레이더 수집신호의 랜덤노이즈 제거방법 및 그 장치 |
| KR10-2019-0143162 | 2019-11-11 | ||
| PCT/KR2019/015255 WO2021095893A1 (ko) | 2019-11-11 | 2019-11-11 | 생체신호 측정용 레이더에서 레이더 수집신호의 랜덤노이즈 제거방법 및 그 장치 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2023502177A JP2023502177A (ja) | 2023-01-20 |
| JP7440131B2 true JP7440131B2 (ja) | 2024-02-28 |
Family
ID=71136373
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2022554166A Active JP7440131B2 (ja) | 2019-11-11 | 2019-11-11 | 生体信号測定用のレーダにおけるレーダ収集信号のランダムノイズの除去方法及びその装置 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US20220381877A1 (ja) |
| JP (1) | JP7440131B2 (ja) |
| KR (1) | KR102122758B1 (ja) |
| CN (1) | CN114667094A (ja) |
| WO (1) | WO2021095893A1 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111812629A (zh) * | 2020-07-21 | 2020-10-23 | 成都阶跃时进科技有限公司 | 一种活体存在检测方法、装置、设备及可读存储介质 |
| KR102381262B1 (ko) * | 2021-12-17 | 2022-04-01 | (주)디지탈엣지 | 레이더 시스템 및 그에서 수행되는 생체 신호 탐지방법 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080077015A1 (en) | 2006-05-17 | 2008-03-27 | Olga Boric-Lubecke | Determining presence and/or physiological motion of one or more subjects with multiple receiver Doppler radar systems |
| JP2012130391A (ja) | 2010-12-20 | 2012-07-12 | Oki Electric Industry Co Ltd | 心拍数検出装置、心拍数検出方法、およびプログラム |
| JP2016057168A (ja) | 2014-09-10 | 2016-04-21 | 三菱電機株式会社 | レーダ装置 |
| JP2016156754A (ja) | 2015-02-25 | 2016-09-01 | 古河電気工業株式会社 | レーダ装置およびレーダ装置の対象物検知方法 |
| JP2016168177A (ja) | 2015-03-12 | 2016-09-23 | 国立大学法人 筑波大学 | 生体情報検出装置および背もたれ付シート |
| US20160336989A1 (en) | 2015-05-14 | 2016-11-17 | University Of Florida Research Foundation, Inc. | Low if architectures for noncontact vital sign detection |
| JP2017129558A (ja) | 2016-01-15 | 2017-07-27 | パナソニックＩｐマネジメント株式会社 | 推定装置および推定方法 |
| CN111812629A (zh) | 2020-07-21 | 2020-10-23 | 成都阶跃时进科技有限公司 | 一种活体存在检测方法、装置、设备及可读存储介质 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009501044A (ja) * | 2005-07-15 | 2009-01-15 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 心臓の活動の検出用装置 |
| CA2720871A1 (en) * | 2008-04-03 | 2009-10-08 | Kai Medical, Inc. | Non-contact physiologic motion sensors and methods for use |
| JP5527277B2 (ja) * | 2011-05-11 | 2014-06-18 | コニカミノルタ株式会社 | 信号処理装置、信号処理方法および生体情報測定装置 |
| CN106999105B (zh) * | 2014-11-28 | 2019-10-18 | 夏普株式会社 | 高频装置 |
| JP6281912B2 (ja) * | 2015-03-23 | 2018-02-21 | 国立大学法人九州工業大学 | 生体信号検出装置 |
| KR102362085B1 (ko) | 2016-07-22 | 2022-02-15 | 한국전자통신연구원 | 다중 캐리어 도플러 레이더 |
| KR101948386B1 (ko) | 2016-10-04 | 2019-02-14 | 재단법인대구경북과학기술원 | 임펄스 레이더를 이용한 타겟의 생체 정보 결정 장치 및 방법 |
-
2019
- 2019-11-11 WO PCT/KR2019/015255 patent/WO2021095893A1/ko active Application Filing
- 2019-11-11 CN CN201980102160.8A patent/CN114667094A/zh active Pending
- 2019-11-11 US US17/775,512 patent/US20220381877A1/en active Pending
- 2019-11-11 JP JP2022554166A patent/JP7440131B2/ja active Active
- 2019-11-11 KR KR1020190143162A patent/KR102122758B1/ko active IP Right Grant
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080077015A1 (en) | 2006-05-17 | 2008-03-27 | Olga Boric-Lubecke | Determining presence and/or physiological motion of one or more subjects with multiple receiver Doppler radar systems |
| JP2012130391A (ja) | 2010-12-20 | 2012-07-12 | Oki Electric Industry Co Ltd | 心拍数検出装置、心拍数検出方法、およびプログラム |
| JP2016057168A (ja) | 2014-09-10 | 2016-04-21 | 三菱電機株式会社 | レーダ装置 |
| JP2016156754A (ja) | 2015-02-25 | 2016-09-01 | 古河電気工業株式会社 | レーダ装置およびレーダ装置の対象物検知方法 |
| JP2016168177A (ja) | 2015-03-12 | 2016-09-23 | 国立大学法人 筑波大学 | 生体情報検出装置および背もたれ付シート |
| US20160336989A1 (en) | 2015-05-14 | 2016-11-17 | University Of Florida Research Foundation, Inc. | Low if architectures for noncontact vital sign detection |
| JP2017129558A (ja) | 2016-01-15 | 2017-07-27 | パナソニックＩｐマネジメント株式会社 | 推定装置および推定方法 |
| CN111812629A (zh) | 2020-07-21 | 2020-10-23 | 成都阶跃时进科技有限公司 | 一种活体存在检测方法、装置、设备及可读存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2023502177A (ja) | 2023-01-20 |
| KR102122758B1 (ko) | 2020-06-26 |
| CN114667094A (zh) | 2022-06-24 |
| US20220381877A1 (en) | 2022-12-01 |
| WO2021095893A1 (ko) | 2021-05-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112754441B (zh) | 一种基于毫米波的非接触式心跳检测方法 | |
| US9993175B2 (en) | Pulsed ultra-wideband sensor and the method thereof | |
| US8979760B2 (en) | Method and apparatus for automatic optimization of doppler imaging parameters | |
| Sekine et al. | Non-contact heart rate detection using periodic variation in Doppler frequency | |
| JP7440131B2 (ja) | 生体信号測定用のレーダにおけるレーダ収集信号のランダムノイズの除去方法及びその装置 | |
| JP5980587B2 (ja) | レーダ装置及び反射信号処理方法 | |
| JP5333427B2 (ja) | 心拍数検出装置、心拍数検出方法、およびプログラム | |
| WO2013118121A1 (en) | A microwave contactless heart rate sensor | |
| US8121222B2 (en) | Systems and methods for construction of time-frequency surfaces and detection of signals | |
| US11311202B2 (en) | Robust real-time heart rate monitoring method based on heartbeat harmonics using small-scale radar | |
| CN112336322A (zh) | 一种无接触式呼吸或心跳检测方法 | |
| KR20210001575A (ko) | 잡음 환경 내 실시간 심장 박동수 탐지 장치 및 그 방법 | |
| US20230136937A1 (en) | Data pre-processing method and exercise vital sign detection radar | |
| Liu et al. | Vital sign extraction in the presence of radar mutual interference | |
| CN113854990A (zh) | 一种心跳检测方法及装置 | |
| CN117838083A (zh) | 一种基于毫米波雷达的体征快速精确检测方法 | |
| JP2022185752A (ja) | 生体情報検出システム、プログラム、及び、生体情報検出方法 | |
| KR20180137690A (ko) | 레이더 신호의 시계열 데이터 분석을 이용한 생체 정보 결정 장치 및 방법 | |
| KR101916591B1 (ko) | 레이더 신호의 주성분 분석을 이용한 생체 정보 결정 장치 및 방법 | |
| KR102381262B1 (ko) | 레이더 시스템 및 그에서 수행되는 생체 신호 탐지방법 | |
| WO2021161551A1 (ja) | レーダ装置 | |
| KR20150128415A (ko) | 고주파 변조 연속파 기반 침입 감지 장치 및 방법 | |
| KR101827509B1 (ko) | 클루터 신호를 이용한 침입탐지와 오탐지 저감 시스템 및 이를 이용한 침입 탐지 및 오탐지 저감 방법 | |
| JP2017191074A (ja) | Ｃｆａｒ処理装置及びレーダ装置 | |
| KR102444685B1 (ko) | 크기와 위상 사이의 코히어런시에 기반한 활력 징후 측정 거리 결정 장치 및 방법 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20220713 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20220714 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20230628 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20230711 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20231011 |
|
| 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: 20240109 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20240207 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7440131 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |