JP4854438B2 - データ処理装置、データ処理方法、医用画像診断装置及び磁気共鳴イメージング装置 - Google Patents
データ処理装置、データ処理方法、医用画像診断装置及び磁気共鳴イメージング装置 Download PDFInfo
- Publication number
- JP4854438B2 JP4854438B2 JP2006251317A JP2006251317A JP4854438B2 JP 4854438 B2 JP4854438 B2 JP 4854438B2 JP 2006251317 A JP2006251317 A JP 2006251317A JP 2006251317 A JP2006251317 A JP 2006251317A JP 4854438 B2 JP4854438 B2 JP 4854438B2
- Authority
- JP
- Japan
- Prior art keywords
- data
- signal power
- data processing
- image
- space
- 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
Links
- 238000012545 processing Methods 0.000 title claims description 135
- 238000002595 magnetic resonance imaging Methods 0.000 title claims description 46
- 238000003672 processing method Methods 0.000 title claims description 28
- 238000002059 diagnostic imaging Methods 0.000 title description 2
- 238000000034 method Methods 0.000 claims description 72
- 230000008569 process Effects 0.000 claims description 20
- 239000002872 contrast media Substances 0.000 claims description 15
- 238000012935 Averaging Methods 0.000 claims description 13
- 238000002347 injection Methods 0.000 claims description 13
- 239000007924 injection Substances 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 12
- 239000000203 mixture Substances 0.000 claims description 12
- 230000014509 gene expression Effects 0.000 claims description 11
- 230000003068 static effect Effects 0.000 claims description 11
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 238000003786 synthesis reaction Methods 0.000 claims description 8
- 230000033001 locomotion Effects 0.000 claims description 7
- 238000003384 imaging method Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 3
- 238000002599 functional magnetic resonance imaging Methods 0.000 claims description 3
- 230000015556 catabolic process Effects 0.000 claims description 2
- 230000003247 decreasing effect Effects 0.000 claims description 2
- 238000006731 degradation reaction Methods 0.000 claims description 2
- 230000003796 beauty Effects 0.000 claims 1
- 230000006872 improvement Effects 0.000 description 28
- 238000002597 diffusion-weighted imaging Methods 0.000 description 24
- 230000006870 function Effects 0.000 description 24
- 238000010586 diagram Methods 0.000 description 14
- 238000007792 addition Methods 0.000 description 13
- 230000008859 change Effects 0.000 description 10
- 238000005481 NMR spectroscopy Methods 0.000 description 8
- 230000006866 deterioration Effects 0.000 description 8
- 230000000694 effects Effects 0.000 description 8
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 8
- 238000013480 data collection Methods 0.000 description 7
- 238000001914 filtration Methods 0.000 description 7
- 238000009792 diffusion process Methods 0.000 description 6
- 230000008439 repair process Effects 0.000 description 5
- 238000002598 diffusion tensor imaging Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000009499 grossing Methods 0.000 description 3
- 238000002603 single-photon emission computed tomography Methods 0.000 description 3
- 206010006322 Breath holding Diseases 0.000 description 2
- 230000003321 amplification Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000002542 deteriorative effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 238000002075 inversion recovery Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 210000005036 nerve Anatomy 0.000 description 1
- 210000004126 nerve fiber Anatomy 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Description
Jain AK. Fundamental Digital Image Processing. Presence Hall Information and System Science Series. p276-283 伊藤聡志、山田芳文 フレネル空間の複式解法を利用したMR映像のSNR改善法 Med Imag Tech 19(5), 355-369 (2001) 伊藤聡志、鈴木藤孝、山田芳文 フレネル変換信号の帯域分割効果を用いたMRI画像の反復的SNR改善法 Vol.20, No.3, 212-216 (2002) 川田聡、南茂夫 科学計測のための画像データ処理 CQ出版社、東京、1994
dynamic_studyでは時系列的に同じ処理対象データを収集するが、造影剤注入等の所要時点を境に、所要時点前のデータを基準データとして高いSNRで収集する。一方、所要時点以降のデータを処理対象データとして、時間分解能を向上させるためにSNRを低下させて収集する。この場合、所要時点前のデータのみを基準データとし、隣接するデータ点を加算平均したPsdをPsと推定すればよい。この場合、Psdの時間軸と基準データとの違いが問題になりうるが、適用においてはその差が無視できる程十分に小さい必要がある。なお、造影効果は低周波数に比べ高周波数成分では相対的に大きな周波数の変化にはならない。
図16は、dynamic_studyで、所要データのみSNR増大させ、所要データのみを基準データとして使用する場合の実施例を示す模式図である。ここで、Psd(S(base))は各スキャンの収集時間が同一ならば、造影剤注入前に収集した複数の基準データ群(S(basei))から、次の式によって加算平均して求められる。又は、MRIの検査やRIの検査ならば始めからNAQを多くし、CTの検査ならば線量を上げたりして基準データ群を収集する。
mean_Ps方式は、造影剤注入前に収集したデータのみを基準データとし、その基準データのPsdをPsと推定する場合に限定されない。mean_Ps方式は、造影剤注入以降に収集した処理対象データも含めて時系列データをすべて加算平均して基準データとし、その基準データのPsdをPsと推定する場合であってもよい。必要なら動きや位置ズレのある範囲内の又は補正後のデータを使用する。
f−MRIでは負荷のON/OFFを一定周期で繰り返しながら時間的に複数のデータを収集するが、基本的にはdynamic_studyの場合と同じであるので基本的には同じ方法が使える。基準データをNAQを多くして収集し、その基準データのPsdを負荷前後の時間軸データに適用してもよい。より実用的にはf−MRIでは負荷前と負荷後で2組のコントラストの異なるデータが収集されるので、Psdを、負荷前のデータのみからなる第1の基準データのPsd[1]と、負荷後のデータのみからなる第2の基準データのPsd[2]との2種類に分けてもよい。負荷前後の全加算平均を用いる場合に比べ、PsdのSNRは1/√2に低減するがコントラストの相違によるPsdの相違は最小化できる。
検査が異なっても同一種類のパラメータデータを使用することが可能となる。理想的なデータとして画像を変形して標準化すればいかなる患者のデータでも標準データを使用可能である。又は、同一患者の過去のデータがあればそれを使用可能であるので理想的なデータの新たな収集が不要となる。
あるまとまった同一の検査単位で収集する同一患者の異種パラメータ画像のデータを使用することもできる。SNRがよくなる方のパラメータをPsに用いる。MRIの検査ならばT2強調(T2W)に対しT1強調(T1W)が、FLAIR(Fluid Attenuated Inversion Recovery)に対しT2W等が適用可能である。
造影を行なう場合は、通常は同一の検査単位毎に造影前後で一回ずつ撮像する。この場合のmean_Ps方式では、造影前のデータについては、造影前のデータのみを基準データとし、その基準データのPsdをPsと推定する一方、造影後のデータについては、造影後のデータに造影前のデータをも加算して基準データとし、その基準データのPsdをPsと推定する。2種類が同一SNRであっても造影後のデータはPs_SNRR=sqrt(2)=1.4倍になり、SNRが改善される。
時間的に変化するデジタル化されたデータ(音声や動画)で、静止した時間内のデータを加算平均した基準データのPsdを用いてPsを推定し、以降発生するPs分布が類似の処理対象データに対して適用する。又は、時間的に変化に乏しい部分を抽出し、その部分にのみ適用する。特にフーリエ空間でのフィルタリングの場合はフーリエ空間での平行移動は、フーリエ空間で位相の変化であるから信号値の二乗に相当するPsの変化はないので多少の動きにはロバストである。ただし、単なるフレーム平均ではSNRが向上するが動きがあるとぶれてしまう。
ここではDTIデータへの適用を示し、各MPG(Motion Probing Gradient)軸のデータを平均化することによりランダムノイズの残差を低減する。
16 コンピュータシステム
51 データ収集手段
52 データ配置手段
53 信号パワー推定手段
55 データ処理手段
56 画像再構成手段
57 合成・差分処理演算手段
58 画像記録手段
59 表示制御手段
Claims (27)
- 処理対象データと類似で、かつ、SNR(Signal to Noise Ratio)が相対的に大きな別のデータである基準データを用いて信号パワーを推定する信号パワー推定手段と、
前記信号パワー推定手段で推定した信号パワーを基にしたWiener_Filterによって、前記処理対象データを処理するデータ処理手段と、を有することを特徴とするデータ処理装置。 - 前記信号パワー推定手段は、フーリエ空間に配置された前記基準データを用いて信号パワーを推定し、前記データ処理手段は、前記信号パワー推定手段で推定した信号パワーを基にしたWiener_Filterによって、前記フーリエ空間に配置された前記処理対象データを処理することを特徴とする請求項1に記載のデータ処理装置。
- 前記信号パワー推定手段は、FREBAS(Fresnel transform Band Splitting)空間に配置された前記基準データを用いて信号パワーを推定し、前記データ処理手段は、前記信号パワー推定手段で推定した信号パワーを基にしたWiener_Filterによって、前記FREBAS空間に配置された前記処理対象データを処理することを特徴とする請求項1に記載のデータ処理装置。
- 前記データ処理手段で処理した処理対象データを1組毎に2次元又は3次元のフーリエ変換に付して実空間の画像に再構成する画像再構成手段を有することを特徴とする請求項1に記載のデータ処理装置。
- 前記画像再構成手段で再構成した画像の合成処理や差分演算処理を行なう合成・差分演算処理手段を有することを特徴とする請求項4に記載のデータ処理装置。
- 前記合成処理を、画素毎に加算する処理又は最大値投影処理とすることを特徴とする請求項5に記載のデータ処理装置。
- 処理対象データと類似で、かつ、SNRが相対的に大きな別のデータである基準データを用いて信号パワーが推定される信号パワー推定工程と、
前記信号パワー推定工程で推定した信号パワーを基にしたWiener_Filterによって、前記処理対象データが処理されるデータ処理工程と、を有することを特徴とするデータ処理方法。 - 前記信号パワー推定工程では、フーリエ空間に配置された前記基準データを用いて信号パワーが推定され、前記データ処理工程では、前記信号パワー推定工程で推定した信号パワーを基にしたWiener_Filterによって、前記フーリエ空間に配置された前記処理対象データが処理されることを特徴とする請求項7に記載のデータ処理方法。
- 前記信号パワー推定工程では、FREBAS空間に配置された前記基準データを用いて信号パワーが推定され、前記データ処理工程では、前記信号パワー推定工程で推定した信号パワーを基にしたWiener_Filterによって、前記FREBAS空間に配置された前記処理対象データが処理されることを特徴とする請求項7に記載のデータ処理方法。
- 前記処理対象データとは別に前記基準データを収集することを特徴とする請求項7に記載のデータ処理方法。
- 複数の基準データを収集する場合、前記信号パワーを、前記複数の基準データを平均して算出することを特徴とする請求項7に記載のデータ処理方法。
- 複数の基準データを収集する場合、前記信号パワーを、前記複数の基準データを重み付き平均して算出することを特徴とする請求項7に記載のデータ処理方法。
- 前記信号パワーを、低周波成分については前記処理対象データの重み付けを大きくする一方、高周波成分については、前記処理対象データと前記処理対象データとは異なるデータとのうち前記信号パワーのSNRが大きい一方のデータの重み付けを大きくして算出することを特徴とする請求項12に記載のデータ処理方法。
- 造影剤注入前のデータを収集する一方、造影剤注入以降のデータを前記造影剤注入前のデータと比較してSNRを低下させて前記処理対象データとして収集し、前記前記造影剤注入前のデータを前記基本データとすることを特徴とする請求項7に記載のデータ処理方法。
- 造影剤注入前のデータを収集する一方、造影剤注入以降のデータを前記造影剤注入前のデータと比較してSNRを低下させて前記処理対象データとして収集し、前記造影剤注入前のデータのみならず、前記造影剤注入後のデータも含めて時系列データをすべて加算平均して前記基準データとすることを特徴とする請求項7に記載のデータ処理方法。
- 所要の刺激前後で分別されるデータ毎に、異なる前記基準データを用いることを特徴とする請求項7に記載のデータ処理方法。
- f−MRI(functional-Magnetic Resonance Imaging)では、負荷前と負荷後でコントラストの異なる2組のデータ毎に、2つの前記基準データを用いることを特徴とする請求項18に記載のデータ処理方法。
- 前記基準データを、異種パラメータの画像データとすることを特徴とする請求項7に記載のデータ処理方法。
- 前記基準データを、複数のデータ群のうち全て又は一部を加算平均したものとすることを特徴とする請求項7に記載のデータ処理方法。
- 前記データ群を、動きを補正した後のデータとすることを特徴とする請求項21に記載のデータ処理方法。
- 前記基準データを、MPG(Motion Probing Gradient)軸のデータ群を平均化したデータとすることを特徴とする請求項7に記載のデータ処理方法。
- 処理対象データと類似で、かつ、SNRが相対的に大きな別のデータである基準データを用いて信号パワーを推定する信号パワー推定手段と、
前記信号パワー推定手段で推定した信号パワーを基にしたWiener_Filterによって、前記処理対象データを処理するデータ処理手段と、
前記データ処理手段で処理した処理対象データを1組毎に2次元又は3次元のフーリエ変換又はFREBAS変換に付して実空間の画像に再構成する画像再構成手段と、
前記画像再構成手段で再構成した画像の合成処理や差分演算処理を行なう合成・差分演算処理手段と、を有することを特徴とする医用画像診断装置。 - 前記信号パワー推定手段は、フーリエ空間に配置された前記基準データを用いて信号パワーを推定し、前記データ処理手段は、前記信号パワー推定手段で推定した信号パワーを基にしたWiener_Filterによって、前記フーリエ空間に配置された前記処理対象データを処理することを特徴とする請求項24に記載の医用画像診断装置。
- 前記信号パワー推定手段は、FREBAS空間に配置された前記基準データを用いて信号パワーを推定し、前記データ処理手段は、前記信号パワー推定手段で推定した信号パワーを基にしたWiener_Filterによって、前記FREBAS空間に配置された前記処理対象データを処理することを特徴とする請求項24に記載の医用画像診断装置。
- 静磁場中の被検体に対して傾斜磁場及び高周波パルスを印加することで磁気共鳴信号を発生させる撮影部と、
前記磁気共鳴信号を検出するための高周波コイルと、
前記高周波コイルにより検出された磁気共鳴信号を周波数帯域空間に配置する収集手段と、
前記収集手段により前記周波数帯域分割空間に配置されたデータを再構成して複数枚の画像を生成する画像生成手段と、
前記周波数帯域分割空間においてSNRが高い帯域ほどゲインを大きくする一方、前記SNRが低い帯域ほどゲインを小さくするものであって、前記複数枚の画像のうち処理対象画像と類似で、かつ、SNRが相対的に大きな別の画像を生成した周波数帯域分割空間におけるデータの信号パワーを、前記処理対象画像の周波数帯域分割空間における信号パワーに適用するフィルタと、を備えたことを特徴とする磁気共鳴イメージング装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006251317A JP4854438B2 (ja) | 2005-11-28 | 2006-09-15 | データ処理装置、データ処理方法、医用画像診断装置及び磁気共鳴イメージング装置 |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005342174 | 2005-11-28 | ||
JP2005342174 | 2005-11-28 | ||
JP2006251317A JP4854438B2 (ja) | 2005-11-28 | 2006-09-15 | データ処理装置、データ処理方法、医用画像診断装置及び磁気共鳴イメージング装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2007167623A JP2007167623A (ja) | 2007-07-05 |
JP4854438B2 true JP4854438B2 (ja) | 2012-01-18 |
Family
ID=38294926
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2006251317A Expired - Fee Related JP4854438B2 (ja) | 2005-11-28 | 2006-09-15 | データ処理装置、データ処理方法、医用画像診断装置及び磁気共鳴イメージング装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4854438B2 (ja) |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4220592B2 (ja) * | 1998-07-13 | 2009-02-04 | 株式会社東芝 | Mri装置 |
FR2864299B1 (fr) * | 2003-12-19 | 2006-06-02 | Ge Med Sys Global Tech Co Llc | Procede optimise de traitement d'images radiologiques |
-
2006
- 2006-09-15 JP JP2006251317A patent/JP4854438B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2007167623A (ja) | 2007-07-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | Deep learning for image enhancement and correction in magnetic resonance imaging—state-of-the-art and challenges | |
JP4861821B2 (ja) | 磁気共鳴イメージング装置および画像データ補正装置 | |
Scott et al. | Beat-to-beat respiratory motion correction with near 100% efficiency: a quantitative assessment using high-resolution coronary artery imaging | |
JP6533571B2 (ja) | 磁気共鳴イメージング装置及び画像作成方法 | |
US11002815B2 (en) | System and method for reducing artifacts in echo planar magnetic resonance imaging | |
CN104067137A (zh) | 具有对流动伪影的抑制的mr成像 | |
US11796617B2 (en) | System and method for reconstruction of magnetic resonance images acquired with partial Fourier acquisition | |
Zoroofi et al. | MRI artifact cancellation due to rigid motion in the imaging plane | |
RU2531969C2 (ru) | Сегментация сердца при формировании изображений сердца в покое и при нагрузке | |
US11835612B2 (en) | System and method for motion correction of magnetic resonance image | |
JP2017529960A (ja) | アーチファクト抑制を有するプロペラmrイメージング | |
EP1790991B1 (en) | Data processing system and data processing method involving a Wiener filter | |
US20240138700A1 (en) | Medical image processing apparatus, method of medical image processing, and nonvolatile computer readable storage medium storing therein medical image processing program | |
Lalush | Magnetic resonance–derived improvements in PET imaging | |
US10302726B2 (en) | Image reconstruction for MRI using multiplexed sensitivity encoding | |
Kathiravan et al. | A review of magnetic resonance imaging techniques | |
JP4854438B2 (ja) | データ処理装置、データ処理方法、医用画像診断装置及び磁気共鳴イメージング装置 | |
JP5575695B2 (ja) | Mri装置 | |
WO2015170394A1 (ja) | 撮像装置、画像処理装置及び画像処理方法 | |
Susanto et al. | Application of Denoising Weighted Bilateral Filter and Curvelet Transform on Brain MR Imaging of Non-cooperative Patients | |
Wang et al. | Multiple B-value model-based residual network (MORN) for accelerated high-resolution diffusion-weighted imaging | |
JP2009219577A (ja) | 医用画像における撮像条件決定装置 | |
Gharabaghi et al. | Super-Resolution Reconstruction of MRI | |
Royuela del Val et al. | Cardio-respiratory motion estimation for compressed sensing reconstruction of free breathing 2D cine MRI | |
JP2023133213A (ja) | 磁気共鳴イメージング方法、磁気共鳴イメージング装置及びプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20090825 |
|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20100418 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20110712 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20110713 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110906 |
|
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: 20111004 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20111025 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20141104 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4854438 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313117 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
LAPS | Cancellation because of no payment of annual fees |