WO2012077694A1 - X線ct装置及び画像再構成方法 - Google Patents
X線ct装置及び画像再構成方法 Download PDFInfo
- Publication number
- WO2012077694A1 WO2012077694A1 PCT/JP2011/078234 JP2011078234W WO2012077694A1 WO 2012077694 A1 WO2012077694 A1 WO 2012077694A1 JP 2011078234 W JP2011078234 W JP 2011078234W WO 2012077694 A1 WO2012077694 A1 WO 2012077694A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- data
- ray
- target data
- window function
- reconstruction
- Prior art date
Links
- 238000003168 reconstitution method Methods 0.000 title 1
- 238000012545 processing Methods 0.000 claims abstract description 61
- 238000006243 chemical reaction Methods 0.000 claims abstract description 20
- 238000000034 method Methods 0.000 claims description 120
- 230000006870 function Effects 0.000 claims description 113
- 230000008569 process Effects 0.000 claims description 43
- 238000012804 iterative process Methods 0.000 claims description 7
- 238000003384 imaging method Methods 0.000 abstract description 19
- 238000003780 insertion Methods 0.000 abstract description 4
- 238000004364 calculation method Methods 0.000 description 55
- 238000012966 insertion method Methods 0.000 description 25
- 238000010586 diagram Methods 0.000 description 9
- 238000005070 sampling Methods 0.000 description 9
- 238000005520 cutting process Methods 0.000 description 6
- 238000003860 storage Methods 0.000 description 5
- 238000012937 correction Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 230000037431 insertion Effects 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 238000003672 processing method Methods 0.000 description 3
- 230000002238 attenuated effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000007480 spreading Effects 0.000 description 2
- 238000003892 spreading Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 210000003027 ear inner Anatomy 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 210000003128 head Anatomy 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
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
-
- 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/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- 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/4078—Fan-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 irradiates a subject with X-rays, measures X-rays transmitted through the subject with an X-ray detector, and reconstructs projection data from multiple directions to obtain a tomographic image of the subject.
- the present invention relates to an image reconstruction process.
- the filtered back projection method is used from the viewpoint of calculation time, calculation accuracy, and memory capacity used.
- the successive approximation method repeatedly performs forward projection processing and backprojection processing, and thus has problems such as noise reduction and artifact reduction although the calculation time and the amount of memory used are problematic. Therefore, in recent years, studies have been made on shortening the calculation time and reducing the amount of memory used for practical application of the successive approximation method.
- a filter-corrected back projection method in multi-slice CT an extended Feldkamp method or a method applying this method is mainly used.
- a reconstruction filter such as a ramp filter or a Shep and Logan filter is applied to projection data, and then a value is embedded (integrated) in an image by back projection processing.
- a ray bundle driven type ray-driven
- a pixel driven type pixel-driven voxel-driven
- a distance driven type distance-driven
- the distance driving type is a method that considers the distance between the pixel boundary and the beam boundary as a reference.
- the distance-driven backprojection process as shown in FIG. 16, the distance between the pixel boundary and the beam boundary is scanned, and the projection values are sequentially embedded in the pixels 104 included in the beam 102.
- the distance-driven back projection process is disclosed in Patent Document 1, for example.
- the filter-corrected back projection method it is possible to increase the speed by using an image generated by the filtered back projection method as an initial image. Therefore, it is desirable that the filter-corrected back projection method can be used in combination even when image reconstruction is performed by the successive approximation method.
- the back projection process described above is also performed in an iterative process of the successive approximation method. Therefore, when the filter-corrected backprojection method and the successive approximation method are used in combination, a tomographic image having no difference in backprojection processing can be obtained by adopting the backprojection processing of the same method, and the development cost can be reduced. it can.
- a 1/4 channel offset (referred to as “quarter offset”) is made in the channel direction of the detector, and a beam of target phase data and counter phase data is obtained. The route is shifted. As a result, the sampling density in the channel direction of the beam can be effectively improved.
- back projection processing may be performed from the closest beam including the facing data.
- high resolution reconstruction high resolution reconstruction
- High-resolution reconstruction is a technique that is useful when photographing the head (particularly when diagnosing a microtissue such as the inner ear).
- JP 2005-522304 gazette JP-A-6-181919 Japanese Patent Laid-Open No. 10-314161
- the processing speed can be increased.
- the present invention has been made in view of the above-described problems, and an object of the present invention is to perform distance-driven back projection processing capable of high-resolution reconstruction and generate a high-resolution tomographic image. Is to provide line CT equipment, etc.
- the present invention is an image reconstruction method executed by an X-ray CT apparatus based on parallel beam data obtained by performing fan-para conversion on fan beam data.
- the parallel beam data of the phase to be the target data, the parallel beam data of the phase opposite to the target data is the counter data, and the presence / absence of the counter data corresponding to the target data is determined.
- Different high-resolution reconstruction is performed for each, and using the target data or the target data and the opposite data, back projection processing is performed based on the distance between the pixel boundary and the beam boundary, and image reconstruction is performed. It is characterized by performing.
- an X-ray CT apparatus and the like that can execute a distance-driven backprojection process capable of high-resolution reconstruction and generate a high-resolution tomographic image at high speed.
- the X-ray CT apparatus of the present invention is an X-ray CT apparatus that performs image reconstruction on the basis of parallel beam data obtained by performing fan-para conversion on fan beam data, the phase of the processing target
- a determination unit configured to determine parallel beam data as target data, the parallel beam data having a phase opposite to the target data as counter data, and determining whether or not the counter data corresponding to the target data exists; and determination by the determination unit Different high-resolution reconstruction is performed for each result, and using the target data or the target data and the opposite data, back projection processing is executed based on the distance between the pixel boundary and the beam boundary, and image reconstruction is performed.
- the image reconstruction unit further includes a determination unit that determines a window function indicating a ratio of the projection value corresponding to the target data or the target data and the opposite data to a pixel value.
- the back projection process is executed by applying the shape of the window function determined by the determination unit.
- the determination unit determines the shape of the window function related to one of the target data and the facing data and the window function related to the other.
- the respective window functions are determined so that the vertically inverted shapes are point-symmetric with each other.
- the determining unit sets the shape of the window function related to the target data and the counter data to a rectangle.
- the determination unit can select a triangle, a trapezoid, and a rectangle as the shape of the window function related to the target data and the facing data.
- the image reconstruction unit generates an initial image of the iterative process in the successive approximation method by a filtered back projection method, and executes the iterative process in the successive approximation method based on the initial image.
- pixel values are weighted according to the distance between the pixel boundary and the beam boundary.
- An image reconstruction method of the present invention is an image reconstruction method executed by an X-ray CT apparatus based on parallel beam data obtained by performing fan-para conversion on fan beam data, and is a phase to be processed And determining the presence or absence of the counter data corresponding to the target data by using the parallel beam data of the target data as the target data, the parallel beam data having a phase opposite to the target data as the counter data, and the determination step.
- a different high resolution reconstruction is performed for each determination result, and using the target data or the target data and the opposite data, back projection processing is performed based on a distance between a pixel boundary and a beam boundary, and an image And an image reconstruction step for performing reconstruction.
- the X-ray CT apparatus 1 is, for example, a multi-slice CT apparatus.
- the scan method is, for example, a rotate-rotate method (third generation).
- the X-ray CT apparatus 1 includes a scanner 2, an operation unit 3, and a bed 7.
- the scanner 2 scans the subject 4 placed on the bed 7 according to the instruction from the operation unit 3.
- Scanner 2 X-ray generator 5, collimator 6, detector 8, central controller 11, X-ray controller 12, scanner controller 13, high-voltage generator 14, collimator controller 15, bed controller 16, bed
- the moving measuring device 17, the driving device 18, the preamplifier 19, the A / D converter 20, and the like are included.
- the operation unit 3 includes an input / output device 31, an arithmetic device 32, and the like.
- the input / output device 31 includes a display device 33, an input device 34, a storage device 35, and the like.
- the arithmetic device 32 includes a reconstruction arithmetic device 36, an image processing device 37, and the like.
- the input device 34 in the operation unit 3 is configured by a mouse, a keyboard, a touch panel, and the like, and inputs measurement parameters and reconstruction parameters such as bed movement speed information and a reconstruction position.
- the display device 33 is configured by a display device such as a liquid crystal display.
- the storage device 35 is configured by a drive device of a hard disk or various storage media.
- the user uses the input device 34 in the operation unit 3 to capture the imaging conditions (bed movement speed, tube current, tube voltage, slice position, etc.) and reconstruction conditions (reconstruction method, high-resolution processing ON / OFF, image slice thickness, back projection) Phase width, region of interest, reconstructed image matrix size, reconstruction filter function, maximum number of iterations of successive approximation, convergence condition, etc.).
- imaging conditions bed movement speed, tube current, tube voltage, slice position, etc.
- reconstruction conditions reconstruction method, high-resolution processing ON / OFF, image slice thickness, back projection
- Phase width region of interest
- reconstructed image matrix size reconstructed image matrix size
- reconstruction filter function maximum number of iterations of successive approximation, convergence condition, etc.
- the central control device 11 Based on the input instruction, the central control device 11 sends control signals necessary for imaging to the X-ray control device 12, the scanner control device 13, and the bed control device 16, and starts imaging in response to the imaging start signal. To do.
- a control signal is sent to the high voltage generator 14 by the X-ray controller 12, a high voltage is applied to the X-ray generator 5, and the X-ray 9 is sent from the X-ray generator 5 to the subject 4 Is irradiated.
- a control signal is sent from the scanner control device 13 to the drive device 18, and the gantry including the X-ray generation device 5, the detector 8, the preamplifier 19, etc. circulates around the subject 4.
- the bed 7 on which the subject 4 is placed moves according to a control signal from the bed control device 16 (stationary scan) or translates in the body axis direction of the subject 4 (spiral scan).
- the X-ray 9 has its irradiation area limited by the collimator 6, is absorbed (attenuated) in each tissue in the subject 4, passes through the subject 4, and is detected by the detector 8.
- the X-ray 9 detected by the detector 8 is converted to a current, amplified by a preamplifier 19, converted to digital data by an A / D converter 20, converted to LOG, calibrated and calculated as a projection data signal Input to device 32.
- the projection data signal input to the computing device 32 becomes input data for image reconstruction processing performed by the reconstruction computing device 36.
- the reconstructed image is stored in the storage device 35 and displayed as a CT image by the display device 33.
- the image is displayed as a CT image by the display device 33.
- a tube voltage and a tube current are applied to the X-ray tube that is the X-ray generation device 5 based on the imaging conditions input from the input device 34 of the operation unit 3 connected to the scanner 2.
- the X-ray CT apparatus electrons having energy corresponding to the applied tube voltage are emitted from the cathode, and the emitted electrons collide with the target (anode), so that X-rays having energy corresponding to the electron energy are emitted.
- 9 is irradiated from the X-ray source of the X-ray tube.
- the irradiated X-ray 9 is arranged at a position facing the X-ray source through the subject 4 and attenuated according to the X-ray attenuation coefficient of the substance (tissue) in the subject 4 that passes through the subject 4.
- the light is received by the detector 8 to obtain projection data.
- the reconstruction calculation device 36 of the X-ray CT apparatus 1 obtains the filter-corrected projection data by superimposing the reconstruction filter on the projection data, and applies the filter-corrected projection data to the imaging conditions.
- X-rays inside the subject 4 are determined by performing back projection (image reconstruction) while applying weights (hereinafter referred to as “view weights”) of substantially the same shape in the view direction regardless of the position of the tomographic image.
- view weights weights
- a tomogram is imaged non-destructively as a distribution map of attenuation coefficients.
- the detector 8 of the X-ray CT apparatus 1 is a one-dimensional detector (also referred to as a “single-row detector” or “single slice”) that is arranged in a one-dimensional direction in the circulation direction for the purpose of imaging a wide range in a short time. . ) May be employed in the direction of the rotation axis (also referred to as “multi-row detector” or “multi-slice detector”).
- the X-ray CT apparatus 1 in which the detector 8 is arranged in a one-dimensional manner in the circumferential direction is called “single slice CT”, and the X-ray CT apparatus 1 arranged in two dimensions is called “multi-slice CT”.
- an X-ray beam spreading in a fan shape is emitted from X-ray generator 5 (X-ray source), and in multi-slice CT, cone-shaped from X-ray generator 5 (X-ray source) according to detector 8 Alternatively, an X-ray beam spreading in a pyramid shape is irradiated.
- X-ray irradiation is performed while circling around the subject 4 placed on the bed 7.
- imaging in which the bed 7 is fixed and the X-ray generator 5 (X-ray source) circulates around the subject 4 in a circular orbit is called “normal scan” or “axial scan”.
- imaging in which the bed 7 moves and the X-ray generator 5 (X-ray source) orbits around the subject 4 in a spiral trajectory is called “helical scan” or “helical scan”.
- the X-ray CT apparatus 1 In the case of the successive approximation method, it is advantageous for the X-ray CT apparatus 1 to use an image generated by the filtered back projection method as an initial image for speeding up. In the iterative processing in the successive approximation method, the X-ray CT apparatus 1 is superior to the application of distance-driven back projection processing for reasons such as high-frequency errors. By performing high-resolution reconstruction based on distance-driven backprojection processing, the same type of backprojection processing can be applied in the iterative processing of the filter-corrected backprojection method and the successive approximation method. As a result, a tomographic image having no difference with respect to the backprojection process can be obtained, and the development cost can be reduced.
- the fan beam has a fan-like spread when viewed from the circumferential axis direction.
- one beam is treated as a triangle in the distance-driven backprojection process.
- the shapes of the target beam and the counter beam are compared, since the shapes of the two are different as shown in FIG. 2, it is not possible to perform high-resolution reconstruction using data related to the counter beam in addition to the target beam.
- the parallel beams are parallel as viewed from the direction of the rotation axis and are equally spaced.
- one beam is treated as a rectangular shape in distance-driven backprojection processing.
- the shapes of the two match as shown in FIG. 3, so the high resolution reconstruction is performed using the data related to the counter beam by the “0 insertion method” described later. Can do.
- the reconstruction calculation device 36 performs fan-para conversion, which is conversion processing from fan beam data to parallel beam data, and performs image reconstruction processing based on the parallel beam data after fan-para conversion.
- the fan-para conversion process may use a known technique.
- parallel beam data having a phase to be processed is set as target data
- parallel beam data having a phase opposite to the target data is set as counter data
- the reconstruction calculation device 36 of the X-ray CT apparatus 1 uses both the filter-corrected back projection method and the successive approximation method, and both employ distance-driven back projection processing.
- the reconstruction calculation device 36 executes a first back projection process shown in FIG.
- the reconstruction calculation device 36 executes the first backprojection process shown in FIG. 4 or the second backprojection process shown in FIG.
- the reconstruction calculation device 36 repeatedly executes the back projection process.
- the reconstruction calculation device 36 determines the presence or absence of opposing data corresponding to each target data based on the imaging conditions, and calculates from where to where the opposing data exists (step 11). That is, the reconstruction calculation device 36 calculates a backprojection phase width that can be used to generate one slice of the reconstructed image from the imaging conditions.
- the reconstruction calculation device 36 calculates the backprojection phase width from the parallel beam data as follows in the reconstruction in which the X-ray beam inclination in the body axis direction is taken into account as in the Feldkamp reconstruction. To do.
- F is the backprojection phase width
- ⁇ d row is the row direction detector element size
- N row is the number of detector rows
- SID is the distance between the source and circuit center
- T is the bed moving speed
- SDD is the source-detection.
- FOV effective field size
- x 0 is reconstruction center position in X direction
- y 0 is reconstruction center position in Y direction
- FOV x is FOV size in X direction
- FOV y is FOV size in Y direction
- the backprojection phase width is calculated from the imaging conditions using the above formula, but the present invention is not limited to this, and the backprojection phase width calculated in advance may be used as the reconstruction condition. You may set directly from. Further, assuming that the projection data is extended by extrapolation in the detector row direction, it may be set larger than the value calculated from the above equation.
- FIG. 6 shows the relationship between the backprojection phase width and the presence or absence of opposing data.
- Reference numerals 41a, 41b, 41c, and 41d indicate back projection phase widths, respectively.
- Reference numerals 42a and 42b denote phases in which no opposing data exists.
- FIG. 6A shows a case where the backprojection phase width is 180 degrees
- FIG. 6B shows a case where the backprojection phase width is 180 degrees to 360 degrees
- FIG. 6C shows a case where the backprojection phase width is 360 degrees.
- FIG. 6 (d) shows an example in which the backprojection phase width is 360 degrees or more.
- the reconstruction calculation device 36 performs a high-resolution transformation of the 0 insertion method (see FIG. 8) described later for the phase range where the opposite data exists, and the data described later for the phase range where the opposite data does not exist. Perform high-resolution conversion using an interpolation method (see FIG. 9). Then, the reconstruction calculation device 36 creates high-resolution projection data sampled twice (step 12). In the case of single-slice CT, high-resolution conversion may be performed by a counter-insertion method (see FIG. 7), which will be described later, for a phase range in which counter-data exists.
- the reconstruction calculation device 36 calculates view weights (step 13), and performs filter correction processing on the high-resolution projection data. Then, the reconstruction calculation device 36 performs distance-driven backprojection processing (see FIG. 11), which will be described later, using the filter-corrected high-resolution projection data subjected to the filter correction, and generates a reconstructed image (step 14).
- the reconstruction calculation device 36 determines the presence / absence of opposing data corresponding to each target data based on the imaging conditions, and calculates the phase range in which the opposing data exists (step 11).
- the first backprojection processing method shown in FIG. 4 is performed in the backprojection processing at the time of iterative processing when an image is generated by the successive approximation reconstruction method after determining the presence or absence of the opposing data (step 11).
- the configuration calculation device 36 performs high-resolution conversion using the zero insertion method for the phase range where the opposite data of the parallel beam projection data exists, and uses the data interpolation method for the phase range where the opposite data does not exist. Performs high-resolution conversion to create double-sampled parallel high-resolution projection data (step 12), calculates view weight (step 13), and generates parallel-high-resolution projection data by distance-driven backprojection (Step 14).
- the first other backprojection processing method shown in FIG. 4 is the backprojection at the time of iterative processing in the case of generating an image by the successive approximation reconstruction method after determining the presence / absence of opposing data (step 11).
- the reconstruction calculation unit 36 performs high-resolution conversion using the data interpolation method both when the parallel data of the parallel beam projection data exists and when it is not (step 12), obtains parallel high-resolution projection data, and obtains the view weight. Is calculated (step 13), and distance-driven backprojection processing (step 14) is performed.
- the reconstruction calculation device 36 determines the presence / absence of opposing data corresponding to each target data based on the imaging conditions, and calculates the phase range in which the opposing data exists (step 21).
- the second backprojection processing method shown in FIG. 5 is performed by performing backprojection processing at the time of iterative processing when an image is generated by the successive approximation reconstruction method after determining the presence or absence of the opposing data (step 21).
- the configuration calculation unit 36 calculates the view weight (step 22) and opposes the window size of the detector element with respect to one of the pixel and the detector element to the phase range where the opposite data of the parallel beam projection data exists.
- An image is generated by performing distance-driven high-resolution backprojection processing that is set to half the size when there is no data (step 23).
- high-resolution conversion may be performed by a counter-insertion method (see FIG. 7), which will be described later, for a phase range in which counter-data exists.
- an initial image that is an input image of the iterative process is preferably generated by the above-described filter-corrected back projection process. As a result, the calculation time of the iterative process can be shortened.
- the three methods are the opposite insertion method, the zero insertion method, and the data interpolation method.
- an offset corresponding to 1/4 channel (“quarter offset”) is provided in the channel direction of the detector 8.
- phase width capable of backprojection is gathered over 360 degrees and the opposite data exists in any phase.
- phase width that can be backprojected is not gathered 360 degrees and there is only a part of the opposing data, a process for obtaining an image with as high a resolution as possible is executed.
- Opposite insertion method prepares projection data for 360 degrees, for example, and embeds (inserts) opposing data in the target data, thereby doubling the number of channels and reducing the channel spacing to half. This is a method of generating projection data for 180 degrees.
- the X-ray CT apparatus 1 performs the filter correction process on the assumption that the double-sampled projection data obtained by the opposite insertion method is obtained by a virtual detector having two times the number of channels. Obtain corrected projection data. Then, the X-ray CT apparatus 1 embeds (backprojects) 180 degree projection data into the pixel while generating a beam passing through the pixel center based on the filter-corrected projection data by interpolation between the nearest beams.
- a high-resolution image can be obtained.
- sampling is effectively doubled, and interpolation processing can be performed between closer projection data including opposing data, so that a high-resolution image can be generated.
- the data can be sampled twice in the channel direction before the filter correction process, and a high resolution can be achieved by using a reconstruction filter with better high frequency characteristics.
- the reconstruction calculation device 36 inserts opposite data between channels corresponding to the target data (hereinafter referred to as “opposite insertion”), thereby sampling in the channel direction. Double the density and then back project 180 degrees. Since the opposite insertion method requires opposite data, it can be applied only to the phase range where the opposite data exists. Therefore, it is necessary to use a relatively slow helical pitch in the case of helical scanning. Since the back projection may be 180 degrees, the reconstruction process can be speeded up.
- the opposing insertion method can be used only in single slice CT.
- the reconstruction calculation device 36 performs back projection processing without using view weights.
- the reconstruction computation device 36 inserts 0 data (hereinafter referred to as “0 insertion”) instead of inserting opposing data between channels as shown in FIG.
- 0 insertion 0 data
- the sampling in the channel direction is densified, and the same backprojection phase width as that in the normal reconstruction process is backprojected. Since the zero insertion method requires opposite data, it can be applied only to the phase range where the opposite data exists.
- the 0 insertion method is different from the opposite insertion method and simply inserts 0 between channels, so it can be applied to projection data obtained by multi-slice CT.
- the view weight is set to 1 during backprojection (“1” is equivalent to the case where the view weight is not used), or the view weight is used. Need to be reconfigured.
- a known view weight can be used. By using view weights, motion artifacts due to subject movement and helical artifacts due to helical scanning can be reduced.
- the data interpolation method is a method of sampling twice by preparing projection data for 360 degrees, for example, and embedding data created by interpolation from projection data of a target phase.
- the X-ray CT apparatus 1 generates projection data for 360 degrees in which the number of channels is doubled and the channel spacing is halved by data interpolation, and is obtained by a virtual detector with twice the number of channels. It is assumed that a high-resolution image can be obtained by embedding (back-projecting) 360-degree projection data into a pixel while generating a beam passing through the pixel center by interpolation between closest beams.
- the reconstruction calculation device 36 performs sampling between channels within the same phase to double the sampling, and backprojects the backprojection phase width.
- the data interpolation method can be applied even when there is no opposite data, and in the case of a helical scan, a relatively high helical pitch can be used.
- the data interpolation method since the data interpolation method only inserts interpolated data between channels as in the zero insertion method, it can also be applied to projection data obtained by multi-slice CT. Regardless, known view weights can be used.
- the spatial resolution is inferior compared to the opposite insertion method and the zero insertion method.
- FIG. 10 shows the trajectory of the X-ray 9 with respect to the image plane 51 showing the set of pixels when the subject 4 is viewed from the side. Specifically, the trajectory until the X-ray 9 starts from the X-ray generator 5, passes through the image plane 51 and reaches the detector 8 is shown.
- the pixel boundary 52 indicates the boundary position between pixels (M1, M2, M3, M4) adjacent in the X-axis direction on the X-Z plane.
- the beam boundary 53 indicates the boundary position of the virtual X-ray beam.
- the beam boundary 53 coincides with the boundary position of detector elements (D1, D2, D3) adjacent in the X-axis direction in the X-Z plane among the detector elements included in the detector 8.
- the beam boundary 53 is a virtual boundary position for dividing the X-ray beam bundle for each detector element for convenience of processing.
- L det is the interval between the pixel boundaries 52, that is, the size of the detector element in the X-axis direction. Usually, L det is constant for all detector elements.
- L1 indicates the width in the X-axis direction where the X-ray beam corresponding to the detector element D1 and the pixel M1 overlap on the X-Z plane.
- L2 indicates the width in the X-axis direction in which the X-ray beam corresponding to the detector element D1 and the pixel M2 overlap in the X-Z plane.
- the window function is generally a function that becomes 0 outside a certain finite interval.
- a window function is used to cut a certain shape into a finite space in an infinitely continuous space. In the embodiment of the present invention, it means the ratio (contribution rate) that the projection contributes to the pixel, or vice versa.
- the reconstruction calculation device 36 determines a window function indicating the ratio of the projection value corresponding to the target data to the pixel value in the first back projection process. In addition, the reconstruction calculation device 36 determines a window function indicating the ratio of the projection value corresponding to the target data and the opposing data to the pixel value in the second backprojection process.
- FIG. 11 illustrates the window function in the first backprojection process.
- the shape of the window function 54a shown in FIG. 11 is a rectangle, and is continuous with the adjacent window function 54a. That is, the horizontal width of the window function 54a is equal to the horizontal width of the beam boundary 53.
- the reconstruction calculation device 36 integrates the projection value of the X-ray beam included in the pixel boundary 52.
- Raw 1 to raw 5 are projection values 55 for each detector element.
- sum 1 and sum 2 are pixel values 56 for each pixel.
- L 1 indicates the width in the X-axis direction where the window function 54a for cutting out the projection value 55 “raw 1 ” and the pixel corresponding to the pixel value 56 “sum 1 ” overlap on the XZ plane.
- L 2 indicates the width in the X-axis direction where the window function 54a for cutting out the projection value 55 and the pixel corresponding to the pixel value 56 “sum 1 ” overlap in the XZ plane.
- L 3 indicates the width in the X-axis direction where the window function 54a for cutting out the projection value 55 “raw 3 ” and the pixel corresponding to the pixel value 56 “sum 1 ” overlap on the XZ plane.
- L 4 indicates the width in the X-axis direction where the window function 54a for cutting out the projection value 55 “raw 3 ” and the pixel corresponding to the pixel value 56 “sum 2 ” overlap on the XZ plane.
- L 5 indicates the width in the X-axis direction in which the window function 54a for cutting out the projection value 55 “raw 4 ” and the pixel corresponding to the pixel value 56 “sum 2 ” overlap on the XZ plane.
- L 6 indicates the width in the X-axis direction where the window function 54a for cutting out the projection value 55 “raw 5 ” and the pixel corresponding to the pixel value 56 “sum 2 ” overlap on the XZ plane.
- the reconstruction calculation device 36 forms the shape of the window function 54a shown in FIG. 11 for the pixels included in the target data based on the distance between the pixel boundary and the beam boundary.
- back projection processing is executed and image reconstruction is performed.
- the back projection process for the pixel value 56 “sum 1 ” is executed by the following expression.
- the back projection process for the pixel value 56 “sum 2 ” is executed by the following expression.
- the reconstruction calculation device 36 may use the area ratio as a weight instead of the width ratio in the X-axis direction.
- the backprojection process for the pixel value 56 “sum 1 ” may be executed by the following equation.
- the back projection process for the pixel value 56 “sum 2 ” may be executed by
- FIG. 12 illustrates a first example of the window function in the second backprojection process.
- the window function 54b shown in FIG. 12 relates to the target data
- the window function 54c relates to the counter data.
- the shapes of the window functions 54b and 54c are rectangular and have a width that is half the width of the window function 54a shown in FIG. Accordingly, the window functions 54b and 54c are not continuous with the adjacent window functions 54b and 54c, respectively.
- the window function 54b and the window function 54c have the same shape.
- the reconstruction calculation device 36 takes the shape of the window functions 54b and 54c as a rectangle having a half width of the window function 54a, and positions the window functions 54b and 54c at the center of the X-ray beam. (The positions of the window functions 54b and 54c are determined so that the midpoint of the window functions 54b and 54c in the width direction coincides with the center of the X-ray beam). Then, the reconstruction calculation device 36 multiplies the projection value 55 as a weight by the area ratio of the window functions 54b and 54c included in the pixel boundary 52 to the whole pixel, and further multiplies the projection weight 55 by the view weight. Integration processing of projection data is performed on the target pixel. The calculation method using the area ratio as a weight is the same as the calculation method described in the example shown in FIG.
- FIG. 13 illustrates a second example of the window function in the second backprojection process.
- the window function 54d shown in FIG. 13 relates to the target data
- the window function 54e relates to the counter data.
- the shape of the window functions 54d and 54e is an isosceles trapezoid and has an area that is half the area of the window function 54a shown in FIG. Accordingly, the window functions 54d and 54e are not continuous with the adjacent window functions 54d and 54e, respectively.
- the window function 54d and the window function 54e have the same shape.
- the reconstruction calculation device 36 makes the shape of the window functions 54d and 54e an isosceles trapezoid having an area that is half the area of the window function 54a, and sets the positions of the window functions 54d and 54e to the X-ray beam. (The positions of the window functions 54d and 54e are determined so that the midpoint in the width direction of the window functions 54d and 54e and the center of the X-ray beam coincide with each other). Then, the reconstruction calculation device 36 multiplies the projection value as a weight by the area ratio of the window functions 54d and 54e included in the pixel boundary 52 to the entire pixel, and further multiplies the projection weight 55 by the view weight. Integration processing of projection data is performed on the pixels. The calculation method using the area ratio as a weight is the same as the calculation method described in the example shown in FIG.
- FIG. 14 illustrates a third example of the window function in the second backprojection process.
- the window function 54f shown in FIG. 14 relates to the target data
- the window function 54g relates to the counter data.
- the shape of the window functions 54f and 54g is an isosceles triangle, and has an area half the area of the window function 54a shown in FIG. Accordingly, the window functions 54f and 54g are continuous only at both apexes of the base and the adjacent window functions 54f and 54g, respectively.
- the window function 54f and the window function 54g have the same shape.
- the reconstruction calculation device 36 sets the shape of the window functions 54f and 54g to an isosceles triangle having an area that is half the area of the window function 54a, and sets the positions of the window functions 54f and 54g to the X-ray beam. (The positions of the window functions 54f and 54g are determined so that the midpoint of the bottom sides of the window functions 54f and 54g and the center of the X-ray beam coincide with each other). Then, the reconstruction calculation unit 36 multiplies the projection value as a weight by the area ratio of the window functions 54f and 54g included in the pixel boundary 52 with respect to the entire pixel, and further multiplies the projection value 55 by the view weight. Integration processing of projection data is performed on the pixels. The calculation method using the area ratio as a weight is the same as the calculation method described in the example shown in FIG.
- the window functions are collectively referred to, the reference numeral “54” is given.
- the adjacent window functions 54 are continuous at all points of the beam boundary 53, and the embedding of the counter data is not considered.
- the adjacent window function 54 is discontinuous at the beam boundary 53 or continuous at only one point, and the embedding of the opposing data is considered. Yes. Accordingly, the examples shown in FIGS. 12 to 14 (second backprojection processing) enable high-resolution reconstruction using the facing data even in the distance-driven backprojection processing.
- the target data window function 54b and the counter data window function 54c which are rectangular alternately exist, and there is no position where both exist. This means that there is no pixel on which back projection processing is performed using both target data and counter data. Therefore, in the example shown in FIG. 12, the spatial resolution can be improved as compared with the examples in FIGS.
- both the window function 54f of the target data that is isosceles trapezoid and the inclined portion of 54g of the opposing data (near the end of the X-ray beam) Exists, and only one of the other portions exists, that is, near the center of the X-ray beam. Since the image noise can be reduced in the part that uses both the target data and the counter data, the example shown in FIG. 13 improves the spatial resolution near the center of the X-ray beam and near the end of the X-ray beam. Image noise can be reduced.
- both the window function 54f of the target data that is an isosceles triangle and 54g of the opposing data always exist. Since the image noise can be reduced in the portion using both the target data and the counter data, the image noise can be reduced in the example shown in FIG. 14 as compared with the examples in FIG. 13 and FIG.
- the image noise and the spatial resolution are in a trade-off relationship, the image noise and the spatial resolution may be adjusted according to the purpose of diagnosis and imaging conditions. Therefore, the X-ray CT apparatus 1 accepts an instruction to improve the spatial resolution or reduce the image noise via the input device. Then, the reconstruction calculation device 36 appropriately selects a triangle, a trapezoid, and a rectangle as the shape of the window function 54 based on an instruction input via the input device 34.
- the reconstruction calculation device 36 selects a rectangle as the shape of the window function 54.
- the reconstruction calculation device 36 selects a triangle as the shape of the window function 54.
- the reconstruction calculation device 36 selects a trapezoid as the shape of the window function 54.
- FIGS. 12 to 14 are common in that the shape of the window function 54 related to one of the target data and the opposite data and the shape obtained by vertically inverting the window function 54 related to the other are point-symmetric with each other. There is a nature of.
- FIG. 15 (a) shows the window function 54b and the window function 54c of the example shown in FIG.
- the shape of the window function 54b and the shape obtained by vertically inverting the window function 54c (the shape of the window function 54c is a rectangle, so there is no change in the shape even if it is turned upside down). , They are point-symmetric with each other.
- FIG. 15B shows the window function 54d and the window function 54e of the example shown in FIG. As shown in FIG. 15 (b), the shape of the window function 54d and the shape obtained by vertically inverting the window function 54e are point-symmetric with each other.
- FIG. 15 (c) shows the window function 54f and the window function 54g of the example shown in FIG. As shown in FIG. 15 (c), the shape of the window function 54f and the shape obtained by vertically inverting the window function 54g are point-symmetric with each other.
- the reconstruction calculation device 36 performs the back projection process by applying the window function 54 having the above-mentioned common properties, so that the opposite data can be obtained even in the distance-driven back projection process. Enables high-resolution reconstruction.
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
Description
距離駆動型は、画素境界とビーム境界との間の距離を基準として考える手法である。距離駆動型の逆投影処理では、図16に示すように、画素境界とビーム境界との間の距離を走査していき、ビーム102内に含まれる画素104に投影値を順次埋め込んでいく。距離駆動型の逆投影処理は、例えば、特許文献1において開示されている。
本発明のX線CT装置は、ファンビームデータに対してファンパラ変換を行うことによって得られるパラレルビームデータを基に画像再構成を行うX線CT装置であって、処理対象とする位相の前記パラレルビームデータを対象データとし、前記対象データと対向する位相の前記パラレルビームデータを対向データとし、前記対象データに対応する前記対向データの存在の有無を判定する判定部と、前記判定部による判定結果ごとに異なるハイレゾ再構成を行い、前記対象データ、又は、前記対象データ及び前記対向データを用いて、画素境界とビーム境界との間の距離に基づいて逆投影処理を実行し、画像再構成を行う画像再構成部と、を具備する。
本発明の画像再構成方法は、ファンビームデータに対してファンパラ変換を行うことによって得られるパラレルビームデータを基にX線CT装置が実行する画像再構成方法であって、処理対象とする位相の前記パラレルビームデータを対象データとし、前記対象データと対向する位相の前記パラレルビームデータを対向データとし、前記対象データに対応する前記対向データの存在の有無を判定する判定ステップと、前記判定ステップによる判定結果ごとに異なるハイレゾ再構成を行い、前記対象データ、又は、前記対象データ及び前記対向データを用いて、画素境界とビーム境界との間の距離に基づいて逆投影処理を実行し、画像再構成を行う画像再構成ステップと、を含む。
なお、以下の説明及び添付図面において、同一の機能を有する構成要素については、同一の符号を付することにより重複説明を省略することにする。
X線CT装置1は、例えば、マルチスライスCT装置である。スキャン方式は、例えば、ローテート-ローテート方式(第3世代)である。X線CT装置1は、スキャナ2と操作ユニット3と寝台7とによって構成される。
記憶装置35は、ハードディスクや各種の記憶媒体のドライブ装置によって構成される。
X線CT装置1では、スキャナ2に接続される操作ユニット3の入力装置34から入力される撮影条件に基づき、X線発生装置5であるX線管に管電圧、管電流が印加される。
)を周回軸方向に拡張した2次元検出器(「多列検出器」、「マルチスライス検出器」とも言う。)を採用しても良い。一般に、検出器8が周回方向の1次元に配置されるX線CT装置1は「シングルスライスCT」、2次元に配置されるX線CT装置1は「マルチスライスCT」と呼ばれる。シングルスライスCTでは、X線発生装置5(X線源)から扇状に広がるX線ビームが照射され、マルチスライスCTでは、検出器8に合わせてX線発生装置5(X線源)から円錐状、もしくは角錐状に広がるX線ビームが照射される。
対向挿入法は、例えば360度分の投影データを準備し、対象データの中に対向データを埋め込む(挿入する)ことによって、チャネル数が2倍となり、チャネル間隔が半分である180度分の投影データを生成する方法である。X線CT装置1は、対向挿入法により得られる2倍サンプリング投影データを、チャネル数が2倍の仮想的な検出器によって得られたものと仮定して、フィルタ補正処理を行うことによって、フィルタ補正投影データを得る。そして、X線CT装置1は、フィルタ補正投影データを基に画素中心を通過するビームを最近接ビーム間の補間により生成しながら、180度分の投影データを画素に埋め込む(逆投影する)ことによって高分解能画像を得ることができる。この方法を用いる場合、サンプリングが実効的に2倍となり、対向データを含めて、より近い投影データ間での補間処理が可能であることから、高分解能な画像を生成できる。フィルタ補正逆投影法では、フィルタ補正処理の前にチャネル方向にデータを2倍サンプリング化しておき、より高周波特性にすぐれた再構成フィルタを用いることによって、高分解能化を行うことができる。
0挿入法では、再構成演算装置36は、図8に示すように、チャネル間に対向データを挿入する代わりに0データを挿入(以下、「0挿入」という。)することにより、チャネル方向のサンプリングを高密度化し、通常の再構成処理と同じ逆投影位相幅分を逆投影する。0挿入法は、対向データが必要であるため、対向データが存在する位相範囲にのみ適用できる。
データ補間法は、例えば360度分の投影データを準備し、対象とする位相の投影データから補間により作成したデータを埋め込むことによって2倍サンプリング化する方法である。X線CT装置1は、データ補間法によってチャネル数が2倍かつチャネル間隔が半分である360度分の投影データを生成し、チャネル数が2倍の仮想的な検出器によって得られたものと仮定して、画素中心を通過するビームを最近接ビーム間の補間により生成しながら、360度分の投影データを画素に埋め込む(逆投影する)ことによって高分解能画像を得ることができる。
検出器素子D1の投影値=(M1の画素値×L1+M2の画素値×L2)/Ldet
の関係が成り立つと言える。
窓関数とは、一般には、ある有限区間以外において0となる関数である。窓関数は、無限に継続している空間の中で、ある有限の空間に、ある特定の形状によって切り出すために用いられる。本発明の実施の形態では、投影が画素に対して寄与する割合(寄与率)、もしくは、その逆を意味する。
sum1=(raw1×L1+raw2×L2+raw3×L3)/Limg
の式によって、画素値56「sum1」に関する逆投影処理を実行する。
sum2=(raw3×L4+raw4×L5+raw5×L6)/Limg
の式によって、画素値56「sum2」に関する逆投影処理を実行する。
sum1=(raw1×S1+raw2×S2+raw3×S3)/(S1+S2+S3)
の式によって、画素値56「sum1」に関する逆投影処理を実行しても良い。
sum2=(raw3×S4+raw4×S5+raw5×S6)/(S4+S5+S6)
の式によって、画素値56「sum2」に関する逆投影処理を実行しても良い。
図11に示す例(第1の逆投影処理)は、隣接する窓関数54がビーム境界53の全ての点において連続であり、対向データの埋め込みが考慮されていない。一方、図12~図14に示す例(第2の逆投影処理)は、隣接する窓関数54がビーム境界53において不連続、または1点のみにおいて連続であり、対向データの埋め込みが考慮されている。従って、図12~図14に示す例(第2の逆投影処理)は、距離駆動型の逆投影処理においても、対向データを用いるハイレゾ再構成が可能となる。
Claims (8)
- ファンビームデータに対してファンパラ変換を行うことによって得られるパラレルビームデータを基に画像再構成を行うX線CT装置であって、
処理対象とする位相の前記パラレルビームデータを対象データとし、前記対象データと対向する位相の前記パラレルビームデータを対向データとし、前記対象データに対応する前記対向データの存在の有無を判定する判定部と、
前記判定部による判定結果ごとに異なるハイレゾ再構成を行い、前記対象データ、又は、前記対象データ及び前記対向データを用いて、画素境界とビーム境界との間の距離に基づいて逆投影処理を実行し、画像再構成を行う画像再構成部と、
を具備するX線CT装置。 - 前記対象データ、又は、前記対象データ及び前記対向データに対応する投影値が画素値に対して寄与する割合を示す窓関数を決定する決定部、
を更に具備し、
前記画像再構成部は、前記決定部によって決定される前記窓関数の形状を適用して逆投影処理を実行する請求項1に記載のX線CT装置。 - 前記決定部は、前記判定部によって前記対向データが存在すると判定された場合には、前記対象データ及び前記対向データのいずれか一方に関する前記窓関数の形状と、もう一方に関する前記窓関数を上下反転させた形状が、互いに点対称となるように、それぞれの前記窓関数を決定する請求項2に記載のX線CT装置。
- 前記決定部は、前記対象データ及び前記対向データに関する前記窓関数の形状を矩形とする請求項3に記載のX線CT装置。
- 前記決定部は、前記対象データ及び前記対向データに関する前記窓関数の形状として、三角形、台形及び矩形を選択可能である請求項3に記載のX線CT装置。
- 前記画像再構成部は、逐次近似法における反復処理の初期画像をフィルタ補正逆投影法によって生成し、前記初期画像に基づいて逐次近似法における反復処理を実行する請求項3に記載のX線CT装置。
- 前記逆投影処理は、前記画素境界とビーム境界との間の距離に応じて画素値の重み付けを行う請求項1に記載のX線CT装置。
- ファンビームデータに対してファンパラ変換を行うことによって得られるパラレルビームデータを基にX線CT装置が実行する画像再構成方法であって、
処理対象とする位相の前記パラレルビームデータを対象データとし、前記対象データと対向する位相の前記パラレルビームデータを対向データとし、前記対象データに対応する前記対向データの存在の有無を判定する判定ステップと、
前記判定ステップによる判定結果ごとに異なるハイレゾ再構成を行い、前記対象データ、又は、前記対象データ及び前記対向データを用いて、画素境界とビーム境界との間の距離に基づいて逆投影処理を実行し、画像再構成を行う画像再構成ステップと、
を含む画像再構成方法。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/989,564 US8971607B2 (en) | 2010-12-10 | 2011-12-07 | X-ray CT apparatus and image reconstruction method |
JP2012547882A JP5858928B2 (ja) | 2010-12-10 | 2011-12-07 | X線ct装置及び画像再構成方法 |
CN201180059648.0A CN103260520B (zh) | 2010-12-10 | 2011-12-07 | X射线ct装置以及图像重构方法 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010-275278 | 2010-12-10 | ||
JP2010275278 | 2010-12-10 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2012077694A1 true WO2012077694A1 (ja) | 2012-06-14 |
Family
ID=46207181
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2011/078234 WO2012077694A1 (ja) | 2010-12-10 | 2011-12-07 | X線ct装置及び画像再構成方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US8971607B2 (ja) |
JP (1) | JP5858928B2 (ja) |
CN (1) | CN103260520B (ja) |
WO (1) | WO2012077694A1 (ja) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014011681A2 (en) * | 2012-07-09 | 2014-01-16 | The Trustees Of The Universtiy Of Pennsylvania | Super-resolution tomosynthesis imaging systems and methods |
WO2014034797A1 (ja) * | 2012-08-31 | 2014-03-06 | 株式会社 東芝 | X線ct装置 |
WO2015016328A1 (ja) * | 2013-07-31 | 2015-02-05 | 株式会社東芝 | X線CT(Computed Tomography)装置、画像処理装置、画像処理方法及び記憶媒体 |
CN104902818A (zh) * | 2013-02-05 | 2015-09-09 | 株式会社日立医疗器械 | X射线ct装置以及图像重构方法 |
WO2016017402A1 (ja) * | 2014-07-30 | 2016-02-04 | 株式会社 日立メディコ | データ処理方法、データ処理装置、及びx線ct装置 |
JP2019005013A (ja) * | 2017-06-22 | 2019-01-17 | コニカミノルタ株式会社 | 解析装置及び解析システム |
US11185294B2 (en) | 2016-02-26 | 2021-11-30 | The Trustees Of The University Of Pennsylvania | Super-resolution tomosynthesis imaging systems and methods |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6218334B2 (ja) * | 2012-11-30 | 2017-10-25 | 株式会社日立製作所 | X線ct装置及びその断層画像撮影方法 |
WO2014192831A1 (ja) * | 2013-05-28 | 2014-12-04 | 株式会社東芝 | 医用画像診断装置及び制御方法 |
JP6521575B2 (ja) * | 2013-06-11 | 2019-05-29 | キヤノンメディカルシステムズ株式会社 | X線コンピュータ断層撮影装置 |
WO2015053787A1 (en) * | 2013-10-11 | 2015-04-16 | Analogic Corporation | Tomosynthesis imaging |
US10241181B2 (en) | 2014-01-13 | 2019-03-26 | Siemens Healthcare Gmbh | Resolution enhancement of diffusion imaging biomarkers in magnetic resonance imaging |
US10410343B2 (en) * | 2015-07-08 | 2019-09-10 | Hitachi, Ltd. | Image computing device, image computing method, and tomograph |
US10115211B2 (en) * | 2016-03-25 | 2018-10-30 | L3 Security & Detection Systems, Inc. | Systems and methods for reconstructing projection images from computed tomography volumes |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH01181847A (ja) * | 1988-01-11 | 1989-07-19 | Toshiba Corp | 画像再構成方法 |
JPH06181919A (ja) * | 1992-12-22 | 1994-07-05 | Yokogawa Medical Syst Ltd | ヘリカルスキャンにおけるスキューアーティファクトの低減方法 |
JPH10286253A (ja) * | 1997-02-17 | 1998-10-27 | Toshiba Corp | X線ct装置 |
JPH10314161A (ja) * | 1997-05-15 | 1998-12-02 | Ge Yokogawa Medical Syst Ltd | 倍密度補間方法およびx線ct装置 |
JP2005224637A (ja) * | 2005-05-12 | 2005-08-25 | Toshiba Corp | X線ct装置及びそのミスアライメント補正方法 |
JP2009089810A (ja) * | 2007-10-05 | 2009-04-30 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1021128A1 (en) * | 1997-10-10 | 2000-07-26 | Analogic Corporation | Computed tomography scanning target detection |
US20050175143A1 (en) * | 2002-06-03 | 2005-08-11 | Osamu Miyazaki | Multi-slice x-ray ct device |
JP2006345979A (ja) * | 2005-06-14 | 2006-12-28 | Ge Medical Systems Global Technology Co Llc | データ処理方法、逆投影方法、再投影方法、画像再構成方法、および放射線ct装置 |
CN101404088B (zh) * | 2008-11-05 | 2011-10-26 | 华中科技大学 | Ct图像重建的方法及系统 |
-
2011
- 2011-12-07 US US13/989,564 patent/US8971607B2/en not_active Expired - Fee Related
- 2011-12-07 WO PCT/JP2011/078234 patent/WO2012077694A1/ja active Application Filing
- 2011-12-07 JP JP2012547882A patent/JP5858928B2/ja not_active Expired - Fee Related
- 2011-12-07 CN CN201180059648.0A patent/CN103260520B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH01181847A (ja) * | 1988-01-11 | 1989-07-19 | Toshiba Corp | 画像再構成方法 |
JPH06181919A (ja) * | 1992-12-22 | 1994-07-05 | Yokogawa Medical Syst Ltd | ヘリカルスキャンにおけるスキューアーティファクトの低減方法 |
JPH10286253A (ja) * | 1997-02-17 | 1998-10-27 | Toshiba Corp | X線ct装置 |
JPH10314161A (ja) * | 1997-05-15 | 1998-12-02 | Ge Yokogawa Medical Syst Ltd | 倍密度補間方法およびx線ct装置 |
JP2005224637A (ja) * | 2005-05-12 | 2005-08-25 | Toshiba Corp | X線ct装置及びそのミスアライメント補正方法 |
JP2009089810A (ja) * | 2007-10-05 | 2009-04-30 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014011681A2 (en) * | 2012-07-09 | 2014-01-16 | The Trustees Of The Universtiy Of Pennsylvania | Super-resolution tomosynthesis imaging systems and methods |
WO2014011681A3 (en) * | 2012-07-09 | 2014-03-06 | The Trustees Of The Universtiy Of Pennsylvania | Super-resolution tomosynthesis imaging systems and methods |
US9743891B2 (en) | 2012-07-09 | 2017-08-29 | The Trustees Of The University Of Pennsylvania | Super-resolution tomosynthesis imaging systems and methods |
US9737279B2 (en) | 2012-08-31 | 2017-08-22 | Toshiba Medical Systems Corporation | X-ray CT apparatus |
JP2014061278A (ja) * | 2012-08-31 | 2014-04-10 | Toshiba Corp | X線ct装置 |
WO2014034797A1 (ja) * | 2012-08-31 | 2014-03-06 | 株式会社 東芝 | X線ct装置 |
CN104902818A (zh) * | 2013-02-05 | 2015-09-09 | 株式会社日立医疗器械 | X射线ct装置以及图像重构方法 |
WO2015016328A1 (ja) * | 2013-07-31 | 2015-02-05 | 株式会社東芝 | X線CT(Computed Tomography)装置、画像処理装置、画像処理方法及び記憶媒体 |
JP2015029913A (ja) * | 2013-07-31 | 2015-02-16 | 株式会社東芝 | X線CT(ComputedTomography)装置、画像処理装置、画像処理方法及び記憶媒体 |
US9224216B2 (en) | 2013-07-31 | 2015-12-29 | Kabushiki Kaisha Toshiba | High density forward projector for spatial resolution improvement for medical imaging systems including computed tomography |
WO2016017402A1 (ja) * | 2014-07-30 | 2016-02-04 | 株式会社 日立メディコ | データ処理方法、データ処理装置、及びx線ct装置 |
US11185294B2 (en) | 2016-02-26 | 2021-11-30 | The Trustees Of The University Of Pennsylvania | Super-resolution tomosynthesis imaging systems and methods |
US11786186B2 (en) | 2016-02-26 | 2023-10-17 | The Trustees Of The University Of Pennsylvania | Super-resolution tomosynthesis imaging systems and methods |
JP2019005013A (ja) * | 2017-06-22 | 2019-01-17 | コニカミノルタ株式会社 | 解析装置及び解析システム |
Also Published As
Publication number | Publication date |
---|---|
CN103260520A (zh) | 2013-08-21 |
US8971607B2 (en) | 2015-03-03 |
US20130243299A1 (en) | 2013-09-19 |
CN103260520B (zh) | 2015-06-17 |
JPWO2012077694A1 (ja) | 2014-05-19 |
JP5858928B2 (ja) | 2016-02-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5858928B2 (ja) | X線ct装置及び画像再構成方法 | |
JP4646810B2 (ja) | 断層撮影像の再構成方法及び断層撮影装置 | |
US6665370B2 (en) | Computed tomography method and apparatus for acquiring images dependent on a time curve of a periodic motion of the subject | |
JP5960048B2 (ja) | 再構成演算装置、再構成演算方法、及びx線ct装置 | |
JP4644785B2 (ja) | コーンビームct画像再構成におけるアーチファクトを低減するための方法及び装置 | |
JP5280450B2 (ja) | X線ct画像形成方法及びそれを用いたx線ct装置 | |
US9123156B2 (en) | X-ray CT apparatus and image reconstruction method | |
US7227923B2 (en) | Method and system for CT imaging using a distributed X-ray source and interpolation based reconstruction | |
JP6571313B2 (ja) | 医用画像診断装置及び制御方法 | |
JP5406063B2 (ja) | 再構成演算装置、再構成演算方法、及びx線ct装置 | |
JP5728304B2 (ja) | X線ct装置及び画像再構成方法 | |
JP4612294B2 (ja) | X線コンピュータ断層装置、x線コンピュータ断層装置制御方法、及びx線コンピュータ断層撮影プログラム | |
JP2007000408A (ja) | X線ct装置 | |
US9592021B2 (en) | X-ray CT device, and method | |
US9406121B2 (en) | X-ray CT apparatus and image reconstruction method | |
JP2004313655A (ja) | 放射線計算断層画像装置および断層画像生成方法 | |
WO2016132880A1 (ja) | 演算装置、x線ct装置、及び画像再構成方法 | |
CN107341836B (zh) | 一种ct螺旋扫描图像重建方法及装置 | |
JP2005137752A (ja) | X線ct装置 | |
JP2024037308A (ja) | 医用画像処理装置および医用画像処理方法 | |
US20060243914A1 (en) | Attenuation map generation from pet scans | |
JPH0620449B2 (ja) | X線断層撮影装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 11847152 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2012547882 Country of ref document: JP Kind code of ref document: A |
|
WWE | Wipo information: entry into national phase |
Ref document number: 13989564 Country of ref document: US |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 11847152 Country of ref document: EP Kind code of ref document: A1 |