超音波診断装置
【課題】血管壁の弾性率計測を行なう超音波診断装置において、解析する心拍に対応する適正なＢモード画像を表示することができる超音波診断装置を提供する。
【解決手段】Ｂ／Ｍモード表示において、Ｍモード画像中で心拍を選択する手段を有し、心拍が選択されたら、心拍を中心とする心拍の１４０％の期間内で予め設定された所定位置に応じて、選択された心拍の所定位置のＢモード画像を表示することにより、前記課題を解決する。
【解決手段】Ｂ／Ｍモード表示において、Ｍモード画像中で心拍を選択する手段を有し、心拍が選択されたら、心拍を中心とする心拍の１４０％の期間内で予め設定された所定位置に応じて、選択された心拍の所定位置のＢモード画像を表示することにより、前記課題を解決する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、血管壁の弾性率計測に好適な超音波診断装置に関し、特に、Ｂモード画像とＭモード画像との同時表示において、解析する心拍に対応する適正なＢモード画像を表示することができる超音波診断装置に関する。
【背景技術】
【０００２】
医療分野において、超音波画像を利用した超音波診断装置が実用化されている。
一般に、この種の超音波診断装置は、超音波プローブ（超音波探触子 以下、プローブとする）と、診断装置本体とを有しており、プローブから被検体に向けて超音波を送信し、被検体からの超音波エコーをプローブで受信して、その受信信号を診断装置本体で電気的に処理することにより超音波画像が形成される。
【０００３】
また、超音波を血管や心臓壁等に送信し、その超音波エコーを受信して、受信信号を解析することで、心臓壁や血管壁等の変位量を求め、この変位量から、血管壁や心臓壁（心筋）等の弾性率を計測することも行なわれている。
【０００４】
例えば、特許文献１には、心拍（心臓拍動）に同期して動く対象物に超音波を送受信して超音波エコーの受信信号を得、この受信信号の振幅および位相を用いて対象物の瞬間的な瞬時的な位置を決定して、心拍に基づく血管壁の大きな振幅変位運動をトラッキングすることにより、血管の弾性率を求めることが記載されている。
具体的には、血管壁の順次の位置に基づき、血管壁の微小振動の運動速度波形を求め、血管壁内部の深さ方向に所定間隔でとられた局所ごとのトラッキング軌跡を求め、局所毎の厚みの時間変化を算出することにより血管の弾性率を求めている。
【０００５】
また、特許文献２にも、同様に、心拍に同期して動く対象物に超音波を送受信して得られた超音波エコーの受信信号から、血管等の変位量を求め、この変位量から弾性率を求める超音波診断装置が記載されている。
この超音波診断装置では、血管等の対象物から得た受信信号を用いてＢモード画像およびＭモード画像を形成し、Ｍモード画像の受信信号から手ブレや体動ブレを検出して、このブレが検出されたＭモード画像の受信信号を用いてプローブと被検体との位置変化量を検出して、この検出結果から受信信号の確度を判断し、確度が高いと判断されたＭモード画像の受信信号を用いて、対象物の変位量を求め、この変位量から血管壁等の弾性率を計測している。
【０００６】
ところで、特許文献２に示されるように、Ｂモード画像およびＭモード画像の同じ表示を行なう超音波診断装置においては、通常、画像を静止画にするフリーズボタンが押圧されると、フリーズボタンが押された時点でＭモード画像を静止し、かつ、フリーズボタンが押された時点のＢモード画像を表示する。
そのため、表示されるＢモード画像は、必ずしも、行なう診断などに適した画像であるとは限らない。
【０００７】
これに対して、特許文献３には、Ｂモード画像とＭモード画像とを同時に表示する際に、心電図等の生体組織の活動情報に基づく波形図を表示すると共に、この波形図の特定時相をＭモード画像中に表示し、かつ、この特定時相に対応するＢモードのサムネイル画像を表示する超音波診断装置が記載されている。また、この特許文献３に記載される超音波診断装置では、Ｍモード画像中で時相を任意に選択することにより、この時相に対応するＢモードのサムネイル画像を表示する。
特許文献３に記載される超音波診断装置においては、このような構成を有することにより、弾性率や歪みなどの生体組織の性状に関する情報等を効率よく取得して、超音波診断を内の診断効率を向上することを可能にしている。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開平１０−５２２６号公報
【特許文献２】特開２０１０−２３３９５６号公報
【特許文献３】特開２００４−８３５０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００９】
心拍に同期して動く対象物の変位量から、血管壁の弾性率等を正確に計測するためには、Ｍモード画像から解析に適した心拍を選択すると共に、Ｂモード画像も、選択した心拍すなわち解析する心拍に対応し、かつ、解析にも好適なＢモード画像である必要が有る。
しかしながら、従来の超音波診断装置では、必ずしも、Ｍモード画像から任意の心拍を選択できるとは限らず、さらに、Ｍモード画像で任意の心拍を選択できても、選択された心拍に対応し、かつ、解析にも好適なＢモード画像が表示されるとは限らない。
【００１０】
本発明の目的は、前記従来技術の問題点を解決することにあり、Ｂモード画像とＭモード画像とを表示して、Ｍモード画像の中から血管弾性率の計測を行なうための解析に適した心拍を選択ことができ、かつ、Ｍモード画像中で選択された心拍に対応し、解析にも好適なＢモード画像を表示することができる超音波診断装置を提供することにある。
【課題を解決するための手段】
【００１１】
前記目的を達成するために、本発明の超音波診断装置は、超音波を送信し、被検体によって反射された超音波エコーを受信して、受信した超音波エコーに応じた受信信号を出力する超音波トランスデューサを有する超音波プローブと、前記超音波トランスデューサが出力した受信信号から、Ｂモード画像およびＭモード画像を形成する、画像形成手段と、表示手段と、前記画像形成手段が形成したＢモード画像およびＭモード画像の少なくとも一方を、前記表示手段に表示させる表示処理手段と、前記Ｂモード画像およびＭモード画像が表示された状態において、表示画像を静止画にするフリーズ手段と、前記Ｂモード画像およびＭモード画像が静止画にされた状態において、Ｍモード画像中で心拍を選択する心拍選択手段とを有し、かつ、前記表示処理手段は、前記Ｍモード画像中で心拍が選択されたら、心拍の長さに対して心拍の前２０％から心拍の後２０％までの期間内で予め設定された所定位置に対応して、前記選択された心拍における前記所定位置のＢモード画像を前記表示手段に表示させることを特徴とする超音波診断装置を提供する。
【００１２】
このような本発明の超音波診断装置において、前記Ｍモード画像中の心拍を検出する拍検出手段を有し、前記表示処理手段は、前記拍検出手段が心拍を検出したら、Ｍモード画像中に、前記拍検出手段が検出した全ての心拍を表示するのが好ましい。
また、前記拍検出手段は、前記Ｍモード画像を解析することにより、前記Ｍモード画像中の心拍を検出するのが好ましい。
【００１３】
また、前記所定位置を選択する選択手段を有するのが好ましい。
また、前記心拍選択手段は、前記拍検出手段が検出した心拍の内、開始から終了までの全てがＭモード画像中で取得されている、最も新しい心拍が選択されたとするのが好ましい。
また、前記所定位置は、心臓の拡張期よりも後で、かつ、心臓の収縮期よりも前の期間以外に設定されるのが好ましい。
【００１４】
前記心拍選択手段によって選択された心拍の位置を調整する位置調整手段を有するのが好ましい。
また、前記位置調整手段によって心拍の位置調整が行なわれたら、前記表示処理手段は、この心拍の位置調整および前記所定位置に応じて、前記表示手段に表示するＢモード画像を対応する位置の画像に変更するのが好ましい。
また、前記Ｂモード画像中において血管壁の境界位置を設定する、血管壁境界設定手段を有するのが好ましい。
【００１５】
また、表示手段に表示されたＢモード画像中で、関心領域を設定するＲＯＩ設定手段を有するのが好ましい。
さらに、前記関心領域の設定指示に応じて、前記超音波トランスデューサによる超音波のフレームレートを、設定指示前よりも高くするのが好ましい。
【発明の効果】
【００１６】
上記構成を有する本発明の超音波診断装置は、Ｂモード画像とＭモード画像とを同時に表示して、血管弾性の計測等を行なうに際し、Ｍモード画像中から目的とする計測を行なうための解析に適した心拍を選択できると共に、心拍に対する所定期間内で予め設定された所定位置に応じて、選択された心拍における所定位置のＢモード画像を表示する。
そのため、本発明の超音波診断装置では、血管弾性率の計測等に適した心拍を選択した上で、この選択した心拍に対応する、解析に適したＢモード画像を表示ずることができ、正確な血管弾性率の計測等を、安定して行なうことができる。
【図面の簡単な説明】
【００１７】
【図１】本発明の超音波診断装置の一例を概念的に示す図である。
【図２】図１に示す超音波診断装置の構成を概念的に示すブロック図である。
【図３】図１に示す超音波診断装置における血管壁の弾性計測の一例を説明するためのフローチャートである。
【図４】血管壁の弾性計測のための超音波診断を説明するための概念図である。
【図５】（Ａ）および（Ｂ）は、図１に示す超音波診断装置での画像表示の一例を示す概念図である。
【図６】（Ａ）および（Ｂ）は、図１に示す超音波診断装置での画像表示の一例を示す概念図、（Ｃ）は、本発明の超音波診断装置の作用を説明するための概念図である。
【図７】（Ａ）〜（Ｃ）は、図１に示す超音波診断装置での画像表示の一例を示す概念図である。
【図８】（Ａ）および（Ｂ）は、図１に示す超音波診断装置での画像表示の一例を示す概念図である。
【図９】図１に示す超音波診断装置での画像表示の一例を示す概念図である。
【図１０】（Ａ）〜（Ｇ）は、図１に示す超音波診断装置での画像表示の一例を示す概念図である。
【図１１】（Ａ）および（Ｂ）は、図１に示す超音波診断装置での画像表示の一例を示す概念図である。
【図１２】図１に示す超音波診断装置での画像表示の一例を示す概念図である。
【発明を実施するための形態】
【００１８】
以下、本発明の超音波診断装置について、添付の図面に示される好適実施例を基に、詳細に説明する。
【００１９】
図１に、本発明の超音波診断装置の一例の外観を概念的に示す。
図１に示すように、超音波診断装置１０は、基本的に、診断装置本体１２と、超音波プローブ１４と、操作パネル１６と、ディスプレイ１８とを有して構成される。また、超音波診断装置１０の下端部には、キャスタ２４が配置されており、人力で容易に装置を移動することが可能になっている。
【００２０】
超音波プローブ１４（以下、プローブ１４とする）は、超音波の送受信を行なって、受信した超音波エコーに応じた受信信号を診断装置本体１０に供給するものである。
このプローブ１４は、被検体に超音波を送信し、被検体に反射された超音波エコーを受信して、受信した超音波エコーに応じた電気信号（受信信号）を出力する、いわゆる超音波トランスデューサ（超音波圧電素子）を１次元的もしくは二次元的に配列してなる、各種の超音波診断装置に利用される公知の超音波プローブである。
【００２１】
本発明において、プローブ１４の種類には、特に限定はなく、コンベックス型、リニア型、セクタ型等の各種の形式が利用可能である。また、体外式プローブでもよいし、ラジアルスキャン方式等の超音波内視鏡用プローブでもよい。さらに、プローブ１４は、ハーモニックイメージングに対応する、送信した超音波の二次以上の高調波を受信するための超音波振動子を有するものであってもよい。
また、図示例においては、プローブ１４と診断装置本体１２とは、ケーブル２０によって接続されている。しかしながら、本発明は、これに限定はされず、プローブ１４内に後述する送信回路２８、受信回路３０、送受信制御部３２等を配置して、無線通信によってプローブ１４と診断装置本体１２とを接続するものであってもよい。
【００２２】
ディスプレイ１８は、公知のディスプレイ（表示装置）である。
超音波診断装置１０において、ディスプレイ１８は、各種の超音波診断装置と同様、プローブ１４が出力した受信信号に応じた超音波画像、被検者の情報、ＧＵＩ（Graphical User Interface)による操作を行なうための選択手段や指示手段、関心領域（Region of Interest 以下、ＲＯＩとする）、後述する血管壁の弾性計測結果等を表示する。
【００２３】
操作パネル１６は、超音波診断装置１０の操作を行なうものである。
図示は省略するが、超音波診断装置１０において、操作パネル１６には、ＢモードやＭモードなどの各種のモードの選択手段、ディスプレイ１８に表示されたカーソルやライン等を移動するためのトラックボール（トラックパッド／タッチパッド）、選択や操作を決定（確定）するためのセットボタン、動画表示と静止画表示との切り換え等を行なうためのフリーズボタン、超音波画像の視野深度の変更手段、ゲイン調整手段、超音波画像を拡大するためのズームボタン等が配置される。
なお、超音波診断装置１０においては、モードとして、ＢモードやＭモードなどの通常の超音波診断装置が有するモードに加え、血管壁の弾性率を計測するためのモードであるＶＥモード（Vascular Elasticity モード）も設定されている。
また、同じく図示は省略するが、操作パネル１６には、ＧＵＩによる操作等を行なうための表示装置であるタッチパネル１６ａも配置される（図６（Ｂ）参照）。
【００２４】
診断装置本体１２は、超音波診断装置１０の全体の動作の制御を行なうと共に、プローブ１４が出力した受信信号に応じた超音波画像を形成してディスプレイ１８に表示させ、さらに、血管弾性率の計測のための各種の処理を行なうものである。
診断装置本体１２は、例えば、コンピュータを利用して構成される。
【００２５】
図２に、超音波診断装置１０の構成をブロック図で概念的に示す。
図２に示すように、診断装置本体１２は、送信回路２８、受信回路３０、送受信制御部３２、画像形成部３４、記憶部３６、境界検出部４０、トラッキング部４２、拍検出部４６、拍選択部４８、弾性率算出部５０、および、表示処理部５２を有する。
また、画像形成部３４は、Ｂモード画像形成部５６およびＭモード画像形成部５８を有する。
【００２６】
前述のプローブ１４は、送信回路２８および受信回路３０に接続される。また、送信回路２８および受信回路３０には、送受信制御部３２が接続される。さらに、受信回路３０は、画像形成部３４に接続される。
画像形成部３４は、表示処理部５２に接続される。また、画像形成部３４のＢモード画像形成部５６およびＭモード画像形成部５８は、記憶部３６に接続される。Ｂモード画像形成部５８は、さらに境界検出部４０にも接続される。
記憶部３６は、トラッキング部４２、拍検出部４６、および表示処理部５２に接続される。また、拍検出部４６および表示処理部５２は、共に、トラッキング部４２および表示処理部５２に接続される。トラッキング部４２は、さらに、弾性率算出部５０にも接続され、この弾性率算出部５０は、表示処理部５２に接続される。また、拍検出部４６は、さらに、拍選択部４８にも接続され、この拍選択部４８も、同じく表示処理部５２に接続される。
【００２７】
送受信制御部３２は、送信回路２８および受信回路３０を介してプローブ１４の超音波ビームの送信方向および超音波エコーの受信方向を、順次、設定する。
また、送受信制御部３２は、設定した送信方向に応じて送信遅延パターンを選択する送信制御機能と、設定した受信方向に応じて受信遅延パターンを選択する受信制御機能とを有している。
【００２８】
送信遅延パターンとは、プローブ１４の複数の超音波トランスデューサから送信される超音波によって所望の方向に超音波ビームを形成するために各超音波トランスデューサの駆動信号に与えられる遅延時間のパターンである。他方、受信遅延パターンとは、複数の超音波トランスデューサによって受信される超音波によって所望の方向からの超音波エコーを抽出するために受信信号に与えられる遅延時間のパターンである。
複数の送信遅延パターンおよび複数の受信遅延パターンが内部メモリ（図示せず）に格納されていて、状況に応じて、適宜、選択して使用する。
【００２９】
送信回路２８は、複数のチャネルを備えており、プローブ１４の複数の超音波トランスデューサにそれぞれ印加する複数の駆動信号を形成する。その際に、送受信制御部３２によって選択された送信遅延パターンに基づいて、複数の駆動信号にそれぞれの遅延時間を与えることができる。
なお、送信回路２８は、複数の超音波トランスデューサから送信される超音波が超音波ビームを形成するように、複数の駆動信号の遅延量を調節して複数の駆動信号をそれぞれプローブ１４の複数の超音波トランスデューサに供給するようにしても良く、複数の超音波トランスデューサから一度に送信される超音波が被検体の撮像領域全体に届くように構成した複数の駆動信号をプローブ１４に供給するようにしても良い。
【００３０】
受信回路３０は、送信回路２８と同様に複数のチャネルを備えており、複数の超音波トランスデューサを介して受信された複数のアナログ信号を増幅し、デジタルの受信信号に変換する。
さらに、送受信制御部３２によって選択された受信遅延パターンに基づいて、複数の受信信号にそれぞれの遅延時間を与え、それらの受信信号を加算することにより、受信フォーカス処理を行う。この受信フォーカス処理によって、超音波エコーの焦点が絞り込まれた音線信号（音線データ）が形成される。
【００３１】
形成された音線データは、画像形成部３４に供給される。
画像形成部３４は、供給された音線データに対して、Ｌｏｇ（対数）圧縮やゲイン調整等のプリプロセス処理を施して超音波画像の画像データを形成し、この画像データを通常のテレビジョン信号の走査方式に従う画像データに変換（ラスター変換）し、さらに階調処理等の必要な画像処理を施した上で表示処理部５２へ出力する。
また、画像形成部３４は、Ｂモード画像を形成するＢモード画像形成部５６と、Ｍモード画像を形成するＭモード画像形成部５８とを有する。Ｂモード画像およびＭモード画像の形成は、公知の方法によればよい。
【００３２】
表示処理部５２は、画像形成部３４から供給された超音波画像の画像データ、記憶部３６から読み出した超音波画像の画像データ、操作パネル１６で行なわれた操作（入力指示）、後述する心拍の検出結果や心拍の選択、血管壁弾性率の計測結果（解析結果）等に応じて、ディスプレイ１８に表示するための表示用データを形成し、ディスプレイ１８に表示させる部位である。
【００３３】
図示例の超音波診断装置１０において、診断装置本体１２の記憶部３６、境界検出部４０、トラッキング部４２、拍検出部４６、拍選択部４８、および、弾性率算出部５０は、主に、血管壁の弾性率を計測するＶＥモードの際に用いられる部位である。
以下、図３のフローチャートおよび図５〜図１２を参照して、ＶＥモードにおける超音波診断装置１０の作用を説明することにより、上記記憶部３６や境界検出部４０等の各部位、ならびに、本発明の超音波診断装置１０について、より詳細に説明する。
なお、以下の説明では、特に記載がなくても、ディスプレイ１８の表示に関しては、表示処理部５２がラインの形成等の必要な処理を行なう。
【００３４】
超音波診断装置１０による超音波診断が開始されると、送受信制御部３２による制御の下、送信回路２８がプローブ１４の超音波トランスデューサから超音波を送信させ、また、受信回路３０はプローブ１４が出力した受信信号を処理して音線信号を形成し、画像形成部３４に出力する。
一例として、Ｂモードが選択され、図４に概念的に示すように、被検者の頸動脈ｃを測定対象として、プローブ１４が首ｎに当てられたとして、画像形成部３４（Ｂモード画像形成部５６）によって形成されたＢモード画像が、表示処理部５２で処理されて、ディスプレイ１８に表示される。
【００３５】
目的とする頸動脈ｃが適正に観察できるようになり、操作パネル１６のモード選択手段（以下の説明では、「操作パネル１６」は省略する）によってＶＥモードが選択されると、表示処理部５２は、図５（Ａ）に概念的に示すように、Ｂモード画像中に、関心領域を示すＲＯＩ６０を表示させる。
【００３６】
この状態では、トラックボールでの操作によってＢモード画像中のＲＯＩ６０の位置を移動できる。また、セットボタンを押すと、ＲＯＩ６０の位置が固定され、トラックボールでの操作によってＲＯＩ６０のサイズを変更できる。
さらに、セットボタンを押すたびに、ＲＯＩ６０の位置変更およびＲＯＩ６０のサイズ調整が、交互に実施可能になる。
【００３７】
この状態からズームボタンが押されると（押下されると）、ＲＯＩ６０の位置やサイズの調整が終了してＲＯＩ６０の設定が指示されたとして、送受信制御部３２は、フレームレートをＲＯＩ６０の設定指示前よりも高くする（例えば、２００Ｈｚ以上あるいはＲＯＩ設定指示前の５倍以上）。また、このズームボタンの押下に応じて、Ｍモード画像形成部５８によって、ＲＯＩ６０のＭモード画像の形成が開始され、図５（Ｂ）に示されるように、ＲＯＩ６０の部分が拡大されたＢモード画像６４、および、ＲＯＩ６０（その選択ライン６２）のＭモード画像６５が、同時に表示される。
なお、Ｂモード画像６４とＭモード画像６５との同時表示（デュアルモード表示）は、公知の超音波診断装置における、いわゆるＢ／Ｍモード表示と同様に行なえば良い。
【００３８】
図５（Ｂ）において、上側がＢモード画像６４で、下側がＭモード画像６５である。
Ｂモード画像６４において、図中横方向はアジマス方向（超音波振動子の配列方向（二次元配列の場合は長手方向））であり、縦方向は深度方向（超音波の送受信方向）で、上方が深度が浅い側（プローブ１４側）である。
また、Ｂモード画像中には、Ｂモード画像中のアジマス方向におけるＭモード画像表示位置（Ｍモード画像の表示ライン）を選択するための、深度方向に延在する選択ライン６２が表示される。この選択ライン６２は、トラックボールによって、アジマス方向（図中左右方向）に移動可能になっている。
【００３９】
また、Ｍモード画像６５において、横方向は時間軸で、時間は左から右に流れており、間隙６５ａの左側が、現在のフレームとなる（すなわち、間隙６５ａの右側は過去のフレーム）。また、Ｂモード画像６４と同様、縦方向は深度方向で、上方が深度が浅い側である。
図５（Ｂ）において、ディスプレイ１８に表示されているＭモード画像６５は、予め位置が設定された選択ライン６２の位置のＭモード画像６５になっている。
【００４０】
ここで、Ｍモード画像形成部５８は、アジマス方向の所定位置（予め設定された所定位置や選択された位置）や、アジマス方向で選択された位置のみではなく、Ｂモード画像６４のアジマス方向の全域について、Ｍモード画像を形成する。
【００４１】
Ｂモード画像形成部５６が形成したＲＯＩ６０のＢモード画像（Ｂモード画像データ）、および、Ｍモード画像形成部５８が形成したＭモード画像（Ｍモード画像データ）は、共に、記憶部３６に記憶される。
なお、記憶部３６が記憶する画像の時間的な量には、特に限定はないが、一般的な心拍が２以上、入る長さであるのが好ましい。従って、記憶部３６は、最新の３秒以上のＢモード画像およびＭモード画像を記憶するのが好ましい。
【００４２】
前述のように、選択ライン６２は、トラックボールによってアジマス方向に移動することができる。
選択ライン６２の位置とＭモード画像とは、連動している。すなわち、トラックボールによって選択ライン６２を左右方向に移動すると、表示処理部５２は、選択ライン６２の位置のＭモード画像をディスプレイ１８に表示する。
【００４３】
操作者が、適正な画像が得られたと判断したら、フリーズボタンが押される。
フリーズボタンが押されると、表示処理部５２は、必要な画像データを記憶部３６から読み出し、図６（Ａ）に示すように、ディスプレイ１８に、フリーズボタンが押された時点が最も右となるように並べ直して、選択ライン６２の位置のＭモード画像６５の静止画を表示させ、また、後述する所定位置のＢモード画像６４の静止画を表示させる。同時に、選択ライン６２が破線になって、移動ができなくなる（非アクティブになる）。
さらに、図６（Ｂ）に示すように、操作パネル１６のタッチパネル１６ａに、後述する血管壁の境界の設定を指示するための「ＡＷ Ｄｅｔ」ボタン、血管壁弾性率の解析開始を指示するため「Elasticity Ana」ボタン、被検者の血圧を入力するための「Ｐｓ」ボタンおよび「Ｐｄ」ボタン、信頼性閾値を入力するための「Quality Factor Threshold」ボタンが表示される。なお、この時点では、「Elasticity Ana」ボタンは、選択できない状態になっている。
【００４４】
また、フリーズボタンが押されると、拍検出部４６が、記憶部３６に記憶されている全てのＭモード画像に対して、心拍の検出（心拍の自動検出）を行なう。心拍の検出結果は、記憶部３６に送られ、対応するＭモード画像に、情報として付加される。
さらに、心拍の検出結果は、表示処理部５２および拍選択部４８にも送られる。
【００４５】
心拍の検出の検出方法には、特に限定はないが、一例として、Ｍモード画像を解析（例えば時間方向にトラッキング）して、横方向の延在する白線（輝線）の深度方向の移動速度（速度の上昇開始時点）や、同白線の深度方向の動きの脈動等を用いて、検出すればよい。あるいは、心拍の検出に心電計（心電図）を利用してもよい。
また、心拍を自動検出するのではなく、操作者がＭモード画像を見て心拍の位置（心拍の開始と終了の位置）を入力するようにしてもよい。
【００４６】
表示処理部５２は、図６（Ａ）に示すように、心拍（心拍の開始／終了位置）の検出結果を、Ｍモード画像６５中に三角マークおよび直線で表示させる。
なお、検出に失敗した心拍が有る場合には、周囲の心拍の間隔等に応じて、適当な位置に心拍を表示してもよい。
また、拍選択部４８は、検出部４６による心拍の検出結果に応じて、心拍の開始から終了までが含まれる完全な心拍（撮り切れた心拍）の内、最も新しい心拍が選択されたとして、その情報を表示処理部５２に送る。すなわち、拍選択部４８は、２本の直線で挟まれた撮り切れている完全な心拍で、かつ、フリーズボタンが押される直前の心拍が選択されたとして、その情報を表示処理部５２に送る。
さらに、拍検出部４８は、記憶部３６に記憶される全てのＭモード画像に対して、同じ心拍が選択されたとする情報を付す。
【００４７】
表示処理部５２は、拍選択部４８による心拍の選択結果に応じて、一例として、Ｍモード画像６５において、選択された心拍を実線で、それ以外の心拍を破線で示す。図６（Ａ）に示す例では、３つの完全な心拍が検出されているので、前述のように、最新の完全な心拍である一番右の心拍が選択されて、この心拍の開始位置および終了位置（開始時点および終了時点）が実線で示され、それ以外の心拍の開始／終了位置は、破線で示される。
なお、このような選択および非選択の区別は、線種に変えて、あるいは加えて、線の色を利用してもよい。
【００４８】
なお、超音波診断装置１０においては、この後、操作者が心拍を選択（選択する心拍を変更）することが可能であり、この工程で確定された心拍が、最終的に選択された心拍として、後の血管弾性計測のための解析に供される。従って、この拍検出部４８による心拍の自動選択は、いわば仮の心拍の選択である。
また、拍検出部４８が自動選択する心拍は、最新の完全な心拍に限定はされず、その前の心拍であってもよく、検出された最も古い心拍であってもよい。さらに、超音波診断装置１０は、拍検出部４８が自動選択する心拍を操作者が選択するための、選択手段を有しても選択できるようにしてもよい。拍検出部４８が自動選択する心拍の選択手段は、ＧＵＩ等を用いた公知の方法で構成すればよい。
【００４９】
また、フリーズが押されると、Ｂモード画像６４は、Ｍモード画像６５で選択された心拍における、予め設定された所定位置（所定位置の時点（所定の時相））の画像となる。
【００５０】
具体的には、超音波診断装置１０においては、心拍の長さ（一拍の時間）に対して、心拍の前２０％の時点から心拍の後２０％の時点までの期間内で、Ｂモード画像６４を表示する所定位置（いわば、心拍（Ｍモード画像）の時間軸上における所定位置）が、予め設定されている。すなわち、図６（Ｃ）に血管の後壁境界のＭモード画像を模して概念的に示すように、心拍の長さをｔとすると、心拍の前０．２ｔの時点から心拍の後０．２ｔの時点までの期間Ｔ（長さｔの心拍を中心とする１．４ｔ）の中から、心拍開始位置や心拍中央のような所定位置が、予め設定されている。
なお、図６（Ｃ）では、図６（Ａ）等と同様、時間の進行方向は左から右であり、上方が深度が浅く、下方が深い。
【００５１】
これに応じて、超音波診断装置１０では、ＲＯＩが設定されてフリーズボタンが押され、さらに、心拍の検出および選択が行なわれたら、選択された心拍における、心拍の開始位置や心拍の中央位置などの、予め設定された所定位置のＢモード画像（この所定位置の時点で撮られたＢモード画像）を表示する。
【００５２】
図示例においては、一例として、心拍の開始位置が所定位置として設定されている。
従って、この時点では、表示処理部５２は、拍選択部４８によって自動選択された、最新の心拍の開始位置のＢモード画像（自動選択された心拍の開始時点に撮られたＢモード画像）を記憶部３６から読み出して、Ｍモード画像６５と共にディスプレイに表示する。
【００５３】
前述のように、通常の超音波診断装置では、Ｂモード画像とＭモード画像との同時表示において、フリーズボタンが押されて表示されるＢモード画像は、フリーズボタンが押された時点の画像である。
他方、超音波診断措置において、血管壁の変位量から血管弾性率を計測する場合のように、心拍に応じた対象部位の変動を用いる計測では、正確な計測を行なうためには、Ｍモード画像から解析に適正な心拍を選択し、かつ、選択した心拍に対応するＢモード画像を表示して、解析を行なう必要がある。
しかしながら、フリーズ時点のＢモード画像を表示する従来の装置では、必ずしも解析に用いる心拍に対応するＢモード画像が表示されるとは、限らない。また、特許文献３に示されるように、心電図等による特定時相のＢモード画像を表示する装置も提案されているが、やはり、Ｍモード画像中で選択された心拍に対応するとは、限らない。
【００５４】
これに対し、本発明の超音波診断装置１０においては、Ｍモード画像中において心拍を選択し、心拍の長さに対して、心拍の前２０％の時点から心拍の後２０％の時点（以下、単純に、『心拍の前２０％から心拍の後２０％』とも言う）までの期間内で、適宜、設定された所定位置に対応して、選択された心拍の所定位置のＢモード画像を表示する。
そのため、本発明によれば、Ｍモード画像中で選択された心拍に対応するＢモード画像を表示して、好適な解析や診断を行なうことができる。例えば、後述するように、Ｂモード画像を用いて血管壁の境界を設定し、Ｂモード画像の時相および設定した血管壁境界を開始点として、Ｍモード画像中で血管壁のトラッキング（追跡）等の解析を行なう際に、余分な時間分のトラッキングを行なうことなく、正確なトラッキングが可能となる。また、選択した心拍における心臓の拡張期や収縮期に対応するＢモード画像など、解析に好適なＢモード画像を表示することもできる。
【００５５】
本発明において、表示するＢモード画像６４が、選択された心拍の前２０％の時点よりも前、および、心拍の後２０％の時点よりも後では、表示するＢモード画像６４が、Ｍモード画像６５で選択された心拍と時間的に離れすぎてしまう。
その結果、選択された心拍に対して、表示されたＢモード画像６４の状態が大きく異なってしまう場合が有り、適正な解析や診断を行なうことが困難になってしまう。また、例えば、後述するように、Ｂモード画像６４を位置および時間的な開示位置としてＭモード画像中でトラッキングを行なう場合には、トラッキングの開始位置から選択された心拍に至るまでの距離（時間）が長くなり過ぎてしまう。その結果、トラッキングにエラーが生じ易くなってしまい、また、トラッキングの結果にスペックルに起因するノイズ等の不要な情報が多く拾われてしまう等の不都合が生じる。
【００５６】
Ｂモード画像６４を表示する所定位置は、心拍の前２０％から心拍の後２０％までの期間の中であれば、心拍の前（心拍開始前）１０％の位置、心拍の開始位置、心臓の収縮期すなわち血管最大径に対応する位置、心拍の中央、心拍の終了位置など、各種の位置が利用可能である。
ここで、心臓の拡張期ｄから収縮期ｓに至る期間は、図６（Ｃ）にも示されるように、血管壁の動きが激しい（血管壁の移動速度が速い）。そのため、この間のＢモード画像は、他の位置に比して、画質が良くない。また、後述するような血管壁のトラッキングを行なう場合に、このように動きが激しい場所を開始点とすると、トラッキングのエラーが発生し易い。
そのため、所定位置は、この心臓の拡張期ｄすなわち心拍の開始位置よりも後で、かつ、収縮期ｓよりも前の期間以外で設定するのが好ましい。
【００５７】
また、図６（Ｃ）にも示されるように、心臓の拡張期の前すなわち心拍の開始位置の前は、血管壁の移動は緩やかで、画質の良いＢモード画像が得られる。加えて、血管壁の弾性の計測など、血管壁の移動量を用いた計測を行なう場合には、多くの場合、解析に重要なのは、血管径が最小から最大に変化する心臓の拡張期から収縮期までの間である。
従って、Ｂモード画像６４を表示する所定位置は、特に、心拍の開始位置から、前２０％までの期間に設定するのが好ましい。すなわち、図６（Ｃ）に示す、心拍（心拍の開始位置を含む）の前０．２ｔの期間に設定するのが好ましい。
【００５８】
本発明の超音波診断装置１０において、Ｂモード画像６４を表示する所定位置は、デフォルトとして予め固定的に設定されていてもよく、もしくは、医師等の操作者が所定位置を選択して設定できるようにしてもよい。
また、操作者が所定位置を選択して設定する場合には、『心拍開始前１０％の位置』、『心拍開始位置』、『心臓の収縮期の位置』、『心拍の中央』、『心拍終了位置』のように、複数の位置を選択肢として設定しておき、操作者が選択して所定位置を設定できるようにしてもよい。あるいは、心拍の前２０％から後２０％までの期間（時間軸）において、任意の位置を選択して所定位置として設定できるようにしてもよい。あるいは、設定された選択肢から選択する方法と、時間軸上の任意の位置を選択する方法とを、選択できるようにしてもよい。
【００５９】
なお、本発明において、所定位置は１つに限定はされず、２以上の所定位置を設定してもよい。
例えば、血管壁の弾性率を簡易的に計測する方法として、１つの心拍における血管の最大径（Ｄｓ）と、最小径（Ｄｄ）との比（Ｄｄ／Ｄｓ）を算出する方法が知られている。これに対応して、血管が最小径となる心臓拡張期（心拍開始位置）と、血管が最大径になる心臓収縮期を所定位置として設定してもよい。
このように、複数の所定位置を設定する場合には、複数のＢモード画像を切り換えて表示するようにすればよい。あるいは、ディスプレイ１８に十分な表示スペースが有る場合には、Ｍモード画像と共に複数のＢモード画像を同時表示してもよい。
【００６０】
Ｍモード画像６５に心拍のラインが表示され、Ｂモード画像６５が選択された心拍に対応する所定位置の画像（前述のように、図示例では心拍開始位置の画像）となると、Ｂモード画像中の選択ライン６２が実線になって、トラックボールによって左右方向に移動可能となる。すなわち、選択ライン６２が、アクティブな状態になる。同時に、Ｍモード画像の心拍を示すラインは、全て、破線になる。なお、各種のラインがアクティブか否かの区別は、先と同様、線種に変えて、あるいは加えて、線の色を利用してもよい。
この状態で、トラックボールによって選択ライン６２を左右方向に移動すると、表示処理部５２は、記憶部３６から、選択ライン６２の位置に対応するＭモード画像を読み出し、心拍の検出結果と共に、その画像をディスプレイ１８に表示させる。すなわち、フリーズ後にも、トラックボールによって選択ライン６２を移動することで、Ｂモード画像６４中のアジマス方向の全域から、Ｂモード画像６４中におけるＭモード画像６５の表示位置（表示ライン）を選択できる。
従って、本例によれば、設定したＲＯＩ６０のアジマス方向の任意の位置のＭモード画像６５を表示して、Ｍモード画像６５、および、Ｍモード画像中の各心拍に対応する画像を、観察／確認することができる。
【００６１】
Ｂモード画像６４の選択ライン６２が移動可能な状態でセットボタンが押されると、Ｍモード画像の表示位置（表示ライン）の選択が終了したとして、図７（Ａ）に示されるように、Ｂモード画像６４の選択ライン６２が破線になって、トラックボールによる移動が不可能な状態となる。同時に、Ｍモード画像６５において、拍選択部４８によって自動選択された最新の心拍を示すラインが、共に実線となる。
【００６２】
Ｍモード画像６５において、最新の心拍を示すラインが、共に実線になると、トラックボールによって、心拍の選択が可能な状態になる。
セットボタンが押された時点では、前述のように最新の心拍が選択（拍検出部４８によって自動剣術）された状態になっている。これに応じて、Ｍモード画像６５の表示では、図７（Ａ）および（Ｂ）に示すように、最新の心拍を示すラインが実線となって選択された状態となっている。
この状態から、例えば、トラックボールが左に回されると、その信号が拍選択部４８に供給される。この回転量に応じて、拍選択部４８は、最新の心拍の１つ前の心拍が選択されたとして、その旨の情報を表示処理部５２に供給する。表示処理部５２は、この情報に応じて、図７（Ｃ）に示すように、最新の心拍の終了に対応するラインを破線とし、次に新しい心拍に対応するラインが実線として、この心拍が選択された状態となる。
さらにトラックボールが左に回されると、拍選択部４８は、その回転量に応じて、３番目に新しい心拍が選択されたとして、その旨の情報を表示処理部５２に供給し、２番目に新しい拍に対応するラインが破線になり、３番目に新しい心拍に対応するラインが実線となって、選択された状態となる。
また、トラックボールを右に回せば、同様に、順次、新しい心拍に対応するラインが選択された状態となる。
【００６３】
さらに、トラックボールの移動によって、選択された心拍が変更された場合には、表示処理部５２は、選択された心拍の情報に応じて、この心拍の前記所定位置に対応するＢモード画像を記憶部３６から読み出して、ディスプレイ１８に表示する。すなわち、Ｂモード画像６４を、新しく選択された心拍の所定位置の画像に変更する。
図示例においては、所定位置として、心拍の開始位置が設定されている。従って、例えば、図７（Ｃ）に示すように、トラックボールによって２番目に新しい心拍が選択された場合には、表示処理部５２は、この心拍の開始位置（開始時点）のＢモード画像を記憶部３６から読み出し、ディスプレイ１８にＢモード画像６４として表示する。
【００６４】
心拍の選択が可能な状態でセットボタンが押されると、心拍の選択が終了したとして、選択された心拍が確定し、選択した心拍の微調整が行なえる状態となる。
また、ディスプレイ１８に表示しているＭモード画像６５中の心拍が選択／確定すると、拍選択部４８は、記憶部３６が記憶している全てのＭモード画像（すなわち、Ｂモード画像６４のアジマス方向の全域のＭモード画像）において、必要に応じて、選択された心拍の情報を変更する。すなわち、図示例においては、最新の心拍以外が選択されて確定された場合には、選択された心拍の情報を新たに選択／確定された心拍に変更する。
また、選択された心拍が確定すると、選択された心拍の情報は、トラッキング部４２に供給される。
【００６５】
一例として、最新の心拍が選択（すなわち、自動選択された心拍からの変更は無し）されたとして、セットボタンが押されると、図８（Ａ）に示すように、まず、選択された心拍の終了に対応する線が細線となり、選択された心拍の開始に対応する線の位置（時間）が、トラックボールによって、矢印ｔで示すように、左右方向（時間方向）に移動可能となり、心拍の開始位置の微調整が可能な状態となる。
必要に応じて、トラックボールによって心拍の開始位置が調整された後、再度、セットボタンが押されると、今度は、図８（Ｂ）に示すように、選択された心拍の終了に対応する線が通常の実線に、同開始に対応する線が細線になって、トラックボールによって、矢印ｔで示すように、選択された心拍の終了に対応する線の位置が、左右方向に移動可能となり、心拍の終了位置の微調整が可能な状態となる。
この心拍の微調整の結果は、微調整を行ったＭモード画像６５のみに反映してもよいが、記憶部３６に記憶されている全てのＭモード画像にも反映されるのが好ましい。
なお、心拍の開始位置が調整された場合には、表示処理部５２は、記憶部３６から、調整された心拍開始位置のＢモード画像を読み出し、ディスプレイ１８に表示するＢモード画像６４を、この画像に変更する。
【００６６】
心拍の微調整の結果は、表示処理部５２にも供給される。表示処理部５２は、心拍の微調整に応じて、選択された心拍中における所定位置が変動した場合には、変動した所定位置に応じたＢモード画像を記憶部３６から読み出し、ディスプレイ１８に表示する。
図示例においては、心拍の開始位置が所定位置として設定されている。従って、必要に応じて心拍の開始位置が微調整されたら、表示処理部５２は、調整された心拍開始位置に対応するＢモード画像を記憶部３６から読み出して、ディスプレイ１８に表示する。
また、心拍の微調整の結果が、記憶部３６に記憶されている全てのＭモード画像にも反映される場合には、この心拍の微調整に応じた所定位置の変動は、記憶部３６に記憶されている全てのＭモード画像にも反映されるのが好ましい。
【００６７】
選択された心拍の終了位置が調整可能な状態で、セットボタンが押されると、前述の図６に示すＢモード画像６４の選択ライン６２が移動可能な状態、すなわち、Ｂモード画像６４中における、Ｍモード画像６５の表示ラインの選択が可能な状態に戻る。
すなわち、図示例の超音波診断装置１０では、「表示ラインの選択」→「心拍の選択」→「心拍の微調整」の各処理を繰り返し行うことができるようになっている。言い換えれば、「表示ラインの選択」→「心拍の選択」→「心拍の微調整」は、ループ状に処理を行なうことが可能になっている。
これにより、より好適に、記憶する全てのＭモード画像から、後述する血管壁の弾性計測のための解析に最適な心拍を選択することが可能になる。
【００６８】
一方、選択された心拍の終了に対応する位置が調整可能な状態で、セットボタンではなく、タッチパネルの「ＡＷ Ｄｅｔ」ボタンが押されると、図９に示すように、Ｂモード画像６４の選択ライン６２、および、Ｍモード画像６５中における心拍を示すラインが、全て、破線となって操作ができない状態となり、血管壁の検出モードとなる。
【００６９】
血管壁の検出モードになると、まず、図１０（Ａ）に示すように、Ｂモード画像６４中に、血管前壁の外膜中膜境界に対応するライン６８が表示される。
このライン６８は、トラックボールによって上下方向（深度方向）に平行移動可能になっている。図１０（Ｂ）に示すように、トラックボールによって移動して、ライン６８を血管前壁の外膜中膜境界の位置に移動したら、セットボタンが押される。
【００７０】
セットボタンが押されると、図１０（Ｃ）に示すように、Ｂモード画像６４において、血管前壁の外膜中膜境界に対応するライン６８が破線になって確定し、血管前壁の内膜内腔境界に対応するライン７０が表示される。
このライン７０も、同様に、トラックボールによって上下方向に移動可能であり、ライン７０を血管前壁の中膜内腔境界の位置に移動したら、セットボタンが押される。
【００７１】
ライン７０が移動可能な状態でセットボタンが押されると、図１０（Ｄ）に示すように、Ｂモード画像６４において、血管前壁の中膜内腔境界に対応するライン７０が破線になって確定し、血管後壁の内膜内腔境界に対応するライン７２が表示される。同様に、トラックボールによってライン７２を血管後壁の中膜内腔境界の位置に移動したら、セットボタンが押される。
さらに、ライン７２が移動可能な状態でセットボタンが押されると、図１０（Ｅ）に示すように、Ｂモード画像６４において、血管後壁の内膜内腔境界に対応するライン７２が破線になって確定し、血管後壁の外膜中膜境界に対応するライン７４が表示される。同様に、トラックボールによってライン７４を血管後壁の外膜中膜境界の位置に移動したら、セットボタンが押される。
【００７２】
血管壁の各境界の情報は、境界検出部４０に供給される。
ライン７４が移動可能な状態でセットボタンが押されると、全ての境界に対応するラインの設定が終了し、境界検出部４０は、設定された内膜内腔境界のライン７２および外膜中膜境界のライン７４を用いて、後壁の内膜内腔境界および外膜中膜境界の自動検出を行なう。両境界の自動検出の結果は、表示処理部５２およびトラッキング部４２に送られ、図１０（Ｆ）に示すように、検出結果が表示される。
なお、これらの境界の自動検出の方法には、特に限定はなく、各種の方法が利用可能である。一例として、Ｂモード画像を解析して、ライン７２およびライン７４の位置において連続する高輝度な部分をトレースして、内膜内腔境界および外膜中膜境界を検出する方法が例示される。
【００７３】
境界検出部４０による、血管後壁の内膜内腔境界および外膜中膜境界の自動検出が終了すると、図１０（Ｆ）に示すように、Ｂモード画像６４にカーソル７８が表示される（血管後壁の自動検出が終了するまでは、このカーソル７８は表示されていない）。
【００７４】
このカーソル７８は、トラックボールによって移動可能になっている。カーソル７８を、自動検出された内膜内腔境界もしくは外膜中膜境界を示すラインに移動して、セットボタンが押されると、カーソル７８に近い側のラインが実線となる。実線となったラインは、修正が可能な状態となる。
一例として、図１０（Ｇ）に示すように、外膜内膜境界を示すライン７４が選択されて実線になったとする。ライン７４に沿って、トラックボールによってカーソル７８を移動して、再度、セットボタンが押されると、境界検出部４０が、カーソルによってなぞられた領域のライン７４を再検出して書き換えが行なわれ、また、その結果がトラッキング部４２に送られる。
【００７５】
後壁の内膜内腔境界および外膜中膜境界の自動検出が終了し、さらに、必要に応じて血管後壁の修正が行なわれると、図１１（Ａ）に示すように、全てのラインが破線になっている状態になり、さらに、図１１（Ｂ）に示すように、タッチパネル１６ａの「Elasticity Ana」ボタンが選択可能な状態となる。
「Elasticity Ana」ボタンが選択可能な状態となったら、「Ｐｓ」ボタンを用いて被検者の心臓収縮期の血圧を、「Ｐｄ」ボタンを用いて被検者の心臓拡張末期の血圧を、それぞれ入力し、さらに、「Quality Factor Threshold」ボタンを用いて、信頼性閾値を入力する。これらの数値の入力は、公知の方法で行なえばよい。
【００７６】
なお、被検者の血圧および信頼性閾値の入力は、血管壁境界の検出が終わった後に行なうのに限定はされず、後述する解析開始の前（後述する「Elasticity Ana」ボタンの押下の前）であれば、どのタイミングで行なってもよい。
また、超音波診断装置１０においては、診断を行なう前に、被検者情報の取得や入力を行なうのが通常であるが、この被検者情報に血圧の情報が有る場合には、これを利用してもよい。
【００７７】
被検者の血圧および信頼性閾値が入力され、「Elasticity Ana」ボタンが押されると、Ｂモード画像の解析が開始され、血管壁の弾性率の計算が行なわれる。
「Elasticity Ana」ボタンが押されると、まず、トラッキング部４２が、Ｍモード画像６５において、選択された心拍における血管前壁（外膜中膜境界および中膜内腔境界）ならびに血管後壁（中膜内腔境界および外膜中膜境界）の動きを追跡する。すなわち、血管前壁および後壁のトラッキングを行なう。
【００７８】
Ｍモード画像６５における血管壁のトラッキングは、先にＢモード画像６４において検出（設定）した、血管前壁の外膜中膜境界、血管前壁の中膜内腔境界、血管後壁の中膜内腔境界、および、血管後壁の外膜中膜境界を、位置的な出発点（焦点深度方向の出発点）として行なう。
他方、Ｍモード画像６５における血管壁のトラッキングの時間的な出発点は、ディスプレイ１８に表示した、選択された心拍の前記所定位置のＢモード画像６４に対応する位置（時相）である。すなわち、Ｍモード画像６５中におけるＢモード画像６４が撮られた時点が、トラッキングの出発点となる。従って、図示例においては、選択された心拍の開始位置が、トラッキングの出発点となる。
【００７９】
前述のように、血管壁の境界を指定したＢモード画像６４は、心拍の前２０％から後２０％の期間内で設定された、選択された心拍の所定位置に対応するＢモード画像である。
従って、図示例においては、トラッキングの出発点となるＢモード画像６４が、選択された心拍あるいは選択された心拍に近いので、トラッキングのエラーや、余分なノイズ等の影響を受けることなく、Ｍモード画像６５において、選択された心拍における血管壁のトラッキングを行なうことができる。
【００８０】
ここで、超音波診断装置１０においては、好ましい態様として、検出（設定）した血管壁の境界のみならず、血管後壁の中において、深さ方向に１以上の測定点が設定されていてもよい。このように、血管後壁の中に、１以上の測定点が設定されている場合には、各測定点毎に、血管壁のトラッキングを行なう。
なお、血管壁中の測定点は、予め設定されていてもよく、あるいは、特定のアルゴリズムに基づいて自動設定してもよく、あるいは、超音波診断装置１０の操作者が画像を見ながら設定してもよく、これらを併用してもよい。
【００８１】
Ｍモード画像６５における血管壁のトラッキングの方法には、特に限定はなく、トラッキングの出発位置からの画像（輝度）の連続性を利用する方法、パターンマッチング法、ゼロクロス法、組織ドプラ法、位相差トラッキング等が例示され、いずれの方法を用いてもよい。
【００８２】
トラッキング部４２によるＭモード画像中での血管壁のトラッキング結果は、弾性率算出部５０および表示処理部５２に供給される。
弾性率算出部５０は、血管壁のトラッキング結果から、まず、血管壁（内膜中膜）の厚みの変化波形、および、血管径（内径）の変化波形を形成する。なお、前述のように、血管壁内に１点以上の測定点が設定されている場合には、血管壁の変化波形は、各測定点の間毎に形成される。
血管壁の厚みの変化波形、および、血管径の変化波形は、表示処理部５２に送られる。
【００８３】
また、弾性率算出部５０は、下記式（１）を用いて、血管の径方向の歪みを算出する。
εi＝Δｈi／ｈdi ・・・（１）
なお、上記式（１）において、εは、各測定点の間における血管の径方向の歪みを、Δｈiは、１つの心拍内で血管壁が最も薄くなる心臓収縮期における各測定点の間の血管壁の厚み変化の最大値を、ｈdiは、血管壁が最も厚くなる心臓拡張末期における各測定点の間の厚みを、それぞれ示す。
【００８４】
さらに、弾性率算出部５０は、先に入力された血圧の最高値および最低値を用いて、下記式（２）によって、血管壁の円周方向の弾性率Ｅθiを算出する。
Ｅθi＝１／２＊［１＋（ｒd／ｈd）］＊［Δｐ／（Δｈi／ｈdi）］） ・・・（２）
もしくは、下記式（３）によって、血管壁の径方向の弾性率Ｅriを算出してもよい。
Ｅri＝Δｐ／（Δｈi／ｈdi） ・・・（３）
なお、上記式（２）および式（３）において、Δｈiおよびｈdiは、先と同様であり、Δｐは、心臓収縮期と心臓拡張末期とにおける血圧差を、ｒdは、心臓拡張末期での血管内腔の半径を、ｈdは、心臓拡張末期における血管壁の厚みを、それぞれ示す。
【００８５】
弾性率を計算したら、弾性率算出部５０は、弾性率の信頼性を計算する。
弾性率の信頼性の計算方法には、特に限定はなく、公知の方法が、各種、利用可能である。一例として、１０００人等の多数の人物の心拍による血管径変化の波形を作成して、これらの多数の波形から、血管径変化のモデル波形を作成し、このモデル波形からのズレの量を用いて、算出した弾性率の信頼性を計算する方法が例示される。
【００８６】
ここで、前述のように、ディスプレイ１８に表示されているＭモード画像６５で心拍が選択／確定すると、記憶部３６が記憶している全てのＭモード画像において、同じ心拍が選択された状態となる。
これに応じて、上述の血管壁のトラッキング、血管壁の厚みおよび血管径の変化波形の作成、血管壁の歪みの計算、血管壁の弾性率および弾性率の信頼性の計算などの処理は、ディスプレイ１８に表示されているＭモード画像６５のみならず、記憶部３６に記憶されている全てのＭモード画像において、選択された心拍に対して行なわれる。すなわち、ディスプレイ１８に表示されているＢモード画像６４のアジマス方向の全領域に対して、対応するＭモード画像を用いて、選択された心拍での上記血管壁の弾性率算出などの処理が行なわれる。
これらの結果は、記憶部３６に記憶されているＭモード画像に情報として付加される。
【００８７】
アジマス方向の全域における演算が終了したら、弾性率算出部５０は、血管壁の弾性率の平均値（Ｅθave）、血管壁の歪みの平均値（Ｓｔｒave）、および、弾性率の信頼性の平均値（ＱＦave）を、演算する。
【００８８】
演算が終了すると、ディスプレイ１８に結果が表示される。
その一例を、図１２に示す。図示例においては、元々、表示されていたＢモード画像６４の図中右側に、このＢモード画像６４に示される血管後壁の弾性率が、Ｂモード画像６４ｅで表示される。さらに、この血管後壁の弾性率を表示するＢモード画像６４ｅの図中右側に、算出した血管壁の弾性率の信頼性が、同じくＢモード画像６４ｑで表示される。
また、Ｂモード画像６４の図中左側に、血管壁の弾性率の平均値（Ｅθave）、血管壁の歪みの平均値（Ｓｔｒave）、および、弾性率の信頼性の平均値（ＱＦave）が、それぞれ、表示される。
【００８９】
血管壁の弾性率は、Ｂモード画像６４ｅにおいて、Ｂモード画像６４で自動検出（あるいは必要に応じて修正）された血管後壁に重ねて、帯状に表示される。また、Ｂモード画像６４ｅの右上側には、弾性率の指標が表示される。図示例においては、画像が高密度であるほど、弾性率が高い。
すなわち、Ｂモード画像６４ｅにおいては、血管後壁に重なる帯の密度が、血管のその位置における血管壁の弾性率を示している。
【００９０】
弾性率の信頼性は、Ｂモード画像６４ｑにおいて、同様に、Ｂモード画像６４中で自動検出された血管後壁に重ねて、帯状に表示される。また、Ｂモード画像６４ｑの右上側には、弾性率の信頼性の指標が表示される。図示例においては、画像が高密度であるほど、弾性率の信頼性が高い。
すなわち、Ｂモード画像６４ｑにおいては、血管後壁に重なる帯の密度が、血管のその位置における血管壁弾性率の信頼性を示している。
【００９１】
なお、この弾性率や弾性率の信頼度の高低は、画像の密度に変えて、あるいは加えて、画像の色によって表現してもよい。
【００９２】
ここで、図１２に示す結果の表示では、先に入力した閾値よりも信頼性が低い、アジマス方向の位置では、結果が自動的に省かれる。
また、結果が省かれた位置に関しては、Ｂモード画像６４ｅにおける弾性率の結果表示の右隅部や、Ｂモード画像６４ｑにおける信頼性の結果表示の右隅部に示されるように、帯の表示が薄くなる。
【００９３】
また、下側のＭモード画像６５においては、選択された心拍に、Ｍモード画像中での血管前壁のトラッキング結果８０および血管後壁のトラッキング結果８２、血管径の変化波形８４、ならびに、血管壁の厚みの変化波形８６が表示される。
なお、前述のように、血管壁の中に深度方向に１以上の測定点が設定されている場合には、血管厚みの変化波形は、各測定点の間毎に、出力するようにしてもよい。
【００９４】
ここで、血管壁の弾性率の計測結果等がディスプレイ１８に表示されると、Ｂモード画像６４における選択ライン６２が実線になって、トラックボールによってアジマス方向に移動可能になる。
Ｂモード画像６４で選択ライン６２が移動されると、表示処理部５２は、記憶部３６から、選択ライン６２の位置対応するＭモード画像を読み出して、ディスプレイ１８に表示する。すなわち、トラックボールによって選択ライン６２が移動されると、Ｍモード画像６５が、この選択ライン６２の位置のＭモード画像に変わり、Ｍモード画像中での血管前壁のトラッキング結果８０および血管後壁のトラッキング結果８２、血管径の変化波形８４、ならびに、血管壁の厚みの変化波形８６が、Ｂモード画像６４の選択ライン６２の位置のデータに変更される。
従って、Ｂモード画像中のアジマス方向の全域において、Ｍモード画像６５および解析結果を表示する表示ラインを、選択することができる。
【００９５】
また、このＢモード画像６４での選択ライン６２の移動に同期して、Ｂモード画像６４ｅの選択ライン６２ｅ、および、Ｂモード画像６４ｑの選択ライン６２ｑも、移動する。
【００９６】
さらに、セットボタンを押した後、Ｂモード画像６４ｅおよびＢモード画像６４ｑにおいて、選択ライン６２ｅおよび選択ライン６２ｑをトラックボールによって移動してアジマス方向の任意の領域を選択し、その後、再度、セットボタンを押すと、選択された領域は前述の信頼性が閾値よりも低かった領域と同様の扱いとなり、データが取り除かれる。
すなわち、検査者が結果を見て、波形等が奇怪しいと感じられる場所が有る場合に、そのデータを取り除くことができ、より正確な解析が可能となる。
【００９７】
なお、このデータの除去は、Ｄｅｌｅｔボタン等の押下によって、１つ前の状態にも土切要にしてもよい。
【００９８】
以上、本発明の超音波診断装置について詳細に説明したが、本発明は、上述の例に限定はされず、本発明の要旨を逸脱しない範囲において、各種の変更や改良を行なってもよいのは、もちろんである。
【産業上の利用可能性】
【００９９】
本発明の超音波診断装置は、心筋梗塞、狭心症、脳疾患等の原因となる動脈硬化の診断を行なう医療現場等において、好適に利用可能である。
【符号の説明】
【０１００】
１０ 超音波診断装置
１２ 診断装置本体
１４ （超音波）プローブ
１６ 操作パネル
１８ ディスプレイ
２０ ケーブル
２４ キャスタ
２８ 送信回路
３０ 受信回路
３２ 送受信制御部
３４ 画像形成部
３６ 記憶部
４０ 境界検出部
４２ トラッキング部
４６ 拍検出部
４８ 拍選択部
５０ 弾性率算出部
５２ 表示処理部
５６ Ｂモード画像形成部
５８ Ｍモード画像形成部
６０ ＲＯＩ
６２ 選択ライン
６４，６４ｅ，６４ｑ，９０ Ｂモード画像
６５，９２ Ｍモード画像
６８，７０，７２，７４ ライン
８０，８２ トラッキング結果
８４ 血管径変化波形
８６ 血管壁厚み変化波形
【技術分野】
【０００１】
本発明は、血管壁の弾性率計測に好適な超音波診断装置に関し、特に、Ｂモード画像とＭモード画像との同時表示において、解析する心拍に対応する適正なＢモード画像を表示することができる超音波診断装置に関する。
【背景技術】
【０００２】
医療分野において、超音波画像を利用した超音波診断装置が実用化されている。
一般に、この種の超音波診断装置は、超音波プローブ（超音波探触子 以下、プローブとする）と、診断装置本体とを有しており、プローブから被検体に向けて超音波を送信し、被検体からの超音波エコーをプローブで受信して、その受信信号を診断装置本体で電気的に処理することにより超音波画像が形成される。
【０００３】
また、超音波を血管や心臓壁等に送信し、その超音波エコーを受信して、受信信号を解析することで、心臓壁や血管壁等の変位量を求め、この変位量から、血管壁や心臓壁（心筋）等の弾性率を計測することも行なわれている。
【０００４】
例えば、特許文献１には、心拍（心臓拍動）に同期して動く対象物に超音波を送受信して超音波エコーの受信信号を得、この受信信号の振幅および位相を用いて対象物の瞬間的な瞬時的な位置を決定して、心拍に基づく血管壁の大きな振幅変位運動をトラッキングすることにより、血管の弾性率を求めることが記載されている。
具体的には、血管壁の順次の位置に基づき、血管壁の微小振動の運動速度波形を求め、血管壁内部の深さ方向に所定間隔でとられた局所ごとのトラッキング軌跡を求め、局所毎の厚みの時間変化を算出することにより血管の弾性率を求めている。
【０００５】
また、特許文献２にも、同様に、心拍に同期して動く対象物に超音波を送受信して得られた超音波エコーの受信信号から、血管等の変位量を求め、この変位量から弾性率を求める超音波診断装置が記載されている。
この超音波診断装置では、血管等の対象物から得た受信信号を用いてＢモード画像およびＭモード画像を形成し、Ｍモード画像の受信信号から手ブレや体動ブレを検出して、このブレが検出されたＭモード画像の受信信号を用いてプローブと被検体との位置変化量を検出して、この検出結果から受信信号の確度を判断し、確度が高いと判断されたＭモード画像の受信信号を用いて、対象物の変位量を求め、この変位量から血管壁等の弾性率を計測している。
【０００６】
ところで、特許文献２に示されるように、Ｂモード画像およびＭモード画像の同じ表示を行なう超音波診断装置においては、通常、画像を静止画にするフリーズボタンが押圧されると、フリーズボタンが押された時点でＭモード画像を静止し、かつ、フリーズボタンが押された時点のＢモード画像を表示する。
そのため、表示されるＢモード画像は、必ずしも、行なう診断などに適した画像であるとは限らない。
【０００７】
これに対して、特許文献３には、Ｂモード画像とＭモード画像とを同時に表示する際に、心電図等の生体組織の活動情報に基づく波形図を表示すると共に、この波形図の特定時相をＭモード画像中に表示し、かつ、この特定時相に対応するＢモードのサムネイル画像を表示する超音波診断装置が記載されている。また、この特許文献３に記載される超音波診断装置では、Ｍモード画像中で時相を任意に選択することにより、この時相に対応するＢモードのサムネイル画像を表示する。
特許文献３に記載される超音波診断装置においては、このような構成を有することにより、弾性率や歪みなどの生体組織の性状に関する情報等を効率よく取得して、超音波診断を内の診断効率を向上することを可能にしている。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開平１０−５２２６号公報
【特許文献２】特開２０１０−２３３９５６号公報
【特許文献３】特開２００４−８３５０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００９】
心拍に同期して動く対象物の変位量から、血管壁の弾性率等を正確に計測するためには、Ｍモード画像から解析に適した心拍を選択すると共に、Ｂモード画像も、選択した心拍すなわち解析する心拍に対応し、かつ、解析にも好適なＢモード画像である必要が有る。
しかしながら、従来の超音波診断装置では、必ずしも、Ｍモード画像から任意の心拍を選択できるとは限らず、さらに、Ｍモード画像で任意の心拍を選択できても、選択された心拍に対応し、かつ、解析にも好適なＢモード画像が表示されるとは限らない。
【００１０】
本発明の目的は、前記従来技術の問題点を解決することにあり、Ｂモード画像とＭモード画像とを表示して、Ｍモード画像の中から血管弾性率の計測を行なうための解析に適した心拍を選択ことができ、かつ、Ｍモード画像中で選択された心拍に対応し、解析にも好適なＢモード画像を表示することができる超音波診断装置を提供することにある。
【課題を解決するための手段】
【００１１】
前記目的を達成するために、本発明の超音波診断装置は、超音波を送信し、被検体によって反射された超音波エコーを受信して、受信した超音波エコーに応じた受信信号を出力する超音波トランスデューサを有する超音波プローブと、前記超音波トランスデューサが出力した受信信号から、Ｂモード画像およびＭモード画像を形成する、画像形成手段と、表示手段と、前記画像形成手段が形成したＢモード画像およびＭモード画像の少なくとも一方を、前記表示手段に表示させる表示処理手段と、前記Ｂモード画像およびＭモード画像が表示された状態において、表示画像を静止画にするフリーズ手段と、前記Ｂモード画像およびＭモード画像が静止画にされた状態において、Ｍモード画像中で心拍を選択する心拍選択手段とを有し、かつ、前記表示処理手段は、前記Ｍモード画像中で心拍が選択されたら、心拍の長さに対して心拍の前２０％から心拍の後２０％までの期間内で予め設定された所定位置に対応して、前記選択された心拍における前記所定位置のＢモード画像を前記表示手段に表示させることを特徴とする超音波診断装置を提供する。
【００１２】
このような本発明の超音波診断装置において、前記Ｍモード画像中の心拍を検出する拍検出手段を有し、前記表示処理手段は、前記拍検出手段が心拍を検出したら、Ｍモード画像中に、前記拍検出手段が検出した全ての心拍を表示するのが好ましい。
また、前記拍検出手段は、前記Ｍモード画像を解析することにより、前記Ｍモード画像中の心拍を検出するのが好ましい。
【００１３】
また、前記所定位置を選択する選択手段を有するのが好ましい。
また、前記心拍選択手段は、前記拍検出手段が検出した心拍の内、開始から終了までの全てがＭモード画像中で取得されている、最も新しい心拍が選択されたとするのが好ましい。
また、前記所定位置は、心臓の拡張期よりも後で、かつ、心臓の収縮期よりも前の期間以外に設定されるのが好ましい。
【００１４】
前記心拍選択手段によって選択された心拍の位置を調整する位置調整手段を有するのが好ましい。
また、前記位置調整手段によって心拍の位置調整が行なわれたら、前記表示処理手段は、この心拍の位置調整および前記所定位置に応じて、前記表示手段に表示するＢモード画像を対応する位置の画像に変更するのが好ましい。
また、前記Ｂモード画像中において血管壁の境界位置を設定する、血管壁境界設定手段を有するのが好ましい。
【００１５】
また、表示手段に表示されたＢモード画像中で、関心領域を設定するＲＯＩ設定手段を有するのが好ましい。
さらに、前記関心領域の設定指示に応じて、前記超音波トランスデューサによる超音波のフレームレートを、設定指示前よりも高くするのが好ましい。
【発明の効果】
【００１６】
上記構成を有する本発明の超音波診断装置は、Ｂモード画像とＭモード画像とを同時に表示して、血管弾性の計測等を行なうに際し、Ｍモード画像中から目的とする計測を行なうための解析に適した心拍を選択できると共に、心拍に対する所定期間内で予め設定された所定位置に応じて、選択された心拍における所定位置のＢモード画像を表示する。
そのため、本発明の超音波診断装置では、血管弾性率の計測等に適した心拍を選択した上で、この選択した心拍に対応する、解析に適したＢモード画像を表示ずることができ、正確な血管弾性率の計測等を、安定して行なうことができる。
【図面の簡単な説明】
【００１７】
【図１】本発明の超音波診断装置の一例を概念的に示す図である。
【図２】図１に示す超音波診断装置の構成を概念的に示すブロック図である。
【図３】図１に示す超音波診断装置における血管壁の弾性計測の一例を説明するためのフローチャートである。
【図４】血管壁の弾性計測のための超音波診断を説明するための概念図である。
【図５】（Ａ）および（Ｂ）は、図１に示す超音波診断装置での画像表示の一例を示す概念図である。
【図６】（Ａ）および（Ｂ）は、図１に示す超音波診断装置での画像表示の一例を示す概念図、（Ｃ）は、本発明の超音波診断装置の作用を説明するための概念図である。
【図７】（Ａ）〜（Ｃ）は、図１に示す超音波診断装置での画像表示の一例を示す概念図である。
【図８】（Ａ）および（Ｂ）は、図１に示す超音波診断装置での画像表示の一例を示す概念図である。
【図９】図１に示す超音波診断装置での画像表示の一例を示す概念図である。
【図１０】（Ａ）〜（Ｇ）は、図１に示す超音波診断装置での画像表示の一例を示す概念図である。
【図１１】（Ａ）および（Ｂ）は、図１に示す超音波診断装置での画像表示の一例を示す概念図である。
【図１２】図１に示す超音波診断装置での画像表示の一例を示す概念図である。
【発明を実施するための形態】
【００１８】
以下、本発明の超音波診断装置について、添付の図面に示される好適実施例を基に、詳細に説明する。
【００１９】
図１に、本発明の超音波診断装置の一例の外観を概念的に示す。
図１に示すように、超音波診断装置１０は、基本的に、診断装置本体１２と、超音波プローブ１４と、操作パネル１６と、ディスプレイ１８とを有して構成される。また、超音波診断装置１０の下端部には、キャスタ２４が配置されており、人力で容易に装置を移動することが可能になっている。
【００２０】
超音波プローブ１４（以下、プローブ１４とする）は、超音波の送受信を行なって、受信した超音波エコーに応じた受信信号を診断装置本体１０に供給するものである。
このプローブ１４は、被検体に超音波を送信し、被検体に反射された超音波エコーを受信して、受信した超音波エコーに応じた電気信号（受信信号）を出力する、いわゆる超音波トランスデューサ（超音波圧電素子）を１次元的もしくは二次元的に配列してなる、各種の超音波診断装置に利用される公知の超音波プローブである。
【００２１】
本発明において、プローブ１４の種類には、特に限定はなく、コンベックス型、リニア型、セクタ型等の各種の形式が利用可能である。また、体外式プローブでもよいし、ラジアルスキャン方式等の超音波内視鏡用プローブでもよい。さらに、プローブ１４は、ハーモニックイメージングに対応する、送信した超音波の二次以上の高調波を受信するための超音波振動子を有するものであってもよい。
また、図示例においては、プローブ１４と診断装置本体１２とは、ケーブル２０によって接続されている。しかしながら、本発明は、これに限定はされず、プローブ１４内に後述する送信回路２８、受信回路３０、送受信制御部３２等を配置して、無線通信によってプローブ１４と診断装置本体１２とを接続するものであってもよい。
【００２２】
ディスプレイ１８は、公知のディスプレイ（表示装置）である。
超音波診断装置１０において、ディスプレイ１８は、各種の超音波診断装置と同様、プローブ１４が出力した受信信号に応じた超音波画像、被検者の情報、ＧＵＩ（Graphical User Interface)による操作を行なうための選択手段や指示手段、関心領域（Region of Interest 以下、ＲＯＩとする）、後述する血管壁の弾性計測結果等を表示する。
【００２３】
操作パネル１６は、超音波診断装置１０の操作を行なうものである。
図示は省略するが、超音波診断装置１０において、操作パネル１６には、ＢモードやＭモードなどの各種のモードの選択手段、ディスプレイ１８に表示されたカーソルやライン等を移動するためのトラックボール（トラックパッド／タッチパッド）、選択や操作を決定（確定）するためのセットボタン、動画表示と静止画表示との切り換え等を行なうためのフリーズボタン、超音波画像の視野深度の変更手段、ゲイン調整手段、超音波画像を拡大するためのズームボタン等が配置される。
なお、超音波診断装置１０においては、モードとして、ＢモードやＭモードなどの通常の超音波診断装置が有するモードに加え、血管壁の弾性率を計測するためのモードであるＶＥモード（Vascular Elasticity モード）も設定されている。
また、同じく図示は省略するが、操作パネル１６には、ＧＵＩによる操作等を行なうための表示装置であるタッチパネル１６ａも配置される（図６（Ｂ）参照）。
【００２４】
診断装置本体１２は、超音波診断装置１０の全体の動作の制御を行なうと共に、プローブ１４が出力した受信信号に応じた超音波画像を形成してディスプレイ１８に表示させ、さらに、血管弾性率の計測のための各種の処理を行なうものである。
診断装置本体１２は、例えば、コンピュータを利用して構成される。
【００２５】
図２に、超音波診断装置１０の構成をブロック図で概念的に示す。
図２に示すように、診断装置本体１２は、送信回路２８、受信回路３０、送受信制御部３２、画像形成部３４、記憶部３６、境界検出部４０、トラッキング部４２、拍検出部４６、拍選択部４８、弾性率算出部５０、および、表示処理部５２を有する。
また、画像形成部３４は、Ｂモード画像形成部５６およびＭモード画像形成部５８を有する。
【００２６】
前述のプローブ１４は、送信回路２８および受信回路３０に接続される。また、送信回路２８および受信回路３０には、送受信制御部３２が接続される。さらに、受信回路３０は、画像形成部３４に接続される。
画像形成部３４は、表示処理部５２に接続される。また、画像形成部３４のＢモード画像形成部５６およびＭモード画像形成部５８は、記憶部３６に接続される。Ｂモード画像形成部５８は、さらに境界検出部４０にも接続される。
記憶部３６は、トラッキング部４２、拍検出部４６、および表示処理部５２に接続される。また、拍検出部４６および表示処理部５２は、共に、トラッキング部４２および表示処理部５２に接続される。トラッキング部４２は、さらに、弾性率算出部５０にも接続され、この弾性率算出部５０は、表示処理部５２に接続される。また、拍検出部４６は、さらに、拍選択部４８にも接続され、この拍選択部４８も、同じく表示処理部５２に接続される。
【００２７】
送受信制御部３２は、送信回路２８および受信回路３０を介してプローブ１４の超音波ビームの送信方向および超音波エコーの受信方向を、順次、設定する。
また、送受信制御部３２は、設定した送信方向に応じて送信遅延パターンを選択する送信制御機能と、設定した受信方向に応じて受信遅延パターンを選択する受信制御機能とを有している。
【００２８】
送信遅延パターンとは、プローブ１４の複数の超音波トランスデューサから送信される超音波によって所望の方向に超音波ビームを形成するために各超音波トランスデューサの駆動信号に与えられる遅延時間のパターンである。他方、受信遅延パターンとは、複数の超音波トランスデューサによって受信される超音波によって所望の方向からの超音波エコーを抽出するために受信信号に与えられる遅延時間のパターンである。
複数の送信遅延パターンおよび複数の受信遅延パターンが内部メモリ（図示せず）に格納されていて、状況に応じて、適宜、選択して使用する。
【００２９】
送信回路２８は、複数のチャネルを備えており、プローブ１４の複数の超音波トランスデューサにそれぞれ印加する複数の駆動信号を形成する。その際に、送受信制御部３２によって選択された送信遅延パターンに基づいて、複数の駆動信号にそれぞれの遅延時間を与えることができる。
なお、送信回路２８は、複数の超音波トランスデューサから送信される超音波が超音波ビームを形成するように、複数の駆動信号の遅延量を調節して複数の駆動信号をそれぞれプローブ１４の複数の超音波トランスデューサに供給するようにしても良く、複数の超音波トランスデューサから一度に送信される超音波が被検体の撮像領域全体に届くように構成した複数の駆動信号をプローブ１４に供給するようにしても良い。
【００３０】
受信回路３０は、送信回路２８と同様に複数のチャネルを備えており、複数の超音波トランスデューサを介して受信された複数のアナログ信号を増幅し、デジタルの受信信号に変換する。
さらに、送受信制御部３２によって選択された受信遅延パターンに基づいて、複数の受信信号にそれぞれの遅延時間を与え、それらの受信信号を加算することにより、受信フォーカス処理を行う。この受信フォーカス処理によって、超音波エコーの焦点が絞り込まれた音線信号（音線データ）が形成される。
【００３１】
形成された音線データは、画像形成部３４に供給される。
画像形成部３４は、供給された音線データに対して、Ｌｏｇ（対数）圧縮やゲイン調整等のプリプロセス処理を施して超音波画像の画像データを形成し、この画像データを通常のテレビジョン信号の走査方式に従う画像データに変換（ラスター変換）し、さらに階調処理等の必要な画像処理を施した上で表示処理部５２へ出力する。
また、画像形成部３４は、Ｂモード画像を形成するＢモード画像形成部５６と、Ｍモード画像を形成するＭモード画像形成部５８とを有する。Ｂモード画像およびＭモード画像の形成は、公知の方法によればよい。
【００３２】
表示処理部５２は、画像形成部３４から供給された超音波画像の画像データ、記憶部３６から読み出した超音波画像の画像データ、操作パネル１６で行なわれた操作（入力指示）、後述する心拍の検出結果や心拍の選択、血管壁弾性率の計測結果（解析結果）等に応じて、ディスプレイ１８に表示するための表示用データを形成し、ディスプレイ１８に表示させる部位である。
【００３３】
図示例の超音波診断装置１０において、診断装置本体１２の記憶部３６、境界検出部４０、トラッキング部４２、拍検出部４６、拍選択部４８、および、弾性率算出部５０は、主に、血管壁の弾性率を計測するＶＥモードの際に用いられる部位である。
以下、図３のフローチャートおよび図５〜図１２を参照して、ＶＥモードにおける超音波診断装置１０の作用を説明することにより、上記記憶部３６や境界検出部４０等の各部位、ならびに、本発明の超音波診断装置１０について、より詳細に説明する。
なお、以下の説明では、特に記載がなくても、ディスプレイ１８の表示に関しては、表示処理部５２がラインの形成等の必要な処理を行なう。
【００３４】
超音波診断装置１０による超音波診断が開始されると、送受信制御部３２による制御の下、送信回路２８がプローブ１４の超音波トランスデューサから超音波を送信させ、また、受信回路３０はプローブ１４が出力した受信信号を処理して音線信号を形成し、画像形成部３４に出力する。
一例として、Ｂモードが選択され、図４に概念的に示すように、被検者の頸動脈ｃを測定対象として、プローブ１４が首ｎに当てられたとして、画像形成部３４（Ｂモード画像形成部５６）によって形成されたＢモード画像が、表示処理部５２で処理されて、ディスプレイ１８に表示される。
【００３５】
目的とする頸動脈ｃが適正に観察できるようになり、操作パネル１６のモード選択手段（以下の説明では、「操作パネル１６」は省略する）によってＶＥモードが選択されると、表示処理部５２は、図５（Ａ）に概念的に示すように、Ｂモード画像中に、関心領域を示すＲＯＩ６０を表示させる。
【００３６】
この状態では、トラックボールでの操作によってＢモード画像中のＲＯＩ６０の位置を移動できる。また、セットボタンを押すと、ＲＯＩ６０の位置が固定され、トラックボールでの操作によってＲＯＩ６０のサイズを変更できる。
さらに、セットボタンを押すたびに、ＲＯＩ６０の位置変更およびＲＯＩ６０のサイズ調整が、交互に実施可能になる。
【００３７】
この状態からズームボタンが押されると（押下されると）、ＲＯＩ６０の位置やサイズの調整が終了してＲＯＩ６０の設定が指示されたとして、送受信制御部３２は、フレームレートをＲＯＩ６０の設定指示前よりも高くする（例えば、２００Ｈｚ以上あるいはＲＯＩ設定指示前の５倍以上）。また、このズームボタンの押下に応じて、Ｍモード画像形成部５８によって、ＲＯＩ６０のＭモード画像の形成が開始され、図５（Ｂ）に示されるように、ＲＯＩ６０の部分が拡大されたＢモード画像６４、および、ＲＯＩ６０（その選択ライン６２）のＭモード画像６５が、同時に表示される。
なお、Ｂモード画像６４とＭモード画像６５との同時表示（デュアルモード表示）は、公知の超音波診断装置における、いわゆるＢ／Ｍモード表示と同様に行なえば良い。
【００３８】
図５（Ｂ）において、上側がＢモード画像６４で、下側がＭモード画像６５である。
Ｂモード画像６４において、図中横方向はアジマス方向（超音波振動子の配列方向（二次元配列の場合は長手方向））であり、縦方向は深度方向（超音波の送受信方向）で、上方が深度が浅い側（プローブ１４側）である。
また、Ｂモード画像中には、Ｂモード画像中のアジマス方向におけるＭモード画像表示位置（Ｍモード画像の表示ライン）を選択するための、深度方向に延在する選択ライン６２が表示される。この選択ライン６２は、トラックボールによって、アジマス方向（図中左右方向）に移動可能になっている。
【００３９】
また、Ｍモード画像６５において、横方向は時間軸で、時間は左から右に流れており、間隙６５ａの左側が、現在のフレームとなる（すなわち、間隙６５ａの右側は過去のフレーム）。また、Ｂモード画像６４と同様、縦方向は深度方向で、上方が深度が浅い側である。
図５（Ｂ）において、ディスプレイ１８に表示されているＭモード画像６５は、予め位置が設定された選択ライン６２の位置のＭモード画像６５になっている。
【００４０】
ここで、Ｍモード画像形成部５８は、アジマス方向の所定位置（予め設定された所定位置や選択された位置）や、アジマス方向で選択された位置のみではなく、Ｂモード画像６４のアジマス方向の全域について、Ｍモード画像を形成する。
【００４１】
Ｂモード画像形成部５６が形成したＲＯＩ６０のＢモード画像（Ｂモード画像データ）、および、Ｍモード画像形成部５８が形成したＭモード画像（Ｍモード画像データ）は、共に、記憶部３６に記憶される。
なお、記憶部３６が記憶する画像の時間的な量には、特に限定はないが、一般的な心拍が２以上、入る長さであるのが好ましい。従って、記憶部３６は、最新の３秒以上のＢモード画像およびＭモード画像を記憶するのが好ましい。
【００４２】
前述のように、選択ライン６２は、トラックボールによってアジマス方向に移動することができる。
選択ライン６２の位置とＭモード画像とは、連動している。すなわち、トラックボールによって選択ライン６２を左右方向に移動すると、表示処理部５２は、選択ライン６２の位置のＭモード画像をディスプレイ１８に表示する。
【００４３】
操作者が、適正な画像が得られたと判断したら、フリーズボタンが押される。
フリーズボタンが押されると、表示処理部５２は、必要な画像データを記憶部３６から読み出し、図６（Ａ）に示すように、ディスプレイ１８に、フリーズボタンが押された時点が最も右となるように並べ直して、選択ライン６２の位置のＭモード画像６５の静止画を表示させ、また、後述する所定位置のＢモード画像６４の静止画を表示させる。同時に、選択ライン６２が破線になって、移動ができなくなる（非アクティブになる）。
さらに、図６（Ｂ）に示すように、操作パネル１６のタッチパネル１６ａに、後述する血管壁の境界の設定を指示するための「ＡＷ Ｄｅｔ」ボタン、血管壁弾性率の解析開始を指示するため「Elasticity Ana」ボタン、被検者の血圧を入力するための「Ｐｓ」ボタンおよび「Ｐｄ」ボタン、信頼性閾値を入力するための「Quality Factor Threshold」ボタンが表示される。なお、この時点では、「Elasticity Ana」ボタンは、選択できない状態になっている。
【００４４】
また、フリーズボタンが押されると、拍検出部４６が、記憶部３６に記憶されている全てのＭモード画像に対して、心拍の検出（心拍の自動検出）を行なう。心拍の検出結果は、記憶部３６に送られ、対応するＭモード画像に、情報として付加される。
さらに、心拍の検出結果は、表示処理部５２および拍選択部４８にも送られる。
【００４５】
心拍の検出の検出方法には、特に限定はないが、一例として、Ｍモード画像を解析（例えば時間方向にトラッキング）して、横方向の延在する白線（輝線）の深度方向の移動速度（速度の上昇開始時点）や、同白線の深度方向の動きの脈動等を用いて、検出すればよい。あるいは、心拍の検出に心電計（心電図）を利用してもよい。
また、心拍を自動検出するのではなく、操作者がＭモード画像を見て心拍の位置（心拍の開始と終了の位置）を入力するようにしてもよい。
【００４６】
表示処理部５２は、図６（Ａ）に示すように、心拍（心拍の開始／終了位置）の検出結果を、Ｍモード画像６５中に三角マークおよび直線で表示させる。
なお、検出に失敗した心拍が有る場合には、周囲の心拍の間隔等に応じて、適当な位置に心拍を表示してもよい。
また、拍選択部４８は、検出部４６による心拍の検出結果に応じて、心拍の開始から終了までが含まれる完全な心拍（撮り切れた心拍）の内、最も新しい心拍が選択されたとして、その情報を表示処理部５２に送る。すなわち、拍選択部４８は、２本の直線で挟まれた撮り切れている完全な心拍で、かつ、フリーズボタンが押される直前の心拍が選択されたとして、その情報を表示処理部５２に送る。
さらに、拍検出部４８は、記憶部３６に記憶される全てのＭモード画像に対して、同じ心拍が選択されたとする情報を付す。
【００４７】
表示処理部５２は、拍選択部４８による心拍の選択結果に応じて、一例として、Ｍモード画像６５において、選択された心拍を実線で、それ以外の心拍を破線で示す。図６（Ａ）に示す例では、３つの完全な心拍が検出されているので、前述のように、最新の完全な心拍である一番右の心拍が選択されて、この心拍の開始位置および終了位置（開始時点および終了時点）が実線で示され、それ以外の心拍の開始／終了位置は、破線で示される。
なお、このような選択および非選択の区別は、線種に変えて、あるいは加えて、線の色を利用してもよい。
【００４８】
なお、超音波診断装置１０においては、この後、操作者が心拍を選択（選択する心拍を変更）することが可能であり、この工程で確定された心拍が、最終的に選択された心拍として、後の血管弾性計測のための解析に供される。従って、この拍検出部４８による心拍の自動選択は、いわば仮の心拍の選択である。
また、拍検出部４８が自動選択する心拍は、最新の完全な心拍に限定はされず、その前の心拍であってもよく、検出された最も古い心拍であってもよい。さらに、超音波診断装置１０は、拍検出部４８が自動選択する心拍を操作者が選択するための、選択手段を有しても選択できるようにしてもよい。拍検出部４８が自動選択する心拍の選択手段は、ＧＵＩ等を用いた公知の方法で構成すればよい。
【００４９】
また、フリーズが押されると、Ｂモード画像６４は、Ｍモード画像６５で選択された心拍における、予め設定された所定位置（所定位置の時点（所定の時相））の画像となる。
【００５０】
具体的には、超音波診断装置１０においては、心拍の長さ（一拍の時間）に対して、心拍の前２０％の時点から心拍の後２０％の時点までの期間内で、Ｂモード画像６４を表示する所定位置（いわば、心拍（Ｍモード画像）の時間軸上における所定位置）が、予め設定されている。すなわち、図６（Ｃ）に血管の後壁境界のＭモード画像を模して概念的に示すように、心拍の長さをｔとすると、心拍の前０．２ｔの時点から心拍の後０．２ｔの時点までの期間Ｔ（長さｔの心拍を中心とする１．４ｔ）の中から、心拍開始位置や心拍中央のような所定位置が、予め設定されている。
なお、図６（Ｃ）では、図６（Ａ）等と同様、時間の進行方向は左から右であり、上方が深度が浅く、下方が深い。
【００５１】
これに応じて、超音波診断装置１０では、ＲＯＩが設定されてフリーズボタンが押され、さらに、心拍の検出および選択が行なわれたら、選択された心拍における、心拍の開始位置や心拍の中央位置などの、予め設定された所定位置のＢモード画像（この所定位置の時点で撮られたＢモード画像）を表示する。
【００５２】
図示例においては、一例として、心拍の開始位置が所定位置として設定されている。
従って、この時点では、表示処理部５２は、拍選択部４８によって自動選択された、最新の心拍の開始位置のＢモード画像（自動選択された心拍の開始時点に撮られたＢモード画像）を記憶部３６から読み出して、Ｍモード画像６５と共にディスプレイに表示する。
【００５３】
前述のように、通常の超音波診断装置では、Ｂモード画像とＭモード画像との同時表示において、フリーズボタンが押されて表示されるＢモード画像は、フリーズボタンが押された時点の画像である。
他方、超音波診断措置において、血管壁の変位量から血管弾性率を計測する場合のように、心拍に応じた対象部位の変動を用いる計測では、正確な計測を行なうためには、Ｍモード画像から解析に適正な心拍を選択し、かつ、選択した心拍に対応するＢモード画像を表示して、解析を行なう必要がある。
しかしながら、フリーズ時点のＢモード画像を表示する従来の装置では、必ずしも解析に用いる心拍に対応するＢモード画像が表示されるとは、限らない。また、特許文献３に示されるように、心電図等による特定時相のＢモード画像を表示する装置も提案されているが、やはり、Ｍモード画像中で選択された心拍に対応するとは、限らない。
【００５４】
これに対し、本発明の超音波診断装置１０においては、Ｍモード画像中において心拍を選択し、心拍の長さに対して、心拍の前２０％の時点から心拍の後２０％の時点（以下、単純に、『心拍の前２０％から心拍の後２０％』とも言う）までの期間内で、適宜、設定された所定位置に対応して、選択された心拍の所定位置のＢモード画像を表示する。
そのため、本発明によれば、Ｍモード画像中で選択された心拍に対応するＢモード画像を表示して、好適な解析や診断を行なうことができる。例えば、後述するように、Ｂモード画像を用いて血管壁の境界を設定し、Ｂモード画像の時相および設定した血管壁境界を開始点として、Ｍモード画像中で血管壁のトラッキング（追跡）等の解析を行なう際に、余分な時間分のトラッキングを行なうことなく、正確なトラッキングが可能となる。また、選択した心拍における心臓の拡張期や収縮期に対応するＢモード画像など、解析に好適なＢモード画像を表示することもできる。
【００５５】
本発明において、表示するＢモード画像６４が、選択された心拍の前２０％の時点よりも前、および、心拍の後２０％の時点よりも後では、表示するＢモード画像６４が、Ｍモード画像６５で選択された心拍と時間的に離れすぎてしまう。
その結果、選択された心拍に対して、表示されたＢモード画像６４の状態が大きく異なってしまう場合が有り、適正な解析や診断を行なうことが困難になってしまう。また、例えば、後述するように、Ｂモード画像６４を位置および時間的な開示位置としてＭモード画像中でトラッキングを行なう場合には、トラッキングの開始位置から選択された心拍に至るまでの距離（時間）が長くなり過ぎてしまう。その結果、トラッキングにエラーが生じ易くなってしまい、また、トラッキングの結果にスペックルに起因するノイズ等の不要な情報が多く拾われてしまう等の不都合が生じる。
【００５６】
Ｂモード画像６４を表示する所定位置は、心拍の前２０％から心拍の後２０％までの期間の中であれば、心拍の前（心拍開始前）１０％の位置、心拍の開始位置、心臓の収縮期すなわち血管最大径に対応する位置、心拍の中央、心拍の終了位置など、各種の位置が利用可能である。
ここで、心臓の拡張期ｄから収縮期ｓに至る期間は、図６（Ｃ）にも示されるように、血管壁の動きが激しい（血管壁の移動速度が速い）。そのため、この間のＢモード画像は、他の位置に比して、画質が良くない。また、後述するような血管壁のトラッキングを行なう場合に、このように動きが激しい場所を開始点とすると、トラッキングのエラーが発生し易い。
そのため、所定位置は、この心臓の拡張期ｄすなわち心拍の開始位置よりも後で、かつ、収縮期ｓよりも前の期間以外で設定するのが好ましい。
【００５７】
また、図６（Ｃ）にも示されるように、心臓の拡張期の前すなわち心拍の開始位置の前は、血管壁の移動は緩やかで、画質の良いＢモード画像が得られる。加えて、血管壁の弾性の計測など、血管壁の移動量を用いた計測を行なう場合には、多くの場合、解析に重要なのは、血管径が最小から最大に変化する心臓の拡張期から収縮期までの間である。
従って、Ｂモード画像６４を表示する所定位置は、特に、心拍の開始位置から、前２０％までの期間に設定するのが好ましい。すなわち、図６（Ｃ）に示す、心拍（心拍の開始位置を含む）の前０．２ｔの期間に設定するのが好ましい。
【００５８】
本発明の超音波診断装置１０において、Ｂモード画像６４を表示する所定位置は、デフォルトとして予め固定的に設定されていてもよく、もしくは、医師等の操作者が所定位置を選択して設定できるようにしてもよい。
また、操作者が所定位置を選択して設定する場合には、『心拍開始前１０％の位置』、『心拍開始位置』、『心臓の収縮期の位置』、『心拍の中央』、『心拍終了位置』のように、複数の位置を選択肢として設定しておき、操作者が選択して所定位置を設定できるようにしてもよい。あるいは、心拍の前２０％から後２０％までの期間（時間軸）において、任意の位置を選択して所定位置として設定できるようにしてもよい。あるいは、設定された選択肢から選択する方法と、時間軸上の任意の位置を選択する方法とを、選択できるようにしてもよい。
【００５９】
なお、本発明において、所定位置は１つに限定はされず、２以上の所定位置を設定してもよい。
例えば、血管壁の弾性率を簡易的に計測する方法として、１つの心拍における血管の最大径（Ｄｓ）と、最小径（Ｄｄ）との比（Ｄｄ／Ｄｓ）を算出する方法が知られている。これに対応して、血管が最小径となる心臓拡張期（心拍開始位置）と、血管が最大径になる心臓収縮期を所定位置として設定してもよい。
このように、複数の所定位置を設定する場合には、複数のＢモード画像を切り換えて表示するようにすればよい。あるいは、ディスプレイ１８に十分な表示スペースが有る場合には、Ｍモード画像と共に複数のＢモード画像を同時表示してもよい。
【００６０】
Ｍモード画像６５に心拍のラインが表示され、Ｂモード画像６５が選択された心拍に対応する所定位置の画像（前述のように、図示例では心拍開始位置の画像）となると、Ｂモード画像中の選択ライン６２が実線になって、トラックボールによって左右方向に移動可能となる。すなわち、選択ライン６２が、アクティブな状態になる。同時に、Ｍモード画像の心拍を示すラインは、全て、破線になる。なお、各種のラインがアクティブか否かの区別は、先と同様、線種に変えて、あるいは加えて、線の色を利用してもよい。
この状態で、トラックボールによって選択ライン６２を左右方向に移動すると、表示処理部５２は、記憶部３６から、選択ライン６２の位置に対応するＭモード画像を読み出し、心拍の検出結果と共に、その画像をディスプレイ１８に表示させる。すなわち、フリーズ後にも、トラックボールによって選択ライン６２を移動することで、Ｂモード画像６４中のアジマス方向の全域から、Ｂモード画像６４中におけるＭモード画像６５の表示位置（表示ライン）を選択できる。
従って、本例によれば、設定したＲＯＩ６０のアジマス方向の任意の位置のＭモード画像６５を表示して、Ｍモード画像６５、および、Ｍモード画像中の各心拍に対応する画像を、観察／確認することができる。
【００６１】
Ｂモード画像６４の選択ライン６２が移動可能な状態でセットボタンが押されると、Ｍモード画像の表示位置（表示ライン）の選択が終了したとして、図７（Ａ）に示されるように、Ｂモード画像６４の選択ライン６２が破線になって、トラックボールによる移動が不可能な状態となる。同時に、Ｍモード画像６５において、拍選択部４８によって自動選択された最新の心拍を示すラインが、共に実線となる。
【００６２】
Ｍモード画像６５において、最新の心拍を示すラインが、共に実線になると、トラックボールによって、心拍の選択が可能な状態になる。
セットボタンが押された時点では、前述のように最新の心拍が選択（拍検出部４８によって自動剣術）された状態になっている。これに応じて、Ｍモード画像６５の表示では、図７（Ａ）および（Ｂ）に示すように、最新の心拍を示すラインが実線となって選択された状態となっている。
この状態から、例えば、トラックボールが左に回されると、その信号が拍選択部４８に供給される。この回転量に応じて、拍選択部４８は、最新の心拍の１つ前の心拍が選択されたとして、その旨の情報を表示処理部５２に供給する。表示処理部５２は、この情報に応じて、図７（Ｃ）に示すように、最新の心拍の終了に対応するラインを破線とし、次に新しい心拍に対応するラインが実線として、この心拍が選択された状態となる。
さらにトラックボールが左に回されると、拍選択部４８は、その回転量に応じて、３番目に新しい心拍が選択されたとして、その旨の情報を表示処理部５２に供給し、２番目に新しい拍に対応するラインが破線になり、３番目に新しい心拍に対応するラインが実線となって、選択された状態となる。
また、トラックボールを右に回せば、同様に、順次、新しい心拍に対応するラインが選択された状態となる。
【００６３】
さらに、トラックボールの移動によって、選択された心拍が変更された場合には、表示処理部５２は、選択された心拍の情報に応じて、この心拍の前記所定位置に対応するＢモード画像を記憶部３６から読み出して、ディスプレイ１８に表示する。すなわち、Ｂモード画像６４を、新しく選択された心拍の所定位置の画像に変更する。
図示例においては、所定位置として、心拍の開始位置が設定されている。従って、例えば、図７（Ｃ）に示すように、トラックボールによって２番目に新しい心拍が選択された場合には、表示処理部５２は、この心拍の開始位置（開始時点）のＢモード画像を記憶部３６から読み出し、ディスプレイ１８にＢモード画像６４として表示する。
【００６４】
心拍の選択が可能な状態でセットボタンが押されると、心拍の選択が終了したとして、選択された心拍が確定し、選択した心拍の微調整が行なえる状態となる。
また、ディスプレイ１８に表示しているＭモード画像６５中の心拍が選択／確定すると、拍選択部４８は、記憶部３６が記憶している全てのＭモード画像（すなわち、Ｂモード画像６４のアジマス方向の全域のＭモード画像）において、必要に応じて、選択された心拍の情報を変更する。すなわち、図示例においては、最新の心拍以外が選択されて確定された場合には、選択された心拍の情報を新たに選択／確定された心拍に変更する。
また、選択された心拍が確定すると、選択された心拍の情報は、トラッキング部４２に供給される。
【００６５】
一例として、最新の心拍が選択（すなわち、自動選択された心拍からの変更は無し）されたとして、セットボタンが押されると、図８（Ａ）に示すように、まず、選択された心拍の終了に対応する線が細線となり、選択された心拍の開始に対応する線の位置（時間）が、トラックボールによって、矢印ｔで示すように、左右方向（時間方向）に移動可能となり、心拍の開始位置の微調整が可能な状態となる。
必要に応じて、トラックボールによって心拍の開始位置が調整された後、再度、セットボタンが押されると、今度は、図８（Ｂ）に示すように、選択された心拍の終了に対応する線が通常の実線に、同開始に対応する線が細線になって、トラックボールによって、矢印ｔで示すように、選択された心拍の終了に対応する線の位置が、左右方向に移動可能となり、心拍の終了位置の微調整が可能な状態となる。
この心拍の微調整の結果は、微調整を行ったＭモード画像６５のみに反映してもよいが、記憶部３６に記憶されている全てのＭモード画像にも反映されるのが好ましい。
なお、心拍の開始位置が調整された場合には、表示処理部５２は、記憶部３６から、調整された心拍開始位置のＢモード画像を読み出し、ディスプレイ１８に表示するＢモード画像６４を、この画像に変更する。
【００６６】
心拍の微調整の結果は、表示処理部５２にも供給される。表示処理部５２は、心拍の微調整に応じて、選択された心拍中における所定位置が変動した場合には、変動した所定位置に応じたＢモード画像を記憶部３６から読み出し、ディスプレイ１８に表示する。
図示例においては、心拍の開始位置が所定位置として設定されている。従って、必要に応じて心拍の開始位置が微調整されたら、表示処理部５２は、調整された心拍開始位置に対応するＢモード画像を記憶部３６から読み出して、ディスプレイ１８に表示する。
また、心拍の微調整の結果が、記憶部３６に記憶されている全てのＭモード画像にも反映される場合には、この心拍の微調整に応じた所定位置の変動は、記憶部３６に記憶されている全てのＭモード画像にも反映されるのが好ましい。
【００６７】
選択された心拍の終了位置が調整可能な状態で、セットボタンが押されると、前述の図６に示すＢモード画像６４の選択ライン６２が移動可能な状態、すなわち、Ｂモード画像６４中における、Ｍモード画像６５の表示ラインの選択が可能な状態に戻る。
すなわち、図示例の超音波診断装置１０では、「表示ラインの選択」→「心拍の選択」→「心拍の微調整」の各処理を繰り返し行うことができるようになっている。言い換えれば、「表示ラインの選択」→「心拍の選択」→「心拍の微調整」は、ループ状に処理を行なうことが可能になっている。
これにより、より好適に、記憶する全てのＭモード画像から、後述する血管壁の弾性計測のための解析に最適な心拍を選択することが可能になる。
【００６８】
一方、選択された心拍の終了に対応する位置が調整可能な状態で、セットボタンではなく、タッチパネルの「ＡＷ Ｄｅｔ」ボタンが押されると、図９に示すように、Ｂモード画像６４の選択ライン６２、および、Ｍモード画像６５中における心拍を示すラインが、全て、破線となって操作ができない状態となり、血管壁の検出モードとなる。
【００６９】
血管壁の検出モードになると、まず、図１０（Ａ）に示すように、Ｂモード画像６４中に、血管前壁の外膜中膜境界に対応するライン６８が表示される。
このライン６８は、トラックボールによって上下方向（深度方向）に平行移動可能になっている。図１０（Ｂ）に示すように、トラックボールによって移動して、ライン６８を血管前壁の外膜中膜境界の位置に移動したら、セットボタンが押される。
【００７０】
セットボタンが押されると、図１０（Ｃ）に示すように、Ｂモード画像６４において、血管前壁の外膜中膜境界に対応するライン６８が破線になって確定し、血管前壁の内膜内腔境界に対応するライン７０が表示される。
このライン７０も、同様に、トラックボールによって上下方向に移動可能であり、ライン７０を血管前壁の中膜内腔境界の位置に移動したら、セットボタンが押される。
【００７１】
ライン７０が移動可能な状態でセットボタンが押されると、図１０（Ｄ）に示すように、Ｂモード画像６４において、血管前壁の中膜内腔境界に対応するライン７０が破線になって確定し、血管後壁の内膜内腔境界に対応するライン７２が表示される。同様に、トラックボールによってライン７２を血管後壁の中膜内腔境界の位置に移動したら、セットボタンが押される。
さらに、ライン７２が移動可能な状態でセットボタンが押されると、図１０（Ｅ）に示すように、Ｂモード画像６４において、血管後壁の内膜内腔境界に対応するライン７２が破線になって確定し、血管後壁の外膜中膜境界に対応するライン７４が表示される。同様に、トラックボールによってライン７４を血管後壁の外膜中膜境界の位置に移動したら、セットボタンが押される。
【００７２】
血管壁の各境界の情報は、境界検出部４０に供給される。
ライン７４が移動可能な状態でセットボタンが押されると、全ての境界に対応するラインの設定が終了し、境界検出部４０は、設定された内膜内腔境界のライン７２および外膜中膜境界のライン７４を用いて、後壁の内膜内腔境界および外膜中膜境界の自動検出を行なう。両境界の自動検出の結果は、表示処理部５２およびトラッキング部４２に送られ、図１０（Ｆ）に示すように、検出結果が表示される。
なお、これらの境界の自動検出の方法には、特に限定はなく、各種の方法が利用可能である。一例として、Ｂモード画像を解析して、ライン７２およびライン７４の位置において連続する高輝度な部分をトレースして、内膜内腔境界および外膜中膜境界を検出する方法が例示される。
【００７３】
境界検出部４０による、血管後壁の内膜内腔境界および外膜中膜境界の自動検出が終了すると、図１０（Ｆ）に示すように、Ｂモード画像６４にカーソル７８が表示される（血管後壁の自動検出が終了するまでは、このカーソル７８は表示されていない）。
【００７４】
このカーソル７８は、トラックボールによって移動可能になっている。カーソル７８を、自動検出された内膜内腔境界もしくは外膜中膜境界を示すラインに移動して、セットボタンが押されると、カーソル７８に近い側のラインが実線となる。実線となったラインは、修正が可能な状態となる。
一例として、図１０（Ｇ）に示すように、外膜内膜境界を示すライン７４が選択されて実線になったとする。ライン７４に沿って、トラックボールによってカーソル７８を移動して、再度、セットボタンが押されると、境界検出部４０が、カーソルによってなぞられた領域のライン７４を再検出して書き換えが行なわれ、また、その結果がトラッキング部４２に送られる。
【００７５】
後壁の内膜内腔境界および外膜中膜境界の自動検出が終了し、さらに、必要に応じて血管後壁の修正が行なわれると、図１１（Ａ）に示すように、全てのラインが破線になっている状態になり、さらに、図１１（Ｂ）に示すように、タッチパネル１６ａの「Elasticity Ana」ボタンが選択可能な状態となる。
「Elasticity Ana」ボタンが選択可能な状態となったら、「Ｐｓ」ボタンを用いて被検者の心臓収縮期の血圧を、「Ｐｄ」ボタンを用いて被検者の心臓拡張末期の血圧を、それぞれ入力し、さらに、「Quality Factor Threshold」ボタンを用いて、信頼性閾値を入力する。これらの数値の入力は、公知の方法で行なえばよい。
【００７６】
なお、被検者の血圧および信頼性閾値の入力は、血管壁境界の検出が終わった後に行なうのに限定はされず、後述する解析開始の前（後述する「Elasticity Ana」ボタンの押下の前）であれば、どのタイミングで行なってもよい。
また、超音波診断装置１０においては、診断を行なう前に、被検者情報の取得や入力を行なうのが通常であるが、この被検者情報に血圧の情報が有る場合には、これを利用してもよい。
【００７７】
被検者の血圧および信頼性閾値が入力され、「Elasticity Ana」ボタンが押されると、Ｂモード画像の解析が開始され、血管壁の弾性率の計算が行なわれる。
「Elasticity Ana」ボタンが押されると、まず、トラッキング部４２が、Ｍモード画像６５において、選択された心拍における血管前壁（外膜中膜境界および中膜内腔境界）ならびに血管後壁（中膜内腔境界および外膜中膜境界）の動きを追跡する。すなわち、血管前壁および後壁のトラッキングを行なう。
【００７８】
Ｍモード画像６５における血管壁のトラッキングは、先にＢモード画像６４において検出（設定）した、血管前壁の外膜中膜境界、血管前壁の中膜内腔境界、血管後壁の中膜内腔境界、および、血管後壁の外膜中膜境界を、位置的な出発点（焦点深度方向の出発点）として行なう。
他方、Ｍモード画像６５における血管壁のトラッキングの時間的な出発点は、ディスプレイ１８に表示した、選択された心拍の前記所定位置のＢモード画像６４に対応する位置（時相）である。すなわち、Ｍモード画像６５中におけるＢモード画像６４が撮られた時点が、トラッキングの出発点となる。従って、図示例においては、選択された心拍の開始位置が、トラッキングの出発点となる。
【００７９】
前述のように、血管壁の境界を指定したＢモード画像６４は、心拍の前２０％から後２０％の期間内で設定された、選択された心拍の所定位置に対応するＢモード画像である。
従って、図示例においては、トラッキングの出発点となるＢモード画像６４が、選択された心拍あるいは選択された心拍に近いので、トラッキングのエラーや、余分なノイズ等の影響を受けることなく、Ｍモード画像６５において、選択された心拍における血管壁のトラッキングを行なうことができる。
【００８０】
ここで、超音波診断装置１０においては、好ましい態様として、検出（設定）した血管壁の境界のみならず、血管後壁の中において、深さ方向に１以上の測定点が設定されていてもよい。このように、血管後壁の中に、１以上の測定点が設定されている場合には、各測定点毎に、血管壁のトラッキングを行なう。
なお、血管壁中の測定点は、予め設定されていてもよく、あるいは、特定のアルゴリズムに基づいて自動設定してもよく、あるいは、超音波診断装置１０の操作者が画像を見ながら設定してもよく、これらを併用してもよい。
【００８１】
Ｍモード画像６５における血管壁のトラッキングの方法には、特に限定はなく、トラッキングの出発位置からの画像（輝度）の連続性を利用する方法、パターンマッチング法、ゼロクロス法、組織ドプラ法、位相差トラッキング等が例示され、いずれの方法を用いてもよい。
【００８２】
トラッキング部４２によるＭモード画像中での血管壁のトラッキング結果は、弾性率算出部５０および表示処理部５２に供給される。
弾性率算出部５０は、血管壁のトラッキング結果から、まず、血管壁（内膜中膜）の厚みの変化波形、および、血管径（内径）の変化波形を形成する。なお、前述のように、血管壁内に１点以上の測定点が設定されている場合には、血管壁の変化波形は、各測定点の間毎に形成される。
血管壁の厚みの変化波形、および、血管径の変化波形は、表示処理部５２に送られる。
【００８３】
また、弾性率算出部５０は、下記式（１）を用いて、血管の径方向の歪みを算出する。
εi＝Δｈi／ｈdi ・・・（１）
なお、上記式（１）において、εは、各測定点の間における血管の径方向の歪みを、Δｈiは、１つの心拍内で血管壁が最も薄くなる心臓収縮期における各測定点の間の血管壁の厚み変化の最大値を、ｈdiは、血管壁が最も厚くなる心臓拡張末期における各測定点の間の厚みを、それぞれ示す。
【００８４】
さらに、弾性率算出部５０は、先に入力された血圧の最高値および最低値を用いて、下記式（２）によって、血管壁の円周方向の弾性率Ｅθiを算出する。
Ｅθi＝１／２＊［１＋（ｒd／ｈd）］＊［Δｐ／（Δｈi／ｈdi）］） ・・・（２）
もしくは、下記式（３）によって、血管壁の径方向の弾性率Ｅriを算出してもよい。
Ｅri＝Δｐ／（Δｈi／ｈdi） ・・・（３）
なお、上記式（２）および式（３）において、Δｈiおよびｈdiは、先と同様であり、Δｐは、心臓収縮期と心臓拡張末期とにおける血圧差を、ｒdは、心臓拡張末期での血管内腔の半径を、ｈdは、心臓拡張末期における血管壁の厚みを、それぞれ示す。
【００８５】
弾性率を計算したら、弾性率算出部５０は、弾性率の信頼性を計算する。
弾性率の信頼性の計算方法には、特に限定はなく、公知の方法が、各種、利用可能である。一例として、１０００人等の多数の人物の心拍による血管径変化の波形を作成して、これらの多数の波形から、血管径変化のモデル波形を作成し、このモデル波形からのズレの量を用いて、算出した弾性率の信頼性を計算する方法が例示される。
【００８６】
ここで、前述のように、ディスプレイ１８に表示されているＭモード画像６５で心拍が選択／確定すると、記憶部３６が記憶している全てのＭモード画像において、同じ心拍が選択された状態となる。
これに応じて、上述の血管壁のトラッキング、血管壁の厚みおよび血管径の変化波形の作成、血管壁の歪みの計算、血管壁の弾性率および弾性率の信頼性の計算などの処理は、ディスプレイ１８に表示されているＭモード画像６５のみならず、記憶部３６に記憶されている全てのＭモード画像において、選択された心拍に対して行なわれる。すなわち、ディスプレイ１８に表示されているＢモード画像６４のアジマス方向の全領域に対して、対応するＭモード画像を用いて、選択された心拍での上記血管壁の弾性率算出などの処理が行なわれる。
これらの結果は、記憶部３６に記憶されているＭモード画像に情報として付加される。
【００８７】
アジマス方向の全域における演算が終了したら、弾性率算出部５０は、血管壁の弾性率の平均値（Ｅθave）、血管壁の歪みの平均値（Ｓｔｒave）、および、弾性率の信頼性の平均値（ＱＦave）を、演算する。
【００８８】
演算が終了すると、ディスプレイ１８に結果が表示される。
その一例を、図１２に示す。図示例においては、元々、表示されていたＢモード画像６４の図中右側に、このＢモード画像６４に示される血管後壁の弾性率が、Ｂモード画像６４ｅで表示される。さらに、この血管後壁の弾性率を表示するＢモード画像６４ｅの図中右側に、算出した血管壁の弾性率の信頼性が、同じくＢモード画像６４ｑで表示される。
また、Ｂモード画像６４の図中左側に、血管壁の弾性率の平均値（Ｅθave）、血管壁の歪みの平均値（Ｓｔｒave）、および、弾性率の信頼性の平均値（ＱＦave）が、それぞれ、表示される。
【００８９】
血管壁の弾性率は、Ｂモード画像６４ｅにおいて、Ｂモード画像６４で自動検出（あるいは必要に応じて修正）された血管後壁に重ねて、帯状に表示される。また、Ｂモード画像６４ｅの右上側には、弾性率の指標が表示される。図示例においては、画像が高密度であるほど、弾性率が高い。
すなわち、Ｂモード画像６４ｅにおいては、血管後壁に重なる帯の密度が、血管のその位置における血管壁の弾性率を示している。
【００９０】
弾性率の信頼性は、Ｂモード画像６４ｑにおいて、同様に、Ｂモード画像６４中で自動検出された血管後壁に重ねて、帯状に表示される。また、Ｂモード画像６４ｑの右上側には、弾性率の信頼性の指標が表示される。図示例においては、画像が高密度であるほど、弾性率の信頼性が高い。
すなわち、Ｂモード画像６４ｑにおいては、血管後壁に重なる帯の密度が、血管のその位置における血管壁弾性率の信頼性を示している。
【００９１】
なお、この弾性率や弾性率の信頼度の高低は、画像の密度に変えて、あるいは加えて、画像の色によって表現してもよい。
【００９２】
ここで、図１２に示す結果の表示では、先に入力した閾値よりも信頼性が低い、アジマス方向の位置では、結果が自動的に省かれる。
また、結果が省かれた位置に関しては、Ｂモード画像６４ｅにおける弾性率の結果表示の右隅部や、Ｂモード画像６４ｑにおける信頼性の結果表示の右隅部に示されるように、帯の表示が薄くなる。
【００９３】
また、下側のＭモード画像６５においては、選択された心拍に、Ｍモード画像中での血管前壁のトラッキング結果８０および血管後壁のトラッキング結果８２、血管径の変化波形８４、ならびに、血管壁の厚みの変化波形８６が表示される。
なお、前述のように、血管壁の中に深度方向に１以上の測定点が設定されている場合には、血管厚みの変化波形は、各測定点の間毎に、出力するようにしてもよい。
【００９４】
ここで、血管壁の弾性率の計測結果等がディスプレイ１８に表示されると、Ｂモード画像６４における選択ライン６２が実線になって、トラックボールによってアジマス方向に移動可能になる。
Ｂモード画像６４で選択ライン６２が移動されると、表示処理部５２は、記憶部３６から、選択ライン６２の位置対応するＭモード画像を読み出して、ディスプレイ１８に表示する。すなわち、トラックボールによって選択ライン６２が移動されると、Ｍモード画像６５が、この選択ライン６２の位置のＭモード画像に変わり、Ｍモード画像中での血管前壁のトラッキング結果８０および血管後壁のトラッキング結果８２、血管径の変化波形８４、ならびに、血管壁の厚みの変化波形８６が、Ｂモード画像６４の選択ライン６２の位置のデータに変更される。
従って、Ｂモード画像中のアジマス方向の全域において、Ｍモード画像６５および解析結果を表示する表示ラインを、選択することができる。
【００９５】
また、このＢモード画像６４での選択ライン６２の移動に同期して、Ｂモード画像６４ｅの選択ライン６２ｅ、および、Ｂモード画像６４ｑの選択ライン６２ｑも、移動する。
【００９６】
さらに、セットボタンを押した後、Ｂモード画像６４ｅおよびＢモード画像６４ｑにおいて、選択ライン６２ｅおよび選択ライン６２ｑをトラックボールによって移動してアジマス方向の任意の領域を選択し、その後、再度、セットボタンを押すと、選択された領域は前述の信頼性が閾値よりも低かった領域と同様の扱いとなり、データが取り除かれる。
すなわち、検査者が結果を見て、波形等が奇怪しいと感じられる場所が有る場合に、そのデータを取り除くことができ、より正確な解析が可能となる。
【００９７】
なお、このデータの除去は、Ｄｅｌｅｔボタン等の押下によって、１つ前の状態にも土切要にしてもよい。
【００９８】
以上、本発明の超音波診断装置について詳細に説明したが、本発明は、上述の例に限定はされず、本発明の要旨を逸脱しない範囲において、各種の変更や改良を行なってもよいのは、もちろんである。
【産業上の利用可能性】
【００９９】
本発明の超音波診断装置は、心筋梗塞、狭心症、脳疾患等の原因となる動脈硬化の診断を行なう医療現場等において、好適に利用可能である。
【符号の説明】
【０１００】
１０ 超音波診断装置
１２ 診断装置本体
１４ （超音波）プローブ
１６ 操作パネル
１８ ディスプレイ
２０ ケーブル
２４ キャスタ
２８ 送信回路
３０ 受信回路
３２ 送受信制御部
３４ 画像形成部
３６ 記憶部
４０ 境界検出部
４２ トラッキング部
４６ 拍検出部
４８ 拍選択部
５０ 弾性率算出部
５２ 表示処理部
５６ Ｂモード画像形成部
５８ Ｍモード画像形成部
６０ ＲＯＩ
６２ 選択ライン
６４，６４ｅ，６４ｑ，９０ Ｂモード画像
６５，９２ Ｍモード画像
６８，７０，７２，７４ ライン
８０，８２ トラッキング結果
８４ 血管径変化波形
８６ 血管壁厚み変化波形
【特許請求の範囲】
【請求項１】
超音波を送信し、被検体によって反射された超音波エコーを受信して、受信した超音波エコーに応じた受信信号を出力する超音波トランスデューサを有する超音波プローブと、
前記超音波トランスデューサが出力した受信信号から、Ｂモード画像およびＭモード画像を形成する、画像形成手段と、
表示手段と、
前記画像形成手段が形成したＢモード画像およびＭモード画像の少なくとも一方を、前記表示手段に表示させる表示処理手段と、
前記Ｂモード画像およびＭモード画像が表示された状態において、表示画像を静止画にするフリーズ手段と、
前記Ｂモード画像およびＭモード画像が静止画にされた状態において、Ｍモード画像中で心拍を選択する心拍選択手段とを有し、
かつ、前記表示処理手段は、前記Ｍモード画像中で心拍が選択されたら、心拍の長さに対して心拍の前２０％から心拍の後２０％までの期間内で予め設定された所定位置に対応して、前記選択された心拍における前記所定位置のＢモード画像を前記表示手段に表示させることを特徴とする超音波診断装置。
【請求項２】
前記Ｍモード画像中の心拍を検出する拍検出手段を有し、
前記表示処理手段は、前記拍検出手段が心拍を検出したら、Ｍモード画像中に、前記拍検出手段が検出した全ての心拍を表示する請求項１に記載の超音波診断装置。
【請求項３】
前記拍検出手段は、前記Ｍモード画像を解析することにより、前記Ｍモード画像中の心拍を検出する請求項２に記載の超音波診断装置。
【請求項４】
前記所定位置を選択する選択手段を有する請求項１〜３のいずれかに記載の超音波診断装置。
【請求項５】
前記心拍選択手段は、前記拍検出手段が検出した心拍の内、開始から終了までの全てがＭモード画像中で取得されている、最も新しい心拍が選択されたとする請求項２〜４のいずれかに記載の超音波診断装置。
【請求項６】
前記所定位置は、心臓の拡張期よりも後で、かつ、心臓の収縮期よりも前の期間以外に設定される請求項１〜５のいずれかに記載の超音波診断装置。
【請求項７】
前記心拍選択手段によって選択された心拍の位置を調整する位置調整手段を有する請求項１〜６のいずれかに記載の超音波診断装置。
【請求項８】
前記位置調整手段によって心拍の位置調整が行なわれたら、前記表示処理手段は、この心拍の位置調整および前記所定位置に応じて、前記表示手段に表示するＢモード画像を対応する位置の画像に変更する請求項７に記載の超音波診断装置。
【請求項９】
前記Ｂモード画像中において血管壁の境界位置を設定する、血管壁境界設定手段を有する請求項１〜８のいずれかに記載の超音波診断装置。
【請求項１０】
表示手段に表示されたＢモード画像中で、関心領域を設定するＲＯＩ設定手段を有する請求項１〜９のいずれかに記載の超音波診断装置。
【請求項１１】
前記関心領域の設定指示に応じて、前記超音波トランスデューサによる超音波のフレームレートを、設定指示前よりも高くする請求項１０に記載の超音波診断装置。
【請求項１】
超音波を送信し、被検体によって反射された超音波エコーを受信して、受信した超音波エコーに応じた受信信号を出力する超音波トランスデューサを有する超音波プローブと、
前記超音波トランスデューサが出力した受信信号から、Ｂモード画像およびＭモード画像を形成する、画像形成手段と、
表示手段と、
前記画像形成手段が形成したＢモード画像およびＭモード画像の少なくとも一方を、前記表示手段に表示させる表示処理手段と、
前記Ｂモード画像およびＭモード画像が表示された状態において、表示画像を静止画にするフリーズ手段と、
前記Ｂモード画像およびＭモード画像が静止画にされた状態において、Ｍモード画像中で心拍を選択する心拍選択手段とを有し、
かつ、前記表示処理手段は、前記Ｍモード画像中で心拍が選択されたら、心拍の長さに対して心拍の前２０％から心拍の後２０％までの期間内で予め設定された所定位置に対応して、前記選択された心拍における前記所定位置のＢモード画像を前記表示手段に表示させることを特徴とする超音波診断装置。
【請求項２】
前記Ｍモード画像中の心拍を検出する拍検出手段を有し、
前記表示処理手段は、前記拍検出手段が心拍を検出したら、Ｍモード画像中に、前記拍検出手段が検出した全ての心拍を表示する請求項１に記載の超音波診断装置。
【請求項３】
前記拍検出手段は、前記Ｍモード画像を解析することにより、前記Ｍモード画像中の心拍を検出する請求項２に記載の超音波診断装置。
【請求項４】
前記所定位置を選択する選択手段を有する請求項１〜３のいずれかに記載の超音波診断装置。
【請求項５】
前記心拍選択手段は、前記拍検出手段が検出した心拍の内、開始から終了までの全てがＭモード画像中で取得されている、最も新しい心拍が選択されたとする請求項２〜４のいずれかに記載の超音波診断装置。
【請求項６】
前記所定位置は、心臓の拡張期よりも後で、かつ、心臓の収縮期よりも前の期間以外に設定される請求項１〜５のいずれかに記載の超音波診断装置。
【請求項７】
前記心拍選択手段によって選択された心拍の位置を調整する位置調整手段を有する請求項１〜６のいずれかに記載の超音波診断装置。
【請求項８】
前記位置調整手段によって心拍の位置調整が行なわれたら、前記表示処理手段は、この心拍の位置調整および前記所定位置に応じて、前記表示手段に表示するＢモード画像を対応する位置の画像に変更する請求項７に記載の超音波診断装置。
【請求項９】
前記Ｂモード画像中において血管壁の境界位置を設定する、血管壁境界設定手段を有する請求項１〜８のいずれかに記載の超音波診断装置。
【請求項１０】
表示手段に表示されたＢモード画像中で、関心領域を設定するＲＯＩ設定手段を有する請求項１〜９のいずれかに記載の超音波診断装置。
【請求項１１】
前記関心領域の設定指示に応じて、前記超音波トランスデューサによる超音波のフレームレートを、設定指示前よりも高くする請求項１０に記載の超音波診断装置。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【公開番号】特開２０１２−２４９８４８（Ｐ２０１２−２４９８４８Ａ）
【公開日】平成２４年１２月２０日（２０１２．１２．２０）
【国際特許分類】
【出願番号】特願２０１１−１２５２８３（Ｐ２０１１−１２５２８３）
【出願日】平成２３年６月３日（２０１１．６．３）
【出願人】（３０６０３７３１１）富士フイルム株式会社 (25,513)
【Ｆターム（参考）】
【公開日】平成２４年１２月２０日（２０１２．１２．２０）
【国際特許分類】
【出願日】平成２３年６月３日（２０１１．６．３）
【出願人】（３０６０３７３１１）富士フイルム株式会社 (25,513)
【Ｆターム（参考）】
[ Back to top ]