被検体情報取得装置
【課題】被検体情報取得装置で超音波を用いて変位を測定するために合成開口法を用いる際の精度を向上させる。
【解決手段】被検体に弾性波を送信し、反射波を受信する装置であって、弾性波と電気信号とを変換可能な素子を少なくとも１方向に複数配列して備える送受信器と、素子に電気信号を入力して弾性波を送信させる素子制御手段と、素子が受信する反射波を検出する検出手段とを有し、素子に入力される電気信号は、複数の素子間において符号化された符号化パルス信号であり、検出手段は、反射波を復号化するとともに、２箇所を焦点とした双曲線および楕円を軸とした両軸の交点について復号化された反射波を合成する合成開口処理を異なる時点で行い、被検体の少なくとも２方向における変位を取得する被検体情報取得装置を用いる。
【解決手段】被検体に弾性波を送信し、反射波を受信する装置であって、弾性波と電気信号とを変換可能な素子を少なくとも１方向に複数配列して備える送受信器と、素子に電気信号を入力して弾性波を送信させる素子制御手段と、素子が受信する反射波を検出する検出手段とを有し、素子に入力される電気信号は、複数の素子間において符号化された符号化パルス信号であり、検出手段は、反射波を復号化するとともに、２箇所を焦点とした双曲線および楕円を軸とした両軸の交点について復号化された反射波を合成する合成開口処理を異なる時点で行い、被検体の少なくとも２方向における変位を取得する被検体情報取得装置を用いる。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波を用いた被検体情報取得装置に関する。
【背景技術】
【０００２】
超音波診断装置は、組織内の音響インピーダンス分布を反映した反射波から生体内組織の構造を映像化する超音波断層像（Ｂ−ｍｏｄｅ像）を表示するだけでなく、血流や組織の運動などを可視化する機能を有していることが一般的である。
【０００３】
さらに最近では超音波により組織部分の硬さを推定して組織診断へ利用することが行われ始めている。体表から静的な圧力を加える、比較的低周波数の超音波により加振を与える、などの手法によって生じる組織内部の歪みを超音波によって計測し、この歪みの大小から組織の硬さを推定する。このような歪みは、超音波によって計測した各点での変位分布を各点間の距離で除することで算出される。
【０００４】
上記のような組織の運動や組織部分の硬さを算出する際には、超音波によって計測した変位を用いることが多く、近年の超音波診断装置では超音波による変位計測が重要な技術課題となっている。
【０００５】
特許文献１には、実効的に双曲線−楕円上に配列したサンプル点でフォーカスしたエコー信号列を用いて変位計測を行う超音波診断装置が開示されている。また、実効的に双曲線−楕円上に配列したサンプル点にフォーカスするために合成開口法を用いる場合、信号ＳＮ比を向上させることを目的として仮想波源を設定し複数の素子で送信を行う手法が開示されている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】国際公開第２０１０／０５３１５６号
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、変位を測定するために仮想波源を設定し複数の素子で送信を行う場合、合成開口法にとって最適な送信条件ではない場合が存在する。
【０００８】
合成開口法を用いる場合、送信波は球面波に近いほど得られる空間分解能が向上する。また、超音波の送受信器（探触子）の素子が配列された面に垂直な方向に対して角度の大きい方向の観察が可能で、例えば２次元断面内の画像化においては円柱状の波面を形成することが好ましい。つまり超音波診断装置において、合成開口法を用いる場合、そのＳＮ比の低下さえ問わなければ、最適な送信方法は１素子だけで送信する手法である。
【０００９】
しかしながら、仮想波源を用いて送信を行う場合、複数素子を用いて送信するため、１素子で送信する場合と比較して理想的な球面波もしくは円柱状の波面を形成することが困難である。そのため、高いＳＮ比の維持と合成開口法にとって最適な送信条件とを両立することには限界があった。
【００１０】
本発明は上記の課題に鑑みてなされたものであり、被検体情報取得装置で超音波を用いて変位を測定するために合成開口法を用いる際の精度を向上させることを目的とする。
【課題を解決するための手段】
【００１１】
本発明は以下の構成を採用する。すなわち、被検体に弾性波を送信し、該送信された弾性波の被検体での反射波を受信して該被検体の情報を取得する被検体情報取得装置であって、弾性波と電気信号とを変換可能な素子を少なくとも１方向に複数配列して備える送受信器と、前記素子に電気信号を入力し、前記素子から前記被検体に弾性波を送信させる素子制御手段と、前記素子が受信する前記送信された弾性波の前記被検体での反射波を検出する検出手段と、を有し、前記素子制御手段が前記素子に入力する電気信号は、前記複数の素子間において符号化された符号化パルス信号であり、前記検出手段は、前記反射波を復号化するとともに、複数の素子が配列された前記１方向における２箇所を焦点とした双曲線に沿った方向と該２箇所を焦点とする楕円に沿った方向とを軸とし、該両軸の交点箇所または該交点に対応する箇所について前記復号化された反射波を合成する合成開口処理を異なる時点において行い、該異なる時点での合成開口処理結果に基づき、前記被検体の少なくとも２方向における変位を取得することを特徴とする被検体情報取得装置である。
【００１２】
本発明はまた、以下の構成を採用する。すなわち、被検体に弾性波を送信し、該送信された弾性波の被検体での反射波を受信して該被検体の情報を取得する被検体情報取得装置であって、弾性波と電気信号とを変換可能な素子を複数配列して備える送受信器と、前記素子に電気信号を入力し、前記素子から前記被検体に弾性波を送信させる素子制御手段と、前記素子が受信する前記送信された弾性波の前記被検体での反射波を検出する検出手段と、を有し、前記素子制御手段が前記素子に入力する電気信号は、前記複数の素子間において符号化された符号化パルス信号であり、前記検出手段は、前記反射波を復号化するとともに、前記復号された反射波を合成する合成開口処理を異なる時点において行い、該異なる時点での合成開口処理結果に基づき、前記被検体の少なくとも２方向における変位を取得することを特徴とする被検体情報取得装置である。
【発明の効果】
【００１３】
本発明によれば、被検体情報取得装置で超音波を用いて変位を測定するために合成開口法を用いる際の精度を向上させることができる。
【図面の簡単な説明】
【００１４】
【図１】被検体情報取得装置のシステム概要を示す図。
【図２】符号列信号を生成する処理を示すフローチャート。
【図３】符号列信号による素子の動作を示す図。
【図４】受信した信号に対する処理を示すフローチャート。
【図５】双曲線と楕円を用いた変調を示す図。
【図６】複数の素子に符号列信号を適用した様子を示す図。
【図７】変位算出の結果の例を示す図。
【発明を実施するための形態】
【００１５】
以下に図面を参照しつつ、本発明の好適な実施の形態を説明する。ただし、以下に記載されている構成部品の寸法、材質、形状及びそれらの相対配置などは、発明が適用される装置の構成や各種条件により適宜変更されるべきものであり、この発明の範囲を以下の記載に限定する趣旨のものではない。
【００１６】
以下の記載においては、本発明の被検体情報取得装置の例として、超音波診断装置について説明する。本発明の超音波診断装置は、被検体に超音波を送信し、被検体内部で反射した超音波（エコー波）を受信して、被検体情報を画像データとして取得する超音波エコー技術を利用した装置である。被検体としては例えば、人や動物等の生体の一部や、生体を模した材料を使用できる。取得される被検体情報は、被検体内部の組織の音響インピーダンスの違いを反映した情報である。本発明でいう超音波とは、一種の弾性波であり、音
波、超音波、音響波などとも呼ばれる。
【００１７】
＜実施例１＞
図１は本発明による被検体情報取得装置のシステム概要を示す図である。本図を用いて本発明の装置の信号の流れと各ブロックの役割を述べ、その後、それぞれの処理について詳細に述べる。
【００１８】
（システム概要と信号の流れ）
システム制御部００１からの制御信号に従って、符号制御ブロック００２は符号列信号を発生させて送信回路系００５へ入力する。送信回路系００５では入力された符号列信号に従って電圧波形を持つ電気信号を生成する。この電圧波形は複数の超音波変換素子００３によって超音波（弾性波）に変換され、探触子（送受信器）００４から被検体０００内部へと送信される。超音波送信の場面において、システム制御部は、本発明の素子制御手段に相当する。
【００１９】
被検体０００内部で反射された超音波（反射波）は、複数の超音波変換素子００３によって検出され、複数の電圧信号（電気信号）に変換され、受信信号として受信回路系００６に入力される。なお、超音波変換素子は、超音波（弾性波）と電気信号を変換可能であればどのようなものでも良い。アポダイゼーション処理を行うために、本発明における超音波の送受信器として、かかる超音波変換素子が少なくとも１方向に複数配列されたものを用いる。
受信回路系００６では複数の電圧信号を増幅し、複数のデジタル信号に変換する。受信回路系００６から出力されたデジタル信号は復号合成ブロック００７に入力される。
【００２０】
復号合成ブロック００７では、入力された複数のデジタル信号と、符号制御ブロック００２から入力された符号列信号に関する情報と、システム制御部から入力された位置情報と変調情報とを用いて復号処理及び合成開口処理を行う。後述する通常アポダイゼーションの結果を第一の合成開口信号として画像生成ブロック００８へ、変調アポダイゼーションの結果を第二の合成開口信号として変位算出ブロック００９へと出力する。
【００２１】
画像生成ブロック００８では入力された第一の合成開口信号を元に包絡線を算出し包絡線データとして画像処理ブロック０１０へと出力する。なお、画像生成ブロック００８では、入力された第一の合成開口信号に対してバンドパスフィルタをかけるなどの各種処理を必要に応じて行ってもよい。
【００２２】
変位算出ブロック００９では異なる時点で取得された受信信号を用いて算出された複数の第二の合成開口信号を用いて変位を算出し、変位情報データとして画像処理ブロック０１０へと出力する。復号合成ブロックと変位算出ブロックが、本発明における検出手段に相当する。
【００２３】
画像処理ブロック０１０では入力された包絡線データに対して強度調整や各種フィルタ処理を行い、組織の音響インピーダンスの違いを反映した分布情報として、インピーダンス輝度データ（いわゆるＢ−ｍｏｄｅ像）を画像表示系０１１に出力する。さらに画像処理ブロック０１０は入力された変位情報データを用いて、変位情報輝度データとして画像表示系０１１に出力する。
【００２４】
画像表示系０１１は入力された変位情報輝度データならびにインピーダンス輝度データを、システム制御部００１からの指示に従い表示する。表示方法は、変位情報輝度データ及びインピーダンス輝度データを重畳して表示しても良いし、並べて表示してもよい。もちろん、夫々の輝度データのみを表示してもよい。表示モードは、システム制御部００１
からの指示によって変更可能である。
【００２５】
なお変位算出ブロック００９は変位を算出するだけでなく、さらにその変位を用いて歪み分布を算出することも可能であり、その歪み分布を歪み情報データとして画像処理ブロック０１０へ出力、さらに画像表示系０１１において表示してもよい。
以上が、本発明による被検体情報取得装置の基本的な構成ならびに信号の流れである。
【００２６】
（符号列信号により制御された超音波送信）
次に符号制御ブロック００２ならびに送信回路系００５の動作について述べる。
符号制御ブロック００２では送信に用いる素子数に適した符号列信号を生成する。図２のフローチャートを用いて符号列信号の生成について説明する。
【００２７】
まず、Ｇｏｌａｙ符号の組｛ａＮ（ｉ）｝、｛ｂＮ（ｉ）｝（ｉ＝１，・・・，２Ｎ）を準備する（ステップＳ２０１）。ここでＧｏｌａｙ符号の組はそれぞれの符号列ａＮ（ｉ）、ｂＮ（ｉ）の自己相関数Ａｊ、Ｂｊが以下の式（１）、式（２）を満たす符号列の組であり、２Ｎは符号列の長さを示す。
【数１】
【００２８】
次に、さらにこれらの符号列に直交するＧｏｌａｙ符号の組を、式（３）のように生成する。（ステップＳ２０２）これは、式（４）で表される要素の組み合わせである。
【数２】
【００２９】
ここで、互いに直交なＧｏｌａｙ符号から構成される、式（５）のような行列Ｇ２を作成する。（ステップＳ２０３）
【数３】
【００３０】
必要な符号の組み合わせ数になるように、Ｈａｄａｍａｒｄ行列とのクロネッカー積を求め、式（６）のように符号列信号ＧＭを生成する。（ステップＳ２０４）
【数４】
【００３１】
なお、Ｈａｄａｍａｒｄ行列は以下の式（７）、式（８）によって算出できる。
【数５】
【００３２】
例として４つの素子に対し符号列信号を用いてどのように動作するかについて、図３を用いて説明する。図中には、各素子に入力される信号の電圧波形が模式的に示されている。
符号制御ブロック００２は先述した手法によって式（９）のような符号列信号（ここではＧ４）を生成する。Ｇ４は、Ｎ＝０のときの符号列信号である。符号列信号は、本発明の符号化パルス信号に相当する。
【数６】
【００３３】
符号制御ブロック００２から送信された符号列信号は送信回路系００５に入力される。送信回路系００５は、４つの素子３０１、素子３０２、素子３０３、素子３０４へ、符号列信号の各列に対応した電圧波形を送信する。
【００３４】
まず、送信回路系００５は、符号列信号の１列目に並んだ（１，１，１，１）を用いてそれぞれの素子に同位相（例えばこの位相を０度の位相とする）の電圧波形を送信する（図３（ａ））。
次に符号列信号の２列目に並んだ（１，−１，１，−１）に従って、素子３０１に位相０度、素子３０２に位相１８０度、素子３０３に位相０度、素子３０４に位相１８０度の電圧波形を送信する（図３（ｂ））。
さらに符号列信号の３列目に並んだ（１，１，−１，−１）に従って、素子３０１に位相０度、素子３０２に位相０度、素子３０３に位相１８０度、素子３０４に位相１８０度の電圧波形を送信する（図３（ｃ））。
最後に符号列信号の４列目に並んだ（１，−１，−１，１）に従って、素子３０１に位相０度、素子３０２に位相１８０度、素子３０３に位相１８０度、素子３０４に位相０度の電圧波形を送信する（図３（ｄ））。
なお図中で黒い素子は位相１８０度の電圧波形を入力されていることを示している。それぞれの送信ごとに、被検体内部で超音波が反射して反射波が発生し、超音波変換素子００３で受信されてアナログの電気信号（受信信号）に変換される。
【００３５】
（受信信号の復号処理）
次に受信回路系００６ならびに復号合成ブロック００７について述べる。
受信回路系００６は、それぞれの送信ごとに複数の超音波変換素子００３で受信された電圧信号を増幅し、複数のデジタル信号に変換する。受信回路系００６から出力されたデジタル信号は復号合成ブロック００７に入力される。
【００３６】
図４を用いて復号合成ブロック００７の動作について述べる。
復号合成ブロック００７では入力された複数のデジタル信号をメモリに保存する（ステ
ップＳ４０１）。所定回数の送信が終了した時点で復号処理を行う（ステップＳ４０２）。
ここで復号化された式（１０）の信号ＲＴｘ，Ｒｘ（ｔ）は、Ｔｘ番目の素子で送信しＲｘ番目の素子で受信したデジタル信号である。このとき、ｇｉｊは符号列信号の（ｉ，ｊ）要素、ｒｉｊ（ｔ）はｉ番目の素子がｊ回目の送信時に受信したデジタル信号である。
【数７】
【００３７】
ここで復号処理について、先述した４素子の送信を例として、より詳細に説明を行う。
符号列信号を用いて４回送信された超音波が被検体内部から反射して戻ってきた受信波形によるデジタル信号がｒｉｊ（ｔ）である。ここでｉは受信した素子を表し（ｉ＝１が素子３０１、ｉ＝２が素子３０２、ｉ＝３が素子３０３、ｉ＝４が素子３０４）、ｊは何回目の送信に対する受信信号であるか（ここではｊ＝１，２，３，４）を意味する。この複数のデジタル信号ｒｉｊ（ｔ）に対して、符号列信号Ｇ４を用いて復号化を行う。
【００３８】
例えば符号列信号の１行目に並んだ（１，１，１，１）を用いて１番目の素子（素子３０１）で送信し、１番目の素子（素子３０１）で受信した信号Ｒ１，１（ｔ）は式（１１）のように表される。
【数８】
【００３９】
この時、２番目の素子（素子３０２）で送信され、被検体内部で反射して戻ってきた受信波形がどのようになっているのかを述べる。２番目の素子からの超音波は（１，−１，１，−１）に従って送信されている。先述した（１，１，１，１）を用いて復号化した場合、１×１＋１×（−１）＋１×１＋１×（−１）＝０となり、送信した超音波が復号化によって０となる。同様に３番目の素子から送信された超音波、４番目の素子から送信された超音波も復号化によって０となる。
つまり、上記のように符号列信号を用いて復号化を行うことで、１番目の素子から送信した超音波の反射波だけを抽出することが出来る。
【００４０】
ここまでの処理の結果をまとめる。符号列信号に従って４つの素子から同時に超音波を送信し、その受信波形を復号化することで、それぞれの素子から送信した超音波の反射波だけを抽出している。これは複数の素子間において符号化された信号を用いることで、反射波を受信後に分離していることになる。復号化した結果、それぞれの素子から送信し、全ての素子で受信した信号（この例では１６種類）を個別に得ることが出来る。さらに復号化したことによって、４回送受信を繰り返したのと同じＳＮ比を実現している。言い換えれば、１素子ずつの送信を実施した場合には、同等のＳＮ比を全ての素子について得るために合計１６回の送信を行う必要があるところ、本発明では４回の送信で実現している。
したがって本発明によれば、合成開口法に適した１素子ずつの送信を実施した場合と同等の信号を、同じ送信回数で、よりＳＮ比を高く得ることができる。
【００４１】
なお、ここではＮ＝０の符号列信号を例にして説明したが、Ｎは正の整数のいずれでも良く、その場合はさらにＳＮ比向上の効果が得られる。
【００４２】
（合成開口処理）
さらに復号合成ブロック００７の動作について説明を続ける。
復号化によって得られた復号化された信号ＲＴｘ，Ｒｘ（ｔ）を用いて、合成開口処理を行う（ステップＳ４０３）。合成開口処理においては、共通する２箇所を焦点とする双曲線群および楕円群の交点をサンプリング点とする。言い換えると、超音波変換素子が配列された１方向における２箇所を焦点とする双曲線に沿った方向と、同じ焦点を持つ楕円に沿った方向とに両軸を取り、その両軸の交点箇所をサンプリング点とする。そして、それらの交点箇所に送受信フォーカスをかけたデータを合成する。
【００４３】
ここでの双曲線群は、２焦点からの超音波パルスの到達時間の差が実質的に同一となるように選ばれた曲線または直線の集合である。この双曲線群内の隣接する該曲線または直線間では、該２焦点からの超音波パルスの到達時間の差の値が、使用超音波波長のｍ分の１だけ異なる。楕円群は、該２焦点からの超音波パルスの到達時間の和が実質的に同一となるように選ばれた曲線の集合である。この楕円群内の隣接する曲線間では、該２焦点からの超音波パルスの到達時間の和の値が、使用超音波の波長のｎ分の１だけ異なる。なお、ｍ、ｎはそれぞれ整数であり、後に述べる変位測定精度の向上のために４以上となることが好ましい。
【００４４】
このようにして定められたサンプリング点は隣接するサンプリング点に対して、双曲線方向、楕円方向共に各々一定の位相差で配列されている。そのため、後に述べる変位測定の際に、位相間隔が揃ったデータとして扱うことが可能で、計算手順が簡略化できる。
このような等間隔の位相差を有するサンプリング点の配置方法の一例を示す。
下式（１２）は双曲線インデックスｉ、楕円インデックスｊ（ｉ＝０，±１，±２，・・・、ｊ＝１，２，３，・・・）に対する各点の座標ｘｉｊ，ｙｉｊを示すものである。
【数９】
ここでｃは音速、Ｔはサンプリング間隔、ｆはプローブの中心からアポダイゼーションのピーク（後述する）までの距離である。
【００４５】
これにより得られる信号はｘ−ｙ平面から楕円−双曲線平面に座標変換したものとなる。図５はその座標変化の概念を説明する図である。
図５（ａ）はアポダイゼーションのピークＰ１、Ｐ２から連続的に超音波を発生し、それらが干渉した際の様子を表している。このような干渉が起きているｘ−ｙ平面において、先述した座標ｘｉｊ，ｙｉｊのデータを取得し、双曲線インデックスｉを横軸に、楕円インデックスｊを縦軸にしてプロットすると、図５（ｂ）のように位相間隔が揃ったデータが得られることが分かる。さらに、直交している２軸（横軸（双曲線インデックス軸）と縦軸（楕円インデックス軸））にそれぞれ独立に変調がかかっているため、各々の軸の波の成分に分離して変位推定を行うことで変位推定の処理規模を抑制することも可能である。
【００４６】
この合成開口処理を行う際に、一般的な超音波装置で行われている通常アポダイゼーシ
ョンと、変位を測定するための変調アポダイゼーションの２種類のアポダイゼーションを用いる。ここでいうアポダイゼーションとは送信や受信の際に複数の素子に対して与える重みであり、通常アポダイゼーションとしては、ガウシアン、ハミング、ハニングなどを用いることが可能である。
【００４７】
ここで変位を測定するための変調アポダイゼーションについて説明する。
複数の素子で超音波を送信する際には、送信に用いた複数の素子に与えた重み（アポダイゼーション）のフーリエ変換が、焦点付近における超音波送受信方向に対して垂直な方向の音圧分布となる。逆に、焦点付近で超音波送受信方向に対して垂直な方向に正弦波状の変調を生じさせたいときは、その変調の逆フーリエ変換のアポダイゼーションを与えれば良い。
【００４８】
例えば、送信時にはｗｔ（ｘ）、受信時にはｗｒ（ｘ）の変調アポダイゼーションを用いることで変調が可能である。それぞれ式（１３）、式（１４）で表す。
【数１０】
ここで、ｘ０＝ｙλｆｘ、σ０＝（ｙλ√２）／σｘ、ｙは焦点を形成する深さ、λは超音波の波長、ｆｘは焦点付近での変調の周波数、σｘは焦点付近で変調された正弦波のガウス包絡線のＦＷＨＭ（半値幅）である。
【００４９】
なお、変調アポダイゼーションのかけ方はこれに限らない。例えば、送受信の内いずれか一方で用いる複数の素子によって形成される開口の中に２箇所の副開口を仮想的に形成し、それらの副開口から可干渉性を有する超音波を送信もしくは受信処理しても、同様に変調を生じさせることが可能である。言い換えれば２つのピークを有する重みを用いることで、変調アポダイゼーションが可能となる。
【００５０】
通常アポダイゼーション（ステップＳ４０４）によって得られた信号は第一の合成開口信号として画像生成ブロック００８へ（ステップＳ４０５）出力される。変調アポダイゼーション（ステップＳ４０６）で得られた信号は第二の合成開口信号として変位算出ブロック００９へ（ステップＳ４０７）と出力される。
なお、通常アポダイゼーションを用いて合成開口する場合は、共通する２焦点を有する双曲線群、楕円群の交点をサンプリング点とせず、プローブ開口から直線的に進む直線上にサンプリング点を設定しても良い。
【００５１】
（変位の算出）
次に変位算出ブロック００９での動作について述べる。
変位算出ブロック００９には、異なる時点で取得された受信信号を用いて算出された合成開口処理結果である、複数の第二の合成開口信号が入力される。
【００５２】
ここでは複数の第二の合成開口信号、すなわち異なる時点での合成開口処理結果間から変位を算出する方法を説明する。
例として双曲線インデックス軸と楕円インデックス軸とにそれぞれ独立にＩＱ信号を求め変位を算出する方法に関して述べる。
双曲線インデックス軸上の位置をｕ、楕円インデックス軸上の位置をｖとした時、入力された第二の合成開口信号は以下の式（１５）のようにモデル化出来る。
【数１１】
Ａは第二の合成開口信号の包絡線、ｆｕ、ｆｖはそれぞれ双曲線インデックス軸、楕円インデックス軸方向の変調周波数、φｕ、φｖはそれぞれ双曲線インデックス軸、楕円インデックス軸方向の位相を表す。
【００５３】
ここで、双曲線インデックス軸方向の波のＩＱ信号Ｉｕ（ｕ，ｖ）、Ｑｕ（ｕ，ｖ）、楕円インデックス軸方向の波のＩＱ信号Ｉｖ（ｕ，ｖ）、Ｑｖ（ｕ，ｖ）は、以下の式（１６）〜式（１９）で算出できる。
【数１２】
【００５４】
ここで、ＩＩ（ｕ，ｖ），ＩＱ（ｕ，ｖ），ＱＩ（ｕ，ｖ），ＱＱ（ｕ，ｖ）は、以下の式（２０）〜式（２３）で表される。
【数１３】
ＬＰＦはカットオフ周波数をｆｕおよびｆｖとするローパスフィルタを意味する。
【００５５】
このとき包絡線Ａ（ｕ，ｖ）と双曲線インデックス軸方向の位相φｕ、楕円インデックス軸方向の位相φｖは以下の式（２４）〜式（２６）で求められる。
【数１４】
【００５６】
異なるタイミングで得られた第二の合成開口信号それぞれの位相をφｕ１、φｕ２、φｖ１、φｖ２とする。このとき、双曲線インデックス軸方向の変位δｕ、楕円インデックス軸方向の変位δｖという２方向での変位は、以下の式（２７）、式（２８）で算出できる。
【数１５】
ここでλｕ、λｖはそれぞれ双曲線インデックス軸方向、楕円インデックス軸方向での波長である。
【００５７】
このようにして必要な観察範囲内の変位分布を算出可能である。このように第二の合成開口信号のＩＱデータを用いて変位を算出できるのは、サンプリングされた点が双曲線方向、楕円方向共に各々一定の位相差で配列されている効果であり、計算規模の低減に寄与する。また、符号列信号を用いたＳＮ比の高い信号に対して変位算出処理を行うため、より高い変位推定精度を実現することが出来る。
【００５８】
（変位算出の変形例１）
上記の例では、双曲線インデックス軸方向、楕円インデックス軸方向に分離したＩＱ信号を用いる手法について述べたが、複数の波形信号から位相情報や変位情報を抽出する手法を用いても構わない。
【００５９】
例えば、以下のような手法でも変位を算出することが出来る。
異なる時間で取得した第二の合成開口信号の包絡線を求め、それらの包絡線の相互相関によって大まかな変位（δｕ’，δｖ’）を推定する。大まかな変位分を補正した後、さらに複素相関によって位相差（φｕ’，φｖ’）を算出する。さらに別途注目位置での瞬時周波数（ｆｕ，ｆｖ）を推定する。
【００６０】
双曲線インデックス軸方向の変位δｕ、楕円インデックス軸方向の変位δｖは以下の式（２９）、式（３０）で算出できる。ここでｃは音速である。
【数１６】
【００６１】
この手法では瞬時周波数を推定し、その推定結果を用いることにより、観察領域内の双曲線群や楕円群のずれによる周波数の変化に対しても安定して変位を求めることが可能である。さらにその瞬時周波数を推定することによって精度向上の効果も得られる。
なお、この手法では観測個所ごとに瞬時周波数を推定するため双曲線群と楕円群との交点でサンプルされた信号を用いなくとも、横方向の変位が測定可能である。その場合であっても、符号化送受信によるＳＮ比向上の効果を得ることが可能である。
【００６２】
（変位算出の変形例２）
また、さらに以下のような手法でも変位を算出することが出来る。
異なる時間で取得した第二の合成開口信号の包絡線を求め、それらの包絡線の相互相関によって大まかな変位（δｕ’，δｖ’）を推定する。大まかな変位を補正した第二の合成開口信号をそれぞれｓ１（ｕ，ｖ）、ｓ２（ｕ，ｖ）として、任意のずれ量Δｕ、Δｖを与えて複素相関の偏角Ｒ（Δｕ，Δｖ：ｕ，ｖ）を算出する。これを式（３１）に示す。
【数１７】
【００６３】
この複素相関の偏角Ｒ（Δｕ，Δｖ：ｕ，ｖ）は、変位（δｕ’＋Δｕ，δｖ’＋Δｖ）が真の変位と一致した時に０になる。局所的には偏角が一次関数で変化すると近似して、偏角Ｒ（Δｕ，Δｖ：ｕ，ｖ）が０付近の値を取る複数の（Δｕ，Δｖ）で線形補間を行い、確からしい変位を算出する。
【００６４】
この手法では変位そのものを推定するため、観察領域内の双曲線群や楕円群のずれによる周波数の変化に対しても安定して変位を求めることが可能であり、変位推定の精度向上の効果が期待できる。なお、この手法では観測個所ごとに変位そのものを推定するため、双曲線群と楕円群との交点でサンプルされた信号を用いなくとも、横方向の変位が測定可能であり、その場合であっても符号化送受信によるＳＮ比向上の効果を得ることが可能である。
【００６５】
（歪みの算出）
次にこの変位データ（δｕ，δｖ）を用いて歪み分布を求める手法について述べる。
変位データを楕円−双曲線平面からｘ−ｙ平面に座標変換し、変位ベクトルとして算出する。この変位ベクトルをそれぞれｘ方向、ｙ方向に微分することで歪み分布を算出することが出来る。なお、楕円−双曲線平面内での歪み分布を先に算出し、その後ｘ−ｙ平面に座標変換しても同様の結果が得られる。
【００６６】
なお、歪み分布を算出する際には微分処理を実施するため空間周波数の高いノイズが発
生する可能性がある。このノイズを低減するために、画像処理ブロック０１０において、算出した歪み分布に対して空間的なローパスフィルタを適用するなどの視認性向上を目的とした処理を行っても良い。
【００６７】
ここまでに述べたように、本発明によれば符号列信号を用いることで合成開口法に適した１素子毎の送信を実施した場合と同等の信号を、よりＳＮ比を高く得ることができる。さらに、そのＳＮ比の高い信号を用いて変位推定するため、より変位推定精度の高い被検体情報処理装置を得ることが出来る。
【００６８】
図７（ａ）は本発明を用いて推定した歪み分布を示す図である。寒天内に周辺よりも硬さの硬い直径９ｍｍの内包物（図中点線で囲まれた領域）を設置し、図示しない探触子（図面上側に接触している）に対して垂直な方向（図面右方向）から０．５％程度の圧縮を加える。圧縮を加える前後での変位を本発明によって算出し、横方向の歪み分布を算出する。図７（ａ）はこの歪み分布をグレイスケールによって示している。
【００６９】
図７（ｂ）は図７（ａ）中の２箇所（２直線上。ｙ＝２９．５ｍｍ、３８．０ｍｍ）における横方向の歪み分布をプロットしたものである。ｙ＝２９．５ｍｍと記された直線上の横方向の歪み分布は実線で、ｙ＝３８．０ｍｍと記された直線上の横方向の歪み分布は点線で示している。点線で示したｙ＝３８．０ｍｍの直線上における横方向の歪み分布はほぼ−０．５％で一定しており、横方向から０．５％程度で圧縮されている状態を計測出来ている。また、実線で示したｙ＝２９．５ｍｍの直線上における横方向の歪み分布においては、中央付近に歪みの小さい領域（−０．２％程度）が存在しており、内包物が周辺の寒天よりも硬いことが計測出来ていることが分かる。
【００７０】
＜実施例２＞
符号列信号に対して位相を変化させた電圧波形を出力する際には、必ずしも一つの符号に対して電圧波形は１波長分である必要はない。すなわち、位相を変化させた電圧波形が、複数の波長分の長さを有していても良い。
【００７１】
このように複数の波長分の長さを有している電圧波形を用いた場合、１波長分の電圧波形を用いた場合に比べて、さらにＳＮ比の向上ならびに変位推定の精度向上の効果が得られる。
なお、電圧波形の長さが長くなるに従って、変位推定の空間分解能が低下するため、先述した空間的なローパスフィルタによって決定される空間分解能以下もしくは同等程度までに電圧波形の長さを抑えることが望ましい。
【００７２】
＜実施例３＞
符号列信号のそれぞれの値を１素子ずつに適用するのではなく複数の素子に対して適用することも可能である。
【００７３】
図６を用いて符号列信号のそれぞれの値を複数の素子に対して適用した例を示す。符号列信号は先述した例と同様のものを用いる。
このとき、送受信器の素子は図示したようにグループ化されている。素子６０１から素子６０３は素子グループＡ（６２１）、素子６０４から素子６０６は素子グループＢ（６２２）、素子６０７から素子６０９は素子グループＣ（６２３）、素子６１０から素子６１２は素子グループＤ（６２４）である。
【００７４】
符号制御ブロック００２から送信された符号列信号は送信回路系００５に入力される。送信回路系００５では１２個の素子６０１から素子６１２へ電圧波形の送信を行う。
まず、符号列信号の１列目に並んだ（１，１，１，１）を用いて、素子グループＡ、素
子グループＢ、素子グループＣ、素子グループＤに同位相（例えばこの位相を０度の位相とする）の電圧波形を送信する。
次に符号列信号の２列目に並んだ（１，−１，１，−１）に従って、素子グループＡに位相０度、素子グループＢに位相１８０度、素子グループＣに位相０度、素子グループＤに位相１８０度の電圧波形を送信する。
【００７５】
このように、符号列信号のそれぞれの値を複数の素子（ここでは３素子）ごとに適用し、超音波を被検体に送信する。さらに被検体内部で反射された超音波を受信し、復号化を行う。復号化した結果は先述した１素子毎の送信を行った場合と同等の信号とは異なり、素子グループＡ、素子グループＢ、素子グループＣ、素子グループＤでそれぞれ送信した信号が分離される。
これらの復号化された信号を用いて合成開口処理を行うステップ以降は先述した内容と同等のため説明を省略する。
【００７６】
このように複数の素子に対して同じ符号列信号を用いて送信を行うことで、送受信に用いる素子数よりも送受信回数を減らすことが可能であり、より高速にデータを取得することが出来る。そのため、より変化の早い変位や歪みを安定に精度良く測定することが可能になる。
【００７７】
なお、符号列信号のそれぞれの値が適用された複数の素子（例えば素子グループＡ）は、仮想的に波源を送受信器（探触子）の内部や被検体内に設定することで、より球面波もしくは円柱状の波面を形成することが可能である。これにより、複数の素子を同時に送信する場合と比較してさらに合成開口による空間分解能を向上させ、変位測定の精度を向上させることが可能である。
【００７８】
なお、上記実施形態１で説明した、変位計測の変形例１、２においては、上記の通り、変位計測の測定誤差を抑える、信号レベルを向上させるという点でより好ましい例である。ただし、これらの例において、計測の測定誤差を抑えるという点においては、必ずしも符号化パルスを用いなくても良い。つまり、送信パルスが符号化されたパルスではない場合においても、変位計測の測定誤差を抑えるという効果を奏するものである。
【符号の説明】
【００７９】
００１：システム制御部，００２：符号制御ブロック，００３：超音波変換素子，００７：復号合成ブロック，００８：画像生成ブロック，００９：変位算出ブロック
【技術分野】
【０００１】
本発明は、超音波を用いた被検体情報取得装置に関する。
【背景技術】
【０００２】
超音波診断装置は、組織内の音響インピーダンス分布を反映した反射波から生体内組織の構造を映像化する超音波断層像（Ｂ−ｍｏｄｅ像）を表示するだけでなく、血流や組織の運動などを可視化する機能を有していることが一般的である。
【０００３】
さらに最近では超音波により組織部分の硬さを推定して組織診断へ利用することが行われ始めている。体表から静的な圧力を加える、比較的低周波数の超音波により加振を与える、などの手法によって生じる組織内部の歪みを超音波によって計測し、この歪みの大小から組織の硬さを推定する。このような歪みは、超音波によって計測した各点での変位分布を各点間の距離で除することで算出される。
【０００４】
上記のような組織の運動や組織部分の硬さを算出する際には、超音波によって計測した変位を用いることが多く、近年の超音波診断装置では超音波による変位計測が重要な技術課題となっている。
【０００５】
特許文献１には、実効的に双曲線−楕円上に配列したサンプル点でフォーカスしたエコー信号列を用いて変位計測を行う超音波診断装置が開示されている。また、実効的に双曲線−楕円上に配列したサンプル点にフォーカスするために合成開口法を用いる場合、信号ＳＮ比を向上させることを目的として仮想波源を設定し複数の素子で送信を行う手法が開示されている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】国際公開第２０１０／０５３１５６号
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、変位を測定するために仮想波源を設定し複数の素子で送信を行う場合、合成開口法にとって最適な送信条件ではない場合が存在する。
【０００８】
合成開口法を用いる場合、送信波は球面波に近いほど得られる空間分解能が向上する。また、超音波の送受信器（探触子）の素子が配列された面に垂直な方向に対して角度の大きい方向の観察が可能で、例えば２次元断面内の画像化においては円柱状の波面を形成することが好ましい。つまり超音波診断装置において、合成開口法を用いる場合、そのＳＮ比の低下さえ問わなければ、最適な送信方法は１素子だけで送信する手法である。
【０００９】
しかしながら、仮想波源を用いて送信を行う場合、複数素子を用いて送信するため、１素子で送信する場合と比較して理想的な球面波もしくは円柱状の波面を形成することが困難である。そのため、高いＳＮ比の維持と合成開口法にとって最適な送信条件とを両立することには限界があった。
【００１０】
本発明は上記の課題に鑑みてなされたものであり、被検体情報取得装置で超音波を用いて変位を測定するために合成開口法を用いる際の精度を向上させることを目的とする。
【課題を解決するための手段】
【００１１】
本発明は以下の構成を採用する。すなわち、被検体に弾性波を送信し、該送信された弾性波の被検体での反射波を受信して該被検体の情報を取得する被検体情報取得装置であって、弾性波と電気信号とを変換可能な素子を少なくとも１方向に複数配列して備える送受信器と、前記素子に電気信号を入力し、前記素子から前記被検体に弾性波を送信させる素子制御手段と、前記素子が受信する前記送信された弾性波の前記被検体での反射波を検出する検出手段と、を有し、前記素子制御手段が前記素子に入力する電気信号は、前記複数の素子間において符号化された符号化パルス信号であり、前記検出手段は、前記反射波を復号化するとともに、複数の素子が配列された前記１方向における２箇所を焦点とした双曲線に沿った方向と該２箇所を焦点とする楕円に沿った方向とを軸とし、該両軸の交点箇所または該交点に対応する箇所について前記復号化された反射波を合成する合成開口処理を異なる時点において行い、該異なる時点での合成開口処理結果に基づき、前記被検体の少なくとも２方向における変位を取得することを特徴とする被検体情報取得装置である。
【００１２】
本発明はまた、以下の構成を採用する。すなわち、被検体に弾性波を送信し、該送信された弾性波の被検体での反射波を受信して該被検体の情報を取得する被検体情報取得装置であって、弾性波と電気信号とを変換可能な素子を複数配列して備える送受信器と、前記素子に電気信号を入力し、前記素子から前記被検体に弾性波を送信させる素子制御手段と、前記素子が受信する前記送信された弾性波の前記被検体での反射波を検出する検出手段と、を有し、前記素子制御手段が前記素子に入力する電気信号は、前記複数の素子間において符号化された符号化パルス信号であり、前記検出手段は、前記反射波を復号化するとともに、前記復号された反射波を合成する合成開口処理を異なる時点において行い、該異なる時点での合成開口処理結果に基づき、前記被検体の少なくとも２方向における変位を取得することを特徴とする被検体情報取得装置である。
【発明の効果】
【００１３】
本発明によれば、被検体情報取得装置で超音波を用いて変位を測定するために合成開口法を用いる際の精度を向上させることができる。
【図面の簡単な説明】
【００１４】
【図１】被検体情報取得装置のシステム概要を示す図。
【図２】符号列信号を生成する処理を示すフローチャート。
【図３】符号列信号による素子の動作を示す図。
【図４】受信した信号に対する処理を示すフローチャート。
【図５】双曲線と楕円を用いた変調を示す図。
【図６】複数の素子に符号列信号を適用した様子を示す図。
【図７】変位算出の結果の例を示す図。
【発明を実施するための形態】
【００１５】
以下に図面を参照しつつ、本発明の好適な実施の形態を説明する。ただし、以下に記載されている構成部品の寸法、材質、形状及びそれらの相対配置などは、発明が適用される装置の構成や各種条件により適宜変更されるべきものであり、この発明の範囲を以下の記載に限定する趣旨のものではない。
【００１６】
以下の記載においては、本発明の被検体情報取得装置の例として、超音波診断装置について説明する。本発明の超音波診断装置は、被検体に超音波を送信し、被検体内部で反射した超音波（エコー波）を受信して、被検体情報を画像データとして取得する超音波エコー技術を利用した装置である。被検体としては例えば、人や動物等の生体の一部や、生体を模した材料を使用できる。取得される被検体情報は、被検体内部の組織の音響インピーダンスの違いを反映した情報である。本発明でいう超音波とは、一種の弾性波であり、音
波、超音波、音響波などとも呼ばれる。
【００１７】
＜実施例１＞
図１は本発明による被検体情報取得装置のシステム概要を示す図である。本図を用いて本発明の装置の信号の流れと各ブロックの役割を述べ、その後、それぞれの処理について詳細に述べる。
【００１８】
（システム概要と信号の流れ）
システム制御部００１からの制御信号に従って、符号制御ブロック００２は符号列信号を発生させて送信回路系００５へ入力する。送信回路系００５では入力された符号列信号に従って電圧波形を持つ電気信号を生成する。この電圧波形は複数の超音波変換素子００３によって超音波（弾性波）に変換され、探触子（送受信器）００４から被検体０００内部へと送信される。超音波送信の場面において、システム制御部は、本発明の素子制御手段に相当する。
【００１９】
被検体０００内部で反射された超音波（反射波）は、複数の超音波変換素子００３によって検出され、複数の電圧信号（電気信号）に変換され、受信信号として受信回路系００６に入力される。なお、超音波変換素子は、超音波（弾性波）と電気信号を変換可能であればどのようなものでも良い。アポダイゼーション処理を行うために、本発明における超音波の送受信器として、かかる超音波変換素子が少なくとも１方向に複数配列されたものを用いる。
受信回路系００６では複数の電圧信号を増幅し、複数のデジタル信号に変換する。受信回路系００６から出力されたデジタル信号は復号合成ブロック００７に入力される。
【００２０】
復号合成ブロック００７では、入力された複数のデジタル信号と、符号制御ブロック００２から入力された符号列信号に関する情報と、システム制御部から入力された位置情報と変調情報とを用いて復号処理及び合成開口処理を行う。後述する通常アポダイゼーションの結果を第一の合成開口信号として画像生成ブロック００８へ、変調アポダイゼーションの結果を第二の合成開口信号として変位算出ブロック００９へと出力する。
【００２１】
画像生成ブロック００８では入力された第一の合成開口信号を元に包絡線を算出し包絡線データとして画像処理ブロック０１０へと出力する。なお、画像生成ブロック００８では、入力された第一の合成開口信号に対してバンドパスフィルタをかけるなどの各種処理を必要に応じて行ってもよい。
【００２２】
変位算出ブロック００９では異なる時点で取得された受信信号を用いて算出された複数の第二の合成開口信号を用いて変位を算出し、変位情報データとして画像処理ブロック０１０へと出力する。復号合成ブロックと変位算出ブロックが、本発明における検出手段に相当する。
【００２３】
画像処理ブロック０１０では入力された包絡線データに対して強度調整や各種フィルタ処理を行い、組織の音響インピーダンスの違いを反映した分布情報として、インピーダンス輝度データ（いわゆるＢ−ｍｏｄｅ像）を画像表示系０１１に出力する。さらに画像処理ブロック０１０は入力された変位情報データを用いて、変位情報輝度データとして画像表示系０１１に出力する。
【００２４】
画像表示系０１１は入力された変位情報輝度データならびにインピーダンス輝度データを、システム制御部００１からの指示に従い表示する。表示方法は、変位情報輝度データ及びインピーダンス輝度データを重畳して表示しても良いし、並べて表示してもよい。もちろん、夫々の輝度データのみを表示してもよい。表示モードは、システム制御部００１
からの指示によって変更可能である。
【００２５】
なお変位算出ブロック００９は変位を算出するだけでなく、さらにその変位を用いて歪み分布を算出することも可能であり、その歪み分布を歪み情報データとして画像処理ブロック０１０へ出力、さらに画像表示系０１１において表示してもよい。
以上が、本発明による被検体情報取得装置の基本的な構成ならびに信号の流れである。
【００２６】
（符号列信号により制御された超音波送信）
次に符号制御ブロック００２ならびに送信回路系００５の動作について述べる。
符号制御ブロック００２では送信に用いる素子数に適した符号列信号を生成する。図２のフローチャートを用いて符号列信号の生成について説明する。
【００２７】
まず、Ｇｏｌａｙ符号の組｛ａＮ（ｉ）｝、｛ｂＮ（ｉ）｝（ｉ＝１，・・・，２Ｎ）を準備する（ステップＳ２０１）。ここでＧｏｌａｙ符号の組はそれぞれの符号列ａＮ（ｉ）、ｂＮ（ｉ）の自己相関数Ａｊ、Ｂｊが以下の式（１）、式（２）を満たす符号列の組であり、２Ｎは符号列の長さを示す。
【数１】
【００２８】
次に、さらにこれらの符号列に直交するＧｏｌａｙ符号の組を、式（３）のように生成する。（ステップＳ２０２）これは、式（４）で表される要素の組み合わせである。
【数２】
【００２９】
ここで、互いに直交なＧｏｌａｙ符号から構成される、式（５）のような行列Ｇ２を作成する。（ステップＳ２０３）
【数３】
【００３０】
必要な符号の組み合わせ数になるように、Ｈａｄａｍａｒｄ行列とのクロネッカー積を求め、式（６）のように符号列信号ＧＭを生成する。（ステップＳ２０４）
【数４】
【００３１】
なお、Ｈａｄａｍａｒｄ行列は以下の式（７）、式（８）によって算出できる。
【数５】
【００３２】
例として４つの素子に対し符号列信号を用いてどのように動作するかについて、図３を用いて説明する。図中には、各素子に入力される信号の電圧波形が模式的に示されている。
符号制御ブロック００２は先述した手法によって式（９）のような符号列信号（ここではＧ４）を生成する。Ｇ４は、Ｎ＝０のときの符号列信号である。符号列信号は、本発明の符号化パルス信号に相当する。
【数６】
【００３３】
符号制御ブロック００２から送信された符号列信号は送信回路系００５に入力される。送信回路系００５は、４つの素子３０１、素子３０２、素子３０３、素子３０４へ、符号列信号の各列に対応した電圧波形を送信する。
【００３４】
まず、送信回路系００５は、符号列信号の１列目に並んだ（１，１，１，１）を用いてそれぞれの素子に同位相（例えばこの位相を０度の位相とする）の電圧波形を送信する（図３（ａ））。
次に符号列信号の２列目に並んだ（１，−１，１，−１）に従って、素子３０１に位相０度、素子３０２に位相１８０度、素子３０３に位相０度、素子３０４に位相１８０度の電圧波形を送信する（図３（ｂ））。
さらに符号列信号の３列目に並んだ（１，１，−１，−１）に従って、素子３０１に位相０度、素子３０２に位相０度、素子３０３に位相１８０度、素子３０４に位相１８０度の電圧波形を送信する（図３（ｃ））。
最後に符号列信号の４列目に並んだ（１，−１，−１，１）に従って、素子３０１に位相０度、素子３０２に位相１８０度、素子３０３に位相１８０度、素子３０４に位相０度の電圧波形を送信する（図３（ｄ））。
なお図中で黒い素子は位相１８０度の電圧波形を入力されていることを示している。それぞれの送信ごとに、被検体内部で超音波が反射して反射波が発生し、超音波変換素子００３で受信されてアナログの電気信号（受信信号）に変換される。
【００３５】
（受信信号の復号処理）
次に受信回路系００６ならびに復号合成ブロック００７について述べる。
受信回路系００６は、それぞれの送信ごとに複数の超音波変換素子００３で受信された電圧信号を増幅し、複数のデジタル信号に変換する。受信回路系００６から出力されたデジタル信号は復号合成ブロック００７に入力される。
【００３６】
図４を用いて復号合成ブロック００７の動作について述べる。
復号合成ブロック００７では入力された複数のデジタル信号をメモリに保存する（ステ
ップＳ４０１）。所定回数の送信が終了した時点で復号処理を行う（ステップＳ４０２）。
ここで復号化された式（１０）の信号ＲＴｘ，Ｒｘ（ｔ）は、Ｔｘ番目の素子で送信しＲｘ番目の素子で受信したデジタル信号である。このとき、ｇｉｊは符号列信号の（ｉ，ｊ）要素、ｒｉｊ（ｔ）はｉ番目の素子がｊ回目の送信時に受信したデジタル信号である。
【数７】
【００３７】
ここで復号処理について、先述した４素子の送信を例として、より詳細に説明を行う。
符号列信号を用いて４回送信された超音波が被検体内部から反射して戻ってきた受信波形によるデジタル信号がｒｉｊ（ｔ）である。ここでｉは受信した素子を表し（ｉ＝１が素子３０１、ｉ＝２が素子３０２、ｉ＝３が素子３０３、ｉ＝４が素子３０４）、ｊは何回目の送信に対する受信信号であるか（ここではｊ＝１，２，３，４）を意味する。この複数のデジタル信号ｒｉｊ（ｔ）に対して、符号列信号Ｇ４を用いて復号化を行う。
【００３８】
例えば符号列信号の１行目に並んだ（１，１，１，１）を用いて１番目の素子（素子３０１）で送信し、１番目の素子（素子３０１）で受信した信号Ｒ１，１（ｔ）は式（１１）のように表される。
【数８】
【００３９】
この時、２番目の素子（素子３０２）で送信され、被検体内部で反射して戻ってきた受信波形がどのようになっているのかを述べる。２番目の素子からの超音波は（１，−１，１，−１）に従って送信されている。先述した（１，１，１，１）を用いて復号化した場合、１×１＋１×（−１）＋１×１＋１×（−１）＝０となり、送信した超音波が復号化によって０となる。同様に３番目の素子から送信された超音波、４番目の素子から送信された超音波も復号化によって０となる。
つまり、上記のように符号列信号を用いて復号化を行うことで、１番目の素子から送信した超音波の反射波だけを抽出することが出来る。
【００４０】
ここまでの処理の結果をまとめる。符号列信号に従って４つの素子から同時に超音波を送信し、その受信波形を復号化することで、それぞれの素子から送信した超音波の反射波だけを抽出している。これは複数の素子間において符号化された信号を用いることで、反射波を受信後に分離していることになる。復号化した結果、それぞれの素子から送信し、全ての素子で受信した信号（この例では１６種類）を個別に得ることが出来る。さらに復号化したことによって、４回送受信を繰り返したのと同じＳＮ比を実現している。言い換えれば、１素子ずつの送信を実施した場合には、同等のＳＮ比を全ての素子について得るために合計１６回の送信を行う必要があるところ、本発明では４回の送信で実現している。
したがって本発明によれば、合成開口法に適した１素子ずつの送信を実施した場合と同等の信号を、同じ送信回数で、よりＳＮ比を高く得ることができる。
【００４１】
なお、ここではＮ＝０の符号列信号を例にして説明したが、Ｎは正の整数のいずれでも良く、その場合はさらにＳＮ比向上の効果が得られる。
【００４２】
（合成開口処理）
さらに復号合成ブロック００７の動作について説明を続ける。
復号化によって得られた復号化された信号ＲＴｘ，Ｒｘ（ｔ）を用いて、合成開口処理を行う（ステップＳ４０３）。合成開口処理においては、共通する２箇所を焦点とする双曲線群および楕円群の交点をサンプリング点とする。言い換えると、超音波変換素子が配列された１方向における２箇所を焦点とする双曲線に沿った方向と、同じ焦点を持つ楕円に沿った方向とに両軸を取り、その両軸の交点箇所をサンプリング点とする。そして、それらの交点箇所に送受信フォーカスをかけたデータを合成する。
【００４３】
ここでの双曲線群は、２焦点からの超音波パルスの到達時間の差が実質的に同一となるように選ばれた曲線または直線の集合である。この双曲線群内の隣接する該曲線または直線間では、該２焦点からの超音波パルスの到達時間の差の値が、使用超音波波長のｍ分の１だけ異なる。楕円群は、該２焦点からの超音波パルスの到達時間の和が実質的に同一となるように選ばれた曲線の集合である。この楕円群内の隣接する曲線間では、該２焦点からの超音波パルスの到達時間の和の値が、使用超音波の波長のｎ分の１だけ異なる。なお、ｍ、ｎはそれぞれ整数であり、後に述べる変位測定精度の向上のために４以上となることが好ましい。
【００４４】
このようにして定められたサンプリング点は隣接するサンプリング点に対して、双曲線方向、楕円方向共に各々一定の位相差で配列されている。そのため、後に述べる変位測定の際に、位相間隔が揃ったデータとして扱うことが可能で、計算手順が簡略化できる。
このような等間隔の位相差を有するサンプリング点の配置方法の一例を示す。
下式（１２）は双曲線インデックスｉ、楕円インデックスｊ（ｉ＝０，±１，±２，・・・、ｊ＝１，２，３，・・・）に対する各点の座標ｘｉｊ，ｙｉｊを示すものである。
【数９】
ここでｃは音速、Ｔはサンプリング間隔、ｆはプローブの中心からアポダイゼーションのピーク（後述する）までの距離である。
【００４５】
これにより得られる信号はｘ−ｙ平面から楕円−双曲線平面に座標変換したものとなる。図５はその座標変化の概念を説明する図である。
図５（ａ）はアポダイゼーションのピークＰ１、Ｐ２から連続的に超音波を発生し、それらが干渉した際の様子を表している。このような干渉が起きているｘ−ｙ平面において、先述した座標ｘｉｊ，ｙｉｊのデータを取得し、双曲線インデックスｉを横軸に、楕円インデックスｊを縦軸にしてプロットすると、図５（ｂ）のように位相間隔が揃ったデータが得られることが分かる。さらに、直交している２軸（横軸（双曲線インデックス軸）と縦軸（楕円インデックス軸））にそれぞれ独立に変調がかかっているため、各々の軸の波の成分に分離して変位推定を行うことで変位推定の処理規模を抑制することも可能である。
【００４６】
この合成開口処理を行う際に、一般的な超音波装置で行われている通常アポダイゼーシ
ョンと、変位を測定するための変調アポダイゼーションの２種類のアポダイゼーションを用いる。ここでいうアポダイゼーションとは送信や受信の際に複数の素子に対して与える重みであり、通常アポダイゼーションとしては、ガウシアン、ハミング、ハニングなどを用いることが可能である。
【００４７】
ここで変位を測定するための変調アポダイゼーションについて説明する。
複数の素子で超音波を送信する際には、送信に用いた複数の素子に与えた重み（アポダイゼーション）のフーリエ変換が、焦点付近における超音波送受信方向に対して垂直な方向の音圧分布となる。逆に、焦点付近で超音波送受信方向に対して垂直な方向に正弦波状の変調を生じさせたいときは、その変調の逆フーリエ変換のアポダイゼーションを与えれば良い。
【００４８】
例えば、送信時にはｗｔ（ｘ）、受信時にはｗｒ（ｘ）の変調アポダイゼーションを用いることで変調が可能である。それぞれ式（１３）、式（１４）で表す。
【数１０】
ここで、ｘ０＝ｙλｆｘ、σ０＝（ｙλ√２）／σｘ、ｙは焦点を形成する深さ、λは超音波の波長、ｆｘは焦点付近での変調の周波数、σｘは焦点付近で変調された正弦波のガウス包絡線のＦＷＨＭ（半値幅）である。
【００４９】
なお、変調アポダイゼーションのかけ方はこれに限らない。例えば、送受信の内いずれか一方で用いる複数の素子によって形成される開口の中に２箇所の副開口を仮想的に形成し、それらの副開口から可干渉性を有する超音波を送信もしくは受信処理しても、同様に変調を生じさせることが可能である。言い換えれば２つのピークを有する重みを用いることで、変調アポダイゼーションが可能となる。
【００５０】
通常アポダイゼーション（ステップＳ４０４）によって得られた信号は第一の合成開口信号として画像生成ブロック００８へ（ステップＳ４０５）出力される。変調アポダイゼーション（ステップＳ４０６）で得られた信号は第二の合成開口信号として変位算出ブロック００９へ（ステップＳ４０７）と出力される。
なお、通常アポダイゼーションを用いて合成開口する場合は、共通する２焦点を有する双曲線群、楕円群の交点をサンプリング点とせず、プローブ開口から直線的に進む直線上にサンプリング点を設定しても良い。
【００５１】
（変位の算出）
次に変位算出ブロック００９での動作について述べる。
変位算出ブロック００９には、異なる時点で取得された受信信号を用いて算出された合成開口処理結果である、複数の第二の合成開口信号が入力される。
【００５２】
ここでは複数の第二の合成開口信号、すなわち異なる時点での合成開口処理結果間から変位を算出する方法を説明する。
例として双曲線インデックス軸と楕円インデックス軸とにそれぞれ独立にＩＱ信号を求め変位を算出する方法に関して述べる。
双曲線インデックス軸上の位置をｕ、楕円インデックス軸上の位置をｖとした時、入力された第二の合成開口信号は以下の式（１５）のようにモデル化出来る。
【数１１】
Ａは第二の合成開口信号の包絡線、ｆｕ、ｆｖはそれぞれ双曲線インデックス軸、楕円インデックス軸方向の変調周波数、φｕ、φｖはそれぞれ双曲線インデックス軸、楕円インデックス軸方向の位相を表す。
【００５３】
ここで、双曲線インデックス軸方向の波のＩＱ信号Ｉｕ（ｕ，ｖ）、Ｑｕ（ｕ，ｖ）、楕円インデックス軸方向の波のＩＱ信号Ｉｖ（ｕ，ｖ）、Ｑｖ（ｕ，ｖ）は、以下の式（１６）〜式（１９）で算出できる。
【数１２】
【００５４】
ここで、ＩＩ（ｕ，ｖ），ＩＱ（ｕ，ｖ），ＱＩ（ｕ，ｖ），ＱＱ（ｕ，ｖ）は、以下の式（２０）〜式（２３）で表される。
【数１３】
ＬＰＦはカットオフ周波数をｆｕおよびｆｖとするローパスフィルタを意味する。
【００５５】
このとき包絡線Ａ（ｕ，ｖ）と双曲線インデックス軸方向の位相φｕ、楕円インデックス軸方向の位相φｖは以下の式（２４）〜式（２６）で求められる。
【数１４】
【００５６】
異なるタイミングで得られた第二の合成開口信号それぞれの位相をφｕ１、φｕ２、φｖ１、φｖ２とする。このとき、双曲線インデックス軸方向の変位δｕ、楕円インデックス軸方向の変位δｖという２方向での変位は、以下の式（２７）、式（２８）で算出できる。
【数１５】
ここでλｕ、λｖはそれぞれ双曲線インデックス軸方向、楕円インデックス軸方向での波長である。
【００５７】
このようにして必要な観察範囲内の変位分布を算出可能である。このように第二の合成開口信号のＩＱデータを用いて変位を算出できるのは、サンプリングされた点が双曲線方向、楕円方向共に各々一定の位相差で配列されている効果であり、計算規模の低減に寄与する。また、符号列信号を用いたＳＮ比の高い信号に対して変位算出処理を行うため、より高い変位推定精度を実現することが出来る。
【００５８】
（変位算出の変形例１）
上記の例では、双曲線インデックス軸方向、楕円インデックス軸方向に分離したＩＱ信号を用いる手法について述べたが、複数の波形信号から位相情報や変位情報を抽出する手法を用いても構わない。
【００５９】
例えば、以下のような手法でも変位を算出することが出来る。
異なる時間で取得した第二の合成開口信号の包絡線を求め、それらの包絡線の相互相関によって大まかな変位（δｕ’，δｖ’）を推定する。大まかな変位分を補正した後、さらに複素相関によって位相差（φｕ’，φｖ’）を算出する。さらに別途注目位置での瞬時周波数（ｆｕ，ｆｖ）を推定する。
【００６０】
双曲線インデックス軸方向の変位δｕ、楕円インデックス軸方向の変位δｖは以下の式（２９）、式（３０）で算出できる。ここでｃは音速である。
【数１６】
【００６１】
この手法では瞬時周波数を推定し、その推定結果を用いることにより、観察領域内の双曲線群や楕円群のずれによる周波数の変化に対しても安定して変位を求めることが可能である。さらにその瞬時周波数を推定することによって精度向上の効果も得られる。
なお、この手法では観測個所ごとに瞬時周波数を推定するため双曲線群と楕円群との交点でサンプルされた信号を用いなくとも、横方向の変位が測定可能である。その場合であっても、符号化送受信によるＳＮ比向上の効果を得ることが可能である。
【００６２】
（変位算出の変形例２）
また、さらに以下のような手法でも変位を算出することが出来る。
異なる時間で取得した第二の合成開口信号の包絡線を求め、それらの包絡線の相互相関によって大まかな変位（δｕ’，δｖ’）を推定する。大まかな変位を補正した第二の合成開口信号をそれぞれｓ１（ｕ，ｖ）、ｓ２（ｕ，ｖ）として、任意のずれ量Δｕ、Δｖを与えて複素相関の偏角Ｒ（Δｕ，Δｖ：ｕ，ｖ）を算出する。これを式（３１）に示す。
【数１７】
【００６３】
この複素相関の偏角Ｒ（Δｕ，Δｖ：ｕ，ｖ）は、変位（δｕ’＋Δｕ，δｖ’＋Δｖ）が真の変位と一致した時に０になる。局所的には偏角が一次関数で変化すると近似して、偏角Ｒ（Δｕ，Δｖ：ｕ，ｖ）が０付近の値を取る複数の（Δｕ，Δｖ）で線形補間を行い、確からしい変位を算出する。
【００６４】
この手法では変位そのものを推定するため、観察領域内の双曲線群や楕円群のずれによる周波数の変化に対しても安定して変位を求めることが可能であり、変位推定の精度向上の効果が期待できる。なお、この手法では観測個所ごとに変位そのものを推定するため、双曲線群と楕円群との交点でサンプルされた信号を用いなくとも、横方向の変位が測定可能であり、その場合であっても符号化送受信によるＳＮ比向上の効果を得ることが可能である。
【００６５】
（歪みの算出）
次にこの変位データ（δｕ，δｖ）を用いて歪み分布を求める手法について述べる。
変位データを楕円−双曲線平面からｘ−ｙ平面に座標変換し、変位ベクトルとして算出する。この変位ベクトルをそれぞれｘ方向、ｙ方向に微分することで歪み分布を算出することが出来る。なお、楕円−双曲線平面内での歪み分布を先に算出し、その後ｘ−ｙ平面に座標変換しても同様の結果が得られる。
【００６６】
なお、歪み分布を算出する際には微分処理を実施するため空間周波数の高いノイズが発
生する可能性がある。このノイズを低減するために、画像処理ブロック０１０において、算出した歪み分布に対して空間的なローパスフィルタを適用するなどの視認性向上を目的とした処理を行っても良い。
【００６７】
ここまでに述べたように、本発明によれば符号列信号を用いることで合成開口法に適した１素子毎の送信を実施した場合と同等の信号を、よりＳＮ比を高く得ることができる。さらに、そのＳＮ比の高い信号を用いて変位推定するため、より変位推定精度の高い被検体情報処理装置を得ることが出来る。
【００６８】
図７（ａ）は本発明を用いて推定した歪み分布を示す図である。寒天内に周辺よりも硬さの硬い直径９ｍｍの内包物（図中点線で囲まれた領域）を設置し、図示しない探触子（図面上側に接触している）に対して垂直な方向（図面右方向）から０．５％程度の圧縮を加える。圧縮を加える前後での変位を本発明によって算出し、横方向の歪み分布を算出する。図７（ａ）はこの歪み分布をグレイスケールによって示している。
【００６９】
図７（ｂ）は図７（ａ）中の２箇所（２直線上。ｙ＝２９．５ｍｍ、３８．０ｍｍ）における横方向の歪み分布をプロットしたものである。ｙ＝２９．５ｍｍと記された直線上の横方向の歪み分布は実線で、ｙ＝３８．０ｍｍと記された直線上の横方向の歪み分布は点線で示している。点線で示したｙ＝３８．０ｍｍの直線上における横方向の歪み分布はほぼ−０．５％で一定しており、横方向から０．５％程度で圧縮されている状態を計測出来ている。また、実線で示したｙ＝２９．５ｍｍの直線上における横方向の歪み分布においては、中央付近に歪みの小さい領域（−０．２％程度）が存在しており、内包物が周辺の寒天よりも硬いことが計測出来ていることが分かる。
【００７０】
＜実施例２＞
符号列信号に対して位相を変化させた電圧波形を出力する際には、必ずしも一つの符号に対して電圧波形は１波長分である必要はない。すなわち、位相を変化させた電圧波形が、複数の波長分の長さを有していても良い。
【００７１】
このように複数の波長分の長さを有している電圧波形を用いた場合、１波長分の電圧波形を用いた場合に比べて、さらにＳＮ比の向上ならびに変位推定の精度向上の効果が得られる。
なお、電圧波形の長さが長くなるに従って、変位推定の空間分解能が低下するため、先述した空間的なローパスフィルタによって決定される空間分解能以下もしくは同等程度までに電圧波形の長さを抑えることが望ましい。
【００７２】
＜実施例３＞
符号列信号のそれぞれの値を１素子ずつに適用するのではなく複数の素子に対して適用することも可能である。
【００７３】
図６を用いて符号列信号のそれぞれの値を複数の素子に対して適用した例を示す。符号列信号は先述した例と同様のものを用いる。
このとき、送受信器の素子は図示したようにグループ化されている。素子６０１から素子６０３は素子グループＡ（６２１）、素子６０４から素子６０６は素子グループＢ（６２２）、素子６０７から素子６０９は素子グループＣ（６２３）、素子６１０から素子６１２は素子グループＤ（６２４）である。
【００７４】
符号制御ブロック００２から送信された符号列信号は送信回路系００５に入力される。送信回路系００５では１２個の素子６０１から素子６１２へ電圧波形の送信を行う。
まず、符号列信号の１列目に並んだ（１，１，１，１）を用いて、素子グループＡ、素
子グループＢ、素子グループＣ、素子グループＤに同位相（例えばこの位相を０度の位相とする）の電圧波形を送信する。
次に符号列信号の２列目に並んだ（１，−１，１，−１）に従って、素子グループＡに位相０度、素子グループＢに位相１８０度、素子グループＣに位相０度、素子グループＤに位相１８０度の電圧波形を送信する。
【００７５】
このように、符号列信号のそれぞれの値を複数の素子（ここでは３素子）ごとに適用し、超音波を被検体に送信する。さらに被検体内部で反射された超音波を受信し、復号化を行う。復号化した結果は先述した１素子毎の送信を行った場合と同等の信号とは異なり、素子グループＡ、素子グループＢ、素子グループＣ、素子グループＤでそれぞれ送信した信号が分離される。
これらの復号化された信号を用いて合成開口処理を行うステップ以降は先述した内容と同等のため説明を省略する。
【００７６】
このように複数の素子に対して同じ符号列信号を用いて送信を行うことで、送受信に用いる素子数よりも送受信回数を減らすことが可能であり、より高速にデータを取得することが出来る。そのため、より変化の早い変位や歪みを安定に精度良く測定することが可能になる。
【００７７】
なお、符号列信号のそれぞれの値が適用された複数の素子（例えば素子グループＡ）は、仮想的に波源を送受信器（探触子）の内部や被検体内に設定することで、より球面波もしくは円柱状の波面を形成することが可能である。これにより、複数の素子を同時に送信する場合と比較してさらに合成開口による空間分解能を向上させ、変位測定の精度を向上させることが可能である。
【００７８】
なお、上記実施形態１で説明した、変位計測の変形例１、２においては、上記の通り、変位計測の測定誤差を抑える、信号レベルを向上させるという点でより好ましい例である。ただし、これらの例において、計測の測定誤差を抑えるという点においては、必ずしも符号化パルスを用いなくても良い。つまり、送信パルスが符号化されたパルスではない場合においても、変位計測の測定誤差を抑えるという効果を奏するものである。
【符号の説明】
【００７９】
００１：システム制御部，００２：符号制御ブロック，００３：超音波変換素子，００７：復号合成ブロック，００８：画像生成ブロック，００９：変位算出ブロック
【特許請求の範囲】
【請求項１】
被検体に弾性波を送信し、該送信された弾性波の被検体での反射波を受信して該被検体の情報を取得する被検体情報取得装置であって、
弾性波と電気信号とを変換可能な素子を少なくとも１方向に複数配列して備える送受信器と、
前記素子に電気信号を入力し、前記素子から前記被検体に弾性波を送信させる素子制御手段と、
前記素子が受信する前記送信された弾性波の前記被検体での反射波を検出する検出手段と、を有し、
前記素子制御手段が前記素子に入力する電気信号は、前記複数の素子間において符号化された符号化パルス信号であり、
前記検出手段は、前記反射波を復号化するとともに、複数の素子が配列された前記１方向における２箇所を焦点とした双曲線に沿った方向と該２箇所を焦点とする楕円に沿った方向とを軸とし、該両軸の交点箇所または該交点に対応する箇所について前記復号化された反射波を合成する合成開口処理を異なる時点において行い、該異なる時点での合成開口処理結果に基づき、前記被検体の少なくとも２方向における変位を取得する
ことを特徴とする被検体情報取得装置。
【請求項２】
前記変位の取得が、前記交点箇所または交点に対応した箇所における前記合成開口処理結果の瞬時周波数の推定結果に基づいて行われる
ことを特徴とする請求項１に記載の被検体情報取得装置。
【請求項３】
前記変位の取得が、前記交点箇所または交点に対応した箇所における前記異なる時点での合成開口処理結果間の複素相関の偏角に基づいて行われる
ことを特徴とする請求項１に記載の被検体情報取得装置。
【請求項４】
前記符号化パルス信号が、複数の素子毎に符号化されている
ことを特徴とする請求項１に記載の被検体情報取得装置。
【請求項５】
被検体に弾性波を送信し、該送信された弾性波の被検体での反射波を受信して該被検体の情報を取得する被検体情報取得装置であって、
弾性波と電気信号とを変換可能な素子を複数配列して備える送受信器と、
前記素子に電気信号を入力し、前記素子から前記被検体に弾性波を送信させる素子制御手段と、
前記素子が受信する前記送信された弾性波の前記被検体での反射波を検出する検出手段と、を有し、
前記素子制御手段が前記素子に入力する電気信号は、前記複数の素子間において符号化された符号化パルス信号であり、
前記検出手段は、前記反射波を復号化するとともに、前記復号された反射波を合成する合成開口処理を異なる時点において行い、該異なる時点での合成開口処理結果に基づき、前記被検体の少なくとも２方向における変位を取得する
ことを特徴とする被検体情報取得装置。
【請求項１】
被検体に弾性波を送信し、該送信された弾性波の被検体での反射波を受信して該被検体の情報を取得する被検体情報取得装置であって、
弾性波と電気信号とを変換可能な素子を少なくとも１方向に複数配列して備える送受信器と、
前記素子に電気信号を入力し、前記素子から前記被検体に弾性波を送信させる素子制御手段と、
前記素子が受信する前記送信された弾性波の前記被検体での反射波を検出する検出手段と、を有し、
前記素子制御手段が前記素子に入力する電気信号は、前記複数の素子間において符号化された符号化パルス信号であり、
前記検出手段は、前記反射波を復号化するとともに、複数の素子が配列された前記１方向における２箇所を焦点とした双曲線に沿った方向と該２箇所を焦点とする楕円に沿った方向とを軸とし、該両軸の交点箇所または該交点に対応する箇所について前記復号化された反射波を合成する合成開口処理を異なる時点において行い、該異なる時点での合成開口処理結果に基づき、前記被検体の少なくとも２方向における変位を取得する
ことを特徴とする被検体情報取得装置。
【請求項２】
前記変位の取得が、前記交点箇所または交点に対応した箇所における前記合成開口処理結果の瞬時周波数の推定結果に基づいて行われる
ことを特徴とする請求項１に記載の被検体情報取得装置。
【請求項３】
前記変位の取得が、前記交点箇所または交点に対応した箇所における前記異なる時点での合成開口処理結果間の複素相関の偏角に基づいて行われる
ことを特徴とする請求項１に記載の被検体情報取得装置。
【請求項４】
前記符号化パルス信号が、複数の素子毎に符号化されている
ことを特徴とする請求項１に記載の被検体情報取得装置。
【請求項５】
被検体に弾性波を送信し、該送信された弾性波の被検体での反射波を受信して該被検体の情報を取得する被検体情報取得装置であって、
弾性波と電気信号とを変換可能な素子を複数配列して備える送受信器と、
前記素子に電気信号を入力し、前記素子から前記被検体に弾性波を送信させる素子制御手段と、
前記素子が受信する前記送信された弾性波の前記被検体での反射波を検出する検出手段と、を有し、
前記素子制御手段が前記素子に入力する電気信号は、前記複数の素子間において符号化された符号化パルス信号であり、
前記検出手段は、前記反射波を復号化するとともに、前記復号された反射波を合成する合成開口処理を異なる時点において行い、該異なる時点での合成開口処理結果に基づき、前記被検体の少なくとも２方向における変位を取得する
ことを特徴とする被検体情報取得装置。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【公開番号】特開２０１３−９９３７０（Ｐ２０１３−９９３７０Ａ）
【公開日】平成２５年５月２３日（２０１３．５．２３）
【国際特許分類】
【出願番号】特願２０１１−２４３４８３（Ｐ２０１１−２４３４８３）
【出願日】平成２３年１１月７日（２０１１．１１．７）
【出願人】（０００００１００７）キヤノン株式会社 (59,756)
【Ｆターム（参考）】
【公開日】平成２５年５月２３日（２０１３．５．２３）
【国際特許分類】
【出願日】平成２３年１１月７日（２０１１．１１．７）
【出願人】（０００００１００７）キヤノン株式会社 (59,756)
【Ｆターム（参考）】
[ Back to top ]