JP2005525205A - 局所的な重み付けされた補間を含む拡散テンソル磁気共鳴イメージング - Google Patents
局所的な重み付けされた補間を含む拡散テンソル磁気共鳴イメージング Download PDFInfo
- Publication number
- JP2005525205A JP2005525205A JP2004505720A JP2004505720A JP2005525205A JP 2005525205 A JP2005525205 A JP 2005525205A JP 2004505720 A JP2004505720 A JP 2004505720A JP 2004505720 A JP2004505720 A JP 2004505720A JP 2005525205 A JP2005525205 A JP 2005525205A
- Authority
- JP
- Japan
- Prior art keywords
- voxel
- representation
- voxels
- diffusion
- fiber
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- 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
- G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
- G01R33/56341—Diffusion imaging
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Signal Processing (AREA)
- Radiology & Medical Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Vascular Medicine (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Theoretical Computer Science (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Generation (AREA)
Abstract
線維構造又は異方性構造(102)を含む被検者(16)をイメージングするためのイメージング方法(150、190)は、ある異方性構造(102)をもつ領域の3次元の見かけの拡散テンソルマップ(162)を取得することを含む。ボクセルの見かけの拡散テンソルは、固有ベクトル(e1、e2、e3)及び固有値(λ1、λ2、λ3)を得るように処理される(164)。3次元線維表現(208)は、固有ベクトル及び固有値を使用して、抽出される(190)。抽出(190)中、ボクセルは、線維表現(208)の近傍の少なくとも選択された次元において局所的に補間される(202)。補間は、局所的な異方性を示すパラメータによって、ボクセルを重み付けすることを含む。補間は、取得されたテンソルマップ(162)より高い追跡精度及び表現解像度をもつ3次元線維表現(208)をもたらす。
Description
本発明は、3次元イメージング技術に関する。本発明は特に、拡散テンソル磁気共鳴イメージング(DT−MRI)による神経線維及び線維束のイメージング、追跡及び表示に関し、これを特に参照して説明される。しかしながら、本発明は更に、他のタイプの線維構造の追跡及びグラフィックレンダリング並びにシングルフォトンエミッションコンピューテッドトモグラフィイメージング(SPECT)、コンピューテッドトモグラフィ(CT)、ポジトロンエミッショントモグラフィ(PET)、などの他のイメージングモダリティと関連するアプリケーションも見出す。
人間及び他の哺乳類の神経組織は、神経線維又は線維束を形成するように配置された細長い軸索部分を有するニューロンを含み、電気化学的信号は、神経線維又は線維束に沿って伝送される。脳において、例えば、非常に高い神経密度によって規定される機能エリアは、一般に、軸索の線維束の構造的に複雑なニューラルネットワークによってリンクされる。軸索の線維束及び他の線維材料は、実質的に他の組織によって囲まれている。
神経疾患の診断、脳外科手術の計画、他の神経学に関連した臨床的アクティビティ及び脳機能に関する研究は、軸索線維及び線維束の非侵襲的なイメージング及び追跡から利益を得ることができる。特に、拡散テンソル磁気共鳴イメージング(DT−MRI)が、軸索の線維束と相関する画像コントラストを提供するために示されている。DT−MRI技法では、磁気共鳴画像が水又は他の流体分子の拡散に関連するコントラストを含むように、拡散鋭敏化(diffusion-sensitizing)磁場勾配が、励起/イメージングシーケンスにおいて印加される。励起/イメージングシーケンス中、選択された方向に拡散勾配を印加することによって、拡散強調画像が取得され、それらの画像から、像空間のそれぞれのボクセル位置ごとに見かけの拡散テンソル係数が得られる。
流体分子は、線維に対し部分的に又は完全に直交する方向と比べて、軸索の線維束の方向に沿ってより速くより簡単に拡散する。従って、見かけの拡散係数の指向性及び異方性は、軸索の線維及び線維束の方向と相関する傾向がある。
DT−MRI画像からの線維構造情報の抽出は、計算集約的であり、処理時間は一般に、臨床的に価値がある画像、ボリューム、モデル又はパラメータのために、数十分から1時間にも及ぶ。合理的な制限内に処理時間を保つとともに信号対雑音比を改善するために、通常、比較的低い解像度の画像が、線維構造の追跡のために取得される。例えば、30−40のスライスのボリューム、約2mmの線形の再構成されたボクセル寸法を有する128x128の各ボクセル、及び2mmのスライス解像度、が一般的である。この比較的粗い解像度は、表示される描画内の追跡された線維表現の外観に影響を及ぼす。線維構造の追跡は一般に3次元的に実施されるので、2のファクタによる取得された画像解像度の増加は、およそ8のファクタによる計算時間及びメモリ使用の増加に対応する。
本発明は、上述した制限その他を克服する改善された装置及び方法を企図する。
本発明の1つの見地によれば、異方性構造を含む被検者をイメージングするイメージング方法が提供される。少なくともいくつかの異方性構造を含む被検者の少なくとも一部の3次元の見かけの拡散テンソルマップが取得される。見かけの拡散テンソルは、固有ベクトル及び固有値を得るために、ボクセルで処理される。3次元線維表現は、固有ベクトル及び固有値を使用して抽出される。抽出中、ボクセルは、線維表現の近傍の少なくとも選択された次元において局所的に補間される。補間は、局所的な異方性を示すパラメータによって、補間されるボクセルを重み付けすることを含む。補間は、取得されたテンソルマップより高い解像度を有する3次元線維表現をもたらす。人間が見ることができる3次元線維表現の表示が生成される。
本発明の別の見地によれば、被検者中の線維構造を追跡する装置が開示される。拡散テンソルマップを取得する手段が提供される。ボクセルに対応する拡散テンソルの順序付けられた(ordered)固有値及び固有ベクトルを決定する手段が提供される。ボクセル固有値及び固有ベクトル並びに少なくとも1つの開始ボクセル選択に基づいて線維構造表現を計算する追跡手段が、提供される。追跡手段と協働し、線維構造表現の計算中に線維構造の近傍で局所的にボクセルを補間することによって線維構造表現の解像度を高める補間手段が提供される。人間が見ることができる媒体に、線維構造表現の少なくとも一部を表示する表示手段が提供される。
本発明の1つの利点は、追跡された線維表現の改善された滑らかさにある。
本発明の別の利点は、改善された計算スピード及び低減されたメモリ使用にある。
本発明の更に別の利点は、より高い解像度を有する改善された追跡精度にある。
本発明の多くの更なる利点及び利益は、以下の好ましい実施例の詳細な説明を読むことにより当業者に明らかになるであろう。
本発明は、さまざまな構成要素及び構成要素の組み合わせ方並びにさまざまなステップ及びステップの組み合わせ方の形をとることができる。図面は、好ましい実施例を説明することのみを目的とし、本発明を制限するものとして解釈されるべきではない。
図1を参照して、磁気共鳴イメージング(MRI)スキャナ10は、一般に、検査領域14を通るz軸に沿って、実質的に均一であり時間的にコンスタントな主磁界B0を生成する超伝導又は常伝導磁石12を有する。ボア(孔)タイプの磁石が図1に示されているが、本発明は、オープンマグネットシステム及び他のタイプのMRIスキャナにも等しく適用可能である。イメージングは、イメージングされている被検者、例えば患者18が少なくとも部分的に検査領域14内に配置されるとともに、一般に関心領域が磁石のアイソセンタにある状態で、磁気共鳴励起及び読み取りシーケンスを実行することによって行われる。脳領域の拡散テンソルMRIイメージングの場合、患者頭部が、図示するようにアイソセンタに配置されることが好ましい。
磁気共鳴シーケンスは、磁気スピンを反転させ又は励起し、磁気共鳴を引き起こし、磁気共鳴を再び集中させ、磁気共鳴を操作し、空間的に符号化し、磁気共鳴を流れるように起こさせ又は拡散符号化し、スピンを飽和させる等のために、被検者16に印加される一連のRF及び磁場勾配パルスを含む。より具体的には、勾配パルス増幅器20が、検査領域14のx軸、y軸及びz軸に沿って磁場勾配を生成するために、電流パルスを全身勾配コイルアセンブリ22に印加する。拡散強調MRIにおいて、選択された磁場勾配が、選択された方向における流体分子の動きに敏感な動き関連コントラストを提供するために印加される。
好適にはデジタルであるRF送信器24は、検査領域にRFパルスを送信するために、RFパルス又はパルスパケットを全身RFコイル26に印加する。一般的なRFパルスは、互いに一緒に利用される短い持続時間のすぐ隣りのパルスセグメントのパケットからなり、印加されるいかなる勾配も、選択された磁気共鳴操作を達成する。RFパルスは、飽和させ、共鳴を励起し、磁化を反転し、共鳴を再び集中させ又は検査領域の選択された部分において共鳴を操作するために使用される。
全身アプリケーションの場合、選択された操作の結果として生成される共鳴信号は、全身RFコイル26によっても拾われる。代わって、被検者の制限された領域にRFパルスを生成する場合、局所的なRFコイルが、選択された領域に隣接して配置される。例えば脳イメージングの場合、挿入可能な頭部コイル28が任意に用いられる。
RFコイル構成及びそのアプリケーションに関係なく、RFコイルの或るもの又は別のものによって拾われる、結果として得られるRF磁気共鳴信号は、RF受信器30、好適にはデジタル受信器によって受け取られ、復調される。シーケンス制御プロセッサ32は、勾配パルス増幅器20、RF送信器24及びRF受信器30を制御して、統合されたMRIパルスシーケンスを生成し、磁気共鳴信号及び任意のエコーを生成する波形を読み出し、結果として得られるMR応答を空間的に拡散符号化するために適当な符号化勾配を印加し、MRピックアップをコーディネートし、オペレーションを受け取る。
MRIシーケンスは、一般に、勾配増幅器22によって生成される複雑な一連の磁場勾配パルス及び/又はスウィープを含み、これらはRFコイル26、28によって送信される選択されたRFパルスとともに、k空間にマップする磁気共鳴エコーを生じさせる。結果として得られる磁気共鳴データは、ソータ34によって分類され、k空間メモリ36に記憶される。拡散テンソル磁気共鳴イメージング(DT−MRI)において、拡散強調なしの0とインデクスを付されたデータ及びN方向に拡散強調をもつ1...Nとインデクスを付されたデータが取得される。好適には、N≧6である。静止の及び拡散強調データが取得されて、対応するk空間メモリ360,361,...,36Nに分類される。静止画像及びそれぞれの拡散強調画像用のk空間データ360,361,...,36Nは、再構成された画像表現S0,S1,...SN、400,401,...40Nを生成するために、一般に逆フーリエ変換プロセッサ又は当技術分野に知られている他の再構成プロセッサである再構成プロセッサ38によって処理される。一般に、画像S0400は、拡散強調なしで得られる。
DT−MRIの場合、選択された方向の見かけの拡散係数(ADC)は、異なる拡散強調をもつ画像の線形回帰によって得られる。選択された方向(i,j)に関するADCは、(i,j)方向に拡散強調された拡散強調画像Si,jを取得することによって決定される。拡散強調画像Si,jは、下式に従って、強調されない画像S0400に関連する。
ここで、Bは、磁界パラメータであり、ADCi,jは、(i,j)方向に関する見かけの拡散係数である。ADCは、方程式(1)から以下のように計算される。
有利には、T1、T2、T2 *のような拡散強調以外のメカニズム及び同様の画像コントラストメカニズムから結果的に得られる画像コントラストは、線形回帰によって実質的に除去される。拡散テンソルプロセッサ42は、拡散テンソルマップを構成するために、方程式(2)によって各ボクセルごとに複数のADC値を計算する。一般に、6つの拡散方向が、各ボクセルにおける拡散テンソルを構成するのに十分な情報を提供する。
図1を引き続き参照し、図2を更に参照して、固有ベクトル/固有値算出プロセッサ46が、ボクセルにおいて拡散テンソル固有ベクトル及び固有値を順序付ける。図2に示すように、拡散テンソルの順序付けられた固有値λ1、λ2、λ3(最大値から最小値まで順序付けられる)及び対応する固有ベクトルe1、e2、e3は、有用な物理的重要性をもつ。最大固有値は、図2にλ1として示される。対応する固有ベクトルe1は、主(major)固有ベクトルと呼ばれ、最も高い拡散係数を有する空間方向とアラインする。残りの固有値λ2、λ3は、図2において、中間の(medium)及び小さい(minor)固有ベクトルと呼ばれる対応する固有ベクトルe2、e3を有する。これらの固有ベクトルe2、e3は、e1と直交し、より低い拡散係数を有する空間方向とアラインする。固有値λ1、λ2、λ3の相対値は、拡散テンソルの異方性の空間的な向き及び大きさを示す。
引き続き図2を参照して、固有ベクトル及び固有値は、楕円体100によって幾何学的に表現可能である。楕円体100の長軸は、固有ベクトルe1、すなわち最も高い拡散係数もつ方向とアラインする。楕円体100の完全な球からの偏りは、拡散テンソルの異方性を表わす。異方性拡散係数テンソルは、例えば固有ベクトルe2、e3の方向のように線維102に対して部分的に又は完全に直交する方向において拡散を呈する傾向がある神経線維束102の影響を反映することができる。対照的に、線維102と平行な拡散、すなわち主固有ベクトルe1の方向に沿った方向の拡散は、e2、e3方向に沿ったものよりも強調されており、大きい。
図1に戻って、当技術分野に知られている異方性比率(fractional anisotropy)マップのような異方性マップ50又は異方性の大きさを重要視する他の異方性画像マップが、順序付けられた固有ベクトル及び固有値から任意に計算される。適切な実施例において、異方性マップが、次式に従って各ボクセルごとに計算される。
ここで、
方程式(3)及び(4)の異方性マップは、拡散線維追跡のための線維の領域を選択する場合に使用するのに特に適していることが分かった。当業者に公知であるように、異方性マップ50は、異方性の大きさの便利な画像表現を提供する。しかしながら、異方性マップ50は方向情報を含まず、従って、例えば脳又は脊椎のような神経学的に複雑な領域において交差し、分岐し、からまり合い又は他の方法で重なり合う多数の線維及び線維束から、選択された軸索の線維束を追跡し又はセグメント化するには適さない。
図1を引き続き参照して、線維追跡/補間プロセッサ52は、線維表現メモリ54に記憶される線維表現を生成するために、拡散テンソルマップ44の中の軸索の線維又は線維束をセグメント化し又は追跡する。線維追跡/補間プロセッサ52は、関連するユーザによって選択される1つ又は複数のシードボクセルから追跡を始める。
適切な選択プロセスにおいて、ユーザは、グラフィックディスプレイ60、キーボード62、マウス又はその他のポインティングデバイス64及び/又は同様の入力又は出力素子を有するパーソナルコンピュータ又はワークステーションのようなユーザインタフェース装置58を操作する。ディスプレイ60は、例えばスライス抽出器68によって抽出される、拡散強調画像、S0画像400又は異方性マップ50のうちの1つからのスライス66又は他の2次元表面のような、被検者16の画像表現を示す。ディスプレイ60は更に、1つ又は複数のシードボクセルを選択するために、ポインティングデバイス64、キーボード62又は他の入力装置を通してユーザによって操作可能である、スーパーインポーズされる対話的なポインタ69又は他のグラフィック選択ツールを示す。
図3は、異方性が楕円体100によって幾何学的に表されている、例示の2次元拡散係数固有ベクトル/固有値マップ110を示す。引き続き図1を参照するとともに、再び図2を参照し、図3を更に参照して、追跡/補間プロセッサ52は、シードボクセル及びそれぞれの後続のボクセルの拡散テンソルの主固有ベクトルe1によって与えられる追跡方向を使用して、線維を反復的に追跡する。図3は、シードボクセル1001で始まる第1の線維表現トラック112と、シードボクセル1002で始まる第2の線維表現トラック114とを示す。
図3には一方向の追跡が示されているが、追跡は、より完全な線維追跡を与えるために、双方向に、すなわちe1及び−e1の両方の方向に任意に行われることができることが分かるであろう。更に、例示としてボクセルの2次元アレイが図3に示されているが、追跡は、一般に、3次元拡散テンソルマップ44を使用して3次元で行われることが分かるであろう。
図1を引き続き参照して、方向の追跡に加えて、線維追跡/補間プロセッサ52は更に、固有ベクトル/固有値算出プロセッサ46と共にテンソルボクセルの局所的な補間を実施する。補間は、計算時間及びメモリローディングを低減するために追跡最前部付近の領域に局所化される。イメージングスライス平面内の二分する補間は、追跡された線維表現の有意な滑らかさを与えることが分かった。追跡精度を改善するために、補間されたボクセルは、例えば近くのボクセルの固有ベクトル及び固有値を平均化し、又は補間されるボクセルを、近くのボクセルの主固有ベクトルの方向の加重平均に対応する拡散方向に割り当てることによって、拡散係数テンソル異方性に関連したパラメータによって重み付けされることが好ましい。
線維表現54は、グラフィック表示装置60又は別の出力装置上のユーザが見ることができる表示のために、レンダリングプロセッサ70によってグラフィカルにフォーマットされる。線維表現54は、異方性マップ50又はその一部、静止画像400、拡散強調画像401,...,40Nの1つ、一様な黒又は白のスクリーン、その他のような背景に、ライン、チューブ、ワイヤフレーム又は他のグラフィック表現としてレンダリングされる。好適には、レンダリングプロセッサ70は、3次元描画を生成し、その3次元描画は、ユーザによってマウス64、キーボード62又は他の入力装置を介して回転され、大きさを変更され又は他の方法で操作されることができる。
図4を参照して、拡散テンソル磁気共鳴イメージング(DT−MRI)データを取得し、関連する空間的マッピング及び画像表現を生成するための適切な方法150が記載されている。被検者は、ステップ152にて、磁気共鳴イメージング(MRI)スキャナ内に位置させられ、固定される。DT−MRIイメージングは、相当な時間期間にわたって多数の画像を取得することを含むので、被検者は、モーションブラー(動きのブレ)及び動きにより引き起こされる画像の位置ずれを最小限にするために固定されることが好ましい。
DT−MRIイメージングデータは、ステップ154にて、例えば選択された拡散強調を生成する付加の磁場勾配パルスを含むスピンエコーシーケンスのようなイメージングシーケンスを使用して取得される。好適には、拡散係数テンソルの選択された見かけの拡散係数(ADC)成分に対応するいくつかの拡散強調156によりデータが取得される複数エコーシーケンスが使用される。概して、6つの見かけの拡散係数が、テンソルを記述するのに十分である。例示の実施例において、強調されない画像(0,0,0)とともに、(x,0,0)、(0,y,0)、(0,0,z)、(x,−y,0)、(x,0,−z)及び(0,y,−z)方向に印加される磁場勾配パルスによる6つの拡散強調156が収集される。しかしながら、代わりに、拡散強調の他の組み合わせが使用されることができる。有利には、複数エコーシーケンスを使用することは、データ取得時間を短縮するとともに、画像間の動きにより引き起こされるブラー又は位置ずれを最小限にする。信号対雑音比を改善するために、複数の画像に関するデータが、それぞれの拡散強調ごとに収集されることが好ましい。イメージングシーケンスは更に、磁場勾配によって引き起こされる渦電流及び他のイメージングアーチファクトを補償するために、付加のRFパルス又は磁場勾配パルス若しくはスイープを任意に含む。
ステップ154にて収集された画像データは、ステップ158にて、静止及び拡散強調画像再構成S0及びSijkを形成するために再構成される。ここで、ijkはさまざまな拡散強調156を示す。当技術分野に知られている逆フーリエ変換の再構成が適切に使用されるが、他の再構成方法が用いられることもできる。
拡散強調画像が取得され再構成されると、それぞれのボクセルの見かけの拡散係数(ADC)が、線形回帰又は別の技法を使用して方程式(2)によって計算され、見かけの拡散係数テンソルマップが、ステップ162にて構成される。固有値及び固有ベクトル166が、ステップ164にて抽出される。
任意に、ステップ168にて、異方性マップが、例えば方程式(3)及び(4)によって計算される。好適には、異方性マップは、関連するユーザに対して表示するためのカラー化された異方性画像172を得るために、例えば異方性値に基づいてボクセルをカラー化することによってステップ170にてレンダリングされる。異方性画像172は、ユーザが線維追跡のための関心領域を選択するのに便利な媒体を提供する。
図5を参照して、DT−MRI画像の線維を追跡するための適切な方法190が記載されている。ステップ192にて、ユーザは、開始の関心領域を選択する。適切な実施例において、選択192は、図4のDT−MRIイメージング方法150において任意に得られる異方性画像172を参照して行われる。好適には、選択された開始の関心領域は、ユーザによって、マウスポインタ又は他のグラフィック選択装置を使用してグラフィカルに示される。選択された関心領域は、単一のボクセル、ボクセルの平面領域又はボクセルの3次元領域でありえる。ユーザは、任意に、終了の関心領域を選択することもできる。このような選択は、一般に、深部白質の2つの機能的脳領域の間に延在する線維を追跡するのに役立つ。
選択された開始の関心領域内の開始ボクセルが、ステップ196にて選択される。このボクセルから始まって、局所的の方向は、ステップ198にて、例えば固有値を順位付け、最大の又は主固有値λ1を識別し、対応する主固有ベクトルe1の方向を識別することによって、主固有ベクトルe1方向(図2を参照)に対応して識別される。ステップ200にて、局所的の方向に沿って現在ボクセルの近傍の次のボクセルが識別される(図3を参照)。好ましい実施例において、正および負の両方(双方向)の追跡が、正及び負の局所的の拡散(e1)方向における次のボクセルを識別することによって実施される。追跡が双方向に進むにつれて、正の線維末端は、正の局所的な方向にボクセルを連続的に識別することによって成長され、負の線維末端は、負の局所的な方向にボクセルを連続的に識別することによって成長される。更に、深部白質の大きく密度の高い領域から伸びる線維を追跡するような特定の状況においては、一方向の線維追跡も企図される。
ステップ202にて、局所化された補間が、現在の1つ又は複数の追跡最前部の近傍で、すなわち次の1つ又は複数のボクセルの近傍で実施される。好適には、補間は、計算時間及びメモリローディングを低減するために追跡最前部(又は双方向追跡の場合は複数の追跡最前部)付近の領域に局所化される。好ましい実施例において、補間は、1つ又は複数の追跡最前部がより高い解像度の局所的な拡散テンソル空間内で進むように、固有ベクトル及び固有値を得る間にテンソルにおいて実施される。結果として得られる追跡された線維表現は、ソースのテンソルマップのより高い補間された局所的な解像度に対応する改善された解像度及び滑らかさを有する。
補間されたより高い解像度の局所的な拡散テンソル空間における追跡精度を維持するために、補間されるボクセルは、局所的な異方性に関連するパラメータによって重み付けされることが好ましい。例えば、補間されるボクセルの固有ベクトル及び固有値は、補間されるボクセルが局所的な異方性を反映するように、近傍のボクセルの固有ベクトル及び固有値の平均として適切に選択される。しかしながら、複数の固有値及び固有ベクトルを重み付けすることは計算的に高価である。線維追跡スピードを改善するために、補間されたボクセルの主固有ベクトルが、隣り合うボクセルの主固有ベクトル方向の加重平均に対応する方向をもつように選択することが十分であることが分かった。選択された主固有ベクトルは、任意に、異方性比率又は別の異方性パラメータに対応する大きさの加重を含む。
現在の好ましい実施例において、補間は、例えば軸取得スライスから構成される画像の場合には軸平面の向きのように、一般には取得された画像スライスの平面の向きに対応する、選択された平面の向きにおいて行われる。しかしながら、原則の矢状面、冠状面及び体軸の平面の向き以外の平面の向きで補間することも企図される。例えば、補間は、任意に、主固有ベクトルe1及び中間固有ベクトルe2(図2を参照)を含む平面の向きにおいて行われることができ、それによって、補間は、線維方向に対応する見込みが最も高い平面に沿って行われる。この企図された実施例において、平面の向きは、成長する線維末端の局所的の異方性の向きに関してシフトし、回転する。3次元補間を実施することもまた企図される。例えば、線形補間は、それぞれの補間された方向の平面解像度を2倍にすることができ、すなわち、128x128ボクセル平面は、256x256ボクセル解像度を与えるように補間される。同様の補間が、スライス間(inter-slice)方向において任意に実施される。
方法190は、一方向に又は双方向に追跡される線維を成長させるために、ステップ198、200、202を反復的に繰り返す。好適には、繰り返しループ内の決定ステップ204が、進行する線維末端の終端について調べる。1つの適切な線維終端基準は、選択された値より低い、例えば方程式(3)において使用されるFA=0.10の閾値以下の異方性比率を含む。低い異方性比率(FA)は、高い等方性の拡散テンソルに対応するので、選択された閾値より低い異方性パラメータを、追跡された線維の終点に関連付けることは合理的である。
適切な別の線維終端基準は、選択された角度より大きい、連続するボクセル間の局所的な方向の変化である。局所的な方向の大きい変化は、線維の分岐点に対応する見込みがある。任意に、このような見かけの分岐点で終端する代わりに、新しい関心領域が、分岐点に又はその周辺に規定され、追跡プロセス170が、分岐する線維又は線維束を追跡するために新しい関心領域について繰り返される。
更に別の適切な線維終端基準は、進行する線維が、ステップ192にてユーザによって選択された終わりの関心領域に入ることである。本願明細書に記載される例示の終端基準のさまざまな組み合わせ及び/又は他の適切な終端基準の付加若しくは置き換えもまた企図される。
決定ステップ204が、線維が終端したことを示す場合、方法は、好適には、関心領域のすべてのボクセルが決定ステップ206において訪れられたかどうかを調べる。訪れられていないボクセルがまだ関心領域に存在する場合、別の開始ボクセルがステップ196にて選択され、追跡ステップ198、200、202、204が反復的に繰り返される。関心領域のそれぞれの開始ボクセルで始まる反復的な追跡は、集合的に、1つ又は複数の追跡された線維又は線維束を表わす3次元の追跡された線維表現208を生成する。好適には、追跡された線維表現は、例えばワイヤメッシュ3次元レンダリング、管状3次元レンダリング、他のタイプのレンダリングを使用して、人間のユーザによる視覚化のために、ステップ210にてグラフィカルにレンダリングされ、表示される。ユーザが、選択された線維を強調表示し、3次元線維表現を回転させ又は他の方法でデータを操ることができるように、レンダリングは対話的であることが好ましい。
本発明は、脳の軸索の線維追跡に関して説明された。しかしながら、当業者であれば、本発明が、実質的な線維又は異方性の組織密度が存在する、脊椎、腎臓及び心臓のような他の解剖学的構造にも適用できることが分かるであろう。
Claims (20)
- 異方性構造を含む被検者をイメージングするためのイメージング方法であって、
少なくともいくつかの異方性構造を含む前記被検者の少なくとも一部の、3次元の見かけの拡散テンソルマップを取得するステップと、
固有ベクトル及び固有値を得るために、ボクセルにおける前記見かけの拡散テンソルを処理するステップと、
前記固有ベクトル及び前記固有値を使用して、3次元線維表現を抽出する抽出ステップと、
前記抽出中、前記線維表現の近傍の少なくとも選択された次元においてボクセルを局所的に補間するステップであって、前記補間が、局所的な異方性を示すパラメータによって、補間されたボクセルを重み付けすることを含み、前記補間が、取得された前記テンソルマップより高い解像度を有する3次元線維表現をもたらす、ステップと、
前記3次元線維表現の、人間が見ることができる表示を生成するステップと、
を含むイメージング方法。 - 前記局所的な異方性を示す前記パラメータが、
主拡散テンソル固有ベクトルの方向、
拡散テンソルの機能的な組み合わせ、
固有値、
固有ベクトル、及び、
異方性比率、
からなるグループから選択される、請求項1に記載のイメージング方法。 - 前記3次元線維表現を抽出する前記ステップは、
前記見かけの拡散テンソルマップ内の開始領域の選択を受け取るステップと、
前記開始領域内のボクセルを選択するステップと、
前記選択されたボクセルの拡散テンソル固有ベクトル及び固有値に基づいて、局所的な方向を決定する決定ステップと、
前記局所的の方向に沿って次のボクセルを識別するステップと、
前記3次元線維表現を抽出するために、前記局所的方向の決定及び前記次のボクセルの識別を複数回繰り返すステップと、
を含む、請求項1又は請求項2に記載のイメージング方法。 - 前記次のボクセルを識別する前記ステップは、補間されるボクセルを識別することを含む、請求項3に記載のイメージング方法。
- 前記線維表現の近傍の少なくとも選択された次元においてボクセルを局所的に補間する前記ステップが、
前記選択されたボクセルと、前記局所的な方向に沿って前記選択されたボクセルと隣り合う別のボクセルと、の間に少なくとも1つのボクセルを補間すること、
を含む、請求項3又は請求項4に記載のイメージング方法。 - 前記3次元線維表現を抽出する前記ステップが、
次のボクセルの異方性尺度が選択された閾値より小さいこと、
前記次のボクセルの局所的な方向の変化が選択された角度より大きいこと、及び、
前記次のボクセルが選択された終わりの関心領域に入ること、
のうちの1つに応じて、前記繰り返すステップを終えることを更に含む、請求項3乃至請求項5のいずれか1項に記載のイメージング方法。 - 前記開始領域の選択を受け取る前記ステップが、
前記3次元の見かけの拡散テンソルマップの少なくとも一部を表わす画像表現の、人間が見ることができる選択表示を生成するステップと、
前記開始領域を選択するために、関連するユーザによって操作可能なグラフィック選択ツールを前記選択表示にスーパーインポーズするステップと、
を含む、請求項3乃至請求項6のいずれか1項に記載のイメージング方法。 - 前記選択表示が、
磁気共鳴画像と、
拡散テンソル磁気共鳴イメージングから得られる異方性画像と、
からなるグループから選択される、請求項7に記載のイメージング方法。 - 前記3次元線維表現を抽出する前記ステップが、
開始ボクセルを選択するステップと、
前記3次元線維表現を構成するために、前記開始ボクセルから始まって、最も大きい固有値に対応する固有ベクトルをボクセルからボクセルへと反復的に追跡するステップと、
を含む、請求項1又は請求項2に記載のイメージング方法。 - 前記ボクセルを局所的に補間する前記ステップが、局所的に補間された前記ボクセルを含むより高い解像度の局所的なテンソル空間を生成することを含み、前記反復的に追跡するステップが、より高い解像度の局所的に補間されたテンソル空間内で行われる、請求項9に記載のイメージング方法。
- 前記3次元の拡散テンソルマップを取得する前記ステップが、
拡散強調磁気共鳴画像表現を取得するステップと、
取得された前記拡散強調磁気共鳴画像表現から、前記拡散テンソルマップを計算するステップと、
を含む、請求項1乃至請求項10のいずれか1項に記載のイメージング方法。 - 前記拡散テンソルマップを計算する前記ステップが、線形回帰によって、前記取得された拡散強調磁気共鳴画像表現から、見かけの拡散テンソル成分を決定することを含む、請求項11に記載のイメージング方法。
- 被検者の線維構造を追跡する装置であって、
拡散テンソルマップを取得する手段と、
ボクセルに対応する前記拡散テンソルの順序付けられた固有値及び固有ベクトルを決定する手段と、
ボクセル固有値及び固有ベクトル並びに少なくとも1つの開始ボクセル選択に基づいて、線維構造表現を計算する追跡手段と、
前記線維構造表現の計算中、前記線維構造の近傍のボクセルを局所的に補間することによって、前記追跡手段と協働して、前記線維構造表現の解像度を高める補間手段と、
人間が見ることができる媒体に、前記線維構造表現の少なくとも一部を表示する表示手段と、
を有する装置。 - 異方性大きさパラメータに対応するボクセルをもつ異方性マップを計算する異方性マッピング手段を更に有し、前記異方性マップが、前記少なくとも1つの開始ボクセル選択のために使用される、請求項13に記載の装置。
- 前記追跡手段が、
前記選択されたボクセルの最も大きい固有値と関連する固有ベクトルの方向に対応する局所的な線維方向を識別する方向計算手段と、
前記局所的な線維方向に沿って前記選択されたボクセルの近傍のボクセルを識別するインクリメント手段と、
前記線維構造表現を反復的に計算するために、前記方向計算手段及び前記インクリメント手段を反復的に起動させるルーピング手段と、
を有する、請求項13又は請求項14に記載の装置。 - 前記補間手段が、
前記選択されたボクセルの近傍の前記拡散テンソルマップの隣り合うボクセルの間に位置する補間されるボクセル位置を識別する位置決め手段と、
補間されるボクセル値を得るために、前記選択されたボクセル及び前記隣り合うボクセルのうち少なくとも1つの見かけの拡散テンソルを、重み付けを伴って組み合わせるバリュー計算手段と、
を有する、請求項15に記載の装置。 - 前記補間手段は、2次元表面及び前記選択されたボクセルを含むボリュームのうちの1つにおいて、局所的にボクセルを補間する、請求項15又は請求項16に記載の装置。
- 前記補間手段が、
前記選択されたボクセルの近傍の前記見かけの拡散テンソルマップの隣り合うボクセルの間で、補間されるボクセル位置を識別する位置決め手段と、
補間されるボクセルについて重み付けされた方向値を得るために、隣り合うボクセルの拡散方向を重み付けを伴って組み合わせる方向値計算手段と、
を有する、請求項15乃至17のいずれか1項に記載の装置。 - 前記拡散テンソルマップを取得する前記手段が、
拡散強調イメージングデータを取得するイメージング手段と、
取得された前記拡散強調イメージングデータを拡散強調画像表現に再構成する再構成手段と、
選択された拡散強調画像表現を選択的に組み合わせることによって、拡散テンソルマップを構成する手段と、
を有する、請求項13乃至請求項18のいずれか1項に記載の装置。 - 前記拡散テンソルマップを構成する前記手段が、前記選択された拡散強調画像から見かけの拡散テンソル係数を計算するために、線形回帰を実施する手段を有する、請求項19に記載の装置。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US10/145,916 US6806705B2 (en) | 2002-05-15 | 2002-05-15 | Diffusion tensor magnetic resonance imaging including local weighted interpolation |
PCT/IB2003/001935 WO2003098251A1 (en) | 2002-05-15 | 2003-05-12 | Diffusion tensor magnetic resonance imaging including local weighted interpolation |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2005525205A true JP2005525205A (ja) | 2005-08-25 |
Family
ID=29418703
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2004505720A Pending JP2005525205A (ja) | 2002-05-15 | 2003-05-12 | 局所的な重み付けされた補間を含む拡散テンソル磁気共鳴イメージング |
Country Status (6)
Country | Link |
---|---|
US (1) | US6806705B2 (ja) |
EP (1) | EP1506425A1 (ja) |
JP (1) | JP2005525205A (ja) |
CN (1) | CN1653349A (ja) |
AU (1) | AU2003224381A1 (ja) |
WO (1) | WO2003098251A1 (ja) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009524070A (ja) * | 2006-01-24 | 2009-06-25 | デスモンド フィッツジェラルド アンド アソシエイツ ピーティーワイ エルティディー | 複数の観察されたテンソル間の改良された補間方法 |
JP2010512174A (ja) * | 2006-12-11 | 2010-04-22 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 巨視的情報に基づくファイバ追跡 |
JP2014068786A (ja) * | 2012-09-28 | 2014-04-21 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴装置 |
JP2019209018A (ja) * | 2018-06-07 | 2019-12-12 | 洋司 西島 | 脳mri検査による認知症の早期診断 |
Families Citing this family (34)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2002338376A1 (en) * | 2001-04-06 | 2002-10-21 | Lawrence R. Frank | Method for analyzing mri diffusion data |
CN1777818B (zh) * | 2003-04-24 | 2010-10-06 | 皇家飞利浦电子股份有限公司 | 纤维跟踪磁共振成像 |
US7355403B2 (en) * | 2005-01-27 | 2008-04-08 | Siemens Medical Solutions Usa, Inc. | Noise reduction in diffusion tensor imaging data using bayesian methods |
US8077938B2 (en) * | 2005-02-09 | 2011-12-13 | Siemens Medical Solutions Usa, Inc. | Diffusion tensor surface visualization |
US7532214B2 (en) * | 2005-05-25 | 2009-05-12 | Spectra Ab | Automated medical image visualization using volume rendering with local histograms |
US20070165989A1 (en) * | 2005-11-30 | 2007-07-19 | Luis Serra Del Molino | Method and systems for diffusion tensor imaging |
US7411393B2 (en) * | 2005-11-30 | 2008-08-12 | Bracco Imaging S.P.A. | Method and system for fiber tracking |
US7889899B2 (en) * | 2006-05-05 | 2011-02-15 | Siemens Medical Solutions Usa, Inc. | Directional statistics visualization tool |
US8041129B2 (en) | 2006-05-16 | 2011-10-18 | Sectra Ab | Image data set compression based on viewing parameters for storing medical image data from multidimensional data sets, related systems, methods and computer products |
JP4972751B2 (ja) * | 2006-11-27 | 2012-07-11 | 株式会社日立製作所 | 神経線維束計測システム及び画像処理システム |
US7355407B1 (en) * | 2006-12-03 | 2008-04-08 | Toshiba Medical Systems Corp. | Methods and apparatus for single-shot magnetic resonance imaging with optimized isotropic diffusion weighting |
US7830381B2 (en) * | 2006-12-21 | 2010-11-09 | Sectra Ab | Systems for visualizing images using explicit quality prioritization of a feature(s) in multidimensional image data sets, related methods and computer products |
JP2009050615A (ja) * | 2007-08-29 | 2009-03-12 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴イメージング装置および磁気共鳴画像表示方法 |
US20090158309A1 (en) * | 2007-12-12 | 2009-06-18 | Hankyu Moon | Method and system for media audience measurement and spatial extrapolation based on site, display, crowd, and viewership characterization |
US8340376B2 (en) | 2008-03-12 | 2012-12-25 | Medtronic Navigation, Inc. | Diffusion tensor imaging confidence analysis |
US8643363B2 (en) * | 2008-11-27 | 2014-02-04 | Syntheticmr Ab | Method of visualizing segmented MR images with absolute-scale values independent of MR scanner settings |
JP5525867B2 (ja) * | 2009-03-04 | 2014-06-18 | 株式会社東芝 | 超音波診断装置、画像処理装置、超音波診断装置の制御方法、及び画像処理方法 |
US9720062B2 (en) | 2010-03-03 | 2017-08-01 | Brain Research Institute Foundation Pty Ltd. | Image processing system |
US8447089B2 (en) * | 2010-04-27 | 2013-05-21 | Chunlei Liu | Systems and methods for susceptibility tensor imaging |
CN102054281B (zh) * | 2010-12-29 | 2012-12-26 | 中国科学院深圳先进技术研究院 | 弥散加权成像方法及系统 |
WO2012125829A2 (en) * | 2011-03-15 | 2012-09-20 | University Of Pittsburgh - Of The Commonwealth System Of Higher Education | Directional diffusion fiber tracking |
US9285449B2 (en) | 2011-06-15 | 2016-03-15 | Chunlei Liu | Systems and methods for imaging and quantifying tissue magnetism with magnetic resonance imaging |
CN103479356B (zh) * | 2012-06-14 | 2015-07-29 | 西门子(深圳)磁共振有限公司 | 一种扩散张量磁共振成像方法 |
US20160061923A1 (en) * | 2013-04-03 | 2016-03-03 | Choukri Mekkaoui | Sheet tractography using diffusion tensor mri |
CN103445780B (zh) * | 2013-07-26 | 2015-10-07 | 浙江工业大学 | 一种扩散加权磁共振成像多纤维重建方法 |
US10712418B2 (en) | 2015-05-15 | 2020-07-14 | The Medical College Of Wisconsin, Inc. | Systems and methods for diffusion-weighted multi-spectral magnetic resonance imaging |
DE102016207910A1 (de) | 2016-05-09 | 2017-11-09 | Siemens Healthcare Gmbh | Aufnahme von diffusionsgewichteten Magnetresonanzbilddaten mit anisotropen Diffusionsrichtungen |
JP7218288B2 (ja) | 2016-11-09 | 2023-02-06 | ランダム ウォーク イメージング アーベー | サンプルに拡散強調磁気共鳴測定を実行する方法 |
FR3060179B1 (fr) * | 2016-12-08 | 2020-04-03 | Safran | Procede et dispositif de determination des orientations d'elements de fibres dans une piece en materiau composite |
CN110021003B (zh) * | 2019-02-14 | 2021-02-02 | 清华大学 | 图像处理方法、图像处理装置和核磁共振成像设备 |
CN109799472B (zh) * | 2019-03-08 | 2022-06-28 | 沈阳工业大学 | 一种基于深度学习的磁共振涡流补偿方法 |
US10593112B1 (en) * | 2019-04-15 | 2020-03-17 | Biosense Webster (Israel) Ltd. | Chamber reconstruction from a partial volume |
CN110473206B (zh) * | 2019-07-24 | 2023-03-21 | 东南大学 | 一种基于超体素与测度学习的弥散张量图像分割方法 |
CN110599560B (zh) * | 2019-08-05 | 2023-07-25 | 上海联影医疗科技股份有限公司 | 磁共振成像方法、装置、存储介质及计算机设备 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1993018415A1 (en) * | 1992-03-09 | 1993-09-16 | University Of Washington | Image neurography and diffusion anisotropy imaging |
WO1995004940A1 (en) | 1993-08-06 | 1995-02-16 | The Government Of The United States Of America, Represented By The Secretary Of The Department Of Health And Human Services | Method and system for measuring the diffusion tensor and for diffusion tension imaging |
US5488297A (en) | 1994-07-19 | 1996-01-30 | The Regents Of The University Of California | Three dimensional anisotropy contrast magnetic resonance imaging |
US5969524A (en) * | 1997-04-14 | 1999-10-19 | The United States Of America As Represented By The Department Of Health And Human Services | Method to significantly reduce bias and variance of diffusion anisotrophy measurements |
JP3538726B2 (ja) * | 1997-06-13 | 2004-06-14 | 株式会社日立メディコ | 拡散テンソルの拡散係数を表示する磁気共鳴装置 |
JPH11100320A (ja) | 1997-09-29 | 1999-04-13 | Risuburan:Kk | 皮膚外用剤 |
US6526305B1 (en) * | 1998-11-25 | 2003-02-25 | The Johns Hopkins University | Method of fiber reconstruction employing data acquired by magnetic resonance imaging |
-
2002
- 2002-05-15 US US10/145,916 patent/US6806705B2/en not_active Expired - Fee Related
-
2003
- 2003-05-12 EP EP03720807A patent/EP1506425A1/en not_active Withdrawn
- 2003-05-12 WO PCT/IB2003/001935 patent/WO2003098251A1/en not_active Application Discontinuation
- 2003-05-12 JP JP2004505720A patent/JP2005525205A/ja active Pending
- 2003-05-12 AU AU2003224381A patent/AU2003224381A1/en not_active Abandoned
- 2003-05-12 CN CNA038108208A patent/CN1653349A/zh active Pending
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009524070A (ja) * | 2006-01-24 | 2009-06-25 | デスモンド フィッツジェラルド アンド アソシエイツ ピーティーワイ エルティディー | 複数の観察されたテンソル間の改良された補間方法 |
JP2010512174A (ja) * | 2006-12-11 | 2010-04-22 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 巨視的情報に基づくファイバ追跡 |
JP2014068786A (ja) * | 2012-09-28 | 2014-04-21 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴装置 |
JP2019209018A (ja) * | 2018-06-07 | 2019-12-12 | 洋司 西島 | 脳mri検査による認知症の早期診断 |
JP7177429B2 (ja) | 2018-06-07 | 2022-11-24 | 洋司 西島 | 脳mri検査による認知症の早期診断 |
Also Published As
Publication number | Publication date |
---|---|
EP1506425A1 (en) | 2005-02-16 |
CN1653349A (zh) | 2005-08-10 |
WO2003098251A1 (en) | 2003-11-27 |
AU2003224381A1 (en) | 2003-12-02 |
US6806705B2 (en) | 2004-10-19 |
US20030214289A1 (en) | 2003-11-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6806705B2 (en) | Diffusion tensor magnetic resonance imaging including local weighted interpolation | |
JP4399353B2 (ja) | リアルタイムトラクトグラフィ | |
US6642716B1 (en) | Diffusion tensor magnetic resonance imaging including fiber rendering using hyperstreamlines | |
JP4414879B2 (ja) | Dtisnrを改善するための画像の遡及的な選択及び様々なタイプの画像アラインメント | |
JP5123191B2 (ja) | 高度に限定された画像再構成法を使用する拡散テンソル・イメージング | |
US6845260B2 (en) | Automatic vessel indentification for angiographic screening | |
US7949170B2 (en) | Image processing method, image processing device, computer aided detection, and method for filtering along the time axis | |
US7034531B1 (en) | Diffusion MRI using spherical shell sampling | |
US20120108946A1 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
JP6979151B2 (ja) | 磁気共鳴イメージング装置及び磁気共鳴画像処理方法 | |
JPWO2002056767A1 (ja) | 高精度コイル感度マップを用いたパラレルmrイメージング | |
JP4070268B2 (ja) | 画像の分解能を向上させる方法、投影画像の一部を拡大させる方法及び画像の分解能を向上させる装置 | |
EP1269212A1 (en) | Magnetic resonance angiography with automated vessel segmentation | |
US20150363951A1 (en) | Systems and methods for objective tractographic processing using features of grid structures of the brain | |
KR101253024B1 (ko) | 자기공명영상을 이용한 삼차원 영상 복원장치 및 방법 | |
US10132902B2 (en) | Intrinsic navigation from velocity-encoding gradients in phase-contrast MRI | |
JP2002000581A (ja) | 磁気共鳴画像化を使用して得られた複素画像の合成 | |
WO2010045003A1 (en) | Method for time-of-arrival mapping in magnetic resonance imaging | |
JP6629247B2 (ja) | 微細構造解析データのデータ品質を評価し向上させる方法 | |
US20210212567A1 (en) | Determination of parametrized characteristics of a tissue | |
Achterberg et al. | Optimal short-time acquisition schemes in high angular resolution diffusion-weighted imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060511 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090226 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20090721 |