(57)【要約】【課題】超音波プローブの移動状況を容易に判別すること。【解決手段】実施例の超音波診断装置は、変化量算出部１５ｃ、推定部１５ｄおよび制御部１７を備える。変化量算出部１５ｃは、超音波プローブ１が送信した超音波の反射波に基づいて時系列に沿って生成された複数の超音波画像間（複数の生データ間、または、複数の表示用超音波画像間）で、複数の局所領域（ＲＯＩ）それぞれのパターンの変化量を算出する。推定部１５ｄは、変化量算出部１５ｃにより算出された複数の局所領域（ＲＯＩ）それぞれのパターンの変化量に基づいて、各超音波画像における走査断面の回転軸の位置を推定する。そして、制御部１７は、推定部１５ｄにより推定された回転軸を各超音波画像に重畳してモニタ２に表示するように制御する。【選択図】図２
2012030053.tif 000002
【特許請求の範囲】
【請求項１】
超音波プローブが送信した超音波の反射波に基づいて時系列に沿って生成された複数の超音波画像間で、複数の局所領域それぞれのパターンの変化量を算出する変化量算出部と、
前記変化量算出部により算出された前記複数の局所領域それぞれのパターンの変化量に基づいて、各超音波画像における走査断面の回転軸の位置を推定する推定部と、
前記推定部により推定された前記回転軸を各超音波画像に重畳して所定の表示部に表示するように制御する表示制御部と、
を備える超音波診断装置。
【請求項２】
前記複数の局所領域それぞれは、前記超音波画像における縦方向に沿って設定される、請求項１の記載の超音波診断装置。
【請求項３】
前記推定部は、前記複数の局所領域それぞれのパターンの変化量を各局所領域の超音波画像内の位置に応じてプロットしたグラフを生成し、当該生成したグラフに基づいて、前記回転軸の位置を推定し、
前記表示制御部は、さらに、前記推定部が生成した前記グラフを前記所定の表示部に表示するように制御する、請求項１に記載の超音波診断装置。
【請求項４】
前記推定部は、前記回転軸の位置を推定するために用いたパターンの変化量が所定の閾値より大きい場合、前記回転軸の位置が推定不可であると判定し、
前記表示制御部は、前記推定部の推定結果が推定不可である場合、前記回転軸を非表示とするように制御する、請求項１に記載の超音波診断装置。
【請求項５】
前記推定部は、前記グラフにおける複数の局所領域それぞれのパターンの変化量の変化率に基づいて、前記走査断面が回転されている速さを推定し、
前記表示制御部は、さらに、前記推定部が推定した前記走査断面が回転されている速さを前記所定の表示部に表示するように制御する、請求項３に記載の超音波診断装置。
【請求項６】
前記変化量算出部は、前記パターンの変化量として相互相関係数を算出する、請求項１に記載の超音波診断装置。
【請求項７】
超音波プローブが送信した超音波の反射波に基づいて時系列に沿って生成された複数の超音波画像間で、複数の局所領域それぞれのパターンの変化量を算出する変化量算出部と、
前記変化量算出部により算出された前記複数の局所領域それぞれのパターンの変化量に基づいて、各超音波画像における走査断面の回転軸の位置を推定する推定部と、
前記推定部により推定された前記回転軸を各超音波画像に重畳して所定の表示部に表示するように制御する表示制御部と、
を備える画像処理装置。
【請求項８】
前記複数の局所領域それぞれは、前記超音波画像における縦方向に沿って設定される、請求項７に記載の画像処理装置。
【請求項９】
前記推定部は、前記複数の局所領域それぞれのパターンの変化量を各局所領域の超音波画像内の位置に応じてプロットしたグラフを生成し、当該生成したグラフに基づいて、前記回転軸の位置を推定し、
前記表示制御部は、さらに、前記推定部が生成した前記グラフを前記所定の表示部に表示するように制御する、請求項７に記載の画像処理装置。
【請求項１０】
前記推定部は、前記回転軸の位置を推定するために用いたパターンの変化量が所定の閾値より大きい場合、前記回転軸の位置が推定不可であると判定し、
前記表示制御部は、前記推定部の推定結果が推定不可である場合、前記回転軸を非表示とするように制御する、請求項７に記載の画像処理装置。
【請求項１１】
前記推定部は、前記グラフにおける複数の局所領域それぞれのパターンの変化量の変化率に基づいて、前記走査断面が回転されている速さを推定し、
前記表示制御部は、さらに、前記推定部が推定した前記走査断面が回転されている速さを前記所定の表示部に表示するように制御する、請求項９に記載の画像処理装置。
【請求項１２】
前記変化量算出部は、前記パターンの変化量として相互相関係数を算出する、請求項７に記載の画像処理装置。
【請求項１３】
変化量算出部が、超音波プローブが送信した超音波の反射波に基づいて時系列に沿って生成された複数の超音波画像間で、複数の局所領域それぞれのパターンの変化量を算出し、
推定部が、前記変化量算出部により算出された前記複数の局所領域それぞれのパターンの変化量に基づいて、各超音波画像における走査断面の回転軸の位置を推定し、
表示制御部が、前記推定部により推定された前記回転軸を各超音波画像に重畳して所定の表示部に表示するように制御する、
ことを含む画像処理方法。
【請求項１４】
前記複数の局所領域それぞれは、前記超音波画像における縦方向に沿って設定される、請求項１３に記載の画像処理方法。
【請求項１５】
前記推定部は、前記複数の局所領域それぞれのパターンの変化量を各局所領域の超音波画像内の位置に応じてプロットしたグラフを生成し、当該生成したグラフに基づいて、前記回転軸の位置を推定し、
前記表示制御部は、さらに、前記推定部が生成した前記グラフを前記所定の表示部に表示するように制御する、請求項１３に記載の画像処理方法。
【請求項１６】
前記推定部は、前記回転軸の位置を推定するために用いたパターンの変化量が所定の閾値より大きい場合、前記回転軸の位置が推定不可であると判定し、
前記表示制御部は、前記推定部の推定結果が推定不可である場合、前記回転軸を非表示とするように制御する、請求項１３の記載の画像処理方法。
【請求項１７】
前記推定部は、前記グラフにおける複数の局所領域それぞれのパターンの変化量の変化率に基づいて、前記走査断面が回転されている速さを推定し、
前記表示制御部は、さらに、前記推定部が推定した前記走査断面が回転されている速さを前記所定の表示部に表示するように制御する、請求項１５に記載の画像処理方法。
【請求項１８】
前記変化量算出部は、前記パターンの変化量として相互相関係数を算出する、請求項１３に記載の画像処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施の形態は、超音波診断装置、画像処理装置および画像処理方法に関する。
【背景技術】
【０００２】
従来、超音波診断装置は、簡便な操作性、被曝のおそれがない非侵襲性、装置規模の小ささなどの様々な利点を備えた医用画像診断装置として、今日の医療において重要な役割を果たしている。
【０００３】
すなわち、超音波診断装置は、超音波プローブを体表から当てるだけの簡便な操作により、例えば、心臓の拍動や胎児の動きといった検査対象の動きの様子をリアルタイムで表示することができる。また、超音波診断装置は、非侵襲性であることから安全性が高く、繰り返して検査を行なうことができる。また、超音波診断装置は、Ｘ線診断装置、Ｘ線ＣＴ（Computed  Tomography）装置、ＭＲＩ（Magnetic  Resonance  Imaging）装置などの他の医用画像診断装置に比べ装置規模が小さく、ベッドサイドへ移動しての検査も容易に行なうことができる。また、被曝のおそれがない超音波診断装置には、片手で持ち運べる程度に小型化された装置も開発されており、かかる超音波診断装置は、産科や在宅医療などの医療現場においても容易に使用することができる。
【０００４】
超音波診断装置のリアルタイム性を利用した診断法としては、例えば、以下のような手技が行われる。断層像である超音波画像に腫瘤性病変が見つかった場合、操作者は、超音波プローブを回転（例えば、９０度回転）させて、腫瘤性病変を異なる断面で観察したり、腫瘤性病変のサイズを計測したりする。かかる観察により、操作者は、腫瘤の形状が円形であるか楕円形であるかといった判定を行なうことができる。また、多くの血管の断面は、円形である。したがって、一つの超音波画像上で円形として見える部位も、例えば、超音波プローブを９０度回転させた状態で生成された超音波画像を参照することで、操作者は、かかる部位が管状であることを確認することができる。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０００−１３２６６４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
ところで、上記した操作者の手技をともなう観察は、リアルタイムで行われるため、診断用の超音波画像は、超音波プローブの回転中も含めた動画像として記録されることが多い。しかしながら２次元の断層像においては、下記のような不便さを伴うことがある。
【０００７】
すなわち、超音波プローブを回転する場合、超音波の走査面は、時系列に沿って変化するために、超音波画像を参照しても、超音波プローブをどのように動かしているのか判別することが困難となる。例えば、超音波画像を参照しても、読影者は、超音波プローブを回転移動させた場合と、超音波プローブを回転させずに垂直方向に移動させた場合とでは共に走査面が変化するために、両者の違いを区別することが困難である。特に、記録された動画像を事後に閲覧する場合や、操作者以外の人間が動画像を参照して読影する場合、超音波画像だけから超音波プローブの移動状況を判断することは、困難である。一方、超音波プローブを断層面に沿って平行移動させた場合は、動画像上では、超音波画像内のパターンが保持されたまま右か左かへ移動するので、読影者は、超音波プローブの移動状況を容易に判断できる。
【課題を解決するための手段】
【０００８】
実施の形態の超音波診断装置は、変化量算出部と、推定部と、表示制御部とを備える。変化量算出部は、超音波プローブが送信した超音波の反射波に基づいて時系列に沿って生成された複数の超音波画像間で、複数の局所領域それぞれのパターンの変化量を算出する。推定部は、前記変化量算出部により算出された前記複数の局所領域それぞれのパターンの変化量に基づいて、各超音波画像における走査断面の回転軸の位置を推定する。表示制御部は、前記推定部により推定された前記回転軸を各超音波画像に重畳して所定の表示部に表示するように制御する。
【図面の簡単な説明】
【０００９】
【図１】図１は、実施例１に係る超音波診断装置の構成を説明するための図である。
【図２】図２は、実施例１に係る画像処理部の構成を説明するための図である。
【図３】図３は、局所領域の一例を説明するための図である。
【図４Ａ】図４Ａは、実施例１に係る推定部の処理の一例を説明するための図（１）である。
【図４Ｂ】図４Ｂは、実施例１に係る推定部の処理の一例を説明するための図（２）である。
【図５Ａ】図５Ａは、実施例１に係る推定部の推定結果に基づいて表示される画像の一例を説明するための図（１）である。
【図５Ｂ】図５Ｂは、実施例１に係る推定部の推定結果に基づいて表示される画像の一例を説明するための図（２）である。
【図６】図６は、実施例１に係る超音波診断装置の処理を説明するためのフローチャートである。
【図７】図７は、実施例２に係る推定部の処理の一例を説明するための図である。
【図８】図８は、実施例２に係る超音波診断装置の処理を説明するためのフローチャートである。
【発明を実施するための形態】
【００１０】
以下、添付図面を参照して、超音波診断装置の実施例を詳細に説明する。
【実施例１】
【００１１】
まず、実施例１に係る超音波診断装置の構成について説明する。図１は、実施例１に係る超音波診断装置の構成を説明するための図である。図１に示すように、実施例１に係る超音波診断装置は、超音波プローブ１と、モニタ２と、入力装置３と、装置本体１０とから構成される。
【００１２】
超音波プローブ１は、例えば、一列に配列された複数の圧電振動子を有し、これら複数の圧電振動子は、後述する装置本体１０が有する送信部１１から供給される駆動信号に基づき超音波を発生し、さらに、被検体Ｐからの反射波を受信して電気信号に変換する。また、超音波プローブ１は、圧電振動子に設けられる整合層と、圧電振動子から後方への超音波の伝播を防止するバッキング材などを有する。
【００１３】
超音波プローブ１から被検体Ｐに超音波が送信されると、送信された超音波は、被検体Ｐの体内組織における音響インピーダンスの不連続面で次々と反射され、反射波信号として超音波プローブ１が有する複数の圧電振動子にて受信される。受信される反射波信号の振幅は、超音波が反射される不連続面における音響インピーダンスの差に依存する。なお、送信された超音波パルスが移動している血流や心臓壁などの表面で反射された場合の反射波信号は、ドプラ効果により、移動体の超音波送信方向に対する速度成分に依存して、周波数偏移を受ける。
【００１４】
モニタ２は、超音波診断装置の操作者が入力装置３を用いて各種設定要求を入力するためのＧＵＩ（Graphical  User  Interface）を表示したり、装置本体１０において生成された超音波画像を表示したりする。
【００１５】
入力装置３は、マウス、キーボード、ボタン、パネルスイッチ、タッチコマンドスクリーン、フットスイッチ、トラックボールなどを有し、超音波診断装置の操作者からの各種設定要求を受け付け、装置本体１０に対して受け付けた各種設定要求を転送する。例えば、本実施例１に係る入力装置３は、後述する画像処理部１５による画像処理の開始を操作者から受け付けるための「回転軸表示モードの開始スイッチ」を有する。
【００１６】
装置本体１０は、超音波プローブ１が受信した反射波に基づいて超音波画像を生成する装置であり、図１に示すように、送信部１１と、受信部１２と、Ｂモード処理部１３と、ドプラ処理部１４と、画像処理部１５と、画像メモリ１６と、制御部１７と、内部記憶部１８とを有する。
【００１７】
送信部１１は、トリガ発生回路、遅延回路およびパルサ(pulser)回路などを有し、超音波プローブ１に駆動信号を供給する。パルサ回路は、所定のレート周波数で、送信超音波を形成するためのレートパルスを繰り返し発生する。また、遅延回路は、超音波プローブ１から発生される超音波をビーム状に集束して送信指向性を決定するために用いられる圧電振動子ごとの遅延時間を、パルサ回路が発生する各レートパルスに対し与える。また、トリガ発生回路は、レートパルスに基づくタイミングで、超音波プローブ１に駆動信号（駆動パルス）を印加する。
【００１８】
受信部１２は、アンプ回路、Ａ／Ｄ変換器、加算器などを有し、超音波プローブ１が受信した反射波信号に対して各種処理を行なって反射波データを生成する。アンプ回路は、反射波信号を増幅してゲイン補正処理を行ない、Ａ／Ｄ変換器は、ゲイン補正された反射波信号をＡ／Ｄ変換して受信指向性を決定するのに必要な遅延時間を与え、加算器は、Ａ／Ｄ変換器によって処理された反射波信号の加算処理を行なって反射波データを生成する。加算器の加算処理により、反射波信号の受信指向性に応じた方向からの反射成分が強調される。
【００１９】
このように、送信部１１および受信部１２は、超音波の送受信における送信指向性と受信指向性とを制御する。
【００２０】
Ｂモード処理部１３は、受信部１２から反射波データを受け取り、対数増幅、包絡線検波処理などを行なって、信号強度が輝度の明るさで表現されるデータ（Ｂモードデータ）を生成する。
【００２１】
ドプラ処理部１４は、受信部１２から受け取った反射波データから速度情報を周波数解析し、ドプラ効果による血流や組織、造影剤エコー(echo  signal  from  contrast  agent)成分を抽出し、平均速度、分散、パワーなどの移動体情報を多点について抽出したデータ（ドプラデータ）を生成する。
【００２２】
画像処理部１５は、Ｂモード処理部１３が生成したＢモードデータや、ドプラ処理部１４が生成したドプラデータから、表示用の超音波画像を生成する。具体的には、画像処理部１５は、超音波スキャンの走査線信号列を、テレビなどに代表されるビデオフォーマットの走査線信号列に変換（スキャンコンバート）することで、Ｂモードデータやドプラデータから表示用画像としての超音波画像（Ｂモード画像やドプラ画像）を生成する。すなわち、Ｂモードデータおよびドプラデータは、スキャンコンバート処理前の超音波画像であり、表示用の超音波画像を生成するための超音波画像である。なお、以下では、Ｂモードデータおよびドプラデータのことを生データと記載する場合がある。
【００２３】
また、画像処理部１５は、種々のパラメータの文字情報、目盛り、ボディーマークなどを、生成した超音波画像に合成した合成画像を生成してモニタ２に出力する。また、画像処理部１５は、表示用の超音波画像や生データに対して後に詳述する画像処理を行なう。
【００２４】
画像メモリ１６は、画像処理部１５が生成した超音波画像や、画像処理部１５が超音波画像を画像処理することで生成した画像を記憶するメモリである。また、画像メモリ１６は、Ｂモード処理部１３やドプラ処理部１４が生成した生データを記憶することも可能である。
【００２５】
制御部１７は、超音波診断装置における処理全体を制御する。具体的には、制御部１７は、入力装置３を介して操作者から入力された各種設定要求や、内部記憶部１８から読込んだ各種制御プログラムおよび各種設定情報に基づき、送信部１１、受信部１２、Ｂモード処理部１３、ドプラ処理部１４および画像処理部１５の処理を制御したり、画像メモリ１６が記憶する超音波画像などをモニタ２にて表示するように制御したりする。
【００２６】
内部記憶部１８は、超音波送受信、画像処理および表示処理を行なうための制御プログラムや、診断情報（例えば、患者ＩＤ、医師の所見など）や、診断プロトコルや各種設定情報などの各種データを記憶する。また、内部記憶部１８は、必要に応じて、画像メモリ１６が記憶する画像の保管などにも使用される。なお、内部記憶部１８が記憶するデータは、図示しないインターフェース回路を経由して、外部の周辺装置へ転送することができる。
【００２７】
以上、実施例１に係る超音波診断装置の全体構成について説明した。かかる構成のもと、実施例１に係る超音波診断装置は、超音波プローブ１が送信した超音波の反射波に基づいて超音波画像を生成する。そして、実施例１に係る超音波診断装置の操作者は、断層像である超音波画像を参照し、当該超音波画像にて着目した部位の形状を観察するために、超音波プローブ１を被検体Ｐの体表にて回転させる手技を行なう。かかる場合、実施例１に係る超音波診断装置は、モニタ２にて、回転中に時系列に沿って順次生成された複数の超音波画像を動画表示する。また、画像メモリ１６は、時系列に沿って順次生成された複数の超音波画像の動画像を記憶する。しかしながら、超音波プローブ１を回転することで生成された２次元の超音波画像の動画像を観察することは、操作者にとって、下記のような不便さを伴うことがある。
【００２８】
すなわち、超音波プローブ１を回転する場合、超音波の走査断面は、時系列に沿って変化するために、超音波画像を参照しても、超音波プローブ１をどのように動かしているのか判別することが困難となる。例えば、超音波画像を参照しても、読影者は、超音波プローブ１を回転移動させた場合と、超音波プローブ１を回転させずに垂直方向に移動させた場合とでは共に走査面が変化するために、両者の違いを区別することが困難である。特に、画像メモリ１６に記録された動画像を事後に閲覧する場合や、操作者以外の人間が動画像を参照して読影する場合、超音波画像だけから超音波プローブ１の移動状況を判断することは、困難である。一方、超音波プローブ１を断層面に沿って平行移動させた場合は、動画像上では、超音波画像内のパターンが保持されたまま右か左かへ移動するので、読影者は、超音波プローブ１の移動状況を容易に判断できる。
【００２９】
このため、実施例１に係る超音波診断装置は、以下、詳細に説明する画像処理部１５の処理により、超音波プローブ１の移動状況を容易に判断することを可能とする。
【００３０】
実施例１に係る画像処理部１５が実行する画像処理について、図２などを用いて詳細に説明する。図２は、実施例１に係る画像処理部の構成を説明するための図である。図２に示すように、実施例１に係る画像処理部１５は、画像生成部１５ａと、画像合成部１５ｂと、変化量算出部１５ｃと、推定部１５ｄとを有する。
【００３１】
画像生成部１５ａは、表示用の超音波画像として、ＢモードデータからＢモード画像を生成し、ドプラデータからドプラ画像を生成し、生成した超音波画像を画像メモリ１６に格納する。例えば、画像生成部１５ａは、操作者である医師または技師が被検体Ｐにあてた超音波プローブ１を回転中に、Ｂモード処理部１３が時系列に沿って順次生成した複数のＢモードデータ（生データ）から、時系列に沿った複数のＢモード画像を生成する。
【００３２】
画像合成部１５ｂは、画像生成部１５ａが生成した超音波画像に、種々のパラメータの文字情報、目盛り、ボディーマークなどを合成した合成画像を生成し、生成した合成画像を制御部１７の制御の元、モニタ２に出力する。
【００３３】
変化量算出部１５ｃは、超音波プローブ１が送信した超音波の反射波に基づいて時系列に沿って生成された複数の超音波画像間で、複数の局所領域それぞれのパターンの変化量を算出する。具体的には、変化量算出部１５ｃは、時系列に沿って生成された複数の表示用の超音波画像間、または、複数の生データ間で、複数の局所領域それぞれのパターンの変化量を算出する。例えば、「回転軸表示モードの開始スイッチ」をＯＮとした状態で操作者が超音波画像の撮影を開始すると、変化量算出部１５ｃは、変化量算出処理を開始する。
【００３４】
ここで、局所領域とは、パターンの変化量を監視するために予め設定された関心領域（Region  Of  Interest）のことである。図３は、局所領域の一例を説明するための図である。
【００３５】
例えば、図３に示すように、変化量算出部１５ｃは、超音波画像に設定される座標を用いて、超音波画像に対して横方向に配列された短冊状の９つのＲＯＩ（ＲＯＩ１〜ＲＯＩ９）を複数の局所領域として各超音波画像内に設定する。すなわち、複数の局所領域である複数のＲＯＩそれぞれは、超音波画像における縦方向に沿って設定される。例えば、ＲＯＩは、超音波の走査線に沿って設定される。
【００３６】
ここで、ＲＯＩ１〜ＲＯＩ９ぞれぞれの各超音波画像内における座標は、内部記憶部１８に予め格納されている。あるいは、ＲＯＩ１〜ＲＯＩ９ぞれぞれの各超音波画像内における座標は、入力装置３が有するマウスを用いて操作者が手動で設定する場合であってもよい。また、ＲＯＩの数は、９つに限定されるものではなく、任意に変更可能である。また、図３に示す一例では、各ＲＯＩが所定の間隔を空けて設定されているが、本実施例１は、各ＲＯＩが一部重複して設定される場合であっても適用可能である。
【００３７】
また、図３に示す一例では、超音波画像の深さ方向を略網羅するように、短冊状の各ＲＯＩが設定されている。しかし、本実施例１では、短冊状の各ＲＯＩの縦方向の長さは、任意に変更することが可能である。例えば、短冊状の各ＲＯＩの縦方向の長さは、超音波画像の深さ方向の一部（例えば、超音波画像の縦方向の中心部分のみ）を網羅するように、設定される場合であってもよい。かかる設定により、変化量算出部１５ｃの処理に要する負荷を軽減することができる。
【００３８】
そして、本実施例１に係る変化量算出部１５ｃは、パターンの変化量として、局所領域間の相互相関係数（Ｒ）を算出する。
【００３９】
以下、変化量算出部１５ｃによる相互相関係数の算出処理が、時系列に沿って連続する２つの超音波画像（表示用の超音波画像）である画像「ｉ」と画像「ｉ＋１」との間で実行される場合について説明する。なお、本実施例１は、近接した画像間の違いが小さい場合、何枚か離れた超音波画像間（例えば、画像「ｉ」と画像「ｉ＋ｍ」（ただし、ｍは、２以上の自然数）との間で、変化量算出部１５ｃによる処理が実行される場合であってもよい。かかる「ｍ」の値は、内部記憶部１８に予め格納されている場合であってもよいし、「回転軸表示モードの開始スイッチ」をＯＮとする前に、操作者により入力装置３を介して入力される場合であってもよい。
【００４０】
例えば、変化量算出部１５ｃは、各局所領域の相互相関係数を、以下に示す式（１）により算出する。
【００４１】
【数１】
2012030053.tif 000003
【００４２】
ここで、式（１）に示す「n」は、相互相関係数の算出対象となる局所領域の信号（画素）の数である。また、式（１）に示す「xi」は、画像「ｉ」の局所領域から抽出された信号値（画素値）であり、式（１）に示す「yi」は、画像「ｉ＋１」の局所領域から抽出された信号値（画素値）である。また、式（１）に示す「x」の上に横棒が付与された記号（バーエックス）は、画像「ｉ」の局所領域から抽出された「xi」全ての値から算出された平均値である。また、式（１）に示す「y」の上に横棒が付与された記号（バーワイ）は、画像「ｉ＋１」の局所領域から抽出された「yi」全ての値から算出された平均値である。
【００４３】
回転軸の近傍に位置する領域間のパターンの変化量は、小さくなることが想定される。すなわち、回転軸の近傍に位置する領域間のパターンの変相互相関係数は、大きくなることが想定される。一方、回転軸より遠い場所に位置する領域間のパターン変化量は、回転によって超音波の走査断面が多く変化することから、大きくなることが想定される。すなわち、回転軸より遠い場所に位置する領域間の相互相関係数は、小さくなることが想定される。
【００４４】
このことから、変化量算出部１５ｃにより算出されるパターンの変化量は、回転軸との位置関係に応じて値が変化する値であるならば、相互相関係数以外にも、既に提案されている種々の手法を用いて算出される場合であってもよい。また、本実施例１は、パターン変化量の算出対象が生データである場合であっても適用可能である。
【００４５】
図２に戻って、推定部１５ｄは、変化量算出部１５ｃにより算出された複数の局所領域それぞれのパターンの変化量に基づいて、各超音波画像における走査断面の回転軸の位置を推定する。ここで、走査断面の回転軸は、超音波プローブ１の回転軸となる。したがって、推定部１５ｄは、複数の局所領域それぞれのパターンの変化量に基づいて、各超音波画像における超音波プローブ１の回転軸の位置を推定する。具体的には、推定部１５ｄは、変化量算出部１５ｃにより算出された各局所領域間の相互相関係数に基づいて、超音波プローブ１（走査断面）の回転軸の位置を推定する。ここで、上述したように、相関係数が大きくなっている場所が超音波画像内で回転軸が位置する場所であると推定することができる。したがって、推定部１５ｄは、変化量算出部１５ｃにより算出された相互相関係数が最大となった局所領域の位置を回転軸の位置と推定することが可能である。しかし、本実施例１に係る推定部１５ｄは、回転軸の位置をより正確に推定するために、以下に示す処理を行なう。
【００４６】
図４Ａおよび図４Ｂは、実施例１に係る推定部の処理の一例を説明するための図である。
【００４７】
例えば、推定部１５ｄは、図４Ａに示すように、複数の局所領域（ＲＯＩ１〜９）それぞれの相互相関係数を各局所領域（ＲＯＩ１〜９）の超音波画像内の位置に応じてプロットしたグラフを生成する。
【００４８】
ここで、推定部１５ｄは、図４Ａのグラフにて相互相関係数の最大値となる場所（ＲＯＩ５の横方向の中心位置）を単純に回転軸の位置として推定することも可能である。しかし、相互相関係数の算出結果には、図４Ａのグラフに示すように、雑音成分も含まれており、単純に最大値を用いて回転軸の位置を推定すると、正確さを欠くことも起こりうる。
【００４９】
そこで、推定部１５ｄは、図４Ｂに示すように、例えば、多項式を用いることで、図４Ａのグラフを近似した近似曲線を生成する。すなわち、推定部１５ｄは、グラフの形状を曲線により近似する。そして、推定部１５ｄは、図４Ｂに示すように、生成した近似曲線における相互相関係数が最大となる位置を超音波画像内の回転軸の位置として推定する。すなわち、推定部１５ｄは、グラフの概型から推定されるピーク値を用いて回転軸の位置を推定する。
【００５０】
そして、推定部１５ｄは、推定した回転軸を、回転軸として判断すべきか否かを判定する処理を行なう。すなわち、推定部１５ｄは、超音波プローブ１（走査断面）の回転軸の位置を推定するために用いたパターンの変化量が所定の閾値より大きい場合、超音波プローブ１（走査断面）の回転軸の位置が推定不可であると判定する。相互相関係数をパターン変化量として算出する本実施例１では、推定部１５ｄは、超音波プローブ１の回転軸の位置を推定するために用いた相互相関係数が所定の閾値より小さい場合、超音波プローブ１の回転軸の位置が推定不可であると判定する。
【００５１】
例えば、推定部１５ｄは、算出された相互相関係数や推定された相互相関係数が予め設定された閾値（例えば、０．８）より小さい場合、回転軸の位置が推定不可であると判定する。
【００５２】
図２に戻って、制御部１７は、推定部１５ｄにより推定された回転軸を各超音波画像に重畳してモニタ２に表示するように制御する。具体的には、図２に示す画像合成部１５ｂは、推定部１５ｄの推定結果に基づいて、超音波プローブ１（走査断面）の回転軸が各超音波画像に重畳された合成画像を生成し、制御部１７は、合成画像をモニタ２に表示するように制御する。すなわち、画像合成部１５ｂは、推定部１５ｄの推定結果に基づいて、超音波プローブ１（走査断面）の回転軸が画像「ｉ＋１」に重畳された合成画像を生成する。
【００５３】
さらに、制御部１７は、推定部１５ｄが生成したグラフをモニタ２に表示するように制御する。すなわち、制御部１７は、超音波プローブ１の回転軸が画像「ｉ＋１」に重畳された合成画像とともに、推定部１５ｄが回転軸を推定するために生成したグラフをモニタ２に表示するように制御する。
【００５４】
図５Ａおよび図５Ｂは、実施例１に係る推定部の推定結果に基づいて表示される画像の一例を説明するための図である。例えば、画像合成部１５ｂは、図５Ａに示すように、推定部１５ｄが推定した画像「ｉ＋１」内の超音波プローブ１の回転軸を、点線の細い矩形により表す合成画像を生成する。そして、モニタ２は、制御部１７の制御により、図５Ａに示す合成画像を表示する。ここで、点線の細い矩形により回転軸を描出することで、読影者は、回転軸が実線の直線で示される場合と比較して、超音波画像の視認性を確保したまま、回転軸の位置を把握することができる。なお、超音波画像の視認性を確保したまま、回転軸の位置を把握することができるならば、回転軸の描出法は、任意の方法により行なうことができる。例えば、画像合成部１５ｂは、透明度を上げた矩形を超音波画像に重畳させた合成画像を生成する場合であってもよい。
【００５５】
また、モニタ２は、制御部１７の制御により、図５Ｂに示すグラフを表示する。例えば、図５Ｂに示す一例では、推定部１５ｄが９つ以上のＲＯＩが設定されることで生成したグラフを示している。なお、モニタ２は、制御部１７の制御により、図４Ｂに示すような近似曲線がグラフ上に重畳された画像を表示する場合であってもよい。
【００５６】
また、制御部１７は、推定部１５ｄの推定結果が推定不可である場合、超音波プローブ１の回転軸を非表示とするように制御する。すなわち、制御部１７は、推定部１５ｄの推定結果が推定不可である場合、画像「ｉ＋１」のみを表示するように制御する。
【００５７】
次に、図６を用いて、実施例１に係る超音波診断装置の処理について説明する。図６は、実施例１に係る超音波診断装置の処理を説明するためのフローチャートである。
【００５８】
図６に示すように、実施例１に係る超音波診断装置は、入力装置３が有する「回転軸表示モードの開始スイッチ」を操作者がＯＮにすることで、回転軸表示モードの超音波画像の撮影要求を受け付けたか否かを判定する（ステップＳ１０１）。ここで、回転軸表示モードの超音波画像の撮影要求を受け付けない場合（ステップＳ１０１否定）、超音波診断装置は、待機状態となる。
【００５９】
一方、回転軸表示モードの超音波画像の撮影要求を受け付けた場合（ステップＳ１０１肯定）、制御部１７は、画像生成部１５ａにより超音波画像が生成されたか否かを判定する（ステップＳ１０２）。ここで、超音波画像が生成されていない場合（ステップＳ１０２否定）、制御部１７は、超音波画像が生成されるまで待機する。
【００６０】
一方、超音波画像が生成された場合（ステップＳ１０２肯定）、制御部１７は、生成された超音波画像をモニタ２にて表示するように制御し（ステップＳ１０３）、表示した画像を画像「ｉ」と設定する（ステップＳ１０４）。
【００６１】
そして、制御部１７は、新規の超音波画像が生成されたか否かを判定し（ステップＳ１０５）、新規の超音波画像が生成されていない場合（ステップＳ１０５否定）、新規の超音波画像が生成されるまで待機する。
【００６２】
一方、新規の超音波画像が生成された場合（ステップＳ１０５肯定）、制御部１７は、新規に生成された超音波画像を画像「ｉ＋１」と設定し（ステップＳ１０６）、変化量算出部１５ｃは、画像「ｉ」と画像「ｉ＋１」との間で、複数のＲＯＩにおける相互相関係数を算出する（ステップＳ１０７）。
【００６３】
そして、推定部１５ｄは、複数のＲＯＩそれぞれの相互相関係数を各ＲＯＩの超音波画像内の位置に応じてプロットしたグラフを生成し、生成したグラフの近似曲線から相互相関係数の最大値を推定する（ステップＳ１０８）。
【００６４】
その後、推定部１５ｄは、推定した最大値が閾値以上であるか否かを判定する（ステップＳ１０９）。ここで、推定した最大値が閾値より小さい場合（ステップＳ１０９否定）、推定部１５ｄは、回転軸推定不可と判定し、制御部１７は、画像「ｉ＋１」をモニタ２に表示するように制御する（ステップＳ１１２）。
【００６５】
一方、推定した最大値が閾値以上である場合（ステップＳ１０９肯定）、推定部１５ｄは、回転軸の位置を推定する（ステップＳ１１０）。すなわち、推定部１５ｄは、グラフの形状から推定した最大値の位置に対応する位置を画像「ｉ＋１」内の回転軸の位置として推定する。
【００６６】
その後、制御部１７は、画像合成部１５ｂにより生成された画像「ｉ＋１」と回転軸との合成画像とともに、グラフをモニタ２に表示するように制御する（ステップＳ１１１）。
【００６７】
ステップＳ１１１またはステップＳ１１２の制御処理の後、制御部１７は、操作者から入力装置３を介して撮影終了要求を受け付けたか否かを判定する（ステップＳ１１３）。ここで、撮影終了要求を受け付けない場合（ステップＳ１１３否定）、制御部１７は、画像「ｉ＋１」を画像「ｉ」と設定し（ステップＳ１１４）、ステップＳ１０５に戻って、新規の超音波画像が生成されたか否かを判定する。すなわち、制御部１７は、ステップＳ１１３にて画像「ｉ」として設定した超音波画像と、ステップＳ１０６にて画像「ｉ＋１」として設定する新規画像との間で変化量算出処理が実行されるように制御する。
【００６８】
一方、撮影終了要求を受け付けた場合（ステップＳ１１３肯定）、制御部１７は、処理を終了する。
【００６９】
上述してきたように、実施例１では、変化量算出部１５ｃは、超音波プローブ１が送信した超音波の反射波に基づいて時系列に沿って生成された複数の超音波画像間（複数の生データ間、または、複数の表示用超音波画像間）で、複数の局所領域（ＲＯＩ）それぞれのパターンの変化量を算出する。推定部１５ｄは、変化量算出部１５ｃにより算出された複数の局所領域（ＲＯＩ）それぞれのパターンの変化量に基づいて、各超音波画像における超音波プローブ１（走査断面）の回転軸の位置を推定する。そして、制御部１７は、推定部１５ｄにより推定された回転軸を各超音波画像に重畳してモニタ２に表示するように制御する。なお、複数の局所領域（ＲＯＩ）それぞれは、超音波画像における縦方向に沿って設定される。
【００７０】
したがって、実施例１では、推定した回転軸が重畳された超音波画像を表示するので、超音波画像を参照するだけで、読影者は、超音波プローブ１の移動状況を容易に判別することが可能となる。また、回転軸が表示されることで、読影者は、表示されている画像が超音波プローブ１の回転中に撮影されたものであることを把握できるので、正確な画像診断を行なうことが可能となる。また、操作者が超音波画像の中心を回転軸として正確に超音波プローブ１を回転させることは、熟練が必要である。したがって、上述した回転軸の表示処理を行なうことで、着目する部位を中心に超音波プローブ１が回転されているか否かを把握することができ、操作者の熟練度を向上させるための支援を行なうことが可能となる。
【００７１】
また、実施例１では、推定部１５ｄは、複数の局所領域（ＲＯＩ）それぞれのパターンの変化量を各局所領域（ＲＯＩ）の超音波画像内の位置に応じてプロットしたグラフを生成し、当該生成したグラフに基づいて、走査断面の回転軸である超音波プローブ１の回転軸の位置を推定する。そして、制御部１７は、さらに、推定部１５ｄが生成したグラフをモニタ２に表示するように制御する。したがって、実施例１では、例えば、グラフの近似曲線を用いてパターンの変化量のピークを推定することができ、回転軸の位置を正確に推定することが可能となる。また、実施例１では、回転軸の位置を推定する根拠となったグラフ（または、グラフおよび近似曲線）を表示するので、読影者は、推定された回転軸の位置の確からしさを確認することが可能となる。
【００７２】
また、実施例１では、推定部１５ｄは、回転軸の位置を推定するために用いたパターンの変化量が所定の閾値より大きい場合、超音波プローブ１の回転軸の位置が推定不可であると判定する。そして、制御部１７は、推定部１５ｄの推定結果が推定不可である場合、回転軸を非表示とするように制御する。したがって、実施例１では、信頼度が低い回転軸の位置を読影者に提示することを回避することができる。
【００７３】
また、実施例１では、変化量算出部１５ｃは、パターンの変化量として相互相関係数を算出する。したがって、実施例１では、既知の演算処理を用いて容易に回転軸の位置を推定することが可能となる。
【実施例２】
【００７４】
実施例２では、超音波プローブ１の回転軸とともに、超音波プローブ１の移動状況を示す指標値が推定される場合について、図７を用いて説明する。図７は、実施例２に係る推定部の処理の一例を説明するための図である。
【００７５】
実施例２に係る超音波診断装置は、図１を用いて説明した実施例１に係る超音波診断装置と同様に構成され、実施例２に係る画像処理部１５は、図２を用いて説明した実施例１に係る画像処理部１５と同様に構成される。しかし、実施例２に係る推定部１５ｄは、さらに、以下の推定処理を行なう。すなわち、実施例２に係る推定部１５ｄは、回転軸とともに、グラフにおける複数の局所領域（ＲＯＩ）それぞれのパターンの変化量の変化率に基づいて、走査断面が回転されている速さを推定する。換言すると、実施例２に係る推定部１５ｄは、超音波プローブ１が回転されている速さを推定する。
【００７６】
まず、実施例２に係る推定部１５ｄは、実施例１と同様に、回転軸の位置を推定するためにグラフの近似曲線を生成する。そして、実施例２に係る推定部１５ｄは、実施例１と同様に、生成した近似曲線を用いて相互相関係数の最大値を推定する。そして、実施例２に係る推定部１５ｄは、図７に示すように、近似曲線において推定した最大値付近の曲率（Ｃ）を算出する。そして、実施例２に係る推定部１５ｄは、図７に示すように、算出した曲率から速さ（Ｓ）を推定する。ここで、超音波プローブ１が回転されている速さが大きい場合、相互相関係数は、推定された最大値を中心にして、速やかに小さくなることが想定される。一方、超音波プローブ１が回転されている速さが小さい場合、相互相関係数は、推定された最大値を中心にして、緩やかに小さくなることが推定される。すなわち、近似曲線における最大値近傍の曲率は、超音波プローブ１が回転されている速さに応じて変化する。そこで、例えば、内部記憶部１８に予め「速さ（Ｓ）と曲率（Ｃ）とを対応付けたテーブル（ＬＵＴ：Lookup  Table）」を格納しておくことで、実施例２に係る推定部１５ｄは、算出した曲率から超音波プローブ１が回転されている速さを推定する。
【００７７】
なお、実施例２に係る推定部１５ｄは、実施例１と同様に、推定した最大値と閾値との比較処理を行ない、最大値が閾値より小さい場合、回転軸とともに速さを推定不可とする。
【００７８】
そして、実施例２に係る制御部１７は、さらに、推定部１５ｄが生成した超音波プローブ１が回転されている速さをモニタ２に表示するように制御する。
【００７９】
次に、図８を用いて、実施例２に係る超音波診断装置の処理について説明する。図８は、実施例２に係る超音波診断装置の処理を説明するためのフローチャートである。
【００８０】
図８に示すように、実施例１に係る超音波診断装置は、入力装置３が有する「回転軸表示モードの開始スイッチ」を操作者がＯＮにすることで、回転軸表示モードの超音波画像の撮影要求を受け付けたか否かを判定する（ステップＳ２０１）。ここで、回転軸表示モードの超音波画像の撮影要求を受け付けない場合（ステップＳ２０１否定）、超音波診断装置は、待機状態となる。
【００８１】
一方、回転軸表示モードの超音波画像の撮影要求を受け付けた場合（ステップＳ２０１肯定）、制御部１７は、画像生成部１５ａにより超音波画像が生成されたか否かを判定する（ステップＳ２０２）。ここで、超音波画像が生成されていない場合（ステップＳ２０２否定）、制御部１７は、超音波画像が生成されるまで待機する。
【００８２】
一方、超音波画像が生成された場合（ステップＳ２０２肯定）、制御部１７は、生成された超音波画像をモニタ２にて表示するように制御し（ステップＳ２０３）、表示した画像を画像「ｉ」と設定する（ステップＳ２０４）。
【００８３】
そして、制御部１７は、新規の超音波画像が生成されたか否かを判定し（ステップＳ２０５）、新規の超音波画像が生成されていない場合（ステップＳ２０５否定）、新規の超音波画像が生成されるまで待機する。
【００８４】
一方、新規の超音波画像が生成された場合（ステップＳ２０５肯定）、制御部１７は、新規に生成された超音波画像を画像「ｉ＋１」と設定し（ステップＳ２０６）、変化量算出部１５ｃは、画像「ｉ」と画像「ｉ＋１」との間で、複数のＲＯＩにおける相互相関係数を算出する（ステップＳ２０７）。
【００８５】
そして、推定部１５ｄは、複数のＲＯＩそれぞれの相互相関係数を各ＲＯＩの超音波画像内の位置に応じてプロットしたグラフを生成し、生成したグラフの近似曲線から相互相関係数の最大値を推定する（ステップＳ２０８）。
【００８６】
その後、推定部１５ｄは、推定した最大値が閾値以上であるか否かを判定する（ステップＳ２０９）。ここで、推定した最大値が閾値より小さい場合（ステップＳ２０９否定）、推定部１５ｄは、回転軸および速さを推定不可と判定し、制御部１７は、画像「ｉ＋１」をモニタ２に表示するように制御する（ステップＳ２１２）。
【００８７】
一方、推定した最大値が閾値以上である場合（ステップＳ２０９肯定）、推定部１５ｄは、回転軸の位置および速さを推定する（ステップＳ２１０）。すなわち、推定部１５ｄは、推定した最大値近傍の近似曲線の曲率を算出し、算出した曲率に対応する速さを、超音波プローブ１の回転されている速さとして推定する。
【００８８】
その後、制御部１７は、画像合成部１５ｂにより生成された画像「ｉ＋１」と回転軸との合成画像とともに、グラフおよび速さをモニタ２に表示するように制御する（ステップＳ２１１）。
【００８９】
ステップＳ２１１またはステップＳ２１２の制御処理の後、制御部１７は、操作者から入力装置３を介して撮影終了要求を受け付けたか否かを判定する（ステップＳ２１３）。ここで、撮影終了要求を受け付けない場合（ステップＳ２１３否定）、制御部１７は、画像「ｉ＋１」を画像「ｉ」と設定し（ステップＳ２１４）、ステップＳ２０５に戻って、新規の超音波画像が生成されたか否かを判定する。
【００９０】
一方、撮影終了要求を受け付けた場合（ステップＳ２１３肯定）、制御部１７は、処理を終了する。
【００９１】
上述してきたように、実施例２では、推定部１５ｄは、グラフにおける複数の局所領域それぞれのパターンの変化量の変化率に基づいて、超音波プローブ１（走査断面）が回転されている速さを推定する。そして、制御部１７は、さらに、推定部１５ｄが推定した超音波プローブ１（走査断面）が回転されている速さをモニタ２に表示するように制御する。したがって、実施例２では、推定した速さも表示するので、超音波画像を参照するだけで、読影者は、超音波プローブ１の移動状況をより容易に判別することが可能となる。また、上述した速さの表示処理を行なうことで、例えば、速さが一定でないことなどを把握することができ、操作者の熟練度をより向上させるための支援を行なうことが可能となる。
【００９２】
なお、上述した実施例１および２では、変化量算出部１５ｃおよび推定部１５ｄの処理が、超音波画像の撮影とともにリアルタイムで実行される場合について説明した。しかし、上述した実施例１および２は、超音波画像の撮影が終了した後に、画像メモリ１６に格納された生データまたは表示用超音波画像を用いて、変化量算出部１５ｃおよび推定部１５ｄの処理が実行される場合であってもよい。
【００９３】
また、上述した実施例１および２では、超音波診断装置において超音波画像に対する画像処理が行なわれる場合について説明した。しかし、超音波診断装置とは独立に設置された画像処理装置により、上述した画像処理が行なわれる場合であってもよい。具体的には、図２に示す画像生成部１５ａ以外の処理部の機能および制御部１７の表示制御機能を有する画像処理装置が、超音波診断装置、または各種の医用画像のデータを管理するシステムであるＰＡＣＳ（Picture  Archiving  and  Communication  Systems）のデータベースや、医用画像が添付された電子カルテを管理する電子カルテシステムのデータベースから時系列に沿った複数の超音波画像を受信して上述した画像処理を行なう場合であってもよい。
【００９４】
また、上述した実施例１および２では、被検体Ｐを２次元で走査する超音波プローブ１を操作者が手動で回転することで走査断面を回転させる際に、走査断面の回転軸を超音波画像に重畳表示させる場合について説明した。かかる場合では、超音波プローブ１の回転軸が走査断面の回転軸となる。しかし、走査断面の回転は、操作者の手技ではなく、機械的な制御により実行される場合がある。
【００９５】
走査断面の回転が機械的な制御により実行される一例としては、超音波プローブ１が、ＴＥＥ（transesophageal  echocardiography）プローブである場合が挙げられる。或いは、走査断面の回転が機械的な制御により実行される一例としては、超音波プローブ１が、複数の振動子がマトリックス状に配置された２Ｄプローブである場合が挙げられる。２Ｄプローブは、マトリックス状に配置された複数の振動子により被検体Ｐを３次元で超音波走査することが可能であるとともに、超音波を集束して送信することで、被検体Ｐを２次元で走査することも可能である。
【００９６】
ＴＥＥプローブや２Ｄプローブにより超音波画像を撮影する際に走査断面を回転させた場合、超音波画像における走査断面の回転軸は、例えば、超音波プローブ１に対して行なわれた超音波送信制御に関するデータを用いて推定することが可能である。しかし、超音波送信制御に関するデータは、必ずしも超音波画像の付帯情報として保存されるものではない。したがって、走査断面の回転が機械的な制御により実行される場合であっても、上述した実施例１および２で説明した変化量算出部１５ｃおよび推定部１５ｄの処理を超音波診断装置や画像処理装置で実行させることで、例えば、読影医は、走査断面の回転軸を超音波画像に重畳表示させることができる。その結果、読影医は、超音波プローブの移動状況を容易に判別することができる。
【００９７】
また、上述した実施例１および２で説明した画像処理方法は、あらかじめ用意された画像処理プログラムをパーソナルコンピュータやワークステーションなどのコンピュータである画像処理装置にて実行することによって実現することができる。この画像処理プログラムは、インターネットなどのネットワークを介して配布することができる。また、この画像処理プログラムは、ハードディスク、フレキシブルディスク（ＦＤ）、ＣＤ−ＲＯＭ、ＭＯ、ＤＶＤなどのコンピュータで読み取り可能な記録媒体に記録され、コンピュータである画像処理装置によって記録媒体から読み出されることによって実行することもできる。
【００９８】
以上、説明したとおり、実施例１および２によれば、超音波プローブの移動状況を容易に判別することが可能となる。
【００９９】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【０１００】
１  超音波プローブ
２  モニタ
３  入力装置
１０  装置本体
１１  送信部
１２  受信部
１３  Ｂモード処理部
１４  ドプラ処理部
１５  画像処理部
１５ａ  画像生成部
１５ｂ  画像合成部
１５ｃ  変化量算出部
１５ｄ  推定部
１６  画像メモリ
１７  制御部
１８  内部記憶部
【図１】
2012030053.tif 000004
【図２】
2012030053.tif 000005
【図４Ａ】
2012030053.tif 000006
【図４Ｂ】
2012030053.tif 000007
【図６】
2012030053.tif 000008
【図７】
2012030053.tif 000009
【図８】
2012030053.tif 000010
【図３】
2012030053.tif 000011
【図５Ａ】
2012030053.tif 000012
【図５Ｂ】
2012030053.tif 000013