JP2022092150A - 画像生成装置、画像生成方法、および画像生成プログラム - Google Patents

画像生成装置、画像生成方法、および画像生成プログラム Download PDF

Info

Publication number
JP2022092150A
JP2022092150A JP2020204769A JP2020204769A JP2022092150A JP 2022092150 A JP2022092150 A JP 2022092150A JP 2020204769 A JP2020204769 A JP 2020204769A JP 2020204769 A JP2020204769 A JP 2020204769A JP 2022092150 A JP2022092150 A JP 2022092150A
Authority
JP
Japan
Prior art keywords
lead
image
out direction
magnetic resonance
coil
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
Application number
JP2020204769A
Other languages
English (en)
Inventor
英明 朽名
Hideaki Kuchina
秀則 竹島
Hidenori Takeshima
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Canon Medical Systems Corp
Original Assignee
Canon Medical Systems Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Canon Medical Systems Corp filed Critical Canon Medical Systems Corp
Priority to JP2020204769A priority Critical patent/JP2022092150A/ja
Priority to US17/496,199 priority patent/US11896361B2/en
Publication of JP2022092150A publication Critical patent/JP2022092150A/ja
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56545Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by finite or discrete sampling, e.g. Gibbs ringing, truncation artefacts, phase aliasing artefacts
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5611Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/4818MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G01R33/4824MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Pathology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Signal Processing (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

Figure 2022092150000001
【課題】画質を向上させた磁気共鳴画像を生成すること。
【解決手段】実施形態に係る画像生成装置は、取得部と、合成部と、画像生成部と、を備える。取得部は、第1のリードアウト方向と第1のリードアウト方向に交差する第2のリードアウト方向とを含む複数のリードアウト方向について収集された磁気共鳴データと、複数のリードアウト方向に対応しローパスフィルタの感度の分布を示す複数のフィルタ感度分布と、磁気共鳴データの収集に用いられる複数のコイルエレメントに対応する複数のコイル感度分布とを取得する。合成部は、フィルタ感度分布とコイル感度分布とをリードアウト方向ごとに合成することにより、リードアウト方向ごとの合成感度分布を生成する。画像生成部は、合成感度分布と磁気共鳴データとに基づいてMR画像を生成する。
【選択図】図5

Description

本明細書及び図面に開示の実施形態は、画像生成装置、画像生成方法、および画像生成プログラムに関する。
磁気共鳴イメージング(Magnetic Resonance Imaging:以下、MRIと呼ぶ)装置におけるイメージングの共通課題として、受信された磁気共鳴(Magnetic Resonance:以下、MRと呼ぶ))信号に対するサンプリングレートの不足(換言するとサンプリング間隔の逆数に相当する撮像視野(Field of view)の不足)に起因する画像の折り返りがある。この課題に対処するための技術として、サンプリングレートと合わせたローパスフィルタの設定や、MR信号に対するオーバーサンプリングなどといった技術がある。
しかしながら、例えば、ラディアルスキャンに代表されるような、リードアウト方向が変化するスキャンにおいては、当該技術におけるローパスフィルタをMR信号に対して用いると、リードアウト方向に応じてローパスフィルタによりカットされる信号の量が異なる。その結果、異なるリードアウト方向間でのk空間に配置されたMRデータの一貫性が損なわれることとなる。これにより、リードアウト方向が変化するスキャンに伴って生成されたMR画像には、例えば、ストリークと呼ばれる筋状のアーチファクトなどのノイズが生じ、画質が低下するという問題がある。画質低下を解消するためにローパスフィルタ(低域通過フィルタ)の透過(通過)帯域の拡大が必要とされるが、MRI装置における性能の限界やS/N(signal-noise ratio:信号雑音比)の低下とのトレードオフ(trade-off:両立性)などから、ローパスフィルタの通過帯域を十分に拡大できない(すべきでない)問題がある。
特開2020-115967号公報
本明細書及び図面に開示の実施形態が解決しようとする課題の一つは、画質を向上させた磁気共鳴画像を生成することにある。ただし、本明細書及び図面に開示の実施形態により解決しようとする課題は上記課題に限られない。後述する実施形態に示す各構成による各効果に対応する課題を他の課題として位置づけることもできる。
実施形態に係る画像生成装置は、取得部と、合成部と、画像生成部と、を備える。取得部は、第1のリードアウト方向と第1のリードアウト方向に交差する第2のリードアウト方向とを含む複数のリードアウト方向をについて収集されたMRデータと、複数のリードアウト方向に対応しローパスフィルタの感度の分布を示す複数のフィルタ感度分布と、MRデータの収集に用いられる複数のコイルエレメントに対応する複数のコイル感度分布とを取得する。合成部は、フィルタ感度分布とコイル感度分布とをリードアウト方向ごとに合成することにより、リードアウト方向ごとの合成感度分布を生成する。画像生成部は、合成感度分布とMRデータとに基づいてMR画像を生成する。
図1は、実施形態に係る画像生成装置の一例を示すブロック図。 図2は、実施形態に係る磁気共鳴イメージング装置の一例を示すブロック図。 図3は、実施形態に係る受信回路の構成の一例を示す図。 図4は、実施形態に係り、本スキャンにおいて、MR信号のリードアウトに関するk空間でのトラジェクトリーの一例を示す図。 図5は、実施形態に係る画像生成処理の手順の一例を示すフローチャート。 図6は、実施形態に係り、k空間におけるMRデータと、画像空間における信号領域、ユーザFOV、再構成領域、阻止帯域およびカットオフ周波数に対応する遮断位置と、フィルタ感度分布との一例を示す図。 図7は、実施形態に係り、k空間におけるMRデータと、画像空間における信号領域、ユーザFOV、再構成領域、阻止帯域およびカットオフ周波数に対応する遮断位置と、フィルタ感度分布との一例を示す図。 図8は、実施形態に係り、k空間におけるMRデータと、画像空間における信号領域、ユーザFOV、再構成領域、阻止帯域およびカットオフ周波数に対応する遮断位置と、フィルタ感度分布との一例を示す図。 図9は、比較例としてMRデータに対して既存の再構成手法により生成されたMR画像と、本実施形態における画像生成処理により生成されたMR画像との一例を示す図。
以下、図面を参照しながら、画像生成装置、画像生成方法、および画像生成プログラムの実施形態について詳細に説明する。図1は、画像生成装置1の一例を示すブロック図である。画像生成装置1は、例えば、本画像生成装置1における各種機能が搭載されたモダリティや、院内などにおけるサーバに搭載される。なお、画像生成装置1における各種機能は、医用画像管理システム(以下、PACS(Picture Archiving and Communication Systems)と呼ぶ)のサーバや、病院情報システム(以下、HIS(Hospital Information System)と呼ぶ)のサーバなどに搭載されてもよい。
また、モダリティは、例えば、磁気共鳴イメージング(Magnetic Resonance Imaging:以下、MRIと呼ぶ)装置、PET(Positron Emission Tomography:陽電子放出コンピュータ断層撮影)-MRI装置、SPECT(single photon emission computed tomography:単一光子放出コンピュータ断層撮影)-MRI装置などのMRIに関する医用画像診断装置である。以下、説明を具体的にするために、画像生成装置1は、MRI装置に搭載されているものとする。このとき、MRI装置は、処理回路15における各種機能を有することとなる。
(実施形態)
図2は、本実施形態に係るMRI装置100の一例を示す図である。図2に示すように、MRI装置100は、静磁場磁石101と、傾斜磁場コイル103と、傾斜磁場電源105と、寝台107と、寝台制御回路109と、送信回路113と、送信コイル115と、受信コイル117と、受信回路119と、撮像制御回路(撮像制御部)121と、システム制御回路(システム制御部)123と、メモリ13と、入力インターフェース127と、ディスプレイ129と、通信インターフェース11と、処理回路15と、を備える。なお、画像生成装置1は、通信インターフェース11と、メモリ13と、処理回路15とに加えて、入力インターフェース127とディスプレイ129とをさらに有していてもよい。
静磁場磁石101は、中空の略円筒状に形成された磁石である。静磁場磁石101は、内部の空間に略一様な静磁場を発生する。静磁場磁石101としては、例えば、超伝導磁石等が使用される。
傾斜磁場コイル103は、中空の略円筒形状に形成されたコイルであり、円筒形の冷却容器の内面側に配置される。傾斜磁場コイル103は、傾斜磁場電源105から個別に電流供給を受けて、互いに直交するX、Y、及びZの各軸に沿って磁場強度が変化する傾斜磁場を発生する。傾斜磁場コイル103によって発生されるX、Y、Z各軸の傾斜磁場は、例えば、スライス選択用傾斜磁場、位相エンコード用傾斜磁場および周波数エンコード用傾斜磁場を形成する。スライス選択用傾斜磁場は、任意に撮像断面を決めるために利用される。位相エンコード用傾斜磁場は、空間的位置に応じて磁気共鳴信号(以下、MR(Magnetic Resonance)信号と呼ぶ)の位相を変化させるために利用される。周波数エンコード用傾斜磁場は、空間的位置に応じてMR信号の周波数を変化させるために利用される。
傾斜磁場電源105は、撮像制御回路121の制御により、傾斜磁場コイル103に電流を供給する電源装置である。
寝台107は、被検体Pが載置される天板1071を備えた装置である。寝台107は、寝台制御回路109による制御のもと、被検体Pが載置された天板1071を、ボア111内へ挿入する。
寝台制御回路109は、寝台107を制御する回路である。寝台制御回路109は、入出力インターフェース17を介した操作者の指示により寝台107を駆動することで、天板1071を長手方向および上下方向、場合によっては左右方向へ移動させる。
送信回路113は、撮像制御回路121の制御により、ラーモア周波数で変調された高周波パルスを送信コイル115に供給する。例えば、送信回路113は、発振部や位相選択部、周波数変換部、振幅変調部、RFアンプなどを有する。発振部は、静磁場中における対象原子核に固有の共鳴周波数のRFパルスを発生する。位相選択部は、発振部によって発生したRFパルスの位相を選択する。周波数変換部は、位相選択部から出力されたRFパルスの周波数を変換する。振幅変調部は、周波数変換部から出力されたRFパルスの振幅を例えばsinc関数に従って変調する。RFアンプは、振幅変調部から出力されたRFパルスを増幅して送信コイル115に供給する。
送信コイル115は、傾斜磁場コイル103の内側に配置されたRF(Radio Frequency)コイルである。送信コイル115は、送信回路113からの出力に応じて、高周波磁場に相当するRFパルスを発生する。
受信コイル117は、傾斜磁場コイル103の内側に配置されたRFコイルである。受信コイル117は、高周波磁場によって被検体Pから放射されるMR信号を受信する。受信コイル117は、受信されたMR信号を受信回路119へ出力する。受信コイル117は、例えば、1以上、典型的には複数のコイルエレメントを有するコイルアレイである。以下、説明を具体的にするために、受信コイル117は、複数のコイルエレメントを有するコイルアレイとして説明する。
なお、受信コイル117は、一つのコイルエレメントにより構成されてもよい。また、図2において送信コイル115と受信コイル117とは別個のRFコイルとして記載されているが、送信コイル115と受信コイル117とは、一体化された送受信コイルとして実施されてもよい。送受信コイルは、被検体Pの撮像部位に対応し、例えば、頭部コイルのような局所的な送受信RFコイルである。
受信回路119は、撮像制御回路121の制御により、受信コイル117から出力されたMR信号に基づいて、デジタルのMR信号(以下、MRデータと呼ぶ)を生成する。具体的には、受信回路119は、受信コイル117から出力されたMR信号に対して、検波、フィルタリングなどの信号処理を施した後、当該信号処理が施されたデータに対してアナログ/デジタル(A/D(Analog to Digital))変換(以下、A/D変換と呼ぶ)して、MRデータを生成する。受信回路119は、生成されたMRデータを、撮像制御回路121に出力する。例えば、MRデータは、コイルエレメントごとに生成され、コイルエレメントを識別するタグとともに、撮像制御回路121に出力される。
図3は、受信回路119の一例を示す図である。受信回路119は、例えば、ローパスフィルタ(Low Pass Filter:低域通過フィルタ)191と、A/D変換器193とを有する。なお、受信回路119には、ローパスフィルタ191およびA/D変換器193以外に、上記信号処理に対応する例えば検波器などの各種回路が搭載されてもよい。
ローパスフィルタ191には、撮像制御回路121を介して、設定機能153により設定されたカットオフ(遮断)周波数が入力される。換言すれば、ローパスフィルタ191における通過帯域は、設定機能153により設定される。なお、カットオフ周波数は、処理回路15から直接ローパスフィルタ191へ入力されてもよい。ローパスフィルタ191は、入力されたカットオフ周波数を用いて、MR信号をフィルタリングする。
A/D変換器193には、撮像制御回路121を介して、設定機能153により設定されたサンプリングレートが入力される。換言すれば、ローパスフィルタ191におけるサンプリング間隔は、設定機能153により設定される。A/D変換器193は、サンプリングレートに応じたサンプリングタイミングで、ローパスフィルタ191を通過したMR信号をサンプリングする。これにより、A/D変換器193は、MRデータを生成する。
撮像制御回路121は、処理回路15から出力された撮像プロトコルに従って、傾斜磁場電源105、送信回路113及び受信回路119等を制御し、被検体Pに対する撮像を行う。撮像プロトコルは、検査の種類に応じたパルスシーケンスを有する。撮像プロトコルには、傾斜磁場電源105により傾斜磁場コイル103に供給される電流の大きさ、傾斜磁場電源105により電流が傾斜磁場コイル103に供給されるタイミング、送信回路113により送信コイル115に供給される高周波パルスの大きさや時間幅、送信回路113により送信コイル115に高周波パルスが供給されるタイミング、受信コイル117によりMR信号が受信されるタイミング等が定義されている。撮像制御回路121は、傾斜磁場電源105、送信回路113及び受信回路119等を駆動して被検体Pを撮像した結果、受信回路119からMRデータを受信すると、受信したMRデータをメモリ13に記憶させる。
撮像制御回路121は、撮像に関するパルスシーケンスを実行することにより、MRデータを収集する。本実施形態において実行される撮像は、第1のリードアウト方向と第1のリードアウト方向に交差する第2のリードアウト方向とを含む複数のリードアウト方向に対応する複数のMR信号を取得するスキャンに相当する。
図4は、本スキャンにおいて、MR信号のリードアウトに関するk空間でのトラジェクトリーkTraの一例を示す図である。本スキャンは、例えば、図3に示すように、2次元的なラディアル収集R2D、3次元的なラディアル収集(スタックオブスターズ(Stack-of-stars)Sos、クッシュボール(koosh-ball)KB)、プロペラ(PROPELLER(periodically rotated overlapping parallel lines with enhanced reconstruction))収集PRP、その他の例えば2つのブレードによる収集BL2などがある。
以下、説明を具体的にするために、本スキャンは、2次元的なラディアル収集R2Dであるものとして説明する。2次元的なラディアル収集R2Dは、以下、簡便に、ラディアル収集と呼ぶ。被検体Pに対する本スキャンとして実施されるラディアル収集におけるリードアウトの回数、すなわちリードアウト方向のトラジェクトリーkTraの本数は、本スキャンに先立って予め設定されているものとする。
撮像制御回路121は、複数のコイル感度分布各々の生成に関するMRデータ(以下、感度データと呼ぶ)を、任意の撮像手法により収集する。複数のコイル感度分布は、MRデータの収集に用いられる複数のコイルエレメントに対応し、当該コイルエレメントの感度の分布を示す。複数のコイル感度分布各々は、複素数のデータで表現される。感度データの収集は、例えば、本スキャンとしてのラディアル収集R2Dに先立って、ロケータスキャンなどを含むプリスキャンにおいて撮像制御回路121により実行される。撮像制御回路121は、例えばプロセッサにより実現される。
「プロセッサ」という文言は、例えば、CPU、GPU(Graphics Processing Unit)、特定用途向け集積回路(Application Specific Integrated Circuit:ASIC)、プログラマブル論理デバイス(例えば、単純プログラマブル論理デバイス(Simple Programmable Logic Device:SPLD)、複合プログラマブル論理デバイス(Complex Programmable Logic Device:CPLD)、及びフィールドプログラマブルゲートアレイ(Field Programmable Gate Array:FPGA))等の回路を意味する。
システム制御回路123は、ハードウェア資源として図示していないプロセッサ、ROM(Read-Only Memory)やRAM(Random Access Memory)等のメモリ等を有し、システム制御機能によりMRI装置100を制御する。具体的には、システム制御回路123は、メモリに記憶されたシステム制御プログラムを読み出してメモリ上に展開し、展開されたシステム制御プログラムに従って本MRI装置100の各回路を制御する。
例えば、システム制御回路123は、入力インターフェース127を介して操作者から入力された撮像条件に基づいて、撮像プロトコルをメモリ13から読み出す。システム制御回路123は、撮像プロトコルを撮像制御回路121に送信し、被検体Pに対する撮像を制御する。システム制御回路123は、例えばプロセッサにより実現される。なお、システム制御回路123は、処理回路15に組み込まれてもよい。このとき、システム制御機能は処理回路15により実行され、処理回路15は、システム制御回路123の代替として機能する。システム制御回路123を実現するプロセッサは、上述と同様な内容なため、説明は省略する。
メモリ13は、システム制御回路123において実行されるシステム制御機能に関する各種プログラム、各種撮像プロトコル、撮像プロトコルを規定する複数の撮像パラメータを含む撮像条件等を記憶する。また、メモリ13は、処理回路15により実現される設定機能153、取得機能155、合成機能157、画像生成機能159を、コンピュータによって実行可能なプログラムの形態で記憶する。
また、メモリ13は、画像生成機能159により生成されたMR画像や、ロケータスキャンなどのプリスキャンにより生成されたプリスキャン画像を記憶する。プリスキャン画像は、例えば、本スキャンにおける撮像視野(以下、FOV(Field of view)と呼ぶ)を設定するための位置決め画像(ロケータ画像とも称される)や、本スキャンにおけるMR画像の生成(再構成)に用いられるコイル感度分布などである。例えば、メモリ13は、受信コイル117における複数のコイルエレメントにそれぞれ対応する複数のコイル感度分布を記憶する。また、メモリ13は、ロケータ画像において設定されたFOVを記憶する。
メモリ13は、プリスキャンの後に実行される本スキャンに関して、A/D変換器193において用いられるサンプリングレート(サンプリング周波数とも称される)またはサンプリング間隔を記憶する。また、メモリ13は、本スキャンに関して、ローパスフィルタ191において用いられるカットオフ周波数を記憶する。メモリ13は、本スキャンに関するMRデータおよび当該MRデータに基づいてMR画像を再構成するアルゴリズムを記憶する。
なお、メモリ13は、通信インターフェース11を介して受信された各種データを記憶してもよい。例えば、メモリ13は、放射線情報システム(RIS:Radiology Information System)等の医療機関内の情報処理システムから受信した被検体Pの検査オーダに関する情報(撮像対象部位、検査目的等)を記憶する。
メモリ13は、例えば、ROM、RAM、フラッシュメモリ等の半導体メモリ素子、HDD(Hard disk Drive)、SSD(Solid State Drive)、光ディスク等により実現される。また、メモリ13は、CD(Compact Disc)-ROMドライブやDVD(Digital Versatile Disc)ドライブ、フラッシュメモリ等の可搬型記憶媒体との間で種々の情報を読み書きする駆動装置等で実現されてもよい。
入力インターフェース127は、操作者からの各種指示(例えば、電源投入指示)や情報入力を受け付ける。入力インターフェース127は、例えば、トラックボール、スイッチボタン、マウス、キーボード、操作面へ触れることで入力操作を行うタッチパッド、表示画面とタッチパッドとが一体化されたタッチスクリーン、光学センサを用いた非接触入力回路、及び音声入力回路等によって実現される。入力インターフェース127は、処理回路15に接続されており、操作者から受け取った入力操作を電気信号へ変換し処理回路15へと出力する。なお、本明細書において入力インターフェース127は、マウス、キーボードなどの物理的な操作部品を備えるものだけに限られない。例えば、MRI装置100とは別体に設けられた外部の入力機器から入力操作に対応する電気信号を受け取り、この電気信号を制御回路へ出力する電気信号の処理回路も入力インターフェース127の例に含まれる。
入力インターフェース127は、ディスプレイ129に表示されたロケータ画像に対して、ユーザの指示によりFOVを入力する。具体的には、入力インターフェース127は、ディスプレイ129に表示されたロケータ画像において、ユーザによる範囲の設定指示によりFOVを入力する。また、入力インターフェース127は、検査オーダに基づくユーザの指示により、本スキャン関する各種撮像パラメータを入力する。
ディスプレイ129は、処理回路15またはシステム制御回路123による制御のもとで、各種のGUI(Graphical User Interface)や、処理回路15によって生成されたMR画像、ロケータ画像などのプリスキャン画像等を表示する。また、ディスプレイ129は、本スキャンやプリスキャンに関する撮像パラメータ撮像、および画像処理に関する各種情報などを表示する。ディスプレイ129は、例えば、CRTディスプレイや液晶ディスプレイ、有機ELディスプレイ、LEDディスプレイ、プラズマディスプレイ、又は当技術分野で知られている他の任意のディスプレイ、モニタ等の表示デバイスにより実現される。
通信インターフェース11は、例えば、HIS、PACSなどとの間でデータ通信を行う。通信インターフェース11および病院情報システムとの通信の規格は、如何なる規格であっても良いが、例えば、HL7(Hearth Level 7)、DICOM、又はその両方等が挙げられる。通信インターフェース11は、RIS等の医療機関内の情報処理システムから受信した被検体Pの検査オーダに関する情報(撮像対象部位、検査目的等)を受信する。また、MRI装置100に画像生成装置1が搭載されていない場合、画像生成装置1における通信インターフェース11は、被検体Pに対する検査において当該被検体Pを撮像するMRI装置100などからMRデータを受信する。このとき、受信されたMRデータは、メモリ13に記憶される。
処理回路15は、例えば、上述のプロセッサなどにより実現される。処理回路15は、設定機能153、取得機能155、合成機能157、画像生成機能159などを備える。設定機能153、取得機能155、合成機能157、画像生成機能159をそれぞれ実現する処理回路15は、設定部、取得部、合成部、および画像生成部に相当する。設定機能153、取得機能155、合成機能157、画像生成機能159などの各機能は、コンピュータによって実行可能なプログラムの形態でメモリ13に記憶されている。例えば、処理回路15は、プログラムをメモリ13から読み出し、実行することで各プログラムに対応する機能を実現する。換言すると、各プログラムを読み出した状態の処理回路15は、設定機能153、取得機能155、合成機能157、画像生成機能159などの各機能を有することとなる。
上記説明では、「プロセッサ」が各機能に対応するプログラムをメモリ13から読み出して実行する例を説明したが、実施形態はこれに限定されない。プロセッサが例えばCPUである場合、プロセッサはメモリ13に保存されたプログラムを読み出して実行することで機能を実現する。一方、プロセッサがASICである場合、メモリ13にプログラムを保存する代わりに、当該機能がプロセッサの回路内に論理回路として直接組み込まれる。なお、本実施形態の各プロセッサは、プロセッサごとに単一の回路として構成される場合に限らず、複数の独立した回路を組み合わせて1つのプロセッサとして構成し、その機能を実現するようにしてもよい。また、単一の記憶回路が各処理機能に対応するプログラムを記憶するものとして説明したが、複数の記憶回路を分散して配置して、処理回路15は個別の記憶回路から対応するプログラムを読み出す構成としても構わない。
処理回路15は、設定機能153により、ユーザにより入力されたFOV(以下ユーザFOVと呼ぶ)に応じて、サンプリングレートを設定する。具体的には、設定機能153はユーザFOVの逆数をサンプリング間隔として、サンプリングレートを設定する。設定機能153は、撮像制御回路121を介して、設定されたサンプリングレートを、A/D変換器193に出力する。これにより、MRデータは、撮像視野に応じて設定されたサンプリングレートでMR信号をサンプリングすることにより生成される。
なお、設定機能153は、ユーザFOVの大きさの逆数を定数倍することにより、オーバーサンプリングに対応するサンプリングレート(以下、オーバーサンプリングレートと呼ぶ)を設定してもよい。このとき、オーバーサンプリングレートの上限(以下、上限レートと呼ぶ)は、MRI装置100の性能限界に基づいて予め設定されているため、設定機能153は、上限レート以下に、オーバーサンプリングレートを設定する。受信回路119により生成されるMRデータは、ユーザFOVに応じたサンプリングレートもしくは、オーバーサンプリングレート(ユーザにより入力された撮像視野に対応するサンプリングレートより高いサインプリングレート)でMR信号をサンプリングしたk空間データに相当する。
オーバーサンプリングは、ユーザFOVに対して、撮像視野を拡張することに対応する。オーバーサンプリングは、例えば、プリスキャン画像において設定されたユーザFOVが撮像空間において水分子などプロトンが存在する領域より小さい場合などの収集レートを高める必要がある場合に設定される。このとき、再構成領域は、拡大されたFOVに対して実行されることとなる。以下、説明を具体的にするために、サンプリングレートは、ユーザFOVの逆数であるものとする。
処理回路15は、設定機能153により、設定されたサンプリングレートに基づいて、ローパスフィルタ191においてMR信号の通過帯域に関するカットオフ周波数を設定する。具体的には、設定機能153は、複数のリードアウト方向による被検体Pのスキャンである本スキャンにおけるユーザFOVと、当該本スキャンにおける傾斜磁場の強度とを用いて、カットオフ周波数を設定する。例えば、設定機能153は、以下の式におけるローパスフィルタ191の帯域を用いて、仮遮断周波数を設定する。
ユーザFOV[cm] = (2 × ローパスフィルタの帯域[Hz]) / 読み出し傾斜磁場強度[Hz/cm]
上式において、読み出し傾斜磁場強度は、本スキャンにおける傾斜磁場の強度に対応する。
設定機能153は、設定されたカットオフ周波数とリードアウト方向とに基づいて、リードアウト方向ごとに、k空間におけるローパスフィルタ191の感度の分布を示すフィルタ感度分布を生成する。設定機能153は、複数のリードアウト方向に対応する複数のフィルタ感度分布を、リードアウト方向と対応付けて、メモリ13に記憶させる。
処理回路15は、取得機能155により、複数のリードアウト方向について収集されたMRデータと、複数のリードアウト方向に対応する複数のフィルタ感度分布と、複数のコイルエレメントに対応する複数のコイル感度分布とを、メモリ13から取得する。取得されるMRデータは、第1のリードアウト方向と第1のリードアウト方向に交差する第2のリードアウト方向とを含む複数のリードアウト方向を用いて収集されたものである。複数のコイル感度分布は、MRデータの収集に関する受信コイル117における複数のコイルエレメントに対応し、当該コイルエレメントの感度の分布を示すものである。複数のフィルタ感度分布は、複数のリードアウト方向に対応しローパスフィルタ191の感度の分布を示すものである。
処理回路15は、合成機能157により、フィルタ感度分布とコイル感度分布とを、リードアウト方向ごとに合成する。これにより、合成機能157は、リードアウト方向ごとの合成感度分布を生成する。合成感度分布は、リードアウト方向ごとに、コイルエレメントの感度とローパスフィルタ191の感度とを合成した合成感度の分布を示すものである。すなわち、合成機能157は、複数のリードアウト方向に対応する複数の合成感度分布を生成する。合成機能157において実行される具体的な処理内容は、第1のリードアウト方向と当該第1のリードアウト方向に交差する第2のリードアウト方向とを含む複数のリードアウト方向を用いて収集されたMRデータに基づいてMR画像を生成する処理(以下、画像生成処理と呼ぶ)において、詳述する。
処理回路15は、画像生成機能159により、被検体Pに対するプリスキャンにより生成されたMRデータ(以下、プリスキャンデータと呼ぶ)を受信回路119から取得してk空間に配置し、k空間に配置されたプリスキャンデータに基づいてプリスキャン画像を生成する。画像生成機能159は、生成されたプリスキャン画像を、メモリ13に記憶させる。
例えば、処理回路15は、画像生成機能159により、ロケータ画像の生成に関するスキャンにより生成されたMRデータを受信回路119から取得してk空間に配置し、k空間に配置されたMRデータに基づいてロケータ画像を生成(再構成)する。また、画像生成機能159は、感度マップの生成に関するスキャンにより生成された感度データを受信回路119から取得してk空間に配置し、k空間に配置された感度データに基づいて感度マップを生成(再構成)する。ロケータ画像や感度マップなどの生成は、既存の再構成手法を用いることできるため、説明は省略する。
処理回路15は、画像生成機能159により、合成感度分布とMRデータとに基づいて、本スキャンに対応するMR画像を生成する。画像生成機能159は、例えば、合成感度分布を用いた感度エンコーディング(以下、SENSE(SENSitivity Encoding))をMRデータに適用することにより、当該MR画像を生成する。なお、画像生成機能159は、当該MR画像の生成に関する正則化項として、圧縮センシング(CS:compressed sensing)や超解像に関する正則化項を用いて、当該MR画像を生成してもよい。画像生成機能159において実行される具体的な処理内容は、画像生成処理と呼ぶ)において、詳述する。
以上のように構成された本実施形態のMRI装置100や画像生成装置1により実行される画像生成処理について、図5乃至図8を用いて説明する。図5は、画像生成処理の手順の一例を示すフローチャートである。まず、画像生成処理の実行に先立って実施される各種処理を説明し、次いで画像生成処理について説明する。
撮像制御回路121は、被検体Pに対してプリスキャンを実行する。処理回路15は、画像生成機能161により、プリスキャンデータに基づいてプリスキャン画像を生成する。画像生成機能161は、生成されたプリスキャン画像、例えば、ロケータ画像およびコイル感度マップを、メモリ13に記憶させる。システム制御回路123は、ロケータ画像をディスプレイ129に表示させる。
入力インターフェース127は、ユーザの指示により、ロケータ画像において、FOV(以下、ユーザFOVと呼ぶ)を入力する。ユーザFOVは、メモリ13に記憶される。処理回路15は、設定機能153により、ユーザFOVに応じて、サンプリングレートを設定し、設定されたサンプリングレートを、本スキャンの実施に先立ってA/D変換器193に出力する。
処理回路15は、設定機能153により、本スキャンにより受信されたMR信号に対して適用されるローパスフィルタ191のカットオフ周波数を設定する。具体的には、設定機能153は、本スキャンにおける傾斜磁場の強度とユーザFOVとに基づいて、カットオフ周波数を設定する。設定機能153は、設定されたカットオフ周波数を、ローパスフィルタ191に出力する。
撮像制御回路121は、被検体Pに対する本スキャンとしてラディアル収集を実行する。具体的には、撮像制御回路121は、本スキャンにおける一つのリードアウト方向に沿って撮像を実行する。これにより、受信コイル117は、当該リードアウト方向におけるMR信号を受信する。受信コイル117により受信されたMR信号は、受信回路119に出力される。このとき、MR信号は、検波器により検波される。
ローパスフィルタ191は、設定されたカットオフ周波数を用いて、検波後のMR信号をフィルタリングする。ローパスフィルタ191によるMR信号のフィルタリングは、ユーザFOVより外側の不必要な信号を低減する。これにより、MR信号などのS/Nが向上するフィルタリングされたMR信号は、A/D変換器193へ出力される。
A/D変換器193は、設定機能153により設定されたサンプリングレートを用いて、ローパスフィルタ191から出力されたMR信号をサンプリングする。これにより、A/D変換器193は、MRデータを生成する。すなわち、A/D変換器193は、ローパスフィルタ191を通過したMR信号に対して、当該サンプリングレートでA/D変換を実行し、MRデータを生成する。MRデータは、撮像制御回路121を介して、メモリ13に記憶される。本スキャンにおいて予め設定された全リードアウト方向(すなわち、リードアウト方向のトラジェクトリーkTraの全本数)に対して、MR信号の収集が完了していれば、画像生成処理が開始される。
(画像生成処理)
(ステップS501)
処理回路15は、取得機能155により、本スキャンにおいて生成されたMRデータと、本スキャンに先だつプリスキャンにより生成された複数のコイル感度分布と、リードアウト方向に応じた複数のフィルタ感度分布とを、メモリ13から取得する。なお、MRI装置100に画像生成装置1が搭載されていない場合、取得機能155は、通信インターフェース11およびネットワークを介して各種MRI装置から、当該MRデータと複数のコイル感度分布と複数のフィルタ感度分布とを取得する。
図6は、k空間ksにおけるリードアウト方向RODおよびサンプリング点SPと、画像空間ISにおける信号領域SA、ユーザFOV UF、再構成領域RR、阻止帯域SB、およびカットオフ周波数に対応する遮断位置CPと、フィルタ感度分布FSDとの一例を示す図である。図6における信号領域SAは、撮像空間において、水分子などプロトンが存在する領域に相当する。阻止帯域SBは、受信コイル117から出力されたMR信号がローパスフィルタ191により阻止される帯域を示している。
図6に示すように、リードアウト方向RODは、k空間ksの中心を通り、kx方向に平行である。図6に示すユーザFOV UFは、サンプリング間隔に対応するため、隣接するサンプリング点SPにより規定されるFOVと同一となる。図6に示すように、カットオフ周波数に対応する遮断位置CPは、フィルタ感度分布FSDにおける感度を示す曲線の変曲点に対応する。
図7は、k空間ksにおけるリードアウト方向RODおよびサンプリング点SPと、画像空間ISにおける信号領域SA、ユーザFOV UF、再構成領域RR、阻止帯域SB、およびカットオフ周波数に対応する遮断位置CPと、フィルタ感度分布FSDとの一例を示す図である。図6と図7との相違は、図7におけるリードアウト方向RODが、k空間ksの中心を通り、ky方向に平行となっていることにある。
図8は、k空間ksにおけるリードアウト方向RODおよびサンプリング点SPと、画像空間ISにおける信号領域SA、各種FOV、再構成領域RR、阻止帯域SB、およびカットオフ周波数に対応する遮断位置CPと、フィルタ感度分布FSDとの一例を示す図である。図6および図7と図8との相違は、図8におけるリードアウト方向RODが、k空間ksの中心を通り、kx方向またはky方向に対して45°傾いていることにある。図6乃至図8に示すように、フィルタ感度分布FSDは、かまぼこ状の感度の分布を有し、リードアウト方向に応じて異なる。
(ステップS502)
処理回路15は、合成機能157により、フィルタ感度分布とコイル感度分布とを、リードアウト方向ごとに合成し、リードアウト方向ごとの合成感度分布を生成する。コイルエレメントiにおけるコイル感度分布をSとし、リードアウト方向
Figure 2022092150000002
におけるフィルタ感度分布を
Figure 2022092150000003
とすると、合成機能157は、本スキャンにおけるk番目のリードアウトに関する合成感度分布Sを、以下の示す式のように、計算する。
Figure 2022092150000004
上式における合成感度分布Sは、コイル感度分布をマスクすることに相当する。合成機能157は、リードアウト方向ごとに、コイルエレメントの総数に応じた合成感度分布を算出する。算出された合成感度分布Sは、リードアウト方向と関連付けられて、メモリ13に記憶される。
(ステップS503)
処理回路15は、画像生成機能159により、合成感度分布SとMRデータとに基づいて、本スキャンに関するMR画像を生成する。例えば、画像生成機能159は、合成感度分布Sを用いたSENSEをMRデータに適用することにより、以下の式を用いて、MR画像を生成する。
Figure 2022092150000005
上式において、xは、MRデータのサンプリンレートに対応する撮像視野または拡大された視野における画像空間のデータ(本スキャン画像)、Fは、k番目のリードアウト方向に対するフーリエ変換、ykは、k番目のリードアウト方向におけるMRデータ、Sは、k番目のリードアウトに関する合成感度分布、R(x)は、本スキャン画像xに関する正則化項を示している。式(1)におけるL2ノルムは、ローパスフィルタ191におけるMR信号のカットオフを考慮しているため、MRデータとMR画像との一貫性を担保している。すなわち、式(1)により、MRデータのサンプリンレートに対応する撮像視野より大きい視野における画像を生成することが可能である。
なお、コイル感度マップを用いない場合は、上式における合成感度分布Sは、k番目のリードアウトに関するフィルタ感度分布に等しくなり、正則化項R(x)に本スキャン画像に関する各種仮定に対応する式が組み込まれる。(1)式を満たすようにxを決めることで、画像生成機能159は、本スキャン画像を生成する。また、画像生成機能159による本スキャン画像の再構成は、上記SENSEに限定されず、例えば、上式(1)のMR画像の生成に関する正則化項R(x)として、圧縮センシング(CS:compressed sensing)や超解像に関する正則化項を用いて、MR画像を生成してもよい。
処理回路15は、画像生成機能159により、生成されたMR画像を、メモリ13に記憶させる。このとき、システム制御回路123は、生成されたMR画像を、ディスプレイ129に表示させてもよい。また、MRI装置100に画像生成装置1が搭載されていない場合、生成されたMR画像は、通信インターフェース11およびネットワークを介して、例えば、PACSのサーバやHISのサーバなどに出力される。
以上に述べた実施形態に係る画像生成装置1によれば、第1のリードアウト方向と第1のリードアウト方向に交差する第2のリードアウト方向とを含む複数のリードアウト方向について収集されたMRデータと、複数のリードアウト方向に対応しローパスフィルタ191の感度の分布を示す複数のフィルタ感度分布と、MRデータの収集に用いられる複数のコイルエレメントに対応する複数のコイル感度分布とを取得し、フィルタ感度分布とコイル感度分布とをリードアウト方向ごとに合成することにより、リードアウト方向ごとの合成感度分布を生成し、合成感度分布とMRデータとに基づいて、MR画像を生成する。本画像生成装置1により生成されたMR画像の大きさは、MRデータのサンプリンレートに対応する撮像視野より大きくすることが可能である。
例えば、本実施形態における画像生成装置1によれば、合成感度分布を用いた感度エンコーディング(SENSE)をMRデータに適用することにより、MR画像を生成する。また、本実施形態における画像生成装置1によれば、MR画像の生成に関する正則化項R(x)として、圧縮センシング(CS)に関する正則化項を用いて、MR画像を生成してもよい。
以上のことから、実施形態に係る画像生成装置1によれば、MRI装置100における性能の限界やS/Nの低下とのトレードオフなどから、ローパスフィルタ191の通過帯域を十分に拡大できない場合であっても、k空間上においてローパスフィルタ191による影響を加味することにより、MRデータとの一貫性を担保して、MR画像を生成することが可能となる。これにより、本画像生成装置1によれば、ストリークの少ない高画質が画像を得ることが可能となる。
図10は、比較例としてローパスフィルタ191による影響を加味していないMRデータに対して既存の再構成手法により生成されたMR画像(1)と、本実施形態における画像生成処理により、ローパスフィルタ191による影響を加味して生成されたMR画像(2)との一例を示す図である。図10に示すように、本実施形態における画像生成処理により生成されたMR画像(2)は、MR画像(1)に比べて、ストリークアーチファクトなどが低減され、画質が向上している。したがって、本実施形態の画像生成装置1によれば、図10に示すように、MRI装置100の性能限界やS/Nの低下とのトレードオフなどによりローパスフィルタ191の通過帯域を十分に拡大できない場合であっても、本スキャンに関して、画質を向上させたMR画像を生成することができる。
本実施形態における技術的思想を画像生成方法で実現する場合、当該画像生成方法は、第1のリードアウト方向と第1のリードアウト方向に交差する第2のリードアウト方向とを含む複数のリードアウト方向について収集されたMRデータと、複数のリードアウト方向に対応しローパスフィルタ191の感度の分布を示す複数のフィルタ感度分布と、MRデータの収集に用いられる複数のコイルエレメントに対応する複数のコイル感度分布とを取得し、フィルタ感度分布とコイル感度分布とをリードアウト方向ごとに合成することにより、リードアウト方向ごとの合成感度分布を生成し、合成感度分布とMRデータとに基づいて、MR画像を生成する。本画像生成方法に関する画生成処理の手順および効果は、実施形態の記載と同様なため、説明は省略する。
本実施形態における技術的思想を画像生成プログラムで実現する場合、当該画像生成プログラムは、コンピュータに、第1のリードアウト方向と第1のリードアウト方向に交差する第2のリードアウト方向とを含む複数のリードアウト方向について収集されたMRデータと、複数のリードアウト方向に対応しローパスフィルタ191の感度の分布を示す複数のフィルタ感度分布と、MRデータの収集に用いられる複数のコイルエレメントに対応する複数のコイル感度分布とを取得し、フィルタ感度分布とコイル感度分布とをリードアウト方向ごとに合成することにより、リードアウト方向ごとの合成感度分布を生成し、合成感度分布とMRデータとに基づいて、MR画像を生成すること、を実現させる。
例えば、MRI装置100などのモダリティやPACSサーバなどにおけるコンピュータに画像生成プログラムをインストールし、これらをメモリ上で展開することによっても、画像生成処理を実現することができる。このとき、コンピュータに当該手法を実行させることのできるプログラムは、磁気ディスク(ハードディスクなど)、光ディスク(CD-ROM、DVDなど)、半導体メモリなどの記憶媒体に格納して頒布することも可能である。画像生成プログラムによる画像生成処理の手順および効果は、本実施形態と同様なため、説明は省略する。
以上説明した少なくとも1つの実施形態等によれば、画質を向上させた磁気共鳴画像を生成することができる。すなわち、少なくとも1つの実施形態等によれば、リードアウト方向が変化するスキャンにおいて、ローパスフィルタ191の通過帯域を十分に拡大できない場合であっても、k空間上で一貫性を保ってMR画像を生成することが可能となり、ストリークの少ない高画質なMR画像を得ることが可能となる。
いくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更、実施形態同士の組み合わせを行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
1 画像生成装置
11 通信インターフェース
13 メモリ
15 処理回路
100 磁気共鳴イメージング装置
101 静磁場磁石
103 傾斜磁場コイル
105 傾斜磁場電源
107 寝台
109 寝台制御回路
111 ボア
113 送信回路
115 送信コイル
117 受信コイル
119 受信回路
121 撮像制御回路
123 システム制御回路
127 入力インターフェース
129 ディスプレイ
153 設定機能
155 取得機能
157 合成機能
159 画像生成機能
191 ローパスフィルタ
193 A/D変換器

Claims (6)

  1. 第1のリードアウト方向と前記第1のリードアウト方向に交差する第2のリードアウト方向とを含む複数のリードアウト方向について収集された磁気共鳴データと、前記複数のリードアウト方向に対応しローパスフィルタの感度の分布を示す複数のフィルタ感度分布と、前記磁気共鳴データの収集に用いられる複数のコイルエレメントに対応する複数のコイル感度分布とを取得する取得部と、
    前記フィルタ感度分布と前記コイル感度分布とを前記リードアウト方向ごとに合成することにより、前記リードアウト方向ごとの合成感度分布を生成する合成部と、
    前記合成感度分布と前記磁気共鳴データとに基づいて、磁気共鳴画像を生成する画像生成部と、
    を備える画像生成装置。
  2. 前記磁気共鳴画像の大きさは、前記磁気共鳴データのサンプリングレートに対応する撮像視野より大きい、
    請求項1に記載の画像生成装置。
  3. 前記画像生成部は、前記合成感度分布を用いた感度エンコーディングを前記磁気共鳴データに適用することにより、前記磁気共鳴画像を生成する、
    請求項1または2に記載の画像生成装置。
  4. 前記画像生成部は、前記磁気共鳴画像の生成に関する正則化項として、圧縮センシングに関する正則化項を用いて、前記磁気共鳴画像を生成する、
    請求項1または2に記載の画像生成装置。
  5. 第1のリードアウト方向と前記第1のリードアウト方向に交差する第2のリードアウト方向とを含む複数のリードアウト方向について収集された磁気共鳴データと、前記複数のリードアウト方向に対応しローパスフィルタの感度の分布を示す複数のフィルタ感度分布と、前記磁気共鳴データの収集に関する複数のコイルエレメントに対応する複数のコイル感度分布とを取得し、
    前記フィルタ感度分布と前記コイル感度分布とを前記リードアウト方向ごとに合成することにより、前記リードアウト方向ごとの合成感度分布を生成し、
    前記合成感度分布と前記磁気共鳴データとに基づいて、磁気共鳴画像を生成する、
    ことを備える画像生成方法。
  6. コンピュータに、
    第1のリードアウト方向と前記第1のリードアウト方向に交差する第2のリードアウト方向とを含む複数のリードアウト方向について収集された磁気共鳴データと、前記複数のリードアウト方向に対応しローパスフィルタの感度の分布を示す複数のフィルタ感度分布と、前記磁気共鳴データの収集に用いられる複数のコイルエレメントに対応する複数のコイル感度分布とを取得し、
    前記フィルタ感度分布と前記コイル感度分布とを前記リードアウト方向ごとに合成することにより、前記リードアウト方向ごとの合成感度分布を生成し、
    前記合成感度分布と前記磁気共鳴データとに基づいて、磁気共鳴画像を生成すること、
    を実現させる画像生成プログラム。
JP2020204769A 2020-12-10 2020-12-10 画像生成装置、画像生成方法、および画像生成プログラム Pending JP2022092150A (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2020204769A JP2022092150A (ja) 2020-12-10 2020-12-10 画像生成装置、画像生成方法、および画像生成プログラム
US17/496,199 US11896361B2 (en) 2020-12-10 2021-10-07 Image generating apparatus, image generation method, and nontransitory computer-readable storage medium storing therein image generation program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2020204769A JP2022092150A (ja) 2020-12-10 2020-12-10 画像生成装置、画像生成方法、および画像生成プログラム

Publications (1)

Publication Number Publication Date
JP2022092150A true JP2022092150A (ja) 2022-06-22

Family

ID=81942193

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020204769A Pending JP2022092150A (ja) 2020-12-10 2020-12-10 画像生成装置、画像生成方法、および画像生成プログラム

Country Status (2)

Country Link
US (1) US11896361B2 (ja)
JP (1) JP2022092150A (ja)

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6777934B2 (en) * 1999-12-08 2004-08-17 Hitachi Medical Corporation Magnetic resonance imaging method and apparatus
JP4152381B2 (ja) * 2002-04-30 2008-09-17 株式会社日立メディコ 磁気共鳴イメージング装置
US7408345B2 (en) * 2006-02-06 2008-08-05 The Board Of Trustees Of The Leland Stanford Junior University Generalized MRI reconstruction with correction for multiple image distortion
JP5048052B2 (ja) * 2007-03-09 2012-10-17 株式会社日立メディコ 磁気共鳴イメージング装置および磁気共鳴イメージング方法
US9594139B2 (en) * 2014-03-28 2017-03-14 General Electric Company Systems and methods for magnetic resonance bone depiction
JP7267752B2 (ja) 2019-01-21 2023-05-02 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置及びk空間軌跡の補正方法

Also Published As

Publication number Publication date
US20220183582A1 (en) 2022-06-16
US11896361B2 (en) 2024-02-13

Similar Documents

Publication Publication Date Title
US9188655B2 (en) Magnetic resonance imaging apparatus and image processing apparatus
JP6073627B2 (ja) 磁気共鳴イメージング装置及び画像処理装置
JP6513336B2 (ja) 磁気共鳴イメージング装置及び画像処理装置
JP6494986B2 (ja) 磁気共鳴イメージング装置及び画像処理装置
JP6548531B2 (ja) 磁気共鳴イメージング装置及び画像処理装置
US9274197B2 (en) Magnetic resonance imaging data sampling methods and systems
JP2023171516A (ja) 医用情報処理装置及び医用情報処理方法
US11016154B2 (en) Magnetic resonance imaging apparatus and image reconstruction program
US20240138700A1 (en) Medical image processing apparatus, method of medical image processing, and nonvolatile computer readable storage medium storing therein medical image processing program
JP7236220B2 (ja) 磁気共鳴イメージング装置
KR101471979B1 (ko) 자기 공명(MR; Magnetic Resonance) 영상의 복셀(Voxel)에 대한 MR 스펙트럼을 획득하는 방법 및 장치
US9618595B2 (en) Magnetic resonance imaging system, data processing apparatus, and method for generating magnetic resonance image
JP2022092150A (ja) 画像生成装置、画像生成方法、および画像生成プログラム
JP7419145B2 (ja) 医用情報処理装置及び医用情報処理方法
JP2022092151A (ja) 画像生成装置、画像生成方法、および画像生成プログラム
CN111009020B (zh) 图像重建方法、装置、计算机设备和存储介质
US10845442B2 (en) Magnetic resonance imaging apparatus, processing apparatus and medical image processing method
WO2017167937A1 (en) Dynamic mr imaging with increased temporal and spatial resolution
JP2022092161A (ja) 磁気共鳴イメージング装置、および画像生成方法
JP5360757B2 (ja) 磁気共鳴画像法および磁気共鳴画像装置
US12019132B2 (en) Magnetic resonance imaging apparatus, image generation method, and computer-readable nonvolatile storage medium storing image generation program
JP7451366B2 (ja) 画像生成装置、画像生成方法、および画像生成プログラム
US20220309669A1 (en) Medical image diagnosis system and non-transitory computer-readable storage medium storing therein image processing program
JP7510840B2 (ja) 情報処理装置、情報処理方法、および情報処理プログラム
JP2005168868A (ja) 磁気共鳴イメージング装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20231030