(57)【特許請求の範囲】
【請求項１】
対象物の周期的運動を表す複数のサイクルを含む運動信号を得るステップと、
前記運動信号の振幅及び傾斜に基づき、前記運動信号に対して運動位相をマッピングするマップを生成するステップと、
を有し、前記マップを生成するステップが、運動信号の各サイクルにおいて下降振幅領域及び上昇振幅領域を特定し、前記複数のサイクル内のそれぞれの下降振幅領域に運動位相の一方のものをマッピングし、前記複数のサイクル内のそれぞれの上昇振幅領域に運動位相の他方のものをマッピングすることを含む、方法。
【請求項２】
前記マップに基づき、所定の関心位相に対する前記複数のサイクルに関するタイムスタンプのセットを特定するステップであって、特定されたタイムスタンプのセットが、共通の信号振幅に対応する、ステップと、
前記特定されたタイムスタンプのセットを示す信号を生成するステップとを更に有する、請求項１に記載の方法。
【請求項３】
反復的な周期的運動と前記反復的な周期的運動に対応する複数のサイクルとを持つ移動する対象物を撮像する間、前記運動信号が検出される、請求項１又は２に記載の方法。
【請求項４】
前記対象物が、被検者であり、前記運動信号は、前記被検者の複数の呼吸又は心周期を示す、請求項１乃至３のいずれかに記載の方法。
【請求項５】
前記マップを生成するステップが、前記下降振幅領域に対して１つのセットの運動位相をマッピングし、前記上昇振幅領域に対して別のセットの運動位相をマッピングするステップを含む、請求項１乃至４のいずれかに記載の方法。
【請求項６】
前記下降振幅領域が、前記対象物の収縮又は拡張のいずれかに対応し、前記上昇振幅領域は、前記対象物の拡張又は収縮の他方に対応する、請求項１乃至５のいずれかに記載の方法。
【請求項７】
前記運動位相が、前記サイクルの開始位相である０％から前記サイクルの終了位相である９９％まで広がり、前記運動位相の１つのセットは、０％から４９％の範囲にわたり、前記運動位相の別のセットは、５０％から９９％の範囲にわたる、請求項５又は６に記載の方法。
【請求項８】
前記複数のサイクルに関する信号最大及び信号最小を特定するステップと、
前記特定された信号最大及び信号最小に基づき、それぞれ平均信号最大及び平均信号最小を決定するステップと、
前記平均信号最大及び前記サイクルに関する信号最大の小さい方から、前記平均信号最小及び前記サイクルに関する信号最小の大きい方まで広がる範囲において、各サイクルの振幅へと前記運動位相をマッピングするステップとを更に有する、請求項１乃至７のいずれかに記載の方法。
【請求項９】
対象物の周期的運動を表す複数のサイクルを含む運動信号を入力し、各サイクルにおける下降振幅領域及び上昇振幅領域を特定する信号傾斜特定部と、
前記運動信号の振幅並びに下降振幅領域及び上昇振幅領域に基づき、前記運動信号に対して運動位相をマッピングする位相対振幅／傾斜マッパーを有し、
前記マッピングが、前記複数のサイクル内のそれぞれの下降振幅領域に運動位相の一方のものをマッピングし、前記複数のサイクル内のそれぞれの上昇振幅領域に運動位相の他方のものをマッピングすることを含む、システム。
【請求項１０】
所定の関心運動位相及び位相対振幅／傾斜マッピングに基づき、運動信号における複数のサイクルに関して、運動信号タイムスタンプのセットを特定するタイムスタンプ特定部を更に有し、前記運動信号タイムスタンプのセットが、共通の信号振幅に対応する、請求項９に記載のシステム。
【請求項１１】
前記運動信号の１つ又は複数のサイクルに関する信号最大及び信号最小を特定する信号特性特定部と、
前記特定された最大及び最小に基づき、平均最大及び平均最小を決定する信号特性プロセッサとを更に有する、請求項９に記載のシステム。
【請求項１２】
前記信号特性特定部が、最大決定閾値を下回る振幅を持つ２つのポイントの間にあるポイントとして、運動サイクルに含まれる評価範囲に関する最大を特定し、
前記最大の前後の前記評価範囲における他のポイントと前記最大とを比較し、及び
前記最大が前記ポイントより大きいとき、前記評価範囲に関する局所最大として前記最大を特定する、請求項１１に記載のシステム。
【請求項１３】
前記位相対振幅／傾斜マッパーが、前記平均最大及び平均最小並びに前記下降及び上昇振幅領域に基づき、前記運動信号に対して運動位相をマッピングする、請求項１１又は１２に記載のシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
本願は一般に、位相対振幅／傾斜マッピングに関し、本書では特に、不規則な周期を持つ周期的信号を参照して説明される。他の用途も、本書において想定される。
【背景技術】
【０００２】
撮像において、患者が息を止めるのではなく自由に呼吸する手順のため、特定の関心呼吸位相（例えば、完全な吸気、完全な呼気、又はこれらの間のどこか）に対応する投影データを特定するのに、患者をスキャンする間に同時に得られる呼吸信号の呼吸位相ベースのゲーティングが用いられる。この手法を用いると、呼吸信号は、（信号最大及び最小に基づき）個別の呼吸サイクルに分けられ、各サイクルは、サイクルの開始（０％）（例えば完全な吸気又は呼気）からサイクルの終了（９９％）までの範囲の位相へと分けられる。
【０００３】
図１は、上記の例を示す。図１において、ｙ軸は、信号振幅（即ち拡張及び収縮を含む肺の機械的な状態）を表し、ｘ軸は、時間を表す。この例に関して、呼吸サイクル１０２は、完全な吸気又は肺拡張で始まるよう規定される。これは、サイクル１０２の信号最大１０４に対応する。この例において、０％の位相は、最大１０４に割り当てられる。サイクル１０２の残りの部分は、例えば、最大１０４に対して、サイクル１０２の２５％、５０％、７５％及び９９％といった複数の位相へと線形に分けられる。
【０００４】
再構成のため、特定の関心位相に関するタイムスタンプのセットが、サイクルのパーセンテージに基づき特定され、再構成器に搬送される。再構成器は、このタイムスタンプを用いて、再構成のための投影データを特定する。例えば、図１において、７５％の位相が関心位相である場合、サイクル１０２に関するタイムスタンプ１１２が、特定され、再構成器に搬送される。図１における他のサイクルに関して、タイムスタンプ１１２、１１４、１１６、１１８及び１２０は、７５％の位相に対応するものとして特定され、再構成器に搬送される。
【発明の概要】
【発明が解決しようとする課題】
【０００５】
残念なことに、この仮定を用いる場合、複数の問題が存在する。（ａ）呼気のレートは一般に、吸気のレートより低い。そこで、この位相が、サイクルにわたり線形に広がらない。（ｂ）呼吸サイクル形状は、（図１で分かるように）サイクル毎に異なる場合がある。（ｃ）吸息の間の肺の機械的な状態（振幅）が、呼息の間の同じであると保証されるわけではない。例えば、図１において、ポイント１１２、１１６及び１２０に対応するタイムスタンプは、肺が収縮し、振幅が最小にある完全な呼気に近いか、又は完全な呼気にある。一方、タイムスタンプ１１４及び１１８は、肺が部分的に拡張され、より大きな振幅を持つ吸気の間の機械的な状態に対応する。結果として、再構成された画像データにおける生体構造が整列することができない。これは、再構成された画像データに運動アーチファクトをもたらす場合がある。
【０００６】
本願の側面は、上述した事項その他に対処する。
【課題を解決するための手段】
【０００７】
１つの側面によれば、ある方法が、複数のサイクルを含む信号を得るステップと、上記信号の振幅及び傾斜に基づき、上記信号に対して運動位相をマッピングするマップを生成するステップとを含む。
【０００８】
別の側面によれば、あるシステムが、所定の関心運動位相及び位相対振幅／傾斜マッピングに基づき、移動する対象物の周期的運動を示す運動信号における複数の運動サイクルに関して、運動信号タイムスタンプのセットを特定するプロセッサを含み、上記運動信号タイムスタンプのセットが、共通の信号振幅に対応する。
【０００９】
別の側面によれば、ある方法が、不規則な周期性を持つノイズが多い周期的信号の運動サイクルにおける複数のピークから一ピークを特定するステップであって、上記ピークが、所定の閾値を下回る振幅を持つ２つのポイントの間にあるポイントに対応する、ステップと、上記ピークの前後のポイントと上記ピークとを比較するステップと、上記ピークが上記周囲のポイントより大きいとき、局所最大として上記ピークを特定するステップとを含む。
【図面の簡単な説明】
【００１０】
【図１】従来技術を示す図である。
【図２】運動位相相関器に関連して例示的な撮像システムを示す図である。
【図３】例示的な運動位相相関器を示す図である。
【図４】運動信号におけるサイクルの振幅及び傾斜に対して運動位相をマッピングする例示的な手法を示す図である。
【図５】運動位相相関器を使用することなく生成された例示的な画像を示す図である。
【図６】運動位相相関器により選択される投影データを用いて生成される例示的な画像を示す図である。
【図７】運動信号における最大を特定する例示的な手法を示す図である。
【図８】サイクルにおける複数の最大の間の局所最大を特定する例示的な手法を示す図である。
【図９】サイクルに関する複数の最大候補の間における最大を拒絶する例示的な手法を示す図である。
【図１０】運動信号を増やす例示的な手法を示す図である。
【図１１】再構成のため関心運動位相に関する投影データを特定する例示的な方法を示す図である。
【図１２】信号における局所最大のセットを決定する例示的な方法を示す図である。
【発明を実施するための形態】
【００１１】
本発明は、様々な要素及び要素の配列の形式並びに様々なステップ及びステップの配列の形式を取ることができる。図面は、好ましい実施形態を説明するためだけにあり、本発明を限定するものとして解釈されるべきものではない。
【００１２】
図２は、例えばコンピュータ断層撮影（ＣＴ）スキャナといった撮像システム２００を示す。しかしながら、他の実施形態では、撮像システム２００は、例えばＸ線、陽電子放出断層撮影（ＰＥＴ）、単光子放出断層撮影（ＳＰＥＣＴ）、磁気共鳴イメージング（ＭＲＩ）、超音波（ＵＳ）及び／又は他の撮像モダリティといった１つ又は複数の他の撮像モダリティを含むことができる。
【００１３】
撮像システム２００は、静止したガントリー２０２及び回転ガントリー２０４を含む。これは、静止したガントリー２０２により支持される。回転ガントリー２０４は、長手方向又はｚ軸の周りで検査領域２０６の周りを回転する。
【００１４】
例えばＸ線管といった放射線源２０８が、回転ガントリー２０４により支持され、回転ガントリー２０４と共に回転し、放射線を放出する。ソース２０８の反対側に、検査領域２０６を横切るよう配置される放射線感知検出器アレイ２１０が、検査領域２０６を横断する放射線を検出し、検出された放射線に対応する投影データを生成する。
【００１５】
例えば寝台といった支持部２１２が、検査領域２０６において対象物又は被検者を支持する。運動モニタ２１４は、スキャンの間、移動する対象物の運動状態をモニタし、運動を示す運動信号を生成する。移動する対象物の例は、以下に限定されるものではないが、患者の肺若しくは心臓、規則的な又は不規則な周期的運動を繰り返す他の対象物を含む。
【００１６】
運動位相相関器２１６は、撮像又はスキャン取得タイムスタンプと運動位相とを相関させる。後述されるように、１つの例において、これは、運動信号におけるサイクルの振幅及び傾斜の両方に基づき、運動信号と運動位相とを相関させることにより実現される。結果として生じる相関は、異なるサイクルに対して、同じ運動位相に対応し、共通の又は同じ機械的な振幅又は状態を持つ投影データのタイムスタンプを特定するために用いられることができる。これは、異なるサイクルにわたる構造の不整列が原因による運動アーチファクトの緩和を容易にすることができる。
【００１７】
再構成器２１８は、投影データを再構成し、検査領域２０６を示すボリュメトリック画像データを生成する。再構成器２１８は、関心位相に関する運動位相相関器２１６を介して特定されるタイムスタンプに基づき、投影データを選択的に再構成するよう構成される。再構成器２１８は、様々な再構成アルゴリズムを使用することができ、単一の又は複数の取得からの投影データが、単一の画像を生成するために再構成されることができる。
【００１８】
汎用のコンピュータシステムが、オペレータ端末２２０として機能する。これは、例えばディスプレイ及び／又はプリンタといった人間が読み取れる出力デバイスと、例えばキーボード及び／又はマウスといった入力デバイスとを含む。端末２２０上のソフトウェアは、オペレータが、システム２００の動作を制御することを可能にする。この制御は、例えば、オペレータが、スキャンプロトコルを選択し、スキャンを開始し、再構成アルゴリズム（例えば、運動ゲート再構成）を選択すること等を可能にすることにより、行われる。
【００１９】
運動位相相関器２１６が、システム１００の一部、例えば、端末２２０等の一部とすることができ、又はシステム１００のリモートにある、例えば、ネットワークにわたり及び／又は他の態様で分散されるワークステーションの一部とすることが出来る点を理解されたい。更に、運動位相相関器２１６は、メモリといったコンピュータ可読媒体に埋め込まれる又はエンコードされるコンピュータ可読命令を処理する１つ又は複数のプロセッサにより実現されることができる。
【００２０】
図３は、例示的な運動位相相関器２１６を示す。
【００２１】
信号特性特定部３０２は、運動信号の１つ又は複数の特性を特定する。これは、運動信号のサイクルの１つ又は複数に関する１つ又は複数の特性を特定することを含むことができる。斯かる特性の例は、最大及び最小並びに／又は他のデータポイントを含むが、これらに限定されるものではない。
【００２２】
信号特性プロセッサ３０４は、特定された信号特性を処理する。これは、特定された信号特性に関する１つ又は複数の統計量（例えば、平均、メディアン、最大、最小等）及び／又は他の情報を含むことができる。例えば、統計量は、１つ又は複数のサイクルに関する信号の振幅に関する。
【００２３】
信号傾斜特定部３０６は、運動信号におけるサイクルにおいて下降する及び上昇する領域を特定する。
【００２４】
位相対振幅／傾斜マッパー３０８は、処理済の特定された信号特性及び特定された傾斜に基づき、運動信号に対して運動位相をマッピングする。
【００２５】
タイムスタンプ特定部３１０は、運動信号におけるタイムスタンプを特定する。この運動信号は、関心運動位相に対して対応し、マッピングに基づき類似する信号振幅を持つ。
【００２６】
運動位相相関器２１６はここで、対象物の撮像と同時に得られる呼吸信号といった不規則な周期性を持つ周期的信号を処理することに関連して説明される。他の適切な信号は、不規則な周期性を持つノイズが多い周期的信号を含む他の周期的信号を含むが、これに限定されるものではない。
【００２７】
図４を参照すると、信号４００は、変化する振幅及び時間幅を持つ複数の呼吸（吸気及び呼気）サイクル４０２１−４０２６を含む。図４において、ｙ軸は、信号振幅（機械的な振幅）を表し、ｘ軸は、時間を表す。
【００２８】
信号特性特定部３０２は、呼吸信号４００の呼吸サイクル４０２１−４０２６の最大４０４１−４０４６及び最小４０６１−４０６６を特定する。この例において、特定のサイクルに関して、最大（又は、最大振幅）は、肺が拡張又は膨張される完全な吸気に対応し、最小（又は、最小振幅）は、肺が収縮又はつぶされる完全な呼気に対応する。
【００２９】
信号特性プロセッサ３０４は、特定された最大４０４１−４０４６及び最小４０６１−４０６６に関してそれぞれ平均最大４０８及び平均最小４１０を決定する。
【００３０】
傾斜特定部３０６は、最大４０４１−４０４６及び最小４０６１−４０６６に基づき、各サイクルに関して、下降する（呼気又は吸気）及び上昇する（吸気又は呼気）領域又は傾斜を特定する。明確さのため、サイクル４０２２に関する傾斜４１２及び４１４だけが、図４において参照符号を介して特定される。図４において、呼気のレート（傾斜４１２）は、吸気のレート（傾斜４１４）よりゆっくりしている点に留意されたい。
【００３１】
位相対振幅／傾斜マッパー３０８は、平均最大４０８及び平均最小４１０の間で、呼気及び吸気傾斜４１２及び４１４に対してそれぞれ運動位相をマッピングする。例えば、サイクル４０２２は、最大４０４２から次の最大４０４３へと広がり、０％の位相は、平均最大４０８へとマッピングされ、４９％の位相は、平均最小４１０へとマッピングされる。
【００３２】
サイクル４０２２の０％及び４９％の位相の間で発生する呼気は、呼気傾斜４１２へとマッピングされ、サイクル４０２２の４９％及び９９％の位相の間で発生する吸気は、吸気傾斜４１４へとマッピングされる。そのようなものとして、呼気及び吸気に関する位相は別々に決定される。これは、呼気の平均レートが吸気の平均レートより一般に低いことを考慮する。そのようなものとして、各サイクルは、０％から９９％の間の位相に対応するデータポイントを持つことが保証される。
【００３３】
いくつかのサイクルは、特定の位相に対応する１つ以上のポイントを持つことができる。これらの場合、第１の又は他の候補ポイントが選択される。サイクル（例えば、サイクル４０２１）が平均最大（例えば、最大４０４１）未満の最大を含む場合、サイクル最大は０％とみなされる。従って、図４において、位相マッピングは、最大４０４１から最小４０６１までの広がりに対して適切にスケール化される。
【００３４】
タイムスタンプ特定部３１０は、運動位相に関する撮像又はスキャン取得タイムスタンプを特定する。例えば、図４において、タイムスタンプ特定部３１０は、運動サイクル４０２１−４０２６の２５％の位相に対応するタイムスタンプ４１６１−４１６６を特定するために用いられる。もちろん、２５％以外の位相に対応するタイムスタンプが特定されることもできる。タイムスタンプ特定部３１０は、タイムスタンプを示す信号を生成し、再構成器２１８（図２）に対して信号を搬送する。再構成器は、この信号に基づき、再構成のための投影データを選択することができる。
【００３５】
これらのタイムスタンプが、移動する対象物の実質的に類似する機械的な振幅に相関するので、再構成された画像データにおける移動する対象物の構造は、うまく整列配置される可能性が高い。これは、図５及び６において見られることができる。
【００３６】
図５は、斯かるタイムスタンプが、再構成のための投影データを選択するのに用いられない場合に再構成される画像５００を示す。投影データにおける生体構造は、整列配置されないか、又は誤って整列配置されており、これは、画像において矢印により示されるように、画像５００にアーチファクト５０２をもたらす。図５におけるアーチファクトは、説明的な目的のために提供され、他の又は異なるアーチファクトが、画像に存在する場合がある。
【００３７】
図６は、特定の運動位相に関して運動位相相関器２１６により特定されるタイムスタンプに対応する投影データに基づき再構成される画像６００を示す。本書において、投影データにおける生体構造は整列配置され、画像において矢印により示されるように、画像５００におけるアーチファクトは存在しない。
【００３８】
上述されたように、信号特性特定部３０２は、運動信号の１つ又は複数の特性、例えば運動信号におけるサイクルに関する最大及び最小を特定する。運動信号がノイズを含むことができる点を理解されたい。これは、信号における偽の最大及び／又は最小、及び／又は不明瞭な真の最大及び／又は最小を提供する場合がある。図７〜９は、ノイズがある場合における最大及び最小を決定する手法を示す。斯かる信号の例は、呼吸、心臓及び／又は他の信号を含む。
【００３９】
最初に、図７を参照して、ノイズが多い信号７０２のサブ部分が示される。ピーク又は最大は、所定の評価範囲に対して、信号７０２において特定される。図示された例において、ピーク７０４は２つのポイント７０６及び７０８の間の単一の最高ポイントとして特定される。これは、所定の検出閾値７１０より小さい。この例において、閾値７１０は、信号のグローバルダイナミックレンジの１／１０、又はグローバル最大及びグローバル最小の間の差であるよう規定される。他の実施形態では、閾値７１０は、他の態様で規定されることができる。
【００４０】
その後、評価範囲における信号７０２の各データポイントが評価される。データポイントが第２の閾値７１２の下に落ちるまで、データポイント７０４の前の評価範囲における信号７０２のデータポイントが検索される。これは、閾値７１０未満の最小から最大への振幅に対応する。同様に、データポイントが第２の閾値７１２の下に落ちるまで、データポイント７０４の後の評価範囲における信号７０２のデータポイントが、検索される。評価範囲における全てのデータポイントが、データポイント７０４に対して比較される。
【００４１】
ピーク７０４が、評価範囲にわたり、最大である場合、それは局所最大として特定される。そうでない場合、検索は、次のポイントへと続く。これは、図８及び９において更に図示される。図８において、閾値８０６に基づき評価範囲８０４において特定されるピーク８０２が、局所最大として特定される。なぜなら、これが、範囲８０４において最も高いポイントだからである。図９において、評価範囲９０４において特定されるピーク９０２は、閾値９０６に基づき、局所最大としては特定されない。なぜなら、これは、範囲９０４において最も高いポイントでないからである。
【００４２】
図１０は、信号を増やす非限定的な技術を示す。
【００４３】
この技術を用いると、運動信号１０００の２つのピーク１００４及び１００６の間の時間距離１００２が、追加のピークの閾値を越える場合、１つ又は複数のピーク１００８及び１０１０が信号１０００に加えられる。ピーク１０１４及び前のピーク１０１６の間の時間距離１０１２並びに最大１０１４及び後続のピーク１０２０の間の時間距離１０１８が共に追加の閾値未満である場合、ピーク１０１４は信号１０００から除去される。閾値は、平均ピークの間隔又はその他に基づかれることができる。
【００４４】
図１１は、再構成のため関心運動位相に関する投影データを特定する方法を示す。
【００４５】
ステップ１１０２において、反復的な周期的運動を持つ移動する対象物をスキャンする間、運動信号が同時に得られる。運動は、呼吸、心拍又は他の周期的運動に関連付けられることができる。
【００４６】
ステップ１１０４において、最大及び最小振幅が、運動信号において特定される。
【００４７】
ステップ１１０６において、グローバル平均最大振幅及びグローバル平均最小振幅が、特定された最大及び最小に基づき決定される。
【００４８】
ステップ１１０８において、グローバル平均最大振幅は、サイクル最大からサイクル最小へと延在する各サイクルの第１の領域に関する開始又は０％のポイントを規定するために用いられる。グローバル最大が、サイクル最大サイクルより大きい場合、開始又は０％のポイントは、サイクル最大で規定される。グローバル最大が、サイクル最大より小さい場合、開始又は０％のポイントは、サイクル最大に続くポイントで規定される。この振幅は、グローバル最大に最も整合する。
【００４９】
ステップ１１１０において、グローバル平均最小振幅は、サイクル最小から次のサイクルの最大へと延在する各サイクルの第２の領域に関する５０％のポイントを規定するために用いられる。グローバル平均最小が、サイクル最小未満の場合、５０％のポイントはサイクル最小で規定される。グローバル最小が、サイクル最小より大きい場合、５０％のポイントは、サイクル最小に続くポイントで規定される。この振幅は、グローバル最小に最も整合する。
【００５０】
ステップ１１１２において、０％から４９％までの運動位相が、第１の領域に対してマッピングされる。
【００５１】
ステップ１１１４において、５０％から９９％までの運動位相が、第２の領域に対してマッピングされる。
【００５２】
ステップ１１１６において、関心位相が特定される。
【００５３】
ステップ１１１８において、特定された位相及び振幅／傾斜の間のマッピングが、マッピングされた振幅に対応するサイクル毎のスキャン取得タイムスタンプを決定するために用いられる。
【００５４】
ステップ１１２０において、特定されたタイムスタンプで取得される投影データが、画像ボリュームへと再構成される。
【００５５】
図１２は、信号最大を決定する方法を示す。
【００５６】
ステップ１２０２において、過去のサイクル最大までの任意の過去のポイントの振幅より振幅が所定の閾値分大きい第１のポイントとして、ピークが決定される。
【００５７】
ステップ１２０４において、ピークに続くポイントがピークを上回るかどうかが決定される。
【００５８】
上回る場合、ステップ１２０２は、次のポイントを用いて繰り返される。
【００５９】
上回らない場合、ステップ１２０６において、ピークに続くポイントが、所定の閾値分ピークより少ないかどうかが決定される。
【００６０】
少ない場合、ステップ１２０８において、このピークは、サイクル最大として特定される。
【００６１】
そうでない場合、ステップ１２０４が、次のポイントを用いて繰り返される。
【００６２】
ステップ１２１０において、最大が、所定の閾値より近いかどうかが決定される。
【００６３】
近い場合、ステップ１２１２において、これらの最大の１つ又は複数が、最大としての特定を外される（除去される）。１つの非限定的な実施形態を用いると、最大間の間隔の平均及び標準偏差が計算され、平均−ｘ＊標準偏差（ここで、ｘは入力パラメータである）より近い最大の１つ又は複数が、最大としての特定を外される。
【００６４】
ステップ１２１４において、最大が、所定の閾値よりはるかに離れているかどうかが決定される。
【００６５】
離れている場合、ステップ１２１６において、これらの最大の１つ又は複数が加えられる。１つの非限定的な実施形態を用いると、最大間の間隔の平均及び標準偏差が計算され、１つ又は複数の最大が、平均＋ｘ＊標準偏差（ここで、ｘは入力パラメータである）よりよりはるか離れている場合、ステップ１２０２から１２０８を用いて、１つ又は複数の最大が、境界ピークの間に加えられる。ここで、適切な最大が特定されるまで、分類閾値が反復的に減弱される。
【００６６】
上記ステップは、コンピュータ可読命令を用いて実現されることができる。これは、コンピュータプロセッサにより実行されるとき、プロセッサに本書に記載されるステップを実行させる。斯かる場合、この命令は、関連コンピュータに関連付けられる及び／又は他の態様でアクセス可能である例えばメモリのようなコンピュータ可読ストレージ媒体に格納される。
【００６７】
上記説明は、ＣＴ及び呼吸信号の文脈において述べられたが、本書において開示される手法は、不規則な周期性を持つ別のノイズが多い周期的信号における実質的に同等のデータポイントを特定するために使用されることができる点を理解されたい。
【００６８】
本発明が、好ましい実施形態を参照して説明されてきた。上記の詳細な説明を読み及び理解すると、第三者は、修正及び変更を思いつくことができる。それらの修正及び変更が添付の特許請求の範囲又はその均等物の範囲内にある限り、本発明は、すべての斯かる修正及び変更を含むものとして構築されることが意図される。
【図１】
0005878920.tif 000002
【図２】
0005878920.tif 000003
【図３】
0005878920.tif 000004
【図４】
0005878920.tif 000005
【図５】
0005878920.tif 000006
【図６】
0005878920.tif 000007
【図７】
0005878920.tif 000008
【図８】
0005878920.tif 000009
【図９】
0005878920.tif 000010
【図１０】
0005878920.tif 000011
【図１１】
0005878920.tif 000012
【図１２】
0005878920.tif 000013