(57)【特許請求の範囲】
【請求項１】
空間に存在している複数の観測対象におけるそれぞれの移動速度として、ライダ装置に対するそれぞれの観測対象の相対速度を算出する信号処理装置であって、
光源から出力されたレーザ光から、光周波数が互いに異なる複数のパルス光が生成されて、それぞれのパルス光が空間に放射され、それぞれの観測対象によって散乱された後のそれぞれのパルス光が散乱光として受信されて、それぞれの散乱光と前記レーザ光との合波光が検波されており、
それぞれの観測対象の移動に伴うドップラー周波数を算出するドップラー周波数算出部と、
前記ドップラー周波数算出部により算出されたそれぞれのドップラー周波数から、それぞれの観測対象の相対速度を算出する速度算出部と、を備え、
前記ドップラー周波数算出部は、
生成された複数のパルス光の中の、１つのパルス光が有している光周波数と残りのパルス光が有している光周波数との周波数差に基づいて、それぞれの合波光が有している光周波数を補正する光周波数補正部と、
前記光周波数補正部による補正後のそれぞれの合波光が有している光周波数と、生成されたそれぞれのパルス光が有している光周波数とから、それぞれの散乱光が有している光周波数に含まれているドップラー周波数を算出する周波数算出処理部と、を含む、
信号処理装置。
【請求項２】
前記ドップラー周波数算出部は、
前記光周波数補正部による光周波数補正後の複数の合波光の周波数スペクトルを積算するスペクトル積算部を備えており、
前記周波数算出処理部は、
前記スペクトル積算部による積算後の周波数スペクトルから、それぞれの合波光が有している光周波数を特定し、それぞれの合波光が有している光周波数と、生成されたそれぞれのパルス光が有している光周波数とから、それぞれの散乱光が有している光周波数に含まれているドップラー周波数を算出することを特徴とする請求項１記載の信号処理装置。
【請求項３】
空間に存在している複数の観測対象におけるそれぞれの移動速度として、ライダ装置に対するそれぞれの観測対象の相対速度を算出する信号処理方法であって、
光源から出力されたレーザ光から、光周波数が互いに異なる複数のパルス光が生成されて、それぞれのパルス光が空間に放射され、それぞれの観測対象によって散乱された後のそれぞれのパルス光が散乱光として受信されて、それぞれの散乱光と前記レーザ光との合波光が検波されており、
ドップラー周波数算出部の光周波数補正部が、生成された複数のパルス光の中の、１つのパルス光が有している光周波数と残りのパルス光が有している光周波数との周波数差に基づいて、それぞれの合波光が有している光周波数を補正し、
ドップラー周波数算出部の周波数算出処理部が、前記光周波数補正部による補正後のそれぞれの合波光が有している光周波数と、生成されたそれぞれのパルス光が有している光周波数とから、それぞれの散乱光が有している光周波数に含まれているドップラー周波数を算出し、
速度算出部が、前記ドップラー周波数算出部により算出されたそれぞれのドップラー周波数から、それぞれの観測対象の相対速度を算出する、
信号処理方法。
【請求項４】
空間に存在している複数の観測対象におけるそれぞれの移動速度として、ライダ装置に対するそれぞれの観測対象の相対速度を算出するライダ装置であって、
レーザ光を出力する光源と、
前記光源から出力されたレーザ光から、光周波数が互いに異なる複数のパルス光を生成するパルス変調部と、
前記パルス変調部により生成されたそれぞれのパルス光を空間に放射し、それぞれの観測対象によって散乱された後のそれぞれのパルス光を散乱光として受信する送受信部と、
前記送受信部により受信されたそれぞれの散乱光と前記光源から出力されたレーザ光との合波光を検波し、それぞれの合波光の検波信号を出力する光検波部と、
それぞれの観測対象の移動に伴うドップラー周波数を算出するドップラー周波数算出部と、
前記ドップラー周波数算出部により算出されたそれぞれのドップラー周波数から、それぞれの観測対象の相対速度を算出する速度算出部と、を備え、
前記ドップラー周波数算出部は、
生成された複数のパルス光の中の、１つのパルス光が有している光周波数と残りのパルス光が有している光周波数との周波数差に基づいて、それぞれの合波光が有している光周波数を補正する光周波数補正部と、
前記光周波数補正部による補正後のそれぞれの合波光が有している光周波数と、生成されたそれぞれのパルス光が有している光周波数とから、それぞれの散乱光が有している光周波数に含まれているドップラー周波数を算出する周波数算出処理部と、を含む、
ライダ装置。
【請求項５】
前記送受信部から放射されるパルス光の放射方向を切り替えるスキャナと、
前記スキャナによる放射方向の切替速度を制御する切替速度制御部と
を備えたことを特徴とする請求項４記載のライダ装置。
【請求項６】
前記切替速度制御部は、
前記送受信部から放射されるパルス光がアイセーフ条件を満足するように、前記送受信部から放射されるパルス光の繰り返し周波数に応じて前記切替速度を制御することを特徴とする請求項５記載のライダ装置。
【請求項７】
前記観測対象が気体であり、
前記光源は、
前記気体の吸収波長帯に含まれる波長を有する第１のレーザ光と、前記第１のレーザ光よりも前記気体による吸収率が低い第２のレーザ光とを順番に前記パルス変調部に出力することを特徴とする請求項４記載のライダ装置。
【請求項８】
前記光源から第１のレーザ光が出力されたときに、前記光検波部から出力されたそれぞれの検波信号の周波数を解析し、前記光源から第２のレーザ光が出力されたときに、前記光検波部から出力されたそれぞれの検波信号の周波数を解析し、それぞれの周波数の解析結果から、前記気体の密度を算出する密度算出部を備えたことを特徴とする請求項７記載のライダ装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、信号処理装置、信号処理方法及びライダ装置に関するものである。
【背景技術】
【０００２】
ライダ装置の中には、複数の観測対象におけるそれぞれの移動速度を算出するライダ装置（以下「従来のライダ装置」という）がある。
従来のライダ装置は、レーザ光を出力する光源と、当該レーザ光からパルス光を生成するパルス光生成部と、当該パルス光を空間に繰り返し放射し、それぞれの観測対象によって散乱された後のそれぞれのパルス光を散乱光として受信する送受信部と、それぞれの散乱光と当該レーザ光との合波光を検波する光合波部とを備えている。また、従来のライダ装置は、光合波部により検波されたそれぞれの合波光の光周波数から、それぞれの観測対象の移動に伴うドップラー周波数を算出し、それぞれのドップラー周波数からそれぞれの観測対象の移動速度を算出する信号処理部を備えている。信号処理部により算出される移動速度は、レーザ光の放射方向の移動速度である。【０００３】
ところで、パルス光を放射してから、観測対象によって散乱された後のパルス光である散乱光を受信するまでの時間に基づいて、測距装置から観測対象までの距離を算出する測距装置が特許文献１に開示されている。
当該測距装置は、複数のパルス光を生成する生成部と、生成部により生成されたそれぞれのパルス光を空間に放射する送信部と、観測対象によって散乱された後のそれぞれのパルス光を散乱光として受信する受信部とを備えている。また、当該測距装置は、送信部からパルス光が放射されてから、受信部により散乱光が受信されるまでの時間に基づいて、観測対象までの距離を算出する距離算出部を備えている。観測対象までの距離が長い場合、又は、パルス光の放射周期が短い場合、送信部が、放射順がＮ（Ｎは、１以上の整数）番目のパルス光を放射してから、受信部が、Ｎ番目のパルス光に対応する散乱光を受信する前に、送信部が、放射順が（Ｎ＋１）番目のパルス光を放射してしまうことがある。このような場合でも、受信部により受信された散乱光が、何番目のパルス光に対応する散乱光であるのかが分かるようにするために、生成部が、光周波数が互いに異なる複数のパルス光を生成し、それぞれのパルス光を送信部に出力している。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】国際公開第２０２０−０７９７７６号
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ライダ装置の遠方に存在している観測対象（以下「遠方の観測対象」という）の移動速度と、当該ライダ装置の近傍に存在している観測対象（以下「近傍の観測対象」という）の移動速度とが同じ速度であることがある。遠方の観測対象の移動速度と近傍の観測対象の移動速度と同じ速度であるときに、遠方の観測対象からの散乱光（以下「遠方散乱光」という）と、近傍の観測対象からの散乱光（以下「近傍散乱光」という）とが、互いに重なり合った状態で、当該ライダ装置に到達することがある。
従来のライダ装置では、送受信部が、互いに重なり合っている、遠方散乱光及び近傍散乱光のそれぞれを受信したとき、遠方散乱光が有している光周波数と近傍散乱光が有している光周波数とが同じ周波数であるため、信号処理部が、遠方散乱光と近傍散乱光とを識別することができない。したがって、当該信号処理部は、近傍の観測対象及び遠方の観測対象のうちのいずれか一方の観測対象の存在を検出して、当該観測対象の移動速度を算出できても、他方の観測対象の存在を検出できないという課題があった。
特許文献１に開示されている測距装置は、光周波数が互いに異なる複数のパルス光を生成する生成部を備えている。しかし、当該測距装置の距離算出部は、遠方散乱光と近傍散乱光とが互いに重なり合っているときに、遠方散乱光と近傍散乱光とを識別するものではない。したがって、仮に当該測距装置を従来のライダ装置に適用したとしても、上記課題を解決することができない。【０００６】
本開示は、上記のような課題を解決するためになされたもので、遠方散乱光と近傍散乱光とが、互いに重なり合った状態であっても、遠方の観測対象の移動速度と近傍の観測対象の移動速度とが同じ速度であるか否かにかかわらず、それぞれの観測対象の移動速度を算出することができる信号処理装置及び信号処理方法を得ることを目的とする。
【課題を解決するための手段】
【０００７】
本開示に係る信号処理装置は、空間に存在している複数の観測対象におけるそれぞれの移動速度として、ライダ装置に対するそれぞれの観測対象の相対速度を算出するものである。
光源から出力されたレーザ光から、光周波数が互いに異なる複数のパルス光が生成されて、それぞれのパルス光が空間に放射され、それぞれの観測対象によって散乱された後のそれぞれのパルス光が散乱光として受信されて、それぞれの散乱光とレーザ光との合波光が検波されている。当該信号処理装置は、それぞれの観測対象の移動に伴うドップラー周波数を算出するドップラー周波数算出部と、ドップラー周波数算出部により算出されたそれぞれのドップラー周波数から、それぞれの観測対象の相対速度を算出する速度算出部とを備える。当該ドップラー周波数算出部は、生成された複数のパルス光の中の、１つのパルス光が有している光周波数と残りのパルス光が有している光周波数との周波数差に基づいて、それぞれの合波光が有している光周波数を補正する光周波数補正部と、光周波数補正部による補正後のそれぞれの合波光が有している光周波数と、生成されたそれぞれのパルス光が有している光周波数とから、それぞれの散乱光が有している光周波数に含まれているドップラー周波数を算出する周波数算出処理部と、を含むものである。
【発明の効果】
【０００８】
本開示によれば、遠方散乱光と近傍散乱光とが、互いに重なり合った状態であっても、遠方の観測対象の移動速度と近傍の観測対象の移動速度とが同じ速度であるか否かにかかわらず、それぞれの観測対象の移動速度を算出することができる。
【図面の簡単な説明】
【０００９】
【図１】実施の形態１に係る信号処理装置１５を含むライダ装置を示す構成図である。
【図２】トリガ生成部４の内部を示す構成図である。
【図３】実施の形態１に係る信号処理装置１５を示す構成図である。
【図４】実施の形態１に係る信号処理装置１５のハードウェアを示すハードウェア構成図である。
【図５】信号処理装置１５が、ソフトウェア又はファームウェア等によって実現される場合のコンピュータのハードウェア構成図である。
【図６】図１に示すライダ装置の処理手順を示すフローチャートである。
【図７】信号処理装置１５の処理手順である信号処理方法を示すフローチャートである。
【図８】ライダ装置から放射されたパルス光Ｐ１と、第１の観測対象からの散乱光Ｒ１と、第２の観測対象からの散乱光Ｒ２とを示す説明図である。
【図９】ライダ装置から放射されたパルス光Ｐ１，Ｐ２の放射時刻Ｔ１，Ｔ２と、第１の観測対象からの散乱光Ｒ１，Ｒ３の受信時刻Ｔ１’，Ｔ３’と、第２の観測対象からの散乱光Ｒ２，Ｒ４の受信時刻Ｔ２’，Ｔ４’とを示す説明図である。
【図１０】図１０Ａは、レンジビン信号（７）が有しているピークスペクトルを示す説明図、図１０Ｂは、パルス光Ｐ１が有している光周波数とパルス光Ｐ２が有している光周波数とが同一周波数である場合の、レンジビン信号（７）が有しているピークスペクトルを示す説明図である。
【図１１】合波光Ｃ１に係る散乱光Ｒ１が存在しているレンジビンと、合波光Ｃ２，Ｃ３に係る散乱光Ｒ２，Ｒ３が存在しているレンジビンと、合波光Ｃ４に係る散乱光Ｒ４が存在しているレンジビンとを示す説明図である。
【図１２】レンジビン（ｎ）の距離とＳＮＲとの対応関係を示す距離特性（ａ−ｓｃｏｐｅ）の説明図である。
【図１３】実施の形態２に係るライダ装置のパルス変調器５を示す構成図である。
【図１４】実施の形態３に係るライダ装置を示す構成図である。
【図１５】実施の形態４に係る信号処理装置１５を含むライダ装置を示す構成図である。
【図１６】実施の形態４に係る信号処理装置１５を示す構成図である。
【図１７】実施の形態４に係る信号処理装置１５のハードウェアを示すハードウェア構成図である。
【図１８】気体の吸収波長帯、光源７１から出力された第１のレーザ光が有する波長及び光源７１から出力された第２のレーザ光が有する波長を示す説明図である。
【発明を実施するための形態】
【００１０】
以下、本開示をより詳細に説明するために、本開示を実施するための形態について、添付の図面に従って説明する。
【００１１】
実施の形態１．
図１は、実施の形態１に係る信号処理装置１５を含むライダ装置を示す構成図である。
図１に示すライダ装置は、光源１、パルス変調部２、送受信部６、光検波部１１及び信号処理装置１５を備えている。
ライダ装置は、空間に存在している複数の観測対象におけるそれぞれの移動速度として、ライダ装置に対するそれぞれの観測対象の相対速度を算出する。観測対象は、固体、液体、又は、気体である。
図１に示すライダ装置では、説明の簡単化のため、観測対象の数が２つであるものとして、一方の観測対象が第１の観測対象、他方の観測対象が第２の観測対象であるとする。ただし、これは一例に過ぎず、観測対象の数が３つ以上であってもよい。
【００１２】
光源１は、例えば、単一周波数のレーザ光を発光するレーザであって、発光スペクトルの線幅が数ＭＨｚ以下の、半導体レーザ、ファイバレーザ、又は、固体レーザによって実現される。あるいは、光源１は、半導体レーザ、ファイバレーザ及び固体レーザのうち、１つ以上のレーザの組み合わせによって実現される。
光源１は、連続光であるレーザ光をパルス変調部２に出力する。光源１から出力されるレーザ光が有している光周波数は、ｆ０である。
【００１３】
パルス変調部２は、光分割部３、トリガ生成部４及びパルス変調器５を備えている。
パルス変調部２は、光源１から出力されたレーザ光から、光周波数が互いに異なる複数のパルス光を生成する。
パルス変調部２により生成される複数のパルス光が、例えば、２つのパルス光Ｐ１，Ｐ２であれば、パルス光Ｐ１が有している光周波数は、ｆ０＋ｆＩＦ１であり、パルス光Ｐ２が有している光周波数は、ｆ０＋ｆＩＦ２である。例えば、ｆＩＦ２＞ｆＩＦ１である。ｆＩＦ２とｆＩＦ１との周波数差Δｆ（＝ｆＩＦ２−ｆＩＦ１）は、第１の観測対象及び第２の観測対象のそれぞれが、想定される最高速度で移動した場合に生じるドップラー周波数の絶対値の２倍よりも大きい。周波数差Δｆが、ドップラー周波数の絶対値の２倍よりも大きければ、第１の観測対象及び第２の観測対象のそれぞれの移動速度が、想定の範囲内の速度である限り、後述する複数の散乱光が有する光周波数は、互いに異なる周波数となる。
パルス変調部２により生成される複数のパルス光が、Ｍ（Ｍは、３つ以上の整数）個のパルス光Ｐ１，Ｐ２，・・・，ＰＭであれば、パルス光Ｐｍ（ｍ＝１，２，・・・，Ｍ）が有している光周波数は、ｆ０＋ｆＩＦｍである。例えば、ｆＩＦＭ＞ｆＩＦ（Ｍ−１）＞・・・＞ｆＩＦ１である。
図１に示すライダ装置では、説明の簡単化のため、パルス変調部２が、２つのパルス光Ｐ１，Ｐ２を生成するものとする。
【００１４】
光分割部３は、例えば、ビームスプリッタ、ファイバ型カプラ、又は、ハーフミラーによって実現される。
光分割部３は、光源１から出力されたレーザ光を所定の比率で２分配する。所定の比率としては、例えば、パルス変調器５側が２、後述する光合波部１２が１の比率である。
光分割部３は、分配後の一方のレーザ光をパルス変調器５に出力し、分配後の他方のレーザ光を参照光として光合波部１２に出力する。
【００１５】
トリガ生成部４は、例えば、パルスジェネレータ、ファンクションジェネレータ、又は、ＦＰＧＡ（Ｆｉｅｌｄ−Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）によって実現される。
トリガ生成部４は、図２に示すように、パルス信号生成部４ａ、基準信号生成部４ｂ、第１のシフト信号生成部４ｃ、第２のシフト信号生成部４ｄ及びスイッチ４ｅを備えている。
トリガ生成部４は、パルス信号をパルス変調器５、後述するアナログデジタル変換器（以下「Ａ／Ｄ変換部」という）１４及び信号処理装置１５のそれぞれに出力する。
また、トリガ生成部４は、変調周波数ｆＩＦ１を示す第１の周波数シフト信号と、変調周波数ｆＩＦ２を示す第２の周波数シフト信号とをパルス変調器５及び信号処理装置１５のそれぞれに出力する。
【００１６】
図２は、トリガ生成部４の内部を示す構成図である。
パルス信号生成部４ａは、パルス幅ΔＴのパルス信号を生成する。
パルス信号生成部４ａは、パルス幅ΔＴのパルス信号をパルス変調器５、Ａ／Ｄ変換部１４及び信号処理装置１５のそれぞれに、周期Ｔｒｅｐで繰り返し出力する。
基準信号生成部４ｂは、周波数ｆの電気信号である基準信号を生成する。
基準信号生成部４ｂは、基準信号を第１のシフト信号生成部４ｃ及び第２のシフト信号生成部４dのそれぞれに出力する。
【００１７】
第１のシフト信号生成部４ｃは、基準信号生成部４ｂから出力された基準信号を取得し、基準信号から、変調周波数ｆＩＦ１を示す第１の周波数シフト信号を生成する。
第１のシフト信号生成部４ｃは、第１の周波数シフト信号をスイッチ４ｅに出力する。
第２のシフト信号生成部４ｄは、基準信号生成部４ｂから出力された基準信号を取得し、基準信号から、変調周波数ｆＩＦ２を示す第２の周波数シフト信号を生成する。
第２のシフト信号生成部４ｄは、第２の周波数シフト信号をスイッチ４ｅに出力する。
スイッチ４ｅは、第１のシフト信号生成部４ｃから出力された第１の周波数シフト信号と、第２のシフト信号生成部４ｄから出力された第２の周波数シフト信号とを順番にパルス変調器５及び信号処理装置１５のそれぞれに出力する。
【００１８】
パルス変調器５は、例えば、音響光学素子を用いている変調素子又はニオブ酸リチウム結晶を用いている変調素子と、半導体光アンプ等の光増幅器とによって実現される。
パルス変調器５は、パルス信号生成部４ａからパルス幅ΔＴのパルス信号を受ける毎に、当該パルス信号に従って、光分割部３から出力されたレーザ光をパルス変調する。即ち、パルス変調器５は、連続光であるレーザ光を、パルス幅ΔＴを有するパルス光に変換する。
また、パルス変調器５は、スイッチ４ｅから第１の周波数シフト信号が出力されれば、パルス幅ΔＴを有するパルス光の光周波数ｆ０を光周波数ｆ０＋ｆＩＦ１にシフトすることによって、光周波数ｆ０＋ｆＩＦ１を有するパルス光Ｐ１を生成する。
パルス変調器５は、パルス光Ｐ１を送信側光学系７に出力する。
パルス変調器５は、スイッチ４ｅから第２の周波数シフト信号が出力されれば、パルス幅ΔＴを有するパルス光の光周波数ｆ０を光周波数ｆ０＋ｆＩＦ２にシフトすることによって、光周波数ｆ０＋ｆＩＦ２を有するパルス光Ｐ２を生成する。
パルス変調器５は、パルス光Ｐ２を送信側光学系７に出力する。
【００１９】
送受信部６は、送信側光学系７、送受分離部８、テレスコープ９及び受信側光学系１０を備えている。
送受信部６は、パルス変調部２により生成されたそれぞれのパルス光Ｐ１，Ｐ２を空間に放射する。
送受信部６は、第１の観測対象によって散乱された後のパルス光Ｐ１を散乱光Ｒ１として受信し、第２の観測対象によって散乱された後のパルス光Ｐ１を散乱光Ｒ２として受信する。
送受信部６は、第１の観測対象によって散乱された後のパルス光Ｐ２を散乱光Ｒ３として受信し、第２の観測対象によって散乱された後のパルス光Ｐ２を散乱光Ｒ４として受信する。
第１の観測対象がパルス光Ｐ１，Ｐ２の放射方向に移動していれば、散乱光Ｒ１，Ｒ３が有している光周波数は、第１の観測対象の移動に伴うドップラー周波数ｆｄｐ１を含んでいる。したがって、散乱光Ｒ１が有している光周波数は、ｆ０＋ｆＩＦ１＋ｆｄｐ１になり、散乱光Ｒ３が有している光周波数は、ｆ０＋ｆＩＦ２＋ｆｄｐ１になる。
第２の観測対象がパルス光Ｐ１，Ｐ２の放射方向に移動していれば、散乱光Ｒ２，Ｒ４が有している光周波数は、第２の観測対象の移動に伴うドップラー周波数ｆｄｐ２を含んでいる。したがって、散乱光Ｒ２が有している光周波数は、ｆ０＋ｆＩＦ１＋ｆｄｐ２になり、散乱光Ｒ４が有している光周波数は、ｆ０＋ｆＩＦ２＋ｆｄｐ２になる。
【００２０】
送信側光学系７は、パルス変調器５から出力されたそれぞれのパルス光Ｐ１，Ｐ２を整形し、整形後のそれぞれのパルス光Ｐ１，Ｐ２を送受分離部８に出力する。パルス光の整形としては、パルス光におけるビーム径の整形のほか、パルス光における広がり角の整形が該当する。
送受分離部８は、例えば、偏光ビームスプリッタと波長板とによって実現される。
送受分離部８は、送信側光学系７から出力された整形後のそれぞれのパルス光Ｐ１，Ｐ２の光軸上に設置されている。
送受分離部８は、送信側光学系７から出力された整形後のそれぞれのパルス光Ｐ１，Ｐ２をテレスコープ９に出力し、テレスコープ９によって集光されたそれぞれの散乱光Ｒ１，Ｒ２，Ｒ３，Ｒ４を受信側光学系１０に出力する。
【００２１】
テレスコープ９は、例えば、複数の屈折レンズ、又は、複数のミラーによって実現される。
テレスコープ９は、送受分離部８から出力された整形後のそれぞれのパルス光Ｐ１，Ｐ２を空間に放射する。
テレスコープ９は、第１の観測対象によって散乱された後のパルス光Ｐ１を散乱光Ｒ１として集光し、第２の観測対象によって散乱された後のパルス光Ｐ１を散乱光Ｒ２として集光する。
テレスコープ９は、第１の観測対象によって散乱された後のパルス光Ｐ２を散乱光Ｒ３として集光し、第２の観測対象によって散乱された後のパルス光Ｐ２を散乱光Ｒ４として集光する。
テレスコープ９は、それぞれの散乱光Ｒ１，Ｒ２，Ｒ３，Ｒ４を送受分離部８に出力する。
受信側光学系１０は、送受分離部８から出力されたそれぞれの散乱光Ｒ１，Ｒ２，Ｒ３，Ｒ４の光軸と、光合波部１２の光軸とが一致するように設置されている。
受信側光学系１０は、送受分離部８から出力されたそれぞれの散乱光Ｒ１，Ｒ２，Ｒ３，Ｒ４を整形し、整形後のそれぞれの散乱光Ｒ１，Ｒ２，Ｒ３，Ｒ４を光合波部１２に出力する。散乱光の整形としては、散乱光におけるビーム径の整形のほか、散乱光における広がり角の整形が該当する。
【００２２】
光検波部１１は、光合波部１２、受光部１３及びＡ／Ｄ変換部１４を備えている。
光検波部１１は、送受信部６により受信されたそれぞれの散乱光Ｒ１，Ｒ２，Ｒ３，Ｒ４と光源１から出力されたレーザ光である参照光との合波光Ｃ１，Ｃ２，Ｃ３，Ｃ４を検波する。
光検波部１１は、それぞれの合波光Ｃ１，Ｃ２，Ｃ３，Ｃ４の検波信号Ｄ１，Ｄ２，Ｄ３，Ｄ４を信号処理装置１５に出力する。
【００２３】
光合波部１２は、例えば、ビームスプリッタ、又は、ファイバ型カプラによって実現される。
光合波部１２は、受信側光学系１０から出力された整形後のそれぞれの散乱光Ｒ１，Ｒ２，Ｒ３，Ｒ４と光分割部３から出力された参照光との合波光Ｃ１，Ｃ２，Ｃ３，Ｃ４を検波する。
即ち、光合波部１２は、それぞれの散乱光Ｒ１，Ｒ２，Ｒ３，Ｒ４と光源１から出力された参照光とを混合することによって、合波光Ｃ１，Ｃ２，Ｃ３，Ｃ４をヘテロダイン検波する。合波光Ｃ１が有している光周波数は、ｆＩＦ１＋ｆｄｐ１であり、合波光Ｃ２が有している光周波数は、ｆＩＦ１＋ｆｄｐ２である。合波光Ｃ３が有している光周波数は、ｆＩＦ２＋ｆｄｐ１であり、合波光Ｃ４が有している光周波数は、ｆＩＦ２＋ｆｄｐ２である。
光合波部１２は、それぞれの合波光Ｃ１，Ｃ２，Ｃ３，Ｃ４を受光部１３に出力する。
【００２４】
受光部１３は、例えば、フォトダイオードによって実現される。
受光部１３は、光合波部１２から出力されたそれぞれの合波光Ｃ１，Ｃ２，Ｃ３，Ｃ４を電気信号に変換する。
受光部１３は、それぞれの電気信号をＡ／Ｄ変換部１４に出力する。
受光部１３は、受信側光学系１０から、散乱光Ｒ１，Ｒ２，Ｒ３，Ｒ４のいずれも出力されていない期間中は、電圧がほぼ０の電気信号をＡ／Ｄ変換部１４に出力する。
【００２５】
Ａ／Ｄ変換部１４は、トリガ生成部４のパルス信号生成部４ａからパルス幅ΔＴのパルス信号が出力されている期間中、受光部１３から出力された電気信号を、アナログ信号からデジタル信号Ｄｉｇ（ｔ）に変換する処理を行う。ｔは、サンプリング時刻を示す変数である。
Ａ／Ｄ変換部１４は、デジタル信号Ｄｉｇ（ｔ）を信号処理装置１５に出力する。受信側光学系１０から散乱光Ｒ１が出力されている期間中、デジタル信号Ｄｉｇ（ｔ）は、合波光Ｃ１の検波信号Ｄ１を示し、受信側光学系１０から散乱光Ｒ２が出力されている期間中、デジタル信号Ｄｉｇ（ｔ）は、合波光Ｃ２の検波信号Ｄ２を示している。受信側光学系１０から散乱光Ｒ３が出力されている期間中、デジタル信号Ｄｉｇ（ｔ）は、合波光Ｃ３の検波信号Ｄ３を示し、受信側光学系１０から散乱光Ｒ４が出力されている期間中、デジタル信号Ｄｉｇ（ｔ）は、合波光Ｃ４の検波信号Ｄ４を示している。
散乱光Ｒ１，Ｒ２，Ｒ３，Ｒ４のいずれも出力されていない期間中、デジタル信号Ｄｉｇ（ｔ）は、ほぼ０の値を示している。
【００２６】
図３は、実施の形態１に係る信号処理装置１５を示す構成図である。
信号処理装置１５は、ドップラー周波数算出部１６、速度算出部１７、ＳＮＲ（Ｓｉｇｎａｌ  ｔｏ  Ｎｏｉｓｅ  Ｒａｔｉｏ）算出部３０及び距離特性算出部３１を備えている。
信号処理装置１５は、Ａ／Ｄ変換部１４から出力されたデジタル信号Ｄｉｇ（ｔ）に基づいて、第１の観測対象及び第２の観測対象におけるそれぞれの移動速度として、ライダ装置に対する第１の観測対象の相対速度Ｖ１と、ライダ装置に対する第２の観測対象の相対速度Ｖ２とを算出する。
図４は、実施の形態１に係る信号処理装置１５のハードウェアを示すハードウェア構成図である。
【００２７】
ドップラー周波数算出部１６は、例えば、図４に示すドップラー周波数算出回路４１によって実現される。
ドップラー周波数算出部１６は、光周波数補正部２１及び周波数算出処理部２８を備えている。
ドップラー周波数算出部１６は、パルス変調部２により生成されたそれぞれのパルス光Ｐ１，Ｐ２が有している光周波数ｆ０＋ｆＩＦ１，ｆ０＋ｆＩＦ２と、光検波部１１から出力されたそれぞれの合波光Ｃ１，Ｃ２，Ｃ３，Ｃ４の検波信号Ｄ１，Ｄ２，Ｄ３，Ｄ４とを取得する。
ドップラー周波数算出部１６は、それぞれのパルス光Ｐ１，Ｐ２が有している光周波数ｆ０＋ｆＩＦ１，ｆ０＋ｆＩＦ２と、それぞれの検波信号Ｄ１，Ｄ２，Ｄ３，Ｄ４とから、それぞれの散乱光Ｒ１，Ｒ２，Ｒ３，Ｒ４が有している光周波数に含まれている、それぞれの観測対象の移動に伴うドップラー周波数を算出する。
即ち、ドップラー周波数算出部１６は、複数のパルス光Ｐ１，Ｐ２の中の、１つのパルス光Ｐ１が有している光周波数ｆ０＋ｆＩＦ１と残りのパルス光Ｐ２が有している光周波数との周波数差ｆ０＋ｆＩＦ２との周波数差Δｆ（＝ｆＩＦ２−ｆＩＦ１）を算出する。
ドップラー周波数算出部１６は、周波数差Δｆと、デジタル信号Ｄｉｇ（ｔ）に含まれているそれぞれの検波信号Ｄ１，Ｄ３とから、第１の観測対象の相対速度Ｖ１に対応するドップラー周波数ｆｄｐ１として、それぞれの散乱光Ｒ１，Ｒ３が有している光周波数に含まれているドップラー周波数ｆｄｐ１を算出する。
ドップラー周波数算出部１６は、周波数差Δｆと、デジタル信号Ｄｉｇ（ｔ）に含まれているそれぞれの検波信号Ｄ２，Ｄ４とから、第２の観測対象の相対速度Ｖ２に対応するドップラー周波数ｆｄｐ２として、それぞれの散乱光Ｒ２，Ｒ４が有している光周波数に含まれているドップラー周波数ｆｄｐ２を算出する。
ドップラー周波数算出部１６は、ドップラー周波数ｆｄｐ１とドップラー周波数ｆｄｐ２とを周波数算出処理部２８に出力する。
【００２８】
光周波数補正部２１は、レンジビン分割部２２、周波数解析部２３、レンジ補正部２４、周波数補正処理部２５、スペクトル積算部２６及びピーク周波数検出部２７を備えている。
光周波数補正部２１は、トリガ生成部４から、第１の周波数シフト信号及び第２の周波数シフト信号のそれぞれを取得する。
光周波数補正部２１は、第２の周波数シフト信号が示す変調周波数ｆＩＦ２から、第１の周波数シフト信号が示す変調周波数ｆＩＦ１を減算することによって、周波数差Δｆを算出する。
光周波数補正部２１は、Ａ／Ｄ変換部１４から、それぞれの検波信号Ｄ１，Ｄ２，Ｄ３，Ｄ４を含んでいるデジタル信号Ｄｉｇ（ｔ）を取得する。
光周波数補正部２１は、周波数差Δｆに基づいて、合波光Ｃ３が有している光周波数ｆＩＦ２＋ｆｄｐ１を補正する。
光周波数補正部２１は、周波数差Δｆに基づいて、合波光Ｃ４が有している光周波数ｆＩＦ２＋ｆｄｐ２を補正する。
【００２９】
レンジビン分割部２２は、Ａ／Ｄ変換部１４から出力されたデジタル信号Ｄｉｇ（ｔ）を時間方向に分割する。デジタル信号Ｄｉｇ（ｔ）における時間方向の分割幅Δｔは、レンジビン幅Ｒｂｗに対応している。したがって、サンプリング時刻ｔ＝１に係る分割後のデジタル信号Ｄｉｇ（１）は、レンジビン（１）に対応し、サンプリング時刻ｔ＝２に係る分割後のデジタル信号Ｄｉｇ（２）は、レンジビン（２）に対応する。また、サンプリング時刻ｔ＝３に係る分割後のデジタル信号Ｄｉｇ（３）は、レンジビン（３）に対応する。
レンジビン分割部２２は、分割後のそれぞれのデジタル信号Ｄｉｇ（ｔ）をレンジビン信号（ｎ）として周波数解析部２３に出力する。ｎは、レンジビンを示す変数であり、ｎ＝１，２，３，・・・である。
【００３０】
周波数解析部２３は、レンジビン分割部２２から出力されたそれぞれのレンジビン信号（ｎ）に対するＦＦＴ（Ｆａｓｔ  Ｆｏｕｒｉｅｒ  Ｔｒａｎｓｆｏｒｍ）処理を実施することで、それぞれのレンジビン信号（ｎ）の周波数スペクトルＦＳ（ｎ）を算出する。
周波数解析部２３は、それぞれの周波数スペクトルＦＳ（ｎ）をレンジ補正部２４に出力する。
【００３１】
レンジ補正部２４は、周波数解析部２３から、それぞれの周波数スペクトルＦＳ（ｎ）を取得する。
レンジ補正部２４は、トリガ生成部４から、第１の周波数シフト信号と第２の周波数シフト信号とを取得する。
レンジ補正部２４は、周波数スペクトルＦＳ（ｎ）から、ピークスペクトルＳｐ１、ピークスペクトルＳｐ２、ピークスペクトルＳｐ３及びピークスペクトルＳｐ４のそれぞれを検出する。
また、レンジ補正部２４は、ピークスペクトルＳｐ１に対応するピーク周波数ｆｐ１を検出し、ピークスペクトルＳｐ２に対応するピーク周波数ｆｐ２を検出する。
レンジ補正部２４は、ピークスペクトルＳｐ３に対応するピーク周波数ｆｐ３を検出し、ピークスペクトルＳｐ４に対応するピーク周波数ｆｐ４を検出する。
【００３２】
レンジ補正部２４は、ピーク周波数ｆｐ１と第１の周波数シフト信号が示す変調周波数ｆＩＦ１との差分の絶対値｜Δｆ１−１｜と、ピーク周波数ｆｐ１と第２の周波数シフト信号が示す変調周波数ｆＩＦ２との差分の絶対値｜Δｆ１−２｜とを算出する。
レンジ補正部２４は、差分の絶対値｜Δｆ１−１｜が差分の絶対値｜Δｆ１−２｜以下であれば、ピークスペクトルＳｐ１を有する合波光が、パルス光Ｐ１に対応する合波光Ｃ１、又は、パルス光Ｐ１に対応する合波光Ｃ２であると判定する。
レンジ補正部２４は、差分の絶対値｜Δｆ１−１｜が差分の絶対値｜Δｆ１−２｜よりも大きければ、ピークスペクトルＳｐ１を有する合波光が、パルス光Ｐ２に対応する合波光Ｃ３、又は、パルス光Ｐ２に対応する合波光Ｃ４であると判定する。
レンジ補正部２４は、ピーク周波数ｆｐ２と変調周波数ｆＩＦ１との差分の絶対値｜Δｆ２−１｜と、ピーク周波数ｆｐ２と変調周波数ｆＩＦ２との差分の絶対値｜Δｆ２−２｜とを算出する。
レンジ補正部２４は、差分の絶対値｜Δｆ２−１｜が差分の絶対値｜Δｆ２−２｜以下であれば、ピークスペクトルＳｐ２を有する合波光が、パルス光Ｐ１に対応する合波光Ｃ１、又は、パルス光Ｐ１に対応する合波光Ｃ２であると判定する。
レンジ補正部２４は、差分の絶対値｜Δｆ２−１｜が差分の絶対値｜Δｆ２−２｜よりも大きければ、ピークスペクトルＳｐ２を有する合波光が、パルス光Ｐ２に対応する合波光Ｃ３、又は、パルス光Ｐ２に対応する合波光Ｃ４であると判定する。
【００３３】
レンジ補正部２４は、ピーク周波数ｆｐ３と変調周波数ｆＩＦ１との差分の絶対値｜Δｆ３−１｜と、ピーク周波数ｆｐ３と変調周波数ｆＩＦ２との差分の絶対値｜Δｆ３−２｜とを算出する。
レンジ補正部２４は、差分の絶対値｜Δｆ３−１｜が差分の絶対値｜Δｆ３−２｜以下であれば、ピークスペクトルＳｐ３を有する合波光が、パルス光Ｐ１に対応する合波光Ｃ１、又は、パルス光Ｐ１に対応する合波光Ｃ２であると判定する。
レンジ補正部２４は、差分の絶対値｜Δｆ３−１｜が差分の絶対値｜Δｆ３−２｜よりも大きければ、ピークスペクトルＳｐ３を有する合波光が、パルス光Ｐ２に対応する合波光Ｃ３、又は、パルス光Ｐ２に対応する合波光Ｃ４であると判定する。
レンジ補正部２４は、ピーク周波数ｆｐ４と変調周波数ｆＩＦ１との差分の絶対値｜Δｆ４−１｜と、ピーク周波数ｆｐ４と変調周波数ｆＩＦ２との差分の絶対値｜Δｆ４−２｜とを算出する。
レンジ補正部２４は、差分の絶対値｜Δｆ４−１｜が差分の絶対値｜Δｆ４−２｜以下であれば、ピークスペクトルＳｐ４を有する合波光が、パルス光Ｐ１に対応する合波光Ｃ１、又は、パルス光Ｐ１に対応する合波光Ｃ２であると判定する。
レンジ補正部２４は、差分の絶対値｜Δｆ４−１｜が差分の絶対値｜Δｆ４−２｜よりも大きければ、ピークスペクトルＳｐ４を有する合波光が、パルス光Ｐ２に対応する合波光Ｃ３、又は、パルス光Ｐ２に対応する合波光Ｃ４であると判定する。
【００３４】
ここでは、説明の便宜上、ピークスペクトルＳｐ１を有する合波光及びピークスペクトルＳｐ２を有する合波光のそれぞれが、パルス光Ｐ１に対応する合波光Ｃ１、又は、パルス光Ｐ１に対応する合波光Ｃ２であるものとする。
また、ピークスペクトルＳｐ３を有する合波光及びピークスペクトルＳｐ４を有する合波光のそれぞれが、パルス光Ｐ２に対応する合波光Ｃ３、又は、パルス光Ｐ２に対応する合波光Ｃ４であるものとする。
この場合、レンジ補正部２４は、ピークスペクトルＳｐ１のレンジビンが、ピークスペクトルＳｐ２のレンジビン以下であれば、ピークスペクトルＳｐ１を有する合波光が、パルス光Ｐ１に対応する合波光Ｃ１であり、ピークスペクトルＳｐ２を有する合波光が、パルス光Ｐ１に対応する合波光Ｃ２であると判定する。レンジ補正部２４は、ピークスペクトルＳｐ１のレンジビンが、ピークスペクトルＳｐ２のレンジビンよりも大きければ、ピークスペクトルＳｐ１を有する合波光が、パルス光Ｐ１に対応する合波光Ｃ２であり、ピークスペクトルＳｐ２を有する合波光が、パルス光Ｐ１に対応する合波光Ｃ１であると判定する。
また、レンジ補正部２４は、ピークスペクトルＳｐ３のレンジビンが、ピークスペクトルＳｐ４のレンジビン以下であれば、ピークスペクトルＳｐ３を有する合波光が、パルス光Ｐ２に対応する合波光Ｃ３であり、ピークスペクトルＳｐ４を有する合波光が、パルス光Ｐ２に対応する合波光Ｃ４であると判定する。レンジ補正部２４は、ピークスペクトルＳｐ３のレンジビンが、ピークスペクトルＳｐ４のレンジビンよりも大きければ、ピークスペクトルＳｐ３を有する合波光が、パルス光Ｐ２に対応する合波光Ｃ４であり、ピークスペクトルＳｐ４を有する合波光が、パルス光Ｐ２に対応する合波光Ｃ３であると判定する。
ここでは、説明の便宜上、ピークスペクトルＳｐ１を有する合波光が合波光Ｃ１であり、ピークスペクトルＳｐ２を有する合波光が合波光Ｃ２であるとする。また、ピークスペクトルＳｐ３を有する合波光が合波光Ｃ３であり、ピークスペクトルＳｐ４を有する合波光が合波光Ｃ４であるとする。
【００３５】
レンジ補正部２４は、合波光Ｃ１，Ｃ２がパルス光Ｐ１に対応する合波光であり、合波光Ｃ３，Ｃ４がパルス光Ｐ２に対応する合波光であると判定すると、後述する図１１に示すように、合波光Ｃ３に係る散乱光Ｒ３が存在しているレンジビン（７）と、合波光Ｃ４に係る散乱光Ｒ４が存在しているレンジビン（１３）とを補正する。
即ち、レンジ補正部２４は、合波光Ｃ３に係る散乱光Ｒ３が存在しているレンジビン（７）から周期Ｔｒｅｐ（＝６）を減算することによって、合波光Ｃ３に係る散乱光Ｒ３が存在しているレンジビン（７）をレンジビン（１）に補正する。
また、レンジ補正部２４は、合波光Ｃ４に係る散乱光Ｒ４が存在しているレンジビン（１３）から周期Ｔｒｅｐ（＝６）を減算することによって、合波光Ｃ４に係る散乱光Ｒ４が存在しているレンジビン（１３）をレンジビン（７）に補正する。
図１１の例では、パルス光Ｐ１の放射時刻Ｔ１とパルス光Ｐ２の放射時刻Ｔ２との時刻差がレンジビン（６）に相当しており、周期Ｔｒｅｐが６である。このため、合波光Ｃ３に係る散乱光Ｒ３が存在している補正後のレンジビンであるｎ’は、１（＝７−６）となり、合波光Ｃ４に係る散乱光Ｒ４が存在している補正後のレンジビンであるｎ’は、７（＝１３−６）となる。
レンジ補正部２４は、複数の周波数スペクトルＦＳ（１）〜ＦＳ（Ｎ）のうち、合波光Ｃ１，Ｃ２に係る散乱光Ｒ１，Ｒ２が存在しているレンジビン（ｎ）についての周波数スペクトルＦＳ（ｎ）を周波数補正処理部２５に出力する。
また、レンジ補正部２４は、複数の周波数スペクトルＦＳ（１）〜ＦＳ（Ｎ）のうち、合波光Ｃ３，Ｃ４に係る散乱光Ｒ３，Ｒ４が存在している補正後のレンジビン（ｎ’）についての周波数スペクトルＦＳ（ｎ’）を周波数補正処理部２５に出力する。
レンジ補正部２４は、合波光Ｃ１が有している光周波数ｆＩＦ１＋ｆｄｐ１（＝ｆｐ１）と、合波光Ｃ２が有している光周波数ｆＩＦ１＋ｆｄｐ２（＝ｆｐ２）と、合波光Ｃ３が有している光周波数ｆＩＦ２＋ｆｄｐ１（＝ｆｐ３）と、合波光Ｃ４が有している光周波数ｆＩＦ２＋ｆｄｐ２（＝ｆｐ４）とを周波数補正処理部２５に出力する。
レンジ補正部２４は、合波光Ｃ１に係る散乱光Ｒ１が存在しているレンジビンと、合波光Ｃ２に係る散乱光Ｒ２が存在しているレンジビンとを距離特性算出部３１に出力する。
【００３６】
周波数補正処理部２５は、トリガ生成部４から、第１の周波数シフト信号と第２の周波数シフト信号とを取得する。
周波数補正処理部２５は、第１の周波数シフト信号が示す変調周波数ｆＩＦ１と第２の周波数シフト信号が示す変調周波数ｆＩＦ２との周波数差Δｆ（＝ｆＩＦ２−ｆＩＦ１）を算出する。
周波数補正処理部２５は、レンジ補正部２４から、合波光Ｃ１，Ｃ２に係る散乱光Ｒ１，Ｒ２が存在しているレンジビン（ｎ）についての周波数スペクトルＦＳ（ｎ）と、合波光Ｃ３，Ｃ４に係る散乱光Ｒ３，Ｒ４が存在している補正後のレンジビン（ｎ’）についての周波数スペクトルＦＳ（ｎ’）とを取得する。
周波数補正処理部２５は、レンジ補正部２４から、合波光Ｃ１が有している光周波数ｆＩＦ１＋ｆｄｐ１と、合波光Ｃ２が有している光周波数ｆＩＦ１＋ｆｄｐ２と、合波光Ｃ３が有している光周波数ｆＩＦ２＋ｆｄｐ１と、レンジビン補正後の合波光Ｃ４が有している光周波数ｆＩＦ２＋ｆｄｐ２とを取得する。
周波数補正処理部２５は、レンジビン補正後の合波光Ｃ３が有している光周波数ｆＩＦ２＋ｆｄｐ１から周波数差Δｆを減算することによって、合波光Ｃ３が有している光周波数を補正する。合波光Ｃ３が有している補正後の光周波数は、ｆＩＦ１＋ｆｄｐ１であり、合波光Ｃ１が有している光周波数ｆＩＦ１＋ｆｄｐ１と同じ周波数である。
周波数補正処理部２５は、レンジビン補正後の合波光Ｃ４が有している光周波数ｆＩＦ２＋ｆｄｐ２から周波数差Δｆを減算することによって、合波光Ｃ４が有している光周波数を補正する。合波光Ｃ４が有している補正後の光周波数は、ｆＩＦ１＋ｆｄｐ２であり、合波光Ｃ２が有している光周波数ｆＩＦ１＋ｆｄｐ２と同じ周波数である。
図３に示す信号処理装置１５では、周波数補正処理部２５が、レンジビン補正後の合波光Ｃ３，Ｃ４が有している光周波数を、レンジビンを補正していない合波光Ｃ１，Ｃ２が有している光周波数と揃えるように補正している。しかし、これは一例に過ぎず、レンジビンを補正していない合波光Ｃ１，Ｃ２が有している光周波数を、レンジビン補正後の合波光Ｃ３，Ｃ４が有している光周波数と揃えるように補正してもよい。
【００３７】
周波数補正処理部２５は、合波光Ｃ１に係る散乱光Ｒ１が存在しているレンジビン（ｎ＝１）についての周波数スペクトルＦＳ（ｎ＝１）をスペクトル積算部２６に出力する。
周波数補正処理部２５は、合波光Ｃ２に係る散乱光Ｒ２が存在しているレンジビン（ｎ＝７）についての周波数スペクトルＦＳ（ｎ＝７）をスペクトル積算部２６に出力する。
周波数補正処理部２５は、合波光Ｃ３に係る散乱光Ｒ３が存在している補正後のレンジビン（ｎ’＝１）についての周波数スペクトルＦＳ（ｎ’＝１）の光周波数を、合波光Ｃ３が有している補正後の光周波数ｆＩＦ１＋ｆｄｐ１に変更する。
周波数補正処理部２５は、光周波数変更後の周波数スペクトルＦＳ（ｎ’＝１）をスペクトル積算部２６に出力する。
周波数補正処理部２５は、合波光Ｃ４に係る散乱光Ｒ４が存在している補正後のレンジビン（ｎ’＝７）についての周波数スペクトルＦＳ（ｎ’＝７）の光周波数を、合波光Ｃ４が有している補正後の光周波数ｆＩＦ１＋ｆｄｐ２に変更する。
周波数補正処理部２５は、光周波数変更後の周波数スペクトルＦＳ（ｎ’＝７）をスペクトル積算部２６に出力する。
【００３８】
スペクトル積算部２６は、周波数補正処理部２５から、周波数スペクトルＦＳ（ｎ＝１）と、周波数スペクトルＦＳ（ｎ＝７）と、光周波数変更後の周波数スペクトルＦＳ（ｎ’＝１）と、光周波数変更後の周波数スペクトルＦＳ（ｎ’＝７）とを取得する。
スペクトル積算部２６は、周波数スペクトルＦＳ（ｎ＝１）と、周波数スペクトルＦＳ（ｎ＝７）と、光周波数変更後の周波数スペクトルＦＳ（ｎ’＝１）と、光周波数変更後の周波数スペクトルＦＳ（ｎ’＝７）とを積算する。スペクトル積算部２６によって積算されることによって、光周波数ｆＩＦ１＋ｆｄｐ１に対応する周波数スペクトルのスペクトル強度と、光周波数ｆＩＦ１＋ｆｄｐ２に対応する周波数スペクトルのスペクトル強度とが大きくなる。
スペクトル積算部２６は、積算後の周波数スペクトルΣＨＦＳをピーク周波数検出部２７に出力する。
【００３９】
ピーク周波数検出部２７は、スペクトル積算部２６から、積算後の周波数スペクトルΣＨＦＳを取得する。
ピーク周波数検出部２７は、積算後の周波数スペクトルΣＨＦＳに含まれている複数のスペクトル強度の中で、閾値以上のスペクトル強度ＦＳｍａｘ１，ＦＳｍａｘ２を特定する。第１の観測対象と第２の観測対象とが空間に存在しているので、２つのスペクトル強度ＦＳｍａｘ１，ＦＳｍａｘ２が特定される。閾値は、ピーク周波数検出部２７の内部メモリに格納されていてもよいし、図１に示すライダ装置の外部から与えられるものであってもよい。
ピーク周波数検出部２７は、それぞれのスペクトル強度ＦＳｍａｘ１，ＦＳｍａｘ２に対応するピーク周波数ｆｐｅａｋ１，ｆｐｅａｋ２を周波数算出処理部２８に出力し、積算後の周波数スペクトルΣＨＦＳをＳＮＲ算出部３０に出力する。
【００４０】
周波数算出処理部２８は、ピーク周波数検出部２７から、それぞれのピーク周波数ｆｐｅａｋ１，ｆｐｅａｋ２を取得する。
周波数算出処理部２８は、トリガ生成部４から、第１の周波数シフト信号と第２の周波数シフト信号とを取得する。
周波数算出処理部２８は、ピーク周波数ｆｐｅａｋ１から第１の周波数シフト信号が示す変調周波数ｆＩＦ１を減算することによって、合波光Ｃ１が有している光周波数ｆＩＦ１＋ｆｄｐ１に含まれているドップラー周波数ｆｄｐ１を算出する。
周波数算出処理部２８は、ピーク周波数ｆｐｅａｋ２から第２の周波数シフト信号が示す変調周波数ｆＩＦ２を減算することによって、合波光Ｃ２が有している光周波数ｆＩＦ２＋ｆｄｐ２に含まれているドップラー周波数ｆｄｐ２を算出する。
周波数算出処理部２８は、それぞれのドップラー周波数ｆｄｐ１，ｆｄｐ２を速度算出処理部２９に出力する。
【００４１】
速度算出部１７は、例えば、図４に示す速度算出回路４２によって実現される。
速度算出部１７は、速度算出処理部２９を備えている。
速度算出処理部２９は、周波数算出処理部２８から、それぞれのドップラー周波数ｆｄｐ１，ｆｄｐ２を取得する。
速度算出処理部２９は、ドップラー周波数ｆｄｐ１から、第１の観測対象の相対速度Ｖ１を算出する。
速度算出処理部２９は、ドップラー周波数ｆｄｐ２から、第２の観測対象の相対速度Ｖ２を算出する。
【００４２】
ＳＮＲ算出部３０は、例えば、図４に示すＳＮＲ算出回路４３によって実現される。
ＳＮＲ算出部３０は、ピーク周波数検出部２７から、積算後の周波数スペクトルΣＨＦＳを取得する。
ＳＮＲ算出部３０は、積算後の周波数スペクトルΣＨＦＳを逆ＦＦＴ処理することによって、それぞれのレンジビン（ｎ）の信号を算出する。
ＳＮＲ算出部３０は、それぞれのレンジビン（ｎ）の信号を帯域外雑音で除算することで、それぞれのレンジビン（ｎ）のＳＮＲを算出する。
ＳＮＲ算出部３０は、それぞれのレンジビン（ｎ）のＳＮＲを距離特性算出部３１に出力する。
【００４３】
距離特性算出部３１は、例えば、図４に示す距離特性算出回路４４によって実現される。
距離特性算出部３１は、レンジ補正部２４から、合波光Ｃ１に係る散乱光Ｒ１が存在しているレンジビン（ｎ＝１）と、合波光Ｃ２に係る散乱光Ｒ２が存在しているレンジビン（ｎ＝７）とを取得する。
距離特性算出部３１は、合波光Ｃ１に係る散乱光Ｒ１が存在しているレンジビン（ｎ＝１）と、Ａ／Ｄ変換部１４のＡ／Ｄ変換レートＲａｔｅと、レンジビン幅Ｒｂｗとから、ライダ装置から第１の観測対象までの距離Ｌ１を算出する。
距離特性算出部３１は、合波光Ｃ２に係る散乱光Ｒ２が存在しているレンジビン（ｎ＝７）と、Ａ／Ｄ変換部１４のＡ／Ｄ変換レートＲａｔｅと、レンジビン幅Ｒｂｗとから、ライダ装置から第２の観測対象までの距離Ｌ２を算出する。
距離特性算出部３１は、レンジビン（ｎ）の距離とＳＮＲ算出部３０により算出されたＳＮＲとの対応関係を示す距離特性（ａ−ｓｃｏｐｅ）を例えば図示せぬ表示装置に表示させる。
【００４４】
図１では、信号処理装置１５の構成要素であるドップラー周波数算出部１６、速度算出部１７、ＳＮＲ算出部３０及び距離特性算出部３１のそれぞれが、図４に示すような専用のハードウェアによって実現されるものを想定している。即ち、信号処理装置１５が、ドップラー周波数算出回路４１、速度算出回路４２、ＳＮＲ算出回路４３及び距離特性算出回路４４によって実現されるものを想定している。
ドップラー周波数算出回路４１、速度算出回路４２、ＳＮＲ算出回路４３及び距離特性算出回路４４のそれぞれは、例えば、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ  Ｓｐｅｃｉｆｉｃ  Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ）、ＦＰＧＡ、又は、これらを組み合わせたものが該当する。
【００４５】
信号処理装置１５の構成要素は、専用のハードウェアによって実現されるものに限るものではなく、信号処理装置１５が、ソフトウェア、ファームウェア、又は、ソフトウェアとファームウェアとの組み合わせによって実現されるものであってもよい。
ソフトウェア又はファームウェアは、プログラムとして、コンピュータのメモリに格納される。コンピュータは、プログラムを実行するハードウェアを意味し、例えば、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）、中央処理装置、処理装置、演算装置、マイクロプロセッサ、マイクロコンピュータ、プロセッサ、あるいは、ＤＳＰ（Ｄｉｇｉｔａｌ  Ｓｉｇｎａｌ  Ｐｒｏｃｅｓｓｏｒ）が該当する。
【００４６】
図５は、信号処理装置１５が、ソフトウェア又はファームウェア等によって実現される場合のコンピュータのハードウェア構成図である。
信号処理装置１５が、ソフトウェア又はファームウェア等によって実現される場合、ドップラー周波数算出部１６、速度算出部１７、ＳＮＲ算出部３０及び距離特性算出部３１におけるそれぞれの処理手順をコンピュータに実行させるためのプログラムがメモリ５１に格納される。そして、コンピュータのプロセッサ５２がメモリ５１に格納されているプログラムを実行する。
【００４７】
また、図４では、信号処理装置１５の構成要素のそれぞれが専用のハードウェアによって実現される例を示し、図５では、信号処理装置１５がソフトウェア又はファームウェア等によって実現される例を示している。しかし、これは一例に過ぎず、信号処理装置１５における一部の構成要素が専用のハードウェアによって実現され、残りの構成要素がソフトウェア又はファームウェア等によって実現されるものであってもよい。
【００４８】
次に、図１に示すライダ装置の動作について説明する。
図６は、図１に示すライダ装置の処理手順を示すフローチャートである。
図７は、信号処理装置１５の処理手順である信号処理方法を示すフローチャートである。
光源１は、光周波数ｆ０の連続光であるレーザ光をパルス変調部２の光分割部３に出力する（図６のステップＳＴ１）。
光分割部３は、光源１からレーザ光を受けると、レーザ光を２分配する（図６のステップＳＴ２）。
光分割部３は、分配後の一方のレーザ光をパルス変調器５に出力し、分配後の他方のレーザ光を参照光として光合波部１２に出力する。
【００４９】
トリガ生成部４のパルス信号生成部４ａは、パルス幅ΔＴのパルス信号を生成する。
パルス信号生成部４ａは、パルス信号をパルス変調器５、Ａ／Ｄ変換部１４及び信号処理装置１５のそれぞれに、周期Ｔｒｅｐで繰り返し出力する。
トリガ生成部４の基準信号生成部４ｂは、周波数ｆの電気信号である基準信号を生成する。
基準信号生成部４ｂは、基準信号を第１のシフト信号生成部４ｃに出力する。
【００５０】
第１のシフト信号生成部４ｃは、基準信号生成部４ｂから基準信号を受けると、基準信号から、変調周波数ｆＩＦ１を示す第１の周波数シフト信号を生成する（図６のステップＳＴ３）。
第１のシフト信号生成部４ｃは、第１の周波数シフト信号をスイッチ４ｅに出力する。
スイッチ４ｅは、第１のシフト信号生成部４ｃから第１の周波数シフト信号を受けると、第１の周波数シフト信号をパルス変調器５及び信号処理装置１５のそれぞれに出力する。
【００５１】
パルス変調器５は、パルス信号生成部４ａからパルス幅ΔＴのパルス信号を受ける毎に、当該パルス信号に従って、光分割部３から出力されたレーザ光をパルス変調する。
即ち、パルス変調器５は、光分割部３から出力されたレーザ光を、パルス幅ΔＴを有するパルス光Ｐに変換する。
パルス変調器５は、スイッチ４ｅから第１の周波数シフト信号が出力されると、パルス光Ｐの光周波数ｆ０を光周波数ｆ０＋ｆＩＦ１にシフトすることによって、光周波数ｆ０＋ｆＩＦ１を有するパルス光Ｐ１を生成する。
パルス変調器５は、パルス光Ｐ１を送信側光学系７に出力する（図６のステップＳＴ４）。
【００５２】
送信側光学系７は、パルス変調器５からパルス光Ｐ１を受けると、パルス光Ｐ１を整形し、整形後のパルス光Ｐ１を送受分離部８に出力する。
送受分離部８は、送信側光学系７から整形後のパルス光Ｐ１を受けると、整形後のパルス光Ｐ１をテレスコープ９に出力する。
テレスコープ９は、送受分離部８から出力された整形後のパルス光Ｐ１を空間に放射する（図６のステップＳＴ５）。
第１の観測対象と第２の観測対象とがパルス光Ｐ１の照射領域内の空間に存在していれば、図８に示すように、パルス光Ｐ１が第１の観測対象及び第２の観測対象のそれぞれによって散乱される。
図８は、ライダ装置から放射されたパルス光Ｐ１と、第１の観測対象からの散乱光Ｒ１と、第２の観測対象からの散乱光Ｒ２とを示す説明図である。
図８の例では、ライダ装置から第１の観測対象までの距離Ｌ１が、ライダ装置から第２の観測対象までの距離Ｌ２よりも短い。このため、第１の観測対象からの散乱光Ｒ１は、第２の観測対象からの散乱光Ｒ２よりも早くライダ装置に戻ってくる。
図９は、ライダ装置から放射されたパルス光Ｐ１，Ｐ２の放射時刻Ｔ１，Ｔ２と、第１の観測対象からの散乱光Ｒ１，Ｒ３の受信時刻Ｔ１’，Ｔ３’と、第２の観測対象からの散乱光Ｒ２，Ｒ４の受信時刻Ｔ２’，Ｔ４’とを示す説明図である。
図９において、横軸は時間を示し、縦軸は光強度を示している。
【００５３】
テレスコープ９は、第１の観測対象によって散乱された後のパルス光Ｐ１を散乱光Ｒ１として集光し、散乱光Ｒ１を送受分離部８に出力する（図６のステップＳＴ６）。
図９の例では、ライダ装置から第１の観測対象までの距離Ｌ１が近いため、テレスコープ９からパルス光Ｐ２が放射される前に、第１の観測対象からの散乱光Ｒ１がテレスコープ９によって集光されている。一方、ライダ装置から第２の観測対象までの距離Ｌ２が遠いため、テレスコープ９からパルス光Ｐ２が放射された後に、第２の観測対象からの散乱光Ｒ２がテレスコープ９によって集光されている。
図９の例では、テレスコープ９からパルス光Ｐ２が放射された後に、第２の観測対象からの散乱光Ｒ２がテレスコープ９によって集光されている。しかし、これは一例に過ぎず、テレスコープ９からパルス光Ｐ２が放射される前に、第２の観測対象からの散乱光Ｒ２がテレスコープ９によって集光されているものであってもよい。
また、図９の例では、テレスコープ９からパルス光Ｐ２が放射される前に、第１の観測対象からの散乱光Ｒ１がテレスコープ９によって集光されている。しかし、これは一例に過ぎず、テレスコープ９からパルス光Ｐ２が放射された後に、第１の観測対象からの散乱光Ｒ１がテレスコープ９によって集光されているものであってもよい。
【００５４】
送受分離部８は、テレスコープ９から第１の観測対象からの散乱光Ｒ１を受けると、散乱光Ｒ１を受信側光学系１０に出力する。
受信側光学系１０は、送受分離部８から散乱光Ｒ１を受けると、散乱光Ｒ１を整形し、整形後の散乱光Ｒ１を光合波部１２に出力する。
【００５５】
光合波部１２は、受信側光学系１０から整形後の散乱光Ｒ１を受けると、整形後の散乱光Ｒ１と光源１から出力された参照光とを混合することによって、合波光Ｃ１をヘテロダイン検波する（図６のステップＳＴ７）。合波光Ｃ１が有している光周波数は、ｆＩＦ１＋ｆｄｐ１である。
光合波部１２は、合波光Ｃ１を受光部１３に出力する。
【００５６】
受光部１３は、光合波部１２から合波光Ｃ１を受けると、合波光Ｃ１を電気信号に変換し、電気信号をＡ／Ｄ変換部１４に出力する。
Ａ／Ｄ変換部１４は、トリガ生成部４のパルス信号生成部４ａからパルス幅ΔＴのパルス信号が出力されている期間中、受光部１３から出力された電気信号を、アナログ信号からデジタル信号Ｄｉｇ（ｔ）に変換する処理を行う。
Ａ／Ｄ変換部１４は、合波光Ｃ１の検波信号Ｄ１を含むデジタル信号Ｄｉｇ（ｔ）を信号処理装置１５のレンジビン分割部２２に出力する。
【００５７】
トリガ生成部４のパルス信号生成部４ａは、パルス信号を出力してから、周期Ｔｒｅｐの時間が経過すると、パルス信号をパルス変調器５、Ａ／Ｄ変換部１４及び信号処理装置１５のそれぞれに出力する。
基準信号生成部４ｂは、基準信号を第１のシフト信号生成部４ｃに出力してから、周期Ｔｒｅｐの時間が経過すると、基準信号を第２のシフト信号生成部４ｄに出力する。
第２のシフト信号生成部４ｄは、基準信号生成部４ｂから基準信号を受けると、基準信号から、変調周波数ｆＩＦ２を示す第２の周波数シフト信号を生成する（図６のステップＳＴ８）。
第２のシフト信号生成部４ｄは、第２の周波数シフト信号をスイッチ４ｅに出力する。
図１に示すライダ装置では、ｆＩＦ２＞ｆＩＦ１であり、変調周波数ｆＩＦ２と変調周波数ｆＩＦ１との周波数差Δｆ（＝ｆＩＦ２−ｆＩＦ１）は、第１の観測対象及び第２の観測対象のそれぞれが、想定される最高速度で移動した場合に生じるドップラー周波数の絶対値の２倍よりも大きい。周波数差Δｆが、ドップラー周波数の絶対値の２倍よりも大きければ、第１の観測対象及び第２の観測対象のそれぞれの移動速度が、想定の範囲内の速度である限り、それぞれの散乱光Ｒ１，Ｒ２，Ｒ３，Ｒ４の光周波数が、互いに異なる周波数になる。
スイッチ４ｅは、第２のシフト信号生成部４ｄから第２の周波数シフト信号を受けると、第２の周波数シフト信号をパルス変調器５及び信号処理装置１５のそれぞれに出力する。
【００５８】
パルス変調器５は、スイッチ４ｅから第２の周波数シフト信号が出力されると、パルス光Ｐの光周波数ｆ０を光周波数ｆ０＋ｆＩＦ２にシフトすることによって、光周波数ｆ０＋ｆＩＦ２を有するパルス光Ｐ２を生成する。
パルス変調器５は、パルス光Ｐ２を送信側光学系７に出力する（図６のステップＳＴ９）。
【００５９】
送信側光学系７は、パルス変調器５からパルス光Ｐ２を受けると、パルス光Ｐ２を整形し、整形後のパルス光Ｐ２を送受分離部８に出力する。
送受分離部８は、送信側光学系７から整形後のパルス光Ｐ２を受けると、整形後のパルス光Ｐ２をテレスコープ９に出力する。
テレスコープ９は、送受分離部８から出力された整形後のパルス光Ｐ２を空間に放射する（図６のステップＳＴ１０）。
【００６０】
テレスコープ９は、第２の観測対象によって散乱された後のパルス光Ｐ１を散乱光Ｒ２として集光し、散乱光Ｒ２を送受分離部８に出力する（図６のステップＳＴ１１）。
また、テレスコープ９は、第１の観測対象によって散乱された後のパルス光Ｐ２を散乱光Ｒ３として集光し、散乱光Ｒ３を送受分離部８に出力する（図６のステップＳＴ１１）。
図９の例では、第２の観測対象によって散乱された後のパルス光Ｐ１である散乱光Ｒ２と、第１の観測対象によって散乱された後のパルス光Ｐ２である散乱光Ｒ３とが同時にテレスコープ９によって集光されている。
ここでの同時は、散乱光Ｒ２の一部と散乱光Ｒ３の一部とが互いに重なり合っていればよく、散乱光Ｒ２の受信時刻と散乱光Ｒ３の受信時刻とが厳密に同じ時刻であるものに限るものではない。したがって、ここでの同時には、散乱光Ｒ２の受信時刻と散乱光Ｒ３の受信時刻とが多少ずれているものも含まれる。
【００６１】
送受分離部８は、テレスコープ９から、第２の観測対象からの散乱光Ｒ２と第１の観測対象からの散乱光Ｒ３とが互いに重なり合っている散乱光Ｒ２，３を受けると、散乱光Ｒ２，３を受信側光学系１０に出力する。
受信側光学系１０は、送受分離部８から散乱光Ｒ２，３を受けると、散乱光Ｒ２，３を整形し、整形後の散乱光Ｒ２，３を光合波部１２に出力する。
【００６２】
光合波部１２は、受信側光学系１０から整形後の散乱光Ｒ２，３を受けると、整形後の散乱光Ｒ２，３と光源１から出力された参照光とを混合することによって、合波光Ｃ２，３をヘテロダイン検波する（図６のステップＳＴ１２）。合波光Ｃ２，３が有している光周波数は、ｆＩＦ１＋ｆｄｐ２，ｆＩＦ２＋ｆｄｐ１である。
光合波部１２は、合波光Ｃ２，３を受光部１３に出力する。
【００６３】
受光部１３は、光合波部１２から合波光Ｃ２，３を受けると、合波光Ｃ２，３を電気信号に変換し、電気信号をＡ／Ｄ変換部１４に出力する。
Ａ／Ｄ変換部１４は、トリガ生成部４のパルス信号生成部４ａからパルス幅ΔＴのパルス信号が出力されている期間中、受光部１３から出力された電気信号を、アナログ信号からデジタル信号Ｄｉｇ（ｔ）に変換する処理を行う。
Ａ／Ｄ変換部１４は、合波光Ｃ２，３の検波信号Ｄ２，３を含むデジタル信号Ｄｉｇ（ｔ）を信号処理装置１５のレンジビン分割部２２に出力する。
【００６４】
その後、テレスコープ９は、第２の観測対象によって散乱された後のパルス光Ｐ２を散乱光Ｒ４として集光し、散乱光Ｒ４を送受分離部８に出力する（図６のステップＳＴ１３）。
送受分離部８は、テレスコープ９から、第２の観測対象からの散乱光Ｒ４を受けると、散乱光Ｒ４を受信側光学系１０に出力する。
受信側光学系１０は、送受分離部８から散乱光Ｒ４を受けると、散乱光Ｒ４を整形し、整形後の散乱光Ｒ４を光合波部１２に出力する。
【００６５】
光合波部１２は、受信側光学系１０から整形後の散乱光Ｒ４を受けると、整形後の散乱光Ｒ４と光源１から出力された参照光とを混合することによって、合波光Ｃ４をヘテロダイン検波する（図６のステップＳＴ１４）。合波光Ｃ４が有している光周波数は、ｆＩＦ２＋ｆｄｐ２である。
光合波部１２は、合波光Ｃ４を受光部１３に出力する。
【００６６】
受光部１３は、光合波部１２から合波光Ｃ４を受けると、合波光Ｃ４を電気信号に変換し、電気信号をＡ／Ｄ変換部１４に出力する。
Ａ／Ｄ変換部１４は、トリガ生成部４のパルス信号生成部４ａからパルス幅ΔＴのパルス信号が出力されている期間中、受光部１３から出力された電気信号を、アナログ信号からデジタル信号Ｄｉｇ（ｔ）に変換する処理を行う。
Ａ／Ｄ変換部１４は、合波光Ｃ４の検波信号Ｄ４を含むデジタル信号Ｄｉｇ（ｔ）を信号処理装置１５のレンジビン分割部２２に出力する。
【００６７】
レンジビン分割部２２は、Ａ／Ｄ変換部１４からデジタル信号Ｄｉｇ（ｔ）を受けると、デジタル信号Ｄｉｇ（ｔ）を時間方向に分割する（図７のステップＳＴ２１）。
デジタル信号Ｄｉｇ（ｔ）における時間方向の分割幅Δｔは、レンジビン幅Ｒｂｗに対応している。
テレスコープ９が、図１１に示すように、レンジビン（０）に相当するサンプリング時刻ｔ＝０のときにパルス光Ｐ１を放射し、レンジビン（６）に相当するサンプリング時刻ｔ＝６のときにパルス光Ｐ２を放射した場合を想定する。
そして、ライダ装置から第１の観測対象までの距離Ｌ１が、レンジビン（１）の距離に対応し、ライダ装置から第２の観測対象までの距離Ｌ２が、レンジビン（７）の距離に対応しているものとする。
この場合、分割後のデジタル信号Ｄｉｇ（１）は、合波光Ｃ１の検波信号Ｄ１を含んでおり、分割後のデジタル信号Ｄｉｇ（７）は、合波光Ｃ２，３の検波信号Ｄ２，３を含んでおり、分割後のデジタル信号Ｄｉｇ（１３）は、合波光Ｃ４の検波信号Ｄ４を含んでいる。
レンジビン分割部２２は、分割後のそれぞれのデジタル信号Ｄｉｇ（ｔ）をレンジビン信号（ｎ）として周波数解析部２３に出力する。
【００６８】
周波数解析部２３は、レンジビン分割部２２から出力されたそれぞれのレンジビン信号（ｎ）に対するＦＦＴ処理を実施することで、それぞれのレンジビン信号（ｎ）の周波数スペクトルＦＳ（ｎ）を算出する（図７のステップＳＴ２２）。
周波数解析部２３は、それぞれの周波数スペクトルＦＳ（ｎ）をレンジ補正部２４に出力する。
レンジビン信号（７）は、合波光Ｃ２，３の検波信号Ｄ２，３を含んでいる。合波光Ｃ２が有している光周波数ｆＩＦ１＋ｆｄｐ２と、合波光Ｃ３が有している光周波数ｆＩＦ２＋ｆｄｐ１とが異なるため、レンジビン信号（７）は、図１０Ａに示すように、２つのピークスペクトルを有している。
レンジビン信号（１）は、合波光Ｃ１の検波信号Ｄ１を含んでいる。このため、レンジビン信号（１）は、１つのピークスペクトルを有している。
レンジビン信号（１３）は、合波光Ｃ４の検波信号Ｄ４を含んでいる。このため、レンジビン信号（１３）は、１つのピークスペクトルを有している。
レンジビン信号（１），（７），（１３）以外のレンジビン信号（ｎ）は、ピークスペクトルを有しておらず、当該レンジビン信号（ｎ）が有している周波数スペクトルは、ほぼ０である。
図１０Ａは、レンジビン信号（７）が有しているピークスペクトルを示す説明図である。
図１０Ｂは、パルス光Ｐ１が有している光周波数とパルス光Ｐ２が有している光周波数とが同一周波数である場合の、レンジビン信号（７）が有しているピークスペクトルを示す説明図である。図１０Ｂの場合、合波光Ｃ２と合波光Ｃ３とを識別することができない。
図１０Ａ及び図１０Ｂにおいて、横軸は時間を示し、縦軸はスペクトル強度を示している。
【００６９】
レンジ補正部２４は、周波数解析部２３から、それぞれの周波数スペクトルＦＳ（ｎ）を取得する。
レンジ補正部２４は、トリガ生成部４から、第１の周波数シフト信号と第２の周波数シフト信号とを取得する。
レンジ補正部２４は、それぞれの周波数スペクトルＦＳ（ｎ）から、複数のピークスペクトルを検出する。
即ち、レンジ補正部２４は、周波数スペクトルＦＳ（１）から、ピークスペクトルＳｐ１を検出する。
レンジ補正部２４は、周波数スペクトルＦＳ（７）から、ピークスペクトルＳｐ２と、ピークスペクトルＳｐ３とを検出する。
また、レンジ補正部２４は、周波数スペクトルＦＳ（１３）から、ピークスペクトルＳｐ４を検出する。
また、レンジ補正部２４は、ピークスペクトルＳｐ１に対応するピーク周波数ｆｐ１を検出し、ピークスペクトルＳｐ２に対応するピーク周波数ｆｐ２を検出する。
レンジ補正部２４は、ピークスペクトルＳｐ３に対応するピーク周波数ｆｐ３を検出し、ピークスペクトルＳｐ４に対応するピーク周波数ｆｐ４を検出する。
【００７０】
レンジ補正部２４は、ピーク周波数ｆｐ１と第１の周波数シフト信号が示す変調周波数ｆＩＦ１との差分の絶対値｜Δｆ１−１｜を算出し、ピーク周波数ｆｐ１と第２の周波数シフト信号が示す変調周波数ｆＩＦ２との差分の絶対値｜Δｆ１−２｜を算出する。
レンジ補正部２４は、差分の絶対値｜Δｆ１−１｜が差分の絶対値｜Δｆ１−２｜以下であれば、ピークスペクトルＳｐ１を有する合波光が、パルス光Ｐ１に対応する合波光Ｃ１、又は、パルス光Ｐ１に対応する合波光Ｃ２であると判定する。
レンジ補正部２４は、差分の絶対値｜Δｆ１−１｜が差分の絶対値｜Δｆ１−２｜よりも大きければ、ピークスペクトルＳｐ１を有する合波光が、パルス光Ｐ２に対応する合波光Ｃ３、又は、パルス光Ｐ２に対応する合波光Ｃ４であると判定する。
レンジ補正部２４は、ピーク周波数ｆｐ２と変調周波数ｆＩＦ１との差分の絶対値｜Δｆ２−１｜を算出し、ピーク周波数ｆｐ２と変調周波数ｆＩＦ２との差分の絶対値｜Δｆ２−２｜を算出する。
レンジ補正部２４は、差分の絶対値｜Δｆ２−１｜が差分の絶対値｜Δｆ２−２｜以下であれば、ピークスペクトルＳｐ２を有する合波光が、パルス光Ｐ１に対応する合波光Ｃ１、又は、パルス光Ｐ１に対応する合波光Ｃ２であると判定する。
レンジ補正部２４は、差分の絶対値｜Δｆ２−１｜が差分の絶対値｜Δｆ２−２｜よりも大きければ、ピークスペクトルＳｐ２を有する合波光が、パルス光Ｐ２に対応する合波光Ｃ３、又は、パルス光Ｐ２に対応する合波光Ｃ４であると判定する。
【００７１】
レンジ補正部２４は、ピーク周波数ｆｐ３と変調周波数ｆＩＦ１との差分の絶対値｜Δｆ３−１｜を算出し、ピーク周波数ｆｐ３と変調周波数ｆＩＦ２との差分の絶対値｜Δｆ３−２｜を算出する。
レンジ補正部２４は、差分の絶対値｜Δｆ３−１｜が差分の絶対値｜Δｆ３−２｜以下であれば、ピークスペクトルＳｐ３を有する合波光が、パルス光Ｐ１に対応する合波光Ｃ１、又は、パルス光Ｐ１に対応する合波光Ｃ２であると判定する。
レンジ補正部２４は、差分の絶対値｜Δｆ３−１｜が差分の絶対値｜Δｆ３−２｜よりも大きければ、ピークスペクトルＳｐ３を有する合波光が、パルス光Ｐ２に対応する合波光Ｃ３、又は、パルス光Ｐ２に対応する合波光Ｃ４であると判定する。
レンジ補正部２４は、ピーク周波数ｆｐ４と変調周波数ｆＩＦ１との差分の絶対値｜Δｆ４−１｜を算出し、ピーク周波数ｆｐ４と変調周波数ｆＩＦ２との差分の絶対値｜Δｆ４−２｜を算出する。
レンジ補正部２４は、差分の絶対値｜Δｆ４−１｜が差分の絶対値｜Δｆ４−２｜以下であれば、ピークスペクトルＳｐ４を有する合波光が、パルス光Ｐ１に対応する合波光Ｃ１、又は、パルス光Ｐ１に対応する合波光Ｃ２であると判定する。
レンジ補正部２４は、差分の絶対値｜Δｆ４−１｜が差分の絶対値｜Δｆ４−２｜よりも大きければ、ピークスペクトルＳｐ４を有する合波光が、パルス光Ｐ２に対応する合波光Ｃ３、又は、パルス光Ｐ２に対応する合波光Ｃ４であると判定する。
【００７２】
ここでは、説明の便宜上、ピークスペクトルＳｐ１を有する合波光及びピークスペクトルＳｐ２を有する合波光のそれぞれが、パルス光Ｐ１に対応する合波光Ｃ１、又は、パルス光Ｐ１に対応する合波光Ｃ２であるものとする。
また、ピークスペクトルＳｐ３を有する合波光及びピークスペクトルＳｐ４を有する合波光のそれぞれが、パルス光Ｐ２に対応する合波光Ｃ３、又は、パルス光Ｐ２に対応する合波光Ｃ４であるものとする。
この場合、レンジ補正部２４は、ピークスペクトルＳｐ１のレンジビンが、ピークスペクトルＳｐ２のレンジビン以下であれば、ピークスペクトルＳｐ１を有する合波光が、パルス光Ｐ１に対応する合波光Ｃ１であり、ピークスペクトルＳｐ２を有する合波光が、パルス光Ｐ１に対応する合波光Ｃ２であると判定する。
レンジ補正部２４は、ピークスペクトルＳｐ１のレンジビンが、ピークスペクトルＳｐ２のレンジビンよりも大きければ、ピークスペクトルＳｐ１を有する合波光が、パルス光Ｐ１に対応する合波光Ｃ２であり、ピークスペクトルＳｐ２を有する合波光が、パルス光Ｐ１に対応する合波光Ｃ１であると判定する。
【００７３】
また、レンジ補正部２４は、ピークスペクトルＳｐ３のレンジビンが、ピークスペクトルＳｐ４のレンジビン以下であれば、ピークスペクトルＳｐ３を有する合波光が、パルス光Ｐ２に対応する合波光Ｃ３であり、ピークスペクトルＳｐ４を有する合波光が、パルス光Ｐ２に対応する合波光Ｃ４であると判定する。
レンジ補正部２４は、ピークスペクトルＳｐ３のレンジビンが、ピークスペクトルＳｐ４のレンジビンよりも大きければ、ピークスペクトルＳｐ３を有する合波光が、パルス光Ｐ２に対応する合波光Ｃ４であり、ピークスペクトルＳｐ４を有する合波光が、パルス光Ｐ２に対応する合波光Ｃ３であると判定する。
ここでは、説明の便宜上、ピークスペクトルＳｐ１を有する合波光が合波光Ｃ１であり、ピークスペクトルＳｐ２を有する合波光が合波光Ｃ２であるとする。また、ピークスペクトルＳｐ３を有する合波光が合波光Ｃ３であり、ピークスペクトルＳｐ４を有する合波光が合波光Ｃ４であるとする。
【００７４】
レンジ補正部２４は、合波光Ｃ１，Ｃ２がパルス光Ｐ１に対応する合波光であり、合波光Ｃ３，Ｃ４がパルス光Ｐ２に対応する合波光であると判定すると、図１１に示すように、合波光Ｃ３が存在しているレンジビン（７）と、合波光Ｃ４が存在しているレンジビン（１３）とを補正する。
即ち、レンジ補正部２４は、合波光Ｃ３が存在しているレンジビン（７）から周期Ｔｒｅｐ（＝６）を減算することによって、合波光Ｃ３が存在しているレンジビン（７）をレンジビン（１）に補正する（図７のステップＳＴ２３）。
また、レンジ補正部２４は、合波光Ｃ４が存在しているレンジビン（１３）から周期Ｔｒｅｐ（＝６）を減算することによって、合波光Ｃ４が存在しているレンジビン（１３）をレンジビン（７）に補正する（図７のステップＳＴ２３）。
【００７５】
図１１は、合波光Ｃ１に係る散乱光Ｒ１が存在しているレンジビンと、合波光Ｃ２，Ｃ３に係る散乱光Ｒ２，Ｒ３が存在しているレンジビンと、合波光Ｃ４に係る散乱光Ｒ４が存在しているレンジビンとを示す説明図である。
図１１において、横軸は時間を示し、縦軸は光強度を示している。図１１では、レンジビンを「レンジ」のように簡略している。
図１１の例では、散乱光Ｒ１が存在しているレンジビンがレンジビン（１）、散乱光Ｒ２，Ｒ３が存在しているレンジビンがレンジビン（７）、散乱光Ｒ４が存在しているレンジビンがレンジビン（１３）である。
パルス光Ｐ１の放射時刻Ｔ１とパルス光Ｐ２の放射時刻Ｔ２との時刻差がレンジビン（６）に相当しており、周期Ｔｒｅｐが６である。このため、合波光Ｃ３に係る散乱光Ｒ３が存在している補正後のレンジビンであるｎ’は、１（＝７−６）となり、合波光Ｃ４に係る散乱光Ｒ４が存在している補正後のレンジビンであるｎ’は、７（＝１３−６）となる。
レンジ補正部２４は、複数の周波数スペクトルＦＳ（１）〜ＦＳ（Ｎ）のうち、合波光Ｃ１，Ｃ２に係る散乱光Ｒ１，Ｒ２が存在しているレンジビン（ｎ）についての周波数スペクトルＦＳ（ｎ）を周波数補正処理部２５に出力する。
また、レンジ補正部２４は、複数の周波数スペクトルＦＳ（１）〜ＦＳ（Ｎ）のうち、合波光Ｃ３，Ｃ４に係る散乱光Ｒ３，Ｒ４が存在している補正後のレンジビン（ｎ’）についての周波数スペクトルＦＳ（ｎ’）を周波数補正処理部２５に出力する。
レンジ補正部２４は、合波光Ｃ１が有している光周波数ｆＩＦ１＋ｆｄｐ１（＝ｆｐ１）と、合波光Ｃ２が有している光周波数ｆＩＦ１＋ｆｄｐ２（＝ｆｐ２）と、合波光Ｃ３が有している光周波数ｆＩＦ２＋ｆｄｐ１（＝ｆｐ３）と、合波光Ｃ４が有している光周波数ｆＩＦ２＋ｆｄｐ２（＝ｆｐ４）とを周波数補正処理部２５に出力する。
レンジ補正部２４は、合波光Ｃ１に係る散乱光Ｒ１が存在しているレンジビンと、合波光Ｃ２に係る散乱光Ｒ２が存在しているレンジビンとを距離特性算出部３１に出力する。
【００７６】
周波数補正処理部２５は、トリガ生成部４から、第１の周波数シフト信号と第２の周波数シフト信号とを取得する。
周波数補正処理部２５は、第１の周波数シフト信号が示す変調周波数ｆＩＦ１と第２の周波数シフト信号が示す変調周波数ｆＩＦ２との周波数差Δｆ（＝ｆＩＦ２−ｆＩＦ１）を算出する。
周波数補正処理部２５は、レンジ補正部２４から、合波光Ｃ１，Ｃ２に係る散乱光Ｒ１，Ｒ２が存在しているレンジビン（ｎ）についての周波数スペクトルＦＳ（ｎ）と、合波光Ｃ３，Ｃ４に係る散乱光Ｒ３，Ｒ４が存在している補正後のレンジビン（ｎ’）についての周波数スペクトルＦＳ（ｎ’）とを取得する。
周波数補正処理部２５は、レンジ補正部２４から、合波光Ｃ１が有している光周波数ｆＩＦ１＋ｆｄｐ１と、合波光Ｃ２が有している光周波数ｆＩＦ１＋ｆｄｐ２と、合波光Ｃ３が有している光周波数ｆＩＦ２＋ｆｄｐ１と、合波光Ｃ４が有している光周波数ｆＩＦ２＋ｆｄｐ２とを取得する。
【００７７】
周波数補正処理部２５は、合波光Ｃ３が有している光周波数ｆＩＦ２＋ｆｄｐ１から周波数差Δｆを減算することによって、合波光Ｃ３が有している光周波数を補正する（図７のステップＳＴ２４）。合波光Ｃ３が有している補正後の光周波数は、ｆＩＦ１＋ｆｄｐ１であり、合波光Ｃ１が有している光周波数ｆＩＦ１＋ｆｄｐ１と同じ周波数である。
周波数補正処理部２５は、合波光Ｃ４が有している光周波数ｆＩＦ２＋ｆｄｐ２から周波数差Δｆを減算することによって、合波光Ｃ４が有している光周波数を補正する（図７のステップＳＴ２４）。合波光Ｃ４が有している補正後の光周波数は、ｆＩＦ１＋ｆｄｐ２であり、合波光Ｃ２が有している光周波数ｆＩＦ１＋ｆｄｐ２と同じ周波数である。
【００７８】
周波数補正処理部２５は、合波光Ｃ１に係る散乱光Ｒ１が存在しているレンジビン（ｎ＝１）についての周波数スペクトルＦＳ（ｎ＝１）をスペクトル積算部２６に出力する。
周波数補正処理部２５は、合波光Ｃ２に係る散乱光Ｒ２が存在しているレンジビン（ｎ＝７）についての周波数スペクトルＦＳ（ｎ＝７）をスペクトル積算部２６に出力する。
周波数補正処理部２５は、合波光Ｃ３に係る散乱光Ｒ３が存在している補正後のレンジビン（ｎ’＝１）についての周波数スペクトルＦＳ（ｎ’＝１）の光周波数を、合波光Ｃ３が有している補正後の光周波数ｆＩＦ１＋ｆｄｐ１に変更する。
周波数補正処理部２５は、光周波数変更後の周波数スペクトルＦＳ（ｎ’＝１）をスペクトル積算部２６に出力する。
周波数補正処理部２５は、合波光Ｃ４に係る散乱光Ｒ４が存在している補正後のレンジビン（ｎ’＝７）についての周波数スペクトルＦＳ（ｎ’＝７）の光周波数を、合波光Ｃ４が有している補正後の光周波数ｆＩＦ１＋ｆｄｐ２に変更する。
周波数補正処理部２５は、光周波数変更後の周波数スペクトルＦＳ（ｎ’＝７）をスペクトル積算部２６に出力する。
【００７９】
スペクトル積算部２６は、周波数補正処理部２５から、周波数スペクトルＦＳ（ｎ＝１）と、周波数スペクトルＦＳ（ｎ＝７）と、光周波数変更後の周波数スペクトルＦＳ（ｎ’＝１）と、光周波数変更後の周波数スペクトルＦＳ（ｎ’＝７）とを取得する。
スペクトル積算部２６は、周波数スペクトルＦＳ（ｎ＝１）と、周波数スペクトルＦＳ（ｎ＝７）と、光周波数変更後の周波数スペクトルＦＳ（ｎ’＝１）と、光周波数変更後の周波数スペクトルＦＳ（ｎ’＝７）とを積算する（図７のステップＳＴ２５）。スペクトル積算部２６によって積算されることによって、光周波数ｆＩＦ１＋ｆｄｐ１に対応する周波数スペクトルのスペクトル強度と、光周波数ｆＩＦ１＋ｆｄｐ２に対応する周波数スペクトルのスペクトル強度とが大きくなる。
スペクトル積算部２６は、積算後の周波数スペクトルΣＨＦＳをピーク周波数検出部２７に出力する。
【００８０】
ピーク周波数検出部２７は、スペクトル積算部２６から、積算後の周波数スペクトルΣＨＦＳを取得する。
ピーク周波数検出部２７は、積算後の周波数スペクトルΣＨＦＳに含まれている複数のスペクトル強度の中で、閾値以上のスペクトル強度ＦＳｍａｘ１，ＦＳｍａｘ２を特定する。第１の観測対象と第２の観測対象とが空間に存在しているので、２つのスペクトル強度ＦＳｍａｘ１，ＦＳｍａｘ２が特定される。
ピーク周波数検出部２７は、スペクトル強度ＦＳｍａｘ１に対応するピーク周波数ｆｐｅａｋ１と、スペクトル強度ＦＳｍａｘ２に対応するピーク周波数ｆｐｅａｋ２とを周波数算出処理部２８に出力する（図７のステップＳＴ２６）。
ピーク周波数ｆｐｅａｋ１は、合波光Ｃ１が有している光周波数ｆＩＦ１＋ｆｄｐ１及び光周波数補正後の合波光Ｃ３が有している光周波数ｆＩＦ１＋ｆｄｐ１のそれぞれに対応している。
ピーク周波数ｆｐｅａｋ２は、合波光Ｃ２が有している光周波数ｆＩＦ１＋ｆｄｐ２及び光周波数補正後の合波光Ｃ４が有している光周波数ｆＩＦ１＋ｆｄｐ２のそれぞれに対応している。
また、ピーク周波数検出部２７は、積算後の周波数スペクトルΣＨＦＳをＳＮＲ算出部３０に出力する。
【００８１】
周波数算出処理部２８は、ピーク周波数検出部２７から、それぞれのピーク周波数ｆｐｅａｋ１，ｆｐｅａｋ２を取得する。
周波数算出処理部２８は、トリガ生成部４から、第１の周波数シフト信号と第２の周波数シフト信号とを取得する。
周波数算出処理部２８は、ピーク周波数ｆｐｅａｋ１から第１の周波数シフト信号が示す変調周波数ｆＩＦ１を減算することによって、合波光Ｃ１及び光周波数補正後の合波光Ｃ３のそれぞれが有している光周波数ｆＩＦ１＋ｆｄｐ１に含まれているドップラー周波数ｆｄｐ１を算出する（図７のステップＳＴ２７）。
周波数算出処理部２８は、ピーク周波数ｆｐｅａｋ２から第２の周波数シフト信号が示す変調周波数ｆＩＦ２を減算することによって、合波光Ｃ２及び光周波数補正後の合波光Ｃ４のそれぞれが有している光周波数ｆＩＦ１＋ｆｄｐ２に含まれているドップラー周波数ｆｄｐ２を算出する（図７のステップＳＴ２７）。
周波数算出処理部２８は、それぞれのドップラー周波数ｆｄｐ１，ｆｄｐ２を速度算出処理部２９に出力する。
【００８２】
速度算出処理部２９は、周波数算出処理部２８から、それぞれのドップラー周波数ｆｄｐ１，ｆｄｐ２を取得する。
速度算出処理部２９は、以下の式（１）に示すように、ドップラー周波数ｆｄｐ１から、第１の観測対象の相対速度Ｖ１を算出する（図７のステップＳＴ２８）。
Ｖ１＝λ×ｆｄｐ１／２    （１）
式（１）において、λは、それぞれのパルス光Ｐ１，Ｐ２の波長である。
【００８３】
速度算出処理部２９は、以下の式（２）に示すように、ドップラー周波数ｆｄｐ２から、第２の観測対象の相対速度Ｖ２を算出する（図７のステップＳＴ２８）。
Ｖ２＝λ×ｆｄｐ２／２    （２）
速度算出処理部２９は、第１の観測対象の相対速度Ｖ１及び第２の観測対象の相対速度Ｖ２のそれぞれを例えば図示せぬ表示装置に表示させる。
【００８４】
ＳＮＲ算出部３０は、ピーク周波数検出部２７から、積算後の周波数スペクトルΣＨＦＳを取得する。
ＳＮＲ算出部３０は、積算後の周波数スペクトルΣＨＦＳを逆ＦＦＴ処理することによって、それぞれのレンジビン（ｎ）の信号を算出する。
ＳＮＲ算出部３０は、それぞれのレンジビン（ｎ）の信号を帯域外雑音で除算することで、それぞれのレンジビン（ｎ）のＳＮＲを算出する。
ＳＮＲ算出部３０は、それぞれのレンジビン（ｎ）のＳＮＲを距離特性算出部３１に出力する。
【００８５】
距離特性算出部３１は、レンジ補正部２４から、合波光Ｃ１に係る散乱光Ｒ１が存在しているレンジビン（ｎ＝１）と、合波光Ｃ２に係る散乱光Ｒ２が存在しているレンジビン（ｎ＝７）とを取得する。
距離特性算出部３１は、以下の式（３）に示すように、合波光Ｃ１に係る散乱光Ｒ１が存在しているレンジビン（ｎ＝１）と、Ａ／Ｄ変換部１４のＡ／Ｄ変換レートＲａｔｅと、レンジビン幅Ｒｂｗとから、ライダ装置から第１の観測対象までの距離Ｌ１を算出する。
Ｌ１＝ｃ×Ｒａｔｅ×ｎ／２
＝ｃ×Ｒａｔｅ×１／２（３）
式（３）において、ｃは、光速である。
【００８６】
距離特性算出部３１は、以下の式（４）に示すように、合波光Ｃ２に係る散乱光Ｒ２が存在しているレンジビン（ｎ＝７）と、Ａ／Ｄ変換部１４のＡ／Ｄ変換レートＲａｔｅと、レンジビン幅Ｒｂｗとから、ライダ装置から第２の観測対象までの距離Ｌ２を算出する。
Ｌ２＝ｃ×Ｒａｔｅ×ｎ／２
＝ｃ×Ｒａｔｅ×７／２（４）
距離特性算出部３１は、図１２に示すように、レンジビン（ｎ）の距離とＳＮＲ算出部３０により算出されたＳＮＲとの対応関係を示す距離特性（ａ−ｓｃｏｐｅ）を例えば図示せぬ表示装置に表示させる。
図１２は、レンジビン（ｎ）の距離とＳＮＲとの対応関係を示す距離特性（ａ−ｓｃｏｐｅ）の説明図である。
図１２において、横軸はレンジビン（ｎ）の距離［ｍ］を示し、縦軸はＳＮＲ［ｄＢ］を示している。
【００８７】
以上の実施の形態１では、信号処理装置１５が、空間に存在している複数の観測対象におけるそれぞれの移動速度として、ライダ装置に対するそれぞれの観測対象の相対速度を算出するように構成した。光源１から出力されたレーザ光から、光周波数が互いに異なる複数のパルス光が生成されて、それぞれのパルス光が空間に放射され、それぞれの観測対象によって散乱された後のそれぞれのパルス光が散乱光として受信されて、それぞれの散乱光とレーザ光との合波光が検波されている。信号処理装置１５は、生成された複数のパルス光が有している光周波数と、それぞれの合波光の検波信号とから、それぞれの散乱光が有している光周波数に含まれている、それぞれの観測対象の移動に伴うドップラー周波数を算出するドップラー周波数算出部１６と、ドップラー周波数算出部１６により算出されたそれぞれのドップラー周波数から、それぞれの観測対象の相対速度を算出する速度算出部１７とを備えている。したがって、信号処理装置１５は、遠方散乱光と近傍散乱光とが、互いに重なり合った状態であっても、遠方の観測対象の移動速度と近傍の観測対象の移動速度とが同じ速度であるか否かにかかわらず、それぞれの観測対象の移動速度を算出することができる。
【００８８】
図１に示すライダ装置は、光を扱う複数の構成要素を備えており、それぞれの構成要素は、他の構成要素と光ファイバによって接続されている。そして、それぞれの構成要素が、他の構成要素と光ファイバを介して光を送受信している。しかし、これは一例に過ぎず、それぞれの構成要素が、他の構成要素と空間伝搬によって光を送受信するようにしてもよい。複数の構成要素は、光源１、光分割部３、パルス変調器５、送信側光学系７、送受分離部８、テレスコープ９、受信側光学系１０、光合波部１２及び受光部１３である。
【００８９】
図３に示す信号処理装置１５では、ドップラー周波数算出部１６が、光周波数補正部２１及び周波数算出処理部２８を備えている。また、光周波数補正部２１は、それぞれの検波信号Ｄ１，Ｄ２，Ｄ３，Ｄ４のＳＮＲを高めるために、レンジビン分割部２２、周波数解析部２３、レンジ補正部２４、周波数補正処理部２５及びスペクトル積算部２６を備えている。
しかし、それぞれの検波信号Ｄ１，Ｄ２のＳＮＲを高めなくても、それぞれの観測対象の移動に伴うドップラー周波数Ｖ１，Ｖ２を算出することが可能であれば、光周波数補正部２１が、レンジビン分割部２２、周波数解析部２３、レンジ補正部２４、周波数補正処理部２５及びスペクトル積算部２６を備えていなくてもよい。
光周波数補正部２１が、レンジビン分割部２２、周波数解析部２３、レンジ補正部２４、周波数補正処理部２５及びスペクトル積算部２６を備えていない場合、ピーク周波数検出部２７は、それぞれの検波信号Ｄ１，Ｄ２をＦＦＴ処理して、それぞれの検波信号Ｄ１，Ｄ２の周波数スペクトルを求める。そして、ピーク周波数検出部２７は、それぞれの周波数スペクトルに含まれている複数のスペクトル強度の中で、閾値以上のスペクトル強度ＦＳｍａｘ１，ＦＳｍａｘ２を特定する。ピーク周波数検出部２７は、スペクトル強度ＦＳｍａｘ１に対応するピーク周波数ｆｐｅａｋ１と、スペクトル強度ＦＳｍａｘ２に対応するピーク周波数ｆｐｅａｋ２とを周波数算出処理部２８に出力する。
【００９０】
実施の形態２．
実施の形態２では、パルス変調器５が、光分岐部５ａ、第１の変調部５ｂ、第２の変調部５ｃ及び光合波部５ｄを備えているライダ装置について説明する。
【００９１】
実施の形態２に係るライダ装置の構成は、実施の形態１に係るライダ装置の構成と同様であり、実施の形態２に係るライダ装置を示す構成図は、図１である。
図１３は、実施の形態２に係るライダ装置のパルス変調器５を示す構成図である。
図１３に示すパルス変調器５は、光分岐部５ａ、第１の変調部５ｂ、第２の変調部５ｃ及び光合波部５ｄを備えている。
光分岐部５ａは、カプラー、又は、光スイッチ等によって実現される。
光分岐部５ａは、光分割部３から出力されたレーザ光を２分岐し、２分岐後の一方のレーザ光を第１の変調部５ｂに出力し、２分岐後の他方のレーザ光を第２の変調部５ｃに出力する。
図１３に示す光分岐部５ａは、光分割部３から出力されたレーザ光を２分岐している。しかし、これは一例に過ぎず、光分岐部５ａは、光分割部３から出力されたレーザ光の出力先を、第１の変調部５ｂ又は第２の変調部５ｃに、交互に切り替えるスイッチであってもよい。
【００９２】
第１の変調部５ｂは、例えば、音響光学素子を用いている変調素子又はニオブ酸リチウム結晶を用いている変調素子と、半導体光アンプ等の光増幅器とによって実現される。
第１の変調部５ｂは、パルス信号生成部４ａからパルス幅ΔＴのパルス信号を受ける毎に、当該パルス信号に従って、光分岐部５ａから出力されたレーザ光をパルス変調することで、当該レーザ光を、パルス幅ΔＴを有するパルス光に変換する。
また、第１の変調部５ｂは、トリガ生成部４のスイッチ４ｅから出力された第１の周波数シフト信号に従って、当該パルス光の光周波数ｆ０を光周波数ｆ０＋ｆＩＦ１にシフトすることによって、光周波数ｆ０＋ｆＩＦ１を有しているパルス光Ｐ１を生成する。
第１の変調部５ｂは、パルス光Ｐ１を光合波部５ｄに出力する。
【００９３】
第２の変調部５ｃは、例えば、音響光学素子を用いている変調素子又はニオブ酸リチウム結晶を用いている変調素子と、半導体光アンプ等の光増幅器とによって実現される。
第２の変調部５ｃは、パルス信号生成部４ａからパルス幅ΔＴのパルス信号を受ける毎に、当該パルス信号に従って、光分岐部５ａから出力されたレーザ光をパルス変調することで、当該レーザ光を、パルス幅ΔＴを有するパルス光に変換する。
また、第２の変調部５ｃは、トリガ生成部４のスイッチ４ｅから出力された第２の周波数シフト信号に従って、当該パルス光の光周波数ｆ０を光周波数ｆ０＋ｆＩＦ２にシフトすることによって、光周波数ｆ０＋ｆＩＦ２を有しているパルス光Ｐ２を生成する。
第２の変調部５ｃは、パルス光Ｐ２を光合波部５ｄに出力する。
【００９４】
光合波部５ｄは、第１の変調部５ｂからパルス光Ｐ１を受けると、パルス光Ｐ１を送信側光学系７に出力する。
光合波部５ｄは、第２の変調部５ｃからパルス光Ｐ２を受けると、パルス光Ｐ２を送信側光学系７に出力する。
【００９５】
次に、実施の形態２に係るライダ装置の動作について説明する。パルス変調器５以外は、図１に示すライダ装置と同様であるため、ここでは、主にパルス変調器５の動作について説明する。
光分岐部５ａは、光分割部３から、連続光であるレーザ光を受けると、レーザ光を２分岐する。
光分岐部５ａは、２分岐後の一方のレーザ光を第１の変調部５ｂに出力し、２分岐後の他方のレーザ光を第２の変調部５ｃに出力する。
【００９６】
第１の変調部５ｂは、パルス信号生成部４ａからパルス幅ΔＴのパルス信号を受ける毎に、当該パルス信号に従って、光分岐部５ａから出力されたレーザ光をパルス変調することで、当該レーザ光を、パルス幅ΔＴを有するパルス光に変換する。
第１の変調部５ｂは、トリガ生成部４のスイッチ４ｅから第１の周波数シフト信号を受けると、第１の周波数シフト信号に従って、当該パルス光の光周波数ｆ０を光周波数ｆ０＋ｆＩＦ１にシフトすることによって、光周波数ｆ０＋ｆＩＦ１を有しているパルス光Ｐ１を生成する。
第１の変調部５ｂは、パルス光Ｐ１を光合波部５ｄに出力する。
【００９７】
第２の変調部５ｃは、パルス信号生成部４ａからパルス幅ΔＴのパルス信号を受ける毎に、当該パルス信号に従って、光分岐部５ａから出力されたレーザ光をパルス変調することで、当該レーザ光を、パルス幅ΔＴを有するパルス光に変換する。
トリガ生成部４のスイッチ４ｅは、第１の周波数シフト信号を第１の変調部５ｂに出力してから、周期Ｔｒｅｐの時間が経過すると、第２の周波数シフト信号を第２の変調部５ｃに出力する。
第２の変調部５ｃは、トリガ生成部４のスイッチ４ｅから第２の周波数シフト信号を受けると、第２の周波数シフト信号に従って、当該パルス光の光周波数ｆ０を光周波数ｆ０＋ｆＩＦ２にシフトすることによって、光周波数ｆ０＋ｆＩＦ２を有しているパルス光Ｐ２を生成する。
第２の変調部５ｃは、パルス光Ｐ２を光合波部５ｄに出力する。
【００９８】
光合波部５ｄは、第１の変調部５ｂからパルス光Ｐ１を受けると、パルス光Ｐ１を送信側光学系７に出力する。
光合波部５ｄは、第１の変調部５ｂからパルス光Ｐ１を受けたのち、周期Ｔｒｅｐの時間が経過することで、第２の変調部５ｃからパルス光Ｐ２を受けると、パルス光Ｐ２を送信側光学系７に出力する。
【００９９】
パルス変調器５が、光分岐部５ａ、第１の変調部５ｂ、第２の変調部５ｃ及び光合波部５ｄを備えているライダ装置は、図１に示すライダ装置と同様に、遠方散乱光と近傍散乱光とが、互いに重なり合った状態であっても、遠方の観測対象の移動速度と近傍の観測対象の移動速度とが同じ速度であるか否かにかかわらず、それぞれの観測対象の移動速度を算出することができる。
【０１００】
実施の形態３．
実施の形態３では、スキャナ６１及び切替速度制御部６２を備えるライダ装置について説明する。
【０１０１】
図１４は、実施の形態３に係るライダ装置を示す構成図である。図１４において、図１と同一符号は同一又は相当部分を示すので説明を省略する。
スキャナ６１は、送受信部６のテレスコープ９から放射されるそれぞれのパルス光Ｐ１，Ｐ２の放射方向を切り替えるものである。
即ち、スキャナ６１は、テレスコープ９から放射されるそれぞれのパルス光Ｐ１，Ｐ２の光軸を時間的に掃引することで、それぞれのパルス光Ｐ１，Ｐ２の走査範囲をスキャンするものである。
切替速度制御部６２は、スキャナ６１による放射方向の切替速度を制御する。
即ち、切替速度制御部６２は、送受信部６により放射されるそれぞれのパルス光Ｐ１，Ｐ２がアイセーフ条件を満足するように、パルス光Ｐ１とパルス光Ｐ２との繰り返し周波数に応じて切替速度を制御する。
【０１０２】
次に、図１４に示すライダ装置の動作について説明する。スキャナ６１及び切替速度制御部６２以外は、図１に示すライダ装置と同様であるため、ここでは、スキャナ６１及び切替速度制御部６２の動作についてのみ説明する。
図１４に示すライダ装置が、アイセーフ条件を満足するように、それぞれのパルス光Ｐ１，Ｐ２を放射する場合、それぞれのパルス光Ｐ１，Ｐ２のパワーと、それぞれのパルス光Ｐ１，Ｐ２における走査範囲のスキャン速度とが制限される。
したがって、図１４に示すライダ装置から放射されるレーザ光の周期Ｔｒｅｐを短くすれば、それぞれのパルス光Ｐ１，Ｐ２における走査範囲のスキャン速度Ｖｓｃａｎを速めることが可能である。一方、図１４に示すライダ装置から放射されるレーザ光の周期Ｔｒｅｐを長くすれば、それぞれのパルス光Ｐ１，Ｐ２における走査範囲のスキャン速度Ｖｓｃａｎを遅くする必要がある。
【０１０３】
ライダ装置が、アイセーフ条件を満足しているか否かを判定する指標として、パルスレーザの被ばく放出限界ＡＥＬが用いられることがある（例えば非特許文献１を参照）。
［非特許文献１］
“レーザ安全ガイドブック”、光産業技術振興協会編、２００６年、新技術コミュニケーションズ出版
【０１０４】
図１４に示すライダ装置から放射されるそれぞれのパルス光Ｐ１，Ｐ２の波長は、例えば、１．５μｍ帯である。
波長が１．５μｍ帯のパルスレーザの場合、単一パルスに対する露光量は、８×１０−３［Ｊ］である。
放出持続時間が３秒以下である、波長が１．５μｍ帯のパルスレーザが、出射するパルスレーザを走査する場合、Φが１［ｍｍ］の円形開口絞りを横切る時間によって、時間基準Ｔが定義される。
時間基準Ｔは、以下の式（５）のように表される。
【０１０５】
0006980163.tif 000002
【０１０６】
時間基準Ｔに含まれるパルスレーザのパルス総数Ｔｏｔａｌは、Ｔｏｔａｌ＝ｆｒｅｐ×Ｔによって表される。ｆｒｅｐ＝１／Ｔｒｅｐである。
また、１パルス当りのエネルギーＥｐｕｌｓｅ＝Ｐｏｗｅｒ×Ｔｒｅｐは、以下の式（６）に示すように、パルス総数Ｔｏｔａｌに依存する補正係数Ｋが乗算された単一パルスに対する被ばく放出限界ＡＥＬｓを超えてはならない。補正係数Ｋは既値である。
Ｅｐｕｌｓｅ  ＜  ＡＥＬｓ×Ｋ            （６）
したがって、図１４に示すライダ装置から放射されるそれぞれのパルス光Ｐ１，Ｐ２がアイセーフ条件を満足するには、繰り返し周波数ｆｒｅｐとスキャン速度Ｖｓｃａｎとが、以下の関係式（７）を満足している必要がある。
【０１０７】
0006980163.tif 000003
【０１０８】
切替速度制御部６２は、繰り返し周波数ｆｒｅｐを示す周波数情報を取得し、周波数情報が示す繰り返し周波数ｆｒｅｐと補正係数Ｋとを関係式（７）に代入することによって、アイセーフ条件を満足するスキャン速度Ｖｓｃａｎを算出する。
切替速度制御部６２は、スキャン速度Ｖｓｃａｎを実現するためのスキャナ駆動信号を生成し、スキャナ駆動信号をスキャナ６１に出力する。
図１４に示すライダ装置から放射されるパルス光の繰り返し周波数ｆｒｅｐが下げられれば、切替速度制御部６２は、スキャン速度Ｖｓｃａｎを速めるためのスキャナ駆動信号を生成する。
図１４に示すライダ装置から放射されるパルス光の繰り返し周波数ｆｒｅｐが上げられれば、切替速度制御部６２は、スキャン速度Ｖｓｃａｎを遅くするためのスキャナ駆動信号を生成する。
【０１０９】
スキャナ６１のスキャン速度Ｖｓｃａｎは、切替速度制御部６２から出力されたスキャナ駆動信号によって制御される。
スキャナ６１は、送受信部６のテレスコープ９から放射されるそれぞれのパルス光Ｐ１，Ｐ２の放射方向をスキャン速度Ｖｓｃａｎで切り替える。
【０１１０】
以上の実施の形態３では、送受信部６から放射されるパルス光の放射方向を切り替えるスキャナ６１と、送受信部６から放射されるパルス光がアイセーフ条件を満足するように、パルス光の繰り返し周波数に応じて、スキャナ６１による放射方向の切替速度を制御する切替速度制御部６２とを備えるように、図１４に示すライダ装置を構成した。したがって、図１４に示すライダ装置は、図１に示すライダ装置と同様に、遠方散乱光と近傍散乱光とが、互いに重なり合った状態であっても、遠方の観測対象の移動速度と近傍の観測対象の移動速度とが同じ速度であるか否かにかかわらず、それぞれの観測対象の移動速度を算出することができる。また、図１４に示すライダ装置は、パルス光がアイセーフ条件を満足する状態で、パルス光の放射方向を切り替えることができる。
【０１１１】
実施の形態４．
実施の形態４では、観測対象が気体であり、気体の吸収波長帯に含まれる波長を有する第１のレーザ光と、第１のレーザ光よりも気体による吸収率が低い第２のレーザ光とを順番に出力する光源７１を備えるライダ装置について説明する。
【０１１２】
実施の形態４では、ライダ装置を設置している気中内に観測対象の気体であるガス（以下「観測対象ガス」という）が存在しているものものとする。ただし、これは一例に過ぎず、ライダ装置を設置している気中内に観測対象ガスがなく、窓等を介して別空間に観測対象ガスが存在しているものであってもよい。
観測対象ガスとしては、例えば、大気中の構成分子が該当する。大気中の構成分子は、窒素、酸素、二酸化炭素、又は、水蒸気である。大気中の構成分子には、窒素酸化物（ＮＯｘ）のような大気汚染物質も含まれる。ライダ装置において、気体の吸収波長帯は既値である。観測対象ガスは、散乱体を含んでいる。散乱体としては、雲、煙、塵、エアロゾル、又は、雨滴等が該当する。
実施の形態４に係る信号処理装置１５は、光源７１が第１のレーザ光を出力したときの検波信号の周波数と、光源７１が第２のレーザ光を出力したときの検波信号の周波数とから、気体の密度を算出する密度算出部７２を備えている。
【０１１３】
図１５は、実施の形態４に係る信号処理装置１５を含むライダ装置を示す構成図である。図１５において、図１と同一符号は同一又は相当部分を示すので説明を省略する。
図１６は、実施の形態４に係る信号処理装置１５を示す構成図である。図１７は、実施の形態４に係る信号処理装置１５のハードウェアを示すハードウェア構成図である。
図１６及び図１７において、図３及び図４と同一符号は同一又は相当部分を示すので説明を省略する。
光源７１は、例えば、単一周波数のレーザ光を発光するレーザであって、発光スペクトルの線幅が数ＭＨｚ以下の、半導体レーザ、ファイバレーザ、又は、固体レーザによって実現される。あるいは、光源７１は、半導体レーザ、ファイバレーザ及び固体レーザのうち、１つ以上のレーザの組み合わせによって実現される。
光源７１は、気体の吸収波長帯に含まれる波長を有する第１のレーザ光と、第１のレーザ光よりも気体による吸収率が低い第２のレーザ光とをパルス変調部２の光分割部３に順番に出力する。
【０１１４】
密度算出部７２は、例えば、図１７に示す密度算出回路８１によって実現される。
密度算出部７２は、光源７１から第１のレーザ光が出力されたときに、光検波部１１から出力されたそれぞれの検波信号Ｄ１〜Ｄ４の周波数を解析し、光源７１から第２のレーザ光が出力されたときに、光検波部１１から出力されたそれぞれの検波信号Ｄ１〜Ｄ４の周波数を解析する。
密度算出部７２は、それぞれの周波数の解析結果から、気体の密度ｄｅｎを算出する。
【０１１５】
図１６では、信号処理装置１５の構成要素であるドップラー周波数算出部１６、速度算出部１７、ＳＮＲ算出部３０、距離特性算出部３１及び密度算出部７２のそれぞれが、図１７に示すような専用のハードウェアによって実現されるものを想定している。即ち、信号処理装置１５が、ドップラー周波数算出回路４１、速度算出回路４２、ＳＮＲ算出回路４３、距離特性算出回路４４及び密度算出回路８１によって実現されるものを想定している。
ドップラー周波数算出回路４１、速度算出回路４２、ＳＮＲ算出回路４３、距離特性算出回路４４及び密度算出回路８１のそれぞれは、例えば、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ＡＳＩＣ、ＦＰＧＡ、又は、これらを組み合わせたものが該当する。
【０１１６】
信号処理装置１５の構成要素は、専用のハードウェアによって実現されるものに限るものではなく、信号処理装置１５が、ソフトウェア、ファームウェア、又は、ソフトウェアとファームウェアとの組み合わせによって実現されるものであってもよい。
信号処理装置１５が、ソフトウェア又はファームウェア等によって実現される場合、ドップラー周波数算出部１６、速度算出部１７、ＳＮＲ算出部３０、距離特性算出部３１及び密度算出部７２におけるそれぞれの処理手順をコンピュータに実行させるためのプログラムが図５に示すメモリ５１に格納される。そして、図５に示すプロセッサ５２がメモリ５１に格納されているプログラムを実行する。
【０１１７】
また、図１７では、信号処理装置１５の構成要素のそれぞれが専用のハードウェアによって実現される例を示し、図５では、信号処理装置１５がソフトウェア又はファームウェア等によって実現される例を示している。しかし、これは一例に過ぎず、信号処理装置１５における一部の構成要素が専用のハードウェアによって実現され、残りの構成要素がソフトウェア又はファームウェア等によって実現されるものであってもよい。
【０１１８】
図１５に示すライダ装置の動作について説明する。光源７１及び密度算出部７２以外は、図１に示すライダ装置と同様であるため、ここでは、主に、光源７１及び密度算出部７２の動作について説明する。
まず、光源７１は、気体の吸収波長帯に含まれる波長を有する第１のレーザ光をパルス変調部２の光分割部３に出力する。光源７１から出力される第１のレーザ光が有する光周波数は、ｆ０Ａである。
光検波部１１は、光源７１から第１のレーザ光が出力されているとき、実施の形態１と同様に、送受信部６により受信されたそれぞれの散乱光Ｒ１，Ｒ２，Ｒ３，Ｒ４と光源１から出力されたレーザ光である参照光との合波光を検波する。
ここでは、説明の便宜上、光検波部１１により検波されるそれぞれの合波光を、Ｃ１Ａ，Ｃ２Ａ，Ｃ３Ａ，Ｃ４Ａのように表記する。
光検波部１１は、それぞれの合波光Ｃ１Ａ，Ｃ２Ａ，Ｃ３Ａ，Ｃ４Ａの検波信号Ｄ１Ａ，Ｄ２Ａ，Ｄ３Ａ，Ｄ４Ａを信号処理装置１５に出力する。
【０１１９】
レンジビン分割部２２、周波数解析部２３、レンジ補正部２４、周波数補正処理部２５及びスペクトル積算部２６のそれぞれは、実施の形態１と同様に動作する。
スペクトル積算部２６は、積算後の周波数スペクトルを密度算出部７２に出力する。
ここでは、説明の便宜上、スペクトル積算部２６から出力される積算後の周波数スペクトルをΣＨＦＳＡのように表記する。
【０１２０】
次に、光源７１は、第１のレーザ光よりも気体による吸収率が低い第２のレーザ光をパルス変調部２の光分割部３に出力する。光源７１から出力される第２のレーザ光が有する光周波数は、ｆ０Ｂである。
図１８は、気体の吸収波長帯、光源７１から出力された第１のレーザ光が有する波長及び光源７１から出力された第２のレーザ光が有する波長を示す説明図である。
図１８において、横軸は波長であり、縦軸は気体に対するレーザ光の透過率である。
破線は、光源７１から出力された第１のレーザ光を示しており、一点鎖線は、光源７１から出力された第２のレーザ光を示している。実線は、気体の吸収波長帯を示している。第１の波長は、気体の吸収波長帯に含まれており、第２の波長は、気体の吸収波長帯に含まれていない。気体による第２の波長の透過率が、気体による第１の波長の透過率よりも大きい。即ち、気体による第２の波長の吸収率が、気体による第１の波長の吸収率よりも小さい。
【０１２１】
光検波部１１は、光源７１から第２のレーザ光が出力されているとき、実施の形態１と同様に、送受信部６により受信されたそれぞれの散乱光Ｒ１，Ｒ２，Ｒ３，Ｒ４と光源１から出力されたレーザ光である参照光との合波光を検波する。
ここでは、説明の便宜上、光検波部１１により検波されるそれぞれの合波光を、Ｃ１Ｂ，Ｃ２Ｂ，Ｃ３Ｂ，Ｃ４Ｂのように表記する。
光検波部１１は、それぞれの合波光Ｃ１Ｂ，Ｃ２Ｂ，Ｃ３Ｂ，Ｃ４Ｂの検波信号Ｄ１Ｂ，Ｄ２Ｂ，Ｄ３Ｂ，Ｄ４Ｂを信号処理装置１５に出力する。
【０１２２】
密度算出部７２は、光源７１から第１のレーザ光が出力されたときに、光検波部１１から出力されたそれぞれの検波信号Ｄ１Ａ〜Ｄ４Ａの周波数を解析し、光源７１から第２のレーザ光が出力されたときに、光検波部１１から出力されたそれぞれの検波信号Ｄ１〜Ｄ４の周波数を解析する。
密度算出部７２は、それぞれの周波数の解析結果から、気体の密度ｄｅｎを算出する。
以下、密度算出部７２による密度ｄｅｎの算出処理を具体的に説明する。
【０１２３】
密度算出部７２は、光周波数補正部２１から、光源７１から第１のレーザ光が出力されたときの、それぞれのスペクトル強度ＦＳｍａｘ１，ＦＳｍａｘ２を取得する。ここでは、説明の便宜上、密度算出部７２により取得されるそれぞれのスペクトル強度を、ＦＳｍａｘ１Ａ（ｎ１Ａ），ＦＳｍａｘ２Ａ（ｎ２Ａ）のように表記する。ｎ１Ａは、スペクトル強度ＦＳｍａｘ１Ａのレンジビンを示し、ｎ２Ａは、スペクトル強度ＦＳｍａｘ２Ａのレンジビンを示している。
また、密度算出部７２は、光周波数補正部２１から、光源７１から第２のレーザ光が出力されたときの、スペクトル強度ＦＳｍａｘ１，ＦＳｍａｘ２を取得する。ここでは、説明の便宜上、密度算出部７２により取得されるそれぞれのスペクトル強度を、ＦＳｍａｘ１Ｂ（ｎ１Ｂ），ＦＳｍａｘ２Ｂ（ｎ２Ｂ）のように表記する。ｎ１Ｂは、スペクトル強度ＦＳｍａｘ１Ｂのレンジビンを示し、ｎ２Ｂは、スペクトル強度ＦＳｍａｘ２Ｂのレンジビンを示している。
【０１２４】
図１６に示す信号処理装置１５では、密度算出部７２が、光周波数補正部２１からスペクトル強度ＦＳｍａｘ１Ａ（ｎ１Ａ），ＦＳｍａｘ２Ａ（ｎ２Ａ），ＦＳｍａｘ１Ｂ（ｎ１Ｂ），ＦＳｍａｘ２Ｂ（ｎ２Ｂ）を取得している。しかし、これは一例に過ぎず、密度算出部７２が、光周波数補正部２１と同様の光周波数補正部を備えていれば、当該光周波数補正部からスペクトル強度ＦＳｍａｘ１Ａ（ｎ１Ａ），ＦＳｍａｘ２Ａ（ｎ２Ａ），ＦＳｍａｘ１Ｂ（ｎ１Ｂ），ＦＳｍａｘ２Ｂ（ｎ２Ｂ）を取得するようにしてもよい。
【０１２５】
密度算出部７２は、スペクトル強度ＦＳｍａｘ１Ａ（ｎ１Ａ），ＦＳｍａｘ２Ａ（ｎ２Ａ），ＦＳｍａｘ１Ｂ（ｎ１Ｂ），ＦＳｍａｘ２Ｂ（ｎ２Ｂ）を以下の式（８）に代入することによって、気体の密度ｄｅｎを算出する。
【０１２６】
0006980163.tif 000004
式（８）において、ｋＯＮは、第１のレーザ光が有している波長の吸収係数であり、既値の係数である。ｋＯＦＦは、第２のレーザ光が有している波長の吸収係数であり、既値の係数である。ｌｎは、底がｅの対数関数を示す数学記号である。
【０１２７】
以上の実施の形態４では、観測対象が気体であり、光源７１が、気体の吸収波長帯に含まれる波長を有する第１のレーザ光と、第１のレーザ光よりも気体による吸収率が低い第２のレーザ光とを順番に出力する。光源７１から第１のレーザ光が出力されたときに、光検波部１１から出力されたそれぞれの検波信号の周波数を解析し、光源７１から第２のレーザ光が出力されたときに、光検波部１１から出力されたそれぞれの検波信号の周波数を解析し、それぞれの周波数の解析結果から、気体の密度を算出する密度算出部７２を備えるように、図１５に示すライダ装置を構成した。したがって、図１５に示すライダ装置は、図１に示すライダ装置と同様に、遠方散乱光と近傍散乱光とが、互いに重なり合った状態であっても、遠方の観測対象の移動速度と近傍の観測対象の移動速度とが同じ速度であるか否かにかかわらず、それぞれの観測対象の移動速度を算出することができる。また、図１５に示すライダ装置は、観測対象である気体の密度を算出することができる。
【０１２８】
なお、本開示は、各実施の形態の自由な組み合わせ、あるいは各実施の形態の任意の構成要素の変形、もしくは各実施の形態において任意の構成要素の省略が可能である。
【産業上の利用可能性】
【０１２９】
本開示は、信号処理装置、信号処理方法及びライダ装置に適している。
【符号の説明】
【０１３０】
１  光源、２  パルス変調部、３  光分割部、４  トリガ生成部、４ａ  パルス信号生成部、４ｂ  基準信号生成部、４ｃ  第１のシフト信号生成部、４ｄ  第２のシフト信号生成部、４ｅ  スイッチ、５  パルス変調器、５ａ  光分岐部、５ｂ  第１の変調部、５ｃ  第２の変調部、５ｄ  光合波部、６  送受信部、７  送信側光学系、８  送受分離部、９  テレスコープ、１０  受信側光学系、１１  光検波部、１２  光合波部、１３  受光部、１４  Ａ／Ｄ変換部、１５  信号処理装置、１６  ドップラー周波数算出部、１７  速度算出部、２１  光周波数補正部、２２  レンジビン分割部、２３  周波数解析部、２４  レンジ補正部、２５  周波数補正処理部、２６  スペクトル積算部、２７  ピーク周波数検出部、２８  周波数算出処理部、２９  速度算出処理部、３０  ＳＮＲ算出部、３１  距離特性算出部、４１  ドップラー周波数算出回路、４２  速度算出回路、４３  ＳＮＲ算出回路、４４  距離特性算出回路、５１  メモリ、５２  プロセッサ、６１  スキャナ、６２  切替速度制御部、７１  光源、７２  密度算出部、８１  密度算出回路。
【要約】
信号処理装置（１５）が、空間に存在している複数の観測対象におけるそれぞれの移動速度として、ライダ装置に対するそれぞれの観測対象の相対速度を算出するように構成した。光源（１）から出力されたレーザ光から、光周波数が互いに異なる複数のパルス光が生成されて、それぞれのパルス光が空間に放射され、それぞれの観測対象によって散乱された後のそれぞれのパルス光が散乱光として受信されて、それぞれの散乱光とレーザ光との合波光が検波されている。信号処理装置（１５）は、生成された複数のパルス光が有している光周波数と、それぞれの合波光の検波信号とから、それぞれの散乱光が有している光周波数に含まれている、それぞれの観測対象の移動に伴うドップラー周波数を算出するドップラー周波数算出部（１６）と、ドップラー周波数算出部（１６）により算出されたそれぞれのドップラー周波数から、それぞれの観測対象の相対速度を算出する速度算出部（１７）とを備えている。【図１】
0006980163.tif 000005
【図２】
0006980163.tif 000006
【図３】
0006980163.tif 000007
【図４】
0006980163.tif 000008
【図５】
0006980163.tif 000009
【図６】
0006980163.tif 000010
【図７】
0006980163.tif 000011
【図８】
0006980163.tif 000012
【図９】
0006980163.tif 000013
【図１０】
0006980163.tif 000014
【図１１】
0006980163.tif 000015
【図１２】
0006980163.tif 000016
【図１３】
0006980163.tif 000017
【図１４】
0006980163.tif 000018
【図１５】
0006980163.tif 000019
【図１６】
0006980163.tif 000020
【図１７】
0006980163.tif 000021
【図１８】
0006980163.tif 000022