JP2017158917A - 医用診断装置、および医用解析プログラム - Google Patents

Abstract

【課題】被検体内における組織性状を解析し、被検体内における非構造物の領域における組織性状を表示可能な医用診断装置の提供。【解決手段】本実施形態に係る医用診断装置は、被検体に対するスキャンにより得られたデータに基づいて、前記被検体内の構造物の位置を推定する推定部と、前記データに基づいて、前記推定部により推定された前記構造物の位置を除く前記被検体内の複数位置の組織性状を解析する解析部と、前記解析部による解析結果を表示する表示部と、を備える。【選択図】 図１

Description

本発明の実施形態は、生体内をスキャンして臓器等の断層像を画像化し、疾患などの診断に関する医用診断装置、および医用解析プログラムに関わるものである。

超音波診断装置は、生体内情報の画像を表示する診断装置である。超音波診断装置は、Ｘ線診断装置やＸ線コンピュータ断層撮影（Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ：ＣＴ）装置などの他の画像診断装置に比べて、安価、無被曝かつ非侵襲性であって、実時間で観測するための有用な医用診断装置として利用されている。超音波診断装置の適用範囲は広い。超音波診断装置は、例えば、心臓などの循環器、肝臓および腎臓などの腹部、抹消血管、産婦人科、乳癌等の診断などに適用されている。

超音波診断装置は、一般的には超音波受信信号（エコー信号）の振幅の大きさを輝度で表現することで、生体組織の形態を映像化している。しかしながら、超音波の受信信号には、他の様々な物理的情報が含まれることが、種々の研究レポートで報告されている。超音波の受信信号に含まれる物理的情報のいくつかは、臨床応用に向けて種々の試みがなされている。

例えば、エコー信号の振幅の統計量を計算し、例えば平均値と分散値との関係を分析することで、視覚的には判断しにくい微小構造物の含有量などを定量化することができる。また、近年、被検体内の臓器の局所的な動き量を分析することで、該臓器の硬さあるいは弾性率といった物理情報を提示するいわゆる超音波エラストグラフィー法も用いられるようになってきている。上記超音波エラストグラフィー法も、画像化される前に超音波信号に含まれる位相情報などを利用する方法である。

また、生体組織は固有の減衰特性を持っている。被検体に照射された超音波は、減衰しながら生体内を伝搬する。このとき、生体内を伝搬した超音波の減衰量が大きい場合、スキャンの途中で十分なエコー信号を受信できなくなる現象が発生する。一方で、エコー信号の減弱の様子を観察することで、生体組織の特徴を観測することもしばしば行われている。例えば、送受信方向のエコー輝度の変化を解析して、対象物体の超音波減衰量を定量する手法が知られており、例えば肝臓を例にとれば、特に脂肪肝の定量診断において有用性が期待されている。具体的には、エコー信号が極端に減少する被検体は、肝臓内に脂肪滴を多く含む脂肪肝であることが推測される。肝硬変の場合においても同様な結果となることがある。

そこで、超音波の減衰量を定量的に診断するための手法が複数提案されている。例えば、中心周波数が異なる複数の超音波パルスを送受信し、取得された複数の信号の強度が深さ方向にどの程度変化するかを比較する。この比較により、被検体に固有の減衰量を推定するという方法がある。生体内における超音波の減衰量は、周波数によって異なることが知られている。このため、複数の周波数信号の強度の変化を比較することによって、対象としている組織に固有の値が求められる。また、広帯域パルスを用いれば、超音波の送受信方向１方向につき１回超音波の送受信することで、上記と同様な効果を得ることができる。

しかしながら、この方法では低周波数領域の信号が組織内を伝播する際に高周波（高調波）成分の信号が発生するため、発生した高周波成分は、被検体に固有の減衰量を推定する際の誤差となってしまうことが考えられる。そこで、一つの超音波送信方向に対して、送信超音波の波形の正負を反転させた２つのパルスを被検体に送信し、得られた受信信号の差分演算を行うことによって、超音波の伝搬中に発生した高周波成分を除去する方法も提案されている。

さらに、上記得られた複数の周波数成分の信号を用いて、減衰量の大きさを反映したカラー表示を行うことも提案されている。例えば、複数の周波数成分にそれぞれ対応する減衰量に異なった色相が割り当てられて、Ｂモード画像に重畳表示される。これにより、深さ方向の色相変化の大きさを見ることによって、ユーザは、視覚的に減衰量の大きさを精度よく把握することができる。また、異なる２周波成分の信号強度差に応じてカラーマッピングすることによっても同様の効果が得られる。さらに、２周波成分の信号強度差を深さ方向に微分することにより断面内の各点における減衰定数が求められるため、この減衰定数の大きさをカラー化して表示することができる。

特開２００１−２３８８８４号公報 特公平３−２４８６８号公報 特公平３−６０４９３号公報 特開平７−５１２７０号公報 特開２０１０−２３３８５９号公報 特許第４３７３６９８号公報 特許第４９２１８２６号公報 特公平７−３２７７６号公報

従来技術に示した方法は、いずれも受信信号の信号強度や周波数特性の深さ方向の変化を計算することによって減衰量を推定している。減衰量の推定は、深さ方向に散乱体の組成や分布が一様であることを仮定している。しかしながら、実際の生体は、いくつもの異なる組織が混在することにより、複雑な構造を有している。このため、上述したような一様な散乱体で成り立っている部分は少ない。たとえば、肝臓においては、数センチメートルの範囲では一様なスペックルのように見える箇所があっても、着目範囲を広げると、血管や腹壁や胆のうなど、目視で構造物が認識できる程度の構造物が、着目領域内に入ってくる。

このような構造物においては、信号の輝度が周囲の一様な実質部分と大きく変わるだけではなく、反射特性も異なる。このため、受信信号に含まれる周波数成分が大きく異なると考えられる。したがって、輝度の変化を解析する手法でも、また、周波数成分を解析する方法でも、上述したような周囲と特性が異なる構造物が解析対象領域内に含まれている場合、減衰量の推定の前提とする仮定が成り立たなくなる。この結果、解析対象領域内において、何らかの数値は算出できるが、計算された値は信頼性が低いものとなってしまう問題がある。

すなわち、上記いずれの方法においても、目立った構造物が解析対象領域内に含まれてしまうと、構造物付近での輝度分布や周波数特性が周囲と異なってしまい，減衰量を正しく算出できない場合がある。

たとえば、図１１に示すように、腹壁や横隔膜のように肝臓の実質とは異なる反射特性および散乱特性を示す部分を解析対象の範囲内に含む領域では、極端に高いあるいは低い減衰量の値が算出される。このとき、算出された減衰量の値は、対象物体の減衰量を表した結果ではなく、構造物によって発生したアーチファクトであると考えられる。

以上のことから、精度のよい減衰量の定量のためには，臨床現場において、できるだけ構造物の目立たない断面を設定し、かつ設定した断面の中でも一様な領域を選んで解析対象領域を設定するようにするという手順が求められる。しかしながら、解析対象領域の設定の判断は、操作者の主観に任されている。このため、解析対象領域の設定の精度や再現性が担保できない問題がある。加えて、解析に適した領域を選択することは、時に困難であるばかりか時間もかかり、操作者にとって負担となる問題がある。すなわち、操作者が構造物を避ける断面を描出したり、定量する際に定量結果が信頼できそうな一様な領域を選択したりすることは、操作者にとって負担であり、精度のよい診断が困難となることがある。

本実施形態の目的は、被検体内における組織性状を解析し、被検体内における非構造物の領域における組織性状を表示可能な医用診断装置および医用解析プログラムを提供することにある。

本実施形態に係る医用診断装置は、被検体に対するスキャンにより得られたデータに基づいて、前記被検体内の構造物の位置を推定する推定部と、前記データに基づいて、前記推定部により推定された前記構造物の位置を除く前記被検体内の複数位置の組織性状を解析する解析部と、前記解析部による解析結果を表示する表示部と、を備える。

図１は、本実施形態に係る医用診断装置としての超音波診断装置の構成を示す構成図である。 図２は、本実施形態に係り、分散比Ｒ_σの値に応じた色相で分散比Ｒ_σが重畳されるＢモード画像である。 図３は、本実施形態に係り、分散比Ｒ_σに関する被走査領域と同一の被走査領域における図２のＢモード画像に、分散比Ｒ_σを、分散比Ｒ_σの値に応じた色相で重畳した分散比重畳画像の一例を示す図である。 図４は、本実施形態に係り、モニタに表示された減衰重畳画像を凡例とともに示す図である。 図５は、本実施形態に係り、構造物推定機能および組織性状解析機能に関する処理手順の一例を示すフローチャートである。 図６は、本実施形態の変形例に係り、構造物推定機能および組織性状解析機能に関する処理手順の一例を示すフローチャートである。 図７は、本実施形態の第１の応用例に係り、計測ＲＯＩ１における代表値の計測が適切であると判定され、計測ＲＯＩ２における代表値の計測が不適であると判定された場合において、計測ＲＯＩの名称に対応する代表値および所定の通知と、計測ＲＯＩを重畳した減衰重畳画像とがモニタに表示される一例を示す図である。 図８は、本実施形態の第１の応用例に係り、代表値計算機能に関する処理手順の一例を示すフローチャートである。 図９は、本実施形態の第２の応用例に係り、計測ＲＯＩの名称に対応する代表値および信頼度のランクと、計測ＲＯＩが重畳された減衰重畳画像とがモニタに表示される一例を示す図である。 図１０は、本実施形態の第２の応用例に係り、信頼度表示機能に関する処理手順の一例を示すフローチャートである。 図１１は、従来技術に関する図面である。

以下、図面を参照しながら本実施形態に係わる医用診断装置について説明する。本実施形態に係る医用診断装置は、例えば、非侵襲性を有する医用診断装置（例えば、磁気共鳴イメージング（ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ ｉｍａｇｉｎｇ：ＭＲＩ）装置など）であれば、いかなる装置であってもよい。説明を具体的にするために、本実施形態に係る医用診断装置は、超音波診断装置として説明する。なお、以下の説明において、略同一の構成を有する構成要素については、同一符号を付し、重複説明は必要な場合にのみ行う。

図１は、本実施形態に係る医用診断装置としての超音波診断装置１の構成を示す構成図である。同図に示すように、超音波診断装置１は、超音波プローブ３、入力装置５、モニタ７、装置本体９を有する。加えて、本超音波診断装置１には、心電計、心音計、脈波計、呼吸センサに代表される図示していない生体信号計測器、不図示の外部記憶装置およびネットワークが、インターフェース回路２７を介して接続されてもよい。

超音波プローブ３は、複数の圧電振動子と、整合層と、複数の圧電振動子の背面側に設けられるバッキング材とを有する。複数の圧電振動子は、圧電セラミックス等の音響／電気可逆的変換素子である。複数の圧電振動子は並列され、超音波プローブ３の先端に装備される。以下、一つの圧電振動子が一チャンネルを構成するものとして説明する。圧電振動子は、後述する超音波送信回路１１から供給される駆動信号に応答して超音波を発生する。超音波プローブ３を介して被検体Ｐに超音波が送信されると、送信された超音波（以下、送信超音波と呼ぶ）は、被検体内の生体組織における音響インピーダンスの不連続面で反射される。

圧電振動子は、反射された超音波を受信し、エコー信号を発生する。エコー信号の振幅は、超音波の反射に関する不連続面を境界とする音響インピーダンスの差に依存する。また、送信超音波が移動している血流、および心臓壁等の表面で反射された場合のエコー信号の周波数は、ドプラ効果により、移動体（血流および心臓壁の表面）の超音波送信方向の速度成分に依存して偏移する。

以下、超音波プローブ３は、１次元的に配列された圧電振動子により構成される１次元アレイにより、被走査領域を２次元的に走査するプローブとして説明する。なお、超音波プローブ３は、１次元アレイを複数の圧電振動子の配列方向と直交する方向に揺動させて３次元走査を実行するメカニカル４次元プローブであってもよい。また、超音波プローブ３は、メカニカル４次元プローブに限定されず、２次元アレイプローブであってもよい。

整合層は、被検体Ｐに対する超音波の送受信を効率よくするために、複数の圧電振動子の超音波放射面側に設けられる。バッキング材は、圧電振動子の後方への超音波の伝搬を防止する。

入力装置５は、インターフェース回路２７を介して装置本体９に接続される。入力装置は、操作者からの各種指示、各種条件、関心領域（ＲＯＩ）の設定指示、種々の画質条件、および設定指示等を装置本体９にとりこむための各種スイッチ、ボタン、トラックボールの他、マウス、キーボード等を有する。なお、入力装置５は、操作面へ触れることで入力操作を行うタッチパッド、および表示画面とタッチパッドとが一体化されたタッチパネルディスプレイ、マイク等を有していてもよい。入力装置５は、入力部や入力インターフェース回路に相当する。

入力装置５は、後述する設定機能２９１、および後述する推定機能２９３を包括的に実行させる機能（以下、構造物推定機能と呼ぶ）に続いて、後述するゲイン逆補正機能２９５、後述する音場特性補正機能２９７、および後述する解析機能２９９等を包括的に実行させる機能（以下、組織性状解析機能と呼ぶ）を実行するための開始指示（以下、減衰定量開始指示と呼ぶ）を入力する。このとき、減衰定量開始指示の入力に関する信号は、後述する制御回路２９に出力される。

なお、入力装置５は、マウス、キーボードなどの物理的な操作部品を備えるものだけに限らない。例えば、超音波診断装置１とは別体に設けられた外部の入力機器から入力操作に対応する電気信号を受け取り、受け取った電気信号を種々の回路へ出力するような電気信号の処理回路も入力装置５の例に含まれる。

モニタ７は、後述する画像生成回路１９、画像合成回路等２３から出力されたビデオ信号に基づいて、生体内の形態学的情報や、血流情報などを画像として表示する。また、モニタ７は、解析機能２９９による解析結果を表示する。モニタとしては、例えば、ＣＲＴディスプレイや液晶ディスプレイ、有機ＥＬディスプレイ、ＬＥＤディスプレイ、プラズマディスプレイ、又は当技術分野で知られている他の任意のディスプレイが適宜利用可能である。モニタ７は、表示部や表示回路に相当する。

装置本体９は、超音波送信回路１１と、超音波受信回路１３と、Ｂモード処理回路１５と、ドプラ処理回路１７と、画像生成回路１９と、画像メモリ２１と、画像合成回路２３と、記憶回路２５と、インターフェース回路２７と、制御回路（中央演算処理装置：Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）２９とを有する。

超音波送信回路１１は、パルス発生器１１１と、送信遅延回路１１３と、パルサ回路１１５とを有する。超音波送信回路１１は、超音波送信部の一例であって、プロセッサを有していてもよい。パルス発生器１１１は、所定のレート周波数ｆｒ Ｈｚ（周期：１／ｆｒ秒）で、送信超音波を形成するためのレートパルスを繰り返し発生する。発生されたレートパルスは、チャンネル数に分配され、送信遅延回路１１３に送られる。

送信遅延回路１１３は、複数のチャンネルごとに、送信超音波をビーム状に収束し、かつ送信指向性を決定するために必要な遅延時間（以下、送信遅延時間と呼ぶ）を、各レートパルスに与える。送信超音波の送信方向または送信方向に関する送信遅延時間（以下、送信遅延パターンと呼ぶ）は、記憶回路２５に記憶される。記憶回路２５に記憶された送信遅延パターンは、制御回路２９により超音波の送信時に参照される。

パルサ回路１１５は、このレートパルスに基づくタイミングで、超音波プローブ３の振動子ごとに電圧パルス（駆動信号）を印加する。これにより、超音波ビームが被検体Ｐに送信される。

超音波送信回路１１は、１フレームに対応する後述のＢモードデータの生成を契機として、減衰定量条件に基づいて、被検体Ｐに超音波（以下、減衰定量用超音波と呼ぶ）を送信する。減衰定量条件とは、減衰定量用超音波を発生させるための送信条件、減衰定量用超音波を受信するための受信条件である。送信条件は、例えば、減衰定量用超音波の送信中心周波数と減衰定量用超音波の送信帯域幅である。受信条件は、例えば、減衰定量用超音波の受信中心周波数と受信帯域幅である。減衰定量条件は、記憶回路２５に記憶される。

なお、送信条件は、上記２種類に限定されず、例えば、周波数が異なる複数の減衰定量用超音波を、一つの走査線に対して送信する条件であってもよい。また、送信条件は、例えば、位相を反転させた２つの減衰定量用超音波を、一つの走査線に対して送信する条件であってもよい。送信条件における送信中心周波数は、例えば、Ｂモード用超音波の送信中心周波数より高い送信中心周波数である。また、送信条件における帯域幅は、Ｂモード用超音波の帯域幅より狭い帯域幅（以下、狭帯域と呼ぶ）である。

入力装置５を介した操作者の減衰定量開始指示により、減衰定量条件は、記憶回路２５から制御回路２９に読み出される。超音波送信回路１１は、減衰定量条件における送信条件にしたがって、制御回路２９により制御される。例えば、超音波送信回路１１は、超音波スキャンにおいて、送信条件にしたがって、周波数が異なる複数の超音波を、超音波プローブ３を介して被検体Ｐに送信する。なお、超音波送信回路１１は、超音波スキャンにおいて、送信条件にしたがって、Ｂモードに関する超音波送信における周波数帯域より狭帯域の超音波を、超音波プローブ３を介して被検体Ｐに送信する。

超音波受信回路１３は、プリアンプ１３１、図示していないアナログディジタル（ａｎａｌｏｇ ｔｏ ｄｉｇｉｔａｌ（以下、Ａ／Ｄと呼ぶ））変換器、受信遅延回路１３３、加算器１３５を有する。超音波受信回路は、超音波受信部の一例であって、プロセッサを有していてもよい。プリアンプ１３１は、超音波プローブ３を介して取り込まれた被検体Ｐからのエコー信号をチャンネル毎に増幅する。Ａ／Ｄ変換器は、増幅された受信エコー信号をディジタル信号に変換する。Ａ／Ｄ変換される前のアナログ信号に対して、アナログゲインが、ＳＴＣ（ｓｅｎｓｉｔｉｖｅ ｔｉｍｅ ｃｏｎｔｒｏｌ）またはＴＧＣ（ｔｉｍｅ ｇａｉｎ ｃｏｎｔｒｏｌ）として与えられる。

受信遅延回路１３３は、ディジタル信号に変換された受信エコー信号に、受信指向性を決定するために必要な遅延時間（以下、受信遅延時間と呼ぶ）を与える。受信遅延回路１３３は、例えば、ディジタルビームフォーマである。受信遅延回路１３３から出力されたディジタル信号に対して、ディジタルゲインが、ＳＴＣまたはＴＧＣとして与えられる。エコー信号の受信方向または受信方向に関する受信遅延時間（以下、受信遅延パターンと呼ぶ）は、後述する記憶回路２５に記憶される。記憶回路２５に記憶された受信遅延パターンは、送信時と同様に、制御回路２９により参照される。

被検体内における超音波の減衰によって、被検体内の深部ほど反射波による信号は、微弱になる。このため、アナログゲインおよびディジタルゲインは、この減衰を補うために、被検体内の深部において反射された超音波に起因する信号ほど、この信号の振幅を増幅するゲインである。

加算器１３５は、遅延時間が与えられた複数のエコー信号を加算する。この加算により、超音波受信回路１３は、受信指向性に応じた方向からの反射成分を強調した受信信号を発生する。この送信指向性と受信指向性とにより超音波送受信の総合的な指向性が決定される。この総合的な指向性により、超音波ビーム（いわゆる「超音波走査線」）が決まる。

超音波受信回路１３は、減衰定量条件における受信条件にしたがって、減衰定量用超音波を受信する。受信条件における受信中心周波数は、減衰定量用超音波の送信中心周波数と略同一の周波数であって、被走査領域の深さ方向に対して変化させずに一定である。また、受信条件における受信帯域幅は、狭帯域と略同一な帯域幅である。超音波受信回路１３は、減衰定量条件における受信条件にしたがって、制御回路２９により制御される。

具体的には、超音波受信回路１３は、１フレームに対応する後述のＢモードデータの生成を契機として、被検体Ｐに送信された減衰定量用超音波の反射波を、受信条件にしたがって受信する。超音波受信回路１３は、減衰定量用超音波の反射波の受信により、減衰定量用受信データを生成する。超音波送信回路１１は、減衰定量用受信データを、Ｂモード処理回路１５に出力する。なお、超音波送信回路１１は、減衰定量用受信データを、制御回路２９および記憶回路２５に出力してもよい。

Ｂモード処理回路１５は、図示していない包絡線検波器、対数変換器などを有する。Ｂモード処理回路１５は、Ｂモード処理部の一例であって、プロセッサを有する。包絡線検波器は、超音波受信回路１３から出力された受信信号に対して包絡線検波を実行する。包絡線検波器は、包絡線検波された信号を、後述する対数変換器に出力する。対数変換器は、包絡線検波された信号に対して対数変換して弱い信号を相対的に強調する。Ｂモード処理回路１５は、対数変換器により強調された信号に基づいて、各走査線および各超音波送受信における深さごとの信号値（Ｂモードデータ）を生成する。

Ｂモードデータは、対数変換器から出力された信号の強度を輝度の明るさとして表現されたデータに相当する。Ｂモード処理回路１５からの出力は、画像生成回路１９に出力される。Ｂモード処理回路１５からの出力は、反射波の強度を輝度で表したＢモード画像として、モニタ７に表示される。Ｂモード処理回路１５は、減衰定量用受信データに基づいて、上記処理手順により、減衰定量用Ｂモードデータを生成する。Ｂモード処理回路１５は、減衰定量用Ｂモードデータを、制御回路２９および記憶回路２５に出力する。

超音波プローブ３がメカニカル４次元プローブである場合や２次元アレイプローブである場合、Ｂモード処理回路１５は、被走査領域におけるアジマス（Ａｚｉｍｕｔｈ）方向、エレベーション（Ｅｌｅｖａｔｉｏｎ）方向、深さ方向（レンジ（Ｒａｎｇｅ）方向）にそれぞれ対応付けて配列された複数の信号値からなる３次元Ｂモードデータを発生してもよい。レンジ方向とは、走査線上の深さ方向である。アジマス方向とは、例えば、１次元アレイにおける圧電振動子の配列方向に沿った電子走査方向である。エレベーション方向とは、例えば、１次元アレイの機械的揺動方向である。

なお、３次元Ｂモードデータは、複数の画素値または複数の輝度値などを、走査線に沿って、アジマス方向、エレベーション方向、レンジ方向にそれぞれ対応付けて配列させたデータであってもよい。また、３次元Ｂモードデータは、被走査領域において予め設定されたＲＯＩに関するデータであってもよい。また、Ｂモード処理回路１５は、３次元Ｂモードデータの代わりにボリュームデータを発生してもよい。以下、Ｂモード処理回路１５で発生されるデータをまとめて、Ｂモードデータと呼ぶ。

ドプラ処理回路１７は、超音波受信回路１３から受け取ったエコー信号に対して速度情報を周波数解析する。ドプラ処理回路１７は、ドプラ処理部の一例であって、プロセッサを有する。ドプラ処理回路１７は、ドプラ効果による血流や組織、造影剤エコー成分を、超音波受信回路１３から受け取ったエコー信号から抽出する。ドプラ処理回路１７は、平均速度、分散、パワー等の血流情報を、走査線上の多点について求める。ドプラ処理回路１７は、得られた血流情報を画像生成回路１９に出力する。ドプラ処理回路１７からの出力は、ドプラ波形画像、平均速度画像、分散画像、パワー画像、およびこれらの組み合わせ画像として、モニタ７にカラー表示される。

例えば、ドプラ処理回路１７は、図示していないミキサー、低域通過フィルタ（Ｌｏｗ Ｐａｓｓ Ｆｉｌｔｅｒ：以下、ＬＰＦと呼ぶ）、速度／分散／Ｐｏｗｅｒ演算回路等を有する。ミキサーは、超音波受信回路１３から出力された受信信号に、送信周波数と同じ周波数ｆ_０を有する基準信号を掛け合わせる。この掛け合わせにより、ドプラ偏移周波数ｆ_ｄの成分の信号と（２ｆ_０＋ｆ_ｄ）の周波数成分を有する信号とが得られる。ＬＰＦは、ミキサーからの２種の周波数成分を有する信号のうち、高い周波数成分（２ｆ_０＋ｆ_ｄ）の信号を取り除く。ドプラ処理回路１７は、高い周波数成分（２ｆ_０＋ｆ_ｄ）の信号を取り除くことにより、ドプラ偏移周波数ｆ_ｄの成分を有するドプラ信号を発生する。

なお、ドプラ処理回路１７は、ドプラ信号を発生するために、直交検波方式を用いてもよい。このとき、受信信号（ＲＦ信号）は、直交検波されＩＱ信号に変換される。ドプラ処理ユニット１４２は、ＩＱ信号を複素フーリエ変換することにより、ドプラ偏移周波数ｆ_ｄの成分を有するドプラ信号を発生する。ドプラ信号は、例えば、血流、組織、造影剤によるドプラ成分である。

速度／分散／Ｐｏｗｅｒ演算回路は、図示していないＭＴＩ（Ｍｏｖｉｎｇ Ｔａｒｇｅｔ Ｉｎｄｉｃａｔｏｒ）フィルタ、ＬＰＦフィルタ、自己相関演算器等を有する。なお、自己相関演算器の代わりに相互相関演算器を有していてもよい。ＭＴＩフィルタは、発生されたドプラ信号に対して、臓器の呼吸性移動や拍動性移動などに起因するドプラ成分（クラッタ成分）を除去する。ＭＴＩフィルタは、ドプラ信号から血流に関するドプラ成分（以下、血流ドプラ成分と呼ぶ）を抽出するために用いられる。ＬＰＦは、ドプラ信号から組織の移動に関するドプラ成分（以下、組織ドプラ成分と呼ぶ）を抽出するために用いられる。

自己相関演算器は、血流ドプラ成分及び組織ドプラ成分に対して自己相関値を算出する。自己相関演算器は、算出された自己相関値に基づいて、血流および組織の平均速度値、分散値、ドプラ信号の反射強度（パワー）等を算出する。速度／分散／Ｐｏｗｅｒ演算回路は、複数のドプラ信号に基づく血流および組織の平均速度値、分散値、ドプラ信号の反射強度等に基づいて、所定領域の各位置におけるカラードプラデータを発生する。以下、ドプラ信号とカラードプラデータとをまとめて、ドプラデータと呼ぶ。

画像生成回路１９は、図示していないディジタルスキャンコンバータ（Ｄｉｇｉｔａｌ Ｓｃａｎ Ｃｏｎｖｅｒｔｅｒ：以下、ＤＳＣと呼ぶ）等を有する。画像生成回路１９は、画像生成部の一例であって、プロセッサを有する。画像生成回路１９は、ＤＳＣに対して、座標変換処理（リサンプリング）を実行する。座標変換処理とは、例えば、Ｂモードデータ、およびドプラデータからなる超音波スキャンの走査線の信号列を、テレビなどに代表される一般的なビデオフォーマットの走査線信号列に変換する処理である。

画像生成回路１９は、座標変換処理により、表示画像としての超音波画像を生成する。具体的には、画像生成回路１９は、Ｂモードデータに基づいてＢモード画像を生成する。画像生成回路１９は、減衰定量用Ｂモードデータに基づいて、減衰Ｂモード画像を生成する。Ｂモード画像および減衰Ｂモード画像は、音波の集束などの超音波プローブの特性や超音波ビーム（例えば、送受信ビーム）の音場特性などが反映された画素値（輝度値）を有する。例えば、Ｂモード画像において、被走査領域において超音波のフォーカス付近では、非フォーカス部分よりも相対的に高輝度となる。

画像生成回路１９は、ドプラデータに基づいて、平均速度画像、分散画像、パワー画像などのドプラ画像を生成する。また、画像生成回路１９は、解析機能２９９により解析された解析結果に基づいて、被走査領域における一部領域の各位置における超音波の減衰量を示す減衰定量画像を生成する。

画像メモリ２１は、発生された超音波画像（Ｂモード画像、平均速度画像、分散画像、パワー画像、減衰定量画像）に対応するデータ（以下、画像データと呼ぶ）を記憶する。画像メモリ２１に記憶された画像データは、入力装置５を介した操作者の指示により、読み出される。画像メモリ２１は、例えば、フリーズする直前の複数のフレームに対応する超音波画像を保存するメモリである。このシネメモリに記憶されている画像を所定のフレームレートで連続表示（シネ表示）することで、超音波動画像が、モニタ７に動画表示される。

画像メモリ２１は、例えば、集積回路記憶装置（ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）、ＲＯＭ（Ｒｅａｄ−Ｏｎｌｙ Ｍｅｍｏｒｙ）等）により実現される。なお、画像メモリ２１の実現は、上記集積回路記憶装置に限定されず、任意の記憶装置であってもよい。

画像合成回路２３は、超音波画像に、種々のパラメータの文字情報および目盛等を合成する。画像合成回路２３は、画像合成部の一例であって、プロセッサを有する。画像合成回路２３は、合成された超音波画像を後述するモニタ７に出力する。画像合成回路２３は、Ｂモード画像に減衰定量画像を位置合わせして重畳した減衰重畳画像を生成する。画像合成回路２３は、発生した減衰重畳画像を、モニタ７に出力する。

記憶回路２５は、種々の情報を記憶するＨＤＤ（ｈａｒｄ ｄｉｓｋ ｄｒｉｖｅ：ハードディスクドライブ）やＳＳＤ（ｓｏｌｉｄ ｓｔａｔｅ ｄｒｉｖｅ：ソリッドステートドライブ）、集積回路記憶装置（ＲＡＭ、ＲＯＭ等）などの記憶装置である。記憶回路２５は、記憶部に相当する。また、記憶回路２５は、ＣＤ−ＲＯＭドライブやＤＶＤドライブ等との間で、種々の情報を読み書きする駆動装置で実現されてもよい。また、記憶回路２５は、磁気ディスク（フロッピー（登録商標）ディスクなど）、光ディスク（ＣＤ−ＲＯＭ、ＤＶＤ、ＭＯなど）、半導体メモリ等の可搬性記憶媒体との間で種々の情報を読み書きする駆動装置で実現されてもよい。

記憶回路２５は、フォーカス深度の異なる複数の受信遅延パターン、および複数の送信遅延パターンを記憶する。記憶回路２５は、本超音波診断装置１の制御プログラム、診断プロトコル、および後述する医用解析プログラムを記憶する。記憶回路２５は、超音波の送受信条件等の各種データ群、診断情報（患者ＩＤ、医師の所見等）を記憶する。記憶回路２５は、超音波受信回路１３により発生された受信信号、Ｂモード処理回路１５により発生されたＢモードデータ、ドプラ処理回路１７により発生されたドプラデータ、および解析機能２９９による解析結果を示す解析データを記憶する。

また、記憶回路２５は、減衰定量用Ｂモードデータの収集に関する被走査領域（減衰Ｂモード画像）に設定される複数の領域の大きさおよび設定位置、構造物推定機能の実行に関するプログラム（以下、構造物推定プログラムと呼ぶ）等を記憶する。記憶回路２５は、構造物推定機能において参照される閾値（以下、構造物判定閾値と呼ぶ）を記憶する。記憶回路２５は、ゲイン逆補正機能２９５で用いられる逆補正データ、音場特性補正機能２９７で用いられる音場特性補正データ、組織性状解析機能の実行に関するプログラム（以下、組織性状解析プログラムと呼ぶ）等を記憶する。以下、構造物推定プログラムと組織性状解析プログラムとをまとめて、医用解析プログラムと呼ぶ。

なお、記憶回路２５は、逆補正データの代わりに、対応表（以下、逆補正対応表と呼ぶ）を記憶してもよい。また、記憶回路２５は、音場特性補正データの代わりに、対応表（以下、音場特性補正対応表と呼ぶ）を記憶してもよい。

逆補正データは、超音波スキャンにより得られたデータ（受信信号、受信データ）に対してＳＴＣまたはＴＧＣにより付与されたアナログゲインおよびディジタルゲインをキャンセルするためのデータである。具体的には、逆補正データは、超音波スキャンにより得られたデータに付与されるゲイン（利得）の深さ方向の応答を示すデータである。すなわち、逆補正データをＢモードデータに適用することにより、ゲイン補正されたＢモードデータは、ゲイン補正前のＢモードデータに変換される。

逆補正対応表とは、ＳＴＣまたはＴＧＣによるゲイン補正をキャンセルするための対応表である。具体的には、逆補正対応表は、超音波スキャンにより得られたデータに付与されるゲイン（利得）の深さ方向の応答を示す対応表であって、ゲイン補正後のＢモードデータを、ゲイン補正前のＢモードデータに変換するための対応表である。

音場特性補正データは、ゲイン補正前のＢモードデータにおける音場特性の依存性をキャンセルするためのデータである。すなわち、音場特性補正データをＢモードデータに適用することにより、音場特性に依存するＢモードデータは、音場特性に依存しないＢモードデータに変換される。音場特性補正対応表とは、Ｂモードデータにおける音場特性の依存性をキャンセルするための対応表である。具体的には、音場特性補正対応表は、音場特性に依存するＢモードデータを、音場特性に依存しないＢモードデータに変換するための対応表である。

音場特性補正データおよび音場特性補正対応表は、例えば、超音波が無減衰であってかつ一様な散乱体を有する物体に対して超音波スキャンを実行したときに得られる深さ方向の画素値（または輝度値）の分布に対応する。音場特性補正データおよび音場特性補正対応表は、超音波が無減衰であってかつ一様なファントムに対して超音波スキャンを実行することにより取得された実測データに基づいて予め生成される。

なお、音場特性補正データおよび音場特性補正対応表は、超音波に対して一定の減衰を有するファントムに対する超音波スキャンにより実測された実測データからこのファントムによる超音波の減衰分を差分することにより生成されてもよい。また、音場特性補正データおよび音場特性補正対応表は、シミュレーションなどの他の手段により生成されてもよい。

記憶回路２５は、Ｂモード画像、平均速度画像、分散画像、パワー画像、減衰定量画像、減衰重畳画像等の各種画像を記憶する。記憶回路２５は、後述する複数の減衰定数にそれぞれ対応する複数の色相を記憶する。記憶回路２５は、減衰定量画像に関する所定の不透明度または透明度を記憶する。記憶回路２５は、後述する分散比の値に応じた色相を記憶する。なお、上述した画像メモリ２１は、記憶回路２５に設けられてもよい。また、構造物推定機能としてＣＦＡＲ（Ｃｏｎｔｒａｓｔ Ｆａｌｓｅ Ａｌａｒｍ Ｒａｔｅ）処理を実行する場合、記憶回路２５は、ＣＦＡＲ処理に関するプログラムを記憶してもよい。

インターフェース回路２７は、入力装置５、不図示の操作パネル、ネットワーク、図示していない外部記憶装置および生体信号計測器に関するインターフェースである。装置本体９によって得られた超音波画像等のデータおよび解析結果等は、インターフェース回路２７とネットワークとを介して他の装置に転送可能である。なお、インターフェース回路２７は、ネットワークを介して、図示していない他の医用画像診断装置で取得された被検体に関する医用画像を、ダウンロードすることも可能である。インターフェース回路２７は、インターフェース部に対応し、プロセッサを有していてもよい。

制御回路２９は、情報処理装置（計算機）としての機能を有し、本超音波診断装置１の装置本体９の動作を制御する制御手段（プロセッサ）である。制御回路２９は、記憶回路２５から画像生成・表示等を実行するための制御プログラムを読み出して各種処理に関する演算・制御等を実行する。制御回路２９は、制御部に相当する。

制御回路２９は、減衰定量開始指示に応答して、記憶回路２５から減衰定量条件を読み出す。制御回路２９は、読み出した減衰定量条件にしたがって、超音波送信回路１１と超音波受信回路１３とを制御する。具体的には、制御回路２９は、読み出した送信条件にしたがって超音波送信回路１１を制御する。これにより、超音波送信回路１１は、１フレームに対応するＢモードデータの生成後に、減衰定量用超音波を、被検体Ｐに送信する。また、制御回路２９は、読み出した受信条件にしたがって超音波受信回路１３を制御する。これにより、超音波受信回路１３は、受信条件にしたがって、被検体Ｐに送信された減衰定量用超音波の反射波を受信する。

本実施形態では、設定機能２９１、推定機能２９３、ゲイン逆補正機能２９５、音場特性補正機能２９７、解析機能２９９にて行われる各処理機能は、コンピュータによって実行可能なプログラムの形態で記憶回路２５へ記憶されている。制御回路２９は、これら機能に対応するプログラムを記憶回路２５から読み出し、実行することで各プログラムに対応する機能を実現するプロセッサである。換言すると、各プログラムを読みだした状態の制御回路２９は、図１の制御回路２９内に示された各機能を有することになる。

なお、上述のそれぞれの機能がプログラムとして構成され、１つの制御回路が各プログラムを実行する場合であってもよいし、推定の機能が専用の独立したプログラム実行回路に実装される場合であってもよい。なお、制御回路２９が有する設定機能２９１、推定機能２９３、ゲイン逆補正機能２９５、音場特性補正機能２９７、解析機能２９９は、それぞれ、設定部、推定部、ゲイン逆補正部、音場特性補正部、解析部の一例である。また、このとき、構造物推定機能を実現するプロセッサは、構造物推定部として機能してもよい。また、組織性状解析機能を実現するプロセッサは、組織性状解析部として機能してもよい。

上記説明において用いた「プロセッサ」という文言は、例えば、ＣＰＵ、ＧＰＵ（Ｇｒａｐｈｉｃａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）或いは、特定用途向け集積回路（Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ：ＡＳＩＣ）、プログラマブル論理デバイス（例えば、単純プログラマブル論理デバイス（Ｓｉｍｐｌｅ Ｐｒｏｇｒａｍｍａｂｌｅ Ｌｏｇｉｃ Ｄｅｖｉｃｅ：ＳＰＬＤ）、複合プログラマブル論理デバイス（Ｃｏｍｐｌｅｘ Ｐｒｏｇｒａｍｍａｂｌｅ Ｌｏｇｉｃ Ｄｅｖｉｃｅ：ＣＰＬＤ）、及びフィールドプログラマブルゲートアレイ（Ｆｉｅｌｄ Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ：ＦＰＧＡ））等の回路を意味する。

プロセッサは、記憶回路１２３に保存されたプログラムを読み出し実行することで機能を実現する。なお、記憶回路１２３にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むよう構成しても構わない。この場合、プロセッサは回路内に組み込まれたプログラムを読み出し実行することで機能を実現する。なお、超音波送信回路１１、超音波受信回路１３、Ｂモード処理回路１５、ドプラ処理回路１７、画像生成回路１９、画像合成回路２３、インターフェース回路２７等の他の回路も同様に、上記のプロセッサなどの電子回路により構成される。

制御回路２９は、減衰定量開始指示に応答して、医用解析プログラムを記憶回路２５から読み出す。制御回路２９は、読み出した医用解析プログラムを実行することにより、構造物推定機能および組織性状解析機能を実現する。具体的には、制御回路２９は、構造物推定プログラムおよび組織性状解析プログラムを、記憶回路２５から読み出す。制御回路２９は、読み出した構造物推定プログラムを実行することにより、減衰定量用超音波を用いた超音波スキャンにより得られた減衰定量用Ｂモードデータに基づいて、被検体内の構造物の位置を推定する。

制御回路２９は、例えば、構造物の位置の推定を契機として、読み出した組織性状解析プログラムを実行することにより、減衰定量用Ｂモードデータに対応する被走査領域において、組織性状を解析する。以下、被検体内の構造物の位置の推定に関する各種機能、および組織性状の解析に関する各種機能について、具体的に説明する。

以下の説明においては、単一の制御回路２９において、設定機能２９１、推定機能２９３、ゲイン逆補正機能２９５、音場特性補正機能２９７、および解析機能２９９等が実行されるものとして説明しているが、複数の独立したプロセッサを組み合わせて処理回路を構成し、各プロセッサがプログラムを実行することにより各種機能を実現するものとしても構わない。また、設定機能２９１、推定機能２９３、ゲイン逆補正機能２９５、音場特性補正機能２９７、および解析機能２９９等は、それぞれ異なる処理回路またはそれぞれ異なるプロセッサで実現されてもよい。

（構造物推定機能）
構造物推定機能とは、例えば、被検体Ｐに対する超音波スキャンにより得られた減衰定量用Ｂモードデータに基づいて、被検体内の複数の位置について不均一さを評価することにより構造物の位置を推定する機能を含む。構造物推定機能として構造物を推定する方法には様々な方法があるが、以下、一例として減衰Ｂモード画像に設定された複数の領域各々に含まれる複数の画素にそれぞれ対応する複数の画素値（または輝度値）の平均値、および分散値を利用した方法を説明する。

この方法は、一般に、例えば肝臓の実質のような一様な散乱体を有する領域において、減衰定量用Ｂモードデータの度数分布（ヒストグラム）がレイリー分布を示すことを利用している。一様な散乱体を有する領域において、減衰定量用Ｂモードデータの度数分布がレイリー分布に従う場合、この領域における複数の画素値による平均値（μ）と分散値（σ^２）とは、以下のような関係となる。

σ^２＝（４／π−１）×μ^２ ・・・（１）
平均値と分散値との計算に関する領域に構造物が含まれている場合、散乱体は非一様であるため、減衰定量用Ｂモードデータの度数分布の分散値は、レイリー分布で計算される分散値より大きくなる。加えて、構造物の特性が周囲と顕著に異なるほど（例えば、構造物の領域が石灰化している領域に対応する場合など）、分散値は大きくなる。このため、被走査領域で設定された複数の領域各々において分差値を算出することにより、領域を代表する画素各々における構造物の有無を判定することができる。

しかしながら、画素値の大きさが大きいほど分散値が大きくなるため、分散値そのものを算出するだけでは、構造物の有無、すなわち一様な散乱体の場合に得られるレイリー分布からの逸脱の有無を判定することはできない。そこで、本実施形態では、設定された領域における画素値の平均値と同じ平均値を有するレイリー分布を、設定された領域に対して仮定する。この仮定のもとで、設定された領域における画素値の平均値と式（１）とを用いて計算された分散値｛（４／π−１）×μ^２｝に対する、設定された領域における画素値の分散値σ^２の比（以下、分散比Ｒ_σと呼ぶ）が計算される。具体的には、分散比Ｒ_σは、以下のような式となる。
Ｒ_σ＝σ^２／｛（４／π−１）×μ^２｝ ・・・（２）
式（２）の右辺における分子（σ^２）は、設定された領域に包含される複数の画素にそれぞれ対応する複数の画素値から計算された実測の分散値である。式（２）の右辺における分母｛（４／π−１）×μ^２｝は、設定された領域に包含される複数の画素にそれぞれ対応する分散の画素値から計算された平均値μと式（１）とを用いて、設定された領域における複数の画素値がレイリー分布を形成すると仮定した場合の分散値である。

分散比Ｒ_σが１に近い場合、設定された領域における画素値の分布は、レイリー分布とみなすことができる。また、分散比Ｒ_σが１より大きな値である場合、設定された領域に含まれる複数の画像地の分布は、レイリー分布から乖離しており、設定された領域には、一様な散乱体から逸脱した構造物が含まれていると推定される。すなわち、分散比Ｒ_σは、設定された領域における構造物の有無の判断の指標に相当する。

設定機能２９１を実現する制御回路２９は、減衰定量用Ｂモードデータが収集された被走査領域（減衰Ｂモード画像）、すなわち減衰定量用超音波による超音波スキャンにおける走査領域において、複数の領域を設定する。設定された複数の領域は、減衰Ｂモード画像における複数の画素各々を中心（以下、中心画素と呼ぶ）または重心として、所定の大きさを有する。減衰Ｂモード画像に設定される領域の大きさは、入力装置５を介した操作者の指示により適宜変更されてもよい。なお、制御回路２９は、所定の画素の幅ごとに一つの領域をスイープさせることで、減衰Ｂモード画像において複数の領域を設定してもよい。

推定機能２９３を実現する制御回路２９は、設定された複数の領域各々に含まれる複数の画素各々に対応する画素値（または輝度値）に基づいて、構造物の位置を推定する。このとき、画素値は、減衰定量用Ｂモードデータにおける画素値に対応する。なお、推定機能２９３により用いられる画素値は、通常のＢモードデータにおける画素値が用いられてもよい。具体的には、制御回路２９は、複数の領域各々において、画素値（または輝度値）の平均値と、分散値とを計算する。制御回路２９は、平均値と分散値とに基づいて、分散比Ｒ_σを計算する。計算された分散比Ｒ_σは、記憶回路２５に記憶されてもよい。また、分散比Ｒ_σは、分散比Ｒ_σの値に応じた色相で、分散比Ｒ_σに関する被走査領域と同一被走査領域におけるＢモード画像に重畳されて、モニタ７に表示されてもよい。

図２は、分散比Ｒ_σの値に応じた色相で分散比Ｒ_σが重畳されるＢモード画像である。図２に示すＢモード画像は、被検体Ｐの肝臓の断面を示している。図３は、分散比Ｒ_σの計算に関する同一被走査領域における図２のＢモード画像に、分散比Ｒ_σを、分散比Ｒ_σの値に応じた色相で重畳した分散比重畳画像の一例を示す図である。図３におけるハッチングの違いは、色相の違いに相当する。なお、図３においては、分散比Ｒ_σの色相は、分散比重畳画像を図示しやすいように３種類としているが、実際には、略連続的な色相として、分散比Ｒ_σの色相を示す凡例とともにモニタ７に表示される。

図３におけるハッチングＬＲは、分散比Ｒ_σが大きい領域を示している。図３に示すように、例えば、血管、腹壁、横隔膜などの構造物および構造物の周辺の領域では、均一なスペックルを有している肝臓の実質と比べて分散比Ｒ_σが大きくなっている。すなわち、図３から明らかなように、構造物および構造物の周辺の領域は、超音波の散乱特性が不均一である。このため、構造物および構造物の周辺の領域は、超音波の減衰量の解析対象としては不適である。

制御回路２９は、計算した分散比Ｒ_σを、設定機能２９１により設定された領域を代表する位置、例えば中心画素の位置に対応付ける。制御回路２９は、構造物判定閾値を、記憶回路２５から読み出す。構造物判定閾値は、１以上の数値である。なお、構造物判定閾値は、入力装置５を介した操作者の指示により適宜変更されてもよい。制御回路２９は、読み出した構造物判定閾値と分散比Ｒ_σとを比較する。制御回路２９は、構造物判定閾値より大きい分散比Ｒ_σを、構造物の位置として推定する。

推定機能２９３を実現する制御回路２９は、減衰Ｂモード画像に対応する被走査領域において推定された構造物の位置、すなわち構造物の領域を、記憶回路２５に出力してもよい。推定された構造物の位置は、後述する組織性状解析機能において用いられる。推定された構造物の領域は、例えば、図３に示すように、例えば、ハッチングＬＲで示す領域に対応する。

なお、受信信号や減衰定量用Ｂモードデータを用いた構造物の推定の手法は、上述した方法に限定されない。例えば、構造物推定機能は、ＣＦＡＲ処理と称される推定的な信号（または画像）の抽出技術を用いてもよい。ＣＦＡＲ処理という用語は、レーダー分野に用いられるものである。本実施形態では、その関連性により説明を具体的にするために、便宜上「ＣＦＡＲ」という語句を用いている。しかしながら、レーダー分野で用いられる方法、あるいは統計量を厳密に使用したものに拘泥されない。

ＣＦＡＲ処理は、例えば次の（１） 〜（３）の手順によって実行される。
（１）設定機能２９１を実現する制御回路２９は、減衰Ｂモード画像において、注目画素Ｐｉ毎に、当該画素Ｐｉの近傍画素を有する領域を設定する。制御回路２９により近傍画素として設定される領域は、減衰Ｂモード画像において十字型に設けられる。しかしながら、設定される領域における近傍画素の配列は、十字型に拘泥されず、例えば演算処理に要する時間が問題とならない場合には、注目画素に隣接する８画素を除く任意の大きさの領域であってもよい。

推定機能２９３を実現する制御回路２９は、設定された領域における輝度平均値（または画素平均値）を計算する。このとき、制御回路２９は、注目画素の輝度値（または画素値）が平均値に影響しないようにするため、注目画素Ｐｉ自体を、輝度平均計算に含めないようにしてもよい。

（２）次に、制御回路２９は、注目画素Ｐｉの画素値から平均値を減算する。制御回路２９は、この減算値を、当該注目画素Ｐｉの位置に対する演算結果Ｋｉとして定義し、記憶回路２５に記憶させる。制御回路２９は、この演算処理を、全ての注目画素Ｐｉについて実行する。

（３）制御回路２９は、予め設定された閾値Ｔを記憶回路２５から読み出す。このとき、閾値Ｔは、構造物判定閾値に対応し、一般的には分散比Ｒ_σに対応する構造物判定閾値とは異なる値である。制御回路２９は、演算結果Ｋｉと閾値Ｔとを比較する。Ｋｉ≧Ｔの場合には，元の輝度を用いて当該注目画素Ｐｉを表示する（構造物の抽出）。一方、Ｋｉ＜Ｔである場合、当該注目画素Ｐｉの輝度値はゼロとすることで表示しない（構造物の除去）。これらの処理を全ての注目画素Ｐｉについて実行することで、当該画像に関するＣＦＡＲ処理を実行することができる。

なお、上記（３）の判定においては、Ｋｉ≧Ｔの場合には輝度をＫｉとして当該注目画素Ｐｉを表示し、Ｋｉ＜Ｔである場合、当該注目画素Ｐｉの輝度値はゼロとすることで表示しない様にしてもよい。

また、構造物推定機能は、より単純には、通常のＢモード画像において設定された複数の領域における輝度値（または画素値）の平均値、および分散値、標準偏差値等と、これらの値各々に対応する構造物判定のための閾値とを比較することにより、構造物の位置を推定してもよい。

（組織性状解析機能）
組織性状解析機能は、組織性状解析プログラムにしたがって制御回路２９により実行される機能である。具体的には、組織性状解析機能は、減衰定量用Ｂモードデータに基づいて、減衰Ｂモード画像において推定された構造物の位置を除く複数の位置に対して組織性状を解析する。組織性状解析機能は、ゲイン逆補正機能２９５と音場特性補正機能２９７と、解析機能２９９と、不図示の表示機能とを有する。以下、組織性状解析機能における各機能について詳述する。組織性状とは、例えば、被検体内を伝搬する減衰定量用超音波の減衰の程度を示す減衰量に関する特徴量である。

なお、組織性状は、減衰量に限定されず、例えば、弾性率（ヤング率）、粘性率、歪みなどの診断対象組織の性状を示す量であってもよい。このとき、組織性状は、例えば、超音波エラストグラフィー法により取得される。解析機能２９９は、静的または動的なエラストグラフィー法に関する各種解析機能を有する。このとき、解析機能２９９を実現する制御回路２９は、組織性状として、被検体内の組織の粘性と弾性とのうち少なくとも一つに関する指標値（粘性パラメータ、弾性パラメータ）を計算する。

代表的な一つの超音波エラストグラフィー法とは、超音波プローブ３により被検体内の組織を体表から圧迫・解放することにより、圧迫・解放時に観測される断面内の各点の歪みの大きさから、被検体内の組織の相対的な硬さを可視化する方法である。

代表的な別の超音波エラストグラフィー法とは、音響放射力や機械的振動を被検体Ｐの体表から被検体内の組織に与えることにより、断面内の各点における組織の変位（歪み）を経時的に観測する方法である。具体的には、動的な超音波エラストグラフィー法は、音響放射力や機械的振動により発生させたせん断波（シア・ウェーブ（ｓｈｅａｒ ｗａｖｅ））の伝播速度を求めることにより、診断対象組織の弾性率や粘性率等を求める方法である。静的または動的な超音波エラストグラフィー法により送受信される超音波は、減衰定量用超音波に対応する。このとき、構造物推定機能において用いられるデータは、例えば、Ｂモード画像の生成に関するＢモードデータが用いられる。また、超音波エラストグラフィー法により生成される画像（弾性画像、粘性画像、歪み画像など）は、減衰定量画像に対応する。

以下、組織性状解析機能は、組織性状として減衰量を計算するものとして説明する。ゲイン逆補正機能２９５の実行に先立って、制御回路２９は、減衰定量用Ｂモードデータにおいて、減衰Ｂモード画像に対応する被走査領域に推定された構造物の位置を除く複数の位置に対応するデータ（以下、減衰定量用部分データと呼ぶ）を推定する。

ゲイン逆補正機能２９５を実現する制御回路２９は、記憶回路２５から逆補正データを読み出す。制御回路２９は、被走査領域において推定された構造物の位置を除く複数の位置に対応する減衰定量用部分データから逆補正データを減算する。これにより、制御回路２９は、ゲイン補正前の減衰定量用部分データを生成する。なお、制御回路２９は、記憶回路２５から逆補正対応表を読み出してもよい。このとき、制御回路２９は、逆補正対応表にしたがって、減衰定量用部分データをゲイン補正前の減衰定量用部分データに変換する。すなわち、ゲイン逆補正機能２９５は、減衰定量用部分データに対するゲイン補正をキャンセルする。

ゲイン逆補正機能２９５により、制御回路２９は、超音波プローブ３が減衰定量用超音波の反射波を受信した時点における受信データ、すなわち深さ方向に沿った深さ方向の超音波信号の純粋な強度を復元する。なお、超音波受信回路１３からゲイン未補正の減衰定量用部分データ（ローデータ）を出力することが可能であれば、上記ゲイン逆補正機能２９５は不要となる。

音場特性補正機能２９７を実現する制御回路２９は、記憶回路２５から音場特性補正データを読み出す。制御回路２９は、ゲイン補正前の減衰定量用部分データから音場特性補正データを減算する。これにより、制御回路２９は、音場特性に依存せず、かつゲイン補正前の減衰定量用部分データ（以下、減衰データと呼ぶ）を生成する。なお、制御回路２９は、記憶回路２５から音場特性補正対応表を読み出してもよい。このとき、制御回路２９は、音場特性補正対応表にしたがって、ゲイン補正前の減衰定量用部分データを、減衰データに変換する。すなわち、制御回路２９は、超音波スキャンにおける音場特性に基づいて、ゲイン補正前の減衰定量用部分データにおける音場特性の依存性をキャンセルする。

音場特性補正機能２９７により、制御回路２９は、超音波ビームの形状、および超音波プローブ３の形状に特有な画素値（または輝度値）の変化量を、ゲイン補正前の減衰定量用Ｂモードデータから排除する。この排除により、制御回路２９は、被検体内の組織により超音波が純粋に減衰する程度を反映した画素値（または輝度値）を有する減衰データを生成する。減衰データは、ゲイン逆補正機能２９５および音場特性補正機能２９７により補正された補正データに対応する。

解析機能２９９を実行する制御回路２９は、被検体Ｐに対する超音波スキャンにより得られた上記減衰データに基づいて、推定された構造物の位置を除く被検体内の複数の位置（以下、非構造物領域と呼ぶ）における組織性状を解析する。すなわち、解析機能２９９を実現する制御回路２９は、被検体内を伝搬する減衰定量用超音波の減衰量を計算する。

具体的には、制御回路２９は、減衰データにおいて、非構造物領域における画素値（または輝度値）において、深さ方向に沿った微分値を計算する。上記微分値は、例えば、非構造物領域における複数の画素各々において、深さ方向に沿って隣接する２つの画素の画素値の差分値を、この２つの画素の間隔（距離）で除した値に相当する。

次いで、制御回路２９は、計算した微分値に対して、超音波の往復分を加味して、１／２を乗算する。これにより、制御回路２９は、減衰定量用超音波の送信中心周波数に関する超音波の減衰量（ｄＢ／ｃｍ）を、非構造物領域に含まれる複数の画素（位置）について計算する。さらに、制御回路２９は、計算した減衰量を減衰定量用超音波の送信中心周波数で除算する。この除算により、制御回路２９は、非構造物領域に含まれる複数の画素（位置）における減衰定数（ｄＢ／ｃｍ／Ｈｚ）を計算する。

減衰定数は、減衰定量用超音波の送信中心周波数に依存しない。また、減衰定数は、非構造物領域に含まれる複数の画素（位置）ごとにほぼ一定の値となる。減衰定数は、診断対象となる組織の性状を表す数値である。上記各種計算により、制御回路２９は、非構造物領域に含まれる複数の位置において減衰定数を示す解析データ（以下、減衰定数データと呼ぶ）を生成する。制御回路２９は、減衰定数データを、画像生成回路１９および記憶回路２５等に出力する。

組織性状は、輝度値または画素値の解析による減衰定数に限定されず、他の解析により減衰量を反映した量であってもよい。例えば、周波数が異なる複数の減衰定量用超音波が一つの走査線に対して送受信された場合、解析機能２９９を実現する制御回路２９は、非構造物領域に含まれる複数の位置におけるドプラデータに基づいて周波数解析を実行することにより、非構造物領域に含まれる複数の位置における組織性状を示すパラメータを計算してもよい。具体的には、制御回路２９は、複数の減衰定量用超音波各々に対応するドプラデータを用いて、周波数の違いに起因する減衰量の差を、組織性状として計算する。

また、Ｂモードに関する超音波送信における周波数帯域より狭帯域の超音波が減衰定量用超音波として被検体Ｐに送信された場合、制御回路２９は、狭帯域の減衰定量用超音波の受信により取得されたデータに対して周波数解析を実行する。次いで、制御回路２９は、周波数解析における周波数特性に基づいて、組織性状として減衰量を計算する。

表示機能を実現する制御回路２９は、解析結果である組織性状を、モニタ７に表示させる。具体的には、制御回路２９は、減衰定数データを画像生成回路１９に出力する。画像生成回路１９は、減衰定数データに基づいて減衰定量画像を生成する。このとき、制御回路２９は、減衰定量画像における減衰定数に応じた色相を付与するために、画像生成回路１９を制御する。この制御により、減衰定量画像は、減衰定数に応じた色相を有する。

なお、構造物の位置に対応する領域において、減衰定量画像は、欠落した状態となる。すなわち、構造物の位置に対応する領域に対応する減衰定量画像は、存在しない。画像生成回路１９は、減衰定量画像を画像合成回路２３に出力する。

画像合成回路２３は、制御回路２９による制御のもとで、減衰定量用Ｂモードデータの収集に関する被走査領域と略同一の被走査領域に関するＢモード画像と、減衰定量画像との位置合わせ（ｒｅｇｉｓｔｒａｔｉｏｎ）を実行する。画像合成回路２３は、制御回路２９による制御のもとで、減衰定量画像を所定の不透明度または透明度に変換する。

画像合成回路２３は、Ｂモード画像に、所定の不透明度または透明度を有する減衰定量画像を重畳させた減衰重畳画像を生成する。画像合成回路２３は、減衰重畳画像に、減衰定数の色相に応じた凡例等を合成する。画像合成回路２３は、凡例等を合成した減衰重畳画像を、モニタ７に出力する。

モニタ７は、構造物の位置を除く被検体内の複数位置各々において、解析結果としての減衰量を表示する。具体的には、モニタ７は、凡例等を合成した減衰重畳画像を表示する。

図４は、モニタ７に表示された減衰重畳画像を凡例とともに示す図である。図４に示すように、減衰重畳画像において、減衰定量画像は、所定の不透明度または透明度でＢモード画像に重畳された状態で表示される。図４に示すように、減衰重畳画像において、推定された構造物の位置を含む領域は、図３に示すように大きい分散比を有し、減衰定数が計算されない領域、すなわち減衰定量画像が存在しない領域である。この領域には背面のＢモード画像が表示される。すなわち、図４において推定された構造物の位置を含む領域は、血管、腹膜、横隔膜等の構造物およびこれらの構造物の付近であって、減衰定数を示す減衰定量画像（カラー画像）が非表示となる領域である。

以下、構造物推定機能および組織性状解析機能に関する処理手順について説明する。図５は、構造物推定機能および組織性状解析機能に関する処理手順の一例を示すフローチャートである。

減衰定量開始指示が入力装置５により入力される（ステップＳａ１）。減衰定量開始指示の入力を契機として、Ｂモード用超音波が、被検体Ｐに対して送受信される（ステップＳａ２）。なお、ステップＳａ１の処理は、ステップＳａ２の処理の後に実行されてもよい。このとき、操作者は、被検体Ｐに対して減衰定数の解析に関するＢモード画像が表示されると、減衰定量開始指示を入力する。加えて、ステップＳａ２の処理に関連して発生された、図２に示すようなＢモード画像は、記憶回路２５に記憶されてもよい。

１フレーム分、すなわち１被走査領域に対応するＢモードデータが収集されると、減衰定量用超音波が、被検体Ｐに対して送受信される（ステップＳａ３）。次いで、減衰定量用Ｂモードデータが発生される。減衰定量用Ｂモードデータに基づいて、構造物の位置が推定される（ステップＳａ４）。なお、構造物の位置の推定は、ステップＳａ２の処理の後に発生されたＢモードデータに基づいて推定されてもよい。このとき、構造物の位置の推定に用いられた分散比をＢモード画像に重畳させた分散比重畳画像（図３参照）が、モニタ７に表示されてもよい。

推定された構造物の位置を除く複数位置の組織性状（例えば、減衰定数）が、減衰定量用Ｂモードデータに基づいて計算される（ステップＳａ５）。ステップＳａ５の処理により、Ｂモード画像に関する被走査領域と略同一の被走査領域における減衰定数データが発生される。構造物の位置を除く複数の位置各々における減衰定数を示す減衰定数データに基づいて、減衰定量画像が生成される（ステップＳａ６）。ステップＳａ２の処理において受信された受信信号に基づいてＢモードデータが発生される。次いで、Ｂモードデータに基づいてＢモード画像が発生される（ステップＳａ７）。

Ｂモード画像と減衰定量画像とが位置合わせされる。加えて、減衰定量画像に対して所定の不透明度または透明度が付与される。所定の不透明度または透明度を有する減衰定量画像をＢモード画像に重畳した減衰重畳画像が生成される。減衰重畳画像がモニタ７に表示される（ステップＳａ８）。

（変形例）
上記実施形態との相違は、減衰定量用Ｂモードデータの収集に関する被走査領域に対応する減衰定量画像を生成し、被走査領域において推定された構造物の位置における減衰定量画像を減衰重畳画像において非表示とすることにある。本変形例では、減衰定量用Ｂモードデータにおいて、減衰定量用部分データを推定する処理は不要となる。

（組織性状解析機能）
ゲイン逆補正機能２９５を実現する制御回路２９は、減衰定量用Ｂモードデータから、逆補正データを減算する。これにより、制御回路２９は、被走査領域全域に亘るゲイン補正前の減衰定量用Ｂモードデータを生成する。

音場特性補正機能２９７を実現する制御回路２９は、ゲイン補正前の減衰定量用Ｂモードデータから音場特性補正データを減算する。これにより、制御回路２９は、音場特性に依存せず、かつゲイン補正前の減衰定量用Ｂモードデータ（以下、減衰Ｂモードデータと呼ぶ）を生成する。

解析機能２９９を実行する制御回路２９は、減衰Ｂモードデータに基づいて、被走査領域全域における組織性状を解析する。具体的には、制御回路２９は、減衰Ｂモードデータにおける画素値（または輝度値）において、深さ方向に沿った微分値を計算する。制御回路２９は、計算した微分値に対して超音波の往復分を加味して１／２を乗算した値（減衰量）を減衰定量用超音波の送信中心周波数で除算する。この除算により、制御回路２９は、収集領域における複数の画素（位置）各々における減衰定数を計算する。

上記各種計算により、制御回路２９は、被走査領域全域において減衰定数を示す解析データ（減衰定数データ）を生成する。制御回路２９は、減衰定数データを、画像生成回路１９および記憶回路２５に出力する。

表示機能を実現する制御回路２９は、減衰定数データを画像生成回路１９に出力する。画像生成回路１９は、減衰定数データに基づいて、減衰定量画像における減衰定数に応じた色相を付与した減衰定量画像を生成する。本変形例における減衰定量画像の大きさは、被走査領域の大きさに略一致する。

画像合成回路２３は、制御回路２９による制御のもとで、減衰定量画像を所定の不透明度または透明度に変換する。画像合成回路２３は、Ｂモード画像に、所定の不透明度または透明度を有する減衰定量画像を重畳させた減衰重畳画像を生成する。制御回路２９は、減衰重畳画像の表示において、被走査領域において推定された構造物の位置に対応する減衰定量画像の一部領域を非表示とするように、画像合成回路２３またはモニタ７を制御する。このとき、モニタ７は、例えば図４に示すような態様で、減衰重畳画像を表示する。すなわち、モニタ７は、推定機能２９３により推定された構造物の位置を除く被検体内の複数位置について、解析機能２９９による解析結果を表示する。

なお、制御回路２９は、上記一部領域を所定の色相でマスクして減衰定量画像を表示するように、画像合成回路２３またはモニタ７を制御してもよい。このとき、モニタ７は、上記一部領域を所定の色相でマスクして、減衰定量画像を表示する。

なお、画像生成回路１９は、制御回路２９による制御のもとで、被走査領域において推定された構造物の位置に対応する減衰定量画像の一部領域を、減衰定量画像から取り除いた部分減衰画像を生成してもよい。このとき、画像合成回路２３は、制御回路２９による制御のもとで、部分減衰画像を所定の不透明度または透明度に変換する。画像合成回路２３は、Ｂモード画像に、所定の不透明度または透明度を有する部分減衰画像を重畳させた減衰重畳画像を生成する。モニタ７は、減衰重畳画像を表示する。

以下、本変形例に係る構造物推定機能および組織性状解析機能に関する処理手順について説明する。図６は、本変形例に係る構造物推定機能および組織性状解析機能に関する処理手順の一例を示すフローチャートである。ステップＳｂ１乃至ステップＳｂ４、およびステップＳｂ７の処理は、図５におけるステップＳａ１乃至ステップＳａ４、およびステップＳａ７の処理とそれぞれ同様なため、説明は省略する。

減衰定数が、減衰定量用Ｂモードデータに基づいて、収集領域にわたって計算される（ステップＳｂ５）。ステップＳｂ５の処理により、Ｂモード画像に関する被走査領域と略同一の被走査領域における減衰定数データが発生される。収集領域における減衰定数データに基づいて、減衰定量画像が生成される（ステップＳｂ６）。

所定の不透明度または透明度を有する減衰定量画像をＢモード画像に重畳した減衰重畳画像が生成される（ステップＳｂ８）。減衰重畳画像における減衰定量画像において構造物の位置に対応する領域を非表示で、または所定の色相でマスクして、減衰重畳画像が、モニタ７に表示される（ステップＳｂ９）。

以上に述べた構成によれば、以下の効果を得ることができる。
本実施形態および本変形例に係る超音波診断装置１によれば、構造物の有無を判定することにより、組織性状（例えば、減衰定数などの減衰量、弾性率、粘性率、硬さなど）を非構造物領域に関して表示することができる。すなわち、本超音波診断装置１によれば、組織性状をカラー画像として表示する際、構造物ありと判定された領域ではカラー非表示となり、背景のＢモード画像を表示することができる。あるいは、本超音波診断装置１によれば、組織性状をカラー画像として表示する際、構造物ありと判定された領域では、カラー画像とは異なる所定の色相でマスクして組織性状を表示することができる。

例えば、本実施形態および本変形例においてモニタ７に表示される図４における減衰重畳画像と従来技術における図１１における画像とを比較すると、図１１における画像においては、構造物の有無によらず全域で減衰定数を算出した結果として、減衰定数が低い値から高い値まで様々な部分が表示される、このため、図１１における画像では、操作者が着目する組織の減衰定数がどの程度なのか判断しにくい状態となっている。一方、本実施形態および本変形例においてモニタ７に表示される図４における減衰重畳画像では、構造物の位置に対応する領域における組織性状は非表示となるため、操作者は、一見して組織性状を示す減衰の程度を把握することが可能となる。

以上のことから、本実施形態および本変形例に係る超音波診断装置１によれば、表示された画像や組織性状の解析結果の数値が妥当であるか否かを操作者が判断することなく、推定された構造物の位置を除く領域において信頼性のある解析結果が表示されるため、得られた解析結果の精度や再現性などが向上させることができる。加えて、本実施形態および本変形例に係る超音波診断装置１によれば、解析結果の表示画像において、表示された画像が妥当であるか否かという判断を操作者が判断する必要がないため、操作者にとって、時間的・精神的な負担を軽減させることでき、診断精度を向上させることができる。

（第１の応用例）
本実施形態および本変形例との相違は、モニタ７に表示された減衰重畳画像に設定されたＲＯＩ（以下、計測ＲＯＩと呼ぶ）における減衰定数の代表値の信頼性を判定し、判定結果に応じて代表値を計算して表示することにある。計測ＲＯＩは、計測ＲＯＩにおける代表値の計算、および代表値に対する信頼性の判定に関する領域である。

設定機能２９１は、モニタ７に表示された減衰重畳画像上に、計測ＲＯＩを設定する機能を有する。なお、減衰定量画像がモニタ７に表示されている場合、設定機能２９１は、減衰定量画像上に計測ＲＯＩを設定してもよい。また、Ｂモード画像がモニタ７に表示されている場合、設定機能２９１は、Ｂモード画像上に計測ＲＯＩを設定してもよい。

解析機能２９９は、計測ＲＯＩに含まれる複数の減衰定数を代表する代表値に対する信頼度を計算する機能と、信頼度と信頼性判定閾値とを比較することにより代表値に対する信頼性を判定する機能と、信頼度が信頼性判定閾値より大きい場合、計測ＲＯＩに関する代表値を計算する機能等とを有する。

表示機能は、入力装置５を介した計測開始の指示（以下、計測開始指示と呼ぶ）の入力に応答して減衰重畳画像上に計測ＲＯＩをモニタ７に表示させる機能と、解析機能２９９による判定結果に応じて代表値または所定の通知をモニタ７に表示させる機能等とを有する。

記憶回路２５は、本応用例に係る各種機能に関するプログラムを記憶する。記憶回路２５は、信頼性判閾値および所定の通知を記憶する。信頼性判定閾値とは、解析機能２９９により計算される代表値の信頼性を判定する閾値である。所定の通知とは、例えば、信頼度が信頼性判定閾値以下である場合、代表値が信頼性を有さないこと、代表値が適切でないこと、および代表値の計測が不可であることなどを、モニタ７を介して操作者等に通知するための文字列等である。

入力装置５は、操作者の指示により計測開始指示を入力する。入力装置５は、トラックボールやパネルボタンなどのインターフェースを介した操作者の操作にしたがって、モニタ７に表示された減衰重畳画像上に計測ＲＯＩの位置の確定指示を入力する。

表示機能を実現する制御回路２９は、入力装置５を介した操作者の計測開始指示に応答して、モニタ７に表示された減衰重畳画像上に、計測ＲＯＩを表示させる。設定機能２９１を実現する制御回路２９は、入力装置５を介した操作者の確定指示に応答して、減衰重畳画像上に、計測ＲＯＩを設定する。

解析機能２９９を実現する制御回路２９は、計測ＲＯＩに含まれる画素値に基づいて、計測ＲＯＩに含まれる複数の減衰定数を代表する代表値に対する信頼度を計算する。代表値は、例えば、減衰定量画像上の計測ＲＯＩに含まれる複数の減衰定数の平均値、中央値、最頻値などである。

具体的には、制御回路２９は、計測ＲＯＩに含まれる画素数（以下、計測画素数と呼ぶ）を、計測ＲＯＩの面積として決定する。制御回路２９は、計測ＲＯＩに含まれる複数の分散比と構造物判定閾値とを比較し、構造物判定閾値より大きい分散比に対応する画素数を、計測ＲＯＩにおける構造物の面積として決定する。制御回路２９は、計測ＲＯＩの面積に対する構造物の面積の割合（以下、構造物面積比と呼ぶ）を計算する。このとき、信頼度は、例えば、１から面積比を差分した差分値となる。

なお、制御回路２９は、計測ＲＯＩの面積に対する非構造物領域の面積の割合（以下、非構造物面積比と呼ぶ）を計算してもよい。信頼度は、非構造物面積比に相当する。このとき、制御回路２９は、計測ＲＯＩに含まれる複数の分散比と構造物判定閾値とを比較し、構造物判定閾値より小さい分散比に対応する画素数を、計測ＲＯＩにおける非構造物の面積として決定する。

また、信頼度は、分散比を用いて計算された構造物面積比または非構造物面積比に基づくことに限定されない。すなわち、信頼度は、計測画素数に対する、減衰定数を示す色相を有する画素数の割合であってもよい。また、信頼度は、式（２）で計算された分散比を、信頼度としてもよい。

なお、信頼度は、計測画素数に対するＢモード画像上の計測ＲＯＩに含まれる領域において所定の輝度値より高い輝度値に関する画素数を、構造物面積比として計算されてもよい。このとき、所定の輝度値は、記憶回路２５に記憶される。所定の輝度値とは、例えば、血管壁などの構造物を代表する輝度値である。

また、信頼度は、計測画素数に対するＢモード画像上の計測ＲＯＩに含まれる領域において所定の輝度値より低い輝度値に関する画素数を、非構造物面積比として計算されてもよい。Ｂモードデータを用いて信頼度を計算した場合、制御回路２９は、信頼度の計算に用いられたＢモードデータを用いて、ゲイン逆補正機能２９５、音場特性補正機能２９７および解析機能２９９により、減衰定数データを生成する。

制御回路２９は、信頼度と信頼性判定閾値とを比較することにより、計測ＲＯＩに関する代表値に対する信頼性を判定する。具体的には、制御回路２９は、信頼度が信頼性判定閾値より大きい場合、計測ＲＯＩにおける代表値の計測が適切であると判定する。制御回路２９は、信頼度が信頼性判定閾値以下である場合、計測ＲＯＩにおける代表値の計測が不適であると判定する。

計測ＲＯＩにおける代表値の計測が適切であると判定された場合、制御回路２９は、計測ＲＯＩ内の複数の位置各々における減衰量に基づいて、計測ＲＯＩ内における減衰量の代表値を計算する。具体的には、制御回路２９は、計測ＲＯＩに含まれる複数の減衰定数に基づいて、減衰定数の平均値、中央値、最頻値等の代表値を計算する。なお、制御回路２９は、減衰データにおける計測ＲＯＩの平均値を深さ方向に沿って微分した微分値（以下、平均微分値と呼ぶ）を用いて、代表値を計算してもよい。このとき、例えば、以下のような計算手順となる。

まず、制御回路２９は、減衰データにおいて計測ＲＯＩに含まれる複数の走査線（以下、計測走査線と呼ぶ）に沿った複数の画素値を平均することにより、計測平均値を生成する。次いで、制御回路２９は、超音波プローブ３が被検体Ｐの体表面に当接された当接面と計測ＲＯＩとの間における計測走査線に沿った複数の画素値を平均することにより、ＲＯＩ上方平均値を生成する。制御回路２９は、ＲＯＩ上方平均値から計測平均値を差分した値を、計測走査線の中心線に沿った計測ＲＯＩの厚みで除算する。この除算による値が、上記平均微分値に相当する。最後に、制御回路２９は、平均微分値に対して超音波の往復分を加味して１／２を乗算した値を減衰定量用超音波の送信中心周波数で除算することにより、代表値を算出する。

なお、制御回路２９は、計測ＲＯＩより深部の領域であって計測走査線に沿った複数の画素値を平均することにより、ＲＯＩ下方平均値を生成してもよい。このとき、制御回路２９は、計測平均値からＲＯＩ下方平均値を差分した値を、計測走査線の中心線に沿った計測ＲＯＩの厚みで除算する。この除算による値が、上記平均微分値に相当する。

制御回路２９は、モニタ７における計測結果表示領域に代表値を表示させるように、画像合成回路２３またはモニタ７を制御する。計測結果表示領域とは、モニタ７における画像の表示領域のうち、計測ＲＯＩが重畳された減衰重畳画像の表示領域とは異なる領域である。モニタ７は、計測ＲＯＩが重畳された減衰重畳画像の表示に加えて、代表値と計測ＲＯＩの名称とを、計測結果表示領域に表示する。

計測ＲＯＩにおける代表値の計測が不適であると判定された場合、制御回路２９は、計測結果表示領域に所定の通知を表示させるように、画像合成回路２３またはモニタ７を制御する。モニタ７は、計測ＲＯＩが重畳された減衰重畳画像の表示に加えて、所定の通知と計測ＲＯＩの名称とを、計測結果表示領域に表示する。

図７は、計測ＲＯＩ１における代表値の計測が適切であると判定され、計測ＲＯＩ２における代表値の計測が不適であると判定された場合において、計測ＲＯＩの名称に対応する代表値および所定の通知と、計測ＲＯＩが重畳された減衰重畳画像とがモニタ７に表示される一例を示す図である。

図７に示すように、計測ＲＯＩ１は、一様な肝臓の実質を含む領域上に設定され、計測ＲＯＩ１内では減衰定量画像が表示されている。このような場合には、計測ＲＯＩ１内には構造物がほとんど存在しないため、計測ＲＯＩ１内において精度よく減衰定数を算出することができる。このため、図７に示すように、計測結果表示領域ＭＤＲには、計測ＲＯＩの名称とともに算出された代表値が表示される。

一方、図７に示すように、計測ＲＯＩ２は、血管領域を含む領域に設定され、計測ＲＯＩ２内では過半数の領域で減衰定量画像が表示されていない。このような場合には、計測ＲＯＩ２において精度よく減衰定数を算出することができないため、代表値は算出されない。このため、図７に示すように、計測結果表示領域ＭＤＲには、計測ＲＯＩの名称とともに、所定の通知「Ｎｏｔ ａｐｐｒｏｐｒｉａｔｅ！」（不適）が表示される。

（代表値計算機能）
代表値計算機能とは、計測ＲＯＩにおける減衰定数の代表値の信頼性を判定し、判定結果に応じて代表値を計算して表示する機能である。代表値計算機能は、上述した設定機能２９１、解析機能２９９および表示機能を包括的に実現する機能を有する。なお、代表値計算機能を実現する制御回路２９（プロセッサ）は、代表値計算部として機能してもよい。

図８は、本応用例に係る代表値計算機能に関する処理手順の一例を示すフローチャートである。図８に示すフローチャートは、図５におけるフローチャートまたは図６におけるフローチャートの続きとして実行される処理手順である。

計測開始指示が入力装置５により入力される（ステップＳｃ１）。計測開始指示の入力に応答して、モニタ７に表示された減衰重畳画像上に、計測ＲＯＩが表示される。確定指示の入力に応答して、減衰重畳画像上に、計測ＲＯＩが設定される（ステップＳｃ２）。設定された計測ＲＯＩ内において構造物が占める面積の割合（構造物面積比）割合が算出される（ステップＳｃ３）。構造物面積比に基づいて信頼度が決定される（ステップＳｃ４）。

信頼度が信頼性判定閾値より大きければ（ステップＳｃ５）、計測ＲＯＩにおける減衰定数の代表値が算出される（ステップＳｃ６）。次いで、計測結果表示領域に代表値が表示される（ステップＳｃ７）。信頼度が信頼性判定閾値より小さければ（ステップＳｃ５）、計測結果表示領域に計測不可であることを示す通知が表示される（ステップＳｃ８）。計測ＲＯＩが設定されるごとに、ステップＳｃ３乃至ステップＳｃ８の処理が繰り返される。なお、複数の計測ＲＯＩが減衰重畳画像上に設定された場合、複数の計測ＲＯＩ各々に対して、ステップＳｃ３乃至ステップＳｃ８の処理が繰り返される。

以上に述べた構成によれば、本実施形態および本変形例に係る効果に加えて、以下の効果を得ることができる。
本応用例に係る超音波診断装置１によれば、操作者により指定された計測ＲＯＩにおける減衰量の代表値を表示する際、計測ＲＯＩに含まれる構造物が多い場合および計測ＲＯＩに含まれる構造物が大きい場合、すなわち計測ＲＯＩにおける組織性状を示す数値の代表値の信頼性を示す信頼度が信頼性判定閾値より小さい場合、代表値を表示せずに、計測不可と表示することができる。すなわち、本応用例に係る超音波診断装置１によれば、計測ＲＯＩ内に含まれる構造物が多い場合や大きい場合は、計測ＲＯＩにおける減衰定数の代表値は信頼できないとして表示せず、計測ＲＯＩを再設定するように、操作者を促すことが可能となる。

以上のことから、本応用例に係る超音波診断装置１によれば、信頼性のある定量結果（組織性状を示す減衰定数などの代表値）が表示されるため、表示された代表値が妥当か否かを操作者が判断する必要がなく、操作者の時間的・精神的な負担が軽減し、加えて、組織性状の計測の再現性や精度を向上させることができる。

（第２の応用例）
第１の応用例との相違は、計測ＲＯＩ内における構造物の有無または構造物の面積等に応じた信頼度によらずに計測ＲＯＩに対応する代表値を計算し、計測ＲＯＩに対応する信頼度と代表値とを表示することにある。

記憶回路２５は、信頼度をランク付けするための少なくとも一つの閾値（以下、ランク閾値と呼ぶ）を記憶する。ランク閾値の数は、ランクの数より１小さい。記憶回路２５は、信頼度のランク付けに関するプログラム（以下、ランク付けプログラムと呼ぶ）を記憶する。なお、信頼度そのものを表示する場合には、ランク閾値およびはランク付けプログラム不要となる。

解析機能２９９を実現する制御回路２９は、信頼度によらずに代表値を計算する。制御回路２９は、計測ＲＯＩに対応する信頼度を算出する。制御回路２９は、ランク閾値とランク付けプログラムを、記憶回路２５から読み出す。制御回路２９は、ランク付けプログラムしたがって、信頼度とランク閾値とを比較することにより、信頼度をランク付けする。以下、説明を簡単にするために、ランク閾値は２つ（以下、高ランク閾値および低ランク閾値と呼ぶ）であるものとする。このとき、ランクの数は、３（例えば、Ａ、Ｂ，Ｃの３ランク）となる。なお、ランク閾値の数は、２に限定されない。すなわち、ランクの数は、３に限定されず、任意に設定可能である。

制御回路２９は、信頼度が高ランク閾値より大きい場合、計測ＲＯＩに関する信頼度のランクを、「Ａ」として決定する。制御回路２９は、信頼度が高ランク閾値以下であって低ランク閾値より大きい場合、計測ＲＯＩに関する信頼度のランクを、「Ｂ」として決定する。制御回路２９は、信頼度が低ランク閾値以下である場合、計測ＲＯＩに関する信頼度のランクを、「Ｃ」として決定する。制御回路２９は、決定したランクを、信頼度に関連する計測ＲＯＩに対応付けて、記憶回路２５に記憶させる。制御回路２９は、モニタ７における計測結果表示領域に代表値とランクとを表示させるように、画像合成回路２３またはモニタ７を制御する。

モニタ７は、計測ＲＯＩが重畳された減衰重畳画像の表示に加えて、代表値とランクと計測ＲＯＩの名称とを、計測結果表示領域に表示する。なお、モニタ７は、信頼度そのものを、代表値および計測ＲＯＩの名称とともに計測結果表示領域に表示してもよい。このとき、モニタ７は、例えば、式（２）で計算された分散比、または非構造物面積比、１から構造物面積比を差分した差分値などを、信頼度として表示する。

図９は、計測ＲＯＩの名称に対応する代表値および信頼度のランクと、計測ＲＯＩが重畳された減衰重畳画像とがモニタ７に表示される一例を示す図である。図９は、信頼度を、Ａ、Ｂ、Ｃと三段階で表示する例を示している。図９に示すように、計測ＲＯＩ１は、一様な肝臓の実質を含む領域上に設定され、計測ＲＯＩ１内では減衰定量画像が表示されている。このような場合には、計測ＲＯＩ１内には構造物がほとんど存在しないため、信頼度のランクは、Ａとして決定され、計測ＲＯＩ１に関する代表値とともに計測結果表示領域ＭＤＲに表示される。

図９に示すように、計測ＲＯＩ２は、血管領域を含む領域に設定され、計測ＲＯＩ２内では過半数の領域で減衰定量画像が表示されていない。このような場合には、計測ＲＯＩ２における過半数の領域において構造物が存在するため、信頼度のランクは、Ｃとして決定され、計測ＲＯＩ２に関する代表値とともに計測結果表示領域ＭＤＲに表示される。

（信頼度表示機能）
信頼度表示機能とは、信頼度もしくは信頼度のランクを、計測ＲＯＩ内における構造物の有無または構造物の面積等に応じた信頼度によらずに計算された代表値とともに、計測結果表示領域ＭＤＲに表示する機能である。信頼度表示機能は、第１の応用例における設定機能２９１、解析機能２９９および表示機能を包括的に実現する機能を有する。

図１０は、本応用例に係る信頼度表示機能に関する処理手順の一例を示すフローチャートである。図１０に示すフローチャートは、図５におけるフローチャートまたは図６におけるフローチャートの続きとして実行される処理手順である。ステップＳｄ１乃至ステップＳｄ４、およびステップＳｄ６の処理は、図８におけるステップＳｃ１乃至ステップＳｃ４、およびステップＳｃ６の処理とそれぞれ同様なため、説明は省略する。

ステップＳｄ４の処理の後、信頼度とランク閾値とを比較することにより、信頼度のランクが決定される（ステップＳｄ５）。なお、信頼度そのものを代表値とともに表示する場合には、ステップＳｄ５の処理は不要となる。

ステップＳｄ６の処理の後、計測結果表示領域に減衰定数の代表値と信頼度のランクとが、計測ＲＯＩの名称とともに表示される（ステップＳｄ７）。計測ＲＯＩが設定されるごとに、ステップＳｄ３乃至ステップＳｄ７の処理が繰り返される。なお、複数の計測ＲＯＩが減衰重畳画像上に設定された場合、複数の計測ＲＯＩ各々に対して、ステップＳｄ３乃至ステップＳｄ６の処理が繰り返され、複数の計測ＲＯＩにそれぞれ対応する複数の代表値および複数のランクが、計測結果表示領域に一覧表示される。

以上に述べた構成によれば、本実施形態および本変形例に係る効果に加えて、以下の効果を得ることができる。
本応用例に係る超音波診断装置１によれば、計測ＲＯＩに含まれる減衰量の代表値を表示する際、計測ＲＯＩに含まれる構造物の大きさや程度に対応する信頼度を併せて表示することができる。すなわち、本応用例に係る超音波診断装置１によれば、計測ＲＯＩにおいて対象組織の組織性状によって信頼性のある結果が得られる場所がほとんど見つからない場合においても、何らかの定量結果（代表値）を得ることができる。例えば、線維化が進行した症例などの肝臓全体が不均一になっている症例などでは、どの場所においても減衰定量結果（代表値）を、信頼度とともに表示することができる。

以上のことから、本応用例に係る超音波診断装置１によれば、定量結果の信頼性の程度を表示することができるため、操作者は、安心して、定量結果を参考情報として参照して、診断に用いることができる。

（第３の応用例）
本実施形態、本変形例、第１の応用例および第２の応用例との相違は、医用診断装置として、ＭＲＩ装置を適用することにある。不図示のＭＲＩ装置における計算機は、画像生成回路１９、画像メモリ２１、画像合成回路２３、記憶回路２５、インターフェース回路２７、制御回路２９、不図示のシーケンサ等を有する。計算機には、入力装置５、モニタ７等が接続される。

計算機における制御回路２９は、上述したように、各種機能を実行する。計算機における記憶回路２５は、制御回路２９で実行される各種機能に関するプログラムなどの各種データを記憶する。

本応用例における推定機能２９３を実現する制御回路２９は、磁気共鳴現象を利用したＭＲスキャンにより得られたデータ（以下、ＭＲデータと呼ぶ）に基づいて、ＭＲスキャンに対応する断面において、上述したように構造物の位置を推定する。

本応用例における解析機能２９９を実現する制御回路２９は、核磁気共鳴現象を利用したＭＲスキャンにより得られたデータに基づいて、断面において推定機能２９３により推定された構造物の位置を除く被検体内の複数位置の組織性状を解析する。本応用例における組織性状とは、例えば、減衰量、弾性率（ヤング率）、粘性率、歪みなどの診断対象組織の性状を示す量である。本応用例において、例えば、組織性状は、ＭＲエラストグラフィー（ＭＲ ｅｌａｓｔｏｇｒａｐｈｙ：ＭＲＥ）により取得される。このとき、組織性状の解析に関するデータは、ＭＲＥスキャンにより得られたＭＲＥデータに対応する。

以下、説明を簡単にするために、ＭＲＥは、被検体Ｐの外部に設けられた非磁性な加振機構による振動に基づき、被検体Ｐの内部を伝播する波を、プロトンの位相を画像化した位相画像を用いて観測するエラストグラフィー法とする。

計算機における画像生成回路１９は、被検体Ｐに対してシーケンサによる制御のもとで実行されたＭＲスキャンにより取得されたＭＲデータに基づいて、ＭＲ画像を生成する。また、画像生成回路１９は、被検体Ｐに対してシーケンサによる制御のもとで実行されたＭＲＥスキャンにより取得されたＭＲＥデータに基づいて、ＭＲＥ画像を生成する。ＭＲ画像またはＭＲＥ画像は、設定機能２９１、推定機能２９３に用いられる。また、ＭＲＥ画像は、解析機能２９９に用いられる。

本応用例においてモニタ７に表示される減衰重畳画像は、図４における背景画像がＭＲ画像であって、減衰定量画像がＭＲＥ画像に対応する。本応用例に係る構造物推定機能および組織性状解析機能に関する処理手順は、ほぼ図５と同様なため、異なる処理について説明する。なお、ＭＲ画像とＭＲＥ画像は、重畳表示されず、個別に表示されても良い。

本応用例において、図５のステップＳａ２に対応する処理は、ＭＲデータを収集するための撮像シーケンスにしたがって被検体Ｐに対してＭＲスキャンを実行する処理となる。本応用例において、図５のステップＳａ３に対応する処理は、ＭＲＥデータを収集するための撮像シーケンスにしたがって被検体に対してＭＲＥスキャンを実行する処理となる。

本応用例において、図５のステップＳａ４に対応する処理は、ＭＲＥデータに基づいて構造物の位置を推定する処理となる。本応用例において、図５のステップＳａ５に対応する処理は、構造物の位置を除く複数の位置におけるＭＲＥデータに基づいて減衰定数を計算する処理となる。

本応用例において、図５のステップＳａ７に対応する処理は、ＭＲ画像を生成する処理となる。本応用例において、図５のステップＳａ８に対応する処理は、ＭＲ画像に減衰定量画像を重畳して表示する処理となる。

なお、本応用例に係るＭＲＩ装置によれば、本実施形態における変形例、第１の応用例、および第２の応用例についても同様に実行することができる。

以上に述べた構成によれば、本実施形態および本変形例に係る効果と同様に、以下の効果を得ることができる。
本応用例に係るＭＲＩ装置によれば、構造物の有無を判定することにより、組織性状を非構造物領域に関して表示することができる。すなわち、本ＭＲＩ装置によれば、組織性状をカラー画像として表示する際、構造物ありと判定された領域ではカラー非表示となり、背景のＭＲ画像を表示することができる。あるいは、本ＭＲＩ装置によれば、組織性状をカラー画像として表示する際、構造物ありと判定された領域では、カラー画像とは異なる所定の色相でマスクして組織性状を表示することができる。

以上のことから、本応用例に係るＭＲＩ装置によれば、表示された画像や組織性状の解析結果の数値が妥当であるか否かを操作者が判断することなく、推定された構造物の位置を除く領域において信頼性のある解析結果が表示されるため、得られた解析結果の精度や再現性などが向上させることができる。加えて、本応用例に係るＭＲＩによれば、解析結果の表示画像において、表示された画像が妥当であるか否かという判断を操作者が判断する必要がないため、操作者にとって、時間的・精神的な負担を軽減させることでき、診断精度を向上させることができる。

加えて、実施形態に係る各機能は、当該処理を実行するプログラムをワークステーション等のコンピュータにインストールし、これらをメモリ上で展開することによっても実現することができる。このとき、コンピュータに当該手法を実行させることのできるプログラムは、磁気ディスク（ハードディスクなど）、光ディスク（ＣＤ−ＲＯＭ、ＤＶＤなど）、半導体メモリなどの記憶媒体に格納して頒布することも可能である。

本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。

１…超音波診断装置、３…超音波プローブ、５…入力装置、７…モニタ、９…装置本体、１１…超音波送信回路、１３…超音波受信回路、１５…Ｂモード処理回路、１７…ドプラ処理回路、１９…画像生成回路、２１…画像メモリ、２３…画像合成回路、２５…記憶回路、２７…インターフェース回路、２９…制御回路（制御プロセッサ（ＣＰＵ））、１１１…パルス発生器、１１３送信遅延回路、１１５…パルサ回路、１３１…プリアンプ、１３３…受信遅延回路、１３５…加算器、２９１…設定機能、２９３…推定機能、２９５…ゲイン逆補正機能、２９７…音場特性補正機能、２９９…解析機能。

Claims (16)

1. 被検体に対するスキャンにより得られたデータに基づいて、前記被検体内の構造物の位置を推定する推定部と、
   前記データに基づいて、前記推定部により推定された前記構造物の位置を除く前記被検体内の複数位置の組織性状を解析する解析部と、
   前記解析部による解析結果を表示する表示部と、
   を備えた医用診断装置。
2. 被検体に対するスキャンにより得られたデータに基づいて、前記被検体内の構造物の位置を推定する推定部と、
   前記データに基づいて、前記被検体内の組織性状を解析する解析部と、
   前記推定部により推定された前記構造物の位置を除く前記被検体内の複数位置について、前記解析部による解析結果を表示する表示部と、
   を備えた医用診断装置。
3. 前記推定部は、前記被検体に対する超音波スキャンにより得られた前記データに基づいて、前記構造物の位置を推定し、
   前記解析部は、前記被検体に対する前記超音波スキャンにより得られた前記データに基づいて、前記複数位置の組織性状を解析する、
   請求項１又は２に記載の医用診断装置。
4. 前記超音波スキャンにおける走査領域において、複数の領域を設定する設定部をさらに備え、
   前記推定部は、
   前記複数の領域各々において、前記データに基づいて分散値と平均値とを計算し、
   前記計算された分散値と平均値とに基づいて、前記走査領域における前記構造物の位置を推定する、
   請求項３に記載の医用診断装置。
5. 前記推定部は、前記超音波スキャンにおける走査領域に対してＣＦＡＲ（Ｃｏｎｔｒａｓｔ Ｆａｌｓｅ Ａｌａｒｍ Ｒａｔｅ）処理を実行することにより、前記走査領域における前記構造物の位置を推定する、
   請求項３に記載の医用診断装置。
6. 前記解析部は、前記被検体内を伝搬する超音波の減衰量に関する特徴量を、前記組織性状として計算し、
   前記表示部は、前記構造物の位置を除く前記被検体内の複数位置各々において、前記解析結果として前記減衰量を表示する、
   請求項３乃至５のうちいずれか一項に記載の医用診断装置。
7. 前記超音波スキャンに対応する超音波画像に対して、関心領域を入力する入力部をさらに備え、
   前記解析部は、前記関心領域内の複数の位置各々における前記減衰量に基づいて、前記関心領域内における前記減衰量を代表する代表値を計算し、
   前記表示部は、前記代表値を前記関心領域とともにさらに表示する、
   請求項６に記載の医用診断装置。
8. 前記解析部は、前記関心領域における画素値に基づいて、前記代表値に関する信頼度を計算し、
   前記表示部は、前記信頼度を前記代表値とともにさらに表示する、
   請求項７に記載の医用診断装置。
9. 前記データに対するゲイン補正をキャンセルするゲイン逆補正部と、
   前記超音波スキャンにおける音場特性に基づいて、前記データにおける前記音場特性の依存性をキャンセルする音場特性補正部と、
   をさらに備え、
   前記解析部は、前記ゲイン逆補正部と前記音場特性補正部とにより補正された補正データに基づいて、前記組織性状として前記減衰量を計算する、
   請求項６乃至８のうちいずれか一項に記載の医用診断装置。
10. 前記超音波スキャンにおいて、周波数が異なる複数の超音波を前記被検体に送信する送信部をさらに備え、
    前記解析部は、前記複数の超音波各々に対応する前記データを用いて、前記周波数の違いに起因する前記減衰量の差を、前記組織性状として計算する、
    請求項６乃至９のうちいずれか一項に記載の医用診断装置。
11. 前記超音波スキャンにおいて、Ｂモードに関する超音波送信における周波数帯域より狭帯域の超音波を送信する送信部をさらに備え、
    前記解析部は、前記狭帯域の超音波の受信により取得されたデータに対して周波数解析を実行し、前記周波数解析における周波数特性に基づいて、前記減衰量を計算する、
    請求項６乃至９のうちいずれか一項に記載の医用診断装置。
12. 前記解析部は、前記組織性状として、前記被検体内の組織の粘性と前記組織の弾性とのうち少なくとも一つを計算する、
    請求項１乃至１１のうちいずれか一項に記載の医用診断装置。
13. 前記表示部は、Ｂモード画像に前記減衰量を、前記減衰量の値に応じた色相で、かつ所定の不透明度または透明度で重畳して表示する、
    請求項６乃至１１のうちいずれか一項に記載の医用診断装置。
14. 前記推定部は、核磁気共鳴現象を前記スキャンとして利用したスキャンにより得られたデータに基づいて、前記構造物の位置を推定し、
    前記解析部は、前記核磁気共鳴現象を利用した前記スキャンにより得られた前記データに基づいて、前記複数位置の組織性状を解析する、
    請求項１又は２に記載の医用診断装置。
15. コンピュータに、
    被検体に対するスキャンにより得られたデータに基づいて、前記被検体内の構造物の位置を推定し、
    前記データに基づいて、前記推定された構造物の位置を除く前記被検体内の複数位置の組織性状を解析し、
    前記解析された組織性状に関する解析結果を、モニタに表示すること、
    を実現させる医用解析プログラム。
16. コンピュータに、
    被検体に対するスキャンにより得られたデータに基づいて、前記被検体内の構造物の位置を推定し、
    前記データに基づいて、前記被検体内の組織性状を解析し、
    前記推定された構造物の位置を除く前記被検体内の複数位置について、前記解析された組織性状に関する解析結果を、モニタに表示すること、
    を実現させる医用解析プログラム。