JP2013523280A - 放射状ｍｒｉデータに対するオシログラム・フィルタリング方法 - Google Patents

Abstract

システムが、データ・レシーバ、サイノグラム・ジェネレータ、プロセッサ、フィルタ・モジュール、出力モジュールを有する。データ・レシーバは、放射状に規則化された磁気共鳴データを受け取るように構成される。サイノグラム・ジェネレータは、磁気共鳴データに対する読出し方向の関数としての視野角に対応する第１のサイノグラムを作り出すように構成される。プロセッサは、角周波数軸を有するオシログラムを作り出すように構成される。オシログラムは、第１のサイノグラムをフーリエ変換したものに対応する。フィルタ・モジュールは、オシログラムの選択された軸に沿って形成された投影内のピークを選択的にフィルタリングするように構成されるが、このピークは、ＲＦ干渉等の干渉信号に関連するものである。この選択された軸は、角周波数軸に直交する。出力モジュールは、フィルタリングを行った投影を変換したものに対応する第２のサイノグラムを形成するように構成される。
【選択図】図１

Description

優先権の主張
本特許出願は、米国特許法第１１９条（ｅ）に基づいて、出願日２０１０年３月３１日、Curtis A. Corumの米国特許仮出願第６１／３１９，７０３号、名称「METHOD OF OSCILLOGRAM FILTERING FOR RADIAL MRI DATA」（代理人整理番号、６００．７４６ＰＲＶ）に対する優先権の利益を主張し、また、同仮出願全体を参照により本明細書に組み込む。

米国政府支援による研究または開発についての申告
本発明は、米国国立衛生研究所、ＮＩＨ ＢＴＲＣ Ｃｏｒｅ３からの授与番号５Ｐ４１ＲＲ００８０７９に基づく政府支援を受けてなされた。米国政府は、本発明の一定の権利を有する。

磁気共鳴画像中のアーチファクトには、断続的で広帯域のノイズが含まれることがある。現在の技術は、こうしたアーチファクトを軽減するのには不十分である。

図面は、必ずしも原寸に比例して描かれておらず、異なる図の中の同様のコンポーネントを同様の数字で説明していることがある。異なる添え字を持つ同様の数字が、同様のコンポーネントの異なる例を表すことがある。図面は、本文書で議論する様々な実施形態を限定ではなく例として概ね例示している。
ソノグラムおよび対応するオシログラムの図である。 オシログラムおよび対応する投影の図である。 オシログラムおよび対応する投影の図である。 １つの例による方法の図である。 １つの例によるシステムの図である。

以下の例は、任意の順序または組合せで組み合わせることができる。

本発明による主題には、放射状ＭＲＩデータ（radial MRI data）から断続的で広帯域のノイズをフィルタリングする方法が含まれる。オシログラム・ドメイン中に示す際の放射状ＭＲＩデータのまばらさ、凸性により、画像データへの影響を最小限に抑えながら、フィルタリングを行って、ノイズ（または他の人為現象的信号）を分離および抑制することが可能になる。

ＲＦ干渉およびジッパーアーチファクト
無線周波数（ＲＦ）アーチファクトが、デカルト・サンプルＭＲＩ画像中に「ジッパー（zipper）」として頻繁に現れる。結果として得られる画像内で、ジッパーアーチファクトは一般に、読出し方向には狭くなり、位相符号化方向に沿って延びる。

放射状（投影（プロジェクション））イメージングにおいては、ジッパーは、最終画像中の特定の線（または面）に局在せず、中間のサイノグラム内に現れる。図１は、視野角について一定であるアーチファクトＡおよび視野角について正弦的に現れたアーチファクトＢを有するサイノグラム（スパイラログラム（spiralogram）またはスパイログラムと称することもある）を示す。

２次元または３次元放射状データセット内には、ジッパーは、視野順序(view ordering)ならびにジッパー干渉の周波数、帯域、位相に依存した構造である「ファン（fan）」アーチファクトまたは「タービン」アーチファクトとして現れることがある。

概ね垂直の（または変動しない）ジッパーアーチファクトに加えて、ジッパーが移動中であること、または時間的に変動する周波数オフセットを有することもある。このことは、放射状取得の視界がアイソセンタに対してオフセットである場合に起こる。この場合、ＲＦ干渉は、実際には固定周波数であってもよいが、ＭＲＩデータは、非アイソセンタ取得の周波数オフセット・テーブルに起因してオフセットである。

移動ジッパーは、一般に、アイソセンタからの距離および視野順序に依存した正弦的経路を描くこととなる。図示の例では、ジッパーのベースは、ｒ-インデックスＢで始まる。

他のフィルタ法
メジアン・フィルタまたは他のフィルタを、（サーチ・アルゴリズムが時間を浪費することがあるが）垂直ジッパーまたは移動ジッパーに適用することもできる。

オシログラム
オシログラムは、サイノグラムの図を変換したものを表す。フーリエ変換が、視野方向に沿って、また、サイノグラムの読出し（ｒ-インデックス）方向に沿って行われる。サイノグラム表現では、ジッパーは、移動するまたは変動しない概ね垂直のフィーチャとして現れ、オシログラム表現では、ジッパーは、水平線として現れる。

本発明の主題の一例によれば、オシログラムの振幅をフィルタリングに使用することができる。

オシログラム・フィルタ
本発明の主題の一例は、オシログラム内のＲＦジッパーアーチファクトに対応する水平線を自動的に認識し、排除するように構成される。１つの方法は、水平軸全体にわたって積分し、１次元投影を作成することを伴う。

次いで、オシログラムの水平線は、投影上で背景の上方に延びるピークに変換される。

実施例
差分画像を作成して、タービン・アーチファクトの除去および所望の画像データの維持を示すことができる。

最終画像内では、タービン・アーチファクトが、復元プロセスによりタービンと同様の放射状の対称形パターンに散らばるジッパーとして現れる。この構造は、ジッパー干渉の特性ならびに放射状の視野順序に依存する。

他の応用例
周期的またはほぼ周期的な動きもオシログラムから検出可能である。これには、拍動流、心臓の動き、呼吸運動が含まれる。オシログラムまたはオシログラムの部分集合から、ナビゲータ型の情報を作成することもできる。

図１は、ＳＷＩＦＴシーケンスを使用することにより作成したサイノグラム内のジッパーアーチファクトを示す。

ＳＷＩＦＴは、空間的情報を符号化するのに使用される磁場勾配の新たな値（または向き）を少なくとも設定するのに必要な時間だけパルス長Ｔ_Ｐを上回る短期の繰返し時間Ｔ_Ｒを有する、周波数変調されたパルスのシーケンスを伴う。画像は、３次元バックプロジェクション・リコンストラクションを使用することにより処理する。一例では、双曲正割ファミリー（hyperbolic secant family）からの周波数変調されたパルス（ＨＳｎパルス）を使用する。ＳＷＩＦＴでは、形成されたパルスが、時間依存の振幅および位相を有するＮ個の異なるサブパルス要素を含むことがある。ＦＭパルス中は、瞬時の共鳴に至るまで、等色線は、有効ＲＦフィールド・ベクトルに追従する。共鳴時には、等色線は、ＲＦパルスの「ハグ（hug）」から解放され、次いで、衰退していく小さな変調を伴ったほぼ自由な歳差運動を示し、これによりスペクトル・コンタミネーションがもたらされる。したがって、こうしたスピン系反応からスペクトル情報を抽出するために、ストカスティック（確率論的）ＮＭＲにおいて位相情報を回復する方法と同様の相互相関手法を使用する処理を行う。ＴＲ＜＜Ｔ１のＳＷＩＦＴについて単位時間毎に理論的に達成可能な信号対雑音比（ＳＮＲ）は、パルスＦＴ（pulsed FT）についてのものと同一である。ＳＷＩＦＴ取得中、適用される画像化勾配は、磁化率または不均質性に起因して、本来備わる固有の勾配の全てを通常は超えるものである。この条件の場合、取得される画像は、横緩和に全く影響されず、信号の強さは、Ｔ１およびスピン密度にのみ依存する。最大Ｔ_１コントラストは、有効フリップ角に依存し、感度とコントラストとの最良の妥協点は、エルンスト角の２倍を越えるフリップ角を有することとなる。フリップ角が非常に小さい場合、Ｔ_１コントラストは、さほど重要でなく、コントラストは、スピン濃度に完全に由来する。画像取得に先立つ、または、画像取得の途中に行われる適当な準備シーケンスによって、他の種類のコントラストを得ることもできる。

ＳＷＩＦＴは、放射状イメージング技法の一例にすぎず、他の放射状イメージング技術を、本発明の主題とともに使用することもできる。

図１では、ｒ-インデックスＡにおいて、高強度のジッパーが、垂直になり、中心の左側に位置する。また、第２のジッパーをｒ-インデックスＢに示す。

図１の左側区画に示したサイノグラム・データは、垂直軸上に視野角を示し、水平軸上に読出し方向を示す。視野角データは順次作り出されるので、視野角を時間軸とみなすこともできる。サイノグラム・データは、２次元または３次元の放射状イメージング・プロトコルで作り出すことができる。

図１に示すサイノグラムは、ＳＷＩＦＴイメージング・プロトコルによって作られる自由誘導減衰（ＦＩＤ）として見ることができる。正弦波帯域は、検査中の対象のマージンを表す。

図示のジッパーアーチファクトは、様々な周波源を有することがある。たとえば、共鳴周波数における干渉が、こうしたアーチファクトを引き起こすことがある。さらに、エイリアシングされた周波数またはその周波数の調波からの干渉が、ジッパーアーチファクトにつながることもある。さらに、ジッパーは、システム内のＲＦ源またはそれに近い周波源により引き起こされることもある。ジッパーアーチファクトは、オフセット・テーブルに関連することもある。

本発明の主題は、フィルタを使用することによりサイノグラムを調整する方法およびシステムを含む。フィルタには、メジアン・フィルタ、ノッチ・フィルタ、または、アーチファクトを除去するようになされた他の種類のフィルタが含まれることがある。メジアン・フィルタは、１次元または複数次元で機能することができ、入力の中央値の出力表示を提供する。メジアン・フィルタは、歪みを低減し、平滑化機能を提供することができる。

ＳＷＩＦＴプロトコルを使用するイメージングは、相対的に高速の切換えおよび１〜４μ秒のオーダの短い持続期間を伴う。こうした短い持続期間により、本明細書で記載するジッパーアーチファクトを作り出すことのある無線周波数干渉源がもたらされることがある。

図１に示すように、サイノグラムをフーリエ変換して、オシログラムを作成する。図１のオシログラムでは、垂直軸（ここでは角周波数を指す）を遅軸（時間）とみなすことができるか、または、これを任意値ｎでインデックス付けすることができる。オシログラムの水平軸も時間軸であり、速軸を指し、ここでは解放前部分と解放後部分とに分けられている。オシログラムのデータは、視野ごとの順序に生のｋ空間データを表す。

図１に示すように、サイノグラム・ドメイン中のＡにおける干渉が、図のオシログラム中のＡ_１’およびＡ_２’に変換する。さらに、サイノグラム・ドメイン中のＢにおける干渉が、オシログラム中のＢ_１’およびＢ_２’に変換する。ＲＦ干渉Ａは、狭帯域もしくは広帯域になること、または、一定周波数を有することがある。ＲＦ干渉Ｂは、移動（moving）周波数に対応する。

画像情報（非移動）は、図に示すようにハッシュ・マークを対角線が横切る三角形領域に限定される。垂直のハッシュ・マークを有する三角形領域で示すように、周期的な動きにより、位置の変わったコンポーネントがもたらされることがある。画像情報の三角形の頂点と周期運動の三角形との垂直距離は、干渉の周波数に相当する。さらに、対象の幅が広いほど、オシログラム表現中の三角形の幅は広くなる。

対象データが視野順序とともに変化するので、本発明の主題により、対象の動きとアーチファクトとを識別することができる。干渉は、位相がコヒーレントであることもある。

サイノグラム内のデータを２つのフーリエ変換によって変換して、オシログラムを作成する。上で述べたように、垂直ジッパーアーチファクトは、時間において相対的に狭く、ｎ軸において、または周波数軸に沿って相対的に広い。

図２は、ｔ_ｆａｓｔ軸に沿ったオシログラム投影の１次元ノッチ・フィルタを示す。

投影Ｉ（Ａ１’）は、オシログラム投影内のアーチファクト・ピークまたは干渉ピークを示し、以下のように表すことができる。

この投影を振幅について積分し、一例では、虚数成分または位相成分について積分することができる。

アーチファクトは、画像情報を越えて延び、自動ピーク発見（検出）アルゴリズムによって識別することができる。図２では、オシログラム投影内の矩形ボックスは位置を示す。

オシログラム内のマスキング
ＲＦ干渉の角周波数座標を使用することにより、２次元オシログラム内にマスク-アウト・ストリップ（mask-out strip）を提供する。さらに、オシログラムの解放前部分（左側）を捨てることもできる。マスク-アウトすることは、フィルタリングを行うことまたは０に設定することを含む。

オシログラム投影内でのマスキングに続いて、変換および復元により画像サイノグラムを作成する。

一例では、オシログラム内の中央の水平線に「キープアウト」・ゾーンを設定することができる。このキープアウト・ゾーンは、ユーザ・インタフェースを使用することによりユーザが選択可能な幅を有することがある。一例では、キープアウト・ゾーンは、投影の中心の周囲約１０％または５〜２０％である。

様々なフィルタを使用することもできる。たとえば、こうしたフィルタを視界の外の画像データに適用する（帯域通過）ことができ、あるいは、こうしたフィルタには、ＲＦノッチ・フィルタが含まれることがある。

さらに、ＳＷＩＦＴを使用する脂肪抑制は、磁化プレパレーション（magnetization preparation）・パルスを含むことがある。磁化プレパレーション・プロトコルでは、複数のノッチ・フィルタを投影上で利用することがある。たとえば、複数のフィルタ、すなわち櫛形フィルタを本発明の主題とともに使用することができる。対象の範囲を超える線形領域を減衰させるために、マスクを使用することもできる。

周期信号に加えて、本発明の主題を使用することにより、（心血流や呼吸等の）動きを修正することもできる。

図３は、オシログラム・ドメイン内で構成され、先の図中に示したようなアーチファクトを除去するようになされたマスク・フィルタを有するオシログラム投影を示す。

図４は、一例による方法４００の流れ図を示す。この方法は、放射状に規則化された磁気共鳴データに対応する第１のサイノグラムを受け取ることを４１０で含むが、この第１のサイノグラムは、磁気共鳴データについての読出し方向の関数として視野角を表す。さらに、本方法は、第１のサイノグラムを変換して、角周波数軸を有するオシログラムを形成することを４２０で含む。また４３０で、本方法は、オシログラムの選択された軸に沿って投影を形成することを含むが、この選択された軸は、角周波数軸に直交する。本方法は、投影内のピークをフィルタリングすることを４４０で含み、このフィルタリングされた投影を変換して第２のサイノグラムを形成することを４５０で含む。

図５は、一例によるシステム５００を示す。このシステムは、データ・レシーバ５１０、サイノグラム・ジェネレータ５２０、プロセッサ５３０、フィルタ・モジュール５４０、出力モジュール５５０を含むことができる。これらの要素は、別々のコンポーネントまたはモジュールとすること、あるいは、任意の組合せで一体化することまたは組み合わせることができる。さらに、これらの要素の２つ以上を単一のモジュールとして組み合わせることができる。データ・レシーバ５１０は、放射状に規則化された磁気共鳴データを受け取るように構成される。サイノグラム・ジェネレータ５２０は、磁気共鳴データに対する読出し方向の関数としての視野角に対応する第１のサイノグラムを作り出すように構成される。プロセッサ５３０は、角周波数軸を有するオシログラムを作り出すように構成される。オシログラムは、第１のサイノグラムを変換したものに相当する。フィルタ・モジュール５４０は、オシログラムの選択された軸に沿って形成された投影を選択的にフィルタリングするように構成される。選択された軸は、角周波数軸に直交する。出力モジュール５５０は、フィルタリングされた投影を変換したものに相当する第２のサイノグラムを形成するように構成される。システム５００は、コンピュータと図示のように構成されたハードウェアとの任意の組合せとして実装することができる。

追記
上の詳細な説明は、その一部分を形成する添付の図面への言及を含んでいる。図面には、本発明を実施することのできる具体的な実施形態が例証として示されている。本明細書では、これらの実施形態を「例」とも称している。こうした例は、図示または説明した要素に加えて他の要素を含むことがある。ただし、本発明者らは、図示または説明した要素のみ提供される例もまた企図している。さらに、本発明者らは、特定の例（もしくはその１つまたは複数の側面）または図示もしくは本明細書で説明した他の例（もしくはその１つまたは複数の側面）について、図示または説明した要素を任意に組み合わせたまたは置き換えた例（またはその１つまたは複数の側面）を使用する例も企図している。

本文書中で参照した全ての刊行物、特許、および特許文献は、個別的に参照により本明細書に組み込むかのごとく、その全体を参照により本明細書に組み込む。本文書と上記のように参照により組み込んだ文献とに矛盾した用法がある場合、組み込んだ参照文献における用法は本文書における用法の補助となるものと考えられたい。すなわち、融和することのできない矛盾の場合、本文書における用法が優先する。

本文書中では、特許文献では一般的であるように、「at least one（少なくとも１）」もしくは「one or more（１または複数）」の使用または何か他の例とは無関係に、「a」または「an」の語を１または複数を含めるように使用している。本文書では、「or」の語を非排他的であることを指すのに用いており、すなわち、特に明記しない限り、「A or B」が「A but not B(ＡだがＢでない）」、「B but not A（ＢだがＡでない）」および「A and B（ＡおよびＢ）」を含むように「or」の語を使用している。添付の特許請求の範囲では、「including」および「in which」の語が、ぞれぞれ「comprising」および「wherein」の語の平易な英語の等価物として使用されている。また、添付の特許請求の範囲では、「including」および「comprising」の語は限定的ではなく、すなわち、請求項中でこうした語の後に挙げる要素とは他の要素を含むシステム、デバイス、物品、または方法が、依然としてその請求項の範囲内のものであるとみなされる。さらに、添付の特許請求の範囲では、「first（第１）」、「second（第２）」、「third（第３）」等の語を単なるラベルとして使用しており、これらは、その対象に数的条件を課すことを意図しない。

本明細書に記載する方法の例は、少なくとも部分的には機械またはコンピュータで実施することが可能である。いくつかの例には、上記の例で説明した方法を実施するための電子デバイスを構成するのに使用可能な命令を用いて符号化したコンピュータ読取り可能な媒体または機械読取り可能な媒体が含まれることがある。こうした方法の一実装例には、マイクロ符号、アセンブリ言語符号、高水準言語符号等の符号が含まれることがある。こうした符号には、様々な方法を実行するためのコンピュータ読取り可能な命令が含まれることがある。この符号は、コンピュータ・プログラム製品の一部分を形成することができる。さらに、この符号は、実行中または他の時間に、１つまたは複数の揮発性または不揮発性の有形的コンピュータ読取り可能媒体に有形的に記憶することができる。これらの有形的コンピュータ読取り可能媒体の例には、ハードディスク、リムーバブル磁気ディスク、リムーバブル光ディスク（たとえば、コンパクト・ディスクおよびデジタル・ビデオ・ディスク）、磁気カセット、メモリカードまたはメモリスティック、ランダム・アクセス・メモリ（ＲＡＭ）、リード・オンリ・メモリ（ＲＯＭ）等が含まれることがあるが、これらに限定されない。

上の説明は、例証となることを意図しており、限定することを意図していない。たとえば、上述の例（またはそれら例の１つまたは複数の側面）を、互いに組み合わせて用いることができる。上の説明を検討することに応じて当業者等により、他の実施形態が用いられる可能性がある。読者が技術的開示内容の特性を迅速に確かめることが可能になるように、米国特許法施行規則（３７Ｃ．Ｆ．Ｒ）１．７２（ｂ）に従って要約書を提供している。この要約書は、それが特許請求の範囲または請求項の意味を解釈または制限するのに用いられないという理解の基に提出したものである。また、上の詳細な説明では、様々な特徴をまとめて、開示内容を合理化することができる。このことは、特許請求されていない開示された特徴が請求項には必須であることを意図するものと解釈されるべきではない。反対に、発明主題は、開示された特定の実施形態の特徴全てには存在しないことがある。したがって、ここで、添付の特許請求の範囲を詳細な説明に組み込み、各請求項は、単独の実施形態として独立し、また、これら実施形態を様々な組合せまたは順序で互いに組み合わせることが可能であることが企図される。本発明の範囲は、添付の特許請求の範囲に関連して、こうした特許請求の範囲の権利が与えられる均等物の完全な範囲とともに確定するべきである。

Claims (22)

1. 放射状に規則化された磁気共鳴データを受け取るように構成されたデータ・レシーバと、
   前記磁気共鳴データに対する読出し方向の関数としての視野角に対応する第１のサイノグラムを作り出すように構成されたサイノグラム・ジェネレータと、
   角周波数軸を有し、前記第１のサイノグラムを変換したものに相当するオシログラムを作り出すように構成されたプロセッサと、
   前記オシログラムの選択された軸に沿って形成された投影を選択的にフィルタリングするように構成されたフィルタ・モジュールであって、前記選択された軸が、前記角周波数軸に直交する、フィルタ・モジュールと、
   前記フィルタリングされた投影を変換したものに相当する第２のサイノグラムを形成するように構成された出力モジュールと、
   を備えるシステム。
2. 前記データ・レシーバが、球状に規則化されたデータを受け取るように構成される、
   請求項１に記載のシステム。
3. 前記データ・レシーバが、２次元データまたは３次元データのうちの少なくとも１つを受け取るように構成される、
   請求項１または２のうちのいずれか１項に記載のシステム。
4. 前記プロセッサが、フーリエ変換または逆フーリエ変換のうちの少なくとも１つを実装するように構成される、
   請求項１〜３のうちのいずれか１項に記載のシステム。
5. 前記プロセッサが、前記視野角に沿って、かつ、前記読出し方向に沿ったフーリエ変換を実装するように構成される、
   請求項１〜４のうちのいずれか１項に記載のシステム。
6. 前記フィルタ・モジュールが、前記オシログラムの振幅成分を使用することにより前記投影を形成するように構成される、
   請求項１〜５のうちのいずれか１項に記載のシステム。
7. 前記フィルタ・モジュールが、前記オシログラムの速軸と一直線にそろった方向に１次元投影を形成するように構成される、
   請求項１〜６のいずれか１項に記載のシステム。
8. 前記フィルタ・モジュールが、ノッチ・フィルタおよびマスク・フィルタのうちの少なくとも１つを実装するように構成される、
   請求項１〜７のいずれか１項に記載のシステム。
9. 前記フィルタ・モジュールが、メジアン・フィルタを実装するように構成される、
   請求項１〜８のいずれか１項に記載のシステム。
10. 前記フィルタ・モジュールが、ピーク検出器を実装することを含む、
    請求項１〜９のいずれか１項に記載のシステム。
11. 前記出力モジュールが、前記第２のサイノグラムを使用することにより画像を作り出すように構成される、
    請求項１〜１０のいずれか１項に記載のシステム。
12. 放射状に規則化された磁気共鳴データに対応し、前記磁気共鳴データに対する読出し方向の関数としての視野角を表す第１のサイノグラムを受け取ることと、
    角周波数軸を有するオシログラムを形成するために、前記第１のサイノグラムを変換することと、
    前記オシログラムの選択された軸に沿って投影を形成することであって、前記選択された軸が、前記角周波数軸に直交することと、
    前記投影内のピークをフィルタリングすることと、
    第２のサイノグラムを形成するために、前記フィルタリングされた投影を変換することと、
    を含む方法。
13. 放射状に規則化された磁気共鳴データに対応する前記第１のサイノグラムを受け取ることが、球状に規則化されたデータを受け取ることを含む、
    請求項１２に記載の方法。
14. 放射状に規則化された磁気共鳴データに対応する前記第１のサイノグラムを受け取ることが、２次元データまたは３次元データのうちの少なくとも１つを受け取ることを含む、
    請求項１２または１３のうちのいずれか１項に記載の方法。
15. 変換することが、フーリエ変換または逆フーリエ変換のうちの少なくとも１つを実行することを含む、
    請求項１２〜１４のうちのいずれか１項に記載の方法。
16. オシログラムを形成するために前記第１のサイノグラムを変換することが、前記視野角に沿ってフーリエ変換を実行することと、前記読出し方向に沿ってフーリエ変換を実行することとを含む、
    請求項１２〜１５のうちのいずれか１項に記載の方法。
17. 前記投影を形成することが、前記オシログラムの振幅成分を使用することにより前記投影を形成することを含む、
    請求項１２〜１６のうちのいずれか１項に記載の方法。
18. 前記選択された軸に沿って前記投影を形成することが、前記オシログラムの速軸と一直線にそろった方向に１次元投影を形成することを含む、
    請求項１２〜１７のうちのいずれか１項に記載の方法。
19. フィルタリングすることが、ノッチ・フィルタリングを行うことまたはマスク・フィルタリングを行うことのうちの少なくとも１つを含む、
    請求項１２〜１８のうちのいずれか１項に記載の方法。
20. フィルタリングすることが、メジアン・フィルタ・アルゴリズムを実装することを含む、
    請求項１２〜１９のうちのいずれか１項に記載の方法。
21. フィルタリングすることが、ピーク検出アルゴリズムを実装することを含む、
    請求項１２〜２０のうちのいずれか１項に記載の方法。
22. 前記第２のサイノグラムを使用することにより画像を生成することを更に含む、
    請求項１２〜２１のうちのいずれか１項に記載の方法。

Images