(57)【特許請求の範囲】
【請求項１】
超音波を被検体に対して送信し、前記被検体から反射される超音波に基づく受信データに基づいて前記被検体における複数の位置での組織の運動情報を検出する運動情報演算手段と、
所定心拍の所定の時間範囲内における前記運動情報の経時変化を検出し、前記所定の時間範囲内にある最初の極点を抽出し、前記最初の極点に到達するまでの時間を運動タイミングとして検出する時間演算手段と、
前記所定心拍より前の心拍における前記運動情報に基づいて、前記所定の心拍における時相を推定する時相推定手段と、
前記推定した時相を基準として、前記運動タイミングと色彩との対応を設定するカラーマップ設定手段と、
前記設定した色彩に基づいて、前記運動タイミングの２次元分布又は３次元分布を表す画像データを生成するカラー像生成手段と
を含むことを特徴とする超音波画像処理装置。
【請求項２】
前記カラーマップ設定手段は、前記推定した時相を基準とした時間範囲をカラーマップ設定範囲として設定し、
前記カラーマップ設定範囲内にある前記運動タイミングに対して前記色彩との対応を設定する
ことを特徴とする請求項１に記載の超音波画像処理装置。
【請求項３】
前記カラー像生成手段は、
前記時間演算手段が前記運動タイミングを検出する度に、前記画像データ中の対応する位置に前記設定した色彩の画素を表示すると共に、
前記時間演算手段が前記検出を行う時相が前記所定心拍と異なる心拍へ移ったことに応じて前記表示した色彩の画素の表示を消去する
ことを特徴とする請求項１または２に記載の超音波画像処理装置。
【請求項４】
前記被検体から反射される超音波に基づく受信データに基づいて前記被検体における複数の位置での組織の速度を検出する速度演算手段を更に含み、
前記時間演算手段は、
前記速度演算手段の検出した速度の時間積分を含む処理によって、前記被検体の運動タイミングを求めることを特徴とする請求項１乃至３のいずれか1項に記載の超音波画像処理装置。
【請求項５】
前記運動情報は、前記被検体における組織の歪み、または変位を表す情報であることを特徴とする請求項１乃至４のいずれか1項に記載の超音波画像処理装置。
【請求項６】
前記カラー像生成手段は、
前記運動タイミングに応じた色彩がＢモード画像上の対応する位置に重畳された画像を生成することを特徴とする請求項１乃至５のいずれか１項に記載の超音波画像処理装置。
【請求項７】
前記時間演算手段は、
閾値を設定するための閾値設定手段を更に含み、
閾値よりも大きな値をもつ前記運動情報に対して、運動タイミングの検出を行わせることを特徴とする請求項１乃至６のいずれか1項に記載の超音波画像処理装置。
【請求項８】
前記時間演算手段は、
前記最初の極点候補が抽出された時点における運動情報の値に、１未満の数値を乗算した値を減衰閾値として設定する減衰閾値設定手段を更に備え、
前記第１ピーク抽出手段は、
前記所定の時間範囲における古い時点から新しい時点へと、順次運動情報の極点を探索することにより最初の極点候補を抽出し、
前記所定の時間範囲内であって前記最初の極点候補よりも新しい時間範囲内において、運動情報の値が減衰閾値を下回った場合に前記最初の極点候補を前記最初の極点として抽出する処理を含むことを特徴とする請求項１乃至６のいずれか1項に記載の超音波画像処理装置。
【請求項９】
前記表示制御手段は、
複数の位置での組織の運動情報のうちの少なくとも一つの運動情報の経時変化を示すための経時グラフ画像を生成し、この経時グラフ画像に前記閾値の範囲を示すための画像を生成する経時グラフ生成手段を含むことを特徴とする請求項７に記載の超音波画像処理装置。
【請求項１０】
前記閾値設定手段は、
操作者からの入力に応じて、前記閾値を変更することを特徴とする請求項９に記載の超音波画像処理装置。
【請求項１１】
前記表示制御手段は、
複数の位置での組織の運動情報のうちの少なくとも一つの運動情報の経時変化を示すための経時グラフ画像を生成し、この経時グラフ画像に前記所定の時間範囲を示すための画像を生成する経時グラフ生成手段を含むことを特徴とする請求項１乃至８のいずれか1項に記載の超音波画像処理装置。
【請求項１２】
前記時間演算手段は、
操作者からの入力に応じて、前記所定の時間範囲を設定するための演算範囲設定手段を含むことを特徴とする請求項１１に記載の超音波画像処理装置。
【請求項１３】
前記カラーマップ設定手段は、
異常に早い運動タイミングと、異常に遅い運動タイミングが分別可能に表示されるように前記運動タイミングと色彩の対応を設定すことを特徴とする請求項１乃至１２のいずれか1項に記載の超音波画像処理装置。
【請求項１４】
前記カラーマップ設定手段は、
異常に早い運動タイミングと、正常の運動タイミング、及び異常に遅い運動タイミングが分別可能に表示されるように前記運動タイミングと色彩の対応を設定することを特徴とする請求項１乃至１２のいずれか1項に記載の超音波画像処理装置。
【請求項１５】
前記時相推定手段は、
前記所定心拍より前の心拍における前記運動情報に基づいて、前記被検体内の心臓の現在の心拍における収縮末期時相又は収縮末期以降の遅れた時相を推定する
ことを特徴とする請求項１乃至１４のいずれか1項に記載の超音波画像処理装置。
【請求項１６】
前記運動情報演算手段は、
前記被検体内の組織の動きを追跡し、前記追跡に従った前記被検体の速度の時間的変化に基づいて前記運動情報を求めることを特徴とする請求項１乃至１５のいずれか1項に記載の超音波画像処理装置。
【請求項１７】
前記運動情報演算手段は、
前記被検体の位置又は形状に応じて、前記運動情報の角度補正を行うことを特徴とする請求項１乃至１６のいずれか1項に記載の超音波画像処理装置。
【請求項１８】
前記表示制御手段は、
複数の前記特徴点に関する解析値又は統計値を表す画像を生成することを特徴とする請求項１乃至１７のいずれか1項に記載の超音波画像処理装置。
【請求項１９】
超音波を被検体に対して送信し、前記被検体から反射される超音波に基づく受信データを取得するスキャン手段と、前記受信データに基づいて前記被検体における複数の位置での組織の運動情報を検出する運動情報演算手段を有する超音波画像処理装置に、
所定心拍の所定の時間範囲内における前記運動情報の経時変化を検出し、前記所定の時間範囲内にある最初の極点を抽出し、前記最初の極点に到達するまでの時間を運動タイミングとして検出する機能と、
前記所定心拍より前の心拍における前記運動情報に基づいて、前記所定の心拍における時相を推定する機能と、
前記推定した時相を基準として、前記運動タイミングと色彩との対応を設定する機能と、
前記設定した色彩に基づいて、前記運動タイミングの２次元分布又は３次元分布を表す画像データを生成する機能と
を実行させること特徴とする超音波画像処理装置の制御プログラム。
【請求項２０】
前記色彩との対応の設定は、前記推定した時相を基準とした時間範囲をカラーマップ設定範囲として設定し、
前記カラーマップ設定範囲内にある前記運動タイミングに対して行われる
ことを特徴とする請求項１９に記載の超音波画像処理装置の制御プログラム。
【請求項２１】
前記画像データの生成は、
前記運動タイミングの検出の度に前記画像データの生成を行うと共に、
前記運動タイミングの検出を行う時相が前記所定心拍と異なる心拍へ移ったことに応じて前記画像データの画素の表示を消去する
ことを特徴とする請求項１９に記載の超音波画像処理装置の制御プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体に対して超音波を送受信することによって得られたエコー情報に基づいて被検体の運動情報を示す画像を生成する超音波画像処理装置及び超音波画像処理装置の制御プログラムに関する。
【背景技術】
【０００２】
被検体内を超音波で走査し、被検体内からの反射波から生成した受信信号を基に、当該被検体の内部を画像化する超音波診断装置が知られている。この超音波診断装置は、超音波探触子から被検体内に超音波を送信し、被検体内部で音響インピーダンスの不整合により生じる反射波を超音波探触子で受信して受信信号を生成する。このような超音波診断装置は、心臓の診断にも有用であることが知られている。なかでも、心筋等の生体組織に関して、その機能を客観的かつ定量的に評価することは、非常に有用である。
【０００３】
ここで、最近注目を集めいている治療法として、重度心不全患者に対する心臓再同期療法（ＣＲＴ）が挙げられる。この心臓再同期療法の適用可否の判定や治療の効果判定のための定量的評価に対しても超音波診断装置を用いることが試みられている。
【０００４】
以下、この心臓再同期療法について、簡単に説明する。重度心不全患者は心壁運動の収縮同期不全を併発することが多い。心臓は電気信号の伝導によって動いているが、重度心疾患の患者には心室内伝導障害がおこることがある。心室内伝導障害とは、心筋を運動させる電気信号の伝わる順序にずれが起きてしまうことがあり、このずれにより、本来はほぼ同時に全体に電気信号が伝わる心室において、早く伝わる部分と遅く伝わる部分ができてしまうことがある。この結果、心壁の収縮が同期せず、血液を充分に送り出せなくなり、心不全の状態となる。このような障害に対し、人工的に電気信号を出して心臓に伝わる電気信号の順序を整え、心臓のポンプ機能を助ける治療が心臓再同期療法であり、具体的には、ペースメーカーを胸の皮膚の下にうめこむことによりなされる。
【０００５】
このような心臓再同期療法は、すでに多数の患者に施されており、劇的な症状の改善が確認されている。しかしながらその一方で、心不全症例であって、この心臓再同期療法を適応しても症状の改善が見られない患者も全体の約３割確認されている。これは、心不全症例の原因が収縮同期不全であるか否かを正確に判断できないためである。現状、心臓再同期療法の適応基準は心電波形のＱＲＳ幅が１３０ｍｓｅｃ未満、左室駆出率（EF）が３５％であることとされているが、この基準では心不全だが収縮同期不全ではない患者も含まれてしまうのである。
【０００６】
そこで、超音波診断装置を用いた定量的評価法により収縮同期不全のみを抽出しようという技術が開発されている。このような技術としては、心筋（心臓壁）の運動速度をドプラ法を用いて検出し、その運動速度を演算・解析する技術（特許文献１）が知られている。この技術によれば、心筋の複数部位における運動速度や変位などの経時変化のピークを自動的に検出することができる。そして、所定心時相からこれらのピークに到達するまでの時間を算出し、その時間に応じて心筋の超音波画像に色づけする。このようにして、心筋全体の運動状態の違いをカラー画像として出力することにより、心筋の各部位における運動タイミングの違いを画像化することができる。
【特許文献１】特開平１０−２６２９７０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００７】
しかしながら、上述の従来技術には下記のような課題がある。特許文献１の技術によれば心筋の運動速度を演算・解析しているが、速度は瞬時的な物理量であるから、心筋速度の経時変化は必ずしも心室の収縮状態に応じてのみピークを有するわけではない。とくに心不全症例では、心筋は異常な動作を行うので、心筋速度は多数のピークを有することが知られている。つまり、運動速度からでは、どのピークが有意な心室収縮を示すピークであるかを安定して評価するのが困難なのである。
【０００８】
以下に、図１２、１３を参照してこの関係を具体的に示す。左心室近辺の超音波画像上に指定されたＲＯＩ１の部位の速度の経時変化を実線、ＲＯＩ２で示される部位の速度の経時変化を点線でグラフに表す。グラフは１心周期の速度変化を示しており、Ｒ波近傍で区切った収縮初期から拡張末期までの速度変化を表している。ここで、図１２では、収縮初期区間まで含めた（ａ）の領域でピーク検出を行った場合を示している。このとき、ＲＯＩ１の運動速度のピークはｔｐ１で示される位置で検出され、ＲＯＩ２の運動速度のピークはｔｐ２で示される位置で検出される。しかしながら、図１３に示されるような、収縮初期区間を除いた（ｂ）の領域でピーク検出を行うとピーク検出の結果は全く異なる。ＲＯＩ２の運動速度のピークは、変わらずｔｐ２の位置で検出されるが、ＲＯＩ１の運動速度のピークはｔｐ３の位置で検出されてしまう。図１２、１３上の超音波画像上には、ピークの到達時間に応じて濃淡が示されているが、図１２の場合には中隔側が早く収縮し、側壁側の方が遅く収縮しているように観察されるが、図１３の場合には全く逆に観察されることがわかる。
【０００９】
したがって、特許文献１に示される技術では、心不全症例の場合の有意な収縮を示すピークの検知を誤る可能性を多く含むため、収縮同期不全を確実に検知することが困難であるという課題があった。
【００１０】
本発明は上記課題に鑑みてなされたものであり、被検体内の組織の運動状態を正しく把握できる医用画像処理装置及び医用画像処理装置の制御プログラムを提供することを目的とする。
【課題を解決するための手段】
【００１１】
上記課題を解決するため、実施形態に係る超音波画像処理装置は、超音波を被検体に対して送信し、前記被検体から反射される超音波に基づく受信データに基づいて前記被検体における複数の位置での組織の運動情報を検出する運動情報演算手段と、所定心拍の所定の時間範囲内における前記運動情報の経時変化を検出し、前記所定の時間範囲内にある最初の極点を抽出し、前記最初の極点に到達するまでの時間を運動タイミングとして検出する時間演算手段と、前記所定心拍より前の心拍における前記運動情報に基づいて、前記所定の心拍における時相を推定する時相推定手段と、
前記推定した時相を基準として、前記運動タイミングと色彩との対応を設定するカラーマップ設定手段と、前記設定した色彩に基づいて、前記運動タイミングの２次元分布又は３次元分布を表す画像データを生成するカラー像生成手段とを含むことを特徴とする。
【００１２】
また、実施形態に係る超音波画像処理装置の制御プログラムは、超音波を被検体に対して送信し、前記被検体から反射される超音波に基づく受信データを取得するスキャン手段と、前記受信データに基づいて前記被検体における複数の位置での組織の運動情報を検出する運動情報演算手段を有する超音波画像処理装置に、所定心拍の所定の時間範囲内における前記運動情報の経時変化を検出し、前記所定の時間範囲内にある最初の極点を抽出し、前記最初の極点に到達するまでの時間を運動タイミングとして検出する機能と、前記所定心拍より前の心拍における前記運動情報に基づいて、前記所定の心拍における時相を推定する機能と、前記推定した時相を基準として、前記運動タイミングと色彩との対応を設定する機能と、前記設定した色彩に基づいて、前記運動タイミングの２次元分布又は３次元分布を表す画像データを生成する機能とを実行させること特徴とする。
【発明の効果】
【００１３】
本発明によれば、被検体内の組織の運動状態を正しく把握できる医用画像処理装置及び医用画像処理装置の制御プログラムを提供することができる。
【発明を実施するための最良の形態】
【００１４】
以下、本発明の一実施形態について説明する。本実施形態では、心臓を対象とし、心筋の組織の局所的な運動について２次元画像を用いて評価する場合について述べる。
【００１５】
（構成）
図１は、本発明の一実施形態に係る超音波画像処理装置の構成図である。超音波プローブ１は、電気信号を超音波に変換する複数の超音波振動子を配列した超音波振動子アレイを備え、この超音波振動子アレイにより被検体に対して超音波の送受波を行う。第1の実施形態において、超音波プローブ１は心臓を対象としたセクタプローブであるものとする。
【００１６】
送信ユニット２は、超音波振動子アレイから超音波を送波するための駆動信号を発生するものであり、所定のスキャンラインに向けて超音波ビームが形成されるよう振動子毎に所定の遅延特性を持たせた駆動信号を発生する。
【００１７】
受信ユニット３は、超音波振動子アレイの各超音波振動子で受波された超音波エコー信号に対して遅延加算処理を行って所定のスキャンラインに対応した超音波エコー信号を生成する。
【００１８】
Ｂモード処理ユニット４は、遅延加算処理された超音波エコー信号に対して包絡線検波処理を施すことにより、超音波エコーの振幅強度に対応したＢモード信号を生成する。又、Ｂモード処理ユニット４は、このＢモード信号の所定断面に係る２次元分布を表したＢモード超音波像を生成する。
【００１９】
組織ドプラ処理ユニット５は、直交検波処理、自己相関処理等を行い、遅延加算処理された超音波エコー信号のドプラ偏移成分に基づいて、被検体内で移動している組織の速度、分散、パワーに対応した組織ドプラ信号を求める。組織ドプラ処理ユニット５は、この速度、分散、パワー値の所定断面に係る２次元分布情報を生成する。
【００２０】
運動情報処理ユニット６は、組織ドプラ処理ユニット５の出力する２次元分布情報に基づいて、運動情報を取得するための各処理を実行する。具体的には、２次元分布情報から得られる局所的な組織速度を用いて時間積分処理を行うことによって局所的な変位や歪みなどの運動情報を演算し、その２次元分布情報を提供する。ここでいう変位（ｄｉｓｐｌａｃｅｍｅｎｔ）とは、被検体内のある組織が任意の時間内にどれだけ移動したかを示すための値であり、歪み（ｓｔｒａｉｎ）とは、被検体内の組織上の２点間の距離が任意の時間内にどれだけ変化したかを示すための値である。
【００２１】
ここで、変位や歪みを得るための時間積分処理は心時相に同期させて行うものとする。具体的には、図示しない心電図ユニットより被検体の心電図を取得し、この心電図から所定の心時相（ＱＲＳ波）を特定し、この心時相にさしかかる度に上記の時間積分処理をリセットさせるものとする。このような運動情報の取得の具体的手法に関しては、特開平２００２−２７２８４５号公報に示される技術が適応可能である。
【００２２】
なお、このような変位又は歪みの情報を得るにあたっては、組織の位置を追跡した上での情報を取得するものとする。本実施形態においては、組織の位置を組織ドプラ情報をもとに追跡するものとするが、これに限られない。Ｂモード情報を用いて組織の位置を追跡するものであってもよい。
【００２３】
さらに、このような運動情報の取得にあたっては角度補正を行ってもよい。通常、ドプラ情報は走査線方向に沿った運動に関する情報である。したがって、例えば左室短軸像が得られている場合には、収縮中心に向かう成分をもとに角度補正を行い、心筋の厚みの変化に関する歪みや変位を抽出してもよい。また、例えば心尖像が得られている場合には心内膜面に沿った方向成分に角度補正を行って心筋の長手方向の伸縮に関する歪みや変位を抽出してもよい。
【００２４】
時間演算ユニット７は、第１ピーク抽出機能７１、閾値設定機能７２、演算範囲設定機能７３、及び減衰閾値設定機能７４を有する。第１ピーク抽出機能７１は、運動情報処理ユニット６から取得した運動情報をもとに、歪みや変位の経時変化が、最初にピークをむかえる時点を検出する。さらに、第１ピーク抽出機能７１は所定の時点からこのピークに達する時点までに要した時間を算出する。このピーク到達時間を指定された範囲内の各点ごとに求めることにより、ピーク到達時間の空間分布情報が得られる。これらの演算は、閾値設定機能７２や演算範囲設定機能７３によって設定された範囲内において行われ、減衰閾値設定機能７４によって設定された条件に基づいて、所定量の減衰をピーク検出の条件としている。これらの機能についての詳しい説明はさらに後述する。
【００２５】
表示制御ユニット８は、カラーマップ設定機能８１、カラー像生成機能８２、経時グラフ生成機能８３、心電図生成機能８４を有する。カラー像生成機能８２は、時間演算ユニット７によって得られたピーク到達時間の空間分布情報に基づいて表示画像を生成する。各点におけるピーク到達時間の情報はカラーマップ設定機能８１の設定に基づいてカラー変換され、Ｂモード像に重畳した画像データを生成する。また、経時グラフ生成機能８３は、画像上に指定されたＲＯＩにおける歪みや変位の運動情報の経時変化をグラフ形式として表示する。心電図生成機能８４は、図示しない心電ユニットを介し、被検体の心電図の画像データを生成する。それぞれの機能が生成する画像データの詳細については後述する。表示装置９は、表示制御ユニットで生成された画像データを表示する。
【００２６】
収縮時相推定ユニット１１は、収縮末期時相や収縮遅延時相を推定する。これをもとに収縮同期不全がもっとも顕著に表示されるよう、ピーク到達時間に対するカラー情報の割り当てが決定される。表示制御ユニット８は、この割り当てに基づいてカラー画像の生成を行う。収縮時相推定ユニットの行う具体的な推定の方法については後述する。
【００２７】
なお、本実施形態における各構成は回路基板等のハードウェアとしてそれぞれが実装されるものであってもよいが、ソフトウェアとして組み込まれてもよい。ソフトウェアとして組み込まれる場合はメモリやハードディスクなどの記憶手段に、それぞれの構成の機能を実現させるためのプログラムが記憶される。そして、内蔵されるＣＰＵにそれらのプログラムが展開されることにより、図１に示されるような概念に基づく動作が行われる。
【００２８】
（動作）
図２は本実施形態に係る超音波画像処理装置の動作を示すフローチャートである。以下、このフローチャートに則して説明を行う。ここではとくに、本実施形態の特徴部である時間演算ユニット７の動作を中心に説明する。
【００２９】
まず、ステップＳ１として、運動情報処理ユニット６の取得する運動情報が設定される。操作者は変位、歪みのどちらに基づいて処理を行うか、及び、心臓の短軸方向、長軸方向のどちらから超音波照射を行うかを選択する。運動情報処理ユニット６はこの設定をうけて、どの方向成分の変位または歪みを取得するかを決定し、この方向成分に従った角度補正を行うことにより、診断に適切な変位や歪みを算出する。
【００３０】
図３に、上述の操作者の設定と、変位や歪みが取得される方向成分との関係を表にして示す。まず、短軸での変位を選択した場合を説明する。収縮始期と収縮末期とでは、心筋の厚みが異なる。短軸像はこの変化を観察するのに適しているため、厚み方向の心筋の変位が算出されるよう角度補正が行われる。具体的には、図４に示されるように、操作者による指定や自動輪郭抽出に基づいて設定された心臓の中心点に向かう方向の変位を、組織の各点において算出する。図４中の○は組織を示し、矢印は算出する方向を示す。図４に示されるように、算出するべき変位の方向と、超音波プローブ１からの超音波照射方向は各点で異なるため、これを補正する。収縮にともなって心筋は中心方向に変位するので、得られる値は心臓の収縮にともなって増加する。収縮にともなって値が増加する場合を図３に（＋）で表し、減少する場合を（−）で表す。短軸像での歪みを選択した場合も同様で、心臓中心に向かう方向の歪みが図５のような角度関係に基づいた補正により求められる。図中の○は組織を示し、２つの○の矢印方向の位置関係をもとに、歪みが求められる。心筋の厚みは収縮にともなって厚くなるので、得られる値は収縮にともなって増加する。
【００３１】
また、長軸での変位を選択した場合を説明する。心臓は、その収縮にともない長手方向（長軸方向）に縮む。長軸像はこの変化を観察するのに適しているため、長手方向の変位が算出されるよう角度補正が行われる。具体的には図６に示されるように、操作者による指定や自動輪郭抽出に基づいて設定された心臓の長軸方向のの変位を、組織の各点において算出する。図３に示されるように、算出するべき変位の方向と、超音波プローブ１からの超音波照射方向は各点で異なるため、これを補正する。収縮にともなって心筋はプローブから向かう方向に変位するので、プローブに向かう方向を正とすると、得られる値は心臓の収縮にともなって増加する。長軸像での歪みを選択した場合も同様で、心臓の長軸方向の歪みが図７に示されるような角度関係に基づいた補正により求められる。
【００３２】
次に、ステップＳ２において、演算範囲設定機能８３が時間演算範囲を設定する。時間演算範囲とは、上述の歪みや変位の運動情報の経時変化のピークを求める時間範囲である。収縮同期不全の診断においては、一心拍内の収縮期間内での心筋各部の収縮タイミングの差を検出することが有用であるので、収縮期間内に時間演算範囲が合わされる。具体的には、操作者が表示装置９に表示される心電波形あるいは指定された領域の歪みや変位の経時変化グラフ上のカーソルを移動することにより、時間演算範囲の始点と終点を定め、これを演算範囲設定機能８２が認知し、所定のメモリに記憶する。なお、具体的な時間演算範囲の設定方法としては、Ｒ波からの何秒の範囲かを数値入力したり、一心拍分の時間を１とした相対的な時間単位での入力でもよい。このようにして、時間演算範囲の始点ｔｓ、終点ｔｅが定まり、時間演算範囲ｔｓ≦ｔ≦ｔｅが設定される。
【００３３】
次に、ステップＳ３として、操作者が、超音波プローブ１を被検体に当接させ、被検体の画像を撮像する。操作者は所望の画像が得られるよう、超音波プローブ１の位置を表示装置９に表示されたリアルタイムの画像を見ながら調節する。これにより撮像かスタートされる。
【００３４】
次に、ステップＳ４として、画像化されている各部の変位や歪みの運動情報が入力値となり、この入力値のピークの探索が開始される。時間演算ユニット７の第１ピーク抽出機能８１により、運動情報処理ユニット６から取得した運動情報のうち、所定の１心拍内でのｔ＝ｔｓの時点における運動情報の空間分布から、ピークの探索を開始する。
【００３５】
次に、ステップＳ５として、第１ピーク抽出機能８１は、時刻ｔにおける入力値が閾値設定機能８２の設定する閾値を超えているか否かを判断する。運動情報は比較的微小なピークを複数有することが多い。このようなピークは本質的でない様々な事情によって現れるものであり、心筋の運動タイミングを示すピークとしては適切ではない。したがって閾値を設けることにより、このようなピークを検出しないようにする。この閾値は操作者により適宜変更可能であり、これを閾値設定機能８２が認知し、所定のメモリに記憶するものとしてもよいし、予め装置に記憶された値を用いてもよい。また、具体的な閾値の適正値は１ｍｍ以下であることが実測により確かめられている。入力値が閾値を超えている場合はステップＳ８へと移り、入力値が閾値を超えていない場合はステップＳ６へと移る。
【００３６】
ステップＳ６ではｔが更新され新たな時点でのピーク探索が継続される。ここでステップＳ７においてｔ＜ｔｅの場合は、有効なピークが検出されないまま探索終了となるが、そうでない場合はステップＳ５に戻る。
【００３７】
ステップＳ８では、第１ピーク抽出機能８１が、時刻ｔにおける入力値がピークであるか否かを判断する。具体的は、入力値を時間微分し、その値が０である時や、微分係数の極性が変化した時点をピークとして検出する。ピークであった場合には、ｔｐ１＝ｔをピーク候補として記憶し、ステップＳ１１にうつる。ピークでなかった場合にはステップＳ９へとうつる。
【００３８】
ステップＳ９ではｔが更新され、新たな時点でのピーク探索が継続される。ここでステップＳ１０においてｔ＜ｔｅの場合は、ｔｓとｔｅの両時点における入力値のうち大きい方を第１ピークとして検出する。ｔｅ≦ｔの場合はステップＳ５に戻る。
【００３９】
ステップＳ１１ではｔが更新される。そして、ステップ１２において、入力値の絶対値が所定の値を下回ったか否かを判断する。ここで所定の値とはステップＳ８におけるｔｐ１の入力値Ｘに減衰閾値設定機能８４の設定する値αを乗じたものである。αは操作者によって設定可能であり、これを減衰閾値設定機能８４が認知し、所定のメモリに記憶するものとしもよいし、予め装置に記憶された値を用いてもよい。また、具体的なαの適正値は０．５〜０．８であることが実測により確かめられている。入力値が所定値を下回っていた場合にはｔｐ１の時点における入力値Ｘが第１ピークとして検知される。そうでない場合には、ステップＳ１３に移る。
【００４０】
ステップＳ１３においては、時点ｔおける入力値がピークであり、かつ時点ｔにおける入力値Ｘ´がｔｐ１における入力値Ｘよりも大きいか否かを判断する。両方の条件を満たす場合にはステップＳ１４にうつり、ピークの候補が変更され、新たにｔｐ１＝ｔとして記憶される。そしてステップＳ１１へ戻る。ステップＳ１３において上述の両方の条件を満たさなかった場合はステップＳ１５へうつり、ｔ＜ｔｅの場合は、ｔｐ１の時点における入力値Ｘが第１ピークとして検知される。そうでない場合はステップＳ１１へ戻る。
【００４１】
このようにして、ｔｓからｔｅまで順次ピークの探索が行われる。なお、ピークの探索は運動情報処理ユニット６から得られた運動情報の２次元分布における各点について行われるものとする。このようにして、第１ピーク到達時間の２次元分布情報が得られる。
【００４２】
ステップＳ１６においては、カラー像生成機能８２が、この第１ピーク到達時間の２次元分布情報を、カラーマップ設定手段８１の設定に従ったカラー情報に変換して、表示装置９に表示させる。このようにして、第１ピーク到達時間の２次元分布情報が一目瞭然に操作者に対して示されることとなる。
【００４３】
なお、上述の説明ではとくに触れていないが、このような動作はリアルタイムで受信されるデータに対して逐一行うものとしてもよいし、過去に撮像した画像を呼び出してきて処理するものとしてもよい。過去に撮像した画像を呼び出す場合はステップＳ１の前に過去画像の呼び出しのステップが付加され、ステップＳ３の超音波プローブによる撮像が捨象される。この場合、表示制御ユニット８は操作者の指定した所望の心拍に対し上述の処理を行い、静止画を生成するものとしてもよい。
【００４４】
また、リアルタイムで行う場合は、各心拍ごとにステップＳ４からステップＳ１６の動作が繰り返される。これによって各心拍ごとに逐一画像が更新される。さらに、表示の処理は、ｔｓからｔｅの全ての範囲のピークの探索が終了した後にステップＳ１６において表示の処理を行ってもよい。また、一つのピークが検出されるたび、対応する位置にカラー情報を表示し、新たな心拍に入ると全てのカラー情報の表示をリセットするものとしてもよい。後者の例によれば、ＱＲＳ波の立ち上がりでカラー情報がリセットされたＢモード像上に、第１ピークを早くむかえた部分から順に彩色が施されるといった態様が繰り返し行われる形で表示される。
【００４５】
ここで、カラーマップ設定機能８１の行う、第１ピーク時間のカラー情報への変換の方法について具体的に説明する。この変換は、カラーマップ設定機能８１が設定するカラーマップに従って行われる。カラーマップとは、どのピーク時間に対してどの色を割り当てるかの設定であり、緑→赤といったようなグラデーションとなるようにカラー情報がそれぞれの時点に対して割り当てられるものである。これは、操作者が手動で定めるものとしてもよいし、装置によって自動設定するものとしてもよい。
【００４６】
また、上記では、緑→赤のグラデーションとなるカラーマップを例示したが、さらに多色を用いたカラーマップを選択できてもよい。まず、２色のカラーマップを用いる場合の色の割り振りは、緑が正常な収縮であるのに対し、赤が遅延している収縮を示すといったように割り振ると好適である。これに対し、３色のカラーマップを用いる場合は、異常な早期収縮を青、正常な収縮を緑、遅延している収縮を赤で表現してもよい。
【００４７】
装置による自動設定の具体的手法としては、例えば、ＱＲＳ波から６００ｍｓｅｃの時点を境として色が変化するようにといったように、心電波形の形状を利用して一定の数値を定めておいてもよいが、本実施形態においては、収縮時相推定ユニット１１が収縮末期時相を自動的に検出しこれを基準に、カラーマップ設定機能８１の設定が行われるものとする。
【００４８】
以下、収縮時相推定ユニット１１の動作について説明する。収縮末期時相の自動検出の具体例としては、本発明の発明者による特許出願の特願平２００４−１６１７９４に示される手法が適応可能である。これは、健常者の例に基づいて心拍数のばらつきを考慮して収縮末期時相を推定する方法である。
【００４９】
しかしながら、重度心不全症例の場合は、健常者の場合と収縮末期時相が異なることが多い。したがって、本実施形態においては、現症例における運動の状態に応じて心拍ごとに収縮末期時相や収縮遅延時相を推定する。具体的な収縮末期時相や収縮遅延時相の推定例としては、観察を行おうとしている直前の心拍における所定領域に対し、以下の４つの処理のいずれかを行う。
【００５０】
１各点で運動情報の重心時間を求め、これを領域内で平均して心筋全体での平均重心時間を求める。２領域内で変位や歪みの運動情報を平均し、この平均運動情報の重心時間を求め、心筋全体での平均重心時間を求める。３各点で運動情報のピークを求め、これを領域内で平均して心筋全体での平均ピーク時間を求める。４領域内で変位や歪みの運動情報を平均し、この平均運動情報のピークを求め、心筋全体での平均ピーク時間を求める。
【００５１】
操作者がこれらのいずれかを選択可能であるとしてもよい。
【００５２】
このようにして推定された収縮末期時相や収縮遅延時相に、観察しようとしている心拍におけるカラーマップに反映させる。具体的には推定された収縮末期時相を基準とした所定範囲をカラーマップの範囲としてもよいし、領域内での運動平均時間情報の分散を求めて収縮の遅延時間量を推定し、これを基準としてカラーマップを決定してもよい。
【００５３】
（表示）
次に、具体的な表示例を参照してさらに本実施形態の機能について説明する。図８は本実施形態における表示例を示す図である。表示制御ユニット８は、カラー像部１１０、全体の心電波形部１２０、指定領域の心電波形部１３０、運動情報の経時グラフ部１４０を表示装置９に表示させる。
【００５４】
カラー像生成機能８２は、被検体のＢモード画像に上述の第１ピーク到達時間をカラー情報として重畳して表示するカラー形態像１１１をあらわすカラー像部１１０を生成する。また、カラー像部１１０には、演算範囲設定機能７３によって設定された時間演算範囲１１２が示される。図８では、ｔｓ及びｔｅは、指定領域の心電波形部１３０及び運動情報の経時グラフ部１４０表示時間領域を１として、ｔｓ＝０．０５、ｔｅ＝０．５８とした場合を示す。また、カラーマップ設定機能８１の設定するカラーマップがカラーバー１１３により示される。図８の場合では、表示時間領域の始点から０．０９秒から０．５６秒の範囲でカラー情報を割り振っている。図８では便宜上、カラー情報を濃淡によって代替しているが、実際には緑→赤といったようなグラデーションのカラー情報が割り当てられる。この表示例では、左室の長軸像が表されており、中隔側１１４で異常な早期収縮がおこっている。したがって、中隔側１１４では、遅延を生じている側壁側１１５よりも濃淡が淡くなっている。
【００５５】
心電図生成機能８４は、図示しない心電図ユニットにより被検体からの心電信号から得られた全体の心電波形部１２０及び指定領域の心電波形部１３０を生成する。表示される範囲は、選択された画像情報の全範囲にわたって表示してもよいし、カラー像部１１０に表示されているＢモード像に対応する時点を基準とした数心拍であってもよい。操作者はこの心電波形１２０を参照し、どの心拍で評価を行うか設定する。設定方法は、具体的にはカーソル１２２を動かすことによって、操作者が手動で定めるものとしてもよいし、一般的な心電波形から規定されるＱＲＳ波の立ち上がりを基準として自動認識するものとしてもよい。また、横波線で示される領域は時間演算範囲の外を示すものである。さらに、全体の心電波形部１２０には収縮末期時相１２３が表示される。操作者はこの収縮末期時相１２３の前後が十分含まれるようにカーソル１２２を動かすことにより、適正な表示時間領域を設定することができる。
【００５６】
さらに、全体の心電波形部１２０上にはカラーマップの基準時点１２４が表示される。この基準時点１２４に従って、緑→赤といったようなグラデーションとなるようにカラーマップがそれぞれの時点に対して割り当てられる。上述したように図８の場合は０．０９秒から０．５６秒にわたってカラー情報が割り振られている。これは、操作者が手動で定めるものとしてもよいし、装置によって自動設定するものとしてもよい。自動設定の具体的手法としては、例えば、ＱＲＳ波から６００ｍｓｅｃなどといったように、心電波形の形状を利用して一定の数値を定めておいてもよいし、上述の収縮末期時相１２３を基準に定めるものとしてもよい。
【００５７】
指定領域の心電波形部１３０は、全体の心電波形部１２０を用いて指定された表示時間領域に対応して心電波形を拡大表示するものである。指定領域の心電波形部１３０は、運動情報の経時グラフ部１４０と対照可能に表示される。
【００５８】
経時グラフ生成機能８３は、カラー像部１１０上に操作者によって指定されたＲＯＩ１１４のにおける運動情報の経時変化を示すグラフである運動情報の経時グラフ部１４０を生成する。操作者は、収縮時相の異常な部位がカラー像部１１０を参照して把握できるので、カラー像部１１０をガイドとしながら任意の部位にＲＯＩを指定し、局所的な運動情報の経時変化の解析を行うことが可能である。図８の場合ではＲＯＩ１の運動情報が実線１４１、ＲＯＩ２の運動情報が破線１４２によって表される。
【００５９】
さらに、表示制御ユニット８は、このようなＲＯＩによって指定された局所的な運動情報の経時変化についての解析値を表示してもよい。図９にこのような解析値をテーブル形式で表示した例を示す。ここでＣＨ１〜ＣＨ４は操作者が指定したＲＯＩに対応するものであり、本表示例においては４つのＲＯＩ（ＣＨ）が指定されている。各ＣＨの数値は各々における第1ピークの値（ｍｍ）を示し、括弧内の数字は第１ピークの到達時間（ｍｓｅｃ）を示す。１行目の第１ピーク値の欄には、ＣＨ１〜４の第１ピーク値の統計値が表示される。μは平均値、σは標準偏差、Δは最大値と最小値の差である。また、２行目の第１ピーク時間の欄も同様で、ＣＨ１〜４の第１ピーク到達時間の統計値が表示される。μは平均値、σは標準偏差、Δは最大値と最小値の差である。
【００６０】
さらに、経時グラフ部１４０においては、演算範囲設定機能８３の設定する時間演算範囲の外が横破線の領域でマスク表示され、閾値設定機能の設定する閾値が点の領域でマスク表示される。操作者はこれらの領域を参照しながら時間演算領域や閾値を設定することができるものとする。なお、図８では正の極性を持つ運動情報を演算の対象としているがこれに限られない。負の極性の運動情報に対して演算を行う場合は負側の領域にマスク表示を行い、運動情報の絶対値に対して演算を行う場合は負側の領域にマスク表示を行うものとする。マスク表示は図８に示されるような領域表示に限られるものではなく、閾値に相当する位置に破線や実線を表示してもよいし、グラフの数値軸の閾値の位置に矢印を表示してもよい。
【００６１】
（効果）
次に本実施形態の効果について以下に列挙する。本実施形態によれば、心筋の変位や歪みなどの運動量を検出し、この運動量の経時変化を演算・解析している。変位や歪みは時間積分を経た物理量であるから、図１２、１３に示されるような速度を用いた場合と異なって経時変化は比較的安定する。したがって、本実施形態によれば運動速度を演算・解析する特許文献１のような技術に比べて、心臓の異常な収縮を示す運動情報のピークを確実に検出できる。
【００６２】
さらに、本実施形態によれば、運動情報の第１ピークの到達時間をカラー情報に変換してＢモード像に重畳している。したがって、異常な収縮を示す運動情報のピークが、その後の収縮による全体動に引きづられて出現するピークよりも小さい場合でも、異常な収縮を示す時点をピークとして検出することができる。したがって、重度心不全症例での収縮初期における異常な収縮を安定的かつ簡便に把握することができる。したがって、収縮同期不全であるか否かを的確に判断することが可能となるため、心臓再同期療法が有効な患者を的確に判定できる。つまり、患者に対して効果の少ない不要な施術を行うことを避けることができるため、患者や術者の負担を大きく軽減できる。  
以下に、図１０を参照してこの関係を具体的に示す。左心室近辺の超音波画像上のＲＯＩ１を実線、ＲＯＩ２を点線に対応させて、１心周期における各部位の変位をグラフに表す。ここで、収縮初期区間まで含めた（ｃ）の領域でピーク検出を行った場合、ＲＯＩ１の変位を示す曲線のピークはｔｐ４で示される収縮初期の位置で異常な挙動を示すピークを有している。しかし、全体動にひきずられて収縮後期になって現れているｔｐ５の位置でもピークを有しているため、機械的にピーク検出を行うとより高い変位を示すｔｐ５の位置をピークとして検出する。このｔｐ５のピークはＲＯＩ２の変位のピークと同じ位置になるため、双方同じ時点でピークを有すると検知される。つまり、ＲＯＩ１の部位とＲＯＩ２の部位はあたかも同期して収縮しているかのような解析結果となってしまうのである。これに対して本実施形態の図８では中隔側で生じている異常な早期収縮を確実に検出できており、遅延を生じている側壁側よりも中隔側が早く収縮し、左室壁間の収縮非同期が生じていることが一目で認識できる。
【００６３】
また、本実施形態においては、所定の閾値を超えた運動情報の入力値のみに対して時間演算を行っており、運動情報の比較的微小なピークを第１のピークとして検出しない。このような微小なピークは本質的でない様々な事情によって現れるものであり、心筋の運動タイミングを示すピークとしては適切ではない。したがって、本実施形態によれば、心臓の異常な収縮を示す運動情報のピークを確実に検出できる。
【００６４】
また、心腔内にも低い値の運動情報が検出されることが一般的であるので、このような閾値を設けない場合には、心腔内にもピークが検出されてカラー情報が重畳されてしまう。しかしながら、収縮同期不全の判定を行う場合にはこのような情報は不要であり、観察の邪魔になる。本実施形態によれば、演算する入力値に閾値を設けることにより、比較的大きな運動情報が検出される心壁の運動情報のピークを中心として検出することができる。したがって、ほぼ心筋のみにカラー情報が重畳された画像を表示することが可能であり、心筋の機能判定が容易となる。
【００６５】
また、運動情報として歪みを用いる場合、歪みは変位よりもスペックルノイズの影響をうけやすいため、その経時変化には微小なピークが存在することが多い。しかしながら、本実施形態によれば、演算する入力値に閾値を設けることにより、比較的大きな運動情報として現れる心筋の収縮に起因する歪みを検出することができる。したがって、本実施形態によれば、心臓の異常な収縮を示す運動情報のピークを確実に検出できる。
【００６６】
さらに、上述のような閾値レベルの最適化のためには図８に示されるような運動情報の経時グラフを参照して、ノイズ除去に最適な閾値設定を行うのが効果的である。本実施形態によれば図８の点の領域に示されるように、閾値によって演算対象外となる範囲がグラフに重畳されてマスク表示される。したがって、操作者は最適な閾値を直感的に把握できる。
【００６７】
一方、本実施形態によれば、減衰閾値を設定することが可能である。心不全症例においては、心壁は不規則な動作をともなうことが多いため、当該心筋部位の収縮末期にのみ必ずピークを有するとは限らない。図１１に示されるように、ピークを３ヶ所持つような場合もある。しかしながら、図１１の場合の収集末期は３番目のｔｐの時点である。これに対して、本実施形態によればピークの検出後、減衰閾値に基づく値まで単調減少するピークのみを抽出する。これにより、不明確なピークを抽出せず、心臓の異常な収縮を示す明確なピークのみを検出できる。
【００６８】
また、本実施形態によれば時間演算範囲を設定することが可能である。短軸像を用いた演算の場合には等容収縮期での回転成分によってＱＲＳ波の立ち上がり直後にピークが現れることが多い。しかしながら、本実施形態によれば、時間演算範囲を設定することにより、このような領域を演算の対象外とすることが可能である。したがって、本実施形態によれば、心臓の異常な収縮を示す運動情報のピークを確実に検出できる。
【００６９】
さらに、上述のような時間演算範囲の最適化のためには図８に示されるような運動情報の経時グラフ部１４０を参照して、最適な設定を行うのが効果的である。本実施形態によれば図８の点の領域に示されるように、演算対象外となる範囲がグラフに重畳されてマスク表示される。したがって、操作者は最適な時間演算範囲を直感的に把握できる。
【００７０】
また、本実施形態によれば、収縮末期時相や収縮遅延時相を自動的に推定し、カラーマップに反映させている。したがって、操作者がカラーマップの設定を行わずとも適切なカラーマップが設定され、操作者の負担が軽減される。さらに、本実施形態においては、過去の心拍における運動情報から次回の心拍における収縮末期時相や収縮遅延時相を推定し、カラーマップの設定に反映している。したがって、１心拍分の撮像の終了を待ってからカラーマップの設定を行う必要がない。つまり、リアルタイムに各心拍において適切なカラーマップが設定できる。これにより、多数の心拍における連続的な観察が容易となり、被検体の運動状態の把握が容易となる。
【００７１】
さらに、本実施形態によれば、カラーマップの設定を２色及び３色以上の色を用いて設定可能である。したがって、２色のカラーマップにより異常な早期収縮のみを浮き立たせるように表示することができると同時に、３色のカラーマップにより異常な早期収縮とともに異常な遅延収縮も一目瞭然に把握することができる。したがって、診断の目的や被検体の状態に応じて、適切な表示態様での画像を取得することができる。したがって、本実施形態によれば、操作者の診断に際しての負担が軽減される。
【００７２】
（変形例）
本発明は、上述のような実施形態に留まらず、その要旨を逸脱しない範囲において様々な変形例が可能であることはいうまでもない。
【００７３】
例えば、被検体に対して運動や薬物により負荷を与え、その前後においてピーク到達時間をカラー表示する技術が本発明者らによって特許出願（特願２００４−２１０１１０）されている。本発明をこれに適応し、第１ピーク抽出機能、閾値設定機能、演算範囲設定機能及び減衰閾値設定機能の少なくともいずれかを、上述のような負荷前後での運動状態の比較のために用いてもよい。このようにしても上述の実施形態で説明したものと同様の効果が得られる。
【００７４】
また、本実施形態では、歪みや変位の運動情報の解析に対して第１ピーク抽出機能、閾値設定機能、演算範囲設定機能及び減衰閾値設定機能を用いたが、これに限られない。速度や分散、パワー値の運動情報に対してこれらの機能を適用可能である。
【００７５】
さらに、本実施形態においては、第１ピーク抽出機能、閾値設定機能、演算範囲設定機能、減衰閾値設定機能、収縮時相推定ユニットを備えているが、必ずしもその全てを備える必要はない。これらのうちの少なくとも一つを含む組み合わせによっても、上述の実施形態に示した効果のうちの対応する効果を奏することは明らかである。
【００７６】
また、上述の実施形態では、受信信号が２次元空間で得られている場合について述べてきたがこれに限られない。受信信号が複数の２次元断面について得られている場合や、３次元空間で得られている場合、さらにはｐｏｌａｒ−ｍａｐのように３次元空間から特徴量を得て２次元平面に演算結果を投影表示したような場合であっても、上述の実施形態の次元の拡張や変換により本発明の適応が可能である。
【００７７】
また、本実施形態では超音波プローブから受信信号を得て表示装置に表示するまでを記しているがこれに限られない。本実施形態のうち特徴的な構成は、運動情報演算ユニット、時間演算ユニット、表示制御ユニット及び収縮時相推定ユニットである。したがって、これらの構成のうち少なくとも一つを有し、外部からの情報を処理する装置としても本発明の適応が可能である。具体的には、ＰＣやワークステーションなどのコンピュータにより、超音波診断装置とは切り離して処理を行う超音波画像処理装置であってもよい。
【図面の簡単な説明】
【００７８】
【図１】本発明に係る実施形態の構成を示すためのブロック図。
【図２】図１の実施形態における動作を示すためのフローチャート。
【図３】図１の実施形態において取得する運動情報の一覧を示す図。
【図４】短軸像における変位を運動情報として取得する場合の角度補正を説明するための模式図。
【図５】短軸像における歪みを運動情報として取得する場合の角度補正を説明するための模式図。
【図６】長軸像における変位を運動情報として取得する場合の角度補正を説明するための模式図。
【図７】長軸像における歪みを運動情報として取得する場合の角度補正を説明するための模式図。
【図８】図１の実施形態において表示制御ユニットが生成する画像の一例を示す図。
【図９】図１の実施形態において表示制御ユニットが生成する演算・解析結果を示すための画像の一例を示す図。
【図１０】図１の実施形態における、第１ピーク抽出機能による効果を説明するための図。
【図１１】図１の実施形態における、減衰閾値設定機能による効果を説明するための図。
【図１２】従来技術を説明するための図。
【図１３】従来技術を説明するための図。
【符号の説明】
【００７９】
１  超音波プローブ
２  送信ユニット
３  受信ユニット
４  Ｂモード処理ユニット
５  組織ドプラ処理ユニット
６  運動情報処理ユニット
７  時間演算ユニット
８  表示制御ユニット
９  表示装置
１１  収縮時相推定ユニット
７１  第１ピーク抽出機能
７２  閾値設定機能
７３  演算範囲設定機能
７４  減衰閾値設定機能
８１  カラーマップ設定機能
８２  カラー像生成機能
８３  経時グラフ生成機能
８４  心電図生成機能
１１０  カラー像部
１１１  カラー形態像
１１２  時間演算範囲
１１３  カラーバー
１２０  全体の心電波形部
１２１  心電波形
１２２  カーソル
１２３  収縮末期時相
１２４  カラーマップの基準位置
１３０  指定領域の心電波形部
１４０  運動情報の経時グラフ部 
【図１】
0004805669.tif 000002
【図２】
0004805669.tif 000003
【図３】
0004805669.tif 000004
【図４】
0004805669.tif 000005
【図５】
0004805669.tif 000006
【図６】
0004805669.tif 000007
【図７】
0004805669.tif 000008
【図８】
0004805669.tif 000009
【図９】
0004805669.tif 000010
【図１０】
0004805669.tif 000011
【図１１】
0004805669.tif 000012
【図１２】
0004805669.tif 000013
【図１３】
0004805669.tif 000014