JP4280018B2 - X-ray computed tomography system - Google Patents
X-ray computed tomography system Download PDFInfo
- Publication number
- JP4280018B2 JP4280018B2 JP2002087047A JP2002087047A JP4280018B2 JP 4280018 B2 JP4280018 B2 JP 4280018B2 JP 2002087047 A JP2002087047 A JP 2002087047A JP 2002087047 A JP2002087047 A JP 2002087047A JP 4280018 B2 JP4280018 B2 JP 4280018B2
- Authority
- JP
- Japan
- Prior art keywords
- data
- ray
- detection element
- field
- view
- 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 - Lifetime
Links
- 238000002591 computed tomography Methods 0.000 title claims description 18
- 238000001514 detection method Methods 0.000 claims description 57
- 238000003384 imaging method Methods 0.000 claims description 31
- 238000013213 extrapolation Methods 0.000 claims description 4
- 230000000007 visual effect Effects 0.000 claims 2
- 238000003325 tomography Methods 0.000 claims 1
- 238000000034 method Methods 0.000 description 35
- 238000010586 diagram Methods 0.000 description 10
- 238000013500 data storage Methods 0.000 description 8
- 230000002093 peripheral effect Effects 0.000 description 7
- 239000011295 pitch Substances 0.000 description 6
- 210000004072 lung Anatomy 0.000 description 5
- 238000007781 pre-processing Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 230000015556 catabolic process Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 239000000470 constituent Substances 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 238000006731 degradation reaction Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000009966 trimming Methods 0.000 description 3
- BUGBHKTXTAQXES-UHFFFAOYSA-N Selenium Chemical compound [Se] BUGBHKTXTAQXES-UHFFFAOYSA-N 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 230000006866 deterioration Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 229910052711 selenium Inorganic materials 0.000 description 2
- 239000011669 selenium Substances 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 239000006096 absorbing agent Substances 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Description
【0001】
【発明の属する技術分野】
本発明は、立体情報を得るために角錐形状のX線ビームで被検体をスキャンするいわゆるコーンビーム型X線コンピュータ断層撮影装置に関する。
【0002】
【従来の技術】
コーンビームスキャン方式では、X線管から放射され、X線コリメータで角錐形状にトリムされたX線で被検体をスキャンする。被検体を透過したX線は、2次元アレイタイプの検出器で検出される。このタイプのX線検出器は、ライン検出器を比較的少数、典型的には、4列並べたものが普及している。しかし、最近では、シンチレータ素子とフォトダイオード素子を組み合わせた固体検出素子、又はX線を電荷に直接的に変換するセレン等の固体検出素子の採用により、32列、さらにそれを越える列数を持ったX線検出器が登場している。このX線検出器の外形は、円筒形、平面形等様々である。
【0003】
コーンビーム画像再構成法には、フェルドカンプ法(Feldkamp method)が一般的に用いられる。フェルドカンプ法は、ファンビームコンボリューション・バックプロジェクション法をもとにした近似的再構成法であり、コンボリューション処理は、コーン角が比較的小さいことを前提として、データをファン投影データと見なして行われる。しかし、バックプロジェクション処理は、実際のレイ(ray)に沿って行われる。
【0004】
つまり、
(1)投影データにZ座標に依存した重みをかけ、
(2)(1)のデータに、ファンビーム再構成と同じ再構成関数をコンボリューションし、
(3)(2)のデータを、コーン角を持つ斜めの実際のレイに沿って逆投影する、
という手順で画像が再構成される。
【0005】
このような画像再構成法は、撮影視野(field of view)の有効長が、撮影視野の径に依存して変化されるという問題を有している。以下に、具体的に説明する。
【0006】
図1には比較的長いRLLに径が設定された撮影視野を側面から示し、図2には、比較的短いRSSに径が設定された撮影視野を示している。撮影視野の径Rは、頭部、肺、胴体等の検査部位が収まる長さに設定される。なお、撮影視野の有効長とは、設定された径Rが維持される撮影視野のスライス方向に関する長さ(被検体の体軸方向に関する撮影視野の長さ)を、“有効長”と定義する。
【0007】
フェルドカンプ法では、1回転分のデータが必要とされるので、画像再構成可能な最大範囲は、円柱形状に限定される。この範囲内で径Rが維持される撮影視野の有効長は、径がRLLの場合には、図3Aに示すように、WLLに限定され、また、径がRSSの場合には、図3Bに示すように、WSSに限定される。このように撮影視野の有効長は、撮影視野の設定される径に対応して変化してしまう。
【0008】
【発明が解決しようとする課題】
本発明の目的は、コーンビーム型X線コンピュータ断層撮影装置において、撮影視野の径に対する有効長の依存性を緩和することにある。
【0009】
【課題を解決するための手段】
本発明のある局面によるX線コンピュータ断層撮影装置は、X線を角錐形状で被検体に照射するように構成されたX線管装置と、前記被検体を透過したX線を検出するスライス方向とチャンネル方向との2方向に関して少なくともスライス方向にアレイされた複数の検出素子を有する検出器と、前記検出素子により検出された実データに基づいて、前記検出素子がアレイされた領域から前記スライス方向に関して外側に位置する拡張領域に対応する仮想データを発生するように構成されたデータ拡張部と、前記実データと、前記仮想データとに基づいて、画像データを再構成するように構成された再構成装置と、撮影視野の径を入力する入力装置と、前記入力された撮影視野の径に基づいて、前記拡張領域の長さを決定するように構成された拡張領域長決定部とを具備する。
【0010】
【発明の実施の形態】
以下、図面を参照して本発明によるX線コンピュータ断層撮影装置を好ましい実施例により説明する。なお、X線コンピュータ断層撮影装置には、X線管とX線検出器とが1体として被検体の周囲を回転するROTATE/ROTATE-TYPE、リング状にアレイされた多数の検出素子が固定され、X線管のみが被検体の周囲を回転するSTATIONARY/ROTATE-TYPE等様々なタイプがあり、いずれのタイプでも本発明を適用可能である。ここでは、ROTATE/ROTATE-TYPEとして説明する。
【0011】
また、画像データ(断層画像データ)を再構成するには、被検体の周囲1周、360°分の投影データが、またハーフスキャン法でも180°+ファン角分の投影データが必要とされる。いずれの再構成方式にも本発明は適用可能である。ここでは、360°法を例に説明する。
【0012】
また、入射X線を電荷に変換するメカニズムは、シンチレータ等の蛍光体でX線を光に変換し更にその光をフォトダイオード等の光電変換素子で電荷に変換する間接変換形と、X線によるセレン等の半導体内での電子正孔対の生成及びその電極への移動すなわち光導電現象を利用した直接変換形とが主流である。検出素子としては、それらのいずれの方式を採用してもよい。
【0013】
また、近年では、X線管とX線検出器との複数のペアを回転リングに搭載したいわゆる多管球型のX線CT装置の製品化が進み、その周辺技術の開発が進んでいる。本発明では、従来からの一管球型のX線CT装置であっても、多管球型のX線CT装置であってもいずれにも適用可能である。ここでは、一管球型として説明する。
【0014】
図4は本実施例に係るX線コンピュータ断層撮影装置の構成を示す図である。ガントリ100には、回転軸Oを中心として回転自在に支持された回転リング102が収容されている。回転リング102には、X線管装置101が搭載される。X線管装置101は、X線管と、トリミング装置とを有する。X線管は、高電圧発生装置109から管電圧の印加及び管電流の供給を受けて、焦点からX線(X-rays)を放射する。トリミング装置は、X線管からのX線を四角形にトリムする。このトリミングにより、X線は、角錐形状に成形される。
【0015】
回転リング102には、X線管装置101とともに、2次元アレイ型のX線検出器103が搭載されている。X線検出器103は、回転軸Oを挟んでX線管装置101に正対する一及びアングルで取り付けられる。X線検出器103は、図2に示すように、複数の検出素子108を有している。複数の検出素子108は、回転軸Oに平行な方向(スライス方向)と、回転軸Oに直交し、且つX線焦点を中心として緩やかにカーブする方向(チャンネル方向)との2方向に関して2次元状に配列される。このような2次元の検出器103は、チャンネル方向に関して1列に配列された検出素子108を、スライス方向に関して複数列並べることで構成するようにしてもよいし、検出素子108がM×Nのマトリクスに配列されてなる複数のモジュールを配列することで構成するようにしてもよい。
【0016】
撮影に際しては、X線管装置101とX線検出器103との間に、被検体が配置される。ヘリカルスキャン時には、被検体とガントリ100の相対的な位置が一定速度で変位される。
【0017】
X線検出器103の出力には、一般的にDAS(Data Acquisition System) と呼ばれるデータ収集システム104が接続されている。このデータ収集システム104には、X線検出器103の各素子の電流信号を電圧に変換するI−V変換器と、この電圧信号をX線の曝射周期に同期して周期的に積分する積分器と、この積分器の出力信号を増幅するアンプと、このプリアンプの出力信号をディジタル信号に変換するアナログ・ディジタル・コンバータとが、チャンネルごとに設けられている。
【0018】
このデータ収集システム104から出力されるデータ(純生データ(pure raw data))は、スリップリング又は非接触信号伝送装置を経由して、前処理装置106に伝送される。前処理装置106は、この純生データに対して前処理を施す。前処理には、例えばチャンネル間の感度不均一補正処理、X線強吸収体、主に金属部による極端な信号強度の低下又は信号脱落を補正する処理等が含まれる。この前処理装置106から出力されるデータ(生データ(raw data))は、磁気ディスク、光磁気ディスク、又は半導体メモリを備えたデータストレージ装置111に記憶される。
【0019】
GUI制御部117は、表示装置116の画面に、グラフィカルユーザインターフェース(GUI)を表示する。グラフィカルユーザインターフェースは、スキャン条件及び再構成条件等の複数の設定項目に対応付けされたアイコン、ボタン、プルダウンメニュー等のグラフィカルな要素を含んでいる。このグラフィカルな要素をマウス等のポインティング・デバイス(入力装置)115により操作することにより、様々な操作を簡易に実現することができるようになっている。なお、再構成条件には、画像データ(ここではボリュームデータ)を再構成する領域を表す略円柱形の撮影視野の大きさが含まれる。撮影視野の大きさは、径Rと長さWとにより定義される。この撮影視野の大きさ、つまり径Rと長さWの入力は、グラフィカルな要素により簡易化されている。
【0020】
拡張領域長決定部112は、入力装置115を介して設定された撮影視野の径Rが維持される撮影視野の有効長を所定長に固定するために、X線検出器103の検出素子108がアレイされている実在のX線有感領域を、スライス方向に関して外側に仮想的に拡張する長さを、入力装置115を介して設定された撮影視野の径Rに基づいて、決定する。なお、実在の有感領域に対して、拡張された仮想的な有感領域を、拡張領域と称する。
【0021】
拡張領域の長さは、撮影視野の径Rに基づいて演算される。または、撮影視野の複数の径に対して、拡張領域の複数の長さをそれぞれ関連付けたテーブルを予め作成し、データストレージ装置111に記憶させておき、このテーブルから、入力装置115を介して設定された撮影視野の径Rに関連付けられている拡張領域の長さを読み出すようにしてもよい。
【0022】
この拡張領域の長さは、詳細は後述するが、操作者により設定された径Rが維持される撮影視野の有効長Wが、任意に設定される撮影視野の径Rの大小に依存しないで、一定又は略一定に固定することを目的として決定される。
【0023】
投影データ拡張部113は、拡張領域長決定部112で決定された長さを有する拡張領域に、実在の検出素子と同じ密度で仮想的にアレイした複数の仮想検出素子に対して設定した仮想的な複数のレイ(射線)のデータ(仮想データ)を、データストレージ装置111に保持されているX線検出器103の検出素子で実測された生データ(実データ)に基づいて、作成する。なお、レイとは、X線管装置101のX線焦点から検出器103の検出素子の中心に引かれる直線として定義され、このレイに沿って逆投影処理が行われる。
【0024】
再構成装置114は、X線管装置101が360°又は180°+ファン角の範囲内の実データと、それから作成された同じ範囲の仮想データとに基づいて、拡張フェルドカンプ再構成法に従って、設定された径Rで上述した一定の長さWを持つ略円柱形の撮影視野の画像データ(ボリュームデータ)を再構成する。表示装置116では、この画像データに基づいて任意断面、3次元のレンダリング画像を作成し、表示する。
【0025】
なお、逆投影処理において、再構成装置114では、厳密には、X線焦点とボクセルの中心とを結ぶ直線(計算上のレイ)に沿って逆投影が行われる。しかし、実際には、上述したように、X線焦点と検出素子中心とを結ぶレイ(実際のレイ)に沿ってX線投影が行われている。計算上のレイと、実際のレイとの空間的な若干のずれは、画質の低下を引き起こしてしまう。
【0026】
再構成装置114では、この計算上のレイと、実際のレイとの空間的な若干のずれを解消して画質低下を軽減するために、実データ及び仮想データに対してずれ補正処理を行う。ずれ補正処理の詳細は、特開平09−19425号公報、米国特許5,825,842に記載されているので、ここでは簡単に説明する。ここであるボクセルに対する逆投影を考える。X線焦点とボクセルの中心とを結んだ直線を延長し、有感領域の平面に交差する点を点Cとする。点Cは(n,m)の検出素子、(n,m+1)の検出素子、(n+1,m)の検出素子、(n+1,m+1)の検出素子各々の中心点間に存在するものとする。この点Cの近傍の複数、ここでは4つの検出素子のデータから距離補間により点Cのデータを推定する。この推定したデータを使って逆投影処理を行うことにより、上述した計算上のレイと実際のレイとの空間的なずれによる画質低下を軽減することができる。
【0027】
図6には、拡張領域長決定部112により撮影視野の比較的長い径RLLに応じて決定される拡張領域の長さ(拡張幅)EDLLを示し、図7には比較的短い径RSSに応じて決定される拡張領域の拡張幅EDSSを示している。また、図8Aは図6の拡張幅EDLLに応じて決まる撮影視野の有効長WLLを示し、図8Bは拡張幅EDSSに応じて決まる撮影視野の有効長WSSを示している。なお、図6、図7では、理解を容易にするために、角度位置が0°にあるときのX線管装置101及び検出器103とともに、それに対向する180°の角度位置にあるX線管装置及び検出器をそれぞれ101´、103´として示している。また、図6、図7において、実データレイを実線で、仮想データレイを破線で示している。また、実データに対応する撮影視野内の領域を斜線で、また仮想データに対応する撮影視野内の領域を網線で示している。
【0028】
仮想データを作成する目的は、撮影視野の径Rに対するその径が維持される有効長Wの依存性を緩和することにある。つまり撮影視野の径Rが様々に設定されたとしても、その径Rが維持される撮影視野の有効長Wが変動しない又はその変動が微小な範囲に抑えられることを実現することにある。この目的達成の為に、検出器103のスライス方向に関して実在有感領域の外側に拡張領域を付加し、換言すると、有感領域をスライス方向に仮想的に拡大し、この拡張領域の仮想データに最外又はその近傍の検出素子の実データをあてる、又は仮想データを、最外の検出素子の実データと近傍の検出素子の実データとから外挿補間により作成し、しかも有効長が所定長に固定されるように、拡張領域の長さWを、撮影視野の径Rに応じて変化させる。
【0029】
図6及び図8Aに示すように、撮影視野の径Rが比較的長い径RLLに設定されたとき、拡張領域の長さはEDLLに決定される。この拡張幅EDLLにより、撮影視野の有効長はWLLに決まる。もちろん、この撮影視野の有効長WLLは、同じ径RLLの条件もとで検出器103の実在の有感領域に依存して決まる有効長WRLLより長い。
【0030】
また、図7及び図8Bに示すように、撮影視野の径Rが比較的短い径RSSに設定されたとき、拡張領域の長さはEDSSに決定される。この拡張幅EDSSにより、撮影視野の有効長はWSSに決まる。
【0031】
拡張領域の長さEDLLは、拡張領域の長さEDSSよりも長く決定されている。それぞれの長さは、最終的な撮影視野の有効長Wが、等しく又は略等しくなる、換言すると最終的な撮影視野の有効長Wが所定長又はほぼ所定長に固定されるように、X線管装置101のX線焦点から検出器103の検出面までの最短距離、撮影視野の径R、撮影視野の所望の有効長W、および検出器103の実在の有感領域の長さに基づいて幾何学的に決定される。
【0032】
拡張領域は検出器103の実在の有感領域外に設定されているのであるから、その拡張領域上の仮想データは、実測されない。従って、作り出す必要がある。ここでは、作成効率と画質劣化抑制とがバランスした2種類の方法を提供する。いずれの方法を採用しても良いし、両方法を装備してユーザ指示に従って選択的に使用しても良いし、両方法を使用して2種類の画像を再構成して最終的にユーザが選択するようにしも良い。
【0033】
その一方の方法は、図9、図10に示すように、仮想データの仮想検出素子とスライス方向に関して同じライン上であって、仮想データの仮想検出素子に最も近い実在の検出素子で検出した実データ、つまり実在の検出素子の中でスライス方向に関して最外(outermost)に位置する検出素子で検出した実データを、そのまま仮想データにあてる方法である。
【0034】
これは実際には、実データを保持しているデータストレージ装置111から再構成装置114へのデータの読み出し制御により実現され得る。つまり、データ拡張部113は、仮想データを発生するに際して、データストレージ装置111を、最外に位置する検出素子で検出した実データと同じアドレスでアクセスし、その最外に位置する検出素子で検出した実データを当該仮想検出素子の仮想データとして再構成装置114へ読み出させる。
【0035】
この方法では、拡張領域上においてスライス方向に並んでいる複数の仮想検出素子の仮想データは、同じ実データで統一される。実際にその仮想データが影響するのは、図6、図7に網線で示したスライスの周辺部分である。実際の検査上では、スライス中心付近に関心部位が配置されるケースが多く、スライスの周辺部分は中心付近よりも重要視されない。また、スライス中心付近に関心部位を配置したとき、スライスの周辺部分は、比較的簡単な組織構造であることが解剖学的には多い。従って、周辺部分の画質低下はそれほど問題にはならない。この変形例としては、最外列の検出素子のデータの代わりに、最外列から一列内側の検出素子の実データを仮想データに割り当てても良い。
【0036】
他の方法は、仮想データを、その仮想検出素子に近傍する複数の実在検出素子で検出された実データから外挿補間により作成する方法である。具体的には、例えば、図9、図10に示すように、仮想データの仮想検出素子に最も近い位置、つまり最外に位置する実在検出素子で検出された実データd1と、その最外の実在検出素子より1つ内側の実在検出素子の実データd2とから、それぞれの実在の検出素子と仮想検出素子との間の検出面上での距離S1、S2に基づいて計算される。もちろん、最外を含めて近傍の3以上の実在県素子の実データを使って補間するようにしてもよい。この方法でも、先の方法より多少改善されるものの、やはり、多少画質(ここでは画像再現性、信頼性を意味する)は低下するが、その影響が強く現れるのがスライスの周辺部分であり、その周辺部分の画質低下は実際上それほど問題にはならないと想定される。この変形列としては、仮想検出素子に近傍する複数の実在検出素子で検出された実データの平均値等を仮想データとして作成しても良い。
【0037】
このように検出器103のスライス方向に関して実在の有感領域の外側に拡張領域を設定し、この拡張領域の投影データを、実データから作成し、しかもその拡張領域の長さWを、撮影視野の径Rに応じて変化させることにより、撮影視野の径Rに対するその有効長Wの依存性が解消又は緩和される。
【0038】
図11に撮影視野の径Rと拡張領域の長さEDとの関係をシュミレーションした結果を示している。なお、図11では拡張領域の長さEDを、拡張する仮想的な検出素子列の列数として表している。撮影視野の直径2Rが約250mm以下では、ほとんど拡張領域を設定しなくてよい。これは実データだけで、所定の有効長Wが達成されることを理由とする。撮影視野の直径2Rが約250mmを超えたあたりから、撮影視野の有効長Wの短縮を抑えるために、拡張領域が設定され、その拡張領域の拡張幅EDは撮影視野の径Rの増大に対してほぼ比例する関係で増大される。
【0039】
また、ヘリカルスキャンでは、撮影視野の径Rに応じて撮影視野の有効長が変動する不具合を抑えるために、1回転あたり寝台の天板が移動する単位距離、つまりヘリカルピッチを変化させる必要があった。図12に太線で示すように、本方法を使うことで、ヘリカルピッチの変動を抑制することができる。
【0040】
図13には、GUI制御部117により表示装置116の画面に表示される再構成条件設定用のグラフィカルユーザインターフェース(GUI)を示している。このGUIには、スキャノグラム201とともに、再構成条件設定項目が表示される。スキャノグラム201上には、撮影視野が四角図形202により示されている。この四角カーソル202の長さ(紙面上下)は撮影視野の有効幅Wに対応し、また、四角カーソル202の幅(紙面左右)は撮影視野の直径2Rに対応している。また、撮影視野の中心線203がスキャノグラム201上に示されている。
【0041】
スキャノグラム201の下側には、再構成関数、フィルタ、スライス厚、スライスピッチ、撮影視野の径R、撮影視野の有効幅W、撮影視野の中心位置(X、Y)、画像枚数それぞれの設定ボタン204乃至212が配列される。スライスピッチは隣接スライスとの間の中心線間距離を表し、このスライスピッチ、スライス厚、そして撮影視野の有効幅Wにより、画像枚数が自動設定される。撮影視野の径Rに対応するボタン208に数値入力すると、それに追従して四角カーソル202の幅が変化し、逆に四角カーソル202の幅をポインタドラッグにより変化させると、それに追従して撮影視野の径Rに対応するボタン208の数値が変化する。同様に、撮影視野の有効幅Wと四角カーソル202の長さとは関連して変化する。また、撮影視野の中心位置(X、Y)と、四角カーソル202の位置及び中心線203とも、関連して変化する。
【0042】
撮影視野の径Rと撮影視野の有効幅Wそれぞれの入力方法としては、数値入力、四角カーソル202の伸縮の他に、図14に示すように、頭部(HEAD)、肺野(LUNGS)、胴体(BODY)等の検査部位の選択肢をリストしたプルダウンメニュー213,214から選択的に指定する方法が、入力操作を支援するために用意されている。頭部(HEAD)、肺野(LUNGS)、胴体(BODY)等の各検査部位には、撮影視野の径Rと撮影視野の有効幅Wそれぞれの標準値が関連付けられている。GUI制御部117は、選択された検査部位に関連付けられている標準値に撮影視野の径Rと撮影視野の有効幅Wそれぞれの値を自動的に設定する。
【0043】
なお、プルダウンメニュー213,214の選択肢の表記方法は、検査部位の名称に限定されることはなく、図15に示すように、S(small)、M(medium)、L(large)というサイズで表記するようにしても良い。
【0044】
(変形例)
本発明は、上述した実施形態に限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で種々変形して実施することが可能である。さらに、上記実施形態には種々の段階が含まれており、開示される複数の構成要件における適宜な組み合わせにより種々の発明が抽出され得る。例えば、実施形態に示される全構成要件から幾つかの構成要件が削除されてもよい。
【0045】
【発明の効果】
本発明によれば、コーンビーム型X線コンピュータ断層撮影装置において、撮影視野の径に対する有効長の依存性を緩和することができる。
【図面の簡単な説明】
【図1】従来において、径がRLLに設定された撮影視野の側面図。
【図2】従来において、径がRSSに設定された撮影視野の側面図。
【図3】従来において図1の撮影視野の有効長WLLと、図2の撮影視野の有効長WSSとを示す図。
【図4】本発明の実施例によるX線コンピュータ断層撮影装置の構成図。
【図5】図4の2次元アレイタイプのX線検出器の斜視図。
【図6】図4の拡張領域決定部により撮影視野の径RLLに応じて決定される拡張領域の拡張幅EDLLを示す図。
【図7】図4の拡張領域決定部により撮影視野の径RSSに応じて決定される拡張領域の拡張幅EDSSを示す図。
【図8】図6の拡張領域により拡大された撮影視野の有効長WLLと、図7の拡張領域により拡大された撮影視野の有効長WSSとを示す図。
【図9】図6の拡張領域上の仮想レイを示す図。
【図10】図7の拡張領域上の仮想レイを示す図。
【図11】図4のデータストレージ装置に記憶される撮影視野の径に対する拡張幅(検出素子列の数)の対応表を示す図。
【図12】図4のデータストレージ装置に記憶される撮影視野の径に対するヘリカルピッチの対応表を示す図。
【図13】図4のGUI制御部により提供される再構成条件設定のためのグラフィカルユーザインタフェースを示す図。
【図14】図4のGUI制御部により提供される再構成条件設定のためのグラフィカルユーザインタフェースの他の例を示す図。
【図15】図4のGUI制御部により提供される再構成条件設定のためのグラフィカルユーザインタフェースのさらに他の例を示す図。
【符号の説明】
100…ガントリ、
102…回転リング、
101…X線管装置、
109…高電圧発生装置、
103…2次元アレイ型X線検出器、
108…検出素子、
104…データ収集システム、
106…前処理装置、
111…データストレージ装置、
117…GUI制御部、
116…表示装置、
115…ポインティング・デバイス(入力装置)、
112…拡張領域長決定部、
113…投影データ拡張部、
114…再構成装置。[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a so-called cone beam X-ray computed tomography apparatus that scans a subject with a pyramid-shaped X-ray beam in order to obtain stereoscopic information.
[0002]
[Prior art]
In the cone beam scanning method, a subject is scanned with X-rays emitted from an X-ray tube and trimmed into a pyramid shape by an X-ray collimator. X-rays transmitted through the subject are detected by a two-dimensional array type detector. As this type of X-ray detector, a relatively small number of line detectors, typically four rows are arranged. Recently, however, the number of columns exceeds 32 by adopting a solid-state detection element that combines a scintillator element and a photodiode element, or a solid-state detection element such as selenium that directly converts X-rays into electric charges. X-ray detectors have appeared. The X-ray detector has various external shapes such as a cylindrical shape and a planar shape.
[0003]
A Feldkamp method is generally used as the cone beam image reconstruction method. The Feldkamp method is an approximate reconstruction method based on the fan beam convolution and back projection method. The convolution process assumes that the cone angle is relatively small and regards the data as fan projection data. Done. However, the back projection process is performed along an actual ray.
[0004]
That means
(1) The projection data is weighted depending on the Z coordinate,
(2) Convolution the same reconstruction function as the fan beam reconstruction into the data of (1),
(3) Back-project the data of (2) along an oblique actual ray with a cone angle,
The image is reconstructed by the procedure.
[0005]
Such an image reconstruction method has a problem that the effective length of the field of view is changed depending on the diameter of the field of view. This will be specifically described below.
[0006]
FIG. 1 shows a field of view with a diameter set to a relatively long R LL from the side, and FIG. 2 shows a field of view with a diameter set to a relatively short R SS . The diameter R of the field of view is set to a length that can accommodate the examination site such as the head, lungs, and trunk. The effective length of the field of view is defined as the “effective length” as the length of the field of view in which the set radius R is maintained (the length of the field of view regarding the body axis direction of the subject). .
[0007]
Since the Feldkamp method requires data for one rotation, the maximum range in which an image can be reconstructed is limited to a cylindrical shape. The effective length of the field of view in which the diameter R is maintained within this range is limited to W LL when the diameter is R LL , as shown in FIG. 3A, and when the diameter is R SS , As shown in FIG. 3B, it is limited to W SS . As described above, the effective length of the photographing field of view changes in accordance with the set diameter of the photographing field of view.
[0008]
[Problems to be solved by the invention]
An object of the present invention is to alleviate the dependence of the effective length on the diameter of the field of view in a cone beam X-ray computed tomography apparatus.
[0009]
[Means for Solving the Problems]
An X-ray computed tomography apparatus according to an aspect of the present invention includes an X-ray tube apparatus configured to irradiate a subject with X-rays in a pyramid shape, and a slice direction for detecting X-rays transmitted through the subject. A detector having a plurality of detection elements arrayed in at least a slice direction with respect to two directions with respect to the channel direction, and the slice direction from the region where the detection elements are arrayed based on actual data detected by the detection elements A data extension unit configured to generate virtual data corresponding to an extension region located outside, a reconstruction unit configured to reconstruct image data based on the real data and the virtual data An apparatus, an input device for inputting a diameter of the photographing field of view, and a length of the extended region are determined based on the input diameter of the photographing field of view. ; And a clad region length determining unit.
[0010]
DETAILED DESCRIPTION OF THE INVENTION
DESCRIPTION OF THE PREFERRED EMBODIMENTS Preferred embodiments of the X-ray computed tomography apparatus according to the present invention will be described below with reference to the drawings. In the X-ray computed tomography system, ROTATE / ROTATE-TYPE, which rotates around the subject as a single X-ray tube and X-ray detector, and a large number of detection elements arrayed in a ring shape are fixed. There are various types such as STATIONARY / ROTATE-TYPE in which only the X-ray tube rotates around the subject, and the present invention can be applied to any type. Here, it is described as ROTATE / ROTATE-TYPE.
[0011]
Further, in order to reconstruct image data (tomographic image data), projection data for one rotation around the subject and 360 ° is required, and projection data for 180 ° + fan angle is also required in the half scan method. . The present invention can be applied to any reconstruction method. Here, the 360 ° method will be described as an example.
[0012]
In addition, the mechanism for converting incident X-rays into electric charges is based on an indirect conversion type in which X-rays are converted into light by a phosphor such as a scintillator and the light is further converted into electric charges by a photoelectric conversion element such as a photodiode. The generation of electron-hole pairs in a semiconductor such as selenium and the transfer to the electrodes, that is, the direct conversion type utilizing a photoconductive phenomenon, are the mainstream. Any of these methods may be adopted as the detection element.
[0013]
In recent years, the so-called multi-tube type X-ray CT apparatus in which a plurality of pairs of X-ray tubes and X-ray detectors are mounted on a rotating ring has been commercialized, and the development of peripheral technologies has been advanced. The present invention can be applied to both a conventional single-tube X-ray CT apparatus and a multi-tube X-ray CT apparatus. Here, a single tube type will be described.
[0014]
FIG. 4 is a diagram showing the configuration of the X-ray computed tomography apparatus according to the present embodiment. The
[0015]
A two-dimensional array
[0016]
When imaging, a subject is placed between the
[0017]
A
[0018]
Data (pure raw data) output from the
[0019]
The
[0020]
In order to fix the effective length of the imaging field in which the diameter R of the imaging field set via the
[0021]
The length of the extended area is calculated based on the diameter R of the imaging field. Alternatively, a table in which a plurality of lengths of the extended area are associated with a plurality of diameters of the imaging field of view is created in advance, stored in the data storage device 111, and set from the table via the
[0022]
Although the details of the length of the extended area will be described later, the effective length W of the imaging field of view in which the diameter R set by the operator is maintained does not depend on the size of the diameter R of the imaging field of view set arbitrarily. It is determined for the purpose of fixing to a constant or substantially constant.
[0023]
The projection
[0024]
The
[0025]
In the backprojection process, strictly speaking, the
[0026]
The reconstructing
[0027]
FIG. 6 shows the length (expansion width) ED LL of the extension region determined by the extension region
[0028]
The purpose of creating the virtual data is to alleviate the dependency of the effective length W on which the diameter of the photographing field of view R is maintained. In other words, even if the photographing field diameter R is variously set, the effective length W of the photographing field in which the diameter R is maintained does not vary or the variation is suppressed to a minute range. In order to achieve this object, an extension region is added outside the actual sensitive region in the slice direction of the
[0029]
As shown in FIGS. 6 and 8A, when the diameter R of the field of view is set to a relatively long diameter RLL , the length of the extended region is determined to be EDLL . The effective length of the field of view is determined as W LL by the extended width ED LL . Of course, the effective length W LL of this field of view is longer than the effective length WR LL determined depending on the actual sensitive area of the
[0030]
As shown in FIGS. 7 and 8B, when the diameter R of the field of view is set to a relatively short diameter R SS , the length of the extended region is determined to be ED SS . The effective length of the field of view is determined as W SS by this extended width ED SS .
[0031]
The extension area length ED LL is determined to be longer than the extension area length ED SS . Each length is equal to or substantially equal to the effective length W of the final imaging field, in other words, the effective length W of the final imaging field is fixed to a predetermined length or almost a predetermined length. Based on the shortest distance from the X-ray focal point of the
[0032]
Since the extended area is set outside the actual sensitive area of the
[0033]
As shown in FIGS. 9 and 10, one of the methods is the same as the virtual data virtual detection element in the slice direction, and is detected by a real detection element closest to the virtual data virtual detection element. This is a method in which data, that is, actual data detected by a detection element located at the outermost position in the slice direction among actual detection elements is directly applied to virtual data.
[0034]
In practice, this can be realized by data read control from the data storage device 111 holding actual data to the
[0035]
In this method, virtual data of a plurality of virtual detection elements arranged in the slice direction on the extended region is unified with the same actual data. The virtual data actually affects the peripheral portion of the slice indicated by the mesh line in FIGS. In actual examination, the region of interest is often placed near the center of the slice, and the peripheral portion of the slice is less important than the vicinity of the center. In addition, when the region of interest is arranged near the center of the slice, the peripheral portion of the slice often has a relatively simple tissue structure. Therefore, the deterioration in image quality in the peripheral part is not so much of a problem. As a modification, instead of the data of the detection elements in the outermost row, the actual data of the detection elements in one row from the outermost row may be assigned to the virtual data.
[0036]
The other method is a method of creating virtual data by extrapolation from real data detected by a plurality of real detection elements near the virtual detection element. Specifically, for example, as shown in FIG. 9 and FIG. 10, the actual data d1 detected by the actual detection element located closest to the virtual detection element of the virtual data, that is, the outermost actual detection element, and the outermost data Calculation is performed based on the distances S1 and S2 on the detection surface between the actual detection element and the virtual detection element from the actual data d2 of the actual detection element one inner side than the actual detection element. Of course, interpolation may be performed using actual data of three or more existing prefecture elements in the vicinity including the outermost part. Although this method is somewhat improved over the previous method, the image quality (which means image reproducibility and reliability in this case) is still somewhat reduced, but the influence of this method appears strongly in the peripheral portion of the slice, It is assumed that the degradation of the image quality in the peripheral portion is not so much a problem in practice. As this modified sequence, an average value of real data detected by a plurality of real detection elements close to the virtual detection element may be created as virtual data.
[0037]
In this way, an extended area is set outside the actual sensitive area with respect to the slice direction of the
[0038]
FIG. 11 shows the result of simulating the relationship between the imaging field diameter R and the extended region length ED. In FIG. 11, the length ED of the extension region is represented as the number of virtual detection element rows to be extended. When the
[0039]
Further, in the helical scan, it is necessary to change the unit distance that the couch top moves, that is, the helical pitch, in order to suppress the problem that the effective length of the field of view varies depending on the diameter R of the field of view. It was. As shown by a thick line in FIG. 12, the fluctuation of the helical pitch can be suppressed by using this method.
[0040]
FIG. 13 shows a graphical user interface (GUI) for setting reconstruction conditions displayed on the screen of the
[0041]
Below the
[0042]
As input methods for the imaging field diameter R and the imaging field effective width W, in addition to numerical input and expansion / contraction of the
[0043]
It should be noted that the method of representing the choices in the pull-down
[0044]
(Modification)
The present invention is not limited to the above-described embodiments, and various modifications can be made without departing from the scope of the invention at the stage of implementation. Furthermore, the above embodiment includes various stages, and various inventions can be extracted by appropriately combining a plurality of disclosed constituent elements. For example, some constituent requirements may be deleted from all the constituent requirements shown in the embodiment.
[0045]
【The invention's effect】
According to the present invention, in the cone beam X-ray computed tomography apparatus, the dependency of the effective length on the diameter of the field of view can be relaxed.
[Brief description of the drawings]
FIG. 1 is a side view of an imaging field of view with a diameter set to R LL in the related art.
FIG. 2 is a side view of an imaging field of view in which the diameter is set to R SS in the related art.
3 is a diagram showing the effective length W LL of the imaging field of view of FIG. 1 and the effective length W SS of the imaging field of view of FIG. 2 in the prior art.
FIG. 4 is a configuration diagram of an X-ray computed tomography apparatus according to an embodiment of the present invention.
5 is a perspective view of the two-dimensional array type X-ray detector of FIG. 4;
6 is a diagram showing an expansion width ED LL of an expansion region determined according to a shooting field diameter R LL by the expansion region determination unit of FIG. 4;
7 is a diagram showing an expansion width ED SS of an expansion area determined according to the diameter R SS of the field of view by the expansion area determination unit in FIG. 4;
8 is a diagram showing an effective length W LL of the photographing field enlarged by the extended region of FIG. 6 and an effective length W SS of the photographing field enlarged by the extended region of FIG. 7;
9 is a diagram showing a virtual ray on the extended area of FIG. 6. FIG.
10 is a view showing a virtual ray on the extended area of FIG. 7;
11 is a view showing a correspondence table of the expansion width (number of detection element arrays) with respect to the diameter of the field of view to be stored in the data storage device of FIG. 4;
12 is a view showing a correspondence table of helical pitches with respect to the diameter of the field of view to be stored stored in the data storage device of FIG. 4;
13 is a diagram showing a graphical user interface for setting reconfiguration conditions provided by the GUI control unit of FIG. 4;
14 is a diagram showing another example of a graphical user interface for setting reconstruction conditions provided by the GUI control unit of FIG. 4;
FIG. 15 is a diagram showing still another example of a graphical user interface for setting reconstruction conditions provided by the GUI control unit of FIG. 4;
[Explanation of symbols]
100 ... Gantry,
102 ... rotating ring,
101 ... X-ray tube device,
109 ... high voltage generator,
103 ... Two-dimensional array type X-ray detector,
108 ... detecting element,
104 ... Data collection system,
106: Pretreatment device,
111 ... data storage device,
117 ... GUI control unit,
116 ... display device,
115 ... pointing device (input device),
112 ... Extended region length determination unit,
113 ... Projection data expansion unit,
114: Reconstruction device.
Claims (10)
前記被検体を透過したX線を検出するスライス方向とチャンネル方向との2方向に関して少なくともスライス方向にアレイされた複数の検出素子を有する検出器と、
前記検出素子により検出された実データに基づいて、前記検出素子がアレイされた領域から前記スライス方向に関して外側に位置する拡張領域に対応する仮想データを発生するように構成されたデータ拡張部と、
前記実データと、前記仮想データとに基づいて、画像データを再構成するように構成された再構成装置と、
撮影視野の径を入力する入力装置と、
前記入力された撮影視野の径に基づいて、前記拡張領域の長さを決定するように構成された拡張領域長決定部とを具備することを特徴とするX線コンピュータ断層撮影装置。An X-ray tube device configured to irradiate a subject with X-rays in a pyramid shape;
A detector having a plurality of detection elements arrayed in at least a slice direction with respect to two directions of a slice direction and a channel direction for detecting X-rays transmitted through the subject;
A data extension unit configured to generate virtual data corresponding to an extension area located outside the area in which the detection elements are arrayed with respect to the slice direction, based on actual data detected by the detection elements;
A reconstruction device configured to reconstruct image data based on the real data and the virtual data;
An input device for inputting the diameter of the field of view;
An X-ray computed tomography apparatus comprising: an extended region length determination unit configured to determine a length of the extended region based on the input imaging visual field diameter.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002087047A JP4280018B2 (en) | 2001-04-03 | 2002-03-26 | X-ray computed tomography system |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001-104915 | 2001-04-03 | ||
JP2001104915 | 2001-04-03 | ||
JP2002087047A JP4280018B2 (en) | 2001-04-03 | 2002-03-26 | X-ray computed tomography system |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2009016849A Division JP2009090139A (en) | 2001-04-03 | 2009-01-28 | X-ray computerized tomographic apparatus |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2002360562A JP2002360562A (en) | 2002-12-17 |
JP4280018B2 true JP4280018B2 (en) | 2009-06-17 |
Family
ID=26613028
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002087047A Expired - Lifetime JP4280018B2 (en) | 2001-04-03 | 2002-03-26 | X-ray computed tomography system |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4280018B2 (en) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4642331B2 (en) * | 2003-08-21 | 2011-03-02 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X-ray CT apparatus and control method thereof |
US7570733B2 (en) * | 2005-06-10 | 2009-08-04 | General Electric Company | Step-and-shoot cardiac CT imaging |
US8340241B2 (en) * | 2006-02-27 | 2012-12-25 | Kabushiki Kaisha Toshiba | Image display apparatus and X-ray computed tomography apparatus |
JP5468190B2 (en) * | 2006-02-27 | 2014-04-09 | 株式会社東芝 | Image display apparatus and X-ray CT apparatus |
JP5196751B2 (en) * | 2006-09-13 | 2013-05-15 | 株式会社東芝 | Computer-aided diagnosis device |
JP2008237886A (en) | 2007-02-28 | 2008-10-09 | Toshiba Corp | X-ray ct device and method of controlling the same |
JP2009006133A (en) | 2007-05-31 | 2009-01-15 | Toshiba Corp | X-ray ct apparatus and method of controlling the same |
JP5390132B2 (en) * | 2007-09-03 | 2014-01-15 | 株式会社東芝 | X-ray CT system |
JP5404767B2 (en) * | 2009-03-04 | 2014-02-05 | 株式会社日立メディコ | X-ray CT apparatus and image correction program |
US8712134B2 (en) * | 2011-10-18 | 2014-04-29 | Kabushiki Kaisha Toshiba | Method and system for expanding axial coverage in iterative reconstruction in computer tomography (CT) |
JP6026132B2 (en) | 2012-04-11 | 2016-11-16 | 東芝メディカルシステムズ株式会社 | X-ray computed tomography apparatus and reconstruction processing program |
-
2002
- 2002-03-26 JP JP2002087047A patent/JP4280018B2/en not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
JP2002360562A (en) | 2002-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2009090139A (en) | X-ray computerized tomographic apparatus | |
US6922457B2 (en) | Computer tomography apparatus | |
JP4282302B2 (en) | X-ray CT system | |
US9271681B2 (en) | X-ray computed tomography apparatus | |
JP4612294B2 (en) | X-ray computed tomography apparatus, x-ray computed tomography control method, and x-ray computed tomography program | |
US8290233B2 (en) | X-ray computed tomography apparatus and image processing method | |
JP4828839B2 (en) | X-ray computed tomography apparatus, image processing apparatus and image processing method | |
JP2004180715A (en) | X-ray computed tomography apparatus | |
JP2005305026A (en) | X-ray computer tomographical apparatus | |
US7529335B2 (en) | Voxel-driven spiral reconstruction for cone-beam computer tomography | |
JP4519434B2 (en) | Super-resolution processing apparatus and medical image diagnostic apparatus | |
JP4280018B2 (en) | X-ray computed tomography system | |
JP4469555B2 (en) | X-ray computed tomography system | |
JPH0661330B2 (en) | Method and apparatus for creating tomographic image | |
US8705688B2 (en) | Medical image processing apparatus, X-ray computed tomography apparatus, and medical image processing method | |
JP5605985B2 (en) | X-ray computed tomography system | |
JP2825352B2 (en) | CT device | |
JP2005095329A (en) | Superresolving apparatus and medical diagnostic imaging apparatus | |
JP4551612B2 (en) | Computed tomography equipment | |
US9784694B2 (en) | X-ray computed tomography apparatus and reconstruction processing method | |
JP2010279427A (en) | X-ray computed tomographic system | |
JP4607364B2 (en) | X-ray CT system | |
JP2825471B2 (en) | CT device | |
JP2010042096A (en) | X-ray computed tomography apparatus | |
JP2007159982A (en) | X-ray computed tomography apparatus |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050325 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20080709 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080715 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080916 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20081202 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090128 |
|
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: 20090303 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20090313 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120319 Year of fee payment: 3 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 4280018 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120319 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130319 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140319 Year of fee payment: 5 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313114 Free format text: JAPANESE INTERMEDIATE CODE: R313111 Free format text: JAPANESE INTERMEDIATE CODE: R313117 |
|
R371 | Transfer withdrawn |
Free format text: JAPANESE INTERMEDIATE CODE: R371 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313117 Free format text: JAPANESE INTERMEDIATE CODE: R313111 Free format text: JAPANESE INTERMEDIATE CODE: R313114 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
EXPY | Cancellation because of completion of term |