JPH0212080A - 磁気共鳴スペクトロメータ - Google Patents
磁気共鳴スペクトロメータInfo
- Publication number
- JPH0212080A JPH0212080A JP1076413A JP7641389A JPH0212080A JP H0212080 A JPH0212080 A JP H0212080A JP 1076413 A JP1076413 A JP 1076413A JP 7641389 A JP7641389 A JP 7641389A JP H0212080 A JPH0212080 A JP H0212080A
- Authority
- JP
- Japan
- Prior art keywords
- frequency
- signal
- circuit
- sampling frequency
- analog
- 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.)
- Pending
Links
- 238000005070 sampling Methods 0.000 claims abstract description 51
- 230000000694 effects Effects 0.000 abstract description 5
- 238000006243 chemical reaction Methods 0.000 abstract description 3
- 238000007689 inspection Methods 0.000 abstract description 2
- 238000001228 spectrum Methods 0.000 description 12
- 238000000034 method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000005415 magnetization Effects 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 241001233037 catfish Species 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000001627 detrimental effect Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- WQGWDDDVZFFDIG-UHFFFAOYSA-N pyrogallol Chemical compound OC1=CC=CC(O)=C1O WQGWDDDVZFFDIG-UHFFFAOYSA-N 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/32—Excitation or detection systems, e.g. using radio frequency signals
- G01R33/36—Electrical details, e.g. matching or coupling of the coil to the receiver
- G01R33/3621—NMR receivers or demodulators, e.g. preamplifiers, means for frequency modulation of the MR signal using a digital down converter, means for analog to digital conversion [ADC] or for filtering or processing of the MR signal such as bandpass filtering, resampling, decimation or interpolation
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
本発明は、受信部におけるアナログディジタル変換器が
スピン共鳴信号又は低周波数帯域に転置されたスピン共
鳴信号をディジタル信号に変換し、フーリエ変換回路が
アナログディジタル変換器に続く信号路に配設される磁
気共鳴スペクトロメータに関する。この種の磁気共鳴ス
ペクトロメータは叶−０３３６０４２８１に記載されて
いる。、 ＯＦ−Ｏ３３６０４２８１では、スピン共鳴
１５号はよヂベースバンドに変換されてからアナログデ
ィジタル変換器によりデジタルデータワードの列に変換
される。ベースバンドへの変換のため直角復調が必要に
なり、スピン共鳴信号は、９０°の位相シフトがある２
つの混合段に印加されるラーマ−周波数を有する発振と
混合される。
スピン共鳴信号又は低周波数帯域に転置されたスピン共
鳴信号をディジタル信号に変換し、フーリエ変換回路が
アナログディジタル変換器に続く信号路に配設される磁
気共鳴スペクトロメータに関する。この種の磁気共鳴ス
ペクトロメータは叶−０３３６０４２８１に記載されて
いる。、 ＯＦ−Ｏ３３６０４２８１では、スピン共鳴
１５号はよヂベースバンドに変換されてからアナログデ
ィジタル変換器によりデジタルデータワードの列に変換
される。ベースバンドへの変換のため直角復調が必要に
なり、スピン共鳴信号は、９０°の位相シフトがある２
つの混合段に印加されるラーマ−周波数を有する発振と
混合される。
公知のスペクトロメータでは２つの混合段の構成１よ、
正確に対称でなければならず、２つの発振の間の位相ジ
ットは正確に９０°で’３ ＋−ｊればならない。アナ
ログ混合段においては、この条件を満足するのはかなり
複雑である。ベースバンドにおける最低有用周波数が低
すぎるため混合段の出力信号を増幅するＤＣ増幅器が必
要であるという別の欠点らある。この増幅器の動作点が
ずれると、信号の処理中に再構成エラーを引き起こす擾
乱ドノフ］−が生じる。またベースバンド及び主周波数
の高周波での低周波数ノイズ（フリッカノイズとし称さ
れる）は特に悪影響をおよぼす。
正確に対称でなければならず、２つの発振の間の位相ジ
ットは正確に９０°で’３ ＋−ｊればならない。アナ
ログ混合段においては、この条件を満足するのはかなり
複雑である。ベースバンドにおける最低有用周波数が低
すぎるため混合段の出力信号を増幅するＤＣ増幅器が必
要であるという別の欠点らある。この増幅器の動作点が
ずれると、信号の処理中に再構成エラーを引き起こす擾
乱ドノフ］−が生じる。またベースバンド及び主周波数
の高周波での低周波数ノイズ（フリッカノイズとし称さ
れる）は特に悪影響をおよぼす。
本発明の目的は、前記の種類であって複雑なアブログ構
成要索４省略しえ、再構成がトリット効果に影響されな
いスペクトロメータを提供するにある。
成要索４省略しえ、再構成がトリット効果に影響されな
いスペクトロメータを提供するにある。
この目的は、アナログディジタル変換器とフーリエ変換
回路との間の信号路に、ディジタルフィルタとサンプリ
ング周波数低減回路とからなる少なくとも１つの回路が
設けられ、その出力信号のサンプリング周波数は入力信
号に対し整数係数だけ低減され、低減されたサンプリン
グ周波数はディジタルフィルタの出力での信号周波数よ
り非整数係数だけ低い本発明により達成される。
回路との間の信号路に、ディジタルフィルタとサンプリ
ング周波数低減回路とからなる少なくとも１つの回路が
設けられ、その出力信号のサンプリング周波数は入力信
号に対し整数係数だけ低減され、低減されたサンプリン
グ周波数はディジタルフィルタの出力での信号周波数よ
り非整数係数だけ低い本発明により達成される。
本発明における信号処理は、木質的にはディジタル式で
ある。従ってドリフト効果とは無縁である。しかし環時
点では、磁気共鳴検査中に発生する周波数帯域（数ＭＨ
ｚ ）のアナログ信号を実時間でフーリエ変換できるフ
ーリエ変換回路は利用できない。いまのところ経済的に
実現しつるのは数百ｋ Ｈｚまでの周波数についてのみ
である。。
ある。従ってドリフト効果とは無縁である。しかし環時
点では、磁気共鳴検査中に発生する周波数帯域（数ＭＨ
ｚ ）のアナログ信号を実時間でフーリエ変換できるフ
ーリエ変換回路は利用できない。いまのところ経済的に
実現しつるのは数百ｋ Ｈｚまでの周波数についてのみ
である。。
従って本発明においては、アナログディジタル変換器と
フーリエ変換回路との間にフーリエ変換回路により処理
可能となる程度までサンプリング周波数を低減する（少
なくとも）１つのサンプリング周波数低減回路が設けら
れる。これはつまり、スピン共鳴信号が最終的にフーリ
エ変換回路に供給される際におけるサンプリング周波数
が、信号自体の周波数より相当低いことを意味する。こ
れは、信号を歪めるスペクトル折返し効果（７ライアシ
ング）を防ぐためシンブリング周波数が信号周波数の少
なくとも２倍に選定される従来のサンプリング方法と対
立するものである。本発明は、スピン共鳴信号の中心周
波数は厄知であり、スピン共鳴信号の帯域がこの中心周
波数に対し狭いという事実を利用する。従って周波数に
ついては、サンプリング周波数は常に、アンプサンプリ
ングによりイｊ用な信号が、ベースバンド外の周波数帯
域であってサンプリング周波数の半分により定められる
周期的ラスク内に転置される。先行するディジタルフィ
ルタは、干渉信号が他の帯域から７ライアシングにより
この帯域に入るのを防ぐ。
フーリエ変換回路との間にフーリエ変換回路により処理
可能となる程度までサンプリング周波数を低減する（少
なくとも）１つのサンプリング周波数低減回路が設けら
れる。これはつまり、スピン共鳴信号が最終的にフーリ
エ変換回路に供給される際におけるサンプリング周波数
が、信号自体の周波数より相当低いことを意味する。こ
れは、信号を歪めるスペクトル折返し効果（７ライアシ
ング）を防ぐためシンブリング周波数が信号周波数の少
なくとも２倍に選定される従来のサンプリング方法と対
立するものである。本発明は、スピン共鳴信号の中心周
波数は厄知であり、スピン共鳴信号の帯域がこの中心周
波数に対し狭いという事実を利用する。従って周波数に
ついては、サンプリング周波数は常に、アンプサンプリ
ングによりイｊ用な信号が、ベースバンド外の周波数帯
域であってサンプリング周波数の半分により定められる
周期的ラスク内に転置される。先行するディジタルフィ
ルタは、干渉信号が他の帯域から７ライアシングにより
この帯域に入るのを防ぐ。
低減されたサンプリング周波数がディジタルフィルタの
出力における信号周波数より非整数係数だけ低い場合に
は、前述の如く常に有用信号がベースバンド、つまり周
波数ゼロ付近の周波数帯域に折込まれないようにできる
。従って、直角復調器からなる公知のスペクトロメータ
と異なり、かかるスペクトロメータでは、スピン共鳴信
号がフーリエ変換に先立って処理される信号路は２つで
なく１つしか必要とされない。
出力における信号周波数より非整数係数だけ低い場合に
は、前述の如く常に有用信号がベースバンド、つまり周
波数ゼロ付近の周波数帯域に折込まれないようにできる
。従って、直角復調器からなる公知のスペクトロメータ
と異なり、かかるスペクトロメータでは、スピン共鳴信
号がフーリエ変換に先立って処理される信号路は２つで
なく１つしか必要とされない。
サンプリング周波数を低減する係数を大きくしすぎると
ディジタルフィルタの装作を満足しにくくなるので、任
意に大きくとることはできない。
ディジタルフィルタの装作を満足しにくくなるので、任
意に大きくとることはできない。
充分大幅な４ノンブリング周波数の低減を行なうには、
信号路にディジタルフィルタとサンプリング周波数低減
回路とからなる回路を２つ続けて接続すればよい。他の
方法としては、アナログ入力信号の周波数とサンプリン
グ周波数との商を０．５より大きい非整数とすることが
ある。
信号路にディジタルフィルタとサンプリング周波数低減
回路とからなる回路を２つ続けて接続すればよい。他の
方法としては、アナログ入力信号の周波数とサンプリン
グ周波数との商を０．５より大きい非整数とすることが
ある。
アナログディジタル変換器のサンプリング周波数は、ア
ンプサンプリング及び関係するアライアシング効果が起
こるよう低くとっであるためアナログ人力信号はより低
い同波数帯域にディジタル的に転置される。信号周波数
と勺ンブリング周波数の商が非整数であるため、信号は
ベースバンドに遠さない。
ンプサンプリング及び関係するアライアシング効果が起
こるよう低くとっであるためアナログ人力信号はより低
い同波数帯域にディジタル的に転置される。信号周波数
と勺ンブリング周波数の商が非整数であるため、信号は
ベースバンドに遠さない。
第１図の参照番号１は、検査ゾーン２に例えば０２１王
の一様かつ定常的で図の平面に対し垂直に延在する磁場
を発生する装置を示す１．装置は、必要ならば勾配磁場
、つまり一様かつ定常的な磁場と同一方向に延在するが
、強度が磁場の方向又はそれと直交するん向に直線的に
変化する磁場ら発生する。この目的のため、装置１は定
常的磁場を発生する発生器３と様々な勾配磁場のための
電流を出力する発生器４とに結合される。送信モードで
は定常的磁場に対し直交して延在する無線周波数磁場を
発生し、受信モードでは検査ゾーンで発生するスピン共
鳴信号を受信する無線周波数」イル５も設けられている
１、無線周波数コイルはスペクトロメータ６の出力に接
続される。スペクトロメータは、送信モードでは無線周
波数１ネルギーを無線周波数コイル５に供給し、受信モ
ードではコイルに誘導されるスピン共１月信号を処叩ツ
る。
の一様かつ定常的で図の平面に対し垂直に延在する磁場
を発生する装置を示す１．装置は、必要ならば勾配磁場
、つまり一様かつ定常的な磁場と同一方向に延在するが
、強度が磁場の方向又はそれと直交するん向に直線的に
変化する磁場ら発生する。この目的のため、装置１は定
常的磁場を発生する発生器３と様々な勾配磁場のための
電流を出力する発生器４とに結合される。送信モードで
は定常的磁場に対し直交して延在する無線周波数磁場を
発生し、受信モードでは検査ゾーンで発生するスピン共
鳴信号を受信する無線周波数」イル５も設けられている
１、無線周波数コイルはスペクトロメータ６の出力に接
続される。スペクトロメータは、送信モードでは無線周
波数１ネルギーを無線周波数コイル５に供給し、受信モ
ードではコイルに誘導されるスピン共１月信号を処叩ツ
る。
発１器３及び４とスペクトロメータ６は、ディジタル中
央コニット７により制御される。そのためにｆイジタル
中央ユニット７は制御部７１を有り−る。ディジタル中
央コニットは、スペクトロメータ６から供給されるデー
タを処理して検査ゾーン２内の核磁化弁孔を再構成する
処理部７２を有する。この処理は、特に実時間で行なわ
れるノーリエ変換からなる。
央コニット７により制御される。そのためにｆイジタル
中央ユニット７は制御部７１を有り−る。ディジタル中
央コニットは、スペクトロメータ６から供給されるデー
タを処理して検査ゾーン２内の核磁化弁孔を再構成する
処理部７２を有する。この処理は、特に実時間で行なわ
れるノーリエ変換からなる。
第２図に示されるスペクトロメータの受信部のブロック
図から分る如く、検査ゾーンから無線周波数」イル５が
受信したＭ号は、前置増幅器６１へ供給される。
図から分る如く、検査ゾーンから無線周波数」イル５が
受信したＭ号は、前置増幅器６１へ供給される。
一様かつ定常的な磁場の誘導が０．２１ Ｔある場合、
スピン共鳴信号の周波数の平均値は約９ＭＨｚになる。
スピン共鳴信号の周波数の平均値は約９ＭＨｚになる。
帯域は、スピン共鳴信号の検出中に活性化される磁気勾
配場の勾配の強度及び検査ゾーンの寸法に依存する。し
かし帯域は一般に６０ ｋ １−１２を越えることはな
い。前置増幅器６１の出力のスペクトルは第３ａ図に示
されている。、第３ａ図においては、概略的に表わされ
たスピン共鳴信号のスペクトルは文字Ｓで示され、それ
より相当広い雑音スペクトルは文字Ｎで示される、。
配場の勾配の強度及び検査ゾーンの寸法に依存する。し
かし帯域は一般に６０ ｋ １−１２を越えることはな
い。前置増幅器６１の出力のスペクトルは第３ａ図に示
されている。、第３ａ図においては、概略的に表わされ
たスピン共鳴信号のスペクトルは文字Ｓで示され、それ
より相当広い雑音スペクトルは文字Ｎで示される、。
前置増幅器の出力信号は、中心周波数９ Ｍ Ｉ−１ｚ
を有するアナログ帯域フィルタ６２を介してアナログデ
ィジタル変換器６３に供給される。アブ・ログ帯域ノイ
ルタロ２の出力すにおける信号のスペクトルは第３ｂ図
に丞されているが、これは本質的には第３ａ図のスペク
トラムに対応するが、ただ第３ｂ図の破線Ｆで概略的に
表わされる伝達関数を有するＷＩ域フィルタにより雑音
帯域がカットされている。
を有するアナログ帯域フィルタ６２を介してアナログデ
ィジタル変換器６３に供給される。アブ・ログ帯域ノイ
ルタロ２の出力すにおける信号のスペクトルは第３ｂ図
に丞されているが、これは本質的には第３ａ図のスペク
トラムに対応するが、ただ第３ｂ図の破線Ｆで概略的に
表わされる伝達関数を有するＷＩ域フィルタにより雑音
帯域がカットされている。
アナログディジタル変換器６３は、アナログ人力信号を
１２ビツトのワード長のデジタルデータワードの列に変
換する。しかし従来の信号ダイナミクスの要件を満たす
には１０ビツトのワード良Ｃも充分である。前置増幅器
６１の利得は、スピン共鳴信号の雑音成分が、変換器６
３の量子化ステップより大きくなるよう選定しなければ
ならない。アナログスピン共鳴信号が４ノンブリングさ
れる周波数ｆａは７．２Ｍ Ｈｚとなる。このサンプリ
ングの結果、信号周波数ｆｓは、９ＭＨｚの値から、周
期ｆａで繰り返す−ｆａ２から十ｆａ２（３，６ＭＨｚ
）までの周波数帯域内にある１、８Ｍ Ｈｚの値ｆｓ
Ｉに変換される。有用周波数帯域外からの干渉信号が周
波数ｆｓｎ及びその近傍に変換されないようにするため
、帯域フィルタ６２は、周波数ｆｓから１．８ＭＨｚよ
り大きい周波数成分ができる限り抑斤されるよう設計さ
れねばならない。
１２ビツトのワード長のデジタルデータワードの列に変
換する。しかし従来の信号ダイナミクスの要件を満たす
には１０ビツトのワード良Ｃも充分である。前置増幅器
６１の利得は、スピン共鳴信号の雑音成分が、変換器６
３の量子化ステップより大きくなるよう選定しなければ
ならない。アナログスピン共鳴信号が４ノンブリングさ
れる周波数ｆａは７．２Ｍ Ｈｚとなる。このサンプリ
ングの結果、信号周波数ｆｓは、９ＭＨｚの値から、周
期ｆａで繰り返す−ｆａ２から十ｆａ２（３，６ＭＨｚ
）までの周波数帯域内にある１、８Ｍ Ｈｚの値ｆｓ
Ｉに変換される。有用周波数帯域外からの干渉信号が周
波数ｆｓｎ及びその近傍に変換されないようにするため
、帯域フィルタ６２は、周波数ｆｓから１．８ＭＨｚよ
り大きい周波数成分ができる限り抑斤されるよう設計さ
れねばならない。
第３Ｃ図は、アナログディジタル変換器６３の出力にお
けるスペクトルを示す。周波数ｒｓ＋ −１，８ＭＨｚ
は、差ｆｓ−ｆａの絶対値に由来する。
けるスペクトルを示す。周波数ｒｓ＋ −１，８ＭＨｚ
は、差ｆｓ−ｆａの絶対値に由来する。
それぞれサンプリング周波数ｆａ及びその高周波数に関
して対称の位置にあり、周波数ピロに対し周波数ｆｓｌ
を有する信号によってサンプリング周波数に対し同一位
置を占める別の成分もある。
して対称の位置にあり、周波数ピロに対し周波数ｆｓｌ
を有する信号によってサンプリング周波数に対し同一位
置を占める別の成分もある。
アナログディジタル変換器６３の出力（１は、アナログ
ディジタル変換器６３と同一のザンブリング信号数で動
作する帯域通過特性を有するデイジタルフィルタ６４に
供給される。この帯域フィルタの中心周波数は、周波数
ｆｓ ＋ （１，８ＭＨｚ ）にあり、このフィルタ
の伝達関数は、有用周波数帯域の他端における雑音が、
信号の接続する信号処理によって雑音成分が信号周波数
帯域に折返されないよう抑圧されるよう選択される。こ
のフィルタの出力ｄにおけるスペクトルは、第３ｄ図に
示されている（第３ｄ図では、第３ａ図乃至第３Ｃ図と
異なる周波数縮尺が用いられている）。
ディジタル変換器６３と同一のザンブリング信号数で動
作する帯域通過特性を有するデイジタルフィルタ６４に
供給される。この帯域フィルタの中心周波数は、周波数
ｆｓ ＋ （１，８ＭＨｚ ）にあり、このフィルタ
の伝達関数は、有用周波数帯域の他端における雑音が、
信号の接続する信号処理によって雑音成分が信号周波数
帯域に折返されないよう抑圧されるよう選択される。こ
のフィルタの出力ｄにおけるスペクトルは、第３ｄ図に
示されている（第３ｄ図では、第３ａ図乃至第３Ｃ図と
異なる周波数縮尺が用いられている）。
ディジタルフィルタの出力信号は、第１のサンプリング
周波数低減回路６５に供給される。この回路は、ｐ番目
毎のデータワードのみが次の処理に進むようにする。こ
こでｐは本実施例では５である整数係数である。従って
、出力ｅでのデータストリームは、この回路の入力での
データストリームに対し５分の１ （＝１／ｐ）に減少
する。つまり回路６５の出力には毎秒１４４万データワ
ードが現われる。
周波数低減回路６５に供給される。この回路は、ｐ番目
毎のデータワードのみが次の処理に進むようにする。こ
こでｐは本実施例では５である整数係数である。従って
、出力ｅでのデータストリームは、この回路の入力での
データストリームに対し５分の１ （＝１／ｐ）に減少
する。つまり回路６５の出力には毎秒１４４万データワ
ードが現われる。
この種のサンプリング及び低減回路は、様々な方法で実
現される。例えば、ゲート回路と、サンプリング周波数
でパルスを３１数し５パルス毎（鯰にリセットされるカ
ウンタとからなり、その時のデータワードが通るようゲ
ート回路を聞くようにすればよい。分周器を介してゲー
トにクロック周波数を供給してもよい。
現される。例えば、ゲート回路と、サンプリング周波数
でパルスを３１数し５パルス毎（鯰にリセットされるカ
ウンタとからなり、その時のデータワードが通るようゲ
ート回路を聞くようにすればよい。分周器を介してゲー
トにクロック周波数を供給してもよい。
第３ｅ図は、回路６５の出力におけるスペクトルを示す
。サンプリング周波数が、値ｆａの５分の１に低減した
結果、サンプリング周波数は新たにｆ ａ＋ −１，４
４ＭＨｚとされるのでスベク］・ルはｆａ、の周期で繰
り返す。折返し後信号周波数はやはり係数５により低減
されるので、回路６５の出力の信号周波数ｆｓｚは３６
０ｋＨｚ （ｆ Ｓ１／ｐ）にはる。
。サンプリング周波数が、値ｆａの５分の１に低減した
結果、サンプリング周波数は新たにｆ ａ＋ −１，４
４ＭＨｚとされるのでスベク］・ルはｆａ、の周期で繰
り返す。折返し後信号周波数はやはり係数５により低減
されるので、回路６５の出力の信号周波数ｆｓｚは３６
０ｋＨｚ （ｆ Ｓ１／ｐ）にはる。
サンプリング周波数の低減は、別のアンダサンプリング
（この場合はディジタル領域で行なわれる）に対応し、
以前は有用周波数から離されていた周波数が有用周波数
帯域へ折返される。例えば、搬送波ｆｓｌから１．０８
ＭＨｚの距離にあるアナログディジタル変換器６３の出
力Ｃの周波数成分は、回路６５によるアンダサンプリン
グにより丁度周波数ｆｓｚ −３６０ｋ）１２に折返さ
れる。アンダサンプリング以前この周波数帯域にまだ雑
音信号があるならば信号雑音比は影響を受ける。これを
避けるため帯域フィルタ６４は、７２０ｋＨｚの帯域（
つまりｆｓＩから３６０ｋＨｚの距離毎に）できる限り
抑圧しなければならない。しかし、抑圧が僅かに離れて
位置する周波数でのみなされるものならば欠点とはなら
ない。その場合成分は、既に０からｆａの帯１ｊｌ（第
３Ｃ図）から０からｆａｔの帯域へ折返されており、折
返された成分は、雑音のみが存在し有用信号は存在しな
い０からｒａｌの帯域の縁に位置することになる。
（この場合はディジタル領域で行なわれる）に対応し、
以前は有用周波数から離されていた周波数が有用周波数
帯域へ折返される。例えば、搬送波ｆｓｌから１．０８
ＭＨｚの距離にあるアナログディジタル変換器６３の出
力Ｃの周波数成分は、回路６５によるアンダサンプリン
グにより丁度周波数ｆｓｚ −３６０ｋ）１２に折返さ
れる。アンダサンプリング以前この周波数帯域にまだ雑
音信号があるならば信号雑音比は影響を受ける。これを
避けるため帯域フィルタ６４は、７２０ｋＨｚの帯域（
つまりｆｓＩから３６０ｋＨｚの距離毎に）できる限り
抑圧しなければならない。しかし、抑圧が僅かに離れて
位置する周波数でのみなされるものならば欠点とはなら
ない。その場合成分は、既に０からｆａの帯１ｊｌ（第
３Ｃ図）から０からｆａｔの帯域へ折返されており、折
返された成分は、雑音のみが存在し有用信号は存在しな
い０からｒａｌの帯域の縁に位置することになる。
以上のことから、周波数ｆｓ１及びｆ−３２は、常に正
確に０からｆａの帯域又はＯからｔ’ａ＋の帯域の中心
に位置するのが特に有利である。信号周波数が一方の側
へずれたとすると、フィルタ６２及び６４への要求はよ
り厳しいものになる。
確に０からｆａの帯域又はＯからｔ’ａ＋の帯域の中心
に位置するのが特に有利である。信号周波数が一方の側
へずれたとすると、フィルタ６２及び６４への要求はよ
り厳しいものになる。
ｎを正整数として次の関数が成り立っているなら、周波
数ｆｓｌが常に好ましい中心位置を占めることが証明さ
れる。
数ｆｓｌが常に好ましい中心位置を占めることが証明さ
れる。
この場合、アナログディジタル変換器の後での信号周波
数は、値ｆｓ＋”ｆｓ−ｆｎｆａを有する。本実施例で
はｎ＝１である。
数は、値ｆｓ＋”ｆｓ−ｆｎｆａを有する。本実施例で
はｎ＝１である。
次の関係が成り立っている場合も同様に好ましい周波数
ｆｓＩの位置が得られる。
ｆｓＩの位置が得られる。
ｔｓ
ｔａ＝ ４ｎ−１■
しかしサンプリング周波数が、式のに従って配分される
場合には、サンプリング後の有用周波数帯域の位置は逆
転する、１つまりサンプリング前にｆｓより古かった周
波数成分は、す゛ンプリング後にはｆ’ｓＩより低くな
る（その逆も成り立つ）。
場合には、サンプリング後の有用周波数帯域の位置は逆
転する、１つまりサンプリング前にｆｓより古かった周
波数成分は、す゛ンプリング後にはｆ’ｓＩより低くな
る（その逆も成り立つ）。
従って核磁化分布が間違って再構成されつる。しかし、
これは適切なアンダサンプリングをさらに行なってこの
周波数帯域の位置を再び逆転してもとの位置を古び占め
るようにづることで防止されるしかし、この像の逆転は
、フーリエ変換後のデータ処理中にソフトウＴ７により
打消すこともできる。
これは適切なアンダサンプリングをさらに行なってこの
周波数帯域の位置を再び逆転してもとの位置を古び占め
るようにづることで防止されるしかし、この像の逆転は
、フーリエ変換後のデータ処理中にソフトウＴ７により
打消すこともできる。
ｔナンプリング周波数ｆａが式（１）及び■により定め
られる値に対し僅かにずれることはかまわないが、商ｆ
ｓ／ｆａが整数となるほど大きくずれてはならない。有
用周波数帯は７ンダサンプリングにより周波数ｆ−Ｑに
対し対称の位置となり、負の周波数は正の周波数帯域に
折返される。かかる場合に再構成を行なうには、スピン
共鳴信号は、直角復調器における如く２つの平行な部分
で位相をシフトさせて処理しなければならない、。
られる値に対し僅かにずれることはかまわないが、商ｆ
ｓ／ｆａが整数となるほど大きくずれてはならない。有
用周波数帯は７ンダサンプリングにより周波数ｆ−Ｑに
対し対称の位置となり、負の周波数は正の周波数帯域に
折返される。かかる場合に再構成を行なうには、スピン
共鳴信号は、直角復調器における如く２つの平行な部分
で位相をシフトさせて処理しなければならない、。
例えばサンプリング周波数が回路６４１Ｃおいて係ａｐ
＝＝４だけ低減される場合も同じ望ましくない事情が生
じる。その場合それは１．８ＭＨｚとなって、丁度（低
減された）信号周波数ｆｓ１に等しくなり、ｆｓ２はゼ
ロとなる。
＝＝４だけ低減される場合も同じ望ましくない事情が生
じる。その場合それは１．８ＭＨｚとなって、丁度（低
減された）信号周波数ｆｓ１に等しくなり、ｆｓ２はゼ
ロとなる。
ｆｓ８ｆｓ＋で置き換え、ｆａをｆａ＋１’置き換えざ
えすれば最も好ましい位置がやはり得られる。容易に証
明される如くｐ＝５の場合に得られるサンプリング周波
数ｆ ａｔ −１，４４ＭＨｚはｎ−１について式（１
）を満足する。。
えすれば最も好ましい位置がやはり得られる。容易に証
明される如くｐ＝５の場合に得られるサンプリング周波
数ｆ ａｔ −１，４４ＭＨｚはｎ−１について式（１
）を満足する。。
サンプリング周波数低域回路６４の出力信号は、帯域通
過特性を有するディジタルフィルタ６６を介して、別の
サンプリング及び低減回路６７へ供給される。フィルタ
６６のりンブリング周波数はｆ ａｔ （１，４４Ｍ
Ｈｚ ＞となる。フィルタの中心周波数は３６０ｋＨｚ
であり、その帯域幅は約７２ｋ Ｈｚである。第３ｆ図
はフィルタの出力の周波数スベク１〜ルを示す。この周
波数スペクトルも、第３ｄ図及び第３ｅ図に対して縮尺
が拡大されている。
過特性を有するディジタルフィルタ６６を介して、別の
サンプリング及び低減回路６７へ供給される。フィルタ
６６のりンブリング周波数はｆ ａｔ （１，４４Ｍ
Ｈｚ ＞となる。フィルタの中心周波数は３６０ｋＨｚ
であり、その帯域幅は約７２ｋ Ｈｚである。第３ｆ図
はフィルタの出力の周波数スベク１〜ルを示す。この周
波数スペクトルも、第３ｄ図及び第３ｅ図に対して縮尺
が拡大されている。
サンプリング周波数低域回路６７は、やはりデータスト
リームを５分の１に低減する。つまり回路７７の入力に
環われる毎秒１４４万データワードの５番目毎のデータ
ワードである毎秒２８８．０００データワードのみが処
理される。従って回路６７の出力のサンプリング周波数
はｆａ＝２８８ｋＨｚの値まで減少する。信号周波数ｆ
Ｓ２については、この値Ｌｔｆ’ｓをｆ’ｓ２で置き換
えｆａｔをｆａ２で置き換えるなら式（１）によって定
められる条件を満足する。サンプリング周波数の低減の
結果、信号周波数ｆ’ｓ２も５分の１になる。つまり、
３６０ｋ ｌ−１ｚから７２ｋＨｚになる。第３ｇ図は
、各回路６７の出力におけるスペクトルを示す。
リームを５分の１に低減する。つまり回路７７の入力に
環われる毎秒１４４万データワードの５番目毎のデータ
ワードである毎秒２８８．０００データワードのみが処
理される。従って回路６７の出力のサンプリング周波数
はｆａ＝２８８ｋＨｚの値まで減少する。信号周波数ｆ
Ｓ２については、この値Ｌｔｆ’ｓをｆ’ｓ２で置き換
えｆａｔをｆａ２で置き換えるなら式（１）によって定
められる条件を満足する。サンプリング周波数の低減の
結果、信号周波数ｆ’ｓ２も５分の１になる。つまり、
３６０ｋ ｌ−１ｚから７２ｋＨｚになる。第３ｇ図は
、各回路６７の出力におけるスペクトルを示す。
回路６７の出力＠号は、処理回路７２の一部をなすフー
リエ変換回路７２０に供給される。入力信号用のクロッ
ク周波数は、周波数ｆａ２と同一であるから、数千デー
タワードが回路７２０において２乃〒３０ｍ５の長さを
有するサンプリング区間中に処理され、その間にスピン
共鳴信号はデータワードの列に変換され、それから回路
７２０がフーリエ変換により半分の大きさの幾つの離散
周波数に対して絶対値及びスペクトル成分の位相をＧ１
専する。この回路は、す”ンプリング値を同時に処理す
る複数のトランスビ」−夕からなるシステムからなるよ
うにできる。このようにして再構成される周波数スペク
トルは、０からｆ ａ （１４４ｋＨｚ ）の周波数帯
域全体をカバーする（第３ｇ図）。しかし、さらなる処
理に本質的な有用信号帯域は、これより相当に狭く、例
えば２０ｋＨｚとなるにすぎ４ｊい。従ってこの周波数
帯域に対応する回路７１０の出力信号のみが、つまり本
実施例では６２ｋＨｚと８２ｋＨｚの間の周波数を伴う
出力信号がさらなる処理を受ける。
リエ変換回路７２０に供給される。入力信号用のクロッ
ク周波数は、周波数ｆａ２と同一であるから、数千デー
タワードが回路７２０において２乃〒３０ｍ５の長さを
有するサンプリング区間中に処理され、その間にスピン
共鳴信号はデータワードの列に変換され、それから回路
７２０がフーリエ変換により半分の大きさの幾つの離散
周波数に対して絶対値及びスペクトル成分の位相をＧ１
専する。この回路は、す”ンプリング値を同時に処理す
る複数のトランスビ」−夕からなるシステムからなるよ
うにできる。このようにして再構成される周波数スペク
トルは、０からｆ ａ （１４４ｋＨｚ ）の周波数帯
域全体をカバーする（第３ｇ図）。しかし、さらなる処
理に本質的な有用信号帯域は、これより相当に狭く、例
えば２０ｋＨｚとなるにすぎ４ｊい。従ってこの周波数
帯域に対応する回路７１０の出力信号のみが、つまり本
実施例では６２ｋＨｚと８２ｋＨｚの間の周波数を伴う
出力信号がさらなる処理を受ける。
以上の記載では、磁気共鳴スペクトロメータは、スピン
共鳴信号がアナログディジタル変換器に直接供給される
ものとして説明された。しかし、これらの信号の周波数
が、アナログディジタル変換器が処理しつる周波数より
高い場合には、ディジタル化の前に混合段によって信号
を低い周波数帯域へ転置する必要がある。
共鳴信号がアナログディジタル変換器に直接供給される
ものとして説明された。しかし、これらの信号の周波数
が、アナログディジタル変換器が処理しつる周波数より
高い場合には、ディジタル化の前に混合段によって信号
を低い周波数帯域へ転置する必要がある。
第１図は、本発明が使用される磁気共１１！Ｓ断層撮影
装置の概略図、第２図は本発明によるスベクｌ−ロメー
タの受信部のブロック図、第３ａ図乃至第３ｇ図は第２
図のｌｉｉ！置の点ａ乃至点Ｑにおけるスペクトルを示
す図である。 １・・・装置、２・・・検査ゾーン、３．４・・・発生
器、５・・・無線周波数コイル、６・・・スペクトロメ
ータ、７・・・ディジタル中央ユニット、６１・・・前
置増幅器、６２・・・アナログ帯域フィルタ、６３・・
・アブ口グディジタル変換器、６４．６６・・・ディジ
タルフィルタ、６５．６７・・・サンプリング周波数低
減回路、７１・・・制御部、７２・・・処即部、７２０
１９．ツーり丁変換回路。 特訂出鎮人 エメ・ベー・フィリップス・フルーイラン
ベンファブリケン 同
装置の概略図、第２図は本発明によるスベクｌ−ロメー
タの受信部のブロック図、第３ａ図乃至第３ｇ図は第２
図のｌｉｉ！置の点ａ乃至点Ｑにおけるスペクトルを示
す図である。 １・・・装置、２・・・検査ゾーン、３．４・・・発生
器、５・・・無線周波数コイル、６・・・スペクトロメ
ータ、７・・・ディジタル中央ユニット、６１・・・前
置増幅器、６２・・・アナログ帯域フィルタ、６３・・
・アブ口グディジタル変換器、６４．６６・・・ディジ
タルフィルタ、６５．６７・・・サンプリング周波数低
減回路、７１・・・制御部、７２・・・処即部、７２０
１９．ツーり丁変換回路。 特訂出鎮人 エメ・ベー・フィリップス・フルーイラン
ベンファブリケン 同
Claims (3)
- （１）受信部におけるアナログディジタル変換器がスピ
ン共鳴信号又は低周波数帯域に転置されたスピン共鳴信
号をディジタル信号に変換し、フーリエ変換回路がアナ
ログディジタル変換器に続く信号路に配設される磁気共
鳴スペクトロメータであって、アナログディジタル変換
器（６３）とフーリエ変換回路（７２０）との間の信号
路に、ディジタルフィルタ（６４）とサンプリング周波
数低減回路（６５）とからなる少なくとも１つの回路が
設けられ、その出力信号のサンプリング周波数は入力信
号に対し整数係数だけ低減され、低減されたサンプリン
グ周波数（ｆａ＿１；ｆａ＿２）はディジタルフィルタ
の出力での信号周波数（ｆｓ＿１；ｆｓ＿２）より非整
数係数だけ低いことを特徴とする磁気共鳴スペクトロメ
ータ。 - （２）アナログディジタル変換器（６３）のアナログ入
力信号の周波数（ｆｓ）とサンプリング周波数（ｆａ）
との商は０．５より大きい非整数であることを特徴とす
る請求項１記載の磁気共鳴スペクトロメータ。 - （３）ｎを正整数とし、ａを＋１又は−１として、サン
プリング周波数低減回路及び／又はディジタルフィルタ
の入力信号の周波数ｆｓとサンプリング周波数ｆａとの
間にｆａ／ｆｓ＝４／（４ｎ＋ａ）の関係が成り立つこ
とを特徴とする請求項１又は２記載の磁気共鳴スペクト
ロメータ。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| DE3811066A DE3811066A1 (de) | 1988-03-31 | 1988-03-31 | Kernresonanz-spektrometer |
| DE3811066.0 | 1988-03-31 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH0212080A true JPH0212080A (ja) | 1990-01-17 |
Family
ID=6351215
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP1076413A Pending JPH0212080A (ja) | 1988-03-31 | 1989-03-28 | 磁気共鳴スペクトロメータ |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US4949040A (ja) |
| EP (1) | EP0336479B1 (ja) |
| JP (1) | JPH0212080A (ja) |
| DE (2) | DE3811066A1 (ja) |
| IL (1) | IL89784A0 (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5533870A (en) * | 1992-11-13 | 1996-07-09 | Kabushiki Kaisha Toyoda Jidoshokki Seisakusho | Piston type compressor |
| JP2011206287A (ja) * | 2010-03-30 | 2011-10-20 | Ge Medical Systems Global Technology Co Llc | 信号処理装置および磁気共鳴イメージング装置 |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4992736A (en) * | 1989-08-04 | 1991-02-12 | General Electric Company | Radio frequency receiver for a NMR instrument |
| US5109854A (en) * | 1990-02-23 | 1992-05-05 | Picker International, Inc. | Roll-over aliasing suppression in undersampled images |
| US5529068A (en) * | 1994-06-16 | 1996-06-25 | The Regents Of The University Of California | Synchronized digital signal processor for MRI reception |
| US5594341A (en) * | 1994-06-27 | 1997-01-14 | Varian Associates, Inc. | Nuclear magnetic resonance receiver, method and system |
| US5739691A (en) * | 1995-11-28 | 1998-04-14 | The Regents Of The University Of California | Multi-frequency digital low pass filter for magnetic resonance imaging |
| WO2003046597A1 (en) * | 2001-11-26 | 2003-06-05 | Koninklijke Philips Electronics N.V. | Magnetic resonance imaging method with reduced acoustic noise |
| US20090285336A1 (en) * | 2006-06-16 | 2009-11-19 | Michael Anthony Pugel | Wideband Out-Of-Band-Receiver |
| US8049505B2 (en) * | 2006-12-19 | 2011-11-01 | Koninlijke Philips Electronics N.V. | MRI system with direct digital receiver using resampling |
| CN102103195B (zh) * | 2009-12-18 | 2014-01-15 | 东软飞利浦医疗设备系统有限责任公司 | 一种宽频带数字磁共振射频接收实现装置及方法 |
| DE102011086288B4 (de) * | 2011-11-14 | 2014-10-09 | Siemens Aktiengesellschaft | Magnetresonanztomographie-System, Empfangsvorrichtung für ein solches System sowie Verfahren zum Gewinnen eines Bildsignals in dem System |
| US9747789B2 (en) * | 2013-03-31 | 2017-08-29 | Case Western Reserve University | Magnetic resonance imaging with switched-mode current-source amplifier having gallium nitride field effect transistors for parallel transmission in MRI |
| FR3034274B1 (fr) * | 2015-03-27 | 2017-03-24 | Stmicroelectronics Rousset | Procede de traitement d'un signal analogique issu d'un canal de transmission, en particulier un signal vehicule par courant porteur en ligne |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US32712A (en) * | 1861-07-02 | Shutter and door fastener | ||
| USRE32712E (en) | 1979-08-20 | 1988-07-12 | General Electric Company | Moving gradient zeugmatography |
| EP0132337A3 (en) * | 1983-07-21 | 1986-12-30 | The Regents Of The University Of California | Apparatus and method for reducing aliasing in sagittal or coronal nmr imaging |
| US4567439A (en) * | 1983-09-23 | 1986-01-28 | Texas Instruments Incorporated | Apparatus for measuring the magnitude of a magnetic field |
| JPS60104272A (ja) * | 1983-11-11 | 1985-06-08 | Jeol Ltd | 積算処理過程において疑似的フイルタ効果を得る方法及びそれを利用した核磁気共鳴測定方法 |
| US4593247A (en) * | 1983-11-25 | 1986-06-03 | General Electric Company | Method of NMR imaging with off-center zoom scan |
| US4673880A (en) * | 1985-08-16 | 1987-06-16 | Technicare Corporation | Phase sensitive detection in multislice magnetic resonance imaging systems |
| US4743851A (en) * | 1985-08-27 | 1988-05-10 | Resonex, Inc. | Apparatus and method for creating non-orthogonal magnetic resonance imaging |
| JPH0620435B2 (ja) * | 1985-09-25 | 1994-03-23 | 株式会社日立メディコ | 核磁気共鳴イメージング装置 |
| DE3604281A1 (de) * | 1986-02-12 | 1987-08-13 | Philips Patentverwaltung | Verfahren zur bestimmung der kernmagnetisierungsverteilung in einer schicht eines untersuchungsbereiches und kernspintomograph zur durchfuehrung des verfahrens |
-
1988
- 1988-03-31 DE DE3811066A patent/DE3811066A1/de not_active Withdrawn
-
1989
- 1989-03-20 US US07/325,630 patent/US4949040A/en not_active Expired - Fee Related
- 1989-03-24 EP EP89200757A patent/EP0336479B1/de not_active Expired - Lifetime
- 1989-03-24 DE DE58909771T patent/DE58909771D1/de not_active Expired - Fee Related
- 1989-03-28 JP JP1076413A patent/JPH0212080A/ja active Pending
- 1989-03-29 IL IL89784A patent/IL89784A0/xx unknown
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5533870A (en) * | 1992-11-13 | 1996-07-09 | Kabushiki Kaisha Toyoda Jidoshokki Seisakusho | Piston type compressor |
| JP2011206287A (ja) * | 2010-03-30 | 2011-10-20 | Ge Medical Systems Global Technology Co Llc | 信号処理装置および磁気共鳴イメージング装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| DE3811066A1 (de) | 1989-10-12 |
| DE58909771D1 (de) | 1997-03-13 |
| US4949040A (en) | 1990-08-14 |
| IL89784A0 (en) | 1989-09-28 |
| EP0336479A2 (de) | 1989-10-11 |
| EP0336479A3 (de) | 1991-01-02 |
| EP0336479B1 (de) | 1997-01-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6653833B2 (en) | Method for sampling a received radio-frequency signal, in particular a radio-frequency signal from a receiving coil in a magnetic resonance apparatus | |
| JPH0212080A (ja) | 磁気共鳴スペクトロメータ | |
| US5170123A (en) | Magnetic resonance imager with digital transmitter/receiver | |
| JPH08194043A (ja) | 核磁気共鳴レシーバ装置及び方法 | |
| DE69129155T2 (de) | Nmr-frequenzverriegelungskreis | |
| KR20140031795A (ko) | 자기 공명 신호들의 송신을 위한 어레인지먼트 | |
| JP6389508B2 (ja) | 無線周波数アンテナ装置、ｍｒイメージングシステム、磁気共鳴信号提供方法及びソフトウェアパッケージ | |
| CN102048538A (zh) | 磁共振成像装置、方法及探头单元控制、成像单元 | |
| US4777438A (en) | Multiple imaging with improved signal-to-noise ratio | |
| US5341136A (en) | Bandpass sigma-delta modulator for analog-to-digital converters | |
| US10371769B2 (en) | MRI apparatus and method using direct A/D of MR signals without frequency down conversion | |
| US4740750A (en) | Reception signal processing apparatus in nuclear magnetic resonance diagnostic apparatus | |
| CA2120746A1 (en) | Nmr receiver with sigma-delta a/d converter | |
| US4737713A (en) | Apparatus and method for processing an electrical signal and increasing a signal-to-noise ratio thereof | |
| EP1486002B1 (en) | Method and apparatus for analog-to-digital conversion | |
| US4873486A (en) | Magnetic resonance spectrometer | |
| Liu et al. | A novel digital magnetic resonance imaging spectrometer | |
| JPH06105824A (ja) | 磁気共鳴信号の処理装置およびその処理方法 | |
| EP1739445B1 (en) | Magnetic resonance involving digital quadrature detection with reference signals of variable frequency | |
| JP2555233B2 (ja) | 核磁気共鳴装置 | |
| JPH08299301A (ja) | Ｍｒイメージング装置 | |
| WO2002048719A2 (en) | Method and apparatus for accurate measurement of communications signals | |
| Redfield et al. | Analog filtering of large solvent signals for improved dynamic range in high-resolution NMR | |
| JPH09117425A (ja) | Ｍｒイメージング装置 | |
| JPH09117424A (ja) | 核磁気共鳴検査装置 |