(57)【要約】  【課題】  従来より短時間で動脈壁の弾性率を求めることができる超音波診断装置を提供する。  【解決手段】  超音波診断装置２０は、生体の動脈で反射した超音波信号に基づいて求められる動脈の壁の歪量のうち動脈を加振する加振器４２による加振に応じた周波数の成分を抽出する帯域制限フィルタ２７と、動脈内の血圧値のうち加振に応じた周波数の成分を抽出する帯域制限フィルタ２９と、帯域制限フィルタ２７によって抽出される周波数を加振に応じた信号である加振信号及び動脈の壁の歪量に基づいて設定するとともに帯域制限フィルタ２９によって抽出される周波数を加振信号及び動脈内の血圧値に基づいて設定する周波数解析部３０と、帯域制限フィルタ２７によって抽出された動脈の壁の歪量の成分及び帯域制限フィルタ２９によって抽出された動脈内の血圧値の成分に基づいて動脈の壁の弾性率を演算する弾性率演算部３１とを備える。【選択図】  図１
2006304965.tif 000002
【特許請求の範囲】
【請求項１】
生体の動脈で反射した超音波信号に基づいて求められる前記動脈の壁の歪量のうち前記動脈を加振する加振器による前記加振に応じた周波数の成分を抽出する歪量加振成分抽出手段と、前記動脈内の血圧値のうち前記加振に応じた周波数の成分を抽出する血圧値加振成分抽出手段と、前記歪量加振成分抽出手段によって抽出される前記周波数を前記加振に応じた信号である加振信号及び前記歪量に基づいて設定する歪量抽出周波数設定手段と、前記血圧値加振成分抽出手段によって抽出される前記周波数を前記加振信号及び前記血圧値に基づいて設定する血圧値抽出周波数設定手段と、前記歪量加振成分抽出手段によって抽出された前記成分及び前記血圧値加振成分抽出手段によって抽出された前記成分に基づいて前記壁の弾性率を演算する弾性率演算手段とを備えたことを特徴とする超音波診断装置。
【請求項２】
前記歪量及び前記血圧値の少なくとも一方に対して前記加振の再現性を周波数毎に取得する再現性取得手段を備えたことを特徴とする請求項１に記載の超音波診断装置。
【請求項３】
前記歪量及び前記血圧値の少なくとも一方と、前記加振との相関性を周波数毎に取得する相関性取得手段を備えたことを特徴とする請求項１に記載の超音波診断装置。
【請求項４】
前記歪量及び前記血圧値の少なくとも一方と、前記加振との間の伝達関数を取得する伝達関数取得手段を備えたことを特徴とする請求項１に記載の超音波診断装置。
【請求項５】
前記弾性率演算手段は、前記歪量加振成分抽出手段によって抽出された前記成分の複数回分の加算平均に基づいて前記弾性率を演算することを特徴とする請求項１に記載の超音波診断装置。
【請求項６】
生体の動脈で反射した超音波信号に基づいて求められる前記動脈の径の変化量のうち前記動脈を加振する加振器による前記加振に応じた周波数の成分を抽出する変化量加振成分抽出手段と、前記動脈内の血圧値のうち前記加振に応じた周波数の成分を抽出する血圧値加振成分抽出手段と、前記変化量加振成分抽出手段によって抽出される前記周波数を前記加振に応じた信号である加振信号及び前記変化量に基づいて設定する変化量抽出周波数設定手段と、前記血圧値加振成分抽出手段によって抽出される前記周波数を前記加振信号及び前記血圧値に基づいて設定する血圧値抽出周波数設定手段と、前記変化量加振成分抽出手段によって抽出された前記成分及び前記血圧値加振成分抽出手段によって抽出された前記成分に基づいて前記動脈の壁の弾性率を演算する弾性率演算手段とを備えたことを特徴とする超音波診断装置。
【請求項７】
前記変化量及び前記血圧値の少なくとも一方に対して前記加振の再現性を周波数毎に取得する再現性取得手段を備えたことを特徴とする請求項６に記載の超音波診断装置。
【請求項８】
前記変化量及び前記血圧値の少なくとも一方と、前記加振との相関性を周波数毎に取得する相関性取得手段を備えたことを特徴とする請求項６に記載の超音波診断装置。
【請求項９】
前記変化量及び前記血圧値の少なくとも一方と、前記加振との間の伝達関数を取得する伝達関数取得手段を備えたことを特徴とする請求項６に記載の超音波診断装置。
【請求項１０】
前記弾性率演算手段は、前記変化量加振成分抽出手段によって抽出された前記成分の複数回分の加算平均に基づいて前記弾性率を演算することを特徴とする請求項６に記載の超音波診断装置。
【請求項１１】
前記弾性率演算手段は、複数回の演算によって求めた複数の前記弾性率を加算平均することを特徴とする請求項１又は請求項２に記載の超音波診断装置。
【請求項１２】
前記弾性率演算手段は、前記生体の１心拍内の各心周期における前記弾性率を演算することを特徴とする請求項１又は請求項２に記載の超音波診断装置。
【請求項１３】
生体の動脈を加振する加振器と、超音波信号に基づいて前記動脈の壁の弾性率を求める超音波診断装置とを備え、
前記超音波診断装置は、前記動脈で反射した前記超音波信号に基づいて求められる前記壁の歪量のうち前記加振に応じた周波数の成分を抽出する歪量加振成分抽出手段と、前記動脈内の血圧値のうち前記加振に応じた周波数の成分を抽出する血圧値加振成分抽出手段と、前記歪量加振成分抽出手段によって抽出される前記周波数を前記加振に応じた信号である加振信号及び前記歪量に基づいて設定する歪量抽出周波数設定手段と、前記血圧値加振成分抽出手段によって抽出される前記周波数を前記加振信号及び前記血圧値に基づいて設定する血圧値抽出周波数設定手段と、前記歪量加振成分抽出手段によって抽出された前記成分及び前記血圧値加振成分抽出手段によって抽出された前記成分に基づいて前記壁の弾性率を演算する弾性率演算手段とを備えたことを特徴とする超音波診断システム。
【請求項１４】
生体の動脈を加振する加振器と、超音波信号に基づいて前記動脈の壁の弾性率を求める超音波診断装置とを備え、
前記超音波診断装置は、前記動脈で反射した前記超音波信号に基づいて求められる前記動脈の径の変化量のうち前記加振に応じた周波数の成分を抽出する変化量加振成分抽出手段と、前記動脈内の血圧値のうち前記加振に応じた周波数の成分を抽出する血圧値加振成分抽出手段と、前記変化量加振成分抽出手段によって抽出される前記周波数を前記加振に応じた信号である加振信号及び前記変化量に基づいて設定する変化量抽出周波数設定手段と、前記血圧値加振成分抽出手段によって抽出される前記周波数を前記加振信号及び前記血圧値に基づいて設定する血圧値抽出周波数設定手段と、前記変化量加振成分抽出手段によって抽出された前記成分及び前記血圧値加振成分抽出手段によって抽出された前記成分に基づいて前記壁の弾性率を演算する弾性率演算手段とを備えたことを特徴とする超音波診断システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波信号に基づいて生体内の動脈の壁（以下「動脈壁」という。）の弾性率を求める超音波診断装置に関するものである。
【背景技術】
【０００２】
生体組織の性状を識別し同定する手段の一つとして、生体組織を構成する弾性繊維、膠原線維、脂肪や血栓等によって弾性率に相違があることを利用して、生体組織に応力が加わったときの生体組織の歪量に基づいて生体組織の弾性率を求める手法が知られている。
【０００３】
例えば、生体組織に応力が加わったときの生体組織の変位を生体組織内に送信した超音波の反射波に基づいて計測し、計測した生体組織の変位に基づいて得られる生体組織の歪量と、生体組織に加わった応力の大きさとに基づいて生体組織の弾性率を求める手法が知られている（例えば、特許文献１、２参照。）。
【０００４】
【特許文献１】特開平１１−１８８０３６号公報
【特許文献２】特開平１０−５２２６号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、従来の手法においては、動脈壁の弾性率を求めるために動脈内の血圧変化値（脈圧）を用いた場合、最低血圧と最高血圧とを計測するために１心拍分の時間（ヒトの場合、約１秒）が必要であるという問題があり、更に加算平均等でＳＮ比を向上させる場合には加算回数に応じて複数心拍分の時間が必要となるという問題がある。
【０００６】
本発明は、従来の問題を解決するためになされたもので、従来より短時間で動脈壁の弾性率を求めることができる超音波診断装置を提供することを目的とする。
【課題を解決するための手段】
【０００７】
本発明の超音波診断装置は、生体の動脈で反射した超音波信号に基づいて求められる前記動脈の壁の歪量のうち前記動脈を加振する加振器による前記加振に応じた周波数の成分を抽出する歪量加振成分抽出手段と、前記動脈内の血圧値のうち前記加振に応じた周波数の成分を抽出する血圧値加振成分抽出手段と、前記歪量加振成分抽出手段によって抽出される前記周波数を前記加振に応じた信号である加振信号及び前記歪量に基づいて設定する歪量抽出周波数設定手段と、前記血圧値加振成分抽出手段によって抽出される前記周波数を前記加振信号及び前記血圧値に基づいて設定する血圧値抽出周波数設定手段と、前記歪量加振成分抽出手段によって抽出された前記成分及び前記血圧値加振成分抽出手段によって抽出された前記成分に基づいて前記壁の弾性率を演算する弾性率演算手段とを備えた構成を有している。
【０００８】
この構成により、本発明の超音波診断装置は、動脈壁の歪量のうち加振器による動脈の加振に応じた周波数の成分と、動脈内の血圧値のうち加振器による動脈の加振に応じた周波数の成分とに基づいて動脈壁の弾性率を演算するので、従来より短時間で動脈壁の弾性率を求めることができる。
【０００９】
また、本発明の超音波診断装置は、前記歪量及び前記血圧値の少なくとも一方に対して前記加振の再現性を周波数毎に取得する再現性取得手段を備えた構成を有している。
【００１０】
この構成により、本発明の超音波診断装置は、動脈壁の歪量及び動脈内の血圧値の少なくとも一方に対して、加振器による動脈の加振が安定して再現されているか否かを周波数毎に評価することが可能になる。
【００１１】
また、本発明の超音波診断装置は、前記歪量及び前記血圧値の少なくとも一方と、前記加振との相関性を周波数毎に取得する相関性取得手段を備えた構成を有している。
【００１２】
この構成により、本発明の超音波診断装置は、動脈壁の歪量及び動脈内の血圧値の少なくとも一方に対して、加振器による動脈の加振が影響している度合いを周波数毎に評価することが可能になる。
【００１３】
また、本発明の超音波診断装置は、前記歪量及び前記血圧値の少なくとも一方と、前記加振との間の伝達関数を取得する伝達関数取得手段を備えた構成を有している。
【００１４】
この構成により、本発明の超音波診断装置は、動脈壁の歪量及び動脈内の血圧値の少なくとも一方に対して、加振器による動脈の加振が伝達する状態を伝達関数の振幅特性及び位相特性に基づいて評価することが可能になる。
【００１５】
また、本発明の超音波診断装置の前記弾性率演算手段は、前記歪量加振成分抽出手段によって抽出された前記成分の複数回分の加算平均に基づいて前記弾性率を演算する構成を有している。
【００１６】
この構成により、本発明の超音波診断装置は、動脈壁の歪量のＳＮ比が低い場合であっても、ノイズの影響を抑制した動脈壁の弾性率を求めることができる。
【００１７】
また、本発明の超音波診断装置は、生体の動脈で反射した超音波信号に基づいて求められる前記動脈の径の変化量のうち前記動脈を加振する加振器による前記加振に応じた周波数の成分を抽出する変化量加振成分抽出手段と、前記動脈内の血圧値のうち前記加振に応じた周波数の成分を抽出する血圧値加振成分抽出手段と、前記変化量加振成分抽出手段によって抽出される前記周波数を前記加振に応じた信号である加振信号及び前記変化量に基づいて設定する変化量抽出周波数設定手段と、前記血圧値加振成分抽出手段によって抽出される前記周波数を前記加振信号及び前記血圧値に基づいて設定する血圧値抽出周波数設定手段と、前記変化量加振成分抽出手段によって抽出された前記成分及び前記血圧値加振成分抽出手段によって抽出された前記成分に基づいて前記動脈の壁の弾性率を演算する弾性率演算手段とを備えた構成を有している。
【００１８】
この構成により、本発明の超音波診断装置は、動脈の径の変化量のうち加振器による動脈の加振に応じた周波数の成分と、動脈内の血圧値のうち加振器による動脈の加振に応じた周波数の成分とに基づいて動脈壁の弾性率を演算するので、従来より短時間で動脈壁の弾性率を求めることができる。
【００１９】
また、本発明の超音波診断装置は、前記変化量及び前記血圧値の少なくとも一方に対して前記加振の再現性を周波数毎に取得する再現性取得手段を備えた構成を有している。
【００２０】
この構成により、本発明の超音波診断装置は、動脈の径の変化量及び動脈内の血圧値の少なくとも一方に対して、加振器による動脈の加振が安定して再現されているか否かを周波数毎に評価することが可能になる。
【００２１】
また、本発明の超音波診断装置は、前記変化量及び前記血圧値の少なくとも一方と、前記加振との相関性を周波数毎に取得する相関性取得手段を備えた構成を有している。
【００２２】
この構成により、本発明の超音波診断装置は、動脈の径の変化量及び動脈内の血圧値の少なくとも一方に対して、加振器による動脈の加振が影響している度合いを周波数毎に評価することが可能になる。
【００２３】
また、本発明の超音波診断装置は、前記変化量及び前記血圧値の少なくとも一方と、前記加振との間の伝達関数を取得する伝達関数取得手段を備えた構成を有している。
【００２４】
この構成により、本発明の超音波診断装置は、動脈の径の変化量及び動脈内の血圧値の少なくとも一方に対して、加振器による動脈の加振が伝達する状態を伝達関数の振幅特性及び位相特性に基づいて評価することが可能になる。
【００２５】
また、本発明の超音波診断装置の前記弾性率演算手段は、前記変化量加振成分抽出手段によって抽出された前記成分の複数回分の加算平均に基づいて前記弾性率を演算する構成を有している。
【００２６】
この構成により、本発明の超音波診断装置は、動脈の径の変化量のＳＮ比が低い場合であっても、ノイズの影響を抑制した動脈壁の弾性率を求めることができる。
【００２７】
また、本発明の超音波診断装置の前記弾性率演算手段は、複数回の演算によって求めた複数の前記弾性率を加算平均する構成を有している。
【００２８】
この構成により、本発明の超音波診断装置は、動脈壁の歪量又は動脈の径の変化量のＳＮ比が低い場合であっても、動脈内の血圧値のＳＮ比が低い場合であっても、ノイズの影響を抑制した動脈壁の弾性率を求めることができる。
【００２９】
また、本発明の超音波診断装置の前記弾性率演算手段は、前記生体の１心拍内の各心周期における前記弾性率を演算する構成を有している。
【００３０】
また、本発明の超音波診断システムは、生体の動脈を加振する加振器と、超音波信号に基づいて前記動脈の壁の弾性率を求める超音波診断装置とを備え、前記超音波診断装置は、前記動脈で反射した前記超音波信号に基づいて求められる前記壁の歪量のうち前記加振に応じた周波数の成分を抽出する歪量加振成分抽出手段と、前記動脈内の血圧値のうち前記加振に応じた周波数の成分を抽出する血圧値加振成分抽出手段と、前記歪量加振成分抽出手段によって抽出される前記周波数を前記加振に応じた信号である加振信号及び前記歪量に基づいて設定する歪量抽出周波数設定手段と、前記血圧値加振成分抽出手段によって抽出される前記周波数を前記加振信号及び前記血圧値に基づいて設定する血圧値抽出周波数設定手段と、前記歪量加振成分抽出手段によって抽出された前記成分及び前記血圧値加振成分抽出手段によって抽出された前記成分に基づいて前記壁の弾性率を演算する弾性率演算手段とを備えた構成を有している。
【００３１】
この構成により、本発明の超音波診断システムは、動脈壁の歪量のうち加振器による動脈の加振に応じた周波数の成分と、動脈内の血圧値のうち加振器による動脈の加振に応じた周波数の成分とに基づいて動脈壁の弾性率を演算するので、従来より短時間で動脈壁の弾性率を求めることができる。
【００３２】
また、本発明の超音波診断システムは、生体の動脈を加振する加振器と、超音波信号に基づいて前記動脈の壁の弾性率を求める超音波診断装置とを備え、前記超音波診断装置は、前記動脈で反射した前記超音波信号に基づいて求められる前記動脈の径の変化量のうち前記加振に応じた周波数の成分を抽出する変化量加振成分抽出手段と、前記動脈内の血圧値のうち前記加振に応じた周波数の成分を抽出する血圧値加振成分抽出手段と、前記変化量加振成分抽出手段によって抽出される前記周波数を前記加振に応じた信号である加振信号及び前記変化量に基づいて設定する変化量抽出周波数設定手段と、前記血圧値加振成分抽出手段によって抽出される前記周波数を前記加振信号及び前記血圧値に基づいて設定する血圧値抽出周波数設定手段と、前記変化量加振成分抽出手段によって抽出された前記成分及び前記血圧値加振成分抽出手段によって抽出された前記成分に基づいて前記壁の弾性率を演算する弾性率演算手段とを備えた構成を有している。
【００３３】
この構成により、本発明の超音波診断システムは、動脈の径の変化量のうち加振器による動脈の加振に応じた周波数の成分と、動脈内の血圧値のうち加振器による動脈の加振に応じた周波数の成分とに基づいて動脈壁の弾性率を演算するので、従来より短時間で動脈壁の弾性率を求めることができる。
【発明の効果】
【００３４】
本発明は、従来より短時間で動脈壁の弾性率を求めることができる超音波診断装置を提供することができるものである。
【発明を実施するための最良の形態】
【００３５】
以下、本発明の実施の形態について、図面を用いて説明する。
【００３６】
（第１の実施の形態）
まず、第１の実施の形態に係る超音波診断システムの構成について説明する。
【００３７】
図１及び図２に示すように、本実施の形態に係る超音波診断システム１０は、超音波信号を用いて生体８０内の動脈９０の壁である動脈壁９１の弾性率（以下「動脈壁弾性率」という。）を求める超音波診断装置２０と、生体８０内の最大血圧値（心臓収縮期の血圧値）及び最小血圧値（心臓拡張期の血圧値）を計測可能な血圧計４１と、生体８０の心拍の周期より短い周期で動脈９０を加振するとともに加振に応じた信号である加振信号を出力する加振器４２とを備えている。
【００３８】
超音波診断装置２０は、超音波振動子２１ａを有し生体８０の体表８１に接触させられる超音波プローブ２１と、超音波振動子２１ａを介して超音波信号の送信及び受信を行う送受信部２２と、送受信部２２によって送信される超音波信号を発生させる送信信号発生部２３と、送受信部２２によって送信される超音波信号及び送受信部２２によって受信される超音波信号に基づいて動脈壁９１の歪量（以下「動脈壁歪量」という。）を演算する歪量演算部２４と、送受信部２２によって送信される超音波信号及び送受信部２２によって受信される超音波信号に基づいて動脈９０の径の変化量（以下「動脈径変化量」という。）を演算する径変化量演算部２５と、血圧計４１によって計測された血圧値及び径変化量演算部２５によって演算された動脈径変化量に基づいて血圧値を演算する血圧値演算部２６と、歪量演算部２４によって演算された動脈壁歪量に対して帯域制限を行う帯域制限フィルタ２７と、血圧値演算部２６によって演算された血圧値に対して帯域制限を行う帯域制限フィルタ２９と、歪量演算部２４によって演算された動脈壁歪量、径変化量演算部２５によって演算された動脈径変化量、血圧値演算部２６によって演算された血圧値及び加振器４２によって出力される加振信号の周波数解析を行って帯域制限フィルタ２７、２９の遮断周波数を設定する周波数解析部３０と、動脈壁弾性率を演算する弾性率演算手段としての弾性率演算部３１と、弾性率演算部３１によって演算された動脈壁弾性率を表示する表示部３２と、送受信部２２、送信信号発生部２３、歪量演算部２４、径変化量演算部２５、血圧値演算部２６、帯域制限フィルタ２７、帯域制限フィルタ２９、周波数解析部３０、弾性率演算部３１及び表示部３２の動作を制御する制御部３３とを備えている。
【００３９】
帯域制限フィルタ２７、２９は、例えばＦＩＲ（Ｆｉｎｉｔｅ  Ｉｍｐｕｌｓｅ  Ｒｅｓｐｏｎｓｅ）フィルタ、ＩＩＲ（Ｉｎｆｉｎｉｔｅ  Ｉｍｐｕｌｓｅ  Ｒｅｓｐｏｎｓｅ）フィルタによって、実現されるのが一般的である。
【００４０】
次に、超音波診断システム１０の動作について説明する。
【００４１】
超音波診断装置２０は、血圧計４１が生体８０内の血圧値を計測し、加振器４２が動脈９０を加振している状態で、以下のように動作する。
【００４２】
制御部３３は、送信信号発生部２３に超音波信号を発生させ、送信信号発生部２３によって発生させられた超音波信号を超音波プローブ２１の超音波振動子２１ａを介して送受信部２２に送信させる。
【００４３】
超音波振動子２１ａから送信された超音波信号は、生体８０内の組織で反射する。制御部３３は、生体８０内の組織で反射した超音波信号を超音波プローブ２１の超音波振動子２１ａを介して送受信部２２に受信させる。
【００４４】
次いで、制御部３３は、動脈壁歪量を演算するよう歪量演算部２４に指示する。歪量演算部２４は、制御部３３から指示を受けると、超音波振動子２１ａによって送信される超音波信号の音響線２１Ａ（図２参照。）上の動脈壁９１のうち超音波振動子２１ａに近い方である動脈前壁９２（図２参照。）に対して動脈９０の内側に設定された内側計測点９２ａ（図２参照。）の位置変位と、動脈前壁９２に対して動脈９０の外側に設定された外側計測点９２ｂ（図２参照。）の位置変位とを、送受信部２２によって送信された超音波信号と、送受信部２２によって受信された超音波信号とに基づいて検出する。そして、歪量演算部２４は、検出した内側計測点９２ａの位置変位と、検出した外側計測点９２ｂの位置変位との差を動脈壁歪量として演算する。ここで、歪量演算部２４によって演算された動脈壁歪量は、加振器４２によって動脈９０が図３（ａ）に示すような周期で加振されているので、生体８０の拍動に起因する成分（以下「拍動成分」という。）に、加振器４２による動脈９０の加振に応じた成分（以下「加振成分」という。）が重畳されて図３（ｂ）に示すようになっている。なお、加振成分は、図３（ｂ）においては部分的にしか示されていないが、加振器４２によって動脈９０が加振されている間は常に拍動成分に重畳されている。
【００４５】
また、制御部３３は、動脈径変化量を演算するよう径変化量演算部２５に指示する。径変化量演算部２５は、制御部３３から指示を受けると、動脈前壁９２の内側計測点９２ａの位置変位と、超音波振動子２１ａによって送信される超音波信号の音響線２１Ａ上の動脈壁９１のうち超音波振動子２１ａから遠い方である動脈後壁９３（図２参照。）に対して動脈９０の内側に設定された内側計測点９３ａ（図２参照。）の位置変位とを、送受信部２２によって送信された超音波信号と、送受信部２２によって受信された超音波信号とに基づいて検出する。そして、径変化量演算部２５は、検出した動脈前壁９２の内側計測点９２ａの位置変位と、検出した動脈後壁９３の内側計測点９３ａの位置変位との差を動脈径変化量として演算する。ここで、径変化量演算部２５によって演算された動脈径変化量は、加振器４２によって動脈９０が図３（ａ）に示すような周期で加振されているので、拍動成分に加振成分が重畳されて図３（ｃ）に示すようになっている。なお、加振成分は、図３（ｃ）においては部分的にしか示されていないが、加振器４２によって動脈９０が加振されている間は常に拍動成分に重畳されている。
【００４６】
次いで、制御部３３は、血圧値を演算するよう血圧値演算部２６に指示する。血圧値演算部２６は、制御部３３から指示を受けると、径変化量演算部２５によって演算された動脈径変化量に基づいて現時点での血圧値を演算する。即ち、血圧値演算部２６は、径変化量演算部２５によって演算された動脈径変化量の立ち上がりを基準として心拍を判断し、動脈９０の径の最大値及び最小値の差に対する血圧計４１によって計測された最大血圧値及び最小血圧値の差（脈圧）の割合（以下「動脈径血圧値変化割合」という。）を何れかの１心拍間において予め求めておき、この動脈径血圧値変化割合と、血圧計４１によって予め計測されていた血圧値と、この血圧値が血圧計４１によって計測された時点から現時点までの動脈径変化量とに基づいて、現時点での血圧値を演算することができる。ここで、血圧値演算部２６によって演算された血圧値は、加振器４２によって動脈９０が図３（ａ）に示すような周期で加振されているので、歪量演算部２４によって演算された動脈壁歪量や径変化量演算部２５によって演算された動脈径変化量と同様に、拍動成分に加振成分が重畳されている。
【００４７】
次いで、制御部３３は、帯域制限フィルタ２７、２９の遮断周波数を設定するよう周波数解析部３０に指示する。
【００４８】
周波数解析部３０は、制御部３３から指示を受けると、加振器４２によって出力される加振信号の図４（ａ）に示すような周波数スペクトルのピークから加振信号の中心周波数を求め、求めた加振信号の中心周波数と、加振信号の周波数スペクトルの半値とに基づいて、中心周波数より低域側の遮断帯域幅と、中心周波数より高域側の遮断帯域幅とを求める。また、周波数解析部３０は、歪量演算部２４によって演算された動脈壁歪量の図４（ｂ）に実線で示すような周波数スペクトルのピークから動脈壁歪量の中心周波数を求め、求めた動脈壁歪量の中心周波数と、上述した低域側の遮断帯域幅とに基づいて、動脈壁歪量の低域側遮断周波数を求めるとともに、動脈壁歪量の中心周波数と、上述した高域側の遮断帯域幅とに基づいて、動脈壁歪量の高域側遮断周波数を求める。そして、周波数解析部３０は、求めた動脈壁歪量の低域側遮断周波数及び高域側遮断周波数を帯域制限フィルタ２７の遮断周波数として設定する。即ち、周波数解析部３０は、歪量抽出周波数設定手段を構成している。なお、図４（ｂ）に破線で示す周波数スペクトルは、加振器４２によって動脈９０が加振されていないときの動脈壁歪量の周波数スペクトル、即ち、動脈壁歪量のうち拍動成分の周波数スペクトルである。
【００４９】
また、周波数解析部３０は、帯域制限フィルタ２７の遮断周波数の設定と同様に、加振器４２によって出力される加振信号と、血圧値演算部２６によって演算された血圧値とに基づいて求めた血圧値の低域側遮断周波数及び高域側遮断周波数を帯域制限フィルタ２９の遮断周波数として設定する。即ち、周波数解析部３０は、血圧値抽出周波数設定手段を構成している。
【００５０】
次いで、制御部３３は、帯域制限を行うよう帯域制限フィルタ２７、２９に指示する。
【００５１】
帯域制限フィルタ２７は、制御部３３から指示を受けると、歪量演算部２４によって演算された動脈壁歪量に対して、周波数解析部３０によって設定された低域側遮断周波数以下の周波数成分と、周波数解析部３０によって設定された高域側遮断周波数以上の周波数成分とを遮断する。したがって、帯域制限フィルタ２７は、例えば図３（ｂ）に示すように拍動成分に加振成分が重畳された動脈壁歪量から、拍動成分をほぼ除去して、図３（ｄ）に示すようにほぼ加振成分のみを抽出することができる。即ち、帯域制限フィルタ２７は、歪量加振成分抽出手段を構成している。
【００５２】
また、帯域制限フィルタ２９は、制御部３３から指示を受けると、血圧値演算部２６によって演算された血圧値に対して、周波数解析部３０によって設定された低域側遮断周波数以下の周波数成分と、周波数解析部３０によって設定された高域側遮断周波数以上の周波数成分とを遮断する。したがって、帯域制限フィルタ２９は、拍動成分に加振成分が重畳された血圧値から、拍動成分を除去して、ほぼ加振成分のみを抽出することができる。即ち、帯域制限フィルタ２９は、血圧値加振成分抽出手段を構成している。
【００５３】
次いで、制御部３３は、動脈壁弾性率を演算するよう弾性率演算部３１に指示する。弾性率演算部３１は、制御部３３から指示を受けると、帯域制限フィルタ２７によって帯域制限を行った動脈壁歪量と、帯域制限フィルタ２９によって帯域制限を行った血圧値とに基づいて、動脈壁弾性率を演算する。
【００５４】
次いで、制御部３３は、動脈壁弾性率を表示するよう表示部３２に指示する。表示部３２は、制御部３３から指示を受けると、弾性率演算部３１によって演算された動脈壁弾性率を表示する。ここで、表示部３２は、一般的な超音波診断装置の基本機能であるＢモード断層画像上に動脈壁弾性率を重ねて表示しても良い。また、表示部３２は、周波数解析部３０で得た周波数スペクトルや、帯域制限フィルタ２７、２９の制御情報等を制御部３３を介して表示しても良い。
【００５５】
以上に説明したように、超音波診断装置２０は、動脈壁歪量のうち加振器４２による動脈９０の加振に応じた周波数の成分と、動脈９０内の血圧値のうち加振器４２による動脈９０の加振に応じた周波数の成分とに基づいて動脈壁９１の動脈壁弾性率を演算するので、従来より短時間で動脈壁９１の弾性率を求めることができる。
【００５６】
なお、歪量演算部２４は、検出した内側計測点９２ａの位置変位と、検出した外側計測点９２ｂの位置変位との差を動脈壁歪量として演算するようになっているが、他の方法によって動脈壁歪量を演算するようになっていても良い。例えば、歪量演算部２４は、動脈前壁９２の内側計測点９２ａと、外側計測点９２ｂとの間に複数の計測点を更に設定し、動脈前壁９２の歪量を層状に演算するようになっていても良い。また、歪量演算部２４は、内側計測点９２ａの運動速度と、外側計測点９２ｂの運動速度とを、送受信部２２によって送信された超音波信号と、送受信部２２によって受信された超音波信号とに基づいて検出し、検出した内側計測点９２ａの運動速度と、検出した外側計測点９２ｂの運動速度との差の時間積分を動脈壁歪量として演算するようになっていても良い。ここで、生体８０内の各計測点の運動速度の検出手法としては、ＦＦＴ（Ｆａｓｔ  Ｆｏｕｒｉｅｒ  Ｔｒａｎｓｆｏｒｍ）ドプラ法、自己相関法等、一般的に用いられている何れの手法が用いられても良い。
【００５７】
また、径変化量演算部２５は、検出した動脈前壁９２の内側計測点９２ａの位置変位と、検出した動脈後壁９３の内側計測点９３ａの位置変位との差を動脈径変化量として演算するようになっているが、他の方法によって動脈径変化量を演算するようになっていても良い。例えば、径変化量演算部２５は、動脈前壁９２の内側計測点９２ａの運動速度と、動脈後壁９３の内側計測点９３ａの運動速度とを、送受信部２２によって送信された超音波信号と、送受信部２２によって受信された超音波信号とに基づいて検出し、検出した動脈前壁９２の内側計測点９２ａの運動速度と、検出した動脈後壁９３の内側計測点９３ａの運動速度との差の時間積分を動脈径変化量として演算するようになっていても良い。
【００５８】
また、径変化量演算部２５は、動脈９０の内径の変化量を動脈径変化量として演算しているが、動脈９０の外径の変化量を動脈径変化量として演算するようになっていても良い。例えば、径変化量演算部２５は、動脈前壁９２の外側計測点９２ｂの位置変位と、超音波振動子２１ａによって送信される超音波信号の音響線２１Ａ上の動脈壁９１のうち超音波振動子２１ａから遠い方である動脈後壁９３に対して動脈９０の外側に設定された外側計測点９３ｂ（図２参照。）の位置変位とを、送受信部２２によって送信された超音波信号と、送受信部２２によって受信された超音波信号とに基づいて検出し、検出した動脈前壁９２の外側計測点９２ｂの位置変位と、検出した動脈後壁９３の外側計測点９３ｂの位置変位との差を動脈９０の外径の変化量として演算するようになっていても良い。
【００５９】
また、血圧値演算部２６は、径変化量演算部２５によって演算された動脈径変化量の立ち上がりを基準として心拍を判断するようになっているが、他の方法によって心拍を判断するようになっていても良い。例えば、血圧値演算部２６は、超音波診断装置に一般的に具備されている心電波形計測機能によって計測された心電波形のＰ波、Ｑ波、Ｒ波、Ｓ波、Ｔ波、Ｕ波（図５参照。）の何れかを基準として心拍を判断するようになっていても良い。また、血圧値演算部２６は、心音波形を基準として心拍を判断するようになっていても良い。
【００６０】
また、周波数解析部３０は、歪量演算部２４によって演算された動脈壁歪量、径変化量演算部２５によって演算された動脈径変化量、血圧値演算部２６によって演算された血圧値のそれぞれに対して加振器４２による加振の再現性を周波数毎に取得する再現性取得手段としても機能するようになっていても良い。超音波診断装置２０は、周波数解析部３０が再現性取得手段としても機能する場合、歪量演算部２４によって演算された動脈壁歪量、径変化量演算部２５によって演算された動脈径変化量、血圧値演算部２６によって演算された血圧値のそれぞれに対して、加振器４２による動脈９０の加振が安定して再現されているか否かを周波数毎に評価することが可能になる。なお、周波数毎の再現性を求める手法として、再現性評価関数を用いる手法が一般的に知られている。
【００６１】
また、周波数解析部３０は、歪量演算部２４によって演算された動脈壁歪量、径変化量演算部２５によって演算された動脈径変化量、血圧値演算部２６によって演算された血圧値のそれぞれと、加振器４２による加振との相関性を周波数毎に取得する相関性取得手段としても機能するようになっていても良い。超音波診断装置２０は、周波数解析部３０が相関性取得手段としても機能する場合、歪量演算部２４によって演算された動脈壁歪量、径変化量演算部２５によって演算された動脈径変化量、血圧値演算部２６によって演算された血圧値のそれぞれに対して、加振器４２による動脈９０の加振が影響している度合いを周波数毎に評価することが可能になる。なお、周波数毎の相関性を求める手法として、コヒーレンス関数を用いる手法が一般的に知られている。
【００６２】
また、周波数解析部３０は、歪量演算部２４によって演算された動脈壁歪量、径変化量演算部２５によって演算された動脈径変化量、血圧値演算部２６によって演算された血圧値のそれぞれと、加振器４２による加振との間の伝達関数を取得する伝達関数取得手段としても機能するようになっていても良い。超音波診断装置２０は、周波数解析部３０が伝達関数取得手段としても機能する場合、歪量演算部２４によって演算された動脈壁歪量、径変化量演算部２５によって演算された動脈径変化量、血圧値演算部２６によって演算された血圧値のそれぞれに対して、加振器４２による動脈９０の加振が伝達する状態を伝達関数の振幅特性及び位相特性に基づいて評価することが可能になる。
【００６３】
また、周波数解析部３０は、周波数スペクトルの半値幅によって遮断帯域幅を設定するようになっているが、他の手法によって遮断帯域幅を設定するようになっていても良い。例えば、周波数解析部３０は、遮断帯域幅を任意に設定するようになっていても良い。
【００６４】
また、弾性率演算部３１は、複数回の演算によって求めた複数の動脈壁弾性率を加算平均するようになっていても良い。弾性率演算部３１は、複数回の演算によって求めた複数の動脈壁弾性率を加算平均するようになっている場合、歪量演算部２４によって演算された動脈壁歪量と、血圧値演算部２６によって演算された血圧値との少なくとも一方のＳＮ比が低い場合であっても、ノイズの影響を抑制した動脈壁弾性率を求めることができる。例えば、弾性率演算部３１は、生体８０の心拍毎に演算した複数の動脈壁弾性率を加算平均することによって、生体８０の各心拍間で再現性のないノイズ成分を除去した動脈壁弾性率を求めることができる。また、弾性率演算部３１は、加振器４２による加振毎に演算した複数の動脈壁弾性率を加算平均することによって、加振器４２による各加振間で再現性のないノイズ成分を除去した動脈壁弾性率を求めることができる。
【００６５】
また、弾性率演算部３１は、歪量演算部２４によって演算された動脈壁歪量のうち帯域制限フィルタ２７によって抽出された成分の複数回分の加算平均に基づいて動脈壁弾性率を演算するようになっていても良い。弾性率演算部３１は、歪量演算部２４によって演算された動脈壁歪量のうち帯域制限フィルタ２７によって抽出された成分の複数回分の加算平均に基づいて動脈壁弾性率を演算するようになっている場合、歪量演算部２４によって演算された動脈壁歪量のＳＮ比が低い場合であっても、ノイズの影響を抑制した動脈壁弾性率を求めることができる。
【００６６】
また、弾性率演算部３１は、生体８０の１心拍内の各心周期（即ち、充満期、緊張期、駆出期及び弛緩期）における動脈壁弾性率を演算するようになっていても良い。
【００６７】
（第２の実施の形態）
まず、第２の実施の形態に係る超音波診断システムの構成について説明する。
【００６８】
なお、本実施の形態に係る超音波診断システムの構成のうち第１の実施の形態に係る超音波診断システム１０（図１参照。）の構成と同様な構成については、超音波診断システム１０の構成と同一の符号を付して詳細な説明を省略する。
【００６９】
図６に示すように、本実施の形態に係る超音波診断システム５０の構成は、超音波信号を用いて動脈壁弾性率を求める超音波診断装置６０を超音波診断装置２０（図１参照。）に代えて超音波診断システム１０が備えた構成と同様である。
【００７０】
超音波診断装置６０の構成は、径変化量演算部２５によって演算された動脈径変化量に対して帯域制限を行う帯域制限フィルタ６１を歪量演算部２４（図１参照。）及び帯域制限フィルタ２７（図１参照。）に代えて超音波診断装置２０が備えた構成と同様である。
【００７１】
帯域制限フィルタ６１は、例えばＦＩＲフィルタ、ＩＩＲフィルタによって、実現されるのが一般的である。
【００７２】
次に、超音波診断システム５０の動作について説明する。
【００７３】
超音波診断装置６０は、血圧計４１が生体８０内の血圧値を計測し、加振器４２が動脈９０を加振している状態で、以下のように動作する。
【００７４】
制御部３３は、送信信号発生部２３に超音波信号を発生させ、送信信号発生部２３によって発生させられた超音波信号を超音波プローブ２１の超音波振動子２１ａを介して送受信部２２に送信させる。
【００７５】
超音波振動子２１ａから送信された超音波信号は、生体８０内の組織で反射する。制御部３３は、生体８０内の組織で反射した超音波信号を超音波プローブ２１の超音波振動子２１ａを介して送受信部２２に受信させる。
【００７６】
次いで、制御部３３は、動脈径変化量を演算するよう径変化量演算部２５に指示する。径変化量演算部２５は、制御部３３から指示を受けると、第１の実施の形態と同様に、動脈径変化量を演算する。
【００７７】
次いで、制御部３３は、血圧値を演算するよう血圧値演算部２６に指示する。血圧値演算部２６は、制御部３３から指示を受けると、第１の実施の形態と同様に、血圧値を演算する。
【００７８】
次いで、制御部３３は、帯域制限フィルタ２９、６１の遮断周波数を設定するよう周波数解析部３０に指示する。
【００７９】
周波数解析部３０は、制御部３３から指示を受けると、第１の実施の形態と同様に、血圧値の低域側遮断周波数及び高域側遮断周波数を帯域制限フィルタ２９の遮断周波数として設定する。
【００８０】
また、周波数解析部３０は、帯域制限フィルタ２９の遮断周波数の設定と同様に、加振器４２によって出力される加振信号と、径変化量演算部２５によって演算された動脈径変化量とに基づいて求めた動脈径変化量の低域側遮断周波数及び高域側遮断周波数を帯域制限フィルタ６１の遮断周波数として設定する。即ち、周波数解析部３０は、変化量抽出周波数設定手段を構成している。
【００８１】
次いで、制御部３３は、帯域制限を行うよう帯域制限フィルタ２９、６１に指示する。
【００８２】
帯域制限フィルタ２９は、制御部３３から指示を受けると、第１の実施の形態と同様に、血圧値演算部２６によって演算された血圧値に対して、周波数解析部３０によって設定された低域側遮断周波数以下の周波数成分と、周波数解析部３０によって設定された高域側遮断周波数以上の周波数成分とを遮断する。したがって、帯域制限フィルタ２９は、拍動成分に加振成分が重畳された血圧値から、拍動成分を除去して、ほぼ加振成分のみを抽出することができる。
【００８３】
また、帯域制限フィルタ６１は、制御部３３から指示を受けると、径変化量演算部２５によって演算された動脈径変化量に対して、周波数解析部３０によって設定された低域側遮断周波数以下の周波数成分と、周波数解析部３０によって設定された高域側遮断周波数以上の周波数成分とを遮断する。したがって、帯域制限フィルタ６１は、例えば図３（ｃ）に示すように拍動成分に加振成分が重畳された動脈径変化量から、拍動成分をほぼ除去して、ほぼ加振成分のみを抽出することができる。即ち、帯域制限フィルタ６１は、変化量加振成分抽出手段を構成している。
【００８４】
次いで、制御部３３は、動脈壁弾性率を演算するよう弾性率演算部３１に指示する。弾性率演算部３１は、制御部３３から指示を受けると、帯域制限フィルタ６１によって帯域制限を行った動脈径変化量と、帯域制限フィルタ２９によって帯域制限を行った血圧値とに基づいて、動脈壁弾性率を演算する。
【００８５】
次いで、制御部３３は、動脈壁弾性率を表示するよう表示部３２に指示する。表示部３２は、制御部３３から指示を受けると、弾性率演算部３１によって演算された動脈壁弾性率を表示する。
【００８６】
以上に説明したように、超音波診断装置６０は、動脈径変化量のうち加振器４２による動脈９０の加振に応じた周波数の成分と、動脈９０内の血圧値のうち加振器４２による動脈９０の加振に応じた周波数の成分とに基づいて動脈壁９１の動脈壁弾性率を演算するので、従来より短時間で動脈壁９１の弾性率を求めることができる。
【００８７】
なお、弾性率演算部３１は、径変化量演算部２５によって演算された動脈径変化量のうち帯域制限フィルタ６１によって抽出された成分の複数回分の加算平均に基づいて動脈壁弾性率を演算するようになっていても良い。弾性率演算部３１は、径変化量演算部２５によって演算された動脈径変化量のうち帯域制限フィルタ６１によって抽出された成分の複数回分の加算平均に基づいて動脈壁弾性率を演算するようになっている場合、径変化量演算部２５によって演算された動脈径変化量のＳＮ比が低い場合であっても、ノイズの影響を抑制した動脈壁弾性率を求めることができる。
【００８８】
（第３の実施の形態）
まず、第３の実施の形態に係る超音波診断システムの構成について説明する。
【００８９】
なお、本実施の形態に係る超音波診断システムの構成のうち第１の実施の形態に係る超音波診断システム１０（図１参照。）や、第２の実施の形態に係る超音波診断システム５０（図６参照。）の構成と同様な構成については、超音波診断システム１０や、超音波診断システム５０の構成と同一の符号を付して詳細な説明を省略する。
【００９０】
本実施の形態に係る超音波診断システムの構成は、超音波信号を用いて動脈壁弾性率を求める超音波診断装置を、超音波診断装置２０（図１参照。）に代えて超音波診断システム１０が備えた構成や、超音波診断装置６０（図６参照。）に代えて超音波診断システム５０が備えた構成と同様である。
【００９１】
本実施の形態に係る超音波診断装置の構成は、図７に示すように、複数の超音波振動子１２１ａを有し生体８０の体表８１に接触させられる超音波プローブ１２１を超音波プローブ２１（図２参照。）に代えて超音波診断装置２０や、超音波診断装置６０が備えた構成と同様である。
【００９２】
次に、本実施の形態に係る超音波診断システムの動作について説明する。
【００９３】
なお、本実施の形態に係る超音波診断システムは、以下の動作を除いて、第１の実施の形態に係る超音波診断システム１０（図１参照。）や、第２の実施の形態に係る超音波診断システム５０（図６参照。）と同様に動作する。
【００９４】
超音波プローブ１２１は、複数の超音波振動子１２１ａが動脈９０の延在方向に配列されるように、生体８０の体表８１に接触させられる。
【００９５】
送受信部２２は、送信信号発生部２３によって発生させられた超音波信号が動脈９０の延在方向に走査するように、超音波プローブ１２１の複数の超音波振動子１２１ａを介して超音波信号を送信し、生体８０内の組織で反射した超音波信号を超音波プローブ１２１の複数の超音波振動子１２１ａを介して受信する。
【００９６】
径変化量演算部２５は、第１の実施の形態や、第２の実施の形態と同様に、動脈９０のうち超音波振動子１２１ａによって送信される超音波信号の音響線１２１Ａ上の箇所の動脈径変化量を、送受信部２２によって送信された超音波信号と、送受信部２２によって受信された超音波信号とに基づいて演算する。即ち、径変化量演算部２５は、複数の超音波振動子１２１ａのそれぞれに対応する動脈９０の延在方向の複数箇所の動脈径変化量を演算する。
【００９７】
また、血圧値演算部２６は、第１の実施の形態や、第２の実施の形態と同様に、動脈９０のうち超音波振動子１２１ａによって送信される超音波信号の音響線１２１Ａ上の箇所の動脈径血圧値変化割合を求める。即ち、血圧値演算部２６は、複数の超音波振動子１２１ａのそれぞれに対応する動脈９０の延在方向の複数箇所の動脈径血圧値変化割合を求める。そして、血圧値演算部２６は、複数箇所の動脈径血圧値変化割合の平均値と、血圧計４１によって予め計測されていた血圧値と、この血圧値が血圧計４１によって計測された時点から現時点までの動脈径変化量とに基づいて、現時点での血圧値を演算する。
【００９８】
以上に説明したように、本実施の形態に係る超音波診断装置は、動脈９０の形状や組織性状等の要因によって動脈９０の延在方向の複数箇所の動脈径変化量にばらつきが生じる場合であっても、血圧値演算部２６によって血圧値を最適に演算することができる。
【００９９】
なお、血圧値演算部２６は、複数箇所の動脈径血圧値変化割合の平均値に基づいて血圧値を演算するのではなく、複数箇所の動脈径血圧値変化割合の最大値及び最小値の何れか一方と、血圧計４１によって予め計測されていた血圧値と、この血圧値が血圧計４１によって計測された時点から現時点までの動脈径変化量とに基づいて、現時点での血圧値を演算するようになっていても、同様に血圧値を最適に演算することができる。
【０１００】
また、血圧値演算部２６は、複数箇所の動脈径血圧値変化割合の平均値に基づいて血圧値を演算するのではなく、複数箇所の動脈径血圧値変化割合の最小二乗値と、血圧計４１によって予め計測されていた血圧値と、この血圧値が血圧計４１によって計測された時点から現時点までの動脈径変化量とに基づいて、現時点での血圧値を演算するようになっていても、同様に血圧値を最適に演算することができる。
【産業上の利用可能性】
【０１０１】
以上のように、本発明に係る超音波診断装置は、従来より短時間で動脈壁の弾性率を求めることができるという効果を有し、超音波信号に基づいて生体内の動脈壁の弾性率を求める超音波診断装置等として有用である。
【図面の簡単な説明】
【０１０２】
【図１】本発明の第１の実施の形態に係る超音波診断システムのブロック図
【図２】図１に示す超音波診断装置の超音波プローブと、生体の動脈との模式図
【図３】（ａ）図１に示す加振器によって出力される加振信号を示す図  （ｂ）図１に示す歪量演算部によって演算される動脈壁歪量を示す図  （ｃ）図１に示す径変化量演算部によって演算される動脈径変化量を示す図  （ｄ）図１に示す帯域制限フィルタによって抽出された動脈壁歪量の成分を示す図
【図４】（ａ）図１に示す加振器によって出力される加振信号の周波数スペクトルを示す図  （ｂ）図１に示す歪量演算部によって演算される動脈壁歪量の周波数スペクトルを示す図
【図５】図１に示す超音波診断装置の心電波形計測機能によって計測された心電波形を示す図
【図６】本発明の第２の実施の形態に係る超音波診断システムのブロック図
【図７】本発明の第３の実施の形態に係る超音波診断システムの超音波診断装置の超音波プローブと、生体の動脈との模式図
【符号の説明】
【０１０３】
１０  超音波診断システム
２０  超音波診断装置
２７  帯域制限フィルタ（歪量加振成分抽出手段）
２９  帯域制限フィルタ（血圧値加振成分抽出手段）
３０  周波数解析部（歪量抽出周波数設定手段、血圧値抽出周波数設定手段、再現性取得手段、相関性取得手段、伝達関数取得手段、変化量抽出周波数設定手段）
３１  弾性率演算部（弾性率演算手段）
４２  加振器
５０  超音波診断システム
６０  超音波診断装置
６１  帯域制限フィルタ（変化量加振成分抽出手段）
８０  生体
９０  動脈
９１  動脈壁（壁）
【図１】
2006304965.tif 000003
【図２】
2006304965.tif 000004
【図３】
2006304965.tif 000005
【図４】
2006304965.tif 000006
【図５】
2006304965.tif 000007
【図６】
2006304965.tif 000008
【図７】
2006304965.tif 000009