KR100964649B1 - Method for correcting magnification from central axis shift on X-ray computed tomography - Google Patents

Abstract

According to the present invention, a virtual model is prepared, the center point of the virtual model is aligned with the center point of the central axis of the X-ray city photographing apparatus, and placed on the X-ray city photographing apparatus, and the X-ray city photographing apparatus scans the virtual model. Obtain a plurality of detector plan data, extract middle plan data of the virtual model from the detector plan data to obtain measurement lengths of middle plan data, and measure a length and reference length of each of the middle plan data of the virtual model. To obtain a magnification ratio of each of the middle plan data of the virtual model, and then use the magnification ratios to propose a method of correcting the magnification caused by the central axis movement in the X-ray city imaging apparatus. To correct the magnification attributable to the movement of the central axis in the A.

computed tomographic, middle plan, 3D image, central axis

Description

Method for correcting magnification from central axis shift on X-ray computed tomography

The present invention relates to a method for correcting the magnification caused by the movement of the central axis in the X-ray Citigraphy apparatus, and more particularly, preparing a virtual model and matching the center point of the virtual model to the center point of the central axis of the X-ray Citigraphy apparatus. Position the image on the X-ray Citigraphy device, scan the virtual model to obtain a plurality of detector plan data, and extract the middle plan data of the virtual model from the detector plan data. After obtaining the measurement length of the plan data, and comparing the measurement length and reference length of each of the middle plan data of the virtual model to obtain a magnification ratio of each of the middle plan data of the virtual model, using the magnification ratios How to calibrate the center point of the central axis of the X-ray Citigraphy Device The X-ray imaging apparatus City present invention relates to a method for correcting the magnification due to the central axis moves.

Today, devices for acquiring medical images have been developed to more accurately diagnose a patient's condition.

Among them, three-dimensional tomographic systems such as X-ray computed tomography have an advantage of obtaining accurate three-dimensional images. For this reason, three-dimensional tomography systems, such as X-ray systems, have become one of the most widely used equipment for acquiring medical images.

In this case, in order to obtain a 3D image, the X-ray device first acquires a plurality of 2D images, and then reconstructs the 2D images through mathematical calculations to obtain a 3D image.

The 3D image reconstruction technique has been technically developed from the 2D image reconstruction technique.

Representative two-dimensional image reconstruction techniques include parallel-beam reconstruction, equal-angle fan-beam reconstruction, and equal-spaced fan-beam reconstruction.

The 3D projection reconstruction technique is a cone-beam reconstruction technique, also called a FDK (Feldkamp) technique, which extends the z-axis to enable 3D processing based on the 2D equal space fan beam reconstruction technique. It is a three-dimensional reconstruction technique.

In this case, the 3D projection reconstruction technique receives 2D projection data from an X-ray Citi apparatus and reconstructs it into a 3D image through mathematical calculation. The X-ray source and detector of the X-ray Citi apparatus are used. The center point of the detector's central axis is not located at the correct position, but the depth changes, causing problems such as blurring in 3D images.

Accordingly, the present invention is to solve the above-mentioned disadvantages and problems of the prior art, to prepare a virtual model, by matching the center point of the virtual model to the center point of the central axis of the X-ray Citigraphy X-ray imaging A plurality of detector plan data are obtained by scanning the virtual model with the X-ray imaging apparatus, extracting middle plan data of the virtual model from the detector plan data, and measuring a measurement length of the middle plan data. After acquiring a magnification ratio of each of the middle plan data of the virtual model by comparing the measured length and the reference length of each of the middle plan data of the virtual model, and using the magnification ratios, X-rays by suggesting how to correct the center point of the central axis X-ray source that eliminates the cause of inability to sharpen three-dimensional images, such as blurring caused by varying depths, rather than the center of the X-ray source of the device and the center axis of the detector It is an object of the present invention to provide a method for correcting the magnification caused by the movement of a central axis in an imaging device.

The object of the present invention is to prepare a virtual model; Positioning a center of the virtual model on the X-ray Citigraphic apparatus by matching a center point of a central axis of the X-ray Citigraphic apparatus; Acquiring a plurality of detector plan data by scanning the virtual model with the X-ray imaging apparatus; Extracting middle plan data of the virtual model from the detector plan data to obtain a measurement length of middle plan data; Obtaining a magnification ratio of each of the middle plan data of the virtual model by comparing the measured length of each of the middle plan data of the virtual model with a reference length; And correcting the center point of the central axis of the X-ray Citi apparatus using the magnification ratios.

In addition, the object of the present invention is obtained when the reference length is the center of the virtual model is located at the center point of the central axis, the center of the virtual model, the center point of the central axis and the center of the X-ray source is located in a straight line, It is also achieved by a method of correcting the magnification caused by the movement of the central axis in the X-ray Citigraphic apparatus, characterized in that the theoretically calculated length of the measurement of the middle plan data.

In addition, the object of the present invention is to correct the center point of the central axis of the X-ray imaging apparatus using the magnification ratios, the back projection in the process of reconstructing the detector plan data in three dimensions using a filter back-projection algorithm The detector plan data may be corrected by multiplying the respective magnification ratios by the respective magnification ratios immediately before the process.

In addition, the object of the present invention is to correct the center point of the central axis of the X-ray imaging apparatus using the magnification ratios to rearrange the position of the source and the detector of the X-ray imaging apparatus with the obtained magnification ratios. It is also achieved by a method for correcting the magnification caused by the movement of the central axis in the X-ray Citigraphy apparatus characterized in that.

In addition, the object of the present invention is to extract the middle plan data of the virtual model from the detector plan data extracts the data of the row located in the middle of the entire row in each of the detector plan data, the extracted intermediate It is also achieved by a method of correcting the magnification caused by the central axis movement in the X-ray Citigraphic apparatus, characterized in that by extracting the data of the virtual model of the data of the row.

In the X-ray Citigraphy apparatus of the present invention, a method of correcting the magnification caused by the movement of the central axis removes blurring phenomenon in which the 3D image is not clear because the depth of the main axis is not located at the center of the object. In the X-ray Citigraphy apparatus capable of acquiring a clear three-dimensional image during X-ray imaging, there is an effect of providing a method of correcting the magnification caused by the movement of the central axis.

Details of the above objects and technical configurations and the effects thereof according to the present invention will be more clearly understood by the following detailed description with reference to the drawings showing preferred embodiments of the present invention. Like numbers refer to like elements throughout.

In the method of correcting the magnification caused by the movement of the central axis in the X-ray imaging apparatus of the present invention, first, a virtual model having a cross section is prepared. In this case, the virtual model is preferably spherical or cylindrical so that the cross section is a circle. Of course, the virtual model is made of a material that increases the difference between the density and the portion other than the virtual model, it is preferable that the X-rays are made of a material having a low transmittance.

Subsequently, the virtual model is placed on the X-ray city photographing apparatus. In this case, the center point of the virtual model should be positioned to coincide with the center point of the central axis of the X-ray city photographing apparatus.

The method of confirming that the center of the virtual model is positioned to coincide with the center point of the central axis of the X-ray city photographing apparatus may be confirmed by scanning the virtual model with the X-ray city photographing apparatus and reconstructing it into a 3D image. have.

Subsequently, the virtual model is scanned by the X-ray imaging apparatus to obtain a plurality of detector plan data.

Then, middle plan data of the virtual model is extracted from the detector plan data.

In this case, middle plan data of the virtual model may be extracted by extracting data of a row located in the middle of all the rows from the detector plan data, and extracting data of the virtual model from the extracted intermediate rows of data. have.

Then, the length of each of the middle plan data of the virtual model, that is, the measurement length is obtained.

Subsequently, a magnification ratio of each of the middle plan data of the virtual model is obtained by comparing the measured length of each of the middle plan data of the virtual model with a reference length.

In this case, the reference length is the center of the virtual model is located at the center point of the central axis, the center of the virtual model, the center point of the central axis and the center of the X-ray source is located in a straight line, the center of the virtual model and the center axis of the When the center points coincide, this is the length calculated theoretically for the measurement length of the middle plan data that can be obtained.

Subsequently, the center points of the central axis of the X-ray imaging apparatus are corrected using the magnification ratios.

In this case, there are two methods of correcting the magnification caused by the movement of the central axis in the X-ray Citigraphy apparatus. First, the obtained magnification ratios are stored, and the detector obtained by the X-ray Citigraphy apparatus is obtained. There is a method of correcting projection data of an object included in the detector plan data, accurately reconstructing plan data, and reconstructing the corrected detector plan data to obtain a three-dimensional image. In consideration of the above, there is a method of rearranging positions of sources and detectors of the X-ray Citigraphy apparatus so that the center point of the central axis of the X-ray Citigraphy apparatus is not moved or moved.

That is, in the first method, the center point of the central axis of the X-ray imaging apparatus does not coincide with the center point of the object or is moved during scanning, and is corrected in software to obtain a clear 3D image by correcting the 3D image reconstruction. The second method is to align the X-ray source and the detector of the X-ray Citigraphy device, and to hardware-correct that the center point of the center axis of the X-ray Citigraphy device does not coincide with the center point of the object or is moved during scanning. That's how.

1A and 1B are diagrams illustrating a method of obtaining detector plan data of an X-ray Citigraphy apparatus.

Referring to FIGS. 1A and 1B, the X-ray city photographing apparatus 100 includes an X-ray source 110 and a detector 120.

The X-ray imaging apparatus 100 irradiates the object 130 having the object point 135 with the X-ray generated by the X-ray source 110 and is received by the detector 120. In this case, the X-rays are irradiated with a cone beam.

An image 140 of the object 130 is formed in the detector 120. When the detector 120 reads the image 140, the object 120 becomes projection data of the object. In this case, the projection data of the object includes the object point data generated by the object point 135. In this case, the data read from the detector 120 is detector plan data.

Therefore, the data read from the detector 120 may be divided into detector plan data, which is data of the entire detector 120, and projection data, which is a part of the detector plan data.

In order to obtain the 3D image, the X-ray imaging apparatus 100 rotates the circular trajectory up to 360 degrees while the X-ray source 110 and the detector 120 move by the same angle about the central axis 150. A plurality of detector plan data, which is data that is a two-dimensional image, is acquired by a scan.

The plurality of detector plan data including the target point 135 data may be reconstructed into a 3D image.

The 3D image reconstruction has various techniques, but in the present invention, the filter back-projection algorithm proposed by Feldkamp, Davis, and Kress (FDK) was used.

First, preprocessing and filtering are performed. The detector plan data acquired by the X-ray imaging apparatus 100 are loaded one by one.

Then, Jacobian weighting processing is performed on the detector plan data.

In this case, the Jacobian weighting process is a process for correcting a change in a value between polar coordinates and rectangular coordinates during a mathematical process.

The detector plan data is filtered using a filter.

At this time, the filtering becomes the core of the filter bag-projection algorithm technique. When the filtering is not performed, the shape of the object part derived as a result of the reconstruction appears blurred and unclear. To solve this problem, filter processing is performed on each projection data so that the shape of the original object can be approximated. The filtering process is performed between the Jacobian weighting process and the back-projection process to be described later.

Then back-projection is performed. The weighting factor of the target point data of the detector plan data according to the distance is calculated.

Then, interpolation is performed using the detector plan data to determine the value of each target point of the three-dimensional volume. And the value of the target point is stored.

The reconstructed slice is then stored.

Each of the detector plan data is sequentially processed through the preprocessing, filtering, back-projection processing, and reconstruction section storage, and then the reconstruction result is stored in the section.

At this time, the X-ray source 110, the center point 155 of the central axis and the center point 125 of the detector are located on one center connecting line 160 to maintain a straight line to the X-ray around the center point 155 of the central axis The source 110 and the center point 125 of the detector are rotated.

As described with reference to FIGS. 1A and 1B, the X-ray source 110, the center point 155 of the central axis, and the center point 125 of the detector may be used to acquire an accurate three-dimensional image in the X-ray imaging apparatus 100. The center connection line 160, which is a straight line, should be positioned. Due to a mechanical defect or vibration during driving of the X-ray imaging apparatus 100, the center point 155 of the central axis may not move while being scanned.

2A and 2B are diagrams showing that the center point of the central axis is moved. In this case, the movement of the center point of the central axis is moved in the depth direction, wherein FIG. 2A shows that the center point of the center axis is located behind the center point of the object, and FIG. 2B shows that the center point of the center axis is the center point of the object. It shows the one located earlier.

2A and 2B, the center points of the X-ray source 110 and the center axis of the detector 120 of the X-ray imaging apparatus 100 coincide with the center point C of the object O. The image D of one object O and the center points C1 and C2 of the central axis are located behind or in front of the center point C of the object O, and the image of the object O photographing the object O ( Comparing D1 and D2, the center point of the central axis is compared to the image D of the object O photographed by being positioned at the center point C of the object O as shown in the drawings. When C1) is located behind the center point C of the object, a phenomenon in which the image D1 of the object O looks larger occurs, and when the center point C2 of the central axis is located in front of the center point C of the object. In this case, a phenomenon in which the image D2 of the object O looks smaller occurs, which is called a magnification problem.

The magnification problem is that the object (O) is fixed, but the center point of the central axis of the X-ray source 110 and the detector 120 of the X-ray imaging apparatus 100 is located and fixed at the center point of the object (O). It is a phenomenon that occurs while moving forward and backward without being.

When the above-described three-dimensional image reconstruction is performed on the detector plan data obtained in the state where the magnification problem occurs, a three-dimensional image in which blurring phenomenon as if the image overlaps is generated is obtained.

3, 4A, 4B, and 5A through 5C are diagrams illustrating a method of correcting a central axis movement of an X-ray Citigraphy apparatus according to an exemplary embodiment.

Referring to FIG. 3, FIG. 3 illustrates an overall X-ray imaging apparatus for illustrating a method of correcting a central axis movement of an X-ray imaging apparatus according to an embodiment of the present invention.

As described above, the X-ray imaging apparatus 200 includes an X-ray source 210 and a detector 220. At this time, the virtual point of the sphere (sphere) or cylinder (cylinder) so that the cross section can be drawn to the circle (circle) in the center point 255 of the central axis 250 that the X-ray source 210 and the detector 220 rotates Model 230 is located.

At this time, the center point of the virtual model 230 is positioned to exactly match the center point 255 of the central axis. Whether the center point of the virtual model 230 coincides with the center point 255 of the central axis may be used as it is. That is, this may be confirmed by using a laser beam provided in the X-ray city photographing apparatus 200. In addition, the center point of the X-ray source 210, the center point 225 of the detector, and the center point of the virtual model 230 are positioned on the center connection line 260 which is a straight line.

In the present invention, the virtual model 230 may be a cylindrical virtual model, the virtual model 230 has a diameter of 10mm or more, the height is preferably 50mm or more, and the material is aluminum Or teflon or the like.

The virtual model 230 is scanned by the X-ray Citigraphic apparatus equipped with the virtual model 230.

In this case, the X-ray source 210 and the detector 220 of the X-ray imaging apparatus are scanned at a predetermined angle to obtain a plurality of detector plan data.

After obtaining the plurality of detector plan data, the middle plan data 245 of the virtual model 230 is extracted from each detector plan data.

Referring to FIGS. 4A and 4B, a method of extracting middle plan data 245 of the virtual model 230 from the detector plan data is illustrated, and FIG. 4A illustrates one detector plan data. Assuming that the detector 220 is composed of 256 * 256 pixels and the detector plan data has obtained such data, the 128th pixel line, which is a line located in the middle of the rows and columns of the detector plan data, is the row. It becomes the centerline in overheating. In this case, the point where the center line of the row and the column meet is the center point 225 of the detector described above and the center point of the detector plan data. When the image generated by the virtual model 230 is read, the projection data 240 of the virtual model 230 is obtained.

In this case, in the middle plan data 245 of the projection data 240 of the virtual model 230, as illustrated in FIG. 4B, the pixel line of the 128th row, which is the middle of the column, is extracted from the detector plan data. If only the pixels having data of the virtual model 230 are selected among the data, the middle plan data 245 of the virtual model 230 can be extracted.

In this case, the middle plan data 245 of the virtual model 230 illustrated in FIG. 4B may appear to be located at the same position as the center point 225 of the detector, that is, the center point of the detector plan data. It is unknown whether the center point 255 of the central axis of the photographing apparatus coincides with the center point of the virtual model 230. That is, with only the middle plan data 245 of the virtual model 230, the center point 255 of the central axis is vertical, that is, the center point of the virtual model 230 coincides with the center point of the virtual model 230, or It is not known whether it is located behind or in front of the center point of the virtual model 230.

Referring to FIGS. 5A through 5C, FIGS. 5A through 5C illustrate middle plan data 245 of the virtual model 230.

At this time, the middle plan data 245 shown in FIG. 5A is the middle plan data 245 having a reference length K calculated theoretically. As described with reference to FIG. 3, the reference length K confirms whether the center point of the virtual model 230 and the center point of the central axis of the X-ray source 210 and the detector 220 are exactly the same, and then the middle plan. The data 245 may be acquired, and the length of the middle plan data 245 (exactly, the number of pixels of the middle plan data 245 may be calculated and obtained) may be obtained. It can be calculated and obtained mathematically using geometry information such as the distance from the X-ray source of the device to the center point of the detector, the distance from the X-ray source to the center point of the virtual model, and the diameter of the virtual model.

On the other hand, the middle plan data 245 shown in FIGS. 5B and 5C is the center point of the central axis of the X-ray source 210 and the detector 220, respectively, behind or in front of the center point of the virtual model 230 (Fig. 2A and FIG. 2B), which are actually obtained middle plan data 245.

In this case, the magnification ratio R, which is an index indicating a deviation from the center point of the central axis and the center point of the virtual model 230 of the middle plan data 245 illustrated in FIGS. 5B and 5C, may be calculated.

That is, it can be obtained by dividing the actual lengths L1 and L2 of the middle plan data 245 actually obtained by the reference length K when the center point of the central axis coincides with the center point of the virtual model 230. , If it is expressed by an equation, is equal to Equation 1 below.

(Where R is the magnification ratio value, K is the reference length, and L is the actual length.)

For example, a virtual model 230 having a 52-pixel number length of the reference length K may be inserted into an X-ray Citigraphy apparatus and scanned to obtain a plurality of detector plan data, and the respective detector plan data may be obtained from the detector plan data. After obtaining the middle plan data, the actual length L of the middle plan data is obtained.

At this time, if the actual length L of any one of the middle plan data has a 56 pixel number length, the magnification ratio of the middle plan data is about 1.077. As shown in FIG. 2A, the center point C1 of the central axis of the X-ray source 110 and the detector 120 is located behind the center point C of the object O, that is, the center point of the virtual model 230. Can be.

The magnification ratios of the remaining middle plan data may be obtained in the same manner as described above.

In this case, the more magnification ratios obtained from the middle plan data and the larger the deviation ratio are from 1, the larger the deviation ratio is, so that the plurality of detector plan data are reconstructed into a 3D image. Ring phenomenon is severely generated.

In order to prevent the blurring, the movement of the center point of the central axis of the X-ray imaging apparatus may be corrected using the obtained magnification ratio. There are two methods.

The first method is a software method of correcting detector plan data obtained by the X-ray Citigraphy apparatus using the magnification ratios.

That is, each magnification ratio is calculated and stored from the detector plan data, the detector plan data are corrected using the magnification ratios, and then reconstructed to obtain a 3D image. In this case, the time for correcting the detector plan data by using the magnification ratios is when the back-projection is performed. The back-projection is performed by replacing the position of the detector plan data to be read with the position changed by the product of the magnification ratios. This is the same as performing the back-projection after all the pixel sizes of the detector plan data are equal to the reference length.

The stored magnification ratios can then be used to continuously correct the detector plan data using the X-ray Citigraphy apparatus before it is corrected by a hardware correction method described below or once again by software. .

Second, there is a hardware method of correcting the movement of the center point of the central axis of the X-ray imaging apparatus using the magnification ratios.

That is, after calculating and obtaining the magnification ratios from the detector plan data obtained by the X-ray imaging apparatus, the X-ray source 210 and the detector 220 of the X-ray imaging apparatus are used by using the magnification ratios. By rearranging the positions of the center axes of the central axis is located in the center of the virtual model (230).

Therefore, the method for correcting the magnification caused by the movement of the central axis in the X-ray Citigraphy apparatus of the present invention obtains detector plan data and projection data by performing a CIT imaging of a virtual model, and extracts middle plan data from the detector plan data. And obtaining a magnification ratio by comparing the theoretical length theoretically calculated with the actual length of the middle plan data, and then presenting a method of correcting the center point of the zoom axis of the X-ray imaging apparatus using the magnification ratio. X-ray source of the device and the center point of the center axis of the detector is not located in the correct position, or blurring caused by the movement can be removed to eliminate the cause of the three-dimensional image.

While the present invention has been particularly shown and described with reference to exemplary embodiments thereof, it is clearly understood that the same is by way of illustration and example only and is not to be taken by way of limitation, Various changes and modifications will be possible.

1A and 1B are diagrams illustrating a method of obtaining detector plan data of an X-ray Citigraphy apparatus.

2A and 2B are diagrams showing that the center point of the central axis is moved. In this case, the movement of the center point of the central axis is moved in the depth direction, wherein FIG. 2A shows that the center point of the center axis is located behind the center point of the object, and FIG. 2B shows that the center point of the center axis is the center point of the object. It shows the one located earlier.

3, 4A, 4B, and 5A through 5C are diagrams illustrating a method of correcting a central axis movement of an X-ray Citigraphy apparatus according to an exemplary embodiment.

<Description of the symbols for the main parts of the drawings>

100,200: X-Ray Citigraphy Device 110,210: X-Ray Source

120,220 Detector 130 Object

230: virtual model

Claims (9)

Preparing a virtual model; Positioning a center point of the virtual model on the X-ray city photographing device by coinciding with a center point of a central axis of the X-ray city photographing device; Acquiring a plurality of detector plan data by scanning the virtual model with the X-ray imaging apparatus; Extracting middle plan data of the virtual model in a one-to-one correspondence with each of the detector plan data to obtain a measurement length of the middle plan data; Obtaining a magnification ratio of each of the middle plan data of the virtual model by comparing the measured length and the reference length of each of the middle plan data of the virtual model corresponding one-to-one with each of the detector plan data; And And correcting the center point of the central axis of the x-ray city photographing apparatus by using the magnification ratios. The method of claim 1, The reference length may be obtained when a center point of the virtual model is located at the center point of the central axis, a center point of the virtual model, a center point of the central axis, and a center point of the X-ray source of the X-ray imaging apparatus are located in a straight line. A method for correcting the magnification caused by the movement of a central axis in an X-ray Citigraphic apparatus, characterized in that the theoretical length of the measured length of the middle plan data. The method of claim 1, Correcting the center point of the central axis of the X-ray Citigraphic apparatus using the magnification ratios In the process of reconstructing the detector plan data in three dimensions by using a filter back-projection algorithm, correcting detector plan data by multiplying each of the detector plan data by respective magnification ratios immediately before a back projection process. A method of correcting magnification attributable to the movement of a central axis in an x-ray imaging apparatus. The method of claim 3, wherein And the filter back-projection algorithm comprises a preprocessing process, a filtering process, a back-projection process, and a reconstructed cross-section storage process. The method of claim 1, Correcting the center point of the central axis of the X-ray Citigraphic apparatus using the magnification ratios And rearranging the positions of the X-ray source and the detector of the X-ray Citigraphic apparatus with the obtained magnification ratios. 6. The method according to any one of claims 1 to 5, Extracting middle plan data of the virtual model corresponding to each of the detector plan data one-to-one may extract data of a row located in the middle of all rows from each of the detector plan data, and extract the data of the extracted intermediate row. And extracting data of the virtual model from among them. The method of claim 6, The virtual model is a method for correcting the magnification due to the movement of the central axis in the x-ray Citigraphic apparatus, characterized in that the cross section is caused. The method of claim 7, wherein And the virtual model is a spherical or cylindrical shape. The method of claim 8, Aligning the center point of the virtual model with the center point of the central axis of the X-ray Citigraphy apparatus and positioning it on the X-ray Citigraphy apparatus Magnification caused by the movement of the central axis in the x-ray city imaging device, wherein the center point of the virtual model is matched with the center point of the central axis of the x-ray city imaging device by using a laser beam provided in the x-ray city imaging device. How to calibrate.

Images