(57)【特許請求の範囲】
【請求項１】
超音波を送受して得られた信号に基づいて血流が含まれる座標系において各座標の血流ベクトルを得ることにより、複数座標の血流ベクトルで構成された各ベクトルフレームからなる複数のベクトルフレームを生成するベクトル演算部と、
互いに隣接する２つのベクトルフレーム間に補間処理を適用することにより各補間フレームを生成し、当該ベクトルフレーム間に１又は複数の補間フレームを追加する補間処理部と、
複数のベクトルフレームとそれらのベクトルフレーム間に追加された複数の補間フレームで構成されるフレーム列に基づいて、血流の仮想的な各粒子の移動先を導出する粒子演算部と、
を有する、
ことを特徴とする超音波診断装置。
【請求項２】
請求項１に記載の超音波診断装置において、
前記補間処理部は、互いに隣接する２つのベクトルフレーム間における血流ベクトルに基づく補間処理により、複数座標の補間ベクトルで構成された前記各補間フレームを生成する、
ことを特徴とする超音波診断装置。
【請求項３】
請求項１または２に記載の超音波診断装置において、
前記粒子演算部は、複数時相に対応した前記フレーム列に含まれる各フレーム内における前記各粒子の血流ベクトルに基づいて、当該フレームよりも後の時相のフレーム内における当該各粒子の移動先の座標を算出する、
ことを特徴とする超音波診断装置。
【請求項４】
請求項３に記載の超音波診断装置において、
複数時相に亘る前記各粒子の移動先の座標を画像内に示した血流表示画像を形成する表示処理部をさらに有する、
ことを特徴とする超音波診断装置。
【請求項５】
請求項４に記載の超音波診断装置において、
前記表示処理部は、複数時相に亘る前記各粒子の移動先の座標を点と線の少なくとも一方による軌跡で示した血流表示画像を形成する、
ことを特徴とする超音波診断装置。
【請求項６】
請求項５に記載の超音波診断装置において、
前記表示処理部は、複数時相に亘る前記各粒子の軌跡を示すにあたり、当該粒子の明示期間に対応した時相の軌跡部分と当該明示期間以外の他期間に対応した時相の軌跡部分とにおいて表示態様を異ならせる、
ことを特徴とする超音波診断装置。
【請求項７】
請求項６に記載の超音波診断装置において、
前記表示処理部は、前記各粒子の軌跡が表示される時相から所定時間前の時相までを当該粒子の明示期間とする、
ことを特徴とする超音波診断装置。
【請求項８】
請求項６に記載の超音波診断装置において、
前記表示処理部は、前記各粒子が生成された時相から所定時間後の時相までを当該粒子の明示期間とする、
ことを特徴とする超音波診断装置。
【請求項９】
請求項１から８のいずれか１項に記載の超音波診断装置において、
超音波を送受して得られた信号に基づいて前記血流を含んだ心臓の超音波画像を形成する画像形成部と、
前記超音波画像内における２つの特徴点を結んだ生成ライン上に前記血流の仮想的な複数の粒子を生成する粒子生成部と、
をさらに有する、
ことを特徴とする超音波診断装置。
【請求項１０】
請求項９に記載の超音波診断装置において、
複数時相に亘って形成される前記超音波画像内における心臓の形状変化に追従するように前記２つの特徴点を移動させることにより、複数時相に亘って心臓の動きに前記生成ラインを追従させる、
ことを特徴とする超音波診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波診断装置に関し、特に、血流の移動情報を得る技術に関する。
【背景技術】
【０００２】
血流に超音波を送受して得られる受信信号から血流の移動情報（運動情報）を得る技術が知られている。例えば、特許文献１には、生体内の血流などの流体に対して超音波を送受して得られる受信信号に基づいて、観測面内の複数点において流体に関する２次元の速度ベクトルを得る技術が記載されている。観測面内の複数点における２次元の速度ベクトルの分布から、流体の流れを示す流線などの診断情報を得ることが可能になり、例えば、心臓等の診断への応用が期待される。
【０００３】
また、特許文献２には、生体内の血流を複数の表示要素の運動として表現した画像を形成する超音波診断装置が開示されている。各表示要素としては、血流の仮想的な粒子が好適であり、特許文献２に記載された装置は、各粒子の現フレームにおける速度ベクトルに基づいて、その粒子の次フレームにおける位置（移動先）を決定し、例えば、複数フレームに亘って各粒子の移動先を追跡することにより得られる各粒子の軌跡等を表示する。これにより、例えば、血流内における渦流、乱流、停滞などの状態を視覚的に直感的に確認できるようになる。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１３−１９２６４３号公報
【特許文献２】特開２００８−７３２７９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
特許文献１，２に記載された技術は、血流の移動情報（運動情報）を得るにあたって実用面での利用価値が極めて高く、これらの技術の更なる応用や改良が期待されている。
【０００６】
本発明は、上述した背景技術に鑑みて成されたものであり、その目的は、超音波を利用して得られる血流の移動情報の精度を高めることにある。
【課題を解決するための手段】
【０００７】
上記目的にかなう好適な超音波診断装置は、超音波を送受して得られた信号に基づいて血流が含まれる座標系において各座標の血流ベクトルを得ることにより、複数座標の血流ベクトルで構成された各ベクトルフレームからなる複数のベクトルフレームを生成するベクトル演算部と、互いに隣接する２つのベクトルフレーム間に補間処理を適用することにより各補間フレームを生成し、当該ベクトルフレーム間に１又は複数の補間フレームを追加する補間処理部と、複数のベクトルフレームとそれらのベクトルフレーム間に追加された複数の補間フレームで構成されるフレーム列に基づいて、血流の仮想的な各粒子の移動先を導出する粒子演算部と、を有することを特徴とする。
【０００８】
上記構成において、血流ベクトルとは、血流の運動に関するベクトル情報であり、例えば、血流内の各座標（各血流部分）における速度と方向を示す速度ベクトルや、血流内の各座標における移動量と方向を示す移動ベクトルなどが好適な具体例である。血流ベクトルは、例えば、特許文献１（特開２０１３−１９２６４３号公報）に記載された技術、つまり２次元の速度ベクトル分布を利用して得ることができるものの、他の公知の技術を利用して血流ベクトルが得られてもよい。
【０００９】
また、上記構成における座標系は、血流ベクトルの位置や大きさや方向などを数学的（幾何学的）に表現するにあたっての基準となる。座標系は、超音波が送受される実空間（血流を含む領域）に対応しており、例えば、超音波を送受するプローブの位置を基準（例えば原点）として定義されるのが一般的であるものの、他の位置が座標系の基準（原点）とされてもよい。
【００１０】
また、各ベクトルフレームは、血流が含まれる座標系の複数座標における血流ベクトルで構成される。例えば、生体内において血流が含まれる断面内で超音波を走査することにより走査面が形成され、複数時相に亘って走査面が繰り返し形成されることにより、複数時相に対応した複数のベクトルフレームが形成される。
【００１１】
また、血流の仮想的な各粒子は、血流（血液の流れ）を解析するにあたって利用される演算上の仮想的な要素であり、例えば、診断対象となる血流内に１つ又は複数の仮想的な粒子が生成される。
【００１２】
そして、上記構成の装置によれば、複数のベクトルフレームに複数の補間フレームが追加されたフレーム列が形成され、そのフレーム列に基づいて血流の仮想的な各粒子の移動先が導出される。複数のベクトルフレームに複数の補間フレームが追加されてフレーム列が形成されるため、複数のベクトルフレームのみでフレーム列が形成される場合に比べ、フレーム列のフレームレートが高められる。そして、フレームレートが高められたフレーム列に基づいて、血流の仮想的な各粒子の移動先が導出される。そのため、フレームレートが低い場合に比べて、移動先の推定精度が高められる。
【００１３】
例えば、複数のベクトルフレームのフレーム間において各粒子の血流ベクトルが比較的大きく変化しているにも関わらず、その変化を考慮せずにそのフレーム間で各粒子の移動先を推定してしまうと、推定された移動先と本来の移動先とは大きくずれてしまう可能性がある。
【００１４】
特に、複数座標における血流ベクトルを得るにあたりカラードプラ法の送受信を利用すると、同じビーム方向に超音波が繰り返し送受されるため、例えばＢモード画像を得る場合に比べてフレームレートが低くなってしまう。その低いフレームレートのままで得られた複数のベクトルフレームのみのフレーム間で各粒子の移動先を推定してしまうと、本来の移動先とは大きくずれてしまう可能性がある。
【００１５】
これに対し、上記装置では、互いに隣接する２つのベクトルフレーム間に、例えば血流ベクトルに基づく補間処理が適用され、そのベクトルフレーム間に１又は複数の補間フレームが追加される。追加される各補間フレームは、例えば複数座標の補間ベクトルで構成されることが望ましい。補間ベクトルは、例えば血流ベクトルに基づく補間処理を適用して得られ、補間ベクトルには、ベクトルフレーム間における血流ベクトルの変化が反映される。
【００１６】
そして、上記装置によれば、複数のベクトルフレームとそれらのベクトルフレーム間に追加された複数の補間フレームで構成されるフレーム列に基づいて、血流の仮想的な各粒子の移動先が導出される。そのフレーム列は、互いに隣接する２つのベクトルフレーム間に１又は複数の補間フレームを追加して構成されており、ベクトルフレーム間における血流ベクトルの変化が反映されたものとなる。
【００１７】
したがって、そのフレーム列に基づいて各粒子の移動先を導出することにより、ベクトルフレーム間における血流ベクトルの変化を考慮しつつ各粒子の移動先を導出することが可能になる。これにより、ベクトルフレーム間における血流ベクトルの変化を無視又は軽視する場合に比べて、移動先の推定精度が格段に高められる。なお、上記装置は、例えば複数のベクトルフレームのフレームレートが比較的低い場合に移動先の推定精度を高めることができることはもちろんであるが、複数のベクトルフレームのフレームレートがそれほど低くない場合においても、移動先の推定精度をさらに高めようとするニーズにも対応できる。
【００１８】
望ましい具体例において、前記補間処理部は、互いに隣接する２つのベクトルフレーム間における血流ベクトルに基づく補間処理により、複数座標の補間ベクトルで構成された前記各補間フレームを生成する、ことを特徴とする。
【００１９】
望ましい具体例において、前記粒子演算部は、複数時相に対応した前記フレーム列に含まれる各フレーム内における前記各粒子の血流ベクトルに基づいて、当該フレームよりも後の時相のフレーム内における当該各粒子の移動先の座標を算出することを特徴とする。
【００２０】
望ましい具体例において、前記超音波診断装置は、複数時相に亘る前記各粒子の移動先の座標を画像内に示した血流表示画像を形成する表示処理部をさらに有する、ことを特徴とする。
【００２１】
望ましい具体例において、前記表示処理部は、複数時相に亘る前記各粒子の移動先の座標を点と線の少なくとも一方による軌跡で示した血流表示画像を形成する、ことを特徴とする。
【００２２】
望ましい具体例において、前記表示処理部は、複数時相に亘る前記各粒子の軌跡を示すにあたり、当該粒子の明示期間に対応した時相の軌跡部分と当該明示期間以外の他期間に対応した時相の軌跡部分とにおいて表示態様を異ならせることを特徴とする。例えば、明示期間に対応した軌跡部分のみを表示する、明示期間に対応した軌跡部分を実線として他期間に対応した軌跡部分を破線とする、又は、明示期間に対応した軌跡部分よりも他期間に対応した軌跡部分の輝度を低下させる等の表示処理により、明示期間に対応した軌跡部分を際だたせて、他期間の軌跡部分により軌跡の表示が混雑することを避ける又は軽減することが望ましい。
【００２３】
望ましい具体例において、前記表示処理部は、前記各粒子の軌跡が表示される時相から所定時間前の時相までを当該粒子の明示期間とすることを特徴とする。この所定時間（第１所定時間）は、固定的に設定された時間でもよいし、ユーザによる調整又は選択により設定されてもよい。
【００２４】
望ましい具体例において、前記表示処理部は、前記各粒子が生成された時相から所定時間後の時相までを当該粒子の明示期間とすることを特徴とする。この所定時間（第２所定時間）は、固定的に設定された時間でもよいし、ユーザによる調整又は選択により設定されてもよい。なお、第２所定時間と第１所定時間は互いに異なっていてもよい。
【００２５】
望ましい具体例において、前記超音波診断装置は、超音波を送受して得られた信号に基づいて前記血流を含んだ心臓の超音波画像を形成する画像形成部と、前記超音波画像内における２つの特徴点を結んだ生成ライン上に前記血流の仮想的な複数の粒子を生成する粒子生成部と、をさらに有することを特徴とする。
【００２６】
望ましい具体例において、前記超音波診断装置は、複数時相に亘って形成される前記超音波画像内における心臓の形状変化に追従するように前記２つの特徴点を移動させることにより、複数時相に亘って心臓の動きに前記生成ラインを追従させることを特徴とする。
【発明の効果】
【００２７】
本発明により、超音波を利用して得られる血流の移動情報の精度が高められる。例えば本発明の好適な態様によれば、複数のベクトルフレームに複数の補間フレームが追加されたフレーム列が形成され、そのフレーム列に基づいて血流の仮想的な各粒子の移動先が導出されるため、ベクトルフレーム間における血流ベクトルの変化を無視又は軽視する場合に比べて、移動先の推定精度が格段に高められる。
【図面の簡単な説明】
【００２８】
【図１】本発明の実施において好適な超音波診断装置の全体構成を示す図である。
【図２】心臓内腔に対応した関心領域の具体例を示す図である。
【図３】流入ラインと流出ラインの他の具体例を示す図である。
【図４】速度ベクトル分布の具体例を示す図である。
【図５】フレーム列を説明するための図である。
【図６】フレーム間の補間処理を説明するための図である。
【図７】複数粒子の生成を説明するための図である。
【図８】各粒子の移動先の演算の具体例を説明するための図である。
【図９】各粒子の移動先の演算の例外処理を説明するための図である。
【図１０】各粒子の移動先の演算の終了条件を説明するための図である。
【図１１】流跡線の具体例を示す図である。
【図１２】血流表示画像の具体例を示す図である。
【発明を実施するための形態】
【００２９】
図１は、本発明の実施において好適な超音波診断装置の全体構成図である。図１の超音波診断装置は、血流（血液の流れ）の移動情報を得る機能を備えており、特に、生体の心臓内における血流の診断に好適である。もちろん、心臓以外における血流が診断対象とされてもよい。
【００３０】
プローブ１０は、例えば生体内の心臓などの診断対象を含む領域に超音波を送受波する超音波探触子である。プローブ１０は、複数の振動素子を備えており、複数の振動素子が電子的に走査制御されて、心臓を含む空間内で超音波ビームが走査される。プローブ１０は、例えば、医師等のユーザ（検査者）に把持されて被検者の体表面上に当接して用いられる。なお、プローブ１０は、被検者の体腔内に挿入して用いられるものであってもよいし、電子的な走査と機械的な走査とを組み合わせた探触子であってもよい。プローブ１０としては例えばコンベックス型が望ましいもののセクタ型やリニア型等であってもよい。
【００３１】
送受信部１２は、送信ビームフォーマーおよび受信ビームフォーマーとしての機能を備えている。つまり、送受信部１２は、プローブ１０が備える複数の振動素子の各々に対して送信信号を出力することにより送信ビームを形成し、さらに、複数の振動素子から得られる複数の受波信号に対して整相加算処理などを施して受信ビームを形成する。これにより、超音波ビーム（送信ビームと受信ビーム）が走査面内において走査され、超音波ビームに対応した受信信号が形成される。なお、超音波の受信信号を得るにあたって、超音波ビームが三次空間内で立体的に走査されてもよいし、送信開口合成等の技術が利用されてもよい。
【００３２】
画像形成部２０は、走査面内から得られる超音波の受信信号に基づいて、超音波画像のデータ（画像データ）を形成する。画像形成部２０は、例えば、超音波の受信信号に対して、検波処理やフィルタ処理やＡＤ変換処理等を施すことにより、Ｂモード画像用のフレームデータを形成する。もちろん、Ｂモード画像以外の公知の超音波画像に係る画像データが形成されてもよい。
【００３３】
ドプラ処理部３０は、超音波ビームに対応した受信信号に含まれるドプラシフト量を計測する。ドプラ処理部３０は、例えば公知のドプラ処理により、血流によって超音波の受信信号内に生じるドプラシフトを計測し、血流についての超音波ビーム方向の速度情報（ドプラ情報）を得る。
【００３４】
速度ベクトル演算部４０は、血流についての超音波ビーム方向の速度情報から、走査面内における２次元の速度ベクトルの分布を形成する。速度ベクトル演算部４０は、例えば特許文献１（特開２０１３−１９２６４３号公報）に説明されるように、血流についての超音波ビーム方向の速度情報に加えて、心臓壁の運動情報を利用して、走査面内の各位置における血流の２次元速度ベクトルを得る。
【００３５】
なお、超音波ビーム方向に沿った１次元の速度情報を利用して、走査面内における２次元の速度ベクトルの分布を形成するにあたっては、公知の様々な手法を利用することができる。もちろん、互いに方向が異なる２本の超音波ビームを形成して、２本の超音波ビームの各々から速度情報を得て、２次元の速度ベクトルを形成するようにしてもよい。
【００３６】
速度ベクトル演算部４０は、超音波が送受される空間に対応した演算用座標系において複数のサンプル点について、各サンプル点ごとに速度ベクトルを得る。例えば、演算用座標系をｘｙｚ直交座標系で表し、超音波の走査面に対応したｘｙ平面内で、各サンプル点ごとに速度ベクトルを得て２次元の速度ベクトルの分布を形成する。なお、超音波の走査に対応した走査座標系、例えばビームの深さ方向ｒとビームの走査方向θによるｒθ座標系において２次元の速度ベクトルの分布が形成されてもよい。
【００３７】
速度ベクトル演算部４０は、複数のサンプル点（複数座標）に対応した速度ベクトルで構成される２次元の速度ベクトルの分布を示す各ベクトルフレームを生成する。また、速度ベクトル演算部４０は、複数時相に亘って複数のベクトルフレームを次々に生成する。
【００３８】
関心領域設定部４２は、画像形成部２０による処理で得られた画像データ内において、心臓内腔に対応した関心領域を設定する。関心領域設定部４２は、画像データ内において心臓内腔の外縁として設定された内腔ラインと、心臓内腔へ流入する血流の流路に設定された流入ラインと、心臓内腔から流出する血流の流路に設定された流出ラインにより囲まれた領域を関心領域とする。
【００３９】
図２は、心臓内腔に対応した関心領域の具体例を示す図である。図２には、画像形成部２０において得られる画像データ２２の具体例が図示されており、図２の画像データ２２内には、周囲を心筋や弁によって囲まれた心臓内腔が含まれている。
【００４０】
関心領域設定部４２は、画像データ２２内において、心臓内腔の外縁として設定された内腔ライン５２と、心臓内腔へ流入する血流の流路に設定された流入ライン５４と、心臓内腔から流出する血流の流路に設定された流出ライン５６により囲まれた領域を関心領域とする。
【００４１】
内腔ライン５２は、心臓内腔の外縁に対応した複数のトレースポイントに基づいて形成される。例えば、画像データ２２に対応した表示画像が表示部８２に表示され、医師等のユーザがその表示画像を見ながら操作デバイス９０を利用して、心臓内腔と心筋の境界上または境界付近にいくつかの（数個程度でよい）トレース基準ポイントを設定する。そしてユーザによって設定されたトレース基準ポイントに基づいて、例えば、補間処理等によりトレース基準ポイント間に複数のトレースポイントが追加される。そして、いくつかのトレース基準ポイントと追加された複数のトレースポイントからなる複数のサンプル点に基づいて内腔ライン５２が形成される。例えば、複数のサンプル点を互いに連結するように内腔ライン５２が形成される。なお、画像データ２２に対する二値化処理等の画像処理により、心臓内腔と心筋の境界が特定され、その境界に沿って内腔ライン５２が形成されてもよい。
【００４２】
流入ライン５４と流出ライン５６は、ユーザからの操作に応じて設定される。例えば、画像データ２２に対応した表示画像を見ながら、医師等のユーザが、始点Ｓと終点Ｅの位置を指定することにより、流入ライン５４と流出ライン５６を設定する。
【００４３】
ユーザにより流入ライン５４と流出ライン５６が初期設定されると、関心領域設定部４２は、内腔ライン５２と流入ライン５４を接続するように流入ライン５４を修正し、内腔ライン５２と流出ライン５６を接続するように流出ライン５６を修正し、さらに、流入ライン５４と流出ライン５６とを接続する。
【００４４】
例えば、関心領域設定部４２は、流入ライン５４の始点Ｓを、その始点Ｓに最も近い内腔ライン５２上のサンプル点（トレースポイント又はトレース基準ポイント）の位置に移動する。また、関心領域設定部４２は、流出ライン５６の終点Ｅを、その終点Ｅに最も近い内腔ライン５２上のサンプル点（トレースポイント又はトレース基準ポイント）の位置に移動する。さらに、関心領域設定部４２は、流入ライン５４の終点Ｅと流出ライン５６の始点Ｓとを結ぶ直線または曲線を形成する。
【００４５】
こうして、関心領域設定部４２は、内腔ライン５２と流入ライン５４と流出ライン５６によって囲まれた領域を形成し、その領域を関心領域として設定する。なお、図２には、流入ライン５４と流出ライン５６を直線とする具体例を示したが、直線以外のラインが利用されてもよい。
【００４６】
図３は、流入ライン５４と流出ライン５６の他の具体例を示す図である。例えば、図３（Ａ）に示す具体例のように、閉曲線の内腔ライン５２が得られている場合には、図３（Ｂ）に示すように、始点Ｓと終点Ｅを内腔ライン５２に沿って結んだ曲線状の流入ライン５４と流出ライン５６が形成されてもよい。
【００４７】
図１に戻り、内腔ライン速度演算部４４は、画像形成部２０において形成された画像データに基づいて、内腔ライン（図２の符号５２）上における心筋（心臓壁）の速度情報を生成する。内腔ライン速度演算部４４は、内腔ライン上の複数サンプル点について、各サンプル点ごとに心筋の速度情報を生成する。
【００４８】
内腔ライン速度演算部４４は、例えば、複数フレームに亘って得られる画像データのフレーム間において、画像データの画素値（輝度値等）に基づく相関演算等を利用したパターンマッチングにより、内腔ライン上の各サンプル点ごとに、複数フレームに亘ってそのサンプル点の移動位置を２次元平面内で追跡する。これにより、各サンプル点ごとに、２次元的な移動情報が得られ、例えば、フレーム間における移動量（移動ベクトル）とフレーム間の時間に基づいて、２次元の速度ベクトルが算出される。画像データがｘｙ直交座標系に対応したデータであれば、ｘｙ直交座標系における速度ベクトルが算出され、画像データがｒθ座標系に対応したデータであれば、ｒθ座標系における速度ベクトルが算出される。
【００４９】
また、内腔ライン速度演算部４４による内腔ライン上の各サンプル点の追跡結果は、関心領域設定部４２に送られ、関心領域設定部４２は、複数サンプル点の移動に追従するように内腔ラインの形状を変化させる、さらに、関心領域設定部４２は、始点Ｓと終点Ｅ（図２参照）をそれらに対応した各サンプル点の移動に追従させる。これにより、内腔ラインの形状の変化に応じて、つまり画像データ内における心臓の運動に追従するように、流入ライン５４と流出ライン５６が設定される。
【００５０】
なお、速度ベクトル演算部４０において、特許文献１に説明される手法により血流の２次元速度ベクトルを得る場合には、心臓壁の運動情報を利用するが、この場合には、心臓壁の運動情報として、内腔ライン速度演算部４４において算出された内腔ライン上の各サンプル点における速度ベクトルが利用される。
【００５１】
補間処理部５０と粒子生成部６０と粒子演算部７０における処理と機能については、後に図面を参照して詳述する。
【００５２】
表示処理部８０は、画像形成部２０から得られる超音波画像の画像データと、粒子演算部７０から得られる演算結果に基づいて、血流表示画像を形成する、表示処理部８０において形成された血流表示画像は表示部８２に表示される。
【００５３】
制御部１００は、図１の超音波診断装置内を全体的に制御する。制御部１００による全体的な制御には、操作デバイス９０を介して、医師や検査技師などのユーザから受け付けた指示も反映される。
【００５４】
図１に示す構成（符号を付された各部）のうち、送受信部１２，画像形成部２０，ドプラ処理部３０，速度ベクトル演算部４０，関心領域設定部４２，内腔ライン速度演算部４４，補間処理部５０，粒子生成部６０，粒子演算部７０，表示処理部８０の各部は、例えば電気電子回路やプロセッサ等のハードウェアを利用して実現することができ、その実現において必要に応じてメモリ等のデバイスが利用されてもよい。また、上記各部に対応した機能の少なくとも一部がコンピュータにより実現されてもよい。つまり、上記各部に対応した機能の少なくとも一部が、ＣＰＵやプロセッサやメモリ等のハードウェアと、ＣＰＵやプロセッサの動作を規定するソフトウェア（プログラム）との協働により実現されてもよい。
【００５５】
表示部８２の好適な具体例は、液晶ディスプレイ等であり、操作デバイス９０は、例えば、マウス、キーボード、トラックボール、タッチパネル、その他のスイッチ類等のうちの少なくとも一つにより実現できる。そして、制御部１００は、例えば、ＣＰＵやプロセッサやメモリ等のハードウェアと、ＣＰＵやプロセッサの動作を規定するソフトウェア（プログラム）との協働により実現することができる。
【００５６】
図１の超音波診断装置の全体構成は以上のとおりである。次に、図１の超音波診断装置により実現される機能の具体例について詳述する。なお、図１に示した構成（符号を付された各部）については、以下の説明において図１の符号を利用する。
【００５７】
図４は、速度ベクトル分布の具体例を示す図である。速度ベクトル演算部４０は、例えば特許文献１（特開２０１３−１９２６４３号公報）に説明されるように、血流についての超音波ビーム方向の速度情報と心臓壁の運動情報を利用して、走査面内の各位置における血流の２次元速度ベクトルを得る。具体的には、ドプラ処理部３０から得られる超音波ビーム方向の速度情報（ドプラ情報）と、内腔ライン速度演算部４４から得られる内腔ライン５２上の各サンプル点における速度情報に基づいて、例えば、図４に示す速度ベクトル分布が形成される。
【００５８】
図４に示す速度ベクトル分布は、心臓の内腔ライン５２（図２参照）を含むｘｙ座標系（直交座標系）で表現されており、ｘｙ座標系内の複数座標において算出された複数の速度ベクトル（血流の速度ベクトル）Ｖで構成されている。速度ベクトル演算部４０は、例えば、超音波の走査に対応した走査座標系、例えばビームの深さ方向ｒとビームの走査方向θによるｒθ座標系において２次元の速度ベクトル分布を形成してから、座標変換処理を行って図４に示すｘｙ座標系の速度ベクトル分布を得る。
【００５９】
速度ベクトル演算部４０は、複数のサンプル点（複数座標）に対応した速度ベクトルＶで構成される各ベクトルフレーム、つまり２次元の速度ベクトル分布を示す各ベクトルフレームを生成する。速度ベクトル演算部４０は、複数時相に亘って複数のベクトルフレームを次々に生成する。これにより、複数のベクトルフレームからなるベクトルフレーム列が得られる。
【００６０】
図５は、フレーム列を説明するための図であり、図５（Ａ）には、ベクトルフレーム列の具体例が図示されている。図５（Ａ）に示す具体例において、ベクトルフレーム列は、複数のベクトルフレーム（１〜５を代表的に図示）で構成される。
【００６１】
複数座標における速度ベクトルで構成された各ベクトルフレームは、ドプラ処理部３０から得られる超音波ビーム方向の速度情報（ドプラ情報）を利用して生成される。ドプラ情報を得るにあたり、例えばカラードプラ法の送受信を利用すると、同じビーム方向に超音波が繰り返し送受されるため、例えばＢモード画像を得る場合等に比べて、ドプラ情報を得る場合の送受信フレームレートは低くなってしまう。各ベクトルフレームを構成する速度ベクトルはドプラ情報を利用して算出されるため、ベクトルフレーム列のフレームレートも比較的低いものとなる。
【００６２】
そこで、補間処理部５０により、ベクトルフレーム列に対してフレーム間補間処理が実行される。図５（Ｂ）は、フレーム間補間処理を施された補間後のフレーム列の具体例が図示されている。図５（Ｂ）に示す具体例において、補正後のフレーム列は、複数のベクトルフレーム（図５（Ａ）に代表的に図示された１〜５）と、それらのベクトルフレーム間に追加された複数の補間フレームで構成される。各補間フレームは、互いに隣接する２つのベクトルフレーム間に速度ベクトルに基づく補間処理を適用して生成される。
【００６３】
なお、補間後のフレーム列のフレームレートは、例えば、表示部８２における表示フレームレートに応じて決定するようにしてもよい。具体的には、例えば、表示フレームレートが６０Ｈｚであれば、補間後のフレーム列のフレームレートも６０Ｈｚとなるように、追加される補間フレームの枚数等が決定される。もちろん、表示フレームレートと補間後のフレーム列のフレームレートは必ずしも一致させる必要はない。
【００６４】
図６は、フレーム間の補間処理を説明するための図である。補間処理部５０は、複数のベクトルフレームからなるベクトルフレーム列（図５参照）内において、互いに隣接する２つのベクトルフレーム間に補間処理を適用し、そのベクトルフレーム間に１又は複数の補間フレームを追加する。各補間フレームは、複数座標における補間ベクトルで構成される。
【００６５】
補間処理部５０は、各座標ごとに、互いに隣接する２つのベクトルフレームから得られる当該座標に対応した２つの速度ベクトルに基づく補間処理により、当該座標に対応した補間ベクトルを算出する。
【００６６】
図６には、座標（ｘ，ｙ）におけるフレーム間補間処理の具体例が図示されている。図６において、ベクトルフレーム（ｎ）とベクトルフレーム（ｎ＋１）は、ベクトルフレーム列内において互いに隣接する２つのベクトルフレームである。図６に示す具体例では、ベクトルフレーム（ｎ）とベクトルフレーム（ｎ＋１）の間に、等間隔で４つの補間フレーム（１）〜（４）が追加される。ベクトルフレーム（ｎ）とベクトルフレーム（ｎ＋１）の時間間隔はΔＴであり、ΔＴ内において４つの補間フレーム（１）〜（４）が等間隔Δｔで追加される。したがって、ΔＴ＝５×Δｔとなる。
【００６７】
補間処理部５０は、例えば、時間間隔に応じた線形補間により、各補間フレームを構成する補間ベクトルを算出する。例えば、ベクトルフレーム（ｎ）内の座標（ｘ，ｙ）における速度ベクトル（ｘ方向成分，ｙ方向成分）が（Ｖｘ０，Ｖｙ０）であり、ベクトルフレーム（ｎ＋１）内の座標（ｘ，ｙ）における速度ベクトル（ｘ方向成分，ｙ方向成分）が（Ｖｘ１，Ｖｙ１）である場合に、補間フレーム（１）〜（４）の各々の座標（ｘ，ｙ）における補間ベクトルのｘ方向成分とｙ方向成分は、それぞれ数１式から数４式により算出される。
【００６８】
【数１】
0005937254.tif 000002
【数２】
0005937254.tif 000003
【数３】
0005937254.tif 000004
【数４】
0005937254.tif 000005
【００６９】
なお、数１式から数４式は、時間間隔に応じた線形補間を利用した場合における一つの具体例に過ぎず、他の数式を利用して線形補間が実現されてもよいし、線形補間以外の補間処理を利用して補間ベクトルが算出されてもよい。さらに、各座標における補間ベクトルを算出する際に、その座標以外の、例えばその座標近傍の速度ベクトルが参照されてもよい。また、例えば、ベクトルフレーム（ｎ）とベクトルフレーム（ｎ＋１）の間における補間処理に、これら２つのベクトルフレーム以外の、例えばこれら２つのベクトルフレームの近傍のベクトルフレーム内における速度ベクトルが利用されてもよい。
【００７０】
補間処理部５０は、速度ベクトル分布が得られた複数座標について、各座標ごとにフレーム間補間処理を実行して補間ベクトルを得ることにより、複数座標における補間ベクトルで構成された各補間フレームを形成する。こうして、複数のベクトルフレームとそれらのベクトルフレーム間に追加された複数の補間フレームで構成される補間後のフレーム列（図５（Ｂ）参照）が得られる。そして、補間後のフレーム列に基づいて、血流に関する仮想的な複数粒子の移動先が算出される。仮想的な複数粒子は、粒子生成部６０により生成される。
【００７１】
図７は、複数粒子の生成を説明するための図である。粒子生成部６０は、血流が含まれる座標系内、つまり超音波画像と２次元速度ベクトル分布が形成された座標系内に、血流に関する複数粒子を設定する。粒子生成部６０は、例えば、超音波画像内において心臓に対して設定された流入ライン５４上に複数粒子を生成する。この場合には、流入ライン５４が複数粒子を生成する生成ラインとなる。
【００７２】
例えば、図７において基本例（１）として示すように、始点Ｓと終点Ｅを直線で結んだ流入ライン５４上に、等間隔で一例に並ぶ複数粒子が生成される。例えば５０個の粒子が流入ライン５４上に等間隔で生成される。なお、流入ライン５４の長さが５０画素（ピクセル）以下の場合には、流入ライン５４上の１画素につき１個の粒子が生成される。もちろん、５０個以外の設定個数で複数粒子が生成されてもよい。ユーザが複数粒子の個数を設定又は変更できるようにしてもよい。
【００７３】
粒子生成部６０は、特定の１フレームのみ（１時相のみ）で複数粒子を発生させてもよいが、複数フレームに亘って定期的に各フレーム内に複数粒子を発生させることが望ましい。例えば、複数のベクトルフレームからなるベクトルフレーム列（図５（Ａ））において、各ベクトルフレームごとに複数粒子を発生させてもよいし、補正後のフレーム列（図５（Ｂ））において、各ベクトルフレーム内と各補間フレーム内に複数粒子を発生させてもよい。もちろん、数フレーム間隔で各フレーム内に複数粒子が生成されてもよい。
【００７４】
関心領域設定部４２により、超音波画像の画像データ内における心臓の運動に追従するように流入ライン５４が設定されていれば、心臓の運動に追従するように位置を修正しつつ複数粒子を発生させることができる。
【００７５】
また、図７において変形例（２）として示すように、流入ライン５４上における速度ベクトルＶの大きさや向きに応じて、複数粒子の発生態様を異ならせてもよい。変形例（２）では、速度ベクトルＶの大きさ（例えば流入ライン５４に対する垂直成分）が大きいほど、その速度ベクトルＶの位置における粒子個数を増大させている。また、変形例（２）では、速度ベクトルＶが正方向の場合にのみ、その速度ベクトルＶの位置に粒子が生成されている。粒子を生成させる向き（正方向または負方向）は、例えばユーザにより設定されてもよい。
【００７６】
また、図７に示す具体例に限らず、粒子生成部６０は、例えば、曲線状の流入ライン５４（図３参照）上に複数粒子を生成してもよい。さらに、流入ライン５４上に限らず、例えばユーザにより指定された生成ライン上または生成領域内に複数粒子が生成されてもよいし、関心領域設定部４２が設定した関心領域内において例えば２次元的に均等に（例えば格子状に）複数粒子が生成されてもよい。
【００７７】
粒子生成部６０により複数粒子が生成されると、粒子演算部７０は、補正後のフレーム列（図５（Ｂ））に基づいて、各粒子ごとにその粒子の移動先を算出する。
【００７８】
図８は、各粒子の移動先の演算の具体例を説明するための図である。図８には、ベクトルフレーム（ｎ）とベクトルフレーム（ｎ＋１）の間に、等間隔で４つの補間フレーム（１）〜（４）が追加された補正後のフレーム列（図６参照）に基づく、一つの粒子Ｐに関する演算の具体例が示されている。
【００７９】
ベクトルフレーム（ｎ）内において粒子Ｐが位置Ｐ０（座標Ｐ０）に存在する場合、例えばベクトルフレーム（ｎ）内の位置Ｐ０に粒子Ｐが生成された場合、まず、位置Ｐ０における血流の速度ベクトルＶ０が利用される。ベクトルフレーム（ｎ）を構成する複数座標の速度ベクトルの中に、位置Ｐ０（座標Ｐ０）に対応する座標の速度ベクトルが存在すれば当該速度ベクトルが速度ベクトルＶ０とされる。もし、位置Ｐ０に対応する座標の速度ベクトルがなければ、位置Ｐ０の近傍にある複数座標の速度ベクトルに基づく線形補間処理（フレーム内補間処理）等により速度ベクトルＶ０が算出される。
【００８０】
そして、速度ベクトルＶ０とフレーム間隔Δｔ（図６参照）を乗算することにより移動ベクトル（大きさが速度ベクトルＶ０のΔｔ倍で速度ベクトルＶ０と同じ方向）が算出され、位置Ｐ０からその移動ベクトルだけ移動した位置Ｐ１（座標Ｐ１）が導出される。これにより得られた位置Ｐ１が、ベクトルフレーム（ｎ）の次のフレーム（次の時相）である補間フレーム（１）における粒子Ｐの位置（移動先の座標）となる。
【００８１】
次に、位置Ｐ１における血流の速度ベクトルＶ１が利用される。補間フレーム（１）を構成する複数座標の補間ベクトルの中に、位置Ｐ１（座標Ｐ１）に対応する座標の補間ベクトルが存在すれば当該補間ベクトルが速度ベクトルＶ１とされ、位置Ｐ１に対応する座標の補間ベクトルがなければ、位置Ｐ１の近傍にある複数座標の補間ベクトルに基づく線形補間処理（フレーム内補間処理）等により速度ベクトルＶ１が算出される。
【００８２】
そして、速度ベクトルＶ１とフレーム間隔Δｔを乗算することにより移動ベクトル（大きさが速度ベクトルＶ１のΔｔ倍で速度ベクトルＶ１と同じ方向）が算出され、位置Ｐ１からその移動ベクトルだけ移動した位置Ｐ２（座標Ｐ２）が導出される。これにより得られた位置Ｐ２が、補間フレーム（１）の次のフレーム（次の時相）である補間フレーム（２）における粒子Ｐの位置（移動先の座標）となる。
【００８３】
補間フレーム（２）に続く補間フレーム（３）と補間フレーム（４）についても、上記と同様な処理が実行される。つまり、粒子Ｐの位置における速度ベクトル（Ｖ２，Ｖ３）とフレーム間隔Δｔを乗算することにより得られる移動ベクトルに基づいて粒子Ｐの移動先の座標が算出される。図８の具体例において、補間フレーム（３）における位置Ｐ３と補間フレーム（４）における位置Ｐ４が粒子Ｐの移動先の座標である。
【００８４】
さらに、位置Ｐ４における血流の速度ベクトルＶ４が利用される。補間フレーム（４）を構成する複数座標の補間ベクトルの中に、位置Ｐ４（座標Ｐ４）に対応する座標の補間ベクトルが存在すれば当該補間ベクトルが速度ベクトルＶ４とされ、位置Ｐ４に対応する座標の補間ベクトルがなければ、位置Ｐ４の近傍にある複数座標の補間ベクトルに基づく線形補間処理（フレーム内補間処理）等により速度ベクトルＶ４が算出される。
【００８５】
そして、速度ベクトルＶ４とフレーム間隔Δｔを乗算することにより移動ベクトル（大きさが速度ベクトルＶ４のΔｔ倍で速度ベクトルＶ４と同じ方向）が算出され、位置Ｐ４からその移動ベクトルだけ移動した位置Ｐ５（座標Ｐ５）が導出される。これにより得られた位置Ｐ５が、補間フレーム（５）の次のフレーム（次の時相）であるベクトルフレーム（ｎ＋１）における粒子Ｐの位置（移動先の座標）となる。
【００８６】
こうして、粒子演算部７０は、ベクトルフレーム（ｎ＋１）の以降に続く複数フレーム（補間フレームまたはベクトルフレーム）においても、上記と同様な処理を実行し、後に説明する終了条件が満たされるまで、次々に粒子Ｐの移動先を導出する。また、粒子演算部７０は、粒子生成部６０が生成した複数粒子の各々について、各粒子が生成されたフレーム（時相）から、その粒子の移動先を導出する。
【００８７】
図８に示す具体例によれば、ベクトルフレーム間に複数の補間フレームが追加された補間後のフレーム列に基づいて各粒子の移動先が導出されるため、複数の補間フレームを追加しない場合に比べて、移動先の推定精度が高められる。
【００８８】
例えば、図８において補間フレーム（１）〜（４）を利用せずに、ベクトルフレーム（ｎ）における粒子Ｐの位置Ｐ０から、速度ベクトルＶ０にベクトルフレーム間隔ΔＴ（図６参照）を乗算して得られる移動ベクトル（大きさが速度ベクトルＶ０のΔＴ倍で速度ベクトルＶ０と同じ方向）だけ移動した位置をベクトルフレーム（ｎ＋１）における粒子Ｐの移動先としてしまうと、図８の具体例で得られる位置Ｐ５（座標Ｐ５）とは大きく異なる移動先が得られてしまう。補間フレーム（１）〜（４）を利用しないと、ベクトルフレーム（ｎ）とベクトルフレーム（ｎ＋１）との間における粒子Ｐの速度ベクトルの変化が反映されないためである。
【００８９】
図９は、各粒子の移動先の演算の例外処理を説明するための図である。粒子演算部７０は、図８を利用して説明した基本処理により各粒子の移動先を導出するが、図９に示すように、各粒子の移動先が内腔ライン５２を超えてしまう場合には、各粒子の移動先を内腔ライン５２上または内腔ライン５２の内側（心腔側）近傍に移動先を修正する。
【００９０】
例えば、図９に示す具体例のように、各フレーム（ベクトルフレーム又は補間フレーム）において粒子Ｐが位置Ｐ６にあり、基本処理により得られる次フレーム（ベクトルフレーム又は補間フレーム）における移動先が位置Ｐ７となる場合、つまり移動ベクトル（破線矢印）と内腔ライン５２が交差する場合には、次フレームにおける移動先が位置Ｐ７から位置Ｐ７´に修正される。図９の具体例において、位置Ｐ７´は、内腔ライン５２と移動ベクトルの交点であるが、当該交点の近傍に、例えば内腔ライン５２の内側（心腔側）の近傍が位置Ｐ７´とされてもよい。
【００９１】
図１０は、各粒子の移動先の演算の終了条件を説明するための図である。粒子演算部７０は、図８を利用して説明した基本処理と図９を利用して説明した例外処理により各粒子の移動先を次々に導出し、図１０に示すように、各粒子の移動先が流出ライン５６を通過した場合に、その粒子の移動先の演算を終了する。
【００９２】
例えば、図１０に示す具体例のように、各フレーム（ベクトルフレーム又は補間フレーム）において粒子Ｐが位置Ｐ８にあり、次フレーム（ベクトルフレーム又は補間フレーム）における移動先が位置Ｐ９となる場合、つまり移動ベクトル（破線矢印）と流出ライン５６が交差する場合に、粒子演算部７０は、次フレームにおける移動先を位置Ｐ９から位置Ｐ９´に修正して、粒子Ｐの移動先の演算を終了する。
【００９３】
なお、図１０の具体例において、位置Ｐ９´は、流出ライン５６と移動ベクトルの交点であるが、当該交点の近傍、例えば流出ライン５６の上側（心腔側）の近傍が位置Ｐ９´とされてもよい。
【００９４】
粒子生成部６０により複数粒子が生成され、粒子演算部７０により各粒子ごとにその粒子の複数時相（複数フレーム）に亘る移動先が次々に算出されると、表示処理部８０は、複数時相に亘る各粒子の移動先の座標を画像内に示した血流表示画像を形成する。表示処理部８０は、例えば、複数時相に亘る各粒子の移動先の座標を点と線の少なくとも一方による軌跡で示した流跡線の画像を形成する。
【００９５】
図１１は、流跡線の具体例を示す図である。図１１には、複数粒子のうちの一つである粒子Ｐに関する１本の流跡線Ｌが図示されている。流跡線Ｌは、補正後のフレーム列（図５（Ｂ），図６）内において複数時相（複数フレーム）に亘って次々に算出された粒子Ｐの移動先（複数時相に対応した複数位置）の軌跡である。例えば、粒子Ｐの複数時相における移動先（例えば図８の位置Ｐ０，Ｐ１，Ｐ２，Ｐ３，Ｐ４，Ｐ５，・・・）を時相順に直線又は曲線（例えばスプライン補間に基づく曲線）で結ぶことにより流跡線Ｌが形成される。なお、流跡線Ｌ上に又は流跡線Ｌに代えて、粒子Ｐの複数時相における移動先（例えば図８の位置Ｐ０，Ｐ１，Ｐ２，Ｐ３，Ｐ４，Ｐ５，・・・）が複数の点による列で表示されてもよい。
【００９６】
また、表示される流跡線Ｌの長さは適宜に調整されることが望ましい。例えば表示処理部８０は、各粒子の軌跡が表示される時相（現時相）から所定時間前の時相（残存時相）までの流跡線Ｌを形成する。つまり、図１１に示す具体例のように、粒子Ｐの流跡線Ｌとして、現時相における粒子Ｐの位置Ｐｎから、現時相よりも所定時間前の残存時相までの流跡線Ｌの部分（実線部分）のみが表示され、残存時相よりも前に形成された流跡線Ｌの部分（破線部分）が表示されない。さらに、ユーザが所定時間を設定できることが望ましい。例えば操作デバイス９０が備えるキーをユーザが操作することにより、１００ｍｓ〜１０００ｍｓの設定可能範囲内において１００ｍｓ単位で所定時間をユーザが決定できるようにしてもよい。
【００９７】
なお、表示処理部８０は、例えば、各粒子が生成された時相から所定時間後の時相までの流跡線Ｌの部分のみを表示するようにしてもよい。
【００９８】
図１２は、血流表示画像の具体例を示す図である。表示処理部８０は、複数粒子について各粒子ごとに流跡線Ｌを形成し、例えば、画像形成部２０から得られる心臓の超音波画像上に複数粒子の流跡線Ｌを示した血流表示画像を形成する。また、ドプラ処理部３０から得られるドプラ情報を利用して形成されるカラードプラ画像上に複数粒子の流跡線Ｌを示した血流表示画像が形成されてもよい。表示処理部８０において形成された血流表示画像は表示部８２に表示される。
【００９９】
また、生成された全粒子に対応した複数の流跡線Ｌを全て表示せずに、流跡線Ｌの表示本数を間引いてもよい。例えば形成された多数の流跡線Ｌの中で１０本につき１本のみを表示させるなどにより、血流表示画像内において多数の流跡線Ｌが混雑することを解消して見易い画像が形成されてもよい。さらに、例えば、操作デバイス９０が備えるキーをユーザが操作することにより、表示させる流跡線Ｌの本数又は間引きの割合をユーザが設定できるようにしてもよい。
【０１００】
表示処理部８０は、複数の表示時相に亘って、各表示時相ごとに血流表示画像を形成する。例えば、各表示時相に対応した心臓の超音波画像又はカラードプラ画像上に、その表示時相に対応した血流表示画像（図１２）が形成される。これにより、複数時相に亘って動的に変化する心臓の動きを超音波画像又はカラードプラ画像で確認しつつ、複数時相に亘って複数粒子の流跡線Ｌが変化する様子を視覚的に且つ動的に確認できるようになる。もちろん、ユーザが望む特定の心拍時相（拡張末期や収縮末期など）における静止画（フリーズ画像）を表示させてもよい。これにより、例えば、心臓内における血液の渦流、乱流、停滞などの状態を視覚的且つ直感的に確認できるようになる。
【０１０１】
以上、本発明の好適な実施形態を説明したが、上述した実施形態は、あらゆる点で単なる例示にすぎず、本発明の範囲を限定するものではない。本発明は、その本質を逸脱しない範囲で各種の変形形態を包含する。
【符号の説明】
【０１０２】
１０  プローブ、１２  送受信部、２０  画像形成部、３０  ドプラ処理部、４０  速度ベクトル演算部、５０  補間処理部、６０  粒子生成部、７０  粒子演算部、８０  表示処理部、１００  制御部。
【要約】
【課題】超音波を利用して得られる血流の移動情報の精度を高める。
【解決手段】速度ベクトル演算部４０は、超音波を送受して得られた信号に基づいて血流の各座標における速度ベクトルを得ることにより、複数座標の速度ベクトルで構成された各ベクトルフレームからなる複数のベクトルフレームを生成する。補間処理部５０は、互いに隣接する２つのベクトルフレーム間に補間処理を適用することにより各補間フレームを生成し、当該ベクトルフレーム間に１又は複数の補間フレームを追加する。補間処理部５０は、互いに隣接する２つのベクトルフレーム間における速度ベクトルに基づく補間処理により、複数座標の補間ベクトルで構成された各補間フレームを生成する。
【選択図】図１【図１】
0005937254.tif 000006
【図２】
0005937254.tif 000007
【図３】
0005937254.tif 000008
【図４】
0005937254.tif 000009
【図５】
0005937254.tif 000010
【図６】
0005937254.tif 000011
【図７】
0005937254.tif 000012
【図８】
0005937254.tif 000013
【図９】
0005937254.tif 000014
【図１０】
0005937254.tif 000015
【図１１】
0005937254.tif 000016
【図１２】
0005937254.tif 000017