(57)【要約】【課題】生体内における対象組織を精度良く抽出する。【解決手段】ボリュームデータに対して複数の参照断面からなる参照断面列５６が設定される。参照断面列５６は一定間隔で配列された複数のマニュアルトレース用参照断面５８とそれ以外の複数の自動トレース用参照断面６０とで構成される。参照断面５８についてはマニュアルトレース処理及び自動的な修正処理が適用され、参照断面６０については補間トレースラインの生成処理及び自動的な修正処理が適用される。最終的に複数のトレースラインに基づいて対象組織が抽出されあるいは体積が演算される。隣接するマニュアルトレースライン間における形状の類似度が小さい場合に、そこに１又は複数マニュアルトレース用参照断面を追加設定するようにしてもよい。【選択図】図８
2008142519.tif 000002
【特許請求の範囲】
【請求項１】
対象組織を含む三次元空間に対して超音波の送受波を行ってボリュームデータを取得するデータ取得手段と、
前記対象組織に対して設定された複数の代表断面上に複数のマニュアルトレースラインが形成された場合に、各マニュアルトレースラインに対して第１修正処理を適用する代表断面処理手段と、
前記第１修正処理後の複数のマニュアルトレースラインに基づく補間処理により前記対象組織に対して設定された複数の非代表断面上に複数の補間トレースラインを形成すると共に、各補間トレースラインに対して第２修正処理を適用する非代表断面処理手段と、
前記第１修正処理後の複数のマニュアルトレースライン及び前記第２修正処理後の複数の補間トレースラインからなるトレースライン列に基づいて、前記ボリュームデータから対象組織データを抽出し、あるいは、前記対象組織の体積を演算する手段と、
を含むことを特徴とする超音波診断装置。
【請求項２】
請求項１記載の装置において、
前記代表断面処理手段は、
前記対象組織に対して設定された断面セットの中から前記複数の代表断面を選択する手段と、
前記選択された各代表断面に対するユーザーのトレース入力を受け付ける手段と、
前記第１修正処理に当たって、前記各マニュアルトレースライン上における各ポイントごとに修正可否を判断し、修正可能なポイントについて実際の組織輪郭に基づいて当該ポイントの位置を修正する第１修正処理手段と、
を含むことを特徴とする超音波診断装置。
【請求項３】
請求項２記載の装置において、
前記第１修正処理手段は、前記各ポイントごとに交差ラインを設定してその交差ライン上でエッジ検出を行うことにより修正の可否を判断する、ことを特徴とする超音波診断装置。
【請求項４】
請求項１記載の装置において、
前記非代表断面処理手段は、
前記第１修正処理後の複数のマニュアルトレースラインに基づいて補間処理により前記各非代表断面上に補間トレースラインを形成する手段と、
前記第２修正処理に当たって、前記各補間トレースライン上における各ポイントごとに修正可否を判断し、修正可能なポイントについて実際の組織輪郭に基づいて当該ポイントの位置を修正する第２修正処理手段と、
を含むことを特徴とする超音波診断装置。
【請求項５】
請求項４記載の装置において、
前記第２修正処理手段は、前記各ポイントごとに交差ラインを設定してその交差ライン上でエッジ検出を行うことにより修正の可否を判断する、ことを特徴とする超音波診断装置。
【請求項６】
請求項１記載の装置において、
前記対象組織を通過する基準ラインの指定を受け付ける手段と、
前記基準ラインの方向に並んで前記基準ラインに直交する断面セットを設定する手段と、
を含み、
前記断面セットの中から前記複数の代表断面及び前記複数の非代表断面が定められる、ことを特徴とする超音波診断装置。
【請求項７】
対象組織を含む三次元空間に対して超音波の送受波を行って得られたボリュームデータを格納し、
前記対象組織に対して設定された複数の代表断面に対する複数のマニュアルトレースを受け付け、
前記各代表断面上に形成された各マニュアルトレースライン、及び、前記対象組織に対して設定された各非代表断面上に形成された各補間トレースライン、における各ポイントについて実際の組織輪郭に基づいて位置修正の可否を判断し、修正可能なポイントの位置を修正する修正処理を実行し、
前記修正処理後の複数のマニュアルトレースライン及び複数の補間トレースラインからなるトレースライン列に基づいて対象組織を抽出し、あるいは、対象組織の体積を演算する、
ことを特徴とするボリュームデータ処理方法。
【請求項８】
対象組織を含む三次元空間に対して超音波の送受波を行ってボリュームデータを取得するデータ取得手段と、
前記対象組織に対して設定された複数の代表断面上に複数のマニュアルトレースラインが形成された場合に、各マニュアルトレースラインに対して第１修正処理を適用する代表断面処理手段と、
前記第１修正処理前又は後の複数のマニュアルトレースラインにおける隣接マニュアルトレースラインペアごとに類似度を演算する類似度演算手段と、
前記隣接マニュアルトレースラインペアごとの類似度に基づいて前記複数の代表断面の他に１又は複数の追加代表断面を加える追加代表断面設定手段と、
前記１又は複数の追加代表断面上に形成された１又は複数のマニュアルトレースラインに対して第１修正処理を適用する追加代表断面処理手段と、
前記第１修正処理後の複数のマニュアルトレースラインに基づく補間処理により前記対象組織に対して設定された複数の非代表断面上に複数の補間トレースラインを形成すると共に、各補間トレースラインに対して第２修正処理を適用する非代表断面処理手段と、
前記第１修正処理後の複数のマニュアルトレースライン及び前記第２修正処理後の複数の補間トレースラインからなるトレースライン列に基づいて、前記ボリュームデータから対象組織データを抽出し、あるいは、前記対象組織の体積を演算する手段と、
を含むことを特徴とする超音波診断装置。
【請求項９】
請求項８記載の装置において、
前記類似度演算手段は、前記隣接マニュアルトレースラインペアごとにそれを構成する２つのマニュアルトレースライン間で相互相関演算を実行することにより前記類似度を演算する、ことを特徴とする超音波診断装置。
【請求項１０】
請求項８記載の装置において、
前記追加代表断面設定手段は、前記類似度が所定の追加条件を満たした隣接代表断面ペア間に前記１又は複数の追加代表断面を加える、ことを特徴とする超音波診断装置。
【請求項１１】
請求項１０記載の装置において、
前記追加断面設定手段は、前記各類似度の大きさに応じて追加代表断面の追加数を決定する、ことを特徴とする超音波診断装置。
【請求項１２】
対象組織を含む三次元空間に対して超音波の送受波を行って得られたボリュームデータを格納し、
前記対象組織に対して設定された複数の代表断面に対する複数のマニュアルトレースを受け付け、
前記複数の代表断面に加えて１又は複数の追加代表断面が設定された場合に、１又は複数の追加代表断面に対する１又は複数のマニュアルトレースを受け付け、
前記各代表断面上に形成された各マニュアルトレースライン、及び、前記対象組織に対して設定された各非代表断面上に形成された各補間トレースライン、における各ポイントについて実際の組織輪郭に基づいて位置修正の可否を判断し、修正可能なポイントの位置を修正する修正処理を実行し、
前記修正処理後の複数のマニュアルトレースライン及び複数の補間トレースラインからなるトレースライン列に基づいて対象組織を抽出し、あるいは、対象組織の体積を演算する、
ことを特徴とするボリュームデータ処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は超音波診断装置及びボリュームデータ処理方法に関し、特に三次元空間内での対象組織の特定又は計測に関する。
【背景技術】
【０００２】
超音波ビームを走査して形成される走査面をそれと直交する方向へ走査することにより三次元空間（三次元エコーデータ取込空間）を形成できる。三次元空間内に存在する対象組織を抽出してその体積を演算する方法として従来からDisk Summation法が周知である。この従来方法では、対象組織を横切って整列する複数のスライスデータのそれぞれが断層画像として表示され、各断層画像上で対象組織の輪郭がマニュアルでトレースされ、これにより各スライスデータごとに断面積が求まる。各スライスデータごとに断面積と厚み（断面間隔）を利用して対象組織部分の要素体積が演算される。それらの要素体積を加算することによって対象組織の体積が求まる。この従来手法ではスライスデータの個数（つまり計測対象となる断面の個数）を増加させないと、精度良く体積を求められない。一方、多くのスライスデータをマニュアルトレースに委ねると、ユーザーの負担が大きい。マニュアルトレースではなく、各断面に対して自動トレース処理を適用することも考えられる。しかし、輪郭が不明瞭な部分があると、トレースが行えないかトレース結果の信頼性が低下してしまう。
【０００３】
特許文献１には三次元空間内に存在する組織（心腔等）の抽出処理が記載されている。特許文献２には輪郭モデルを利用して組織を抽出することが記載されている。特許文献３にはマニュアルトレース結果を補正する技術が記載されている。特許文献４−７には組織の体積を演算する技術が記載されている。
【０００４】
【特許文献１】特開２００４−１５９９９７号公報
【特許文献２】特開２００２−２２４１１６号公報
【特許文献３】特開平１１−１６４８３４号公報
【特許文献４】特開平１０−３０５０３３号公報
【特許文献５】特開２００２−３３０９６８号公報
【特許文献６】特開２０００−１０７１８３号公報
【特許文献７】特開２０００−２９６１２９号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
対象組織の特定又は計測に際し、高い精度を維持しつつユーザーの負担を軽減することが望まれるが、それに十分に応えるものは未だ実現されていない。上記したいずれの文献にも、組織輪郭が不明瞭な部位と明瞭な部位とが混在する場合において、マニュアルトレースと自動トレースのそれぞれの利点を発揮させるための構成は記載されていない。例えば、妊娠後期における子宮内の胎盤に関し、胎盤周囲において子宮に接している領域部分については輝度差が乏しく、羊水に接している領域部分については輝度差が大きい。そのような状況において胎盤のトレースに当たってユーザーの負担を軽減し、且つ、トレース精度を高めたいとの要請がある。
【０００６】
本発明の目的は、ユーザーにあまり負担をかけずに対象組織の特定精度あるいは計測精度を高められるようにすることにある。
【課題を解決するための手段】
【０００７】
本発明は、対象組織を含む三次元空間に対して超音波の送受波を行ってボリュームデータを取得するデータ取得手段と、前記対象組織に対して設定された複数の代表断面上に複数のマニュアルトレースラインが形成された場合に、各マニュアルトレースラインに対して第１修正処理を適用する代表断面処理手段と、前記第１修正処理後の複数のマニュアルトレースラインに基づく補間処理により前記対象組織に対して設定された複数の非代表断面上に複数の補間トレースラインを形成すると共に、各補間トレースラインに対して第２修正処理を適用する非代表断面処理手段と、前記第１修正処理後の複数のマニュアルトレースライン及び前記第２修正処理後の複数の補間トレースラインからなるトレースライン列に基づいて、前記ボリュームデータから対象組織データを抽出し、あるいは、前記対象組織の体積を演算する手段と、を含むことを特徴とする超音波診断装置に関するものである。
【０００８】
望ましくは、ボリュームデータ内における対象組織に対して複数の断面からなる断面セットが設定される。この場合には対象組織の全体がカバーされるように断面セットが設定されるのが望ましい。断面セットは走査面セットに相当していてもよいし、走査面セットとは別に設定されてもよい。断面セットにおいて個々の断面間隔は均等であってもよいし、非均等であってもよい。形態が大きく変化するような部位に対しては断面密度を高めるようにしてもよい。いずれにしても、断面セットは、複数の代表断面及び複数の非代表断面を含む。対象組織に対して複数の代表断面が定められると、各代表断面上においてマニュアルトレースラインが形成される。本発明によれば、全断面に対してマニュアルトレースを行うのではなく、代表断面に対してマニュアルトレースを行えばよいので、ユーザーの負担を軽減できる。各マニュアルトレースラインに対して第１修正処理が適用される。第１修正処理は、修正要否の判断と修正の実行とを含む概念である。各マニュアルトレースラインについては、少なくとも修正要否の判断がなされる。このように、各マニュアルトレースラインは必要に応じて修正される。例えば、対象組織の輪郭（境界）が明瞭な部分については自動トレースが可能であり、しかも自動トレースの方がマニュアルトレースよりも一般に精度が良いので（細かい凹凸に対応し易い）、マニュアルトレースの結果が全体的に又は部分的に自動的に修正される（修正済みマニュアルトレースラインの形成）。この場合、例えば対象組織の輪郭の明瞭性を修正処理適用の要否を判断する基準とすることができる。複数の代表断面上にそれぞれ第１修正処理が適用されたマニュアルトレースラインが形成されると、それらを基準として、個々の非代表断面上に補間トレースラインを形成することが容易となる。例えば、隣接する２つの代表断面（隣接代表断面ペア）を構成する２つのマニュアルトレースライン（隣接マニュアルトレースライン）の間に、つまり１又は複数の非代表断面上に、補間処理によって１又は複数の補間トレースラインが形成される。各補間トレースラインに対して第２修正処理が適用される。第２修正処理は、上記第１修正処理と同様に、修正要否の判断と修正の実行とを含む概念である。つまり、各補間トレースラインは必要に応じて実際の組織輪郭に基づいて修正され、これにより修正済み補間トレースラインが形成される。複数の断面上に形成されたトレースライン列（第１修正処理後の複数のマニュアルトレースライン及び第２修正処理後の複数の補間トレースラインを含む）に基づいて、対象組織の表面形状を特定可能であり、これにより対象組織の表面を表す三次元画像が形成され、あるいは、対象組織の体積が演算される。個々のマニュアルトレースライン及び個々の補間トレースラインを実際の組織輪郭に基づいて（且つ、必要に応じて）修正可能であるので演算計測の精度を高められる。換言すれば、上記方法によれば、マニュアルトレースと自動トレースを有機的に組み合わせて両者の利点を効果的に発揮させることができる。
【０００９】
望ましくは、前記代表断面処理手段は、前記対象組織に対して設定された断面セットの中から前記複数の代表断面を選択する手段と、前記選択された各代表断面に対するユーザーのトレース入力を受け付ける手段と、前記第１修正処理に当たって、前記各マニュアルトレースライン上における各ポイントごとに修正可否を判断し、修正可能なポイントについて実際の組織輪郭に基づいて当該ポイントの位置を修正する第１修正処理手段と、を含む。代表断面は自動的に指定され、あるいは、ユーザーにより指定される。断面セットを構成する断面の個数あるいは代表断面の個数は、自動的に設定され、あるいはユーザーにより指定される。ユーザーによる組織輪郭のトレースは通常断層画像上で行われる。マニュアルトレースラインを描く場合、どうしても細かい凹凸を十分にトレースできない場合があるが、組織輪郭が明瞭である（あるいは識別できる）限りにおいて、マニュアルトレース結果を自動的に修正すれば（オートトレースの重畳適用）、上記のような細かい凹凸を高精度にトレースできる。あるいは、そのような細かい凹凸に対して忠実にマニュアルトレースを行う必要がないので、ユーザーの負担を軽減できる。
【００１０】
望ましくは、前記第１修正処理手段は、前記各ポイントごとに交差ラインを設定してその交差ライン上でエッジ検出を行うことにより修正の可否を判断する。エッジ検出に当たっては輝度値の変化（微分値）等を参照するのが望ましい。
【００１１】
望ましくは、前記非代表断面処理手段は、前記第１修正処理後の複数のマニュアルトレースラインに基づいて補間処理により前記各非代表断面上に補間トレースラインを形成する手段と、前記第２修正処理に当たって、前記各補間トレースライン上における各ポイントごとに修正可否を判断し、修正可能なポイントについて実際の組織輪郭に基づいて当該ポイントの位置を修正する第２修正処理手段と、を含む。この構成によれば、個々の非代表断面上に自動的に補間トレースラインを生成した上で、各補間トレースラインを実際の組織輪郭に基づいて修正して高精度のトレース結果を得られる。非代表断面についてはマニュアルでの作業が不要となるので、ユーザーの負担を軽減できる。また、単に補間トレースラインを生成するのではなく、それを実際の組織輪郭に基づいて修正し得るので補間処理結果の信頼性を高められる。
【００１２】
望ましくは、前記第２修正処理手段は、前記各ポイントごとに交差ラインを設定してその交差ライン上でエッジ検出を行うことにより修正の可否を判断する。望ましくは、前記対象組織を通過する基準ラインの指定を受け付ける手段と、前記基準ラインの方向に並んで前記基準ラインに直交する断面セットを設定する手段と、を含み、前記断面セットの中から前記複数の代表断面及び前記複数の非代表断面が定められる。この構成によれば、断面セットはトレース前に予め設定され、それを構成する一部の断面がそれぞれ代表断面になり、他の一部の断面がそれぞれ非代表断面になる。但し、複数の代表断面及び複数の非代表断面を段階的に設定することによって結果として断面セットが構成されてもよい。
【００１３】
また、本発明は、対象組織を含む三次元空間に対して超音波の送受波を行って得られたボリュームデータを格納し、前記対象組織に対して設定された複数の代表断面に対する複数のマニュアルトレースを受け付け、前記各代表断面上に形成された各マニュアルトレースライン、及び、前記対象組織に対して設定された各非代表断面上に形成された各補間トレースライン、における各ポイントについて実際の組織輪郭に基づいて位置修正の可否を判断し、修正可能なポイントの位置を修正する修正処理を実行し、前記修正処理後の複数のマニュアルトレースライン及び複数の補間トレースラインからなるトレースライン列に基づいて対象組織を抽出し、あるいは、対象組織の体積を演算する、ことを特徴とするボリュームデータ処理方法に関する。
【００１４】
また、本発明は、対象組織を含む三次元空間に対して超音波の送受波を行ってボリュームデータを取得するデータ取得手段と、前記対象組織に対して設定された複数の代表断面上に複数のマニュアルトレースラインが形成された場合に、各マニュアルトレースラインに対して第１修正処理を適用する代表断面処理手段と、前記第１修正処理前又は後の複数のマニュアルトレースラインにおける隣接マニュアルトレースラインペアごとに類似度を演算する類似度演算手段と、前記隣接マニュアルトレースラインペアごとの類似度に基づいて前記複数の代表断面の他に１又は複数の追加代表断面を加える追加代表断面設定手段と、前記１又は複数の追加代表断面上に形成された１又は複数のマニュアルトレースラインに対して第１修正処理を適用する追加代表断面処理手段と、前記第１修正処理後の複数のマニュアルトレースラインに基づく補間処理により前記対象組織に対して設定された複数の非代表断面上に複数の補間トレースラインを形成すると共に、各補間トレースラインに対して第２修正処理を適用する非代表断面処理手段と、前記第１修正処理後の複数のマニュアルトレースライン及び前記第２修正処理後の複数の補間トレースラインからなるトレースライン列に基づいて、前記ボリュームデータから対象組織データを抽出し、あるいは、前記対象組織の体積を演算する手段と、を含むことを特徴とする超音波診断装置に関するものである。
【００１５】
上記構成によれば、マニュアルトレース、補間処理及び自動トレースを組み合わせて対象組織の計測を行える。特に、隣接代表断面ペア間において必要に応じて１又は複数の追加代表断面が事後的に設定され、つまり、隣接マニュアルトレースラインペア間における類似度が小さい場合に（組織形状の変化が大となる部分に対して）、その隣接マニュアルトレースラインペア間に１又は複数のマニュアルトレースラインを追加的に生成可能である。これにより、例えば、対象組織において形状変化が大きい部分には代表断面が密に設定されるので（つまり、そこにマニュアルトレースラインが密に形成されるので）、補間処理と相俟って、対象組織の計測精度を向上できる。勿論、そのような追加設定の要否は類似度の大きさに応じて決定されるので、代表断面（当初の代表断面及び追加の代表断面）の個数を必要最小限にして、マニュアルトレースの負担が著しく増大することを防止できる。上記構成では、複数の当初マニュアルトレースラインと１又は複数の追加マニュアルトレースラインとが生成される。類似度の演算は、第１修正処理前の複数の当初マニュアルトレースラインに対して行われてもよいし、第１修正処理後の複数の当初マニュアルトレースラインに対して行われてもよい。前者の場合には、追加処理後において、すべてのマニュアルトレースラインに対して第１修正処理が適用されることになる。前者の場合には、複数の当初マニュアルトレースラインに対する１回目の第１修正処理と、複数の追加マニュアルトレースラインに対する２回目の第１修正処理と、が段階的に適用される。
【００１６】
望ましくは、前記類似度演算手段は、前記隣接マニュアルトレースラインペアごとにそれを構成する２つのマニュアルトレースライン間で相互相関演算を実行することにより前記類似度を演算する。相互相関演算によれば、隣り合うマニュアルトレースライン間において形状がどの程度似ているのかを表す情報を得られる。
【００１７】
望ましくは、前記追加代表断面設定手段は、前記類似度が所定の追加条件を満たした隣接代表断面ペア間に前記１又は複数の追加代表断面を加える。望ましくは、前記追加断面設定手段は、前記各類似度の大きさに応じて追加代表断面の追加数を決定する。
【００１８】
また、本発明は、対象組織を含む三次元空間に対して超音波の送受波を行って得られたボリュームデータを格納し、前記対象組織に対して設定された複数の代表断面に対する複数のマニュアルトレースを受け付け、前記複数の代表断面に加えて１又は複数の追加代表断面が設定された場合に、１又は複数の追加代表断面に対する１又は複数のマニュアルトレースを受け付け、前記各代表断面上に形成された各マニュアルトレースライン、及び、前記対象組織に対して設定された各非代表断面上に形成された各補間トレースライン、における各ポイントについて実際の組織輪郭に基づいて位置修正の可否を判断し、修正可能なポイントの位置を修正する修正処理を実行し、前記修正処理後の複数のマニュアルトレースライン及び複数の補間トレースラインからなるトレースライン列に基づいて対象組織を抽出し、あるいは、対象組織の体積を演算する、ことを特徴とするボリュームデータ処理方法に関する。
【発明の効果】
【００１９】
以上説明したように、本発明によれば、ユーザーにあまり負担をかけずに対象組織の特定精度あるいは計測精度を高められる。
【発明を実施するための最良の形態】
【００２０】
以下、本発明の好適な実施形態を図面に基づいて説明する。
【００２１】
図１には、本発明に係る超音波診断装置の全体構成がブロック図として示されている。この超音波診断装置は医療の分野において用いられ、特に生体内における対象組織を抽出し、その体積を演算する機能を有している。対象組織としては、胎盤、悪性腫瘍、胆嚢、甲状腺等を挙げることができる。以下の説明においては、胎盤を対象組織としてそれを抽出する処理が示されている。
【００２２】
図１において、３Ｄプローブ１０は、例えば体表面上に当接して用いられあるいは体腔内に挿入して用いられる超音波送受波器である。本実施形態において、３Ｄプローブ１０は２Ｄアレイ振動子を有している。２Ｄアレイ振動子は第一方向及び第二方向に整列した複数の振動素子によって構成されるものである。２Ｄアレイ振動子により超音波ビームが形成され、その超音波ビームは二次元走査される。これにより三次元空間としての三次元エコーデータ取込空間が形成される。具体的には、その三次元空間は複数の走査面の集合体として構成され、各走査面は超音波ビームを一次元走査することによって構成される。２Ｄアレイ振動子に代えて、１Ｄアレイ振動子を機械的に走査することによって、上記同様の三次元空間を形成することが可能である。
【００２３】
送信部１２は送信ビームフォーマーとして機能する。送信部１２は所定の遅延関係をもって複数の送信信号を上記２Ｄアレイ振動子に対して供給する。これによって送信ビームが形成される。生体内からの反射波は２Ｄアレイ振動子によって受波され、これによってその２Ｄアレイ振動子から複数の受信信号が受信部１４に対して出力される。受信部１４は複数の受信信号に対して整相加算処理を実行し、これによって整相加算後の受信信号（ビームデ−タ）を出力する。その受信信号に対しては所定の信号処理が施され、例えば検波、対数変換処理などが施され、信号処理後の受信信号であるビームデータが３Ｄデータメモリ１６に格納される。
【００２４】
３Ｄデータメモリ１６は、生体内の送受波空間である三次元空間に相当する三次元記憶空間を有している。３Ｄデータメモリ１６の書き込みあるいは読み出しに際して各データに対する座標変換が実行される。本実施形態においては、３Ｄデータメモリ１６の書き込み時に、送受波座標系から記憶空間座標系への座標変換が実行される。これによって、後に説明するようにボリュームデータが構成される。ボリュームデータは、複数の走査面に相当する複数のフレームデータ（スライスデータ）の集合体であり、各フレームデータは複数のビームデータからなるものである。各ビームデータは深さ方向に整列した複数のエコーデータからなるものである。ちなみに、この３Ｄデータメモリ１６以降の各構成については専用ハードウェアとして構成することも可能であるし、ソフトウェアの機能として実現することも可能である。
【００２５】
三次元画像形成部１８は、３Ｄデータメモリ１６に格納されたボリュームデータに基づいて、例えばボリュームレンダリング法に基づく画像処理を実行し、これによって三次元超音波画像を構成する。その画像データは表示処理部２６に送られる。任意断層画像形成部２０は、ユーザーによって設定された三次元空間内の任意断面に対応する断層画像を形成するものである。その場合において３Ｄデータメモリ１６内における任意断面に相当するデータアレイが読み出され、それらに基づいて任意断層画像としてのＢモード画像が形成される。その画像データは表示処理部２６に送られる。
【００２６】
組織抽出部２２は、以下に詳述する本実施形態特有の画像処理を実行するモジュールであり、三次元空間内あるいはボリュームデータ内に存在する対象組織（あるいは対象組織データ）をトレース処理によって抽出するものである。この場合においては、マニュアルトレース処理と補間処理とが併用されており、また、それぞれの処理結果に対する自動修正処理が利用されている。それらについては後に詳述する。抽出された対象組織データは表示処理部２６に送られ、対象組織の画像表示に利用される他、本実施形態においては体積演算部２４に送られている。
【００２７】
この体積演算部２４は上述した体積演算法（Disk Summation法）により対象組織の体積を求めるモジュールである。すなわち、組織抽出部２２により、対象組織の全体にわたって複数の閉ループとしてのトレースライン列が構成されるため、それらのトレースラインに基づいて対象組織の体積を近似的に求めるものである。その場合においては、各閉ループ間すなわち各断面間の距離も利用される。演算された体積値のデータは表示処理部２６へ送られる。体積演算法としては上記の他にAverage Rotation法等がある。
【００２８】
上記の三次元画像形成部１８、任意断層画像形成部２０及び組織抽出部２２等の各モジュールは、ユーザーによって選択された動作モードに応じて機能し、表示処理部２６にはそれぞれのモードに対応するデータが入力される。表示処理部２６は、入力されるデータに対して画像合成処理や色付け処理などの処理を施し、その結果であるデータを表示部２８に出力する。表示部２８には、選択された動作モードに応じて、三次元超音波画像、任意断層画像、抽出された組織の三次元画像、体積値等が表示される。ちなみに、三次元空間全体の三次元画像と対象組織の三次元画像とを合成して表示することも可能である。
【００２９】
制御部３０は図１に示される各構成の動作制御を行っており、特に入力部３２によってユーザーにより設定されたパラメータに基づいて、上述した組織抽出処理及び体積演算の動作を制御している。また、制御部３０は３Ｄデータメモリ１６へのデータ書き込み制御等を担っている。入力部３２は、キーボードやトラックボールなどを有する操作パネルによって構成されるものである。制御部３０はＣＰＵ及び動作プログラムなどによって構成されるものである。単一のＣＰＵが三次元画像処理、任意断層画像形成処理、組織抽出処理及び体積演算を実行するようにしてもよい。
【００３０】
次に、図２以降の各図を参照しつつ対象組織の抽出処理を具体的に説明する。図２にはそのプロセスがフローチャートとして示されている。
【００３１】
Ｓ１０１では、上述した３Ｄプローブを用いてデータが収集される。この場合においては、複数の走査面が形成される。図３には、一つの走査面に対応する断面データ４０が示されている。本実施形態では、観測対象は胎盤であり、図３にはその断面が符号４２Ａで示されている。ちなみに、符号４１は胎児を示しており、符号４３は子宮を示しており、符号４５は羊水を示している。図示されるように胎盤４２Ａは子宮４３の内壁に接合されており、断面データ４０上においてその接合部分、すなわち境界部分については輝度差があまりなく不明瞭なものとなっている。
【００３２】
図２のＳ１０２では、図１に示した３Ｄデータメモリ上にボリュームデータが構築される。すなわち、図３に示したような断面データ４０が３Ｄデータメモリ１６に順次格納され、当該メモリ上にボリュームデータが構築されることになる。その様子が図４に示されている。ボリュームデータ（三次元空間）４４内には胎盤（胎盤データ）４２が含まれている。なお、他の組織については図示省略されている。
【００３３】
図２のＳ１０３においては、表示部に任意断層画像を表示させながら、その断面の位置を適宜調整することにより、図５に示されるような基準断面４６が設定される。この場合においては、対象組織である胎盤４２の全体が断面として表れるように、特に、この断面が最大となるように基準断面４６の位置決めを行うのが望ましい。ただし、後に説明するように断面セットとしての参照断面列（断面セット）が設定されるため、基準断面４６はそのような参照断面列が対象組織である胎盤４２の全体をカバーするように設定されれば充分である。図６には基準断面４６に相当する断層画像４６Ａが拡大図として示されている。胎盤（の断面）４２Ｂは子宮（の断面）４３に貼り付いている。
【００３４】
図２のＳ１０４では、図７に示されるように基準断面に相当する断層画像４６Ａ上において、対象組織である胎盤４２Ｂについて両端の点がユーザーにより設定される。両端の点は図７において符号５０，５２で表されている。それらを結ぶ直線が基準線５４である。またこのＳ１０４では以下に説明する参照断面列を構成する断面の個数ｍが設定される。
【００３５】
基準線５４を設定すると、図８に示されるように、三次元空間に相当するボリュームデータ４４に対して参照断面列５６が自動的に生成される。参照断面列５６は、図７に示した基準線に直交する複数の断面として構成され、すなわち、それは一方の点５０から他方の点５２までにわたって等間隔あるいは非等間隔で複数の断面を並べたものに相当する。ここで、参照断面列５６は、複数のマニュアルトレース用参照断面５８と、複数の自動トレース用参照断面６０とで構成される。ここで、マニュアルトレース用参照断面５８は予め定められた個数だけ形成されており、その個数はｎである。例えば５〜１０の範囲内にｎの値が定められる。後述する事後的な追加処理が適用されると、ｎ個を越えて、マニュアルトレース用参照断面が追加設定されることもある。マニュアルトレース用参照断面は代表断面に相当するものであり、その代表断面についてだけマニュアルトレースが行われるので、ユーザーの負担を大幅に軽減できる。一方、後に説明するように個々の自動トレース用参照断面６０については補間処理により自動トレースが実行される。
【００３６】
具体的に説明すると、図２のＳ１０５においては選択されたｎ個の参照断面が画面上に断層画像として表示される。図９には、ｎ個の参照断面５８が示されている。図１０には複数の参照断面に対応する複数の断層画像５８Ａが示されている。それらの断層画像５８Ａは画像リスト（画像集）６４を構成するものである。この場合において、一枚ずつ断層画像を表示させてもよいし、複数の断層画像を並べて同時に表示するようにしてもよい。
【００３７】
図２のＳ１０６では、個々の断層画像５８Ａ（図１０）に対してマニュアルトレース処理が実行され、すなわち入力部が用いられて画像観察を行いながらユーザーにより胎盤の輪郭すなわち境界がトレースされる。そのトレース結果が図１１に示されている。符号４２Ｂは胎盤の断面を示し、符号４２Ｃは胎盤の輪郭を表している。その輪郭は子宮に接触して不明瞭となっている部分と羊水４５に接触して比較的明瞭になっている部分とを有している。不明瞭な部分であっても人間の視覚判断によりある程度のトレースを行うことが可能である。符号６６はマニュアルトレースラインが表されており、閉ループとして輪郭４２Ｃがなぞられている。マニュアルトレースライン６６は子宮４３側の部分６６Ｂと羊水４５側の部分６６Ａを有しているが、いずれも細かい凹凸までは対応できておらず概形がラインとして模擬されている。
【００３８】
図２のＳ１０７では、図１２に示されるようにマニュアルトレースラインの自動的な修正処理が実行される。この修正処理は各マニュアルトレースラインごとに行われ、すなわち代表断面としてのマニュアルトレース用参照断面ごとに実行される。図１２に示されるように、マニュアルトレースライン上の個々のポイントごとにその周辺に対してエッジ検出処理が実行され、エッジが検出されたポイントについてはエッジ上の位置に当該ポイントをシフトさせる処理が実行される。ただし、エッジが検出されないような場合にはマニュアルトレース結果がそのまま保存される。つまり、各マニュアルトレースラインに第１修正処理が実行される。その結果、図１２に示されるように羊水４５側のライン部分については輪郭４２Ｃに沿って忠実にトレースした自動トレースライン部分６８Ａが構成され、その一方、子宮４３側の部分については輝度差があまりないため自動トレースを適正に行うことができない結果、マニュアルトレースライン部分６８Ｂとして元の形状が保存される。その部分については細かな凹凸形状をトレースして結果を得ることはできないが、自動トレース処理による誤差等を考慮すればマニュアルトレース結果を尊重することが優先される。従って、このような処理によれば、マニュアルトレースと、自動トレースあるいは自動的トレースライン修正と、のそれぞれの利点及び欠点に着目し、両者の利点を巧みに引き出し、また両者の欠点を互いに補い合うことが可能となる。また自動的なトレースライン修正を前提とすれば、ユーザーによってマニュアルトレースを行う場合における負担を軽減できるという利点もある。さらに抽出精度及び体積演算精度を高められるという利点がある。
【００３９】
図１３には、ポイントの修正の具体的な手法が示されている。トレースライン６６上のあるポイント７２に着目すると、当該ポイント７２からトレースライン６６に直交する方向７４が定義され、当該方向７４に沿ってポイント７２からエッジ探索が両方向に実行される。そのエッジ探索の範囲が符号７６で示されている。エッジ探索は個々の位置において微分値を演算し、その微分値が所定値以上となった場合にエッジとみなすものであり、図示の例においては、輪郭４２Ｃへ探索が到達した時点で微分値が増大した結果、その位置７８が新たなトレースラインのポイントとして更新されることになる。このような処理を各トレースライン上の個々のポイントについて適用すれば、結果として輪郭が鮮明な部分についてはその輪郭に合致したトレースライン部分を描くことが可能となる。
【００４０】
本実施形態においては微分処理が適用されていたが、エッジを検出できる限りにおいて各種の処理を適用することが可能である。また本実施形態においては輝度差があまりないような不鮮明な部分については強制的なポイント修正が行われないため、すなわちマニュアルでのトレース結果が尊重されるため、トレース修正が過剰に行われてしまう問題を未然に防止できる。ちなみに探索する方向については直交方向以外にも各種の方向が考えられ、例えば組織の中心点が既知であるような場合にはその方向とポイントとを結ぶ直線として探索方向を定義するようにしてもよい。
【００４１】
図２のＳ１０８では、図８に示した複数の自動トレース用参照断面について自動トレース処理が実行される。この場合においては、図１４に示されるように、複数のマニュアルトレース用参照断面５８上に形成された、第１修正処理後の複数のマニュアルトレースラインとしての複数の複合トレースラインが基準とされ、それらを基礎として補間処理を行うことにより、複数の閉ループを面状に繋ぎ合わせたものとしてトレース面７０が構成される。この場合、完全なる三次元的な曲面を定義する必要はないが、少なくとも個々の自動トレース用参照断面６０について補間トレースラインが定義できるように補間処理を実行する。
【００４２】
図１５には非代表断面に相当する自動トレース用参照断面の断面データ６０が示されている。対象組織である胎盤４２Ｂの輪郭４２Ｃをおおよそ取り囲む閉ループとして補間トレースライン７４が自動的に設定されている。ただし、それは、実際のトレース結果ではなく、それを挟んで存在する二つあるいは複数の複合トレースラインに基づく補間処理によって生成されたものであるので、補間トレースライン７４は輪郭４２Ｃから全体的にあるいは部分的に外れている。
【００４３】
そこで、図２のＳ１０９においては、図１２及び図１３を用いて説明したように、自動的にトレースラインを修正する処理が実行される。具体的には、図１６に示されるように、羊水４５に接しているライン部分について自動的なトレースの修正が行われ、その一方、胎盤４３に接している部分については補間によるトレース結果がそのまま保存、尊重される。これにより部分的に修正された補間トレースライン７６が構成される。つまり、第２修正処理が実行される。このような処理は個々の自動トレース用参照断面について実行される。
【００４４】
以上の処理を行うと、図１７に示されるように、対象組織としての胎盤全体をその形態に沿って包み込んだ自動修正処理後のトレース面を定義することができ、当該トレース面７２は具体的には複数の閉ループすなわち複数のトレースラインからなるものである。それらのトレースラインは、上述した複数の複合トレースライン（第１修正処理後のマニュアルトレースライン）及び複数の第２修正処理後の補間トレースラインからなるものである。そこで、Ｓ１１１においては、抽出された組織の三次元画像が表示され、また演算された体積値が画面上に表示される。
【００４５】
従来においては、高精度の体積演算のためには、例えば２００枚あるいは３００枚といった多数の断面に対してマニュアルトレースを正確に行う必要があり、極めてユーザー負担が生じていたが、あるいは事実上そのような処理を行えないために体積演算値の信頼性が低下していたが、上記実施形態によれば、複数の代表断面についてだけマニュアルトレースを行えばよいので、ユーザーの負担を大幅に軽減しつつ、しかも高精度の演算結果を得ることが可能である。すなわち、マニュアルトレース結果に対しては自動的なトレース修正処理が適用されるため、ユーザーとしては組織の輪郭が明瞭な部分についてはあまり細かな凹凸にとらわれずにトレース処理を行え、その面においてもユーザー負担を軽減できる。また、非代表断面についてはマニュアルトレースを行う必要がなく、自動トレース処理を適用し、さらにその自動トレース処理結果に対して自動的な修正処理を適用して信頼性の高いものに加工できるので、その意味においてもユーザー負担を軽減でき、また高精度の演算を行えるという利点を得られる。
【００４６】
上記実施形態においては、基準断面の設定、両端点の指定、代表断面数の指定、また断面セット全体の個数等についてユーザーによって設定が行われていたが、それらの全部あるいは一部を自動的に設定するようにしてもよい。また状況に応じて自動的にそれらの値を可変するようにしてもよい。上記実施形態において、複数の断面は基本的には互いに均等間隔に設定されるが、例えば組織の形態変化が激しい部位などについては断面の間隔を小さくし、すなわち密に断面を設定し、それ以外の部位においては断面の間隔を大きくし、すなわち粗に断面を設定するようにしてもよい。また、上記実施形態においては、自動的なトレース修正の際に修正の可否がポイント毎に自動的に判断されているため、過剰に修正が行われてしまうような問題を未然に防止できる。すなわち境界あるいは輪郭の明瞭性に応じてマニュアルトレース結果あるいは補間トレース結果を尊重するかあるいは修正を施すのかを適応的に判断することが可能であり、様々な状況に対応した処理を実現できるという利点がある。ちなみに、自動トレース及びトレースラインの修正方法としては従来から提案されている各種の手法を適用することが可能である。
【００４７】
次に、図１８〜図２３を用いて変形例について説明する。図１８には、対象組織の抽出処理プロセスがフローチャートとして示されている。なお、図２に示した構成と同様の工程には同一符号を付し、その説明を省略する。
【００４８】
図１８に示されるプロセスにおいては、図２におけるＳ１０７とＳ１０８との間に、必要に応じて、１又は複数の追加のマニュアルトレース用参照断面（つまり、１又は複数の追加のマニュアルトレースライン）を生成する工程Ｓ２００を含んでいる。この工程Ｓ２００は、形状の変化が大きい部位について、マニュアルトレース遂行密度を高めて組織抽出精度を向上するためのものである。この処理内容が図１９に概念的に示されている。以下に説明する隣接マニュアルトレースラインペア１００−１，１００−２，１００−３，…の中で、対象組織部位の形状変化が大きく、このために類似度が小さくなっている隣接マニュアルトレースラインペアが特定された場合に、当該隣接マニュアルトレースラインペア（図１９において、１００−１，１００−２参照）の間に、追加のマニュアルトレース用参照断面１０２が事後的に追加設定される。なお、図１８に示したプロセスにおいて、追加処理工程Ｓ２００に相当する工程をＳ１０６とＳ１０７との間に設けることも可能である。その場合には、１又は複数の追加のマニュアルトレースラインが生成された後に、それらに対して自動修正処理（Ｓ１０７）が適用される。
【００４９】
図１８に示されるプロセスを具体的に説明する。Ｓ２０１以降の一連の工程が隣接マニュアルトレースラインペアを１つずつシフトさせながら、隣接マニュアルトレースラインペアごとに実行される。隣接マニュアルトレースラインペアは、隣り合う２つのマニュアルトレースライン（自動修正を経ている場合には２つの複合トレースライン）で構成される。追加処理前のマニュアルトレース用参照断面の個数をｎとした場合、ｎ−１個の隣接マニュアルトレースラインペアが構成される。
【００５０】
Ｓ２０１では、現在注目している隣接マニュアルトレースラインペア間で類似度が演算される。類似度の演算に当たっては例えば相互相関演算が実行される。すなわち、２つのマニュアルトレースラインの形状の類似度を表す情報が求められる。当該情報は、マニュアルトレース用参照断面を事後的に追加するか否か、換言すれば、形状の変化が大きい部位であるか否かを判断する指標となるものである。
【００５１】
図２０には、類似度演算の具体例が示されている。まず、隣接関係にある２つのマニュアルトレースライン１０２Ａ，１０２Ｂのそれぞれがサンプリングされて、それぞれに対して、後述するペアリングを適正に行う観点から、ポイント列１０４Ａ，１０４Ｂが設定される。そして、所定のルールに従って、それらのポイント列１０４Ａ，１０４Ｂ間で、複数のポイントペアが定義される。各ポイントペアは、対応関係（符号１０６参照）にある２つのポイントで構成される。そのルールとしては、状況に応じて相応しい各種のルールを選択的に採用することが可能である。例えば、各ポイント列１０４Ａ，１０４Ｂについて、それぞれ、最も近接している関係にある１番目のポイントペアを特定し、そこから順次、所定方向に回りながらペアリングを行って、最終的にＮ番目のポイントペアを特定するようにしてもよい。
【００５２】
図２１には、そのような処理の結果が示されており、隣接する２つのマニュアルトレース用参照断面５８Ａ，５８Ｂ上の２つのマニュアルトレースライン１０２Ａ，１０２Ｂ上には２つのポイント列１０４Ａ，１０４Ｂが設定されている。その後、それぞれのマニュアルトレースライン１０２Ａ，１０２Ｂについて重心点が求められ、各ポイントごとに重心点と当該ポイントとの間の距離が演算される。符号１０８に示されるように、各ポイント列１０４Ａ，１０４Ｂごとに、サンプル点（ポイント）番号ｔごとに、重心までの距離を表すと、グラフ（関数）ｆA(t)，ｆB(t)が得られる。そして、以下の（１）式の相互相関演算を実行することにより、相互相関係数として類似度αが得られる（類似度α：０．０〜１．０）。なお、類似度αを求める方法は上記のものには限られない。いずれにしても、２つのトレースライン間における形状の類似度を指標する情報が求められる。
【００５３】
【数１】
2008142519.tif 000003
【００５４】
図１８に戻って、Ｓ２０２では、上記類似度αが所定値以下か否かが判断される。所定値よりも類似度αが大きければ類似性ありと判断され、処理がＳ２０５へ移行し、断面追加は行われない。一方、所定値よりも類似度αが小さければ類似性なし、つまり対象組織の形状変化が大きいと判断されて、Ｓ２０３において、隣接する２つのマニュアルトレース用参照断面の間に、１又は複数の追加のマニュアルトレース用参照断面が設定される。追加する個数については１つであってもよいし、状況に応じて可変設定されてもよい。その後、追加された個々のマニュアルトレース用参照断面上において、上記Ｓ１０５〜Ｓ１０７と同様の処理が適用される。すなわち、個々の追加参照断面が画像として表示され、その画像を見ながらマニュアルトレースが実行され、その後、それによる個々のマニュアルトレースラインが自動的に修正される。なお、このような追加のマニュアルトレース作業は、すべての隣接マニュアルトレースラインペアに対して追加要否の判断がなされた後に行われてもよい。Ｓ２０５においては、すべての隣接マニュアルトレースラインペアについて上記の追加要否の判断及び必要な追加処理が行われたか否かが判断されて、残りの隣接マニュアルトレースラインペアがあればそれに対して上記処理が繰り返し適用される。
【００５５】
上記のＳ２０３に関し、マニュアルトレース用参照断面の追加方法については幾つかの手法が考えられる。それについて図２２及び図２３を用いて説明する。なお、選択する手法に応じて、図１８に示したプロセスの修正が必要となる場合がある。
【００５６】
図２２には、３つのタイプ（追加方法）が示されている。上段のタイプが適用された場合には、上記の所定値（閾値）が１つ利用されて、その所定値よりも類似度αが小さくなった場合に、１つだけマニュアルトレース用参照断面が追加される。上段には、更に、その条件と具体例が示されている。具体例においては、参照断面ＢＣ間、参照断面ＣＤ間にそれぞれ１つ参照断面が追加されている（ハッチング参照）。中段のタイプでは、それぞれ異なる複数の閾値が利用される。例えば、条件として図示されているように複数の閾値が利用され、各閾値区間ごとに追加断面数として１，２又は３が選択される。勿論、類似度が最高閾値よりも大きければ参照断面の追加は行われない。その具体例においては、参照断面ＢＣ間に１つの参照断面が追加されており、参照断面ＣＤ間には２つの参照断面が追加されている。下段のタイプにおいては、追加数として予めｋが定められおり、ｋ個の追加断面が動的に配分される。これを図２３を用いて説明する。最初に、増分Δαが（αMAX−αMIN）／ｋで求まっているものとする。αMAXは初期状態におけるすべてのαの中での最大値を示し、αMINは初期状態におけるすべてのαの中の最小値を示す。図２３の（１）に示すように、最初に、すべての断面間において最も小さなαを有する断面間が特定され、そこに１つの追加断面が挿入される。ここではＣＤ間に１つの追加断面の挿入が決定される。この処理に伴い、当該最も小さなα（ＣＤ間の当初α）に上記Δαが加算される（その時点での最小値のかさ上げ）。次に、図２３の（２）に示すように、先の加算処理後の状態において、最も小さなαが特定される。ここでは、やはりＣＤ間のαがなお最も小さいために、ＣＤ間に更に１つの追加断面の挿入が決定される。そして、上記同様に、当該加算修正後のαに更にΔαが加算される。次に図２３の（３）に示すように、２番目の加算処理後の状態において、最も小さなαが特定される。ここでは、ＢＣ間のαが最も小さいため、ＢＣ間に追加断面の挿入が決定される。このようにして、追加断面の挿入がｋ回繰り返されると、追加断面の挿入処理が終了する。この処理条件及び具体例が図２２の下段にも示されている。上記においては、幾つかの追加処理方法を示したが、状況に応じて適当な手法を選択すればよい。図２２の下段に示した方法によれば、追加上限の範囲内でより追加すべき区間により多くの断面を追加できるという利点がある。図２２の上段及び中段に示した方法によれば簡便な処理を期待できる。
【図面の簡単な説明】
【００５７】
【図１】本発明に係る超音波診断装置の実施形態を示すブロック図である。
【図２】対象組織の抽出処理を説明するためのフローチャートである。
【図３】子宮についての断面データを示す図である。
【図４】胎盤を含むボリュームデータを示す図である。
【図５】ボリュームデータに対する基準断面の設定を説明するための図である。
【図６】基準断面に相当する断層画像を示す図である。
【図７】断層画像上における基準線の設定を説明するための図である。
【図８】ボリュームデータに対して基準線を基準に自動的に形成された参照断面列を示す図である。
【図９】参照断面列に含まれる複数のマニュアルトレース用参照断面を示す図である。
【図１０】複数のマニュアルトレース用参照断面に対応する複数の断層画像を示す図である。
【図１１】マニュアルトレースラインを説明するための図である。
【図１２】マニュアルトレースラインの自動的な修正を説明するための図である。
【図１３】トレースライン上の点の位置の自動的な修正を説明するための図である。
【図１４】補間処理による複数の補間トレースラインの生成方法を説明するための図である。
【図１５】補間処理によって生成された補間トレースラインを説明するための図である。
【図１６】自動的に修正された補間トレースラインを示す図である。
【図１７】対象組織の全体に沿って設定された複数のトレースラインからなるトレース面を示す図である。
【図１８】他の実施形態に係る抽出処理を示すフローチャートである。
【図１９】マニュアルトレース用参照断面の追加を概念的に示す図である。
【図２０】隣接マニュアルトレースラインペアの形状関係を示す図である。
【図２１】隣接マニュアルトレースラインペア間における類似度の演算方法を説明するための図である。
【図２２】幾つかの断面追加処理方法を説明するための図である。
【図２３】図２２の下段に示す方法を具体的に説明するための図である。
【符号の説明】
【００５８】
１０  ３Ｄプローブ、１２  送信部、１４  受信部、１６  ３Ｄデータメモリ、１８  三次元画像形成部、２０  任意断層画像形成部、２２  組織抽出部、２４  体積演算部、２６  表示処理部。
【図１】
2008142519.tif 000004
【図２】
2008142519.tif 000005
【図３】
2008142519.tif 000006
【図４】
2008142519.tif 000007
【図５】
2008142519.tif 000008
【図６】
2008142519.tif 000009
【図７】
2008142519.tif 000010
【図８】
2008142519.tif 000011
【図９】
2008142519.tif 000012
【図１０】
2008142519.tif 000013
【図１１】
2008142519.tif 000014
【図１２】
2008142519.tif 000015
【図１３】
2008142519.tif 000016
【図１４】
2008142519.tif 000017
【図１５】
2008142519.tif 000018
【図１６】
2008142519.tif 000019
【図１７】
2008142519.tif 000020
【図１８】
2008142519.tif 000021
【図１９】
2008142519.tif 000022
【図２０】
2008142519.tif 000023
【図２１】
2008142519.tif 000024
【図２２】
2008142519.tif 000025
【図２３】
2008142519.tif 000026