WO2018167355A1 - Patient movement correction method for cone-beam computed tomography - Google Patents

Patient movement correction method for cone-beam computed tomography Download PDF

Info

Publication number
WO2018167355A1
WO2018167355A1 PCT/FI2018/000005 FI2018000005W WO2018167355A1 WO 2018167355 A1 WO2018167355 A1 WO 2018167355A1 FI 2018000005 W FI2018000005 W FI 2018000005W WO 2018167355 A1 WO2018167355 A1 WO 2018167355A1
Authority
WO
WIPO (PCT)
Prior art keywords
ray
projection
coordinate system
image
estimate
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
PCT/FI2018/000005
Other languages
English (en)
French (fr)
Inventor
Mikko LILJA
Kalle KARHU
Jaakko Lähelmä
Kustaa Nyholm
Ari Hietanen
Timo Müller
Sakari Kettunen
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Planmeca Oy
Original Assignee
Planmeca Oy
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Planmeca Oy filed Critical Planmeca Oy
Priority to EP18723569.2A priority Critical patent/EP3595528B1/en
Priority to JP2019551297A priority patent/JP7116078B2/ja
Priority to RU2019131095A priority patent/RU2766743C1/ru
Priority to CN201880028718.8A priority patent/CN110602990B/zh
Priority to BR112019019206-9A priority patent/BR112019019206B1/pt
Priority to KR1020197030330A priority patent/KR102369653B1/ko
Publication of WO2018167355A1 publication Critical patent/WO2018167355A1/en
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/008Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/40Arrangements for generating radiation specially adapted for radiation diagnosis
    • A61B6/4064Arrangements for generating radiation specially adapted for radiation diagnosis specially adapted for producing a particular type of beam
    • A61B6/4085Cone-beams
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/46Arrangements for interfacing with the operator or the patient
    • A61B6/461Displaying means of special interest
    • A61B6/466Displaying means of special interest adapted to display 3D data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • A61B6/5264Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to motion
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/005Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/006Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/44Constructional features of apparatus for radiation diagnosis
    • A61B6/4429Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
    • A61B6/4435Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure
    • A61B6/4441Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure the rigid structure being a C-arm or U-arm
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/501Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the head, e.g. neuroimaging or craniography
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10116X-ray image
    • G06T2207/10124Digitally reconstructed radiograph [DRR]
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00Image generation
    • G06T2211/40Computed tomography
    • G06T2211/424Iterative

Definitions

  • the invention concerns patient movement correction methods in the field of cone-beam computed tomography.
  • the following disclosure relates to three-dimensional X-ray cone-beam computed tomography for medical applications, wherein a plurality of X-ray projection images acquired from different angles are used to reconstruct 3-D cross-sectional images of an anatomy of a patient.
  • the duration of the X-ray projection image acquisition is typically of the order of 10—30 seconds because the X-ray tube (X-ray source) and sensor (X-ray detector) must physically travel the spatial trajectory corresponding to the acquisition angles.
  • the imaging trajectory is typically realized by a rotation and translation mechanism.
  • the imaging trajectory should be known with a sufficient accuracy and the imaged object should remain sufficiently stationary during the X- ray projection image acquisition in order for the reconstructed CBCT image to be sharp and true to the anatomy, because the projection image measurements are assumed to represent co-registered integrated views of a stationary object. This results in a set of geometrically consistent measurements that can be used to reconstruct the attenuation distribution representing the studied anatomy. Whenever these assumptions are violated, the accuracy of the reconstructed image is degraded as a consequence of the projection measurements becoming mutually inconsistent.
  • the accuracy of the final image typically depends on how accurately the assumptions made in the reconstruction process correspond to the actual, physical image acquisition process.
  • the estimated spatial positions of the X-ray source and detector corresponding to each acquired X-ray image affect the computation of the ray paths during the CBCT reconstruction process. Due to inherent manufacturing and operating tolerances and potential deformation of the imaging device the realized rotation angles and positions tend to deviate from the ideal values according to the assumed form of the imaging trajectory. A systematic deviation, however, can be addressed by using different calibration methods that are repeated after certain period of time or operation cycles.
  • CBCT imaging applying a fixed coordinate system for the modelling and compensation of patient motion is not ideal, as the intrinsic geometric degrees-of-freedom in CBCT imaging are not separated by the coordinate system.
  • a CBCT imaging device it is particular that the X-ray beams diverge and form a pyramid-shaped cone.
  • a shift along the isoray adjoining the X-ray source and the center of the X-ray detector will only affect the magnification factor, whereas an in-plane shift along the X-ray detector' s pixel array will result in a maximal shift of the imaged object within its projection image.
  • preventing a net transformation arising as a result of an applied geometric correction by known means of a rigid registration of the resulting corrected CBCT reconstruction and the uncorrected CBCT reconstruction is computationally expensive, especially if applied repeatedly during the geometric correction process.
  • CBCT medical cone-beam computed tomography
  • the purpose of the correction process is to improve the resulting CBCT reconstruction image quality by improving the geometric consistency of the X-ray image measurements and, in turn, better to better satisfy the requirements of tomographic image reconstruction.
  • the intrinsic coordinate system enables defining the applied geometric degrees-of-freedom in a manner that corresponds to their relative importance to both the correction process as well as the resulting image quality.
  • the disclosed patient movement correction process takes as its input the data that is normally required for computing a CBCT image reconstruction: a set of X-ray projection images and an estimate of the 3-D projection geometry corresponding to the spatial positions of the X-ray source and X-ray detector during the acquisition of the X- ray projection images.
  • an intermediate reconstruction is first computed using the estimated projection geometry.
  • a corrective geometric transformation improving the geometric correspondence of each accessed X-ray projection image with the rest of the X-ray projection images is established with the transformation corresponding to a virtual movement of the X-ray source and detector during the image acquisition.
  • a projection image-specific rotating coordinate system is applied for determining the corrective geometric transformation.
  • the correction process including the computation of an intermediate reconstruction using the current estimate of the projection geometry and the subsequent optimization of the corrective transformations, may be iterated for a number of times.
  • a final CBCT reconstruction is computed using the X-ray projection images and the final estimate of the projection geometry corresponding to the corrective transformations.
  • the benefit of the disclosed patient movement correction process before reconstructing the final CBCT image is an improvement in the projection geometric consistency, which results in a higher image quality in terms of sharpness, level of detail and contrast.
  • a further benefit of the correction process is that by potentially preventing a re-scan due to projection geometry-related image quality degradation, the radiation dose incurred to the patient may be decreased.
  • Fig. 1 shows an example of one typical CBCT imaging apparatus.
  • Fig. 2 shows a flow chart of the disclosed patient movement correction method.
  • Fig. 3 shows a flow chart detailing the step 203 of Fig. 2.
  • Fig. 4 shows in the context of a CBCT imaging apparatus a rotating coordinate system applied in the disclosed patient movement correction method.
  • a number of 2-D X-ray projection images are measured using an X-ray source and X-ray detector that are rotated around the imaged anatomy.
  • the exposure of each X-ray projection image takes typically in the range of milliseconds and when using typical scanning velocities, continuous movement of the X-ray source and detector does not any significant motion blurring in the projection images.
  • the 3-D spatial positions of the X-ray source and detector in a suitable reference coordinate system comprising the imaging or projection geometry must be available to the reconstruction algorithm.
  • a typical description of the imaging geometry consists of the 3-D positions of the focus of the X-ray source and the center of the X- ray detector, and information to sufficiently uniquely determine the orientation of the X-ray detector. Such information may consist of e.g. rotation angles in a reference coordinate system that, when applied according to a predefined convention, will determine the directions of the horizontal and vertical axes of the detector's pixel array.
  • the projection geometry description is typically based on knowledge of the imaging device's physical measures as well as the ideal exposure trajectories corresponding to given imaging program. Furthermore, a periodic calibration process is typically performed to ensure a sufficient accuracy of the projection geometry.
  • the spatial propagation of X-radiation through the imaged anatomy is modeled.
  • a rectilinear propagation is assumed for simplicity and the X-ray beam paths from the source to the detector are modeled as line integrals, which are computed based on the information contained in the projection geometry description.
  • the spatial overlaps of the X-ray paths and the elements of the 3-D image voxel array used for the reconstruction are solved using a suitable projector algorithm and the projection geometry description.
  • the total X-ray attenuation distribution corresponding to the imaged anatomy can be reconstructed using a suitable method such as the well-known FDK algorithm.
  • a fundamental assumption employed in the reconstruction process is that the imaged anatomy has remained sufficiently stationary during the acquisition of the X-ray projection images.
  • the rationale is that the X-ray projections should represent co- registered measurements of a stationary object that can then be consistently combined to reconstruct the 3-D structure of the object.
  • a systematic error in the projection geometry can be compensated by a periodic calibration process, which eliminates the effects arising from any deviations from the assumed, ideal X-ray projection image acquisition trajectory.
  • a more difficult problem arises from patient movement during the imaging, which is unpredictable and random in its nature. Although it prevented to a degree by supporting the patient during the imaging, the relatively long duration of X-ray projection image acquisition in CBCT imaging, of the order of 10 seconds, makes it unfeasible to completely eliminate patient movement.
  • an intermediate CBCT reconstruction is first computed using the 2-D X-ray projection images and the estimated imaging geometry. It is sufficient to compute the intermediate CBCT reconstruction at a coarser resolution than is typically used when making reconstructions for diagnostic purposes.
  • the intermediate reconstruction serves to aggregate the information from all X-ray projection images with the appearance of the reconstruction reflecting the mutual geometric consistency of the measured 2-D X-ray projection images.
  • the optimization of the projection-image-specific geometry is based on measuring the similarity of the physical X-ray projection images and the corresponding re-projected data of the intermediate CBCT reconstruction that forms a digitally reconstructed radiograph (DRR) .
  • DRR digitally reconstructed radiograph
  • the intermediate reconstruction serves as an aggregate of all projection images
  • its reprojection reflects the sum of all projection images and the maximum similarity of the measured X-ray projection image and the forward projected images can be expected to be maximized when the corresponding reprojection geometry matches the average correct projection geometry in the sense of the intermediate reconstruction.
  • the emitted X-ray beams diverge and form a pyramid-shaped cone.
  • a shift along the isoray adjoining the X- ray source and the center of the X-ray detector will only affect the magnification factor, whereas a shift along the X-ray detector plane will result in a maximal shift of the imaged object within its projection image.
  • the disclosed approach adopts a rotating coordinate system that is attached to the physical positions of the X-ray source and detector during the image acquisition. Namely, two of the coordinate axes are attached to the rectangular X-ray detector pixel array and the remaining perpendicular axis to the normal of the detector' s pixel array.
  • the net transformation can be estimated by mapping the transformation corresponding to each transformed projection image from the rotating coordinate system to the fixed coordinate system.
  • the corresponding displacements in the fixed coordinate system can be computed based on the known horizontal axes of the projection-specific rotated coordinate systems, and the average value can be taken to represent the net displacement in the fixed coordinate system.
  • the inverse of the net transformation can be mapped back to the rotating coordinate systems and subtracted from the projection image- specific transformations. As a result, the net transformation in the fixed coordinate system is eliminated.
  • the goodness of a corrective geometric transformation of a given X-ray projection image is measured by the similarity of the forward-projected image corresponding to the transformed projection geometry and the original X-ray projection image.
  • the similarity of the reprojection and the X-ray projection can be measured, e.g. by the mean squared difference of the images, correlation coefficient, or gradient correlation coefficient.
  • the optimal (in the sense of the similarity measure) geometric transformation for each projection image given an intermediate CBCT reconstruction can then be determined by finding the extremum of the similarity measure between the forward-projected image and the X-ray projection image as a function of the parameters of the geometric transformation.
  • the optimization process then includes computing an intermediate CBCT reconstruction using the initial estimate of the X-ray projection geometry; accessing all or a subset of the measured X-ray projection images; (for each accessed projection image) establishing a projection image-specific corrective transformation by finding the maximum similarity between the measured X-ray projection image and the corresponding forward projection of the intermediate CBCT reconstruction as a function of the parameters of the geometric transformation performed in the rotating coordinate system; estimating the net transformation in a fixed reference coordinate system and subtracting the corresponding transformation from the transformation parameters in the rotating coordinate system; computing a final CBCT reconstruction when a sufficient correction result is estimated to have been obtained.
  • a medical CBCT imaging apparatus 100 which includes a vertical base construction 101 from which horizontally extended a support structure 102, a patient support means 107 and an arm part 103 which supports a structure supporting the imaging means, an arm part 104.
  • an X-ray source 105 and a receiver means of X-ray image information (X- ray detector) 109 that are arranged with respect to the patient support means 107 such that an imaging station 108 positioned between the X-ray source 105 and the receiver means of X-ray image information 109 is formed such that a beam generated by the X-ray source 105 is alignable to go through the imaging station 108 towards the receiver means of X-ray image information 109.
  • the arm part 104 supporting the imaging means is arranged to be rotatable, and also its location with respect to the structure supporting it 103 and/or the patient support station 108 may be arranged changeable.
  • the arrangement includes a control means, of which Fig. 1 shows a control panel 106 placed in connection with the support structure 102 supporting the patient support means 107.
  • the imaging apparatus 100 can be arranged to be connected to a controller 110 via a cable, the controller including a computer arranged with a means for processing image information produced by the imaging apparatus, and a display 111 on which images can be shown.
  • the controller 110 further comprises at least one processor 112 and at least one memory 113.
  • the at least one processor 112 may be configured to execute computer programs and the at least one memory 113 is configured to store computer programs and related data.
  • the controller 110 may be a general purpose computer or a specifically manufactured device for implementing the process described below.
  • Figure 2 describes the steps of the disclosed patient movement correction method that may be used, for example, for processing images acquire by the imaging arrangement of figure 1.
  • the method is based on finding corrective geometric transformations for the initially estimated projection geometry in a rotating coordinate system attached to the positions of the X-ray source 105 and X-ray detector 109.
  • the input data is acquired comprising the measured X-ray projection images and an initial estimate of the projection geometry corresponding to the acquisition process.
  • the projection geometry defines the physical trajectory of the X-ray source 105 and X-ray detector 109 when they are rotated and translated around the imaging station 108 by means of the arms 103 and 104.
  • the projection geometry also determines the estimated orientation of the X-ray detector 109.
  • the form of the projection geometry is typically based on a geometry calibration procedure utilizing e.g. a known reference phantom with radiopaque markers.
  • an intermediate CBCT reconstruction for the purpose of the projection optimization method is computed using the input data acquired in step 200.
  • the intermediate reconstruction is understood to aggregate all available physical and geometric information acquired during the measurement process in step 200. In the event of geometric inconsistency, this is reflected by the intermediate reconstruction by e.g. blurriness of the reconstructed details.
  • step 202 corrective geometric transformations for the projection images are established in the rotating coordinate system.
  • the purpose of the geometric transformations is to compensate for the intrinsic geometric inconsistency in the initial estimate of the projection geometry.
  • the details of step 202 are explained below in reference to Fig. 3 but, more generally, the corrective geometric transformation is sought by finding the optimal geometric transformation.
  • the goodness of a transformation is defined by assigning a similarity value to it.
  • the similarity value is computed by comparing a reprojected digitally reconstructed radiograph (DRR) of the intermediate CBCT image reconstruction to the corresponding measured X-ray projection image, with the applied projection geometry corresponding to the evaluated geometric transformation.
  • a higher similarity value is taken as an indication of a better corrective geometric transformation.
  • step 203 the net geometric transformation in a fixed coordinate system is subtracted.
  • the fixed coordinate system is defined typically in reference to static components of the imaging apparatus 100, such as components of the imaging arrangement fixed in connection with the imaging station 108.
  • the net geometric transformation is computed by linearly transforming the geometric transformations established in step 202 from the rotating coordinate system to the fixed coordinate system.
  • the linear transformation is readily obtained based on the known coordinate axes of the rotated and fixed coordinate systems.
  • its inverse is linearly transformed from the fixed coordinate system to its rotated coordinate system.
  • the inverse of the net transformation then corresponds to a set of geometric transformation parameters for each projection image in the rotating coordinate system. Adding these values to the geometric parameter values established in step 202 will result in the net transformation being canceled out in the fixed coordinate system.
  • step 204 the transformation parameter values resulting from step 203 are applied to the initial projection geometry estimate to obtain the corrected projection geometry estimate.
  • steps 201- 203 may be repeated iteratively a plurality of times before advancing to step 204.
  • step 205 the final CBCT reconstruction is computed using the corrected projection geometry estimate.
  • the final CBCT reconstruction is computed in a normal manner with the exception to the situation where steps 201-204 were not applied that the initial projection geometry estimate is replaced by the corrected projection geometry estimate obtained in steps 201-20 .
  • Figure 3 describes as an example details of step 202 according to Fig 2.
  • step 300 the assessed geometric transformation is applied to the initial geometry of a projection image in the rotating coordinate system. Again, the coordinate system is defined to coincide with the spatial positions and orientation of the X-ray source 105 and X-ray detector 109 during the physical acquisition of the considered projection image.
  • Applying the transformation in the rotating frame of reference involves a linear mapping from the fixed coordinate system, where the imaging geometry is typically defined, to the rotating coordinate system, where the transformation is performed, followed by an inverse linear mapping from the rotating coordinate system to the fixed coordinate system.
  • the initial projection geometry of the projection image is mapped to the rotating coordinate system, translated along this axis of the rotating coordinate system by the given amount, which corresponds to a virtual movement of the X-ray source and detector, and then mapped back to the fixed coordinate system.
  • the transformed projection geometry is expressed in the fixed coordinate system and the obtained transformed projection geometry is used as input for the subsequent step 301.
  • a reprojected DRR image of the intermediate CBCT reconstruction computed in step 201 is computed using the transformed projection geometry obtained in step 300.
  • the computation of the DRR image can be performed using a standard algorithm such as the Siddon raycasting method.
  • the input of the algorithm consists of the end points of the 3-D X-ray representation and the source image of which the DRR image is computed including knowledge of its spatial position and orientation expressed in the same coordinate system. In the described setting, this coordinate system corresponds to the fixed coordinate system.
  • the effect of the geometric transformation applied in step 300 is to change the end points of each virtual X-ray path through the intermediate CBCT reconstruction, which propagates the effect of the geometric transformation to the obtained DRR image.
  • step 302 the similarity between the DRR image obtained in step 301 and the X-ray projection image acquired in step 200 is evaluated.
  • the similarity is based on a pointwise comparison of the images using established approaches such as the average squared difference of the images or their cross-correlation.
  • the specific measure used for evaluating the similarity is not significant to the described method.
  • the obtained similarity value is assigned to the geometric parameters given as input to step 300. A higher similarity is taken as an indication of more suitable geometric transformation parameters.
  • step 303 the optimal geometric transformation parameters are established by finding the parameters corresponding to the highest similarity value obtained by applying steps 300-302.
  • steps 300-302 are repeatedly evaluated by a suitable minimization algorithm such as the well-known Nelder- ead simplex algorithm to establish the optimal geometric transformation parameters.
  • a suitable minimization algorithm such as the well-known Nelder- ead simplex algorithm to establish the optimal geometric transformation parameters.
  • Figure 4 illustrates the rotating coordinate system applied in step 202.
  • a C-arm part 403 typically supports the X-ray source 404 and X-ray detector 405.
  • the C-arm is supported by the vertical base construction 401 and a shoulder arm part 402.
  • the X-ray beams diverge and form a pyramid-shaped cone.
  • a shift along the isoray adjoining the X-ray source 404 and the center of the X-ray detector 405 will only affect the magnification factor, whereas a shift in the plane of the X-ray detector's pixel array will result in a maximal shift of the imaged object within its projection image.
  • the disclosed approach adopts a rotating uvw coordinate system 406 that is attached to the physical positions and orientations of the X-ray source and detector pixel array during the image acquisition.
  • a rotating uvw coordinate system 406 that is attached to the physical positions and orientations of the X-ray source and detector pixel array during the image acquisition.
  • the result of the disclosed patient movement correction method is an improved estimate of the projection geometry corresponding to the physical acquisition of the X-ray projection images and correspondingly a CBCT reconstruction image, where the effect of geometric inconsistency has been reduced.
  • the disclosed method may be implemented as computer software executed in a computing device.
  • the software is embodied on a computer readable medium so that it can be provided to the computing device, such as the controller 110 of figure 1.
  • the components of the exemplary embodiments can include computer readable medium or memories for holding instructions programmed according to the teachings of the present embodiments and for holding data structures, tables, records, and/or other data described herein.
  • Computer readable medium can include any suitable medium that participates in providing instructions to a processor for execution.
  • Computer-readable media can include, for example, a floppy disk, a flexible disk, hard disk, magnetic tape, any other suitable magnetic medium, a CD-ROM, CD ⁇ R, CD+RW, DVD, DVD- RAM, DVD ⁇ RW, DVD+R, HD DVD, HD DVD-R, HD DVD-R , HD DVD-RAM, Blu-ray Disc, any other suitable optical medium, a RAM, a PROM, an EPROM, a FLASH-EPROM, any other suitable memory chip or cartridge, a carrier wave or any other suitable medium from which a computer can read.

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Pathology (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Pulmonology (AREA)
  • Human Computer Interaction (AREA)
  • Quality & Reliability (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Neurology (AREA)
  • Neurosurgery (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
PCT/FI2018/000005 2017-03-17 2018-03-19 Patient movement correction method for cone-beam computed tomography Ceased WO2018167355A1 (en)

Priority Applications (6)

Application Number Priority Date Filing Date Title
EP18723569.2A EP3595528B1 (en) 2017-03-17 2018-03-19 Patient movement correction for cone-beam computed tomography
JP2019551297A JP7116078B2 (ja) 2017-03-17 2018-03-19 コーンビーム型コンピュータ断層撮影における患者の動き補正方法
RU2019131095A RU2766743C1 (ru) 2017-03-17 2018-03-19 Способ коррекции перемещения пациента для конусно-лучевой компьютерной томографии
CN201880028718.8A CN110602990B (zh) 2017-03-17 2018-03-19 用于锥形束计算机断层扫描的患者移动校正方法
BR112019019206-9A BR112019019206B1 (pt) 2017-03-17 2018-03-19 Método de correção de movimento de paciente para tomografia computadorizada de feixe cônico, meio de armazenamento legível por computador, e, aparelho
KR1020197030330A KR102369653B1 (ko) 2017-03-17 2018-03-19 콘-빔 컴퓨터 단층촬영을 위한 환자 움직임 정정 방법

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FI20175244A FI20175244A7 (fi) 2017-03-17 2017-03-17 Itsekalibroiva lääketieteellinen kuvannuslaite
FI20175244 2017-03-17

Publications (1)

Publication Number Publication Date
WO2018167355A1 true WO2018167355A1 (en) 2018-09-20

Family

ID=62143215

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/FI2018/000005 Ceased WO2018167355A1 (en) 2017-03-17 2018-03-19 Patient movement correction method for cone-beam computed tomography

Country Status (8)

Country Link
US (1) US10636179B2 (enExample)
EP (1) EP3595528B1 (enExample)
JP (1) JP7116078B2 (enExample)
KR (1) KR102369653B1 (enExample)
CN (1) CN110602990B (enExample)
FI (1) FI20175244A7 (enExample)
RU (1) RU2766743C1 (enExample)
WO (1) WO2018167355A1 (enExample)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3639752B1 (en) 2018-10-03 2025-05-21 Canon Medical Systems Corporation Analyzing apparatus and analyzing program
CN110766629B (zh) * 2019-10-17 2022-03-01 广州华端科技有限公司 Cbct系统几何校正方法、装置、计算机设备和存储介质
CN111553849B (zh) * 2020-03-28 2023-03-24 中国人民解放军战略支援部队信息工程大学 基于局部特征匹配的锥束ct几何伪影去除方法及装置
CN113520421A (zh) * 2020-05-06 2021-10-22 北京友通上昊科技有限公司 X射线切面成像方法、存储介质及成像系统
CN112529820B (zh) * 2020-12-29 2024-06-14 广州开云影像科技有限公司 Cbct系统半扇束几何校正方法、装置、设备和介质
CN112837246A (zh) * 2021-03-25 2021-05-25 中科超精(南京)科技有限公司 一种图像引导系统几何校正方法
US11712215B2 (en) 2021-04-13 2023-08-01 Canon Medical Systems Corporation Devices, systems, and methods for motion-corrected medical imaging
JP7563304B2 (ja) * 2021-06-16 2024-10-08 トヨタ自動車株式会社 X線ct装置及び画像生成方法
EP4120182B1 (de) * 2021-07-12 2025-08-27 DENTSPLY SIRONA Inc. Selbstkalibrierendes dentales dvt unterstützt durch maschinelles lernen
CN113963056B (zh) * 2021-09-07 2022-08-26 于留青 Ct图像重建方法、装置、电子设备以及存储介质
DE102022200806A1 (de) * 2022-01-25 2023-07-27 Siemens Healthcare Gmbh Bereitstellen eines Ergebnisdatensatzes
CN114612569A (zh) * 2022-02-23 2022-06-10 广州华端科技有限公司 一种几何标定中标记物重叠的处理方法与装置
CN116894798B (zh) * 2023-09-11 2023-12-05 金华飞光科技有限公司 光固化3d打印机的投影畸形矫正方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060058620A1 (en) * 2004-08-19 2006-03-16 Georgios Sakas Setting a relative position of a radiation device and a patient
US20110116703A1 (en) * 2006-11-02 2011-05-19 Accuray Incorporated Target tracking using direct target registration
US20120301004A1 (en) * 2010-01-13 2012-11-29 The Australian National University Computed Tomography Imaging Process And System
US20150103972A1 (en) * 2012-06-05 2015-04-16 Koninklijke Philips N.V. Motion layer decomposition calibration of x-ray ct imagers

Family Cites Families (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5251128A (en) * 1990-11-19 1993-10-05 General Electric Company Motion artifact reduction in projection imaging
EP1374949B1 (en) * 2002-06-17 2005-12-14 Nucletron B.V. Real time radiation treatment planning system
US6888924B2 (en) * 2003-06-25 2005-05-03 General Electric Company Method, apparatus, and medium for calibration of tomosynthesis system geometry using fiducial markers with non-determined position
US7653229B2 (en) * 2003-12-23 2010-01-26 General Electric Company Methods and apparatus for reconstruction of volume data from projection data
JP4438053B2 (ja) * 2004-05-11 2010-03-24 キヤノン株式会社 放射線撮像装置、画像処理方法及びコンピュータプログラム
WO2006082557A2 (en) * 2005-02-01 2006-08-10 Koninklijke Philips Electronics N.V. Apparatus and method for correction or extension of x-ray projections
DE102005005087B4 (de) * 2005-02-03 2012-09-27 Siemens Ag Verfahren zur Korrektur nicht reproduzierbarer Geometriefehler eines Röntgen-C-Bogengerätes
US8199168B2 (en) * 2005-11-15 2012-06-12 General Electric Company System and method for 3D graphical prescription of a medical imaging volume
EP2086502B1 (en) 2006-11-09 2013-04-17 Alcon Research, Ltd. Punctal plug comprising a water-insoluble polymeric matrix
FR2919096A1 (fr) * 2007-07-19 2009-01-23 Gen Electric Procede de correction de recalage d'images radiographiques
WO2010103623A1 (ja) 2009-03-10 2010-09-16 三菱重工業株式会社 放射線治療装置制御装置および特定部位位置計測方法
JP6214128B2 (ja) * 2010-11-22 2017-10-18 キヤノン株式会社 画像処理装置、画像処理方法、及び記憶媒体
CN102488528B (zh) * 2011-12-07 2013-04-24 华中科技大学 一种层析成像几何参数的校准方法
KR20150079560A (ko) * 2012-08-20 2015-07-08 오란게덴탈 게엠베하 운트 코카게 콘빔 컴퓨터 단층 촬영 장치의 기하 특성화 및 교정
CN103729827A (zh) * 2012-10-11 2014-04-16 南京普爱射线影像设备有限公司 影像增强器c形臂x射线机三维重建重叠伪影校正方法
CN103519833B (zh) * 2013-06-05 2015-10-07 东南大学 一种旋转c型臂x射线机的三维校正重建方法
US10512441B2 (en) * 2013-08-29 2019-12-24 Siemens Healthcare Gmbh Computed tomography having motion compensation
EP3175791B1 (en) * 2013-11-04 2021-09-08 Ecential Robotics Method for reconstructing a 3d image from 2d x-ray images
CN103735282B (zh) * 2014-01-06 2015-09-23 北京大学 一种锥束ct系统探测器几何校正装置及其校正方法
US10478148B2 (en) * 2016-02-19 2019-11-19 The Johns Hopkins University Self-calibrating projection geometry for volumetric image reconstruction

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060058620A1 (en) * 2004-08-19 2006-03-16 Georgios Sakas Setting a relative position of a radiation device and a patient
US20110116703A1 (en) * 2006-11-02 2011-05-19 Accuray Incorporated Target tracking using direct target registration
US20120301004A1 (en) * 2010-01-13 2012-11-29 The Australian National University Computed Tomography Imaging Process And System
US20150103972A1 (en) * 2012-06-05 2015-04-16 Koninklijke Philips N.V. Motion layer decomposition calibration of x-ray ct imagers

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
KAROLCZAK MAREK ET AL: "Implementation of a cone-beam reconstruction algorithm for the single-circle source orbit with embedded misalignment correction using homogeneous coordinates", MEDICAL PHYSICS, AIP, MELVILLE, NY, US, vol. 28, no. 10, 1 October 2001 (2001-10-01), pages 2050 - 2069, XP012011303, ISSN: 0094-2405, DOI: 10.1118/1.1406514 *

Also Published As

Publication number Publication date
KR102369653B1 (ko) 2022-03-07
EP3595528B1 (en) 2023-07-19
FI20175244L (fi) 2018-09-18
CN110602990B (zh) 2023-09-01
JP2020509887A (ja) 2020-04-02
EP3595528A1 (en) 2020-01-22
FI20175244A7 (fi) 2018-09-18
US20180268574A1 (en) 2018-09-20
CN110602990A (zh) 2019-12-20
BR112019019206A2 (pt) 2020-04-14
RU2766743C1 (ru) 2022-03-15
JP7116078B2 (ja) 2022-08-09
EP3595528C0 (en) 2023-07-19
KR20190126882A (ko) 2019-11-12
US10636179B2 (en) 2020-04-28

Similar Documents

Publication Publication Date Title
EP3595528B1 (en) Patient movement correction for cone-beam computed tomography
US12056794B2 (en) Self-calibration procedure for digital breast tomosynthesis imaging apparatus
US6888924B2 (en) Method, apparatus, and medium for calibration of tomosynthesis system geometry using fiducial markers with non-determined position
US9082182B2 (en) Extracting patient motion vectors from marker positions in x-ray images
JP4904048B2 (ja) イメージング・システムの取得時幾何学的構成のための方法及び装置
US9538974B2 (en) Methods and systems for correcting table deflection
US20100266220A1 (en) Features-based 2d-3d image registration
US20110123081A1 (en) Correcting and reconstructing x-ray images using patient motion vectors extracted from marker positions in x-ray images
Niebler et al. Projection‐based improvement of 3D reconstructions from motion‐impaired dental cone beam CT data
US20170086767A1 (en) Image quality in computed tomography using redundant information in production data sets
WO2017071956A1 (en) Computed tomography image generation apparatus
RU2789105C2 (ru) Способ автокалибровки устройства для формирования изображений цифровой томографической реконструкции груди
HK40018965A (en) Patient movement correction method for cone-beam computed tomography
HK40018965B (zh) 用於锥形束计算机断层扫描的患者移动校正方法
CN114119801A (zh) 三维数字减影血管造影方法、装置、电子设备及存储介质
CN113963056A (zh) Ct图像重建方法、装置、电子设备以及存储介质
BR112019019206B1 (pt) Método de correção de movimento de paciente para tomografia computadorizada de feixe cônico, meio de armazenamento legível por computador, e, aparelho
US20230237716A1 (en) Providing a result dataset
Wicklein et al. An object-independent measure for improving misalignment correction in c-arm ct
RU2021108170A (ru) Способ автокалибровки устройства для формирования изображений цифровой томографической реконструкции груди

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: 18723569

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2019551297

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

REG Reference to national code

Ref country code: BR

Ref legal event code: B01A

Ref document number: 112019019206

Country of ref document: BR

ENP Entry into the national phase

Ref document number: 20197030330

Country of ref document: KR

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 2018723569

Country of ref document: EP

ENP Entry into the national phase

Ref document number: 2018723569

Country of ref document: EP

Effective date: 20191017

ENP Entry into the national phase

Ref document number: 112019019206

Country of ref document: BR

Kind code of ref document: A2

Effective date: 20190916