(57)【特許請求の範囲】
【請求項１】
測定対象が放射する電波の周波数帯域に対応して当該電波を受信する単一のアンテナと、
前記アンテナで電圧として受信されたマイクロ波電力を所定のサンプリング周波数でデジタル化するＡＤ変換器と、
デジタル信号に変換された前記マイクロ波電力に対してフーリエ変換を行い、当該マイクロ波電力を周波数成分に分解すると共に周波数成分毎に複素フーリエ係数を積算して輝度温度ＴBを導き出す信号処理部と、
を備えたことを特徴とするマイクロ波放射計。
【請求項２】
前記信号処理部はリップル補正部を備え、
前記リップル補正部は計測された輝度温度ＴBのスペクトルを補正するためのリップルパターンを備え、当該リップルパターンを用いて計測された輝度温度ＴBを補正することを特徴とする請求項１記載のマイクロ波放射計。
【請求項３】
前記ＡＤ変換器及び前記信号処理部からなる処理系統を複数備えると共に、前記アンテナで受信したマイクロ波電力を複数の周波数帯に分配する分配器を備え、前記分配器の出力を各処理系統に入力したことを特徴とする請求項１記載のマイクロ波放射計。
【請求項４】
前記ＡＤ変換器でのサンプリング周波数を２のべき乗に相当する周波数で行うことを特徴とする請求項１記載のマイクロ波放射計。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、環境中のあらゆる物質から放射される自然由来のマイクロ波電力を計測するマイクロ波放射計に関する。
【背景技術】
【０００２】
例えば、人工衛星や航空機などから地球の表面付近を観測するリモートセンシング分野では、マイクロ波放射計を用いて環境中のマイクロ波帯の電磁波を受信し、受信した電磁波の電力から物体の輝度温度を観測した後に、当該輝度温度から計測対象の温度や総量を推定することが行われている。マイクロ波は、赤外線よりもさらに波長の長い電磁波であり、微小な水滴によって減衰しにくい。このため、前記マイクロ波放射計を用いれば、適切な波長を選択することにより、赤外線では観測できなかった雲下の海面水温や土壌水分量等の物理量を観測可能である。
【０００３】
図５は従来公知の全電力マイクロ波放射計の受信器の構成を示すブロック図である。所望の周波数帯域Ｂ（数１００ＭＨｚ〜１ＧＨｚ程度）のマイクロ波電力は狭帯域アンテナ１０で受信され、低雑音受信機（ＬＮＡ）で増幅された後に、バンドパスフィルタ（ＢＰＦ）で特定の帯域幅に切り出されて、更にダイオード検波器（Ｄｉｏｄｅ  Ｄｅｔｅｃｔｏｒ）で検出される。その後、検出された電力は時間積分（∫）され、求められたカウント値（Ｃ）が輝度温度（ＴB）に校正される。図示はされてないが、前記低雑音受信機（ＬＮＡ）とバンドパスフィルタ（ＢＰＦ）の間に周波数変換部が設けられ、必要に応じて周波数変換がなされる場合もある。
【０００４】
環境中に存在する自然物由来のマイクロ波電力は微弱である一方、近年では無線通信分野におけるマイクロ波帯の活用が進んでいる。このため、計測対象である周波数帯域Ｂの近傍の周波数帯域で無線通信が行われた場合に、人工電波による干渉（ＲＦＩ）が混入し、輝度温度を正しく計測できないといったリスクが懸念される。計測された輝度温度がＲＦＩで汚染されている場合、そこから汚染前の輝度温度を復元することは困難であり、マイクロ波を利用したリモートセンシングに深刻な問題である。
【０００５】
一方、人工的に放射されるマイクロ波の帯域幅はせいぜい数１０ＭＨｚである。このため、輝度温度を計測する帯域幅を狭くし、可及的に広い周波数範囲にわたって連続的に多数の輝度温度を計測すれば、前記ＲＦＩで汚染された周波数帯域の識別が可能になる。そして、図６に示すように、前記ＲＦＩで汚染された周波数帯域の輝度温度を取り除き、当該周波数帯域の前後に存在する汚染されていない周波数帯域の輝度温度で補間することにより、汚染された周波数帯域についてもＲＦＩ汚染のない正しい輝度温度を得ることができる。
【０００６】
しかし、図５に示された従来公知のマイクロ波放射計では、アンテナ１０、バンドパスフィルタ（ＢＰＦ）及びダイオード検波器（Ｄｉｏｄｅ  Ｄｅｔｅｃｔｏｒ）を含む受信機一式につき、周波数帯域Ｂについて一つの輝度温度を計測できるのみである。このため、複数の周波数帯域についてそれぞれの輝度温度を計測するためには、周波数帯域毎に受信機が必要となるが、マイクロ波放射計のシステム全体で搭載可能な受信機の数は搭載スペースの制限により十数程度であり、計測可能な周波数帯域の個数には制約があった。
【０００７】
単一のアンテナを用いながら複数の周波数帯域の夫々について輝度温度を計測する手段としては、例えば非特許文献１や非特許文献２に示されるように、アンテナによって広い周波数帯域のマイクロ波電力を受信する一方、この受信周波数帯域を連続的にカバーするように帯域幅の狭いバンドパスフィルタ（ＢＰＦ）を複数設け、各バンドパスフィルタ（ＢＰＦ）が切り出した周波数帯域毎に輝度温度を計測することが考えられる。しかし、マイクロ波放射計のシステム全体で計測可能な輝度温度の周波数帯域の個数は数１０個が限界であり、計測する輝度温度の周波数帯域の個数に比例して装置構成が大型化してしまうといった課題がある。
【０００８】
また、図５に示した従来のマイクロ波放射計のように高周波を受信する装置では、アンテナからダイオード検波器に向かう信号と、機器のアナログ信号経路内で反射して逆方向に向かう信号とが干渉して定在波が発生してしまい（図２（ａ）中の実線）、輝度温度のスペクトルに系統的誤差が発生する現象（リップル）が不可避である。この現象は輝度温度を周波数方向に連続的に取得することで初めて認識できるため、広い周波数帯域について輝度温度を一つ又は数個しか計測できない従来のマイクロ波放射計ではリップル現象の存在を識別することができず、計測された輝度温度に潜在的な系統誤差が含まれてしまう課題がある。
【先行技術文献】
【非特許文献】
【０００９】
【非特許文献１】Microwave radiometer for cloud carthography: A 22-channel ground-based microwave radiometer for atmospheric research
【非特許文献２】A multichannel radiometric profiler of temperature, humidity, and cloud liquid
【発明の概要】
【発明が解決しようとする課題】
【００１０】
本発明はこのような課題に鑑みなされたものであり、その目的とするところは、人工電波の干渉を除外して自然由来のマイクロ波を正しく計測することができると共に、リップルによる輝度温度の系統誤差を補正することができ、広い周波数帯域について正確に輝度温度を計測することが可能なマイクロ波放射計を提供することにある。
【課題を解決するための手段】
【００１１】
すなわち、本発明のマイクロ波放射計は、測定対象が放射する電波の周波数帯域に対応して当該電波を受信する単一のアンテナと、前記アンテナで電圧として受信されたマイクロ波電力を所定のサンプリング周波数でデジタル化するＡＤ変換器と、デジタル信号に変換された前記マイクロ波電力に対してフーリエ変換を行い、当該マイクロ波電力を周波数成分に分解すると共に周波数成分毎に複素フーリエ係数を積算して輝度温度ＴBを導き出す信号処理部と、を備えたことを特徴とするものである。
【発明の効果】
【００１２】
本発明によれば、輝度温度を計測する帯域幅を狭くしつつも、可及的に広い周波数範囲にわたって連続的に輝度温度を計測することが可能であり、人工電波の干渉を除外して自然由来のマイクロ波を正しく計測することができると共に、リップルによる輝度温度の系統誤差を補正することができ、広い周波数帯域について正確に輝度温度を計測することが可能となる。
【図面の簡単な説明】
【００１３】
【図１】本発明方法が適用されるマイクロ波放射計の第一実施形態を示すブロック図である。
【図２】リップル補正の前後における輝度温度のスペクトルを示す図である。
【図３】本発明方法が適用されるマイクロ波放射計の第二実施形態を示すブロック図である。
【図４】本発明方法が適用されるマイクロ波放射計の第三実施形態を示すブロック図である。
【図５】従来公知のマイクロ波放射計を示すブロック図である。
【図６】人工電波干渉（ＲＦＩ）によって汚染された輝度温度と周辺の輝度温度から補間したＲＦＩ汚染前の輝度温度との関係を示す図である。
【発明を実施するための形態】
【００１４】
以下、添付図面を用いて本発明のマイクロ波放射計を詳細に説明する。
【００１５】
図１は本発明を適用したマイクロ波放射計の第一実施形態を示すブロック図であり、本発明の基本的構成を示している。
【００１６】
このマイクロ波放射計は、測定対象が放射する電波を受信するアンテナ１と、前記アンテナ１で電圧として受信されたマイクロ波電力（以下、「受信信号」という）を増幅する低雑音受信機（ＬＮＡ）と、受信信号の周波数成分を特定の帯域幅で抽出するバンドパスフィルタ（ＢＰＦ）と、受信信号を所定のサンプリング周波数でデジタル化するＡＤ変換器（ＡＤＣ）と、ＡＤ変換器によってデジタル信号に変換された受信信号を周波数成分に分解すると共に周波数成分毎に輝度温度を導き出す信号処理部２と、を備えている。
【００１７】
前記アンテナ１としては観測する周波数帯域に最適な単一のアンテナを選択して使用することが可能であり、例えばマイクロ波のＣバンドからＫａバンドを受信できる超広帯域アンテナを選択することができる。具体的には、マイクロ波電力の水平・垂直偏波成分を同時に受信できるクアッドリッジフィードホーンを使用した。
【００１８】
前記バンドパスフィルタ（ＢＰＦ）としては、受信信号の周波数成分を帯域幅８１９２ＭＨｚで抽出するものを使用した。このバンドパスフィルタ（ＢＰＦ）は後段のＡＤ変換器（ＡＤＣ）での処理の際に生じる雑音成分を抑制する目的に特化したアンチエイリアシングフィルタ（ＡＡＦ）としても機能する。尚、前記バンドパスフィルタ（ＢＰＦ）の帯域幅はＡＤ変換機のサンプリング周波数に応じて設計変更することが可能であり、例えば帯域幅１６３８４ＭＨｚのバンドパスフィルタ（ＢＰＦ）を使用しても良い。
【００１９】
従って、前記アンテナ１から入力された受信信号は前記低雑音受信機（ＬＮＡ）で増幅された後、前記バンドパスフィルタ（ＢＰＦ）によって帯域制限され、更に前記ＡＤ変換器（ＡＤＣ）によってアナログ信号からデジタル信号に変換される。変換されたデジタル信号はアンテナで受信されたマイクロ波電力の時間波形を示しており、これが後段の信号処理部に入力される。尚、図示はされてないが、前記低雑音受信機（ＬＮＡ）とバンドパスフィルタ（ＢＰＦ）の間に周波数変換部を設け、必要に応じて周波数変換を行ってもよい。
【００２０】
一方、前記信号処理部２はＰＬＤの一種であるＦＰＧＡ（ｆｉｅｌｄ−ｐｒｏｇｒａｍｍａｂｌｅ  ｇａｔｅ  ａｒｒａｙ）によって構成されている。このＦＰＧＡには高速フーリエ変換（以下、「ＦＦＴ」という）が実装されており、前段のＡＤ変換機によってデジタル化された受信信号を各周波数成分に分解すると共に、各周波数成分の複素フーリエ係数を計算する。求められた複素フーリエ係数は周波数毎に積分部（∫）において積算され、各周波数成分の帯域幅に対応するカウント値（Ｃ1〜Ｃm）が出力された後、当該カウント値が周波数毎の輝度温度（ＴB1〜ＴBm）に校正されて外部に出力される。尚、前記信号処理部２は前述のＦＰＧＡ以外に、別な種類の集積回路（ＩＣ）や、ソフトウェアで実現することも可能である。
【００２１】
ここで、ｍは分解された周波数成分の数を表しており、ＦＦＴ点数の１／２の数値である。前記ＡＤ変換器のサンプリング周波数をｘ、ＦＦＴ点数を２mとすれば、ｘ／２の周波数範囲をｍ区間に分割した各々帯域幅についてマイクロ波の輝度温度ＴBを計測できる。例えば典型的な例として、前記ＡＤ変換器のサンプリング速度を１６３８４ＭＨｚ、ＦＦＴ点数を１０２４点とした場合、８１９２ＭＨｚの広い周波数範囲を１６ＭＨｚの狭い帯域幅で５１２個に分割し、５１２個の帯域幅のそれぞれについて輝度温度ＴBを同時に計測することが可能となる。
【００２２】
このようにして取得された５１２個の輝度温度ＴB（輝度温度ＴBのスペクトル）を分析すれば、人工電波の干渉（ＲＦＩ）で汚染された周波数帯域を識別し、当該帯域に対応する輝度温度ＴBを補正することで、自然由来のマイクロ波放射による輝度温度を正しく抽出することが可能となる。
【００２３】
以上説明してきたように、図１に示した本実施形態のマイクロ波放射計は単一のアンテナ１、単一のＡＤ変換器（ＡＤＣ）及び単一の信号処理用ＬＳＩ（例えばＦＰＧＡ）で実現可能であり、複数の周波数帯域の夫々について輝度温度を計測しつつも装置の小型化を図ることが可能である。
【００２４】
次に、本発明を適用したマイクロ波放射計の第二実施形態について説明する。この第二実施形態では第一実施形態のマイクロ波放射計と基本的な考え方は同じであるが、リップル現象の補正を行っている点において第一実施形態のマイクロ波放射計と異なる。
【００２５】
マイクロ波放射計におけるアナログ信号の経路内、すなわち前記アンテナ１から前記ＡＤ変換機（ＡＤＣ）までの経路内においては、受信信号の特定の周波数成分が反射して干渉を生じる場合がある。このような干渉が生じると、本来は受信信号から計測された輝度温度ＴBのスペクトルが図２（ｂ）中に実線で示されるように周波数方向に緩やかに変化するのに対し、図２（ａ）中に実線で示すように、計測された輝度温度ＴBのスペクトルが周波数方向に乱れてしまう。この現象がリップルである。
【００２６】
このリップルを補正するため、図３に示すように、第二実施形態のマイクロ波放射計では前記信号処理部２にリップル補正部を設けている。前記リップル補正部２には各マイクロ波放射計に固有のリップルパターンが予め格納されており、当該リップル補正部は前記リップルパターンを用いて計測によって得られた輝度温度ＴBのスペクトルを補正する。
【００２７】
前記リップルパターンは前記実施形態１で取得された輝度温度ＴBのスペクトルから生成することができる。具体的には、受信信号から得られた輝度温度の周波数方向のスペクトルを用いて低次の多項式近似式を導き出し、図２（ａ）に点線で示される回帰曲線を得る。既に説明したように、第１実施形態のマイクロ波放射計では、広い周波数範囲について狭い帯域幅毎に輝度温度ＴBを連続的に計測可能なので、これら輝度温度の実測値から回帰曲線を容易に求めることが可能である。
【００２８】
この回帰曲線を自然由来の輝度温度のスペクトルの近似値とし、計測されたリップルを含む輝度温度ＴBのスペクトルから差し引くことで、放射計の内部で発生したリップルのパターンを得ることが出来る。すなわち、基準時刻における輝度温度ＴBのスペクトルをＴB0（ｆ）、回帰曲線の関数をＭ（ｆ）とすると、リップルパターンＲ（ｆ）は以下のように表すことができる。
Ｒ（ｆ）＝ＴB0（ｆ）−Ｍ（ｆ）
【００２９】
そして、以降の時刻において計測された輝度温度TBのスペクトルから前記リップルパターンＲ（ｆ）を差し引くことで、図２（ｂ）に示すように、リップルの影響を取り除いた自然由来のマイクロ波による輝度温度ＴBのスペクトルを得ることができ、輝度温度をより正確に計測することが可能となる。
【００３０】
前記リップルは放射計内部のアナログ経路の反射による定在波に由来しているため、当該アナログ経路の構造に依存しており、前記リップルパターンは個々のマイクロ波放射計に固有のものである。このため、前記リップルパターンが経時的に大きく変動することはなく、前記リップル補正部によるリップルパターンの取得及び格納は、例えばマイクロ波の計測の開始時または一定の時間間隔で行えば十分である。
【００３１】
尚、前記リップル補正部に格納されるリップルパターンを得るための回帰曲線としては、前述の如く受信信号から得られた多数の計測値を利用して多項近似式を導き出しても良いが、マイクロ波放射計の測定対象の放射伝達モデルが解っている場合、その放射伝達モデルを利用して前記回帰曲線を生成してもよい。
【００３２】
次に、本発明を適用したマイクロ波放射計の第三実施形態について説明する。
【００３３】
アナログ信号をデジタル信号に変換する際、前記ＡＤ変換器（ＡＤＣ）のサンプリング周波数としてはアナログ信号の帯域幅の２倍よりも高い周波数を用いることが必要である。このため、前記アンテナ１で受信したマイクロ波の周波数範囲が著しく広く、前記ＡＤ変換器（ＡＤＣ）のサンプリング周波数ｘの１／２を超える場合は、受信信号を単一のＡＤ変換器（ＡＤＣ）でデジタル信号に変換することが困難となる。
【００３４】
広帯域の受信信号に対応するため、この第三実施形態のマイクロ波放射計では、図４に示すように、前記アンテナ１から入力された受信信号は低雑音受信機（ＬＮＡ）で増幅後に５系統に分割された後に、バンドパスフィルタ（ＢＰＦ１〜ＢＰＦ５）を介してＡＤ変換機（ＡＤＣ１〜ＡＤＣ５）に入力されており、この点において第一実施形態又は第二実施形態のマイクロ波放射計と相違している。各ＡＤ変換器（ＡＤＣ１〜ＡＤＣ５）の後段に設けられた信号処理部（ＤＭ１〜ＤＭ５）は前記第一実施形態又は第二実施形態の信号処理部と同一である。
【００３５】
図４に示すマイクロ波放射計では、アンテナ１から入力される受信信号の周波数帯は０〜４０ＧＨｚであり、８ＧＨｚの帯域幅をデジタル化可能なＡＤ変換器および信号処理部を５セット使用することによって０〜４０ＧＨｚの帯域幅を広げて有する受信信号の処理を可能としている。使用するＡＤ変換器は、０〜２４ＧＨｚの周波数帯のアナログ信号を入力することができ、このアナログ信号を１６ＧＨｚのサンプリング周波数でデジタル化することが可能である。広帯域アンテナ１から入力された０〜４０ＧＨｚの受信信号は、低雑音受信機（ＬＮＡ）で増幅された後に、先ずは分配器で０〜８ＧＨｚ、８〜１６ＧＨｚ、１６〜２４ＧＨｚ及び２４〜４０ＧＨｚの四つの周波数帯に分配される。
【００３６】
分配された四つの周波数帯のうち、帯域幅が８ＧＨｚの最初の3系統では、入力された受信信号がアンチエイリアシングフィルタ（ＡＡＦ）を兼ねるバンドパスフィルタ（ＢＰＦ１〜ＢＰＦ３）に入力され、更にＡＤ変換器（ＡＤＣ１〜ＡＤＣ３）によってデジタル化され、後段の信号処理部（ＤＭ１〜ＤＭ３）に入力される。
【００３７】
一方、４番目の系統が処理すべき周波数帯は２４〜４０ＧＨｚであって、前記ＡＤ変換器に入力可能な最高周波数を超えている。このため、４番目の系統に入力された受信信号は局部発振器及びミキサからなる周波数変換部を介して０〜１６ＧＨｚの周波数帯に周波数変換され、その後に分配器によって０〜８ＧＨｚ及び８〜１６ＧＨｚの２系統に分配される。この後、これらの受信信号はバンドパスフィルタ（ＢＰＦ４〜ＢＰＦ５）に入力され、更にＡＤ変換器（ＡＤＣ４〜ＡＤＣ５）によってデジタル化され、後段の信号処理部（ＤＭ４〜ＤＭ５）に入力される。
【００３８】
そして、各信号処理部（ＤＭ１〜ＤＭ５）では、前述の第一実施形態で説明したように、１０２４点のＦＦＴが実行されて５１２個の輝度温度ＴBに変換される。すなわち、この第三実施形態のマイクロ波放射計によれば、０〜４０ＧＨｚの広範囲な周波数帯におけるマイクロ波から約１６ＭＨｚの狭い帯域幅で合計２５６０個の輝度温度ＴBを連続的に計測することが可能となる。
【００３９】
尚、前述の第一実施形態乃至第三実施形態においては、デジタル化された受信信号に対して前記信号処理部で高速フーリエ変換（ＦＦＴ）を行っているが、ＦＦＴにおける周波数の分割点数は２のべき乗に制限される。このため、前記ＡＤ変換器から出力された受信信号の帯域全体を２のべき乗で分割することになり、一つ一つの輝度温度ＴBに対応する周波数は一般的に端数を含んで解りにくくなる。例えば、ＡＤ変換器のサンプリング周波数を１６０００ＭＨｚとした場合、取得される信号の周波数は第一ナイキスト領域において０〜８０００ＭＨｚとなり、これを１０２４点について高速フーリエ変換した場合、各輝度温度ＴBの周波数は以下のようになる。
８０００／５１２×ｎ（ＭＨｚ）  （但し、ｎ=０，１，２，…，５１１）
【００４０】
これを実際の周波数として表現すると、順番に０ＭＨｚ、１５．６２５ＭＨｚ、３１．２５ＭＨｚ、…、７９８４.３７５ＭＨｚとなり、直感的に把握しにくい結果となる。
【００４１】
そこで、前記ＡＤ変換器のサンプリング周波数を２のべき乗に相当する周波数、例えば１６３８４MHzとすることで、取得される信号の周波数は０〜８１９２ＭＨｚとなる。これを１０２４点について高速フーリエ変換した場合、各輝度温度ＴBの周波数は０、１６．０、３２．０、…、８１７６．０ＭＨｚとなり、直感的に把握し易い結果となる。このように前述の第一実施形態乃至第三実施形態のマイクロ波放射計を構築する際には、前記アンテナ１が受信するマイクロ波の周波数帯をカバーする範囲で、前記ＡＤ変換器（ＡＤＣ）のサンプリング周波数を２のべき乗に相当する周波数に設定するのが好ましい。
【００４２】
以上説明してきたように、本発明のマイクロ波放射計によれば、輝度温度を計測する帯域幅を狭くしつつも、可及的に広い周波数範囲にわたって連続的に多数の輝度温度を計測することが可能であり、ＲＦＩで汚染された周波数帯域を識別することが可能になると共に、当該周波数帯域の前後に存在する汚染されていない周波数帯域の輝度温度を用いることで、汚染された周波数帯域についてもＲＦＩ汚染のない正しい輝度温度を得ることが可能となる。
【００４３】
また、本発明によれば、可及的に広い周波数範囲にわたって連続的に多数の輝度温度を計測しつつも、単一のアンテナ及び単一の受信機でこれを実現することができ、人工衛星や航空機を用いたリモートセンシングに好適な小型のマイクロ波放射計を提供することが可能となる。
【００４４】
更に、本発明によれば、可及的に広い周波数範囲にわたって連続的に多数の輝度温度を計測することが可能なので、計測された輝度温度ＴBのスペクトルを補正するためのリップルパターンを容易に取得することができ、リップルの影響を取り除いて輝度温度をより正確に計測することが可能となる。
【要約】      （修正有）【課題】人工電波の干渉を除外して自然由来のマイクロ波を正しく計測することができると共に、リップルによる輝度温度の計測誤差を補正することができ、広い周波数帯域について正確に輝度温度を計測することが可能なマイクロ波放射計を提供する。【解決手段】測定対象が放射する電波を受信する単一のアンテナ１と、アンテナで電圧として受信されたマイクロ波電力を所定のサンプリング周波数でデジタル化するＡＤ変換器と、デジタル信号に変換されたマイクロ波電力に対してフーリエ変換を行い、マイクロ波電力を周波数成分に分解すると共に周波数成分毎に複素フーリエ係数を積算して輝度温度ＴBを導き出す信号処理部２と、を備えている。【選択図】図１【図１】
0006727600.tif 000002
【図２】
0006727600.tif 000003
【図３】
0006727600.tif 000004
【図４】
0006727600.tif 000005
【図５】
0006727600.tif 000006
【図６】
0006727600.tif 000007