▶ メイヨ フオンデーシヨン フオー メデイカル エジユケーシヨン アンド リサーチの特許一覧
特表2024-504635超音波振動測定を用いた膀胱過活動の非侵襲的決定のためのシステム及び方法
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公表特許公報(A)
(11)【公表番号】
(43)【公表日】2024-02-01
(54)【発明の名称】超音波振動測定を用いた膀胱過活動の非侵襲的決定のためのシステム及び方法
(51)【国際特許分類】
A61B 8/08 20060101AFI20240125BHJP
A61B 5/20 20060101ALI20240125BHJP
【ＦＩ】
A61B8/08
A61B5/20
【審査請求】未請求
【予備審査請求】未請求
(21)【出願番号】P 2023542881
(86)(22)【出願日】2022-01-17
(85)【翻訳文提出日】2023-09-11
(86)【国際出願番号】 US2022012675
(87)【国際公開番号】W WO2022155573
(87)【国際公開日】2022-07-21
(31)【優先権主張番号】63/138,116
(32)【優先日】2021-01-15
(33)【優先権主張国・地域又は機関】US
(81)【指定国・地域】
(71)【出願人】
【識別番号】500480919
【氏名又は名称】メイヨ フオンデーシヨン フオー メデイカル エジユケーシヨン アンド リサーチ
(74)【代理人】
【識別番号】100134832
【弁理士】
【氏名又は名称】瀧野 文雄
(74)【代理人】
【識別番号】100165308
【弁理士】
【氏名又は名称】津田 俊明
(74)【代理人】
【識別番号】100115048
【弁理士】
【氏名又は名称】福田 康弘
(72)【発明者】
【氏名】ファテミ モスタファ
(72)【発明者】
【氏名】アリザッド アズラ
(72)【発明者】
【氏名】ローゼン デイヴィッド ピー．
【テーマコード（参考）】
4C038
4C601
【Ｆターム（参考）】
4C038DD06
4C601DD18
4C601DD19
4C601DD23
4C601DD30
4C601EE11
4C601EE20
4C601JB34
4C601JB50
(57)【要約】
排尿筋過活動を決定するために、超音波膀胱振動測定（ＵＢＶ）を用いて湾曲した組織壁の粘弾性を決定するためのシステム及び方法が提供される。ＵＢＶは、特定の場合には、湾曲した膀胱壁にラム波を励起する集束超音波放射力、及び、そのような湾曲した壁を通って伝播する組織変形を追跡するパルスエコー技術を利用する、非侵襲的な技術である。
【選択図】図１
【選択図】図１
【特許請求の範囲】
【請求項1】
組織体積の特性を表すパラメータを特徴付けるための方法であって、トランスデューサを用いて、超音波エコーデータを形成するために応力を受ける前記組織体積に沿った複数の位置によって反射された超音波エネルギを検出することであって、前記組織体積は流体材料と剛体材料とを空間的に分離する組織壁によって形成される、検出すること、
前記超音波エコーデータから波速度時系列データを決定すること、
前記決定された波速度時系列データから過渡的ピークのアンサンブルを生成すること、及び
前記過渡的ピークのアンサンブルから前記組織体積の特性を表すパラメータを決定すること、
を備える方法。
【請求項2】
前記組織体積の特性を表すパラメータは、群速度（ＧＶ）、群速度の二乗（ＧＶ２）、群速度のべき乗、誘発ラム波の到着時間、前記組織体積の推定弾性特性、排尿筋圧（Ｐｄｅｔ）、又は前記組織体積の排尿筋過活動指数の少なくとも一つである、請求項１に記載の方法。
【請求項3】
前記過渡的ピークのアンサンブルを生成することは、下側包絡線データ系列を形成すること、及び前記波速度時系列データから前記下側包絡線データ系列を減算することを含む、請求項１に記載の方法。
【請求項4】
前記組織体積に沿って伝搬する機械的変形の位相速度を計算すること、及び前記組織体積の幾何学的境界によって引き起こされる前記組織体積における前記機械的変形の減衰を考慮するために、反対称ラム波又はせん断波関数を適合させることを更に備える、請求項１に記載の方法。
【請求項5】
前記波速度時系列データから外れ値データ点を除去すること、及び前記除去された外れ値データ点上に前記波速度時系列データを補間することを更に備える、請求項１に記載の方法。
【請求項6】
前記組織壁の表面に沿った前記組織壁の運動に対応する機械的変形を表す分散データを決定するためにスペクトル分析を行うことを更に備える、請求項１に記載の方法。
【請求項7】
前記トランスデューサを用いて前記超音波エネルギを検出することは、（ｉ）音響放射力、（ｉｉ）電気機械入力、及び（ｉｉｉ）機械入力の少なくとも一つを前記組織体積に印加することによって前記応力を生じさせることを含む、請求項１に記載の方法。
【請求項8】
分散データに基づくラム波計算又はせん断波計算の一つを用いて前記組織体積の粘弾性パラメータを推定することを更に備える、請求項１に記載の方法。
【請求項9】
分散データに基づいて、前記組織壁によって生成される圧力に関連するパラメータを決定することを更に備え、前記パラメータは、（ｉ）前記組織体積に沿って伝搬する機械的変形に関連する波の位相速度、（ｉｉ）前記機械的変形に関連する波の群速度、（ｉｉｉ）前記機械的変形に関連する波ピーク速度、（ｉｖ）前記組織壁の弾性、（ｖ）前記組織壁の粘性、（ｖｉ）前記組織体積のコンプライアンス、及び（ｖｉｉ）前記波速度の関数の一つ又は複数を含む、請求項１に記載の方法。
【請求項10】
前記組織体積は被験者の膀胱を含む、請求項１に記載の方法。
【請求項11】
前記組織体積の特性を表すパラメータは、データ系列ｘｅｔｐ及びｘｌｅの数学関数として決定される排尿筋過活動（ＤＯ）指数Ｉであり、ｘｅｔｐ（ｔ）は、前記過渡的ピークのアンサンブル（ＥＴＰ）データ系列の連続信号を表し、ｘｌｅ（ｔ）は、前記下側包絡線（ＬＥ）データ系列の連続信号を表し、各々は時間の関数である、請求項１に記載の方法。
【請求項12】
前記排尿筋過活動（ＤＯ）指数Ｉは、【数1】
【請求項13】
組織体積の特性を表すパラメータを特徴付けるためのシステムであって、超音波エコーデータを形成するために応力を受ける前記組織体積に沿った複数の位置によって反射された超音波エネルギを検出するために構成されたトランスデューサであって、前記組織体積は流体材料と剛体材料とを空間的に分離する組織壁によって形成される、トランスデューサ、及び
前記超音波エコーデータから波速度時系列データを決定し、
前記決定された波速度時系列データから過渡的ピークのアンサンブルを生成し、
前記過渡的ピークのアンサンブルから前記組織体積の特性を表すパラメータを決定する
ように構成されたプロセッサ、
を備えるシステム。
【請求項14】
前記組織体積の特性を表すパラメータは、群速度（ＧＶ）、群速度の二乗（ＧＶ２）、群速度のべき乗、誘発ラム波の到着時間、前記組織体積の推定弾性特性、排尿筋圧（Ｐｄｅｔ）、又は前記組織体積の排尿筋過活動指数の少なくとも一つである、請求項１３に記載のシステム。
【請求項15】
前記プロセッサは、下側包絡線データ系列を形成し、及び前記波速度時系列データから前記下側包絡線データ系列を減算することによって、前記過渡的ピークのアンサンブルを生成するように更に構成される、請求項１３に記載のシステム。
【請求項16】
前記プロセッサは、前記組織体積に沿って伝搬する機械的変形の位相速度を計算し、及び前記組織体積の幾何学的境界によって引き起こされる前記組織体積における前記機械的変形の減衰を考慮するために、反対称ラム波又はせん断波関数を適合させるように更に更に構成される、請求項１３に記載のシステム。
【請求項17】
前記プロセッサは、前記波速度時系列データから外れ値データ点を除去し、及び前記除去された外れ値データ点上に前記波速度時系列データを補間するように更に構成される、請求項１３に記載のシステム。
【請求項18】
前記プロセッサは、前記組織壁の表面に沿った前記組織壁の運動に対応する機械的変形を表す分散データを決定するためにスペクトル分析を行うように更に構成される、請求項１３に記載のシステム。
【請求項19】
前記トランスデューサは、（ｉ）音響放射力、（ｉｉ）電気機械入力、及び（ｉｉｉ）機械入力の少なくとも一つを前記組織体積に印加することによって前記応力を生じさせることによって超音波エネルギを検出するように構成される、請求項１３に記載のシステム。
【請求項20】
前記プロセッサは、分散データに基づくラム波計算又はせん断波計算の一つを用いて前記組織体積の粘弾性パラメータを推定するように更に構成される、請求項１３に記載のシステム。
【請求項21】
前記プロセッサは、分散データに基づいて、前記組織壁によって生成される圧力に関連するパラメータを決定するように更に構成され、前記パラメータは、（ｉ）機械的変形に関連する波の位相速度、（ｉｉ）前記機械的変形に関連する波の群速度、（ｉｉｉ）前記機械的変形に関連する波ピーク速度、（ｉｖ）前記組織壁の弾性、（ｖ）前記組織壁の粘性、（ｖｉ）前記組織体積のコンプライアンス、及び（ｖｉｉ）前記波速度の関数の一つ又は複数を含む、請求項１３に記載のシステム。
【請求項22】
前記組織体積は被験者の膀胱を含む、請求項１３に記載のシステム。
【請求項23】
前記組織体積の特性を表すパラメータは、データ系列ｘｅｔｐ及びｘｌｅの数学関数として決定される排尿筋過活動（ＤＯ）指数Ｉであり、ｘＥＴＰ（ｔ）は、前記過渡的ピークのアンサンブル（ＥＴＰ）データ系列の連続信号を表し、ｘＬＥ（ｔ）は、前記下側包絡線（ＬＥ）データ系列の連続信号を表し、各々は時間の関数である、請求項１３に記載のシステム。
【請求項24】
前記排尿筋過活動（ＤＯ）指数Ｉは、【数2】
【発明の詳細な説明】
【技術分野】
【0001】
この出願は、“Ｓｙｓｔｅｍ ａｎｄ Ｍｅｔｈｏｄ ｆｏｒ Ｎｏｎ－Ｉｎｖａｓｉｖｅ Ｄｅｔｅｒｍｉｎａｔｉｏｎ ｏｆ Ｂｌａｄｄｅｒ Ｏｖｅｒａｃｔｉｖｉｔｙ Ｕｓｉｎｇ Ｕｌｔｒａｓｏｕｎｄ Ｖｉｂｒｏｍｅｔｒｙ”と題された２０２１年１月１５日に出願された米国仮特許出願第６３／１３８，１１６号の利益を主張するものであり、当該仮出願は全ての目的のために参照によってその全体が記載されているのと同様に本明細書に組み込まれる。
連邦政府の支援を受けた研究に関する声明
連邦政府の支援を受けた研究に関する声明
【0002】
この発明は、米国国立衛生研究所から授与されたＤＫ０９９２３１の下での政府支援を受けて行われた。米国政府は、本発明における一定の権利を有する。
【背景技術】
【0003】
本開示は、膀胱過活動を分析するためのシステム及び方法に関する。より具体的には、本開示は、膀胱壁に機械波を励起し、超音波パルスエコー技術を用いてその運動を追跡する、超音波振動測定を用いて膀胱の活動を決定するための方法に関する。
【0004】
膀胱コンプライアンス（排尿筋圧の変化による膀胱容積の変化として定義される）の低下は、神経因性膀胱障害を有する患者の亜集団において特に重要であり、膀胱圧が正常なパラメータの範囲内に収まることを確実にするために１年ごとに尿流動態検査（ＵＤＳ）で注意深く監視されるべきである。尿流動態検査は、現在、膀胱コンプライアンスを測定するための臨床評価における至適基準と考えられているが、不快感と感染のリスクを伴う。典型的な処置は約４５分を要し、膀胱及び膣又は直腸にカテーテルを載置し、所定の速度で膀胱を満たし、膀胱が満たされる際の排尿筋圧として定義される排尿筋圧（Ｐｄｅｔ）の変化を測定することが必要である。コンプライアンス（弾性）膀胱は、充填量に適応して拡張し、その結果、排尿筋圧が低くなるが、非コンプライアンス（非弾性）膀胱は、容易に拡張せず、膀胱充填中に排尿筋圧が上昇する。尿流動態検査は、処置中及び処置後に痛みを伴う場合が多く、尿道カテーテル処置による外傷及び感染症が男女ともに報告されている。不快感及び臨床的リスクに加えて、尿流動態検査は、これらの患者の治療及び経過観察を行う施設にとって労力及び資源を消耗するものである。
【0005】
排尿筋過活動（ＤＯ）とは、排尿の充填期における膀胱の排尿筋の不随意収縮の尿動態的所見である。これらの不随意収縮は、ＵＤＳによって測定されるＰｄｅｔに変動をもたらす。ＤＯは神経因性膀胱患者では一般的であり、上部尿路問題のリスクを伴うため、この集団において診断が重要である。神経因性膀胱患者集団及び非神経因性膀胱患者集団の両方において、ＤＯは、頻尿及び夜間頻尿とともに尿意切迫感を特徴とする過活動膀胱症候群（ＯＡＢ）も伴う場合が多い。ＯＡＢは、この集団の１０％に発生すると考えられており、生活の質の大幅な低下と関連する。
【0006】
下部尿路症状を有する患者におけるＤＯの評価は、臨床的な判断に有用な情報を提供するが、一般に、実質的な臨床上の必要性がある場合にのみ評価される。たとえば、一般に、ＤＯは、非神経因性ＯＡＢの治療において保存療法が失敗した場合にのみ評価されることが推奨される。これは、膀胱圧の測定に従来必要とされる尿路へのカテーテル処置による、ＵＤＳに伴う不快感及び罹患の可能性が理由である。たとえば、９６３人の非神経因性患者について報告された結果を考慮したレビューでは、ＵＤＳ後の尿路感染率は２８％であることが示された。
【0007】
ＤＯの評価に超音波を用いるアプローチは、主に、膀胱の幾何学的特性を測定することに焦点を当てている。たとえば、超音波は、膀胱壁の厚さを測定するために用いられており、膀胱壁の厚さは、排尿筋肥大に起因するＤＯと相関することが提案されている。同様に、収縮中の膀胱形状の変化が超音波撮像によって測定される。そのようなアプローチは、膀胱形状と膀胱力学との相関性に依拠する。これらの技術及び他の利用可能な技術も、非侵襲的で比較的痛みが少ないＤＯの定量的測定を提供しない。それよりも、臨床診療において、臨床医は、しばしば、上述したリスクや欠点にかかわらず、尿流動態検査に頼る。
【0008】
したがって、膀胱機能障害を患う患者において、膀胱の特性に関する有益な情報を非侵襲的かつ無痛で提供するシステム及び方法が必要とされている。
【発明の概要】
【0009】
本開示は、排尿筋過活動（ＤＯ）を決定するために超音波膀胱振動測定（ＵＢＶ）を用いて比較的痛みが少なく非侵襲的に膀胱の粘弾性を決定するシステム及び方法を提供することによって、上述した欠点に対処するものである。本開示は、失禁を経験する患者の早期評価及び高頻度の経過観察のために用いられ得るシステム及び方法を提供することによって、上述した欠点を更に克服するものである。ＵＢＶは、壁にラム波を励起するために集束される超音波放射力と、組織変形を追跡するためのパルスエコー技術とを用いる非侵襲的技術である。この点に関して、“膀胱”という用語は、人体の“泌尿器の”膀胱を指すとは限らず、より一般に、ある程度の弾性を有し得る組織壁によって形成された組織体積を指し得る。したがって、泌尿器膀胱は膀胱の一例である。クロススペクトル分析は、膀胱壁の弾性特性に関連する波速度を計算するために用いられ得る。
【0010】
厳密に形状に基づく膀胱の評価とは異なり、ＵＢＶは、膀胱の弾性及び荷重を特徴付けるより直接的な手段を提供するエラストグラフィ技術である。エラストグラフィ技術は、大きな変形や活発な筋肉収縮によって生じるものを含む、組織荷重によって生じる組織弾性の変化に対する感度を有する。この感度は、心臓の心室、横隔膜筋及び肺、及び大きな血管を含む様々な生物学的導管において組織が維持する圧力を評価及び追跡することが可能な手段として活用され得る。膀胱に関して、ＵＢＶは、ＵＢＶで測定された群速度の二乗（ＧＶ２）、及び弾性と、ＵＤＳで決定されたＰｄｅｔとの中度から強度の相関性を示している。
【0011】
ＵＢＶは、膀胱壁の弾性及び荷重状態に対し高感度であるため、排尿筋圧の時間的変動がＵＢＶ測定によって有用に検出され得る。いくつかの構成において、ＤＯは、ＵＢＶ測定値の時間分解系列における過渡的ピークの分析によって検出され得る。非限定的な例において、定期的なＵＤＳを受ける神経因性膀胱患者の集団がＤＯ測定に関して評価された。タイムスタンプ付きのＧＶ２データ系列が、過渡的ピークデータ系列のアンサンブルと下側包絡線（ＬＥ）データ系列（ＥＴＰ）とに分解された。ＤＯ指数Ｉは、これらのＤＯ関連ピークの特徴付けに用いられた。同じ分析が、比較のために同時に測定されたＰｄｅｔデータ系列に適用された。
【0012】
本開示で説明される一態様によると、超音波振動測定を用いてＤＯを非侵襲的に決定するためのシステム及び方法が開示される。このシステムは、超音波放射力を用いて膀胱壁を“タップ”し、ラム波を励起する。システムはパルスエコー超音波を用いて壁の運動を追跡することを更に含む。波速度分散（周波数の関数としての波速度の変化）を得るために、組織運動のフーリエ空間分析が用いられる。分析的ラム波モデルは、弾性組織体積の粘弾性を推定するために分散データに適合される。その後、被験者に関して、たとえば排尿筋過活動などの組織体積の特性を表すパラメータが決定され得る。
【0013】
一態様において、組織体積の特性を表すパラメータを特徴付けるための方法が提供される。この方法は、トランスデューサを用いて、超音波エコーデータを形成するために応力を受ける組織体積に沿った複数の位置によって反射された超音波エネルギを検出することを含む。組織体積は、流体材料と剛体材料とを空間的に分離する組織壁によって形成される。この方法は、超音波エコーデータから過渡的ピークのアンサンブルを生成すること、及び過渡的ピークのアンサンブルから組織体積の特性を表すパラメータを決定することも含む。
【0014】
この方法のいくつかの態様において、組織体積の特性を表すパラメータは、群速度（ＧＶ）、群速度の二乗（ＧＶ２）、群速度のべき乗、誘発ラム波の到着時間、組織体積の推定弾性特性、排尿筋圧（Ｐｄｅｔ）、又は組織体積の排尿筋過活動指数の少なくとも一つである。
【0015】
一態様において、組織体積の特性を表すパラメータを特徴付けるためのシステムが提供される。このシステムは、超音波エコーデータを形成するために応力を受ける組織体積に沿った複数の位置によって反射された超音波エネルギを検出するために構成されたトランスデューサを含む。組織体積は、流体材料と剛体材料とを空間的に分離する組織壁によって形成される。当該システムは、超音波エコーデータから過渡的ピークのアンサンブルを生成し、及び過渡的ピークのアンサンブルから組織体積の特性を表すパラメータを決定するように構成されたプロセッサも含む。
【0016】
このシステムのいくつかの態様において、組織体積の特性を表すパラメータは、群速度（ＧＶ）、群速度の二乗（ＧＶ２）、群速度のべき乗、誘発ラム波の到着時間、組織体積の推定弾性特性、排尿筋圧（Ｐｄｅｔ）、又は組織体積の排尿筋過活動指数の少なくとも一つである。
【0017】
本発明の上記及び他の態様及び利点は、以下の説明から明らかになる。本明細書において、本明細書の一部を成す添付図面が参照され、図面には、本発明の好適な実施形態が例示的に示される。そのような実施形態は、必ずしも本発明の全範囲を表すものではないので、本発明の範囲を解釈するために、本明細書及び特許請求の範囲が参照される。
【図面の簡単な説明】
【0018】
【図1】超音波放射力システムのブロック図である。
【図5B】排尿筋過活動を決定するための非限定的な方法例の別のフローチャートである。
【図5C】排尿筋過活動を決定するための非限定的な方法例のまた別のフローチャートである。
【図8】本開示の一態様に係る、半無限弾性媒質上の薄層及び流体負荷の図である。
【図9】本開示の一態様に係る、組織弾性を推定するために得られた分散データに適合されたラム波分散方程式のプロットである。
【図11A-D】図１１ＡＤＯ膀胱に関するＧＶ２データ系列の非限定的な例のグラフである。図１１ＢはＤＯ膀胱に関するＧＶ２データ系列の非限定的な例の別のグラフである。図１１Ｃは非ＤＯ膀胱に関するＧＶ２データ系列の非限定的な例のグラフである。図１１Ｄは非ＤＯ膀胱に関するＧＶ２データ系列の非限定的な例の別のグラフである。
【図12A-D】図１２ＡはＤＯ膀胱に関するＰｄｅｔデータ系列の非限定的な例のグラフである。図１２ＢはＤＯ膀胱に関するＰｄｅｔデータ系列の非限定的な例の別のグラフである。図１２Ｃは非ＤＯ膀胱に関するＰｄｅｔデータ系列の非限定的な例のグラフである。図１２Ｄは非ＤＯ膀胱に関するＰｄｅｔデータ系列の非限定的な例の別のグラフである。
【図13A-C】図１３ＡはＰｄｅｔ信号指数に対するＧＶ２信号指数の散布図の非限定的な例である。図１３ＢはＧＶ２データに適用されるＤＯの分類子として信号指数（Ｉ）を用いて生成されたＲＯＣ曲線の非限定的な例である。図１３ＣはＰｄｅｔデータに適用されるＤＯの分類子として信号指数（Ｉ）を用いて生成されたＲＯＣ曲線の非限定的な例である。
【図14A-C】図１４Ａは膀胱例のサンプリング間隔統計値のヒストグラムの非限定的な例である。図１４Ｂは膀胱例のサンプリング間隔統計値のヒストグラムの別の非限定的な例である。図１４Ｃは非限定的な膀胱例の各々に関する実験期間のヒストグラムの非限定的な例である。
【図15】ハンペルフィルタリング前及びハンペルフィルタリング後のデータ系列のグラフの非限定的な例である。
【図16A】ＥＴＰデータ系列を生成するための非限定的な方法例のフローチャートである。
【図16B】生データ系列のグラフ及び処理後のグラフの非限定的な例を示す。
【発明を実施するための形態】
【0019】
排尿筋過活動を決定するために超音波膀胱振動測定（ＵＢＶ）を用いて膀胱の粘弾性を決定するためのシステム及び方法が提供される。ＵＢＶは、組織壁にラム波を励起するために集束される超音波放射力と、組織変形を追跡するためのパルスエコー技術とを用いる非侵襲的技術である。この点に関して、“膀胱”という用語は、人体の“泌尿器の”膀胱を指すとは限らず、より一般に、ある程度の弾性を有し得る組織壁によって形成された組織体積を指し得る。したがって、泌尿器膀胱は、膀胱の一例である。クロススペクトル分析は、膀胱壁の弾性特性に関連する波速度を計算するために用いられ得る。
【0020】
飛行時間法は、壁に沿った組織運動から組織の群速度（ＧＶ）を推定するために用いられてよく、この測定値は、保存及びタイムスタンプされ得る。この方法は、保存及びタイムスタンプされたＧＶ測定値の二乗（ＧＶ２）から過渡的ピークのアンサンブルを生成することも含む。
【0021】
超音波エラストグラフィ法は、様々な軟組織の弾性を測定するために用いられてきた。過去１０年間で、音響放射力（ＡＲＦ）に基づくいくつかの技術が、組織弾性を調査するために開発されてきた。超音速せん断撮像は、ＡＲＦを用いて組織内にせん断波を誘発し、組織のせん断弾性率のマップを生成する。最近では、組織の弾性及び粘性を測定するためのせん断波分散超音波振動計（ＳＤＵＶ）が開発されている。
【0022】
上記技術は一般に、たとえば乳房、肝臓、腎臓、筋肉、及び前立腺などの器官に適用するために、（理論的には境界がない）バルク組織内の波伝搬を仮定して設計されている。また上記技術は、検査される組織が無限の厚さであることも仮定している。しかしながら、この仮定は膀胱には適切でなく、予想される膀胱壁の厚さは最大でも数ミリメートルである。また、膀胱壁は（内部の）液体と（外部の）その他の体腔部分との境界面として機能するため、上述した方法における現行のモデルは、膀胱の独自の特性に適さない。
【0023】
最近の研究により、たとえば心臓壁などのプレート状器官の粘弾性を非侵襲的に定量化するためにラム波分散超音波振動計（ＬＤＵＶ）を用いることが調査されている。ＬＤＵＶは、ＳＤＵＶをベースとした技術であり、機械アクチュエータ又は集束される超音波放射力を用いて関心対象の媒質にラム波を励起し、パルスエコートランスデューサを用いて変形を追跡し、位相速度を推定する。媒質の弾性及び粘性を推定するために、測定されたラム波速度分散（周波数の関数としての位相速度の変化）にラム波分散方程式が適合される。しかしながら、ＬＤＵＶにおいて用いられるラム波分散方程式は、平坦プレートモデルを仮定しており、プレートが両面において水状の流体（尿、血液、心膜液など）に包囲されていることを仮定している。一方、膀胱は、空間的に湾曲した壁（湾曲プレート）を有し、片側には尿を含むが、表層側では結合組織、脂肪、及び腹膜液に包囲されている。したがって、従来用いられているＬＤＵＶ技術を膀胱に容易に適用することはできない。
【0024】
たとえば剛性などの膀胱の機械的特性を測定するための他の方法は、膀胱壁の特性を直接的に評価するバイオセンサを用いることを含む。バイオセンサを用いて行われた剛性測定は、膀胱壁の即時的な生体力学的特性を表すと考えられ得る。この測定値は、膀胱が空であるか、充填しているか、自発的に収縮しているかにかかわらず得られ得る。当然、測定は、バイオセンサが接触している膀胱激領域のみで行われるので、その結果生じる剛性データは、その領域に特有であり、領域の不均質性を検査するために用いられ得る。バイオセンサは膀胱の剛性の直接的な測定を提供するが、バイオセンサは、剛性が測定される膀胱組織と接触するカテーテルに取り付けられるので、患者に対する侵襲的な処置となる。
【0025】
特に図１を参照すると、本開示において説明される方法を実施することができる超音波システム１０００は、刺激点１００４で組織１００６に放射力を伝達する、又はラム波１０１０を付与するために用いられ得る集束ビーム１００２を生成し、そして期間Δτの後にたとえば受信点１００８などからの受信器ビーム１０１２で組織変形を検出する複数の素子を有する超音波トランスデューサ１０１４を含む。それらの素子は、所望のビート周波数によって異なる超音波周波数で、それぞれの連続波シンセサイザによって駆動され得る。集束ビーム１００２は、測定される標的組織１００６を目標としてよく、それに応答して、標的組織は、決定された周波数のラム波１０１０で振動又は発振する。したがってこれらの素子は、所定のビート周波数で標的組織を発振させる力発生器として機能する。
【0026】
標的組織の振動は、超音波システム１０００によって測定され得る。以下で更に詳しく説明するように、超音波システム１０００は、超音波トランスデューサ１０１４を駆動して、標的組織１００６に集束又は非集束超音波ビーム１００２を印加し、標的組織１００６によって反射された受信器ビーム１０１２を有するエコー信号を受信する。これらのエコー信号の位相及び振幅は、たとえば群速度（ＧＶ）、群速度の二乗（ＧＶ２）、群速度のべき乗、誘発されたラム波の到着時間、推定弾性特性、排尿筋圧（Ｐｄｅｔ）、たとえば標的組織１００６の排尿筋過活動指数などの排尿筋活動など、機械的特性や組織特性といった組織体積の特性を表すパラメータを測定するために、後述のように処理され得る。
【0027】
図２を参照すると、トランスデューサアレイ２３は、送信器１３によって生成されたパルスによって通電されると各々が超音波エネルギのバーストを生成する、複数の個別に駆動される素子１１を含む。被検体からトランスデューサアレイ２３に戻り反射する超音波エネルギは、各トランスデューサ素子１１によって電気信号に変換され、スイッチ１５のセットを介して個別に受信器９に印加される。送信器１３、受信器９、及びスイッチ１５は、人間のオペレータによって入力されたコマンドに応答するデジタルコントローラ１９の制御下で動作する。スイッチ１５がそれらの送信位置に設定され、送信器１３が各トランスデューサ素子１１を通電するために瞬間的にゲートオンされ、その後スイッチ１５がそれらの受信位置に設定され、そして各トランスデューサ素子１１によって生成されたエコー信号が受信器９に印加される一連のエコーを取得することによって、スキャンが完了する。各トランスデューサ素子１１からの個別のエコー信号は、受信器９において結合され、表示システム１７上の画像に直線を生成するために用いられる単一のエコー信号が生成される。
【0028】
送信器１３は、生成された超音波エネルギがビーム又はパルスとして方向付け又はステアリングされるようにトランスデューサアレイ２３を駆動する。したがって、トランスデューサアレイ２３を物理的に動かすのではなく、このビームを角度のセットを介して点から点へと動かすことによってＢスキャンが行われ得る。これを達成するために、送信器１３は、連続したトランスデューサ素子１１に印加されるそれぞれのパルス２０に時間遅延（Ｔｉ）を付与する。時間遅延がゼロ（Ｔｉ＝０）である場合、トランスデューサ素子１１は同時に通電され、その結果生じる超音波ビームは、トランスデューサアレイ２３の中心を起点としてトランスデューサ面に垂直な軸３２１に沿って方向付けられる。時間遅延（Ｔｉ）の増加とともに、超音波ビームは、中心軸３２１から角度θだけ下方に方向付けられる。
【0029】
セクタスキャンは、連続する励起において時間遅延Ｔｉを漸次変更することによって行われる。したがって、角度θは、送信ビームを一続きの方向にステアリングするために増分的に変更される。ビームの方向が中心軸３２１より上にある場合、パルス７のタイミングは逆になる。
【0030】
更に図２を参照すると、超音波エネルギの各バーストによって生成されたエコー信号は、超音波ビームに沿って連続した位置（Ｒ）に位置する反射物体から発出する。これらは、トランスデューサアレイ２３の各トランスデューサ素子１１によって個別に感知され、特定の時点におけるエコー信号の大きさのサンプルは、特定の範囲（Ｒ）で発生する反射の量を表す。しかしながら、焦点Ｐと各トランスデューサ素子１１との間の伝搬経路の違いにより、これらのエコー信号は同時に発生せず、それらの振幅は等しくない。受信器９の機能は、これらの個別のエコー信号を増幅及び復調し、各々に適切な時間遅延を付与し、それらを合計して、角度θに向けられた超音波ビームに沿って範囲Ｒに位置する各焦点Ｐから反射した合計超音波エネルギを正確に示す単一のエコー信号を提供することである。
【0031】
各トランスデューサ素子１１からのエコーによって生成された電気信号を同時に合計するために、受信器９の各個別のトランスデューサ素子チャネルに時間遅延が導入される。リニアアレイ２３の場合、各チャネルに導入される遅延は２つの成分に分割することができ、一方の成分はビームステアリング時間遅延と称され、他方の成分はビーム集束時間遅延と称される。受信のためのビームステアリング及びビーム集束時間遅延は、上述した送信遅延と全く同じ遅延（Ｔｉ）である。しかしながら、各受信器チャネルに導入される集束時間遅延成分は、エコー信号が発出する範囲Ｒにおける受信ビームの動的集束を提供するために、エコーの受信中に連続的に変化する。
【0032】
デジタルコントローラ１９の指示の下、受信器９は、受信器９のステアリングが送信器１３によってステアリングされたビームの方向を追跡し、連続した範囲でエコー信号をサンプリングし、ビームに沿って点Ｐに動的に集束するための適切な遅延を提供するように、スキャン中に遅延を提供する。したがって、超音波パルスの各放出は、超音波ビームに沿って位置する対応する一連の点Ｐからの反射音の量を表す一連のデータ点の取得をもたらす。
【0033】
適切な時間遅延を選択することにより、組織のいくつかの点からの振動情報を測定するために複数の集束位置からのエコーが受信され得る。近接位置にある２つの点に関するトランスデューサの横方向分解能の限界は、異なる位置に異なる送信符号を割り当てることによって改善され得る。
【0034】
表示システム１７は、受信器９によって生成された一連のデータ点を受信し、所望の画像を生成する形式にデータを変換する。たとえば、Ａスキャンが所望される場合、一連のデータ点の大きさは、時間の関数として単にグラフ化される。Ｂスキャンが所望される場合、系列内の各データ点は、画像内の画素の輝度を制御するために用いられ、そして画像の表示のために用いられ得るデータを提供するために、連続的なステアリング角度（θ）における一連の測定から成るスキャンが行われる。
【0035】
特に図３を参照すると、送信器１３は、集合的に９００で示されるチャネルパルス符号メモリのセットを含む。各パルス符号メモリ９００は、生成される超音波パルス９０４の周波数を決定するビットパターン９０２を格納する。このビットパターンは、マスタクロックによって各パルス符号メモリ９００から読み出され、信号をトランスデューサ２３の駆動に適した電力レベルに増幅するドライバ９０６に印加される。図３に示す例において、ビットパターンは、５メガヘルツ（“ＭＨｚ”）の超音波パルス９０４を生成するための４つの“１”ビットと４つの“０”ビットとが交互になったシーケンスである。これらの超音波パルス９０４が印加されるトランスデューサ素子２３は、超音波エネルギを生成することによって応答する。
【0036】
上述したように、所望の方法で超音波エネルギの送信ビームをステアリングするために、Ｎ個のチャネルの各々に関してパルス９０４が生成され、そして適切な量だけ遅延される。これらの遅延は、デジタルコントローラ１９からの制御信号を受信する送信制御装置９０８によって提供される。制御信号が受信されると、送信制御装置９０８は、クロック信号を第１の送信チャネルにゲートする。その後、連続した遅延時間間隔の各々において、クロック信号は、電通される全てのチャネルが超音波パルス９０４を生成するまで、次のチャネルパルス符号メモリ９００にゲートされる。各送信チャネルは、その全ビットパターン９０２が送信された後、リセットされてよく、その後、送信器１３は、デジタルコントローラ１９からの次の制御信号を待機する。
【0037】
特に図４を参照すると、受信器９は、時間利得制御部１００、ビーム形成部１０１、及び中間プロセッサ１０２の３つの部分で構成される。時間利得制御部１００は、Ｎ＝１２８個の受信器チャネルごとの増幅器１０５と、時間利得制御回路１０６とを含む。各増幅器１０５の入力は、受信するエコー信号を受信及び増幅するためにトランスデューサ素子１１のそれぞれ一つに接続される。増幅器１０５によって提供される増幅の量は、時間利得制御回路１０６によって駆動される制御ライン１０７を介して制御される。当該技術において周知であるように、エコー信号の範囲が増大すると、その振幅は減少する。その結果、より遠くの反射体から発出するエコー信号が近くの反射体からのエコー信号よりも多く増幅されない限り、画像の輝度は、範囲（Ｒ）の関数として急速に減少する。この増幅は、ＴＧＣリニアポテンショメータ１０８をセクタスキャンの全範囲にわたり比較的均一な輝度を提供する値に手動で設定するオペレータによって制御される。エコー信号が取得される時間間隔により、エコー信号が発出する範囲が決定され、この時間間隔は、ＴＧＣ制御回路１０６によって分割される。ポテンショメータの設定は、エコー信号が取得時間間隔にわたり増加し続ける量で増幅されるように、それぞれの時間間隔の間に増幅器１０５の利得を設定するために用いられる。
【0038】
受信器９のビーム形成部１０１は、Ｎ＝１２８個の個別の受信器チャネル１１０を含む。各受信器チャネル１１０は、入力１１１においてＴＧＣ増幅器１０５の一つからアナログエコー信号を受信し、そしてＩバス１１２及びＱバス１１３においてデジタル化出力値のストリームを生成する。これらのＩ値及びＱ値の各々は、特定の範囲（Ｒ）におけるエコー信号エンベロープのサンプルを表す。これらのサンプルは、合計点１１４及び１１５において他の受信器チャネル１１０の各々からのＩ及びＱサンプルと合計されると、ステアドビーム（θ）上の範囲Ｒに位置する点Ｐから反射したエコー信号の大きさ及び位相を示すように、上述した方法で遅延されている。
【0039】
図４を更に参照すると、中間プロセッサ部１０２は、合計点１１４及び１１５からのビームサンプルを受信する。各ビームサンプルのＩ値及びＱ値は、点（Ｒ、θ）からの反射音の大きさの同相成分及び直交成分を表す１６ビットのデジタル数値である。中間プロセッサ１０２は、これらのビームサンプルに様々な計算を行うことができ、再構成される画像の種類によって選択肢が決定される。
【0040】
たとえば、従来の超音波画像は、エコー信号の大きさをそのＩ成分及びＱ成分
から計算する検出プロセッサ１２０によって生成され得る。
【数1】
【0041】
その結果生じる大きさ値は１２１において表示システム１７に出力され、それによって、各画素において反射したエコーの大きさが示された画像が生じる。
【0042】
本発明の方法は、中間プロセッサ１０２の一部を形成する組織量プロセッサ又は特性プロセッサ１２２の特性を表すパラメータによって実施され得る。以下で詳しく説明するように、このプロセッサ１０２は、対象組織の一連の測定中に取得されたＩ及びＱビームサンプルを受信し、組織の機械的特性（すなわち厚さ、活動など）を計算する。
【0043】
本開示によると、ＵＢＶ測定における過渡的ピークの分析によってＤＯを検出するための方法が提供される。ピークは、これらの測定において、ピークと下側包絡線（ＬＥ）との間で異なる信号エネルギを特徴付けるＤＯ指数（Ｉ）によって特徴付けられ得る。提供される非限定的な例において、ＤＯ膀胱と非ＤＯ膀胱との間で統計的に有意なＩの差が観察された。
【0044】
いくつかの構成において、ＧＶ２データ系列は、たとえば神経因性膀胱症状を有する被験者などの被験者から収集され得る。取得されたデータは、データ系列をＥＴＳ系列とＬＥ系列とに分解することによってＤＯ関連の過渡的ピークを識別するために処理され得る。分解された系列を単一の数値によって特徴付けるために、ＤＯ指数Ｉが決定され得る。
【0045】
ＤＯに関連するＧＶ２測定系列における変動を特徴付けるために、そのような変動は、膀胱充填に関連する任意のより大きな全体傾向から分離され得る。膀胱壁滑動の結果である可能性が低く、たとえば大きな外れ測定値のように特徴付けを左右し得る変動は除去され得る。後者を考慮する場合、統計的フィルタリングアプローチ（たとえばハンペルフィルタリング）が用いられ得る。系列の主な傾向からＤＯ関連の変動を分離するために、系列は、それぞれＬＥ系列とＥＴＰ系列とに分解され得る。
【0046】
ピーク識別の前に、膀胱の漏れが報告されたＧＶ２の値は、所与のデータ系列から除去され得る。排尿の開始に伴う正常な排尿筋収縮が含まれることを避けるために、実験の最終量での測定値もピーク識別から除外され得る。そのような外れ値は、たとえば、測定中の患者の動きや異常に弱いＵＢＶ信号に起因し得る。これらの外れ値を除去するために、たとえばハンペルフィルタなどの統計的外れ値フィルタがＧＶ２データ系列に適用され得る。このフィルタは、外れ値であると決定された点をサンプルウィンドウの中央値に置き換えることによって機能する。フィルタは、
である場合、中央値ｍｉを有するスライディングウィンドウの中心点であるデータ系列ｘｉからの所与のサンプルなどの外れ値を識別してよい。
【数2】
【0047】
ここで、ｎσは、標準偏差の閾値数であり、σｉは、サンプリングウィンドウの推定標準偏差である。サンプリングウィンドウの標準偏差は、以下のように中央絶対偏差から推定され得る。
【数3】
【0048】
ここで、ｗは、ウィンドウの中心点の両側にある画素の数であり、パラメータｋは、中央絶対偏差を標準偏差にスケーリングする。
【0049】
外れ値フィルタリングの後、結果として生じるデータ系列は、ＥＴＰ系列とＬＥ系列とに分解され得る。膀胱充填に関連するデータ系列内の任意の傾向を捕捉するために、生データ系列（ｘ）に平滑化フィルタが適用され得る。平滑化は、たとえば２次Ｓａｖｉｔｚｋｙ－Ｇｏｌａｙフィルタによって系列に適用され得る。Ｓａｖｉｔｚｋｙ－Ｇｏｌａｙフィルタは、スライディングウィンドウ内のデータに（２次）多項式回帰を適用する。時間間隔は不規則な間隔を有し得るため、離散ウィンドウサイズではなく平滑化係数が用いられ得る。平滑化係数は、０～１の間の値をとるように選択されてよく、平滑化によって減衰する信号エネルギの割合に対応する（たとえば、平滑化係数が０．５の場合、信号エネルギの５０％が減衰する）。平滑化係数は、データの主な傾向を除くデータ系列内の全ての変動をフィルタリングするために、たとえば０．９５などの強い値に設定され得る。
【0050】
平滑化の後、ｘは、フィルタリングされたデータ系列（ｘｓ）と差分されてよく、そして差分系列の負の値に関連する時点（すなわち、（ｘ－ｘｓ）＜０である時点）が識別され得る。時間境界における増加が過渡的ピークとして扱われないように、開始時点及び終了時点も識別された時点に含まれ得る。
【0051】
その後、識別された時点におけるｘからの値が収集され得る。識別されなかった時点に値を割り当てるために、線形補間が用いられ得る。その結果生じるＬＥ系列（ｘＬＥ）は、ｘと同じ要素数を有してよい。その後、ＥＴＰ系列（ｘＥＴＰ）は、ｘからｘＬＥを差分することによって決定され得る。その後、ＥＴＰ系列及びＬＥ系列は、更なる分析のために用いられ得る。
【0052】
診断メトリック
【0053】
ピークに対するＤＯ膀胱と非ＤＯ膀胱との比較を行うために、識別された過渡的ピークを特徴付ける単一のメトリックが用いられ得る。また、そのようなメトリックは、原理的に、上述したＵＢＶ測定の可変の実験期間及び範囲を考慮すると、実験ごとに比較可能でなくてはならない。排尿筋過活動（ＤＯ）指数と呼ばれるこのメトリックは、Ｉで示される。
【数4】
【0054】
ここで、ｘＥＴＰ（ｔ）は、ＥＴＰデータ系列の連続信号を示し、ｘＬＥ（ｔ）は、ＬＥデータ系列の連続信号を示し、各々は、所与の実験に関する時間の関数である。基本的に、Ｉは、（それぞれＥＥＴＰ及びＥＬＥと示される）過渡的ピーク及び下側包絡線からの信号エネルギの比の平方根である。よって、ＤＯを有する膀胱からのデータ系列は、ｘＬＥ（ｔ）に比べてｘＥＴＰ（ｔ）における信号エネルギが大きいので、ＤＯを有さない膀胱よりもＩが大きいことが予想される。Ｉは単位がないため、ＲＭＳ（ｘＥＴＰ）をＲＭＳ（ｘＬＥ）で正規化することによって、より長い実験期間、及び測定されるＧＶ２又はＰｄｅｔにおけるより広い範囲が考慮されるので、実験間の比較を行うことが合理的である。
【0055】
図５Ａを参照すると、排尿筋過活動を決定するための非限定的な方法例４００のフローチャートが示される。ステップ４０２において、システムは、たとえば指定された期間、次の測定値の取得を待機してよい。ステップ４０４において、音響放射力励起信号が形成されてよく、そしてステップ４０６において、関心対象の組織に振動パルスが印加され得る。ステップ４０８において、運動信号が決定されてよく、ここで、連続的なＢモードによって組織運動を追跡するために、たとえば膀胱壁である伝搬線に沿ったいくつかの点で運動を測定するためにパルスエコー超音波（検出ビーム）が用いられる。
【0056】
ステップ４１０において、時間及び距離の関数としてインパルスの伝搬距離を測定することによって群速度が決定される。その直線の傾きが群速度である。群速度を測定するための方法の非限定的な例は、パルスのピークが既知の距離を移動するために要する時間を測定することを含む。これは、組織内のパルス伝搬を監視することによって行われ得る。監視は、関心対象領域のＢモード超音波における連続的なＡライン（無線周波数エコー信号）を相関させることによって行われ得る。パルスが、その原点（プッシュビームが作用した点）から、原点から既知の距離にある点まで移動するのに要する時間が測定される。別の非限定的な例において、パルスが物体内の一点から別の点まで通過するのに要する時間を測定することができる。より正確な測定値を得るために、複数の点でプロセスを繰り返すことができ、そして、群速度は、距離対時間のプロットにおける回帰直線の傾きとして推定され得る。ステップ４１２において、群速度はＧＶ２データ系列を形成するために二乗され、タイムスタンプされ、更なる処理のために保存される。
【0057】
生ＧＶ２データ系列のタイムスタンプされたデータ系列は、ステップ４１４において最後の時点に到達した後、ステップ４１６において、印加された励起から収集又は集計され得る。ステップ４１８において、生データをフィルタリングすることによって、平滑化されたデータ系列が収集され得る。ステップ４２０において、選択時点が識別され、系列から除外され得る。ステップ４２２において、除外された時点において系列が補間され得る。ステップ４２４において、下側包絡線（ＬＥ）データ系列が形成され得る。ステップ４２６において、生データ系列からＬＥ系列が減算され得る。その後、ステップ４２８において、ＥＴＰデータ系列からのエネルギが生成され得る。
【0058】
図５Ｂを参照すると、排尿筋過活動を決定するための非限定的な例示的方法４３０のフローチャートが示される。ステップ４３２において、システムは、たとえば指定された期間、次の測定値の取得を待機してよい。ステップ４３４において、音響放射力励起信号が形成されてよく、そして、ステップ４３６において、関心対象の組織に振動パルスが印加され得る。ステップ４３８において、運動信号が決定されてよく、ここで、連続的なＢモードによって組織運動を追跡するために、たとえば膀胱壁である伝搬ラインに沿ったいくつかの点で運動を測定するためにパルスエコー超音波（検出ビーム）が用いられる。ステップ４４０において、群速度が決定される。ステップ４４２において、群速度はＧＶ２データ系列を形成するために二乗され、タイムスタンプされ、そして、ステップ４４４において最後の時点に到達すると、更なる処理のために保存される。最後の時点に到達していない場合、方法はステップ４３４で再開し、音響放射力励起信号を再形成してよい。
【0059】
図５Ｃを参照すると、排尿筋過活動を決定するための非限定的な例示的方法４５０が示される。ステップ４５２において、システムは、たとえば指定された期間、次のＵＢＶ測定値の取得を待機してよい。ステップ４５４において、音響放射力励起信号が形成されてよく、そして、ステップ４５６において、関心対象の組織に振動パルスが印加され得る。ステップ４５８において、運動信号が決定されてよく、ここで、たとえば連続的なＢモードによって組織運動を追跡するために、たとえば膀胱壁である伝搬ラインに沿ったいくつかの点で運動を測定するためにパルスエコー超音波（検出ビーム）が用いられる。ステップ４６０において、群速度が決定される。ステップ４６２において、位相速度分散データが決定され得り、そして反対称ラム波モデルに適合されてよく、これは、ステップ４６４において最後の周波数に到達すると、更なる処理のために保存され得る。最後の周波数に到達していない場合、方法は、ステップ４７０において後続の振動周波数を決定すること、及びステップ４５８において新たな運動信号を決定するために再開することを含んでよい。ステップ４６４において最後の周波数に到達する場合、ステップ４６６において機械的特性が決定及び保存され得る。ステップ４６８において最後の時点に到達する場合、方法は終了してよく、そうでない場合、方法はステップ４５２において再開し、次のＵＢＶ信号を待機及び取得してよい。
【0060】
ここで図５Ｄを参照すると、組織２１の一つ又は複数の機械的特性を測定するための方法の例示的ステップ５００を記載するフローチャートが提供される。プロセスはプロセスブロック５０２で音響放射力励起信号の形成から始まり、プロセスブロック５０４で、被験者に印加され、それによって対象の組織に振動パルスが印加される。図６を参照すると、超音波膀胱振動測定（ＵＢＶ）の基本的原理の図が示される。ここでも、“膀胱”という用語は、人体における“泌尿器の”膀胱を指すのではなく、より一般に、ある程度の弾性を有し得る組織壁で形成された組織体積を指し得る。ＵＢＶは、放射力６０４（プッシュビーム）を生成して、関心対象の媒質にインパルスラム波（長さ２００～６００μｓ）を励起するために、集束超音波を用いる。放射力励起６０４は、たとえば６００μｓのトーンバーストであってよい。
【0061】
たとえば膀胱壁などの膜状媒質に印加される集束超音波放射力は、円筒波を励起する。尿で充填した膀胱は、たとえば非圧縮性で非粘性の流体に浸漬された非圧縮性で均質な等方性の固体としてモデル化され得る。そのような幾何学に関する反対称ラム波分散方程式（２）は以下の通りである。
【数5】
【0062】
ここで、
は流体中の圧縮波の波数であり、
は組織中の圧縮波の波数であり、
はラム波数であり、ωは角周波数であり、
は各周波数であり、
は周波数依存性ラム波速度であり、
はせん断波数であり、
は密度であり、ｈは膀胱壁の半分の厚さ（Ｈ＝２ｈ）である。組織、尿、及び血液の場合、圧力波数（ｋｐ及びｋｆ）が同様であり、せん断波及びラム波（ｋｓ）及び（ｋＬ）より大幅に小さいので、
である。組織、尿、及び血液の密度が類似していることにより、分散方程式（２）は、以下のように簡略化される。
【数6】
【数7】
【数8】
【数9】
【数10】
【数11】
【数12】
【数13】
【数14】
【0063】
あるいは、膀胱内でせん断波を励起するために、外部機械振動器（不図示）が用いられ得る。外部機械振動器を用いると、超音波放射力は印加されず、これらの例では、超音波は、膀胱内にせん断波を発生させるためではなく膀胱壁における波を検出するために用いられる。機械振動器は、皮膚に機械エネルギを伝達するプラスティックシャフトを有するスピーカのような機構を含んでよい。機械エネルギは、腹部及び膀胱を振動させる。
【0064】
再び図５Ｄ及び図６を参照すると、プロセスブロック５０６において、連続的なＢモードによって組織運動を追跡するために、たとえば膀胱壁である伝搬ライン６０８に沿ったいくつかの点で運動を測定するために検出ビーム６０６（たとえばパルスエコー超音波）が用いられる。以下の説明は、膀胱壁に関して記載されるが、他の弾性組織体積も同様に考慮され得る。Ｂモードスキャンは、膀胱壁に４００～６００μｓのインパルスを励起し、機械運動を追跡するために、たとえば図６Ａに示すようなＣ４－２リニアアレイトランスデューサ６０２を備えたＶｅｒａｓｏｎｉｃｓ（登録商標）超音波撮像プラットフォームによって行われ得る。Ｖｅｒａｓｏｎｉｃｓは、ワシントン州Ｒｅｄｍｏｎｄ所在のＶｅｒａｓｏｎｉｃｓ社の登録商標である。検出ビーム６０６は、約３．０ＭＨｚ（特定の例において３．３ＭＨｚ）に周波数の中心を有する、たとえば約２．０ｋＨｚ（別の例では２．５ｋＨｚ）のパルス繰返し周波数（ＰＲＦ）で送信され得る。図６Ｂは、ＵＢＶ励起ビーム６３０のスケッチによってヒトの膀胱の生体内Ｂモードを示す。
【0065】
次に、膀胱壁運動を時間の関数として計算するために、受信エコーのクロススペクトル分析が用いられ得る。インパルスは主にインパルスの時間長の逆数までの周波数成分を含むので、単一のインパルスプッシュから多数の位相速度が抽出され得る。プロセスブロック５０８に示すように、膀胱壁運動の２次元高速フーリエ変換（２Ｄ－ＦＦＴ）は、周波数の関数としてのラム波速度の変化、又はラム波分散を計算するために用いられ得る。言い換えると、運動のフーリエ空間分析は、座標が周波数ｆ及び波数ｋであるｋ空間が得られる時間の関数として用いられ得る。波速度はｃ＝ｋ／ｆであるため、各周波数における位相速度は、プロセスブロック５１０において、所与の周波数におけるピークを探索し、所与のピークの周波数座標による波数座標で割ることによって計算され得る。プロセスブロック５１８に示すように膀胱粘弾性を推定するために、プロセスブロック５１２において、ラム波分散方程式が分散データに適合され得る。プロセスブロック５１４によって示すように、デジタルコントローラは、高調波励起アプローチが用いられる場合に最後の周波数が測定されたかを決定する。否である場合、プロセスブロック５１６において、別の周波数が選択され、各所望の所定周波数でプロセスブロック５０６～５１２が繰り返される。
【0066】
上述したように、ラム波モデルは、膀胱が両側を流体で包囲されていることを仮定しており、状況によっては、膀胱の機械的特性を定量化するために十分ではない場合がある。本開示の別の態様によると、膀胱壁が片側を尿に、他方側を軟質結合組織及び筋肉に包囲されているという点を考慮するためのモデルが導出され得る。
【0067】
図８を参照すると、片側が半無限流体８０４に、他方の側が半無限固体８０６に包囲されたプレート又は薄層８０２を含むラム波伝搬の分析モデルが作成され得る。このモデルは、膀胱壁の自然な環境を模倣しているために有利であり得る。このモデルは、図８に示すように、片側が剛性半空間と接触し、上面に流体が結合した、厚さｈの薄層８０２を含む。波運動のポテンシャル関数は、以下のように記述され得る。
【数15】
【0068】
ここで、
及びｋは波数であり、下付き文字ｐは一次波又は圧縮波を示し、下付き文字ｓはせん断波を示し、ｆは流体を示す。以下の境界条件は、ｘ２＝ｈで適用され得る。
【数16】
【数17】
【0069】
及び層とｘ２＝０における剛性半空間との間の境界での以下の境界条件：
【数18】
【0070】
（４）のポテンシャル関数を用いて上記の境界条件を適用すると、ｘ２＝ｈにおける硬質基板での平面応力は、
とすると、次のように記述され得る。
【数19】
【数20】
【0071】
ｘ２＝０におけるせん断応力は、
である。
【数21】
【0072】
ｘ２＝ｈにおけるせん断応力は、
である。
【数22】
【0073】
ｘ２＝ｈにおける境界条件によるｘ１方向に沿った変位：
【数23】
【0074】
ｘ２＝０における流体負荷による基板の平面応力は、
である。
【数24】
【0075】
ｘ２＝０における流体負荷によるｘ２方向への変位は、
である。
【数25】
【0076】
ｘ２＝ｈにおける流体負荷によるｘ２方向への変位は、
である。
【数26】
【0077】
これらの方程式を方程式系に組み立てると、下の行列Ａとなる。非自明の結果を得るために、Ａの行列式はゼロに設定される（ｄｅｔ（Ａ）＝０）。これらの結果を用いて、特定の周波数範囲に関する分散曲線がプロット化され得る。
【数27】
【0078】
分散方程式（１４）は、方程式（２）に関して説明した上述のラム波分散と同様に、膀胱壁の機械的特性を定量化するために用いられ得る。行列Ａの行列式は、ゼロに設定され、数値的に解くと、所望の分散曲線が得られる。ＵＢＶにおける分散行列Ａの使用により、膀胱壁を包囲する媒質が考慮され、放射力に起因する膀胱壁での波伝搬に影響を及ぼす応力及び歪みをより完全に記述することができる。
【0079】
また、膀胱組織は、幾何学的減衰（この減衰は幾何学的形状の境界によって生じる）と、材料粘性による減衰（この減衰は媒質が無限であっても生じる）との両方を経験し得る。幾何学的減衰は、粘性がない場合にもプレートにおける大きな分散を示す反対称ラム波モデルを用いることによって考慮され得る。粘性による減衰は、組織の弾性及び粘性を二つの個別の因子として記述するＶｏｉｇｔモデルを用いることによって考慮される。膀胱壁の場合、減衰が大きく、励起点から数センチメートル以内に波が消滅し得るので、波の反射及び巻き返しの可能性を無視し、上述したように、膀胱を無限プレートとしてモデル化する。膀胱組織はＶｏｉｇｔ材料として挙動すると仮定され得るので、せん断弾性率はμ＝μ１＋ｉωμ２であり、ここでμ１及びμ２は、それぞれ組織の弾性及び粘性である。Ｖｏｉｇｔモデルの代わりに、組織の機械的特性を特徴付けるための他のレオロジーモデル及びモデルのないアプローチが用いられてよい。組織は、弾性又は粘弾性のいずれかであると仮定され得る。
【0080】
ただし、上述したように、Ｖｏｉｇｔモデルの代わりに、たとえばマックスウェル、一般化マックスウェル、ツェナーなどの他のレオロジーモデルが用いられ得る。方程式３は、プロセスブロック５１０に示すラム波分散速度を得るために適合されてよく、この例では、入力は、ラム波速度対周波数と超音波システムのＢモードスキャンから推定された壁厚さ（ｈ）である。また、Ｂモード画像の限られた空間分解能により、ラム波分散方程式（２）で分散データを適合させるために用いられる壁厚さ（ｈ）推定値には、ある程度の誤差がある。しかしながら、壁厚さ（ｈ）測定値の誤差によるラム波分散速度の誤差は、０．２ｍ／ｓ以下、又は約１０％である。あるいは、壁厚さ（ｈ）測定値の精度を高めるために、より高い周波数の超音波が用いられてよい。
【0081】
図５Ｄのプロセスブロック５１２を特に参照すると、図７Ａ及び図７Ｂに示す平坦な粘弾性プレート７００及び湾曲した粘弾性プレート７０２によってモデル化された膀胱壁におけるラム波伝搬の有限要素分析は、ラム波分散に対する曲率の影響を調査するために設計され得る。固体粘弾性プレート７００、７０２の２次元軸対称有限要素モデル（ＦＥＭ）は、ＡＢＡＱＵＳ６．８－３（ＳＩＭＵＬＩＡ、Ｐｒｏｖｉｄｅｎｃｅ、ＲＩ）を用いて設計され得る非圧縮性の尿及び血液を模倣した流体に浸漬され得る。プレート７００、７０２を包囲する流体は、たとえば２．２ＧＰａの体積弾性率及びたとえば１ｇ／ｃｍ３の密度を有する音響要素（ＡＣＡＸ８）によって表され得る。小さな運動の影響下にある非粘性流体を左右する平衡方程式は、
である。
【数28】
【0082】
ここで、ｐは動圧、ｘは空間位置、
及び
は、それぞれ流体粒子の速度及び加速度、ｐは流体の密度である。γは、速度ごとの体積当たりの力単位での“体積抗力”と定義される。動圧ｐ、体積弾性率Ｋ、及び体積ひずみεを結び付ける流体の構成的挙動は、以下のように公式化される。
【数29】
【数30】
【数31】
【0083】
音響媒質及び固体部分は、“タイ”境界条件と結合され、たとえば１．５×１０６Ｎ・ｓ・ｍ－３の音響インピーダンスが境界面に定められる。
【0084】
平坦プレート７００及び湾曲プレート７０２の両者は粘弾性であり、たとえば１．０ｇ・ｃｍ３の密度、０．４９９のポアソン比、Ｅ＝２５ｋＰａのヤング率、及び、たとえばｇ１＝０．９９０１及びτ１＝２．４ｘ１０－６ｓである２パラメータＰｒｏｎｙ級数に関して定義される機械的特性を有してよい。Ｐｒｏｎｙ級数は、一般化マックスウェルモデルに容易に関連付けることができ、Ｖｏｉｇｔモデルに更に関連付けられるので、選択された機械的特性により、たとえばμ１＝８．３３ｋＰａ及びμ２＝２Ｐａ・ｓのＶｏｉｇｔ材料が推定される。プレート７００、７０２は、ＡＢＡＱＵＳにおける８ノード四次軸対称要素（ＣＡＸ８Ｈ）でメッシュ化され得る。
【0085】
平坦プレート７００は、たとえば厚さ２ｍｍ、長さ６０ｍｍであってよい。湾曲プレート７０２は、たとえば厚さ２ｍｍ、４ｃｍの内径ｒ、及び曲率９０度であってよい。図７Ａ及び図７Ｂに示す両方のモデルにおいて、プレート７００、７０２の遠位端７０４は、剛体運動を避けるために一定の境界条件を有してよく、プレート７００、７０２は、反射が生じないように十分大きくてよい。両方のモデルにおいて、運動を励起するために、たとえばプレートの厚さＨを通る線源から５００μｓにわたり１００μｍの振幅を有するステップ関数として実施されるインパルス変位が用いられ得る。励起７０６の位置は、図７Ａ及び図７Ｂに示す。上述したように、モデルは、ＡＢＡＱＵＳの動的陽解法ソルバで解くことができる。
【0086】
変位は、たとえば合計２０ｍｓにわたり１０ｋＨｚのサンプリングレートで、たとえば図７Ａ及び図７Ｂに示すようにｘ軸に沿って０．５ｍｍごとに測定され得る。平坦プレート７００において、たとえばサンプルの１１の等間隔の深さでｚ軸に沿って示される、励起７０６に平行な変異が記録され得る。湾曲プレート７０２において、変位は、たとえばｚ軸に沿って示すように励起７０６に平行、かつｒ軸に沿って示すようにサンプルの弧に垂直な１１の等間隔の深さで記録され得る。中間深さ線７０８は、分析のために用いられ得る。両方のプレート７００、７０２に関する変位対時間の、図５Ｄのプロセスブロック５０８に示すフーリエ空間分析が、ラム波分散曲線を得るために用いられ得る。ＦＥＭシミュレーションに機械的特性が挿入されたラム波分散方程式（２）が比較のためにプロット化され得る。ＵＢＶ技術は、曲面における軸方向の運動を測定するので、励起７０６及びリニアトランスデューサ要素（不図示）に平行な運動成分が検出され得る。変位ベクトル（ｄｘ＝ｒｄθ）における曲率の補正がプレート７０２の曲率を考慮する上で十分であるかを決定するために、励起に平行な運動と弧に垂直な運動との比較が行われ得る。
【0087】
ここで、図５Ｄのプロセスブロック５０８及び５１０を特に参照すると、上述したように、インパルスは周波数成分を含む。短い有限長のトーンバーストは、トーンバーストの持続時間の逆数までの周波数成分を生成し、これにより、単一のプッシュから多数の位相速度を抽出することが容易になり得る。インパルスは、個々の周波数位相速度の速度とは異なる群速度で伝搬することができ、これについては後述される。分散ダイアグラムは位相速度に関連するので、いくつかの実装において、インパルス変位データから各周波数で位相速度が抽出され得る。よって（６）を用いて変位対時間データに２次元高速フーリエ変換（２ＤＦＦＴ）が行われ得る。
【数32】
【0088】
励起ビームからの距離（ｘ）及び時間（ｔ）の関数として、図７Ａ及び図７Ｂに示すように、ｕｚ（ｘ，ｔ）は、励起ビーム７０６に垂直な膀胱壁の運動である。ここでｋは波数であり、ｆは波の時間周波数である。ｋ空間の座標は、波数（１／λ）及び周波数（ｆ）である。多くのエネルギを運ぶ最も低次の伝搬モードは、ゼロ次反対称（Ａ０）ラム波モードである。Ａ０モードのラム波位相速度分散は、各周波数に関して最大の大きさを有する画素を発見し、そして、波速度ｃ＝λｆであるため、ｋ空間の波数座標（１／ｍ）で周波数座標（Ｈｚ）を割ることによって計算され得る。このように得られた分散データが図９に示される。その後、平坦プレート７００に関するラム波分散方程式（２）は、膀胱の弾性及び粘性並びにＢモードによる膀胱の推定厚さを測定するために、分散データに適合され得る。湾曲プレート７０２の場合、変位ベクトルは、弧の長さｘ＝ｒθを認識することによって構成されてよく、ｒは曲率半径であり、θは角度である。この関係により、湾曲プレート７０２の曲率が補正される。
【0089】
あるいは、群速度は、膀胱の特性を評価するために用いられてもよい。この方法は、位相速度を用いて膀胱の特性を評価するほど正確ではない場合があるが、上述したように、この方法の利点は、測定及び計算の単純さである。たとえば、膀胱の特性を評価するための群速度アプローチは、上述した調和励起及び放射力インパルスアプローチによって用いられるような、図５Ｄのプロセスブロック５０８に示すフーリエドメイン分析を必要としない。群速度は、時間及び距離の関数としてインパルスの伝搬距離を測定することによって計算される。その直線の傾斜が群速度である。より具体的には、群速度を測定する方法の一つは、パルスのピークが既知の距離を移動するのに要する時間を測定することである。これは、組織におけるパルス伝搬を監視することによって行われ得る。監視は、関心領域のＢモード超音波における連続的なＡライン（無線周波数エコー信号）を相関付けることによって行われ得る。パルスが、その原点（プッシュビームが作用した点）から、原点から既知の距離にある点まで移動するのに要する時間が測定される。あるいは、パルスが物体内の一つの点から別の点まで通過するのに要する時間を測定してもよい。より正確な測定値を得るために、プロセスは複数の点で繰り返されてよく、群速度は、距離対時間のプロットにおける回帰直線の傾斜として推定され得る。
【0090】
群速度は、帯域幅の中心周波数の位相速度に等しいので、膀胱壁の弾性に関連する。群速度（ｃｇ）は、以下の式によって弾性（μ）に関連付けられ得る。
【数33】
【0091】
ここで、ｐは組織密度であり、たとえば水の密度であることが仮定される。せん断波速度とラム波速度との不一致を補正するために、１．２の補正係数が足され得る。
【0092】
非限定的な例
【0093】
ＵＢＶ測定の分析によってＤＯを検出することの実現性を調査するために、本明細書で説明される方法が用いられる。ＧＶ２及びＰｄｅｔデータ系列は、神経因性膀胱患者が同時のＵＢＶ及びＵＤＳ測定を受けた実験から収集された。これらのデータ系列に処理が適用され、データ系列をＥＴＳ系列とＬＥ系列とに分解することによってＤＯ関連の過渡的ピークが識別された。最後に、分解された系列を単一の数値で特徴付けるために、ＤＯ指数Ｉが決定された。
【0094】
患者集団
【0095】
非限定的な例は、神経因性膀胱を有する７６人の成人患者からのデータを含む。その結果であるＵＤＳチャートからの泌尿器科医の診断により、患者のうち２３人がＵＤＳ測定中にＤＯを示したことが決定された。
【0096】
ＵＤＳ及びＵＢＶ測定
【0097】
ＵＢＶ測定値は、患者の定期的な膀胱内圧測定及びＵＤＳ中に収集された。膀胱内圧測定中、圧力センサを備えたカテーテルを通して膀胱内に生理食塩水が注入された。膀胱の排尿筋によって維持される圧力（すなわちＰｄｅｔ）は、膀胱内のカテーテルからの圧力測定値と、患者の膣又は直腸内に載置された腹圧を測定する追加の圧力センサからの圧力測定値とを差分することによって推定された。ＵＤＳの充填期において、各増分量に複数の測定値が得られた。ＵＢＶ測定の取得時、ＵＢＶ及びＵＤＳ測定に関するタイムスタンプが記録され、ＵＢＶ測定に関連する及びＰｄｅｔ測定がマークされた。
【0098】
膀胱内圧測定中、各充填体積で複数のＵＢＶ測定値（２～３）が収集された。タイムスタンプされたＵＢＶ／ＵＤＳ測定間の時間間隔の各実験に関して、平均及び標準偏差のヒストグラムが計算された。同じ充填体積で行われた測定及び充填体積間で行われた測定に関して個別の平均及び標準偏差が収集された。
【0099】
ＵＢＶ測定値は、プログラマブル超音波システム（Ｖｅｒａｓｏｎｉｃｓ、Ｒｅｄｍｏｎｄ、ＷＡ）及び２．５ＭＨｚの中心周波数を有する曲線アレイトランスデューサ（Ｃ４－２、ＡＴＬ／Ｐｈｉｌｉｐｓ、Ｂｏｔｈｅｗｅｌｌ、ＷＡ）を用いて収集された。音響放射力のメカニズムによって組織に波を誘導するために６００～９００マイクロ秒のトーンバースト音響パルスが膀胱壁に印加され、そしてコヒーレントな配合のために三つのステアリング角度を有する超高速平面波撮像シーケンス（毎秒２５００フレーム）が、結果として生じる膀胱壁に沿ったラム波の伝搬を捕捉した。その結果生じる同位相及び直交（ＩＱ）データから、自己相関技術を用いて粒子運動速度が推定された。その後、膀胱壁に沿って伝搬するラム波の群速度（ＧＶ）が推定された。弾性は、バルク波の場合により厳密であるように、波速度の二乗に相関することが予想されるため、分析のために群速度が二乗された（ＧＶ２）。
【0100】
データ準備及びピーク識別
【0101】
ＤＯに関連するＧＶ２及びＰｄｅｔ測定系列の変動を特徴付けるために、そのような変動は、膀胱充填に関連する任意のより大きな全体傾向から分離された。たとえば大きな外れ測定値のように、膀胱壁活動の結果である可能性が低く、特徴を左右し得る変動は除去された。後者を考慮するために、統計的フィルタリングアプローチ（たとえばハンペルフィルタリング）が用いられた。系列の主な傾向からＤＯ関連の変動を分離するために、系列は、それぞれＬＥ系列とＥＴＰ系列とに分解された。
【0102】
ピーク識別の前に、膀胱の漏れが報告されたＧＶ２及びＰｄｅｔの値は、所与のデータ系列から除去された。この例では、排尿の開始に伴う正常な排尿筋収縮が含まれることを避けるために、実験の最終量での測定値もピーク識別から除外された。そのような外れ値は、たとえば、測定中の患者の動きや異常に弱いＵＢＶ信号に起因し得る。これらの外れ値を除去するために、ハンペルフィルタと呼ばれる統計的外れ値フィルタがＧＶ２及びＰｄｅｔデータ系列に適用された。このフィルタは、外れ値であると決定された点をサンプルウィンドウの中央値に置き換えることによって機能する。フィルタは、外れ値を識別するために以下を用いる。中央値ｍｉを有するスライディングウィンドウの中心点であるデータ系列ｘｉからの所与のサンプルは、
の場合、外れ値と見なされる。
【数34】
【0103】
ここで、ｎσは、標準偏差の閾値数であり、σｉは、サンプリングウィンドウの推定標準偏差である。サンプリングウィンドウの標準偏差は、
として中央絶対偏差から推定され得る。
【数35】
【0104】
ここで、ｗは、ウィンドウの中心点の両側にある画素の数であり、パラメータｋは、中央絶対偏差を標準偏差にスケーリングする。標準正規分布の場合、
である。ＵＢＶデータ分析のために、以下の理由から、ｎσは６に設定され、ｗは５に設定された。フィルタは、不均一にサンプリングされたデータに関して固定された時間ウィンドウを許可しないが、外れ値が検出された場合を除き、信号を修正することはない。同様に、ウィンドウは、信号の長さに対して大きく設定された（中央値のサンプリング間隔及び中央値の持続時間に基づく～１６点に対し、ｗ＝５のウィンドウ内に１１点）。識別する可能性のある異常な外れ値の数に関して外れ値閾値を控えめに設定することにより、フィルタは、系列内の極端な外れ値点のみを変更し、それ以外の残りの点を変更せずに留めるように構成され得る。
【数36】
【0105】
フィルタによって修正された点が少なく、サンプルウィンドウが大きいことにより、固定されていない時間ウィンドウの影響は最小限であることが予想された。ｎσを６に設定することは、ガウス分布の場合、データの９９．９９９７％が変更されずに残ることを意味する。ガウス分布を有するフィルタウィンドウは、データ系列内の異常な外れ値を識別するための信頼性が高く客観的な発見的教授法を提供し、人間の判断に依存しないと考えられる。
【0106】
外れ値フィルタリングの後、結果として生じるデータ系列は、ＥＴＰ系列とＬＥ系列とに分解された。膀胱充填に関連するデータ系列内の任意の傾向を捕捉するために、生データ系列（ｘ）に平滑化フィルタが適用された。平滑化は、２次Ｓａｖｉｔｚｋｙ－Ｇｏｌａｙフィルタを介して系列に適用された。Ｓａｖｉｔｚｋｙ－Ｇｏｌａｙフィルタは、スライディングウィンドウ内のデータに（２次）多項式回帰を適用する。時間間隔は不規則な間隔を有したため、離散ウィンドウサイズではなく平滑化係数が用いられた。平滑化係数は、平滑化によって減衰する信号エネルギの割合に対応する０と１との間の値をとる（たとえば、０．５の平滑化係数は、信号エネルギの５０％を減衰する）。平滑化係数は、データの主な傾向を除いてデータ系列内の全ての変動をフィルタするために、０．９５に設定された。
【0107】
平滑化の後、ｘは、フィルタリングされたデータ系列（ｘｓ）と差分され、差分系列の負の値に関連する時点（すなわち、（ｘ－ｘｓ）＜０である時点）が識別された。時間境界における増加が過渡的ピークとして扱われないように、開始時点及び終了時点も識別された時点に含まれた。
【0108】
その後、識別された時点におけるｘからの値が収集された。識別されなかった時点に値を割り当てるために、線形補間が用いられた。その結果生じる系列は、ｘと同じ要素数を有した。その結果生じる系列を、ＬＥ系列（ｘＬＥ）と称する。
【0109】
その後、ｘからｘＬＥを差分することによってＥＴＰ系列（ｘＥＴＰ）が決定された。ＥＴＰ系列及びＬＥ系列は、更なる分析のために用いられ得る。
【0110】
診断メトリック
【0111】
前項で識別されたピークに対するＤＯ膀胱と非ＤＯ膀胱との比較を行うために、識別された過渡的ピークを特徴付ける単一のメトリックが用いられた。また、そのようなメトリックは、原理的に、上述したＵＢＶ測定の可変の実験期間及び範囲を考慮すると、実験ごとに比較可能でなくてはならない。ＤＯ指数と呼ばれるこのメトリックは、Ｉで示される。
【数37】
【0112】
ここで、ｘＥＴＰ（ｔ）は、ＥＴＰデータ系列の連続信号を示し、ｘＬＥ（ｔ）は、ＬＥデータ系列の連続信号を示し、各々は、所与の実験に関する時間の関数である。基本的に、Ｉは、過渡的ピーク及び下側包絡線からの（それぞれＥＥＴＰ及びＥＬＥと示される）信号エネルギの比の平方根である。よって、ＤＯを有する膀胱からのデータ系列は、ｘＬＥ（ｔ）に比べてｘＥＴＰ（ｔ）における信号エネルギが大きいので、ＤＯを有さない膀胱よりもＩが大きいことが予想される。Ｉは単位がないため、ＲＭＳ（ｘＥＴＰ）をＲＭＳ（ｘＬＥ）で正規化することによって、より長い実験期間、及び測定されるＧＶ２又はＰｄｅｔにおけるより広い範囲が考慮されるので、実験間で比較を行うことが合理的である。
【0113】
ＲＯＣ分析によると、この分析は、ＤＯに関して０．７０の感度及び０．７５の特異度をもたらすことが示された。これは、ＵＤＳから記録された同時Ｐｄｅｔ測定値に同じ分析を適用することによって生成された感度及び特異度（それぞれ０．７０及び０．８３）に匹敵する。
【0114】
ＵＢＶ測定は、充填体積を増加してＵＤＳ測定と同時に行われた。このアプローチの利点は、ＵＢＶ及び同時ＵＤＳ測定値が、泌尿器科医のＤＯ診断に用いられるＵＤＳチャート（グランドトゥルース）と同時に収集された点である。このアプローチは、臨床検査とＵＢＶ研究とが異なる時間に行われた場合に生じるばらつきを除去する。このアプローチの別の利点は、同時に取得されたＧＶ２及びＰｄｅｔの両方に対するピーク分析を比較することが可能な点である。ＵＤＳデータからのＰｄｅｔ及びＵＢＶデータからのＧＶ２にこの方法を適用することにより、同様の診断性能が提供される。ＵＢＶ及びＵＤＳは、提案された診断指標に対して同様の情報を提供した。したがってＵＢＶは、ＤＯ診断のためのＵＤＳに代替する非侵襲的方法を提供する。
【0115】
図１０Ａ～Ｂを参照すると、膀胱活動を測定するための上述したＵＢＶ法の検証が示される。ＤＯ（図１０Ａ）及び非ＤＯ（図１０Ｂ）に関するＧＶ２及びＰｄｅｔデータ系列の非限定的な例である。図１０Ａ～Ｂは、Ｐｄｅｔ及びＧＶ２の両データ系列の二軸プロットを示す。図１０Ａにおける矢印は、ＤＯに関連するピークを示す。
【0116】
これらのピークは、系列をそれぞれ過渡的ピーク（ＥＴＰ）からのエネルギと下側包絡線（ＬＥ）信号からのエネルギとに分解した場合、データ系列の大きな傾向から分離可能であり得る。非ＤＯ膀胱において過渡的な変動が観察可能であったが、これらの変動に関連する信号エネルギは比較的低かった。このパターンは、信号ＥＴＰ及びＬＥそれぞれの比の平方根であるＤＯ指数Ｉで捕捉された。特に、この例では、ＤＯを有さない膀胱に対し、ＤＯを有する膀胱の例において、Ｐｄｅｔ及びＧＶ２の両方に関してより高い値のＩが観察された。
【0117】
この特徴付けアプローチの統計及び診断分析は、非パラメトリック仮説試験（Ｗｉｌｃｏｘｏｎ試験）及び受信者動作特性曲線（ＲＯＣ）分析を用いて適用された。その分析結果は、表１にまとめられる。例となる膀胱と一致して、ＧＶ２及びＰｄｅｔに関する中央値ＤＯ指数は、ＤＯを有さない膀胱（０．２５及び０．３１）よりもＤＯを有する膀胱（それぞれ０．４３及び０．７０）で大きかった。Ｗｉｌｃｏｘｏｎ試験により、ＧＶ２及びＰｄｅｔの両方における統計的有意な確率差（ｐ＜０．０１）が示され、統計的有意なＲＯＣ曲線下面積（ＡＵＣ）が示唆される。
【0118】
図１１Ａ～Ｄを参照すると、ＥＴＰ系列及びＬＥ系列への生系列の分解を示すために、生系列に下側包絡線（ＬＥ）が重ね合わせられたＧＶ２値を示す、図１０Ａ～Ｂの非限定的な例が示される。分離されたＥＴＰ系列及びＬＥ系列は、信号指数（Ｉ）を計算するために用いられた。データ系列に関して識別されたピークのアンサンブルは、生系列と下側包絡線との差である。
【0119】
図１２Ａ～Ｄを参照すると、ＥＴＰ系列及びＬＥ系列への生系列の分解を示すために、生系列に下側包絡線（ＬＥ）が重ね合わせられたＰｄｅｔ値を示す、図１０Ａ～Ｂ及び図１１Ａ～Ｄの非限定的な例が示される。分離されたＥＴＰ系列及びＬＥ系列は、信号指数（Ｉ）を計算するために用いられた。データ系列に関して識別されたピークのアンサンブルは、生系列と下側包絡線との差である。
【0120】
図１３Ａ～Ｃを参照すると、ＧＶ２及びＰｄｅｔそれぞれについて計算されたＩの関連ＲＯＣ曲線とともに、統計及び診断分析のグラフ結果の非限定的な例が示される。図１３Ａは、Ｐｄｅｔ信号指数に対するＧＶ２信号指数の散布図である。破線は、それぞれＧＶ２及びＰｄｅｔのカットオフを示す。図１３Ｂ及び図１３Ｃにおいて、ＧＶ２（図１３Ｂ）及びＰｄｅｔ（図１３Ｃ）に適用されたＤＯの分類子として信号指数（Ｉ）を用いて生成されたＲＯＣ曲線が示される。円形マーカは、選択された最適カットオフを示す。最適カットオフは、基準最近接（０，１）（すなわち、ＲＯＣ曲線の左上角部に最も近い曲線上の点）を用いて推定された。対応するカットオフ値は、ＧＶ２について計算されたＩに関して０．３２７、Ｐｄｅｔから計算されたＩに関して０．５６７であった。これらに関連する感度及び特異度は、表１に報告される。図１３Ａは、７６人全ての患者によるＵＢＶ／ＵＤＳ実験に関して計算されたＤＯ指数の散布図によってこれらの統計結果をグラフで示し、ＧＶ２及びＰｄｅｔのＲＯＣ曲線（図１３Ｂ、図１３Ｃ）を伴う。ピーク特徴付けの分類特性は、ＧＶ２とＰｄｅｔとの間で同等（すなわち、同等の感度、特異度、及びＡＵＣ）であった。
【表1】
【0121】
図１４Ａ～Ｃを参照すると、非限定的な膀胱実験の各々に関するサンプリング間隔統計値（図１４Ａ及び図１４Ｂ）及び実験期間（図１４Ｃ）のヒストグラムの非限定的な例が示される。各実験のサンプリング間隔の標準偏差（ＳＤ）（図１４Ａ）及び平均（図１４Ｂ）は、同じ充填体積の間の間隔及び充填体積間の間隔について計算された。テキストボックスで報告された統計は、四分位間（ＩＱ）範囲を有するそれぞれのヒストグラムに関する。
【0122】
図１５を参照すると、ハンペルフィルタリング前及びハンペルフィルタリング後のデータ系列の非限定的な例が示される。
【0123】
図１６Ａを参照すると、ＥＴＰデータ系列を生成するための非限定的な方法例１６００のフローチャートが示される。ステップ１６０２において、一般に（ｘ，ｔ）のフォーマットであってよい生データ系列がアクセス又は取得され得る。データは、データアーカイブシステム又はデータベースからアクセスされてよく、あるいは、たとえば図１～４に示すようなシステムを用いて取得され得る。ステップ１６０４において、データ系列が平滑化され得る。非限定的な例において、生データ系列は、２次Ｓａｖｉｔｚｋｙ－Ｇｏｌａｙフィルタを用いて平滑化され得る。ステップ１６０６において、平滑化されたデータ系列は、生データ系列から減算され得る。ステップ１６０８において、負の値を有する時点、又は結果として生じる系列の開始又は終了の時点が識別され得る。ステップ１６１０において、識別された時点における生データ系列“ｘ”からの対応する値が収集され得る。ステップ１６１２において、除外された時点に生データ系列が補間され得る。ステップ１６１４において、上述したように、下側包絡線（ＬＥ）データ系列が生成され得る。ステップ１６１６において、ＬＥデータ系列が生データ系列から減算されてよく、ステップ１６１８において、ＥＴＰデータ系列が生成され得る。
【0124】
図１６Ｂを参照すると、生データ系列及び処理後のグラフの非限定的な例が示される。生データ系列は、２次Ｓａｖｉｔｚｋｙ－Ｇｏｌａｙフィルタを用いて平滑化されて示される。生データは、平滑化データによって差分され、そして、値が平滑化曲線を下回る（すなわち負の値を有する）時点が識別された。識別された時点は、点の補間によって生データ系列からＬＥ系列を画定するために用いられた。ＬＥ系列と生データ系列との差により、ＥＴＰ系列が決定された。
【0125】
組織体積の特性を表すパラメータを評価するための方法及びシステムは、有形非一時的メモリに格納された命令によって制御されるデータ処理電子回路（プロセッサ）の使用を必要とし得ることが理解される。メモリは、ランダムアクセスメモリ（ＲＡＭ）、読取専用メモリ（ＲＯＭ）、フラッシュメモリ、又は制御ソフトウェア又は他の命令及びデータを格納するために適した他の任意のメモリやそれらの組合せであってよい。当業者は、本開示において説明された方法の機能（たとえば実験データの収集及び／又は格納、求められるパラメータを定義するためのデータ処理など）を定義する命令又はプログラムが、非書入れ可能な記憶媒体（たとえばＲＯＭなどのコンピュータ内の読取専用メモリデバイス、又はＣＤ－ＲＯＭやＤＶＤディスクなどのコンピュータＩ／Ｏアタッチメントによって読取り可能なデバイス）に永久的に格納された情報、書入れ可能な媒体（たとえばフロッピーディスク、取外し可能フラッシュメモリ、及びハードドライブ）に変更可能に格納された情報、又は有線又は無線コンピュータネットワークを含む通信媒体を介してコンピュータに伝達される情報を含むがこれらに限定されない多数の形式でプロセッサに供給され得ることも容易に認識すべきである。また、本開示で説明される方法は、ソフトウェアにおいて具体化され得るが、方法を実施するために用いられる機能は、任意選択的又は代替的に、一部又は全体が、たとえば組合せ論理、特定用途向け集積回路（ＡＳＩＣ）、フィールドプログラマブルゲートアレイ（ＦＰＧＡ）、又は他のハードウェアや、ハードウェア、ソフトウェア、及び／又はファームウェア構成要素の何らかの組み合わせなどのファームウェア及び／又はハードウェア構成要素において具体化されてもよい。したがって、本開示において説明されるステップを実行するために特別にプログラムされたプロセッサを用いるシステム及び方法、及びそのようなステップを実行するようにプロセッサを統治する命令のセットを含むコンピュータプログラム製品は、本開示の範囲内である。
【0126】
本発明は、一つ又は複数の好適な実施形態に関して説明されており、明確に記載された実施形態以外にも多数の等価物、代替物、変形物、及び修正物が可能であり、本発明の範囲内であることを理解すべきである。たとえば、本開示で説明される方法は、超音波エコーデータのスペクトル分析に基づくが、代替として分散データを決定するために時間領域分析が用いられ得ることが理解される。一例において、群速度又は波ピークの速度が膀胱筋圧に相関付けされ得る。同様に、軟組織における機械波の複雑な伝搬は、（物体粒子の変位が波伝搬方向に垂直であり、定義上、そのような伝搬に影響を及ぼす媒質内の境界が存在しないという波の考えに基づく）せん断波モデル、又は（波がプレートに沿って、又は表面と平行に伝搬し、粒子がプレート表面に垂直に、したがって波伝搬方向に対し垂直に動く場合に定義される）ラム波に基づく波モデルによって近似され得る。ラム波ベースのモデルにおいて、（有限）プレートの２つの表面は、波を前方に導く境界である。膀胱壁を伝搬する波の場合、本開示で説明されるいくつかの実施形態において、ラム波モデルは、実用的かつ良好な近似を提供する。ただし、（たとえば上述したせん断波モデルなどの）他の波モデルも、膀胱での波動を近似するために一般的に用いられ得る。
【図1】
【図2】
【図3】
【図4】
【図5A】
【図5B】
【図5C】
【図5D】
【図6A-B】
【図7A-B】
【図8】
【図9】
【図10A-B】
【図11A-D】
【図12A-D】
【図13A-C】
【図14A-C】
【図15】
【図16A】
【図16B】
【国際調査報告】