JP6505343B2 - レーザレーダ装置 - Google Patents
レーザレーダ装置 Download PDFInfo
- Publication number
- JP6505343B2 JP6505343B2 JP2018563960A JP2018563960A JP6505343B2 JP 6505343 B2 JP6505343 B2 JP 6505343B2 JP 2018563960 A JP2018563960 A JP 2018563960A JP 2018563960 A JP2018563960 A JP 2018563960A JP 6505343 B2 JP6505343 B2 JP 6505343B2
- Authority
- JP
- Japan
- Prior art keywords
- signal
- spectrum
- frequency
- light
- unit
- 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
- 238000001228 spectrum Methods 0.000 claims description 411
- 230000003287 optical effect Effects 0.000 claims description 133
- 230000010354 integration Effects 0.000 claims description 116
- 238000012937 correction Methods 0.000 claims description 106
- 238000012545 processing Methods 0.000 claims description 103
- 238000000034 method Methods 0.000 claims description 30
- 238000010521 absorption reaction Methods 0.000 description 45
- 239000000443 aerosol Substances 0.000 description 32
- 238000010586 diagram Methods 0.000 description 23
- 230000003321 amplification Effects 0.000 description 10
- 238000003199 nucleic acid amplification method Methods 0.000 description 10
- 230000006866 deterioration Effects 0.000 description 9
- 230000015654 memory Effects 0.000 description 9
- 239000013307 optical fiber Substances 0.000 description 8
- 230000003595 spectral effect Effects 0.000 description 7
- 238000000862 absorption spectrum Methods 0.000 description 6
- 238000001514 detection method Methods 0.000 description 5
- 238000009825 accumulation Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 238000002617 apheresis Methods 0.000 description 1
- 239000005427 atmospheric aerosol Substances 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- GQYHUHYESMUTHG-UHFFFAOYSA-N lithium niobate Chemical compound [Li+].[O-][Nb](=O)=O GQYHUHYESMUTHG-UHFFFAOYSA-N 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/4802—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- 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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/06—Systems determining position data of a target
- G01S17/08—Systems determining position data of a target for measuring distance only
- G01S17/10—Systems determining position data of a target for measuring distance only using transmission of interrupted, pulse-modulated waves
-
- 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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
- G01S17/95—Lidar systems specially adapted for specific applications for meteorological use
-
- 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/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/4804—Auxiliary means for detecting or identifying lidar signals or the like, e.g. laser illuminators
-
- 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/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/481—Constructional features, e.g. arrangements of optical elements
- G01S7/4811—Constructional features, e.g. arrangements of optical elements common to transmitter and receiver
- G01S7/4813—Housing arrangements
-
- 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/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/481—Constructional features, e.g. arrangements of optical elements
- G01S7/4818—Constructional features, e.g. arrangements of optical elements using optical fibres
-
- 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/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/483—Details of pulse systems
- G01S7/486—Receivers
- G01S7/487—Extracting wanted echo signals, e.g. pulse 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/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/497—Means for monitoring or calibrating
-
- 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/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/481—Constructional features, e.g. arrangements of optical elements
- G01S7/4814—Constructional features, e.g. arrangements of optical elements of transmitters alone
- G01S7/4815—Constructional features, e.g. arrangements of optical elements of transmitters alone using multiple transmitters
-
- 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)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Optics & Photonics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Optical Radar Systems And Details Thereof (AREA)
Description
このレーザレーダ装置は、大気中の分子に吸収される波長である吸収波長のレーザ光と、大気中の分子に吸収されない波長である非吸収波長のレーザ光とを発生させる光源を有している。大気中の分子は、固有の吸収スペクトルを有しており、吸収波長のレーザ光の波長は、大気中の分子の吸収スペクトルの波長に対応している。
また、このレーザレーダ装置は、光源から発生された吸収波長のレーザ光及び非吸収波長のレーザ光を大気中に放射し、大気中のエアロゾルにより散乱された吸収波長のレーザ光の散乱光及び非吸収波長のレーザ光の散乱光を受信する光アンテナを有している。
以下、吸収波長のレーザ光の散乱光を「吸収波長散乱光」と称し、非吸収波長のレーザ光の散乱光を「非吸収波長散乱光」と称する。
この光カプラは、吸収波長散乱光と吸収波長のレーザ光との合波光を出力するとともに、非吸収波長散乱光と非吸収波長のレーザ光との合波光を出力する。
以下、吸収波長散乱光と吸収波長のレーザ光との合波光を「吸収波長合波光」と称し、非吸収波長散乱光と非吸収波長のレーザ光との合波光を「非吸収波長合波光」と称する。
このレーザレーダ装置は、光受信機から出力された吸収波長合波光の受信信号と非吸収波長合波光の受信信号とから大気中の分子の濃度を算出する信号処理器を有している。
信号処理器は、信号対雑音比であるSN比を高めるために、同一時間レンジ内の複数の信号スペクトルにおける同一周波数の信号成分同士を積算することで、複数の信号スペクトルを積算する。
信号処理器は、複数の信号スペクトルを積算すると、積算後の信号スペクトルにおける複数の信号成分を比較し、複数の信号成分の中で、最も大きな信号成分を吸収波長合波光のピーク強度として特定する。
信号処理器は、SN比を高めるために、同一時間レンジ内の複数の信号スペクトルにおける同一周波数の信号成分同士を積算することで、複数の信号スペクトルを積算する。
信号処理器は、複数の信号スペクトルを積算すると、積算後の信号スペクトルにおける複数の信号成分を比較し、複数の信号成分の中で、最も大きな信号成分を非吸収波長合波光のピーク強度として特定する。
光学的厚さは、各々の時間レンジにおける吸収波長合波光のピーク強度と非吸収波長合波光のピーク強度との比である。
しかし、風速あるいは風向の変化に伴って大気中の散乱体の位置が変化している環境下では、レーザ光の散乱光にドップラーシフトが生じる。
ドップラーシフトが生じると、ピーク強度に対応する周波数であるピーク位置が変化するため、複数の信号スペクトルにおける同一周波数の信号成分同士を積算しても、SN比が高まらなくなる。その結果、吸収波長合波光のピーク強度の特定精度及び非吸収波長合波光のピーク強度の特定精度が劣化するため、分子濃度の算出精度が劣化してしまうという課題があった。
図1はこの発明の実施の形態1によるレーザレーダ装置を示す構成図である。
図1において、基準光源1は波長切替器2から大気中の分子に吸収される波長のレーザ光である第1のレーザ光Lb1を出力する旨を示す波長切替信号を受けると、第1のレーザ光Lb1を光分配器3に出力する光源である。
また、基準光源1は波長切替器2から大気中の分子に吸収されない波長のレーザ光である第2のレーザ光Lb2を出力する旨を示す波長切替信号を受けると、第2のレーザ光Lb2を光分配器3に出力する。
例えば、観測対象の分子が二酸化炭素であれば、吸収波長は1572.192nm、非吸収波長は1573.193nmなどに設定される。
なお、基準光源1としては、例えば、線幅が数MHz以下のレーザ光を出力する半導体レーザ、ファイバレーザ、固体レーザ、あるいは、これらのレーザの組み合わせが用いられる。
ここでは、波長切替器2から出力される波長切替信号に従って基準光源1が第1のレーザ光Lb1又は第2のレーザ光Lb2を出力する例を示しているが、基準光源1として、第1のレーザ光Lb1及び第2のレーザ光Lb2を同時に出力可能な光源を用いるようにしてもよい。
光分配器3は基準光源1から出力された第1のレーザ光Lb1及び第2のレーザ光Lb2を所定の比率でパルス変調器4と光カプラ8に分配する。
パルス変調器4は光分配器3と光ファイバケーブルで接続されており、例えば音響光学素子又はニオブ酸リチウム結晶を用いた変調素子で構成されている。
パルス変調器4は光分配器3により分配された第1のレーザ光Lb1及び第2のレーザ光Lb2のパルス変調を行う。
光増幅器5はパルス変調器4から出力された第1のレーザ光Lb1及び第2のレーザ光Lb2の強度を増幅し、強度増幅後の第1のレーザ光Lb1及び強度増幅後の第2のレーザ光Lb2を光サーキュレータ6に出力する。
光サーキュレータ6は光増幅器5と光ファイバケーブルで接続されており、例えば偏光ビームスプリッタ及び波長板で構成されている。
光サーキュレータ6は光増幅器5から出力された強度増幅後の第1のレーザ光Lb1及び強度増幅後の第2のレーザ光Lb2を光アンテナ7に出力し、光アンテナ7から出力された第1のレーザ光Lb1の散乱光である第1の散乱光Sb1及び第2のレーザ光Lb2の散乱光である第2の散乱光Sb2を光カプラ8に出力する。
光アンテナ7は光サーキュレータ6から出力された第1のレーザ光Lb1及び第2のレーザ光Lb2のビーム径を大きくして第1のレーザ光Lb1及び第2のレーザ光Lb2を大気中に放射する。
また、光アンテナ7は大気中のエアロゾル(散乱体)又はハードターゲットにより散乱された第1のレーザ光Lb1の散乱光である第1の散乱光Sb1及び第2のレーザ光Lb2の散乱光である第2の散乱光Sb2を受信して、第1の散乱光Sb1及び第2の散乱光Sb2を光サーキュレータ6に出力する。
光カプラ8は光サーキュレータ6から第1の散乱光Sb1が出力されると、第1の散乱光Sb1と光分配器3により分配された第1のレーザ光Lb1とを合波し、第1の散乱光Sb1と第1のレーザ光Lb1との合波光である第1の合波光Cb1を50対50の比率で2つに分岐して光受信機9に出力する。
光カプラ8は光サーキュレータ6から第2の散乱光Sb2が出力されると、第2の散乱光Sb2と光分配器3により分配された第2のレーザ光Lb2とを合波し、第2の散乱光Sb2と第2のレーザ光Lb2との合波光である第2の合波光Cb2を50対50の比率で2つに分岐して光受信機9に出力する。
光受信機9は光カプラ8から出力された2つの第1の合波光Cb1をヘテロダイン検波して、第1の合波光Cb1の受信信号である第1の受信信号R1をアナログデジタル変換器であるA/D変換器10に出力する。
また、光受信機9は光カプラ8から出力された2つの第2の合波光Cb2をヘテロダイン検波して、第2の合波光Cb2の受信信号である第2の受信信号R2をA/D変換器10に出力する。
第1の受信信号R1及び第2の受信信号R2は電気信号である。
信号処理器11はA/D変換器10と電気信号ケーブルで接続されており、波長切替器2から出力される波長切替信号を制御する。
また、信号処理器11はA/D変換器10から出力された受信データDR1,DR2に基づいて大気中の分子の濃度Cを算出する。
図3はこの発明の実施の形態1によるレーザレーダ装置の信号処理器11を示すハードウェア構成図である。
図2及び図3において、スペクトル算出部21はフーリエ変換部22、スペクトル積算処理部23及び雑音除去部24を備えている。
スペクトル算出部21はA/D変換器10から出力された受信データDR1の信号スペクトルである第1の信号スペクトルSP1を算出するとともに、A/D変換器10から出力された受信データDR2の信号スペクトルである第2の信号スペクトルSP2を算出する。
フーリエ変換部22はA/D変換器10から出力された受信データDR1に対する高速フーリエ変換処理であるFFTを実施することで、受信データDR1の周波数領域の信号として、時間レンジ毎の複数の第1の信号スペクトルSP1を算出する。
また、フーリエ変換部22はA/D変換器10から出力された受信データDR2に対するFFTを実施することで、受信データDR2の周波数領域の信号として、時間レンジ毎の複数の第2の信号スペクトルSP2を算出する。
スペクトル積算処理部23はフーリエ変換部22により算出された同一時間レンジ内の複数の第1の信号スペクトルSP1を積算して、積算後の第1の信号スペクトルSP1を雑音除去部24に出力する処理を実施する。
また、スペクトル積算処理部23はフーリエ変換部22により算出された同一時間レンジ内の複数の第2の信号スペクトルSP2を積算して、積算後の第2の信号スペクトルSP2を雑音除去部24に出力する処理を実施する。
雑音除去部24はスペクトル積算処理部23による積算後の第1の信号スペクトルSP1に含まれている雑音成分を除去して、雑音成分除去後の第1の信号スペクトルSP1を周波数シフト補正部25に出力する処理を実施する。
また、雑音除去部24はスペクトル積算処理部23による積算後の第2の信号スペクトルSP2に含まれている雑音成分を除去して、雑音成分除去後の第2の信号スペクトルSP2を周波数シフト補正部25に出力する処理を実施する。
周波数シフト補正部25は光アンテナ7から放射された第1のレーザ光Lb1の周波数に対する雑音除去部24による雑音成分除去後の第1の信号スペクトルSP1の周波数のシフトを補正する処理を実施する。
また、周波数シフト補正部25は光アンテナ7から放射された第2のレーザ光Lb2の周波数に対する雑音除去部24による雑音成分除去後の第2の信号スペクトルSP2の周波数のシフトを補正する処理を実施する。
スペクトル積算部26は周波数シフト補正部25により周波数のシフトが補正された複数の第1の信号スペクトルSP1を積算して、積算後の第1の信号スペクトルSP1をピーク強度算出部28に出力する処理を実施する。
周波数シフト補正部25により周波数のシフトが補正された複数の第1の信号スペクトルSP1は、同一時間レンジ内の複数の時刻の信号スペクトルが積算された補正後の信号スペクトルである。
なお、スペクトル積算部26は、周波数シフト補正部25により周波数のシフトが補正された複数の第1の信号スペクトルSP1における同一周波数の信号成分同士を積算することで、複数の第1の信号スペクトルSP1を積算する。
周波数シフト補正部25により周波数のシフトが補正された複数の第2の信号スペクトルSP2は、同一時間レンジ内の複数の時刻の信号スペクトルが積算された補正後の信号スペクトルである。
なお、スペクトル積算部26は、周波数シフト補正部25により周波数のシフトが補正された複数の第2の信号スペクトルSP2における同一周波数の信号成分同士を積算することで、複数の第2の信号スペクトルSP2を積算する。
分子濃度算出部27はスペクトル積算部26により積算された第1の信号スペクトルSP1とスペクトル積算部26により積算された第2の信号スペクトルSP2とから大気中の分子の濃度Cを算出する。
ピーク強度算出部28は例えば図3に示すピーク強度算出回路46によって実現されるものである。
ピーク強度算出部28はスペクトル積算部26による積算後の第1の信号スペクトルSP1のピーク強度である第1のピーク強度Ponを算出する処理を実施する。
また、ピーク強度算出部28はスペクトル積算部26による積算後の第2の信号スペクトルSP2のピーク強度である第2のピーク強度Poffを算出する処理を実施する。
分子濃度算出処理部29はピーク強度算出部28により算出された第1のピーク強度Ponとピーク強度算出部28により算出された第2のピーク強度Poffとから大気中の分子の濃度Cを算出する処理を実施する。
ただし、信号処理器11の構成要素は、専用のハードウェアで実現されるものに限るものではなく、信号処理器11がソフトウェア、ファームウェア、または、ソフトウェアとファームウェアとの組み合わせで実現されるものであってもよい。
コンピュータのメモリは、例えば、RAM(Random Access Memory)、ROM(Read Only Memory)、フラッシュメモリ、EPROM(Erasable Programmable Read Only Memory)、EEPROM(Electrically Erasable Programmable Read Only Memory)などの不揮発性又は揮発性の半導体メモリ、磁気ディスク、フレキシブルディスク、光ディスク、コンパクトディスク、ミニディスク、DVD(Digital Versatile Disc)などが該当する。
信号処理器11がソフトウェア又はファームウェアで実現される場合、フーリエ変換部22、スペクトル積算処理部23、雑音除去部24、周波数シフト補正部25、スペクトル積算部26、ピーク強度算出部28及び分子濃度算出処理部29の処理手順をコンピュータに実行させるためのプログラムをメモリ51に格納し、コンピュータのプロセッサ52がメモリ51に格納されているプログラムを実行するようにすればよい。
図5及び図6は信号処理器11がソフトウェア又はファームウェアで実現される場合の処理手順を示すフローチャートである。
また、図3では、信号処理器11の構成要素のそれぞれが専用のハードウェアで実現される例を示し、図4では、信号処理器11がソフトウェア又はファームウェアで実現される例を示しているが、信号処理器11における一部の構成要素が専用のハードウェアで実現され、残りの構成要素がソフトウェア又はファームウェアで実現されるものであってもよい。
基準光源1は、波長切替器2から大気中の分子に吸収される波長のレーザ光である第1のレーザ光Lb1を出力する旨を示す波長切替信号を受けると、第1のレーザ光Lb1を光分配器3に出力する。
また、基準光源1は、波長切替器2から大気中の分子に吸収されない波長のレーザ光である第2のレーザ光Lb2を出力する旨を示す波長切替信号を受けると、第2のレーザ光Lb2を光分配器3に出力する。
この実施の形態1では、基準光源1が、第1のレーザ光Lb1と第2のレーザ光Lb2を交互に出力するものとするが、基準光源1が、第1のレーザ光Lb1と第2のレーザ光Lb2を同時に出力するものであってもよい。
光分配器3は、基準光源1から第2のレーザ光Lb2を受けると、第2のレーザ光Lb2を所定の比率でパルス変調器4と光カプラ8に分配する。
パルス変調器4は、光分配器3から第2のレーザ光Lb2を受けると、第2のレーザ光Lb2のパルス変調を行う。
例えば、パルス変調器4は、第1のレーザ光Lb1及び第2のレーザ光Lb2のパルス変調として、第1のレーザ光Lb1及び第2のレーザ光Lb2の光強度の形状がガウス形状又は矩形形状になるような第1のレーザ光Lb1及び第2のレーザ光Lb2の強度変調を行う。
また、パルス変調器4は、第1のレーザ光Lb1及び第2のレーザ光Lb2のパルス変調として、第1のレーザ光Lb1及び第2のレーザ光Lb2の周波数をシフトさせる周波数シフトを行う。第1のレーザ光Lb1及び第2のレーザ光Lb2の周波数シフト量としては、例えば、20MHz〜200MHzの範囲である。
これにより、パルス変調器4からパルス変調後の第1のレーザ光Lb1及びパルス変調後の第2のレーザ光Lb2が光増幅器5に出力される。
光増幅器5は、パルス変調器4からパルス変調後の第2のレーザ光Lb2を受けると、パルス変調後の第2のレーザ光Lb2の強度を増幅し、強度増幅後の第2のレーザ光Lb2を光サーキュレータ6に出力する。
光サーキュレータ6は、光増幅器5から強度増幅後の第1のレーザ光Lb1を受けると、強度増幅後の第1のレーザ光Lb1を光アンテナ7に出力する。
光サーキュレータ6は、光増幅器5から強度増幅後の第2のレーザ光Lb2を受けると、強度増幅後の第2のレーザ光Lb2を光アンテナ7に出力する。
光アンテナ7は、光サーキュレータ6から強度増幅後の第2のレーザ光Lb2を受けると、強度増幅後の第2のレーザ光Lb2のビーム径を大きくして、強度増幅後の第2のレーザ光Lb2を大気中に放射する。
光アンテナ7は、強度増幅後の第1のレーザ光Lb1を大気中に放射した後、大気中のエアロゾル又はハードターゲットにより散乱された第1のレーザ光Lb1の散乱光である第1の散乱光Sb1を受信して、第1の散乱光Sb1を光サーキュレータ6に出力する。
光アンテナ7は、強度増幅後の第2のレーザ光Lb2を大気中に放射した後、大気中のエアロゾル又はハードターゲットにより散乱された第2のレーザ光Lb2の散乱光である第2の散乱光Sb2を受信して、第2の散乱光Sb2を光サーキュレータ6に出力する。
光カプラ8は、光サーキュレータ6から第1の散乱光Sb1を受けると、第1の散乱光Sb1と光分配器3から出力された第1のレーザ光Lb1とを合波し、第1の散乱光Sb1と第1のレーザ光Lb1との合波光である第1の合波光Cb1を50対50の比率で2つに分岐して光受信機9に出力する。
光カプラ8は、光サーキュレータ6から第2の散乱光Sb2を受けると、第2の散乱光Sb2と光分配器3から出力された第2のレーザ光Lb2とを合波し、第2の散乱光Sb2と第2のレーザ光Lb2との合波光である第2の合波光Cb2を50対50の比率で2つに分岐して光受信機9に出力する。
光受信機9は、光カプラ8から50対50の比率で分岐された2つの第2の合波光Cb2を受けると、2つの第2の合波光Cb2をヘテロダイン検波して、第2の合波光Cb2の受信信号である第2の受信信号R2をA/D変換器10に出力する。
光受信機9による第1の合波光Cb1及び第2の合波光Cb2のヘテロダイン検波は、同相雑音を軽減することが可能なバランスド検出である。
ここで、図7は光受信機9から出力される第1の受信信号R1及び第2の受信信号R2の波形例を示す説明図である。
図7Aは第1の受信信号R1の波形例を示し、図7Bは第2の受信信号R2の波形例を示している。
図7において、横軸は時間、縦軸は受信信号の信号強度である。
図7では、時間レンジをNに分けて、第1時間レンジビン〜第N時間レンジビンの信号強度を示している。
信号処理器11は、波長切替器2から出力される波長切替信号を制御する。
この実施の形態1では、上述したように、基準光源1が、第1のレーザ光Lb1と第2のレーザ光Lb2を交互に出力するように、信号処理器11が波長切替器2から出力される波長切替信号を制御している。
信号処理器11は、A/D変換器10からデジタル信号である受信データDR1,DR2を受けると、受信データDR1,DR2に基づいて大気中の分子の濃度Cを算出する。
以下、信号処理器11による分子濃度Cの算出処理を具体的に説明する。
また、フーリエ変換部22は、A/D変換器10から受信データDR2が出力されると(図5のステップST3:YESの場合)、当該受信データDR2に対するFFTを実施することで、当該受信データDR2の周波数領域の信号として、時間レンジ毎の複数の第2の信号スペクトルSP2を算出する(図5のステップST4)。
この実施の形態1では、フーリエ変換部22は、第1時間レンジビン〜第N時間レンジビンの受信データDR1,DR2に対するFFTを実施することで、第1時間レンジビン〜第N時間レンジビンにおける複数の第1の信号スペクトルSP1及び複数の第2の信号スペクトルSP2を算出している。
なお、フーリエ変換部22は、A/D変換器10から受信データDR1,DR2が出力されていないときは(図5のステップST1,ST3:NOの場合)、A/D変換器10から受信データDR1又は受信データDR2が出力されるまで待機する。
図8Aはフーリエ変換部22により算出された第1時間レンジビン及び第2時間レンジビンの第1の信号スペクトルSP1を示し、図8Bはフーリエ変換部22により算出された第1時間レンジビン及び第2時間レンジビンの第2の信号スペクトルSP2を示している。
図8において、横軸は周波数、縦軸は信号強度である。
図8は、大気中のエアロゾルの位置に変化がない例を示している。
大気中のエアロゾルの位置に変化がない場合、同一時間レンジ内の複数の第1の信号スペクトルSP1のピーク強度Ponに対応する周波数であるピーク位置は同じ位置になる。また、時間レンジが異なっていても、図8Aに示すように、第1の信号スペクトルSP1のピーク強度Ponに対応する周波数であるピーク位置は同じ位置となる。
また、大気中のエアロゾルの位置に変化がない場合、同一時間レンジ内の複数の第2の信号スペクトルSP2のピーク強度Poffに対応する周波数であるピーク位置は同じ位置になる。また、時間レンジが異なっていても、図8Bに示すように、第2の信号スペクトルSP2のピーク強度Poffに対応する周波数であるピーク位置は同じ位置となる。
図8A及び図8Bの例では、第1時間レンジビン及び第2時間レンジビンのいずれであっても、第1の信号スペクトルSP1のピーク強度Pon及び第2の信号スペクトルSP2のピーク強度Poffが、全て第6周波数レンジビンに存在している。
図9に示すように、第1時間レンジビン内の複数の第1の信号スペクトルSP1のピーク強度Ponに対応する周波数が異なっており、第1時間レンジビン内の複数の第2の信号スペクトルSP2のピーク強度Poffに対応する周波数とが異なっている。
即ち、第1時間レンジビン内の時刻t1の第1の信号スペクトルSP1(t1)のピーク強度Ponに対応する周波数と、時刻t2の第1の信号スペクトルSP1(t2)のピーク強度Ponに対応する周波数とが異なっている。また、第1時間レンジビン内の時刻t1の第2の信号スペクトルSP2(t1)のピーク強度Poffに対応する周波数と、時刻t2の第2の信号スペクトルSP2(t2)のピーク強度Poffに対応する周波数とが異なっている。
図9では、第1時間レンジビンだけを例示しているが、実際には、第2時間レンジビン〜第N時間レンジビンも存在している。Nは3以上の整数である。
一方、大気中のエアロゾルの位置に変化がある場合、図9に示すように、同一時間レンジ内であっても、ピーク強度に対応する周波数であるピーク位置が異なる位置となる。このため、後段のスペクトル積算処理部23が、同一時間レンジ内の複数の第1の信号スペクトルSP1及び複数の第2の信号スペクトルSP2をそれぞれ積算しても、所望の信号成分が積み上がらず、SN比が高まらない。
ただし、大気中のエアロゾルの位置に変化がある場合でも、複数の第1の信号スペクトルSP1及び複数の第2の信号スペクトルSP2をそれぞれ短時間積算(例えば4秒)することで、後段の周波数シフト補正部25において、第1の信号スペクトルSP1及び第2の信号スペクトルSP2の周波数のシフト量fd1,fd2を算出することが可能になる程度にはSN比が高まる。
1つの時間レンジの信号スペクトルのみでは、信号強度が小さいため、後段の周波数シフト補正部25において、第1の信号スペクトルSP1の周波数のシフト量fd1及び第2の信号スペクトルSP2の周波数のシフト量fd2を算出することは困難である。
また、スペクトル積算処理部23は、フーリエ変換部22が時間レンジ毎の複数の第2の信号スペクトルSP2を算出すると、同一時間レンジ内の複数の第2の信号スペクトルSP2を積算し、積算後の第2の信号スペクトルSP2を雑音除去部24に出力する(図5のステップST6)。
同一時間レンジ内の複数の第1の信号スペクトルSP1の積算処理は、同一時間レンジ内の複数の第1の信号スペクトルSP1における同一周波数の信号成分同士を積算するものである。
また、同一時間レンジ内の複数の第2の信号スペクトルSP2の積算処理は、同一時間レンジ内の複数の第2の信号スペクトルSP2における同一周波数の信号成分同士を積算するものである。
なお、大気中のエアロゾルの位置が風の影響で変化している場合、スペクトル積算処理部23が、同一周波数の信号成分同士を積算しても、SN比が大きく高まることはないが、後段の周波数シフト補正部25が、周波数のシフト量fd1,fd2を算出することが可能になる程度にはSN比が高まる。
また、雑音除去部24は、スペクトル積算処理部23から積算後の第2の信号スペクトルSP2を受けると、積算後の第2の信号スペクトルSP2に含まれている雑音成分を除去し、雑音成分除去後の第2の信号スペクトルSP2を周波数シフト補正部25に出力する(図5のステップST8)。
図10Aは雑音除去部24により雑音成分の除去処理が行われる前の信号スペクトル及び雑音スペクトルを示し、図10Bは雑音除去部24により雑音成分の除去処理が行われた後の信号スペクトルを示している。
以下、雑音除去部24による雑音成分の除去処理を具体的に説明する。
雑音除去部24は、雑音成分の除去処理を行う前に、光受信機9の熱雑音及び光増幅器5から放射される増幅光雑音などの雑音成分を観測して、その雑音成分を示す雑音スペクトルを記録する。
雑音除去部24は、スペクトル積算処理部23による積算後の第1の信号スペクトルSP1から記録している雑音スペクトルを減算することで、第1の信号スペクトルSP1に含まれている雑音成分を除去し、図10Bに示すような第1の信号スペクトルSP1を周波数シフト補正部25に出力する。
また、雑音除去部24は、スペクトル積算処理部23による積算後の第2の信号スペクトルSP2から記録している雑音スペクトルを減算することで、第2の信号スペクトルSP2に含まれている雑音成分を除去し、図10Bに示すような第2の信号スペクトルSP2を周波数シフト補正部25に出力する。
この周波数のシフト量fd1は、大気中のエアロゾルの位置の変化に伴う第1の散乱光Sb1のドップラーシフト量に相当する。
ドップラーシフト量の算出処理自体は公知の技術であるため詳細な説明を省略するが、光アンテナ7から放射された第1のレーザ光Lb1の周波数と、第1の散乱光Sb1の周波数との差分からドップラーシフト量を求めることができる。
周波数シフト補正部25は、周波数のシフト量fd1を算出すると、雑音成分除去後の第1の信号スペクトルSP1の周波数をシフト量fd1だけシフトすることで、雑音成分除去後の第1の信号スペクトルSP1の周波数をドップラーシフトのないときの周波数と一致させる補正を行う(図5のステップST9)。
周波数シフト補正部25は、周波数のシフト量fd2を算出すると、雑音成分除去後の第2の信号スペクトルSP2の周波数をシフト量fd2だけシフトすることで、雑音成分除去後の第2の信号スペクトルSP2の周波数をドップラーシフトのないときの周波数と一致させる補正を行う(図5のステップST10)。
図11は、周波数シフト補正部25によって同一時間レンジ内の複数の第1の信号スペクトルSP1の周波数がシフトされるとともに、同一時間レンジ内の複数の第2の信号スペクトルSP2の周波数がシフトされることで、エアロゾルに対応する信号成分の周波数が、同じ周波数レンジビンに存在するようになることを表している。
図11の例では、シフト補正後の複数の第1の信号スペクトルSP1に含まれているエアロゾルに対応する信号成分の周波数はf1になっている。また、シフト補正後の複数の第2の信号スペクトルSP2に含まれているエアロゾルに対応する信号成分の周波数はf2になっている。
ただし、前段のスペクトル積算処理部23によって積算処理が行われているため、図11における1つの第1の信号スペクトルSP1は、図9における複数の第1の信号スペクトルSP1と対応している。また、図11における1つの第2の信号スペクトルSP2は、図9における複数の第2の信号スペクトルSP2と対応している。
例えば、スペクトル積算処理部23における第1の信号スペクトルSP1及び第2の信号スペクトルSP2の積算回数が5回であれば、図11における第1の信号スペクトルSP1(t1)は、図9における同一時間レンジ内の5つの時刻t=1〜5の第1の信号スペクトルSP1(t1)〜SP1(t5)と対応している。図11における第2の信号スペクトルSP2(t1)は、図9における同一時間レンジ内の5つの時刻t=1〜5の第2の信号スペクトルSP2(t1)〜SP2(t5)と対応している。
また、図11における第1の信号スペクトルSP1(t2)は、図9における同一時間レンジ内の5つの時刻t=6〜10の第1の信号スペクトルSP(t6)〜SP1(t10)1と対応している。図11における第2の信号スペクトルSP2(t2)は、図9における同一時間レンジ内の5つの時刻t=6〜10の第2の信号スペクトルSP2(t6)〜SP2(t10)と対応している。
第1の信号スペクトルSP1及び第2の信号スペクトルSP2の周波数の補正方法は、上記の例に限るものではない。
例えば、周波数シフト補正部25が、同一時間レンジ内の複数の第1の信号スペクトルSP1の周波数が同じドップラーシフト量となるように、複数の第1の信号スペクトルSP1の周波数を補正するようにしてもよい。また、周波数シフト補正部25が、同一時間レンジ内の複数の第2の信号スペクトルSP2の周波数が同じドップラーシフト量となるように、複数の第2の信号スペクトルSP2の周波数を補正するようにしてもよい。
具体的には、同一時間レンジ内の複数の第1の信号スペクトルSP1のドップラーシフト量が、同一時間レンジ内の基準とする第1の信号スペクトルSP1のドップラーシフト量と同じになるように、周波数シフト補正部25が、同一時間レンジ内の複数の第1の信号スペクトルSP1の周波数を補正する例が考えられる。同一時間レンジ内の基準とする第1の信号スペクトルSP1としては、例えば、時刻t1の第1の信号スペクトルSP1が考えられる。
また、同一時間レンジ内の複数の第2の信号スペクトルSP2のドップラーシフト量が、同一時間レンジ内の基準とする第2の信号スペクトルSP2のドップラーシフト量と同じになるように、周波数シフト補正部25が、同一時間レンジ内の複数の第2の信号スペクトルSP2の周波数を補正する例が考えられる。
信号処理器11が処理を開始してからの時間が観測時間を経過すると(ステップST11:YESの場合)、ステップST1〜ST10の繰り返し処理は終了する。
なお、観測時間中、周波数シフト補正部25によって第1時間レンジビン〜第M時間レンジビンにおける複数の第1の信号スペクトルSP1及び複数の第2の信号スペクトルSP2の周波数のシフトが補正された場合、周波数のシフトが補正された第1時間レンジビン〜第M時間レンジビンにおける複数の第1の信号スペクトルSP1及び複数の第2の信号スペクトルSP2が、スペクトル積算部26に記憶される。
また、スペクトル積算部26は、記憶している同一時間レンジ内の複数の第2の信号スペクトルSP2を積算し、積算後の第2の信号スペクトルSP2をピーク強度算出部28に出力する(図6のステップST22)。
同一時間レンジ内の複数の第1の信号スペクトルSP1の積算処理は、各々の第1の信号スペクトルSP1における同一周波数の信号成分同士を積算するものである。
また、同一時間レンジ内の複数の第2の信号スペクトルSP2の積算処理は、各々の第2の信号スペクトルSP2における同一周波数の信号成分同士を積算するものである。
周波数シフト補正部25によって周波数がシフトされている複数の第1の信号スペクトルSP1が積算されるとともに、周波数シフト補正部25によって周波数がシフトされている複数の第2の信号スペクトルSP2が積算されているので、大気中のエアロゾルの位置に変化がある場合でも、所望の信号成分であるエアロゾルに対応する信号成分が積み上がり、SN比が高まる。
また、ピーク強度算出部28は、スペクトル積算部26による積算後の第2の信号スペクトルSP2における各々の信号成分を比較して、最も大きな信号成分である第2のピーク強度Poffを特定する(図6のステップST24)。
以下、分子濃度算出処理部29による分子濃度Cの算出処理を具体的に説明する。
また、分子濃度算出処理部29は、第(m−1)時間レンジビンの第1のピーク強度Ponであるピーク強度Pon(m−1)と、第(m−1)時間レンジビンの第2のピーク強度Poffであるピーク強度Poff(m−1)とを取得する。
光学的厚さODは、第m時間レンジビンにおけるピーク強度Pon(m)とピーク強度Poff(m)との比である。
式(1)において、Poff(m−1)/Pon(m−1)は、光学的厚さODを規格化するために、Pon(m)/Poff(m)に乗算している。
特に、大気分子が水蒸気で、波長1.5um帯での吸収スペクトル強度が小さい場合において分子濃度計測の精度を上げる場合には、吸収波長と非吸収波長のピーク強度を高いSN比で計測して、差分を算出する必要がある。この場合、数分から数十分といった長時間の積算が必要となるが、この実施の形態1によって、大気中のエアロゾルの位置に変化が生じる場合でも、分子濃度の算出精度の劣化を防止することができる効果を奏する。
上記実施の形態1では、同一時間レンジ内の複数の第1の信号スペクトルSP1及び同一時間レンジ内の複数の第2の信号スペクトルSP2の積算処理として、スペクトル積算処理部23が、複数の第1の信号スペクトルSP1及び複数の第2の信号スペクトルSP2における同一周波数の信号成分同士を積算する例を示している。
この実施の形態2では、スペクトル積算処理部23が、複数の第1の信号スペクトルSP1及び複数の第2の信号スペクトルSP2における異なる周波数の信号成分同士を積算する例を説明する。
具体的には、以下の通りである。
スペクトル積算処理部23には、予め、同一時間レンジ内の時刻t1〜t6の第1の信号スペクトルSP1及び第2の信号スペクトルSP2における信号成分の積算を行う周波数の組み合わせが複数用意されている。
例えば、同一時間レンジ内の時刻t1〜t6の信号スペクトルの積算を行う際の周波数の組み合わせについてのパターンとして、組み合わせパターンAと、組み合わせパターンBと、組み合わせパターンCとを例示する。
積算後の信号スペクトルの第m周波数レンジビンにおける信号成分
=時刻t1の第m周波数レンジビンにおける信号成分
+時刻t2の第m周波数レンジビンにおける信号成分
+時刻t3の第m周波数レンジビンにおける信号成分
+時刻t4の第(m+1)周波数レンジビンにおける信号成分
+時刻t5の第(m+1)周波数レンジビンにおける信号成分
+時刻t6の第(m+1)周波数レンジビンにおける信号成分
積算後の信号スペクトルの第m周波数レンジビンにおける信号成分
=時刻t1の第m周波数レンジビンにおける信号成分
+時刻t2の第m周波数レンジビンにおける信号成分
+時刻t3の第(m+1)周波数レンジビンにおける信号成分
+時刻t4の第(m+1)周波数レンジビンにおける信号成分
+時刻t5の第(m+2)周波数レンジビンにおける信号成分
+時刻t6の第(m+2)周波数レンジビンにおける信号成分
積算後の信号スペクトルの第m周波数レンジビンにおける信号成分
=時刻t1の第m周波数レンジビンにおける信号成分
+時刻t2の第(m+1)周波数レンジビンにおける信号成分
+時刻t3の第(m+2)周波数レンジビンにおける信号成分
+時刻t4の第(m+3)周波数レンジビンにおける信号成分
+時刻t5の第(m+4)周波数レンジビンにおける信号成分
+時刻t6の第(m+5)周波数レンジビンにおける信号成分
また、スペクトル積算処理部23は、組み合わせパターンCに従って時刻t1〜時刻t6における第1の信号スペクトルSP1及び第2の信号スペクトルSP2をそれぞれ積算する。
スペクトル積算処理部23は、選択した第1の信号スペクトルSP1を積算後の第1の信号スペクトルSP1として、ピーク強度算出部28に出力する。
また、スペクトル積算処理部23は、組み合わせパターンA,B,Cに従ってそれぞれ積算した第2の信号スペクトルSP2における最大の信号成分を比較し、積算後の3つの第2の信号スペクトルSP2の中から、最大の信号成分が最も大きい第2の信号スペクトルSP2を選択する。
スペクトル積算処理部23は、選択した第2の信号スペクトルSP2を積算後の第2の信号スペクトルSP2として、ピーク強度算出部28に出力する。
上記実施の形態1では、周波数シフト補正部25が、第1のレーザ光Lb1の周波数に対する第1の信号スペクトルSP1の周波数のシフトを補正するとともに、第2のレーザ光Lb2の周波数に対する第2の信号スペクトルSP2の周波数のシフトを補正する例を説明している。
これにより、大気中のエアロゾルの位置に変化が生じる場合でも、分子濃度Cの算出精度の劣化を防止することができる。
しかし、乱流のように時間的に風速又は風向が変化する環境下では、フーリエ変換部22により算出される第1の信号スペクトルSP1及び第2の信号スペクトルSP2に含まれているエアロゾルに対応する信号成分の周波数帯域が広がることがある。
エアロゾルに対応する信号成分の周波数帯域が広がることで、エアロゾルに対応する信号成分が複数の周波数レンジビンに存在し、第1のピーク強度Ponの特定精度及び第2のピーク強度Poffの特定精度が劣化することがある。
この実施の形態3では、時間的に風速又は風向が変化する環境下でも、第1のピーク強度Ponの特定精度及び第2のピーク強度Poffの特定精度の劣化を防止して、分子濃度の算出精度の劣化を防止する例を説明する。
図12はこの発明の実施の形態3によるレーザレーダ装置の信号処理器11を示す構成図である。図12において、図2と同一符号は同一または相当部分を示すので説明を省略する。
周波数シフト補正部30は、図2の周波数シフト補正部25と同様に、例えば図3に示す周波数シフト補正回路44によって実現されるものである。
周波数シフト補正部30は、図2の周波数シフト補正部25と同様に、第1のレーザ光Lb1の周波数に対する雑音成分除去後の複数の第1の信号スペクトルSP1の周波数のシフトを補正するとともに、第2のレーザ光Lb2の周波数に対する雑音成分除去後の複数の第2の信号スペクトルSP2の周波数のシフトを補正する処理を実施する。
周波数シフト補正部30は雑音除去部24による雑音成分除去後の第1の信号スペクトルSP1及び第2の信号スペクトルSP2のそれぞれにおいて、閾値Th以上の信号成分を有する1つ以上の周波数の中から、いずれか1つの周波数を選択し、選択していない残りの周波数の信号成分を選択した周波数に移動させることで、選択した周波数の信号成分を閾値Th以上の信号成分の総和とする波形の補正処理を行う。
以下、周波数シフト補正部30による波形の補正処理を具体的に説明する。
図13では、第1時間レンジビン及び第2時間レンジビンのそれぞれにおいて、1つの第1の信号スペクトルSP1だけを代表的に図示している。また、第1時間レンジビン及び第2時間レンジビンのそれぞれにおいて、1つの第2の信号スペクトルSP2だけを代表的に図示している。
周波数シフト補正部30は、雑音除去部24による雑音成分除去後の第1の信号スペクトルSP1及び雑音成分除去後の第2の信号スペクトルSP2における各々の信号成分と事前に設定された閾値Thとを比較して、閾値Th以上の信号成分を有する周波数を特定する。事前に設定された閾値Thとして、エアロゾルに対応する信号成分の想定される信号強度の半分程度の値が周波数シフト補正部30に設定されている。
一方、第1時間レンジビンにおける第1の信号スペクトルSP1については、閾値Th以上の信号成分を有する周波数は4つである。このため、第1時間レンジビンにおける第1の信号スペクトルSP1については、周波数シフト補正部30によって波形の補正処理が行われる。
また、第1時間レンジビンにおける第2の信号スペクトルSP2については、閾値Th以上の信号成分を有する周波数は3つである。このため、第1時間レンジビンにおける第2の信号スペクトルSP2については、周波数シフト補正部30によって波形の補正処理が行われる。
周波数シフト補正部30は、例えば、1つの周波数fcを選択すると、残り3つの周波数fa,fb,fdの信号成分を周波数fcに移動させることで、周波数fcの信号成分を4つの周波数fa,fb,fc,fdの信号成分の総和とする。
ここでは、周波数fcを選択しているが、これは一例に過ぎず、周波数fa,fb,fdのいずれかを選択してもよい。
また、周波数fa〜fdの範囲内の周波数を選択し、その選択した周波数に4つの周波数fa,fb,fc,fdの信号成分を移動させることで、その選択した周波数の信号成分を4つの周波数fa,fb,fc,fdの信号成分の総和とするようにしてもよい。
周波数シフト補正部30は、例えば、1つの周波数ffを選択すると、残り2つの周波数fe,fgの信号成分を周波数ffに移動させることで、周波数ffの信号成分を3つの周波数fe,ff,fgの信号成分の総和とする。
ここでは、周波数ffを選択しているが、これは一例に過ぎず、周波数fe又は周波数fgのいずれかを選択してもよい。
また、周波数fe〜fgの範囲内の周波数を選択し、その選択した周波数に3つの周波数fe,ff,fgの信号成分を移動させることで、その選択した周波数の信号成分を3つの周波数fe,ff,fgの信号成分の総和とするようにしてもよい。
周波数シフト補正部30は、波形の補正処理を実施した後、図2の周波数シフト補正部25と同様に、複数の第1の信号スペクトルSP1の周波数のシフトを補正する処理及び複数の第2の信号スペクトルSP2の周波数のシフトを補正する処理を実施する。
この実施の形態3では、乱流のように時間的に風速又は風向が変化する環境下でも、第1のピーク強度Ponの特定精度及び第2のピーク強度Poffの特定精度が向上するため、分子濃度の算出精度の劣化を防止することができる。
上記実施の形態1〜3では、分子濃度算出部27が、大気中の分子の濃度Cを算出する例を示しているが、この実施の形態4では、分子濃度算出部27が、分子濃度Cの分布を求める例を説明する。
方向切替器12は例えばミラーなどを備えており、信号処理器11から出力される制御信号に従って光アンテナ7からレーザ光が放射される方向を切り替える切替器である。
分子濃度算出部31はピーク強度算出部28及び分子濃度算出処理部32を備えている。
分子濃度算出部31は、図2の分子濃度算出部27と同様に、スペクトル積算部26により積算された第1の信号スペクトルSP1とスペクトル積算部26により積算された第2の信号スペクトルSP2とから大気中の分子の濃度Cを算出する処理を実施する。
分子濃度算出部31は、図2の分子濃度算出部27と異なり、光アンテナ7からレーザ光が放射される方向を切り替える旨を示す制御信号を方向切替器12に出力し、方向切替器12によりレーザ光が放射される方向が切り替えられる毎に、大気中の分子の濃度Cを算出することで、分子濃度Cの分布を求める処理を実施する。
分子濃度算出処理部32は、図2の分子濃度算出処理部29と同様に、ピーク強度算出部28により算出された第1のピーク強度Ponとピーク強度算出部28により算出された第2のピーク強度Poffとから大気中の分子の濃度Cを算出する処理を実施する。
分子濃度算出処理部32は、図2の分子濃度算出処理部29と異なり、光アンテナ7からレーザ光が放射される方向を切り替える旨を示す制御信号を方向切替器12に出力する。
また、分子濃度算出処理部32は、方向切替器12によりレーザ光が放射される方向が切り替えられる毎に、大気中の分子の濃度Cを算出することで、分子濃度Cの分布を求める処理を実施する。
方向切替器12及び分子濃度算出処理部32以外は、上記実施の形態1と同様であるため、ここでは方向切替器12及び分子濃度算出処理部32の動作だけを説明する。
分子濃度算出処理部32は、光アンテナ7からレーザ光が放射される方向を示す制御信号を方向切替器12に出力する。
方向切替器12は、光アンテナ7から放射されたレーザ光を分子濃度算出処理部32から出力された制御信号が示す方向に設定する。
分子濃度算出処理部32は、大気中の分子の濃度Cを算出すると、光アンテナ7からレーザ光が放射される方向を切り替える旨を示す制御信号を方向切替器12に出力する。
方向切替器12は、分子濃度算出処理部32から出力された制御信号に従って光アンテナ7からレーザ光が放射される方向を切り替える。
方向切替器12と分子濃度算出処理部32が、上記の処理を繰り返し実施することで、分子濃度算出処理部32は、複数の放射方向に対応する分子の濃度Cを算出する。これにより、大気中の分子濃度Cの分布が求められる。
上記実施の形態1〜4では、基準光源1が、第1のレーザ光Lb1と第2のレーザ光Lb2とを出力する機能を有するものについて示したが、第1のレーザ光Lb1を出力する基準光源1aと、第2のレーザ光Lb2を出力する基準光源1bとを別個に設けるようにしてもよい。
基準光源1aは第1のレーザ光Lb1を出力する光源である。
基準光源1bは第2のレーザ光Lb2を出力する光源である。
光スイッチドライバ13は波長切替器2から第1のレーザ光Lb1を出力する旨を示す波長切替信号を受けると、基準光源1aと光分配器3を接続する旨を示す制御信号を光スイッチ14に出力する。
また、光スイッチドライバ13は波長切替器2から第2のレーザ光Lb2を出力する旨を示す波長切替信号を受けると、基準光源1bと光分配器3を接続する旨を示す制御信号を光スイッチ14に出力する。
また、光スイッチ14は光スイッチドライバ13から基準光源1bと光分配器3を接続する旨を示す制御信号を受けると、基準光源1bと光分配器3を接続する。
図16は、基準光源1a,1b、光スイッチドライバ13及び光スイッチ14を図1のレーザレーダ装置に適用している例を示しているが、基準光源1a,1b、光スイッチドライバ13及び光スイッチ14を図14のレーザレーダ装置に適用するようにしてもよい。
その他の点では、上記実施の形態1〜4と同様であるため詳細な説明を省略する。
この実施の形態5でも、上記実施の形態1と同様に、大気中のエアロゾルの位置に変化が生じる場合でも、分子濃度Cの算出精度の劣化を防止することができる効果を奏する。
Claims (8)
- 大気中の分子に吸収される波長の第1のレーザ光が大気中の散乱体によって散乱された第1の散乱光と前記第1のレーザ光との合波光である第1の合波光の信号スペクトルとして、同一時間レンジ内の複数の第1の信号スペクトルを算出するとともに、前記分子に吸収されない波長の第2のレーザ光が前記散乱体によって散乱された第2の散乱光と前記第2のレーザ光との合波光である第2の合波光の信号スペクトルとして、同一時間レンジ内の複数の第2の信号スペクトルを算出するスペクトル算出部と、
前記第1のレーザ光の周波数に対する前記複数の第1の信号スペクトルの周波数のシフトを補正するとともに、前記第2のレーザ光の周波数に対する前記複数の第2の信号スペクトルの周波数のシフトを補正する周波数シフト補正部と、
前記周波数シフト補正部により補正された複数の第1の信号スペクトルを積算するとともに、前記周波数シフト補正部により補正された複数の第2の信号スペクトルを積算するスペクトル積算部と、
前記スペクトル積算部により積算された第1及び第2の信号スペクトルから大気中の分子の濃度を算出する分子濃度算出部と、を備え、
前記スペクトル算出部は、
前記第1及び第2の合波光を検波する光受信機から前記第1の合波光の受信信号である第1の受信信号をフーリエ変換することで、前記第1の受信信号における時間レンジ毎の複数の第1の信号スペクトルを算出するとともに、前記光受信機から前記第2の合波光の受信信号である第2の受信信号をフーリエ変換することで、前記第2の受信信号における時間レンジ毎の複数の第2の信号スペクトルを算出するフーリエ変換部と、
前記フーリエ変換部により算出された同一時間レンジ内の複数の第1の信号スペクトルを積算するとともに、前記フーリエ変換部により算出された同一時間レンジ内の複数の第2の信号スペクトルを積算するスペクトル積算処理部とを備え、
前記周波数シフト補正部は、前記第1のレーザ光の周波数に対する前記スペクトル積算処理部により積算された第1の信号スペクトルの周波数のシフトを補正するとともに、前記第2のレーザ光の周波数に対する前記スペクトル積算処理部により積算された第2の信号スペクトルの周波数のシフトを補正し、
複数の信号スペクトルにおける信号成分の積算を行う周波数の組み合わせが複数用意されており、
前記スペクトル積算処理部は、用意されている周波数の組み合わせ毎に、前記フーリエ変換部により算出された同一時間レンジ内の複数の第1の信号スペクトルにおける当該組み合わせが示す周波数の信号成分を積算し、複数の組み合わせに対応する積算後の信号スペクトルの中から、いずれか1つの信号スペクトルを積算後の第1の信号スペクトルとして選択し、用意されている周波数の組み合わせ毎に、前記フーリエ変換部により算出された同一時間レンジ内の複数の第2の信号スペクトルにおける当該組み合わせが示す周波数の信号成分を積算し、複数の組み合わせに対応する積算後の信号スペクトルの中から、いずれか1つの信号スペクトルを積算後の第2の信号スペクトルとして選択する、レーザレーダ装置。 - 大気中の分子に吸収される波長の第1のレーザ光が大気中の散乱体によって散乱された第1の散乱光と前記第1のレーザ光との合波光である第1の合波光の信号スペクトルとして、同一時間レンジ内の複数の第1の信号スペクトルを算出するとともに、前記分子に吸収されない波長の第2のレーザ光が前記散乱体によって散乱された第2の散乱光と前記第2のレーザ光との合波光である第2の合波光の信号スペクトルとして、同一時間レンジ内の複数の第2の信号スペクトルを算出するスペクトル算出部と、
前記第1のレーザ光の周波数に対する前記複数の第1の信号スペクトルの周波数のシフトを補正するとともに、前記第2のレーザ光の周波数に対する前記複数の第2の信号スペクトルの周波数のシフトを補正する周波数シフト補正部と、
前記周波数シフト補正部により補正された複数の第1の信号スペクトルを積算するとともに、前記周波数シフト補正部により補正された複数の第2の信号スペクトルを積算するスペクトル積算部と、
前記スペクトル積算部により積算された第1及び第2の信号スペクトルから大気中の分子の濃度を算出する分子濃度算出部と
を備え、
前記周波数シフト補正部は、
前記第1及び第2の信号スペクトルの周波数のシフトを補正する前に、
前記第1及び第2の信号スペクトルのそれぞれにおいて、信号成分が閾値以上である1つ以上の周波数の中から、いずれか1つの周波数を選択し、前記選択した周波数の信号成分を、前記1つ以上の周波数の信号成分の総和とする補正処理を行う、レーザレーダ装置。 - 前記スペクトル算出部は、
前記スペクトル積算処理部により積算された第1の信号スペクトルに含まれている雑音成分を除去して、雑音成分除去後の第1の信号スペクトルを前記周波数シフト補正部に出力するとともに、前記スペクトル積算処理部により積算された第2の信号スペクトルに含まれている雑音成分を除去して、雑音成分除去後の第2の信号スペクトルを前記周波数シフト補正部に出力する雑音除去部を備えたことを特徴とする請求項1記載のレーザレーダ装置。 - 前記周波数シフト補正部は、
選択していない周波数の信号成分を前記選択した周波数に移動させることで、前記選択した周波数の信号成分を、前記1つ以上の周波数の信号成分の総和とする補正処理を行う請求項2に記載のレーザレーダ装置。 - 前記分子濃度算出部は、
前記スペクトル積算部により積算された第1の信号スペクトルのピーク強度である第1のピーク強度を算出するとともに、前記スペクトル積算部により積算された第2の信号スペクトルのピーク強度である第2のピーク強度を算出するピーク強度算出部と、
前記ピーク強度算出部により算出された第1のピーク強度と前記ピーク強度算出部により算出された第2のピーク強度とから大気中の分子の濃度を算出する分子濃度算出処理部とを備えたことを特徴とする請求項1から請求項4までの何れか1項に記載のレーザレーダ装置。 - 前記分子濃度算出処理部は、前記ピーク強度算出部により算出された第1のピーク強度と前記ピーク強度算出部により算出された第2のピーク強度との比を算出し、前記第1のピーク強度と前記第2のピーク強度との比から大気中の分子の濃度を算出することを特徴とする請求項5記載のレーザレーダ装置。
- 前記第1のレーザ光及び前記第2のレーザ光を出力する光源と、
前記光源から出力された第1及び第2のレーザ光を大気中にそれぞれ放射し、前記散乱体により散乱された前記第1及び第2のレーザ光の散乱光である第1及び第2の散乱光をそれぞれ受信する光アンテナと、
前記光アンテナにより受信された第1の散乱光と前記第1のレーザ光を合波して、前記第1の散乱光と前記第1のレーザ光との合波光である第1の合波光を出力するとともに、前記光アンテナにより受信された第2の散乱光と前記第2のレーザ光を合波して、前記第2の散乱光と前記第2のレーザ光との合波光である第2の合波光を出力する光カプラと、
前記光カプラから出力された第1の合波光を検波して、前記第1の合波光を前記スペクトル算出部に出力するとともに、前記光カプラから出力された第2の合波光を検波して、前記第2の合波光を前記スペクトル算出部に出力する光受信機とを備えたことを特徴とする請求項1から請求項6までの何れか1項に記載のレーザレーダ装置。 - 前記光アンテナから前記第1及び第2のレーザ光が放射される方向を切り替える方向切替器を備え、
前記分子濃度算出部は、前記方向切替器により第1及び第2のレーザ光が放射される方向が切り替えられる毎に、大気中の分子の濃度を算出することで、分子濃度の分布を求めることを特徴とする請求項7記載のレーザレーダ装置。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2017/002282 WO2018138766A1 (ja) | 2017-01-24 | 2017-01-24 | レーザレーダ装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP6505343B2 true JP6505343B2 (ja) | 2019-04-24 |
JPWO2018138766A1 JPWO2018138766A1 (ja) | 2019-04-25 |
Family
ID=62978123
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2018563960A Active JP6505343B2 (ja) | 2017-01-24 | 2017-01-24 | レーザレーダ装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US11543524B2 (ja) |
EP (1) | EP3564652B1 (ja) |
JP (1) | JP6505343B2 (ja) |
CN (1) | CN110234979B (ja) |
WO (1) | WO2018138766A1 (ja) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3428685B1 (en) * | 2016-04-05 | 2020-09-09 | Mitsubishi Electric Corporation | Laser radar device |
CN116964488A (zh) * | 2021-02-26 | 2023-10-27 | 三菱电机株式会社 | 激光雷达装置 |
JP6980163B1 (ja) * | 2021-02-26 | 2021-12-15 | 三菱電機株式会社 | 信号処理装置、信号処理方法及びライダ装置 |
JP7034398B1 (ja) * | 2021-06-30 | 2022-03-11 | 三菱電機株式会社 | レーザレーダ装置 |
DE102021121211A1 (de) * | 2021-08-16 | 2023-02-16 | OSRAM Opto Semiconductors Gesellschaft mit beschränkter Haftung | Lasersystem zur Entfernungsmessung und Verfahren zur Entfernungsmessung |
CN116879227B (zh) * | 2023-07-07 | 2024-06-11 | 安徽农业大学 | 一种基于激光光谱的快速气体反演方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1183629A (ja) * | 1997-09-02 | 1999-03-26 | Sumitomo Heavy Ind Ltd | 差分吸収ライダ計測装置 |
JP2009503486A (ja) * | 2005-07-29 | 2009-01-29 | キネテイツク・リミテツド | レーザ測定のデバイスおよび方法 |
JP2013083467A (ja) * | 2011-10-06 | 2013-05-09 | Japan Aerospace Exploration Agency | 光学式遠隔気流計測装置の有色ノイズ低減方法とその装置 |
WO2014024508A1 (ja) * | 2012-08-08 | 2014-02-13 | 三菱電機株式会社 | レーダ装置 |
WO2014041852A1 (ja) * | 2012-09-14 | 2014-03-20 | 三菱電機株式会社 | レーザーレーダ装置及び測定対象物の速度算出方法 |
JP2015108557A (ja) * | 2013-12-05 | 2015-06-11 | 三菱電機株式会社 | レーザレーダ装置 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0411097D0 (en) * | 2004-05-19 | 2004-06-23 | Qinetiq Ltd | Laser radar device and method |
US7260507B2 (en) * | 2005-09-09 | 2007-08-21 | Itt Manufacturing Enterprises, Inc. | Method for improving the performance accuracy in differential absorption lidar for oil and gas pipeline leak detection and quantification |
CN102288973B (zh) * | 2011-06-28 | 2013-03-20 | 华中科技大学 | 基于频谱特征信号的多气象参数同步测量方法及激光雷达 |
CN103946716B (zh) * | 2011-12-21 | 2016-05-11 | 三菱电机株式会社 | 激光雷达装置 |
WO2014045627A1 (ja) * | 2012-09-24 | 2014-03-27 | 三菱電機株式会社 | コヒーレントライダ装置及びレーザレーダ装置 |
CN103852595B (zh) * | 2012-12-07 | 2016-06-29 | 中国电子科技集团公司第三研究所 | 一种低空风探测的方法 |
US8541744B1 (en) | 2013-03-09 | 2013-09-24 | Jian Liu | Method and apparatus for wavelength locking free optical frequency comb based differential absorption Lidar |
CN103472255B (zh) * | 2013-09-16 | 2017-01-18 | 南京牧镭激光科技有限公司 | 全光纤多普勒相干激光雷达风速测量装置 |
WO2015087380A1 (ja) * | 2013-12-09 | 2015-06-18 | 三菱電機株式会社 | レーザレーダ装置 |
-
2017
- 2017-01-24 US US16/471,430 patent/US11543524B2/en active Active
- 2017-01-24 EP EP17894031.8A patent/EP3564652B1/en active Active
- 2017-01-24 CN CN201780083987.XA patent/CN110234979B/zh active Active
- 2017-01-24 JP JP2018563960A patent/JP6505343B2/ja active Active
- 2017-01-24 WO PCT/JP2017/002282 patent/WO2018138766A1/ja unknown
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1183629A (ja) * | 1997-09-02 | 1999-03-26 | Sumitomo Heavy Ind Ltd | 差分吸収ライダ計測装置 |
JP2009503486A (ja) * | 2005-07-29 | 2009-01-29 | キネテイツク・リミテツド | レーザ測定のデバイスおよび方法 |
JP2013083467A (ja) * | 2011-10-06 | 2013-05-09 | Japan Aerospace Exploration Agency | 光学式遠隔気流計測装置の有色ノイズ低減方法とその装置 |
WO2014024508A1 (ja) * | 2012-08-08 | 2014-02-13 | 三菱電機株式会社 | レーダ装置 |
WO2014041852A1 (ja) * | 2012-09-14 | 2014-03-20 | 三菱電機株式会社 | レーザーレーダ装置及び測定対象物の速度算出方法 |
JP2015108557A (ja) * | 2013-12-05 | 2015-06-11 | 三菱電機株式会社 | レーザレーダ装置 |
Also Published As
Publication number | Publication date |
---|---|
CN110234979A (zh) | 2019-09-13 |
CN110234979B (zh) | 2022-01-11 |
US11543524B2 (en) | 2023-01-03 |
JPWO2018138766A1 (ja) | 2019-04-25 |
EP3564652A4 (en) | 2020-01-01 |
US20200309950A1 (en) | 2020-10-01 |
EP3564652A1 (en) | 2019-11-06 |
EP3564652B1 (en) | 2021-09-01 |
WO2018138766A1 (ja) | 2018-08-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6505343B2 (ja) | レーザレーダ装置 | |
CN110073240B (zh) | 激光雷达装置 | |
US11327176B2 (en) | Laser radar device | |
JP6301040B2 (ja) | コヒーレントライダー装置 | |
WO2017187510A1 (ja) | 距離計測装置、距離計測方法、及び形状計測装置 | |
JP6177338B2 (ja) | 風計測ライダ装置 | |
CN109891263B (zh) | 激光雷达装置 | |
US11506521B2 (en) | Optical fiber characteristics measurement apparatus and optical fiber characteristics measurement method | |
EP3605140B1 (en) | Laser radar device | |
EP2805179B1 (en) | Optical system for range finding | |
JP2023539144A (ja) | Lidarシステムにおけるノイズ校正およびターゲット検出 | |
JP2007333518A (ja) | 差分吸収ライダ装置 | |
EP2801163A1 (en) | Optical phased array | |
US20230032919A1 (en) | Techniques for time domain automatic gain control for a frequency modulated continuous wave (fmcw) light detection and ranging (lidar) system | |
JP2015524621A (ja) | 調節可能レーザを動的に掃引するためのシステムおよび方法 | |
JP2020079789A (ja) | 複数の並置された伝送経路を有するレーダーの伝送位相偏移の決定 | |
JP6847325B1 (ja) | コヒーレントライダ装置 | |
JP5941877B2 (ja) | 光パルス試験装置及び光パルス試験方法 | |
CN116601529A (zh) | 激光雷达装置和风计测方法 | |
JP7130175B2 (ja) | 光空間通信装置 | |
JP7515754B1 (ja) | レーザレーダ装置および信号処理方法 | |
WO2023228805A1 (ja) | 周波数掃引特性測定装置、LiDAR装置および周波数掃引特性測定方法 | |
WO2017009939A1 (ja) | レーザレーダ装置 | |
WO2022175998A1 (ja) | 光学装置及び信号処理方法 | |
WO2019181696A1 (ja) | 測距装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20190130 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20190130 |
|
A871 | Explanation of circumstances concerning accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A871 Effective date: 20190130 |
|
A975 | Report on accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A971005 Effective date: 20190218 |
|
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: 20190226 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20190326 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6505343 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |