(57)【要約】  二つの連続する心拍の間の時間間隔（δｔｉ）又は前記時間間隔の逆数（１／δｔｉ）を表わす複数の（ＲＲ）サンプルを含むＲＲシリーズを処理するための本発明の方法は、（Ｎ）サンプル（ＲＲｉ）が予め決められた時間長さ（ｎ）を有するメイン時間ウインドゥにおいて選択されること、前記メインウインドゥが（ｍ）サブウインドゥ（Ｆｊ）に切り分けられること、サブウインドゥ（Ｆｊ）に含まれるサンプル（ＲＲｉ）に基づいて中間パラメータ（Ａｊ）が各サブウインドゥ（Ｆｊ）について計算されること、及び最終パラメータが中間パラメータ（Ａｊ）の関数として計算されることにある。本発明は、心調律の変動性を分析するために、特に患者の痛み又はストレスを評価するために使用される。【選択図】  図５
2008513073.tif 000002
【特許請求の範囲】
【請求項１】
二つの連続する心拍を分離する時間間隔（δｔｉ）又は前記時間間隔の逆数（１／δｔｉ）を表わす複数のサンプル（ＲＲｉ）からなるＲＲシリーズを処理するための方法において、（Ｎ）サンプル（ＲＲｉ）が予め決められた時間長さ（ｎ）を有するメイン時間ウインドゥにおいて選択されること、前記メインウインドゥが（ｍ）サブウインドゥ（Ｆｊ）に切り分けられること、サブウインドゥ（Ｆｊ）に含まれるサンプル（ＲＲｉ）に基づいて中間パラメータ（Ａｊ）が各サブウインドゥ（Ｆｊ）について計算されること、及び最終パラメータが中間パラメータ（Ａｊ）の関数として計算されることを特徴とする方法。
【請求項２】
最終パラメータが、メインウインドゥの時間長さ（ｎ）未満の予め決められた時間間隔（ｐ）だけメイン時間ウインドゥをシフトさせることによって、反復によって計算されることを特徴とする請求項１に記載の方法。
【請求項３】
シリーズ（ＲＲ）のサンプル（ＲＲｉ）が二つの連続する心拍を分離する時間間隔（δｔｉ）を表わす場合、各中間パラメータ（Ａｊ）は少なくともサブウインドゥ（Ｆｊ）におけるＲＲシリーズの最小値（Ｐｉ）の関数であり、シリーズ（ＲＲ）のサンプル（ＲＲｉ）が二つの連続する心拍を分離する時間間隔（δｔｉ）の逆数（１／δｔｉ）を表わす場合、各中間パラメータ（Ａｊ）は少なくともサブウインドゥ（Ｆｊ）におけるＲＲシリーズの最大値（Ｐｉ）の関数であることを特徴とする請求項１又は２に記載の方法。
【請求項４】
最終パラメータのための計算された値（ＡＵＣｍａｘ）が中間パラメータ（Ａｊ）の最大値に比例することを特徴とする請求項１〜３のいずれか一項に記載の方法。
【請求項５】
最終パラメータのための計算された値（ＡＵＣｍｏｙ）が中間パラメータ（Ａｊ）の平均値に比例することを特徴とする請求項１〜３のいずれか一項に記載の方法。
【請求項６】
中間パラメータ（Ａｊ）の計算のために包絡線（Ｃｊ）が規定され、シリーズ（ＲＲ）のサンプル（ＲＲｉ）が二つの連続する心拍を分離する時間間隔（δｔｉ）を表わす場合、前記包絡線は測定された最小点（Ｐｉ）に連結し、シリーズ（ＲＲ）のサンプル（ＲＲｉ）が二つの連続する心拍を分離する時間間隔（δｔｉ）の逆数（１／δｔｉ）を表わす場合、前記包絡線は測定された最大点（Ｐｉ）を連結することを特徴とする請求項１〜５のいずれか一項に記載の方法。
【請求項７】
各中間パラメータ（Ａｊ）は包絡線（Ｃｊ）によって画定される面積の関数であることを特徴とする請求項６に記載の方法。
【請求項８】
パラメータの計算に先立ち、０．１Ｈｚより大きいか又は０．１Ｈｚと等しい、好ましくは０．１Ｈｚ〜０．１５Ｈｚの間にある遮断周波数を有する高域フィルタによってＲＲシリーズがろ波されることを特徴とする請求項１〜７のいずれか一項に記載の方法。
【請求項９】
中間パラメータ（Ａｊ）の計算に先立ち、各サブウインドゥ（Ａｊ）においてＲＲシリーズのサンプル（ＲＲｉ）がメインウインドゥの全幅にわたって正規化されることを特徴とする請求項１〜８のいずれか一項に記載の方法。
【請求項１０】
心調律の変動性を分析するためのシステムであって、アナログ心臓信号を獲得する手段（１，２）、前記心臓信号をサンプリングする手段（４）、及び二つの連続する心拍を分離する時間間隔（δｔｉ）又は前記時間間隔の逆数（１／δｔｉ）を表わす複数のサンプル（ＲＲｉ）からなるＲＲシリーズを構築するように設計されている、サンプリングされた信号を処理する手段（５）を含むものにおいて、前記処理手段（５）が、請求項１〜９のいずれか一項に記載の方法に従って少なくとも一つの最終パラメータをシリーズ（ＲＲ）から自動的に計算するようにさらに設計されていることを特徴とするシステム。
【請求項１１】
痛みの評価のための請求項１０に記載のシステムであって、最終的に計算されたパラメータが痛みのレベルを特徴付けることを特徴とするシステム。
【請求項１２】
ストレスの評価のための請求項１０に記載のシステムであって、最終的に計算されたパラメータがストレスのレベルを特徴付けることを特徴とするシステム。
【請求項１３】
生物によって感じられている痛みを評価するための請求項１０に記載の分析用システムの使用。
【請求項１４】
生物によって感じられているストレスを評価するための請求項１０に記載の分析用システムの使用。
【請求項１５】
生物の心調律の変動性を分析するための方法であって、以下の主要段階を含むことを特徴とする方法：
−  生物のアナログ心臓信号を獲得する、
−  前記心臓信号をサンプリングし、二つの連続する心拍を分離する時間間隔（δｔｉ）又は前記時間間隔の逆数（１／δｔｉ）を表わす複数のサンプル（ＲＲｉ）からなるＲＲシリーズを構築する、
−  請求項１〜９のいずれか一項に記載の方法に従ってＲＲシリーズを処理する。
【請求項１６】
最終的に計算されたパラメータが痛みのレベルを特徴付けることを特徴とする、生物によって感じられている痛みの評価のための請求項１５に記載の方法の使用。
【請求項１７】
最終的に計算されたパラメータがストレスのレベルを特徴付けることを特徴とする、生物によって感じられているストレスの評価のための請求項１５に記載の方法の使用。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、生物における心調律の変動性を分析するための一般的な方法に関し、及び生物、特に硬膜外型の局所麻酔下で意識がある生物又は全身麻酔下で意識がない生物における痛みを評価するための及び生物によって感じられているストレスを評価するためのその使用に関する。
【背景技術】
【０００２】
生理学的な見地からは、生物の心臓は全ての外部の影響から隔離された場合、洞房結節の作用のためメトロノームのように極めて規則的な速度で自動的に収縮する。洞房結節は独立した神経インパルスを生成し、この作用を通して心筋の自発的な収縮を引き起す。それにもかかわらず、心臓は隔離されておらず、副交感神経系及び交感神経系を介して自律神経系（ＡＮＳ）に接続されている。この自律神経系は心臓の活動に影響を及ぼす：交感系は心調律を加速させ、一方、副交感系は心調律を低下させる。
【０００３】
従って、ある程度は自律的であるが、心臓は自律神経系の影響に供されており、このことは特に生物の体にとって心調律をその必要とする機能に適合させながらも合理的な限界内にとどまらせることを可能にする。結果として、心調律における、特に心調律の変動性（二つの心拍を分離する時間間隔の変動性）における時間にわたる変化の分析は、心臓系の活動に関する、特に自律神経系の活動に関する重要な情報が得られることを可能にする。ＡＮＳの活動に関する知識は多数の臨床状況において診断に到達するための価値ある助けであり得る。この主題において、以下の文献が例えば参照されることができる：Ｌａｃｒｏｉｘ  Ｄ，Ｌｏｇｉｅｒ  Ｒ，Ｋａｃｅｔ  Ｓ，Ｈａｚａｒｄ  Ｊ−Ｒ，Ｄａｇａｎｏ  Ｊ（１９９２）：“Ｅｆｆｅｃｔｓ  ｏｆ  ｃｏｎｓｅｃｕｔｉｖｅ  ａｄｍｉｎｉｓｔｒａｔｉｏｎ  ｏｆ  ｃｅｎｔｒａｌ  ａｎｄ  ｐｅｒｉｐｈｅｒａｌ  ａｎｔｉｃｈｏｌｉｎｅｒｇｉｃ  ａｇｅｎｔｓ  ｏｎ  ｒｅｓｐｉｒａｔｏｒｙ  ｓｉｎｕｓ  ａｒｒｈｙｔｈｍｉａ  ｉｎ  ｎｏｒｍａｌ  ｓｕｂｊｅｃｔｓ”，Ｊｏｕｒｎａｌ  ｏｆ  ｔｈｅ  Ａｕｔｏｎｏｍｉｃ  Ｎｅｒｖｏｕｓ  Ｓｙｓｔｅｍ，Ｖｏｌ  ３９，ｐｐ２１１−２１８。
【０００４】
これらの心調律の変動及びその結果としてのＡＮＳの活動を研究するため、瞬間的な心調律（又は周波数）における時間にわたる変化を表わす信号のスペクトル分析について様々な技術が１９７０年から開発されてきた。この信号はアナログ生物電気信号のサンプリング後に得られるものであり、生物の心調律を特徴付け、以下では心臓信号と称する。
【０００５】
心調律の変動性を分析する現在公知の主要な方法は以下のものにある：
−  いかなる侵襲的又は非侵襲的手段によって心臓信号を獲得する（例えば、心電計による心電計（ＥＣＧ）信号の獲得、又は動脈中に挿入されるカテーテルに接続された血圧センサの使用、又は赤外パルスセンサの使用であり、これらの例は非限定的である）、
−  この信号から、二つの連続する心拍の間の時間間隔を表わす複数のサンプル（ＲＲｉ）からなるいわゆるＲＲシリーズを構築する、
−  そしてＲＲシリーズのスペクトル分析を実行する。
【０００６】
特に、心臓信号から由来するＲＲシリーズのスペクトル分析は二つの主要段階で通常実行される。
【０００７】
第一段階では、ＲＲシリーズのスペクトル密度についての曲線が様々な公知手段を使用して例えば０Ｈｚ〜２Ｈｚで計算される。最も一般的に使用される方法は、予め規定された重みづけウインドゥによって重みづけられた予め規定された時間ウインドゥにおけるＲＲシリーズの離散高速フーリエ変換を計算することである。予測された実行に依存して、これは長方形の重みづけウインドゥであることができるか、又は例えばカイザー、ハミング、ハニング又はバーレット重みづけウインドゥであることができる。同様に、計算時間ウインドゥは予め規定されかつ定数であることもでき、又は計算時間ウインドゥは予め決められた寸法であるが時間にわたってスライドさせることもできる。例えば、フーリエ変換は、２５６秒のスライドする時間ウインドゥにおいて実行され、ＲＲシリーズに適用され、カイザー重みづけに供されてウインドゥイングによるエッジ効果を制限する。
【０００８】
第二段階では、第一段階の最後で得られたスペクトル密度曲線から、予め決められかつ所望によりユーザが調整可能な周波数境界の間で自動的にスペクトルパワー（スペクトル密度曲線の下の面積）が計算される。かかるスペクトルパワー計算は自律神経系（ＡＮＳ）の活動を特徴付ける定量的な情報が得られることを可能とし、従ってＡＮＳによる心臓の調節を調査して分析する手段を構成する。
【０００９】
しかしながら、上述のスペクトル分析はいくつかの欠点を有する。
【００１０】
ＲＲシリーズは非静的である。結果として、この種類の非静的シリーズに対するフーリエ変換の適用は、不明確でたぶん不正確ですらあり、しかも確実性を持って解釈されることができない結果を生ずる。高速フーリエ変換（又は等価物）によってスペクトル密度カーブを計算することは、コンピュータ処理力及び／又はコンピュータ処理時間の見地から比較的求められており、これは現在ではスペクトル分析のこの方法を不適切にし、特に携帯型システムにおいて実時間でのこの方法の実行を困難にする。
【００１１】
同様に、許容可能な周波数解像度を得るために、高速フーリエ変換は、ＲＲシリーズの有意な数のサンプルに対応する比較的大きな時間ウインドゥ（例えば２５６秒）にわたって計算されなければならない。結果として、スペクトル分析のこの方法は記憶効果を伴ない、これはＲＲシリーズに生じるいかなる変化の認識を遅延させる。
【発明の開示】
【発明が解決しようとする課題】
【００１２】
その主目的として、本発明はＲＲシリーズを自動的に処理する新規の方法を提案し、この方法は、スペクトル分析に基づく先行技術の方法に内在する上述の欠点を緩和し、高度の感度を持ってＡＮＳの活動を特徴付ける定量的値（パラメータ）を計算することを可能にする。
【課題を解決するための手段】
【００１３】
本発明によるＲＲシリーズを処理するための方法は、二つの連続する心拍を分離する時間間隔（δｔｉ）又は前記時間間隔の逆数（１／δｔｉ）を表わすサンプル（ＲＲｉ）の時間分析方法に基づいている。
【００１４】
この方法は、予め決められた長さ（ｎ）を有するメイン時間ウインドゥにおいて（Ｎ）サンプル（ＲＲｉ）を選択すること、前記メインウインドゥを（ｍ）サブウインドゥ（Ｆｊ）に切り分けること、サブウインドゥ（Ｆｊ）に含まれるサンプル（ＲＲｉ）に基づいて中間パラメータ（Ａｊ）を各サブウインドゥ（Ｆｊ）について計算すること、及び最終パラメータを中間パラメータ（Ａｊ）の関数として計算することにある。
【００１５】
本発明は、その第二の目的として、心調律の変動性を分析するためのシステムを有し、この装置はアナログ心臓信号を獲得する手段、前記心臓信号をサンプリングする手段、及び二つの連続する心拍を分離する時間間隔（δｔｉ）又は前記時間間隔の逆数（１／δｔｉ）を表わす複数のサンプル（ＲＲｉ）からなるＲＲシリーズを構築するように設計されている、サンプリングされた信号を処理する手段を含み、前記処理手段が、ＲＲシリーズを処理するための上述の方法に従って少なくとも一つのパラメータをシリーズ（ＲＲ）から自動的に計算するようにさらに設計されている。
【００１６】
本発明の第三の目的は、生物の心調律の変動性を分析するための方法である。この方法は以下の主要段階を含む：
−  生物のアナログ心臓信号を獲得する、
−  前記心臓信号をサンプリングし、二つの連続する心拍を分離する時間間隔（δｔｉ）又は前記時間間隔の逆数（１／δｔｉ）を表わす複数のサンプル（ＲＲｉ）からなるＲＲシリーズを構築する、
−  上述の方法に従ってＲＲシリーズを処理する。
【００１７】
本明細書において、及び特に特許請求の範囲において、用語「心臓信号」は生物の瞬間的な心調律（又は周波数）を特徴付けるいかなる物理的信号に言及する。本発明を実行するために、様々な侵襲的又は非侵襲的技術が前記心臓信号を獲得するために使用されることができる。公知の侵襲的技術は、例えば、動脈中に挿入されるカテーテルに接続された血圧センサを使用することにある。公知の非侵襲的方法（こちらの方が好ましい）は、例えば赤外パルスセンサの使用、心臓周期の測定を可能にする超音波センサの使用、心電計で使用される種類のセンサの使用、又は心電計（ＥＣＧ）信号の獲得を含む。心電計（ＥＣＧ）信号の獲得は実用的に最も一般的に使用される方法である。何故なら、その非侵襲性の特徴に加えて、この方法は例えば赤外パルスセンサによって得られる心臓信号よりもずっと正確な心臓信号の獲得を可能にするからである。
【００１８】
本明細書において、及び特に特許請求の範囲において、用語「ＲＲシリーズ（ｓｅｒｉｅｓ）」は、生物の心調律を特徴付けるアナログ心臓信号のサンプリング後に得られた幾つかの連続サンプル（ＲＲｉ）のシリーズを一般的に意味し、各サンプル（ＲＲｉ）は、二つの連続する心拍の間の時間間隔（δｔｉ）又は前記時間間隔の逆数（１／δｔｉ）を一般的な態様で特徴とする。
【００１９】
添付の図面を参照して以下に既述される好ましい実施態様では、このＲＲシリーズはＥＣＧ信号のＲ波から一層特異的に構築される。しかし、これは本発明の範囲を限定しない。ＥＣＧタイプの心臓信号の場合、いわゆるＲＲシリーズは、ＲＲシリーズを構築するためのＥＣＧ信号の他の復極波（Ｐ，Ｑ，Ｓ又はＴ）を使用することによって構築されることもできるが、しかしながら、これはＥＣＧ信号のＲ波を使用した場合より正確さが劣る。同様に、心臓信号がＥＣＧ信号でない場合、ＲＲシリーズのサンプルは、ＥＣＧ信号の二つの連続するＲ波を分離する時間間隔を決定することによって計算されるのではなく、代わりに、心臓信号において二つの連続する心拍間の時間間隔を測定することによって一層一般的な態様で決定される。
【００２０】
本発明の方法又はシステムによって計算される最終パラメータは、いかなる刺激に影響を与えるＡＮＳ活動であって心調律（又は周波数）の変動性を生じるＡＮＳ活動の特徴付けを一般的な態様で可能とする。
【００２１】
本発明の第一の重要な用途は医療又は外科手術の分野に見出され、生物の痛みのレベルを評価することである。
【００２２】
医療又は外科手術の分野では、患者によって感じられている痛みのレベルの知識を有することが、特に前記痛みに対する最良の可能な治療を可能にしそして提供することができるようにするために実際に重要である。
【００２３】
現在までに知られている最も一般的な痛み評価方法は、患者が感じている痛みのレベルを予め確立された痛みスケール上でカーソル又はその等価物によって患者が指すことを可能にするスライディングスケールの使用に基づく主観的な方法である。この方法は、一方では完全に主観的である故全く信頼できないという欠点を有し、他方では意識のある患者に対してしか使用するのに好適でないという欠点を有する。
【００２４】
さらに、局所麻酔下で意識のある患者の場合（例えば硬膜外型の局所麻酔下の患者の場合）、上述の主観的評価方法は麻酔剤投与の実時間での自動的制御のために使用されることができない。
【００２５】
全身麻酔の分野では、意識がない患者によって感じられている痛みのレベルを知ることができることも関心がある。前記痛みのレベルの知識は鎮痛剤投与の最良の可能な適合を可能にする。
【００２６】
上述の理由のため、生物によって感じられている痛みを評価する客観的な方法を提案することにおいてかなりの関心がある。前記方法は、痛みのレベルの尺度であるパラメータの自動的な計算を生じる必要がある。
【００２７】
従って、本発明は、他の目的として、生物によって感じられている痛みを評価するための、上述の分析システムの使用並びに心調律の変動性を分析する上述の方法の応用を有し、この場合、最終的に計算されたパラメータが痛みのレベルを特徴付ける。
【００２８】
文脈に応じて、本発明によれば、痛みの評価は、機械的に換気されていない意識のある生物（即ち、全身麻酔下の患者とは特に対照的に呼吸頻度が正常で変動性であり、制御された換気装置によって支配されていない生物）に対して並びに意識のない生物（特に、全身麻酔下の生物）に対して実行されることができる。後者（全身麻酔下の生物）の場合、痛みのレベルの評価は全身麻酔中の痛覚脱失のレベルに関する間接的な知識を提供する。
【００２９】
本発明の第二の重要な用途は、生物のストレスのレベルを評価するための準医療活動の分野に見出されることができる。
【００３０】
図面の記述
本発明の他の特徴及び利点は、非限定的実施例示の詳細な記述を読むことによって、及び添付の図面を参照することによって一層明確になるであろう。図中、
−  図１は、本発明による分析システムの一例の主要要素の概略を示す。
−  図２は、ＥＣＧ信号における心拍を特徴付ける一組の波（ＰＱＲＳＴ）を示す。
−  図３は、アナログＥＣＧ信号のサンプリング後に得られたディジタルＥＣＧ信号の一例を示す。
−  図４は、図３の信号から構築されたＲＲシリーズ（ＲＲ信号とも称する）の一例を示す。
−  図５は、ろ波及び正規化後のＲＲシリーズの一例を示す。
【発明を実施するための最良の形態】
【００３１】
心調律の変動性の分析システム
図１は、本発明による、生物（以下、「患者」と称する）の心臓信号の獲得及び処理を可能にする心調律の変動性の分析システムの一例を示す。
【００３２】
このシステムは以下のものを含む：
−  心電計（ＥＣＧ）モニタ２の入力に接続されたいくつかの測定電極１を含む、ＥＣＧ信号を獲得するための従来の手段、
−  ＥＣＧモニタ２の出力によって与えられるＥＣＧ信号を処理するための手段３。
【００３３】
ＥＣＧ信号の処理手段３は、アナログ−ディジタル変換器４及びプログラムされた処理ユニット５を含む。変換器４の入力はＥＣＧモニタ２の出力に接続されており、変換器４の出力は処理ユニット５の入力ポートに接続されている。本発明の特定の非限定的実施態様では、処理ユニット５はマイクロコンピュータからなり、変換器４は前記マイクロコンピュータの通信ポート（例えばＲＳ２３２シリアルポート又はＵＳＢポート）に接続されている。
【００３４】
作動モードでは、電極１は、患者の体に付与され、ＥＣＧモニタは、ＥＣＧ信号と呼ばれるアナログ電気信号を従来公知の態様で出力し、このＥＣＧ信号は各心拍についてのものであり、例えば図２に示されるような信号の形態を有する。
【００３５】
図２を参照すると、各心拍について、この心電計（ＥＣＧ）信号は一組の電気波からなる：
−  Ｐ波。これは房の復極に対応し、低い振幅及びドーム形状を表す；
−  ＰＱ空間。これは房心間伝導時間を表す；
−  Ｒ波。これは実際には室の心収縮又は心拍のマーカと考えられており、ＱＲＳ複合体は室の収縮を反映する；そして
−  Ｔ波。これは室の再分極を反映する。
【００３６】
このアナログＥＣＧ信号は変換器４によって、例えば２５６Ｈｚの値の予め規定されたサンプリング周波数（ｆｃ）でディジタル化される。
【００３７】
変換器４の出力によって与えられたサンプリング信号（図３に示される信号）は、以下に詳述する特定の処理ソフトウェア（痛み評価ソフトウェア）を使用して処理ユニット５によって処理される。この痛み評価ソフトウェアは、処理ユニット５のメモリに記憶されており、実行されたとき、アナログ−ディジタル変換器４によって送達されるディジタル信号から患者の痛みのレベルの尺度である二つの別個のパラメータ（ＡＵＣｍａｘ及びＡＵＣｍｏｙ）を自動的に計算することを可能にする。
【００３８】
この痛み評価ソフトウェアの好ましい変形例は以下に詳述される。
【００３９】
痛み評価ソフトウェアのアルゴリズムの一例
本発明の好ましい実施態様では、痛み評価ソフトウェアのアルゴリズムの主要連続段階は以下のようである：
１．ＲＲｉサンプルの獲得及び予め決めされた周波数ｆでの再サンプリング
２．ｎ秒（ｎ＞１／ｆ）のメイン時間ウインドゥに含まれるＲＲｉサンプルの選択
３．ろ波
４．信号の正規化
５．メインウインドゥにおける最小値（Ｐｉ）の測定
６．メイン時間ウインドゥのｍ（ｍ≧２）個のサブウインドゥ（Ｆｊ）への切り分け、及び各サブウインドゥ（Ｆｊ）について以前の段階で測定された最小値（Ｐｉ）を連結する下方の包絡線（Ｃｊ）によって画定される面積（Ａｊ）の計算
７．ＡＵＣｍｏｙ及びＡＵＣｍａｘパラメータの計算
８．ｐ秒（ｐ≦ｎ）の長さの時間間隔によってメイン時間ウインドゥをｎ秒だけシフトさせ、段階２の計算を繰り返す。
【００４０】
実際には、システムは実時間で又は遅延時間で使用されるようにプログラムされることができる。
【００４１】
システムが遅延時間で使用される場合、段階１は第一ステップでは全体の望ましい分析期間にわたって全てのＲＲｉサンプルを得るために実時間で実行される；次に、これらの連続的なＲＲｉサンプルはそれらの実体が処理ユニットのメモリ中の獲得ファイルに記憶される。第二ステップでは、獲得ファイルに記憶されているＲＲｉ間隔値を使用して遅延時間で段階２〜８が連続ループで実行される。
【００４２】
システムが実時間で作動する場合、一方では獲得段階１が、そして他方では続く処理段階２〜８が、平行に実行される二つの別個のソフトウェアモジュールによって実行され、これらは、第二処理に供給する第一獲得モジュール（段階１）及びバッファファイル又は等価物による計算モジュール（段階２〜８）である。
【００４３】
段階１〜８は以下に詳述される：
段階１：ＲＲｉサンプルの獲得及び予め決められた周波数ｆでの再サンプリング
ＲＲｉサンプルの獲得は、アナログ−ディジタル変換器４によって与えられるディジタルＥＣＧ信号（図３の信号）を構成する連続ディジタルデータを入力として受取る第一ソフトウェアサブモジュールによって実行される。ＥＣＧ信号の各データ（又は点）は、ＥＣＧ信号の瞬間的振幅ＥＣＧｉによって及びサンプリング瞬間ｔｉ（ｔｉ＝ｎｉ／ｆｃであり、ここでｎｉはサンプル数であり、ｆｃは変換器４のサンプリング周波数を表わす）によって規定される。
【００４４】
ＲＲｉサンプルの獲得のための第一サブモジュールは、変換器４によって与えられるディジタル信号におけるそれぞれの連続Ｒｉピークを自動的に測定してＲＲｊサンプルの連続からなるＲＲシリーズ（図４）を自動的に構築するように設計されている。各ＲＲｊサンプルは以下の一対の同格物によって規定される：ｔｉ（サンプリング瞬間（又は数））、ピークＲｉを続くピークＲｊ＋１から分離する（別の変形例ではこれは以前のピークＲｊ−１にも当てはまることができる）時間間隔δｔｉ（例えばサンプリング期間ｆｃの倍数で表される）。
【００４５】
Ｒ波は従来ＱＲＳ複合体の最高振幅を有する最も狭い部分であるので、Ｒ波は極めて高い正確さで心拍を測定するのに好ましく使用される波であり、時間間隔δｔｉは実際には二つの連続する心拍を分離する時間に対応する。それにもかかわらず、別の変形例では、ＲＲシリーズを測定して構築するためにＥＣＧ信号の心拍の他の波（例えばＱ波又はＳ波）を使用することを予見することができる。
【００４６】
上述の第一サブモジュールによって与えられるＲＲシリーズ（図４）は、予め決められた周波数ｆで第二ソフトウェアサブモジュールによって自動的に再サンプリングされる。この再サンプリング周波数ｆは、患者の最大生理学的心周波数の二倍より大きくなければならない。すべての起こりうる生理学的状況を考慮に入れるため、再サンプリング周波数ｆは６Ｈｚより大きい値に固定されることが好ましい。
【００４７】
この再サンプリングの目的は、ＲＲｉサンプルが時間の観点から見て等距離にあるＲＲシリーズ（即ち、換言すればサンプリング瞬間が規則的であるＲＲシリーズ）を出力することである。この再サンプリングは、例えば線形補間法のような補間法によって公知の態様で行われる。
【００４８】
段階２：ｎ秒（ｎ＞１／ｆ）のメイン時間ウインドゥに含まれるＲＲｉサンプルの選択
この段階はＮ個の連続的ＲＲｉサンプル（Ｎ＝ｎ．ｆ）を単離することにある。概略としては、例えば８Ｈｚの再サンプリング周波数を有する５１２個の連続的ＲＲｉサンプル（Ｎ＝５１２）に対応する６４秒（ｎ＝６４）のメインウインドゥが選択される。
【００４９】
以下の段階３〜７はこのメインウインドゥに含まれるサンプルに適用される。
【００５０】
段階３：ろ波／［ｆ１；ｆ２］
この段階は、予め規定された周波数帯域幅［ｆ１；ｆ２］内に存在する周波数のみを維持するために、メインウインドゥ内に含まれるＲＲシリーズのサンプルに帯域通過フィルタ（ろ波器）を適用することにある。
【００５１】
特に、周波数帯域幅［ｆ１；ｆ２］は、周波数帯域幅［０．０５Ｈｚ；５Ｈｚ］に等しいか又はこの周波数帯域幅内に存在する。周波数帯域幅［ｆ１；ｆ２］は［０．１Ｈｚ；１Ｈｚ］に等しいことが好ましい。
【００５２】
帯域通過ろ波のこの段階を実行するために、例えば周波数ｆ１で遮断周波数を有するディジタル高域フィルタであって、周波数ｆ２で遮断周波数を有するディジタル低域フィルタと直列に接続されたディジタル高域フィルタが使用される。周波数帯域幅［ｆ１；ｆ２］にセンターリングされた無限インパルス応答を有する再帰型選択的フィルタ（ＩＩＲフィルタ）を使用することも可能である。
【００５３】
高域フィルタ（遮断周波数ｆ１）は、０．１Ｈｚより下の低周波数のろ波及び同時に信号中の多数のアーチファクトの除去をその目的として有する。従って、実際、遮断周波数は０．１Ｈｚを超えるか又は０．１Ｈｚに等しく、０．１Ｈｚ〜０．１５Ｈｚの間にあることが好ましい。それは信号の平均値を除去することができるという利点も有する。高域フィルタの使用を省略することも可能である。この場合、中間パラメータ（Ａｊ）を計算する前に信号をその平均に再センターリングすることが好ましい。
【００５４】
低域フィルタ（遮断周波数ｆ２）は、高周波数、典型的には１Ｈｚを超える高周波数のろ波をその目的とする。何故なら、実際にはこれらの高周波数はいかなる関心のある情報も含んでいないからである。この低域フィルタの実行は好ましいが、それにもかかわらず任意であり、本発明の実行には必須でない。
【００５５】
段階４：信号の正規化
この段階は、以下の一般的に使用される式に従って、段階３から生じる信号のノルム（ｎｏｒｍ）Ｓを第一ステップに計算するソフトウェアサブモジュールによって行われる：
2008513073.tif 000003
式中、ｐは１より大きい整数であり、Ｓｉは信号の離散値である。
【００５６】
別の変形例では、ノルムは以下の式によって計算されることもできる：
2008513073.tif 000004
【００５７】
好ましくは、本発明の実行のためには、ｐは２に等しいように選択される。
【００５８】
第二ステップでは、ソフトウェアサブモジュールは、以前に計算されたノルムＳで信号の各値Ｓｉを割ることによって信号の正規化を実行する。
【００５９】
この段階４は、最終結果に関するより良好な感度（痛みのレベルの尺度となるＡＵＣｍｏｙ及びＡＵＣｍａｘパラメータの感度）の獲得を可能にする。
【００６０】
この正規化はメインウインドゥの全幅に対して、即ち、正規化中、メインウインドゥの（Ｎ）ＲＲｉサンプルを考慮に入れて行われることが必須である。対照的に、各サブウインドゥ（Ａｊ）に対して行われる正規化は全く無駄である。
【００６１】
段階５：メインウインドゥにおける最小値の測定
図５は、上述の段階４で得られたＲＲ信号の一例（ろ波及び正規化後の）を示す。この図では、横軸は、ろ波されて正規化された信号のゼロ値に対応する。
【００６２】
段階５は、ソフトウェアサブモジュールによって実行され、このサブモジュールは、例えば勾配の反転（又は等価な態様においては信号導関数の符号の変化）を測定するアルゴリズムによって、この信号の最小値（図５の点Ｐｉ）を測定する。
【００６３】
段階６：メイン時間ウインドゥのｍ（ｍ≧２）個のサブウインドゥ（Ｆｊ）への切り分け、及び各サブウインドゥ（Ｆｊ）について以前の段階で測定された最小値（Ｐｉ）を連結する下方の包絡線（Ｃｊ）によって画定される面積（Ａｊ）の計算
段階６に対応するソフトウェアサブモジュールは、各サブウインドゥ（Ｆｊ）について段階５で測定された最小値（点Ｐｉ）を連結しかつサブウインドゥ（Ｆｊ）に含まれる下方の包絡線（図５の曲線Ｃｊ）を自動的に計算する。この下方包絡線（Ｃｊ）は例えば線形補間法によって、即ち最小点（Ｐｉ）を連結する直線セグメントを計算することによって決定される。
【００６４】
これに続いて、段階６に対応するソフトウェアサブモジュールは、下方包絡線（Ｃｊ）及び横軸（信号ゼロ値）によって画定される面積（Ａｊ）を各サブウインドゥ（Ｆｊ）について計算する。この面積（Ａｊ）は図５の細い平行線を引かれた面積に相当する。
【００６５】
下方包絡線（Ｃｊ）が直線セグメントによって近似される場合、面積（Ａｊ）のこの計算は、二つの連続する最小点（Ｐｉ）と横軸とによって画定される台形面積の合計を計算することを必要とする。
【００６６】
代替的実施態様では、面積（Ａｊ）は、横軸とは異なりかつ好ましくは（しかし必須でない）横軸より下に位置された水平直線、例えば図５の一点直線Ｄ及び下方包絡線（Ｃｊ）によって画定されることができる。同様に、別の実施態様では、面積（Ａｊ）の計算のためのこの直線Ｄは横軸より上に位置させることができ、特に最高縦座標値の点Ｍｉ（図５の特定の例では点Ｍ１０）を通過するように位置させることができる。
【００６７】
非限定的例として、方法は６４秒のメインウインドゥ時間でテストされ、段階６の実行のためにこのメインウインドゥは１６秒の四つ（ｍ＝４）のサブウインドゥＦ１〜Ｆ４に切り分けられた。
【００６８】
段階７：ＡＵＣｍｏｙ及びＡＵＣｍａｘパラメータの計算
二つのＡＵＣｍｏｙ及びＡＵＣｍａｘパラメータは痛みレベルの測定を可能にする。
【００６９】
ＡＵＣｍｏｙは（ｍ）サブウインドゥＦｊの面積Ａｊの平均に等しい（又はより一般的には前記平均に比例する)。
【００７０】
ＡＵＣｍａｘは（ｍ）サブウインドゥＦｊについて計算された面積Ａｊの最大値に等しい（又はより一般的には前記最大値に比例する)。
【００７１】
段階８：ｐ秒（ｐ≦ｎ）の長さの時間間隔によって時間ウインドゥをｎ秒だけシフトさせ、段階２の計算を繰り返す。
メイン計算ウインドゥのスライディング及び段階２〜８の繰り返しのためのステップ（ｐ）は、ＡＵＣｍｏｙ及びＡＵＣｍａｘパラメータの感度に影響を与え、従って感度後の探索（ｔｈｅ  ｓｏｕｇｈｔ  ａｆｔｅｒ  ｓｅｎｓｉｔｉｖｉｔｙ）に依存する。概略的には、及び実際には、２秒の長さのステップ（ｐ）を選択した場合に極めて満足できる結果が得られている。
【００７２】
さらなる実施態様では、上述のろ波段階（段階３）は、例えばＲＲｉサンプルが獲得されるにつれ及びＲＲｉサンプルが獲得されたときに連続的態様で選択段階２の前に実行されることができた。同様に、別の変形例の実行では、このろ波段階は正規化段階（上述の段階４）の後に実行されることができた。
【実施例】
【００７３】
実験的確認
テストは、ＡＵＣｍｏｙパラメータ及びＡＵＣｍａｘパラメータがそれぞれ生物によって感じられている痛い刺激の強度に相関することを示した。従って、これらのパラメータは生物の痛みの強度の測定を可能にする。
【００７４】
上述のＡＵＣｍｏｙ及びＡＵＣｍａｘパラメータと感じられた痛みのレベルとの間のこの相関を確認するため、これらのＡＵＣｍｏｙ及びＡＵＣｍａｘパラメータは、心臓血管病、糖尿病又は神経学的疾患の経歴のない、ＡＳＡ  Ｉ又はＩＩに分類される成人分娩婦の集団に対して、硬膜外痛覚脱失（ＥＡ）下で膣分娩によって臨月で出産する状況についてテストされた。
【００７５】
計算されたＡＵＣｍｏｙ及びＡＵＣｍａｘ値は、硬膜外麻酔の処方前及び処方後の両方でアナログ可視的痛みスケール（ＥＶＡ）を使用して分娩婦によって与えられた痛みレベルの指示と比較された。
【００７６】
この比較は、ＡＵＣｍｏｙ及びＡＵＣｍａｘパラメータが、ＥＶＡを使用して分娩婦によって与えられた痛みレベルの指示と相関することを証明することを可能にした。
【００７７】
ＡＵＣｍａｘパラメータが最も信頼性があり、従って好ましいパラメータであることは特に注目すべきである。
【００７８】
二つのＡＵＣｍｏｙ及びＡＵＣｍａｘパラメータは意識のある患者でテストされた。そして、本発明による有利な態様ではこれらのパラメータは生物自身の呼吸頻度から独立していることを決定することが可能であった。
【００７９】
二つのＡＵＣｍｏｙ及びＡＵＣｍａｘパラメータは、呼吸頻度が制御された換気装置によって支配されている全身麻酔下にある意識のない患者の痛みの評価のためにも成功裏にテストされた。全身麻酔中のこれらのパラメータを知ることは全身麻酔中の鎮痛成分の間接的な測定を可能にする。
【００８０】
しかし、本発明は、二つの特異的な上述のＡＵＣｍｏｙ及びＡＵＣｍａｘパラメータに限定されない。これらのパラメータは好ましい例として与えられたにすぎない。
【００８１】
より一般的には、目標とされた結果（痛みの評価）に影響を与える本発明の本質的要素は、最終パラメータがメインウインドゥのＲＲｉサンプルの組からそれらの実体において直接計算される解決策（これは本発明の一部でない）とは対照的に、時間サブウインドゥに含まれるＲＲｉサンプルに対して計算された中間パラメータから最終パラメータを計算することであるということが見出された。従って、一般的に、本発明の方法は、時間サブウインドゥにおける中間パラメータＡｊの計算を可能にし、最終パラメータがこれらの中間パラメータの関数において痛み刺激の強度と相関する、ＲＲｉサンプルのいかなる処理方法を使用しても実行されることができる。
【００８２】
同様に、予期せぬことに、痛み刺激は呼吸ピーク（最小点Ｐｉ）に極めて顕著に影響を与えることが証明された；痛み刺激が強いほど、最小値（Ｐｉ）に相当するピークの振幅は弱くなり、いくつかの場合には実際にはゼロであることがありうる。対照的に、曲線の最大値（Ｍｉ）は痛み刺激によって顕著に影響を与えられない。結果として、これらのピーク（Ｐｉ）の振幅の変動性を考慮にいれることを可能にする、各サブウインドゥ（Ｆｊ）におけるＲＲｉサンプルのいかなる処理方法も各サブウインドゥにおける中間パラメータを計算するために好適であろう。従って、本発明は上述の面積計算に限定されず、上述の面積計算は一例として与えられたにすぎない。非限定的例として、面積（Ａｊ）の計算は、最小点（Ｐｉ）に対応するピークの振幅の関数における中間パラメータの計算によって置換されることができる。同様に、さらなる実施態様では、好適な形態のウェーブレット（例えばドゥーベチィーのウェーブレット）に基づいてサブウインドゥ（Ｆｊ）のＲＲｉサンプルにウェーブレット変換を適用することによって各中間パラメータ（Ａｊ）を計算することが可能である。
【００８３】
添付の図面を参照して記述された実施態様では、計算アルゴリズムは、二つの連続する心拍を分離する時間間隔（δｔｉ）を特徴とするサンプル（ＲＲｉ）のＲＲシリーズに対して適用される。別の実施態様では、本発明は、二つの連続する心拍の間の時間間隔（δｔｉ）の逆数（１／δｔｉ）を特徴とするサンプル（ＲＲｉ）のＲＲシリーズに基づいて実行されることができる。この場合、上述の段階１〜４を適用することによって、図５の曲線とは逆の曲線が得られ、点Ｐｉは最大値となり点Ｍｉは最小値となり、上述（段階６）の面請（Ａｊ）の計算は最大点（Ｐｉ）を連結する上方包絡線（Ｃｊ）に基づいて実行される。
【００８４】
上述の計算アルゴリズム並びにそこから得られるパラメータ、即ちＡＵＣｍｏｙ及びＡＵＣｍａｘパラメータは、生物の痛みの評価に対する用途の観点に限定されず、これらは一層一般的な態様で生物（意識のある又は意識のない）における心調律の変動性を測定するために使用されることができ、その結果としてＡＮＳ（自律神経系）の活動を変更する可能性のあるいかなる刺激のＡＮＳに対する影響を研究するために使用されることができる。
【００８５】
従って、痛みの評価に対して適用される上述の指示は、生物（意識のある又は意識のない）によって感じられているストレスのレベルを評価するために、他の適用分野及び特に準医療の分野にそのまま使用されて転用されることができる。
【図面の簡単な説明】
【００８６】
【図１】本発明による分析システムの一例の主要要素の概略を示す。
【図２】ＥＣＧ信号における心拍を特徴付ける一組の波（ＰＱＲＳＴ）を示す。
【図３】アナログＥＣＧ信号のサンプリング後に得られたディジタルＥＣＧ信号の一例を示す。
【図４】図３の信号から構築されたＲＲシリーズ（ＲＲ信号とも称する）の一例を示す。
【図５】ろ波及び正規化後のＲＲシリーズの一例を示す。
【図１】
2008513073.tif 000005
【図２】
2008513073.tif 000006
【図３】
2008513073.tif 000007
【図４】
2008513073.tif 000008
【図５】
2008513073.tif 000009
【国際調査報告】
2008513073.tif 000010
2008513073.tif 000011
2008513073.tif 000012
2008513073.tif 000013
2008513073.tif 000014
2008513073.tif 000015
2008513073.tif 000016
2008513073.tif 000017