JP7046277B1 - 物体透視装置、制御回路、記憶媒体、および物体透視方法 - Google Patents

物体透視装置、制御回路、記憶媒体、および物体透視方法 Download PDF

Info

Publication number
JP7046277B1
JP7046277B1 JP2021538716A JP2021538716A JP7046277B1 JP 7046277 B1 JP7046277 B1 JP 7046277B1 JP 2021538716 A JP2021538716 A JP 2021538716A JP 2021538716 A JP2021538716 A JP 2021538716A JP 7046277 B1 JP7046277 B1 JP 7046277B1
Authority
JP
Japan
Prior art keywords
phase
frequency
composite image
data
fluoroscope
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
JP2021538716A
Other languages
English (en)
Other versions
JPWO2022162869A1 (ja
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.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Application granted granted Critical
Publication of JP7046277B1 publication Critical patent/JP7046277B1/ja
Publication of JPWO2022162869A1 publication Critical patent/JPWO2022162869A1/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • 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/06—Systems determining position data of a target
    • G01S13/42—Simultaneous measurement of distance and other co-ordinates
    • G01S13/426—Scanning radar, e.g. 3D radar
    • 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/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904—SAR modes
    • G01S13/9088—Circular SAR [CSAR, C-SAR]
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N22/00—Investigating or analysing materials by the use of microwaves or radio waves, i.e. electromagnetic waves with a wavelength of one millimetre or more
    • 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/06—Systems determining position data of a target
    • G01S13/08—Systems for measuring distance only
    • G01S13/32—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated
    • G01S13/36—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated with phase comparison between the received signal and the contemporaneously transmitted signal
    • G01S13/38—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated with phase comparison between the received signal and the contemporaneously transmitted signal wherein more than one modulation frequency is used
    • 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/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9004—SAR image acquisition techniques
    • G01S13/9005—SAR image acquisition techniques with optical processing of the SAR signals
    • 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/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021—SAR image post-processing techniques
    • 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/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904—SAR modes
    • G01S13/9064—Inverse SAR [ISAR]
    • 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/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904—SAR modes
    • G01S13/9089—SAR having an irregular aperture
    • 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/024—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using polarisation effects
    • 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
    • G01S7/298—Scan converters
    • 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/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/887—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Electromagnetism (AREA)
  • Signal Processing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Health & Medical Sciences (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

測定エリア内に配置された被測定対象である電波散乱体（３０）に対して照射された複数の周波数を含んだ電波の反射波（Ｗ）に基づいて被測定対象である電波散乱体（３０）の画像を生成する物体透視装置（１０Ａ）であって、反射波（Ｗ）に基づいたイメージングを行うことによって得られた複数の画像に対し、複数の画像の画素毎に複素加算を行うことで複数の画像を合成した位相合成画像を生成する位相合成画像生成部（２８）を備える。

Description

本開示は、物体を透視する物体透視装置、制御回路、記憶媒体、および物体透視方法に関する。
被測定対象である物体にミリ波、テラヘルツ波などを照射し、物体からの反射波を用いて物体を透視する装置の１つに、ＭＦ(Matched Filter、マッチドフィルタ)を用いて物体のイメージングを行う物体透視装置がある。この物体透視装置は、測定用の送受信アンテナを決められたパスに従って移動させながら物体に電波を照射し、物体の散乱点からの反射波を受信し、反射波の時系列データを記録する。
物体透視装置は、種々の散乱点から受信した種々の波形を有した反射波を重ね合わせることで受信信号の波形である受信波形を生成する。この時、物体透視装置は、測定エリアと送受信アンテナとの位置関係が既知であるならば、測定エリア上の任意の座標にある散乱点からの反射波の波形を推定することが可能である。物体透視装置は、ＭＦを用いたイメージングを行う場合、測定エリア内に網羅的に観測点を設けておき、各々の観測点に対して反射波の推定波形を生成し、推定波形と受信波形との相関をとることによって、各観測点からの反射強度をマッピングする。この場合において、推定波形と受信波形との相関値は、位相情報を持つ相関ベクトルとして得られる。すなわち、相関ベクトルは、推定波形と受信波形との相関を示すとともに位相の情報を有している。物体透視装置は、観測点と散乱点とが近い場合には高い相関値を得ることができ、観測点の近くに散乱点が存在しない場合には低い相関値を得ることができる。
物体透視装置は、ＭＦを用いたイメージングを行う際に、複数の周波数の電波を用いて同じ測定を繰り返すことによって、送受信アンテナの配置、および測定エリア周辺の散乱物に起因するサイドローブを低減できる。複数の周波数の電波で測定した画像を合成する方法としては、画素毎に相関ベクトルの大きさを積算する電力合成方式と、画素毎に位相情報を考慮して複素加算を行う位相合成方式とがある。
位相合成方式は、電力合成方式と比較して高いサイドローブ低減効果が得られるものの、異なる周波数の電波で測定した同一観測点の相関ベクトルの位相を正確に一致させなければ複素加算時に相関ベクトルが打ち消されてしまう。このため、位相合成方式には、観測点と送受信アンテナとの位置関係の情報、または測定系の周波数特性の影響を考慮した信号処理が必要となる。
位相合成を実現する方法としては、ある画素を異なる周波数の電波で測定した際に発生する位相オフセットを何らかの手段により推定ないしは事前に測定し、測定結果から差し引くことで位相を補正する方法がある。特許文献１に記載のミリ波画像処理装置は、アンテナを用いて受信した信号の信号処理を行う信号処理装置と、信号処理装置の外部に配置された校正信号発生装置とを備えている。校正信号発生装置は、信号処理装置に適用する位相オフセット情報を取得し、信号処理装置は、位相オフセット情報を用いてアンテナが受信する電波の位相を補正している。
特開２００７－２５６１７１号公報
しかしながら、上記特許文献１の技術では、校正信号発生装置が必要となるので、装置構成が複雑になるという問題があった。
本開示は、上記に鑑みてなされたものであって、物体までの正確な距離を測定できない環境下であっても、簡易な構成で正確に物体を透視することができる物体透視装置を得ることを目的とする。
上述した課題を解決し、目的を達成するために、本開示は、測定エリア内に配置された被測定対象である物体に対して照射された複数の周波数を含んだ電波の反射波に基づいて物体の画像を生成する物体透視装置であって、照射された電波に含まれる各周波数と同一の周波数の反射波のみに基づいたイメージングを行うことによって得られた複数の画像に対し、複数の画像の画素毎に複素加算を行うことで複数の画像を合成した位相合成画像を生成する位相合成画像生成部を備える。また、物体透視装置は、送信アンテナから放射され物体で反射されて受信アンテナで受信された反射波の波形のデータである受信波形データ、物体に対する送信アンテナおよび受信アンテナの位置を示す位置データ、および送信アンテナから放射される電波の周波数のデータである周波数データに基づき、反射波の推定波形と受信波形データとの相関を示す相関ベクトルを画素毎に積算する電力合成方式のマッチドフィルタを用いたイメージング方法によって周波数毎の物体の画像を合成した電力合成画像を生成する電力合成画像生成部を備える。また、物体透視装置は、電力合成画像、受信波形データ、位置データ、および周波数データに基づいて、電力合成画像で特定値よりも大きな反射強度を示す観測座標を探索する処理と、観測座標での相関ベクトルの位相を、画像を位相合成する際の位相の補正に用いる補正位相として算出する処理とを実行する補正位相計算部を備える。位相合成画像生成部は、補正位相、受信波形データ、位置データ、および周波数データに基づき、相関ベクトルを画素毎に複素加算する位相合成方式のマッチドフィルタを用いたイメージング方法によって周波数毎の物体の画像を合成した位相合成画像を生成する。
本開示にかかる物体透視装置は、物体までの正確な距離を測定できない環境下であっても、簡易な構成で正確に物体を透視することができるという効果を奏する。
実施の形態１にかかる物体透視装置の構成を示す図 実施の形態１にかかる物体透視装置が適用するアンテナ移動パスの他の例を説明するための図 実施の形態１にかかる物体透視装置による受信信号の生成処理を説明するための図 実施の形態１にかかる物体透視装置が設定する観測点を説明するための図 実施の形態１にかかる物体透視装置が算出する、反射波の推定波形と受信波形との相関関係を説明するための図 実施の形態１にかかる物体透視装置が用いる高周波信号の周波数の更新パターンを説明するための図 実施の形態１にかかる物体透視装置が備える直交検波部の構成を説明するための図 実施の形態１にかかる物体透視装置が備える波形記録部の構成を示す図 実施の形態１にかかる物体透視装置が備える電力合成画像生成部による電力合成画像の生成処理を説明するための図 実施の形態１にかかる物体透視装置が備える補正位相計算部による周波数毎の補正位相の算出処理を説明するための図 実施の形態１にかかる物体透視装置が備える位相合成画像生成部による位相合成画像の生成処理を説明するための図 実施の形態１にかかる物体透視装置による位相合成画像の生成処理手順を示すフローチャート 実施の形態２にかかる物体透視装置の構成を示す図 実施の形態３にかかる物体透視装置による位相合成画像の生成処理手順を示すフローチャート 実施の形態３にかかる物体透視装置が算出する観測座標集合の組み合わせと、各組み合わせにおけるユークリッド距離の総和を説明するための図 実施の形態１に係る物体透視装置が備える処理回路をプロセッサおよびメモリで実現する場合の処理回路の構成例を示す図 実施の形態１に係る物体透視装置が備える処理回路を専用のハードウェアで構成する場合の処理回路の例を示す図
以下に、本開示の実施の形態にかかる物体透視装置、制御回路、記憶媒体、および物体透視方法を図面に基づいて詳細に説明する。
実施の形態１．
図１は、実施の形態１にかかる物体透視装置の構成を示す図である。物体透視装置１０Ａは、被測定対象である電波散乱体（物体）３０にミリ波、テラヘルツ波などの高周波信号の電波を照射し、電波散乱体３０からの反射波Ｗを用いて電波散乱体３０を透視する装置である。物体透視装置１０Ａは、ＭＦを用いたイメージングを行うことで、電波散乱体３０を透視する。物体透視装置１０Ａは、複数の周波数の高周波信号を用いた位相合成方式でＭＦを実行する。なお、実施の形態１では、物体透視装置１０Ａは、電波散乱体３０までの正確な距離を測定できない環境下で電波散乱体３０を透視するものとする。
物体透視装置１０Ａは、送信アンテナ１１と、受信アンテナ１２と、直交検波部２１と、高周波信号発生部２４と、波形記録部２２と、位置制御部２５と、周波数制御部２３とを備えている。また、物体透視装置１０Ａは、電力合成画像生成部２６と、補正位相計算部２７と、位相合成画像生成部２８とを備えている。なお、送信アンテナ１１は、物体透視装置１０Ａと別構成であってもよい。また、受信アンテナ１２は、物体透視装置１０Ａと別構成であってもよい。
周波数制御部２３は、高周波信号の周波数を指定した周波数指令を高周波信号発生部２４に出力することで、高周波信号発生部２４から出力される高周波信号の周波数を制御する。周波数制御部２３は、高周波信号発生部２４に対して指定する高周波信号の周波数を、予め決められたパターンで更新する。周波数制御部２３は、周波数指令で示される高周波信号の周波数を周波数データとして波形記録部２２に出力する。
高周波信号発生部２４は、周波数制御部２３からの周波数指令に従って、測定に用いる高周波信号を生成し、送信アンテナ１１および直交検波部２１に出力する。
高周波信号発生部２４の出力は、直交検波部２１のローカル入力および送信アンテナ１１のローカル入力に結合されている。高周波信号発生部２４は、送信アンテナ１１および直交検波部２１に、種々の周波数の高周波信号を供給する。高周波信号発生部２４の例は、高周波信号発生器である。
送信アンテナ１１は、高周波信号発生部２４から出力された高周波信号を電波として空間に放射する部品である。送信アンテナ１１は、測定対象である電波散乱体３０に高周波信号の電波を照射する。送信アンテナ１１としては、例えば、ホーンアンテナ、基板上に形成されたパターンアンテナ、複数のアンテナにより構成されるアレーアンテナなどが用いられる。
受信アンテナ１２は、電波散乱体３０からの反射波Ｗを受信する部品である。受信アンテナ１２は、散乱点１Ａ～１Ｃなどの複数の電波散乱点によって反射された反射波Ｗを受信し、直交検波部２１に出力する。受信アンテナ１２としては、例えば、ホーンアンテナ、基板上に形成されたパターンアンテナ、複数のアンテナにより構成されるアレーアンテナなどが用いられる。なお、受信アンテナ１２は、必ずしも送信アンテナ１１と同一構造である必要はない。
位置制御部２５は、送信アンテナ１１の位置および受信アンテナ１２の位置の少なくとも一方を制御する。位置制御部２５は、送信アンテナ１１および受信アンテナ１２の両方の位置を制御してもよいし、送信アンテナ１１のみの位置を制御してもよいし、受信アンテナ１２のみの位置を制御してもよい。
受信アンテナ１２の位置が固定されている場合に、位置制御部２５は、送信アンテナ１１の位置を制御する。送信アンテナ１１の位置が固定されている場合に、位置制御部２５は、受信アンテナ１２の位置を制御する。
送信アンテナ１１の位置は、送信アンテナ１１による電波の放射方向に対応している。以下では、位置制御部２５が、送信アンテナ１１および受信アンテナ１２の両方の位置を制御する場合について説明する。また、送信アンテナ１１と受信アンテナ１２との相対位置は不変であり、位置制御部２５は、送信アンテナ１１および受信アンテナ１２をまとめて移動させるものとする。
位置制御部２５は、送信アンテナ１１および受信アンテナ１２を搭載した搬送機構（図示せず）に接続されており、搬送機構の位置を制御することで送信アンテナ１１および受信アンテナ１２の位置を制御する。
位置制御部２５は、送信アンテナ１１および受信アンテナ１２を、予め決められたアンテナ移動パス７１に従って移動させるとともに、送信アンテナ１１および受信アンテナ１２の位置を示す位置データを波形記録部２２に出力する。送信アンテナ１１の位置が固定されている場合、送信アンテナ１１の位置データは、固定値である。また、受信アンテナ１２の位置が固定されている場合、受信アンテナ１２の位置データは、固定値である。位置制御部２５が波形記録部２２に出力する位置データは、位置制御部２５が搬送機構に出力する移動指令に対応している。
アンテナ移動パス７１は、例えば、図１に示すように、測定対象である電波散乱体３０を取り囲むように移動させる円に沿った経路である。なお、アンテナ移動パス７１は、円に沿った経路に限らない。
図２は、実施の形態１にかかる物体透視装置が適用するアンテナ移動パスの他の例を説明するための図である。アンテナ移動パス７１の他の例であるアンテナ移動パス７２は、特定の平面内で送信アンテナ１１および受信アンテナ１２を上下左右に移動させる経路である。アンテナ移動パス７２に設定される平面をＸＹ平面とすると、位置制御部２５は、Ｘ方向の移動とＹ方向の移動とを種々組み合わせることで、送信アンテナ１１および受信アンテナ１２をＸＹ平面内で種々の方向に移動させる。
物体透視装置１０Ａは、アンテナ移動パス７２を適用する場合、ＸＹ平面内で送信アンテナ１１および受信アンテナ１２を移動させるＸＹステージを用いる。ＸＹステージは、特定平面内の２つの軸であって互いに直交する２つの軸をＸ軸およびＹ軸とした場合に、Ｘ軸方向およびＹ軸方向に移動可能なステージである。
直交検波部２１は、受信アンテナ１２から得られる受信信号を、高周波信号発生部２４から供給される高周波信号によってダウンコンバージョンすることで複素信号であるベースバンド信号を得る。このベースバンド信号には、送信アンテナ１１から出力される高周波信号と受信アンテナ１２で受信した反射波との振幅差および位相差の情報である振幅位相差情報が含まれている。直交検波部２１は、振幅位相差情報を含んだベースバンド信号を波形記録部２２に出力する。直交検波部２１の例は、直交検波回路である。
波形記録部２２は、直交検波部２１から出力されるベースバンド信号をアナログ信号からデジタル信号に変換して得られる受信波形データを記録する。波形記録部２２は、振幅位相差情報を含んだ受信波形データと、周波数データと、位置データとを対応付けして記録する。
電力合成画像生成部２６は、波形記録部２２に記録されている、受信波形データ、周波数データ、および位置データを用いて、電力合成方式のＭＦを用いたイメージング方法によって測定対象である電波散乱体３０の電力合成画像を生成する。電力合成画像は、電力合成方式のＭＦを用いたイメージング方法によって生成された画像である。
電力合成方式は、全ての観測点に対して周波数毎に、推定波形と受信波形との間の相関値、すなわち位相情報を有した相関ベクトルを算出し、得られた相関ベクトルの大きさのみを積算する手法である。電力合成方式は、各々の相関ベクトルの位相がばらばらであっても、位相を考慮しないためロバストにイメージングが行えるというメリットがある。
電力合成画像生成部２６が生成する電力合成画像の座標のうち、反射強度が最大である観測点の位置を示す座標が、最大反射強度座標である。電力合成画像生成部２６は、最大反射強度座標を含む電力合成画像のデータを補正位相計算部２７に出力する。
補正位相計算部２７は、電力合成画像生成部２６から出力された電力合成画像と、波形記録部２２に記録されている、受信波形データ、周波数データ、および位置データとに基づいて、周波数毎の補正位相Arg(cn(λ))を算出する。この場合において、補正位相計算部２７は、電力合成画像生成部２６から出力された電力合成画像から最大反射強度座標を探索し、最大反射強度座標における補正位相Arg(cn(λ))を周波数毎に算出する。補正位相計算部２７は、周波数毎の補正位相Arg(cn(λ))を位相合成画像生成部２８に出力する。
位相合成画像生成部２８は、補正位相計算部２７から出力された周波数毎の補正位相Arg(cn(λ))と、波形記録部２２に記録されている、受信波形データ、周波数データ、および位置データとから、位相合成画像を生成する。この場合において、位相合成画像生成部２８は、位相合成方式のイメージング手法によって位相合成画像を生成する。位相合成画像は、位相合成方式のＭＦを用いたイメージング方法によって生成された画像である。位相合成画像生成部２８は、イメージング結果である位相合成画像を表示装置等の外部装置に出力する。
ここで、物体透視装置１０Ａの測定系で得られる反射波の時系列データと、受信信号との関係について説明する。図３は、実施の形態１にかかる物体透視装置による受信信号の生成処理を説明するための図である。
図３に示すように、異なる位置の散乱点１Ａ～１Ｃからの反射波は、異なる波形を有している。受信アンテナ１２は、反射波の時系列データを重ね合わせることによって、受信信号を生成する。ここでの受信アンテナ１２は、散乱点１Ａからの反射波と、散乱点１Ｂからの反射波と、散乱点１Ｃからの反射波とを重ね合わせることによって、受信信号の波形である受信波形を生成する。このように、受信信号は、全ての散乱点からの反射波を重ね合わせた信号である。
この場合において、散乱点が配置されている測定エリアの位置と、送信アンテナ１１および受信アンテナ１２の位置である送受信アンテナ位置との位置関係が既知であるならば、測定エリア上の任意の座標にある散乱点からの反射波の波形を推定することが可能である。測定エリアの位置は、予め物体透視装置１０Ａに登録しておく。送受信アンテナ位置は、位置制御部２５が波形記録部２２に記録させる位置データに対応している。
物体透視装置１０Ａは、ＭＦを用いたイメージングを行う際に、測定エリア内に網羅的に観測点を設けておき、各々の観測点に対して反射波の推定波形を生成する。図４は、実施の形態１にかかる物体透視装置が設定する観測点を説明するための図である。
物体透視装置１０Ａは、測定エリア内の種々の位置に網羅的に観測点を設定する。図４では、物体透視装置１０Ａが、観測点２Ａ～２Ｃを設定した場合を示している。
物体透視装置１０Ａは、反射波の推定波形と受信波形との相関を取ることによって、各観測点からの反射強度をマッピングする。図５は、実施の形態１にかかる物体透視装置が算出する、反射波の推定波形と受信波形との相関関係を説明するための図である。
物体透視装置１０Ａの電力合成画像生成部２６は、反射波の推定波形と受信波形との相関を取る。ここでの電力合成画像生成部２６は、観測点２Ａにおける反射波の推定波形と受信波形との相関を取る。また、電力合成画像生成部２６は、観測点２Ｂにおける反射波の推定波形と受信波形との相関を取り、観測点２Ｃにおける反射波の推定波形と受信波形との相関を取る。
電力合成画像生成部２６は、反射波の推定波形と受信波形との間の相関値を示す相関結果に基づいて、各観測点からの反射強度をマッピングする。この場合における相関値は、位相情報を持つ相関ベクトルとして得られる。観測点と散乱点との間の距離が近い場合には高い相関値が得られ、観測点の近くに散乱点が存在しない場合には低い相関値が得られる。
物体透視装置１０Ａが実行するＭＦを用いたイメージングでは、複数の周波数を用いて同じ測定を繰り返すことによって、送信アンテナ１１および受信アンテナ１２の配置、および電波散乱体３０の周辺にある散乱物に起因するサイドローブを低減することができる。
物体透視装置１０Ａは、複数の周波数で測定した画像を合成する際に、画素毎に相関ベクトルの大きさを積算する電力合成方式と、画素毎に位相情報を考慮して複素加算を行う位相合成方式との両方を用いる。
位相合成方式では、異なる周波数で測定した同一観測点の相関ベクトルの位相を正確に一致させなければ複素加算時に相関ベクトルが打ち消されてしまうので、物体透視装置１０Ａは、観測点の位置と送受信アンテナ位置との間の位置関係、および測定系の周波数特性の影響を考慮した信号処理を実行する。
ここで、物体透視装置１０Ａが実行する位相合成画像の生成処理の概要について説明する。物体透視装置１０Ａは、複数の周波数で測定した画像を位相合成方式によって合成する場合、異なる周波数で測定した同一観測点の相関ベクトルの位相を一致させるように、観測系に含まれる位相オフセットを除去する。すなわち、物体透視装置１０Ａは、位相合成を実現するために、ある画素を異なる周波数で測定した際に発生する、相関ベクトルの位相オフセットを、測定結果の位相から差し引くことで電波の位相を補正する。
図４に示した測定系において、時刻ｔにおける送受信アンテナ位置と、ある観測点ｎとの間の往復の伝搬距離をｄn（ｔ）とすると、受信信号の位相θ（ｔ，λ）は、以下の式（１）で表される。測定に用いた電波の波長が波長λであり、波長λにおける測定系に含まれる固定位相オフセット量がφ（λ）である。
Figure 0007046277000001
また、測定系に含まれる固定位相オフセットを考慮しない時刻ｔにおける観測点ｎからの反射波位相の推定値θ∧（ｔ,λ）は以下の式（２）で表される。θ∧は、「θ」の真上にハット記号が乗っていることを示している。送受信アンテナ位置と観測点との間の実際の距離と推定値との固定誤差がＥである。
Figure 0007046277000002
時刻ｔにおける観測点ｎからの反射波ｙn（ｔ,λ）と、その推定値ｙn ∧（ｔ,λ）は、式（１）および式（２）を用いて、以下の式（３）および式（４）でそれぞれ表される。ｙn ∧は、「ｙn」の真上にハット記号が乗っていることを示している。
Figure 0007046277000003
Figure 0007046277000004
なお、条件を簡略化するため、反射波の振幅変動は無視できるものとする。この場合、波長λ、観測点ｎにおける相関値ｃn（λ）は、以下の式（５）で表される。
Figure 0007046277000005
式（５）に示した通り、各観測点における相関ベクトルの位相成分は、測定に使用する電波の波長λのみに依存し、観測点の場所に依存しないことが分かる。そこで、本実施の形態の物体透視装置１０Ａは、各観測点における相関ベクトルの位相成分が、測定に使用する電波の波長λのみに依存し、観測点の場所に依存しないことを利用する。物体透視装置１０Ａは、測定エリア内で、特定値よりも大きな反射強度が得られる座標（キャリブレーション座標）での相関ベクトルの位相Arg(cn(λ))を周波数毎に算出し、この位相Arg(cn(λ))を、周波数を合成する際の補正位相として用いる。
物体透視装置１０Ａは、補正位相を用いて周波数毎の画像を補正したうえで、画像を合成することで位相合成画像を生成する。これにより、物体透視装置１０Ａは、異なる周波数で取得した画像を位相合成する際に相関ベクトルが打ち消されることを回避する。
以下、物体透視装置１０Ａが備える各構成要素の詳細な構成および動作について説明する。図６は、実施の形態１にかかる物体透視装置が用いる高周波信号の周波数の更新パターンを説明するための図である。図６に示すグラフの横軸は時間である。図６の上段に示すグラフの縦軸は、送受信アンテナ位置であり、図６の下段に示すグラフの縦軸は、周波数制御部２３が高周波信号発生部２４に対して指定する周波数である。
位置制御部２５は、一定の速度で送信アンテナ１１および受信アンテナ１２を移動させる。高周波信号の周波数の更新パターンＰＴとしては、図６に示すように、位置制御部２５によって設定される送受信アンテナ位置に対して、一定範囲の周波数をステップ上に繰り返すようなパターンが考えられる。
図６の上段に示すグラフでは、送受信アンテナ位置が、順番に位置Ｐ１、位置Ｐ２、位置Ｐ３となるように、位置制御部２５が、送信アンテナ１１および受信アンテナ１２の位置を制御する場合を示している。
また、図６の下段に示すグラフでは、周波数が、周波数Ｆ１、周波数Ｆ２、周波数Ｆ３とステップ上に上昇するように、周波数制御部２３が、周波数を制御する場合を示している。周波数Ｆ１は、送受信アンテナ位置が、位置Ｐ１である場合の周波数である。また、周波数Ｆ２は、送受信アンテナ位置が、位置Ｐ２である場合の周波数であり、周波数Ｆ３は、送受信アンテナ位置が、位置Ｐ３である場合の周波数である。
周波数制御部２３は、周波数を特定の大きさまで到達させると、再び、周波数Ｆ１、周波数Ｆ２、周波数Ｆ３とステップ上に上昇するように周波数を制御する。周波数制御部２３は、これらの処理を繰り返す。
なお、図６に示した更新パターンＰＴは、一例であり、周波数と送受信アンテナ位置との組み合わせは、組み換えられてもよいし、また、測定の順番が入れ替えられてもよい。
ここで、直交検波部２１の具体例について説明する。図７は、実施の形態１にかかる物体透視装置が備える直交検波部の構成を説明するための図である。直交検波部２１は、ミキサ３１，３２と、９０度位相部３３とを有している。９０度位相部３３の例は、９０度位相器である。
直交検波部２１へは、高周波信号発生部２４から供給される高周波信号がローカル入力され、受信アンテナ１２から受信信号が入力される。高周波信号発生部２４からの高周波信号は、９０度位相部３３およびミキサ３２に入力される。受信アンテナ１２からの受信信号は、ミキサ３１およびミキサ３２に入力される。
９０度位相部３３は、高周波信号から同一周波数で９０度位相差を持つ高周波信号を生成してミキサ３１に出力する。これにより、ミキサ３２には、高周波信号発生部２４からの高周波信号がそのまま入力され、ミキサ３１には、高周波信号発生部２４からの高周波信号に対して同一周波数で９０度位相差を持つ高周波信号が入力される。
ミキサ３２は、高周波信号発生部２４からの高周波信号と受信信号とを混合して出力する。ミキサ３１は、９０度位相差を持つ高周波信号と受信信号とを混合して出力する。これにより、直交検波部２１は、受信アンテナ１２から出力された受信信号を、ダウンコンバージョンすることで、複素信号であるベースバンド信号（受信波形データ）を算出し、波形記録部２２に出力する。
つぎに、波形記録部２２の詳細な構成および動作について説明する。波形記録部２２は、直交検波部２１から出力される複素信号である受信波形データと、位置制御部２５から出力される位置データと、周波数制御部２３から出力される周波数データとを記録するメモリ部を有している。
図８は、実施の形態１にかかる物体透視装置が備える波形記録部の構成を示す図である。波形記録部２２は、複素信号で示される受信波形データと、周波数データと、位置データとが対応付けされた対応情報４４を格納するメモリ部４１を有している。
また、波形記録部２２は、直交検波部２１から出力される複素信号をデジタル信号に変換してメモリ部４１に記録させるＡＤＣ（Analog to Digital Converter、アナログデジタル変換器）４２，４３を有している。
ＡＤＣ４２は、直交検波部２１のミキサ３１から出力される複素信号をデジタル信号に変換してメモリ部４１に記録させる。ＡＤＣ４３は、直交検波部２１のミキサ３２から出力される複素信号をデジタル信号に変換してメモリ部４１に記録させる。
図８に示す対応情報４４は、図６に示す送受信アンテナ位置および周波数が設定された場合にメモリ部４１に格納される情報である。
対応情報４４に格納されている周波数データ、位置データ、および受信波形データは、電力合成画像生成部２６、補正位相計算部２７、および位相合成画像生成部２８によって読み出される。例えば、受信波形データの受信波形ｒ１は、周波数データの周波数Ｆ１および位置データの位置Ｐ１に対応している。
つぎに、電力合成画像生成部２６の詳細な構成および動作について説明する。電力合成画像生成部２６は、波形記録部２２に記録されている受信波形データ、周波数データ、および位置データを読み出す。電力合成画像生成部２６は、波形記録部２２に記録されている受信波形データ、周波数データ、および位置データを用いて電力合成画像を生成する。
図９は、実施の形態１にかかる物体透視装置が備える電力合成画像生成部による電力合成画像の生成処理を説明するための図である。図９では、電力合成画像生成部２６が、波形記録部２２から読み出すデータである対応情報４４と、読み出した対応情報４４を用いて生成する電力合成画像との関係を示している。
電力合成画像生成部２６が、波形記録部２２から読み出す対応情報４４は、受信波形データと、周波数データと、位置データとが対応付けされた情報である。電力合成画像生成部２６は、読み出した対応情報４４を同一周波数のデータにグルーピングし、周波数毎にＭＦを用いた画像を生成する。
ここでは、電力合成画像生成部２６が、周波数Ｆ１のデータから画像５１を生成し、周波数Ｆ２のデータから画像５２を生成し、周波数Ｆ３のデータから画像５３を生成する場合を示している。
その後、電力合成画像生成部２６は、周波数毎に生成した画像を、電力合成方式によって合成し、１枚の電力合成画像を生成する。図９では、電力合成画像生成部２６が、画像５１～５３を電力合成方式によって合成し、１枚の電力合成画像５５を生成した場合を示している。
つぎに、補正位相計算部２７の詳細な構成および動作について説明する。補正位相計算部２７は、電力合成画像生成部２６から出力される電力合成画像と、波形記録部２２に記録されている受信波形データと、周波数データと、位置データとから、周波数毎の補正位相Arg(cn(λ))を算出する。
図１０は、実施の形態１にかかる物体透視装置が備える補正位相計算部による周波数毎の補正位相の算出処理を説明するための図である。図１０では、補正位相計算部２７が、電力合成画像５５と、対応情報４４とを用いて周波数毎の補正位相Arg(cn(λ))を算出する場合について説明する。
補正位相計算部２７は、読み出した対応情報４４を同一周波数のデータにグルーピングする。また、補正位相計算部２７は、電力合成画像生成部２６から出力される電力合成画像５５の最大反射強度座標を探索する。ここでは、最大反射強度座標が、観測点２Ａである場合について説明する。
補正位相計算部２７は、最大反射強度座標である観測点２Ａの補正位相Arg(cn(λ))を周波数毎に算出する。具体的には、補正位相計算部２７は、観測点２Ａにおける周波数Ｆ１～Ｆ３の位相を算出する。観測点２Ａにおける周波数Ｆ１～Ｆ３の位相は、それぞれ電力合成画像５５を生成する際の画像５１～５３の補正位相Arg(cn(λ))に用いられる。
すなわち、観測点２Ａにおける周波数Ｆ１の位相Arg(cn(λ))は、電力合成画像５５を生成する際の画像５１の補正位相Arg(cn(λ))に用いられる。同様に、観測点２Ａにおける周波数Ｆ２の位相Arg(cn(λ))は、電力合成画像５５を生成する際の画像５２の補正位相Arg(cn(λ))に用いられる。観測点２Ａにおける周波数Ｆ３の位相Arg(cn(λ))は、電力合成画像５５を生成する際の画像５３の補正位相Arg(cn(λ))に用いられる。図１０では、周波数Ｆ１～Ｆ３の補正位相を、補正位相Arg(cA(λ))で示している。
つぎに、位相合成画像生成部２８の詳細な構成および動作について説明する。位相合成画像生成部２８は、補正位相計算部２７から出力される周波数毎の補正位相Arg(cn(λ))と、波形記録部２２に記録されている受信波形データと、周波数データと、位置データとから、位相合成画像を生成する。
図１１は、実施の形態１にかかる物体透視装置が備える位相合成画像生成部による位相合成画像の生成処理を説明するための図である。図１１では、位相合成画像生成部２８が、補正位相Arg(cn(λ))と、対応情報４４とを用いて位相合成画像５６を生成する場合について説明する。
位相合成画像生成部２８は、波形記録部２２から読み出した対応情報４４を同一周波数のデータにグルーピングし、周波数毎にＭＦを用いた画像５１～５３を生成する。なお、位相合成画像生成部２８は、電力合成画像生成部２６から、ＭＦを用いた周波数毎の画像５１～５３を取得してもよい。
位相合成画像生成部２８は、周波数毎に作成した画像５１～５３に対し、各画素の相関ベクトルを、補正位相計算部２７が計算した周波数毎の補正位相Arg(cn(λ))分だけ逆回転させてから、複素加算することで、位相合成画像５６を生成する。この位相合成画像５６が測定エリアでの画像に対応している。
このように、物体透視装置１０Ａは、測定エリア内で特定値よりも大きな反射強度が得られる座標、例えば、最大反射強度座標で周波数毎に相関ベクトルの位相Arg(cn(λ))を算出する。物体透視装置１０Ａは、この位相Arg(cn(λ))を補正位相Arg(cn(λ))として用いることで、異なる周波数で取得した画像５１～５３を位相合成する際に相関ベクトルが打ち消されることを回避する。
また、物体透視装置１０Ａは、補正位相Arg(cn(λ))を求めるためのキャリブレーション座標として、電力合成方式により合成した画像５１～５３から、特定値よりも反射強度の強い座標、例えば、最大反射強度座標を選択する。これにより、物体透視装置１０Ａは、測定エリアに校正信号発生装置を配置することなく、また既知の座標に校正用の散乱体を設置することなく補正位相Arg(cn(λ))を用いて測定エリアでの位相合成画像５６を生成することができる。
また、物体透視装置１０Ａが、周波数毎に作成した画像５１～５３の位相を校正する際に、校正信号発生装置が定期的に校正信号を発生させる必要がなく、また位相の校正中に反射波の測定を停止する必要がない。また、物体透視装置１０Ａは、測定環境内の正確な位置に校正用の散乱体の一例であるＲＦＩＤ（Radio Frequency Identification、高周波識別子）タグを設置する必要がないので、測定環境が限定されることはない。
つぎに、物体透視装置１０Ａによる位相合成画像の生成処理手順について説明する。図１２は、実施の形態１にかかる物体透視装置による位相合成画像の生成処理手順を示すフローチャートである。物体透視装置１０Ａは、周波数および送受信アンテナ位置を変化させながら反射波を受信し、反射波の情報である受信波形データを記録する（ステップＳ１０）。
具体的には、位置制御部２５が、送受信アンテナ位置を制御し、周波数制御部２３が周波数を制御する。受信アンテナ１２は、電波散乱体３０からの反射波を受信すると、受信信号を、直交検波部２１に出力する。直交検波部２１は、受信アンテナ１２からの受信信号および高周波信号発生部２４からの高周波信号を用いて、複素信号であるベースバンド信号を生成する。
波形記録部２２は、ベースバンド信号から受信波形データを生成して記録する。また、波形記録部２２は、位置制御部２５が出力した位置データと、周波数制御部２３が出力した周波数データと、を受信波形データに対応付けして記録する。
電力合成画像生成部２６は、波形記録部２２に記録されている、受信波形データ、周波数データ、および位置データを用いて、電力合成画像を生成する（ステップＳ２０）。
補正位相計算部２７は、電力合成画像内で最大反射強度座標を探索する（ステップＳ３０）。補正位相計算部２７は、周波数毎に、最大反射強度座標における相関ベクトルの位相を補正位相として算出する（ステップＳ４０）。位相合成画像生成部２８は、補正位相を用いて位相合成画像を生成する（ステップＳ５０）。
このように実施の形態１では、物体透視装置１０Ａが、電波散乱体３０からの反射波に基づいたイメージングを行うことによって得られた複数の画像に対し、複数の画像の画素毎に複素加算を行うことで複数の画像を合成した位相合成画像を生成している。これにより、物体透視装置１０Ａは、複数の周波数を用いて測定した画像を位相合成する際に必要となる、各周波数における補正位相量を、測定エリアに校正信号発生装置を設置することなく、または既知の座標に校正用の散乱体を設置することなく求めることが可能となる。したがって、物体透視装置１０Ａは、電波散乱体３０までの正確な距離を測定できない環境下であっても、複数の周波数の電波を用いて簡易な構成で正確に電波散乱体３０を透視することができる。
実施の形態２．
つぎに、図１３を用いて実施の形態２について説明する。実施の形態２では、送信アンテナ１１および受信アンテナ１２を固定し、電波散乱体３０を移動させる。
図１３は、実施の形態２にかかる物体透視装置の構成を示す図である。図１３の各構成要素のうち図１に示す実施の形態１の物体透視装置１０Ａと同一機能を達成する構成要素については同一符号を付しており、重複する説明は省略する。
物体透視装置１０Ｂは、物体透視装置１０Ａと比較して、測定系が異なる。物体透視装置１０Ｂは、送信アンテナ１１および受信アンテナ１２を固定し、代わりに被測定対象である電波散乱体３０を移動させる。物体透視装置１０Ｂは、被測定対象を載せて回転する回転台５０を備えている。
物体透視装置１０Ｂの位置制御部２５は、回転台５０の回転位置を制御する。回転台５０が回転することによって、電波散乱体３０は、種々の位置で電波を反射する。これにより、物体透視装置１０Ｂの受信アンテナ１２は、物体透視装置１０Ａと同様の反射波Ｗを受信する。したがって、物体透視装置１０Ｂは、物体透視装置１０Ａと同様の処理によって、位相合成画像を生成することができる。
なお、図１３では、電波散乱体３０を回転台５０によって回転させる構成を示しているが、物体透視装置１０Ｂは、回転台５０の代わりに、実施の形態１の図２で説明したＸＹステージのような機構を用いて電波散乱体３０を移動させてもよい。
このように、実施の形態２では、物体透視装置１０Ｂが、回転台５０によって電波散乱体３０を回転させるとともに、物体透視装置１０Ａと同様の処理によって位相合成画像を生成している。これにより、物体透視装置１０Ｂは、物体透視装置１０Ａと同様に、電波散乱体３０までの正確な距離を測定できない環境下であっても、複数の周波数の電波を用いて簡易な構成で正確に電波散乱体３０を透視することができる。
また、物体透視装置１０Ｂは、送受信アンテナ位置を固定することができるので、送信アンテナ１１と高周波信号発生部２４とをつなぐ配線の移動、受信アンテナ１２と直交検波部２１とをつなぐ配線の移動を考慮する必要がない。
実施の形態３．
つぎに、図１４を用いて実施の形態３について説明する。実施の形態３では、補正位相計算部２７における補正位相の計算手順が、実施の形態１，２と異なる。
実施の形態３の物体透視装置は、物体透視装置１０Ａであってもよいし、物体透視装置１０Ｂであってもよい。以下の説明では、実施の形態３の物体透視装置が物体透視装置１０Ａである場合について説明する。
図１４は、実施の形態３にかかる物体透視装置による位相合成画像の生成処理手順を示すフローチャートである。物体透視装置１０Ａは、実施の形態１で説明したステップＳ１０の処理と同様の処理によって、周波数および送受信アンテナ位置を変化させながら反射波を受信し、反射波の情報である受信波形データを記録する（ステップＳ１１０）。
電力合成画像生成部２６は、実施の形態１で説明したステップＳ２０の処理と同様の処理によって、電力合成画像を生成する（ステップＳ１２０）。
補正位相計算部２７は、電力合成画像生成部２６が出力した電力合成画像内の観測点の座標、すなわち観測座標から、反射強度の高い上位Ｍ（Ｍは２以上の自然数）個の観測座標を選定する（ステップＳ１３０）。補正位相計算部２７は、選定したＭ個の観測座標に対し、波形記録部２２に記録されている、受信波形データ、周波数データ、および位置データから、周波数毎の補正位相Arg(cn(λk))（ｎ＝１，・・・Ｍ、ｋ＝１・・・Ｋ）を算出する。ただし、Ｋは観測に使用した周波数の数を表す。
補正位相計算部２７は、以下の式（６）で表されるような、周波数毎の補正位相Arg(cn(λk))を要素に持つ補正位相ベクトルＰn（ｎ＝１・・・Ｍ）を定義する。
Figure 0007046277000006
補正位相計算部２７は、選定したＭ個の観測座標からｘ（ｘは、Ｍ以下の自然数）個の観測座標を抽出する。補正位相計算部２７は、抽出したｘ個の観測座標集合Ｌから２つの観測座標を選び出したlＣ2個の組み合わせについて、補正位相ベクトルのユークリッド距離を求める。補正位相計算部２７は、ユークリッド距離の総和を計算する。
補正位相計算部２７は、ユークリッド距離の総和をMＣx個の組み合わせ全てに対して計算し、補正位相ベクトルのユークリッド距離の総和が最も小さくなる組み合わせを選択する。すなわち、補正位相計算部２７は、上位Ｍ個の観測座標から周波数毎の位相が最も近いｘ個の観測点を選定する（ステップＳ１４０）。
補正位相計算部２７は、選定した観測点の、補正位相ベクトルの平均値を算出する。すなわち、補正位相計算部２７は、選定したｘ個の観測点、それぞれについて周波数毎の補正位相ベクトルを求め、周波数毎に補正位相ベクトルの平均値を算出する（ステップＳ１５０）。
このように、補正位相計算部２７は、上位Ｍ個の観測座標を抽出し、上位Ｍ個の観測座標に含まれるｘ個の観測座標の組み合わせのうち、周波数毎の位相差が最も小さくなる組み合わせとなるｘ個の観測座標を抽出する。さらに、補正位相計算部２７は、抽出したｘ個の観測座標に対し、周波数毎の位相の平均値を補正位相として算出する。
補正位相計算部２７は、算出した補正位相ベクトルの平均値を、位相合成画像生成部２８に出力する周波数毎の補正位相とする。補正位相計算部２７は、ｘ個の観測点について、周波数毎に求めた補正位相（補正位相ベクトルの平均値）を用いて、位相合成画像を生成する（ステップＳ１６０）。なお、ｘはＭと同じ値でもよい。
図１５は、実施の形態３にかかる物体透視装置が算出する観測座標集合の組み合わせと、各組み合わせにおけるユークリッド距離の総和を説明するための図である。図１５では、Ｍ＝４、ｘ＝３とした場合における、観測座標集合Ｌの組み合わせと、各組み合わせにおけるユークリッド距離の総和を式で表している。
このように、実施の形態３によれば、物体透視装置１０Ａが、補正位相の計算に複数の観測点の情報を用いることができるので、実施の形態１，２と比較して、補正位相を安定的に求めることが可能となる。
ここで、物体透視装置１０Ａ，１０Ｂのハードウェア構成について説明する。物体透視装置１０Ａ，１０Ｂは、同様のハードウェア構成を有しているので、以下では実施の形態１にかかる物体透視装置１０Ａのハードウェア構成について説明する。
実施の形態１にかかる物体透視装置１０Ａにおいて、直交検波部２１、波形記録部２２、周波数制御部２３、高周波信号発生部２４、位置制御部２５、電力合成画像生成部２６、補正位相計算部２７、および位相合成画像生成部２８は、処理回路により実現される。処理回路は、メモリに格納されるプログラムを実行するプロセッサおよびメモリであってもよいし、専用のハードウェアであってもよい。処理回路は制御回路とも呼ばれる。
図１６は、実施の形態１に係る物体透視装置が備える処理回路をプロセッサおよびメモリで実現する場合の処理回路の構成例を示す図である。図１６に示す処理回路９０は制御回路であり、プロセッサ９１およびメモリ９２を備える。処理回路９０がプロセッサ９１およびメモリ９２で構成される場合、処理回路９０の各機能は、ソフトウェア、ファームウェア、またはソフトウェアとファームウェアとの組み合わせにより実現される。ソフトウェアまたはファームウェアはプログラムとして記述され、メモリ９２に格納される。処理回路９０では、メモリ９２に記憶されたプログラムをプロセッサ９１が読み出して実行することにより、各機能を実現する。すなわち、処理回路９０は、物体透視装置１０Ａの処理が結果的に実行されることになるプログラムを格納するためのメモリ９２を備える。このプログラムは、処理回路９０により実現される各機能を物体透視装置１０Ａに実行させるためのプログラムであるともいえる。このプログラムは、プログラムが記憶された記憶媒体により提供されてもよいし、通信媒体など他の手段により提供されてもよい。
上記プログラムは、図１２のステップＳ１０からＳ５０の処理を物体透視装置１０Ａに実行させるプログラムであるとも言える。すなわち、上記プログラムは、受信波形データを記録するステップと、電力合成画像を生成するステップと、最大反射強度座標を探索するステップと、最大反射強度座標における相関ベクトルの位相を補正位相として算出するステップと、補正位相を用いて位相合成画像を生成するステップとを物体透視装置１０Ａに実行させるプログラムであるとも言える。
ここで、プロセッサ９１は、例えば、ＣＰＵ（Central Processing Unit）、処理装置、演算装置、マイクロプロセッサ、マイクロコンピュータ、またはＤＳＰ（Digital Signal Processor）などである。また、メモリ９２は、例えば、ＲＡＭ（Random Access Memory）、ＲＯＭ（Read Only Memory）、フラッシュメモリ、ＥＰＲＯＭ（Erasable Programmable ＲＯＭ）、ＥＥＰＲＯＭ（登録商標）（Electrically ＥＰＲＯＭ）などの、不揮発性または揮発性の半導体メモリ、磁気ディスク、フレキシブルディスク、光ディスク、コンパクトディスク、ミニディスク、またはＤＶＤ（Digital Versatile Disc）などが該当する。
図１７は、実施の形態１に係る物体透視装置が備える処理回路を専用のハードウェアで構成する場合の処理回路の例を示す図である。図１７に示す処理回路９３は、例えば、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ＡＳＩＣ（Application Specific Integrated Circuit）、ＦＰＧＡ（Field Programmable Gate Array）、またはこれらを組み合わせたものが該当する。処理回路９３については、一部を専用のハードウェアで実現し、一部をソフトウェアまたはファームウェアで実現するようにしてもよい。このように、処理回路９３は、専用のハードウェア、ソフトウェア、ファームウェア、またはこれらの組み合わせによって、上述の各機能を実現することができる。
以上の実施の形態に示した構成は、一例を示すものであり、別の公知の技術と組み合わせることも可能であるし、実施の形態同士を組み合わせることも可能であるし、要旨を逸脱しない範囲で、構成の一部を省略、変更することも可能である。
１Ａ～１Ｃ 散乱点、２Ａ～２Ｃ 観測点、１０Ａ，１０Ｂ 物体透視装置、１１ 送信アンテナ、１２ 受信アンテナ、２１ 直交検波部、２２ 波形記録部、２３ 周波数制御部、２４ 高周波信号発生部、２５ 位置制御部、２６ 電力合成画像生成部、２７ 補正位相計算部、２８ 位相合成画像生成部、３０ 電波散乱体、３１，３２ ミキサ、３３ ９０度位相部、４１ メモリ部、４２，４３ ＡＤＣ、４４ 対応情報、５０ 回転台、５１～５３ 画像、５５ 電力合成画像、５６ 位相合成画像、７１，７２ アンテナ移動パス、９０，９３ 処理回路、９１ プロセッサ、９２ メモリ、Ｗ 反射波。

Claims (9)

  1. 測定エリア内に配置された被測定対象である物体に対して照射された複数の周波数を含んだ電波の反射波に基づいて前記物体の画像を生成する物体透視装置であって、
    照射された前記電波に含まれる各周波数と同一の周波数の前記反射波のみに基づいたイメージングを行うことによって得られた複数の画像に対し、前記複数の画像の画素毎に複素加算を行うことで前記複数の画像を合成した位相合成画像を生成する位相合成画像生成部と、
    送信アンテナから放射され前記物体で反射されて受信アンテナで受信された前記反射波の波形のデータである受信波形データ、前記物体に対する前記送信アンテナおよび前記受信アンテナの位置を示す位置データ、および前記送信アンテナから放射される電波の周波数のデータである周波数データに基づき、前記反射波の推定波形と前記受信波形データとの相関を示す相関ベクトルを画素毎に積算する電力合成方式のマッチドフィルタを用いたイメージング方法によって前記周波数毎の前記物体の画像を合成した電力合成画像を生成する電力合成画像生成部と、
    前記電力合成画像、前記受信波形データ、前記位置データ、および前記周波数データに基づいて、前記電力合成画像で特定値よりも大きな反射強度を示す観測座標を探索する処理と、前記観測座標での前記相関ベクトルの位相を、前記画像を位相合成する際の位相の補正に用いる補正位相として算出する処理とを実行する補正位相計算部と、
    を備え、
    前記位相合成画像生成部は、前記補正位相、前記受信波形データ、前記位置データ、および前記周波数データに基づき、前記相関ベクトルを画素毎に複素加算する位相合成方式のマッチドフィルタを用いたイメージング方法によって前記周波数毎の前記物体の画像を合成した前記位相合成画像を生成する、
    ことを特徴とする物体透視装置。
  2. 前記位相合成画像生成部は、
    前記位相合成画像を生成する際に、前記相関ベクトルの位相オフセットから前記周波数毎に前記補正位相を差し引いて前記相関ベクトルを画素毎に複素加算する、
    ことを特徴とする請求項１に記載の物体透視装置。
  3. 前記補正位相計算部は、
    前記電力合成画像に含まれる前記観測座標のうち、最大の反射強度を示す最大反射強度座標を探索し、前記最大反射強度座標における前記周波数毎の位相を、前記補正位相として算出する、
    ことを特徴とする請求項１または２に記載の物体透視装置。
  4. 前記補正位相計算部は、
    前記電力合成画像に含まれる前記観測座標のうち、反射強度の上位複数個の観測座標を選定し、前記上位複数個の観測座標における前記周波数毎の位相の平均値を、前記補正位相として算出する、
    ことを特徴とする請求項１または２に記載の物体透視装置。
  5. 前記補正位相計算部は、
    前記上位複数個の観測座標に含まれる複数個の観測座標のうち、前記周波数毎の位相差が最も小さくなる組み合わせとなる複数個の観測座標を抽出し、抽出した複数個の観測座標における前記周波数毎の位相の平均値を、前記補正位相として算出する、
    ことを特徴とする請求項４に記載の物体透視装置。
  6. 前記電波は、サブテラヘルツ帯からテラヘルツ帯の信号である、
    ことを特徴とする請求項１から５の何れか１つに記載の物体透視装置。
  7. 測定エリア内に配置された被測定対象である物体に対して照射された複数の周波数を含んだ電波の反射波に基づいて前記物体の画像を生成する制御回路であって、
    照射された前記電波に含まれる各周波数と同一の周波数の前記反射波のみに基づいたイメージングを行うことによって得られた複数の画像に対し、前記複数の画像の画素毎に複素加算を行うことで前記複数の画像を合成した位相合成画像を生成する処理と、
    送信アンテナから放射され前記物体で反射されて受信アンテナで受信された前記反射波の波形のデータである受信波形データ、前記物体に対する前記送信アンテナおよび前記受信アンテナの位置を示す位置データ、および前記送信アンテナから放射される電波の周波数のデータである周波数データに基づき、前記反射波の推定波形と前記受信波形データとの相関を示す相関ベクトルを画素毎に積算する電力合成方式のマッチドフィルタを用いたイメージング方法によって前記周波数毎の前記物体の画像を合成した電力合成画像を生成する処理と、
    前記電力合成画像、前記受信波形データ、前記位置データ、および前記周波数データに基づいて、前記電力合成画像で特定値よりも大きな反射強度を示す観測座標を探索する処理と、前記観測座標での前記相関ベクトルの位相を、前記画像を位相合成する際の位相の補正に用いる補正位相として算出する処理と、
    を、前記物体を透視する物体透視装置に実施させ、
    前記位相合成画像を生成する処理では、前記補正位相、前記受信波形データ、前記位置データ、および前記周波数データに基づき、前記相関ベクトルを画素毎に複素加算する位相合成方式のマッチドフィルタを用いたイメージング方法によって前記周波数毎の前記物体の画像を合成した前記位相合成画像を生成する、
    ことを特徴とする制御回路。
  8. 測定エリア内に配置された被測定対象である物体に対して照射された複数の周波数を含んだ電波の反射波に基づいて前記物体の画像を生成するプログラムを記憶した記憶媒体であって、
    前記プログラムは、照射された前記電波に含まれる各周波数と同一の周波数の前記反射波のみに基づいたイメージングを行うことによって得られた複数の画像に対し、前記複数の画像の画素毎に複素加算を行うことで前記複数の画像を合成した位相合成画像を生成する処理と、
    送信アンテナから放射され前記物体で反射されて受信アンテナで受信された前記反射波の波形のデータである受信波形データ、前記物体に対する前記送信アンテナおよび前記受信アンテナの位置を示す位置データ、および前記送信アンテナから放射される電波の周波数のデータである周波数データに基づき、前記反射波の推定波形と前記受信波形データとの相関を示す相関ベクトルを画素毎に積算する電力合成方式のマッチドフィルタを用いたイメージング方法によって前記周波数毎の前記物体の画像を合成した電力合成画像を生成する処理と、
    前記電力合成画像、前記受信波形データ、前記位置データ、および前記周波数データに基づいて、前記電力合成画像で特定値よりも大きな反射強度を示す観測座標を探索する処理と、前記観測座標での前記相関ベクトルの位相を、前記画像を位相合成する際の位相の補正に用いる補正位相として算出する処理と、
    を、前記物体を透視する物体透視装置に実施させ、
    前記位相合成画像を生成する処理では、前記補正位相、前記受信波形データ、前記位置データ、および前記周波数データに基づき、前記相関ベクトルを画素毎に複素加算する位相合成方式のマッチドフィルタを用いたイメージング方法によって前記周波数毎の前記物体の画像を合成した前記位相合成画像を生成する、
    ことを特徴とする記憶媒体。
  9. 測定エリア内に配置された被測定対象である物体に対して照射された複数の周波数を含んだ電波の反射波に基づいて前記物体の画像を生成する物体透視方法であって、
    制御回路が、照射された前記電波に含まれる各周波数と同一の周波数の前記反射波のみに基づいたイメージングを行うことによって得られた複数の画像に対し、前記複数の画像の画素毎に複素加算を行うことで前記複数の画像を合成した位相合成画像を生成する第１の画像生成ステップと、
    送信アンテナから放射され前記物体で反射されて受信アンテナで受信された前記反射波の波形のデータである受信波形データ、前記物体に対する前記送信アンテナおよび前記受信アンテナの位置を示す位置データ、および前記送信アンテナから放射される電波の周波数のデータである周波数データに基づき、前記反射波の推定波形と前記受信波形データとの相関を示す相関ベクトルを画素毎に積算する電力合成方式のマッチドフィルタを用いたイメージング方法によって前記周波数毎の前記物体の画像を合成した電力合成画像を生成する第２の画像生成ステップと、
    前記電力合成画像、前記受信波形データ、前記位置データ、および前記周波数データに基づいて、前記電力合成画像で特定値よりも大きな反射強度を示す観測座標を探索する処理と、前記観測座標での前記相関ベクトルの位相を、前記画像を位相合成する際の位相の補正に用いる補正位相として算出する算出ステップと、
    を含み、
    前記第１の画像生成ステップでは、前記補正位相、前記受信波形データ、前記位置データ、および前記周波数データに基づき、前記相関ベクトルを画素毎に複素加算する位相合成方式のマッチドフィルタを用いたイメージング方法によって前記周波数毎の前記物体の画像を合成した前記位相合成画像を生成する、
    ことを特徴とする物体透視方法。
JP2021538716A 2021-01-29 2021-01-29 物体透視装置、制御回路、記憶媒体、および物体透視方法 Active JP7046277B1 (ja)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2021/003238 WO2022162869A1 (ja) 2021-01-29 2021-01-29 物体透視装置、制御回路、記憶媒体、および物体透視方法

Publications (2)

Publication Number Publication Date
JP7046277B1 true JP7046277B1 (ja) 2022-04-01
JPWO2022162869A1 JPWO2022162869A1 (ja) 2022-08-04

Family

ID=81255846

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021538716A Active JP7046277B1 (ja) 2021-01-29 2021-01-29 物体透視装置、制御回路、記憶媒体、および物体透視方法

Country Status (5)

Country Link
US (1) US20230358875A1 (ja)
EP (1) EP4266030A4 (ja)
JP (1) JP7046277B1 (ja)
CN (1) CN116745605A (ja)
WO (1) WO2022162869A1 (ja)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07120548A (ja) * 1993-09-03 1995-05-12 Mitsubishi Electric Corp レーダ装置
JP2005328853A (ja) * 2004-04-23 2005-12-02 Ge Medical Systems Global Technology Co Llc 画像生成方法およびｍｒｉ装置
JP2008145230A (ja) * 2006-12-08 2008-06-26 Nippon Hoso Kyokai <Nhk> イメージング装置
JP2012002797A (ja) * 2010-05-17 2012-01-05 Japan Radio Co Ltd レーダ装置
JP2018059895A (ja) * 2016-09-29 2018-04-12 パナソニック株式会社 マルチレーダシステム

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007256171A (ja) 2006-03-24 2007-10-04 Nec Corp ミリ波画像処理装置及びミリ波画像処理方法
US10690767B2 (en) * 2017-12-11 2020-06-23 Raytheon Company System and method for synthetic aperture radar image formation

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07120548A (ja) * 1993-09-03 1995-05-12 Mitsubishi Electric Corp レーダ装置
JP2005328853A (ja) * 2004-04-23 2005-12-02 Ge Medical Systems Global Technology Co Llc 画像生成方法およびｍｒｉ装置
JP2008145230A (ja) * 2006-12-08 2008-06-26 Nippon Hoso Kyokai <Nhk> イメージング装置
JP2012002797A (ja) * 2010-05-17 2012-01-05 Japan Radio Co Ltd レーダ装置
JP2018059895A (ja) * 2016-09-29 2018-04-12 パナソニック株式会社 マルチレーダシステム

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Yulei Qian，SAR Image Formation From Azimuth Periodically Gapped Raw Data Via Complex ISTA，IEEE Conference Proceedings of Asian and Pacific Conference on Synthetic Aperture Radar (APSAR)，2019年，Vol.2019 No.APSAR，pp.1-5
YULEI QIAN: "SAR Image Formation From Azimuth Periodically Gapped Raw Data Via Complex ISTA", IEEE CONFERENCE PROCEEDINGS OF ASIAN AND PACIFIC CONFERENCE ON SYNTHETIC APERTURE RADAR (APSAR), vol. 2019, JPN6021013308, 2019, pages 1 - 5, XP033751297, ISSN: 0004613863, DOI: 10.1109/APSAR46974.2019.9048466 *
若山俊夫，合成開口レーダの技術動向，電子情報通信学会誌，2012年，Vol.95 No.2，pp.142-144
若山俊夫: "合成開口レーダの技術動向", 電子情報通信学会誌, vol. 95, no. 2, JPN6021013309, 2012, pages 142 - 144, ISSN: 0004613864 *

Also Published As

Publication number Publication date
US20230358875A1 (en) 2023-11-09
CN116745605A (zh) 2023-09-12
WO2022162869A1 (ja) 2022-08-04
EP4266030A4 (en) 2024-02-21
JPWO2022162869A1 (ja) 2022-08-04
EP4266030A1 (en) 2023-10-25

Similar Documents

Publication Publication Date Title
US11860262B2 (en) System and methods for three dimensional modeling of an object using a radio frequency device
US10830805B2 (en) Method, system and computer program for correcting a radiation pattern
JP2019128235A (ja) レーダ装置
JP2011237257A (ja) レーダ装置及び目標探知方法
JPWO2018147025A1 (ja) 物体検知装置、物体検知方法、及びプログラム
US20200209377A1 (en) Position measurement device, position measurement method, and program recording medium
JP2015094760A5 (ja)
JP6678554B2 (ja) アンテナ測定装置
López-Pérez et al. Surface accuracy improvement of the Yebes 40 meter radiotelescope using microwave holography
US11789116B2 (en) Multi-direction phased array calibration
JP7046277B1 (ja) 物体透視装置、制御回路、記憶媒体、および物体透視方法
WO2019180767A1 (ja) 物体検知装置、物体検知方法、及びコンピュータ読み取り可能な記録媒体
WO2020189198A1 (ja) 電磁界分布生成プログラム、フィルタ生成プログラム及び電磁界分布生成装置
JP6289672B2 (ja) 合成開口レーダ信号処理装置及び合成開口レーダ信号処理プログラム
JP6491263B2 (ja) 到来電波を測定する測定装置及びその測定方法
JP5458124B2 (ja) 電磁波イメージング装置及び電磁波イメージング方法
JPWO2020136871A1 (ja) レーダ信号処理装置及びレーダ信号処理方法
JP4519478B2 (ja) 目標距離測定装置
JP2000321314A (ja) アンテナ位相中心位置測定法およびアンテナ装置
Ferrero et al. Impact of mechanical accuracy in mmW spherical measurements
JP6415392B2 (ja) 信号処理装置
JP7367876B2 (ja) 物体検知装置及び物体検知方法及びプログラム
JP2005083814A (ja) レーダ装置
JP5972316B2 (ja) 干渉測定システム及び干渉測定方法
Zhao et al. Sparse equivalent source method for radiated emission analysis of shielded circuits

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210630

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210630

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20210630

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20211012

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20211210

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220322

R150 Certificate of patent or registration of utility model

Ref document number: 7046277

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150