WO2005077278A1 - 断層撮影像の再構成方法及び断層撮影装置 - Google Patents
断層撮影像の再構成方法及び断層撮影装置 Download PDFInfo
- Publication number
- WO2005077278A1 WO2005077278A1 PCT/JP2005/002264 JP2005002264W WO2005077278A1 WO 2005077278 A1 WO2005077278 A1 WO 2005077278A1 JP 2005002264 W JP2005002264 W JP 2005002264W WO 2005077278 A1 WO2005077278 A1 WO 2005077278A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- reconstructing
- width
- tomographic image
- data
- detector
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims description 105
- 238000012937 correction Methods 0.000 claims abstract description 76
- 238000003384 imaging method Methods 0.000 claims description 38
- 238000012545 processing Methods 0.000 claims description 36
- 230000005855 radiation Effects 0.000 claims description 17
- 238000003325 tomography Methods 0.000 claims description 10
- 230000003247 decreasing effect Effects 0.000 claims description 2
- 230000005540 biological transmission Effects 0.000 claims 1
- 230000002123 temporal effect Effects 0.000 claims 1
- 230000006870 function Effects 0.000 description 85
- 238000010586 diagram Methods 0.000 description 17
- 230000008569 process Effects 0.000 description 16
- 230000008707 rearrangement Effects 0.000 description 9
- 230000007423 decrease Effects 0.000 description 8
- 238000001514 detection method Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 6
- 238000004088 simulation Methods 0.000 description 5
- 238000013459 approach Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000006866 deterioration Effects 0.000 description 2
- 238000004904 shortening Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000010894 electron beam technology Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 230000015654 memory Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 238000007669 thermal treatment Methods 0.000 description 1
- 230000004584 weight gain Effects 0.000 description 1
- 235000019786 weight gain Nutrition 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/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/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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/412—Dynamic
-
- 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 method for reconstructing a tomographic image from projection data in a computer tomographic apparatus (hereinafter, referred to as a CT apparatus) using a fan beam or a cone beam, and in particular,
- the present invention relates to a method for reconstructing a tomographic image from projection data having a backprojection phase width of 7T [ rad ] or more, and a CT apparatus for realizing the method.
- An X-ray CT apparatus receives X-rays emitted from an X-ray source and transmitted through a subject with an X-ray detector arranged at a position facing the X-ray source to acquire projection data.
- the X-ray source and detector facing each other around the subject are rotated around the rotation axis, and projection data at different rotation angles (phases) of the X-ray source and detector are collected.
- the inside of the subject is non-destructively imaged.
- Such an X-ray CT apparatus uses a single-row detector in which detection elements are arranged one-dimensionally (linearly), and an X-ray CT apparatus and detection elements are two-dimensionally arranged.
- an X-ray CT device that uses a multi-row detector.
- the simplest X-ray CT imaging method is a normal scan method in which an X-ray source and a detector are rotated around a rotation axis in a range of 2 ⁇ to perform imaging, and the normal scan method is used.
- the scan range of the projected data is 2 ⁇ [rad].
- the X-ray source and the detector can be rotated during one rotation.
- the same projection data (line integral) will be measured twice.
- an imaging method in which the smaller the data redundancy is, the smaller the scan range is 2 ⁇ is also adopted.
- the beam when the position of the X-ray source and the position of one detection element are switched as shown in FIGS. 14 (a) and (b) is equivalent, and as shown in FIG. 14 (c), the fan beam Assuming that the maximum fan angle is 2 ym, when the X-ray source moves by ⁇ + 2 ⁇ m, all beams necessary for image reconstruction are projected.
- the shadow data can be measured. That is, this range is the minimum scan range.
- phase range of data that can be back-projected is different for each pixel. That is, as shown in Fig. 15, for example, as shown in Fig. 15, for the pixel pi, the data whose phase range is ⁇ or more around the pixel pi is used for reconstruction.For the pixel ⁇ 2, the phase range is ⁇ or less around the pixel ⁇ 2. This means that only the data of this type is used.
- FIG. 16 shows this as a sinogram showing projection data with the horizontal axis representing the fan angle (the angle between the center beam and each beam) and the vertical axis representing the orbiting phase angle ⁇ .
- FIG. 16 is a sinogram showing the minimum complete data set, and two areas indicated by hatched triangles above and below are mutually redundant data.
- the following document 1 proposes to assign a weight to a predetermined area of projection data.
- Patent Document 1 JP 2001-299738 A
- Ref. 1 proposes a weighting function applied to an intermediate data set.
- a virtual fan that does not depend on the actual physical maximum fan angle
- the shape of the weight function is trapezoidal, triangular, or nonlinearly deformed in the phase direction.
- the scan range approaches 2 ⁇ it will approach a triangle from a trapezoid. This means that as the scanning range approaches 2 ⁇ , the area where the weighting factor becomes 1 or less increases, and the data contribution rate decreases significantly compared to the scanning range 2 ⁇ where all the weighting factors are 1, and the noise Can be significantly increased (ie, the SNR decreases).
- an object of the present invention is to provide a reconstruction method using a weight function that can reduce body motion artifacts and can be applied to reconstruction of projection data in all scan ranges.
- any correction angle width (body motion correction range) and back projection phase width (reconstruction A weighting function is created based on the width in the view direction used for image reconstruction, and image reconstruction is performed using this weighting function.
- the image reconstruction method of the present invention is configured as follows.
- the radiation source and the detector which are arranged opposite to each other with the object in between, are circulated around a predetermined circling axis, and the detector detects the transmitted radiation emitted from the radiation source and transmitted through the object.
- the correction angle width and the backprojection phase width is set, and the setting of the correction angle width and the backprojection phase width is not set.
- the method may include a step of setting the other value based on the set one value.
- the correction angle width is ⁇ and the backprojection phase width is 2F ⁇ , 0 ⁇ ( 2F-1) (however, ⁇ ⁇ 2F_2 eeil (1 ° g2F ) ).
- ⁇ is called a correction angle width index
- F is called a back projection phase width index.
- the correction angle width is used to remove discontinuity between data due to body movement occurring between a set of projection data and projection data measured subsequently, and to correct data redundancy.
- the correction angle width index ⁇ expresses the guaranteed width of the slope portion of the weight function.
- the correction angle width is set, for example, corresponding to a range of an area for correcting data discontinuity at an end of the projection data (hereinafter, referred to as a data discontinuity area).
- the correction angle width may be changed according to the magnitude of the noise amount in the reconstructed image or the magnitude of the body motion artifact.
- the correction angle width can be reduced in direct proportion to the back projection phase width if.
- the back projection phase width (data width) is determined in consideration of data redundancy, SN, and the time width (time resolution) that contributes to the image, and the data width of the minimum complete data set [ ⁇ + the maximum value of the fan angle 2 times] or more.
- the weight function created based on the corrected angle width index ⁇ and the back projection phase width index F is such that the weight (coefficient of weight) in the data discontinuous region is equal to that in another region equivalent to the data discontinuous region.
- the weights are set to be smaller than the weights.For example, it is assumed that the first sub-weight function and the second sub-weight function obtained by shifting the second sub-weight function by a predetermined phase are added and normalized. be able to.
- the second sub-weight function is shifted by the phase width to obtain a second sub-weight function.
- the second sub-weight function is obtained by adding the first sub-weight function and the second sub-weight function and normalizing (multiplying by the sub-weight gain).
- a function having a trapezoidal shape with a top side of ⁇ _ ⁇ and a bottom side of ⁇ + ⁇ is given.
- discontinuity between data is corrected according to a user request indicated by a correction angle width index, and there is no image distortion due to redundancy. It will be.
- the scan range is ⁇ or more, it can be applied to projection data with an arbitrary backprojection phase width of 2 ⁇ or less and 2 ⁇ or more, and the problem of discontinuity in processing bounded by the scan range 2 ⁇ is solved. .
- the tomographic image reconstruction method of the present invention further includes a step of rearranging a fan beam emitted from the radiation source into a parallel beam, and a weight function for reconstructing the parallel beam.
- w ( ⁇ ) is the correction angle width ⁇ ⁇ [rad] and 2 (N — “ ⁇ F_ ⁇ / 2 ⁇ 2 ⁇ 0 or more, where ⁇ is the rotation phase (view phase) at the time of projection data detection. ⁇ obtained from the integer of
- Vl ⁇ -one F + 2 ( N- one)
- W2 (2 * (-F) + ⁇ ) * W1 / + W1 if [ ⁇ > 0, M ⁇ F]
- the above-mentioned weight function can be represented by the following equation (3) using the following equation (4) and equation (8).
- Ws is the sub-weight for the parallel beam
- ⁇ is the view phase
- ⁇ and ⁇ are the central view phases of the sub-weights
- ⁇ is the sub-weight reference width
- G is the sub-weight gain
- ⁇ is 2 ( N — D ⁇ F_ ⁇ / 2 ⁇ 2 ⁇ is an integer greater than or equal to 0.
- At least one of the correction angle width index ⁇ and the back projection width index F is input by the user, and the others are automatically determined.
- Equation (3) it is possible to modify Equation (3) so as to make the sub-weights non-linear, and use a non-linear weight function. That is, the nonlinear weight function is expressed by Expression (3 ′) using Expression (4), Expression (8), and Expression (9).
- the weight function w ( ⁇ ,) for reconstructing the fan beam when the projection phase of the fan beam is ⁇ and the fan angle is ⁇ is the correction angle width ⁇ ⁇ [rad] and 2 ( N_1) ⁇ F- ⁇ / 2 ⁇ 2 ⁇ ( ⁇ is an integer of 0 or more) [Ss
- W2 (2 * (MF) + £) * Wl / f + W1 if [ ⁇ > 0, M ⁇ F]
- the weight function we, ⁇ ) for the fan beam can be expressed by the equation do) by using the same expression method as the weight function wp (e) for the parallel beam described above. It can be represented by the used equation ( ⁇ '). Any of these may be used for the fan beam. ( Ten)
- the projection data is data detected while moving the object in the rotation axis direction along with the rotation of the radiation source and the detector.
- the method includes a step of interpolating the projection data and creating projection data of a plane orthogonal to the rotation axis (two-dimensional back projection processing).
- the tomography apparatus of the present invention provides a tomographic image of a region of interest of the object from a radiation source and a detector arranged opposite to each other with the object interposed therebetween and projection data detected by the detector.
- the reconstruction unit is equipped with the above-described tomographic image reconstruction method. It is characterized by the following.
- the imaging control means can change the back projection phase width depending on the imaging site.
- the imaging control unit increases the SNR by increasing the back projection phase width to perform imaging.
- imaging is performed by improving the time resolution by reducing the back projection phase width.
- the tomography apparatus of the present invention has means for moving the object in the orbital axis direction, and the imaging control means corresponds to the moving speed of the object in the orbital axis direction and adjusts the correction angle used for the reconstruction.
- the width and / or backprojection phase width can be changed.
- the detector may be a single-row detector or a multi-row detector.
- the reconstruction means may use the same weighting factor for each row of detectors, or use a different weighting factor for at least one detector row than the other rows.
- the tomography apparatus further includes an input unit that receives information from the user regarding the correction angle width and the back projection phase width.
- an input unit that receives information from the user regarding the correction angle width and the back projection phase width.
- the minimum correction angle width is maintained by setting the correction angle width index ⁇ to 0 or more ( ⁇ > 0), so that a rectangular shape is not obtained.
- the weight of the present invention has a two-step trapezoidal shape as shown in FIG.
- FIG. 1 is a diagram showing an outline of a CT apparatus for realizing the tomographic image reconstruction method of the present invention.
- the CT apparatus mainly includes a scanner 40, an operation unit 50, and a bed 60 on which a subject is mounted and moved.
- the scanner 40 includes a central controller 400, an X-ray controller 401, a high-voltage generator 402, a high-voltage switching unit 403, an X-ray generator 404, an X-ray detector 405, a preamplifier 406, a scanner controller 407, and a driving device. 408, a collimator control device 409, a bed control device 410, a bed movement measurement device 411, and the like.
- the operation unit 50 includes an input / output device 51 that also includes a display device, an input device, and a storage device, and an arithmetic device 52 that includes a reconstruction arithmetic device, an image processing device, and the like.
- the input device consists of a mouse, keyboard, etc., and inputs measurement and reconstruction parameters such as bed moving speed information and reconstruction position.
- measurement and reconstruction parameters such as bed moving speed information and reconstruction position.
- the user can set a correction angle width index ⁇ and a back projection phase width index F as necessary parameters, and these indexes are also input from the input device.
- the corrected angle width index ⁇ and the back projection phase width index F will be described later in detail.
- the storage device stores the information input by the input device, the processing result in the arithmetic device 52, and the like.
- the display device displays such information and various data such as reconstructed images.
- the reconstruction operation device processes data obtained from the X-ray detector, and the image processing device performs various processes on the reconstructed image and displays the processed image on a display device.
- the central processing control device 400 controls the imaging conditions (bed moving speed, tube current, tube voltage, slice position, etc.) and reconstruction parameters (region of interest, reconstructed image size, Based on the back projection phase width, reconstruction filter function, etc.), control signals necessary for imaging are transmitted to the X-ray controller 401, the bed controller 410, and the scanner controller 407. Start.
- a control signal is sent from the X-ray controller 401 to the high-voltage generator 402, and the high voltage passes through the high-voltage switching unit 403 and is applied to the X-ray generator 404.
- the emitted X-ray is applied to the subject, and the transmitted light is incident on the X-ray detector 405.
- a control signal is sent from the scanner control device 407 to the drive device 408, and the X-ray generator 404, the X-ray detector 405, and the preamplifier 406 are controlled to orbit around the subject.
- the couch 60 on which the subject is placed is stationary by the couch controller 410 during this orbit.
- the spiral orbit scan 42 as shown in FIG. 2B, the bed 60 is moved at a predetermined pitch in the direction of the orbital axis of the X-ray generator 404 or the like.
- the distance ⁇ that the bed moves relative to the imaging form during one rotation of the imaging system is the ratio (D / ⁇ ) to the width D of the detection element in the orbital direction.
- the beam pitch is defined as the helical pitch, and the ratio of the detector to the total length in the orbital direction is defined as the beam pitch.
- the helical pitch is used up to “2”, which can cover the entire imaging area almost completely in consideration of the facing data.
- the irradiation area of the X-rays emitted from the X-ray generator 404 is limited by the collimator 412 controlled by the collimator controller 409, and is absorbed (attenuated) by each tissue in the subject. While passing through the subject, it is detected by the X-ray detector 405.
- the X-rays detected by the X-ray detector 405 are converted into a current there, amplified by the preamplifier 406, and input to the arithmetic unit 52 of the operation unit 50 as a projection data signal.
- the projection data signal input to the arithmetic unit 52 is subjected to image reconstruction processing by a reconstruction arithmetic unit in the arithmetic unit 52. This reconstructed image is stored in a storage device in the input / output device 51, and is displayed on the display device as a CT image.
- the X-ray detector 405 is a single-array detector 11 in which detection elements are arranged one-dimensionally as shown in FIG. 3 (a), or a single-row detector as shown in FIG. 3 (b). May be any of the multi-row detectors 12 arranged in a plurality of rows in the direction of the rotation axis (the direction of the arrow in the figure). In Fig. 3 (a), the detector 11 is shown as a straight line. Typically, the detectors are set so that the distance from the X-ray source 10 to each detector element and the angle between adjacent X-ray beams are equal. Those arranged in an arc shape are used. In a single-row detector, the X-ray beam is orthogonal to the orbital axis. In the case of a multi-row detector, a larger area can be imaged at a time than a single-array detector. It has a tilt angle (cone angle).
- the processing executed by the reconstruction arithmetic device includes, for example, rearrangement processing for associating fan beam projection data obtained from a fan beam with parallel beam parallel data, fan beam projection data or parallel processing. Calculates the weighting function to be applied to the beam projection data and applies a weighting function to the projection data. Data correction processing, filter processing by superimposing a reconstruction filter on the parallel projection data. Filter correction to generate parallel beam projection data. Processing and filtering Back-projection processing for back-projecting parallel beam projection data to the back-projection area corresponding to the region of interest, and data for creating data such as a circular orbit by interpolating data obtained during helical trajectory scanning Interpolation processing.
- the detector 405 is a multi-row detector as shown in FIG. 3B, cone angle correction processing for further multiplying each row of projection data by a coefficient depending on the tilt angle of radiation is performed. Do.
- Figure 4 (a) shows the procedure of the reconstruction method.
- the shooting method is circular orbit scan, spiral scan Either one may be used, but the case of circular orbit scanning will be described here.
- the projection data is obtained by irradiating the subject with the X-rays emitted from the X-ray generator 404 and transmitting the X-rays to the X-ray detector 405 at the same time as the X-ray generator 404.
- X-ray detector 405 and preamplifier 406 orbit around the subject (step 101).
- the scan range is 2F TT determined by the back projection phase width index F input from the input device, and is ⁇ or more when considered as a parallel beam. That is, a numerical value satisfying F ⁇ 0.5 is input as the back projection phase width index F.
- the user determines an appropriate value in consideration of the amount of X-ray exposure, time resolution, noise, etc.
- the obtained fan beam projection data is rearranged into parallel beam projection data (Step 102).
- This rearrangement process is performed by combining data of different phases and different fan angles from a fan beam radiated radially as viewed from the orbital axis direction as shown in Fig. 5 (a).
- This is a process of converting the beam into a parallel beam as seen from the direction of the orbital axis as shown in FIG.
- Sl and S2 indicate the positions of the radiation source and the detector.
- the address is calculated by using a parallel beam, which calculates the address of the X-ray beam passing through the reconstructed pixel on the detector at each projection phase.
- the inverse trigonometric function that requires a long processing time and has a high computational load, and the distance calculation between the X-ray source and the reconstructed pixel are not required, and the inverse trigonometric function is replaced with the product-sum operation. Therefore, there is an advantage of speeding up despite the increase in the operation time required for the rearrangement process.
- a correction for removing the image distortion due to the body motion correction and the redundancy of the projection data is performed (step 103).
- a weighting function w ( ⁇ ) to be applied to the parallel beam is obtained based on the correction angle width index ⁇ and the back projection phase width index F input from the input device.
- the correction angle width index ⁇ is set by the user in consideration of the degree of body movement according to the imaging region of the subject in order to eliminate discontinuity in data due to body movement. Is set to ⁇ , the correction angle width is set so as not to exceed the [scan width (2F7i) _7r], that is, 0 ⁇ (2F-1).
- the weighting function w ( ⁇ ) ( ⁇ is the projection phase) is, for example, as shown in FIG. Equation (1) (for a fan beam) or Equation (2) (for a parallel beam) is satisfied for the determined correction area of the projection data,
- a trapezoidal sub-weighting function Ws ( ⁇ ) with a base of ⁇ + ⁇ and an upper side of ⁇ _ ⁇ is determined, and this sub-weighting function Ws ( ⁇ ) is determined by the scan width 2FTT for the center of the data width. It is obtained by adding and shifting the phase shifted (2 ⁇ _ ( ⁇ + ⁇ )). Normalization is a process of obtaining the power by averaging the weights of each phase.
- Vl ⁇ -F + 2 (N_1) if one F + 2 (N- ⁇ 0]
- W2 (2 * (M-F) + ⁇ ) * W1 / ⁇ + W1 if [ ⁇ > 0, ⁇ F]
- the weight function Wp for the parellel beam is expressed by the following equation (4) and the equation (8). ⁇
- G 2 ' N (8)
- Ws is the sub weight for the parallel beam
- ⁇ is the view phase
- ⁇ and ⁇ are the center view phases of the sub weights
- ⁇ is the sub weight reference width
- G is the sub weight gain
- ⁇ is 2 (N -.
- D ⁇ F_ ⁇ / 2 ⁇ is an integer of 0 or more to be 2
- New correction angle width index epsilon at least one of backprojection width index F is entered by the user, others are determined automatically.
- Equation (3) it is also possible to modify Equation (3) so as to make the sub-weights non-linear, and use a non-linear weight function. That is, the nonlinear weight function is expressed by Expression (3 ′) using Expression (4), Expression (8), and Expression (9).
- the same weight function can be used for all channel positions of the parallel beam (positions corresponding to the fan angle position of the fan beam), and the weight can be stored in a small number of memories.
- the CT reconstruction can include redundant data without using a weight having a limited phase width.
- the parallel beam projection data to which the weight function thus obtained is applied is subjected to reconstruction filter processing by superimposing a reconstruction filter function in the channel direction (step 104), and then performing back projection processing.
- Known methods can be used for the reconstruction filter and the back projection process.
- the CT apparatus using the single-row detector has been described.
- a plurality of CT apparatuses arranged in the rotation axis direction are used. Since projection data is obtained for the detection elements in a column, the above-described rearrangement process (step 102) and correction process using a weighting function (step 103) are performed on the projection data for each column.
- the parallel beam projection data is multiplied by a cosine of several cone angles corresponding to each column, and reconstruction filter processing for superimposing a reconstruction filter function in the channel direction of the parallel beam in each column is performed.
- a CT image can be obtained in the same manner as a CT device using a single-row detector.
- the weighting function applied to the projection data of each column may be the same or different.
- the correction angle width index can be set in consideration of the movement of the part to be imaged. For example, when photographing a plurality of parts by normal scanning, the back projection phase width 2FTT is widened for a part where the amount of noise is important, and narrow for a part where time resolution is important.
- the data redundancy differs according to the position of the pixel to be reconstructed, and if the correction angle width differs according to the redundancy in order to use data efficiently. You may let it. For example, the correction angle width is increased at a position where redundancy is high, and the correction angle width is reduced at a position where redundancy is low.
- the reconstruction method proposed by the present applicant Japanese Patent Laid-Open No. 2004-188163
- This method simplifies the conventional arcsin calculation by calculating an approximate straight line of a curve indicating the position of the radiation source with respect to the position of the parallel beam projection data in the channel direction.By adopting this method, Significant speedup can be realized compared with the conventional reconstruction method.
- projection data obtained by circular orbit scanning has been described as an example.
- the present invention can be similarly applied to spiral orbit scanning.
- the streak-like artefat at the corresponding position can be obtained only with the filter-corrected 2D back projection due to the discontinuity of the data at the imaging end phase. Occurs. Therefore, as shown in Fig. 2 (b), data obtained from the spiral orbit is subjected to data interpolation to correct the data into circular orbit data, and then a filter correction process is performed.
- Such data interpolation is performed prior to rearrangement into parallel beam projection data (step 102) in the reconstruction procedure shown in FIG.
- the projection data obtained by the spiral trajectory scan is interpolated in the moving direction of the bed to create data such as a circular trajectory.
- the projection data thus interpolated is subjected to the above-described rearrangement processing (step 102), correction processing using a weighting function (step 103), reconstruction filter processing (step 104), and back projection processing (step 105).
- CT images can be obtained in the same way as circular orbit data.
- the degree of the artifact generated at the imaging end in the spiral trajectory scan is determined by the degree of discontinuity in the movement trajectory of the X-ray source. That is, the degree of the artifact changes depending on the moving speed of the subject (the moving speed of the bed). Therefore, when determining the correction angle width index ⁇ and the back projection phase width index F, it is necessary to consider the bed moving speed in addition to the X-ray exposure dose, time resolution, noise, etc., which are considered for normal scan. is necessary. For example, when the moving speed is fast, decrease F. By decreasing F, the maximum value of ⁇ that can be set also decreases, but to obtain a correction effect, ⁇ It is desirable to make it as large as possible.
- Fig. 4 (b) shows an example of the procedure of the reconstruction method.
- body movement correction of the projection data without rearranging the obtained fan beam projection data and correction for removing image distortion due to redundancy are performed (step 112).
- the data correction processing is based on the correction angle width index ⁇ and the back projection phase width index F input from the input device, and a weight function w ( ⁇ , ⁇ ) applied to the fan beam Ask for.
- the weight function Wf ( ⁇ , ⁇ ) for the fan beam can be expressed by the expression do) by using the same expression method as the weight function Wp (e) for the parallel beam described above, and the sub-weighted nonlinear weight Can be represented by the following equation (10 ′).
- FIG. 10 shows an example of a weight function for fan beam back projection according to the present embodiment.
- the weight shape can be changed as shown in (c).
- the result is as shown in Fig. 7 (c), where ⁇ / 2 + ⁇ ⁇ F ⁇ M + ⁇ / If it becomes 2, it will be as shown in Fig. 7 (a).
- the fan beam projection data to which the weight function thus obtained is applied is subjected to reconstruction filter processing by superimposing a reconstruction filter function in the channel direction (step 113), and then performing an inverse trigonometric function operation and X-ray Backprojection processing involving the calculation of the distance between the source and the reconstructed pixels is performed to obtain a CT image (step 114).
- reconstruction filter processing by superimposing a reconstruction filter function in the channel direction (step 113), and then performing an inverse trigonometric function operation and X-ray Backprojection processing involving the calculation of the distance between the source and the reconstructed pixels is performed to obtain a CT image (step 114).
- Known methods can be used for the reconstruction filter and the back projection process.
- discontinuity at the data end can be corrected with a constant correction effect regardless of the value of the backprojection phase width.
- the rearrangement process from the fan beam to the parallel beam is not performed, and the inverse trigonometric function operation and the X-ray source Configuration It may be applied to a reconstruction method for calculating the distance between pixels, or after correcting fan beam projection data with a weight function, rearrangement processing to parallel beams, reconstruction filter processing, back projection processing May be applied.
- weighting of fan beam projection data is performed by cutting projection data in continuous scans and spiral trajectory scans. It can also be used as a bay window. Specifically, it can be applied to interventional imaging that performs surgery or thermal treatment while using contrast imaging or CT equipment as a monitor.
- the correction angle width and the back projection phase width may use values set as defaults in advance, or at least one of them may be set arbitrarily by the user. May be.
- a weighting function is obtained before imaging based on a correction angle width or back projection phase width set in advance as a default or by a user. It is possible. In the ex-post case, the correction angle width and the back projection phase width are changed and applied.
- a tomographic apparatus using X-rays has been described, but the present invention is not limited to this, and the present invention is also applicable to a tomographic apparatus using neutrons, positrons, gamma rays, and light. It is.
- the scanning method is not limited to any of the first, second, third, and fourth generations.
- the detectors are also arranged on a cylindrical surface centered on the X-ray source, a flat detector, a detector arranged on a spherical surface centered on the X-ray source, and on a cylindrical surface centered on the orbit axis. It can be applied to any type of detector, such as a detector.
- the generation of motion artifact was confirmed by fixing the corrected angle width index ⁇ to 0.4 and changing the back projection phase width index F to 0.8, 0.9, 1.0, and 1.1. Also in the conventional reconstruction method, the generation of motion artifact was confirmed using a weighting function in which the backprojection phase width index was similarly changed. In the conventional weight function, the slope of the weight is automatically determined when the back projection phase width index is determined.
- the image SD value was measured as an evaluation of noise.
- the noise amount slightly increases as compared with the conventional reconstruction method.
- the image SD and motion artifact are in a trade-off relationship, and the conventional reconstruction method could not suppress motion artifacts because the image SD was emphasized, whereas the reconstruction method of the present invention Now, it is possible to significantly reduce motion artifacts without causing a significant increase in noise.
- FIG. 1 is a diagram showing the overall configuration of an X-ray CT apparatus to which the present invention is applied
- FIG. 2 is a view for explaining an imaging method adopted by the X-ray CT apparatus of the present invention.
- FIG. 3 is a view showing a detector of the X-ray CT apparatus of the present invention, wherein (a) shows a single-row detector and (b) shows a multi-row detector.
- FIG. 4 is a diagram showing a procedure of image reconstruction according to the present invention, wherein (a) shows the first embodiment and (b) shows the second embodiment.
- FIG. 6 is a view for explaining the concept of a weight function used in the tomographic image reconstruction method of the present invention.
- FIG. 7 is a diagram showing an example of a parallel beam backprojection weighting function employed in the tomographic image reconstruction method of the present invention.
- FIG. 8 is a diagram showing an example of a parallel beam backprojection weighting function employed in the tomographic image reconstruction method of the present invention.
- FIG. 10 is a diagram showing an example of a fan beam back projection weight function employed in the tomographic image reconstruction method of the present invention.
- FIG. 11 is a diagram for explaining the shape of a weight function used in the present invention.
- FIG. 12 A diagram showing a relationship between a weight shape and a motion artifact in a reconstructed image.
- FIG. 13 A diagram showing a comparison between a reconstructed image by a reconstruction method of the present invention and a reconstructed image by a conventional reconstruction method.
- FIG. 15 A diagram for explaining a difference in redundancy depending on pixels when backprojection is performed from data having a minimum backprojection phase width of ⁇ + 2 ⁇ m.
- FIG. 16 is a diagram showing an example of a sinogram showing a minimum complete data set
- FIG. 17 is a diagram showing a conventional weight function in which a sinogram is divided into three or less regions including two triangular regions on a virtual sinogram, and different weights are given to the two triangular regions in a circular phase direction.
- Fig.18 Diagram for explaining the relationship between image noise and imaging time.
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)
- Image Analysis (AREA)
Abstract
Description
Claims
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005518018A JP4646810B2 (ja) | 2004-02-16 | 2005-02-15 | 断層撮影像の再構成方法及び断層撮影装置 |
US10/588,257 US7653224B2 (en) | 2004-02-16 | 2005-02-15 | Image reconstruction method and tomograph |
EP05710215.4A EP1716809B1 (en) | 2004-02-16 | 2005-02-15 | Tomogram reconstruction method and tomograph |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004-038833 | 2004-02-16 | ||
JP2004038833 | 2004-02-16 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2005077278A1 true WO2005077278A1 (ja) | 2005-08-25 |
Family
ID=34857821
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2005/002264 WO2005077278A1 (ja) | 2004-02-16 | 2005-02-15 | 断層撮影像の再構成方法及び断層撮影装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US7653224B2 (ja) |
EP (1) | EP1716809B1 (ja) |
JP (1) | JP4646810B2 (ja) |
CN (1) | CN100457043C (ja) |
WO (1) | WO2005077278A1 (ja) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007089673A (ja) * | 2005-09-27 | 2007-04-12 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
JP2007252898A (ja) * | 2006-02-27 | 2007-10-04 | Toshiba Corp | 画像表示装置及びx線ct装置 |
JP2007325853A (ja) * | 2006-06-09 | 2007-12-20 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
JP2012034972A (ja) * | 2010-08-10 | 2012-02-23 | Toshiba Corp | X線ct装置及び、x線ct装置の信号処理プログラム |
JP2012081279A (ja) * | 2010-10-13 | 2012-04-26 | Toshiba Corp | Ctイメージングの投影データ重み付け方法 |
JP2013166033A (ja) * | 2006-02-27 | 2013-08-29 | Toshiba Corp | 画像表示装置及びx線ct装置 |
WO2013161443A1 (ja) * | 2012-04-24 | 2013-10-31 | 株式会社 日立メディコ | X線ct装置及び画像再構成方法 |
WO2016129433A1 (ja) * | 2015-02-12 | 2016-08-18 | 株式会社日立製作所 | X線ct装置、画像処理装置及び画像再構成方法 |
JP5978429B2 (ja) * | 2011-04-28 | 2016-08-24 | 株式会社日立製作所 | 医用画像処理装置、医用画像処理方法 |
JP2020078462A (ja) * | 2018-11-13 | 2020-05-28 | キヤノンメディカルシステムズ株式会社 | X線診断システム及び再構成処理システム |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8938104B2 (en) * | 2008-08-29 | 2015-01-20 | Varian Medical Systems International Ag | Systems and methods for adaptive filtering |
JP2011005042A (ja) * | 2009-06-26 | 2011-01-13 | Canon Inc | 光音響イメージング装置及び光音響イメージング方法 |
JP5455694B2 (ja) * | 2010-02-09 | 2014-03-26 | 株式会社日立ハイテクノロジーズ | 荷電粒子線装置 |
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 |
RU2595766C2 (ru) * | 2010-10-27 | 2016-08-27 | Конинклейке Филипс Электроникс Н.В. | Идентификация и уменьшение искажения изображения |
JP5728304B2 (ja) * | 2011-06-21 | 2015-06-03 | 株式会社日立メディコ | X線ct装置及び画像再構成方法 |
JP5968316B2 (ja) * | 2011-07-08 | 2016-08-10 | 株式会社日立製作所 | 画像再構成装置及び画像再構成方法 |
CN104144650B (zh) * | 2012-03-27 | 2016-08-31 | 株式会社日立制作所 | 放射线拍摄装置和图像处理方法 |
JP5840070B2 (ja) * | 2012-05-08 | 2016-01-06 | 富士フイルム株式会社 | 光音響計測装置および光音響計測装置用プローブ |
DE102012217163B4 (de) * | 2012-09-24 | 2022-06-02 | Siemens Healthcare Gmbh | Verfahren zur Rekonstruktion von CT-Bilddaten mit gewichteter Rückprojektion, einschließlich Recheneinheit und CT-System für dieses Verfahren |
JP2014113444A (ja) * | 2012-11-16 | 2014-06-26 | Sony Corp | 画像処理装置、画像処理方法、およびプログラム |
WO2014199995A1 (ja) * | 2013-06-11 | 2014-12-18 | 株式会社 東芝 | X線コンピュータ断層撮影装置、およびスキャン開始タイミング決定方法 |
JP6640527B2 (ja) * | 2014-10-31 | 2020-02-05 | キヤノンメディカルシステムズ株式会社 | X線ct装置 |
JP6824133B2 (ja) * | 2017-09-28 | 2021-02-03 | 富士フイルム株式会社 | 画像処理装置、画像処理方法、及び画像処理プログラム |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS605127A (ja) * | 1983-05-02 | 1985-01-11 | ゼネラル・エレクトリツク・カンパニイ | 物体スライスの像を構成する装置 |
JPH08280664A (ja) * | 1995-04-13 | 1996-10-29 | Hitachi Medical Corp | 医療用診断ct装置 |
JP2001299738A (ja) * | 2000-04-14 | 2001-10-30 | Toshiba Corp | X線コンピュータ断層撮影装置 |
JP2002078702A (ja) * | 2000-09-06 | 2002-03-19 | Hitachi Medical Corp | 断層写真像作成方法及び断層写真像作成装置。 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NL8204086A (nl) * | 1982-10-22 | 1984-05-16 | Philips Nv | Werkwijze en inrichting voor het bepalen van een stralingsverzwakkingsverdeling in een vlak van een lichaam. |
US5606585A (en) * | 1995-12-21 | 1997-02-25 | General Electric Company | Methods and apparatus for multislice helical image reconstruction in a computer tomography system |
-
2005
- 2005-02-15 CN CNB2005800050669A patent/CN100457043C/zh not_active Expired - Fee Related
- 2005-02-15 US US10/588,257 patent/US7653224B2/en active Active
- 2005-02-15 EP EP05710215.4A patent/EP1716809B1/en not_active Not-in-force
- 2005-02-15 WO PCT/JP2005/002264 patent/WO2005077278A1/ja not_active Application Discontinuation
- 2005-02-15 JP JP2005518018A patent/JP4646810B2/ja not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS605127A (ja) * | 1983-05-02 | 1985-01-11 | ゼネラル・エレクトリツク・カンパニイ | 物体スライスの像を構成する装置 |
JPH08280664A (ja) * | 1995-04-13 | 1996-10-29 | Hitachi Medical Corp | 医療用診断ct装置 |
JP2001299738A (ja) * | 2000-04-14 | 2001-10-30 | Toshiba Corp | X線コンピュータ断層撮影装置 |
JP2002078702A (ja) * | 2000-09-06 | 2002-03-19 | Hitachi Medical Corp | 断層写真像作成方法及び断層写真像作成装置。 |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4509903B2 (ja) * | 2005-09-27 | 2010-07-21 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置 |
JP2007089673A (ja) * | 2005-09-27 | 2007-04-12 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
JP2016041387A (ja) * | 2006-02-27 | 2016-03-31 | 株式会社東芝 | 画像表示装置 |
JP2007252898A (ja) * | 2006-02-27 | 2007-10-04 | Toshiba Corp | 画像表示装置及びx線ct装置 |
JP2013166033A (ja) * | 2006-02-27 | 2013-08-29 | Toshiba Corp | 画像表示装置及びx線ct装置 |
JP2007325853A (ja) * | 2006-06-09 | 2007-12-20 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
JP4509971B2 (ja) * | 2006-06-09 | 2010-07-21 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置 |
JP2012034972A (ja) * | 2010-08-10 | 2012-02-23 | Toshiba Corp | X線ct装置及び、x線ct装置の信号処理プログラム |
JP2012081279A (ja) * | 2010-10-13 | 2012-04-26 | Toshiba Corp | Ctイメージングの投影データ重み付け方法 |
JP5978429B2 (ja) * | 2011-04-28 | 2016-08-24 | 株式会社日立製作所 | 医用画像処理装置、医用画像処理方法 |
WO2013161443A1 (ja) * | 2012-04-24 | 2013-10-31 | 株式会社 日立メディコ | X線ct装置及び画像再構成方法 |
US9406121B2 (en) | 2012-04-24 | 2016-08-02 | Hitachi Medical Corporation | X-ray CT apparatus and image reconstruction method |
JPWO2013161443A1 (ja) * | 2012-04-24 | 2015-12-24 | 株式会社日立メディコ | X線ct装置及び画像再構成方法 |
WO2016129433A1 (ja) * | 2015-02-12 | 2016-08-18 | 株式会社日立製作所 | X線ct装置、画像処理装置及び画像再構成方法 |
JPWO2016129433A1 (ja) * | 2015-02-12 | 2017-11-24 | 株式会社日立製作所 | X線ct装置、画像処理装置及び画像再構成方法 |
US10398392B2 (en) | 2015-02-12 | 2019-09-03 | Hitachi, Ltd. | X-ray CT apparatus, image processing device and image reconstruction method |
JP2020078462A (ja) * | 2018-11-13 | 2020-05-28 | キヤノンメディカルシステムズ株式会社 | X線診断システム及び再構成処理システム |
JP7305334B2 (ja) | 2018-11-13 | 2023-07-10 | キヤノンメディカルシステムズ株式会社 | X線診断システム及び再構成処理システム |
Also Published As
Publication number | Publication date |
---|---|
US7653224B2 (en) | 2010-01-26 |
CN100457043C (zh) | 2009-02-04 |
JPWO2005077278A1 (ja) | 2007-10-18 |
EP1716809B1 (en) | 2013-11-06 |
CN1917811A (zh) | 2007-02-21 |
US20080273778A1 (en) | 2008-11-06 |
EP1716809A1 (en) | 2006-11-02 |
EP1716809A4 (en) | 2011-04-06 |
JP4646810B2 (ja) | 2011-03-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2005077278A1 (ja) | 断層撮影像の再構成方法及び断層撮影装置 | |
JP5106978B2 (ja) | X線ct装置 | |
US8031831B2 (en) | Voltage and or current modulation in dual energy computed tomography | |
JP6169558B2 (ja) | コントラスト依存の解像度をもつ画像 | |
JP5194095B2 (ja) | 半撮影域のみをカバーする縮小サイズ検出器を利用するコンピュータ断層撮影システムに用いる装置及び方法 | |
JP5220309B2 (ja) | X線断層撮影装置 | |
US20050249431A1 (en) | Method for post- reconstructive correction of images of a computer tomograph | |
US10789738B2 (en) | Method and apparatus to reduce artifacts in a computed-tomography (CT) image by iterative reconstruction (IR) using a cost function with a de-emphasis operator | |
JP6482934B2 (ja) | 画像処理装置、放射線検出装置および画像処理方法 | |
JP5726288B2 (ja) | X線ct装置、および方法 | |
JP2008279153A (ja) | X線ct装置 | |
JP2008006032A (ja) | X線ct装置およびx線ct撮影方法 | |
US7792238B2 (en) | Method and system for reconstructing cone-beam projection data with reduced artifacts | |
WO2016017402A1 (ja) | データ処理方法、データ処理装置、及びx線ct装置 | |
JP5606667B2 (ja) | X線ct装置 | |
US9330443B1 (en) | Noise reduction in image domain for spectral computed tomography | |
EP3221848B1 (en) | X-ray imaging device for an object | |
US10070841B2 (en) | Arithmetic device, X-ray CT apparatus, and image reconstruction method | |
US10383589B2 (en) | Direct monochromatic image generation for spectral computed tomography | |
US10950012B2 (en) | Radiation image generation method and image processing device | |
CN110073412B (zh) | 使用交替否定的图像噪声估计 | |
JP4887132B2 (ja) | X線ct装置 | |
JP7483654B2 (ja) | 医用画像処理装置および医用画像処理方法 | |
WO2017130657A1 (ja) | X線ct装置、撮影条件設定方法及び撮影条件設定プログラム | |
CN107341836B (zh) | 一种ct螺旋扫描图像重建方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AK | Designated states |
Kind code of ref document: A1 Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BW BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE EG ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NA NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW |
|
AL | Designated countries for regional patents |
Kind code of ref document: A1 Designated state(s): GM KE LS MW MZ NA SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IS IT LT LU MC NL PL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application | ||
WWE | Wipo information: entry into national phase |
Ref document number: 2005518018 Country of ref document: JP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2005710215 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 10588257 Country of ref document: US |
|
WWE | Wipo information: entry into national phase |
Ref document number: 200580005066.9 Country of ref document: CN |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
WWW | Wipo information: withdrawn in national office |
Country of ref document: DE |
|
WWP | Wipo information: published in national office |
Ref document number: 2005710215 Country of ref document: EP |