(57)【要約】【課題】三次元空間内に関心領域を容易に設定できるようにする。【解決手段】楕円体軸設定部３２は、ユーザがディスプレイ２２を見ながら入力する切断面位置情報および楕円体パラメータに基づいて楕円体の長軸と短軸を設定する。長軸端点検出部３４、正面図短軸端点検出部３６および側面図短軸端点検出部３８は、楕円体軸設定部３２において設定された長軸・短軸、および、エッジ抽出部２４から出力される各時相ごとの対象組織表面画像に基づいて長軸・短軸それぞれの端点を検出する。長軸・短軸設定部４０は、入力される長軸・短軸の端点の中から適切な長軸長・短軸長を選定して関心領域発生部４２に出力し、関心領域発生部４２において三次元の関心領域が生成される。【選択図】        図１
2004202132.tif 000002
【特許請求の範囲】
【請求項１】
対象組織を含む三次元空間に対して超音波を送受波して複数のボクセルデータから成るボリュームデータを形成するボリュームデータ形成手段と、
前記三次元空間内に設定される互いに交差関係にある複数のプレーンに基づいて、前記対象組織を貫く第一基準線を設定する基準線設定手段と、
前記第一基準線、および、それに交差する第二基準線に沿って、ボクセルデータを参照して、前記対象組織の境界を探索する境界探索手段と、
前記探索された境界に基づいて、前記対象組織に関する三次元の関心領域を特定するための複数の形状パラメータを演算する形状パラメータ演算手段と、
前記複数の形状パラメータに基づいて前記三次元空間内に前記関心領域を設定する関心領域設定手段と、
を有することを特徴とする超音波診断装置。
【請求項２】
請求項１記載の超音波診断装置であって、
前記対象組織の構造に基づいて決定される基準点を設定する基準点設定手段をさらに有し、
前記基準線設定手段は、前記基準点を通る前記第一基準線の傾きを設定する、ことを特徴とする超音波診断装置。
【請求項３】
請求項２記載の超音波診断装置であって、
前記ボリュームデータに基づいて、前記各プレーン上において前記対象組織を含むプレーン画像を形成するプレーン画像形成手段をさらに有し、
前記基準点設定手段は、ディスプレイに表示される前記各プレーン画像を利用してユーザが入力する基準点設定情報に基づいて前記基準点を設定する、
ことを特徴とする超音波診断装置。
【請求項４】
請求項３記載の超音波診断装置であって、
前記基準線設定手段は、ディスプレイに表示される前記各プレーン画像を利用してユーザが入力する基準線傾き情報であって、前記第一基準線の前記各プレーン画像への投影線の傾きを表す基準線傾き情報に基づいて前記第一基準線の傾きを設定する、
ことを特徴とする超音波診断装置。
【請求項５】
請求項２記載の超音波診断装置であって、
前記基準点設定手段は、前記基準点として、前記ボリュームデータに基づいて演算される前記対象組織の重心点を設定する、
ことを特徴とする超音波診断装置。
【請求項６】
請求項５記載の超音波診断装置であって、
前記基準線設定手段は、前記ボリュームデータを利用して得られる前記対象組織の境界情報に基づいて前記第一基準線の傾きを設定する、
ことを特徴とする超音波診断装置。
【請求項７】
請求項１から６いずれか１項記載の超音波診断装置であって、
前記複数のプレーンは、互いに直交関係にあることを特徴とする超音波診断装置。
【請求項８】
請求項１から７いずれか１項記載の超音波診断装置であって、
前記第二基準線は、前記第一基準線に直交することを特徴とする超音波診断装置。
【請求項９】
請求項２から７いずれか１項記載の超音波診断装置であって、
前記第二基準線は、前記第一基準線に直交し且つ前記基準点を通り複数設定されることを特徴とする超音波診断装置。
【請求項１０】
対象組織を含む三次元空間に対して超音波を送受波して複数のボクセルデータから成るボリュームデータを形成するボリュームデータ形成手段と、
前記三次元空間内において互いに直交関係にある三つのプレーンを設定するプレーン設定手段と、
前記対象組織を貫く第一基準線を設定する第一基準線設定手段であって、前記三つのプレーンのうちの二つのプレーンの各々に投影される前記第一基準線の投影線に基づいて、前記第一基準線を設定する第一基準線設定手段と、
前記第一基準線が投影された一方のプレーン内において、前記第一基準線に直交する第二基準線を設定する第二基準線設定手段と、
前記第一基準線が投影された他方のプレーン内において、前記第一基準線に直交する第三基準線を設定する第三基準線設定手段と、
前記ボクセルデータに基づいて、前記第一基準線と前記対象組織の境界との交点である第一基準線端点を演算し、前記第二基準線と前記対象組織の境界との交点である第二基準線端点を演算し、前記第三基準線と前記対象組織の境界との交点である第三基準線端点を演算する基準線端点演算手段と、
前記第一基準線端点に基づいて前記第一基準線上において所定の長さに決定される第一軸を設定し、前記第二基準線端点に基づいて前記第二基準線上において所定の長さに決定される第二軸を設定し、前記第三基準線端点に基づいて前記第三基準線上において所定の長さに決定される第三軸を設定する軸設定手段と、
前記設定された第一軸、第二軸および第三軸に基づいて、前記三次元空間内に楕円体形状の関心領域を設定する関心領域設定手段と、
を有することを特徴とする超音波診断装置。
【請求項１１】
請求項１０記載の超音波診断装置であって、
前記ボリュームデータに基づいて、前記各プレーン上において前記対象組織を含むプレーン画像を形成するプレーン画像形成手段と、
前記対象組織の構造に基づいて決定される基準点を設定する基準点設定手段であって、ディスプレイに表示される前記各プレーン画像を利用してユーザが入力する基準点設定情報に基づいて前記基準点を設定する基準点設定手段と、をさらに有し、
前記第一基準線設定手段は、ディスプレイに表示される前記各プレーン画像を利用してユーザが入力する前記投影線の傾き角度に基づいて前記基準点を通る前記第一基準線を設定し、
前記第二基準線および前記第三基準線はそれぞれ前記基準点を通る、
ことを特徴とする超音波診断装置。
【請求項１２】
請求項１１記載の超音波診断装置であって、
前記プレーン設定手段は、ユーザが入力する各プレーンの位置情報に基づいて前記三つのプレーンの位置を決定し、
前記プレーン画像形成手段は、前記プレーン画像として、前記各プレーン上における前記対象組織の切断面画像を形成する、
ことを特徴とする超音波診断装置。
【請求項１３】
請求項１２記載の超音波診断装置であって、
前記軸設定手段は、前記第一軸、第二軸および第三軸の各軸の長さとして、各軸が前記対象組織の境界を超えて外部に突き出る長さに設定し、
前記関心領域設定手段は、前記対象組織の境界内部を取り囲む楕円体形状の関心領域を設定する、
ことを特徴とする超音波診断装置。
【請求項１４】
請求項１３記載の超音波診断装置であって、
前記ボリュームデータおよび前記関心領域に基づいて、前記対象組織の境界内部が前記関心領域をはみ出したことを判定する関心領域判定手段をさらに有することを特徴とする超音波診断装置。
【請求項１５】
請求項１４記載の超音波診断装置であって、
前記関心領域判定手段は、前記関心領域の表面上に前記対象組織の境界内部のボクセルが含まれる場合、前記対象組織の境界内部が前記関心領域をはみ出したと判定することを特徴とする超音波診断装置。
【請求項１６】
請求項１５記載の超音波診断装置であって、
前記対象組織の境界内部が前記関心領域をはみ出した場合、前記第一軸、第二軸および第三軸を再設定する関心領域修正手段をさらに有することを特徴とする超音波診断装置。
【請求項１７】
請求項１６記載の超音波診断装置であって、
前記関心領域修正手段は、前記関心領域の表面上における前記対象組織の境界内部のはみ出し部分の位置に基づいて、前記第一軸、第二軸および第三軸の長さを再設定することを特徴とする超音波診断装置。
【請求項１８】
請求項１７記載の超音波診断装置であって、
前記関心領域の表面上における前記対象組織の境界内部のはみ出し部分に所定の表示処理を施したはみ出し画像を形成するはみ出し画像形成手段をさらに有することを特徴とする超音波診断装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、超音波診断装置に関し、特に関心領域を設定できる三次元超音波診断装置に関する。
【０００２】
【従来の技術】
二次元超音波診断装置では、関心領域を設定することで注目部位の診断機能を向上させていた。つまり、関心領域内においてカラードップラ画像を形成するためのドップラ演算を行うことで演算処理時間を短縮させ、また、対象組織の面積演算の際に関心領域内のみを演算対象とすることで演算精度を高めていた。関心領域は、三次元超音波診断装置においても利用され、演算処理時間の短縮や、面積あるいは体積演算の精度向上が図られている（例えば、特許文献１参照）。
【０００３】
【特許文献１】
特開２００１−３７７５６号公報
【０００４】
【発明が解決しようとする課題】
ところが、三次元超音波診断装置における関心領域の設定には、関心領域の三次元形状の設定や、三次元空間内における関心領域の位置の設定が必要であり、その設定操作は操作者にとって複雑で煩わしいものである。
【０００５】
そこで、本発明は、三次元空間内に関心領域を容易に設定できるようにすることを目的とする。
【０００６】
【課題を解決するための手段】
（１）上記目的を達成するために、本発明に係る超音波診断装置は、対象組織を含む三次元空間に対して超音波を送受波して複数のボクセルデータから成るボリュームデータを形成するボリュームデータ形成手段と、前記三次元空間内に設定される互いに交差関係にある複数のプレーンに基づいて、前記対象組織を貫く第一基準線を設定する基準線設定手段と、前記第一基準線、および、それに交差する第二基準線に沿って、ボクセルデータを参照して、前記対象組織の境界を探索する境界探索手段と、前記探索された境界に基づいて、前記対象組織に関する三次元の関心領域を特定するための複数の形状パラメータを演算する形状パラメータ演算手段と、前記複数の形状パラメータに基づいて前記三次元空間内に前記関心領域を設定する関心領域設定手段と、を有するものとする。
【０００７】
上記構成において、対象組織の境界とは診断を必要とする部位を規定するものであり、例えば、対象組織の表面部に相当する。
【０００８】
上記構成によれば、超音波診断装置が形状パラメータを演算してこの形状パラメータに基づいた関心領域を設定するため、ユーザによる関心領域の設定操作が簡略化され、三次元空間内への関心領域の設定が容易になる。
【０００９】
望ましくは、前記対象組織の構造に基づいて決定される基準点を設定する基準点設定手段をさらに有し、前記基準線設定手段は、前記基準点を通る前記第一基準線の傾きを設定するものとする。
【００１０】
望ましくは、前記ボリュームデータに基づいて、前記各プレーン上において前記対象組織を含むプレーン画像を形成するプレーン画像形成手段をさらに有し、前記基準点設定手段は、ディスプレイに表示される前記各プレーン画像を利用してユーザが入力する基準点設定情報に基づいて前記基準点を設定するものとする。
【００１１】
上記構成によれば、ユーザは対象組織に関する基準点の位置を視覚的に確認しながら基準点情報を設定することができるため、基準点の設定が容易になり、延いては関心領域の設定が容易になる。
【００１２】
望ましくは、前記基準線設定手段は、ディスプレイに表示される前記各プレーン画像を利用してユーザが入力する基準線傾き情報であって、前記第一基準線の前記各プレーン画像への投影線の傾きを表す基準線傾き情報に基づいて前記第一基準線の傾きを設定するものとする。
【００１３】
上記構成によれば、ユーザはプレーン画像を利用して第一基準線を視覚的に確認しながらその傾きを設定することができるため、第一基準線の設定が容易になり、延いては関心領域の設定が容易になる。
【００１４】
望ましくは、前記基準点設定手段は、前記基準点として、前記ボリュームデータに基づいて演算される前記対象組織の重心点を設定するものとする。
【００１５】
望ましくは、前記基準線設定手段は、前記ボリュームデータを利用して得られる前記対象組織の境界情報に基づいて前記第一基準線の傾きを設定するものとする。
【００１６】
上記構成において、第一基準線の傾きは、例えば、第一基準線が対象組織の境界内部を最も長く貫くように設定されてもよく、また、対象組織境界の特定部位を通るように設定されてもよい。
【００１７】
上記構成によれば、対象組織の構造に基づいて第一基準線を設定することが可能になる。つまり、対象組織の構造を反映した関心領域の設定が可能になる。
【００１８】
望ましくは、前記複数のプレーンは、互いに直交関係にあるものとする。
【００１９】
上記構成によれば、複数のプレーンが互いに直交関係にあるため、各プレーンの位置関係の把握が容易になる。
【００２０】
望ましくは、前記第二基準線は、前記第一基準線に直交するものとする。
【００２１】
望ましくは、前記第二基準線は、前記第一基準線に直交し且つ前記基準点を通り複数設定されるものとする。
【００２２】
（２）また、上記目的を達成するために、本発明に係る超音波診断装置は、対象組織を含む三次元空間に対して超音波を送受波して複数のボクセルデータから成るボリュームデータを形成するボリュームデータ形成手段と、前記三次元空間内において互いに直交関係にある三つのプレーンを設定するプレーン設定手段と、前記対象組織を貫く第一基準線を設定する第一基準線設定手段であって、前記三つのプレーンのうちの二つのプレーンの各々に投影される前記第一基準線の投影線に基づいて、前記第一基準線を設定する第一基準線設定手段と、前記第一基準線が投影された一方のプレーン内において、前記第一基準線に直交する第二基準線を設定する第二基準線設定手段と、前記第一基準線が投影された他方のプレーン内において、前記第一基準線に直交する第三基準線を設定する第三基準線設定手段と、前記ボクセルデータに基づいて、前記第一基準線と前記対象組織の境界との交点である第一基準線端点を演算し、前記第二基準線と前記対象組織の境界との交点である第二基準線端点を演算し、前記第三基準線と前記対象組織の境界との交点である第三基準線端点を演算する基準線端点演算手段と、前記第一基準線端点に基づいて前記第一基準線上において所定の長さに決定される第一軸を設定し、前記第二基準線端点に基づいて前記第二基準線上において所定の長さに決定される第二軸を設定し、前記第三基準線端点に基づいて前記第三基準線上において所定の長さに決定される第三軸を設定する軸設定手段と、前記設定された第一軸、第二軸および第三軸に基づいて、前記三次元空間内に楕円体形状の関心領域を設定する関心領域設定手段と、を有するものとする。
【００２３】
上記構成において、三つのプレーンは任意の位置に設定可能であり、例えば、対象組織に関する上面、正面および側面に設定される。また、第一基準線端点、第二基準線端点および第三基準線端点は、対象組織の境界の形状に応じて、それぞれ複数であってもよい。対象組織としては、例えば、心臓左室内の心腔が挙げられる。この場合、左室内壁が対象組織の境界に相当する。
【００２４】
上記構成によれば、第一軸、第二軸および第三軸に基づいて楕円体形状およびその位置が決定されるため、三次元空間内に楕円体形状の関心領域が容易に設定できる。
【００２５】
望ましくは、前記ボリュームデータに基づいて、前記各プレーン上において前記対象組織を含むプレーン画像を形成するプレーン画像形成手段と、前記対象組織の構造に基づいて決定される基準点を設定する基準点設定手段であって、ディスプレイに表示される前記各プレーン画像を利用してユーザが入力する基準点設定情報に基づいて前記基準点を設定する基準点設定手段と、をさらに有し、前記第一基準線設定手段は、ディスプレイに表示される前記各プレーン画像を利用してユーザが入力する前記投影線の傾き角度に基づいて前記基準点を通る前記第一基準線を設定し、前記第二基準線および前記第三基準線はそれぞれ前記基準点を通るものとする。
【００２６】
上記構成により、ユーザは、所望の位置に基準点および第一基準線を設定することができ、また、これら設定情報に基づいて超音波診断装置が楕円体形状の関心領域を設定するため、ユーザ所望の位置に、且つ、容易に関心領域を設定することができる。
【００２７】
望ましくは、前記プレーン設定手段は、ユーザが入力する各プレーンの位置情報に基づいて前記三つのプレーンの位置を決定し、前記プレーン画像形成手段は、前記プレーン画像として、前記各プレーン上における前記対象組織の切断面画像を形成するものとする。
【００２８】
上記構成によれば、ユーザは所望の位置に三つの切断面を設定し、この切断面を利用して第一基準線を設定することができる。
【００２９】
望ましくは、前記軸設定手段は、前記第一軸、第二軸および第三軸の各軸の長さとして、各軸が前記対象組織の境界を超えて外部に突き出る長さに設定し、前記関心領域設定手段は、前記対象組織の境界内部を取り囲む楕円体形状の関心領域を設定するものとする。
【００３０】
上記構成によれば、関心領域となる楕円体形状が、診断対象となる対象組織の境界内部よりも大きく設定される。このため、診断対象となる部位のほぼ全部、望ましくは診断対象となる部位の全てを取り囲む関心領域の設定が可能になる。
【００３１】
望ましくは、前記ボリュームデータおよび前記関心領域に基づいて、前記対象組織の境界内部が前記関心領域をはみ出したことを判定する関心領域判定手段をさらに有するものとする。
【００３２】
望ましくは、前記関心領域判定手段は、前記関心領域の表面上に前記対象組織の境界内部のボクセルが含まれる場合、前記対象組織の境界内部が前記関心領域をはみ出したと判定するものとする。
【００３３】
望ましくは、前記対象組織の境界内部が前記関心領域をはみ出した場合、前記第一軸、第二軸および第三軸を再設定する関心領域修正手段をさらに有するものとする。
【００３４】
上記構成によれば、関心領域が修正されるため、関心領域の設定精度が高まり、延いては超音波診断装置による診断精度が高まる。
【００３５】
望ましくは、前記関心領域修正手段は、前記関心領域の表面上における前記対象組織の境界内部のはみ出し部分の位置に基づいて、前記第一軸、第二軸および第三軸の長さを再設定するものとする。
【００３６】
望ましくは、前記関心領域の表面上における前記対象組織の境界内部のはみ出し部分に所定の表示処理を施したはみ出し画像を形成するはみ出し画像形成手段をさらに有するものとする。
【００３７】
上記構成によれば、ユーザはディスプレイに表示されるはみ出し画像を見ることで、診断対象となる部位のはみ出し部分を容易に確認することができる。したがって、はみ出し部分の診断への影響を推し量ることが可能になり、例えば、診断への影響が少ないと判断した場合、関心領域の再設定を行うことなく診断を続行することが可能になる。
【００３８】
【発明の実施の形態】
以下、本発明の好適な実施の形態を図面に基づいて説明する。
【００３９】
図１には、本発明に係る超音波診断装置の好適な実施形態が示されており、図１はその全体構成を示すブロック図である。
【００４０】
送受信部１２は、プローブ１０を介して対象組織を含む空間内に超音波を送受波することで各時相ごとの三次元超音波画像を、三次元データメモリ１４へ出力する。各時相ごとに得られた三次元超音波画像は、座標変換部１６において表示座標系に変換された後、各時相ごとに表示画像形成部２０およびエッジ抽出部２４に出力される。
【００４１】
表示画像形成部２０は、座標変換部１６から出力される座標変換された三次元超音波画像に基づいて、この三次元超音波画像を二次元表面上に投影した投影画像を形成する。投影画像の形成は、例えばボリュームレンダリング法を用いて行われる。また、表示画像形成部２０は、三次元超音波画像内における対象組織に関する切断面画像を形成する。切断面画像としては、対象組織の上面からみた切断面（上面図）、正面からみた切断面（正面図）および側面からみた切断面（側面図）の三つの切断面画像が形成される。各切断面の位置はユーザからの切断面位置情報に基づいて設定される。表示画像形成部２０で形成された投影画像と三つの切断面画像はディスプレイ２２に表示される（ディスプレイ２２に表示される画像は、例えば、図２から図５に示すものである）。
【００４２】
エッジ抽出部２４は、座標変換部１６から出力される座標変換された三次元超音波画像に基づいて、各時相ごとに対象組織の表面画像を抽出する。三次元超音波画像は複数のボクセルデータが規則的に配列されたものである。エッジ抽出部２４は、三次元超音波画像に対して二値化処理を施すことで、複数のボクセルを、対象組織に対応するボクセル（対象組織ボクセル）と対象組織以外の部位に対応するボクセル（非対象組織ボクセル）とに分別する。そして、三次元超音波画像を構成する全てのボクセルの中で、隣接するボクセルとして少なくとも一つの非対象組織を有する対象組織ボクセルの集合を抽出し、抽出したボクセル集合による対象組織表面画像を抽出する。
【００４３】
関心領域発生ブロック３０は、エッジ抽出部２４から出力される対象組織表面画像、並びに、ユーザから入力される切断面位置情報および楕円体パラメータに基づいて、三次元超音波画像内に楕円体形状の関心領域を形成する。ここで、関心領域発生ブロック３０内各部の動作について図２から図５を利用して説明する。なお、以下の説明において図１の記載部分については図１における符号を付すものとする。
【００４４】
図２から図５は、楕円体形状の三次元関心領域の形成手法を説明するための図であり、図２は楕円体形状の三次元関心領域７０を含む三次元超音波画像の斜視図であり、三次元関心領域７０の内部には図示しない対象組織画像が存在する。図３は三次元関心領域７０を含む三次元超音波画像の上面図、図４は三次元関心領域７０を含む三次元超音波画像の正面図、図５は三次元関心領域７０を含む三次元超音波画像の側面図をそれぞれ示している。ここで上面図とは対象組織の上面からみた切断面を、正面図とは対象組織の正面からみた切断面を、側面図とは対象組織の側面からみた切断面を、それぞれ示している。
【００４５】
楕円体軸設定部３２は、ユーザから入力される切断面位置情報および楕円体パラメータに基づいて楕円体の長軸と短軸を設定する。長軸と短軸の設定にあたり、ユーザは、図３から図５に示す各切断面の位置を決定する。例えば、図３の上面図の位置は、ユーザが切断面（図２に示すｙｚ平面）をｘ軸方向の所望の位置に移動することで設定される。同様に図４の正面図および図５の側面図についても、ユーザが各面をその面に垂直な座標軸方向に移動させることで所望の位置に設定される。ユーザが設定する切断面位置情報は、表示画像形成部２０にも入力されており、表示画像形成部２０は切断面位置情報に基づいて座標変換部１６から出力される三次元超音波画像の各切断面における切断面画像を形成する。このように、図３から図５に示す上面図、正面図および側面図に対応する各画像は表示画像形成部２０において画像形成されてディスプレイ２２に表示されるため、ユーザはディスプレイ２２に表示される切断面画像を見ながら、各切断面位置を設定することができる。
【００４６】
各切断面画像は、三次元超音波画像内に対象組織を取り囲む楕円体形状の関心領域を形成するために利用される。このため、各切断面は切断面内の対象組織画像がなるべく大きくなるように、望ましくは最大になるように、その切断位置が設定される。つまり、ユーザはディスプレイ２２に表示される各切断面における対象組織の切断画像を見ながら、各切断面における対象組織の切断画像がなるべく大きくなるような切断面位置を設定する。
【００４７】
以上のようにして各切断面が設定された後、楕円体軸設定部３２はユーザから入力される楕円体パラメータ情報（Ｘ0,Ｙ0,Ｚ0,θ, φ）に基づいて、楕円体の長軸と短軸を設定する。楕円体の中心点位置（Ｘ0,Ｙ0,Ｚ0）は、各切断面において対象組織断面図のなるべく中心部分に位置するように設定されるものであり、ユーザが三つの切断面を見ながら設定する。例えば、上面図に基づいて対象組織の中心点になるようにＹ0,Ｚ0を設定した後、正面図と側面図を参照しながらＸ0を設定する。もちろん正面図に基づいてＸ0,Ｙ0を設定した後に上面図および側面図を参照してＺ0を設定してもよく、あるいは側面図に基づいてＺ0, Ｘ0を設定した後に上面図および正面図を参照してＹ0を設定してもよい。このようにして、楕円体の中心点位置、すなわち、座標（Ｘ0,Ｙ0,Ｚ0）に対応した三次元空間内における一つの点が設定される。
【００４８】
次に、ユーザは設定した楕円体の中心点位置に基づいて、楕円体長軸の傾きであるθ, φを設定する。傾きθは正面図におけるｘ軸に対する長軸の傾きである。ユーザは正面図の画像を見ながら正面図内において、楕円体の中心点を通り対象組織内でなるべく長い直線、望ましくは最長の直線を求めてその傾きθを設定する。同様にしてユーザは側面図の画像を見ながら側面図内において、楕円体の中心点を通り対象組織内でなるべく長い直線、望ましくは最長の直線を求めてその傾きφを設定する。このようにして、楕円体の長軸、すなわち、座標（Ｘ0,Ｙ0,Ｚ0）および傾きθ,φに対応した三次元空間内における一本の直線が設定される。
【００４９】
そして、楕円体軸設定部３２は、設定された楕円体の中心点位置および長軸に基づいて楕円体の短軸を設定する。短軸は、三次元空間内において、楕円体中心点を通り長軸に直交する平面内において二本設定される。つまり、正面図内において長軸に直交する直線（正面図短軸）と、側面図内において長軸に直交する直線（側面図短軸）の二本が設定される。
【００５０】
以上のようにして、楕円体軸設定部３２において、三次元関心領域である楕円体の長軸、正面図短軸および側面図短軸が設定される。
【００５１】
図１に戻り、長軸端点検出部３４は、楕円体軸設定部３２において設定された長軸、および、エッジ抽出部２４から出力される各時相ごとの対象組織表面画像に基づいて長軸の端点を検出する。つまり、対象組織表面画像と長軸との複数の交点座標を各時相ごとに算出する。例えば、対象組織が楕円体に近い形状の場合、対象組織表面画像と長軸との交点は二点となる。各時相ごとに算出された複数の交点座標は長軸・短軸設定部４０に出力される。
【００５２】
正面図短軸端点検出部３６は、楕円体軸設定部３２において設定された正面図短軸、および、エッジ抽出部２４から出力される各時相ごとの対象組織表面画像に基づいて正面図短軸の端点を検出する。つまり、対象組織表面画像と正面図短軸との複数の交点座標を各時相ごとに算出する。例えば、対象組織が楕円体に近い形状の場合、対象組織表面画像と正面図短軸との交点は二点となる。各時相ごとに算出された複数の交点座標は長軸・短軸設定部４０に出力される。
【００５３】
側面図短軸端点検出部３８は、楕円体軸設定部３２において設定された側面図短軸、および、エッジ抽出部２４から出力される各時相ごとの対象組織表面画像に基づいて側面図短軸の端点を検出する。つまり、対象組織表面画像と側面図短軸との複数の交点座標を各時相ごとに算出する。例えば、対象組織が楕円体に近い形状の場合、対象組織表面画像と側面図短軸との交点は二点となる。各時相ごとに算出された複数の交点座標は長軸・短軸設定部４０に出力される。
【００５４】
長軸・短軸設定部４０は、長軸端点検出部３４、正面図短軸端点検出部３６および側面図短軸端点検出部３８から出力される交点座標に基づいて長軸長および短軸長を設定する。つまり、長軸端点検出部３４から各時相ごと出力される複数の交点の中から、楕円体中心点からの距離が最も大きい交点（長軸交点）を各時相ごとに選定する。そして、各時相ごとに選定した長軸交点の中で、楕円体中心点からの距離が最も大きい交点を選定し、選定した長軸交点と楕円体中心点との距離を長軸長とする。
【００５５】
また、長軸・短軸設定部４０は、正面図短軸端点検出部３６から各時相ごと出力される複数の交点の中から、楕円体中心点からの距離が最も大きい交点（正面図短軸交点）を各時相ごとに選定する。そして、各時相ごとに選定した正面図短軸交点の中で、楕円体中心点からの距離が最も大きい交点を選定し、選定した正面図短軸交点と楕円体中心点との距離を正面図短軸長とする。
【００５６】
さらに、長軸・短軸設定部４０は、側面図短軸端点検出部３８から各時相ごと出力される複数の交点の中から、楕円体中心点からの距離が最も大きい交点（側面図短軸交点）を各時相ごとに選定する。そして、各時相ごとに選定した側面図短軸交点の中で、楕円体中心点からの距離が最も大きい交点を選定し、選定した側面図短軸交点と楕円体中心点との距離を側面図短軸長とする。
【００５７】
長軸・短軸設定部４０は、複数の時相の中から最大の軸長を選定するが、例えば、対象組織が心臓心室の場合、心臓心室は心拍運動により収縮拡張を繰り返しているため、拡張末期時点において長軸・短軸を設定することが望ましい。拡張末期時点は心電波形のＲ波によって確認されるため、Ｒ波発生時点の時相において長軸長、正面図短軸長および側面図短軸長が選定される。
【００５８】
関心領域発生部４２は、長軸・短軸設定部４０から出力される長軸長、正面図短軸長および側面図短軸長に基づいて楕円体形状の関心領域を発生する。関心領域発生部４２で発生する楕円体の一般式は次式で表現される。
【数１】
2004202132.tif 000003
数１は、中心点が原点にあり、長軸半径をａ´、二つの短軸半径をｂ´,ｃ´とする楕円体領域の一般式である。つまり、ａ´,ｂ´,ｃ´を設定することで、数１に基づいて三次元空間内における原点を中心点とする楕円体形状の三次元関心領域が決定される。つまり、長軸・短軸設定部４０から出力される長軸長、正面図短軸長および側面図短軸長を、それぞれ数１のａ´,ｂ´,ｃ´に代入することで、楕円体形状の中心点が原点にある対象組織の外縁が設定できる。ただし、関心領域としては対象組織よりもやや広めに設定されることが望ましいため、例えば、長軸長、正面図短軸長および側面図短軸長よりも１〜２割り程度大きい値をａ´,ｂ´,ｃ´に代入してもよい。数１で表現される楕円体は中心点が原点に位置する。そこで数１に基づいて、この楕円体をＺ軸回りにθ、Ｙ軸周りにφ回転させ、さらに中心点を位置（Ｘ0,Ｙ0,Ｚ0）に平行移動させる一般式を次のように作成する。
【００５９】
点Ｐ（ｘ,ｙ,ｚ）がＺ軸回りにθ回転する場合の変換は次の行列で表現できる。
【数２】
2004202132.tif 000004
さらに、点Ｐ（ｘ,ｙ,ｚ）がＹ軸回りにφ回転する場合の変換は次の行列で表現できる。
【数３】
2004202132.tif 000005
したがって、点Ｐ（ｘ,ｙ,ｚ）をＺ軸回りにθ、Ｙ軸回りにφ回転させた上で位置（Ｘ0,Ｙ0,Ｚ0）に平行移動させる場合の変換は次の行列で表現できる。
【数４】
2004202132.tif 000006
数４の変換式を数１に適用することで、数１で表現される楕円体を、Ｚ軸回りにθ、Ｙ軸回りにφ回転させ、且つ、中心点を位置（Ｘ0,Ｙ0,Ｚ0）に平行移動させることができる。関心領域発生部４２は、長軸・短軸設定部４０から出力される長軸長（ａ´）、正面図短軸長（ｂ´）および側面図短軸長（ｃ´）に基づいて、また、ユーザにより設定された楕円体パラメータ情報（Ｘ0,Ｙ0,Ｚ0,θ, φ）に基づいて、数１および数４を利用して、三次元超音波画像内に楕円体形状の三次元関心領域を発生する。
【００６０】
以上のようにして、関心領域発生ブロック３０から出力される楕円体形状の三次元関心領域は、関心領域エッジ抽出部５０に出力され、ここで関心領域の表面画像が抽出される。関心領域の表面は、数１における不等号を等号に置き換えた式に、数４の変換式を適用することで得られる。
【００６１】
関心領域適合判定部５２は、関心領域エッジ抽出部５０から出力される関心領域の表面画像、および、座標変換部１６から出力される三次元超音波画像が二値化処理部５４により二値化処理された二値化画像に基づいて、関心領域内に対象組織が収まっているか否かを判定する。二値化処理部５４から出力される二値化画像は、対象組織に対応するボクセル（対象組織ボクセル）と対象組織以外の部位に対応するボクセル（非対象組織ボクセル）とに分別された画像である。そこで、関心領域適合判定部５２は、関心領域の表面画像と二値化画像との比較を行い、関心領域の表面画像上に対象組織ボクセルが存在する場合、対象組織画像が関心領域をはみ出したものと判定する。関心領域適合判定部５２は、関心領域の表面画像上における対象組織ボクセル部分、つまり、はみ出し表面領域に例えば色付け処理を施して対象組織が関心領域をはみ出していることを明示した画像を形成し、表示画像形成部２０に出力する。また、関心領域適合判定部５２は、はみ出し表面領域の画像を関心領域修正手段６０へも出力する。
【００６２】
関心領域修正手段６０は、対象組織が関心領域をはみ出した場合に、関心領域を修正する手段である。次に、図６を利用して関心領域修正手段６０について説明する。
【００６３】
図６は、関心領域修正手段６０の内部構成を示すブロック図である。関心領域修正手段６０は、重心検出部７２、長軸成分検出部７４、正面図短軸成分検出部７６、側面図短軸成分検出部７８および長軸短軸拡大率設定部８０で構成される。重心検出部７２には、関心領域適合判定部５２から出力される、はみ出し表面領域の画像が入力される。はみ出し表面領域の画像は各時相ごとに形成されており、重心検出部７２は各時相ごとのはみ出し表面領域を所定時間内で積算する。つまり、所定時間内において一時相でもはみ出しを経験した表面部分を全て抽出する。そして、積算した表面領域の重心点を算出する。なお、積算を行う所定時間としては、例えば対象組織が心臓左室の場合、一心拍が好適である。
【００６４】
長軸成分検出部７４は、重心検出部７２が出力する重心点の位置に基づいて、長軸成分を算出する。長軸成分の算出手法を図７を利用して説明する。長軸成分検出部７４は、重心検出部７２が出力するはみ出し表面領域の重心点から関心領域の長軸へ垂線を引いてこの垂線と長軸との交点を抽出し、抽出した交点と重心点側に位置する長軸端点との距離ｇ１を算出し、長軸短軸拡大率設定部８０へ出力する。
【００６５】
図６に戻り、正面図短軸成分検出部７６は、重心検出部７２が出力する重心点の位置に基づいて、この重心点から関心領域の正面図短軸へ垂線を引いてこの垂線と正面図短軸との交点を抽出し、抽出した交点と重心点側に位置する正面図短軸端点との距離ｇ２を算出し、長軸短軸拡大率設定部８０へ出力する。また、側面図短軸成分検出部７８は、重心検出部７２が出力する重心点の位置に基づいて、この重心点から関心領域の側面図短軸へ垂線を引いてこの垂線と側面図短軸との交点を抽出し、抽出した交点と重心点側に位置する側面図短軸端点との距離ｇ３を算出し、長軸短軸拡大率設定部８０へ出力する。
【００６６】
長軸短軸拡大率設定部８０は、長軸成分検出部７４、正面図短軸成分検出部７６および側面図短軸成分検出部７８からそれぞれ出力されるｇ１，ｇ２，ｇ３に基づいて、長軸および二つの短軸の拡大率を決定する。拡大率はｇ１，ｇ２，ｇ３のそれぞれの大きさに応じて予め設定されているものとする。例えば、ｇ１が長軸半径ａの１／４以下の場合には、はみ出し表面領域が長軸端点近くに存在すると判断して、長軸半径を２割大きくする。ｇ１が長軸半径ａの１／４より大きく１／２以下の場合には、はみ出し表面領域が長軸端点にそれほど近くないものと判断して、長軸半径を１割大きくする。同様に、ｇ２，ｇ３の大きさに基づいて、正面図短軸半径、側面図短軸半径の拡大率が決定される。長軸短軸拡大率設定部８０で設定した各軸の拡大率は長軸・短軸設定部（図１の符号４０）に出力され、各軸が拡大される。
【００６７】
図１に示した長軸・短軸設定部４０は、拡大した長軸長、正面図短軸長および側面図短軸長を関心領域発生部４２に出力し、関心領域発生部４２は拡大した関心領域を生成して関心領域適合判定部５２に出力する。このようにして、長軸・短軸を繰り返し拡大修正することで、対象組織がはみ出さない三次元関心領域が形成される。なお、関心領域修正手段６０にはユーザからの指示に基づいて、自動修正を停止させる機能があると好適である。これにより、対象組織のある部位が関心領域をはみ出しても、その部位が診断に差し支えないとユーザが判断して関心領域を固定することができる。
【００６８】
関心領域発生部４２で生成された三次元関心領域は、体積計測部６２にも出力される。体積計測部６２は、三次元関心領域内において対象組織の各時相における体積を計測するものであり、二値化処理部５４から出力される二値化画像に基づいて、三次元関心領域内における対象組織ボクセル数をカウントして対象組織の体積を計測する。各時相ごとに得られた対象組織の体積の計測結果はグラフ形成部６４に出力され、各時相ごとの対象組織の体積、つまり計測時刻と体積値を示したグラフが生成され表示画像形成部２０を介してディスプレイ２２に出力される。なお、グラフ形成部６４は、各時相ごとの体積計測結果に基づいて、例えば、体積変化率を演算してグラフ化してもよい。
【００６９】
前述したように、関心領域適合判定部５２は、はみ出し表面領域に色付け処理を施して、対象組織が関心領域をはみ出していることを明示した画像を形成し、表示画像形成部２０に出力する。そして、表示画像形成部２０は、座標変換部１６から出力される三次元超音波画像と、はみ出し表面領域に色付け処理を施した関心領域表面を合成して、ディスプレイ２２に表示する。
【００７０】
図８は、図１の超音波診断装置のディスプレイ２２に表示される対象組織画像を示す図である。対象組織画像８２は、三次元関心領域８４をはみ出した部分８６に色付け処理が施されて表示されている。ユーザはこの画像を見ることで、対象組織のどの部分が関心領域からはみ出しているのかを知ることができる。
【００７１】
図１に示した実施形態において、座標変換部１６の直後に並進回転移動キャンセル処理部が挿入されるとさらに好適である。
【００７２】
図９は、並進回転移動キャンセル処理部８７の内部構成を示すブロック図である。並進回転移動キャンセル処理部８７には座標変換部１６から出力される三次元超音波画像が入力され、二値化処理部８８において、対象組織に対応するボクセル（対象組織ボクセル）と対象組織以外の部位に対応するボクセル（非対象組織ボクセル）とに分別された二値化画像が形成される。
【００７３】
心室用ＲＯＩ（関心領域）発生器９０は、対象組織である心臓心室の外縁を成すＲＯＩの座標を発生する。心室用ＲＯＩは、例えば楕円形状であり、ユーザは超音波画像を見ながら楕円の長軸や短軸の長さ、中心点の位置、楕円の傾きなどの初期値を、ＲＯＩの中に心室の画像が収まるように設定する。この際、ユーザは超音波画像を見ながら一心拍分の運動を観察した上で、全ての時相においてＲＯＩが心臓左室を含むように、トラックボールなどを操作して初期値を決定する。ＲＯＩの設定はユーザによるマニュアル設定に限られるものではなく、心室の動きに応じて装置設定されるものでもよい。
【００７４】
心室用ゲート回路９２は、心室用ＲＯＩ内のボクセルデータのみを通過させる回路である。つまり、心室用ゲート回路９２の一方の入力端子には心室用ＲＯＩ発生器９０から出力されるＲＯＩの座標が入力され、他方の入力端子に入力される二値化画像において心室用ＲＯＩ内に属する座標のボクセルデータのみを抽出し、心腔抽出部９４に出力する。心腔抽出部９４は、ＲＯＩ内の二値化画像から心室内部の心腔画像を抽出する。心室重心演算部９６は、心腔抽出部９４から出力される心室内部画像における重心点の座標を各時相毎に算出する。算出された心室重心点の座標は読み出しアドレス発生器１１２および心室重心点メモリ９８に出力される。
【００７５】
弁輪部用ＲＯＩ（関心領域）発生器１００は、心室端に位置する弁輪部の外縁を成すＲＯＩの座標を発生する。弁輪部用ＲＯＩは、例えば楕円形状であり、ユーザは超音波画像を見ながら楕円の長軸や短軸の長さ、中心点の位置、楕円の傾きなどの初期値を、ＲＯＩの中に弁輪部の画像が収まるように設定する。この際、ユーザは超音波画像を見ながら一心拍分の運動を観察した上で、全ての時相においてＲＯＩが弁輪部を含むように、トラックボールなどを操作して初期値を決定する。ＲＯＩの設定はユーザによるマニュアル設定に限られるものではなく、弁輪部の動きに応じて装置設定されるものでもよい。
【００７６】
弁輪部ゲート回路１０２は、弁輪部用ＲＯＩ内のボクセルデータのみを通過させる回路である。つまり、弁輪部用ゲート回路１０２の一方の入力端子には弁輪部用ＲＯＩ発生器１００から出力されるＲＯＩの座標が入力され、他方の入力端子に入力される二値化画像において弁輪部用ＲＯＩ内に属する座標のボクセルデータのみを抽出し、弁輪部抽出部１０４に出力する。弁輪部抽出部１０４は、ＲＯＩ内の二値化画像から弁輪部画像を抽出する。弁輪部重心演算部１０６は、弁輪部抽出部１０４から出力される弁輪部画像に対して弁輪部の重心点の座標を各時相毎に算出する。算出された弁輪部重心点の座標は読み出しアドレス発生器１１２よび弁輪部重心点メモリ１０８に出力される。
【００７７】
心室重心点メモリ９８には心室の拡張末期時の心室重心点の座標が記憶される。拡張末期を知らせるトリガとしては心電波形のＲ波を利用する。つまり、拡張末期時に得られるＲ波をトリガとして、心室重心演算部９６から出力される心室重心点の座標を拡張末期時の心室重心点の座標として記憶しておく。同様にＲ波をトリガとして拡張末期時の弁輪部重心点の座標が、弁輪部重心演算部１０６から弁輪部重心点メモリ１０８に記憶される。
【００７８】
読み出し制御部１１０は、読み出しアドレス発生器１１２及びメモリ制御部１１４で構成されており、各時相のボリューム間における心室の並進移動量及び回転移動量をキャンセルした超音波画像を形成すべく、ボリュームメモリ１１６からボクセルデータを読み出す。つまり、読み出しアドレス発生器１１２は、心室重心点メモリ９８から拡張末期時の心室重心点の座標を取得し、また、弁輪部重心点メモリ１０８から拡張末期時の弁輪部重心点の座標を取得する。さらに、心室重心演算部９６から現ボリュームにおける心室重心点の座標を取得し、また、弁輪部重心演算部１０６から現ボリュームにおける弁輪部重心点の座標を取得する。
【００７９】
読み出しアドレス発生器１１２は、現ボリュームの心室重心点が拡張末期時の心室重心点に重なるように、かつ、現ボリュームにおける心室重心点と弁輪部重心点を通る直線が、拡張末期時における心室重心点と弁輪部重心点を通る直線に重なるような読み出しアドレスを演算する。
【００８０】
ボリュームメモリ１１６には二値化処理部８８から出力されたボクセルデータが原画像のアドレスのままボリューム毎にコピーされており、メモリ制御部１１４は、読み出しアドレス発生器１１２が算出した読み出しアドレスに従ってボリュームメモリ１１６からボクセルデータを読み出して、表示画像形成部（図１の符号２０）およびエッジ抽出部（図１の符号２４）に出力する。この結果、ボリュームメモリ１１６から出力されるボクセルデータに基づいた超音波画像は、対象組織である心臓左室の身体全体に対する相対的な並進移動および回転移動がキャンセルされた超音波画像となる。
【００８１】
【発明の効果】
以上説明したように、本発明に係る超音波診断装置により三次元空間内に関心領域を容易に設定することができる。
【図面の簡単な説明】
【図１】本発明に係る超音波診断装置の好適な実施形態を示すブロック図である。
【図２】楕円体形状の三次元関心領域を含む三次元超音波画像の斜視図である。
【図３】三次元関心領域を含む三次元超音波画像の上面図である。
【図４】三次元関心領域を含む三次元超音波画像の正面図である。
【図５】三次元関心領域を含む三次元超音波画像の側面図である。
【図６】関心領域修正手段の内部構成を示すブロック図である。
【図７】長軸成分の算出手法を説明するための図である。
【図８】ディスプレイに表示される対象組織画像を示す図である。
【図９】並進回転移動キャンセル処理部の内部構成を示すブロック図である。
【符号の説明】
２４  エッジ抽出部、３２  楕円体軸設定部、３４  長軸端点検出部、３６  正面図短軸端点検出部、３８  側面図短軸端点検出部、４０  長軸・短軸設定部、４２  関心領域発生部、５０  関心領域エッジ抽出部、５２  関心領域適合判定部、６０  関心領域修正手段。
【図１】
2004202132.tif 000007
【図２】
2004202132.tif 000008
【図３】
2004202132.tif 000009
【図４】
2004202132.tif 000010
【図５】
2004202132.tif 000011
【図６】
2004202132.tif 000012
【図７】
2004202132.tif 000013
【図８】
2004202132.tif 000014
【図９】
2004202132.tif 000015