WO2004034908A1 - X線断層撮影装置 - Google Patents
X線断層撮影装置Info
- Publication number
- WO2004034908A1 WO2004034908A1 PCT/JP2003/010971 JP0310971W WO2004034908A1 WO 2004034908 A1 WO2004034908 A1 WO 2004034908A1 JP 0310971 W JP0310971 W JP 0310971W WO 2004034908 A1 WO2004034908 A1 WO 2004034908A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- data
- projection data
- projection
- phase
- radiation
- Prior art date
Links
- 230000005855 radiation Effects 0.000 claims abstract description 193
- 238000012937 correction Methods 0.000 claims abstract description 56
- 238000004364 calculation method Methods 0.000 claims abstract description 52
- 238000000034 method Methods 0.000 claims description 124
- 238000003325 tomography Methods 0.000 claims description 36
- 230000008569 process Effects 0.000 claims description 35
- 230000008707 rearrangement Effects 0.000 claims description 20
- 238000001914 filtration Methods 0.000 claims description 7
- 238000001514 detection method Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 60
- 240000008042 Zea mays Species 0.000 abstract 1
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 abstract 1
- 235000002017 Zea mays subsp mays Nutrition 0.000 abstract 1
- 235000005822 corn Nutrition 0.000 abstract 1
- 230000006870 function Effects 0.000 description 29
- 230000000875 corresponding effect Effects 0.000 description 22
- 238000003384 imaging method Methods 0.000 description 17
- 238000005259 measurement Methods 0.000 description 13
- 238000010586 diagram Methods 0.000 description 12
- 230000007423 decrease Effects 0.000 description 7
- 230000001276 controlling effect Effects 0.000 description 5
- 238000011161 development Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- LKJPSUCKSLORMF-UHFFFAOYSA-N Monolinuron Chemical compound CON(C)C(=O)NC1=CC=C(Cl)C=C1 LKJPSUCKSLORMF-UHFFFAOYSA-N 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000006866 deterioration Effects 0.000 description 2
- 230000018199 S phase Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000013481 data capture Methods 0.000 description 1
- 230000002542 deteriorative effect Effects 0.000 description 1
- 238000010894 electron beam technology Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000001926 trapping method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/027—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4064—Arrangements for generating radiation specially adapted for radiation diagnosis specially adapted for producing a particular type of beam
- A61B6/4085—Cone-beams
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/421—Filtered back projection [FBP]
Definitions
- the present invention relates to a tomographic apparatus that generates a tomographic image of a subject using projection data obtained through a radiation detector from a radiation source that moves in a body axis direction relative to the subject.
- the conventional three-dimensional back projection method will be described.
- the F e 1 dkamp method, the Wang method, the IHCB method, and the PI-Method which have been proposed as three-dimensional back projection methods, use the slice (body axis) direction and channel ( This is a three-dimensional backprojection method in which a cone beam having a spread (tilt angle) in both directions is captured as a collection of multiple rows of fan beams having spread only in the channel direction.
- For the obtained fan beam projection data or the parallel beam projection data obtained by replacing the fans in parallel by the rearrangement process. Perform the same filter correction processing as in the two-dimensional back projection method, and follow the beam trajectory. To obtain a reconstructed image.
- Fig. 7 shows the reconfigurable conditions of the Wang method
- Fig. 8 shows the reconfigurable conditions of PI-Method.
- FOV indicates an effective visual field
- 3 ⁇ 0 indicates the distance between the tube and the orbital axis of the CT device
- S ID indicates the distance between the X-ray tube and the CT device.
- the Wang method is a method in which the Fe1dkamp method is adapted to imaging of a spiral trajectory, and the backprojection phase width is ⁇ to 2 ⁇ .
- PI-Method is disclosed in Japanese Patent Application Laid-Open No. H11-12553334. This is also a backprojection method corresponding to imaging of a spiral trajectory, but it is a reconstruction method that backprojects a ⁇ range with a different phase for each poxel in order to improve the bed moving speed in the Wang method.
- the backprojection phase range for each poxel can be set to ⁇ by restricting the X-ray beam to be backprojected in the vertical direction by a spiral opposite to the X-ray focal position.
- An example of the IHCB method is disclosed in JP-A-11-11823. This method consists of an algorithm that backprojects a different backprojection phase range for each poxel, whose backprojection phase width is either ⁇ or the total possible data range that varies for each poxel.
- the Feldkamp method is an image reconstruction method for capturing circular orbits, and does not support spiral orbital imaging.
- the Wang method is an image reconstruction method corresponding to imaging of a spiral orbit, and is performed by the conventional two-dimensional backprojection method by increasing the backprojection phase width beyond ⁇ (to increase data redundancy). Such effects of the movement of the subject can be corrected, but the efficiency of data utilization is low, and the helical pitch at the time of imaging (hereinafter referred to as measurement throughput) must be very small.
- the PI-Me thod and I HCB methods significantly improve measurement throughput compared to the Wang method by improving the backprojection phase range in the Wang method, but complete data redundancy. Because of the back projection method in the ⁇ range, which is excluded from the above, data discontinuity occurs at the start phase and end phase of the back projection phase range due to the movement of the subject, resulting in strong artifacts in the image. Cheap.
- Data redundancy refers to the extent of the phase range over which data of a certain phase is acquired in addition to the data of the opposite phase.
- data redundancy differs for each bath cell. For example, as shown in Fig. 22, when back projection is performed from data obtained by rotating the radiation source by 180 degrees, the data phase range contributing to each reconstructed pixel differs, and pixel a has a phase range of 180 degrees or more. Although the data is within the phase range, pixel b cannot have data below 180 degrees.
- the beam width in the body axis direction must be taken into account, and data redundancy varies in a complicated manner for each pixel. For this reason, complicated redundancy correction processing is required.
- the increased computation is as follows: (1) The increase in the number of executions of the one-dimensional rearrangement process; 2) increase the number of executions of reconstruction filter processing,
- SOD is the distance between the radiation source and the orbital axis
- ⁇ is the phase angle of the parallel beam
- J is the relative distance of the radiation source relative to the subject per rotation of the scanner on the radiation detector 13
- tj is the position in the channel direction the z-direction position of the radiation source 1 1 zs
- orbiting phase of the radiation source is a z s when 0 and z so. Therefore, if these operations can be simplified, a large increase in speed can be realized.
- an object of the present invention is to provide a tomographic apparatus capable of obtaining a high-quality tomographic image by suppressing distortion due to data discontinuity by using the data in the three-dimensional backprojection operation without losing data redundancy.
- an object of the present invention is to simplify the aresin operation relating to the fan-to-parallel beam conversion and the backprojection processing in the three-dimensional backprojection operation according to the set FOV range, without deteriorating the image quality. It is an object of the present invention to provide a tomography apparatus that achieves a significantly higher speed.
- the present invention provides a method in which a subject on a bed is arranged between a radiation source and a radiation detector which are arranged to face each other, and the radiation source and the radiation detector are rotated.
- the bed is configured to be movable with respect to the axis, the radiation emitted from the radiation source and transmitted through the subject is detected using the radiation source detector, and the target data is detected from the detected projection data.
- the reconstruction means can perform back projection in which the used projection data phase width is 180 degrees or more for each poxel.
- the projection data range is determined, a reconstruction filter is superimposed, the same or opposite phase data is weighted for each phase in this projection data range, and the filtered projection data is used to irradiate the radiation beam irradiation trajectory.
- the three-dimensional back-projection is performed along the trace over the determined back-projectable data range.
- the tomographic apparatus determines the projection data phase range to be used for each poxel, this phase range is set so that the absolute value of the radiation beam inclination angle is the same at both ends of the projection data.
- the projection data phase range can be determined for each pixel, thereby making it possible to use a small cone angle and projection data. By correcting the distortion, distortion due to discontinuity in the data phase direction can be suppressed and a high-quality image can be obtained. In addition, an image can be generated at high speed because redundancy processing that requires a complicated operation is not required.
- the data range is determined such that a difference between absolute values of cone angles at both ends of the projection data range to be used becomes small. It is characterized by doing.
- phase width of the projection data to be used is determined so as to be equal to the phase width for each poxel.
- the absolute value of the cone angle at both ends of the projection data range to be specifically used is determined by the determining means for determining the projection data phase range to be used for each of the vota cells.
- the characteristic is that the projection data range is determined so that the difference between the projection data is small, or the phase width of the projection data to be used is determined to be the same as that of each poxel. Data can be used.
- the absolute values of the radiation beam inclination angles at both ends of the projection data exactly equal, the position of the detector row direction can be calculated simultaneously from the start direction or end direction of the data, and the inverse of each reconstructed poxel can be calculated. Since the same phase range is used at the time of projection, the weighting function for redundancy correction can be determined by a single formula, and high-speed calculations can be performed.
- the present invention is further characterized in that, in the item described in the item 1, the projection data range in which the back projection can be performed is any one of 270 degrees and 360 degrees.
- the tomography apparatus uses a projection data range that can be back-projected using either 270 degrees or 360 degrees. By performing weighting using 270 °, discontinuity at the data edge can be reduced to the maximum. This 270 degree data corrects the discontinuity at the 180 degree data end with the least discontinuous data phase with a 90 degree phase difference to maximize the data discontinuity. Therefore, it is possible to realize a higher-quality image reconstruction.
- the present invention provides the method according to any one of Items 1 to 4, wherein the number of images per rotation is a multiple of the number C of sides of the polygon display pixel from the obtained projection data.
- the number of images per rotation is a multiple of the number C of sides of the polygon display pixel from the obtained projection data.
- Is converted to data containing fan beam data and parallel beam data, and a filter is superimposed on the projection data, and the projection phase in the same channel position and in the circling direction is ⁇ 2 (N l, 2, 3, ... ) [Rad]
- Different data are grouped in groups, and this group is back-projected to a rectangular image array.
- the tomography apparatus described in the above item 5 and item 6 is a method for speeding up the back projection which requires the most computation time in image generation.
- the fact that the reconstructed image array shape is a polygon and that imaging is performed while performing a circular motion around the reconstructed image is used, and is described in Item 5.
- fan beam reconstruction was performed by taking images with multiple views of the number of sides of the display pixels, and in the method described in Section 6, the number of views was changed to data that was a multiple of the number of sides of the display pixels by rearrangement processing. Convert and perform parallel beam reconstruction.
- the image can be generated by calculating the channel position within 1/4 of one round ( ⁇ 2 [rad]).
- the amount of computation is 1Z4 when considered in the case of full reconstruction, and a result close to the case of performing parallel computation using four computing units is obtained even though computation is performed with one computing unit, and low cost High performance can be achieved in a short time.
- any one of Items 1 to 6 a pixel interval in a body axis direction of an image by a polygon display pixel, and a relative movement of an object and the radiation source in a rotation axis direction.
- An association means for associating the distances with each other is provided.
- association means is configured such that a relationship between a pixel interval rpitch in the body axis direction of the rectangular image and a relative movement speed J of a target object and the radiation source in the rotation axis direction is satisfied.
- the present invention provides a method in which a subject on a bed is arranged between a radiation source arranged opposite to the radiation source and a radiation detector constituted by two-dimensionally arranging detection elements, and the radiation source and the radiation detector are circulated.
- the bed is movable with respect to this orbit.
- the radiation emitted from the radiation source and transmitted through the subject is detected by the radiation detector.
- a phase range is determined, and a curve representing the radiation source position with respect to the channel direction position of the parallel beam projection data corresponding to the region of interest obtained by the parallel beam having a parallel shape viewed from the direction of the circling axis emitted from the radiation source is determined.
- Approximate straight lines are calculated and corrected by multiplying each row of projection data by a coefficient that depends on the tilt angle of the radiation from the radiation source.
- the fan shape is viewed from the direction of the orbit axis emitted from the radiation source.
- the parallel projection data is superimposed with a reconstruction filter to generate filter-processed parallel projection data.
- the parallel beam projection data that has undergone the filtering process is generated.
- the three-dimensional back-projection is performed on the back-projection region corresponding to the region of interest along the near-field irradiation locus using the approximate straight line.
- the tomographic apparatus converts filtered parallel beam projection data based on the projection data range that can be back-projected determined by the use data phase range calculation means into a parallel shape as viewed from the direction of the orbit axis emitted from the radiation source. Approximate irradiation of the radiation beam calculated using the approximate straight line by the approximate straight line calculation means for calculating the approximate straight line for the curve indicating the radiation source position with respect to the channel direction position corresponding to the region of interest of the parallel beam projection data obtained by the parallel beam Since the three-dimensional backprojection was performed along the trajectory onto the backprojection area corresponding to the region of interest, the calculation of the focal position of the conventional parallel beam included the arcsin calculation, and the calculation addition was large. Is replaced with an approximation line, simplifying the computational complexity of the parallel beam 3D back projection method. It is possible to significantly speed up turned into.
- the reconstructing means is arranged to correct the data redundancy in each phase in accordance with the determined phase width of the projection data in the phase direction.
- Weighting for redundancy correction for generating a weighting coefficient from a weighting function is performed, and the parallel beam three-dimensional backprojection unit obtains the projection data in the determined projection data phase range by the redundancy correction weighting unit.
- Weighted weighting factor and back projection area It is characterized in that three-dimensional back projection is performed on the area along the approximate trajectory.
- the tomographic apparatus described in Item 11, further comprising: determining the phase range of ⁇ [rad] in the view direction in determining the phase range of the projection data; and determining the weight by the weighting means for redundancy correction.
- Redundancy correction can be performed by a function. For this reason, the data is made redundant (backprojection phase width is made larger than 180 degrees), and the discontinuity at the end of the data (at the end of shooting) is reduced by weighting with a weight function. It is possible to obtain an image in which the influence of the movement of the subject is reduced as much as possible.
- the use data phase range calculating means includes: setting a projection data range that can be back-projected for each reconstructed poxel to a maximum of a beam back-projected for each poxel; The cone angle is determined to be the narrowest.
- the backprojection phase range for each Vota cell is determined by the use data phase range calculation means so that the maximum cone angle is minimized, and therefore, the image quality based on the cone angle is determined. It is possible to reduce the influence of deterioration and obtain good image quality, and improve the relative moving speed of the subject and the focus in the Z direction (so-called measurement throughput).
- the projection data range that can be back-projected for each reconstructed Vota cell is changed in the phase direction range of the beam back-projected for each poxel. Is determined to be as narrow as possible.
- the back-projection phase range for each poxel is determined by the use data phase range calculation means so as to reduce the number of views, so that the time resolution for each poxel is improved. it can.
- good image quality can be obtained in a portion where the subject moves greatly.
- the back projection phase range for each poxel is set to the time range captured at the same time as possible so that the time position of each poxel in the display image is close, the time width contributing to the reconstructed image is reduced. It can be shortened to improve the time resolution.
- FIG. 1 is a schematic view showing a tomography apparatus according to the present invention.
- FIG. 2 is a block diagram of the tomography apparatus shown in FIG.
- FIGS. 3A and 3B are conceptual diagrams showing focal trajectories of a circular orbit scan and a spiral orbit scan.
- Figures 4A and 4B are the waist side views of a single radiation detector and a multi-row radiation detector.
- Figures 5A and 5B show the collimation thickness of the X-ray beam per row for a single radiation detector and a multi-row radiation detector.
- FIG. 6 is a flowchart showing a processing operation by a general reconstruction means.
- FIG. 7 is a plan view showing the positional relationship between the radiation detector 13 and the beam as reconfigurable conditions by the Wang method.
- FIG. 8 is a plan view showing the positional relationship between the radiation detector 13 and the beam as a reconfigurable condition by the PI-Method method.
- FIG. 9 is a flowchart showing the processing operation of one embodiment of the reconfiguration means of the present invention.
- FIG. 10 is a flowchart showing the operation of the use data position range determination processing shown in FIG.
- Figures 11A and 11B are perspective views showing the spiral trajectory of the radiation source and radiation detector.
- 12A and 12B are a perspective view and a development view illustrating the operation of the rearrangement processing shown in FIG.
- FIGS. 13A and 13B are a perspective view and a development view illustrating the operation of another sorting process.
- Figure 14A shows the spiral measurement when using 180 degree data.
- FIG. 14B is a characteristic diagram showing a weighting function corresponding to spiral measurement when using 180-degree data.
- Fig. 15A shows the spiral measurement when using the 270 degree data.
- FIG. 15B is a characteristic diagram showing a weighting function corresponding to another spiral measurement when using 270 ° data.
- FIG. 16 is a flowchart showing the processing operation of another embodiment of the reconstructing means of the present invention.
- FIG. 18 is a perspective view showing an example of the back projection data phase range when the phase range is between 180 degrees and 360 degrees (1 ⁇ f ⁇ 2).
- FIGS. 20A and 20B are characteristic diagrams of a weighting function for explaining the redundancy weighting process shown in FIG.
- 21A to 21C are characteristic diagrams of a weighting function for each phase for explaining the redundancy weighting process shown in FIG.
- FIG. 22 is a plan view showing a phase range in which back projection is possible.
- FIG. 23 is a flowchart showing the processing operation of another embodiment of the reconstructing means of the present invention.
- FIG. 24 is a flowchart showing the processing operation of another embodiment of the reconstructing means of the present invention.
- FIG. 25 is a view for explaining backprojection for each group in the grouping process shown in FIG.
- FIG. 26 is a view for explaining back projection in another group unit in the grouping process shown in FIG.
- FIG. 27 is a flowchart showing the processing operation of another embodiment of the reconstructing means of the present invention.
- FIGS. 28A, 28B, and 28C are diagrams illustrating the calculation process of the approximate straight line illustrated in FIG. 27, and a diagram illustrating the relationship between the cone angle, the X-ray source, and the reconstruction.
- FIG. 29 is a diagram for explaining a process of determining a projection data phase range.
- FIG. 30 is a flowchart showing the processing operation of another embodiment of the reconstructing means of the present invention.
- FIG. 31 is a flowchart showing the operation of the phase range calculation processing of the use data of FIG. 30;
- FIGS. 32A to 32D show the relationship between the fan beam and the parallel beam.
- FIG. 1 is an external configuration diagram showing a tomography apparatus according to an embodiment of the present invention.
- Scanner 1 used for imaging, couch 2 for moving with object to be tested placed on it, measurement consisting of mouse and keyboard, etc.Input for inputting reconstruction parameters such as couch movement distance information and reconstruction position It has a device 3, an arithmetic device 4 for processing data obtained from the multi-row detector, and a display device 5 for displaying a reconstructed image.
- FIG. 2 is a block diagram showing a main part of the tomography apparatus shown in FIG.
- the scanner 1 includes a bed 2, a radiation source 11 such as a radiation generator having a high-voltage switching unit 8, a high-voltage generator 9, and a radiation controller 10, and a radiation source 11 via a subject 12.
- Radiation detectors 13 arranged opposite to each other, orbiting drive device 14 for driving this radiation detector 13 and radiation source 11 in the orbital direction, collimator 15 for controlling the radiation area irradiated from radiation source 11 etc. have.
- a control signal required for imaging is sent from the central controller 20 to the radiation controller 10, the bed controller 18 and the scanner controller 17, and the imaging is started in response to the imaging start signal.
- a control signal is sent from the radiation control device 10 to the high voltage generator 9, a high voltage is applied to the radiation source 11, and the radiation is emitted from the radiation source 11 to the object 12. Is done.
- a control signal is sent from the scanner control device 17 to the orbiting drive device 14, and the radiation source 11, the radiation detector 13 and the preamplifier 21 orbit relative to the object 12.
- the couch 2 on which the subject 12 is placed is stopped by the couch controller 18 during circular orbit scanning, and is translated in the orbital direction of the radiation source 11 and the like during spiral orbit scanning.
- the orbiting drive 14, the scanner controller 17 and the bed controller 18 make the radiation source 11 and the radiation detector 13 orbit relative to the subject 12 and the subject 1.
- a drive device that is relatively movable in the axial direction is configured.
- the irradiation area of the radiation emitted from the radiation source 11 is limited by the collimator 15, absorbed and attenuated by each tissue in the subject 12, passed through the subject 12, and detected by the radiation detector 13. You.
- the radiation detected by the radiation detector 13 is converted into a current, amplified by the preamplifier 21, and input to the arithmetic unit 4 as a projection data signal.
- the projection data signal input to the arithmetic unit 4 reconstructs the image in the arithmetic unit 4. It is processed by the reconstruction means 22.
- the reconstructed image is stored in the storage device 23 in the input / output device 50, and displayed as a tomographic image on the display device 5 by the image processing device 26.
- FIGS. 3A and 3B are conceptual diagrams showing focal orbits of a circular orbit scan and a spiral orbit scan.
- Figure 3A shows the trajectory 24a of the radiation source (focal point) during a circular orbit scan
- Figure 3B shows the trajectory 24b of the radiation source (focal point) during a helical orbit scan.
- the detector row was one row
- the image was captured in a circular orbit as shown in the movement trajectory 24 b.
- the image of the radiation source position was accurately reproduced by performing filter-corrected two-dimensional back projection.
- Can be When the image is taken in a spiral trajectory like the movement trajectory 24b, streak-like artifacts are found only at the position of the imaging end due to the discontinuity of the data due to the discontinuity of the data only in the filtered 2D back projection. Will happen.
- the degree of artifact in this case is determined by the degree of discontinuity in the X-ray source trajectory, that is, the degree of artifact varies depending on the moving speed of the subject.
- the helical pitch the ratio of the moving speed of the object to the thickness of the X-ray beam in the orbital direction
- 4A and 4B are schematic side views showing a single-row radiation detector 13a and a multi-row radiation detector 13b, respectively.
- the multi-row radiation detectors 13b which are narrower per row than the single-row radiation detectors 13a in Fig. 4A, are arranged in multiple rows in the orbital axis direction. Detector wider than detector 13a is realized.
- FIG. 3 is a schematic side view showing a thickness (hereinafter, referred to as a detector collimation thickness DCT).
- a detector collimation thickness DCT a thickness
- the detector collimation thickness DCT is thinner than that of the single-row radiation detector 13a shown in Fig. 5A.
- the spatial resolution (body axis resolution) of the obtained tomographic image in the rotation axis direction improves as the detector collimation thickness decreases.
- FIG. 6 is an example of a flowchart showing a process in the Fe1dkamp reconstruction method. '
- FIG. 9 is a flowchart showing the processing operation of the reconstructing means 22 in the tomography apparatus according to one embodiment of the present invention. Note that this flowchart assumes processing performed for each slice.
- reconstruction means 22 is used data phase range calculation means for determining a projection data phase range that can be back-projected for each reconstruction voxel, and a tilt angle of radiation from a radiation source for each column of projection data.
- Cone angle correction means for multiplying by a coefficient depending on the temperature, and one-dimensional rearrangement processing means for obtaining parallel beam projection data from fan beam projection data obtained from a fan beam having a fan shape viewed from the orbital axis emitted from the radiation source
- Filter correction means for generating a filtered parallel beam projection data by superimposing a reconstruction filter on the parallel beam projection data; and a filtered parallel beam projection data based on the determined backprojectable projection data range.
- a parallel beam three-dimensional backprojection means for performing three-dimensional backprojection on the backprojection region corresponding to the region of interest Have.
- step S4 a data range to be used for each poxel is determined by the used data phase range calculating means.
- step S5 each row of the projection data is multiplied by a coefficient depending on the tilt angle of the radiation by the cone angle correction means, and in step S6, the orbit emitted from the radiation source by the one-dimensional rearrangement processing means is multiplied.
- the fan beam projection data obtained from the fan beam having a fan shape viewed from the axial direction is associated with the parallel beam projection data.
- step S7 the reconstruction data for the parallel beam projection data is A filter is superimposed to generate filtered parallel beam projection data.
- step 8 the filtered parallel beam projection data is three-dimensionally back-projected to the back projection area corresponding to the region of interest based on the projection data range that can be back projected determined by the parallel beam three dimensional back projection means.
- step S4 the data range to be used is determined by the used data phase range calculating means for each poxel for all the poxels in the slice.
- the distance between the radiation and radiation source 11 and the center of rotation is SOD, and the scanner on the radiation detector 13 3
- the radiation source for the subject per rotation 1 The relative movement distance in the body axis direction of 1 (for example, the table feed amount) is J, the orbital phase of the fan beam source is] 3, the beam toward the reconstructed poxel I (x, y, z) and the beam relative to the center beam
- the traveling direction of the parallel beam is W and the direction perpendicular to the traveling direction W (channel direction of the parallel beam) is T
- the W coordinate is represented by Equation 4 and Equation 5, respectively.
- T, y, ⁇ ) x * c os 0 + y * s i n ⁇ i)
- a parallel beam of ⁇ phase forms a reconstructed water cell I (x, y, z).
- V-axis of the radiation detector 13 when it crosses the radiation detector 13 whose distance from the radiation source 11 is SID (in the same direction as the z-axis, its origin is at the center of the detector) Let H (x, y, ⁇ ) be the coordinate of the system formed by X and Y—Y axes. Note that the ⁇ ⁇ axis and the V axis match, but the ⁇ axis differs from the scan start position as the origin position, and the V axis has the detector center as the origin position.
- H (x, y, ⁇ ) (z-J ( ⁇ + ⁇ ) / 2 ⁇ ) ⁇ SID / (s_tz_dist (, y, 0) + W (x, y, ⁇ ))
- phase range determined in step S4 based on FIG. 28C will be described. Since this phase range differs for each vota cell as described above, it is determined for each voxel, and is represented by f here.
- the phase range is determined by selecting a plane beam opening angle of the cone beam with respect to the vertical line in the detector direction, that is, the one that minimizes the cone angle. F is usually between 1 and 2.
- FIG. 10 shows an example of an algorithm for calculating the data range with the smallest cone angle.
- the initial value of ⁇ is 1 f ⁇ 2, and the calculated phase accuracy is Q (Q is usually the phase angle per view, but if processing time is prioritized, one view is used.
- Q is usually the phase angle per view, but if processing time is prioritized, one view is used.
- (X, y, ⁇ ) and the sum of ⁇ (x, y,) and err (, y, ⁇ ) (hereafter referred to as err), and this sum err
- the phase range (Bs ⁇ 0 and Be) is expressed by the following equations (1) and (2).
- the force for which the phase range was determined using the simplest method as described above is the minimum of the function err (0) in the phase range (one f ⁇ / 2 ⁇ ⁇ and f ⁇ / 2).
- This is a value calculation problem, using existing methods, such as the Brents method (B rents Method) or the golden section method (goldensection search), or a combination of these methods to minimize err ( ⁇ ) and ⁇ + ⁇ + f ⁇ May be calculated.
- Speed can be increased.
- phase range f ⁇ used to backproject the reconstructed poxel ( ⁇ ⁇ , y, z) is such that the absolute value of the beam tilt angle (cone angle) of the X-ray beam at the end of the phase range ⁇ is small.
- B s and ⁇ B e may be determined by determining the orbital phase ⁇ of the parallel beam and extending the data range to both ends of the data range as shown in Expressions 14 and 15 below.
- the filter correction in the reconstruction is a filter according to the distance from the revolving axis in the reconstructed image, and it is necessary to apply a filter according to the cone angle to correct the effect of the beam tilt.
- the data before filter correction is P para ( ⁇ , t, v)
- the data after filter correction is f ⁇ ⁇ ara ( ⁇ , t, v)
- the reconstruction filter function is g (t)
- the filter processing can be expressed as shown in Equation 16 above, and the cone angle correction is the part shown in Equation 17 among them.
- the cone angle correction term is a coefficient corresponding to the detector row position V (cone angle), and therefore can be performed before and after filter correction.
- a known technique in a three-dimensional back projection method such as the Fe 1 dk amp method can be applied.
- Equation 16 t 'is an integral variable in Equation 16.
- FIG. 9 the reordering process (rebi-nning) by the one-dimensional reordering unit in step S6 shown in FIG. 9 will be described.
- Figures 32A and 32B show the relationship between the fan beam and the parallel beam.
- 32A to 32C show the 180 ° reconstruction in the fan beam
- FIG. 32D shows the 180 ° reconstruction in the parallel beam.
- the fan beam illuminated like a fan as viewed from the orbital axis as shown in Fig. 11A and Fig. 11B is applied to Fig. 12A (Fig.
- the one-dimensional rearrangement process of rearranging into parallel beams parallel to the orbital axis as shown in Fig. 12B is performed.
- Figures 13A and 13B show parallel beams that can be arranged in the orbital axis direction, and will be described later.
- Equation 18 The reordering process can be expressed as shown in Equation 18.
- filter correction There are two types of filter correction, a convolution method that performs convolution in real space and a Fourier method that performs multiplication in Fourier space.
- the former method of convolution is the convolution of the inverse Fourier-transformed filter function in real space.
- the latter Fourier method is a process that uses the Fourier transform to convert to Fourier space, multiplies it by a filter function (spatial frequency filter), and then applies the inverse Fourier transform. '
- Equation 29 real space filtering by the convolution method is as shown in Equation 29 when the inverse Fourier transform g (t) of G ( ⁇ ) is expressed as shown in Equation 28.
- the application direction of the filter is set to the T direction for simplicity, but the filter can be applied to a high-dimensional direction that combines the V direction, the T direction, and the ⁇ direction.
- the projection data is treated as continuous data.
- the projection data is actually discrete data, it has to be calculated discretely using a known interpolation method. This discrete calculation method has been performed so far, and is the same as the filter correction used in the weighted spiral correction reconstruction.
- the reconstructed poxel is denoted by I (X, y, z), and the V-axis position corresponding to the orbital axis on the cylindrical detector centered on radiation source 11 is denoted by v, This V axis Assuming that the position on the T axis that is almost orthogonal to is t, the reconstructed poxel I (x, y, z) is shown in Equation 22.
- V (z—J0 / 2 ⁇ ) -SID / (SODcosa—xsin0 + ycos0) in order to calculate at a higher speed without sacrificing accuracy.
- each pixel is back-projected from the data with the least error (data with a small cone angle) for each pixel, and used for each pixel to achieve this.
- Determining means for determining a data phase range specifically, the phase range is determined by determining a data phase range for each bath cell so that the absolute value of the radiation beam inclination angle is the same at both ends of the data; It is possible to use projection data with a small cone angle, and to obtain an image with reduced discontinuity in the data phase direction by performing correction using a weighting function for each poxel while maintaining redundancy. it can.
- a weighting function for each poxel while maintaining redundancy. it can.
- 270 degrees of data in the phase direction and performing weighting as shown in Fig. 2OA it is possible to reduce discontinuity at the end of data to the maximum.
- This 270-degree data consists of the discontinuity points at the 180-degree data end as shown in Figs.
- the weight function for redundancy correction is determined by only one equation, so that the calculation can be performed at high speed.
- the errors are mutually corrected at the positions of ⁇ ⁇ 2 and ⁇ in Fig. 15B.
- FIG. 16 is a flowchart showing the processing operation of the reconfiguration means 22 according to another embodiment of the present invention.
- step S9 weighting for redundancy correction is performed in step S9, and then three-dimensional back projection is performed in step S8. ing.
- the reconstruction means further changes the shape according to the phase width with respect to the filtered projection data over the projection data range f ⁇ obtained by the use data phase range means.
- Fig. 20 ⁇ and Fig. 21A- 21C weighting that changes linearly in the view direction is performed, but weighting that changes nonlinearly as shown in Fig. 20B may be performed.
- the nonlinear weight function W "(") shown in FIG. 20B can be calculated from the weight function W ( ⁇ ) as shown in Equations 26 to 28 as an example.
- B ⁇ 2 ⁇ is described, but the case where ⁇ > 2 ⁇ can be easily calculated using the same concept.
- Such a tomographic apparatus is a three-dimensional reconstruction method that determines the phase range of f ⁇ [rad] in the view direction in the above-described phase range calculation processing for each Vota cell, and performs redundancy correction using a weight function.
- a weight function By making the data redundant (backprojection phase width greater than 180 degrees) and by weighting it with a weight function, it is possible to reduce discontinuities at the data edges (at the start / end of imaging). , Minimizing the effects of subject movement Image can be obtained.
- the redundancy (projection phase range) of the obtained projection data differs for each poxel.
- the data phase range that can be backprojected differs for each reconstructed pixel. Can obtain data in the phase range of 180 degrees or more, but pixel b can only obtain data of 180 degrees or less.
- the data redundancy differs for each pixel, when backprojecting from projection data of 360 degrees or less, complicated redundancy correction processing is required at the time of backprojection.
- redundancy correction process also relates to the measurement throughput (the relative movement speed between the focus and the subject).
- redundancy is used, and data with a phase range of 180 ° or more is used for each poxel to prevent discontinuity due to movement and the like, thereby improving data efficiency. Also leads to.
- FIG. 23 is a flowchart showing a processing operation of the reconfiguration means 22 according to still another embodiment of the present invention.
- step S11 rearrangement processing is performed on projection data shot in a multiple view of four. I have to.
- step S7 the filter is corrected.
- step S9 weighting processing for redundancy correction is performed in step S9, and the projection data grouped in step S8 is back-projected into a square image in group units.
- a means for acquiring projection data in which the number of shots per rotation is a multiple of 4 is provided, and the reconstructing means 22 is a means for superimposing a filter on the projection data.
- grouping means for grouping data having the same channel position and different projection phases in the circumferential direction by N ⁇ ⁇ 2 (N 1, 2, 3, ......) [rad].
- the reconstructed image array shape is rectangular and that shooting is performed while moving circularly around the reconstructed image.
- the image is photographed in a multiple view of 4 and the fan beam reconstruction is performed.
- the view number is converted to data of a multiple of 4 by the rearrangement process and the parallel beam reconstruction is performed. I have to.
- N 1, 2, 3, ......) [r a d].
- the amount of computation is 1 Z 4 in the case of full reconstruction, and a result close to that obtained by performing parallel computation using four computing units despite obtaining one computing unit is obtained. . In other words, high performance can be realized at low cost.
- the number of views at the time of shooting can be made a multiple of 4, and it is possible to reconstruct directly from the fan beam without performing reordering (rebinning).
- the display pixel is polygonal and has a C-side, the above-mentioned circumferential phase is 2 ⁇ / C [r ad]
- the beam is radiated from the focal position S ( ⁇ ), which is the phase j8, and passes through the poxel (X, y).
- the backprojection processing for each group is expressed by Equations 29 to 32. Is done.
- Reference numeral 13 1 indicates a reconstruction area.
- the beam radiated from phase; 3 + 2 ⁇ and passing through the poxel (x, y, z + J) is applied to the radiation detector in the same manner as when the poxel I (x, y, z) is radiated from phase j3. Irradiation is performed at the position V in the rotation axis direction. Similarly, the beam emitted from the phase] 3 + ⁇ passes through the Vota cell I (- ⁇ , 1 y, z + J / 2) and is emitted to the orbital position V on the radiation detector.
- the MDCT correlates the pixel spacing of the poxel in the body axis direction and the relative movement speed of the object and the radiation source in the orbital axis direction, thereby obtaining the body axis direction.
- the position can be calculated at a high speed, and in the case of generating an image from data of a plurality of rounds obtained by photographing with a spiral run, back projection that requires the longest time for image generation can be speeded up.
- the helical period in the body axis direction and the period in the body axis direction of the reconstructed voxel are synchronized.
- the beam passing through the poxel I (_x, 1 y, ⁇ ⁇ J / 2 + Z) whose axial position is (N ⁇ J / 2) + Z [mm] are the positions on the radiation detector where they intersect.
- the use data phase range calculation for each pixel cell can be similarly performed by, for example, changing H (X, y, ⁇ ) in the above equation to equation 33.
- the tomographic apparatus can be applied to apparatuses using X-rays, gamma rays, neutron rays, positrons, electromagnetic energy, and light.
- the scanning method is not limited to the first to fourth generations, and is not limited to the deviation method, and is also used for multi-tube CT and donut-shaped tube CT equipped with a plurality of radiation sources.
- the shape of the radiation detector is also a detector arranged on the surface of the cylinder centering on the radiation source, a flat detector, a detector arranged on a spherical surface centering on the radiation source, and a cylinder centering on the orbit axis It can be applied to any radiation detector such as a detector arranged on the surface.
- the position of the radiation detector corresponding to the reconstructed poxel was calculated each time.
- N / 4 circumference is an example in the case where the shape of the display pixels is a quadrangle, and if the display pixels are hexagonal, they may be grouped by N / 6 circumferences.
- FIG. 27 is a flowchart showing the processing operation of the reconstruction means 22 in the tomography apparatus according to one embodiment of the present invention.
- the reconstruction means 22 is a use data phase range calculation means for determining a projection data phase range that can be back-projected for each reconstruction poxel, and a parallel shape parallel when viewed from the direction of the orbit axis emitted from the radiation source.
- Approximation straight line calculation means for calculating an approximation straight line for a curve indicating a radiation source position with respect to a channel direction position corresponding to a region of interest of parallel beam projection data obtained by a beam; and a radiation source for each column of projection data.
- Cone angle correction means for multiplying by a coefficient depending on the inclination angle of the radiation, and parallel beam projection data obtained from fan beam projection data obtained from a fan beam having a fan shape viewed from the direction of the orbit emitted from the radiation source 1
- Filter correction means for generating parallel beam projection data, and interest in filtered parallel beam projection data based on the determined back-projectable projection data range along the approximate irradiation trajectory of the radiation beam calculated using an approximate straight line
- a parallel beam three-dimensional backprojection unit for performing three-dimensional backprojection on a backprojection area corresponding to the area is provided.
- step S4 the data range to be used for each poxel is determined by the use data phase range calculation unit, and in step S14, the approximate straight line calculation unit determines the data range to be used from the orbit axis emitted from the radiation source An approximate straight line is calculated for the curve indicating the radiation source position with respect to the channel direction position of the parallel beam projection data obtained from the parallel shaped parallel beam.
- step S5 each row of the projection data is multiplied by a coefficient depending on the inclination angle of the radiation by the cone angle correction means, and in step S6, the rays are emitted from the radiation 3 ⁇ 4g by the one-dimensional rearrangement processing means.
- the beam projection data is associated with the parallel beam projection data.
- a reconstruction filter is superimposed on the parallel beam projection data by the filter correction means to generate filtered parallel beam projection data.
- filtered parallel beam projection data is calculated using an approximate straight line. Three-dimensional back projection is performed on the back projection area corresponding to the region of interest along the approximate irradiation trajectory.
- Steps S4 to S7 are the same as those already described in the other embodiments.
- the z-coordinate Z i of the focal point at the channel i position of the parallel beam is given by the following equation 34, and the approximated straight line z A excluding the arcsin operation is given by equation 35 Is performed.
- the position of the normal beam in the t-axis direction in the channel i is denoted by ti.
- Equation 36 is minimized to determine A and B. At the minimum, the derivatives of A and B in Eq. 44 are zero, as shown in Eqs. 37 and 38.
- Equation 39 Equation 36 and Equation 37
- Equation 4 5 Can be obtained.
- step S shown in FIG. 2 7 15 Based on the determined projection data range that can be back-projected, the parallel beam projection data filtered is converted to the inverse projection corresponding to the region of interest along the approximate irradiation trajectory of the radiation beam calculated using the approximate straight line.
- the parallel beam three-dimensional back projection means for performing three-dimensional back projection on the projection area will be described.
- the reconstructed poxel is represented by I (x, y, z), the relative movement distance of the radiation source 11 to the subject per one rotation of the scanner on the radiation detector], and the radiation source
- V be the position along the orbital axis on the cylindrical radiation detector 13 centered at 1 1
- t be the position on the T axis that is almost perpendicular to this
- T (x, y, ⁇ ) be the coordinates. 46 to number 50 are obtained.
- ⁇ indicates an image array in which I (x, y, z) are gathered
- 111 indicates an X-ray beam.
- Equation 46 shows the weighted three-dimensional back projection along the beam trajectory over the back projection data range determined by the data phase range calculation means.
- Equation 50 shows the radiation beam trajectory calculated using the approximate straight line.
- V (ZI—Zs) ⁇ SID / L ( ⁇ , x, y) [Equation 5 0
- Equation 50 shows that the calculation of arcsin is included in the conventional calculation of the focal point z position of the parallel beam, and the calculation load is large, as can be seen from the comparison with Equation 1.
- this arcsin operation is replaced with an approximate straight line, the amount of operation of the parallel beam three-dimensional backprojection method can be simplified and a large increase in speed can be achieved.
- the phase range of f ⁇ [rad] is determined in the view direction, and a three-dimensional reconstruction method is performed in which redundancy correction is performed using a weight function.
- the same backprojection phase range is used for all the bottom cells.
- the position of the focal point in the z direction in the channel direction in the parallel beam is not the same, and the maximum cone angle that is back-projected at each poxel becomes large.
- Increasing the maximum cone angle used means that a wider detector is required in the z-direction of the revolving axis, that is, the relative movement speed of the subject and the focal point in the z-direction.
- the measurement throughput decreases.
- the maximum cone angle of the beam used for back projection is reduced as described above, it is possible to improve the measurement throughput by reconstructing with a narrow detector in the Z direction.
- FIG. 30 is a flowchart showing the processing operation of the reconstruction means 22 in the tomography apparatus according to another embodiment of the present invention.
- the reconstruction means 22 are used data phase range calculation means for determining the projection data phase range that can be back-projected for each reconstruction Vota cell, and a parallel shape parallel when viewed from the direction of the orbit axis emitted from the radiation source.
- Approximation line calculation means for calculating an approximation line for a curve indicating the radiation source position with respect to the channel direction position of the parallel beam projection data obtained by the beam; and a coefficient depending on the inclination angle of the radiation for each column of the projection data.
- Cone angle correction means for multiplying, one-dimensional rearrangement processing means for associating fan beam projection data obtained from a fan beam having a fan shape viewed from the direction of the circling axis emitted from the radiation source with parallel beam projection data, A filter complement that generates a filtered projection data by superimposing a reconstruction filter on the corrected projection data Means and a redundancy function for performing redundancy correction using a weighting function whose shape changes according to the phase width for the filtered projection data in one of the projection data ranges f ⁇ determined by the used data phase range calculation means.
- the reconstructing means 22 generates filter-processed parallel beam projection data by superimposing a reconstruction filter on the parallel beam projection data by the filter correcting means in step S7, as in the case of FIG. 27.
- the redundancy correction weighting means performs redundancy correction using a weighting function on the filtered projection data generated by the filter correction means over the data range f ⁇ determined by the use data phase range calculation means. .
- the filtered parallel beam projection data is filtered based on the projection data range that can be back-projected determined by the parallel beam three-dimensional back projection means.
- three-dimensional back projection is performed on the back projection area corresponding to the region of interest along the approximate irradiation trajectory of the radiation beam calculated using the approximate straight line.
- step S4 the data range is determined for each poxel so that the difference in cone angle (corresponding to the absolute value of V) at the end of the back projection data range (data start ⁇ end position) is minimized. Similar reconstruction processing can be performed based on the determined data range.
- phase range can be determined by the following method.
- the distance between the radiation source and the rotation center is S OD
- the scanner on the radiation detector is J
- the relative movement distance of the radiation source with respect to the subject per rotation is 3
- the rotation phase of the fan beam source is 3
- the beam and center toward the reconstructed poxel Assuming that the beam opening angle with respect to the beam is ⁇ and the orbital phase of the parallel beam is ⁇ , the radiation source position S ( ⁇ ) ⁇ S (x s , y s , z s ) is given by the above equation (2).
- Equation (12) Considering this as a parallel beam by the rearrangement process, the above equation (12) is obtained.
- the traveling direction of the parallel beam is w
- the direction perpendicular to w is t
- the t coordinate and the w coordinate when the ⁇ phase parallel beam passes through the coordinates (x, y) are as described above.
- Equations (13) and (14) are obtained, and the distance between the radiation source and the tz plane (the plane perpendicular to the parallel beam passing through the orbital axis) is represented by Equation (6) described above.
- the reconstructed voxel I (x, y, z) is back-projected in the phase range where the difference in the backprojection phase range for voxels whose z-position is in the same plane is small.
- the end position B s phase range f [pi to be used, z-direction position of the radiation detector beam from B e irradiation Isa crosses the transmitted radiation detector reconfiguration Pokuseru within the detector range since it is necessary to position, when theta 0 orbiting phase when the z-direction position of the focal point is the reconstruction voxel position, the number 3 Oyopi number 5 4 And ⁇ should be chosen to be as close as possible to ⁇ 0 — ⁇ ⁇ / 2. [Equation 53]
- the time resolution for each poxel can be improved by determining the backprojection phase range for each poxel so that the number of views is small, and by combining with the weighted backprojection described earlier, Good image quality can be obtained in a part where the movement of the subject is large. Furthermore, by setting the backprojection phase range for each poxel to the time range taken at the same time as possible so that the time position of each poxel in the display image becomes closer, the time width contributing to the reconstructed image is reduced. Shorten it Thus, the time resolution can be improved.
- the backprojection phase range is ideally the same for all bottle cells, but the relative movement speed of the subject and the focus is high, and data acquisition in the same phase range for all poxels is not possible.
- the backprojection phase range can be determined for each poxel as close as possible.
- a body axis (circumferential axis) direction filter that can change the spatial frequency characteristic in the column direction with respect to the projection data is preferably used.
- the superposition of the body axis direction filter may be performed before the one-dimensional rearrangement processing or may be performed after the one-dimensional rearrangement processing. Also, it may be included in the filter correction processing.
- the tomographic apparatus using X-rays is used.
- the present invention is not limited to this, and the present invention is also applicable to a tomographic apparatus using neutron rays, positron-gamma rays, or light.
- the scanning method is not limited to any of the first, second, third, and fourth generation methods, and includes a multi-tube CT, a force-sword scan CT, and an electron beam equipped with a plurality of radiation sources. It can also be used for CT.
- the radiation detectors are shaped like a radiation detector, a flat detector, a radiation detector placed on a spherical surface around the radiation source, and a circular axis around the radiation source. It can be applied to any radiation detector such as a radiation detector placed on the surface of a cylinder.
- the tomography apparatus is not limited to a spiral orbital scan, but can be applied to a circular orbital scan.
- interpolation methods such as Lagrange interpolation are actually used, and the phase direction, the column direction of the radiation detector, It is desirable to calculate discretely by interpolation in three directions in the channel direction.
- arcsin is approximated by one approximation line, but is approximated by using a plurality of approximation lines (using different approximation lines according to the distance from the orbital axis). Is also good.
- non-linear function value in the present invention may be calculated using a pre-calculation (tabulation) of calculated values and interpolation for speeding up.
- the process (S 4) of FIG. 9 for determining the phase range of the projection data used for each poxel can be applied to other embodiments of the reconstruction means 22.
- the redundancy correction weighting process (S 9) of FIG. 16 can be applied to other embodiments of the reconfiguration unit 22.
- the process (S11) of rearranging the data into a multiple view data of 4 (S11) and the process of grouping (S12) in FIG. 23 can be applied to other embodiments of the reconstruction means 22.
- the reordering process (S13) in FIG. 24 can be applied to other embodiments of the reconstructing means 22.
- the approximate straight line calculation processing (S14) and the parallel beam three-dimensional back projection processing (S15) in FIG. 27 are applicable to other embodiments of the reconstruction means 22.
- the tomography apparatus of the present invention when reconstruction is performed from data obtained by scanning, distortion due to data discontinuity can be reduced as much as possible, and streak is included in the reconstructed image. High-quality images can be obtained without causing artefacts.
- a curve indicating the position of the radiation source with respect to the position in the channel direction of parallel beam projection data obtained by a parallel beam having a parallel shape viewed from the direction of the orbit axis emitted from the radiation source By calculating an approximation straight line to, the arcsin operation used so far can be simplified, the speed can be greatly increased, and a high-quality image can be obtained in a short time.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Theoretical Computer Science (AREA)
- Pulmonology (AREA)
- General Physics & Mathematics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Description
Claims
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN03824376.8A CN1688254B (zh) | 2002-10-18 | 2003-08-28 | X线断层摄影装置 |
US10/524,341 US7684539B2 (en) | 2002-10-18 | 2003-08-28 | X-ray tomograph |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002304463 | 2002-10-18 | ||
JP2002-304463 | 2002-10-18 | ||
JP2003078125A JP4360817B2 (ja) | 2002-10-18 | 2003-03-20 | 放射線断層撮影装置 |
JP2003-78125 | 2003-03-20 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2004034908A1 true WO2004034908A1 (ja) | 2004-04-29 |
Family
ID=32109481
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2003/010971 WO2004034908A1 (ja) | 2002-10-18 | 2003-08-28 | X線断層撮影装置 |
Country Status (4)
Country | Link |
---|---|
US (1) | US7684539B2 (ja) |
JP (1) | JP4360817B2 (ja) |
CN (1) | CN1688254B (ja) |
WO (1) | WO2004034908A1 (ja) |
Families Citing this family (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6917663B2 (en) * | 2003-06-16 | 2005-07-12 | Kabushiki Kaisha Toshiba | Cone-beam reconstruction apparatus and computed tomography apparatus |
JP4785441B2 (ja) * | 2005-06-23 | 2011-10-05 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置 |
JP4509971B2 (ja) * | 2006-06-09 | 2010-07-21 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置 |
CN101529472B (zh) * | 2006-10-18 | 2012-12-05 | 皇家飞利浦电子股份有限公司 | 用于对感兴趣区域成像的成像系统、成像方法和计算机程序 |
US8116426B2 (en) * | 2008-11-11 | 2012-02-14 | Kabushiki Kaisha Toshiba | Computed tomography device and method using circular-pixel position-adaptive interpolation |
US8159493B2 (en) | 2009-02-03 | 2012-04-17 | Xerox Corporation | Adaptive grand tour |
US9025848B2 (en) * | 2010-06-17 | 2015-05-05 | Hitachi Medical Corporation | X-ray CT apparatus and control method for reconstructing X-ray CT images from filtered projection data |
JP5601683B2 (ja) * | 2010-06-28 | 2014-10-08 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 画像生成装置およびプログラム並びにx線ct装置 |
JP5728304B2 (ja) * | 2011-06-21 | 2015-06-03 | 株式会社日立メディコ | X線ct装置及び画像再構成方法 |
JP5968316B2 (ja) * | 2011-07-08 | 2016-08-10 | 株式会社日立製作所 | 画像再構成装置及び画像再構成方法 |
CN104039232B (zh) * | 2012-04-24 | 2016-09-28 | 株式会社日立制作所 | X射线ct装置及图像重构方法 |
US10652444B2 (en) | 2012-10-30 | 2020-05-12 | California Institute Of Technology | Multiplexed Fourier ptychography imaging systems and methods |
US9864184B2 (en) | 2012-10-30 | 2018-01-09 | California Institute Of Technology | Embedded pupil function recovery for fourier ptychographic imaging devices |
CN108761752A (zh) | 2012-10-30 | 2018-11-06 | 加州理工学院 | 傅立叶重叠关联成像系统、设备和方法 |
CN103961122B (zh) * | 2013-01-31 | 2018-07-31 | 通用电气公司 | 用于非等γ角CT系统中数据转换的方法和装置 |
CA2919985A1 (en) | 2013-07-31 | 2015-02-05 | California Institute Of Technology | Aperture scanning fourier ptychographic imaging |
WO2015027188A1 (en) | 2013-08-22 | 2015-02-26 | California Institute Of Technoloby | Variable-illumination fourier ptychographic imaging devices, systems, and methods |
CN106029171B (zh) * | 2014-02-24 | 2019-01-22 | 国立研究开发法人量子科学技术研究开发机构 | 放射线治疗用运动物体跟踪装置、放射线治疗用照射区域决定装置以及放射线治疗装置 |
US11468557B2 (en) | 2014-03-13 | 2022-10-11 | California Institute Of Technology | Free orientation fourier camera |
CN107111118B (zh) | 2014-12-22 | 2019-12-10 | 加州理工学院 | 用于厚样本的epi照明傅立叶重叠关联成像 |
CN107209362B (zh) * | 2015-01-21 | 2020-11-06 | 加州理工学院 | 傅立叶重叠关联断层摄影 |
US9829695B2 (en) | 2015-01-26 | 2017-11-28 | California Institute Of Technology | Array level Fourier ptychographic imaging |
CN107106112B (zh) * | 2015-02-12 | 2020-06-19 | 株式会社日立制作所 | X射线ct装置、图像处理装置及图像重构方法 |
AU2016233588A1 (en) | 2015-03-13 | 2017-10-26 | California Institute Of Technology | Correcting for aberrations in incoherent imaging system using fourier ptychographic techniques |
JP6654102B2 (ja) * | 2016-05-31 | 2020-02-26 | 株式会社日立製作所 | 粒子線治療システム |
US10568507B2 (en) | 2016-06-10 | 2020-02-25 | California Institute Of Technology | Pupil ptychography methods and systems |
US11092795B2 (en) | 2016-06-10 | 2021-08-17 | California Institute Of Technology | Systems and methods for coded-aperture-based correction of aberration obtained from Fourier ptychography |
WO2019090149A1 (en) | 2017-11-03 | 2019-05-09 | California Institute Of Technology | Parallel digital imaging acquisition and restoration methods and systems |
EP4111179A4 (en) * | 2020-02-27 | 2023-11-15 | Shenzhen Xpectvision Technology Co., Ltd. | IMAGING SYSTEM |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08187240A (ja) * | 1995-01-11 | 1996-07-23 | Toshiba Corp | コンピュータ断層撮影装置 |
JPH10290798A (ja) * | 1997-04-17 | 1998-11-04 | Ge Yokogawa Medical Syst Ltd | 投影データ測定方法および装置並びにx線ct装置 |
JPH114823A (ja) * | 1997-06-17 | 1999-01-12 | Toshiba Corp | X線コンピュータ断層撮影装置 |
JPH11253434A (ja) * | 1998-01-13 | 1999-09-21 | Koninkl Philips Electronics Nv | 断層x線撮像方法及び装置 |
JP2000225114A (ja) * | 1999-02-04 | 2000-08-15 | Toshiba Corp | X線ct装置 |
JP2002291732A (ja) * | 2001-04-04 | 2002-10-08 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
JP2003026080A (ja) * | 2001-07-10 | 2003-01-29 | Kayaba Ind Co Ltd | バウランプ装置 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5047931A (en) * | 1989-05-30 | 1991-09-10 | General Electric Company | Parallel computation of fan-beam back-projection reconstruction algorithm in computed tomography |
US5390112A (en) * | 1993-10-04 | 1995-02-14 | General Electric Company | Three-dimensional computerized tomography scanning method and system for imaging large objects with smaller area detectors |
JP2914891B2 (ja) * | 1995-07-05 | 1999-07-05 | 株式会社東芝 | X線コンピュータ断層撮影装置 |
CN1107290C (zh) * | 1996-05-02 | 2003-04-30 | 西门子公司 | 计算机层析x射线摄影机的图像再现方法 |
US6947584B1 (en) * | 1998-08-25 | 2005-09-20 | General Electric Company | Volume imaging system |
US6097784A (en) * | 1998-09-30 | 2000-08-01 | Picker International, Inc. | 3D image reconstruction for helical partial cone beam data |
JP2002095655A (ja) * | 2000-09-26 | 2002-04-02 | Shimadzu Corp | Ct装置 |
US6678346B2 (en) * | 2001-10-11 | 2004-01-13 | Ge Medical Systems Global Technology Company Llc | Cone-beam CT scanner with image reconstruction using multiple sub-images |
US6490333B1 (en) * | 2001-12-28 | 2002-12-03 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for cone-tilted parallel sampling and reconstruction |
US20030128801A1 (en) * | 2002-01-07 | 2003-07-10 | Multi-Dimensional Imaging, Inc. | Multi-modality apparatus for dynamic anatomical, physiological and molecular imaging |
-
2003
- 2003-03-20 JP JP2003078125A patent/JP4360817B2/ja not_active Expired - Lifetime
- 2003-08-28 US US10/524,341 patent/US7684539B2/en active Active
- 2003-08-28 WO PCT/JP2003/010971 patent/WO2004034908A1/ja active Application Filing
- 2003-08-28 CN CN03824376.8A patent/CN1688254B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08187240A (ja) * | 1995-01-11 | 1996-07-23 | Toshiba Corp | コンピュータ断層撮影装置 |
JPH10290798A (ja) * | 1997-04-17 | 1998-11-04 | Ge Yokogawa Medical Syst Ltd | 投影データ測定方法および装置並びにx線ct装置 |
JPH114823A (ja) * | 1997-06-17 | 1999-01-12 | Toshiba Corp | X線コンピュータ断層撮影装置 |
JPH11253434A (ja) * | 1998-01-13 | 1999-09-21 | Koninkl Philips Electronics Nv | 断層x線撮像方法及び装置 |
JP2000225114A (ja) * | 1999-02-04 | 2000-08-15 | Toshiba Corp | X線ct装置 |
JP2002291732A (ja) * | 2001-04-04 | 2002-10-08 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
JP2003026080A (ja) * | 2001-07-10 | 2003-01-29 | Kayaba Ind Co Ltd | バウランプ装置 |
Also Published As
Publication number | Publication date |
---|---|
US20060165211A1 (en) | 2006-07-27 |
CN1688254A (zh) | 2005-10-26 |
US7684539B2 (en) | 2010-03-23 |
CN1688254B (zh) | 2010-09-29 |
JP2004188163A (ja) | 2004-07-08 |
JP4360817B2 (ja) | 2009-11-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2004034908A1 (ja) | X線断層撮影装置 | |
US7903779B2 (en) | Apparatus and method for reconstruction of volumetric images in a divergent scanning computed tomography system | |
US9662084B2 (en) | Method and apparatus for iteratively reconstructing tomographic images from electrocardiographic-gated projection data | |
JP5133690B2 (ja) | ボクセルに依存する補間を用いる画像再構成 | |
US6285733B1 (en) | Computed tomography method utilizing a conical radiation beam | |
JP2003199740A (ja) | X線ct装置、x線ctの3次元再構成方法、及びx線ctの重み設定方法 | |
US9805448B2 (en) | Method and apparatus for computed tomography using asymetric filter for volume half reconstruction | |
JP2007098025A (ja) | 画像再構成方法、およびx線ct装置 | |
JP2000107168A (ja) | 被検体内の関心領域の3次元コンピュ―タ断層撮影方法及び装置 | |
JP2000081318A (ja) | 3次元コンピュ―タトモグラフィイメ―ジングのためのスキャニングおよびデ―タ収集方法およびイメ―ジング装置 | |
US6426989B2 (en) | Computed tomography method | |
JP5641811B2 (ja) | コンピュータ断層撮影装置 | |
JP2001095794A (ja) | 斜めのスライスを使用するコンピュータ断層撮影法において画像を再構成するための装置および方法 | |
JP2008501465A (ja) | 螺旋相対運動及び円錐ビームを用いたコンピュータ断層撮影方法 | |
JP3980696B2 (ja) | 画像再構成処理装置 | |
JP4025530B2 (ja) | X線ct装置 | |
US20240153160A1 (en) | Method for use in ct reconstruction | |
US6542572B2 (en) | Computed tomography method involving a helical relative motion | |
JP2006187453A (ja) | X線ct装置 | |
JP2011502683A (ja) | 関心領域の画像を決定するイメージング装置、イメージング方法及びコンピュータプログラム | |
JP2004113271A (ja) | 断層撮影装置 | |
US7292717B2 (en) | Computer tomography method using redundant measured values | |
JP3950101B2 (ja) | X線コンピュータ断層撮影装置 | |
JP4639331B2 (ja) | コーンビームct装置 | |
CN117979905A (zh) | 扫描仪和图像重建的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AK | Designated states |
Kind code of ref document: A1 Designated state(s): CN US |
|
AL | Designated countries for regional patents |
Kind code of ref document: A1 Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LU MC NL PT RO SE SI SK TR |
|
DFPE | Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101) | ||
121 | Ep: the epo has been informed by wipo that ep was designated in this application | ||
ENP | Entry into the national phase |
Ref document number: 2006165211 Country of ref document: US Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 10524341 Country of ref document: US |
|
WWE | Wipo information: entry into national phase |
Ref document number: 20038243768 Country of ref document: CN |
|
122 | Ep: pct application non-entry in european phase | ||
WWP | Wipo information: published in national office |
Ref document number: 10524341 Country of ref document: US |