音響電気相互作用を用いて物体中の物理的なパラメータを決定する装置および方法
【課題】音響電気相互作用を用いて物体中の物理的なパラメータを決定する装置および方法を提供する。
【解決手段】本発明は物体中の誘電関数を決定するための装置に関する。装置は、前記物体を通るマイクロ波放射を送信する１つの送信アンテナ４２と、前記送信されたマイクロ波放射を受信する１つの受信アンテナ４３と、前記物体中の密度変化を生成するために前記物体を通る超音波放射を射出する１つの超音波送信機と、前記物体中の前記音響電気相互作用δを決定するために前記密度変化を介して送信されるマイクロ波放射を解析する手段と、前記音響電気相互作用から前記物体中の誘電関数を計算する手段と、を備える。発明は、物体中の誘電関数を決定する方法にも関する。
【解決手段】本発明は物体中の誘電関数を決定するための装置に関する。装置は、前記物体を通るマイクロ波放射を送信する１つの送信アンテナ４２と、前記送信されたマイクロ波放射を受信する１つの受信アンテナ４３と、前記物体中の密度変化を生成するために前記物体を通る超音波放射を射出する１つの超音波送信機と、前記物体中の前記音響電気相互作用δを決定するために前記密度変化を介して送信されるマイクロ波放射を解析する手段と、前記音響電気相互作用から前記物体中の誘電関数を計算する手段と、を備える。発明は、物体中の誘電関数を決定する方法にも関する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、請求項１および請求項１９のプリアンブルに記載の物体の誘電関数を決定することによって、物質の内部の、温度または密度のような複数の物理的なパラメータを決定するための装置に関する。発明はまた、請求項１２のプリアンブルに記載の物体の内部の誘電関数を決定するための方法、および、請求項１８に記載の食品の中での温度の局所分布を決定するための装置に関する。
【背景技術】
【０００２】
物体を破壊することなく、物体に侵入せず、物体を切り裂くことなく、任意の物体の、温度、密度および他の複数の内部パラメータに関する情報を得るために、様々な型の放射は、所望の複数のパラメータを再構成することを可能にする情報を提供するのに利用可能である。
【０００３】
放射の特殊型式を選んで、分析法の選択において適切な予想結果を組込む４つの異なる場合がある。これらは、２つの問題領域によって分類される。
【０００４】
−選ばれた放射に対する物体の透明度
−選ばれた放射の波長に関して必要とされる物体での解像度
ケース１Ａ
（物体は、測定に用いられた放射に対して、透明であるか弱く吸収する。また、達成される解像度は、放射波長に等しいか、放射波長よりも小さい。）
ただ一つの情報源は、例えば、以下に述べるような手法を使用して近い場をプローブすることにより得られる。
【０００５】
−原子間力顕微鏡検査（ＡＦＭ）は、
検査中の物体の表面上の構造を読み出すために、物体の表面上に高精度で位置決めされるサブ波長サイズのステンシル上の力を読み出す。
【０００６】
−ラスタートンネル顕微鏡検査（ＲＴＭ）は、
物体の表面の電子状態の情報を生成するために、力の代わりに、検査中の物体の表面に接近して位置決めされるサブ波長サイズのプローブからトンネル電流を測定する。
【０００７】
−光学的近接場顕微鏡検査は、
電磁放射が、薄いプローブ上でのサブ波長解像度で光学的特性の表面画像を生成することを用いた放射の波長よりもはるかに小さい穴を必要とする微視的に小さな複数の穴を介して通り抜ける。
【０００８】
−インピーダンス断層撮影法は、
１セットの電極は検査中の物体に着けられ、また、すべてのプローブ間のインピーダンスが測定される。この方法は検査中の物体の内部のいくつかの特性を計算することを可能にする。しかし、解像度は一般に劣っている。この方法は、例えば抗凝固薬の影響を評価する前後で心臓の領域のインピーダンスを測定するなどの異なる複数の手法での成功で使用されてきている。
【０００９】
一般的な特徴として、上述した方法の高解像度は、選ばれた放射の固有波長に起因しているのではなく、サブ波長解像度を提供する別の拘束条件（横隔膜、ステンシルのようにほとんど機械的）に起因する。一般的な欠点は検査中の物体の厚さ要求によって与えられる−上記の複数の方法は、解像度を失わずに、非常に制限のある深さでの単に表面情報または内部情報のいずれかを生成する。
【００１０】
ケース１Ｂ
（物体は、測定に用いられる放射に対して、透明であるか、または弱く吸収することである。また、解像度は放射波長よりもはるかに大きい。）
このケースはすべての直接撮像および光伝送方法によってカバーされる。このレジーム（ｒｅｇｉｍｅ）で電磁放射を用いるとすると、以下のものがある。
【００１１】
−ライダー
−Ｘ線
分析の手段として、散乱が役割を果たさないので、レイトレーシングおよび１対１写像の方法が適切である。それは、放射源と受信機との間に位置する物体体積によって、ある与えられた位置で採用されたピクセル情報のそれぞれが影響するに過ぎないことを、解像度の損失なしで仮定することが可能である。
【００１２】
この領域の最近の進展はパッシブレーダーであり、受信機のまわりの環境ですべての身体に固有の熱放射が測定され撮像される。このレーダー方法は、どんな送信された信号をも必要とせず、したがって、跡をたどることができない。
【００１３】
非電磁気的な方法では、商業上利用可能な以下のものがある。
【００１４】
−超音波断層法
−核磁気共鳴（ＮＭＲ）
ケース２Ａ
（物体は、測定に用いられた放射に対して適度に吸収している。また、解像度は放射波長に等しいかそれよりも小さい。）
物体が測定に用いられた放射に対して適度に吸収しているという事実は、調査することが可能であるプローブに厚さ限界を設けることになる。
【００１５】
このケースについては、今日、最先端技術に関して利用可能な実現可能な方法はない。
【００１６】
ケース２Ｂ
（物体は、測定に用いられた放射に対して適度に吸収している。また、解像度は放射波長よりもはるかに大きい。）
この場合、ほとんどの無線周波数およびマイクロ波周波数の応用は見つかる（特に、検査中の物体が、損失があり、それが損失のない環境に組み込まれた場合）。また、マイクロ波断層撮影法が利用可能である。これらの方法で、最もよくあるものが次のものである。
【００１７】
−（アクティブ）レーダー（ＲＡＤＡＲ）
ソースとターゲットの間および受信機へ戻る信号のランタイムは、送信機（モノスタティックレーダー）と同じ場所に受信機を配置することによってか、もしくは、送信機（バイスタティックレーダー）と、ソースとターゲットの相対速度による周波数変化とが評価される（ドップラーレーダー）よりも異なる位置に受信機を配置することによって測定される。
【００１８】
このように、測定に用いられた放射に対して適度に吸収する物体にとって、温度、密度、組成のような物理的なパラメータを決定するための装置を開発する必要がある。所望の解像度は放射波長よりもはるかに大きい。
【発明の開示】
【課題を解決するための手段】
【００１９】
本発明の目的は、任意に形成された物体の誘電関数を決定するための装置を提供することである。
【００２０】
目的は、製品密度における制御可能な変化を作成するための超音波を用いて、請求項１および請求項１９の特徴部分で定義されるような装置、および請求項１２の特徴部分に記載の方法によって達成される。その後、装置は、密度変化を読出し、かつ誘電関数の空間分布にそれを関連づけるために、マイクロ波放射を用いる。これは、請求項１８の特徴部分に定義されるように、物体の温度、含水量および密度を決定するために一方で用いられる。
【００２１】
本発明での利点は、空間分布の解像度が、第１型放射（例えばマイクロ波放射)の波長に限定されず、むしろ第２型放射（例えば超音波、Ｘ線)の波長によって決定されることである。
【００２２】
本発明での別の利点は、温度、含水量などのような物理的特性の非接触測定が、複数の仮想プローブとして発明を適用して確立されてもよいということである。
【００２３】
他の複数の物体および複数の利点は、発明の詳細な記述からの当業者にとって明白になる。
【発明を実施するための最良の形態】
【００２４】
本発明に先立って、材料のみの内部特性を再構成するためのツールとして以下のものが存在する（ただし、回析および散乱は優勢である）。
【００２５】
−マイクロ波断層撮影法
−超音波断層法
両方のケースでは、解像度は使用された放射の波長によって決定される。
【００２６】
本発明では、超音波およびマイクロ波の方法は組み合わせられる。物体再構成は純粋なマイクロ波逆散乱方法、およびそれぞれの限定を有する純粋な超音波断層法方法によって行うことが可能である。ここで、超音波は、物体再構成ツールとして用いられず、調査される物体での密度変化を生成するツールとして用いられる。この密度変化は、物体再構成に用いられる送信されたマイクロ波放射線に、位相および周波数の変化を生じる。したがって、この方法の利用可能な解像度は、（典型的な超音波治療周波数用の１ミリメートルよりも小さな）超音波の解像度によって決定される。密度読み出しは、マイクロ波放射（ある周波数では、減衰はさらに合理的な複数の浸入度（例えばＳ、ＩＳＭ５．８、Ｘバンド）を許可する）を用いて行われる。この方法は、ミリメートル解像度がミリメートル波長を必要とするマイクロ波断層撮影法の複数のアプローチの基本的な困難を回避する。不運なことに、ミリメートル放射は、いくつかの波長内で興味深いほとんどの物体に吸収され、したがって、いかなる内部のパラメータも抽出されることは不可能である。上記の分類では、本発明は領域１Ｂ、２Ａおよび２Ｂをカバーする。そのような方法は本発明に先立って知られていない。
【００２７】
本発明によって記述されたシステムは、食品産業で用いられることが好適である。食品産業では、正確に食品の温度を制御することは多くの場合、重要である。例えば、食品が冷凍されることになっている場合、製品全体が冷凍されることが重要である。製品全体（例えば鶏フィレ）が冷凍されたことを保証することが可能でない場合、製品を廃棄しなければならないか、または短い保存寿命を備えた製品を配達しなければならない。したがって、製品の冷凍の非破壊および非接触の制御が必要である。この問題は、誘電関数を測定し、下記に記述されるようにそれを温度の分布に変換することによって解決されてもよい。
【００２８】
しかしながら、システムは決してこの種の産業に制限されているわけではない。見込みのある他の応用は以下がある。
−コンクリート硬化（建設業）
−接着剤硬化（飛行機製造）
−医療用画像（機能的な脳断層撮影法、脊髄の断層撮影法）
−地盤調査、トラッキングパイプおよび地中管
−救命および救助道具（瓦礫下の人を検知）
−掃海（大きくなる領域での特にプラスチックの機雷）
以下では、好ましい実施形態が要約される。上記の他の応用領域でこの方法を適応させするために形状に必要とされる修正は小さい。
【００２９】
以下では、持続波（ＣＷ）マイクロ波およびパルス波列超音波に基づいたシステムは、単純性の目的のために記述される。記述された方法はこのケースに制限されない。振幅変調（ＡＭ）、周波数変調（ＦＭ）、周波数変調持続波（ＦＭＣＷ）、パルス符号変調（ＰＣＭ）、位相変調（ＰＭ）およびウェーブレットに基づいた変調方式（ＷＭ）のような、両者のための他の変調スキーム、電磁波、および、超音波は適用可能であり、ある他の応用には最適である。
【００３０】
図１は発明による装置４０を記述する。システムはコンベア手段１１に接近して配置される。それは、センサ測定ギャップ１３を介して検査下の製品１２を運搬する。システム４０は、マイクロ波部５０、超音波部７０および評価ユニット６０から成る。システムは、本実施形態では、２固定周波数マイクロ波発生器５１および５２、ならびに固定周波数超音波発生器７１を具備する。第１マイクロ波発生器５１は第１の固定マイクロ波振動数ｆ１（例えば５．８１８ＧＨｚ）を持っており、少なくとも１つの送信アンテナ４２につながれる。また、第２マイクロ波発生器５２は第２の固定マイクロ波振動数ｆ２（例えば５．８ＧＨｚ）を持っており、好ましくはミキサのようなダウンコンバータ５４につながれる。ダウンコンバータは、送信される（少なくとも１つの受信アンテナ４３によって集められる）マイクロ波信号、および、第２マイクロ波発生器５２からの受信マイクロ波信号を、低い中間周波数にシフトする。これは、検査下の製品１２を通って送信されるマイクロ波信号が振幅および位相で評価されることを可能にする。それは更に、フィルタユニット５９、アナログ／デジタル変換器ＡＤＣ５５、１セットの信号処理部５６、および、システムを制御し、かつデータを評価するために必要なアルゴリズムを含んでいる評価プロセッサ６０を含む。結果は表示ユニット６５に提出される。システム４０はまた、送信アンテナ４２および受信アンテナ４３に加えて、測定ギャップ１３のまわりにすべてグループ化して、（明瞭さの目的のために１つだけ示されている）１セットのトランスデューサー７２を含む。トランスデューサーは、検査下の製品１２を通る、超音波周波数ｆＵＳ（例えば４．５ＭＨｚ）を持っている超音波信号を発する。これは、超音波速度で移動する密度変位を引き起こす。同時に、第１マイクロ波発生器５１からのマイクロ波信号が、送信アンテナ４２から射出される。この信号はまた、検査下の製品１２を通って移動する。マイクロ波信号は、マイクロ波周波数を保ったまま製品を通って移動することにより減衰および位相遅れを示す。検査下の製品１２の超音波が密度変位を生じるそれらの体積では、マイクロ波信号の一部分が周波数シフトされ、そして、上側波帯と下側波帯が生じる。透過したマイクロ波信号はマイクロ波受信アンテナ４３を用いて集められる。受信信号はダウンコンバータユニット５４を用いてダウンコンバートされる。その後、低周波信号は、フィルタユニット５９を用いてフィルタリングされ、ＡＤＣ５５を用いてアナログデジタル変換される。デジタル信号は受信信号処理部５６を用いて評価される。受信信号処理部は、標準の最新技術のデジタルフィルターを用いて、入ってくるデジタル信号をゼロ周波数に変換する。
【００３１】
このフィルタリングの結果は、当業者によく知られているように送信４２および受信４３アンテナとの間のＳ２１パラメータに相当する。これは周波数がシフトされていない。上記では、我々は、受信アンテナ４３をマイクロ波ポート２と称し、送信アンテナ４２をマイクロ波ポート１と称する。
【００３２】
本発明によって記述されるシステムでは、第１信号パス５９、５５、５６に平行して、第２セットのバンドパスフィルター５８、別のＡＤＣ５５および第２のデジタル信号処理部５７がある。
【００３３】
バンドパスフィルター５９は、両方のマイクロ波発生器５１および５２の間の差周波数に合わせられる。それは本実施形態では５．８１８ＧＨｚ−５．８ＧＨｚ＝１８ＭＨｚである。第２バンドパスフィルター５７は、マイクロ波発生器間の差周波数（例えば１８ＭＨｚ）に超音波信号発生器７１の中心周波数（例えば４．５ＭＨｚ）を加えた周波数に合わせられる。したがって、５８、５５および５７を収納しているこの第２のデジタル信号処理部のパスは、入ってくる信号を、超音波周波数によって周波数シフトされたゼロ周波数に変換する。したがって、測定結果は、超音波とマイクロ波信号との間の断面に制限されている。
【００３４】
第１の５９，５５，５６と第２の５８，５５，５７のデジタル受信機のＩＦ帯域幅は、超音波発生器７１によって生成された、超音波周波数ｆＵＳの半分に選ばれる。これは、超音波振動子の複数の位相を変えることによって、周波数シフトを最適化するために要求される。
【００３５】
製品１２の超音波メトリックを得る第１段階の間に、トランスデューサー７２から射出される超音波放射を収集する超音波受信機７３が存在して、より詳細に下に記述されるような減衰（Ｔ５６）およびランタイムを評価しなければならない。上記で、我々は、超音波受信機７３をマイクロ波ポート６と称し、トランスデューサー７２をマイクロ波ポート５と称する。減衰およびランタイムが、超音波評価ユニット７４で評価されるが、これは、評価ユニット６０に統合されても当然よい。
【００３６】
図２は、検査中の製品へ射出された放射を示す。トランスデューサー７２は、この例において、検査下の製品１２を通る超音波パルス９１を射出する。これは、超音波速度で移動する密度変位を引き起こす。同時に、マイクロ波信号９０は、送信アンテナ４２から射出され、製品１２を通って移動し、領域９５以外において不変のマイクロ波振動数で、減衰および位相遅れを示す。ここで、超音波は密度変位を引き起こす。この領域では、マイクロ波信号の一部分は周波数シフトされ、上述したように、上側波帯および下側波帯が生じる。送信されたマイクロ波信号９０は受信アンテナ４３を用いて集められる。超音波９１は、誘電関数の空間分布を決定する次の段階中に用いられる超音波メトリックを得る処理の間に受信機７３に集められる。
【００３７】
図３は、図１に関連して記述されるようなシステムを用いて、発明による測定原理を記述するフローチャートを示す。
【００３８】
基本的に、本発明の方法は、物質の誘電性および音響電気特性のマイクロ波超音波組合せ測定法である。ここで、解像度は超音波波長から受け継ぐ。
【００３９】
測定手順は後述されるように３つのフェーズから成る。
【００４０】
（フェーズ１）
（超音波メトリックの取得）
このフェーズでは、今後、超音波メトリックと称する、局所超音波のランタイムおよび減衰の特性のマップが確立される。
【００４１】
位相プログラミングロジックを用いて、複数の超音波振動子７２間の複数の位相を変えることによって、超音波場の任意の所望の位相形式が生成されることが可能である。超音波の半波長のオーダーの幾何学的な大きさを備えた点に超音波パワーを集中させて、すべての超音波振動子の複数の位相を制御することは可能である。できるだけ小さな体積の媒体中に超音波を集中させると、送信されたマイクロ波信号の周波数変位は最大に達する。したがって、複数の超音波振動子の位相はマイクロ波信号を最適化するために変えられる。超音波パルスと到達した最大周波数シフトとの間の遅延時間を評価することは、焦点がアンテナからのどの距離で検査下の製品２の内部に位置するかを決めることを可能にする。この測定は、所定の解像度で検査中の製品全体を覆う１セットの点で繰り返される。
【００４２】
その結果、独立した焦点ごとに選択される複数の位相とアンテナに関する位置とを含むテーブルが得られる。同時に、最大の信号の強度が、局所超音波減衰をマップすることを可能にする測定物体の至る所からのこれらの複数の測定点の各点から得られる。
【００４３】
超音波信号の局所強度は、すべての超音波振動子間のランタイムおよび減衰値の測定により計算される。（もちろん、位相のいかなる選択も、この層での各点のマイクロ波信号を最大化することにより最適化される。）複数のトランスデューサーに近い製品の層に対するこれらの遅延時間および減衰値を仮定することで、最も接近している焦点の位相が得られる。
【００４４】
１つの焦点に超音波パワーを集中させるために、伝送の複数の位相を同調させ、別の焦点上に集中させるために受信の複数の位相を同調させることによって、第１の層の２つの焦点の間のランタイムが得られる。
【００４５】
複数の焦点の周囲かつ点の次の層に接近しているところで有効なこれらの値を仮定することで、次の層の交互の点の位相および振幅値が得られる。（もちろん、位相のいかなる選択も、任意の層の各点のマイクロ波信号を最大限にすることにより最適化される）。
【００４６】
検査中の製品全体が走査されるまで、この処理が繰り返される。
【００４７】
結果は、すべての焦点にとってマイクロ波信号の強度に加えて、すべての走査された焦点間での、超音波信号の局所減衰と、超音波信号の局所位相遅れとのテーブル（「超音波メトリック」）である。
【００４８】
超音波メトリックは参照物体上で得られてもよい。それは、分析されることになっている複数の物体の代表である。その後、測定は、物体ごとの超音波メトリックを得る必要なしに、そのような物体上でなされてもよい。
【００４９】
メトリック自体も、発明の本質的な結果と見なすことが可能であり、自主的な複数の応用として用いることが可能である。更に、複数の参照物体上で得られる複数のメトリックは、フェーズ１による複数の測定を促進する手段として用いられてもよい。
【００５０】
（フェーズ２）
（マイクロ波相互作用の評価）
上記の生成された超音波メトリックとマイクロ波応答に基づいて、音響電気相互作用が、マイクロ波アンテナに最も近い層から開始して１層また１層と得られる。１層また１層の方法でこの解析が進むことは要求されないが、そのようにすることは、後の３Ｄ画像処理に便利であることがわかる。
【００５１】
各焦点で測定されたマイクロ波信号の強度は、以下の積によって決定される。すなわち、
（ａ） 超音波信号の局所強度と、
（ｂ） 圧縮率と、
（ｃ） 焦点での材料の誘電関数
との積である。
【００５２】
すべての焦点での超音波信号の局所強度が、メトリックから分かるので、マイクロ波アンテナに最も近い層上の入射および周波数シフトされた送信マイクロ波信号の間の相互作用は、この焦点での誘電関数にもたらすグリーンの関数定理を適用することによって得られる。この特定の焦点の相互作用以外の点相互作用は可能ではない。なぜならば、マイクロ波側波帯応答が、超音波焦点が測定の間に延長されている領域から始まらなければならないからである。したがって、方法の解像度は、阻害しない方法で（数センチメートルのオーダーの）マイクロ波波長によってではなく、（２５０マイクロメーターまでの）超音波信号の波束解像度から与えられる。しかしながら、入射マイクロ波信号は、送信アンテナから焦点までの途中、およびまた受信アンテナへ行く途中で、近隣の要素によって影響を受ける。焦点のマイクロ波信号は、検査中の製品のすべての誘電点に依存し、差異および入射場振幅で一次形式によって表わされる。受信アンテナに集められた場も、全ての未知の差異を含んでいる一次形式によって記述される。各測定について、すべての未知の差異を含んでいる双一次形式は得られる。各測定について、新しい方程式が生成される。各焦点の方程式があるので、方程式システムは、反復のない１対１の方法で解決することが可能である。
【００５３】
結果は、超音波メトリックと同様な基礎的な特別な構造で検査中の製品の音響電気特性と誘電特性のマップである。
【００５４】
（フェーズ３）
（音響電気誘電特性の計算）
超音波減衰はそれほど温度に依存はしない。対照的に、製品の圧縮率とともに超音波ランタイムおよび誘電関数は、強い温度依存を示す。
【００５５】
圧縮率と誘電関数との比率は、温度の関数を与える。誘電および音響電気マップを用いて、測定物体の温度は得られる。
【００５６】
第３フェーズのさらなる詳細は、図６および７ａ−７ｄに関連して記述される。
【００５７】
３つのフェーズを詳細に記述して、その測定は、図３を参照して記述される。
【００５８】
フローはステップ１００からスタートする。それは、第１の周波数ωｔｒａｎｓｍｉｔ＝２πｆ１でのマイクロ波信号が送信アンテナ４２から発信され、ωｔｒａｎｓｍｉｔとωｒｅｃｅｉｖｅの混合周波数でのマイクロ波信号が受信アンテナ４３で受信される。減衰Ｓ２１および周波数オフセットδで、２つの信号間のオフセット周波数Ｓ’２１での信号生成は、ステップ１０１で測定され、次のステップ１０２で測定された減衰Ｓ２１は以前に記録された参照減衰Ｓ２１，０に対して比較される。それは空の測定ギャップ１３（すなわち、検査中の物体がない場合）で測定された減衰に対応する。測定された減衰が、ギャップに検査中の物体がないときの減衰に等しい場合には、フローは点１０３に戻り、減衰はステップ１０１に戻り減衰が測定される。
【００５９】
物体が測定ギャップ１３に導入される場合、フローはステップ１０４に続く。ここでは、超音波メトリックが得られる。このステップは図４に関連してより詳細に記述される。
【００６０】
物体の空間誘電特性は、その後、ステップ１０４で得られたメトリックを使用して測定され計算される。この手順は、図５に関連してより詳細に記述される。
【００６１】
物体の誘電特性が決定された場合は、例えば、温度、含水量、密度などの他の物理的特性が、（所定のε（Ｔ）モデルに基づいた）誘電特性の空間分布を用いて、ステップ１０６で、決定されてもよい。そのようなモデルは現在の出願人に譲渡された、公表されたＰＣＴ出願ＷＯ０２／１８９２０に記述されるように、先行技術の中で知られている。
【００６２】
図４は、超音波メトリックを得る処理を開示したフローチャートを示す。フローはステップ１２０からスタートする。そこでは、超音波放射は物体中の点に集中する。超音波は、側波帯パスで信号を生成する。それは、マイクロ波信号によって測定された周波数変位に対応し、δと記載される。そして、ステップ１２１で測られる音響電気効率信号、ステップ１２２ではチェックが、音響電気効率信号が最大であるかどうかを決定するためになされる。もし最大でない場合にはフローはステップ１２３を介して戻り、超音波信号の位相の値が更新され、ステップ１２０に戻る。最大周波数変位が得られるまで、その処理が繰り返される。フローがステップ１２４に続いた場合、上に記述されるような焦点の位置に関する情報と共に超音波信号の位相は、メモリに格納される。ステップ１２５では、検査中の製品１２の超音波メトリックを得るために測定されるべき別の点があるかどうかが決定される。そうでなければ、ステップ１２７でメトリック得るための処理が終了する。または、フローは線分１２６を介してステップ１２０へ戻る。
【００６３】
（既知の超音波メトリックに基づいた誘電関数の測定（図４を参照））
図５ａは、物体の準備の間に、物理的に物体を精査せず、内部温度のような、物体中の物理的特性を決定するために、食品のような物体の中で誘電関数を決定するための第１の実施形態を示す。
【００６４】
フローはステップ１１０で開始する。ここで物体での点が選択される。超音波メトリックを得る処理の間に用いられた点を選択することは有利である。選択された点は方程式１−１７での点３に対応する。
【００６５】
超音波放射は、その後ステップ１１１およびステップ１１２でこの点に集中される。ＳパラメータＳ３１およびＳ２３は、図６に関連してより詳細に記述されるように、測定される。
【００６６】
ステップ１１３で、別の点が選択されるべきかどうかの決定が下される。別の点が選択されるべきである場合、フローは、ステップ１１０に戻る。ここで、新しい点はステップ１１１１および１１２の前で選択され繰り返される。選択されない場合には、フローは、ステップ１１４に続き、測定されたＳパラメータを備えた行列が逆行列化され、仮想受信機のＳ３１または仮想送信機のＳ３２を解く。
【００６７】
各選択された点ｘの誘電関数ε（ｘ）は、その後、先行技術のアルゴリズムを用いて、ステップ１１５で計算される。選択された点での温度は、その後、図３のステップ１０６によって示されるように計算される。
【００６８】
図５ｂは、物理的に物体を精査せずに、物質特性（例えば脳腫瘍の存在）のような、物体中の２つの位置間の物理的特性を決定するために、食品のような物体中の誘電関数を決定するための第２の実施形態を示す。
【００６９】
フローはステップ２１０で開始する。ここで、物体中の１組の点が選択されている。超音波メトリックを得る処理の間に用いられた点を選択することは有利である。選択された点は方程式１−１７に３および４に対応する。
【００７０】
超音波放射は、その後、ステップ２１１で両方の点に集中する。ステップ２１２で、Ｓ３１、Ｓ２３、Ｓ４１、Ｓ２４、Ｓ４’１、Ｓ２４’、Ｓ３’１およびＳ２３’のＳパラメータが、図７に関連して詳細に記述されるように、測定される。
【００７１】
ＳパラメータＳ４３（すなわち選択された点の間の減衰）は、ステップ２１３で計算される。点３は、この実施形態で、仮想送信機として作動し、点４は仮想受信機として機能する。
【００７２】
選択された点ｘ、ｙ（すなわち、方程式１−７の点３，４）の間での誘電関数の平均値
【数１】
【００７３】
は、その後、ステップ２１４で計算される。
【００７４】
ステップ２１５で、別の組の点が選択されるべきかどうかの決定が下される。別の組の点が選択されるべきである場合には、フローは、ステップ２１０に戻る。ここで、ステップ２１１から２１４までが繰り返される前に新しい組が選択される。選択されるべきでない場合には、フローは図３のステップ１０６に続く。ここで所望の物理的特性が計算される。
【００７５】
（発明の第１の使用）
図６は、本発明の第１の使用の機能を概略的に示す。１つの超音波メトリックｕ（ｘ、ｔ）が製品内のすべての点ｘに対して得られる場合、以下のステップの適用によりすべての点で誘電率を計算することが可能である。
【００７６】
１） 点３のうちの１つに超音波を集中させる。超音波は、送信アンテナ１から受信アンテナ２に送られたマイクロ波信号の振動数シフトに関する焦点にだけ影響することが知られており、これにより、側波帯（すなわちマイクロ波ベース周波数（ｆ１）±超音波周波数（ｆＵＳ））で信号を生成する。
【００７７】
２） 複数のサイドバンドの少なくとも１つで信号強度を測定する。両側波帯の中で信号強度が測定される場合、測定からのより信頼できる結果が得られる。受信アンテナ２で測定された信号強度は、次のように表現されてもよい：
【数２】
【００７８】
ここで、Ｓ２１は、測定ギャップにある製品１２によって引き起こされた減衰であり、Ｖ２（ｔ）は側波帯で測定された信号強度であり、Ｖ１（ｔ）は送信アンテナ１から送られた信号の信号強度である。Ｓ２３は受信アンテナ２への点３の間の減衰であり、α３は、超音波がマイクロ波側波帯信号（音響電気利得と称する）に変換される点３での効率を決定する係数である。ｕ３（ｘ、ｔ）は点３での超音波メトリックであり、Ｓ３１は送信アンテナ１および点３の間の減衰である。
【００７９】
第１の近似では、効率αは次のように表現することが可能である。
【数３】
【００８０】
ここで、Δεは、超音波放射ｙによって引き起こされる圧力波による誘電率の変化である。圧縮モジュールκに対して、
【数４】
【００８１】
の関係が確立する。Κの値は、当業者に知られていて、より詳細に検討しない。
【００８２】
３） 製品１２中の図６に３を表示する全ての所望の点に対して処理を繰り返す。
【００８３】
４） 逆散乱アルゴリズムでのすべての測定データを使用し、製品中の誘電関数の空間分布を計算する。
【００８４】
物体が、測定装置に関して、相対的に低速で移動し、かつ、下の関係を満たす場合、射出された超音波とマイクロ波放射の補償は考慮する必要がない。
【数５】
【００８５】
ｖｏｂｊは、測定ギャップ１３での物体移動の速度であり、ｔｍｅａｓは完全な処理ための測定時間であり、ｖＵＳは物体中の超音波の速度であり、ｆＵＳは超音波周波数であり、ｄＦｏｃａｌは焦点の直径である。
【００８６】
相対速度が高い場合、超音波の集束は移動を補償するために、測定ステップの間、物体の中に焦点を維持するために、超音波放射の調整を含んでいなければならない。さらに、ドップラー偏移を回避するために、
【数６】
【００８７】
である。
【００８８】
（発明の第２の使用）
図７ａ−７ｄは、製品中の点３および４の２点間の誘電率を計算する際の、本発明の第２の使用の主関数を示す。第１の点３はソースであると見なしてもよく、第２の点４は受信機であると見なしてもよい。
【００８９】
主関数は、図６に関連して記述したものと非常に似ているが、２つの焦点３および４が超音波放射によって同時に生成されたので、２つの上側波帯および、２つの下側波帯が生成されるという違いがある。第１の上側波帯および下側波帯は、図６に関連して記述されたものと同様であり、第２の上側波帯および下側波帯は、２倍の超音波周波数（すなわち、マイクロ波ベース周波数（ｆ１）±２＊超音波周波数（２ｆＵＳ））を持っている。同じ超音波周波数がこの目的に用いられる場合、第２のオーダーの側波帯を生成するために、２つの異なる超音波周波数を選ぶことが可能である。図１に関連して記述された装置は、この例で、第２の上側波帯および下側波帯に応じて調整される余分な側波帯パスが付加される必要がある。
【００９０】
次の関係は、それぞれ単一の仮想ソースとしての点３および４にとって確立されることが可能である。
【数７】
【００９１】
図７ｂにしたがって、焦点を３から３’へ移動させ、焦点を４から４’へ移動することによって、新しい関係が表現することができる。
【数８】
【００９２】
図７ａから、点３および４の間の求められた減衰を含む関係は次のように表現されてもよい。
【数９】
【００９３】
方程式６は、７ｘ７問題を解くのに用いられず、適切な近似と取り替えられる。方程式１６および１７を参照。
【００９４】
図７ｃは、３および４に対応する２重のソースの関係を示す。
【数１０】
【００９５】
点３’および４’の関係が次のように表現されてもよい。
【数１１】
【００９６】
方程式１０は７ｘ７問題と８ｘ８問題を解くのに用いられず、適切な近似と取り替えられる。８ｘ８問題のためには方程式１５を参照し、７ｘ７問題のためには方程式１６および１７を参照する。
【００９７】
次の関係は図７ａ−７ｃから明白である。
【数１２】
【００９８】
方程式１１−１４は複数のＳパラメータを除去するために用いられる。それは図７ｄに示されるようなＳパラメータになる。求められたＳ４３であるＳパラメータがあり、問題を解くための１０の方程式（すなわち、方程式１−１０）を必要とするいくつかの未知のＳパラメータと共に、全く興味のないＳ３’４’である１つのＳパラメータがある。
【００９９】
有限要素についてのツィエンキービッツ（Ｚｉｅｎｋｉｅｗｉｃｚ）によって導入されたトリックを適用することによって、点３と点４との間の減衰を見つけるのに必要とされる方程式の数を減らすことは可能である。
【０１００】
方程式１０は用いられず、次の近似が代わりに用いられる。
【数１３】
【０１０１】
２度ツィエンキービッツのトリック（それは方程式６および１０の必要を除去する）を適用することによって、必要とされる方程式の数を８つの方程式のみに減らすことは可能でさえある。複数の方程式の代わりに用いられる近似は次のとおりである。
【数１４】
【０１０２】
点３と点４との間および点３’と点４’との間の減衰Ｓ４３は、必要とされる方程式を対数にすることにより計算することが可能であり、複数の解析点が適切に選ばれる限り解は常に利用可能である方程式と同数の未知を備えた方程式１から１０が方程式の非同次線形システムになる。これらの点の役割を「仮想プローブ」として示す点４と点３の間のマイクロ波ランタイムを得るためにＳ４３のシステムを解かねばならない。
【０１０３】
上に記述されたシステムは、「仮想送信機」（すなわち、点３）および「仮想受信機」（すなわち、点４）を使用する。発明の第１の使用法に到達するそれぞれ本物の送信または受信アンテナと一致するために、これらの点の１つを容易に置くことが可能である。複数の物理的なプローブアンテナの場所に両仮想プローブを置くことは、本発明に先立って既知の従来からのマイクロ波測定技術になる。
【０１０４】
解決される物理的問題に依存して、単一仮想プローブ概念（仮想受信機または仮想送信機）または２つの仮想プローブ概念を利用する。複数の仮想プローブによって生成／受信された特定のビームパターンを生じるために、複数のプローブのセット（例えば、仮想プローブのアレイ）を用いることは可能である。
【０１０５】
異なるプローブ構成は、掃海、物質解析、鉱物探査、医学的応用などとして応用に用いられてもよい。
【０１０６】
（方法の短縮した数学的導出）
電磁放射は、マクスウェル方程式によって統率され、ベクトルの電場Ｅは、次に示すように座標に依存する３次元の空間ｘおよび時間ｔで書かれているヘルムホルツ形式に容易に変形される。
【数１５】
【０１０７】
ここで、Δはラプラス演算子であり、ε０は真空の誘電率であり、εｒ（３ｘ３テンソルである）は与えられた位置での材料の局所的で相対的な誘電関数であり、μ０は真空透磁率を表し、μｒは検査中の材料の局所比透磁率を表す。この短縮した導出では、μｒは単位テンソル１（３ｘ３）に設定する。当業者にとっては、類似の方法が同時にεｒとμｒを解くことによって導出されることは明白である。
【０１０８】
同時に、テンソルの３ｘ３応力振幅ｙおよび媒体の局所音速ｖを備えた超音波も、次に示す同様の形式にすることが可能である。
【数１６】
【０１０９】
両方の微分方程式の解は、放射ソースの位置を考慮して実行される。処理の要点に集中させて、ゼロにならない振幅を備えたいかなる超音波も（圧縮または剪断型である）材料に応力を生じる。この応力は材料の局所圧縮によって反映される。この圧縮によって、分極電荷の密度は影響を受ける。既知の事実として、誘電性の物体のいかなる圧縮も相対的な誘電関数テンソルεｒを次のように変更する。
【数１７】
【０１１０】
この関係は、本発明で利用された超音波伝播と電磁波との間のカップリングを生じる。相互作用の強度は、音響電気光学の相互作用、３ｘ３ｘ３テンソルであるαによって決定される。含まれる物理学の全体像のために、上記の関係は（例えば、キャビテーションおよび他の非線形効果が無視される）比較的小さな超音波に対してだけ保たれることを言及しなければならない。
【０１１１】
その後、電磁気的に解決される完全なシステムは、次によって与えられる。
【数１８】
【０１１２】
当業者にとって、時間ｔでフーリエ変換が適用される場合、この種の微分方程式が周波数空間ωでの畳み込みになることは明白である。
【数１９】
【０１１３】
ここで、円で囲まれた積算演算子Ｅ（ｘ，ω）は、（畳み込み積分の前の、結果として生じる規格化定数を除いて）次の完全形になる（例えば、１９８８年にG.Doetschによる“Anleitung zum praktischen gebrauch der Laplace transformation”に見出される）周波数畳み込み積分を示す。
【数２０】
【０１１４】
したがって、単一の周波数超音波励起および物体に入射する単一の周波数マイクロ波信号を仮定すると、複数の受信マイクロ波信号は、入射するマイクロ波周波数に一部分を含むだけでなく、畳み込み積分によって生じた超音波およびマイクロ波周波数の差および和での複数の側波帯にも含む。
【０１１５】
上記の関係は、マイクロ波場から情報を抽出するために、超音波を適切に位相制御し、パルス波列を用いることによって、全く新しい世界を提示する。
【０１１６】
（単一の仮想プローブ）
単一の仮想プローブを含むパスに沿って解決するべき方法を適用する。これは、（点３または４のいずれかへの全ての関係が消える）上述した到来する線形方程式系を解くどんな送信パラメータにも依存する、仮想送信機または仮想受信機のどちらかに相当する。波の伝播メカニズムはこのケースには同一である。（一様な、境界条件無しの）理想的なケースについては、次の伝播関係に達する。
【数２１】
【０１１７】
（２つの仮想プローブ）
さらに、２つの仮想プローブを通るパスに沿って解決するために方法を適用することが可能である。これは、（すべて方程式が存在する）上述された到来する９ｘ９線形方程式系を解決するどんな伝送パラメータにも依存する、仮想送信機または仮想受信機のいずれかに相当する。（一様な、境界条件無しの）理想的なケースについては、次の伝播関係に達する。
【数２２】
【０１１８】
最初の２つの方程式は、仮想送信機の役割を務める解析点Ｘで側波帯の生成を示す。第３の方程式は、仮想受信機の役割を務める別の解析点Ｙで集中させることによって、最初の式の上部の第２の側波帯の生成を示す。周波数オフセットは、集中させることを遂行するために使用した超音波の周波数によって決定された、点Ｘでのηおよび点Ｙでのηで示される。これらが、ある応用では、両方の点Ｘ、Ｙにとって同じ周波数ではないことがあることに注意する。
【０１１９】
第１方程式は、側波帯オフセットＸを備えた所定の位置ξで側波帯の生成を述べる。第２の方程式は、強度ｑのソースが位置Ｘに置かれる場合、検査中の物体全体を通る側波帯の伝播を述べる。したがって、方法は、物体の内部の任意の複数の位置でのマイクロ波ソースを合成することによって、物体を「調査する」ことを可能にする。そして、このソースを動かす場合、このソースから生成された放射を測定する。
【０１２０】
発明は、マイクロ波発生器および超音波発生器に関連して記述されたが、物体内の密度変位を生じるために、他の型の放射が用いられてもよいことは明白である。しかしながら、放射は同時に射出されなければならず、さらに、変位を生じるために射出された放射間の周波数に差がなければならない。解像度は、物体で最短の波長を有する放射によって決定される。
【０１２１】
したがって、密度変位を生じ、かつ、それによって、発明を用いて材料の誘電関数を決定するために、異なる（例えばわずかに０．５Ｈｚだけ異なる）周波数を持っている２つのマイクロ波信号で物体を同時に照らすことは可能である。射出された放射の可能な組合せは、これに制限するわけではないが、マイクロ波、超音波およびＸ線の任意の組合せ含んでいる。
【図面の簡単な説明】
【０１２２】
【図１】発明のシステムを示す。
【図２】検査中の製品へ射出される放射線を示す。
【図３】検査中の製品の内部で、温度のような物理的特性を決定するためにフローチャートを示す。
【図４】超音波メトリックを得る処理を示すフローチャートを示す。
【図５ａ】検査中の製品内の誘電関数の空間分布を決定する処理の２つの実施形態を示すフローチャートを示す。
【図５ｂ】検査中の製品内の誘電関数の空間分布を決定する処理の２つの実施形態を示すフローチャートを示す。
【図６】本発明の第１の使用の主関数を示す。
【図７ａ】本発明の第２の使用の主関数を示す。
【図７ｂ】本発明の第２の使用の主関数を示す。
【図７ｃ】本発明の第２の使用の主関数を示す。
【図７ｄ】本発明の第２の使用の主関数を示す。
【技術分野】
【０００１】
本発明は、請求項１および請求項１９のプリアンブルに記載の物体の誘電関数を決定することによって、物質の内部の、温度または密度のような複数の物理的なパラメータを決定するための装置に関する。発明はまた、請求項１２のプリアンブルに記載の物体の内部の誘電関数を決定するための方法、および、請求項１８に記載の食品の中での温度の局所分布を決定するための装置に関する。
【背景技術】
【０００２】
物体を破壊することなく、物体に侵入せず、物体を切り裂くことなく、任意の物体の、温度、密度および他の複数の内部パラメータに関する情報を得るために、様々な型の放射は、所望の複数のパラメータを再構成することを可能にする情報を提供するのに利用可能である。
【０００３】
放射の特殊型式を選んで、分析法の選択において適切な予想結果を組込む４つの異なる場合がある。これらは、２つの問題領域によって分類される。
【０００４】
−選ばれた放射に対する物体の透明度
−選ばれた放射の波長に関して必要とされる物体での解像度
ケース１Ａ
（物体は、測定に用いられた放射に対して、透明であるか弱く吸収する。また、達成される解像度は、放射波長に等しいか、放射波長よりも小さい。）
ただ一つの情報源は、例えば、以下に述べるような手法を使用して近い場をプローブすることにより得られる。
【０００５】
−原子間力顕微鏡検査（ＡＦＭ）は、
検査中の物体の表面上の構造を読み出すために、物体の表面上に高精度で位置決めされるサブ波長サイズのステンシル上の力を読み出す。
【０００６】
−ラスタートンネル顕微鏡検査（ＲＴＭ）は、
物体の表面の電子状態の情報を生成するために、力の代わりに、検査中の物体の表面に接近して位置決めされるサブ波長サイズのプローブからトンネル電流を測定する。
【０００７】
−光学的近接場顕微鏡検査は、
電磁放射が、薄いプローブ上でのサブ波長解像度で光学的特性の表面画像を生成することを用いた放射の波長よりもはるかに小さい穴を必要とする微視的に小さな複数の穴を介して通り抜ける。
【０００８】
−インピーダンス断層撮影法は、
１セットの電極は検査中の物体に着けられ、また、すべてのプローブ間のインピーダンスが測定される。この方法は検査中の物体の内部のいくつかの特性を計算することを可能にする。しかし、解像度は一般に劣っている。この方法は、例えば抗凝固薬の影響を評価する前後で心臓の領域のインピーダンスを測定するなどの異なる複数の手法での成功で使用されてきている。
【０００９】
一般的な特徴として、上述した方法の高解像度は、選ばれた放射の固有波長に起因しているのではなく、サブ波長解像度を提供する別の拘束条件（横隔膜、ステンシルのようにほとんど機械的）に起因する。一般的な欠点は検査中の物体の厚さ要求によって与えられる−上記の複数の方法は、解像度を失わずに、非常に制限のある深さでの単に表面情報または内部情報のいずれかを生成する。
【００１０】
ケース１Ｂ
（物体は、測定に用いられる放射に対して、透明であるか、または弱く吸収することである。また、解像度は放射波長よりもはるかに大きい。）
このケースはすべての直接撮像および光伝送方法によってカバーされる。このレジーム（ｒｅｇｉｍｅ）で電磁放射を用いるとすると、以下のものがある。
【００１１】
−ライダー
−Ｘ線
分析の手段として、散乱が役割を果たさないので、レイトレーシングおよび１対１写像の方法が適切である。それは、放射源と受信機との間に位置する物体体積によって、ある与えられた位置で採用されたピクセル情報のそれぞれが影響するに過ぎないことを、解像度の損失なしで仮定することが可能である。
【００１２】
この領域の最近の進展はパッシブレーダーであり、受信機のまわりの環境ですべての身体に固有の熱放射が測定され撮像される。このレーダー方法は、どんな送信された信号をも必要とせず、したがって、跡をたどることができない。
【００１３】
非電磁気的な方法では、商業上利用可能な以下のものがある。
【００１４】
−超音波断層法
−核磁気共鳴（ＮＭＲ）
ケース２Ａ
（物体は、測定に用いられた放射に対して適度に吸収している。また、解像度は放射波長に等しいかそれよりも小さい。）
物体が測定に用いられた放射に対して適度に吸収しているという事実は、調査することが可能であるプローブに厚さ限界を設けることになる。
【００１５】
このケースについては、今日、最先端技術に関して利用可能な実現可能な方法はない。
【００１６】
ケース２Ｂ
（物体は、測定に用いられた放射に対して適度に吸収している。また、解像度は放射波長よりもはるかに大きい。）
この場合、ほとんどの無線周波数およびマイクロ波周波数の応用は見つかる（特に、検査中の物体が、損失があり、それが損失のない環境に組み込まれた場合）。また、マイクロ波断層撮影法が利用可能である。これらの方法で、最もよくあるものが次のものである。
【００１７】
−（アクティブ）レーダー（ＲＡＤＡＲ）
ソースとターゲットの間および受信機へ戻る信号のランタイムは、送信機（モノスタティックレーダー）と同じ場所に受信機を配置することによってか、もしくは、送信機（バイスタティックレーダー）と、ソースとターゲットの相対速度による周波数変化とが評価される（ドップラーレーダー）よりも異なる位置に受信機を配置することによって測定される。
【００１８】
このように、測定に用いられた放射に対して適度に吸収する物体にとって、温度、密度、組成のような物理的なパラメータを決定するための装置を開発する必要がある。所望の解像度は放射波長よりもはるかに大きい。
【発明の開示】
【課題を解決するための手段】
【００１９】
本発明の目的は、任意に形成された物体の誘電関数を決定するための装置を提供することである。
【００２０】
目的は、製品密度における制御可能な変化を作成するための超音波を用いて、請求項１および請求項１９の特徴部分で定義されるような装置、および請求項１２の特徴部分に記載の方法によって達成される。その後、装置は、密度変化を読出し、かつ誘電関数の空間分布にそれを関連づけるために、マイクロ波放射を用いる。これは、請求項１８の特徴部分に定義されるように、物体の温度、含水量および密度を決定するために一方で用いられる。
【００２１】
本発明での利点は、空間分布の解像度が、第１型放射（例えばマイクロ波放射)の波長に限定されず、むしろ第２型放射（例えば超音波、Ｘ線)の波長によって決定されることである。
【００２２】
本発明での別の利点は、温度、含水量などのような物理的特性の非接触測定が、複数の仮想プローブとして発明を適用して確立されてもよいということである。
【００２３】
他の複数の物体および複数の利点は、発明の詳細な記述からの当業者にとって明白になる。
【発明を実施するための最良の形態】
【００２４】
本発明に先立って、材料のみの内部特性を再構成するためのツールとして以下のものが存在する（ただし、回析および散乱は優勢である）。
【００２５】
−マイクロ波断層撮影法
−超音波断層法
両方のケースでは、解像度は使用された放射の波長によって決定される。
【００２６】
本発明では、超音波およびマイクロ波の方法は組み合わせられる。物体再構成は純粋なマイクロ波逆散乱方法、およびそれぞれの限定を有する純粋な超音波断層法方法によって行うことが可能である。ここで、超音波は、物体再構成ツールとして用いられず、調査される物体での密度変化を生成するツールとして用いられる。この密度変化は、物体再構成に用いられる送信されたマイクロ波放射線に、位相および周波数の変化を生じる。したがって、この方法の利用可能な解像度は、（典型的な超音波治療周波数用の１ミリメートルよりも小さな）超音波の解像度によって決定される。密度読み出しは、マイクロ波放射（ある周波数では、減衰はさらに合理的な複数の浸入度（例えばＳ、ＩＳＭ５．８、Ｘバンド）を許可する）を用いて行われる。この方法は、ミリメートル解像度がミリメートル波長を必要とするマイクロ波断層撮影法の複数のアプローチの基本的な困難を回避する。不運なことに、ミリメートル放射は、いくつかの波長内で興味深いほとんどの物体に吸収され、したがって、いかなる内部のパラメータも抽出されることは不可能である。上記の分類では、本発明は領域１Ｂ、２Ａおよび２Ｂをカバーする。そのような方法は本発明に先立って知られていない。
【００２７】
本発明によって記述されたシステムは、食品産業で用いられることが好適である。食品産業では、正確に食品の温度を制御することは多くの場合、重要である。例えば、食品が冷凍されることになっている場合、製品全体が冷凍されることが重要である。製品全体（例えば鶏フィレ）が冷凍されたことを保証することが可能でない場合、製品を廃棄しなければならないか、または短い保存寿命を備えた製品を配達しなければならない。したがって、製品の冷凍の非破壊および非接触の制御が必要である。この問題は、誘電関数を測定し、下記に記述されるようにそれを温度の分布に変換することによって解決されてもよい。
【００２８】
しかしながら、システムは決してこの種の産業に制限されているわけではない。見込みのある他の応用は以下がある。
−コンクリート硬化（建設業）
−接着剤硬化（飛行機製造）
−医療用画像（機能的な脳断層撮影法、脊髄の断層撮影法）
−地盤調査、トラッキングパイプおよび地中管
−救命および救助道具（瓦礫下の人を検知）
−掃海（大きくなる領域での特にプラスチックの機雷）
以下では、好ましい実施形態が要約される。上記の他の応用領域でこの方法を適応させするために形状に必要とされる修正は小さい。
【００２９】
以下では、持続波（ＣＷ）マイクロ波およびパルス波列超音波に基づいたシステムは、単純性の目的のために記述される。記述された方法はこのケースに制限されない。振幅変調（ＡＭ）、周波数変調（ＦＭ）、周波数変調持続波（ＦＭＣＷ）、パルス符号変調（ＰＣＭ）、位相変調（ＰＭ）およびウェーブレットに基づいた変調方式（ＷＭ）のような、両者のための他の変調スキーム、電磁波、および、超音波は適用可能であり、ある他の応用には最適である。
【００３０】
図１は発明による装置４０を記述する。システムはコンベア手段１１に接近して配置される。それは、センサ測定ギャップ１３を介して検査下の製品１２を運搬する。システム４０は、マイクロ波部５０、超音波部７０および評価ユニット６０から成る。システムは、本実施形態では、２固定周波数マイクロ波発生器５１および５２、ならびに固定周波数超音波発生器７１を具備する。第１マイクロ波発生器５１は第１の固定マイクロ波振動数ｆ１（例えば５．８１８ＧＨｚ）を持っており、少なくとも１つの送信アンテナ４２につながれる。また、第２マイクロ波発生器５２は第２の固定マイクロ波振動数ｆ２（例えば５．８ＧＨｚ）を持っており、好ましくはミキサのようなダウンコンバータ５４につながれる。ダウンコンバータは、送信される（少なくとも１つの受信アンテナ４３によって集められる）マイクロ波信号、および、第２マイクロ波発生器５２からの受信マイクロ波信号を、低い中間周波数にシフトする。これは、検査下の製品１２を通って送信されるマイクロ波信号が振幅および位相で評価されることを可能にする。それは更に、フィルタユニット５９、アナログ／デジタル変換器ＡＤＣ５５、１セットの信号処理部５６、および、システムを制御し、かつデータを評価するために必要なアルゴリズムを含んでいる評価プロセッサ６０を含む。結果は表示ユニット６５に提出される。システム４０はまた、送信アンテナ４２および受信アンテナ４３に加えて、測定ギャップ１３のまわりにすべてグループ化して、（明瞭さの目的のために１つだけ示されている）１セットのトランスデューサー７２を含む。トランスデューサーは、検査下の製品１２を通る、超音波周波数ｆＵＳ（例えば４．５ＭＨｚ）を持っている超音波信号を発する。これは、超音波速度で移動する密度変位を引き起こす。同時に、第１マイクロ波発生器５１からのマイクロ波信号が、送信アンテナ４２から射出される。この信号はまた、検査下の製品１２を通って移動する。マイクロ波信号は、マイクロ波周波数を保ったまま製品を通って移動することにより減衰および位相遅れを示す。検査下の製品１２の超音波が密度変位を生じるそれらの体積では、マイクロ波信号の一部分が周波数シフトされ、そして、上側波帯と下側波帯が生じる。透過したマイクロ波信号はマイクロ波受信アンテナ４３を用いて集められる。受信信号はダウンコンバータユニット５４を用いてダウンコンバートされる。その後、低周波信号は、フィルタユニット５９を用いてフィルタリングされ、ＡＤＣ５５を用いてアナログデジタル変換される。デジタル信号は受信信号処理部５６を用いて評価される。受信信号処理部は、標準の最新技術のデジタルフィルターを用いて、入ってくるデジタル信号をゼロ周波数に変換する。
【００３１】
このフィルタリングの結果は、当業者によく知られているように送信４２および受信４３アンテナとの間のＳ２１パラメータに相当する。これは周波数がシフトされていない。上記では、我々は、受信アンテナ４３をマイクロ波ポート２と称し、送信アンテナ４２をマイクロ波ポート１と称する。
【００３２】
本発明によって記述されるシステムでは、第１信号パス５９、５５、５６に平行して、第２セットのバンドパスフィルター５８、別のＡＤＣ５５および第２のデジタル信号処理部５７がある。
【００３３】
バンドパスフィルター５９は、両方のマイクロ波発生器５１および５２の間の差周波数に合わせられる。それは本実施形態では５．８１８ＧＨｚ−５．８ＧＨｚ＝１８ＭＨｚである。第２バンドパスフィルター５７は、マイクロ波発生器間の差周波数（例えば１８ＭＨｚ）に超音波信号発生器７１の中心周波数（例えば４．５ＭＨｚ）を加えた周波数に合わせられる。したがって、５８、５５および５７を収納しているこの第２のデジタル信号処理部のパスは、入ってくる信号を、超音波周波数によって周波数シフトされたゼロ周波数に変換する。したがって、測定結果は、超音波とマイクロ波信号との間の断面に制限されている。
【００３４】
第１の５９，５５，５６と第２の５８，５５，５７のデジタル受信機のＩＦ帯域幅は、超音波発生器７１によって生成された、超音波周波数ｆＵＳの半分に選ばれる。これは、超音波振動子の複数の位相を変えることによって、周波数シフトを最適化するために要求される。
【００３５】
製品１２の超音波メトリックを得る第１段階の間に、トランスデューサー７２から射出される超音波放射を収集する超音波受信機７３が存在して、より詳細に下に記述されるような減衰（Ｔ５６）およびランタイムを評価しなければならない。上記で、我々は、超音波受信機７３をマイクロ波ポート６と称し、トランスデューサー７２をマイクロ波ポート５と称する。減衰およびランタイムが、超音波評価ユニット７４で評価されるが、これは、評価ユニット６０に統合されても当然よい。
【００３６】
図２は、検査中の製品へ射出された放射を示す。トランスデューサー７２は、この例において、検査下の製品１２を通る超音波パルス９１を射出する。これは、超音波速度で移動する密度変位を引き起こす。同時に、マイクロ波信号９０は、送信アンテナ４２から射出され、製品１２を通って移動し、領域９５以外において不変のマイクロ波振動数で、減衰および位相遅れを示す。ここで、超音波は密度変位を引き起こす。この領域では、マイクロ波信号の一部分は周波数シフトされ、上述したように、上側波帯および下側波帯が生じる。送信されたマイクロ波信号９０は受信アンテナ４３を用いて集められる。超音波９１は、誘電関数の空間分布を決定する次の段階中に用いられる超音波メトリックを得る処理の間に受信機７３に集められる。
【００３７】
図３は、図１に関連して記述されるようなシステムを用いて、発明による測定原理を記述するフローチャートを示す。
【００３８】
基本的に、本発明の方法は、物質の誘電性および音響電気特性のマイクロ波超音波組合せ測定法である。ここで、解像度は超音波波長から受け継ぐ。
【００３９】
測定手順は後述されるように３つのフェーズから成る。
【００４０】
（フェーズ１）
（超音波メトリックの取得）
このフェーズでは、今後、超音波メトリックと称する、局所超音波のランタイムおよび減衰の特性のマップが確立される。
【００４１】
位相プログラミングロジックを用いて、複数の超音波振動子７２間の複数の位相を変えることによって、超音波場の任意の所望の位相形式が生成されることが可能である。超音波の半波長のオーダーの幾何学的な大きさを備えた点に超音波パワーを集中させて、すべての超音波振動子の複数の位相を制御することは可能である。できるだけ小さな体積の媒体中に超音波を集中させると、送信されたマイクロ波信号の周波数変位は最大に達する。したがって、複数の超音波振動子の位相はマイクロ波信号を最適化するために変えられる。超音波パルスと到達した最大周波数シフトとの間の遅延時間を評価することは、焦点がアンテナからのどの距離で検査下の製品２の内部に位置するかを決めることを可能にする。この測定は、所定の解像度で検査中の製品全体を覆う１セットの点で繰り返される。
【００４２】
その結果、独立した焦点ごとに選択される複数の位相とアンテナに関する位置とを含むテーブルが得られる。同時に、最大の信号の強度が、局所超音波減衰をマップすることを可能にする測定物体の至る所からのこれらの複数の測定点の各点から得られる。
【００４３】
超音波信号の局所強度は、すべての超音波振動子間のランタイムおよび減衰値の測定により計算される。（もちろん、位相のいかなる選択も、この層での各点のマイクロ波信号を最大化することにより最適化される。）複数のトランスデューサーに近い製品の層に対するこれらの遅延時間および減衰値を仮定することで、最も接近している焦点の位相が得られる。
【００４４】
１つの焦点に超音波パワーを集中させるために、伝送の複数の位相を同調させ、別の焦点上に集中させるために受信の複数の位相を同調させることによって、第１の層の２つの焦点の間のランタイムが得られる。
【００４５】
複数の焦点の周囲かつ点の次の層に接近しているところで有効なこれらの値を仮定することで、次の層の交互の点の位相および振幅値が得られる。（もちろん、位相のいかなる選択も、任意の層の各点のマイクロ波信号を最大限にすることにより最適化される）。
【００４６】
検査中の製品全体が走査されるまで、この処理が繰り返される。
【００４７】
結果は、すべての焦点にとってマイクロ波信号の強度に加えて、すべての走査された焦点間での、超音波信号の局所減衰と、超音波信号の局所位相遅れとのテーブル（「超音波メトリック」）である。
【００４８】
超音波メトリックは参照物体上で得られてもよい。それは、分析されることになっている複数の物体の代表である。その後、測定は、物体ごとの超音波メトリックを得る必要なしに、そのような物体上でなされてもよい。
【００４９】
メトリック自体も、発明の本質的な結果と見なすことが可能であり、自主的な複数の応用として用いることが可能である。更に、複数の参照物体上で得られる複数のメトリックは、フェーズ１による複数の測定を促進する手段として用いられてもよい。
【００５０】
（フェーズ２）
（マイクロ波相互作用の評価）
上記の生成された超音波メトリックとマイクロ波応答に基づいて、音響電気相互作用が、マイクロ波アンテナに最も近い層から開始して１層また１層と得られる。１層また１層の方法でこの解析が進むことは要求されないが、そのようにすることは、後の３Ｄ画像処理に便利であることがわかる。
【００５１】
各焦点で測定されたマイクロ波信号の強度は、以下の積によって決定される。すなわち、
（ａ） 超音波信号の局所強度と、
（ｂ） 圧縮率と、
（ｃ） 焦点での材料の誘電関数
との積である。
【００５２】
すべての焦点での超音波信号の局所強度が、メトリックから分かるので、マイクロ波アンテナに最も近い層上の入射および周波数シフトされた送信マイクロ波信号の間の相互作用は、この焦点での誘電関数にもたらすグリーンの関数定理を適用することによって得られる。この特定の焦点の相互作用以外の点相互作用は可能ではない。なぜならば、マイクロ波側波帯応答が、超音波焦点が測定の間に延長されている領域から始まらなければならないからである。したがって、方法の解像度は、阻害しない方法で（数センチメートルのオーダーの）マイクロ波波長によってではなく、（２５０マイクロメーターまでの）超音波信号の波束解像度から与えられる。しかしながら、入射マイクロ波信号は、送信アンテナから焦点までの途中、およびまた受信アンテナへ行く途中で、近隣の要素によって影響を受ける。焦点のマイクロ波信号は、検査中の製品のすべての誘電点に依存し、差異および入射場振幅で一次形式によって表わされる。受信アンテナに集められた場も、全ての未知の差異を含んでいる一次形式によって記述される。各測定について、すべての未知の差異を含んでいる双一次形式は得られる。各測定について、新しい方程式が生成される。各焦点の方程式があるので、方程式システムは、反復のない１対１の方法で解決することが可能である。
【００５３】
結果は、超音波メトリックと同様な基礎的な特別な構造で検査中の製品の音響電気特性と誘電特性のマップである。
【００５４】
（フェーズ３）
（音響電気誘電特性の計算）
超音波減衰はそれほど温度に依存はしない。対照的に、製品の圧縮率とともに超音波ランタイムおよび誘電関数は、強い温度依存を示す。
【００５５】
圧縮率と誘電関数との比率は、温度の関数を与える。誘電および音響電気マップを用いて、測定物体の温度は得られる。
【００５６】
第３フェーズのさらなる詳細は、図６および７ａ−７ｄに関連して記述される。
【００５７】
３つのフェーズを詳細に記述して、その測定は、図３を参照して記述される。
【００５８】
フローはステップ１００からスタートする。それは、第１の周波数ωｔｒａｎｓｍｉｔ＝２πｆ１でのマイクロ波信号が送信アンテナ４２から発信され、ωｔｒａｎｓｍｉｔとωｒｅｃｅｉｖｅの混合周波数でのマイクロ波信号が受信アンテナ４３で受信される。減衰Ｓ２１および周波数オフセットδで、２つの信号間のオフセット周波数Ｓ’２１での信号生成は、ステップ１０１で測定され、次のステップ１０２で測定された減衰Ｓ２１は以前に記録された参照減衰Ｓ２１，０に対して比較される。それは空の測定ギャップ１３（すなわち、検査中の物体がない場合）で測定された減衰に対応する。測定された減衰が、ギャップに検査中の物体がないときの減衰に等しい場合には、フローは点１０３に戻り、減衰はステップ１０１に戻り減衰が測定される。
【００５９】
物体が測定ギャップ１３に導入される場合、フローはステップ１０４に続く。ここでは、超音波メトリックが得られる。このステップは図４に関連してより詳細に記述される。
【００６０】
物体の空間誘電特性は、その後、ステップ１０４で得られたメトリックを使用して測定され計算される。この手順は、図５に関連してより詳細に記述される。
【００６１】
物体の誘電特性が決定された場合は、例えば、温度、含水量、密度などの他の物理的特性が、（所定のε（Ｔ）モデルに基づいた）誘電特性の空間分布を用いて、ステップ１０６で、決定されてもよい。そのようなモデルは現在の出願人に譲渡された、公表されたＰＣＴ出願ＷＯ０２／１８９２０に記述されるように、先行技術の中で知られている。
【００６２】
図４は、超音波メトリックを得る処理を開示したフローチャートを示す。フローはステップ１２０からスタートする。そこでは、超音波放射は物体中の点に集中する。超音波は、側波帯パスで信号を生成する。それは、マイクロ波信号によって測定された周波数変位に対応し、δと記載される。そして、ステップ１２１で測られる音響電気効率信号、ステップ１２２ではチェックが、音響電気効率信号が最大であるかどうかを決定するためになされる。もし最大でない場合にはフローはステップ１２３を介して戻り、超音波信号の位相の値が更新され、ステップ１２０に戻る。最大周波数変位が得られるまで、その処理が繰り返される。フローがステップ１２４に続いた場合、上に記述されるような焦点の位置に関する情報と共に超音波信号の位相は、メモリに格納される。ステップ１２５では、検査中の製品１２の超音波メトリックを得るために測定されるべき別の点があるかどうかが決定される。そうでなければ、ステップ１２７でメトリック得るための処理が終了する。または、フローは線分１２６を介してステップ１２０へ戻る。
【００６３】
（既知の超音波メトリックに基づいた誘電関数の測定（図４を参照））
図５ａは、物体の準備の間に、物理的に物体を精査せず、内部温度のような、物体中の物理的特性を決定するために、食品のような物体の中で誘電関数を決定するための第１の実施形態を示す。
【００６４】
フローはステップ１１０で開始する。ここで物体での点が選択される。超音波メトリックを得る処理の間に用いられた点を選択することは有利である。選択された点は方程式１−１７での点３に対応する。
【００６５】
超音波放射は、その後ステップ１１１およびステップ１１２でこの点に集中される。ＳパラメータＳ３１およびＳ２３は、図６に関連してより詳細に記述されるように、測定される。
【００６６】
ステップ１１３で、別の点が選択されるべきかどうかの決定が下される。別の点が選択されるべきである場合、フローは、ステップ１１０に戻る。ここで、新しい点はステップ１１１１および１１２の前で選択され繰り返される。選択されない場合には、フローは、ステップ１１４に続き、測定されたＳパラメータを備えた行列が逆行列化され、仮想受信機のＳ３１または仮想送信機のＳ３２を解く。
【００６７】
各選択された点ｘの誘電関数ε（ｘ）は、その後、先行技術のアルゴリズムを用いて、ステップ１１５で計算される。選択された点での温度は、その後、図３のステップ１０６によって示されるように計算される。
【００６８】
図５ｂは、物理的に物体を精査せずに、物質特性（例えば脳腫瘍の存在）のような、物体中の２つの位置間の物理的特性を決定するために、食品のような物体中の誘電関数を決定するための第２の実施形態を示す。
【００６９】
フローはステップ２１０で開始する。ここで、物体中の１組の点が選択されている。超音波メトリックを得る処理の間に用いられた点を選択することは有利である。選択された点は方程式１−１７に３および４に対応する。
【００７０】
超音波放射は、その後、ステップ２１１で両方の点に集中する。ステップ２１２で、Ｓ３１、Ｓ２３、Ｓ４１、Ｓ２４、Ｓ４’１、Ｓ２４’、Ｓ３’１およびＳ２３’のＳパラメータが、図７に関連して詳細に記述されるように、測定される。
【００７１】
ＳパラメータＳ４３（すなわち選択された点の間の減衰）は、ステップ２１３で計算される。点３は、この実施形態で、仮想送信機として作動し、点４は仮想受信機として機能する。
【００７２】
選択された点ｘ、ｙ（すなわち、方程式１−７の点３，４）の間での誘電関数の平均値
【数１】
【００７３】
は、その後、ステップ２１４で計算される。
【００７４】
ステップ２１５で、別の組の点が選択されるべきかどうかの決定が下される。別の組の点が選択されるべきである場合には、フローは、ステップ２１０に戻る。ここで、ステップ２１１から２１４までが繰り返される前に新しい組が選択される。選択されるべきでない場合には、フローは図３のステップ１０６に続く。ここで所望の物理的特性が計算される。
【００７５】
（発明の第１の使用）
図６は、本発明の第１の使用の機能を概略的に示す。１つの超音波メトリックｕ（ｘ、ｔ）が製品内のすべての点ｘに対して得られる場合、以下のステップの適用によりすべての点で誘電率を計算することが可能である。
【００７６】
１） 点３のうちの１つに超音波を集中させる。超音波は、送信アンテナ１から受信アンテナ２に送られたマイクロ波信号の振動数シフトに関する焦点にだけ影響することが知られており、これにより、側波帯（すなわちマイクロ波ベース周波数（ｆ１）±超音波周波数（ｆＵＳ））で信号を生成する。
【００７７】
２） 複数のサイドバンドの少なくとも１つで信号強度を測定する。両側波帯の中で信号強度が測定される場合、測定からのより信頼できる結果が得られる。受信アンテナ２で測定された信号強度は、次のように表現されてもよい：
【数２】
【００７８】
ここで、Ｓ２１は、測定ギャップにある製品１２によって引き起こされた減衰であり、Ｖ２（ｔ）は側波帯で測定された信号強度であり、Ｖ１（ｔ）は送信アンテナ１から送られた信号の信号強度である。Ｓ２３は受信アンテナ２への点３の間の減衰であり、α３は、超音波がマイクロ波側波帯信号（音響電気利得と称する）に変換される点３での効率を決定する係数である。ｕ３（ｘ、ｔ）は点３での超音波メトリックであり、Ｓ３１は送信アンテナ１および点３の間の減衰である。
【００７９】
第１の近似では、効率αは次のように表現することが可能である。
【数３】
【００８０】
ここで、Δεは、超音波放射ｙによって引き起こされる圧力波による誘電率の変化である。圧縮モジュールκに対して、
【数４】
【００８１】
の関係が確立する。Κの値は、当業者に知られていて、より詳細に検討しない。
【００８２】
３） 製品１２中の図６に３を表示する全ての所望の点に対して処理を繰り返す。
【００８３】
４） 逆散乱アルゴリズムでのすべての測定データを使用し、製品中の誘電関数の空間分布を計算する。
【００８４】
物体が、測定装置に関して、相対的に低速で移動し、かつ、下の関係を満たす場合、射出された超音波とマイクロ波放射の補償は考慮する必要がない。
【数５】
【００８５】
ｖｏｂｊは、測定ギャップ１３での物体移動の速度であり、ｔｍｅａｓは完全な処理ための測定時間であり、ｖＵＳは物体中の超音波の速度であり、ｆＵＳは超音波周波数であり、ｄＦｏｃａｌは焦点の直径である。
【００８６】
相対速度が高い場合、超音波の集束は移動を補償するために、測定ステップの間、物体の中に焦点を維持するために、超音波放射の調整を含んでいなければならない。さらに、ドップラー偏移を回避するために、
【数６】
【００８７】
である。
【００８８】
（発明の第２の使用）
図７ａ−７ｄは、製品中の点３および４の２点間の誘電率を計算する際の、本発明の第２の使用の主関数を示す。第１の点３はソースであると見なしてもよく、第２の点４は受信機であると見なしてもよい。
【００８９】
主関数は、図６に関連して記述したものと非常に似ているが、２つの焦点３および４が超音波放射によって同時に生成されたので、２つの上側波帯および、２つの下側波帯が生成されるという違いがある。第１の上側波帯および下側波帯は、図６に関連して記述されたものと同様であり、第２の上側波帯および下側波帯は、２倍の超音波周波数（すなわち、マイクロ波ベース周波数（ｆ１）±２＊超音波周波数（２ｆＵＳ））を持っている。同じ超音波周波数がこの目的に用いられる場合、第２のオーダーの側波帯を生成するために、２つの異なる超音波周波数を選ぶことが可能である。図１に関連して記述された装置は、この例で、第２の上側波帯および下側波帯に応じて調整される余分な側波帯パスが付加される必要がある。
【００９０】
次の関係は、それぞれ単一の仮想ソースとしての点３および４にとって確立されることが可能である。
【数７】
【００９１】
図７ｂにしたがって、焦点を３から３’へ移動させ、焦点を４から４’へ移動することによって、新しい関係が表現することができる。
【数８】
【００９２】
図７ａから、点３および４の間の求められた減衰を含む関係は次のように表現されてもよい。
【数９】
【００９３】
方程式６は、７ｘ７問題を解くのに用いられず、適切な近似と取り替えられる。方程式１６および１７を参照。
【００９４】
図７ｃは、３および４に対応する２重のソースの関係を示す。
【数１０】
【００９５】
点３’および４’の関係が次のように表現されてもよい。
【数１１】
【００９６】
方程式１０は７ｘ７問題と８ｘ８問題を解くのに用いられず、適切な近似と取り替えられる。８ｘ８問題のためには方程式１５を参照し、７ｘ７問題のためには方程式１６および１７を参照する。
【００９７】
次の関係は図７ａ−７ｃから明白である。
【数１２】
【００９８】
方程式１１−１４は複数のＳパラメータを除去するために用いられる。それは図７ｄに示されるようなＳパラメータになる。求められたＳ４３であるＳパラメータがあり、問題を解くための１０の方程式（すなわち、方程式１−１０）を必要とするいくつかの未知のＳパラメータと共に、全く興味のないＳ３’４’である１つのＳパラメータがある。
【００９９】
有限要素についてのツィエンキービッツ（Ｚｉｅｎｋｉｅｗｉｃｚ）によって導入されたトリックを適用することによって、点３と点４との間の減衰を見つけるのに必要とされる方程式の数を減らすことは可能である。
【０１００】
方程式１０は用いられず、次の近似が代わりに用いられる。
【数１３】
【０１０１】
２度ツィエンキービッツのトリック（それは方程式６および１０の必要を除去する）を適用することによって、必要とされる方程式の数を８つの方程式のみに減らすことは可能でさえある。複数の方程式の代わりに用いられる近似は次のとおりである。
【数１４】
【０１０２】
点３と点４との間および点３’と点４’との間の減衰Ｓ４３は、必要とされる方程式を対数にすることにより計算することが可能であり、複数の解析点が適切に選ばれる限り解は常に利用可能である方程式と同数の未知を備えた方程式１から１０が方程式の非同次線形システムになる。これらの点の役割を「仮想プローブ」として示す点４と点３の間のマイクロ波ランタイムを得るためにＳ４３のシステムを解かねばならない。
【０１０３】
上に記述されたシステムは、「仮想送信機」（すなわち、点３）および「仮想受信機」（すなわち、点４）を使用する。発明の第１の使用法に到達するそれぞれ本物の送信または受信アンテナと一致するために、これらの点の１つを容易に置くことが可能である。複数の物理的なプローブアンテナの場所に両仮想プローブを置くことは、本発明に先立って既知の従来からのマイクロ波測定技術になる。
【０１０４】
解決される物理的問題に依存して、単一仮想プローブ概念（仮想受信機または仮想送信機）または２つの仮想プローブ概念を利用する。複数の仮想プローブによって生成／受信された特定のビームパターンを生じるために、複数のプローブのセット（例えば、仮想プローブのアレイ）を用いることは可能である。
【０１０５】
異なるプローブ構成は、掃海、物質解析、鉱物探査、医学的応用などとして応用に用いられてもよい。
【０１０６】
（方法の短縮した数学的導出）
電磁放射は、マクスウェル方程式によって統率され、ベクトルの電場Ｅは、次に示すように座標に依存する３次元の空間ｘおよび時間ｔで書かれているヘルムホルツ形式に容易に変形される。
【数１５】
【０１０７】
ここで、Δはラプラス演算子であり、ε０は真空の誘電率であり、εｒ（３ｘ３テンソルである）は与えられた位置での材料の局所的で相対的な誘電関数であり、μ０は真空透磁率を表し、μｒは検査中の材料の局所比透磁率を表す。この短縮した導出では、μｒは単位テンソル１（３ｘ３）に設定する。当業者にとっては、類似の方法が同時にεｒとμｒを解くことによって導出されることは明白である。
【０１０８】
同時に、テンソルの３ｘ３応力振幅ｙおよび媒体の局所音速ｖを備えた超音波も、次に示す同様の形式にすることが可能である。
【数１６】
【０１０９】
両方の微分方程式の解は、放射ソースの位置を考慮して実行される。処理の要点に集中させて、ゼロにならない振幅を備えたいかなる超音波も（圧縮または剪断型である）材料に応力を生じる。この応力は材料の局所圧縮によって反映される。この圧縮によって、分極電荷の密度は影響を受ける。既知の事実として、誘電性の物体のいかなる圧縮も相対的な誘電関数テンソルεｒを次のように変更する。
【数１７】
【０１１０】
この関係は、本発明で利用された超音波伝播と電磁波との間のカップリングを生じる。相互作用の強度は、音響電気光学の相互作用、３ｘ３ｘ３テンソルであるαによって決定される。含まれる物理学の全体像のために、上記の関係は（例えば、キャビテーションおよび他の非線形効果が無視される）比較的小さな超音波に対してだけ保たれることを言及しなければならない。
【０１１１】
その後、電磁気的に解決される完全なシステムは、次によって与えられる。
【数１８】
【０１１２】
当業者にとって、時間ｔでフーリエ変換が適用される場合、この種の微分方程式が周波数空間ωでの畳み込みになることは明白である。
【数１９】
【０１１３】
ここで、円で囲まれた積算演算子Ｅ（ｘ，ω）は、（畳み込み積分の前の、結果として生じる規格化定数を除いて）次の完全形になる（例えば、１９８８年にG.Doetschによる“Anleitung zum praktischen gebrauch der Laplace transformation”に見出される）周波数畳み込み積分を示す。
【数２０】
【０１１４】
したがって、単一の周波数超音波励起および物体に入射する単一の周波数マイクロ波信号を仮定すると、複数の受信マイクロ波信号は、入射するマイクロ波周波数に一部分を含むだけでなく、畳み込み積分によって生じた超音波およびマイクロ波周波数の差および和での複数の側波帯にも含む。
【０１１５】
上記の関係は、マイクロ波場から情報を抽出するために、超音波を適切に位相制御し、パルス波列を用いることによって、全く新しい世界を提示する。
【０１１６】
（単一の仮想プローブ）
単一の仮想プローブを含むパスに沿って解決するべき方法を適用する。これは、（点３または４のいずれかへの全ての関係が消える）上述した到来する線形方程式系を解くどんな送信パラメータにも依存する、仮想送信機または仮想受信機のどちらかに相当する。波の伝播メカニズムはこのケースには同一である。（一様な、境界条件無しの）理想的なケースについては、次の伝播関係に達する。
【数２１】
【０１１７】
（２つの仮想プローブ）
さらに、２つの仮想プローブを通るパスに沿って解決するために方法を適用することが可能である。これは、（すべて方程式が存在する）上述された到来する９ｘ９線形方程式系を解決するどんな伝送パラメータにも依存する、仮想送信機または仮想受信機のいずれかに相当する。（一様な、境界条件無しの）理想的なケースについては、次の伝播関係に達する。
【数２２】
【０１１８】
最初の２つの方程式は、仮想送信機の役割を務める解析点Ｘで側波帯の生成を示す。第３の方程式は、仮想受信機の役割を務める別の解析点Ｙで集中させることによって、最初の式の上部の第２の側波帯の生成を示す。周波数オフセットは、集中させることを遂行するために使用した超音波の周波数によって決定された、点Ｘでのηおよび点Ｙでのηで示される。これらが、ある応用では、両方の点Ｘ、Ｙにとって同じ周波数ではないことがあることに注意する。
【０１１９】
第１方程式は、側波帯オフセットＸを備えた所定の位置ξで側波帯の生成を述べる。第２の方程式は、強度ｑのソースが位置Ｘに置かれる場合、検査中の物体全体を通る側波帯の伝播を述べる。したがって、方法は、物体の内部の任意の複数の位置でのマイクロ波ソースを合成することによって、物体を「調査する」ことを可能にする。そして、このソースを動かす場合、このソースから生成された放射を測定する。
【０１２０】
発明は、マイクロ波発生器および超音波発生器に関連して記述されたが、物体内の密度変位を生じるために、他の型の放射が用いられてもよいことは明白である。しかしながら、放射は同時に射出されなければならず、さらに、変位を生じるために射出された放射間の周波数に差がなければならない。解像度は、物体で最短の波長を有する放射によって決定される。
【０１２１】
したがって、密度変位を生じ、かつ、それによって、発明を用いて材料の誘電関数を決定するために、異なる（例えばわずかに０．５Ｈｚだけ異なる）周波数を持っている２つのマイクロ波信号で物体を同時に照らすことは可能である。射出された放射の可能な組合せは、これに制限するわけではないが、マイクロ波、超音波およびＸ線の任意の組合せ含んでいる。
【図面の簡単な説明】
【０１２２】
【図１】発明のシステムを示す。
【図２】検査中の製品へ射出される放射線を示す。
【図３】検査中の製品の内部で、温度のような物理的特性を決定するためにフローチャートを示す。
【図４】超音波メトリックを得る処理を示すフローチャートを示す。
【図５ａ】検査中の製品内の誘電関数の空間分布を決定する処理の２つの実施形態を示すフローチャートを示す。
【図５ｂ】検査中の製品内の誘電関数の空間分布を決定する処理の２つの実施形態を示すフローチャートを示す。
【図６】本発明の第１の使用の主関数を示す。
【図７ａ】本発明の第２の使用の主関数を示す。
【図７ｂ】本発明の第２の使用の主関数を示す。
【図７ｃ】本発明の第２の使用の主関数を示す。
【図７ｄ】本発明の第２の使用の主関数を示す。
【特許請求の範囲】
【請求項１】
物体を通る第１型の放射を送信する少なくとも１つの第１の送信機（４２）と、前記送信された第１型の放射を受信する少なくとも１つの受信機（４３）と、を備える、前記物体中の誘電関数を決定するための装置であって、
前記物体を通る第２型の放射を射出し、前記第１および第２型の放射は、異なる周波数内容を有し前記物体での密度変化を生成する、少なくとも１つの第２の送信機と、
前記物体中の前記音響電気相互作用（δ）を決定するために、前記密度変化を介して送信される前記第１型の放射を解析する手段と、
前記決定された音響電気相互作用に基づいて、前記物体中の前記誘電関数を計算する計算手段と、を備えることを特徴とする装置。
【請求項２】
前記少なくとも１つの送信機（４２）に接続し、第１の固定周波数（ｆ１）を有する送信信号を生成し送信する第１の発生器（５１）をさらに備える請求項１に記載の装置。
【請求項３】
減衰を決定する手段をさらに備え、
該決定手段は、
前記少なくとも１つの受信機（４３）からの前記受信した第１型の放射と、第２の固定周波数（ｆ２）を有する局部発振器信号とを混合することによって中間周波数（ＩＦ）信号を生成し、前記局部発振器信号は第２の発生器（５２）によって生成される、ミキサと、
前記ＩＦ信号の位相と振幅を評価することによって、前記音響電気相互作用を決定する評価ユニットと、を備える請求項１または請求項２に記載の装置。
【請求項４】
前記射出された第１および第２型の放射は、前記物体に関して移動されるように配置される請求項１から請求項３のいずれか１項に記載の装置。
【請求項５】
前記装置を通って前記物体を運ぶコンベア（１１）をさらに備え、前記装置は静止している請求項４に記載の装置。
【請求項６】
前記装置は、静止物体に関して移動される請求項４に記載の装置。
【請求項７】
前記第２型の放射は、第３の固定周波数（ｆＵＳ）を有する信号で、第３の発生器（７１）によって生成された信号である請求項１から請求項７のいずれか１項に記載の装置。
【請求項８】
前記物体中の前記音響電気相互作用を決定するために使用される前記物体のメトリックに対応する、ランタイムおよび減衰のマッピングを決定するために、前記物体を通って射出される前記第２型の放射を受信する少なくとも１つの受信機（７３）をさらに備える請求項１から請求項７のいずれか１項に記載の装置。
【請求項９】
前記メトリックの１部分である焦点ごとに前記受信した第２型の放射の位相を決定する手段をさらに備える請求項８に記載の装置。
【請求項１０】
前記第１および第２型の放射は、マイクロ波放射、超音波、および、Ｘ線の任意の組合せを備える請求項１から請求項９のいずれか１項に記載の装置。
【請求項１１】
前記物体は食品であり、前記計算された誘電性関数に基づいて食品の中での温度の局所分布を計算する手段をさらに備える請求項１から請求項１０のいずれか１項に記載の装置。
【請求項１２】
少なくとも１つの第１の送信機（４２）から物体を通る第１型の放射を送信し、少なくとも１つの第１の受信機（４３）で前記送信された第１型の放射を受信するステップを備える、前記物体中の誘電関数を決定する方法であって、
少なくとも１つの第２の送信機から、前記物体を通る、第２型の放射を射出し、前記第１および第２型の放射は異なる周波数内容を有し、前記物体中での密度変化を生成していて、
前記物体での音響電気相互作用を決定するために、前記密度変化を介して送信される前記第１型の放射を解析し、
前記音響電気相互作用から、前記物体中の前記誘電関数を計算するステップを備えることを特徴とする方法。
【請求項１３】
前記物体中の前記音響電気相互作用を決定するために、前記第１型の放射を解析するステップは、前記物体のメトリックを得るステップを備える請求項１２に記載の方法。
【請求項１４】
前記メトリックを得るステップは、
ａ）射出された第２型の放射を物体中の点へ集中させ、
ｂ）最大の音響電気効率信号を得るために前記音響電気効率信号を測定している間に、前記第２型の放射の位相を調整し、
ｃ）焦点の位置と共に前記位相の値をメモリに格納し、
ｄ）前記物体のメトリックが完成するまで、ステップａ）からｃ）を繰り返す請求項１３に記載の方法。
【請求項１５】
前記物体中の前記誘電関数を計算するステップは、
前記物体の内部の少なくとも１点を選択し、
前記第２型の放射を前記少なくとも１点に集中させ、
前記受信した第１型の放射の減衰を決定し、
前記メトリックを使用して前記誘電関数を決定するステップを備える請求項１３または請求項１４に記載の方法。
【請求項１６】
前記物体中の誘電関数を計算するステップは、
物体の内部の少なくとも１組の点の選択し、
前記第２型の放射を前記少なくとも１組の点に集中させ、
前記少なくとも１組の点とっての第１型の放射の減衰を決定し、
メトリックを使用して前記少なくとも１組の点の間の前記減衰および前記誘電関数を決定するステップを備える請求項１３または請求項１４に記載の方法。
【請求項１７】
前記第１および第２型の放射は、マイクロ波放射、超音波、および、Ｘ線の任意の組合せであるように選択される請求項１２から請求項１６のいずれか１項に記載の方法。
【請求項１８】
食品を通る第１型の放射を送信する少なくとも１つの第１の送信機（４２）と、
前記送信された第１型の放射を受信する少なくとも１つの第１の受信機（４３）と、を備える、前記食品中の温度の局所分布を決定するための装置であって、
前記食品を通る超音波を射出し、前記第１および第２型の放射は、異なる周波数内容を有し、前記食品中の密度変化を生成する、少なくとも１つの第２の送信機と、
前記食品中の音響電気相互作用（δ）を決定するために、前記密度変化を介して送信される前記第１型の放射を解析する手段と、
前記音響電気相互作用に基づいて前記食品中の前記誘電関数を計算し、前記計算された誘電関数に基づいて前記食品中の温度の局所分布を計算する手段と、を備えることを特徴とする装置。
【請求項１９】
物体の特性を決定するための装置であって、
物体を通る第１型の放射と第２型の放射を送信し、前記第１および第２型の放射は、異なる周波数内容を有する送信ユニットと、
前記物体の特性を決定するために、前記送信された第２型の放射によって引き起こされた前記物体中の密度変化を介して送信される前記第１型の放射を解析する評価ユニットと、を備えることを特徴とする装置。
【請求項２０】
前記評価ユニットは、前記物体中の音響電気相互作用を決定し、前記決定された音響電気相互作用に基づいて前記物体中の誘電関数を計算する請求項１９に記載の装置。
【請求項２１】
前記物体の特性は前記物体の温度分布を含み、前記評価ユニットは前記計算された誘電関数に基づいて前記物体の前記温度分布を計算する請求項２０に記載の装置。
【請求項２２】
前記送信ユニットは、前記物体を通る前記第１型の放射を送信する少なくとも１つの第１の送信アンテナと、前記物体を通る前記第２型の放射を送信する少なくとも１つの第２の送信アンテナとを含む請求項１９から請求項２１のいずれか１項に記載の装置。
【請求項２３】
前記少なくとも第１の送信アンテナに接続し、第１の固定周波数を有する送信信号を生成し送信する第１の発生器をさらに備える請求項２２に記載の装置。
【請求項２４】
前記物体を通って送信される前記第１型の放射を受信する受信機と、
前記受信された第１型の放射と、第２の固定周波数を有する局部発振器信号とを混合することによって中間周波数（ＩＦ）信号を生成し、前記局部発振器信号は第２の発生器によって生成される、ミキサと、をさらに備え、
前記評価ユニットは、前記ＩＦ信号の位相および振幅を評価することによって前記音響電気相互作用を決定する請求項２２または請求項２３に記載の装置。
【請求項２５】
前記装置を通って前記物体を運ぶコンベヤをさらに備え、前記装置は静止している請求項１９から請求項２４のいずれか１項に記載の装置。
【請求項２６】
前記装置は、前記送信された第１型の放射と第２型の放射が前記物体に関して移動するように、静止物体に関して移動される請求項１９から請求項２４のいずれか１項に記載の装置。
【請求項２７】
前記第２型の放射は、第３の発生器によって生成される第３の固定周波数を有する信号である請求項１９から請求項２６のいずれか１項に記載の装置。
【請求項２８】
前記物体中の音響電気相互作用を決定するために前記評価ユニットによって使用される前記物体のメトリックに対応する、ランタイムおよび減衰のマッピングを決定するために、前記物体を通って送信される前記第２型の放射を受信する少なくとも１つの受信アンテナをさらに備える請求項１９から請求項２７のいずれか１項に記載の装置。
【請求項２９】
前記評価ユニットは、前記メトリックの一部分である焦点ごとに、第２型の放射の位相を決定する請求項２８に記載の装置。
【請求項３０】
前記第１および第２型の放射は、マイクロ波放射、超音波、および、Ｘ線の任意の組合せを備える請求項１９から請求項２９のいずれか１項に記載の装置。
【請求項１】
物体を通る第１型の放射を送信する少なくとも１つの第１の送信機（４２）と、前記送信された第１型の放射を受信する少なくとも１つの受信機（４３）と、を備える、前記物体中の誘電関数を決定するための装置であって、
前記物体を通る第２型の放射を射出し、前記第１および第２型の放射は、異なる周波数内容を有し前記物体での密度変化を生成する、少なくとも１つの第２の送信機と、
前記物体中の前記音響電気相互作用（δ）を決定するために、前記密度変化を介して送信される前記第１型の放射を解析する手段と、
前記決定された音響電気相互作用に基づいて、前記物体中の前記誘電関数を計算する計算手段と、を備えることを特徴とする装置。
【請求項２】
前記少なくとも１つの送信機（４２）に接続し、第１の固定周波数（ｆ１）を有する送信信号を生成し送信する第１の発生器（５１）をさらに備える請求項１に記載の装置。
【請求項３】
減衰を決定する手段をさらに備え、
該決定手段は、
前記少なくとも１つの受信機（４３）からの前記受信した第１型の放射と、第２の固定周波数（ｆ２）を有する局部発振器信号とを混合することによって中間周波数（ＩＦ）信号を生成し、前記局部発振器信号は第２の発生器（５２）によって生成される、ミキサと、
前記ＩＦ信号の位相と振幅を評価することによって、前記音響電気相互作用を決定する評価ユニットと、を備える請求項１または請求項２に記載の装置。
【請求項４】
前記射出された第１および第２型の放射は、前記物体に関して移動されるように配置される請求項１から請求項３のいずれか１項に記載の装置。
【請求項５】
前記装置を通って前記物体を運ぶコンベア（１１）をさらに備え、前記装置は静止している請求項４に記載の装置。
【請求項６】
前記装置は、静止物体に関して移動される請求項４に記載の装置。
【請求項７】
前記第２型の放射は、第３の固定周波数（ｆＵＳ）を有する信号で、第３の発生器（７１）によって生成された信号である請求項１から請求項７のいずれか１項に記載の装置。
【請求項８】
前記物体中の前記音響電気相互作用を決定するために使用される前記物体のメトリックに対応する、ランタイムおよび減衰のマッピングを決定するために、前記物体を通って射出される前記第２型の放射を受信する少なくとも１つの受信機（７３）をさらに備える請求項１から請求項７のいずれか１項に記載の装置。
【請求項９】
前記メトリックの１部分である焦点ごとに前記受信した第２型の放射の位相を決定する手段をさらに備える請求項８に記載の装置。
【請求項１０】
前記第１および第２型の放射は、マイクロ波放射、超音波、および、Ｘ線の任意の組合せを備える請求項１から請求項９のいずれか１項に記載の装置。
【請求項１１】
前記物体は食品であり、前記計算された誘電性関数に基づいて食品の中での温度の局所分布を計算する手段をさらに備える請求項１から請求項１０のいずれか１項に記載の装置。
【請求項１２】
少なくとも１つの第１の送信機（４２）から物体を通る第１型の放射を送信し、少なくとも１つの第１の受信機（４３）で前記送信された第１型の放射を受信するステップを備える、前記物体中の誘電関数を決定する方法であって、
少なくとも１つの第２の送信機から、前記物体を通る、第２型の放射を射出し、前記第１および第２型の放射は異なる周波数内容を有し、前記物体中での密度変化を生成していて、
前記物体での音響電気相互作用を決定するために、前記密度変化を介して送信される前記第１型の放射を解析し、
前記音響電気相互作用から、前記物体中の前記誘電関数を計算するステップを備えることを特徴とする方法。
【請求項１３】
前記物体中の前記音響電気相互作用を決定するために、前記第１型の放射を解析するステップは、前記物体のメトリックを得るステップを備える請求項１２に記載の方法。
【請求項１４】
前記メトリックを得るステップは、
ａ）射出された第２型の放射を物体中の点へ集中させ、
ｂ）最大の音響電気効率信号を得るために前記音響電気効率信号を測定している間に、前記第２型の放射の位相を調整し、
ｃ）焦点の位置と共に前記位相の値をメモリに格納し、
ｄ）前記物体のメトリックが完成するまで、ステップａ）からｃ）を繰り返す請求項１３に記載の方法。
【請求項１５】
前記物体中の前記誘電関数を計算するステップは、
前記物体の内部の少なくとも１点を選択し、
前記第２型の放射を前記少なくとも１点に集中させ、
前記受信した第１型の放射の減衰を決定し、
前記メトリックを使用して前記誘電関数を決定するステップを備える請求項１３または請求項１４に記載の方法。
【請求項１６】
前記物体中の誘電関数を計算するステップは、
物体の内部の少なくとも１組の点の選択し、
前記第２型の放射を前記少なくとも１組の点に集中させ、
前記少なくとも１組の点とっての第１型の放射の減衰を決定し、
メトリックを使用して前記少なくとも１組の点の間の前記減衰および前記誘電関数を決定するステップを備える請求項１３または請求項１４に記載の方法。
【請求項１７】
前記第１および第２型の放射は、マイクロ波放射、超音波、および、Ｘ線の任意の組合せであるように選択される請求項１２から請求項１６のいずれか１項に記載の方法。
【請求項１８】
食品を通る第１型の放射を送信する少なくとも１つの第１の送信機（４２）と、
前記送信された第１型の放射を受信する少なくとも１つの第１の受信機（４３）と、を備える、前記食品中の温度の局所分布を決定するための装置であって、
前記食品を通る超音波を射出し、前記第１および第２型の放射は、異なる周波数内容を有し、前記食品中の密度変化を生成する、少なくとも１つの第２の送信機と、
前記食品中の音響電気相互作用（δ）を決定するために、前記密度変化を介して送信される前記第１型の放射を解析する手段と、
前記音響電気相互作用に基づいて前記食品中の前記誘電関数を計算し、前記計算された誘電関数に基づいて前記食品中の温度の局所分布を計算する手段と、を備えることを特徴とする装置。
【請求項１９】
物体の特性を決定するための装置であって、
物体を通る第１型の放射と第２型の放射を送信し、前記第１および第２型の放射は、異なる周波数内容を有する送信ユニットと、
前記物体の特性を決定するために、前記送信された第２型の放射によって引き起こされた前記物体中の密度変化を介して送信される前記第１型の放射を解析する評価ユニットと、を備えることを特徴とする装置。
【請求項２０】
前記評価ユニットは、前記物体中の音響電気相互作用を決定し、前記決定された音響電気相互作用に基づいて前記物体中の誘電関数を計算する請求項１９に記載の装置。
【請求項２１】
前記物体の特性は前記物体の温度分布を含み、前記評価ユニットは前記計算された誘電関数に基づいて前記物体の前記温度分布を計算する請求項２０に記載の装置。
【請求項２２】
前記送信ユニットは、前記物体を通る前記第１型の放射を送信する少なくとも１つの第１の送信アンテナと、前記物体を通る前記第２型の放射を送信する少なくとも１つの第２の送信アンテナとを含む請求項１９から請求項２１のいずれか１項に記載の装置。
【請求項２３】
前記少なくとも第１の送信アンテナに接続し、第１の固定周波数を有する送信信号を生成し送信する第１の発生器をさらに備える請求項２２に記載の装置。
【請求項２４】
前記物体を通って送信される前記第１型の放射を受信する受信機と、
前記受信された第１型の放射と、第２の固定周波数を有する局部発振器信号とを混合することによって中間周波数（ＩＦ）信号を生成し、前記局部発振器信号は第２の発生器によって生成される、ミキサと、をさらに備え、
前記評価ユニットは、前記ＩＦ信号の位相および振幅を評価することによって前記音響電気相互作用を決定する請求項２２または請求項２３に記載の装置。
【請求項２５】
前記装置を通って前記物体を運ぶコンベヤをさらに備え、前記装置は静止している請求項１９から請求項２４のいずれか１項に記載の装置。
【請求項２６】
前記装置は、前記送信された第１型の放射と第２型の放射が前記物体に関して移動するように、静止物体に関して移動される請求項１９から請求項２４のいずれか１項に記載の装置。
【請求項２７】
前記第２型の放射は、第３の発生器によって生成される第３の固定周波数を有する信号である請求項１９から請求項２６のいずれか１項に記載の装置。
【請求項２８】
前記物体中の音響電気相互作用を決定するために前記評価ユニットによって使用される前記物体のメトリックに対応する、ランタイムおよび減衰のマッピングを決定するために、前記物体を通って送信される前記第２型の放射を受信する少なくとも１つの受信アンテナをさらに備える請求項１９から請求項２７のいずれか１項に記載の装置。
【請求項２９】
前記評価ユニットは、前記メトリックの一部分である焦点ごとに、第２型の放射の位相を決定する請求項２８に記載の装置。
【請求項３０】
前記第１および第２型の放射は、マイクロ波放射、超音波、および、Ｘ線の任意の組合せを備える請求項１９から請求項２９のいずれか１項に記載の装置。
【図１】
【図２】
【図３】
【図４】
【図５ａ】
【図５ｂ】
【図６】
【図７ａ】
【図７ｂ】
【図７ｃ】
【図７ｄ】
【図２】
【図３】
【図４】
【図５ａ】
【図５ｂ】
【図６】
【図７ａ】
【図７ｂ】
【図７ｃ】
【図７ｄ】
【公表番号】特表２００８−５１９９７９（Ｐ２００８−５１９９７９Ａ）
【公表日】平成２０年６月１２日（２００８．６．１２）
【国際特許分類】
【出願番号】特願２００７−５４１１４２（Ｐ２００７−５４１１４２）
【出願日】平成１７年１１月１０日（２００５．１１．１０）
【国際出願番号】ＰＣＴ／ＳＥ２００５／００１６９４
【国際公開番号】ＷＯ２００６／０５２２０２
【国際公開日】平成１８年５月１８日（２００６．５．１８）
【出願人】（５００４６１３８７）エフエムシー・フードテク・エービー (4)
【Ｆターム（参考）】
【公表日】平成２０年６月１２日（２００８．６．１２）
【国際特許分類】
【出願日】平成１７年１１月１０日（２００５．１１．１０）
【国際出願番号】ＰＣＴ／ＳＥ２００５／００１６９４
【国際公開番号】ＷＯ２００６／０５２２０２
【国際公開日】平成１８年５月１８日（２００６．５．１８）
【出願人】（５００４６１３８７）エフエムシー・フードテク・エービー (4)
【Ｆターム（参考）】
[ Back to top ]