弾性画像生成方法及び超音波診断装置
【課題】検者の経験や熟練度に依存することなく、客観的で確定的な組織鑑別を可能とする弾性画像を提供する。
【解決手段】被検体の生体組織に圧迫を加えて超音波を送信し（２）、被検体から発生する反射エコー信号を計測して得られたフレームデータを用いて、取得時刻が異なる一対のフレームデータに基づいて複数の計測点における生体組織の弾性情報と各計測点に加わる圧迫の程度に関する計測圧迫条件を求め（１３）、求めた計測圧迫条件が予め設定された基準圧迫条件を満たすときの弾性情報に基づいて弾性画像を生成して（１８）、表示画面（１０）に表示することにより、検者の経験や熟練度に依存することなく、客観的で確定的な組織鑑別を可能とする弾性画像を提供する。
【解決手段】被検体の生体組織に圧迫を加えて超音波を送信し（２）、被検体から発生する反射エコー信号を計測して得られたフレームデータを用いて、取得時刻が異なる一対のフレームデータに基づいて複数の計測点における生体組織の弾性情報と各計測点に加わる圧迫の程度に関する計測圧迫条件を求め（１３）、求めた計測圧迫条件が予め設定された基準圧迫条件を満たすときの弾性情報に基づいて弾性画像を生成して（１８）、表示画面（１０）に表示することにより、検者の経験や熟練度に依存することなく、客観的で確定的な組織鑑別を可能とする弾性画像を提供する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、弾性画像生成方法及び超音波診断装置に係り、特に、被検体の生体組織に圧迫を加えて計測される組織弾性を画像化して診断に供する弾性画像生成方法及び超音波診断装置に関する。
【背景技術】
【０００２】
超音波診断装置は、被検体の表面に超音波探触子を当て、その探触子から被検体に超音波を送信し、被検体内部からの超音波の反射波を受信し、その受信信号である反射エコー信号に基づいて被検体の各部の生体情報の断層像などの画像により表示して診断に供するものである。
【０００３】
一方、特許文献１、２等に記載されているように、被検体に圧迫力を加えて内部に生じた組織の歪みや弾性率等の弾性に関する情報（以下、弾性情報という。）を求めて弾性画像を生成し、その弾性画像を表示して検者の病変部等の診断に供することが行われている。例えば、検者は、圧迫力を加減しながら計測される弾性画像を観察することにより、正常組織、ガン細胞、腫瘍などの病変を鑑別する。
【０００４】
ところで、生体組織は、圧迫に対して非線形の弾性応答を示し、圧縮量などの圧迫の程度（以下、圧迫条件という。）に応じて硬さが変化することが知られている。例えば、非特許文献１によれば、図２に示すような乳腺組織の乳腺部位と脂肪部位を圧縮した際の応力（ｋＰａ）と歪み（％）の関係の実測値が報告されている。同図によれば、脂肪部位の応力―歪みは直線的な関係であるが、乳腺部位の応力―歪みの関係は非線形である。したがって、図２の応力−歪み曲線の傾きで与えられる弾性率（ヤング率）は、図３に模式的に示すように、脂肪部位は一定値であるが、乳腺部位は加えられた歪みに応じて変化することが分かる。
【０００５】
一方、非特許文献２には、歪みと弾性率の関係のデータを取得した後、歪み−弾性率曲線を関数で近似し、最小二乗法などにより最もよく近似する曲線からその非線形性を表す非線形パラメタを抽出し、その組織の非線形性を表す情報として評価する試みが提案されている。
【０００６】
このように、組織の硬さの情報を用いて組織の良悪性を鑑別するにあたって、例えば弾性率がどのような圧迫条件の下で計測されたかが非常に重要な情報であり、圧縮量などの圧迫条件が大きく異なる条件で計測された診断対象については、良悪性を誤って診断してしまうことも考えられる。
【０００７】
【特許文献１】特開平５−３１７３１３号号公報
【特許文献２】特開２０００−６０８５３号公報
【非特許文献１】Krouskop T, et al : Elastic Moduli of Breast and Prostate Tissues Under Compression. Ultrasonic Imaging 20 : 260-274, 1998.
【非特許文献２】「超音波による組織の非線形弾性特性の画像化」新田他、電子情報通信学会論文誌 2001/12 Vol.J84-A No.1
【発明の開示】
【発明が解決しようとする課題】
【０００８】
しかしながら、特許文献１、２等に記載された従来の技術では、圧迫の程度に応じて硬さが変化する生体組織の非線形性を考慮していない。そのため、検者は様々な圧迫条件下において複数の弾性画像を取得し、その中から検者が主観的に適切と考える１又は複数の弾性画像に基づいて診断することを余儀なくされている。その結果、従来技術によれば、適切な弾性画像を選択するために経験や熟練を要し、診断の信頼性及び客観性に問題が残る。また、圧迫条件を変化させながら適切と思われる弾性画像を繰返し取得しなければならず、検者に負担がかかるとともに、検査効率が悪いという問題がある。
【０００９】
本発明は、超音波診断装置において、検者の経験や熟練度に依存することなく、客観的で確定的な組織鑑別を可能とする弾性画像を提供することを課題とする。
【課題を解決するための手段】
【００１０】
上記の課題を解決するため、本発明は、被検体の生体組織に圧迫を加えて超音波を送信し、前記被検体から発生する反射エコー信号を計測して得られたフレームデータを用いて、取得時刻が異なる一対の前記フレームデータに基づいて複数の計測点における生体組織の弾性情報と前記各計測点に加わる圧迫の程度に関する計測圧迫条件を求め、求めた前記計測圧迫条件が予め設定された基準圧迫条件を満たすときの前記弾性情報に基づいて弾性画像を生成して表示画面に表示することを特徴とする。
【００１１】
このように、各計測点に加わる計測圧迫条件が予め設定された基準圧迫条件を満たして計測及び生成された弾性画像を用いているから、検者の主観を廃することができ、検者の経験や熟練度に依存することなく、客観的で確定的な組織鑑別を可能とすることができ、かつ検査の効率を向上できる。また、診断に適した弾性画像の選定作業が不要となるから検者の負担が軽減される。
【００１２】
本発明の超音波診断装置は、被検体の生体組織に圧迫を加え超音波を送信して得られる反射エコー信号のフレームデータを取得するフレームデータ取得手段と、該フレームデータ取得手段により取得された取得時刻が異なる一対の前記フレームデータに基づいて複数の計測点における生体組織の弾性情報と前記各計測点に加わる圧迫の程度を示す計測圧迫条件を求める演算手段と、該演算手段により求められた前記弾性情報に基づいて弾性画像を生成する弾性画像生成手段と、前記演算手段により求められた前記計測圧迫条件が予め設定された基準圧迫条件を満たすか否か評価し、前記基準圧迫条件を満たす前記弾性画像を表示手段に表示させる圧迫条件評価手段とを有して構成することができる。
【００１３】
本発明の圧迫条件には、計測点の歪み、計測点の変位量、計測点の応力、被検体の体表に加えられた圧力のいずれか一つを採用することができる。なお、本発明でいう歪み及び変位量は、圧迫による生体組織の歪み変化又は変位のゼロから計測時点までの積分量をいい、その時々の歪み変化又は変位と区別する。また、基準圧迫条件は、生体組織の弾性を評価するのに適した圧迫条件を設定する。
【００１４】
また、弾性情報には、弾性率、粘弾性率、弾性率の非線形性に係る非線形パラメタを採用することができる。
【００１５】
また、フレームデータに対応する画像領域内に設定された関心領域内の計測点について、計測圧迫条件と計測された弾性情報との関係を線図（グラフ）にして表示画面に表示することができる。この場合、線図の基準圧迫条件の位置にマークを表示することが好ましい。
【００１６】
また、弾性画像に関心領域の輪郭を表示するとともに、関心領域に対応付けて関心領域内の複数の計測点の弾性情報の平均値を数値で表示することができる。弾性画像は、各画素の色相又は輝度を各弾性情報に応じて割り付けて生成される。
【発明の効果】
【００１７】
本発明によれば、検者の経験や熟練度に依存することなく、客観的で確定的な組織鑑別を可能とする弾性画像を提供することができる。
【発明を実施するための最良の形態】
【００１８】
以下、本発明の弾性画像生成方法及び装置を実施形態に基づいて説明する。図１に、本発明の弾性画像生成方法を実施する一実施形態の超音波診断装置のブロック構成図を示す。図に示すように、被検体１に当接して用いられる超音波の探触子２は、被検体１との間で超音波を送信及び受信する複数の振動子を有して形成されている。探触子２は、送信回路３から供給される超音波パルスにより駆動される。送受信制御回路４は、探触子２の複数の振動子を駆動する超音波パルスの送信タイミングを制御して、被検体１内に設定される焦点に向けて超音波ビームを形成するようになっている。また、送受信制御回路４は、探触子２の振動子の配列方向に電子的に超音波ビームを走査するようになっている。
【００１９】
一方、探触子２は、被検体１内から発生する反射エコー信号を受信して受信回路５に出力する。受信回路５は、送受信制御回路４から入力されるタイミング信号に従って、反射エコー信号を取り込んで増幅などの受信処理を行う。受信回路５により処理された反射エコー信号は、整相加算回路６において複数の振動子により受信された反射エコー信号の位相を合わせて加算することにより増幅される。整相加算回路６において整相加算された反射エコー信号は、信号処理部７に入力され、ゲイン補正、ログ圧縮、検波、輪郭強調、フィルタ処理等の信号処理がなされる。
【００２０】
信号処理部７により処理された反射エコー信号は、白黒スキャンコンバータ８に導かれて超音波ビームの走査面に対応した２次元の断層像データ（ディジタルデータ）に変換される。これらの信号処理部７と白黒スキャンコンバータ８によって断層像（Ｂモード像）の画像再構成手段が構成される。白黒スキャンコンバータ８から出力される断層像データは、切替加算部９を介して画像表示器１０に供給されてＢモード像が表示されるようになっている。
【００２１】
一方、整相加算回路６から出力される反射エコー信号は、フレームデータ取得部１１に導かれる。フレームデータ取得部１１は、超音波ビームの走査面（断層面）に対応する反射エコー信号群を、フレームデータとして複数フレーム分を取得してメモリなどに格納する。変位計測部１２は、フレームデータ取得部１１に格納されている取得時刻が異なる複数対のフレームデータを順次取り込み、取り込んだ一対のフレームデータに基づいて断層面における複数の計測点の変位ベクトルを求め、変位フレームデータとして弾性情報演算部１３に出力するようになっている。
【００２２】
本実施例の弾性情報演算部１３は、変位フレームデータに基づいて各計測点の生体組織の歪み変化を求める歪み演算部と、歪みフレームデータに基づいて各計測点の生体組織の弾性率を求める弾性率演算部と、歪み演算部で求めた歪み変化を積分して計測圧迫条件の一例である歪みを求める歪み演算部を有して構成されている。弾性情報演算部１３にて求められた弾性率のフレームデータは弾性情報処理部１４に出力され、歪みは圧迫条件評価部１８に出力されるようになっている。
【００２３】
弾性情報処理部１４は、弾性情報演算部１３から入力される各弾性情報のフレームデータに対して、座標平面内におけるスムージング処理、コントラスト最適化処理、フレーム間における時間軸方向のスムージング処理などの様々な画像処理を施して、カラースキャンコンバータ１５に送出するようになっている。
【００２４】
カラースキャンコンバータ１５は、弾性情報処理部１４により処理された弾性率のフレームデータを取り込み、設定された弾性率のカラーマップに従って、フレームデータの画素ごとに色調コードを付与してカラー弾性画像を生成するようになっている。
【００２５】
カラースキャンコンバータ１５により生成されたカラー弾性画像は、切替加算部９を介して画像表示器１０に表示されるようになっている。また、切替加算部９は、白黒スキャンコンバータ８から出力される白黒の断層像と、カラースキャンコンバータ１５から出力されるカラー弾性画像とを入力し、両画像を切り替えていずれか一方を表示させる機能と、両画像の一方を半透明にして加算合成して画像表示器１０に重ねて表示させる機能と、両画像を並べて表示させる機能を有して形成されている。また、切替加算部９から出力される画像データは、装置制御インターフェイス部１９の制御に従ってシネメモリ２０に格納されるようになっている。シネメモリ２０に格納された画像データは、装置制御インターフェイス部１９の制御に従って画像表示器１０に表示されるようになっている。
【００２６】
本実施形態の特徴に係る圧迫条件評価部１８は、弾性情報演算部１３により求められた計測圧迫条件である歪みを取り込み、その歪みが装置制御インターフェイス部１９から入力される基準圧迫条件を満たすか否か判断するようになっている。その判断の結果、取り込んだ歪みが基準圧迫条件を満たす場合は、切替加算部９に指令を送って、カラースキャンコンバータ１５から出力されるカラー弾性画像を画像表示器１０に出力表示させるようになっている。
【００２７】
このように構成される本実施形態の基本的な動作について説明する。まず、探触子２により被検体１に圧迫を加えて被検体１に超音波ビームを走査し、走査面からの反射エコー信号を連続的に受信する。そして、整相加算回路６から出力される反射エコー信号に基づいて、信号処理部７及び白黒スキャンコンバータ８により断層像が再構成され、切替加算器９を介して画像表示器１０に表示される。
【００２８】
一方、フレームデータ取得部１１は、反射エコー信号を取り込んでフレームレートに同期させてフレームデータを繰り返し取得し、内蔵されたフレームメモリ内に時系列順に保存する。そして、取得時刻が異なる一対のフレームデータを単位として、連続的に複数対のフレームデータを選択して変位計測部１２に出力する。変位計測部１２は、選択された一対のフレームデータを1次元もしくは２次元相関処理し、走査面における複数の計測点の変位を計測して変位フレームデータを生成する。この変位ベクトルの検出法としては、例えば特開平５−３１７３１３号公報等に記載されているブロックマッチング法又はグラジェント法が知られている。ブロックマッチング法は、画像を例えばＮ×Ｎ画素からなるブロックに分け、現フレーム中の着目しているブロックに最も近似しているブロックを前フレームから探索し、これに基づいて計測点の変位を求める。また、一対のＲＦ信号フレームデータの同一領域における自己相関を計算して変位を算出することができる。
【００２９】
弾性情報演算部１３は、変位フレームデータを取り込んで、各計測点の歪み変化を求め、求めた歪み変化に基づいて弾性情報である弾性率を演算し、弾性率フレームデータを弾性情報処理部１４に出力する。歪み変化の演算は、公知のように変位を空間微分することによって計算される。また、求めた歪み変化に基づいて各計測点の弾性率を演算する。弾性率を求める場合は、圧力計測部１７により計測された圧力の計測値を取り込み、これに基づいて各計測点における応力を演算する。圧力計側部１７は、探触子2の超音波送受信面と被検体1との間に設けられた圧力センサ１６により検出された圧力に基づいて、被検体１内部の計測点における応力を演算する。つまり、弾性情報演算部１３の弾性率演算部は、各計測点における応力と、弾性情報演算部１３で求めた歪みフレームデータから走査面上の各計測点の弾性率Ｅ（例えば、ヤング率）を演算し、弾性情報処理部１４に出力する。
【００３０】
弾性情報処理部１４は、入力される弾性率をスムージング処理などの処理を施してカラースキャンコンバータ１５に出力する。カラースキャンコンバータ１５は、弾性情報に基づいてカラー弾性画像を生成する。カラー画像は、例えば２５６階調化による色調のグラデーションで画素単位ごとにフレームデータの弾性率に応じて色付けされる。なお、カラースキャンコンバータ１５に代えて、白黒スキャンコンバータを用いることができる。この場合は、弾性率が大きい領域は輝度を明るく、逆に小さい領域は輝度を暗くするなどにより、良性又は悪性を鑑別できるようにすることができる。
【００３１】
以下に、本実施形態を用いて本発明の弾性画像生成方法を実施する具体的な実施例を説明する。なお、各実施例は、本実施形態の特徴部である弾性情報演算部１３、圧迫条件評価部１８、装置制御インターフェイス部１９及びシネメモリ２０等により実施される。
【実施例１】
【００３２】
本発明の弾性画像生成方法の一実施例として、図４に示す生体組織１〜５を含む部位に圧迫を加えて弾性画像を生成する場合を説明する。本実施例では、弾性情報演算部１３において弾性フレームデータとして生体組織１〜５の各部位における計測点の弾性率を弾性情報として求めるものとする。なお、図４は、弾性フレームデータと弾性画像との位置関係を対応付けるために、１フレームの画像の体裁で示している。
【００３３】
まず、圧迫条件評価部１８には、装置制御インターフェイス部１９から検者によって関心領域（以下、ＲＯＩという。）としてＲＯＩ１、ＲＯＩ２が入力設定される。圧迫条件評価部１８は、弾性情報演算部１３からＲＯＩ１、ＲＯＩ２の中に分布する各計測点の弾性率データ（弾性率の数値データ）と、ＲＯＩ１、ＲＯＩ２の中に分布する各計測点の歪みの数値データを取り込む。本実施例では、圧迫の程度を評価するための圧迫条件として歪みを用いている。また、歪み（ε＝ΣΔε）は、圧縮ゼロから現時点までの歪み変化の積分値であり、一対のフレームデータ間ごとに得られる歪み変化（Δε）とは区別する。
【００３４】
圧迫条件評価部１８は、取り込んだ歪みと弾性率に基づいて、ＲＯＩ１、ＲＯＩ２のそれぞれについて、例えば、歪みε（％）の平均値と弾性率Ｅ（ｋＰａ）の平均値を演算し、それらの関係を図５に示すようにグラフにプロットする。そして、歪み（圧迫条件）−弾性情報（弾性率）曲線のグラフ画像を構築し、圧迫条件評価部１８に記憶するとともに、切替加算部９に送出して画像表示器１０に表示する。これにより、画像表示器１０には、ＲＯＩ１、ＲＯＩ２内の組織に関する圧迫条件−弾性情報曲線がグラフ画像として表示される。
【００３５】
例えば、基準圧迫条件を歪み４〜６％に設定すると、図５に示すように、基準圧迫条件の範囲が圧迫条件−弾性情報曲線に表示される。一方、圧迫条件評価部１８は、基準圧迫条件を満たした時点において計測されたＲＯＩ１、ＲＯＩ２の弾性率Ｅ１、Ｅ２を解析する。そして、図４に示したＲＯＩ１、ＲＯＩ２に対応する組織の弾性診断結果として、図５に示すように弾性率Ｅ１、Ｅ２の数値データをグラフ画像に表示して、ＲＯＩ１、ＲＯＩ２の基準圧迫条件を満たすときの弾性率を弾性診断結果として確定する。なお、一般に、ＲＯＩ１とＲＯＩ２の計測圧迫条件は異なることがあるから、基準圧迫条件が同じであっても、同じ計測時刻においてＲＯＩ１とＲＯＩ２が基準圧迫条件を満たすというわけではない。
【００３６】
ここで、基準圧迫条件は、範囲として設定する方法に限らない。例えば、歪み５％を基準圧迫条件として設定し、この歪み５％を跨ぐ前後において計測された弾性率の値から線形補間することにより歪み５％における弾性率の値を求めるようにすることができる。
【００３７】
また、例えば、圧迫操作を繰り返して行う弾性情報を計測する過程で、基準圧迫条件を満たして計測された弾性率の値により、基準圧迫条件を満たすたびに図５の弾性率を更新するようにすることができる。また、基準圧迫条件を満たして計測された現在までの弾性率の値との加算平均値を表示させるようにすることができる。
【００３８】
また、本実施例では、弾性情報としてヤング率に代表される弾性率を例に示したが、粘弾性率、非線形性などの組織の弾性特性に基づいた指標であれば用いることができる。
【００３９】
また、計測圧迫条件は、ＲＯＩ内の歪みに限らず、変位計測部１２で計測された変位量（変位ゼロからの積算を含む）、圧力計測部１２で計測された圧迫の圧力に基づいて導出されたＲＯＩ内に伝わった応力など、圧迫の状態を反映した情報を用いることができる。
【００４０】
上述したように、本実施例によれば、生体組織の画像にＲＯＩを設定して弾性情報の計測操作を実施するだけで、そのＲＯＩに対応する組織に加えられる計測圧迫条件が評価に適切な基準圧迫条件を満たしているときに取得された弾性情報である弾性率が数値により表示される。したがって、経験や熟練がない検者であっても、ＲＯＩに係る生体組織の適切な弾性情報を得ることができるから、信頼性及び客観性を有する診断を行うことができる。また、圧迫条件を変化させながら適切と思われる弾性画像を繰返し取得する操作が不要となり、検者の負担を軽減して、検査効率を向上できる。また、本実施例によれば、超音波によるリアルタイム性を損なうことなく、確定的な診断情報を臨床現場において効率よく取得することができる。
【実施例２】
【００４１】
実施例１では、ＲＯＩの弾性情報を数値データとして画像に出力する例を示した。しかし、本発明により生成する弾性画像は、これに限られるものではなく、図６に示す実施例２のように、基準圧迫条件を満たした時点において計測されたＲＯＩ内の弾性率に基づき、弾性率に応じてＲＯＩ内の画像を階調化した弾性画像を生成して表示することができる。
【００４２】
例えば、実施例１で設定したＲＯＩ内の計測圧迫条件が基準圧迫条件を満たしたときに、ＲＯＩ内において計測された弾性率の大きさに応じて、白黒輝度又は色相に階調を与えた画像を生成して表示する。このとき、弾性率の数値データとひずみの数値データも画像上に表示するようにしてもよい。
【実施例３】
【００４３】
実施例２では、各ＲＯＩについて、各計測点に対応する画素単位で弾性率を階調化して表示する例を説明した。しかし、図７の実施例３に示すように、実施例２のＲＯＩに代えて、計測領域の全体に画素を構築できる程度の小さな局所的ＲＯＩを格子状に複数配置することができる。この場合は、各局所的ＲＯＩ毎に独立に、計測圧迫条件が基準圧迫条件を満たしたときの各局所的ＲＯＩの弾性率の大きさを求め、局所的ＲＯＩごとに白黒輝度又は色相に階調を与えて、観測領域全域にわたって弾性画像を形成して表示する。なお、局所的ＲＯＩは、１又は複数の計測点を含む領域として設定することができる。
【００４４】
また、基準圧迫条件を満たして割当てられた各局所的ＲＯＩの弾性画像の階調は、繰り返し計測される弾性率の計測圧迫条件が基準圧迫条件を満たすたびに更新して表示されるようにすることができる。また、基準圧迫条件を満たして取得されたこれまでの弾性率の値の加算平均値に基づいて、弾性画像の階調を更新して表示すれば、弾性画像の精度を向上させることができる。
【００４５】
また、計測中に計測圧迫条件が基準圧迫条件を満たさなかったＲＯＩについては、弾性率の大きさを反映した階調が付与されないようにすれば、弾性評価不可能であったことを認識できる。
【実施例４】
【００４６】
実施例１〜３では、弾性情報として弾性率を用いる例を説明した。本実施例４では、図３に示した歪み−弾性率曲線の非線形性を表す非線形パラメタを弾性情報として画像化する例を説明する。つまり、計測圧迫条件が基準圧迫条件を満たした場合における各計測点の非線形パラメタの値を求め、求めた非線形パラメタに基づいて階調化した弾性画像を生成する。
【００４７】
前述したとおり、非特許文献１には、乳腺組織の弾性特性の非線形性に関して、脂肪組織は大きな歪みの範囲までほぼ一定の弾性率の線形な応答を示すのに対し、繊維組織や浸潤性腫瘍は歪みの増加に伴い弾性率も顕著に増加する歪み硬化が認められることが報告されている。また、繊維組織よりも浸潤性の悪性腫瘍の方が、歪み硬化の程度が大きく、非線形の程度も大きくなる。また、一般に生体組織は非線形の弾性応答を示し、例えば、歪み−弾性率の関係は、図３に示した曲線に従うことが知られている。非特許文献２には、この歪みと弾性率の関係のデータを取得した後、所定の関数でその曲線を近似し、最小二乗法などにより最もよく近似する曲線からその非線形パラメタを抽出し、その組織の非線形性を表す情報として評価する試みが提案されている。
【００４８】
つまり、非特許文献２によれば、Ｅをヤング率と歪みε（＝ΣΔε）の非線形性を解析する目的において、圧迫に対する組織の応答の曲線を、次式（１）のように１次関数として仮定している。
【００４９】
Ｅ＝Ｅ０＋α×ε （１）
ここで、Ｅ０は一定値、αは非線形パラメタであり、ヤング率−歪み曲線の傾きである。しかし、図８に示すように、大きな歪みになるほど、実際の組織の非線形性が顕著になり、非線形パラメタαは圧迫条件に大きく依存する。したがって、（１）式のようにＥがεの１次関数であるとする仮定する非線形パラメタαは、評価区間を定めなければ定量情報とは成り得ず、確定診断には採用できない。
【００５０】
そこで、本実施例では、実施例１〜３と同様、基準圧迫条件を設定し、計測圧迫条件が基準圧迫条件を満たして計測された弾性率と歪みとのデータに基づいて非線形パラメタαを求め、求めた非線形パラメタαを数値情報や画像情報として検者に提供するようにすることを特徴とする。具体的には、実施例１〜３における「弾性率」を「非線形パラメタ」に置き換えればよい。
【００５１】
さらに、非線形パラメタαを一層忠実に求めるためには、広範囲の計測圧迫条件において、組織の非線形応答を忠実に近似する曲線として、次式（２）に示す高次の関数により近似することが好ましい。γは、２以上の自然数である。
【００５２】
Ｅ＝Ｅ０＋α×εγ （２）
（２）式による非線形パラメタαも、組織の非線形性が大きいほど大きな値となり、線形に近づくほど限りなくゼロに近くなる。さらに、一般化すると、次式（３）の指数関数で近似できる。
【００５３】
Ｅ＝（Ｅ０−１）＋ ｅｘｐ（α×ε) （３）
これらの式（２）、（３）で示される関係式を仮定すれば、計測圧迫条件の評価区間を任意に設定しても、高精度に確定した非線形パラメタαを得ることができ、定量的な評価が実現できる。
【００５４】
例えば、従来は、図８のようにαの値が圧迫条件に依存するため、歪み−弾性率の関係のデータを蓄積した後、オフライン処理により計測圧迫条件について適当な評価区間を確認して解析する必要がある。この点、本実施例によれば、計測圧迫条件に依存しない確定的な非線形パラメタαが得られる。つまり、時刻ｔ（０）において計測された歪み−弾性率の関係の（ε（０）、Ｅ（０)）と、過去の時刻ｔ（−１）、ｔ（−２）、…において得られた歪み−弾性率の関係（ε（−１）、Ｅ（−１））、（ε（−２）、Ｙ（−２））、 …を歪み−弾性率の関係の蓄積データとして適用し、時刻ｔ（０）においてリアルタイムに非線形パラメタαを評価しても、その非線形パラメタαの数値情報、画像情報は圧迫条件に依存しない客観的な情報として得ることができる。
【００５５】
したがって、本実施例４によれば、図７で説明した場合と同様、非線形パラメタαを弾性情報とし、リアルタイムの弾性画像を生成して表示することが可能である。その結果、例えば、Ｂモード像や弾性率画像によって認識できない悪性組織であっても、非線形パラメタの評価によってのみ異常な組織を検出することができる場合があるので、このような場合に検査において異常組織領域を見逃してしまうことはなく、検者の主観に依らない診断を行うことができる。
【実施例５】
【００５６】
実施例１〜３においては、ＲＯＩに属する組織が圧迫操作により移動することについて考慮していない。本実施例５では、圧迫操作により移動する組織に追従させてＲＯＩを移動させて、弾性情報の評価の精度を向上させる。つまり、画像に対するＲＯＩの位置を固定すると、図９（ａ）〜（ｃ）に示すように、圧迫操作の経過に応じてＲＯＩ（ｔ０）、ＲＯＩ（ｔ１）、ＲＯＩ（ｔ２）内に入ってくる組織が変化してしまうことがあり、計測値の精度が悪くなる。
【００５７】
そこで、本実施例５では、図９（ａ’）〜（ｃ’）に示すように、圧迫の程度に応じて、ＲＯＩ内に属する組織の行き先を追従して、ＲＯＩ（ｔ０）、ＲＯＩ（ｔ１）、ＲＯＩ（ｔ２）を移動させることにより、任意の圧迫条件の下において、同一組織領域の情報のみを抽出することが可能となる。その結果、評価精度を向上させることが可能になる。なお、組織の動きに追従させてＲＯＩを移動させる処理は、前述した組織の変位ベクトルを求めるブロックマッチング法又はグラジェント法を適用することができる。また、局所的ＲＯＩや計測点についても、組織の移動に応じて追従させることができる。
【実施例６】
【００５８】
ここで、上記実施例１〜５等により生成された歪み−弾性率曲線のグラフ画像、基準圧迫条件を満たしたときの弾性率画像、非線形パラメタα画像等の弾性画像を、画像表示器１０に表示する際の表示形態の実施例について説明する。
【００５９】
切替加算部９と画像表示器１０の機能により、２画面を同時に表示可能にすれば、例えば、右画面に現時刻におけるＢモード像や弾性画像をリアルタイムで表示し、左画面に計測圧迫条件−弾性情報のグラフ、基準圧迫条件を満たして計測された弾性画像、もしくは非線形パラメタ画像を表示することができる。また、複数画面の内のいずれかに圧迫条件に依存しない定量の弾性情報を表示する画像を表示することができる。
【００６０】
このように、表示画像の組み合わせは、その他にもカラードプラなど従来の診断画像を含め、任意に選択できるようにすることができる。また、一つの画面に左右に並べて表示する方法に限らず、上下に並べて表示してもよい。さらに、一つの画面内に設けた小窓に別の画像を表示するような形式を採用することもできる。
【実施例７】
【００６１】
一般に、リアルタイム弾性画像の計測においては、超音波診断装置にフリーズ指令を与えると、フリーズまでに取得された時間的に連続した所定フレーム数の弾性画像がシネメモリ２０に確保される。そして、シネメモリ２０に確保された弾性画像を連続的に読み出して画像表示器１０に動画として表示することができるようになっている。
【００６２】
そこで、本実施例７では、圧迫操作を行いながら計測した複数の弾性画像をシネメモリ２０に格納するとともに、基準圧迫条件を満たした計測圧迫条件で得られた弾性画像のフレームに、基準圧迫条件を満たした弾性画像であることを示すマークを表示するようにする。例えば、基準圧迫条件を満たした弾性画像には、丸印（○）を画像のどこかに表示する。また、基準圧迫条件に対して計測圧迫条件が小さければ上向きの矢印（↑）を、大きければ下向きの矢印（↓）を、弾性画像に表示する。また、それらのマーク表示に代えて又はマーク表示と共に、音声によって検者に報知するようにしてもよい。
【００６３】
本実施例によれば、検者は、超音波診断装置からリアルタイムにフィードバックされるガイド情報により、計測圧迫条件が基準圧迫条件を満たす方向に調整することができるから、適切な弾性画像を取得することが可能となる。また、同時に、診断に適用する画像や弾性に関する数値情報を検者の主観によらずに確定的に決定することができる。
【実施例８】
【００６４】
実施例７では、本発明の弾性画像生成方法によりリアルタイムで生成された弾性画像を、画像表示器１０に連続的に表示する例を示したが、本発明はこれに限らず、計測圧迫条件が基準圧迫条件を満たしたときの弾性情報に対応する弾性画像のみを生成して画像表示器１０に表示するようにする。つまり、計測圧迫条件が基準圧迫条件を満たさない弾性画像は表示しないようにする。
【実施例９】
【００６５】
実施例７では、基準圧迫条件に対する計測圧迫条件の大小に応じて弾性画像にマークを表示するようにしたが、本発明はこれに限らず、計測圧迫条件が基準圧迫条件に近づくほど、カラー弾性画像の明るさが強くなるなど、弾性画像の見易さが増加するなどの工夫をすることができる。例えば、Ｂモード像に半透明のカラー弾性画像を重ね合わせた画像を生成する場合は、カラー弾性画像のブレンディング比率を計測圧迫条件が基準圧迫条件に近づくほど大きくするようにすることもできる。
【実施例１０】
【００６６】
また、計測圧迫条件が基準圧迫条件を満たしたときの弾性画像、関心部位の弾性率の数値情報、その他診断の判断材料となる情報は、弾性診断の結果として、超音波診断装置に備えられている各患者の診断レポート機能に自動的に反映するようにすることができる。
【００６７】
また、上記の各実施例においては、対象組織を圧迫することにより弾性情報を取得する例について説明したが、本発明の弾性画像生成方法は、例えば、加振映像法（Y. Yamakoshi et al : Ultrasonic Imaging of Internal Vibration of Soft Tissue under Forced Vibration, IEEE Trans UFFC 1990; 37; 45-53.）などの手法によって得られた弾性情報に対しても適用することができる。
【図面の簡単な説明】
【００６８】
【図１】本発明の弾性画像生成方法を実施する一実施形態の超音波診断装置のブロック構成図である。
【図２】乳腺組織の乳腺部位と脂肪部位を圧迫した際の歪みと応力の関係を示す線図である。
【図３】乳腺組織の乳腺部位と脂肪部位を圧迫した際の歪みと弾性率（ヤング率）の関係を模式的に示す線図である。
【図４】生体組織に圧迫を加えて弾性画像を生成する方法を説明する図である。
【図５】実施例１の弾性画像生成方法により生成される弾性画像の歪み―弾性率曲線のグラフ画像である。
【図６】実施例２の弾性画像生成方法により生成される弾性画像の一例である。
【図７】実施例３の弾性画像生成方法により生成される弾性画像の一例である。
【図８】実施例４の弾性画像生成方法により生成する非線形パラメタが圧迫条件に大きく依存することを説明する線図である。
【図９】圧迫操作によりＲＯＩ内に属する組織の移動に追従させてＲＯＩを移動させる実施例を説明する図である。
【符号の説明】
【００６９】
１ 被検体
２ 探触子
３ 送信回路
４ 送受信制御回路
５ 受信回路
６ 整相加算回路
７ 信号処理部
８ 白黒スキャンコンバータ
９ 切替加算部
１０ 画像表示器
１１ フレームデータ取得部
１２ 変位計測部
１３ 弾性情報演算部
１５ カラースキャンコンバータ
１６ 圧力センサ
１７ 圧力計測部
１８ 圧迫条件評価部
１９ 装置制御インターフェイス部
２０ シネメモリ
【技術分野】
【０００１】
本発明は、弾性画像生成方法及び超音波診断装置に係り、特に、被検体の生体組織に圧迫を加えて計測される組織弾性を画像化して診断に供する弾性画像生成方法及び超音波診断装置に関する。
【背景技術】
【０００２】
超音波診断装置は、被検体の表面に超音波探触子を当て、その探触子から被検体に超音波を送信し、被検体内部からの超音波の反射波を受信し、その受信信号である反射エコー信号に基づいて被検体の各部の生体情報の断層像などの画像により表示して診断に供するものである。
【０００３】
一方、特許文献１、２等に記載されているように、被検体に圧迫力を加えて内部に生じた組織の歪みや弾性率等の弾性に関する情報（以下、弾性情報という。）を求めて弾性画像を生成し、その弾性画像を表示して検者の病変部等の診断に供することが行われている。例えば、検者は、圧迫力を加減しながら計測される弾性画像を観察することにより、正常組織、ガン細胞、腫瘍などの病変を鑑別する。
【０００４】
ところで、生体組織は、圧迫に対して非線形の弾性応答を示し、圧縮量などの圧迫の程度（以下、圧迫条件という。）に応じて硬さが変化することが知られている。例えば、非特許文献１によれば、図２に示すような乳腺組織の乳腺部位と脂肪部位を圧縮した際の応力（ｋＰａ）と歪み（％）の関係の実測値が報告されている。同図によれば、脂肪部位の応力―歪みは直線的な関係であるが、乳腺部位の応力―歪みの関係は非線形である。したがって、図２の応力−歪み曲線の傾きで与えられる弾性率（ヤング率）は、図３に模式的に示すように、脂肪部位は一定値であるが、乳腺部位は加えられた歪みに応じて変化することが分かる。
【０００５】
一方、非特許文献２には、歪みと弾性率の関係のデータを取得した後、歪み−弾性率曲線を関数で近似し、最小二乗法などにより最もよく近似する曲線からその非線形性を表す非線形パラメタを抽出し、その組織の非線形性を表す情報として評価する試みが提案されている。
【０００６】
このように、組織の硬さの情報を用いて組織の良悪性を鑑別するにあたって、例えば弾性率がどのような圧迫条件の下で計測されたかが非常に重要な情報であり、圧縮量などの圧迫条件が大きく異なる条件で計測された診断対象については、良悪性を誤って診断してしまうことも考えられる。
【０００７】
【特許文献１】特開平５−３１７３１３号号公報
【特許文献２】特開２０００−６０８５３号公報
【非特許文献１】Krouskop T, et al : Elastic Moduli of Breast and Prostate Tissues Under Compression. Ultrasonic Imaging 20 : 260-274, 1998.
【非特許文献２】「超音波による組織の非線形弾性特性の画像化」新田他、電子情報通信学会論文誌 2001/12 Vol.J84-A No.1
【発明の開示】
【発明が解決しようとする課題】
【０００８】
しかしながら、特許文献１、２等に記載された従来の技術では、圧迫の程度に応じて硬さが変化する生体組織の非線形性を考慮していない。そのため、検者は様々な圧迫条件下において複数の弾性画像を取得し、その中から検者が主観的に適切と考える１又は複数の弾性画像に基づいて診断することを余儀なくされている。その結果、従来技術によれば、適切な弾性画像を選択するために経験や熟練を要し、診断の信頼性及び客観性に問題が残る。また、圧迫条件を変化させながら適切と思われる弾性画像を繰返し取得しなければならず、検者に負担がかかるとともに、検査効率が悪いという問題がある。
【０００９】
本発明は、超音波診断装置において、検者の経験や熟練度に依存することなく、客観的で確定的な組織鑑別を可能とする弾性画像を提供することを課題とする。
【課題を解決するための手段】
【００１０】
上記の課題を解決するため、本発明は、被検体の生体組織に圧迫を加えて超音波を送信し、前記被検体から発生する反射エコー信号を計測して得られたフレームデータを用いて、取得時刻が異なる一対の前記フレームデータに基づいて複数の計測点における生体組織の弾性情報と前記各計測点に加わる圧迫の程度に関する計測圧迫条件を求め、求めた前記計測圧迫条件が予め設定された基準圧迫条件を満たすときの前記弾性情報に基づいて弾性画像を生成して表示画面に表示することを特徴とする。
【００１１】
このように、各計測点に加わる計測圧迫条件が予め設定された基準圧迫条件を満たして計測及び生成された弾性画像を用いているから、検者の主観を廃することができ、検者の経験や熟練度に依存することなく、客観的で確定的な組織鑑別を可能とすることができ、かつ検査の効率を向上できる。また、診断に適した弾性画像の選定作業が不要となるから検者の負担が軽減される。
【００１２】
本発明の超音波診断装置は、被検体の生体組織に圧迫を加え超音波を送信して得られる反射エコー信号のフレームデータを取得するフレームデータ取得手段と、該フレームデータ取得手段により取得された取得時刻が異なる一対の前記フレームデータに基づいて複数の計測点における生体組織の弾性情報と前記各計測点に加わる圧迫の程度を示す計測圧迫条件を求める演算手段と、該演算手段により求められた前記弾性情報に基づいて弾性画像を生成する弾性画像生成手段と、前記演算手段により求められた前記計測圧迫条件が予め設定された基準圧迫条件を満たすか否か評価し、前記基準圧迫条件を満たす前記弾性画像を表示手段に表示させる圧迫条件評価手段とを有して構成することができる。
【００１３】
本発明の圧迫条件には、計測点の歪み、計測点の変位量、計測点の応力、被検体の体表に加えられた圧力のいずれか一つを採用することができる。なお、本発明でいう歪み及び変位量は、圧迫による生体組織の歪み変化又は変位のゼロから計測時点までの積分量をいい、その時々の歪み変化又は変位と区別する。また、基準圧迫条件は、生体組織の弾性を評価するのに適した圧迫条件を設定する。
【００１４】
また、弾性情報には、弾性率、粘弾性率、弾性率の非線形性に係る非線形パラメタを採用することができる。
【００１５】
また、フレームデータに対応する画像領域内に設定された関心領域内の計測点について、計測圧迫条件と計測された弾性情報との関係を線図（グラフ）にして表示画面に表示することができる。この場合、線図の基準圧迫条件の位置にマークを表示することが好ましい。
【００１６】
また、弾性画像に関心領域の輪郭を表示するとともに、関心領域に対応付けて関心領域内の複数の計測点の弾性情報の平均値を数値で表示することができる。弾性画像は、各画素の色相又は輝度を各弾性情報に応じて割り付けて生成される。
【発明の効果】
【００１７】
本発明によれば、検者の経験や熟練度に依存することなく、客観的で確定的な組織鑑別を可能とする弾性画像を提供することができる。
【発明を実施するための最良の形態】
【００１８】
以下、本発明の弾性画像生成方法及び装置を実施形態に基づいて説明する。図１に、本発明の弾性画像生成方法を実施する一実施形態の超音波診断装置のブロック構成図を示す。図に示すように、被検体１に当接して用いられる超音波の探触子２は、被検体１との間で超音波を送信及び受信する複数の振動子を有して形成されている。探触子２は、送信回路３から供給される超音波パルスにより駆動される。送受信制御回路４は、探触子２の複数の振動子を駆動する超音波パルスの送信タイミングを制御して、被検体１内に設定される焦点に向けて超音波ビームを形成するようになっている。また、送受信制御回路４は、探触子２の振動子の配列方向に電子的に超音波ビームを走査するようになっている。
【００１９】
一方、探触子２は、被検体１内から発生する反射エコー信号を受信して受信回路５に出力する。受信回路５は、送受信制御回路４から入力されるタイミング信号に従って、反射エコー信号を取り込んで増幅などの受信処理を行う。受信回路５により処理された反射エコー信号は、整相加算回路６において複数の振動子により受信された反射エコー信号の位相を合わせて加算することにより増幅される。整相加算回路６において整相加算された反射エコー信号は、信号処理部７に入力され、ゲイン補正、ログ圧縮、検波、輪郭強調、フィルタ処理等の信号処理がなされる。
【００２０】
信号処理部７により処理された反射エコー信号は、白黒スキャンコンバータ８に導かれて超音波ビームの走査面に対応した２次元の断層像データ（ディジタルデータ）に変換される。これらの信号処理部７と白黒スキャンコンバータ８によって断層像（Ｂモード像）の画像再構成手段が構成される。白黒スキャンコンバータ８から出力される断層像データは、切替加算部９を介して画像表示器１０に供給されてＢモード像が表示されるようになっている。
【００２１】
一方、整相加算回路６から出力される反射エコー信号は、フレームデータ取得部１１に導かれる。フレームデータ取得部１１は、超音波ビームの走査面（断層面）に対応する反射エコー信号群を、フレームデータとして複数フレーム分を取得してメモリなどに格納する。変位計測部１２は、フレームデータ取得部１１に格納されている取得時刻が異なる複数対のフレームデータを順次取り込み、取り込んだ一対のフレームデータに基づいて断層面における複数の計測点の変位ベクトルを求め、変位フレームデータとして弾性情報演算部１３に出力するようになっている。
【００２２】
本実施例の弾性情報演算部１３は、変位フレームデータに基づいて各計測点の生体組織の歪み変化を求める歪み演算部と、歪みフレームデータに基づいて各計測点の生体組織の弾性率を求める弾性率演算部と、歪み演算部で求めた歪み変化を積分して計測圧迫条件の一例である歪みを求める歪み演算部を有して構成されている。弾性情報演算部１３にて求められた弾性率のフレームデータは弾性情報処理部１４に出力され、歪みは圧迫条件評価部１８に出力されるようになっている。
【００２３】
弾性情報処理部１４は、弾性情報演算部１３から入力される各弾性情報のフレームデータに対して、座標平面内におけるスムージング処理、コントラスト最適化処理、フレーム間における時間軸方向のスムージング処理などの様々な画像処理を施して、カラースキャンコンバータ１５に送出するようになっている。
【００２４】
カラースキャンコンバータ１５は、弾性情報処理部１４により処理された弾性率のフレームデータを取り込み、設定された弾性率のカラーマップに従って、フレームデータの画素ごとに色調コードを付与してカラー弾性画像を生成するようになっている。
【００２５】
カラースキャンコンバータ１５により生成されたカラー弾性画像は、切替加算部９を介して画像表示器１０に表示されるようになっている。また、切替加算部９は、白黒スキャンコンバータ８から出力される白黒の断層像と、カラースキャンコンバータ１５から出力されるカラー弾性画像とを入力し、両画像を切り替えていずれか一方を表示させる機能と、両画像の一方を半透明にして加算合成して画像表示器１０に重ねて表示させる機能と、両画像を並べて表示させる機能を有して形成されている。また、切替加算部９から出力される画像データは、装置制御インターフェイス部１９の制御に従ってシネメモリ２０に格納されるようになっている。シネメモリ２０に格納された画像データは、装置制御インターフェイス部１９の制御に従って画像表示器１０に表示されるようになっている。
【００２６】
本実施形態の特徴に係る圧迫条件評価部１８は、弾性情報演算部１３により求められた計測圧迫条件である歪みを取り込み、その歪みが装置制御インターフェイス部１９から入力される基準圧迫条件を満たすか否か判断するようになっている。その判断の結果、取り込んだ歪みが基準圧迫条件を満たす場合は、切替加算部９に指令を送って、カラースキャンコンバータ１５から出力されるカラー弾性画像を画像表示器１０に出力表示させるようになっている。
【００２７】
このように構成される本実施形態の基本的な動作について説明する。まず、探触子２により被検体１に圧迫を加えて被検体１に超音波ビームを走査し、走査面からの反射エコー信号を連続的に受信する。そして、整相加算回路６から出力される反射エコー信号に基づいて、信号処理部７及び白黒スキャンコンバータ８により断層像が再構成され、切替加算器９を介して画像表示器１０に表示される。
【００２８】
一方、フレームデータ取得部１１は、反射エコー信号を取り込んでフレームレートに同期させてフレームデータを繰り返し取得し、内蔵されたフレームメモリ内に時系列順に保存する。そして、取得時刻が異なる一対のフレームデータを単位として、連続的に複数対のフレームデータを選択して変位計測部１２に出力する。変位計測部１２は、選択された一対のフレームデータを1次元もしくは２次元相関処理し、走査面における複数の計測点の変位を計測して変位フレームデータを生成する。この変位ベクトルの検出法としては、例えば特開平５−３１７３１３号公報等に記載されているブロックマッチング法又はグラジェント法が知られている。ブロックマッチング法は、画像を例えばＮ×Ｎ画素からなるブロックに分け、現フレーム中の着目しているブロックに最も近似しているブロックを前フレームから探索し、これに基づいて計測点の変位を求める。また、一対のＲＦ信号フレームデータの同一領域における自己相関を計算して変位を算出することができる。
【００２９】
弾性情報演算部１３は、変位フレームデータを取り込んで、各計測点の歪み変化を求め、求めた歪み変化に基づいて弾性情報である弾性率を演算し、弾性率フレームデータを弾性情報処理部１４に出力する。歪み変化の演算は、公知のように変位を空間微分することによって計算される。また、求めた歪み変化に基づいて各計測点の弾性率を演算する。弾性率を求める場合は、圧力計測部１７により計測された圧力の計測値を取り込み、これに基づいて各計測点における応力を演算する。圧力計側部１７は、探触子2の超音波送受信面と被検体1との間に設けられた圧力センサ１６により検出された圧力に基づいて、被検体１内部の計測点における応力を演算する。つまり、弾性情報演算部１３の弾性率演算部は、各計測点における応力と、弾性情報演算部１３で求めた歪みフレームデータから走査面上の各計測点の弾性率Ｅ（例えば、ヤング率）を演算し、弾性情報処理部１４に出力する。
【００３０】
弾性情報処理部１４は、入力される弾性率をスムージング処理などの処理を施してカラースキャンコンバータ１５に出力する。カラースキャンコンバータ１５は、弾性情報に基づいてカラー弾性画像を生成する。カラー画像は、例えば２５６階調化による色調のグラデーションで画素単位ごとにフレームデータの弾性率に応じて色付けされる。なお、カラースキャンコンバータ１５に代えて、白黒スキャンコンバータを用いることができる。この場合は、弾性率が大きい領域は輝度を明るく、逆に小さい領域は輝度を暗くするなどにより、良性又は悪性を鑑別できるようにすることができる。
【００３１】
以下に、本実施形態を用いて本発明の弾性画像生成方法を実施する具体的な実施例を説明する。なお、各実施例は、本実施形態の特徴部である弾性情報演算部１３、圧迫条件評価部１８、装置制御インターフェイス部１９及びシネメモリ２０等により実施される。
【実施例１】
【００３２】
本発明の弾性画像生成方法の一実施例として、図４に示す生体組織１〜５を含む部位に圧迫を加えて弾性画像を生成する場合を説明する。本実施例では、弾性情報演算部１３において弾性フレームデータとして生体組織１〜５の各部位における計測点の弾性率を弾性情報として求めるものとする。なお、図４は、弾性フレームデータと弾性画像との位置関係を対応付けるために、１フレームの画像の体裁で示している。
【００３３】
まず、圧迫条件評価部１８には、装置制御インターフェイス部１９から検者によって関心領域（以下、ＲＯＩという。）としてＲＯＩ１、ＲＯＩ２が入力設定される。圧迫条件評価部１８は、弾性情報演算部１３からＲＯＩ１、ＲＯＩ２の中に分布する各計測点の弾性率データ（弾性率の数値データ）と、ＲＯＩ１、ＲＯＩ２の中に分布する各計測点の歪みの数値データを取り込む。本実施例では、圧迫の程度を評価するための圧迫条件として歪みを用いている。また、歪み（ε＝ΣΔε）は、圧縮ゼロから現時点までの歪み変化の積分値であり、一対のフレームデータ間ごとに得られる歪み変化（Δε）とは区別する。
【００３４】
圧迫条件評価部１８は、取り込んだ歪みと弾性率に基づいて、ＲＯＩ１、ＲＯＩ２のそれぞれについて、例えば、歪みε（％）の平均値と弾性率Ｅ（ｋＰａ）の平均値を演算し、それらの関係を図５に示すようにグラフにプロットする。そして、歪み（圧迫条件）−弾性情報（弾性率）曲線のグラフ画像を構築し、圧迫条件評価部１８に記憶するとともに、切替加算部９に送出して画像表示器１０に表示する。これにより、画像表示器１０には、ＲＯＩ１、ＲＯＩ２内の組織に関する圧迫条件−弾性情報曲線がグラフ画像として表示される。
【００３５】
例えば、基準圧迫条件を歪み４〜６％に設定すると、図５に示すように、基準圧迫条件の範囲が圧迫条件−弾性情報曲線に表示される。一方、圧迫条件評価部１８は、基準圧迫条件を満たした時点において計測されたＲＯＩ１、ＲＯＩ２の弾性率Ｅ１、Ｅ２を解析する。そして、図４に示したＲＯＩ１、ＲＯＩ２に対応する組織の弾性診断結果として、図５に示すように弾性率Ｅ１、Ｅ２の数値データをグラフ画像に表示して、ＲＯＩ１、ＲＯＩ２の基準圧迫条件を満たすときの弾性率を弾性診断結果として確定する。なお、一般に、ＲＯＩ１とＲＯＩ２の計測圧迫条件は異なることがあるから、基準圧迫条件が同じであっても、同じ計測時刻においてＲＯＩ１とＲＯＩ２が基準圧迫条件を満たすというわけではない。
【００３６】
ここで、基準圧迫条件は、範囲として設定する方法に限らない。例えば、歪み５％を基準圧迫条件として設定し、この歪み５％を跨ぐ前後において計測された弾性率の値から線形補間することにより歪み５％における弾性率の値を求めるようにすることができる。
【００３７】
また、例えば、圧迫操作を繰り返して行う弾性情報を計測する過程で、基準圧迫条件を満たして計測された弾性率の値により、基準圧迫条件を満たすたびに図５の弾性率を更新するようにすることができる。また、基準圧迫条件を満たして計測された現在までの弾性率の値との加算平均値を表示させるようにすることができる。
【００３８】
また、本実施例では、弾性情報としてヤング率に代表される弾性率を例に示したが、粘弾性率、非線形性などの組織の弾性特性に基づいた指標であれば用いることができる。
【００３９】
また、計測圧迫条件は、ＲＯＩ内の歪みに限らず、変位計測部１２で計測された変位量（変位ゼロからの積算を含む）、圧力計測部１２で計測された圧迫の圧力に基づいて導出されたＲＯＩ内に伝わった応力など、圧迫の状態を反映した情報を用いることができる。
【００４０】
上述したように、本実施例によれば、生体組織の画像にＲＯＩを設定して弾性情報の計測操作を実施するだけで、そのＲＯＩに対応する組織に加えられる計測圧迫条件が評価に適切な基準圧迫条件を満たしているときに取得された弾性情報である弾性率が数値により表示される。したがって、経験や熟練がない検者であっても、ＲＯＩに係る生体組織の適切な弾性情報を得ることができるから、信頼性及び客観性を有する診断を行うことができる。また、圧迫条件を変化させながら適切と思われる弾性画像を繰返し取得する操作が不要となり、検者の負担を軽減して、検査効率を向上できる。また、本実施例によれば、超音波によるリアルタイム性を損なうことなく、確定的な診断情報を臨床現場において効率よく取得することができる。
【実施例２】
【００４１】
実施例１では、ＲＯＩの弾性情報を数値データとして画像に出力する例を示した。しかし、本発明により生成する弾性画像は、これに限られるものではなく、図６に示す実施例２のように、基準圧迫条件を満たした時点において計測されたＲＯＩ内の弾性率に基づき、弾性率に応じてＲＯＩ内の画像を階調化した弾性画像を生成して表示することができる。
【００４２】
例えば、実施例１で設定したＲＯＩ内の計測圧迫条件が基準圧迫条件を満たしたときに、ＲＯＩ内において計測された弾性率の大きさに応じて、白黒輝度又は色相に階調を与えた画像を生成して表示する。このとき、弾性率の数値データとひずみの数値データも画像上に表示するようにしてもよい。
【実施例３】
【００４３】
実施例２では、各ＲＯＩについて、各計測点に対応する画素単位で弾性率を階調化して表示する例を説明した。しかし、図７の実施例３に示すように、実施例２のＲＯＩに代えて、計測領域の全体に画素を構築できる程度の小さな局所的ＲＯＩを格子状に複数配置することができる。この場合は、各局所的ＲＯＩ毎に独立に、計測圧迫条件が基準圧迫条件を満たしたときの各局所的ＲＯＩの弾性率の大きさを求め、局所的ＲＯＩごとに白黒輝度又は色相に階調を与えて、観測領域全域にわたって弾性画像を形成して表示する。なお、局所的ＲＯＩは、１又は複数の計測点を含む領域として設定することができる。
【００４４】
また、基準圧迫条件を満たして割当てられた各局所的ＲＯＩの弾性画像の階調は、繰り返し計測される弾性率の計測圧迫条件が基準圧迫条件を満たすたびに更新して表示されるようにすることができる。また、基準圧迫条件を満たして取得されたこれまでの弾性率の値の加算平均値に基づいて、弾性画像の階調を更新して表示すれば、弾性画像の精度を向上させることができる。
【００４５】
また、計測中に計測圧迫条件が基準圧迫条件を満たさなかったＲＯＩについては、弾性率の大きさを反映した階調が付与されないようにすれば、弾性評価不可能であったことを認識できる。
【実施例４】
【００４６】
実施例１〜３では、弾性情報として弾性率を用いる例を説明した。本実施例４では、図３に示した歪み−弾性率曲線の非線形性を表す非線形パラメタを弾性情報として画像化する例を説明する。つまり、計測圧迫条件が基準圧迫条件を満たした場合における各計測点の非線形パラメタの値を求め、求めた非線形パラメタに基づいて階調化した弾性画像を生成する。
【００４７】
前述したとおり、非特許文献１には、乳腺組織の弾性特性の非線形性に関して、脂肪組織は大きな歪みの範囲までほぼ一定の弾性率の線形な応答を示すのに対し、繊維組織や浸潤性腫瘍は歪みの増加に伴い弾性率も顕著に増加する歪み硬化が認められることが報告されている。また、繊維組織よりも浸潤性の悪性腫瘍の方が、歪み硬化の程度が大きく、非線形の程度も大きくなる。また、一般に生体組織は非線形の弾性応答を示し、例えば、歪み−弾性率の関係は、図３に示した曲線に従うことが知られている。非特許文献２には、この歪みと弾性率の関係のデータを取得した後、所定の関数でその曲線を近似し、最小二乗法などにより最もよく近似する曲線からその非線形パラメタを抽出し、その組織の非線形性を表す情報として評価する試みが提案されている。
【００４８】
つまり、非特許文献２によれば、Ｅをヤング率と歪みε（＝ΣΔε）の非線形性を解析する目的において、圧迫に対する組織の応答の曲線を、次式（１）のように１次関数として仮定している。
【００４９】
Ｅ＝Ｅ０＋α×ε （１）
ここで、Ｅ０は一定値、αは非線形パラメタであり、ヤング率−歪み曲線の傾きである。しかし、図８に示すように、大きな歪みになるほど、実際の組織の非線形性が顕著になり、非線形パラメタαは圧迫条件に大きく依存する。したがって、（１）式のようにＥがεの１次関数であるとする仮定する非線形パラメタαは、評価区間を定めなければ定量情報とは成り得ず、確定診断には採用できない。
【００５０】
そこで、本実施例では、実施例１〜３と同様、基準圧迫条件を設定し、計測圧迫条件が基準圧迫条件を満たして計測された弾性率と歪みとのデータに基づいて非線形パラメタαを求め、求めた非線形パラメタαを数値情報や画像情報として検者に提供するようにすることを特徴とする。具体的には、実施例１〜３における「弾性率」を「非線形パラメタ」に置き換えればよい。
【００５１】
さらに、非線形パラメタαを一層忠実に求めるためには、広範囲の計測圧迫条件において、組織の非線形応答を忠実に近似する曲線として、次式（２）に示す高次の関数により近似することが好ましい。γは、２以上の自然数である。
【００５２】
Ｅ＝Ｅ０＋α×εγ （２）
（２）式による非線形パラメタαも、組織の非線形性が大きいほど大きな値となり、線形に近づくほど限りなくゼロに近くなる。さらに、一般化すると、次式（３）の指数関数で近似できる。
【００５３】
Ｅ＝（Ｅ０−１）＋ ｅｘｐ（α×ε) （３）
これらの式（２）、（３）で示される関係式を仮定すれば、計測圧迫条件の評価区間を任意に設定しても、高精度に確定した非線形パラメタαを得ることができ、定量的な評価が実現できる。
【００５４】
例えば、従来は、図８のようにαの値が圧迫条件に依存するため、歪み−弾性率の関係のデータを蓄積した後、オフライン処理により計測圧迫条件について適当な評価区間を確認して解析する必要がある。この点、本実施例によれば、計測圧迫条件に依存しない確定的な非線形パラメタαが得られる。つまり、時刻ｔ（０）において計測された歪み−弾性率の関係の（ε（０）、Ｅ（０)）と、過去の時刻ｔ（−１）、ｔ（−２）、…において得られた歪み−弾性率の関係（ε（−１）、Ｅ（−１））、（ε（−２）、Ｙ（−２））、 …を歪み−弾性率の関係の蓄積データとして適用し、時刻ｔ（０）においてリアルタイムに非線形パラメタαを評価しても、その非線形パラメタαの数値情報、画像情報は圧迫条件に依存しない客観的な情報として得ることができる。
【００５５】
したがって、本実施例４によれば、図７で説明した場合と同様、非線形パラメタαを弾性情報とし、リアルタイムの弾性画像を生成して表示することが可能である。その結果、例えば、Ｂモード像や弾性率画像によって認識できない悪性組織であっても、非線形パラメタの評価によってのみ異常な組織を検出することができる場合があるので、このような場合に検査において異常組織領域を見逃してしまうことはなく、検者の主観に依らない診断を行うことができる。
【実施例５】
【００５６】
実施例１〜３においては、ＲＯＩに属する組織が圧迫操作により移動することについて考慮していない。本実施例５では、圧迫操作により移動する組織に追従させてＲＯＩを移動させて、弾性情報の評価の精度を向上させる。つまり、画像に対するＲＯＩの位置を固定すると、図９（ａ）〜（ｃ）に示すように、圧迫操作の経過に応じてＲＯＩ（ｔ０）、ＲＯＩ（ｔ１）、ＲＯＩ（ｔ２）内に入ってくる組織が変化してしまうことがあり、計測値の精度が悪くなる。
【００５７】
そこで、本実施例５では、図９（ａ’）〜（ｃ’）に示すように、圧迫の程度に応じて、ＲＯＩ内に属する組織の行き先を追従して、ＲＯＩ（ｔ０）、ＲＯＩ（ｔ１）、ＲＯＩ（ｔ２）を移動させることにより、任意の圧迫条件の下において、同一組織領域の情報のみを抽出することが可能となる。その結果、評価精度を向上させることが可能になる。なお、組織の動きに追従させてＲＯＩを移動させる処理は、前述した組織の変位ベクトルを求めるブロックマッチング法又はグラジェント法を適用することができる。また、局所的ＲＯＩや計測点についても、組織の移動に応じて追従させることができる。
【実施例６】
【００５８】
ここで、上記実施例１〜５等により生成された歪み−弾性率曲線のグラフ画像、基準圧迫条件を満たしたときの弾性率画像、非線形パラメタα画像等の弾性画像を、画像表示器１０に表示する際の表示形態の実施例について説明する。
【００５９】
切替加算部９と画像表示器１０の機能により、２画面を同時に表示可能にすれば、例えば、右画面に現時刻におけるＢモード像や弾性画像をリアルタイムで表示し、左画面に計測圧迫条件−弾性情報のグラフ、基準圧迫条件を満たして計測された弾性画像、もしくは非線形パラメタ画像を表示することができる。また、複数画面の内のいずれかに圧迫条件に依存しない定量の弾性情報を表示する画像を表示することができる。
【００６０】
このように、表示画像の組み合わせは、その他にもカラードプラなど従来の診断画像を含め、任意に選択できるようにすることができる。また、一つの画面に左右に並べて表示する方法に限らず、上下に並べて表示してもよい。さらに、一つの画面内に設けた小窓に別の画像を表示するような形式を採用することもできる。
【実施例７】
【００６１】
一般に、リアルタイム弾性画像の計測においては、超音波診断装置にフリーズ指令を与えると、フリーズまでに取得された時間的に連続した所定フレーム数の弾性画像がシネメモリ２０に確保される。そして、シネメモリ２０に確保された弾性画像を連続的に読み出して画像表示器１０に動画として表示することができるようになっている。
【００６２】
そこで、本実施例７では、圧迫操作を行いながら計測した複数の弾性画像をシネメモリ２０に格納するとともに、基準圧迫条件を満たした計測圧迫条件で得られた弾性画像のフレームに、基準圧迫条件を満たした弾性画像であることを示すマークを表示するようにする。例えば、基準圧迫条件を満たした弾性画像には、丸印（○）を画像のどこかに表示する。また、基準圧迫条件に対して計測圧迫条件が小さければ上向きの矢印（↑）を、大きければ下向きの矢印（↓）を、弾性画像に表示する。また、それらのマーク表示に代えて又はマーク表示と共に、音声によって検者に報知するようにしてもよい。
【００６３】
本実施例によれば、検者は、超音波診断装置からリアルタイムにフィードバックされるガイド情報により、計測圧迫条件が基準圧迫条件を満たす方向に調整することができるから、適切な弾性画像を取得することが可能となる。また、同時に、診断に適用する画像や弾性に関する数値情報を検者の主観によらずに確定的に決定することができる。
【実施例８】
【００６４】
実施例７では、本発明の弾性画像生成方法によりリアルタイムで生成された弾性画像を、画像表示器１０に連続的に表示する例を示したが、本発明はこれに限らず、計測圧迫条件が基準圧迫条件を満たしたときの弾性情報に対応する弾性画像のみを生成して画像表示器１０に表示するようにする。つまり、計測圧迫条件が基準圧迫条件を満たさない弾性画像は表示しないようにする。
【実施例９】
【００６５】
実施例７では、基準圧迫条件に対する計測圧迫条件の大小に応じて弾性画像にマークを表示するようにしたが、本発明はこれに限らず、計測圧迫条件が基準圧迫条件に近づくほど、カラー弾性画像の明るさが強くなるなど、弾性画像の見易さが増加するなどの工夫をすることができる。例えば、Ｂモード像に半透明のカラー弾性画像を重ね合わせた画像を生成する場合は、カラー弾性画像のブレンディング比率を計測圧迫条件が基準圧迫条件に近づくほど大きくするようにすることもできる。
【実施例１０】
【００６６】
また、計測圧迫条件が基準圧迫条件を満たしたときの弾性画像、関心部位の弾性率の数値情報、その他診断の判断材料となる情報は、弾性診断の結果として、超音波診断装置に備えられている各患者の診断レポート機能に自動的に反映するようにすることができる。
【００６７】
また、上記の各実施例においては、対象組織を圧迫することにより弾性情報を取得する例について説明したが、本発明の弾性画像生成方法は、例えば、加振映像法（Y. Yamakoshi et al : Ultrasonic Imaging of Internal Vibration of Soft Tissue under Forced Vibration, IEEE Trans UFFC 1990; 37; 45-53.）などの手法によって得られた弾性情報に対しても適用することができる。
【図面の簡単な説明】
【００６８】
【図１】本発明の弾性画像生成方法を実施する一実施形態の超音波診断装置のブロック構成図である。
【図２】乳腺組織の乳腺部位と脂肪部位を圧迫した際の歪みと応力の関係を示す線図である。
【図３】乳腺組織の乳腺部位と脂肪部位を圧迫した際の歪みと弾性率（ヤング率）の関係を模式的に示す線図である。
【図４】生体組織に圧迫を加えて弾性画像を生成する方法を説明する図である。
【図５】実施例１の弾性画像生成方法により生成される弾性画像の歪み―弾性率曲線のグラフ画像である。
【図６】実施例２の弾性画像生成方法により生成される弾性画像の一例である。
【図７】実施例３の弾性画像生成方法により生成される弾性画像の一例である。
【図８】実施例４の弾性画像生成方法により生成する非線形パラメタが圧迫条件に大きく依存することを説明する線図である。
【図９】圧迫操作によりＲＯＩ内に属する組織の移動に追従させてＲＯＩを移動させる実施例を説明する図である。
【符号の説明】
【００６９】
１ 被検体
２ 探触子
３ 送信回路
４ 送受信制御回路
５ 受信回路
６ 整相加算回路
７ 信号処理部
８ 白黒スキャンコンバータ
９ 切替加算部
１０ 画像表示器
１１ フレームデータ取得部
１２ 変位計測部
１３ 弾性情報演算部
１５ カラースキャンコンバータ
１６ 圧力センサ
１７ 圧力計測部
１８ 圧迫条件評価部
１９ 装置制御インターフェイス部
２０ シネメモリ
【特許請求の範囲】
【請求項１】
被検体の生体組織に圧迫を加えて超音波を送信し、前記被検体から発生する反射エコー信号を計測して得られたフレームデータを用いて、取得時刻が異なる一対の前記フレームデータに基づいて複数の計測点における生体組織の弾性情報と前記各計測点に加わる圧迫の程度に関する計測圧迫条件を求め、求めた前記計測圧迫条件が予め設定された基準圧迫条件を満たすときの前記弾性情報に基づいて表示画面に表示する弾性画像を生成する弾性画像生成方法。
【請求項２】
被検体の生体組織に圧迫を加え超音波を送信して得られる反射エコー信号のフレームデータを取得するフレームデータ取得手段と、該フレームデータ取得手段により取得された取得時刻が異なる一対の前記フレームデータに基づいて複数の計測点における生体組織の弾性情報と前記各計測点に加わる圧迫の程度を示す計測圧迫条件を求める演算手段と、該演算手段により求められた前記弾性情報に基づいて弾性画像を生成する弾性画像生成手段と、前記演算手段により求められた前記計測圧迫条件が予め設定された基準圧迫条件を満たすか否か評価し、前記基準圧迫条件を満たす前記弾性画像を表示手段に表示させる圧迫条件評価手段とを有してなる超音波診断装置。
【請求項３】
請求項２に記載の超音波診断装置において、
前記圧迫条件評価手段は、前記フレームデータに対応する画像領域内に設定された関心領域内の前記計測点における前記計測圧迫条件と該計測圧迫条件に対応する前記弾性情報との関係のグラフ画像を生成し、該グラフ画像の前記基準圧迫条件に対応する位置にマークを表示して前記表示画面に表示させることを特徴とする超音波診断装置。
【請求項４】
請求項３に記載の超音波診断装置において、
前記圧迫条件評価手段は、前記関心領域内の複数の計測点の前記弾性情報の平均値を前記グラフ画像に数値で表示することを特徴とする超音波診断装置。
【請求項５】
請求項２に記載の超音波診断装置において、
前記圧迫条件評価手段は、前記弾性画像に前記関心領域を表示し、かつ該関心領域に対応付けて該関心領域内の複数の計測点の前記弾性情報の平均値を数値で表示することを特徴とする超音波診断装置。
【請求項６】
請求項２に記載の超音波診断装置において、
前記計測圧迫条件は、前記計測点の歪み、前記計測点の変位量、前記計測点の応力、前記被検体の体表に加えられた圧力のいずれか一つであり、
前記弾性情報は、弾性率、粘弾性率、弾性率の非線形性に係る非線形パラメタ、ヒステリシスのいずれか一つであることを特徴とする超音波診断装置。
【請求項７】
請求項２に記載の超音波診断装置において、
前記圧迫条件評価手段は、前記演算手段により求められた前記計測圧迫条件が予め設定された基準圧迫条件を満たすか否か評価し、前記基準圧迫条件を満たす前記弾性画像にその旨を表示するマークを付するとともに、前後の複数の前記弾性画像に前記計測圧迫条件が前記基準圧迫条件よりも小さいか大きいかを表示するマークを付して前記表示手段に表示させることを特徴とする超音波診断装置。
【請求項８】
請求項２に記載の超音波診断装置において、
前記圧迫条件評価手段は、前記演算手段により求められた前記計測圧迫条件が予め設定された基準圧迫条件を満たすか否か評価し、前記基準圧迫条件を満たす前記弾性画像の輝度又は色相と、前後の複数の前記弾性画像の輝度又は色相を異ならせて前記表示手段に表示させることを特徴とする超音波診断装置。
【請求項１】
被検体の生体組織に圧迫を加えて超音波を送信し、前記被検体から発生する反射エコー信号を計測して得られたフレームデータを用いて、取得時刻が異なる一対の前記フレームデータに基づいて複数の計測点における生体組織の弾性情報と前記各計測点に加わる圧迫の程度に関する計測圧迫条件を求め、求めた前記計測圧迫条件が予め設定された基準圧迫条件を満たすときの前記弾性情報に基づいて表示画面に表示する弾性画像を生成する弾性画像生成方法。
【請求項２】
被検体の生体組織に圧迫を加え超音波を送信して得られる反射エコー信号のフレームデータを取得するフレームデータ取得手段と、該フレームデータ取得手段により取得された取得時刻が異なる一対の前記フレームデータに基づいて複数の計測点における生体組織の弾性情報と前記各計測点に加わる圧迫の程度を示す計測圧迫条件を求める演算手段と、該演算手段により求められた前記弾性情報に基づいて弾性画像を生成する弾性画像生成手段と、前記演算手段により求められた前記計測圧迫条件が予め設定された基準圧迫条件を満たすか否か評価し、前記基準圧迫条件を満たす前記弾性画像を表示手段に表示させる圧迫条件評価手段とを有してなる超音波診断装置。
【請求項３】
請求項２に記載の超音波診断装置において、
前記圧迫条件評価手段は、前記フレームデータに対応する画像領域内に設定された関心領域内の前記計測点における前記計測圧迫条件と該計測圧迫条件に対応する前記弾性情報との関係のグラフ画像を生成し、該グラフ画像の前記基準圧迫条件に対応する位置にマークを表示して前記表示画面に表示させることを特徴とする超音波診断装置。
【請求項４】
請求項３に記載の超音波診断装置において、
前記圧迫条件評価手段は、前記関心領域内の複数の計測点の前記弾性情報の平均値を前記グラフ画像に数値で表示することを特徴とする超音波診断装置。
【請求項５】
請求項２に記載の超音波診断装置において、
前記圧迫条件評価手段は、前記弾性画像に前記関心領域を表示し、かつ該関心領域に対応付けて該関心領域内の複数の計測点の前記弾性情報の平均値を数値で表示することを特徴とする超音波診断装置。
【請求項６】
請求項２に記載の超音波診断装置において、
前記計測圧迫条件は、前記計測点の歪み、前記計測点の変位量、前記計測点の応力、前記被検体の体表に加えられた圧力のいずれか一つであり、
前記弾性情報は、弾性率、粘弾性率、弾性率の非線形性に係る非線形パラメタ、ヒステリシスのいずれか一つであることを特徴とする超音波診断装置。
【請求項７】
請求項２に記載の超音波診断装置において、
前記圧迫条件評価手段は、前記演算手段により求められた前記計測圧迫条件が予め設定された基準圧迫条件を満たすか否か評価し、前記基準圧迫条件を満たす前記弾性画像にその旨を表示するマークを付するとともに、前後の複数の前記弾性画像に前記計測圧迫条件が前記基準圧迫条件よりも小さいか大きいかを表示するマークを付して前記表示手段に表示させることを特徴とする超音波診断装置。
【請求項８】
請求項２に記載の超音波診断装置において、
前記圧迫条件評価手段は、前記演算手段により求められた前記計測圧迫条件が予め設定された基準圧迫条件を満たすか否か評価し、前記基準圧迫条件を満たす前記弾性画像の輝度又は色相と、前後の複数の前記弾性画像の輝度又は色相を異ならせて前記表示手段に表示させることを特徴とする超音波診断装置。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【公開番号】特開２００７−２８２９３２（Ｐ２００７−２８２９３２Ａ）
【公開日】平成１９年１１月１日（２００７．１１．１）
【国際特許分類】
【出願番号】特願２００６−１１５１２０（Ｐ２００６−１１５１２０）
【出願日】平成１８年４月１９日（２００６．４．１９）
【出願人】（０００１５３４９８）株式会社日立メディコ (1,613)
【Ｆターム（参考）】
【公開日】平成１９年１１月１日（２００７．１１．１）
【国際特許分類】
【出願日】平成１８年４月１９日（２００６．４．１９）
【出願人】（０００１５３４９８）株式会社日立メディコ (1,613)
【Ｆターム（参考）】
[ Back to top ]