JPH04241839A - Ｍｒｉ画像処理方法 - Google Patents
Ｍｒｉ画像処理方法Info
- Publication number
- JPH04241839A JPH04241839A JP3012813A JP1281391A JPH04241839A JP H04241839 A JPH04241839 A JP H04241839A JP 3012813 A JP3012813 A JP 3012813A JP 1281391 A JP1281391 A JP 1281391A JP H04241839 A JPH04241839 A JP H04241839A
- Authority
- JP
- Japan
- Prior art keywords
- image
- images
- time
- composite image
- mri
- 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
- 238000003672 processing method Methods 0.000 title claims description 7
- 239000002131 composite material Substances 0.000 claims description 29
- 238000002595 magnetic resonance imaging Methods 0.000 claims description 19
- 238000003384 imaging method Methods 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 13
- 230000005415 magnetization Effects 0.000 claims description 6
- 238000000034 method Methods 0.000 description 16
- 238000002075 inversion recovery Methods 0.000 description 11
- 238000006243 chemical reaction Methods 0.000 description 7
- 230000002159 abnormal effect Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【０００１】
【産業上の利用分野】本発明は、医療用の磁気共鳴イメ
ージング（ＭＲＩ）装置（以下、必要に応じて単にＭＲ
Ｉという）において、新たに画像を撮像（スキャン）す
ることなく、既に撮像された一定条件下での複数枚の画
像群から、計算機上で新たな合成画像を生成するための
画像処理方法に関する。
ージング（ＭＲＩ）装置（以下、必要に応じて単にＭＲ
Ｉという）において、新たに画像を撮像（スキャン）す
ることなく、既に撮像された一定条件下での複数枚の画
像群から、計算機上で新たな合成画像を生成するための
画像処理方法に関する。
【０００２】
【従来の技術】従来、ＭＲＩにおいて画像に含まれるパ
ラメータとしては、単位面積当たりの共鳴にあずかる水
素原子核の数に相当するＰｄ（水素原子核（プロトン）
密度）や、スピンの動き易さないし置かれた環境を反映
するＴ１（縦磁化緩和時間）及びＴ２（横磁化緩和時間
），Ｖ（流速），ＴＥ（エコー時間），ＴＲ（繰返し時
間）等、数多く存在し、ＭＲＩで撮像した画像は上記パ
ラメータの組合せ方により臨床的に持つ意味が変わって
くる。このため、実際のルーチン検査においても各パラ
メータを様々に変えながら多くの画像を撮像せざるを得
ない。更に、ＭＲＩは一般的に他のイメージング装置（
例えばＸ線ＣＴ等）に比べて撮像時間が長いため、上記
のように臨床検査において複数のパラメータを組合せて
多くの撮像を繰り返すことは、検査時間の増大を招くば
かりでなく、患者にとっての負担も増すことになる。 一方、ＭＲＩにより撮像する場合、装置の物理的な制約
、例えばマグネット、電源、高周波装置、コンピュータ
等の切替時間の制約等により、各パラメータの範囲には
、実際に撮像できる限界が存在する。従って、各パラメ
ータの持つ物理的な意味から有効と思われる撮像パラメ
ータの組合せ（例えばエコー時間ＴＥや繰返し時間ＴＲ
の組合せ）によっては、実際上、撮像できないものが出
てくる場合があり、このような場合は実際の撮像手法で
は画像化することができなかった。
ラメータとしては、単位面積当たりの共鳴にあずかる水
素原子核の数に相当するＰｄ（水素原子核（プロトン）
密度）や、スピンの動き易さないし置かれた環境を反映
するＴ１（縦磁化緩和時間）及びＴ２（横磁化緩和時間
），Ｖ（流速），ＴＥ（エコー時間），ＴＲ（繰返し時
間）等、数多く存在し、ＭＲＩで撮像した画像は上記パ
ラメータの組合せ方により臨床的に持つ意味が変わって
くる。このため、実際のルーチン検査においても各パラ
メータを様々に変えながら多くの画像を撮像せざるを得
ない。更に、ＭＲＩは一般的に他のイメージング装置（
例えばＸ線ＣＴ等）に比べて撮像時間が長いため、上記
のように臨床検査において複数のパラメータを組合せて
多くの撮像を繰り返すことは、検査時間の増大を招くば
かりでなく、患者にとっての負担も増すことになる。 一方、ＭＲＩにより撮像する場合、装置の物理的な制約
、例えばマグネット、電源、高周波装置、コンピュータ
等の切替時間の制約等により、各パラメータの範囲には
、実際に撮像できる限界が存在する。従って、各パラメ
ータの持つ物理的な意味から有効と思われる撮像パラメ
ータの組合せ（例えばエコー時間ＴＥや繰返し時間ＴＲ
の組合せ）によっては、実際上、撮像できないものが出
てくる場合があり、このような場合は実際の撮像手法で
は画像化することができなかった。
【０００３】このような場合、従来のＭＲＩでは「合成
画像処理（Ｓｙｎｔｈｅｓｉｃ Ｉｍａｇｉｎｇ）」と
呼ばれる手法により、実際に撮像することなく計算機上
で人工的に新しい画像を生成することが行なわれている
。すなわち、上記合成画像処理は、実際に撮像した複数
の画像群の画素濃度値からＭＲＩ信号強度式に従って各
パラメータの空間分布画像（Ｐｄ画像，Ｔ１画像，Ｔ２
画像）を計算する（これらは計算画像と呼ばれる）と共
に、これらの計算画像に対しエコー時間ＴＥ及び繰返し
時間ＴＲを任意に与え、ＭＲＩ信号強度式に従って再度
、新しい画像を生成する方法である。
画像処理（Ｓｙｎｔｈｅｓｉｃ Ｉｍａｇｉｎｇ）」と
呼ばれる手法により、実際に撮像することなく計算機上
で人工的に新しい画像を生成することが行なわれている
。すなわち、上記合成画像処理は、実際に撮像した複数
の画像群の画素濃度値からＭＲＩ信号強度式に従って各
パラメータの空間分布画像（Ｐｄ画像，Ｔ１画像，Ｔ２
画像）を計算する（これらは計算画像と呼ばれる）と共
に、これらの計算画像に対しエコー時間ＴＥ及び繰返し
時間ＴＲを任意に与え、ＭＲＩ信号強度式に従って再度
、新しい画像を生成する方法である。
【０００４】
【発明が解決しようとする課題】しかし、一般にＭＲＩ
画像は元々画像に含まれるノイズが多いため、単にＰｄ
画像，Ｔ１画像，Ｔ２画像からＭＲＩ信号強度式に従っ
て生成した合成画像は一般に実際に撮像した画像に比べ
てノイズが一層増大した画像になることが多い。また、
合成画像を生成するための演算時間もかなりかかるため
、実際のルーチン検査に用いるには問題が多かった。更
に、Ｐｄ画像，Ｔ１画像，Ｔ２画像に基づき合成する画
像の計算式としてＭＲＩ信号強度式をそのまま使用して
いるため、Ｔ１強調，Ｔ２強調の重み付けを任意に変え
ることができないという問題があった。
画像は元々画像に含まれるノイズが多いため、単にＰｄ
画像，Ｔ１画像，Ｔ２画像からＭＲＩ信号強度式に従っ
て生成した合成画像は一般に実際に撮像した画像に比べ
てノイズが一層増大した画像になることが多い。また、
合成画像を生成するための演算時間もかなりかかるため
、実際のルーチン検査に用いるには問題が多かった。更
に、Ｐｄ画像，Ｔ１画像，Ｔ２画像に基づき合成する画
像の計算式としてＭＲＩ信号強度式をそのまま使用して
いるため、Ｔ１強調，Ｔ２強調の重み付けを任意に変え
ることができないという問題があった。
【０００５】本発明は、上記問題点を解決するためにな
されたもので、その目的とするところは、上記従来の合
成画像処理手法を改善し、ノイズがより少なく臨床的に
有効な合成画像を短時間で得ることができるＭＲＩ画像
処理方法を提供することにある。
されたもので、その目的とするところは、上記従来の合
成画像処理手法を改善し、ノイズがより少なく臨床的に
有効な合成画像を短時間で得ることができるＭＲＩ画像
処理方法を提供することにある。
【０００６】
【課題を解決するための手段】上記目的を達成するため
に、本発明は、磁気共鳴イメージング装置により一定条
件下で撮像した複数枚の画像群から、画像を構成するパ
ラメータとしての縦磁化緩和時間Ｔ１、横磁化緩和時間
Ｔ２及びプロトン密度Ｐｄを求め、これらのＴ１，Ｔ２
，Ｐｄと別のパラメータであるエコー時間ＴＥ及び繰返
し時間ＴＲとに基づいて計算により合成画像を生成する
ＭＲＩ画像処理方法において、前記各パラメータＴ１，
Ｔ２，Ｐｄ，ＴＥ，ＴＲの合成画像に対する重み付けを
可変として任意に撮像条件を変更することにより合成画
像を生成するものである。
に、本発明は、磁気共鳴イメージング装置により一定条
件下で撮像した複数枚の画像群から、画像を構成するパ
ラメータとしての縦磁化緩和時間Ｔ１、横磁化緩和時間
Ｔ２及びプロトン密度Ｐｄを求め、これらのＴ１，Ｔ２
，Ｐｄと別のパラメータであるエコー時間ＴＥ及び繰返
し時間ＴＲとに基づいて計算により合成画像を生成する
ＭＲＩ画像処理方法において、前記各パラメータＴ１，
Ｔ２，Ｐｄ，ＴＥ，ＴＲの合成画像に対する重み付けを
可変として任意に撮像条件を変更することにより合成画
像を生成するものである。
【０００７】
【作用】本発明によれば、ＭＲＩ画像の各パラメータの
相関関係に自由度を持たせ、任意にＴ１強調、Ｔ２強調
等の重み付けが可能な合成画像を作成することができる
。 すなわち、画像パラメータＰｄ，Ｔ１，Ｔ２，ＴＥ，Ｔ
Ｒに基づく合成画像の演算式における重み付けパラメー
タを多数持ち、これらを所定の範囲で可変とすることに
より、実際の撮像で制約されるＰｄ，Ｔ１，Ｔ２，ＴＥ
，ＴＲの値に縛られない任意の合成画像が作成可能とな
る。 例えば実際の撮像では、Ｔ１とＴ２とは負の相関にあり
、ＴＥ，ＴＲは正の値でなければならないが、本発明に
よれば、これらの制約に関わりなく対象組織の特徴に応
じた合成画像を人工的に作成することができる。また、
Ｔ１画像を、変換表の折線近似により信号強度比に対応
するＴ１を探索することによって作成すれば、合成画像
を演算する過程で最も時間がかかるＴ１演算の短縮化が
可能になり、全体として合成画像の作成時間を大幅に短
縮することができる。更に、異常画素の補間や背景ノイ
ズの除去等の補正処理を行うことにより、画像演算過程
におけるノイズの発生を極力抑え、合成画像のＳ／Ｎ比
を向上させることができる。
相関関係に自由度を持たせ、任意にＴ１強調、Ｔ２強調
等の重み付けが可能な合成画像を作成することができる
。 すなわち、画像パラメータＰｄ，Ｔ１，Ｔ２，ＴＥ，Ｔ
Ｒに基づく合成画像の演算式における重み付けパラメー
タを多数持ち、これらを所定の範囲で可変とすることに
より、実際の撮像で制約されるＰｄ，Ｔ１，Ｔ２，ＴＥ
，ＴＲの値に縛られない任意の合成画像が作成可能とな
る。 例えば実際の撮像では、Ｔ１とＴ２とは負の相関にあり
、ＴＥ，ＴＲは正の値でなければならないが、本発明に
よれば、これらの制約に関わりなく対象組織の特徴に応
じた合成画像を人工的に作成することができる。また、
Ｔ１画像を、変換表の折線近似により信号強度比に対応
するＴ１を探索することによって作成すれば、合成画像
を演算する過程で最も時間がかかるＴ１演算の短縮化が
可能になり、全体として合成画像の作成時間を大幅に短
縮することができる。更に、異常画素の補間や背景ノイ
ズの除去等の補正処理を行うことにより、画像演算過程
におけるノイズの発生を極力抑え、合成画像のＳ／Ｎ比
を向上させることができる。
【０００８】
【実施例】以下、図に沿って本発明の実施例を説明する
。まず、図１はこの実施例における画像処理の流れを示
している。本実施例において、まず、合成画像を作成す
るために予め撮像しておく画像群とその条件は以下に示
すとおりである。■繰返し時間ＴＲが等しくエコー時間
ＴＥが異なる２つのスピンエコー（ＳＥ）画像１１，１
２の組。ここでＳＥ画像とは、スピンエコー法により原
子核の核スピンを９０°傾けるように共鳴させたときの
画像であり、スピン密度及びＴ２を反映している。■エ
コー時間ＴＥが等しいＳＥ画像１１とインバージョンリ
カバリ（ＩＲ）画像２１との組。ここでＩＲ画像とは、
反転回復法により原子核の核スピンを１８０°反転させ
るように共鳴させたときの画像であり、Ｔ１が強調され
ている。■スピンの回転角度が異なる２つのフィールド
エコー（ＦＥ）画像３１，３２の組。ここで、ＦＥ画像
とはスピン密度及びＴ１が反映されている画像である。
。まず、図１はこの実施例における画像処理の流れを示
している。本実施例において、まず、合成画像を作成す
るために予め撮像しておく画像群とその条件は以下に示
すとおりである。■繰返し時間ＴＲが等しくエコー時間
ＴＥが異なる２つのスピンエコー（ＳＥ）画像１１，１
２の組。ここでＳＥ画像とは、スピンエコー法により原
子核の核スピンを９０°傾けるように共鳴させたときの
画像であり、スピン密度及びＴ２を反映している。■エ
コー時間ＴＥが等しいＳＥ画像１１とインバージョンリ
カバリ（ＩＲ）画像２１との組。ここでＩＲ画像とは、
反転回復法により原子核の核スピンを１８０°反転させ
るように共鳴させたときの画像であり、Ｔ１が強調され
ている。■スピンの回転角度が異なる２つのフィールド
エコー（ＦＥ）画像３１，３２の組。ここで、ＦＥ画像
とはスピン密度及びＴ１が反映されている画像である。
【０００９】この実施例においては、上記■における２
つのＳＥ画像１１，１２の組について各画像の画像信号
強度ＩＳＥ，ＩＳＥ′から信号強度比Ｒ２を算出し、こ
れに基づいてＴ２画像を演算する。また、上記■におけ
るＳＥ画像１１とＩＲ画像２１の組について各画像の画
像信号強度ＩＳＥ，ＩＩＲから信号強度比Ｒ１Ａを算出
し、これに基づいてＴ１画像を演算するか、あるいは、
上記■におけるＦＥ画像３１，３２の組について各画像
の画像信号強度ＩＦＥ，ＩＦＥ′から信号強度比Ｒ１Ｂ
を算出し、これに基づいてＴ１画像を演算する。ここで
各画像の信号強度を、各々のパラメータに応じてＩＳＥ
（ＴＥ，ＴＲ），ＩＳＥ′（ＴＥ′，ＴＲ），ＩＩＲ（
ＴＥ，ＴＲ′，ＴＩ），ＩＦＥ（ＴＥ，ＴＲ，θ），Ｉ
ＦＥ′（ＴＥ，ＴＲ，θ′）と表すものとする。なお、
各パラメータにおいてＴＥ′はＳＥ画像１２のエコー時
間、ＴＲ′はＩＲ画像２１の繰返し時間、θ′はＦＥ画
像３２のスピン回転角度、ＴＩは反転時間である。
つのＳＥ画像１１，１２の組について各画像の画像信号
強度ＩＳＥ，ＩＳＥ′から信号強度比Ｒ２を算出し、こ
れに基づいてＴ２画像を演算する。また、上記■におけ
るＳＥ画像１１とＩＲ画像２１の組について各画像の画
像信号強度ＩＳＥ，ＩＩＲから信号強度比Ｒ１Ａを算出
し、これに基づいてＴ１画像を演算するか、あるいは、
上記■におけるＦＥ画像３１，３２の組について各画像
の画像信号強度ＩＦＥ，ＩＦＥ′から信号強度比Ｒ１Ｂ
を算出し、これに基づいてＴ１画像を演算する。ここで
各画像の信号強度を、各々のパラメータに応じてＩＳＥ
（ＴＥ，ＴＲ），ＩＳＥ′（ＴＥ′，ＴＲ），ＩＩＲ（
ＴＥ，ＴＲ′，ＴＩ），ＩＦＥ（ＴＥ，ＴＲ，θ），Ｉ
ＦＥ′（ＴＥ，ＴＲ，θ′）と表すものとする。なお、
各パラメータにおいてＴＥ′はＳＥ画像１２のエコー時
間、ＴＲ′はＩＲ画像２１の繰返し時間、θ′はＦＥ画
像３２のスピン回転角度、ＴＩは反転時間である。
【００１０】次に、Ｔ１画像演算のための前記信号強度
比Ｒ１Ａは、次の数式１により算出するものとする。但
し、Ｒ１Ａが負の値になった場合（異常画素の場合）は
、後述する画素の補間を行なう。
比Ｒ１Ａは、次の数式１により算出するものとする。但
し、Ｒ１Ａが負の値になった場合（異常画素の場合）は
、後述する画素の補間を行なう。
【００１１】
【数１】Ｒ１Ａ＝（ＩＳＥ−ＩＩＲ）／ＩＳＥ
【００１
２】また、同じくＴ１計算画像のための前記信号強度比
Ｒ１Ｂは、次の数式２により算出するものとする。
２】また、同じくＴ１計算画像のための前記信号強度比
Ｒ１Ｂは、次の数式２により算出するものとする。
【００１３】
【数２】
Ｒ１Ｂ＝（ｓｉｎθ／ｓｉｎθ′）・（ＩＦＥ′／ＩＦ
Ｅ）
Ｅ）
【００１４】次に、Ｔ２画像演算のための前記信号
強度比Ｒ２は、次の数式３により算出するものとする。 但し、Ｒ２が負の値になった場合（異常画素の場合）は
、後述する画素の補間を行なう。
強度比Ｒ２は、次の数式３により算出するものとする。 但し、Ｒ２が負の値になった場合（異常画素の場合）は
、後述する画素の補間を行なう。
【００１５】
【数３】Ｒ２＝ＩＳＥ／ＩＳＥ′
【００１６】次いで、画素の補間方法について説明する
。基本的には、図２に示すように対象画素（異常画素）
ｅに隣接する８近傍画素からの線形補間により、対象と
する異常画素の値を求める。ここで、対象画素ｅ、８近
傍画素ａ〜ｄ，ｆ〜ｉの画素値をそのままｅ，ａ〜ｄ，
ｆ〜ｉとすると、対象画素の画素値ｅは次の数式４で表
される。
。基本的には、図２に示すように対象画素（異常画素）
ｅに隣接する８近傍画素からの線形補間により、対象と
する異常画素の値を求める。ここで、対象画素ｅ、８近
傍画素ａ〜ｄ，ｆ〜ｉの画素値をそのままｅ，ａ〜ｄ，
ｆ〜ｉとすると、対象画素の画素値ｅは次の数式４で表
される。
【００１７】
【数４】ｅ＝（ａ＋ｂ＋ｃ＋ｄ＋ｆ＋ｇ＋ｈ＋ｉ）／８
【００１８】また、図３のように対象画素ｅの８近傍画
素中に異常画素（図中、●にて示す）がある場合や、図
４のように対象画素ｅが縁にあったり図５のように角に
あって８近傍画素の一部が欠ける場合には、それぞれ残
りの画素により以下の数式５、数式６、数式７を各々用
いて補間する。
素中に異常画素（図中、●にて示す）がある場合や、図
４のように対象画素ｅが縁にあったり図５のように角に
あって８近傍画素の一部が欠ける場合には、それぞれ残
りの画素により以下の数式５、数式６、数式７を各々用
いて補間する。
【００１９】
【数５】
ｅ＝（ａ＋ｂ＋ｆ＋ｇ＋ｈ＋ｉ）／６
【００２０】
【数６】ｅ＝（ｄ＋ｆ＋ｇ＋ｈ＋ｉ）／５
【００２１】
【数７】ｅ＝（ｆ＋ｈ＋ｉ）／３
【００２２】次に、Ｔ１画像の計算方法について説明す
る。まず、ＳＥ画像１１及びＩＲ画像２１に基づく信号
強度比Ｒ１Ａを用いる場合には、次の数式８に従って計
算する。
る。まず、ＳＥ画像１１及びＩＲ画像２１に基づく信号
強度比Ｒ１Ａを用いる場合には、次の数式８に従って計
算する。
【００２３】
【数８】
【００２４】また、２つのＦＥ画像３１，３２に基づく
信号強度比Ｒ１Ｂを用いる場合には、次の数式９に従っ
て計算する。
信号強度比Ｒ１Ｂを用いる場合には、次の数式９に従っ
て計算する。
【００２５】
【数９】
【００２６】但し、上記の数式８又は数式９からＴ１を
一義的に求めることはできないので、上記数式に基づい
てＲ１Ａ，Ｒ１ＢとＴ１との変換表を作成し、この変換
表からＴ１を探索する。この変換表は、例えば引数に対
応させたＴ１をｎ等分し、それぞれのＴ１と信号強度比
との対応関係を表にしたものであり、この変換表のＴ１
及び信号強度比からなる座標を折線近似したグラフに基
づき信号強度比に対応するＴ１を求めるものである。な
お、図６は変換表に基づくＲ１ＡまたはＲ１ＢとＴ１と
の関係を折線近似したグラフであり、小さなプロット（
・）は変換値、大きなプロット（●）は折線近似式上の
値、また、ｓは変換表のサイズを示している。このよう
な手法により、図１に示すように信号強度比Ｒ１Ａ，Ｒ
１Ｂに基づいてＴ１画像４１Ａ，４１Ｂを計算する。
一義的に求めることはできないので、上記数式に基づい
てＲ１Ａ，Ｒ１ＢとＴ１との変換表を作成し、この変換
表からＴ１を探索する。この変換表は、例えば引数に対
応させたＴ１をｎ等分し、それぞれのＴ１と信号強度比
との対応関係を表にしたものであり、この変換表のＴ１
及び信号強度比からなる座標を折線近似したグラフに基
づき信号強度比に対応するＴ１を求めるものである。な
お、図６は変換表に基づくＲ１ＡまたはＲ１ＢとＴ１と
の関係を折線近似したグラフであり、小さなプロット（
・）は変換値、大きなプロット（●）は折線近似式上の
値、また、ｓは変換表のサイズを示している。このよう
な手法により、図１に示すように信号強度比Ｒ１Ａ，Ｒ
１Ｂに基づいてＴ１画像４１Ａ，４１Ｂを計算する。
【００２７】次に、２つのＳＥ画像１１，１２に基づく
信号強度比Ｒ２に従ってＴ２画像４２を計算するには、
以下の数式１０を用いる。なお、この数式１０において
ｌｎは自然対数である。また、ＴＥ′＞ＴＥを条件とす
る。
信号強度比Ｒ２に従ってＴ２画像４２を計算するには、
以下の数式１０を用いる。なお、この数式１０において
ｌｎは自然対数である。また、ＴＥ′＞ＴＥを条件とす
る。
【００２８】
【数１０】Ｔ２＝（ＴＥ′−ＴＥ）／ｌｎ（Ｒ２）
【０
０２９】なお、これらのＴ１画像及びＴ２画像はＳＥ画
像、ＩＲ画像及びＦＥ画像からなる複数枚の画像群に基
づいて計算し、各々加重平均を求めることで図１に示す
ように最終的なＴ１画像４１ＡＸ，４１ＢＸ及びＴ２画
像４２Ｘを算出する。
０２９】なお、これらのＴ１画像及びＴ２画像はＳＥ画
像、ＩＲ画像及びＦＥ画像からなる複数枚の画像群に基
づいて計算し、各々加重平均を求めることで図１に示す
ように最終的なＴ１画像４１ＡＸ，４１ＢＸ及びＴ２画
像４２Ｘを算出する。
【００３０】次に、Ｐｄ画像５１を、ＳＥ画像１１と、
２つのＴ１画像４１ＡＸ，４１ＢＸのうちの何れか一方
と、Ｔ２画像４２Ｘとを用いて、以下の数式１１により
計算する。
２つのＴ１画像４１ＡＸ，４１ＢＸのうちの何れか一方
と、Ｔ２画像４２Ｘとを用いて、以下の数式１１により
計算する。
【００３１】
【数１１】
Ｐｄ＝ＩＳＥ／｛１−ｅｘｐ（−ＴＲ／Ｔ１）｝・
ｅｘｐ（ＴＥ／Ｔ２）
ｅｘｐ（ＴＥ／Ｔ２）
【００３２】そして、このＰｄ画
像５１と２つのＴ１画像４１ＡＸ，４１ＢＸのうちの何
れか一方とＴ２画像４２Ｘとを用いて、以下の数式１２
，数式１３，数式１４により中間合成画像６１としての
ＳＥ合成画像、ＩＲ合成画像、ＦＥ合成画像をそれぞれ
計算する。なお、数式１２において合成画像の重み付け
パラメータであるｐ，ｑ，ｊ，ｍ，ｎは、０＜ｐ，ｑ＜
１、並びに、０＜ｊ，ｍ，ｎ＜１の関係にある。更に、
数式１３，数式１４において合成画像の重み付けパラメ
ータであるｐ，ｑ，ｒ，ｊ，ｍ，ｎは、０＜ｐ，ｑ，ｒ
＜１、並びに、０＜ｊ，ｍ，ｎ＜１の関係にある。これ
らの重み付けパラメータは上述した範囲内で任意に変更
することができ、これによって任意の撮像条件を得るこ
とができる。
像５１と２つのＴ１画像４１ＡＸ，４１ＢＸのうちの何
れか一方とＴ２画像４２Ｘとを用いて、以下の数式１２
，数式１３，数式１４により中間合成画像６１としての
ＳＥ合成画像、ＩＲ合成画像、ＦＥ合成画像をそれぞれ
計算する。なお、数式１２において合成画像の重み付け
パラメータであるｐ，ｑ，ｊ，ｍ，ｎは、０＜ｐ，ｑ＜
１、並びに、０＜ｊ，ｍ，ｎ＜１の関係にある。更に、
数式１３，数式１４において合成画像の重み付けパラメ
ータであるｐ，ｑ，ｒ，ｊ，ｍ，ｎは、０＜ｐ，ｑ，ｒ
＜１、並びに、０＜ｊ，ｍ，ｎ＜１の関係にある。これ
らの重み付けパラメータは上述した範囲内で任意に変更
することができ、これによって任意の撮像条件を得るこ
とができる。
【００３３】
【数１２】
【００３４】
【数１３】
【００３５】
【数１４】
【００３６】こうして計算された中間合成画像６１に対
し、次にマスク処理を行って背景ノイズを除去する。す
なわち、人体の断層画像以外の背景部分（周囲の空気の
部分であるからＡｉｒ Ｓｃａｎともいう）は、上記合
成画像の計算過程により背景としての画素値から外れる
部分が多くなり、診断上好ましくないため、所定のしき
い値によって画像を２値化したマスクパターン８１によ
り背景ノイズを除去する。なお、この背景ノイズの除去
方法は本発明の要旨ではないため、詳述を省略するが、
本出願人による特願平２−３０７８７６号に詳しく説明
されている。以上のようにして、ＳＥ画像、ＩＲ画像及
びＦＥ画像等の複数枚の画像群から、最終的な合成画像
７１を得ることができる。
し、次にマスク処理を行って背景ノイズを除去する。す
なわち、人体の断層画像以外の背景部分（周囲の空気の
部分であるからＡｉｒ Ｓｃａｎともいう）は、上記合
成画像の計算過程により背景としての画素値から外れる
部分が多くなり、診断上好ましくないため、所定のしき
い値によって画像を２値化したマスクパターン８１によ
り背景ノイズを除去する。なお、この背景ノイズの除去
方法は本発明の要旨ではないため、詳述を省略するが、
本出願人による特願平２−３０７８７６号に詳しく説明
されている。以上のようにして、ＳＥ画像、ＩＲ画像及
びＦＥ画像等の複数枚の画像群から、最終的な合成画像
７１を得ることができる。
【００３７】
【発明の効果】以上のように本発明によれば、画像パラ
メータＴ１，Ｔ２，Ｐｄ，ＴＥ，ＴＲの合成画像に対す
る重み付けを可変として任意に撮像条件を変更するもの
であるため、いわゆるｓｈｏｒｔＳＥ画像のようなＴ１
強調画像からｌｏｎｇＳＥ画像のようなＴ２強調画像を
精度覧よく作成することができる。これにより、スキャ
ン時間が比較的長くかかるＴ２強調画像を人工的に合成
することができる。また、実際の撮像では物理的にスキ
ャンできない条件、例えばエコー時間ＴＥが非常に短い
場合や繰返し時間ＴＲが非常に長い場合等における画像
での組織の変化を見たい場合には、撮像条件を計算機上
で変えることで目的とする画像を人工的かつ簡単に合成
することができる。
メータＴ１，Ｔ２，Ｐｄ，ＴＥ，ＴＲの合成画像に対す
る重み付けを可変として任意に撮像条件を変更するもの
であるため、いわゆるｓｈｏｒｔＳＥ画像のようなＴ１
強調画像からｌｏｎｇＳＥ画像のようなＴ２強調画像を
精度覧よく作成することができる。これにより、スキャ
ン時間が比較的長くかかるＴ２強調画像を人工的に合成
することができる。また、実際の撮像では物理的にスキ
ャンできない条件、例えばエコー時間ＴＥが非常に短い
場合や繰返し時間ＴＲが非常に長い場合等における画像
での組織の変化を見たい場合には、撮像条件を計算機上
で変えることで目的とする画像を人工的かつ簡単に合成
することができる。
【００３８】また、一般にＭＲＩ画像は、例えばＳＥ画
像の場合、Ｔ１が増加することによる画素濃度値の増減
と、Ｔ２が増加することによる画素濃度値の増減とが反
比例しており、従って濃度値から見ると、Ｔ１とＴ２の
効果が相殺されていることになる。すなわち、∂Ｄ／∂
Ｔ１と∂Ｄ／∂Ｔ２（Ｄは画像濃度値）とが負の相関を
もっているところを、本発明によって正の相関に変更し
た合成画像を作成することにより、Ｔ１とＴ２との相関
効果が画像の濃度値に正確に現れるようにすることがで
きる。
像の場合、Ｔ１が増加することによる画素濃度値の増減
と、Ｔ２が増加することによる画素濃度値の増減とが反
比例しており、従って濃度値から見ると、Ｔ１とＴ２の
効果が相殺されていることになる。すなわち、∂Ｄ／∂
Ｔ１と∂Ｄ／∂Ｔ２（Ｄは画像濃度値）とが負の相関を
もっているところを、本発明によって正の相関に変更し
た合成画像を作成することにより、Ｔ１とＴ２との相関
効果が画像の濃度値に正確に現れるようにすることがで
きる。
【００３９】このように本発明によれば、ノイズが少な
く所望の撮像条件による臨床的に有効な合成画像を、極
めて短時間のうちに人工的に作成することができるもの
である。
く所望の撮像条件による臨床的に有効な合成画像を、極
めて短時間のうちに人工的に作成することができるもの
である。
【図１】本発明の一実施例における全体の処理の流れを
示す図である。
示す図である。
【図２】画素の補間方法の説明図である。
【図３】画素の補間方法の説明図である。
【図４】画素の補間方法の説明図である。
【図５】画素の補間方法の説明図である。
【図６】Ｔ１計算方法の説明図である。
１１，１２ ＳＥ画像
２１ ＩＲ画像
３１，３２ ＦＥ画像
４１Ａ，４１Ｂ，４１ＡＸ，４１ＢＸ Ｔ１画像４２
，４２Ｘ Ｔ２画像 ５１ Ｐｄ画像 ６１ 中間合成画像 ７１ 合成画像 ８１ マスクパターン
，４２Ｘ Ｔ２画像 ５１ Ｐｄ画像 ６１ 中間合成画像 ７１ 合成画像 ８１ マスクパターン
Claims (1)
- 【請求項１】 磁気共鳴イメージング装置により一定
条件下で撮像した複数枚の画像群から、画像を構成する
パラメータとしての縦磁化緩和時間Ｔ１、横磁化緩和時
間Ｔ２及びプロトン密度Ｐｄを求め、これらのＴ１，Ｔ
２，Ｐｄと別のパラメータとしてのエコー時間ＴＥ及び
繰返し時間ＴＲとに基づいて計算により合成画像を生成
するＭＲＩ画像処理方法において、前記各パラメータＴ
１，Ｔ２，Ｐｄ，ＴＥ，ＴＲの合成画像に対する重み付
けを可変として任意に撮像条件を変更することにより合
成画像を生成することを特徴とするＭＲＩ画像処理方法
。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP3012813A JPH04241839A (ja) | 1991-01-08 | 1991-01-08 | Ｍｒｉ画像処理方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP3012813A JPH04241839A (ja) | 1991-01-08 | 1991-01-08 | Ｍｒｉ画像処理方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH04241839A true JPH04241839A (ja) | 1992-08-28 |
Family
ID=11815828
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP3012813A Pending JPH04241839A (ja) | 1991-01-08 | 1991-01-08 | Ｍｒｉ画像処理方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH04241839A (ja) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002086821A1 (fr) * | 2001-04-19 | 2002-10-31 | Kabushiki Kaisha Toshiba | Procede et dispositif de traitement d'image |
| JP2002323567A (ja) * | 2001-04-26 | 2002-11-08 | Shimadzu Corp | Ｅｃｔ装置 |
| EP1408340A2 (en) * | 2002-10-11 | 2004-04-14 | Communications Research Laboratory, Independent Administrative Institution | Method for processing magnetic resonance imaging image information |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6264348A (ja) * | 1985-09-05 | 1987-03-23 | 横河電機株式会社 | 核磁気共鳴撮像装置 |
| JPS6264347A (ja) * | 1985-09-05 | 1987-03-23 | 横河電機株式会社 | 核磁気共鳴撮像装置 |
-
1991
- 1991-01-08 JP JP3012813A patent/JPH04241839A/ja active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6264348A (ja) * | 1985-09-05 | 1987-03-23 | 横河電機株式会社 | 核磁気共鳴撮像装置 |
| JPS6264347A (ja) * | 1985-09-05 | 1987-03-23 | 横河電機株式会社 | 核磁気共鳴撮像装置 |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002086821A1 (fr) * | 2001-04-19 | 2002-10-31 | Kabushiki Kaisha Toshiba | Procede et dispositif de traitement d'image |
| US7492947B2 (en) | 2001-04-19 | 2009-02-17 | Kabushiki Kaisha Toshiba | Image processing method and image processing apparatus |
| US7974473B2 (en) | 2001-04-19 | 2011-07-05 | Kabushiki Kaisha Toshiba | Image processing method and image processing apparatus |
| JP2002323567A (ja) * | 2001-04-26 | 2002-11-08 | Shimadzu Corp | Ｅｃｔ装置 |
| EP1408340A2 (en) * | 2002-10-11 | 2004-04-14 | Communications Research Laboratory, Independent Administrative Institution | Method for processing magnetic resonance imaging image information |
| EP1408340A3 (en) * | 2002-10-11 | 2005-10-05 | National Institute of Information and Communications Technology Incorporated Administrative Agency | Method for processing magnetic resonance imaging image information |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kociołek et al. | Does image normalization and intensity resolution impact texture classification? | |
| US11341616B2 (en) | Methods and system for selective removal of streak artifacts and noise from images using deep neural networks | |
| US7949170B2 (en) | Image processing method, image processing device, computer aided detection, and method for filtering along the time axis | |
| JP3976684B2 (ja) | 画像における動きの影響を低減する方法および装置 | |
| RU2216751C2 (ru) | Способ и устройство для формирования изображений магнитного резонанса | |
| EP1506427B1 (en) | Real-time tractography | |
| WO2019026407A1 (ja) | 医用撮像装置及び医用画像処理方法 | |
| US8213715B1 (en) | System and method for digital image intensity correction | |
| US4989142A (en) | Three-dimensional images obtained from tomographic slices with gantry tilt | |
| Moratal et al. | k-Space tutorial: an MRI educational tool for a better understanding of k-space | |
| JP2000513978A (ja) | アーチファクト減少像生成方法及びその装置 | |
| JP2005525208A (ja) | Ｄｔｉｓｎｒを改善するための画像の遡及的な選択及び様々なタイプの画像アラインメント | |
| JP4070268B2 (ja) | 画像の分解能を向上させる方法、投影画像の一部を拡大させる方法及び画像の分解能を向上させる装置 | |
| US20190033409A1 (en) | Reconstructing magnetic resonance images with different contrasts | |
| KR20010113535A (ko) | 이미지 처리 장치 및 그 방법, 촬상 장치 및 기록 매체 | |
| JP5154751B2 (ja) | 医用画像処理装置 | |
| US6192265B1 (en) | Diagnostic image processing method | |
| US6862334B2 (en) | Method of reducing noise in volume imaging | |
| US20060066911A1 (en) | Edge detection and correcting system and method | |
| JPH04241839A (ja) | Ｍｒｉ画像処理方法 | |
| JP2985934B2 (ja) | Ｍｒイメージング装置 | |
| JP2008054738A (ja) | 磁気共鳴イメージング装置 | |
| Karim et al. | Medical image zooming by using rational bicubic ball function | |
| Aliaksandrauna | Adjusting videoendoscopic 3D reconstruction results using tomographic data | |
| JP7228386B2 (ja) | 医用画像処理装置、磁気共鳴イメージング装置、および医用画像処理方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 19980721 |