JP5740547B2 - 頭蓋内の弾性を測定する方法 - Google Patents

頭蓋内の弾性を測定する方法 Download PDF

Info

Publication number
JP5740547B2
JP5740547B2 JP2011521326A JP2011521326A JP5740547B2 JP 5740547 B2 JP5740547 B2 JP 5740547B2 JP 2011521326 A JP2011521326 A JP 2011521326A JP 2011521326 A JP2011521326 A JP 2011521326A JP 5740547 B2 JP5740547 B2 JP 5740547B2
Authority
JP
Japan
Prior art keywords
ultrasound
icp
pulse
ultrasonic
signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP2011521326A
Other languages
English (en)
Other versions
JP2011529744A (ja
Inventor
ヴィテルボ,ヴィター,ダニエル デ
ヴィテルボ,ヴィター,ダニエル デ
Original Assignee
ラパポート,アーサー
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by ラパポート,アーサー filed Critical ラパポート,アーサー
Publication of JP2011529744A publication Critical patent/JP2011529744A/ja
Application granted granted Critical
Publication of JP5740547B2 publication Critical patent/JP5740547B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0808Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the brain
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • A61B8/14Echo-tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/03Detecting, measuring or recording fluid pressure within the body other than blood pressure, e.g. cerebral pressure; Measuring pressure in body tissues or organs
    • A61B5/031Intracranial pressure
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/40Detecting, measuring or recording for evaluating the nervous system
    • A61B5/4058Detecting, measuring or recording for evaluating the nervous system for evaluating the central nervous system
    • A61B5/4064Evaluating the brain
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/485Diagnostic techniques involving measuring strain or elastic properties
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5223Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Surgery (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Neurology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Neurosurgery (AREA)
  • Physiology (AREA)
  • Hematology (AREA)
  • Psychology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Description

本発明は、頭蓋内圧を非侵襲的に測定する方法に関する。
一般的に、ヒトのような哺乳類は、一定の頭蓋内血液量を有し、それ故、一定の頭蓋内圧(「ICP」)を有する。しかし、種々の正常および病的な症状は、頭蓋内圧に変化を生じさせるかもしれない。頭蓋内圧の上昇は、脳への血流を減少させるかもしれず、脳が自動的に圧縮され、ついにはヘルニアになるかもしれない。頭蓋内圧の上昇のうち最も多い原因は、頭部外傷である。頭蓋内圧の上昇の次に多い原因は、限定されないが、揺さぶられっ子症候群、硬膜外血腫、硬膜下血腫、脳出血、髄膜炎、脳炎、鉛中毒、ライ症候群、ビタミンA過剰症、糖尿病性ケトアシドーシス、水中毒、脳腫瘍、頭蓋腔における他の腫瘤または血栓、脳膿瘍、脳卒中、ADEM(「急性散在性脳脊髄炎」)代謝疾患、水頭症、並びに硬膜静脈洞および静脈血栓症を含む。頭蓋内圧の変化は、常にモニターし、外科的介入ができることが必要とされるため、医療において、頭蓋内圧をモニターする技術の開発が、重要な目標として残っている。米国特許第6,875,176号。
従来の頭蓋内圧モニター装置は、硬膜外カテーテル、クモ膜下用ボルト/ねじ、脳室造瘻術カテーテル、および光ファイバーカテーテルを含む。これらの方法およびシステムの全ては侵襲的であり、かつ、高度に訓練された神経外科医による侵襲的な外科処置が必要とされる。また、これらの技術には、頭蓋内圧を迅速的または定期的にモニターすることを満足させる技術が無い。さらに、これら全ての従来技術はICPを局所的に測定し、局所的なICPが脳全体のICPを反映するという推測が成されるものである。米国特許第6,875,176号の技術では、既存する方法のこうした限界を示している。
広く受け入れられる、ICPを非侵襲的に測定する方法がない。しかし、臨床的には、ICPを測定する有効な手段の開発は、ICPにより臨床転帰が予測でき、より有効な治療へと変わるかもしれないため、非常に重要である。例えば、外傷性脳損傷後では、頭蓋内圧は上昇傾向があり、速やかに認識されかつ治療する必要がある。Zanierら、Critical Care 11:R7(「2007」)。ICPを測定する既存の標準器は、硬膜外変換器、または脳室内もしくは脳実質内カテーテルの配置を含む直接的で侵襲的な測定が必要とされる。Frankら、Zentralbl Neurochir 61 (「4」):177−80(「2000」)。侵襲的な方法の使用では、感染、出血、外科手術のミスによる損傷度が高まる。Czosnykaら、J.Neurol.Neurosurg.Psychiatry 75:813−821(「2004」)。
耳音響放射(「Frankら、Zentralbl Neurochir 61(「4」):177−80(「2000」)」)を測定すること、および経頭蓋ドップラーを用いて超音波(Ragauskasら、Innovative non−invasive method for absolute intracranial pressure measurement[オンライン]、[2008年7月30日に検索]、インターネットから検索される<URL: http://www.neurosonology.org/bern2002/abs_12.html>)を測定することを含む、ICPを非侵襲的に測定するための種々の異なる技術が探究されてきた。
例えば、米国特許第6,702,743号(「‘743の特許」)は、ICPを測定する非侵襲的な手段を開示している。超音波プローブは、患者の頭部に配置され、次いで患者の頭蓋骨および脳を通って伝播する超音波を発生させるために使用される。

超音波パルスは、超音波プローブと垂直する経路にある頭蓋骨および軟組織から反射される。次いで、発生したエコーEG信号の一部分が選択され、エコーEG信号はエコーパルスグラフ(「EPG」)信号を発生させるために選択した部分にわたって積分される。しかし、‘743の特許の方法を使用してICPを決定するためには、操作者がEPGの一部分を手動で選択するか、または「門付け」(gate)し、次いで、どれが脳の目的部位について最適なEPG波形を提供しているかを決定するため、各門においてEPG波形を見直さなければならない。
出願人は、専門家によるEPG波の手動による見直しを必要とせず、ICP、より一般的には、脳の弾性を非侵襲的に測定する新規の方法を開発した。ICPは、複数の超音波パルスから発生した一式の相互作用した超音波信号(「IUS」)の決定に基づきICPを算出する、模擬人工ニューラルネットワーク(SANN」)に連結されたアルゴリズムを使用して決定される。本発明の方法およびシステムは、EPG波を手動で見直しすることなく、ICPを迅速に決定することができる。
本発明の目的は、エコーパルスグラフ(「EPG」)の形で表される、反射された超音波信号の解析に基づきICPを非侵襲的に測定する技術を提供することである。
ICPは、まず、患者の頭蓋内に少なくとも約1MHzの超音波パルスを送信することにより測定される。次いで、この超音波パルスは、第3脳室の壁を含む頭蓋のさまざまな構造により反射される。反射信号は変換器により受信され、一式の情報が生成される。
本発明は、患者から複数の超音波信号を取得する。第3脳室の壁の状態は、脳に出入り(「収縮および拡張」)する血流に起因して常に変化しているので、コンピュータは、それぞれの波形にある偏差に基づき、第3脳室の領域を突き止めるために各信号を比較する。
第3脳室が突き止められると、第3脳室内の波の部分に沿ったデータポイントは、ICPを算出するために使用される。ICP値は、ICP値が既知の患者から得たICPデータと抽出値とを相関させるアルゴリズムにより算出される。システムが既知のICP値との参照により比較または学習されると、コンピュータにより算出は自動的に完了する。
超音波を送受信し、人工ニューラルネットワークを学習させるのに好ましい装置のブロック図を示す。 頭蓋内圧モニターシステムのフローチャートである。 1つの超音波反射信号(「USRS」)全体を示す。 USRSの部分と脳の部分との間の相関を示す。 ICP値を取得するために、ニューラルネットワークおよびアルゴリズムにUSRSのデータポイントを入力するフローチャートを示す。 入力として使用された特定のEPGポイントの例を示す。 QRS群および超音波信号の隣り合った比較図を示す。 学習処理で使用される逆伝播法のフローチャートを示す。 学習処理がICP値の測定可能範囲を作り出す方法の一例を示す。 学習処理のためのハードウェアの一実施形態を示す。
本発明は、身体内のICP、および、より一般的には、さまざまな器官または腔内またはこれらに最も隣接する組織の弾性を、非侵襲的に測定するための方法を目的とする。一実施形態では、ICPは、経頭蓋ドップラー信号を使用して頭蓋腔を超音波照射(isonate)することにより決定される。まず、第3脳室の前部壁および後部壁の位置が特定され、ICP波のグラフが作成される。次いで、ICPは、ニューラルネットワークを使用してICP波から算出される。より一般的には、本発明の方法およびシステムは、種々の異なる組織の組織弾性を測定するために使用され得る。
本発明一実施形態では、本発明の方法およびシステムは、超音波プローブを使用する。そのようなプローブは、例えば、チタン酸ジルコン酸鉛(「PZT」)、ポリ二フッ化ピリニデン(「PVDF」)、PZTセラミック/ポリマー合成物などでできた電極により活性化される1つ以上の圧電素子から構成されてよい。電極は電圧源に接続され、電圧波形は印加され、そして圧電素子は印加された電圧の周波数と対応する周波数で大きさが変化する。電圧波形が印加されると、圧電素子は、励起波形に含まれる周波数で連結された媒体内に超音波を放出する。逆に、超音波が圧電素子を打ちつけると、素子は電極間に対応する電圧を生成する。本発明は、当技術分野で周知されている数々の超音波プローブのうち任意のものを使用して実施されてよい。
図2により、本発明の方法の概要を提供する。超音波プローブ1は、頭蓋腔内に超音波を送信する。超音波プローブは、患者の頭部に配置され、次いで患者の頭蓋骨および脳を通って伝播する超音波パルスを発生させるために使用される。超音波パルスは、頭蓋9の後頭部分から反射されるとともに、脳組織8を横切る間に遭遇する半剛性または剛性の構造から反射される。横切る間に超音波パルスが遭遇するそのような構造の1つは、脳室の前部壁および後部壁を含む第3脳室である。超音波パルスは、反射されて超音波プローブ1に戻り、超音波獲得システム7に入力される。本発明の方法およびシステムとともに、市販で入手できる超音波装置を使用してよい(AdvancedTransducer Services,Inc、[オンライン]、[2008年7月30日に検索]を参照されたい。インターネットから検索される<URL: www.atsultrasound.com/>)。信号は、ロードアルゴリズムを備えたマイクロプロセッサ6により解釈され得、当該アルゴリズムが第3脳室壁の位置を特定し、EPGポイントをICP値に相関させる。
図1は、ICPを測定するために使用されるシステムの一実施形態を示す。システムはラップトップコンピュータのような任意の中央処理装置(「CPU」)またはマイクロプロセッサ6、ユニバーサル・シリアル・バス(「USB」)インターフェース5、デジタル信号プロセッサ(「DSP」)4、増幅器、アナログ−デジタル変換器(「ADC」)3、超音波回路2、および、送信器、受信器、および超音波パルスを発生させる超音波プローブ1を備える。システムは、心拍を測定する手段と統合される。図1で示される実施形態は、CPU6、アナログ−デジタル変換器3、および超音波プローブ1を有する、本発明の一例のシステムの見本形態のみを示すことを理解されたい。これらの部品全ては、一般的な電子機器供給業者から市販されている。
一般的な、市販部品が、本発明のシステムで使用されてよい。以下に記載する特定の部品は例示に過ぎず、本発明のシステムはこれらの部品に限定されない。例えば、DSP4は、Texas InstrumentsのC2000 DSCおよびTMS320C20x、Canberraの2060モデル、CEVA−Xl 641、CEVA−X1622、CEVA−X1620、またはCEVA−TeakLite−IIIでよい。DSP4は、プローブ1を介し、少なくとも1MHzの周波数を用いた電気パルスまたは信号を発生させ、プローブ1を通る反射波またはエコーを検出し、検出されたデジタル信号を処理する責任がある。レンジは、調査する信号に準じてDSP4のファームウェアにおいて変更することができる。
測定周期は、コンピュータ6から開始信号がDSP4に受信されると開始される。
これに応えてDSP4は、一連の超音波パルスを発生させるようにプローブ1に指示する。市販で入手できる超音波プローブが、本発明の方法およびシステムとともに使用されてよい(Advanced Transducer Services,Inc、[オンライン]、[2008年7月30日に検索]を参照されたい。インターネットから検索される<URL: www.atsultrasound.com/>)。超音波プローブ1は、少なくとも約1MHz、最大で約10MHzの周波数の超音波を送信できるべきである。
超音波源および検出器は、送信モード、または種々の反射モードもしくは散乱モードで用いられ得、圧力波がせん断波に移転することを試験するモードおよびその逆のモードを含む。超音波検出技術は、超音波照射された組織からの音響放出(「s」)をモニターするために使用されてもよい。検出技術は、後方散乱のような音響散乱、または音響放出の変化の測定に関与する。組織特性に関連する音響散乱または音響放出のデータの例は、音響信号の振幅の変化、音響信号のフェーズの変化、音響信号の周波数の変化、問い合わせ信号に対して散乱または射出された信号の長さの変化、心周期および/または呼吸周期の音響信号の初期および/または他の時期の最大振幅および/または最小振幅、並びに心周期の最大振幅および/または最小振幅とそれに続く振幅の平均値または分散もしくは分布の比率の変化、同じ位置における異なる時間および/または異なる位置における同じ時間での散乱または射出された信号の時間的分散または空間的分散の変化、移転の速度または加速度のような内因性による脳組織の移転または弛緩の変化に関する可能性のある全ての比率を含む。複数の音響問い合わせ信号は、同じまたは異なる周波数、パルス長、パルス繰り返し周波数、および強度において使用されてよく、複数の問い合わせ信号は、同時的および/または連続的に、同じ位置または複数位置から送られてよい。単一または複数の問い合わせ信号からの散乱または放出は、単一または複数の周波数、単一または複数の時間、および単一または複数の位置において検出され得る。
図3は、単一の超音波による反射信号(「USRS」)を示す。図表では、この超音波信号は、エコーパルスグラフすなわちEPG10として参照される。信号は、超音波パルスの経路における前部頭蓋および後部頭蓋の円蓋、および頭蓋内容物の解剖学的位置を指す、相互作用する信号である。第3脳室を含む脳を超音波照射する超音波信号は、特定の周波数特性を有する。リターン信号が変わっていない場合、EPGは解剖学的構造を測定し、同じ波形を反映して戻すだけである。しかし、超音波照射された超音波信号がその経路において全てのものと相互作用する場合、詳細には第3脳室の運動性だが、結果生じる波形すなわちEPGは相互作用的であり、一式の反射信号を取得するためにフィルターにかけられ、ICPを算出することができる。例えば、図4は、分類された相互作用的EPGである。波形の識別可能部分は、プローブ11内の反射信号(「」)、前部頭蓋の円蓋、硬膜、髄膜12の反射信号(「」)、脳13の反射信号(「」)、第3脳室14の反射信号(「」)、並びに硬膜および後部頭蓋の円蓋15の反射信号(「」)に一致する。
任意の心周期(「収縮期および拡張期」)の間、複数のEPG測定を行うことができ、図7は、心周期と脳との間の関連性を示す、EPGおよびQRS群の隣り合った比較図である。第3脳室の壁は、心周期(「収縮期および拡張期」)の間に拡張し、収縮する。それ故、心周期の間、超音波プローブに対する第3脳室の壁の位置は、変化する。
本発明の一実施形態では、少なくとも10のEPGの測定がされる。別の実施形態では、少なくとも25のEPGの測定がされる。第3実施形態では、少なくとも50のEPGの測定がされる。第4実施形態では、少なくとも100のEPGの測定がされる。EPG信号は、強度および時間の関数として、各々デジタル化され表示スクリーン上に表示される。図5に示されるように、全EPGのうち第3脳室領域から作り出された16から19のポイントは、ICP値を算出するためにアルゴリズムに入力される。これらのポイントは、図6により明白に表わされており、これはどのようにEPGの第3脳室領域が経時的(「t」)に21から35の島ごとのポイントに分割されるのかを示す。
これらのポイントは、EPGの単離部分からデジタル化されたデータポイントの個別の束を表しており、これは次いで方程式に基づきICPを算出するために使用される。
ICP=Σtanh(ΣIxW+b)W+b
であり、ここで、Iはエコーパルスグラフ21から35のうち選択した部分の全データポイントからなる入力行列を表し、Wは学習処理を通して取得される重み行列であり、bはコンピュータ6により割り当てられたランダム・バイアス定数である。
入力行列は(「n×k」)の数学行列であり、ここでn行はサンプルの数に等しく、本発明の一実施形態では、この値は少なくとも10である。k列は脳室壁間に認められた各EPGが伴うデータポイントに等しい。行列は、公知の数学的手段により算出される。
W値または重み行列は、学習処理または相関処理を通して取得され、これは一回で行われなければならない。SANNを学習させる方法は、V.D.De ViterboおよびJ.C.Belchior、Artificial Neural Networks Applied for Studying Metallic Complexes,Journal of Computational Chemistry,vol.22,no.14,1691−1701(「2001」)に説明されている。学習処理は、ネットワークに入力および所望の出力を繰り返し提示することから構成される逆伝播アルゴリズムである。重みは、ネットワークにおいて所望の誤差が達成されるまで段階的に補正される。この方法は図8に示される。本発明の一実施形態では、逆伝播方法は次の式により実行され、
Figure 0005740547
ここで△W jiはl番目の層のj番目の要素とその前の層のi番目の要素との間の重みに対する補正を表す。数値out l−1はl−1番目の層についての出力結果を含む。パラメータηおよびμはそれぞれ、割り当てられた学習率およびモメンタム定数である。これらは、学習処理の収束率を決定する。通常、これらのパラメータは、最良の収束率を取得するように、動的に調整される。学習状態の間に導入される誤差は、
Figure 0005740547
および、
Figure 0005740547
として算出され、ここでyはl番目の層の出力結果out と比較される出力ターゲットである。次いでネットワーク誤差は、
Figure 0005740547
として算出することができる。学習処理のために、ニューロンの振舞いは、中間層についてはシグモイド関数により、出力層においては線形関数により算出された。
関数を最小化するために、本発明の一実施形態は、Levenbergにより提案され、Marquardtにより実行されたロバスト法(Marquardtら、J Soc Ind Appl Math 11:431(「1963」)を使用する。ロバスト法は、最急降下法およびニュートン法の動的調整を通して機能する。これは、最小値を探索する点で一層迅速であることに有益性がある。レーベンバーグ・マルカート法(LMM)によれば、重みの更新行列は、
Figure 0005740547
として算出することができ、ここでHはヘッセ行列であり、βは変数パラメータであり、通常はβ=0.01として開始する。後者は最小値検索の間に局所的な誤差の推定に従って変更される。Iは恒等行列である。LMMが使用されるときに最も困難な課題は、Hの算出に起因し得、これは
Figure 0005740547
により近似され、ここでJはヤコビ行列であり、
Figure 0005740547
により得られ、ここでlは全重みの相対誤差(方程式(4))である。ヘッセ行列を解くためのこの近似式は、二次導関数の計算を回避し、これにより算出を単純にする。前述の近似を方程式(5)に代用すると、
Figure 0005740547
が取得され、方程式(8)は、βが0に近づく場合は純粋にガウス・ニュートン法に、またはβが無限に近づくときは最急降下法に近似する。
本発明に従うと、この手段では、最初に、ICP値はランダムに割り当てられたW値を用いた方程式により算出される。結果生じるICP値はテスト値である。参照ICP値は、ICPを測定する公知の侵襲的手段により測定される。次いで学習は、公知の侵襲的方法から取得された参照ICP値に対して、そのテストICP値を比較することに関与する。ICP値の差が許容できる誤差よりも大きい場合、ランダムなW値は調整される。W値を調整すると、方程式を使用して新しいテストICP値が算出され、この値は参照ICP値と再び比較される。重みの値を調整し、新しいICP値を算出し、参照ポイントに対して当該ICP値を比較するというこの学習処理は、この処理から算出されたICP値が、参照値に対して許容できる誤差の範囲に収まるまで繰り返される。これが生じると、W値はコンピュータ6に記憶され、テストICP値として取得された、その特定のICP値と自動的に相関される。本発明の一実施形態では、ニューラルネットワークを学習させるアルゴリズムは以下の通りである。
●BEGIN
○WHILE START=ON
■ADCからデジタル化されたエコーのサンプルを獲得する
■サンプルをファイルに保存する
■サンプルのグラフを描く
■正しい波を選択する(手動処理)
■波が正しい場合
●START=OFF(手動)
■END IF
○END WHILE
○ニューラルネットワークの入力数=306
○隠れたニューロン数=2
○Wl(2X306)=ランダム値
○W2(1X2)=ランダム値
■WHILE(誤差>0.001)
●非侵襲的ICP値=
W2*(TANH(W1*デジタル化されたエコー)))
●誤差=侵襲的ICP値−非侵襲的ICP値
●レーベンバーグ・マルカート法を使用して
WlおよびW2を算出する
●Wl=Wl+DWl
●W2=W2+DW2
■END WHILE
●END BEGIN
コンピュータが重みの値および対応するICP値のインデックスまたはデータベースを作成するために、この学習処理は全ての可能なICP値について完了されなければならない。学習後、コンピュータ6は、侵襲的な処置無しに、各一式の入力およびICP値について適切なW値を一致させることにより自動的にICP値を算出することができる。図9は、学習処理によりICP値の測定可能範囲が広くなる様子を例示する。9人の患者のICP値を、3つの異なるICP値で3人ずつ3グループ分取得し、彼らの超音波データを初期データとして本発明に入力することにより、本発明に比較のベースラインが与えられる。また、本発明の動作範囲は、学習された既知のICP値の範囲に等しいだろう。
それ故、ニュートラルネットワークは、頭部の動的特性を相関させるアルゴリズム(「ACDPH」)20である。アルゴリズムは入力されたデータを使用してICP波を作り出す。次いで、EPG波に沿った時間(「t」)における各ポイントが、複数のEPG波にわたって描画される。理解されるように、最大「n」個までのサンプルを単一のEPG波から作ることができる。次いで、複数のEPG波について、各時間(「t」)におけるEPG波の振幅を示すグラフが、各時間(「t」)ごとに作成される。頭蓋の後頭部分のように心周期の間中変化しない構造については、時間(「t」)における複数のEPG波からの抽出サンプルを示すグラフは、直線である。第3脳室に沿ったポイントについてはそうではない。グラフ上で、これは、周期中のEPG波の振幅の変化に反映されている。より具体的には、この変化は正弦波パターンを持つICP波として表され、心周期を通して壁が拡張および収縮していることを反映する。学習中、ADCPH20は、入力されたポイントの上方および下方の境界を取得し、侵襲的な装置から取得された患者のICP値と、そのデータを相関させる。学習後、ADCPHは、侵襲的な方法を使用すること無く自動的に患者のICPを算出することができる。本発明の一実施形態では、ICP値を取得するアルゴリズムは以下の通りである。
●BEGIN
○WHILE START=ON
■学習済みニューラルネットワークWlおよびW2をロードする
■ADCからデジタル化されたエコーのサンプルを取得する
■サンプルをファイルに保存する
■サンプルのグラフを描く
■正しい波を選択する(手動処理)
■波が正しい場合(手動による工程)
●非侵襲的ICP値=
W2*(TANH(W1*デジタル化されたエコー))
■END IF
■END WHILE
●END BEGIN
図10は、即時学習処理のための、ハードウェアの一実施形態を示す。侵襲的にICPモニターするハードウェアから得たデータ、非侵襲的にICPをモニターするハードウェアから得たデータ、および心電図(EKG)のデータは、DSPに入力され、次いで当該DSPはUSBインターフェースを介してラップトップコンピュータに接続される。データ出力のために、本発明の一実施形態では、ラップトップコンピュータは、そのモニターにEPG、EKG、および算出された非侵襲的なICP値を表示する。
‘743の特許と対照的に、本発明は、第3脳室における経時的変化を考慮するため、より精密なICP測定を提供する。‘743の特許は、脳組織を通る血流が最初に脳を出る時点に依存する。さらに、エコーEG信号およびECGからEPGを発生させた後、先行技術の特許は、ICP値に一致するEPGの部分を選択することを操作者に依存する。本発明では、コンピュータプログラムがグラフの関連部分である第3脳室を特定する。最後に‘743の特許は、方程式ICP=ρ(「t/T」)*[t/T]−βに基づきICPを算出し、当該方程式は、p(「t/T」)を定義するために4つの異なる方程式に依存する。
本発の明範囲は、本明細書において先に具体的に示され、説明されたことに限定されない。特許およびさまざまな公表文献を含む数多くの参考文献が、本発明の説明において引用され、考察されている。そのような参考文献の引用および考察は、単に本発明の説明を明白にするものであり、どの参考文献もが本明細書で説明した本発明に対する先行技術であることを認めるものではない。本明細書で引用され、考察された全ての参考文献は、参照によりその内容全体が本明細書に組み込まれる。本明細書で説明したことについて、変化、修正、および他の実行が、本発明の精神および範囲から逸脱すること無く、当技術分野で通常の技術を有する当業者には考え出されるであろう。本発明の特定の実施形態を示し、説明してきたが、変化および修正は本発明の精神および範囲から逸脱すること無く、当技術分野で通常の技術を有する当業者には理解されるであろう。前述の説明および添付の図面に掲げた事柄は、例示のみとして提供され、限定されるものではない。本発明の実際の範囲は、以下の請求の範囲に定義されることを目的とする。

Claims (15)

  1. 哺乳類の身体内の組織の弾性を測定するための方法であって、
    a、記身体内の組織に、少なくとも1つの超音波パルスを送信するステップと、
    b、前記送信した超音波パルスの反射信号を検出するステップと、
    以下の数式を使用して前記組織の弾性を算出するステップとを含み、
    Σtanh(“ΣIxW+b”)W+b
    上記数式において、Iが前記反射信号に基づくデータポイントの入力行列であり、Wが調整用の重み行列であり、bが割り当てられたランダムバイアス定数である方法。
  2. 前記超音波パルスが、少なくとも1MHzの周波数を有する請求項1の方法。
  3. 前記超音波パルスが、少なくとも5MHzの周波数を有する請求項1の方法。
  4. 前記超音波パルスが、少なくとも10MHzの周波数を有する請求項1の方法。
  5. 少なくとも10の超音波パルスが前記身体内に送信される、請求項1の方法。
  6. 前記超音波パルスが、0から200ボルトの間の振幅を有する請求項1の方法。
  7. 前記超音波パルスが、14ビートの間、10メガ−サンプル/秒の速度で超音波受信器により受信される、請求項1の方法。
  8. 前記超音波パルスが、14ビートの間、1000メガ−サンプル/秒の速度で超音波受信器により受信される、請求項1の方法。
  9. 前記超音波パルスが、14ビートの間、1万メガ−サンプル/秒の速度で超音波受信器により受信される、請求項1の方法。
  10. 前記超音波パルスが、14ビートの間、10万メガ−サンプル/秒の速度で超音波受信器により受信される、請求項1の方法。
  11. 哺乳類の被験体の頭蓋内圧(「ICP」)を測定するための方法であって、
    a、前記被験体の頭蓋の標的に少なくとも1つの超音波パルスを送信するステップと、
    b、前記送信した超音波信号の反射信号を検出するステップと、
    前記反射信号を使用し、以下の数式によって前記頭蓋内圧を算出するステップとを含み、
    Σtanh(“ΣIxW+b”)W+b
    上記数式において、Iが前記反射信号に基づくデータポイントの入力行列であり、Wが調整用の重み行列であり、bが割り当てられたランダムバイアス定数である方法。
  12. 前記標的が、前記被験体の頭蓋腔内にあるまたは隣接する少なくとも1つの半剛性または剛性の構造である、請求項11の方法。
  13. 前記標的が、前記頭蓋の後頭部分を含む請求項11の方法。
  14. 前記標的が、第3脳室を含む請求項11の方法。
  15. a、超音波プローブと、
    b、前記超音波プローブに接続された超音波獲得システムと、
    c、前記超音波獲得システムからの信号を受け取り取り、請求項1乃至14の何れか1項に記載の方法を実行するように構成されたマイクロプロセッサシステム、を備える頭蓋内圧モニターシステム。
JP2011521326A 2008-07-30 2009-07-30 頭蓋内の弾性を測定する方法 Expired - Fee Related JP5740547B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US8482708P 2008-07-30 2008-07-30
US61/084,827 2008-07-30
PCT/US2009/052263 WO2010014815A1 (en) 2008-07-30 2009-07-30 Method for measuring intracranial elasticity

Publications (2)

Publication Number Publication Date
JP2011529744A JP2011529744A (ja) 2011-12-15
JP5740547B2 true JP5740547B2 (ja) 2015-06-24

Family

ID=41610728

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2011521326A Expired - Fee Related JP5740547B2 (ja) 2008-07-30 2009-07-30 頭蓋内の弾性を測定する方法

Country Status (11)

Country Link
US (2) US9028416B2 (ja)
EP (1) EP2303137B1 (ja)
JP (1) JP5740547B2 (ja)
KR (1) KR20110036820A (ja)
CN (1) CN102112061B (ja)
AU (1) AU2009276486B2 (ja)
BR (1) BRPI0916409A2 (ja)
CA (1) CA2732411C (ja)
IL (1) IL210839A (ja)
RU (1) RU2011102338A (ja)
WO (1) WO2010014815A1 (ja)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ITGE20110090A1 (it) * 2011-08-10 2013-02-11 Esaote Spa Dispositivo e metodo per la misurazione di parametri di elasticita' di un corpo in esame mediante ultrasuoni
CN103654760B (zh) * 2012-09-10 2016-08-03 焦文华 无创颅内压测量方法及应用该方法的无创颅内压分析仪
CN103190930B (zh) * 2013-04-19 2014-10-29 重庆大学 基于超声波声弹性效应的颅内压监测仪
EP4260803A3 (en) * 2014-06-11 2024-01-17 Nihon Kohden Corporation Apparatus for detecting increase in intracranial pressure
US10863910B2 (en) 2015-05-29 2020-12-15 Nisonic As Method for detecting pulsatile dynamics of the optic nerve sheath, diagnostic methods, medical uses, non-invasive markers, systems and transducer devices
WO2016205824A1 (en) 2015-06-19 2016-12-22 Neural Analytics, Inc. Transcranial doppler probe
US11589836B2 (en) 2016-01-05 2023-02-28 Novasignal Corp. Systems and methods for detecting neurological conditions
JP2019504670A (ja) 2016-01-05 2019-02-21 ニューラル アナリティクス、インコーポレイテッド 臨床指標を決定するためのシステム及び方法
JP2019500155A (ja) 2016-01-05 2019-01-10 ニューラル アナリティクス、インコーポレイテッド 一体型プローブ構造
CN106055853A (zh) * 2016-07-29 2016-10-26 四川大学 无初值依赖和角度依赖的正交各向异性弹性系数反演方法
CN108241840A (zh) * 2016-12-26 2018-07-03 河南农业大学 一种刺吸电位图谱(epg)波形自动识别方法
CN108986909B (zh) * 2018-06-29 2020-06-12 清华大学 基于超声弹性成像的软组织弹性和粘弹性表征方法及装置
CN115861295A (zh) * 2023-02-09 2023-03-28 南京左右脑医疗科技集团有限公司 大脑中线结构识别方法、装置和存储介质

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5388583A (en) 1993-09-01 1995-02-14 Uab Vittamed Method and apparatus for non-invasively deriving and indicating of dynamic characteristics of the human and animal intracranial media
IL119623A0 (en) * 1996-11-15 1997-02-18 Inta Medics Ltd Non-invasive real time diagnosis of migraine
US6231509B1 (en) 1997-12-05 2001-05-15 Royce Johnson Apparatus and method for monitoring intracranial pressure
US6475147B1 (en) * 1999-01-27 2002-11-05 The United States Of America As Represented By The United States National Aeronautics And Space Administration Ultrasonic apparatus and technique to measure changes in intracranial pressure
JP2003521284A (ja) * 1999-05-10 2003-07-15 インタ−メディクス リミテッド 頭蓋内圧の非侵襲性モニタリング
US6702743B2 (en) * 2000-05-26 2004-03-09 Inta-Medics, Ltd. Ultrasound apparatus and method for tissue resonance analysis
US6328694B1 (en) * 2000-05-26 2001-12-11 Inta-Medics, Ltd Ultrasound apparatus and method for tissue resonance analysis
US20060079773A1 (en) * 2000-11-28 2006-04-13 Allez Physionix Limited Systems and methods for making non-invasive physiological assessments by detecting induced acoustic emissions
AU2002239360A1 (en) 2000-11-28 2002-06-11 Allez Physionix Limited Systems and methods for making non-invasive physiological assessments
DE60329617D1 (de) * 2002-10-18 2009-11-19 Global Monitors Inc Dichte/solut-monitor von multimodalitäten und signalverarbeitungsschema
EP1633234A4 (en) * 2003-06-03 2009-05-13 Physiosonics Inc SYSTEMS AND METHODS FOR DETERMINING INTRACRANIAL FA NON-INVASIVE PRESSURE AND ACOUSTIC TRANSDUCER ASSEMBLIES FOR USE THEREIN
US7250882B2 (en) 2005-06-25 2007-07-31 Georgia Tech Research Corporation High speed data converter testing devices, methods, & systems
US8926515B2 (en) * 2008-05-15 2015-01-06 Uab Vittamed Method and apparatus for continuously monitoring intracranial pressure

Also Published As

Publication number Publication date
AU2009276486B2 (en) 2015-08-06
EP2303137B1 (en) 2015-09-09
CN102112061A (zh) 2011-06-29
IL210839A0 (en) 2011-04-28
CN102112061B (zh) 2014-06-18
US9028416B2 (en) 2015-05-12
US20120136249A1 (en) 2012-05-31
RU2011102338A (ru) 2012-09-10
KR20110036820A (ko) 2011-04-11
US9724068B2 (en) 2017-08-08
CA2732411C (en) 2017-02-28
CA2732411A1 (en) 2010-02-04
WO2010014815A1 (en) 2010-02-04
IL210839A (en) 2017-08-31
JP2011529744A (ja) 2011-12-15
US20160095574A1 (en) 2016-04-07
AU2009276486A1 (en) 2010-02-04
BRPI0916409A2 (pt) 2016-02-16
EP2303137A1 (en) 2011-04-06
EP2303137A4 (en) 2012-09-12

Similar Documents

Publication Publication Date Title
JP5740547B2 (ja) 頭蓋内の弾性を測定する方法
JP3875581B2 (ja) 超音波診断システム
CA2428872C (en) Systems and methods for making non-invasive physiological assessments
US10172527B2 (en) Method and apparatus for measuring a physical parameter in mammal soft tissues by propagating shear waves
US7022077B2 (en) Systems and methods for making noninvasive assessments of cardiac tissue and parameters
US7547283B2 (en) Methods for determining intracranial pressure non-invasively
US20060079773A1 (en) Systems and methods for making non-invasive physiological assessments by detecting induced acoustic emissions
US6387051B1 (en) Method and apparatus for non-invasively deriving and indicating of dynamic characteristics of the human and animal intracranial media
EP1806099A1 (en) Ultrasonographic device and ultrasonographic device control method
JP3842285B2 (ja) 超音波診断装置および超音波診断装置の制御方法
JP2005532097A (ja) 心組織およびパラメータの非侵襲的評価システムおよび方法
JP5014132B2 (ja) 超音波診断装置
JP7393435B2 (ja) 頭蓋内コンプライアンスの非侵襲的な評価及び/又は監視方法及び装置
Corazza Shear wave elastography to assess the effect of botulinum toxin in muscle hypertonia following stroke
Hongo et al. Estimation of the location of intracranial vascular diseases using several sensors

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20120709

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20130314

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20130315

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130419

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20130829

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130902

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20131127

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20131204

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20131227

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20140801

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140901

R150 Certificate of patent or registration of utility model

Ref document number: 5740547

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees