KR101704390B1 - Spiral ct systems and reconstruction methods - Google Patents

Abstract

The present invention relates to a spiral CT system and a reconstruction method thereof. In some embodiments, a method is provided for compensating for data loss due to large pitch using weighted processing of the complementary projection data of the projection data obtained using the spiral CT system. After the data is supplemented, the projection data is rearranged into parallel beam data having a cone angle, the cone weight cosine weighting process and one-dimensional filtering are performed, and the parallel beam backprojection is performed later to obtain a reconstructed image. In some embodiments, by using the above-described method, it is possible to improve the passing rate of the hydrate by improving the belt speed by at least 1 time under the premise that the area of the conventional probe and the speed of the slip ring are not changed, It is possible to keep the quality unchanged.

Description

[0001] SPIRAL CT SYSTEMS AND RECONSTRUCTION METHODS [0002]

The present invention relates to radiation imaging, and more particularly, to a large pitch CT system and its reconstruction method.

The CT-type baggage security system is already widely used in public places such as airports and terminals as an important means of detecting explosions. The CT detection system obtains projection data of the hydrate through scanning, acquires a tomographic image by a reconstruction algorithm, identifies the explosion through an identification algorithm, and sends an alarm. Such a scanning method generally uses a method in which an X-ray light source and a slip ring equipped with a probe are rotated and an object advances through a belt, so that the scanning orbit belongs to a spiral orbit. In order to satisfy the requirement for the passage rate of the baggage at the security search location, the advancement speed of the belt must always be a certain limit value, for example 0.3 m / s to 0.5 m / s. In order to acquire the complete projection data necessary to reconstruct the image, the slip ring speed must be increased or the row number of the probe must be increased. Considering mechanical strength and stability of the slip ring, the number of revolutions can not be increased to infinity. On the other hand, considering the limit value of the cone angle of the X-ray source and the unit price of the probe hardware, the number of rows of the probe can not be increased infinitely. Taking into consideration the various factors mentioned above, increasing the scanning pitch is the most effective way to improve the passage rate of the hydrate.

Recently, in the field of spiral CT reconstruction algorithms, for example, Katsevich algorithm, PI algorithm, CB-FBP algorithm and so on have achieved significant results. However, since all of these algorithms must satisfy a certain pitch, projection data is lost when the pitch and the cone angle are large, and the error of the reconstruction result is also large, and further, artifacts are caused . In order to meet the demand for image quality, the pitch factor is generally not greater than 1.5.

For secure CT systems, the realization of reconfiguration algorithms is also an important evaluation factor. In a reconstruction algorithm, a non-one-dimensional shift-invariant form of filtering, a cone beam back projection with a distance weighting factor, a large number of nonlinear equation calculations , The use of a relatively large backprojection angular range will all degrade the efficiency of the reconstruction algorithm, so preferably the use of such an algorithm should be avoided. Therefore, each of the above-described reconstruction algorithms includes a part that affects algorithm efficiency and is mainly used when the projection data is perfect or excessive.

As noted above, current prior art techniques can not be used directly in large pitch CT systems because they do not take into account the problem of data loss under large pitches.

The present invention provides a spiral CT system and its reconstruction method that can satisfy image reconstruction when a pitch is large based on one or a plurality of problems existing in the prior art.

According to an aspect of the present invention, there is provided a method of calculating a minimum row quantity of a probe required to cover a Tam window by a pitch of a cone-beam spiral CT system and an interval between rows of probes of plural rows ; Complementing the lost projection data by performing weighting processing on the complementary projection data when the number of rows of the probe of the cone beam helical CT system is smaller than the minimum number of rows; Rearranging the complementary projection data into parallel beam data having a cone angle; Performing one-dimensional filtering according to a row direction of a virtual probe defined when rearranging parallel beam data after rearranging into parallel beam data after proceeding with cone angle cosine weighting processing on rearranged parallel beam data; And performing a parallel beam backprojection with a cone angle that is not weighted for the filtered data to obtain a reconstructed image; Beam CT reconstruction method.

According to another aspect of the present invention there is provided an apparatus for calculating a minimum row count of a probe required to cover a Tam window by the pitch of a cone beam helical CT system and the spacing between rows of probes of a plurality of rows; An apparatus for compensating lost projection data by progressing a weighting process on complementary projection data when the number of rows of probes of the cone beam helical CT system is smaller than the minimum number of rows; A device for rearranging the complementary projection data into parallel beam data having a cone angle; A device for performing one-dimensional filtering according to a row direction of a virtual probe defined when rearranging parallel beam data after rearranging into parallel beam data after proceeding with a cone angle cosine weighting process on rearranged parallel beam data; And a device for obtaining a reconstructed image by advancing a parallel beam reverse projection having a cone angle without weighting processing on the filtered data.

In some embodiments, by improving the belt speed by at least 1 time under the premise that the area of the conventional probe and the speed of the slip ring are not changed by using the above-described technique, the passing rate of the hydrate can be improved and the quality of the reconstructed image Can be kept unchanged.

BRIEF DESCRIPTION OF THE DRAWINGS For a better understanding of the present invention, the present invention will be described in more detail with reference to the following drawings.
FIG. 1 is a view showing a helical orbital scan in a CT system according to an embodiment of the present invention, wherein the gray area of the probe indicates a full window.
Figure 2 shows a definition of a pair of conjugate rays in a sector scan according to an embodiment of the present invention.
3 is a structural view of a CT facility according to an embodiment of the present invention.
4 is a structural block diagram of a computer data processor as shown in FIG.
5 is a structural block diagram of a controller according to an embodiment of the present invention.
6 is a flowchart of a reconfiguration method according to an embodiment of the present invention.
7 is a diagram showing a conjugate projection interpolation value when scanning according to the helical trajectory according to the embodiment of the present invention.
Figure 8 shows a diagram rearranged with a parallel beam and a virtual probe having cone angles.
9 shows a reconstruction result when the pitch is 3.9 cm in one specific embodiment of the present invention.
Fig. 10 shows a reconstruction result when the pitch is 6 cm in one specific embodiment of the present invention.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS It should be noted that the embodiments described herein are illustrative and not intended to limit the present invention. In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention, but it will be understood by those skilled in the art that the present invention may be practiced without these specific details. In other instances, well-known structures, materials, or methods are not described in detail so as not to obscure the spirit of the present invention.

Reference in the specification to "one embodiment", "an embodiment", "an example", or "an instance" means that a particular feature, structure, or characteristic described in connection with the embodiment Quot; is included in at least one embodiment of < / RTI > Accordingly, the appearances of the phrases " in one embodiment, "" an embodiment," " an example, "or" an example, " In addition, certain features, structures, or characteristics may be used in combination with one or more embodiments or examples in any suitable combination and / or subcombination. Also, those skilled in the art will appreciate that the term "and / or" as used herein includes any and all combinations of one or more related features.

In contrast, while the prior art does not satisfy the reconstruction requirement when the pitch is large, especially when the pitch coefficient is greater than 1.5, the embodiment of the present invention uses weighted processing of the complementary projection data of the projection data obtained using the spiral CT system To compensate for data loss due to large pitches. After all the data has been supplemented, the projection data is rearranged into parallel beam data having a cone angle, the cone weight cosine weighting process and the one-dimensional filtering process are performed, and finally the parallel beam backward projection is performed to reconstruct the reconstructed image . In some embodiments, by using the above-described method to improve the belt rate by more than one-fold under the premise that the conventional probe area and the slip ring speed are not changed, the passing rate of the hydrate is improved and the quality of the reconstructed image is improved It can be kept unchanged.

, X레이 광원(110)으로부터 탐측기(120)까지의 거리를 D라 한다. 그리고 탐측기(120)의 행 수량을

, 행 사이의 간격을

라고 한다. 그러면 피치계수를 아래와 같이 정의내릴 수 있다.BRIEF DESCRIPTION OF THE DRAWINGS FIG. 1 is a diagram showing a helical orbital scan of a CT system according to an embodiment of the present invention; FIG. 1, the rotation radius of the X-ray light source 110 is denoted by R, the distance (referred to as a pitch) by which the belt advances when the slip ring rotates once is h, cylindrical equiangular detectors 120 ) Angle of view of the sector beam

, And the distance from the X-ray light source 110 to the probe 120 is D And the row quantity of the probe 120

, The spacing between rows

. The pitch coefficient can then be defined as:

The rotation of the slip ring and the translation of the belt form a helical orbit 131 that moves relatively. In a coordinate system in which an object is scanned, the motion trajectory of the X-ray light source can be expressed by the following equation.

과 피치 h 등 파라미터가 정해졌을 때, Tam window를 피복하는데 필요한 최소 탐측기의 행 수량을 계산할 수 있다.In the helical scan trajectory, the projection data necessary for correct reconstruction is obtained by projecting data on the projection data on the probe, i.e., the projection data, which is covered by the projection on the probe, i.e., the Tam-Danielsson window As shown in Fig. From this,

And pitch h, parameters can be calculated to calculate the minimum number of rows of probes required to cover the Tam window.

From the above equation (3), it is possible to derive the maximum pitch and the maximum pitch coefficient allowed by the system in the CT system in which the row number of the probe, the spacing between the rows, and the sector angle are fixed.

일 때, 최대 피치계수는 1.3이다.For example,

, The maximum pitch coefficient is 1.3.

, X레이 광원(210)으로부터 탐측기(220) 까지의 거리를 D 라 한다. 2차원 부채꼴 주사에서, 도 2에 도시된 바와 같이, 일 직선을 지나는 두갈래의 방사선, 즉 방사원의 위치(210)으로부터 위치(210')까지의 방사선과 위치(210')로부터 위치(210)까지의 방사선을 한쌍의 공액 광선

이라 한다. 콘 빔 나선형 주사(cone beam spiral scanning)에서는 엄밀한 한쌍의 공액 광선이 존재하지 않으나, 유사한 관계가 있는 상호보완투영을 정의할 수 있다.2 shows a definition of a pair of conjugate rays in a sector scan according to an embodiment of the present invention. 2, the rotation radius of the X-ray light source 210 is R, and the distance (referred to as a pitch) in which the belt advances when the slip ring rotates once is defined as an angle of view of the sector beam of the mirror surface of

, And a distance from the X-ray light source 210 to the probe 220 is D In a two-dimensional sector scan, as shown in FIG. 2, a bifurcated radiation passing through a straight line, from the position 210 to the position 210 'of the radiation source and from the position 210'Lt; RTI ID = 0.0 > rays < / RTI &

Quot; In cone beam spiral scanning, there is no rigid pair of conjugate rays, but a mutually complementary projection with a similar relationship can be defined.

Thus, in the back projection, the spiral reconstruction algorithm can either use the complementary projection data to advance the excess weighting process or to supplement the lost data.

3 is a structural view of a CT facility according to an embodiment of the present invention. 3, the CT apparatus according to the present embodiment includes a frame 20, a mounting unit 40, a controller 50, a computer data processor 60, and the like. The frame 20 includes a radiation source 10 for emitting detection X-rays such as an X-ray machine; And a probe and collecting device (30). The mounting unit 40 mounts the inspected hydrate 70 and passes through the scanning area between the radiation source 10 of the frame 20 and the probe and collecting device 30. [ At the same time, the frame 20 rotates about the advancing direction of the inspected hydrate 70, allowing the cone beam emitted from the radiation source 10 to pass through the inspected hydrate 70, 70).

The probing and collecting device 30 may be, for example, a probe and a data collecting device having an integral module structure. For example, a plurality of rows of probes are used to detect an X-ray transmitted through an object to obtain an analog signal and to convert it into a digital signal, and output projection data of the inspected object 70 by X-ray. The controller 50 is for controlling the synchronous operation of each part of the overall system. The computer data processor 60 is for processing the data collected by the data collector, processing and reconstructing the data, and outputting the results.

3, the radiation source 10 can be placed on one side of the subject, and the probe and collecting device 30 is disposed on the other side of the inspected hydrate 70, and the penetration of the inspected hydrate 70 And a data collector for obtaining data and / or polygonal projection data. The data collector includes a data amplification and shaping circuit, which can operate in (current) integration or pulse (count) mode. The data output cable of the probe and collecting device 30 is connected to the controller 50 and the computer data processor 60 and stores the data collected by the trigger command in the computer data processor 60.

4 is a structural block diagram of a computer data processor as shown in FIG. As shown in Fig. 4, the data collected in the data collector is stored in the memory 61 via the interface unit 68 and the bus (BUS) A read only memory (ROM) 62 stores setting information and a program of a computer data processor. The random access memory (RAM) 63 is for temporarily storing various data of the operation process of the processor 66. The memory 61 also stores a computer program for data processing. The internal bus 64 includes the memory 61, the read only memory 62, the random access memory 63, the input device 65, the processor 66, the display device 67, and the interface unit 68, Lt; / RTI >

When the user inputs an operation command via an input device 65 such as a keyboard, a mouse, etc., the instruction code of the computer program instructs the processor 66 to execute a predetermined data reconstruction algorithm, For example, on a display device 67 such as an LCD display, or outputs a processing result in hard copy format such as direct printing.

5 is a structural block diagram of a controller according to an embodiment of the present invention. 5, the controller 50 includes a control unit 51 for controlling the radiation source 10, the mounting unit 40, and the probe and collecting apparatus 30 in accordance with an instruction from the computer 60, ; A trigger signal generating unit (52) for generating a trigger command for triggering the operation of the radiation source (10), the probe and collecting apparatus (30) and the mount unit (40) under the control of the control unit; A first drive facility 53 for driving the loading unit 40 to transport the inspected hydrate 70 in accordance with a trigger command generated by the trigger signal generating unit 52 under the control of the control unit 51; And a second drive facility (54) for causing the trigger signal generation unit (52) to rotate the frame (20) in accordance with a trigger command generated under the control of the control unit (51). The projection data obtained by the probe and collecting device 30 may be stored in the computer 60 and applied to CT tomographic image reconstruction to obtain tomographic image data of the inspected hydrate 70. According to another embodiment, the CT imaging system described above may be a dual-energy CT system. That is, the X-ray radiation source 10 of the frame 20 can emit high-energy rays and low-energy rays, and the probe and collecting device 30 can emit different energy The computer data processor 60 proceeds to dual energy CT reconstruction to obtain the equivalent atomic number and equivalent electron density data of each single layer of the inspected hydrate 70. [

을 계산한다.6 is a flowchart of a reconfiguration method according to an embodiment of the present invention. As shown in Fig. 6, in step S61, the minimum row quantity of the probe required to cover the Tam window is calculated according to the pitch of the cone beam helical CT system and the spacing between the rows of the probes of plural rows. For example, the minimum number of probe rows required to cover the Tam window by parameters such as pitch, spacing between rows

.

보다 작으면, 상호보완투영을 통하여 손실된 데이터를 보완한다. 도 7에 도시된 바와 같이, 두가지 경우로 나뉜다.In step S62, when the number of rows of the probes of the cone beam helical CT system is smaller than the minimum number of rows of the probes, the weighted processing is performed on the complemented projection data to compensate the lost projection data. For example, if the actual

, The complementary projection compensates for the missing data. As shown in Fig. 7, there are two cases.

일때,One)

when,

일때,2) Similarly,

when,

Where s is the coordinate in the direction of the probe (cone angle), α is the coordinate in the direction of the probe column (sector angle), λ is the projection angle, smin is the minimum value in the direction of the probe row (cone angle), smax is the probe row R is the radius of rotation of the X-ray source, h is the distance that the belt advances when the slip ring rotates once, that is, pitch, and D is the distance from the X-ray source to the column surface probe.

를 원추각을 가지는 평행 빔 데이터

로 재배열한다. 도 8에 도시된 바와 같이, 재배열된 가상 탐측기는, 회전 중심을 지나는 z축에 평행되는 구형 영역이며, 그 폭이 실제 기둥면 탐측기의 부채꼴 각이 덮는 시야 직경

와 같다.In step S63, all the complementary projection data

Parallel beam data having a cone angle

. As shown in Fig. 8, the rearranged virtual probe is a spherical region parallel to the z-axis passing through the center of rotation and the width of which is the field diameter of the actual circular-

.

In step S64, the data after the rearrangement is weighted first by the cone angle cosine, and then the one-dimensional filtering is performed. The filtering direction follows the row direction of the virtual probe. Since the filtering process is a one-dimensional shift-invariant form, filtering can be performed using fast Fourier transform. The filter kernel selects the R-L convolution kernel in the parallel beam filtering reverse projection algorithm.

In step S65, the parallel beam reverse projection having the cone angle without weighting processing proceeds to obtain a reconstructed image.

According to the method of the present embodiment, it is possible to improve the belt speed by one time in the right and left direction without changing conditions such as the area of the conventional probe and the number of revolutions of the slip ring, and further improve the passing rate of the baggage, do. On the other hand, by using the method of the present application, it is possible to reduce the system cost by appropriately reducing the number of rows of the probe by changing the conventional CT system design under a specific use environment.

For example, when the rotation radius of the X-ray source is 50 cm, the distance from the radiation source to the probe is 80 cm, the sector angle of the probe is 60, the number of rows of the probe is 32, , The maximum pitch coefficient should be 1.3 and the corresponding maximum pitch should be 3.9 cm. Using the above-described method of the present application, the maximum pitch can be increased to 6 cm, and the pitch coefficient is 2.

9 shows a reconstruction result when the pitch is 3.9 cm in one specific embodiment of the present invention. Fig. 10 shows a reconstruction result when the pitch is 6 cm in one specific embodiment of the present invention. As shown, the present invention not only improves scanning speed, but also maintains the quality of the reconstructed image.

According to some embodiments, the belt speed can be further improved by using a conventional probe, a slip ring, etc., and the passing rate can be improved, so that the present invention can be applied to security fields such as detonation detection and drug detection. In addition, the method of compensating lost data using the complementary projection described above can be used when the projection data is not perfect. Even if the pitch coefficient is larger than 1.5, the quality of the reconstructed image can still be secured. At the same time, the reconstruction algorithm provided by this embodiment has a shift-invariant filtering form, using a parallel beam back projection method without a distance weighting factor, using a minimum of 180 degrees back projection angular range , It is unnecessary to calculate the nonlinear equations. Therefore, it is much simpler than the conventional technique, and the security real time requirement can be satisfied.

In the foregoing detailed description, various embodiments of the reconstruction method and the spiral CT system have been described by way of structural diagrams, flowcharts and / or examples. It will be appreciated by those skilled in the art that such a structure, a flowchart, and / or an example may include one or more functions and / or operations, (Firmware), or any combination thereof. In one embodiment, some of the topics described in the embodiments of the present invention may be implemented through application specific integrated circuits (ASICs), field programmable gate arrays (FPGAs), digital signal processors (DSPs) Those skilled in the art will recognize that some aspects of the embodiments disclosed herein may be embodied in whole or in parts in their entirely equivalent circuitry and that one or more computer programs running on one or more computers (e.g., One or a plurality of programs executed in one or a plurality of microprocesses) and executed by one or a plurality of processes (for example, one or a plurality of programs executed in one or a plurality of microprocesses) ), May be implemented in firmware, or may be It will also be appreciated by those skilled in the art that the present disclosure has the ability to write to and / or write firmware and / or firmware to circuit design and / or software. Those skilled in the art will also appreciate that, It will be appreciated that the principles of the present disclosure may be distributed in various types of program products and that all exemplary embodiments of the present disclosure may be applied, regardless of the specific type of signal bearing medium actually being deployed. Examples of media include storage media such as floppy diskettes, hard drives, compact discs (CD), digital versatile discs (DVD), digital magnetic tape, computer memory, and the like, as well as digital and / or analog communication media Optical cable, waveguide, wired communication link, wireless communication link, etc.) It is, but is not limited to this.

Although the invention has been described with reference to certain exemplary embodiments, it is to be understood that the terminology used herein is for the purpose of description and is not intended to be limiting of the invention. The present invention may be embodied in various ways without departing from the spirit and scope of the invention. It is therefore to be understood that the above-described embodiments are not limited by any of the above-described specific portions, but can be construed broadly within the spirit and scope defined by the claims. It is therefore to be understood that the appended claims are intended to cover all modifications and equivalents falling within the scope of the appended claims and their equivalents.

Claims (5)

Cone Beam Spiral CT reconstruction with cone beam helical CT system,
Calculating a minimum row count of a probe required to cover the Tam window by the pitch of the cone beam helical CT system and the spacing between the rows of the probes of the plurality of rows;
Complementing the lost projection data by performing weighting processing on the complementary projection data when the number of rows of the probe of the cone beam helical CT system is smaller than the minimum number of rows;
Rearranging the complementary projection data into parallel beam data having a cone angle;
Performing one-dimensional filtering according to a row direction of a virtual probe defined when rearranging parallel beam data after rearranging into parallel beam data after proceeding with cone angle cosine weighting processing on rearranged parallel beam data; And
Proceeding a parallel beam backprojection having a cone angle that is not weighted for the filtered data to obtain a reconstructed image; Containing
Cone beam helical CT reconstruction method. The method according to claim 1,
The step of compensating lost projection data by progressing a weighting process on the mutually complementary projection data includes:

when,

when,

The method comprising:
Where s is the coordinate of the probe row (cone angle), a is the probe column (sector angle) direction coordinate, lambda is the projection angle, smin is the minimum value of the direction of the probe row (cone angle) Where R is the radius of rotation of the X-ray source, h is the distance or pitch at which the belt advances when the slip ring rotates once, D is the distance from the X-ray source to the column surface probe And P (.) Represents the projection data. The method according to claim 1,
Wherein the one-dimensional filtering uses an RL convolution kernel. An apparatus for calculating a minimum row count of a probe required to cover a Tam window by a pitch of a cone beam helical CT system and a spacing between rows of probes of a plurality of rows;
An apparatus for compensating lost projection data by progressing a weighting process on complementary projection data when the number of rows of probes of the cone beam helical CT system is smaller than the minimum number of rows;
A device for rearranging the complementary projection data into parallel beam data having a cone angle;
A device for performing one-dimensional filtering according to a row direction of a virtual probe defined when rearranging parallel beam data after rearranging into parallel beam data after proceeding with a cone angle cosine weighting process on rearranged parallel beam data; And
And a device for obtaining a reconstructed image by proceeding with a parallel beam backprojection having a cone angle that is not weighted for the filtered data
Cone beam spiral CT system. 5. The method of claim 4,
An apparatus for compensating lost projection data by performing a weighting process on mutually complementary projection data,

when,

when,

Lt; RTI ID = 0.0 >
Where s is the coordinate of the probe row (cone angle), a is the probe column (sector angle) direction coordinate, lambda is the projection angle, smin is the minimum value of the direction of the probe row (cone angle) Where R is the radius of rotation of the X-ray source, h is the distance or pitch at which the belt advances when the slip ring rotates once, D is the distance from the X-ray source to the column surface probe Distance, and P (.) Represents the projection data.

Images