JP4467522B2 - Tomographic image constructing apparatus and method - Google Patents

Tomographic image constructing apparatus and method Download PDF

Info

Publication number
JP4467522B2
JP4467522B2 JP2005512615A JP2005512615A JP4467522B2 JP 4467522 B2 JP4467522 B2 JP 4467522B2 JP 2005512615 A JP2005512615 A JP 2005512615A JP 2005512615 A JP2005512615 A JP 2005512615A JP 4467522 B2 JP4467522 B2 JP 4467522B2
Authority
JP
Japan
Prior art keywords
image
tomographic
images
dimensional
scanning line
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
JP2005512615A
Other languages
Japanese (ja)
Other versions
JPWO2005011499A1 (en
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.)
Hitachi Healthcare Manufacturing Ltd
Original Assignee
Hitachi Medical Corp
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 Hitachi Medical Corp filed Critical Hitachi Medical Corp
Publication of JPWO2005011499A1 publication Critical patent/JPWO2005011499A1/en
Application granted granted Critical
Publication of JP4467522B2 publication Critical patent/JP4467522B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment

Description

本発明は、2枚程度の少ない枚数のレントゲン画像から断層像を構成する断層像構成装置及び方法に関する。   The present invention relates to a tomographic image forming apparatus and method for forming a tomographic image from as few as two X-ray images.

被検体の断層像は、通常X線CT装置などの断層像生成専用装置を用いて撮影する。近年、医療施設に最も普及しているX線撮影装置を使い、低いX線被曝で撮影されたレントゲン画像を例えば2枚用いて、それらの2枚のレントゲン画像から前記被検体の断層像やその断層像から三次元画像を構成したいという要求がある。   The tomographic image of the subject is usually captured using a tomographic image generation apparatus such as an X-ray CT apparatus. In recent years, for example, two X-ray images taken with low X-ray exposure using an X-ray imaging apparatus most widely used in medical facilities, and the tomographic image of the subject and There is a demand to construct a three-dimensional image from a tomogram.

上記2枚のレントゲン画像を用いて画像構成する技術は、“IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 22, NO. 6, JUNE 2003”、710〜721頁、“Kinematic and Deformation Analysis of 4-D Coronary Arterial Trees Reconstructed From Cine Angiograms”に記載された方法(「従来方法」と称する)がある。この従来方法は、直交するレントゲン画像からの投影により心臓に走行する樹状の血管の中心線を算出し、その算出された中心線に沿ってそれぞれのレントゲン画像から投影される血管の太さ情報を付加することで、三次元画像を構成する。   The technique of image composition using the above two X-ray images is “IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 22, NO. 6, JUNE 2003”, 710-721, “Kinematic and Deformation Analysis of 4-D Coronary Arterial. There is a method described in “Trees Reconstructed From Cine Angiograms” (referred to as “conventional method”). This conventional method calculates the center line of a dendritic blood vessel that travels to the heart by projection from orthogonal X-ray images, and information on the thickness of the blood vessels projected from each X-ray image along the calculated center line Is added to construct a three-dimensional image.

しかしながら、上記従来方法では、直交する位置関係の2枚のレントゲン画像を用いた三次元画像の構成に限定されているので、任意の異なる角度に配置された2枚のレントゲン画像を逆投影して断層像を構成することは一切言及されていない。   However, since the conventional method is limited to the configuration of a three-dimensional image using two X-ray images that are orthogonal to each other, two X-ray images arranged at arbitrary different angles are back-projected. There is no mention of constructing a tomographic image.

本発明の断層像構成装置は、レントゲン装置によって撮影された被検体の複数の異なる方向からのレントゲン画像のそれぞれを二値化して複数の二値化画像を作成する二値化画像作成手段と、前記二値化画像作成手段によって作成された前記複数の二値化画像の互いに対応する位置に所定方向に順次移動する走査線を設定する走査線設定手段と、前記走査線設定手段によって設定された走査線が順次移動する位置における前記複数の二値化画像の画素点を逆投影することにより断層像を構成する断層像構成手段と、を備え、前記走査線設定手段は、2次元の走査線方向についてランレングスの出現頻度を演算し、前記演算された出現頻度に基づいて、2次元の走査線方向のうちの一方を設定する走査線方向設定手段を含む。 The tomographic image construction apparatus of the present invention is a binarized image creating means for binarizing each of the X-ray images from a plurality of different directions of the subject imaged by the X-ray device to create a plurality of binarized images, A scanning line setting unit that sets scanning lines that sequentially move in predetermined directions to positions corresponding to each other of the plurality of binarized images created by the binarized image creating unit, and set by the scanning line setting unit A tomographic image forming unit configured to form a tomographic image by back projecting pixel points of the plurality of binarized images at a position where the scanning line sequentially moves, and the scanning line setting unit includes a two-dimensional scanning line Scan line direction setting means for calculating the appearance frequency of the run length for the direction and setting one of the two-dimensional scan line directions based on the calculated appearance frequency is included.

これによって、被検体へのX線被曝の低減に配慮し、より少ない枚数、つまり2枚のレントゲン画像を用いて被検体の断層像を得ることができる。ここでのレントゲン画像は、平面検出器等を用いて撮影されたデジタル画像や、フィルム画像をデジタル化した画像を含む。また、上記より少ない枚数は、3枚以上のレントゲン画像であってもよい。この場合は、上記X線被曝の低減の効果は減少するものの、結果として構成される断層像の精度が向上される。   Thus, considering to reduce X-ray exposure to the subject, a tomographic image of the subject can be obtained using a smaller number of X-ray images, that is, two X-ray images. The X-ray image here includes a digital image taken using a flat detector or the like, or an image obtained by digitizing a film image. Further, the smaller number than the above may be three or more X-ray images. In this case, although the effect of reducing the X-ray exposure is reduced, the accuracy of the resulting tomographic image is improved.

また、本発明の望ましい実施形態によれば、前記走査線設定手段は、前記断層像構成手段によって構成された断層像のそれぞれについて、前記走査線の移動の正方向と逆方向における断層像との相関に基づき偽領域を検出する偽領域検出手段と、前記偽領域検出手段によって検出された偽領域を削除する偽領域削除手段と、を含む。   Also, according to a preferred embodiment of the present invention, the scanning line setting means calculates a tomographic image in each of the tomographic images formed by the tomographic image forming means in a direction opposite to the forward direction of the scanning line. A false area detecting unit for detecting a false area based on the correlation; and a false area deleting unit for deleting the false area detected by the false area detecting unit.

これにより、走査線の移動方向が正方向の断層像と走査線の移動方向が逆方向の断層像との相関に基づいて偽領域を削除するので、より真正な領域だけを用いた正確な断層像が得られる。前記相関は所定の値に満たない領域を検出すれば、その領域が偽領域であるので、効率的に偽領域を削除することができる。   As a result, the false region is deleted based on the correlation between the tomographic image in which the scanning line movement direction is the positive direction and the tomographic image in which the scanning line movement direction is the reverse direction. An image is obtained. If an area where the correlation is less than a predetermined value is detected, the area is a false area, so that the false area can be efficiently deleted.

また、本発明の望ましい実施形態によれば、前記走査線設定手段は、前記走査線上で画素値が所定の値をとる領域が一つである位置、及び操作者により指定された位置のいずれかの位置を走査線の初期位置として設定する。   Also, according to a preferred embodiment of the present invention, the scanning line setting means is either one of a position where a pixel value takes a predetermined value on the scanning line and a position designated by an operator. Is set as the initial position of the scanning line.

これにより、走査線の初期位置は任意に設定できることになるから、抽出したい領域の近傍に設定すれば、より効率的に断層像を構成のための2値化領域を移動させることができるから、結果として断層像の構成のための演算を高速化できる。   As a result, the initial position of the scanning line can be arbitrarily set, so if it is set in the vicinity of the region to be extracted, the binarized region for constructing the tomographic image can be moved more efficiently. As a result, the computation for constructing the tomographic image can be speeded up.

また、本発明の望ましい実施形態によれば、前記断層像構成装置は、前記断層像構成手段によって構成された断層像を隣接して複数構成し、前記複数構成された断層像から所定のレンダリング法により三次元画像を構成する三次元画像構成手段と、前記三次元画像構成手段によって構成された三次元画像を表示する画像表示手段と、を更に備える。   Also, according to a preferred embodiment of the present invention, the tomographic image constructing device comprises a plurality of adjacent tomographic images formed by the tomographic image forming means, and a predetermined rendering method from the plurality of configured tomographic images. 3D image constructing means for constructing a 3D image, and image display means for displaying the 3D image constructed by the 3D image constructing means.

これによれば、前記断層像から三次元画像を構成することで、被検体の形状を容易に把握できるようになる。   According to this, it is possible to easily grasp the shape of the subject by constructing a three-dimensional image from the tomographic image.

本発明の断層像構成方法は、レントゲン装置によって撮影された被検体の複数の異なる方向からのレントゲン画像のそれぞれを二値化して複数の二値化画像を作成する二値化画像作成工程と、前記二値化画像作成工程によって作成された前記複数の二値化画像の互いに対応する位置に所定方向に順次移動する走査線を設定する走査線設定工程と、前記走査線設定工程によって設定された走査線が順次移動する位置における前記複数の二値化画像の画素点を逆投影することにより断層像を構成する断層像構成工程と、を備え、前記走査線設定工程は、2次元の走査線方向についてランレングスの出現頻度を演算し、前記演算された出現頻度に基づいて、2次元の走査線方向のうちの一方を設定する走査線方向設定工程を含む。 The tomographic image constructing method of the present invention includes a binarized image creating step of binarizing each of the X-ray images from a plurality of different directions of the subject imaged by the X-ray device to create a plurality of binarized images, A scanning line setting step for setting scanning lines that sequentially move in predetermined directions to positions corresponding to each other of the plurality of binarized images created by the binarized image creating step, and the scanning line setting step A tomographic image forming step of forming a tomographic image by back projecting pixel points of the plurality of binarized images at positions where the scanning line sequentially moves, and the scanning line setting step includes a two-dimensional scanning line A scan line direction setting step of calculating a run length appearance frequency for the direction and setting one of the two-dimensional scan line directions based on the calculated appearance frequency is included.

これによって、被検体へのX線被曝の低減に配慮し、より少ない枚数、つまり2枚のレントゲン画像を用いて被検体の断層像を得ることができる。   Thus, considering to reduce X-ray exposure to the subject, a tomographic image of the subject can be obtained using a smaller number of X-ray images, that is, two X-ray images.

以下、添付図面に従って、本発明に係るレントゲン画像を用いた断層像構成装置及び方法の望ましい実施の形態について詳説する。   Hereinafter, preferred embodiments of a tomographic image constructing apparatus and method using X-ray images according to the present invention will be described in detail with reference to the accompanying drawings.

図1に、本実施の形態に係る断層像構成装置10の構成を示す。断層像構成装置10は、中央処理装置(以下、CPUという)11、主メモリ12、磁気ディスク13、表示メモリ14、CRT15、コントローラ16、マウス17、キーボード18を備えており、これらは共通バス19を介して電気的に接続されている。   FIG. 1 shows a configuration of a tomographic image forming apparatus 10 according to the present embodiment. The tomographic image forming apparatus 10 includes a central processing unit (hereinafter referred to as CPU) 11, a main memory 12, a magnetic disk 13, a display memory 14, a CRT 15, a controller 16, a mouse 17, and a keyboard 18, and these are a common bus 19. It is electrically connected via.

主メモリ12は画像やデータの一時記憶、処理用の領域として用いられ、処理結果は表示メモリ14を介してCRT15に表示されるとともに磁気ディクス13に格納され、再表示や結果参照に利用される。磁気ディスク13には、レントゲン画像が格納されている。レントゲン画像は、LAN1Aを介して接続されたレントゲン装置1Bから取得することも可能である。レントゲン装置1BはX線透視撮影台、デジタルラジオグラフィ装置などあらゆるレントゲン画像が撮影可能な装置のすべてを含む。なお、上記レントゲン画像は、平面検出器等を用いて撮影されたデジタル画像や、フィルム画像をデジタル化した画像も含む。   The main memory 12 is used as an area for temporary storage and processing of images and data, and processing results are displayed on the CRT 15 via the display memory 14 and stored in the magnetic disk 13 for redisplay and result reference. . An X-ray image is stored on the magnetic disk 13. The X-ray image can also be acquired from the X-ray apparatus 1B connected via the LAN 1A. The X-ray apparatus 1B includes all apparatuses capable of capturing all X-ray images, such as an X-ray fluoroscopic table and a digital radiography apparatus. The X-ray image includes a digital image taken using a flat detector or the like, and an image obtained by digitizing a film image.

(第一の実施形態)
次に、図2のフローチャートに従って断層像構成処理を説明する。ここでは2枚のレントゲン画像から断層像を構成する場合について説明する。これらのレントゲン画像31、32、断層像33、三次元画像34は、図3に示す位置関係となっている。断層像33は、レントゲン画像31,32を二値化し、その二値化されたレントゲン画像31、32の一部を逆投影して構成し、三次元画像34は断層像33からレンダリングアルゴリズムを用いて構成する。
(First embodiment)
Next, tomographic image construction processing will be described with reference to the flowchart of FIG. Here, a case where a tomographic image is constructed from two X-ray images will be described. These X-ray images 31, 32, the tomographic image 33, and the three-dimensional image 34 have the positional relationship shown in FIG. The tomographic image 33 is formed by binarizing the X-ray images 31 and 32 and partially projecting the binarized X-ray images 31 and 32, and the three-dimensional image 34 is generated from the tomographic image 33 using a rendering algorithm. Configure.

(ステップ20)
図4に示すように、X線源41によって撮影されたレントゲン画像31は、目的臓器a、その他臓器b、cを1、その他の領域を0とする二値化処理が行われる。同様に、X線源42によって撮影されたレントゲン画像32は、目的臓器a、その他臓器b、cを1、その他の領域を0とする二値化処理が行われる。これらの2枚のレントゲン画像が撮影されたX線源の位置にそれぞれ投影中心51、投影中心52を配置する。それら投影中心からの投影線に囲まれる領域をメモリに記憶すると、図5のようになる。つまり、目的臓器aとその目的臓器以外の領域b、c(「偽領域」と称する)とはその区別がつかない状態でメモリに記憶される。ところが、図6のように、目的臓器が一つだけ撮影された各レントゲン画像からは、投影中心61、投影中心62に偽領域は現れない。つまり、目的臓器aだけがメモリに記憶される。
(Step 20)
As shown in FIG. 4, the X-ray image 31 captured by the X-ray source 41 is subjected to binarization processing in which the target organ a, the other organs b and c are 1, and the other regions are 0. Similarly, the X-ray image 32 taken by the X-ray source 42 is subjected to binarization processing in which the target organ a, the other organs b and c are set to 1, and the other regions are set to 0. A projection center 51 and a projection center 52 are arranged at the position of the X-ray source where these two X-ray images are taken, respectively. When the area surrounded by the projection line from the projection center is stored in the memory, it is as shown in FIG. That is, the target organ a and the regions b and c (referred to as “false regions”) other than the target organ are stored in the memory in a state where they cannot be distinguished from each other. However, as shown in FIG. 6, no false region appears at the projection center 61 and the projection center 62 from each X-ray image obtained by photographing only one target organ. That is, only the target organ a is stored in the memory.

そこで、断層像構成装置10では、この性質を利用して、目的臓器が一個の所を最初に探すために、例えば図7(a)の番号1の直線を最初として二値化処理後、逆投影処理し、順次、隣の2,3,・・・又は2',3',・・・へ移行して二値化処理及び逆投影処理する。二回目以降の逆投影時には、隣接する領域との連続性により選別する。例えば、図7(a)の領域71や72のように、途中で消失するような血管にあっては抽出しないものとする。   Therefore, in the tomographic image constructing apparatus 10, in order to search for a target organ first using this property, for example, the straight line of number 1 in FIG. Projection processing is performed, and the processing proceeds to adjacent 2,3,... Or 2 ′, 3 ′,. In the second and subsequent back projections, selection is performed based on continuity with adjacent areas. For example, it is assumed that a blood vessel that disappears on the way like the region 71 or 72 in FIG. 7A is not extracted.

ここでの二値化処理は、一般にしきい値処理によるが手動でトレースして目的臓器の内部を1、外部を0としてもよい。また、ここでは画像のX方向の一線値上の1領域を探しているが、図8のように、ランレングス(連続する画素数)が最大の領域81を見つけて、その最大領域81以外の領域に0値を代入して、一線値上の1領域を作り、上記一線値上の1領域を探す方法に代えてもよい。   The binarization processing here is generally based on threshold processing, but may be manually traced to set the inside of the target organ to 1 and the outside to 0. In this example, one area on the one-line value in the X direction of the image is searched. As shown in FIG. 8, an area 81 having the maximum run length (the number of consecutive pixels) is found, and other areas than the maximum area 81 are detected. A method of substituting a 0 value into a region to create one region on one line value and searching for one region on the one line value may be used.

例えば、図9の横方向の点線上には1領域が91、92、93の3個所あるから目的臓器を探すための開始位置にはならないが、図7(b)の点線箇所は1領域が一個所で、開始位置となり得る。逆投影開始位置の設定は、図10に示されるように、ROI101で囲んだ領域から開始してもよい。   For example, since there are three areas 91, 92, and 93 on the horizontal dotted line in FIG. 9, it does not become a starting position for searching for the target organ, but the dotted line in FIG. 7B has one area. It can be the starting position at one place. The setting of the backprojection start position may be started from an area surrounded by the ROI 101 as shown in FIG.

また、逆投影の開始位置は、二値化処理された2枚のレントゲン画像(「二値化画像」と称する)同士で対応する位置に設定する。   Also, the back projection start position is set to a position corresponding to the two binarized X-ray images (referred to as “binarized images”).

(ステップ21)
二値化画像は、図6のように逆投影処理が行われる。ここで、逆投影した領域の形状は、血管を抽出することとしているので、多くはその断面が円形であることから、図11に示すように再構成領域は円形で近似する。また、図12に示すように、逆投影の方向が直交する場合は楕円で近似することとする。このように目的臓器の解剖学的な知見から、任意の形状に設定できる。
(Step 21)
The binarized image is subjected to back projection processing as shown in FIG. Here, since the shape of the back-projected region is to extract a blood vessel, in many cases, the cross-section is circular, so that the reconstruction region is approximated by a circle as shown in FIG. Further, as shown in FIG. 12, when the backprojection directions are orthogonal, the approximation is made by an ellipse. Thus, it can be set to an arbitrary shape from the anatomical knowledge of the target organ.

例えば、前記任意形状の設定は、画面に円、楕円などの任意の形状とメニュー情報を表示し、その表示されたメニューを目的臓器に応じてマウスで選択して入力できるようになっている。   For example, the setting of the arbitrary shape is such that an arbitrary shape such as a circle or an ellipse and menu information are displayed on the screen, and the displayed menu can be selected and input with the mouse according to the target organ.

(ステップ22)
逆投影画像は上記逆投影された結果から再構成され、その再構成された断層像は磁気ディスク等に保存する(画像a)。
(Step 22)
The backprojected image is reconstructed from the result of the backprojection, and the reconstructed tomographic image is stored on a magnetic disk or the like (image a).

(ステップ23)
断層像の再構成ラインは、図7(a)に示すように、その再構成位置を、走査線設定方向(図7(a)、図7(b)では上向き、図7(c)では下向き)に、n(n=1,2等)更新される。
(Step 23)
As shown in FIG. 7 (a), the reconstruction line of the tomographic image has its reconstruction position facing upward in the scanning line setting direction (FIGS. 7 (a) and 7 (b) upward, and downward in FIG. 7 (c)). ) Is updated by n (n = 1, 2, etc.).

(ステップ24)
断層像は、上記更新ライン毎に二値化、逆投影を繰り返し、それぞれの更新ライン毎に得られる(画像b)。
(Step 24)
A tomogram is obtained for each update line by repeating binarization and backprojection for each update line (image b).

(ステップ25)
ステップ22で保存した逆投影画像(画像a)とステップ24で逆投影した断層像(画像b)との類似性を相関演算で求め、両者の断層像の略同一アドレスが画素“1”であれば連続性があるとし、その連続性がない領域は偽領域として画像(b)から削除する。
(Step 25)
Similarity between the backprojected image (image a) stored in step 22 and the tomogram (image b) backprojected in step 24 is obtained by correlation calculation, and the substantially same address of both tomographic images is the pixel “1”. If there is continuity, an area having no continuity is deleted from the image (b) as a false area.

また、連続性があるとは、画像間の引き算をして引き算の総和が所定の値より小であることを意味する。その相関値として、二値化画像のときは、双方の画像領域における相対応する画素の画素値がどちらも画素“1”である割合をいうこともある。   Further, “continuity” means that the subtraction between the images is performed and the sum of the subtraction is smaller than a predetermined value. As the correlation value, in the case of a binarized image, it may be the ratio that the pixel values of the corresponding pixels in both image regions are both pixels “1”.

なお、本発明において、「相関」、「連続性」には、実データ同士の直接的な相関だけでなく、外挿補間法により求めたデータ(図示省略)との相関も含まれる。   In the present invention, “correlation” and “continuity” include not only direct correlation between actual data but also correlation with data (not shown) obtained by extrapolation.

このようにして、断層像構成装置10では、偽領域が削除された正確な断層像を得ることができる。   In this way, the tomographic image construction apparatus 10 can obtain an accurate tomographic image from which the false region is deleted.

(ステップ26)
偽領域を削除した断層像は保存される。
(Step 26)
The tomographic image from which the false region is deleted is stored.

(ステップ27)
図13に示されるコーンビームX線源131、132による2枚のレントゲン画像又は図14に示されるようなファンビームX線源141、142を平行移動して得られるレントゲン画像について、三次元画像再構成に必要な数の断層像が得られるまで、上述の処理(断層像の再構成と偽領域の削除)を繰り返す。
(Step 27)
For the two X-ray images by the cone beam X-ray sources 131 and 132 shown in FIG. 13 or the X-ray image obtained by translating the fan beam X-ray sources 141 and 142 as shown in FIG. The above processing (reconstruction of tomographic images and deletion of false regions) is repeated until the number of tomographic images necessary for the configuration is obtained.

なお、上述の説明ではX方向に走査、一線値上1領域を探しているが、Y方向走査の方が良好に抽出できる場合には次のようにする。まず、図15のように、図10上のX方向、Y方向についてのランレングス(連続する画素数)の分布を求める。X方向走査の方が分布のピーク値が低ければ、X方向走査が正しいことになる。違う場合はY方向を選択することになるが、間違ってX方向又はY方向が選択されたときは、その走査線方向が誤りである旨のエラーメッセージを画面に表示する(図示省略)。   In the above description, scanning in the X direction and searching for one region above the line value are performed. However, when the Y direction scanning can be extracted more favorably, the following is performed. First, as shown in FIG. 15, the distribution of run lengths (number of continuous pixels) in the X direction and Y direction in FIG. 10 is obtained. If the peak value of the distribution is lower in the X direction scanning, the X direction scanning is correct. If they are different, the Y direction is selected. If the X direction or the Y direction is selected by mistake, an error message indicating that the scanning line direction is incorrect is displayed on the screen (not shown).

なお、ランレングスの分布にピークが複数ある場合は、被検体の形状や大きさ等を考慮して走査線方向を選択する。   When there are a plurality of peaks in the run length distribution, the scanning line direction is selected in consideration of the shape and size of the subject.

また、図16のフローチャートに示すように、前記ランレングスの分布から走査線方向を自動的に決定するようにしてもよい。   Also, as shown in the flowchart of FIG. 16, the scanning line direction may be automatically determined from the run length distribution.

(ステップ161)
図15に示されるように、X方向とY方向のランレングスの関係を求める。そして、求められたX方向ランレングスの平均値はY方向ランレングスの平均値より小か否かを比較判定する。この比較判定の対象となるX方向ランレングスの平均値はその平均値の所定の定数倍であってもよい。その判定が小であればステップ163に進み、前記判定が小でなければステップ162に進む。
(Step 161)
As shown in FIG. 15, the relationship between the run lengths in the X direction and the Y direction is obtained. Then, it is determined whether or not the obtained average value of the X-direction run length is smaller than the average value of the Y-direction run length. The average value of the X-direction run lengths to be compared and determined may be a predetermined constant multiple of the average value. If the determination is small, the process proceeds to step 163, and if the determination is not small, the process proceeds to step 162.

(ステップ162)
走査方向をY方向に設定して、二値化及び逆投影処理を行い、ステップ164に進む。
(Step 162)
The scanning direction is set to the Y direction, binarization and back projection processing are performed, and the process proceeds to step 164.

(ステップ163)
走査方向をX方向に設定して、二値化及び逆投影処理を行う。
(Step 163)
The scanning direction is set to the X direction, and binarization and back projection processing are performed.

(ステップ164)
ステップ25で説明した連続性により、1領域を選別する。
(Step 164)
One region is selected based on the continuity described in step 25.

(ステップ28)
ステップ27で得られたレントゲン画像についてステップ22乃至ステップ26の繰り返しが終了したか否かを判断する。
(Step 28)
It is determined whether or not the repetition of Step 22 to Step 26 is completed for the X-ray image obtained in Step 27.

(ステップ29)
単純な管状の血管などでは上述のような処理を1回行うことで本物の領域を抽出できるが、図7(b)に示すように、一方向(図10の場合は、下から上へ)の連続性だけでは本物の領域を全て抽出できない場合がある。このような場合は、反対方向(図7(a)、7(b)及び7(c)の場合は、上から下へ)の連続性を用いることで抽出できる。そこで、走査位置が処理範囲の端に達したら、走査線設定方向を反転して逆投影と偽領域の削除を繰り返す。
(Step 29)
In a simple tubular blood vessel or the like, a real region can be extracted by performing the above-described process once. However, as shown in FIG. 7B, one direction (in the case of FIG. 10, from bottom to top) In some cases, it is not possible to extract all the real regions only by the continuity of. In such a case, it can be extracted by using the continuity in the opposite direction (in the case of FIGS. 7 (a), 7 (b) and 7 (c), from top to bottom). Therefore, when the scanning position reaches the end of the processing range, the scanning line setting direction is reversed, and back projection and false region deletion are repeated.

このような処理の繰り返しは所定の回数行うようにしてもよいし、再構成した断層像を表示して操作者に指示入力させるようにしてもよい。   Such a process may be repeated a predetermined number of times, or a reconstructed tomographic image may be displayed and an operator may input an instruction.

上述のように、走査線設定方向を反転して処理を繰り返すことで、枝分かれした複雑な血管などでも細部まで抽出することができる。例えば、図7(b)に示すように1回の処理(この場合下から上方向)では抽出できない部分でも、図7(c)に示すように処理を繰り返すことで抽出することができる。   As described above, by reversing the scanning line setting direction and repeating the process, it is possible to extract even the complicated blood vessels that are branched. For example, even a portion that cannot be extracted by one process (in this case, from the bottom to the top) as shown in FIG. 7B can be extracted by repeating the process as shown in FIG. 7C.

(ステップ2A)
上述のように、断層像構成装置10では、レントゲン画像から断層像を得ることができるが、断層像から三次元画像を構成することで、被検体の形状を容易に把握できるようになる。三次元画像の構成は、サーフェイス法やボリュームレンダリング法などで行うことができる。三次元画像は、例えば、被検体の隣接する断層像を当該被検体の体軸方向に複数枚得て、一つの目的臓器を抽出し、該抽出された目的臓器を仮想的にメモリ上に前記被検体の体軸方向に積み上げて、その積み上げられた目的臓器を投影面に投影して構成される。
(Step 2A)
As described above, the tomographic image constructing apparatus 10 can obtain a tomographic image from an X-ray image. However, by forming a three-dimensional image from a tomographic image, the shape of the subject can be easily grasped. A three-dimensional image can be constructed by a surface method or a volume rendering method. The three-dimensional image is obtained by, for example, obtaining a plurality of adjacent tomographic images of the subject in the body axis direction of the subject, extracting one target organ, and virtually extracting the extracted target organ on the memory. It is configured by stacking the subject in the body axis direction and projecting the stacked target organs onto the projection plane.

また、ステップ21で説明したように、再構成領域は円形や楕円形で近似しているため、図17に示されるように、レントゲン画像171とレントゲン画像172とから三次元画像173を構成したとき、その分岐部分の大きさが不連続的に変化し現れる。そこで、断層像構成装置10では、図18(a)乃至18(f)に示すような結合処理を行う。   Further, as described in step 21, since the reconstruction area is approximated by a circle or an ellipse, when a three-dimensional image 173 is formed from the X-ray image 171 and the X-ray image 172 as shown in FIG. The size of the branching portion appears discontinuously changing. Therefore, the tomographic image forming apparatus 10 performs a combining process as shown in FIGS. 18 (a) to 18 (f).

まず、図18(a)に示されるように血管二本の領域181、182と血管一本の領域183に分け、図18(b)に示されるように領域181、182、183を分離する。そして、領域181、182、183の接続部分の断面像を図18(c)と図18(d)に示す。例えば、血管一本(図18(a)の183)が血管二本(図18(a)の181、182)に分岐する場合、血管一本だけが写っている断層像(d)(図18(d))を修正する。この場合、図18(e)に示す断層像(e)において、領域181の同心円184と断層像(d)の1領域の論理積の領域、及び断層像182の同心円185と断層像(d)の論理積の領域を残し、他は修正する。修正した結果が図18(f)である。これによって、領域181、182、183が滑らかに結合する。なお、同心円の設定方法の例を図19に示す。領域181と領域182間の辺縁部分の最短距離dRは(1)式によって求められる:
dR=RO−rc1−rc2 (1)
ここで、ROは領域181と領域182の各中心間の距離、rc1は領域181の半径、rc2は領域182の半径である。
First, as shown in FIG. 18A, the region is divided into two blood vessel regions 181 and 182 and one blood vessel region 183, and the regions 181, 182 and 183 are separated as shown in FIG. 18B. 18C and 18D show cross-sectional images of connection portions of the regions 181, 182, and 183. For example, when one blood vessel (183 in FIG. 18A) branches into two blood vessels (181 and 182 in FIG. 18A), a tomographic image (d) showing only one blood vessel (FIG. 18). (D)) is corrected. In this case, in the tomographic image (e) shown in FIG. 18 (e), the logical product region of the concentric circle 184 of the region 181 and one region of the tomographic image (d) and the concentric circle 185 of the tomographic image 182 and the tomographic image (d) Leave the logical product area, and modify others. The corrected result is shown in FIG. As a result, the regions 181, 182, and 183 are smoothly joined. An example of a concentric circle setting method is shown in FIG. The shortest distance dR of the edge portion between the region 181 and the region 182 is obtained by the equation (1):
dR = RO−rc1−rc2 (1)
Here, RO is the distance between the centers of the region 181 and the region 182, rc 1 is the radius of the region 181, and rc 2 is the radius of the region 182.

このようにして求めたdRを用いて、領域181に対応する同心円の半径R1及び領域182に対応する同心円の半径R2は、それぞれ(2)式、(3)式によって求められる:
R1=rc1+0.5×dR (2)
R2=rc2+0.5×dR (3)
(ステップ2B)
断層像や三次元画像は、CRT15に表示される。
Using the thus obtained dR, the radius R1 of the concentric circle corresponding to the region 181 and the radius R2 of the concentric circle corresponding to the region 182 are obtained by equations (2) and (3), respectively:
R1 = rc1 + 0.5 x dR (2)
R2 = rc2 + 0.5 x dR (3)
(Step 2B)
A tomographic image or a three-dimensional image is displayed on the CRT 15.

次に、三次元画像に合わせて表示するカテーテルやステントなどの応用例について説明する。   Next, application examples such as a catheter and a stent displayed in accordance with a three-dimensional image will be described.

図20では、撮影方向の異なるレントゲン画像(二次元画像)201、202と、このレントゲン画像を基に構成した三次元画像203を並べて表示している。二次元画像201、202にはカテーテルを重畳されて表示しており、これらの二次元画像201、202を用いて三次元画像を構成すると、カテーテルの血管内の位置関係が3次元画像203に重畳されて、例えば点線のように表示される。   In FIG. 20, X-ray images (two-dimensional images) 201 and 202 having different imaging directions and a three-dimensional image 203 configured based on these X-ray images are displayed side by side. A catheter is superimposed on the two-dimensional images 201 and 202, and when these two-dimensional images 201 and 202 are used to form a three-dimensional image, the positional relationship of the catheter in the blood vessel is superimposed on the three-dimensional image 203. For example, it is displayed like a dotted line.

図21では血管の三次元画像211にステント三次元画像212を重畳して表示している。ステント3次元画像212は、複数種類のステント画像が磁気ディスクに記憶されている。この重畳表示の起動法は、例えば「ステント3D」という画面に表示されたソフトスイッチ213をマウスで制御すると、血管とステントの位置関係が血管三次元画像211に重畳して表示される。また、ステント像との重畳の対象は、2次元のレントゲン画像であって、前記三次元画像との重畳画像と同時表示してもよい。これによって、ステントの立体的かつ平面的な形状の把握が可能となる。   In FIG. 21, a stent three-dimensional image 212 is superimposed on a three-dimensional blood vessel image 211 and displayed. In the stent three-dimensional image 212, a plurality of types of stent images are stored on the magnetic disk. For example, when the soft switch 213 displayed on the screen “Stent 3D” is controlled with a mouse, the positional relationship between the blood vessel and the stent is displayed superimposed on the blood vessel three-dimensional image 211. The target of superimposition with the stent image is a two-dimensional X-ray image, and may be displayed simultaneously with the superimposition image with the three-dimensional image. As a result, the three-dimensional and planar shape of the stent can be grasped.

本実施形態によれば、CPU11にレントゲン装置1Bによって撮影された被検体の複数の異なる角度のレントゲン画像のそれぞれを二値化して複数の二値化画像を作成する二値化画像作成機能と、この二値化画像作成機能によって作成された前記複数の二値化画像の互いに対応する位置に所定方向に順次移動する走査線を設定する走査線設定機能と、この走査線設定機能によって設定された走査線が順次移動する位置における前記複数の二値化画像の画素点を逆投影することにより断層像を構成する断層像構成機能を備えるので、被検体へのX線被曝の低減に配慮し、より少ない枚数、つまり2枚のレントゲン画像を用いて被検体の断層像を得ることができる。ここでのレントゲン画像は、平面検出器等を用いて撮影されたデジタル画像や、フィルム画像をデジタル化した画像を含む。また、上記より少ない枚数は、3枚以上のレントゲン画像であってもよい。この場合は、上記X線被曝の低減の効果は減少するものの、結果として構成される断層像の精度が向上される。また、本実施形態は、二値化画像を作成してから断層像を得るための逆投影処理を行う手順を行っていたが、逆投影処理を行ってから二値化画像を構成する手順でも、求める断層像は構成可能である。   According to the present embodiment, the binarized image creating function for creating a plurality of binarized images by binarizing each of a plurality of X-ray images of the subject taken by the X-ray apparatus 1B on the CPU 11; A scanning line setting function for setting scanning lines that sequentially move in predetermined directions to positions corresponding to each other of the plurality of binarized images created by the binarized image creating function, and the scanning line setting function Since it has a tomographic image forming function for forming a tomographic image by back projecting the pixel points of the plurality of binarized images at the position where the scanning line sequentially moves, considering reduction of X-ray exposure to the subject, A tomographic image of the subject can be obtained using a smaller number of X-ray images, that is, two X-ray images. The X-ray image here includes a digital image taken using a flat detector or the like, or an image obtained by digitizing a film image. Further, the smaller number than the above may be three or more X-ray images. In this case, although the effect of reducing the X-ray exposure is reduced, the accuracy of the resulting tomographic image is improved. Further, in the present embodiment, the procedure for performing the back projection processing for obtaining the tomographic image after creating the binarized image has been performed, but the procedure for constructing the binarized image after performing the back projection processing is also performed. The desired tomographic image can be constructed.

また、本実施形態によれば、前記走査線設定機能は、前記断層像構成機能によって構成された断層像のそれぞれについて、前記走査線の移動の正方向と逆方向における断層像との相関に基づき偽領域を検出する偽領域検出機能と、この偽領域検出機能によって検出された偽領域を削除する。   Further, according to the present embodiment, the scanning line setting function is based on the correlation between the tomographic image in the forward direction and the reverse direction of the scanning line movement for each of the tomographic images formed by the tomographic image forming function. A false area detection function for detecting a false area and a false area detected by the false area detection function are deleted.

これにより、走査線の移動方向が正方向の断層像と走査線の移動方向が逆方向の断層像との相関に基づいて偽領域を削除するので、より真正な領域だけを用いた正確な断層像が得られる。前記相関は所定の値に満たない領域を検出すれば、その領域が偽領域であるので、効率的に偽領域を削除することができる。   As a result, the false region is deleted based on the correlation between the tomographic image in which the scanning line movement direction is the positive direction and the tomographic image in which the scanning line movement direction is the reverse direction. An image is obtained. If an area where the correlation is less than a predetermined value is detected, the area is a false area, so that the false area can be efficiently deleted.

また、本実施形態によれば、前記走査線設定機能は、前記走査線上で画素値が所定の値をとる領域が1個の位置、及び操作者により指定された位置のいずれかの位置を走査線の初期位置として設定する。   According to the present embodiment, the scanning line setting function scans one of the positions where the pixel value has a predetermined value on the scanning line and the position designated by the operator. Set as the initial position of the line.

これにより、走査線の初期位置は任意に設定できることになるから、抽出したい領域の近傍に設定すれば、より効率的に断層像を構成のための2値化領域を移動させることができるから、結果として断層像の構成のための演算を高速化できる。   As a result, the initial position of the scanning line can be arbitrarily set, so if it is set in the vicinity of the region to be extracted, the binarized region for constructing the tomographic image can be moved more efficiently. As a result, the computation for constructing the tomographic image can be speeded up.

更にまた、本形態によれば、前記設定手段は、前記走査線が所定の位置に達したときは前記走査線方向を反転させて走査線の位置を設定する。このように処理を繰り返すことで、より正確な断層像を得ることができる。   Furthermore, according to this embodiment, when the scanning line reaches a predetermined position, the setting means reverses the scanning line direction and sets the position of the scanning line. By repeating the process in this way, a more accurate tomographic image can be obtained.

所定の位置としては、レントゲン画像での断層像構成範囲の端を用いることができる。   As the predetermined position, the end of the tomographic image construction range in the X-ray image can be used.

(第二の実施形態)
上記第一の実施形態では、ライン毎に断層像の再構成と偽領域の削除を行う場合について説明しているが、まず所定の範囲について断層像を構成し、その後断層像同士の相関に基づいて偽領域を削除するようにしてもよい。以下に、この場合の処理について説明する。
(Second embodiment)
In the first embodiment, the case where tomographic image reconstruction and false region deletion are performed for each line has been described. First, a tomographic image is formed for a predetermined range, and then based on the correlation between the tomographic images. The false area may be deleted. The process in this case will be described below.

図22に、本実施形態に係る断層像と三次元画像の再構成手順を示す。図23に、再構成した断層像の相関処理の原理を示す。図24乃至28は、図23中の逆投影像の例である。   FIG. 22 shows a reconstruction procedure of a tomographic image and a three-dimensional image according to this embodiment. FIG. 23 shows the principle of correlation processing of the reconstructed tomographic image. 24 to 28 are examples of backprojected images in FIG.

(ステップ220)
図23に示すように、レントゲン画像31は、目的臓器を1、その他の領域を0とする二値化処理が行われる。同様に、レントゲン画像32は、目的臓器を1、その他の領域を0とする二値化処理が行われる。
(Step 220)
As shown in FIG. 23, the X-ray image 31 is subjected to binarization processing in which the target organ is 1 and the other areas are 0. Similarly, the X-ray image 32 is subjected to binarization processing in which the target organ is 1 and the other areas are 0.

(ステップ221)
二値化画像は、図23のように、逆投影処理が行われる。逆投影された領域の形状は、図23(a)のように、投影線に囲まれる多角形でも、図23(b)のように、円形近似してもよい。この円形近似では、円形近似後に領域拡張法を使って、絞ってもよい。なお、投影結果が1領域の場合、拡張開始点としてもよい。
(Step 221)
The binarized image is subjected to back projection processing as shown in FIG. The shape of the back-projected region may be a polygon surrounded by a projection line as shown in FIG. 23A or a circular approximation as shown in FIG. In this circular approximation, the region expansion method may be used to narrow down after circular approximation. When the projection result is one area, the expansion start point may be used.

具体的には、図24乃至図28のように、レントゲン画像31については各投影中心241乃至281で投影し、レントゲン画像32については各投影中心242乃至282で投影する。   Specifically, as shown in FIGS. 24 to 28, the X-ray image 31 is projected at the projection centers 241 to 281 and the X-ray image 32 is projected at the projection centers 242 to 282.

(ステップ222)
逆投影画像は上記逆投影された結果から再構成され、その再構成された断層像は磁気ディスク等に保存する(画像a)。
(Step 222)
The backprojected image is reconstructed from the result of the backprojection, and the reconstructed tomographic image is stored on a magnetic disk or the like (image a).

(ステップ223)
断層像の再構成ラインは、n(n=1,2等)更新される。
(Step 223)
The tomographic reconstruction line is updated by n (n = 1, 2, etc.).

(ステップ224)
断層像の再構成ラインは全ライン終了しているか否かを判定する。終了していればステップ225に進み、終了していなければステップ221に戻る。
(Step 224)
It is determined whether or not all the reconstruction lines of the tomographic image have been completed. If completed, the process proceeds to step 225. If not completed, the process returns to step 221.

断層像は、上記更新ライン毎に二値化、逆投影を繰り返し、それぞれの更新ライン毎に得られる。   A tomographic image is obtained for each update line by repeating binarization and back projection for each update line.

(ステップ225)
ステップ221からステップ224までを例えば、図24、図25、図26、図27、図28まで繰り返した結果、逆投影により生じた領域が一箇所のところを探す。この場合では、図28の結果が探されることになる。
(Step 225)
Steps 221 to 224 are repeated until, for example, FIG. 24, FIG. 25, FIG. 26, FIG. 27, and FIG. In this case, the result of FIG. 28 is searched.

(ステップ226)
ステップ225で探した領域が一箇所である図28の領域283から図27の領域273、図26の領域263、図25の領域253、図24の領域243と連続性のある領域を隣接画像について探す。そして、連続性がない領域は偽領域として図24、図25、図26から削除する。また、この削除の際に、逆投影後最後に一定の割合以上相関のある画像は残しても良い。
(Step 226)
The region searched in Step 225 is a single region from the region 283 in FIG. 28 to the region 273 in FIG. 27, the region 263 in FIG. 26, the region 253 in FIG. 25, and the region 243 in FIG. look for. And the area | region which does not have continuity is deleted from FIG.24, FIG.25, FIG.26 as a false area | region. Further, at the time of this deletion, an image having a correlation of a certain ratio or more may be left at the end after back projection.

このようにして、断層像構成装置10では、偽領域が削除された正確な断層像を得ることができる。   In this way, the tomographic image construction apparatus 10 can obtain an accurate tomographic image from which the false region is deleted.

(ステップ227)
偽領域を削除した断層像は保存される。
(Step 227)
The tomographic image from which the false region is deleted is stored.

(ステップ228)
二値化画像のラインは三次元画像再構成に必要な数の断層像が得られるまで、更新される。
(Step 228)
The lines of the binarized image are updated until the number of tomographic images necessary for three-dimensional image reconstruction is obtained.

(ステップ229)
ステップ228のラインの更新は、全ラインが終了したか否かを判定する。全ラインが終了していればステップ22Aに進み、全ラインが終了していなければステップ226に戻る。
(Step 229)
The line update in step 228 determines whether all lines have been completed. If all lines have been completed, the process proceeds to step 22A. If all lines have not been completed, the process returns to step 226.

(ステップ22A)
断層像の構成処理は走査線方向を反転して繰り返す処理が終了したか否かを判定する。繰り返し処理が終了すればステップ22Cに進み、繰り返し処理が終了していなければステップ22Bに進む。
(Step 22A)
In the tomographic image construction process, it is determined whether or not the process of reversing the scanning line direction has been completed. If the repetition process is completed, the process proceeds to step 22C. If the repetition process is not completed, the process proceeds to step 22B.

(ステップ22B)
断層像の構成処理は走査線方向を反転して行う。この走査線方向の反転とは、二値化画像に対し、例えば最初に画像上部から下部へ走査線を更新した場合、次に画像下部から上部へ走査線を更新する。
(Step 22B)
The tomographic image construction process is performed by reversing the scanning line direction. The reversal of the scanning line direction means that when a scanning line is first updated from the upper part to the lower part of the binarized image, for example, the scanning line is updated from the lower part to the upper part of the image.

(ステップ22C)
三次元画像は、ステップ227に保存された断層像から構成される
(ステップ22D)
断層像や三次元画像は、CRT15に表示される。
(Step 22C)
The three-dimensional image is composed of tomographic images stored in step 227 (step 22D).
A tomographic image or a three-dimensional image is displayed on the CRT 15.

本実施形態に示す手順は、第一の実施形態の場合と同様にレントゲン画像から断層像を再構成及び偽領域の削除を行うことができる。   The procedure shown in this embodiment can reconstruct a tomographic image from an X-ray image and delete a false region as in the case of the first embodiment.

また、本実施形態によれば、図29に示されるように、3次元画像291、3次元画像291の投影中心をX方向に移動した三次元画像292、三次元画像292を反転させた三次元画像293や、図30に示されるように、3次元画像291をY方向に投影中心を移動した三次元画像301、三次元画像301の投影中心をX方向に移動した三次元画像302、三次元画像302を反転させた三次元画像303など、多様な三次元画像を構成することができる。また、枝分かれした複雑な血管などでも細部まで抽出し、断層像や三次元画像を構成することができる。   In addition, according to the present embodiment, as shown in FIG. 29, a three-dimensional image 291 obtained by moving the projection center of the three-dimensional image 291 and the three-dimensional image 291 in the X direction, and a three-dimensional image obtained by inverting the three-dimensional image 292. 30, a three-dimensional image 301 obtained by moving the projection center of the three-dimensional image 291 in the Y direction, a three-dimensional image 302 obtained by moving the projection center of the three-dimensional image 301 in the X direction, as illustrated in FIG. Various three-dimensional images such as a three-dimensional image 303 obtained by inverting the image 302 can be formed. In addition, it is possible to extract tomograms and three-dimensional images by extracting even the details of branched complex blood vessels.

なお、断層像の再構成に用いるレントゲン画像の枚数の設定や、再構成した断層像の三次元画像への応用は、第一の実施形態と同様に行うことができる。   The setting of the number of X-ray images used for the reconstruction of the tomographic image and the application of the reconstructed tomographic image to the three-dimensional image can be performed as in the first embodiment.

以上説明したように本発明によれば、レントゲン画像から断層像を再構成することができる。   As described above, according to the present invention, a tomographic image can be reconstructed from an X-ray image.

図1は、本発明の第一の実施形態に係る断層像構成装置のハードウェア構成を示す図であり;FIG. 1 is a diagram showing a hardware configuration of a tomographic image forming apparatus according to the first embodiment of the present invention; 図2は、断層像構成の処理手順を示す図であり;FIG. 2 is a diagram showing a processing procedure for tomographic image construction; 図3は、レントゲン画像、断層像、三次元画像の位置関係を示す図であり;FIG. 3 is a diagram showing the positional relationship between an X-ray image, a tomographic image, and a three-dimensional image; 図4は、二値化画像の走査線方向と1領域の関係の例を示す図であり;FIG. 4 is a diagram showing an example of the relationship between the scanning line direction of a binarized image and one region; 図5は、ランレングス分布の例を示す図であり;FIG. 5 is a diagram showing an example of a run length distribution; 図6は、二値化画像で1領域が一個所の場合の逆投影像を示す図であり;FIG. 6 is a diagram showing a back-projected image in the case of a binarized image with one area. 図7(a)、7(b)及び7(c)は、二値化画像における逆投影開始位置とその移動を示す図であり;7 (a), 7 (b) and 7 (c) are diagrams showing the back projection start position and its movement in the binarized image; 図8は、目的臓器を一個所にする場合にランレングスを利用した例を示す図であり;FIG. 8 is a diagram showing an example of using run length when the target organ is one; 図9は、3個の被検体とレントゲン画像の関係の例を示す図であり;FIG. 9 is a diagram showing an example of the relationship between three subjects and X-ray images; 図10は、逆投影開始位置をROIで囲んだ場合の例を示す図であり;FIG. 10 is a diagram showing an example of a case where the back projection start position is surrounded by an ROI; 図11は、逆投影像を円形近似した例を示す図であり;FIG. 11 is a diagram showing an example of circular approximation of a backprojected image; 図12は、再構成領域を楕円で近似できる場合を示す図であり;FIG. 12 shows a case where the reconstruction area can be approximated by an ellipse; 図13は、断層像のコーンビーム再構成の例を示す図であり;FIG. 13 is a diagram showing an example of cone beam reconstruction of a tomographic image; 図14は、断層像の近似的再構成を示す図であり;FIG. 14 is a diagram showing an approximate reconstruction of a tomogram; 図15は、図10のレントゲン画像を二値化して逆投影した場合を示す図であり;FIG. 15 is a diagram showing a case where the X-ray image of FIG. 10 is binarized and back-projected; 図16は、走査線方向を決めるフローチャートであり;FIG. 16 is a flowchart for determining the scan line direction; 図17は、分岐のある三次元画像の例を示す図であり;FIG. 17 is a diagram showing an example of a branched three-dimensional image; 図18(a)乃至18(f)は、分岐部の補正処理の例を示す図であり;18 (a) to 18 (f) are diagrams showing an example of correction processing of a branch portion; 図19は、分岐部の補正処理における同心円の設定方法を示す図であり;FIG. 19 is a diagram showing a concentric circle setting method in the correction processing of the bifurcation; 図20は、三次元画像とカテーテル像との重畳表示の例であり;FIG. 20 is an example of superimposed display of a three-dimensional image and a catheter image; 図21は、三次元画像とステント像との重畳表示の例であり;FIG. 21 is an example of a superimposed display of a three-dimensional image and a stent image; 図22は、断層像構成の処理手順を示す図であり;FIG. 22 is a diagram showing a processing procedure for constructing a tomographic image; 図23は、全ラインについての逆投影後の相関処理を示す図であり;FIG. 23 is a diagram showing correlation processing after back projection for all lines; 図24は、二値化画像に1領域が複数ある場合の逆投影像の例を示す図であり;FIG. 24 is a diagram showing an example of a backprojected image when there are a plurality of one regions in the binarized image; 図25は、二値化画像に1領域が複数ある場合の逆投影像の他の例を示す図であり;FIG. 25 is a diagram showing another example of the backprojected image when there are a plurality of one regions in the binarized image; 図26は、二値化画像に1領域が複数ある場合の逆投影像の他の例を示す図であり;FIG. 26 is a diagram showing another example of a backprojected image when there are a plurality of one regions in the binarized image; 図27は、二値化画像に1領域が一個所ある場合の逆投影像の例を示す図であり;FIG. 27 is a diagram showing an example of a backprojected image when one region is present in the binarized image; 図28は、二値化画像に1領域が一個所ある場合の逆投影像の他の例を示す図であり;FIG. 28 is a diagram showing another example of a backprojected image in the case where there is one area in the binarized image; 図29は、断層像から構成した三次元画像を横方向から見た場合の例を示す図であり;FIG. 29 is a diagram showing an example when a three-dimensional image composed of tomographic images is viewed from the lateral direction; 図30は、断層像から構成した三次元画像を斜め方向から見た場合の例を示す図である。FIG. 30 is a diagram illustrating an example when a three-dimensional image configured from tomographic images is viewed from an oblique direction.

Claims (11)

レントゲン装置によって撮影された被検体の複数の異なる方向からのレントゲン画像のそれぞれを二値化して複数の二値化画像を作成する二値化画像作成手段と、
前記二値化画像作成手段によって作成された前記複数の二値化画像の互いに対応する位置に所定方向に順次移動する走査線を設定する走査線設定手段と、
前記走査線設定手段によって設定された走査線が順次移動する位置における前記複数の二値化画像の画素点を逆投影することにより断層像を構成する断層像構成手段と、を備え、
前記走査線設定手段は、2次元の走査線方向についてランレングスの出現頻度を演算し、前記演算された出現頻度に基づいて、2次元の走査線方向のうちの一方を設定する走査線方向設定手段を含むことを特徴とする断層像構成装置。
Binarized image creating means for binarizing each of a plurality of X-ray images from a plurality of different directions of the subject imaged by the X-ray apparatus to create a plurality of binarized images;
Scanning line setting means for setting scanning lines that sequentially move in a predetermined direction to positions corresponding to each other of the plurality of binarized images created by the binarized image creating means;
A tomographic image forming unit configured to form a tomographic image by back projecting pixel points of the plurality of binarized images at positions where the scanning lines set by the scanning line setting unit sequentially move,
The scanning line setting means calculates a run length appearance frequency in a two-dimensional scanning line direction, and sets one of the two-dimensional scanning line directions based on the calculated appearance frequency. A tomographic image forming apparatus comprising means.
前記走査線設定手段は、前記断層像構成手段によって構成された断層像のそれぞれについて、前記走査線の移動の正方向と逆方向における断層像との相関に基づき偽領域を検出する偽領域検出手段と、前記偽領域検出手段によって検出された偽領域を削除する偽領域削除手段と、を含むことを特徴とする請求項1に記載の断層像構成装置。  The scanning line setting means detects a false area for each of the tomographic images formed by the tomographic image forming means based on the correlation between the tomographic image in the forward direction and the reverse direction of the movement of the scanning line. The tomographic image forming apparatus according to claim 1, further comprising: a false area deleting unit that deletes the false area detected by the false area detecting unit. 前記偽領域削除手段は、所定の範囲について複数の断層像を構成し、それらの後断層像同士の相関に基づいて偽領域を削除することを特徴とする請求項2に記載の断層像構成装置。  The tomographic image constructing apparatus according to claim 2, wherein the false region deleting unit forms a plurality of tomographic images for a predetermined range and deletes the false regions based on a correlation between the subsequent tomographic images. . 前記走査線設定手段は、前記走査線上で画素値が所定の値をとる領域が一つである位置、及び操作者により指定された位置のいずれかの位置を走査線の初期位置として設定することを特徴とする請求項1に記載の断層像構成装置。  The scanning line setting unit sets, as an initial position of the scanning line, one of a position where the pixel value takes a predetermined value on the scanning line and a position designated by the operator. The tomographic image forming apparatus according to claim 1, wherein: 前記断層像構成手段は、前記断層像を構成する際に、その目的臓器の形状に応じて前記目的臓器の近似形状を設定する近似形状設定手段を含むことを特徴とする請求項1に記載の断層像構成装置。  2. The tomographic image constructing unit includes an approximate shape setting unit configured to set an approximate shape of the target organ according to the shape of the target organ when the tomographic image is configured. Tomographic image construction device. 前記断層像構成手段によって構成された断層像を隣接して複数構成し、前記複数構成された断層像から所定のレンダリング法により三次元画像を構成する三次元画像構成手段と、
前記三次元画像構成手段によって構成された三次元画像を表示する画像表示手段と、を更に備えることを特徴とする請求項1に記載の断層像構成装置。
A plurality of tomographic images formed by the tomographic image forming means adjacent to each other, and a three-dimensional image forming means for forming a three-dimensional image from the plurality of formed tomographic images by a predetermined rendering method;
2. The tomographic image construction apparatus according to claim 1, further comprising image display means for displaying a three-dimensional image constituted by the three-dimensional image construction means.
前記三次元画像構成手段は、前記複数枚のレントゲン画像に含まれるカテーテル部分から三次元カテーテル像を構成し、前記画像表示手段は、前記構成された三次元カテーテル像と前記三次元画像構成手段によって構成された三次元画像とを重畳して表示することを特徴とする請求項6に記載の断層像構成装置。  The three-dimensional image constructing means composes a three-dimensional catheter image from the catheter portions included in the plurality of X-ray images, and the image display means comprises the constructed three-dimensional catheter image and the three-dimensional image constructing means. The tomographic image forming apparatus according to claim 6, wherein the tomographic image forming apparatus displays the superimposed three-dimensional image in a superimposed manner. レントゲン装置によって撮影された被検体の複数の異なる方向からのレントゲン画像のそれぞれを二値化して複数の二値化画像を作成する二値化画像作成手段と、
前記二値化画像作成手段によって作成された前記複数の二値化画像の互いに対応する位置に所定方向に順次移動する走査線を設定する走査線設定手段であって、2次元の走査線方向についてランレングスの出現頻度を演算し、前記演算された出現頻度に基づいて、2次元の走査線方向のうちの一方を設定する走査線方向設定手段を含む走査線設定手段と、
前記走査線設定手段によって設定された走査線が順次移動する位置における前記複数の二値化画像の画素点を逆投影することにより断層像を構成する断層像構成手段と、
前記断層像構成手段によって構成された断層像を隣接して複数構成し、前記複数構成された断層像から所定のレンダリング法により三次元画像を構成する三次元画像構成手段と、
前記三次元画像構成手段によって構成された三次元画像を表示する画像表示手段と、
前記三次元画像構成手段によって構成された三次元画像の分岐部を分離演算する分離演算手段と、
前記分離演算手段によって分離演算された分岐部の同心円を求め、前記同心円に基づいて前記分岐部を滑らかに結合演算する結合演算手段と、を備える断層像構成装置。
Binarized image creating means for binarizing each of a plurality of X-ray images from a plurality of different directions of the subject imaged by the X-ray apparatus to create a plurality of binarized images;
Scanning line setting means for setting scanning lines that sequentially move in a predetermined direction to positions corresponding to each other of the plurality of binarized images created by the binarized image creating means, the scanning line setting direction being two-dimensional A scanning line setting unit including a scanning line direction setting unit configured to calculate a run length appearance frequency and set one of the two-dimensional scanning line directions based on the calculated appearance frequency ;
A tomographic image constructing means for constructing a tomographic image by back projecting pixel points of the plurality of binarized images at positions where the scanning lines set by the scanning line setting means sequentially move;
A plurality of tomographic images formed by the tomographic image forming means adjacent to each other, and a three-dimensional image forming means for forming a three-dimensional image from the plurality of formed tomographic images by a predetermined rendering method;
Image display means for displaying a three-dimensional image constituted by the three-dimensional image construction means;
Separation calculation means for separating and calculating a branch part of the three-dimensional image constituted by the three-dimensional image construction means;
A tomographic image construction apparatus comprising: a concatenation computing unit that obtains concentric circles of the branch portion separated and computed by the separation computation unit, and smoothly joins and computes the branch portions based on the concentric circles.
レントゲン装置によって撮影された被検体の複数の異なる方向からのレントゲン画像のそれぞれを二値化して複数の二値化画像を作成する二値化画像作成手段と、
前記二値化画像作成手段によって作成された前記複数の二値化画像の互いに対応する位置に所定方向に順次移動する走査線を設定する走査線設定手段であって、2次元の走査線方向についてランレングスの出現頻度を演算し、前記演算された出現頻度に基づいて、2次元の走査線方向のうちの一方を設定する走査線方向設定手段を含む走査線設定手段と、
前記走査線設定手段によって設定された走査線が順次移動する位置における前記複数の二値化画像の画素点を逆投影することにより断層像を構成する断層像構成手段と、
前記断層像構成手段によって構成された断層像を隣接して複数構成し、前記複数構成された断層像から所定のレンダリング法により三次元画像を構成する三次元画像構成手段と、
前記三次元画像構成手段によって構成された三次元画像を表示する画像表示手段と、
複数のステント像を記憶するステント像記憶手段と、
前記ステント像記憶手段によって記憶された複数のステントから所望のステント像を読出すステント像読出手段と、を備え、
前記画像表示手段は、前記ステント像読出手段によって読み出されたステント像と前記三次元画像構成手段によって構成された三次元画像とを重畳して表示することを特徴とする断層像構成装置。
Binarized image creating means for binarizing each of a plurality of X-ray images from a plurality of different directions of the subject imaged by the X-ray apparatus to create a plurality of binarized images;
Scanning line setting means for setting scanning lines that sequentially move in a predetermined direction to positions corresponding to each other of the plurality of binarized images created by the binarized image creating means, the scanning line setting direction being two-dimensional A scanning line setting unit including a scanning line direction setting unit configured to calculate a run length appearance frequency and set one of the two-dimensional scanning line directions based on the calculated appearance frequency ;
A tomographic image constructing means for constructing a tomographic image by back projecting pixel points of the plurality of binarized images at positions where the scanning lines set by the scanning line setting means sequentially move;
A plurality of tomographic images formed by the tomographic image forming means adjacent to each other, and a three-dimensional image forming means for forming a three-dimensional image from the plurality of formed tomographic images by a predetermined rendering method;
Image display means for displaying a three-dimensional image constituted by the three-dimensional image construction means;
A stent image storage means for storing a plurality of stent images;
A stent image reading means for reading a desired stent image from a plurality of stents stored by the stent image storage means,
The tomographic image construction apparatus, wherein the image display means superimposes and displays the stent image read by the stent image readout means and the 3D image constituted by the 3D image construction means.
レントゲン装置によって撮影された被検体の複数の異なる方向からのレントゲン画像のそれぞれを二値化して複数の二値化画像を作成する二値化画像作成工程と、
前記二値化画像作成工程によって作成された前記複数の二値化画像の互いに対応する位置に所定方向に順次移動する走査線を設定し、2次元の走査線方向についてランレングスの出現頻度を演算し、前記演算された出現頻度に基づいて、2次元の走査線方向のうちの一方を設定する走査線設定工程と、
前記走査線設定工程によって設定された走査線が順次移動する位置における前記複数の二値化画像の画素点を逆投影することにより断層像を構成する断層像構成工程と、
前記断層像構成工程によって構成された断層像を隣接して複数構成し、前記複数構成された断層像から所定のレンダリング法により三次元画像を構成する三次元画像構成工程と、
前記三次元画像構成工程によって構成された三次元画像を表示する画像表示工程と、
前記三次元画像構成工程によって構成された三次元画像の分岐部を分離演算する分離演算工程と、
前記分離演算工程によって分離演算された分岐部の同心円を求め、前記同心円に基づいて前記分岐部を滑らかに結合演算する結合演算工程と、を備える断層像構成方法。
A binarized image creating step for binarizing each of the X-ray images from a plurality of different directions of the subject imaged by the X-ray device to create a plurality of binarized images;
A scan line that sequentially moves in a predetermined direction is set at a position corresponding to each other of the plurality of binarized images created by the binarized image creating step, and the appearance frequency of run length is calculated in the two-dimensional scan line direction. A scanning line setting step for setting one of the two-dimensional scanning line directions based on the calculated appearance frequency ;
A tomographic image forming step of forming a tomographic image by back projecting pixel points of the plurality of binarized images at positions where the scanning lines set by the scanning line setting step sequentially move;
A plurality of tomographic images configured by the tomographic image configuration step are configured adjacently, and a three-dimensional image configuration step of configuring a three-dimensional image from the plurality of configured tomographic images by a predetermined rendering method;
An image display step for displaying the three-dimensional image constructed by the three-dimensional image construction step;
A separation calculation step of separating and calculating a branch portion of the three-dimensional image constituted by the three-dimensional image construction step;
A tomographic image constructing method comprising: a concatenation calculation step of obtaining concentric circles of the branch portions subjected to the separation operation in the separation operation step and smoothly combining the branch portions based on the concentric circles.
レントゲン装置によって撮影された被検体の複数の異なる方向からのレントゲン画像のそれぞれを二値化して複数の二値化画像を作成する二値化画像作成工程と、
前記二値化画像作成工程によって作成された前記複数の二値化画像の互いに対応する位置に所定方向に順次移動する走査線を設定し、2次元の走査線方向についてランレングスの出現頻度を演算し、前記演算された出現頻度に基づいて、2次元の走査線方向のうちの一方を設定する走査線設定工程と、
前記走査線設定工程によって設定された走査線が順次移動する位置における前記複数の二値化画像の画素点を逆投影することにより断層像を構成する断層像構成工程と、
前記断層像構成工程によって構成された断層像を隣接して複数構成し、前記複数構成された断層像から所定のレンダリング法により三次元画像を構成する三次元画像構成工程と、
前記三次元画像構成工程によって構成された三次元画像を表示する画像表示工程と、
複数のステント像を記憶するステント像記憶工程と、
前記ステント像記憶工程によって記憶された複数のステントから所望のステント像を読出すステント像読出工程と、を備え、
前記画像表示工程は、前記ステント像読出工程によって読み出されたステント像と前記三次元画像構成工程によって構成された三次元画像とを重畳して表示することを特徴とする断層像構成方法。
A binarized image creating step for binarizing each of the X-ray images from a plurality of different directions of the subject imaged by the X-ray device to create a plurality of binarized images;
A scan line that sequentially moves in a predetermined direction is set at a position corresponding to each other of the plurality of binarized images created by the binarized image creating step, and the appearance frequency of run length is calculated in the two-dimensional scan line direction. A scanning line setting step for setting one of the two-dimensional scanning line directions based on the calculated appearance frequency ;
A tomographic image forming step of forming a tomographic image by back projecting pixel points of the plurality of binarized images at positions where the scanning lines set by the scanning line setting step sequentially move;
A plurality of tomographic images configured by the tomographic image configuration step are configured adjacently, and a three-dimensional image configuration step of configuring a three-dimensional image from the plurality of configured tomographic images by a predetermined rendering method;
An image display step for displaying the three-dimensional image constructed by the three-dimensional image construction step;
A stent image storing step of storing a plurality of stent images;
A stent image reading step of reading a desired stent image from a plurality of stents stored by the stent image storage step,
The tomographic image construction method characterized in that the image display step superimposes and displays the stent image read by the stent image readout step and the three-dimensional image constructed by the three-dimensional image construction step.
JP2005512615A 2003-08-05 2004-08-04 Tomographic image constructing apparatus and method Expired - Fee Related JP4467522B2 (en)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2003287013 2003-08-05
JP2003287013 2003-08-05
PCT/JP2004/011494 WO2005011499A1 (en) 2003-08-05 2004-08-04 Tomogram constituting system and method

Publications (2)

Publication Number Publication Date
JPWO2005011499A1 JPWO2005011499A1 (en) 2006-09-14
JP4467522B2 true JP4467522B2 (en) 2010-05-26

Family

ID=34113991

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2005512615A Expired - Fee Related JP4467522B2 (en) 2003-08-05 2004-08-04 Tomographic image constructing apparatus and method

Country Status (2)

Country Link
JP (1) JP4467522B2 (en)
WO (1) WO2005011499A1 (en)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7892177B2 (en) 2005-02-28 2011-02-22 Scimed Life Systems, Inc. Systems and methods for estimating the length and position of a stent to be applied within a patient
CN101711681B (en) * 2008-10-07 2012-06-27 株式会社东芝 Three-dimensional image processing apparatus
JP5324880B2 (en) * 2008-10-21 2013-10-23 アズビル株式会社 CT apparatus and metal shape extraction method
JP5324883B2 (en) * 2008-10-24 2013-10-23 アズビル株式会社 CT apparatus and metal shape extraction method
JP5597399B2 (en) * 2010-01-08 2014-10-01 株式会社東芝 Medical diagnostic imaging equipment
US20140031676A1 (en) * 2011-04-12 2014-01-30 Koninklijke Philips N.V. Embedded 3d modelling
JP5989312B2 (en) * 2011-08-18 2016-09-07 東芝メディカルシステムズ株式会社 Image processing display device and image processing display program
CN103714513B (en) * 2012-09-29 2017-05-31 清华大学 Artifact correction method and equipment in CT imagings
CN103713329B (en) * 2012-09-29 2016-12-21 清华大学 CT imaging positions method and the equipment of object
JP6425906B2 (en) * 2014-03-19 2018-11-21 キヤノンメディカルシステムズ株式会社 Medical image diagnostic apparatus and medical image processing apparatus
EP3247301B1 (en) * 2015-01-22 2020-10-28 Koninklijke Philips N.V. Endograft visualization with optical shape sensing
US11515031B2 (en) 2018-04-16 2022-11-29 Canon Medical Systems Corporation Image processing apparatus, X-ray diagnostic apparatus, and image processing method

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03240169A (en) * 1990-02-16 1991-10-25 Hitachi Medical Corp Center positioning method for reformed image of medical image diagnosis device
JP3166232B2 (en) * 1991-08-31 2001-05-14 株式会社島津製作所 X-ray equipment
JPH10234717A (en) * 1997-02-28 1998-09-08 Shimadzu Corp X-ray photographing device
DE19919907C2 (en) * 1999-04-30 2003-10-16 Siemens Ag Method and device for catheter navigation in three-dimensional vascular tree images
JP4405002B2 (en) * 1999-09-10 2010-01-27 阿部 慎一 Stent graft design device

Also Published As

Publication number Publication date
JPWO2005011499A1 (en) 2006-09-14
WO2005011499A1 (en) 2005-02-10

Similar Documents

Publication Publication Date Title
US7742629B2 (en) System and method for three-dimensional reconstruction of a tubular organ
US7574026B2 (en) Method for the 3d modeling of a tubular structure
US9659375B2 (en) Methods and systems for transforming luminal images
JP4467522B2 (en) Tomographic image constructing apparatus and method
US6996205B2 (en) Methods and apparatus to facilitate review of CT colonography exams
US20080205722A1 (en) Method and Apparatus for Automatic 4D Coronary Modeling and Motion Vector Field Estimation
JP2006246941A (en) Image processing apparatus and vessel tracking method
US9582934B2 (en) Method and system for efficient extraction of a silhouette of a 3D mesh
US7689018B2 (en) Anomaly detection in volume data structure information
WO2022109906A1 (en) Method and system for rendering and synthesizing three-dimensional blood vessel, and storage medium
JP2004008304A (en) Method for generating and displaying three-dimensional shape using multidirectional projection image
JP4493436B2 (en) Image interpretation support method, apparatus and program
US8165375B2 (en) Method and system for registering CT data sets
JP4336083B2 (en) Image diagnosis support apparatus and image diagnosis support method
JP2019010382A (en) Image positioning apparatus, method, and program
JPH1031753A (en) Method for preparing three-dimensional image and medical image diagnostic device
JP2011115546A (en) Ultrasonic three-dimensional reconstruction method and device, and ultrasonic system
JP2011182946A (en) Medical image display and medical image display method
NL1044271B1 (en) A method for direct manipulation and visualization of the 3d internal structures of a tubular object as they are in reality without any noticeable distortion
JP2006087827A (en) Diagnostic imaging apparatus and image processing system, program and method
JPS6113374A (en) Diagnostic device of image
JP6878035B2 (en) Medical image processing equipment, medical diagnostic imaging equipment and medical image processing programs
EP4334885A1 (en) Stitching multiple images to create a panoramic image
WO2022171539A1 (en) Displacement estimation of interventional devices
CN114267433A (en) Medical image processing method and medical image processing device

Legal Events

Date Code Title Description
A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20061206

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20061206

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20090713

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20090721

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20091208

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20100115

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

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

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130305

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140305

Year of fee payment: 4

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

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

LAPS Cancellation because of no payment of annual fees