JP2023009819A - 撮影装置 - Google Patents
撮影装置 Download PDFInfo
- Publication number
- JP2023009819A JP2023009819A JP2021113419A JP2021113419A JP2023009819A JP 2023009819 A JP2023009819 A JP 2023009819A JP 2021113419 A JP2021113419 A JP 2021113419A JP 2021113419 A JP2021113419 A JP 2021113419A JP 2023009819 A JP2023009819 A JP 2023009819A
- Authority
- JP
- Japan
- Prior art keywords
- vibration
- magnetic field
- image
- region
- size
- 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
- 238000003384 imaging method Methods 0.000 title abstract description 12
- 238000002091 elastography Methods 0.000 claims abstract description 9
- 230000001902 propagating effect Effects 0.000 claims description 20
- 230000003068 static effect Effects 0.000 claims description 14
- 230000000644 propagated effect Effects 0.000 abstract 2
- 230000006866 deterioration Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 52
- 230000001427 coherent effect Effects 0.000 description 35
- 238000006073 displacement reaction Methods 0.000 description 30
- 238000002595 magnetic resonance imaging Methods 0.000 description 20
- 238000012045 magnetic resonance elastography Methods 0.000 description 19
- 238000000034 method Methods 0.000 description 15
- 239000013598 vector Substances 0.000 description 10
- 238000012545 processing Methods 0.000 description 9
- 238000004364 calculation method Methods 0.000 description 7
- 235000019589 hardness Nutrition 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 230000000052 comparative effect Effects 0.000 description 6
- 238000005259 measurement Methods 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 210000003205 muscle Anatomy 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 206010073306 Exposure to radiation Diseases 0.000 description 1
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000000090 biomarker Substances 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
-
- 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/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- General Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Biomedical Technology (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Signal Processing (AREA)
- Biophysics (AREA)
- Chemical & Material Sciences (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
Description
本発明は、撮影装置に関し、特に、組織の硬さを測定するMRエラストグラフィ(Magnetic Resonance Elastography : MRE)や超音波エラストグラフィの撮影が可能な撮影装置に関する。
医療現場において、患者に対して放射線被曝等の影響がないMRI装置(Magnetic Resonance Imaging装置:磁気共鳴画像装置)が使用されている。MRI装置では、人体の各細胞に含まれる水素原子核(プロトン)に対して、プロトンのスピンに応じた高周波磁場を印加して励起し、励起したプロトンが元の状態に戻る(緩和)際に発する電磁波に基づいて、プロトンの密度が異なる部分(例えば、水と脂肪)を濃淡で表した画像を得ることが可能である。
MRE(Magnetic Resonance Elastography:磁気共鳴エラストグラフィ)は、対象物に振動(体幹部の場合には、50Hz~100Hz程度)を加えながら、MRI装置で撮像することで、対象部内部の「硬さ」の違いによる振動波の伝播の違いを利用し、硬さを画像化する撮像法である(特許文献1,2参照)。
MRE(Magnetic Resonance Elastography:磁気共鳴エラストグラフィ)は、対象物に振動(体幹部の場合には、50Hz~100Hz程度)を加えながら、MRI装置で撮像することで、対象部内部の「硬さ」の違いによる振動波の伝播の違いを利用し、硬さを画像化する撮像法である(特許文献1,2参照)。
MRエラストグラフィにおいて、体内を振動波が伝播する場合、振動源から付与される振動の進行方向に対して、硬さの異なる臓器や骨、筋肉、体表等で伝播波が反射したり散乱したりして進行方向が変わることがある。したがって、計測される伝播波は、さまざまな方向の成分を含むこととなる。主方向の伝播波と、反射波、散乱波が干渉すると、伝播波画像(Wave Image:WI)の局所波長が変化する、すなわち、誤差が発生する可能性がある。この変化(誤差)が弾性率マップ(elastogram)の精度を低下させる恐れがある。
このような反射波、散乱波に対して、これまでにも、Directional Filter(DF)と呼ばれる画像フィルタを使用して、反射等で進行方向が変わった伝播波を除去することが行われている(非特許文献1参照)。非特許文献1に記載の技術では、主方向の伝播波(非特許文献1では上部から下部に向かう伝播波)に対してDFを適用することで、反射波(非特許文献1では下部から上部に向かう波、もしくは波紋状に上方向に広がる波)を除去している。
A.Manduca et al.,"Spatio-temporal directional filtering forimproved inversion of MR elastography images", Medical Image Analysis 7 (2003)465-473
(従来技術の問題点)
非特許文献1に記載のDF(Directional Filter)では、特定の一方向の波以外を除去する技術であるため、波が特定の一方向に伝播することが前提となっている。ファントム(性能評価用の試料)のような理想的な状態では、非特許文献1に記載の技術は適用可能であるが、生体内は、臓器や骨等のため、伝播する波の主方向が変化することがほとんどである。逆に言えば、生体内で伝播波が一方向だけに進むことは、ほぼあり得ない。
したがって、非特許文献1に記載の技術では、生体内を伝播する波に対して、途中で主方向が変化した波まで除去してしまう恐れが高い問題がある。主方向の波が除去されると、MRE画像の精度が低下したり、MRE画像が形成できない問題がある。
非特許文献1に記載のDF(Directional Filter)では、特定の一方向の波以外を除去する技術であるため、波が特定の一方向に伝播することが前提となっている。ファントム(性能評価用の試料)のような理想的な状態では、非特許文献1に記載の技術は適用可能であるが、生体内は、臓器や骨等のため、伝播する波の主方向が変化することがほとんどである。逆に言えば、生体内で伝播波が一方向だけに進むことは、ほぼあり得ない。
したがって、非特許文献1に記載の技術では、生体内を伝播する波に対して、途中で主方向が変化した波まで除去してしまう恐れが高い問題がある。主方向の波が除去されると、MRE画像の精度が低下したり、MRE画像が形成できない問題がある。
なお、MREではなくMRI画像やCT画像において、撮像パラメータや画像処理パラメータなどに国際的な基準がある(QIBA:Quantitative Imaging Biomarkers Alliance)。MREにおいても同様な枠組みが形成されつつあり、「伝播波に反射や散乱が含まれていない一方向に伝播する領域(コヒーレントな伝播波領域)を選択して硬さを測定する」というものがテーマとして挙げられている。現状では、コヒーレントな伝播波領域の判断は、観察者の主観によって決定されており、客観的な領域設定がなされていない。
本発明は、利用者の主観で領域設定を行う場合に比べて、領域設定の精度および再現性の低下を解決することを技術的課題とする。
前記技術的課題を解決するために、請求項1に記載の発明の撮影装置は、
静磁場と、位置に応じて磁場が変化する傾斜磁場と、プロトンの磁気共鳴条件に基づいて予め設定された交番磁場と、を被検者の被検査部に対して発生させる磁場発生装置と、
被検査部に振動を付与する振動付与部材と、
前記傾斜磁場を印加する時期に応じて取得された電磁波に基づいて、電磁波信号の位相に応じたMR位相画像を取得するMR位相画像の取得手段であって、予め定められた時間間隔をあけて複数の前記MR位相画像を取得する前記MR位相画像の取得手段と、
MR位相画像の全画像領域の中から区分けされた判別対象領域において、伝播する振動の波長が長くなると、前記判別対象領域のサイズを大きくする領域サイズ調整手段と、
前記領域サイズ調整手段で調整された判別対象領域に基づいて、伝播する振動の主方向が揃っているか否かを判別する判別手段と、
前記判別手段で振動の主方向が揃っていない判別対象領域の波が除去されたMR位相画像に基づいて、MRエラストグラフィ画像を作成するMRE画像作成手段と、
を備えたことを特徴とする。
静磁場と、位置に応じて磁場が変化する傾斜磁場と、プロトンの磁気共鳴条件に基づいて予め設定された交番磁場と、を被検者の被検査部に対して発生させる磁場発生装置と、
被検査部に振動を付与する振動付与部材と、
前記傾斜磁場を印加する時期に応じて取得された電磁波に基づいて、電磁波信号の位相に応じたMR位相画像を取得するMR位相画像の取得手段であって、予め定められた時間間隔をあけて複数の前記MR位相画像を取得する前記MR位相画像の取得手段と、
MR位相画像の全画像領域の中から区分けされた判別対象領域において、伝播する振動の波長が長くなると、前記判別対象領域のサイズを大きくする領域サイズ調整手段と、
前記領域サイズ調整手段で調整された判別対象領域に基づいて、伝播する振動の主方向が揃っているか否かを判別する判別手段と、
前記判別手段で振動の主方向が揃っていない判別対象領域の波が除去されたMR位相画像に基づいて、MRエラストグラフィ画像を作成するMRE画像作成手段と、
を備えたことを特徴とする。
請求項2に記載の発明は、請求項1に記載の撮影装置において、
前記振動付与部材で付与される振動の伝播波長が大きくなるにつれて、前記判別対象領域の大きさを大きくする前記領域サイズ調整手段、
を備えたことを特徴とする。
前記振動付与部材で付与される振動の伝播波長が大きくなるにつれて、前記判別対象領域の大きさを大きくする前記領域サイズ調整手段、
を備えたことを特徴とする。
請求項1に記載の発明によれば、利用者の主観でコヒーレント領域を判別する場合に比べて精度よく判別可能であり、弾性率計測の精度および再現性の向上を可能にすることができる。
請求項2に記載の発明によれば、振動付与部材で付与される振動の伝播波長に応じて、判別対象領域の大きさを自動的に調整できる。
請求項2に記載の発明によれば、振動付与部材で付与される振動の伝播波長に応じて、判別対象領域の大きさを自動的に調整できる。
次に図面を参照しながら、本発明の実施の形態の具体例(以下、実施例と記載する)を説明するが、本発明は以下の実施例に限定されるものではない。
なお、以下の図面を使用した説明において、理解の容易のために説明に必要な部材以外の図示は適宜省略されている。
なお、以下の図面を使用した説明において、理解の容易のために説明に必要な部材以外の図示は適宜省略されている。
図1は本発明の実施例1の磁気共鳴撮影装置の説明図である。
図1において、本発明の撮影装置の一例としての実施例1の磁気共鳴撮影装置1は、磁場発生装置の一例としての磁石部2を有する。磁石部2には、内部を水平方向に貫通する貫通孔3が形成されている。貫通孔3には、寝た状態の被検者4が支持される寝台6が貫通可能である。
磁石部2は、静磁場印加部材の一例としての静磁場発生磁石11を有する。なお、静磁場発生磁石として、超電導電磁石や永久磁石を使用することが可能である。静磁場発生磁石11の内側には、傾斜磁場印加部材の一例としての傾斜磁場発生コイル12が配置されている。傾斜磁場発生コイル12の内側には、励起磁場印加部材の一例としての高周波磁場発生コイル13が配置されている。高周波磁場発生コイル13の内側には、受信部の一例として、電磁波を受信する受信コイル14が配置されている。
図1において、本発明の撮影装置の一例としての実施例1の磁気共鳴撮影装置1は、磁場発生装置の一例としての磁石部2を有する。磁石部2には、内部を水平方向に貫通する貫通孔3が形成されている。貫通孔3には、寝た状態の被検者4が支持される寝台6が貫通可能である。
磁石部2は、静磁場印加部材の一例としての静磁場発生磁石11を有する。なお、静磁場発生磁石として、超電導電磁石や永久磁石を使用することが可能である。静磁場発生磁石11の内側には、傾斜磁場印加部材の一例としての傾斜磁場発生コイル12が配置されている。傾斜磁場発生コイル12の内側には、励起磁場印加部材の一例としての高周波磁場発生コイル13が配置されている。高周波磁場発生コイル13の内側には、受信部の一例として、電磁波を受信する受信コイル14が配置されている。
実施例1では、被検者4には、被検査部の一例としての大腰筋のMRE測定をするために、背骨側の体表面の部分に、振動付与部材の一例としての振動パッド16が支持されている。
前記磁石部2には、情報処理装置の一例としてのコンピュータ装置21がケーブルCbを介して電気的に接続されている。したがって、コンピュータ装置21は、磁石部2との間で、静磁場発生磁石11等の制御信号や受信コイル14での検知信号等が送受信可能に構成されている。コンピュータ装置21は、コンピュータ本体22と、表示部の一例としてのディスプレイ23と、入力部の一例としてのキーボード24およびマウス25と、を有する。なお、実施例1では、コンピュータ装置21と磁石部2とをケーブルCbで接続する構成を例示したが、これに限定されず、携帯電話回線やBluetooth(登録商標)、無線LAN等、任意の無線通信方式で情報の送受信を行うことも可能である。
(実施例1のコンピュータ本体22の制御部の説明)
図2は実施例1の磁気共鳴撮影装置におけるコンピュータ本体の機能ブロック図である。
図2において、実施例1のコンピュータ本体22の制御部41は、外部との信号の入出力および入出力信号レベルの調節等を行うI/O(入出力インターフェース)、必要な起動処理を行うためのプログラムおよびデータ等が記憶されたROM(リードオンリーメモリ)、必要なデータ及びプログラムを一時的に記憶するためのRAM(ランダムアクセスメモリ)、ROM等に記憶された起動プログラムに応じた処理を行うCPU(中央演算処理装置)ならびにクロック発振器等を有するコンピュータ装置により構成されており、前記ROM及びRAM等に記憶されたプログラムを実行することにより種々の機能を実現することができる。
制御部41には、基本動作を制御する基本ソフト、いわゆる、オペレーティングシステムOS、アプリケーションプログラムの一例としての撮影装置制御プログラムAP1、その他の図示しないソフトウェアが記憶されている。
図2は実施例1の磁気共鳴撮影装置におけるコンピュータ本体の機能ブロック図である。
図2において、実施例1のコンピュータ本体22の制御部41は、外部との信号の入出力および入出力信号レベルの調節等を行うI/O(入出力インターフェース)、必要な起動処理を行うためのプログラムおよびデータ等が記憶されたROM(リードオンリーメモリ)、必要なデータ及びプログラムを一時的に記憶するためのRAM(ランダムアクセスメモリ)、ROM等に記憶された起動プログラムに応じた処理を行うCPU(中央演算処理装置)ならびにクロック発振器等を有するコンピュータ装置により構成されており、前記ROM及びRAM等に記憶されたプログラムを実行することにより種々の機能を実現することができる。
制御部41には、基本動作を制御する基本ソフト、いわゆる、オペレーティングシステムOS、アプリケーションプログラムの一例としての撮影装置制御プログラムAP1、その他の図示しないソフトウェアが記憶されている。
(実施例1の制御部41に接続された要素)
制御部41には、キーボード24やマウス25、受信コイル14等の信号出力要素からの出力信号が入力されている。
また、実施例1の制御部41は、ディスプレイ23、静磁場発生磁石11、傾斜磁場発生コイル12、高周波磁場発生コイル13等の被制御要素へ制御信号を出力している。
制御部41には、キーボード24やマウス25、受信コイル14等の信号出力要素からの出力信号が入力されている。
また、実施例1の制御部41は、ディスプレイ23、静磁場発生磁石11、傾斜磁場発生コイル12、高周波磁場発生コイル13等の被制御要素へ制御信号を出力している。
(制御部41の機能)
実施例1の制御部41の撮影装置制御プログラムAP1は、下記の機能手段(プログラムモジュール)51~55を有する。
実施例1の制御部41の撮影装置制御プログラムAP1は、下記の機能手段(プログラムモジュール)51~55を有する。
磁場制御手段51は、磁石部2を制御して、被検者4の被検査部をMR撮影するための磁場を制御する。実施例1の磁場制御手段51は、繰り返し時間記憶手段51aと、エコー時間記憶手段51bと、静磁場印加手段51cと、傾斜磁場印加手段51dと、交番磁場の印加手段の一例としての高周波磁場印加手段51eと、を有する。
図3は実施例1の磁場の印加および振動の付与の説明図であり、横軸に時間を取ったグラフである。
繰り返し時間記憶手段51aは、被検者4の被検査部に含まれるプロトンを励起するために印加される交番磁場の一例としての高周波磁場を印加する間隔である繰り返し時間TRを記憶する。
繰り返し時間記憶手段51aは、被検者4の被検査部に含まれるプロトンを励起するために印加される交番磁場の一例としての高周波磁場を印加する間隔である繰り返し時間TRを記憶する。
エコー時間記憶手段51bは、高周波磁場が印加されてから、励起されたプロトンが元の状態に戻る(緩和する)際に発する電磁波を取得するまでの間隔であるエコー時間TEを記憶する。なお、実施例1では、繰り返し時間TRおよびエコー時間TEは、予め設定されているが、磁気共鳴撮影装置1の利用者が手動で入力して、設定、変更が可能に構成することも可能である。
静磁場印加手段51cは、静磁場発生磁石11を制御して、静磁場を発生させる。実施例1の静磁場印加手段51cは、一例として、3[T]の静磁場を発生させる。
傾斜磁場印加手段51dは、傾斜磁場発生コイル12を制御して、位置に応じて磁場が変化するMRE用の傾斜磁場(勾配磁場)を発生させる。従って、傾斜磁場が振動検出傾斜磁場MEG(motion encoding gradient)と呼ばれる磁場である。実施例1の傾斜磁場印加手段51dは、図3に示すように、互いに直交するスライス(slice)方向、リードアウト(read out)方向およびフェーズ(phase)方向の3軸方向において、スライス方向(スライス軸)に傾斜磁場を発生させる。
高周波磁場印加手段51eは、高周波磁場発生コイル13を制御して、プロトンを励起する周波数に対応する交番磁場である高周波磁場を発生させる。
傾斜磁場印加手段51dは、傾斜磁場発生コイル12を制御して、位置に応じて磁場が変化するMRE用の傾斜磁場(勾配磁場)を発生させる。従って、傾斜磁場が振動検出傾斜磁場MEG(motion encoding gradient)と呼ばれる磁場である。実施例1の傾斜磁場印加手段51dは、図3に示すように、互いに直交するスライス(slice)方向、リードアウト(read out)方向およびフェーズ(phase)方向の3軸方向において、スライス方向(スライス軸)に傾斜磁場を発生させる。
高周波磁場印加手段51eは、高周波磁場発生コイル13を制御して、プロトンを励起する周波数に対応する交番磁場である高周波磁場を発生させる。
振動付与制御手段52は、振動パッド16により被検査部に振動を付与する。
受信手段の一例としての信号取得手段53は、エコー時間TEの時期に、受信コイル14を介して被検者4のプロトンが緩和する際に発生する電磁波信号を取得する。したがって、実施例1では、図3に示す振動が振動付与制御手段52で付与された状態で、エコー時間TEにおいて受信コイル14で信号を測定することで、振動(Vibration)を付与しながらのMR画像を撮像する。
受信手段の一例としての信号取得手段53は、エコー時間TEの時期に、受信コイル14を介して被検者4のプロトンが緩和する際に発生する電磁波信号を取得する。したがって、実施例1では、図3に示す振動が振動付与制御手段52で付与された状態で、エコー時間TEにおいて受信コイル14で信号を測定することで、振動(Vibration)を付与しながらのMR画像を撮像する。
MR画像取得手段(MRE画像作成手段)54は、信号取得手段53が取得した電磁波信号に基づいて、MR画像を作成する。
信号処理手段54aは、受信した電磁波信号において信号処理をする。実施例1の信号処理手段54aは、信号取得手段53が取得した信号を実数部r(real part)とし、信号取得手段53が取得した信号をπ/2位相を遅らせた信号を虚数部i(imaginary part)とする。すなわち、受信した電磁波信号に基づいた複素数、いわゆる、MRIの技術分野におけるk空間(周波数空間)の信号を生成する。そして、実数部rと虚数部iに対して、フーリエ逆変換(実施例では高速フーリエ逆変換)を行って、実空間の信号R,Iに変換する。そして、実空間における実数部Rと虚数部Iとに基づいて、複素平面における強度M=(R2+I2)1/2と、位相φ=tan-1(I/R)とを演算する。なお、この演算は、従来のMRI装置において導入されており、公知であるため、これ以上の詳細な説明は省略する。
信号処理手段54aは、受信した電磁波信号において信号処理をする。実施例1の信号処理手段54aは、信号取得手段53が取得した信号を実数部r(real part)とし、信号取得手段53が取得した信号をπ/2位相を遅らせた信号を虚数部i(imaginary part)とする。すなわち、受信した電磁波信号に基づいた複素数、いわゆる、MRIの技術分野におけるk空間(周波数空間)の信号を生成する。そして、実数部rと虚数部iに対して、フーリエ逆変換(実施例では高速フーリエ逆変換)を行って、実空間の信号R,Iに変換する。そして、実空間における実数部Rと虚数部Iとに基づいて、複素平面における強度M=(R2+I2)1/2と、位相φ=tan-1(I/R)とを演算する。なお、この演算は、従来のMRI装置において導入されており、公知であるため、これ以上の詳細な説明は省略する。
MR強度画像の作成手段(MR強度画像の取得手段)54bは、信号処理手段54aで算出された強度Mに基づいて、MR強度画像を作成する。なお、MR強度画像は、一般的にMRI画像として、診断に使用される画像である。
MR位相画像の作成手段(MR位相画像の取得手段)54cは、電磁波信号から算出された位相φに基づいて、MR位相画像(MREではWave Imageとして利用される)を作成する。ここで、実施例1のMR位相画像の作成手段54cでは、エコー時間TEで測定された電磁波信号から算出された位相φに基づいて、MR位相画像を作成する。
判別方向設定手段54dは、各ボクセル(最小画素)における主方向を判別するために、主方向の候補としての判別方向を設定する。実施例1では、判別方向は、一例として、8方向(図4Bにおける上から下、下から上、右から左、左から右、右上から左下、左下から右上、左上から右下、右下から左上)が順に設定される。したがって、まず、上から下の方向が設定されて、主方向の判別が行われると、次の方向が設定されて判別が行われるという処理が8つの方向に対して順に行われる。なお、実施例1では、一例として判別方向が8つの場合を例示したが、これに限定されず、1度ずつずらした360度分の方向とする等、要求される精度や計算処理速度等に応じて任意に変更可能である。
図4は判別方向と変位の読取方向の説明図であり、図4Aは判別方向と読取方向が同一の場合の説明図、図4Bは判別方向と読取方向が逆方向の場合の説明図である。
読取方向設定手段54eは、各ボクセルにおいて設定された判別方向111に沿って並んだ複数の位置での振動の変位を取得する際に、変位を読み取る読取方向112(図4B参照)を設定する。実施例1では、図4Bに示すように、判別方向111とは逆方向に読取方向112が設定されている。一例として、判別方向111が、図4Bにおける「上から下」の場合、判別方向に沿って並んだ3つの点103a~103cの各位置での振動の変位を取得する際に、読取方向112は「下から上」となり、103c,103b,103aの順に振動の変位が読み取られるように、読取方向112が設定される。
読取方向設定手段54eは、各ボクセルにおいて設定された判別方向111に沿って並んだ複数の位置での振動の変位を取得する際に、変位を読み取る読取方向112(図4B参照)を設定する。実施例1では、図4Bに示すように、判別方向111とは逆方向に読取方向112が設定されている。一例として、判別方向111が、図4Bにおける「上から下」の場合、判別方向に沿って並んだ3つの点103a~103cの各位置での振動の変位を取得する際に、読取方向112は「下から上」となり、103c,103b,103aの順に振動の変位が読み取られるように、読取方向112が設定される。
変位取得手段54fは、MR位相画像の各ボクセルにおいて、設定された判別方向に沿って3つの点103a~103cを抽出して、各点(各位置)における振動の変位を取得する。この時、変位の取得は、読取方向に沿って順に取得する。また、実施例1の変位取得手段54fでは、時間ごとに変化していく各時期のMR位相画像に対して、振動パッド16の振動の周期Tの1/4の時間間隔をあけて、各点103a~103cの変位を取得するように設定されている。具体的には、実施例1では、時刻(位相)0,T/4,T/2,3T/4における各点103a~103cの変位について、点103cの時刻0の変位、点103bの時刻0の変位、点103aの時刻0の変位、点103cの時刻T/4の変位、点103bの時刻T/4の変位、点103aの時刻T/4の変位、点103cの時刻T/2の変位、点103bの時刻T/2の変位、点103aの時刻T/2の変位、点103cの時刻3T/4の変位、点103bの時刻3T/4の変位、点103aの時刻3T/4の変位の順に取得する。なお、実施例1では、各点103a~103cを等間隔の点としたが、これに限定されない。すなわち、後述する正弦波への近さの判定では、各点103a~103cの間隔が等間隔である方が誤差は少なくなりやすいが、判定自体については等間隔でなくても可能であり、等間隔に制限されない。また、3点に限定されず、4点以上とすることも可能である。
波形生成手段54gは、変位取得手段54fで取得された変位に基づいて、判別波形113(図5参照)を生成する。
図5は残差の演算方法の説明図である。
残差算出手段54hは、判別波形113が正弦波に対する近さを導出するために、振動パッド16で付与される振動である正弦波との差である残差を算出する。なお、実施例1では、振動パッド16で正弦波状の振動が付与されるが、付与される振動が正弦波状ではない場合は、付与される振動(f(t))との残差を導出することが可能である。
実施例1では、正弦波との残差Rを算出するために、波の変位量をuとし、振幅(最大の変位量)をu0とし、波数ベクトルをkとし、位置ベクトルをrとし、振動の各振動数をωとし、時間をtとし、初期位相をθとし、磁場の不均一性に起因する成分をCとした場合に、以下の式(1)を考える。
u(r,t)=u0・sin(kr+ωt+θ)+C …式(1)
残差算出手段54hは、判別波形113が正弦波に対する近さを導出するために、振動パッド16で付与される振動である正弦波との差である残差を算出する。なお、実施例1では、振動パッド16で正弦波状の振動が付与されるが、付与される振動が正弦波状ではない場合は、付与される振動(f(t))との残差を導出することが可能である。
実施例1では、正弦波との残差Rを算出するために、波の変位量をuとし、振幅(最大の変位量)をu0とし、波数ベクトルをkとし、位置ベクトルをrとし、振動の各振動数をωとし、時間をtとし、初期位相をθとし、磁場の不均一性に起因する成分をCとした場合に、以下の式(1)を考える。
u(r,t)=u0・sin(kr+ωt+θ)+C …式(1)
この式において、未知の定数は、u0、k、θ、Cである。kの値を決めると、他の変数は正弦波の線形最小二乗近似で推定可能であり、u(r,t)の近似解u′(r,t)が求まる。u(r,t)とu′(r,t)との残差が最小となるkを推定するために以下の手順(1)~(3)の処理を行う。
(1)適当なkの初期値を与える。
(2)最小二乗近似により近似解u′(r,t)を算出し、実測値u(r,t)、すなわち、判別波形113との残差Rを算出する。
(3)最急降下法により残差Rが最小になるまでkを更新する。
(1)~(3)の手順で得られた最小の残差Rを、8つの判別方向111毎に算出する。
(1)適当なkの初期値を与える。
(2)最小二乗近似により近似解u′(r,t)を算出し、実測値u(r,t)、すなわち、判別波形113との残差Rを算出する。
(3)最急降下法により残差Rが最小になるまでkを更新する。
(1)~(3)の手順で得られた最小の残差Rを、8つの判別方向111毎に算出する。
主方向判別手段54iは、8つの判別方向111における判別波形113の中で最も正弦波に近い判別波形113が得られた判別方向111を主方向と判別する。実施例1では、残差算出手段54hで得られた残差Rを、プロファイル方向(判別方向)を180度変えた残差Rと比較し、残差Rが小さい方向を記録する。これを、全判別方向(360度分、実施例1では8方向分)行う。よって、8つの方向で行ったとすると、180度変える前と後の2方向の比較が4セット繰り返されることとなる。それぞれで残差Rが小さい方向が記録されるので、4方向のベクトルとなる。それぞれのベクトルの大きさを同じ(例えば「1」)に正規化して、ベクトルの和を算出する。この合成ベクトルが各ボクセルの波の進行方向(主方向)となる。
図6は実施例1のコヒーレント領域を判別するために領域を区分けする説明図であり、図6Aは小領域(カーネル)に区分けする説明図、図6Bは画素(ボクセル)とカーネルの説明図、図6Cは図6Aの全域において各画素の局所波長を画素値で示した一例の図である。
局所波長推定手段54jは、MR位相画像の作成手段54cで作成されたMR位相画像に基づいて、MR位相画像内の各所における伝播波長の分布を導出する。前記主方向判別手段54iで判別された各ボクセルの主方向にDFを適用し、反射や散乱の影響を取り除いた状態で、伝播波の局所位相変化から局所波長を推定する。この、DFの適用と局所波長の推定を、全ボクセル分繰り返す。
ただし、局所波長の推定方法は上記に限定されず、推定精度が担保されている方法であれば代用することが可能である。
局所波長推定手段54jは、MR位相画像の作成手段54cで作成されたMR位相画像に基づいて、MR位相画像内の各所における伝播波長の分布を導出する。前記主方向判別手段54iで判別された各ボクセルの主方向にDFを適用し、反射や散乱の影響を取り除いた状態で、伝播波の局所位相変化から局所波長を推定する。この、DFの適用と局所波長の推定を、全ボクセル分繰り返す。
ただし、局所波長の推定方法は上記に限定されず、推定精度が担保されている方法であれば代用することが可能である。
領域区分け手段54kは、MR位相画像の全域101を、予め定められた広さの領域の一例として小領域(カーネル)103に区分けする。図6Aにおいて、小領域(カーネル)103のサイズは、MR位相画像の全域101の大きさや分解能等の、仕様に応じて任意に変更可能であるが、少なくとも3画素×3画素以上の広さが必要である。
領域サイズ調整手段54mは、各ボクセル102において、コヒーレント領域の判別を行うための判別対象領域(カーネル)103の大きさを調整する。図6B、図6Cにおいて、実施例1の領域サイズ調整手段54mは、カーネル103のサイズを、ボクセル102を伝播する振動の波長に応じて調整する。実施例1では、波長が大きくなるほど、カーネル103のサイズが大きくなるように調整する。具体的には、振動付与制御手段52で付与される振動の伝播波長(局所波長)が、予め定められた閾値よりも高い場合(図6Cで黒色の領域101b、弾性率の高い領域)は、カーネル103のサイズを大きくし、前記閾値よりも低い場合(図6Cで白色の領域101a、弾性率の低い領域)はカーネル103のサイズを小さくする。
一例として、長波長の場合はカーネル103を8画素×8画素とし、低波長の場合はカーネル103を4画素×4画素とすることが可能であるが、具体的な数値は例示した数値に限定されず、要求される精度や仕様、利用者の設定等に応じて任意に変更可能である。
また、伝播波長と閾値を使用して判別を行う構成を例示したがこれに限定されない。例えば、伝播波長に限らず、弾性率のように「伝播波長に関連するパラメータ」を使用することも可能である。
一例として、長波長の場合はカーネル103を8画素×8画素とし、低波長の場合はカーネル103を4画素×4画素とすることが可能であるが、具体的な数値は例示した数値に限定されず、要求される精度や仕様、利用者の設定等に応じて任意に変更可能である。
また、伝播波長と閾値を使用して判別を行う構成を例示したがこれに限定されない。例えば、伝播波長に限らず、弾性率のように「伝播波長に関連するパラメータ」を使用することも可能である。
図7はWave-direction MapとCoherence Mapであり、図7AはWave-direction Map、図7Bは図6Cに対応する局所波長Map、図7CはCoherence Mapである。
コヒーレント領域判別手段54nは、領域サイズ調整手段54mで調整されたカーネル103に基づいて、各カーネル103がコヒーレントな領域であるか否かを判別する。実施例1のコヒーレント領域判別手段54nは、まず、MR位相画像(Wave Image画像)から主方向を算出した画像(Wave-direction map:図7A参照)に対し、ボクセル毎に局所波長を算出して(図7B参照)、カーネル103のサイズを調整し、カーネル103内のベクトル平均(主方向のベクトルの大きさを1と仮定)を算出する(Coherence map:図7C参照)。そして、コヒーレント領域判別手段54nは、予め定められた閾値に基づいて、各カーネル103のベクトル平均の値が1に近ければ(閾値以上であれば)コヒーレント、0に近ければ(閾値未満であれば)インコヒーレントと判定する。すなわち、1に近ければ、その領域を伝播する振動はコヒーレント、すなわち、振動の位相(振動の進行方向)が揃っていて、0に近ければ、その領域を伝播する振動はインコヒーレント、すなわち、振動の進行方向が揃っていない、と判別される。
コヒーレント領域判別手段54nは、領域サイズ調整手段54mで調整されたカーネル103に基づいて、各カーネル103がコヒーレントな領域であるか否かを判別する。実施例1のコヒーレント領域判別手段54nは、まず、MR位相画像(Wave Image画像)から主方向を算出した画像(Wave-direction map:図7A参照)に対し、ボクセル毎に局所波長を算出して(図7B参照)、カーネル103のサイズを調整し、カーネル103内のベクトル平均(主方向のベクトルの大きさを1と仮定)を算出する(Coherence map:図7C参照)。そして、コヒーレント領域判別手段54nは、予め定められた閾値に基づいて、各カーネル103のベクトル平均の値が1に近ければ(閾値以上であれば)コヒーレント、0に近ければ(閾値未満であれば)インコヒーレントと判定する。すなわち、1に近ければ、その領域を伝播する振動はコヒーレント、すなわち、振動の位相(振動の進行方向)が揃っていて、0に近ければ、その領域を伝播する振動はインコヒーレント、すなわち、振動の進行方向が揃っていない、と判別される。
除去手段54pは、主方向以外の振動を除去する。実施例1では、コヒーレント領域判別手段54nで判別されたコヒーレントな領域内に関して、ボクセル毎に、主方向判別手段54iで判別された主方向以外の方向についてDF:Directional Filterで振動の成分を除去する。これにより、弾性率算出に誤差を引き起こす散乱波や反射波の除去が可能となる。なお、DF自体は、公知であるため、詳細な説明は省略する。
したがって、実施例1のMR画像取得手段(MRE画像作成手段)54は、MR位相画像の作成手段54cで作成されたMR位相画像から、除去手段54pでボクセル毎に主方向以外の振動が除去されて弾性率が導出された画像を、最終的なMRE画像として出力(作成)する。
したがって、実施例1のMR画像取得手段(MRE画像作成手段)54は、MR位相画像の作成手段54cで作成されたMR位相画像から、除去手段54pでボクセル毎に主方向以外の振動が除去されて弾性率が導出された画像を、最終的なMRE画像として出力(作成)する。
画像表示手段55は、MR画像取得手段54で作成された画像を、ディスプレイ23に表示する。すなわち、被検査部の断面画像であるMR強度画像(通常の診断で利用する画像)と、MRE画像とがディスプレイ23に表示される。
(実施例1の作用)
前記構成を備えた実施例1の磁気共鳴撮影装置1では、MR位相画像において、領域101をボクセル102に区分けして、各ボクセル102において伝播する振動の波長に応じてサイズが調整されたカーネル103毎に、そのカーネル内がコヒーレントな領域であるかが判別される。生体内では音響インピーダンスが異なる組織が含まれるため反射や散乱が生じやすく、たとえコヒーレントな領域であっても硬さの算出に誤差が生じる場合がある。この問題に対し、判別された主方向にあわせたDFを適用することで、反射や散乱を抑制すると共に、主方向以外の振動成分を除去できる。したがって、反射波や散乱波等の主方向以外の波が除去されない場合に比べて、MRE画像において各部位の硬さ等の算出を精度よく行うことができ、診断の精度向上にも寄与する。
前記構成を備えた実施例1の磁気共鳴撮影装置1では、MR位相画像において、領域101をボクセル102に区分けして、各ボクセル102において伝播する振動の波長に応じてサイズが調整されたカーネル103毎に、そのカーネル内がコヒーレントな領域であるかが判別される。生体内では音響インピーダンスが異なる組織が含まれるため反射や散乱が生じやすく、たとえコヒーレントな領域であっても硬さの算出に誤差が生じる場合がある。この問題に対し、判別された主方向にあわせたDFを適用することで、反射や散乱を抑制すると共に、主方向以外の振動成分を除去できる。したがって、反射波や散乱波等の主方向以外の波が除去されない場合に比べて、MRE画像において各部位の硬さ等の算出を精度よく行うことができ、診断の精度向上にも寄与する。
また、カーネル103のサイズが大きくなると、対象となる波の分布領域が大きく(さまざまな主方向をもつ波がカーネル内に含まれるように)なり、潜在的にコヒーレンス度が低くなる傾向にある。一方、カーネルサイズが大きくなると、対象となる波の分布領域が小さくなるため、コヒーレンス度が高くなる傾向にある。MREではこれと同様の問題が弾性率マップの空間分解能に影響する。具体的には、波長が短いと、ある領域内に含まれる波の数(ピーク、波の極大と極小の数)が多くなるため、弾性率マップの空間分解能が向上する。一方、波長が長いと、弾性率マップの空間分解能は低下する。反射や散乱などの波の干渉は弾性率算出にエラーを引き起こすが、空間分解能が低い場合、干渉領域だけでなくその周囲においても弾性率算出にエラーが生じてしまう。
したがって、カーネル103のサイズが波長に関わらず固定の場合は、付与される振動の波長や、振動が伝播する部位の弾性率によっては、十分な精度でコヒーレント領域を判定できなかったり、エラー領域が含まれたりして、弾性率計測の精度が低下する恐れがある。
これに対して、実施例1では、振動の波長に応じてカーネル103のサイズが調整されている。したがって、伝播する振動の波長が長い場合に、1つのカーネル103に含まれる波の数を増やすことが可能であり、カーネル103のサイズが固定の場合に比べて、コヒーレント領域の抽出精度を向上し、弾性率計測の精度を向上させることができる。
これに対して、実施例1では、振動の波長に応じてカーネル103のサイズが調整されている。したがって、伝播する振動の波長が長い場合に、1つのカーネル103に含まれる波の数を増やすことが可能であり、カーネル103のサイズが固定の場合に比べて、コヒーレント領域の抽出精度を向上し、弾性率計測の精度を向上させることができる。
(実験例)
図8は比較例の説明図であり、固定したカーネルサイズを用いて抽出したコヒーレント領域をMR位相画像上に重ね合わせて表示した図である。図8Aは解像度が低い場合のMR位相画像に重ね合わせた図、図8Bは図8Aの場合において解像度が高い場合にMR位相画像に重ね合わせた図、図8Cは図8Bと同じ解像度の場合において付与される振動の周波数が低い(波長が長い)場合のMR位相画像に重ね合わせた図である。
図9は比較例の説明図であり、固定したカーネルサイズを用いて抽出したコヒーレント領域を弾性率マップ上に重ね合わせて表示した図である。図9Aは図8Aと同条件の場合で弾性率マップ上に重ね合わせた図、図9Bは図8Bと同条件の場合で弾性率マップ上に重ね合わせた図、図9Cは図8Cと同条件の場合で弾性率マップ上に重ね合わせた図である。
図10は実験例の説明図であり、カーネルサイズを可変にした実施例1で抽出したコヒーレント領域をMR位相画像上に重ね合わせて表示している。図10Aは解像度が図7Aと同条件の場合でMR位相画像上に重ね合わせた図、図10Bは図7Bと同条件の場合でMR位相画像上に重ね合わせた図、図10Cは図7Cと同条件の場合でMR位相画像上に重ね合わせた図である。
図11は実験例の説明図であり、カーネルサイズを可変にした実施例1で抽出したコヒーレント領域を弾性率マップ上に重ね合わせて表示した図である。図11Aは図10Aと同条件の場合で弾性率マップ上に重ね合わせた図、図11Bは図10Bと同条件の場合で弾性率マップ上に重ね合わせた図、図11Cは図10Cと同条件の場合で弾性率マップ上に重ね合わせた図である。
図8は比較例の説明図であり、固定したカーネルサイズを用いて抽出したコヒーレント領域をMR位相画像上に重ね合わせて表示した図である。図8Aは解像度が低い場合のMR位相画像に重ね合わせた図、図8Bは図8Aの場合において解像度が高い場合にMR位相画像に重ね合わせた図、図8Cは図8Bと同じ解像度の場合において付与される振動の周波数が低い(波長が長い)場合のMR位相画像に重ね合わせた図である。
図9は比較例の説明図であり、固定したカーネルサイズを用いて抽出したコヒーレント領域を弾性率マップ上に重ね合わせて表示した図である。図9Aは図8Aと同条件の場合で弾性率マップ上に重ね合わせた図、図9Bは図8Bと同条件の場合で弾性率マップ上に重ね合わせた図、図9Cは図8Cと同条件の場合で弾性率マップ上に重ね合わせた図である。
図10は実験例の説明図であり、カーネルサイズを可変にした実施例1で抽出したコヒーレント領域をMR位相画像上に重ね合わせて表示している。図10Aは解像度が図7Aと同条件の場合でMR位相画像上に重ね合わせた図、図10Bは図7Bと同条件の場合でMR位相画像上に重ね合わせた図、図10Cは図7Cと同条件の場合でMR位相画像上に重ね合わせた図である。
図11は実験例の説明図であり、カーネルサイズを可変にした実施例1で抽出したコヒーレント領域を弾性率マップ上に重ね合わせて表示した図である。図11Aは図10Aと同条件の場合で弾性率マップ上に重ね合わせた図、図11Bは図10Bと同条件の場合で弾性率マップ上に重ね合わせた図、図11Cは図10Cと同条件の場合で弾性率マップ上に重ね合わせた図である。
次に、本発明の効果を確認するために実験を行った。
なお、実施例の除去手段54pでDFを適用すると、弾性率の誤差が小さくなり、比較例と実験例との差が分かりにくくなる(弾性率の誤差がごまかされてしまう)。よって、図9、図11では、比較例と実験例との差をわかりやすくするために、除去手段54pでのDFを適用する前の弾性率マップを示している。
図8Aでは、付与された振動の周波数が100Hz、マトリクスサイズが128×128画素のMR位相画像に対し、13×13画素に固定したカーネルサイズを用いてコヒーレント領域を抽出した。図9Aでは、図8AのMR位相画像から作成した弾性率マップ上に抽出したコヒーレント領域を重ね合わせて表示している。図10Aでは、図8AのMR位相画像に対し、カーネル103のサイズが可変の実施例1の構成でコヒーレント領域の抽出を行った。図11Aでは、弾性率マップ上に可変のカーネルサイズで抽出したコヒーレント領域を重ね合わせて表示している。図9A(図8A)、図11A(図10A)からわかるように、図9Aでは、波長に対してカーネル103のサイズが大きすぎるために一部(特に画像の左上部分)のコヒーレント領域が抽出できていなかった部分についても、図11Aに示す実施例1では十分な精度で抽出できることがわかる。
なお、実施例の除去手段54pでDFを適用すると、弾性率の誤差が小さくなり、比較例と実験例との差が分かりにくくなる(弾性率の誤差がごまかされてしまう)。よって、図9、図11では、比較例と実験例との差をわかりやすくするために、除去手段54pでのDFを適用する前の弾性率マップを示している。
図8Aでは、付与された振動の周波数が100Hz、マトリクスサイズが128×128画素のMR位相画像に対し、13×13画素に固定したカーネルサイズを用いてコヒーレント領域を抽出した。図9Aでは、図8AのMR位相画像から作成した弾性率マップ上に抽出したコヒーレント領域を重ね合わせて表示している。図10Aでは、図8AのMR位相画像に対し、カーネル103のサイズが可変の実施例1の構成でコヒーレント領域の抽出を行った。図11Aでは、弾性率マップ上に可変のカーネルサイズで抽出したコヒーレント領域を重ね合わせて表示している。図9A(図8A)、図11A(図10A)からわかるように、図9Aでは、波長に対してカーネル103のサイズが大きすぎるために一部(特に画像の左上部分)のコヒーレント領域が抽出できていなかった部分についても、図11Aに示す実施例1では十分な精度で抽出できることがわかる。
図8Bでは、付与された振動の周波数が100Hz、マトリクスサイズが256×256画素で実験を行った。なお、図8Bおよび図9Bで使用されたカーネル103のサイズは、13×13画素で固定されている。図9B(図8B)、図11B(図10B)において、図8Bの条件では、波長に対してカーネル103のサイズが適切であったため、実験結果に差異はほとんどなかった。
図8Cでは、付与された振動の周波数が50Hz、マトリクスサイズが256×256画素で実験を行った。なお、図8Cおよび図9Cで使用されたカーネル103のサイズは、13×13画素で固定されている。図9C(図8C)、図11C(図10C)において、図8Cの条件では、波長に対してカーネル103のサイズが小さすぎたため、エラー領域(特に、図9Cの右上の部分)がコヒーレント領域として認識されてしまったが、図11Cではエラー領域がコヒーレント領域内に含まれなくなっている。
図8Cでは、付与された振動の周波数が50Hz、マトリクスサイズが256×256画素で実験を行った。なお、図8Cおよび図9Cで使用されたカーネル103のサイズは、13×13画素で固定されている。図9C(図8C)、図11C(図10C)において、図8Cの条件では、波長に対してカーネル103のサイズが小さすぎたため、エラー領域(特に、図9Cの右上の部分)がコヒーレント領域として認識されてしまったが、図11Cではエラー領域がコヒーレント領域内に含まれなくなっている。
したがって、図8~図11の実験結果から、図8B(図9B)のようにたまたまカーネル103のサイズが適切であれば、精度が十分なコヒーレント領域の抽出を行うことも可能であるが、カーネル103の適切なサイズは、利用者の経験や技量に依存するところがあり、安定して高精度の弾性率計測を行うことが困難である。これに対して、実施例1では、伝播する振動の波長に応じて自動的にカーネル103のサイズが適切なサイズに調整されており、利用者の経験に頼ることなく、安定して高精度な弾性率計測を行うことが可能である。
また、実施例1では、振動波が伝播中に伝播の主方向が変化していっても、ボクセル毎に主方向が判別されて、主方向以外の振動波が除去される。よって、特定の一方向以外の振動波が除去される従来の構成に比べて、振動の主方向の波を誤って除去したり減衰させたりすることが低減され、MRE画像の精度を向上させることができる。
また、振動波の主方向を作業者(測定者)が手動で選択する従来構成では、作業者の技量や経験によって主方向以外の波を除去する精度が大きく異なってくる。これに対して、実施例1では、自動的に各ボクセルの主方向が判定されており、作業者の技量等に依存せず、安定して振動の主方向以外の波を除去することができる。
また、振動波の主方向を作業者(測定者)が手動で選択する従来構成では、作業者の技量や経験によって主方向以外の波を除去する精度が大きく異なってくる。これに対して、実施例1では、自動的に各ボクセルの主方向が判定されており、作業者の技量等に依存せず、安定して振動の主方向以外の波を除去することができる。
(変更例)
以上、本発明の実施例を詳述したが、本発明は、前記実施例に限定されるものではなく、特許請求の範囲に記載された本発明の要旨の範囲内で、種々の変更を行うことが可能である。本発明の変更例(H01)~(H07)を下記に例示する。
(H01)前記実施例において、磁石部2がリング状、いわゆる、トンネル型の磁気共鳴撮影装置を例示したが、これに限定されない。例えば、磁石部2がコの字型、いわゆる、オープン型の磁気共鳴撮影装置にも適用可能である。
(H02)前記実施例において、例示した具体的な数値は、設計や仕様等に応じて、任意に変更可能である。
以上、本発明の実施例を詳述したが、本発明は、前記実施例に限定されるものではなく、特許請求の範囲に記載された本発明の要旨の範囲内で、種々の変更を行うことが可能である。本発明の変更例(H01)~(H07)を下記に例示する。
(H01)前記実施例において、磁石部2がリング状、いわゆる、トンネル型の磁気共鳴撮影装置を例示したが、これに限定されない。例えば、磁石部2がコの字型、いわゆる、オープン型の磁気共鳴撮影装置にも適用可能である。
(H02)前記実施例において、例示した具体的な数値は、設計や仕様等に応じて、任意に変更可能である。
(H03)前記実施例において、振動パッド16の具体的な形状等は、設計や仕様等に応じて変更可能である。
(H04)前記実施例において、波長の判別は、閾値を超える/超えない、の2段階で判別を行う場合を例示したがこれに限定されない。2つ以上の閾値を使用して、3段階以上の判別を行うことも可能である。これに伴って、カーネル103のサイズも3段階以上とすることも可能である。
(H05)前記実施例において、カーネル103の形状が四角形状の場合を例示したがこれに限定されない。円形状や三角形状、六角形状等、任意の形状とすることが可能である。
(H06)前記実施例において、磁場を利用するMRエラストグラフィ技術に適用する場合を例示したが、これに限定されない。例えば、超音波を利用して弾性率を導出する超音波エラストグラフィにも適用可能である。
(H04)前記実施例において、波長の判別は、閾値を超える/超えない、の2段階で判別を行う場合を例示したがこれに限定されない。2つ以上の閾値を使用して、3段階以上の判別を行うことも可能である。これに伴って、カーネル103のサイズも3段階以上とすることも可能である。
(H05)前記実施例において、カーネル103の形状が四角形状の場合を例示したがこれに限定されない。円形状や三角形状、六角形状等、任意の形状とすることが可能である。
(H06)前記実施例において、磁場を利用するMRエラストグラフィ技術に適用する場合を例示したが、これに限定されない。例えば、超音波を利用して弾性率を導出する超音波エラストグラフィにも適用可能である。
(H07)前記実施例において、各手段54d~54iで説明したように3つの点103a~103cの各位置での振動の変位を使用して主方向の算出を行う場合を例示したがこれに限定されず、主方向が導出可能な任意の方法を採用可能である。例えば、主方向の導出方法の他の例として、以下のProfile Gradient法を使用することが可能である。
まず、各ボクセルにおいて離散フーリエ変換を行って、位相画像を導出し、位相画像で判別方向(PA:Profile Angle、プロファイル角度φ)に沿った位相の勾配(PG: Profile Gradient、振動初期位相の空間勾配、PG(φ))を取得する。この処理を、各PAに対して行う。
そして、PG(φ)に対して、PGの最大値(振幅)をPG0とし、プロファイル角度に対するPG(φ)の初期位相をθとした場合、PG(φ)=PG0・sin(φ+θ)が成立する。
この式において、未知の定数はPG0とθであり、実測値PG(φ)と推定値PG′(φ)との残差が最小となるPG0、θを最小二乗法で算出できる。
求まったPG0、θにおける推定値PG′(φ)が最小となるプロファイル角度φが、振動の伝播の主方向となる。
まず、各ボクセルにおいて離散フーリエ変換を行って、位相画像を導出し、位相画像で判別方向(PA:Profile Angle、プロファイル角度φ)に沿った位相の勾配(PG: Profile Gradient、振動初期位相の空間勾配、PG(φ))を取得する。この処理を、各PAに対して行う。
そして、PG(φ)に対して、PGの最大値(振幅)をPG0とし、プロファイル角度に対するPG(φ)の初期位相をθとした場合、PG(φ)=PG0・sin(φ+θ)が成立する。
この式において、未知の定数はPG0とθであり、実測値PG(φ)と推定値PG′(φ)との残差が最小となるPG0、θを最小二乗法で算出できる。
求まったPG0、θにおける推定値PG′(φ)が最小となるプロファイル角度φが、振動の伝播の主方向となる。
1…撮影装置、
2…磁場発生装置、
4…被検者、
16…振動付与部材、
54…MRE画像作成手段、
54c…MR位相画像の取得手段、
54f…変位取得手段、
54i…主方向判別手段、
54m…領域サイズ調整手段、
54p…除去手段、
101…MR位相画像の全画像領域、
103…判別対象領域、
111…判別方向、
113…判別波形、
u(r、t)…変位。
2…磁場発生装置、
4…被検者、
16…振動付与部材、
54…MRE画像作成手段、
54c…MR位相画像の取得手段、
54f…変位取得手段、
54i…主方向判別手段、
54m…領域サイズ調整手段、
54p…除去手段、
101…MR位相画像の全画像領域、
103…判別対象領域、
111…判別方向、
113…判別波形、
u(r、t)…変位。
Claims (2)
- 静磁場と、位置に応じて磁場が変化する傾斜磁場と、プロトンの磁気共鳴条件に基づいて予め設定された交番磁場と、を被検者の被検査部に対して発生させる磁場発生装置と、
被検査部に振動を付与する振動付与部材と、
前記傾斜磁場を印加する時期に応じて取得された電磁波に基づいて、電磁波信号の位相に応じたMR位相画像を取得するMR位相画像の取得手段であって、予め定められた時間間隔をあけて複数の前記MR位相画像を取得する前記MR位相画像の取得手段と、
MR位相画像の全画像領域の中から区分けされた判別対象領域において、伝播する振動の波長が長くなると、前記判別対象領域のサイズを大きくする領域サイズ調整手段と、
前記領域サイズ調整手段で調整された判別対象領域に基づいて、伝播する振動の主方向が揃っているか否かを判別する判別手段と、
前記判別手段で振動の主方向が揃っていない判別対象領域の波が除去されたMR位相画像に基づいて、MRエラストグラフィ画像を作成するMRE画像作成手段と、
を備えたことを特徴とする撮影装置。 - 前記振動付与部材で付与される振動の伝播波長が大きくなるにつれて、前記判別対象領域の大きさを大きくする前記領域サイズ調整手段、
を備えたことを特徴とする請求項1に記載の撮影装置。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2021113419A JP2023009819A (ja) | 2021-07-08 | 2021-07-08 | 撮影装置 |
PCT/JP2022/024691 WO2023282051A1 (ja) | 2021-07-08 | 2022-06-21 | 撮影装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2021113419A JP2023009819A (ja) | 2021-07-08 | 2021-07-08 | 撮影装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2023009819A true JP2023009819A (ja) | 2023-01-20 |
Family
ID=84801515
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2021113419A Pending JP2023009819A (ja) | 2021-07-08 | 2021-07-08 | 撮影装置 |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP2023009819A (ja) |
WO (1) | WO2023282051A1 (ja) |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5285244B2 (ja) * | 2007-07-13 | 2013-09-11 | 株式会社日立メディコ | 磁気共鳴イメージング装置 |
JP6985737B2 (ja) * | 2018-01-16 | 2021-12-22 | 東京都公立大学法人 | 撮影装置 |
JP7289522B2 (ja) * | 2019-12-20 | 2023-06-12 | 東京都公立大学法人 | 撮影装置 |
CN112244813B (zh) * | 2020-10-22 | 2022-08-26 | 无锡鸣石峻致医疗科技有限公司 | 一种低场核磁共振弹性测量方法及系统 |
-
2021
- 2021-07-08 JP JP2021113419A patent/JP2023009819A/ja active Pending
-
2022
- 2022-06-21 WO PCT/JP2022/024691 patent/WO2023282051A1/ja active Application Filing
Also Published As
Publication number | Publication date |
---|---|
WO2023282051A1 (ja) | 2023-01-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5843876B2 (ja) | 磁気共鳴イメージング装置および磁化率強調画像生成方法 | |
JP5443531B2 (ja) | 骨構造をキャラクタライズするための拡散型磁気共鳴法の手順 | |
JP5865262B2 (ja) | 電気特性断層画像化方法及びシステム | |
US10942239B2 (en) | MRI method for determining signature indices of an observed tissue from signal patterns obtained by motion-probing pulsed gradient MRI | |
RU2616984C2 (ru) | Магнитно-резонансная (mr) томография электрических свойств | |
US9709641B2 (en) | Magnetic resonance imaging apparatus, image processing apparatus, and susceptibility map calculation method | |
US9069998B2 (en) | Determining electrical properties of tissue using magnetic resonance imaging and least squared estimate | |
US10203387B2 (en) | MR imaging with enhanced susceptibility contrast | |
US10698060B2 (en) | Magnetic resonance imaging apparatus and medical image processing apparatus | |
JP3976515B2 (ja) | 磁気共鳴イメージング装置及び画像処理方法 | |
CN102103197A (zh) | 相位图像数据组中背景相位的确定方法 | |
US20220180575A1 (en) | Method and system for generating magnetic resonance image, and computer readable storage medium | |
JP7115889B2 (ja) | 医用画像診断支援装置、および、磁気共鳴イメージング装置 | |
JP7289522B2 (ja) | 撮影装置 | |
JP6928077B2 (ja) | ディクソン型の水/脂肪分離を用いるmr撮像 | |
US20140103925A1 (en) | Determining Electrical Properties of Tissue Using Complex Magnetic Resonance Images | |
WO2023282051A1 (ja) | 撮影装置 | |
JP7230149B2 (ja) | 磁気共鳴イメージングで得た画像の処理方法、画像処理プロブラム、及び、計算機 | |
US11519987B2 (en) | Magnetic resonance fingerprinting thermometry | |
US8952693B2 (en) | Method for principal frequency magnetic resonance elastography inversion | |
JP7171047B2 (ja) | 撮影装置 | |
JPWO2017134830A1 (ja) | 医用画像診断支援装置、および、磁気共鳴イメージング装置 | |
JP7151996B2 (ja) | 撮影装置 | |
JP2015019813A (ja) | 磁気共鳴イメージング装置及び補正用b1マップを計算する方法 | |
US20170011526A1 (en) | Method and magnetic resonance system for segmenting a balloon-type volume |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20240405 |