JP5385394B2 - ハイブリッド型の核/mr画像形成における透過データを使用したmrセグメンテーション - Google Patents

ハイブリッド型の核/mr画像形成における透過データを使用したmrセグメンテーション Download PDF

Info

Publication number
JP5385394B2
JP5385394B2 JP2011526607A JP2011526607A JP5385394B2 JP 5385394 B2 JP5385394 B2 JP 5385394B2 JP 2011526607 A JP2011526607 A JP 2011526607A JP 2011526607 A JP2011526607 A JP 2011526607A JP 5385394 B2 JP5385394 B2 JP 5385394B2
Authority
JP
Japan
Prior art keywords
map
data
transmission data
pet
radiation
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
Application number
JP2011526607A
Other languages
English (en)
Other versions
JP2012503177A (ja
Inventor
ガニョン,ダニエル
フゥ,ジチュイアン
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
Koninklijke Philips Electronics NV
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 Koninklijke Philips NV, Koninklijke Philips Electronics NV filed Critical Koninklijke Philips NV
Publication of JP2012503177A publication Critical patent/JP2012503177A/ja
Application granted granted Critical
Publication of JP5385394B2 publication Critical patent/JP5385394B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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/4808Multimodal MR, e.g. MR combined with positron emission tomography [PET], MR combined with ultrasound or MR combined with computed tomography [CT]
    • G01R33/481MR combined with positron emission tomography [PET] or single photon emission computed tomography [SPECT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • 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/5608Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10088Magnetic resonance imaging [MRI]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10104Positron emission tomography [PET]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing

Description

本発明は、解剖学的な画像形成システムにおける特定の応用を見い出すものであり、特に、結合されたPET-MR画像形成を含むものであるが、他の核画像形成システム等における応用をも見い出すものである。
しかし、記載される技術は、他の画像形成システム、他の画像形成シナリオ、他の画像分析技術等における応用を見い出す場合があることを理解されたい。
陽電子放出型断層撮影(PET: Position Emission Tomography)による画像形成では、PETデータにおける減衰は、典型的に減衰マップを使用して補正される。結合されたPETによるコンピュータ断層撮影(CT)では、CTデータは、PETデータについて正確な減衰マップを容易に生成することができる。CTデータ及び画像は、画像形成された組織の放射線の減衰特性に基づいている。しかし、組み合わせたPET-磁気共鳴(MR: magnetic Resonance)画像形成では、MR画像は、典型的には水ダイポール(H1)からなる共振特性を示し、この特性は、患者の放射線の減衰特性の関数ではない。
PETスキャナは、CTスキャナと通常は結合されるか、最近ではMRスキャナと結合される。MR画像は、解剖学的臓器又は類似の組織のタイプからなる分割された領域を示すマップを提供する。たとえば、MR画像は、通常、暗い領域として骨を示し、これは、暗い領域として示される空気から骨を区別することにおいて問題を生じさせる可能性がある。空気及び骨は、非常に異なる減衰特性を有する。骨の構造は、典型的に、表面において密度の高い皮質骨又はシェル、及び内部において弾力のある骨梁を含んでいる。解剖学の情報は、典型的に、骨と空気の領域とを区別することが必要とされる。皮質のシェルがエアポケットに隣接する場合、MR画像は、両者を1つの暗い領域として示す場合がある。従って、皮質骨には、減衰マップにおいて誤った値(大気に一致する値)が割り当てられる場合があり、逆に、空気には、減衰マップにおいて誤った値(皮質骨に一致する値)が割り当てられる場合がある。皮質骨の領域を分割する技術は、存在するものの、それらは常に信頼できるものではない。
本発明は、上述された課題及び他の課題を克服する、核画像における減衰を補正するMRに基づいた減衰補正マップを洗練“refine”させる新たな洗練されたシステム及び方法を提供するものである。
1態様によれば、解剖学的な画像補正システムは、MR取得スキャンの間に被検体のMR画像データを取得する磁気共鳴(MR)イメージャ、及び、核取得スキャンの間に被検体の核画像データを取得し、同時に、核スキャナの検査領域に位置される放射線源からの透過データを測定する核スキャナを含む。本システムは、MR画像データからの減衰補正(AC)マップを生成し、測定された透過データを使用して、生成されたACマップを繰り返し洗練させ、洗練されたACマップを生成し、洗練されたACマップを使用して、減衰について核画像データを補正する。
別の態様によれば、磁気共鳴(MR)減衰補正(AC)マップを洗練させる方法は、被検体のMR画像から減衰補正(AC)マップを生成するステップと、被検体を通して、被検体から離れて配置されている放射線源から放射線を放出するステップと、被検体を通して透過された放射線から透過データを測定するステップとを含む。本方法は、ACマップから推定された透過データを生成するステップと、測定された透過データと推定された透過データとの比較に基づいてACマップを調整することで、洗練されたACマップを生成するステップを更に含む。
別の態様によれば、PETスキャナは、放射線検出器のリング、放射線検出器のリング内に固定して設けられる放射線源、請求項1に係る方法を実行して、洗練された減衰補正(AC)マップを生成するプロセッサを含む。さらに、PETスキャナは、洗練されたACマップにより放射線検出器からのPETデータを補正し、補正されたPETデータを減衰が補正されたPET画像に再構成する再構成アルゴリズムを更に含む。核イメージにおける減衰を補正する減衰補正マップを洗練させる装置は、被検体のMR画像から減衰補正(AC)マップを生成する手段、被検体から離れて配置され、被検体を通して放射線を透過させる手段、被検体を通して透過された放射線から透過データを測定する手段を含む。さらに、本装置は、ACマップから推定される透過データを生成する手段、測定された透過データと推定された透過データとの比較に基づいてACマップを調整することで、洗練されたACマップを生成する手段を更に含む。
本発明の1つの利点は、減衰補正マップの精度が改善されることである。
本発明の別の利点は、MRに基づく減衰補正マップにおける骨のボクセルと空気のボクセルとの間の不明確さを解決することにある。
本発明の更なる利点は、以下の詳細な説明を読んで理解することで当業者により理解されるであろう。
本発明は、様々なコンポーネント及びコンポーネントのアレンジメントの形態、様々なステップ及びステップのアレンジメントの形態をとる場合がある。図面は、様々な態様を説明するためのものであって、本発明を限定するものとして解釈されるべきものではない。
本発明の様々な態様に係る、MR減衰マップを洗練させるポイント又はラインソースの使用を容易にするシステムを例示する図である。 本発明の様々な態様に係る、核画像形成装置におけるポイントソースの位置を例示する図である。 測定されたポイントソースの透過データと、それぞれのデータが含まれる複数のヒストグラムのビンとの間のグラフ上の関係を例示する図である。 レイトレーシングアルゴリズムを使用してプロセッサにより生成されたシミュレートされた又は推定されたポイントソースの透過データと、それぞれのシミューレートされたデータが含まれる複数のヒストグラムのビンとの間のグラフ上の関係を例示する図である。 測定されたポイントソースの透過データの値と、それぞれの差分データが含まれる複数のヒストグラムのビンとの間の違いに関するグラフ上の関係を例示する図である。 本発明の様々な態様に係る、核イメージにおける減衰を補正する洗練されたACマップを生成する方法を例示する図である。 3次元画像表現を生成するために個々のプロセッサ又は共有される再構成プロセッサにより再構成される画像データを生成する、MR画像形成装置と核スキャナのような複数の画像形成装置を含む例示的なホスピタルシステムを例示する図である。
図1は、本発明の様々な態様に係る、MR減衰マップを洗練させるポイントソースの使用を容易にするシステム10を例示する図である。MR画像形成又はスキャニング装置12は、MR減衰補正(AC)マップ16を生成するために身体の輪郭、肺、軟組織、骨構造等が分割される被検体の組織又は他の関心のあるボリュームのMRに基づく画像14を生成する。たとえばポイントソースである1以上の放射線源18は、被検体(たとえば関心のあるボリューム又は患者)の境界の外側であって、核画像形成装置20の検査領域内に配置される。ポイントソースの減衰又は透過データ22は、PETスキャナ又は単電子放出型断層撮影(SPECT: Single Photon Emission Computed Tomography)スキャナのような核画像形成装置20により測定される。なお、減衰及び透過は、一方が、吸収されない(透過される)放射線の全体量の一部を暗示する意味で、互いに逆のことを意味する。数学的に、何れかを使用することができる。本実施の形態では「透過」データが使用されるが、透過データは、透過情報と減衰情報の両者を含むことを理解されたい。1つの固定されたポイントソースから、被検体又は患者の2次元(2D)の投影核イメージ23が生成される。その間、推定される透過データ24は、MR減衰マップ16を使用して被検体のボリュームを通して投影するレイトレーシング技術又はアルゴリズム26を使用してシミュレートされる。測定されたポイントソースの透過データ22及びシミュレートされた透過データ24は、コンパレータ28により比較され、減衰マップ上で実行される反復的な微調整アルゴリズム30により差が最小にされる。従って、準最適なMR減衰マップ16は、測定されたポイントソースの透過データを使用して最適化され、洗練された(たとえば最適な)MR減衰マップ32が生成される。
MR画像が異なる組織のタイプに分割されるとき、それぞれの組織のタイプは、対応する放射線の減衰特性について推定される。それぞれのセグメントのMR画像及び推定された透過値から、被検体を通して任意の放射線に沿った推定される透過を計算することができる。実際に収集されるポイントソースの透過のそれぞれの放射線に沿って、推定される透過が計算される。対応する放射線に沿って推定された透過と実際に測定された透過とを比較することで、推定された透過値における誤差が決定される。多数の放射線が、異なる分割された領域を通過するのと同様に、共通の分割された領域であって、それぞれの領域において異なる経路長をもつ領域をも通過する。空気であると推定された領域及び骨であると推定された領域を容易に確かめることができる。さらに、放射線に沿って、推定された透過と実際に測定された透過との間の適合が最適化されるまで、それぞれ分割された領域の減衰/透過値は、減衰ファクタを繰り返し調整することで改善される。
同様に、予め選択された閾値にないデータを最適化することが不可能なことは、患者の動きを示す場合がある。身体の輪郭及び全体として内部領域の位置を繰り返し調節すること、又は身体の輪郭内にある内部領域の位置を繰り返し調節することで、患者の動きに適合すること、及び/又は患者の動きを検出することを最適化することができる。
MR画像形成装置12及び核イメージング又はスキャニング装置20は、様々な実施の形態によれば、個別の画像形成装置であるか、組み合わせ又はデュアルモダリティ画像形成装置である場合がある。たとえば、デュアルモダリティイメージャの場合、被検体は、MR画像形成のモダリティを使用してスキャニングされ、次いで、ポイントソースが画像形成装置に位置された後、核イメージング(たとえばPET又はSPECT)モダリティを使用して再びスキャンされる。このようにして、被検体の動きは、スキャンの間で最小にすることができる。さらに、この実施の形態及び他の実施の形態では、ポイントソースの位置は、ポイントソースからの透過データの検出を容易にするために、核イメージングスキャンの間で一貫性があるように、一定に維持されるか、及び/又は予め決定される場合がある。
システム10は、MR画像データ14、準最適なACマップ16、測定されたポイントソースデータ22、シミュレートされた又は推定されたポイントソースの透過データ24、及び洗練されたACマップ32を記憶するメモリ34を含む。メモリは、シミュレートされた透過データ24を生成するために使用されるレイトレーシングアルゴリズム26、コンパレータ28(たとえばシミュレートされたデータを測定されたポイントソースデータに比較するコンピュータ実行可能な命令)、洗練されたACマップ32を生成するとき、ACマップ16を洗練させるために使用される反復的な微調整アルゴリズム(iterative fine-tuning algorithm)30を更に記憶する。
さらに、システム10は、メモリ34に記憶されたデータを分析し、メモリに記憶されるアルゴリズムを実行して、メモリ34における記憶のために新たなデータを生成するプロセッサ36を含む。たとえば、プロセッサは、1以上のMR再構成アルゴリズム38を実行して、被検体のMRスキャンの間に取得されたrawデータであるMRデータからMR画像14を再構成する。同様に、プロセッサは、1以上の核イメージ再構成アルゴリズム39を実行して、被検体の核スキャンの間に取得されたrawデータである核スキャンデータ41から投影核イメージ23を再構成する。さらに、プロセッサは、推定された透過データ24を生成し、コンパレータ28を実行して、推定されたデータ24と測定されたポイントソースの透過データ22との間の差を決定し、反復的な微調整アルゴリズム30を実行して、オリジナルのACマップ16から洗練されたACマップ32を生成する。本明細書で使用される「アルゴリズム」は、メモリ34に持続して記憶され、プロセッサ36により実行される1以上のコンピュータにより実行可能な命令を意味するように解釈される場合がある。
従って、システム10は、被検体と共に核画像形成の検査領域においてポイントソース(又はラインソース)を配置するように利用される。1例では、放射線のポイントソース又はラインソースは、イメージングアイソトープとポイント又はラインソースとを識別するのを容易にするため、被検体のPET画像を生成するために使用されるPETアイソトープとは異なるエネルギーをもつ放射性同位体である。ポイント又はラインソースからの放射線は、PET検出器により検出され、(たとえばメモリ34に記憶され、プロセッサ36により実行される)エネルギー識別アルゴリズム42を使用してPETデータから分離され、特に、被検体を通した投影である透過の放射線データを生成するために使用される。2以上の放射線源を使用するか、患者とは相対的に放射線源を回転することで、3次元の放射線の減衰マップが容易に生成される。ポイントソース又はラインソースは、ガントリー又は核イメージング装置20又は患者支持体の構造に一時的に挿入されるか、又は永続的に搭載することができる。たとえば、レシービング構造は、患者支持体又は穴を画定する構造における予め選択された既知の位置で定義される。ソースの既知の位置は、投影データの分析を容易にする。代替的に、ソースの位置は、投影データを分析することで決定することができる。
別の実施の形態によれば、減衰マップは、MRスキャナ12のような解剖のイメージングモダリティから導出される。良好な性能のため、多数のポイントソース18が実際に利用される場合がある。しかし、以下の例では、1つのポイントソースとのPET/MRの組み合わせが記載される。この例によれば、MRスキャナ12を使用して、被検体又は患者のMR画像が生成される。身体の輪郭及び他の内部の臓器の境界が確定される。実質的又は完全に軟組織である知られている領域が識別され、ラベル付けされる。空気及び骨の組織を含む、肺組織のような特徴付けすることができない領域は、ラベル付けされないままにされる。
同じ被検体又は患者の(たとえばPET又はSPECTといった)核イメージが生成される(たとえば核イメージデータが取得される)。同時に、異なるエネルギーのウィンドウが開き、患者又は被検体の境界の外であって、核イメージャのガントリにおける検査領域内に適切に意図される放射線のポイントソース18からのカウントを受ける。たとえば、SPECT画像形成のため、標準的なエネルギーウィンドウが使用される。PET画像形成のため、エネルギーウィンドウ及び「シングル“single”」モードが使用される。ポイントソースは、患者又は被検体を画像形成するために使用される放射線同位体から発生される放射データよりも高いエネルギーである(たとえば、511kevでPETデータを画像形成するとき、662kevで137Csである)。身体の輪郭又は外形は、散乱からポイントソースの放射線を識別するのを更に容易にする。さらに、PET放射線よりも高いエネルギーをもつポイントソースを使用することで、PET放射線の散乱からポイントソースの放射線を識別することが容易となる。これは、ポイントソースの放射線は、低エネルギーのPET放射線の散乱とオーバラップしないためである。ソースからのデータが取得され、メモリ34に記憶されるヒストグラム又は投影に蓄積される。減衰係数は、分割されたMR画像の既知の領域に(たとえばACマップ16を生成するときに、プロセッサにより)割り当てられる。外部のポイントソース(図2参照)と同じ幾何学的形状をもつMRに基づいた減衰マップ16を通して推定された投影データ24を生成するため、レイトレーシングアルゴリズムを使用して、数学モデルが構築される。減衰マップ16を通して推定された又はシミュレートされた数学投影データ24が、ポイントソースの測定された透過データ22に厳密に一致するまで、未知の領域には、空気又は骨の値に代替的に割り当てられる。適切に分割され、空気及び骨の組織の様々な不明確な領域について値が適切に割り当てられた洗練されたACマップ32は、核イメージ23の減衰の補正を実行するために使用することができる。ひとたび空気及び骨の領域が決定されると、様々な軟組織の領域の減衰ファクタが任意に最適化される。
図2は、本発明の1以上の態様に係る、核画像形成装置20に位置されるポイントソース18を例示する。(たとえば被検体、患者である)関心のあるボリューム60は、核画像形成装置20のガントリ64により定義される画像形成領域におけるテーブル又は患者支持体62に配置される。関心のあるボリュームは、関心のあるボリュームのMR画像を使用して容易に曖昧さをなくすことができない、(たとえば肺組織、大気又はガスポケット等といった)骨の組織又は空気の場合がある不明瞭な領域66を有する。したがって、ポイントソース18は、関心のある領域に隣接するテーブル62の近くに位置され、ガントリ64に搭載される核検出器70により検出される放射線透過ライン68を放出する。1実施の形態では、ガントリ64は、回転可能なガントリであり、核検出器70は、透過ライン68を検出するのと同様に、関心のあるボリュームに投与される放射性同位体からの核スキャンデータを取得して核イメージを生成するため、検査領域に関して回転される。
別の実施の形態では、異なる透過エネルギーを有する2以上のポイントソース18,18’は、完全な視野を形成するために検査領域の周りに配置される。このように、MR画像データにおける骨と空気のボクセルとの間の不明確さを解決するのを容易にするため、2以上のラインの統合が生成される。
別の実施の形態では、核検出器70は、PET検出器であり、取得されたPETデータにおける減衰は、PET画像を再構成するとき、上述された洗練されたACマップを使用して補正される。
別の実施の形態では、核検出器70は、SPECT(Single Position Emission Computed Tomography)検出器であり、取得されたSPECTデータにおける減衰は、SPECT画像を再構成するとき、上述された洗練されたACマップを使用して補正される。
図3Aは、測定されたポイントソースの透過データ22と、それぞれのデータが含まれる複数のヒストグラムビンとのグラフ上の関係71を例示する。PET減衰の補正のために容易に使用されるようになるまでACマップを繰り返し洗練させるために、測定されたポイントソースの透過データ22は、推定されたポイントソース透過データ24(図3B)と共に使用される。
図3Bは、レイトレーシングアルゴリズム26(図1)を使用してプロセッサ36により生成されたシミュレートされた又は推定されたポイントソースの透過データ24と、それぞれのシミュレートされたデータが含まれる複数のヒストグラムビンとの間のグラフ上の関係を例示する。推定されたポイントソースの透過データ24は、測定されたポイントソースデータ22と比較され、PET画像を再構成するとき、ACマップがPETスキャンデータにおける減衰の補正における使用の準備があるかを判定する。
図3Cは、測定されたポイントソースの透過データの値とシミュレートされた透過データの値との差と、それぞれの差のデータが含まれる複数のヒストグラムのビンとのグラフ上の関係74を例示する。差の値は、プロセッサ36が図1のコンパレータ28を実行するときに決定される。それぞれのACマップでの繰り返しの調節により、差の値が低減される。ひとたび、推定された透過データとシミュレートされた透過データとの間の差の値が十分に小さくなると(たとえば予め決定された閾値以下)、洗練されたACマップは、PET減衰の補正において使用されるために記憶される。
図4は、本発明の様々な態様に係る、核イメージにおける減衰を補正する洗練されたACマップを生成する方法を例示する。ステップ78で、被検体のMR画像が生成される。ステップ80で、ACマップは、被検体の分割されたMR画像から生成される。ステップ82で、被検体は、核スキャナ(たとえばPET又はSPECT)の検査領域に配置され、ポイントソースは、被検体の外側に位置され、且つ検査領域内に位置される。ポイントソースは、(たとえば約662kevである)137Csのような被検体に投与される放射線同位体(たとえば511kev)よりも高いエネルギーをもつガンマ線を放出する。代替的に、133Ba(たとえば約360kev)のような、(たとえば511kevよりも低い)低エネルギーがポイントソースについて使用される。
ステップ84で、核スキャナを使用して、ポイントソースからの透過データが測定され、被検体の核イメージデータが生成される。たとえば、シングルモードで動作しているPETは、図3Aに示されるように、核同位体の放出データと同時にポイントソースの透過データを測定する。ステップ86は、図3Bに示されるように、ACマップに関するレイトレーシング技術を使用して透過データが導出又は推定される。ステップ88で、図3Cに示されるように、ステップ84からの透過データとステップ86からの投影データとの間の差が決定される。ステップ90で、二乗平均平方根誤差(rms)等のような、ステップ88で決定された差を表すエラーファクタが計算される。
ステップ92で、ステップ90で計算されたファクタを使用して減衰マップが変更され、本方法は、次の投影データの生成、測定された透過データへの比較、差の決定及びrmsファクタの生成のためにステップ86に戻る。ステップ96で、最小のrmsエラーが発見されたかに関する判定が行われる。最小のrmsエラーが発見されていないと判定された場合、本方法は、ステップ92に戻る。最小のrmsエラーが発見されたと判定された場合、ステップ97で、ステップ84で生成された被検体の核イメージにおける減衰の補正で使用するため、改善された減衰マップが保存される。ステップ98で、ステップ84で生成された核イメージデータは、減衰が補正された核イメージ(たとえばPET又はSPECT)を生成するため、ステップ97で記憶されたACマップを使用して再構成される。
図5を参照して、例示的なホスピタルシステム100は、MR画像形成装置12、核スキャナ20(たとえばPET又はSPECT)等のような複数の画像形成装置を含み、これら複数の画像形成装置は、3D画像の表現を生成するため、個別の再構成プロセッサ又は共有される再構成プロセッサ102により再構成される画像形成データを生成する。画像の表現は、ネットワーク104から中央のメモリ106又はローカルのメモリ108に伝達される。
ネットワークに接続されるステーション110で、オペレータは、ユーザインタフェース124を使用して、選択された3Dの患者のMR減衰マップを中央のメモリ106及びローカルのメモリ108に移動するか、中央のメモリ106とローカルのメモリ108との間に移動する。ビデオプロセッサ116は、ディスプレイ122の第一の窓(view port)1181において選択された減衰マップ(又はMR画像)を表示する。核イメージは、第二の窓1182において表示される。第三の窓1183は、減衰マップと核イメージの重ね合わせを表示する。たとえば、ユーザは、MR減衰マップにおけるランドマークを対応する構造又は核イメージにおけるランドマークに位置決めすることができる。たとえば、オペレータは、インタフェース124を通して、減衰マップの画像におけるランドマークに対応する核イメージのランドマークを(たとえばマウス、スタイラス、又は他の適切なユーザ入力装置を使用して)選択する。代替的に、減衰マップは、プロセッサ116におけるプログラムにより自動的に位置決めすることができる。ユーザインタフェース124におけるプロセッサ36(図1)は、次いで補正アルゴリズムを実行し、減衰マップにおける不明確な領域において補うときに採用すべき適切な組織のタイプ(たとえば骨又は空気)を推測する。
次いで、他の応用で使用される場合がある、減衰が補正された核イメージを再構成するために洗練された減衰マップが使用される。たとえば、治療計画ステーション130は、治療のセッションを計画するために減衰が補正されたPET画像を使用する。ひとたびオペレータの満足に対して計画されると、計画された治療は、自動化された手順に対して必要に応じて、計画されたセッションを実現する治療装置132に転送することができる。他のステーションは、様々な他の計画の手順において減衰が補正されたPET画像を使用する場合がある。
別の実施の形態では、窓1183で表示されたオーバレイは、核イメージに対してMR画像に重み付けし、逆に、MR画像に対して核イメージを重み付けするために調節可能である。たとえば、メカニカル又はディスプレイ122に提示され、入力装置により操作される場合があるスライダバー又はノブ(図示せず)は、MR画像又は核イメージの重みを変化させるために調節される場合がある。1実施の形態では、オペレータは、(窓1181で示される)MR画像データから、MR及び核イメージデータの多数の及び/又は連続した組み合わせを通して、(窓1182に示される)核イメージデータへの画像を窓1183において調節する。たとえば、MR画像データの核画像データに対する比率は、0:1から1:0に不連続に又は連続して調節される。別の選択肢として、MR画像は、グレイスケールで表示され、核イメージは、カラー化することができる。MR画像における解剖学的ランドマークは、核イメージを被検体に関連付けるのに役立つ。
本発明は、幾つかの実施の形態を参照して記載された。先の詳細な説明を読んで理解することで、他の変更及び代替が行われる場合がある。本発明は、特許請求の範囲又はその等価なものの範囲に含まれる限り、全ての係る変更及び代替を含むとして解釈されるべきことが意図される。

Claims (15)

  1. 磁気共鳴(MR)取得スキャンの間に被検体のMR画像データを取得するMRイメージャと、
    核取得スキャンの間に前記被検体の核画像データを取得する核スキャナであって、前記核スキャナの検査領域に位置される放射線源からの透過データを測定する核スキャナと、
    前記MR画像データから減衰補正(AC)マップを生成し、前記測定された透過データを使用して、前記ACマップを繰り返し洗練させ、洗練されたACマップを生成し、前記洗練されたACマップを使用して、減衰について前記核画像データを補正するプロセッサと、
    を備える解剖学的な画像補正システム。
  2. 前記プロセッサは、
    前記ACマップにレイトレーシングアルゴリズムを実行して、推定された透過データを生成し、
    前記推定された透過データと前記測定された透過データとの間の差を決定する比較を行い、
    前記ACマップにおける減衰の値を調節して、前記推定された透過データと前記測定された透過データとの間の差を低減する、
    請求項1記載のシステム。
  3. 前記プロセッサは、最小自乗平均の誤差値を計算して、前記推定された透過データと前記測定された透過データとの間の差を最小にする、
    請求項2記載のシステム。
  4. 前記プロセッサは、前記差が最小にされたことの判定に応じて、前記洗練されたACマップをメモリに記憶する、
    請求項2記載のシステム。
  5. 前記核スキャナは、陽電子放出型断層撮影(PET)用のスキャナであり、
    前記PET用のスキャナは、前記被検体に投与される放射線同位体から放出される放射線であって、約511kevのエネルギーレベルを有する放射線からデータを取得し
    前記放射線源は、約511kevとは異なるエネルギーレベルを有する放射性物質である
    請求項記載のシステム。
  6. 前記放射線源は、固定された位置に設けられるか、又は、前記核スキャナのガントリの検査領域において前記被検体が位置される患者の支持体及び前記ガントリの少なくとも1つに設けられる、
    請求項1記載のシステム。
  7. 磁気共鳴(MR)減衰補正(AC)マップを洗練させる方法であって、
    被検体のMR画像データからACマップを生成するステップと、
    前記被検体を通して、前記被検体の放射線源から放射線を放出するステップと、
    前記被検体を通して透過された放射線から透過データを測定するステップと、
    前記ACマップから推定された透過データを生成するステップと、
    前記測定された透過データと前記推定された透過データとの比較に基づいて前記ACマップを調整することで、洗練されたACマップを生成するステップと、
    を含む方法。
  8. 陽電子放出型断層撮影(PET)用のデータを生成するステップと、
    洗練されたACマップにより前記PETデータを補正するステップと、
    補正されたPETデータをPET画像に再構成するステップと、
    を更に含む請求項記載の方法。
  9. 前記測定された透過データと前記推定された透過データとの間の第一の誤差を決定するステップと、
    前記ACマップを調整して前記第一の誤差を低減するステップと、
    調整されたACマップを使用して第二の推定された透過データを決定するステップと
    前記測定された透過データと前記第二の推定された透過データとの間の第二の誤差を決定するステップと
    前記調整されたACマップを更に調整して前記第二の誤差を低減するステップと
    を更に含む請求項記載の方法。
  10. 前記第二の誤差を決定するステップは、
    二乗平均平方根(rms)の値を繰り返し決定し、決定されたrms値が最適化されるまで前記ACマップを調整するステップを含む、
    請求項記載の方法。
  11. 1以上のコンピュータを制御して請求項10記載の方法を実行させるソフトウェアを記憶したコンピュータ読み取り可能な記録媒体。
  12. 陽電子放出型断層撮影(PET)用のスキャナは、前記被検体に投与される放射性同位体から放出される放射線であって、約511kevのエネルギーレベルを有する放射線からのデータを取得し、
    放射線源は、約511kevとは異なるエネルギーレベルを有する放射性物質である、
    請求項記載の方法。
  13. 陽電子放出型断層撮影(PET)用のスキャナであって、
    放射線検出器のリングと、
    前記放射線検出器のリング内の固定して設けられた放射線源と、
    洗練された減衰補正(AC)マップを生成するために請求項記載の方法を実行するプロセッサと、
    前記洗練されたACマップにより前記放射線検出器からのPETデータを補正し、補正されたPETデータを減衰が補正されたPET画像に再構成する再構成アルゴリズムと、
    を備えるPCT用のスキャナ。
  14. 単電子放出型コンピュータ断層撮影(SPECT)用のスキャナであって、
    複数の放射線検出器と、
    放射線検出器の検査領域内に固定して設けられた放射線源と、
    洗練された減衰補正(AC)マップを生成するために請求項記載の方法を実行するプログラムと、
    前記洗練されたACマップにより前記放射線検出器からのSPECTデータを補正し、補正されたSPECTデータを減衰が補正されたSPECT画像に再構成する再構成アルゴリズムと、
    を備えるSPECT用のスキャナ。
  15. 減衰補正(AC)マップを洗練させる方法であって、
    取得された磁気共鳴(MR)画像データから前記ACマップを生成するステップと、
    放射線のポイントソースから検出された透過データを使用して前記MR画像データにおける骨のボクセルと空気のボクセルとを識別するステップと、
    検出されたポイントソースの透過データに基づいて前記ACマップを更新するステップと、
    を含む方法。
JP2011526607A 2008-09-17 2009-09-09 ハイブリッド型の核/mr画像形成における透過データを使用したmrセグメンテーション Expired - Fee Related JP5385394B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US9755508P 2008-09-17 2008-09-17
US61/097,555 2008-09-17
PCT/IB2009/053944 WO2010032167A2 (en) 2008-09-17 2009-09-09 Mr segmentation using transmission data in hybrid nuclear/mr imaging

Publications (2)

Publication Number Publication Date
JP2012503177A JP2012503177A (ja) 2012-02-02
JP5385394B2 true JP5385394B2 (ja) 2014-01-08

Family

ID=42039963

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2011526607A Expired - Fee Related JP5385394B2 (ja) 2008-09-17 2009-09-09 ハイブリッド型の核/mr画像形成における透過データを使用したmrセグメンテーション

Country Status (6)

Country Link
US (1) US8406495B2 (ja)
EP (1) EP2338142A2 (ja)
JP (1) JP5385394B2 (ja)
CN (1) CN102265307B (ja)
RU (1) RU2504841C2 (ja)
WO (1) WO2010032167A2 (ja)

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8218848B2 (en) * 2008-07-23 2012-07-10 Siemens Aktiengesellschaft System and method for the generation of attenuation correction maps from MR images
DE102008047840B4 (de) * 2008-09-18 2018-08-30 Siemens Healthcare Gmbh Verfahren zur wenigstens teilweisen Ermittlung und/oder Anpassung einer zur Schwächungskorrektur von Positronenemissionstomographie-Bilddatensätzen genutzten Schwächungskarte in einem kombinierten Magnetresonanz-Positronenemissionstomographie-Gerät
WO2012120422A1 (en) * 2011-03-07 2012-09-13 Koninklijke Philips Electronics N.V. Mr segmentation using nuclear emission data in hybrid nuclear imaging/mr
CN103458790A (zh) * 2011-03-17 2013-12-18 皇家飞利浦有限公司 多模态心脏成像
CN103562960B (zh) * 2011-05-24 2017-06-09 皇家飞利浦有限公司 用于生成图像的图像区域与元素类之间的分配的装置
RU2014137323A (ru) 2012-02-16 2016-04-10 Конинклейке Филипс Н.В. Реконструкция радионуклидных изображений с пространственной коррекцией
US9204817B2 (en) * 2012-04-19 2015-12-08 General Electric Company Attenuation correction in positron emission tomography using magnetic resonance imaging
EP2847617B1 (en) * 2012-05-08 2019-06-12 Koninklijke Philips N.V. Spect/pet imaging system
KR101350496B1 (ko) * 2012-05-22 2014-01-16 서울대학교산학협력단 광자감쇠지도 생성 방법 및 mri 융합영상시스템
US8923592B2 (en) * 2012-05-29 2014-12-30 General Electric Company Methods and systems for performing attenuation correction
US8867814B2 (en) * 2012-10-04 2014-10-21 General Electric Company Methods and systems for generating a positron emission tomography attenuation correction map
US9002082B2 (en) 2012-12-27 2015-04-07 General Electric Company Axially varying truncation completion for MR-based attenuation correction for PET/MR
US9364192B2 (en) * 2013-03-15 2016-06-14 Siemens Medical Solutions Usa, Inc. Error estimates in quantitative functional imaging
RU2523929C1 (ru) * 2013-03-18 2014-07-27 Корпорация "САМСУНГ ЭЛЕКТРОНИКС Ко., Лтд." Система и способ для автоматического планирования видов в объемных изображениях мозга
US9311707B2 (en) * 2013-09-27 2016-04-12 General Electric Company System and method for attenuation correction of phantom images
WO2015189815A1 (en) * 2014-06-13 2015-12-17 Siemens Medical Solutions Usa, Inc. Multiple emission energies in single photon emission computed tomography
US20160066874A1 (en) * 2014-09-10 2016-03-10 The General Hospital Corporation Attenuation correction of positron emission tomography data using magnetic resonance images depicting bone density variations
JP6467221B2 (ja) * 2014-12-22 2019-02-06 キヤノン株式会社 画像処理装置および方法
US9508165B1 (en) * 2015-06-30 2016-11-29 General Electric Company Systems and methods for peak tracking and gain adjustment
CN108351396B (zh) * 2015-11-02 2021-09-14 皇家飞利浦有限公司 用于组织分类的方法、计算机程序产品和磁共振成像系统
US10210634B2 (en) 2016-07-20 2019-02-19 Shanghai United Imaging Healthcare Co., Ltd. System and method for segmenting medical image
KR102419590B1 (ko) * 2018-03-07 2022-07-08 지멘스 메디컬 솔루션즈 유에스에이, 인크. 가압 가스 이온 챔버-기반 선량 교정기에서의 교정 바이어스 감소
CN109308728B (zh) * 2018-10-25 2023-01-03 上海联影医疗科技股份有限公司 正电子发射型计算机断层扫描图像处理方法及装置

Family Cites Families (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5376795A (en) * 1990-07-09 1994-12-27 Regents Of The University Of California Emission-transmission imaging system using single energy and dual energy transmission and radionuclide emission data
US5449913A (en) 1993-11-03 1995-09-12 Chang; Wei Apparatus for producing attenuation scan
US5672877A (en) * 1996-03-27 1997-09-30 Adac Laboratories Coregistration of multi-modality data in a medical imaging system
EP1099123A1 (en) 1998-04-27 2001-05-16 Duke University Transmission scanning technique for gamma camera coincidence imaging
US6310968B1 (en) * 1998-11-24 2001-10-30 Picker International, Inc. Source-assisted attenuation correction for emission computed tomography
US7652259B2 (en) * 2000-08-21 2010-01-26 Spectrum Dynamics Llc Apparatus and methods for imaging and attenuation correction
US6628983B1 (en) * 2000-10-25 2003-09-30 Koninklijke Philips Electronics N.V. Nuclear imaging systems and methods with feature-enhanced transmission imaging
TW542993B (en) 2001-07-12 2003-07-21 Inst Information Industry Multi-dimension and multi-algorithm document classifying method and system
DE10231061A1 (de) * 2002-07-10 2004-01-22 Philips Intellectual Property & Standards Gmbh Verfahren und System zur Verbesserung des Informationsgehaltes in einem Bild
US7197171B2 (en) * 2003-02-19 2007-03-27 Elgems Ltd. Nuclear imaging
GB0414277D0 (en) * 2004-06-25 2004-07-28 Leuven K U Res & Dev Orthognatic surgery
US20060004274A1 (en) * 2004-06-30 2006-01-05 Hawman Eric G Fusing nuclear medical images with a second imaging modality
US8108024B2 (en) * 2004-12-15 2012-01-31 Koninklijke Philips Electronics N.V. Registration of multi-modality images
US9468777B2 (en) * 2005-10-17 2016-10-18 Alberta Health Services Integrated external beam radiotherapy and MRI system
US7402807B2 (en) * 2005-11-02 2008-07-22 Siemens Medical Solutions Usa, Inc. Method for reducing an electronic time coincidence window in positron emission tomography
US7348564B2 (en) 2005-12-12 2008-03-25 General Electric Company Multi modality imaging methods and apparatus
US7465927B2 (en) * 2006-03-23 2008-12-16 Siemens Medical Solutions Usa, Inc. Attenuation correction for nuclear medical imaging scanners with simultaneous transmission and emission acquisition
DE102006046287A1 (de) * 2006-09-29 2008-04-03 Siemens Ag Felderzeugungseinheit einer kombinierten MR-PET-Anlage
WO2008064319A2 (en) 2006-11-22 2008-05-29 The General Hospital Corporation Attenuation correction of pet image using image data acquired with an mri system
JP4584232B2 (ja) * 2006-11-27 2010-11-17 富士フイルムRiファーマ株式会社 脳模型
DE102007034953B4 (de) * 2007-07-26 2016-09-22 Siemens Healthcare Gmbh Verfahren zur Bewegungsvorgänge berücksichtigenden Aufnahme von Messdaten eines Patienten und zugehörige medizinische Einrichtung
US8098916B2 (en) * 2007-10-30 2012-01-17 General Electric Company System and method for image-based attenuation correction of PET/SPECT images

Also Published As

Publication number Publication date
US20110164801A1 (en) 2011-07-07
RU2011115056A (ru) 2012-10-27
CN102265307A (zh) 2011-11-30
WO2010032167A3 (en) 2011-04-21
US8406495B2 (en) 2013-03-26
CN102265307B (zh) 2014-10-22
RU2504841C2 (ru) 2014-01-20
EP2338142A2 (en) 2011-06-29
WO2010032167A2 (en) 2010-03-25
JP2012503177A (ja) 2012-02-02

Similar Documents

Publication Publication Date Title
JP5385394B2 (ja) ハイブリッド型の核/mr画像形成における透過データを使用したmrセグメンテーション
CN109741284B (zh) 用于校正pet成像中由呼吸运动引起的失配的系统和方法
CN102711617B (zh) 扫描规划视场调整器、确定器和/或质量评估器
JP5462865B2 (ja) 不完全な解剖学的画像を補償するための非減衰補正pet放出画像の使用
CN103315760B (zh) 用于核医疗成像中基于发射数据进行衰减补偿的系统和方法
CN106456098B (zh) 衰减图的生成方法及系统
CN107635469B (zh) 基于pet系统中的散射符合的衰减图的估计
CN111080706B (zh) 重建图像数据集的方法、装置、计算机程序和数据载体
EP2245592B1 (en) Image registration alignment metric
CN106999135B (zh) 辐射发射成像系统和方法
US20220313176A1 (en) Artificial Intelligence Training with Multiple Pulsed X-ray Source-in-motion Tomosynthesis Imaging System
CN109124666A (zh) 一种确定辐射剂量调制线的方法、系统和装置
JP6855476B2 (ja) 組織分類の方法、コンピュータプログラム、及び磁気共鳴イメージングシステム
KR102283934B1 (ko) Pet 영상에 기반한 가상 ct 영상 및 감쇠보정 pet 영상 생성 방법 및 시스템
CN109077746A (zh) 一种确定辐射剂量调制线的方法、系统和装置
KR101350496B1 (ko) 광자감쇠지도 생성 방법 및 mri 융합영상시스템
CN111050652A (zh) 关于图像引导的应用采用谱(多能量)图像数据
CN111316327A (zh) 对移动对象的pet数据的衰减校正
US20230056685A1 (en) Methods and apparatus for deep learning based image attenuation correction
CN117853703A (zh) 介入物识别方法、成像系统及非暂态计算机可读介质

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20120906

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: 20130910

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20131003

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees