WO2014153561A1 - Computer accessible medium for modulating x-ray beam intensity - Google Patents
Computer accessible medium for modulating x-ray beam intensity Download PDFInfo
- Publication number
- WO2014153561A1 WO2014153561A1 PCT/US2014/031609 US2014031609W WO2014153561A1 WO 2014153561 A1 WO2014153561 A1 WO 2014153561A1 US 2014031609 W US2014031609 W US 2014031609W WO 2014153561 A1 WO2014153561 A1 WO 2014153561A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- arrangement
- radiation
- ray
- resultant
- ray radiation
- Prior art date
Links
- 230000005855 radiation Effects 0.000 claims abstract description 89
- 238000003384 imaging method Methods 0.000 claims abstract description 21
- 210000003484 anatomy Anatomy 0.000 claims abstract description 16
- 238000000034 method Methods 0.000 claims description 118
- 238000002591 computed tomography Methods 0.000 claims description 56
- 230000009467 reduction Effects 0.000 description 26
- 238000013459 approach Methods 0.000 description 13
- 238000012545 processing Methods 0.000 description 9
- 230000008901 benefit Effects 0.000 description 7
- 230000000694 effects Effects 0.000 description 7
- 230000006870 function Effects 0.000 description 5
- 208000011231 Crohn disease Diseases 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 4
- 230000033001 locomotion Effects 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000002595 magnetic resonance imaging Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000009206 nuclear medicine Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000001303 quality assessment method Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000007476 Maximum Likelihood Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000000747 cardiac effect Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000004321 preservation Methods 0.000 description 2
- 230000004224 protection Effects 0.000 description 2
- 230000000306 recurrent effect Effects 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 230000002792 vascular Effects 0.000 description 2
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 1
- 208000017667 Chronic Disease Diseases 0.000 description 1
- 208000012671 Gastrointestinal haemorrhages Diseases 0.000 description 1
- 208000022559 Inflammatory bowel disease Diseases 0.000 description 1
- 101100328463 Mus musculus Cmya5 gene Proteins 0.000 description 1
- 101100043621 Mus musculus Ssxa1 gene Proteins 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 208000003386 Radiation-Induced Neoplasms Diseases 0.000 description 1
- 235000019892 Stellar Nutrition 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000003187 abdominal effect Effects 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000003292 diminished effect Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 210000005069 ears Anatomy 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000002594 fluoroscopy Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000011478 gradient descent method Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 235000003642 hunger Nutrition 0.000 description 1
- 229910052740 iodine Inorganic materials 0.000 description 1
- 239000011630 iodine Substances 0.000 description 1
- 230000005865 ionizing radiation Effects 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000008155 medical solution Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 210000004197 pelvis Anatomy 0.000 description 1
- 230000010412 perfusion Effects 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000000264 spin echo pulse sequence Methods 0.000 description 1
- 230000037351 starvation Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- JLYXXMFPNIAWKQ-UHFFFAOYSA-N γ Benzene hexachloride Chemical compound ClC1C(Cl)C(Cl)C(Cl)C(Cl)C1Cl JLYXXMFPNIAWKQ-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/06—Diaphragms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4035—Arrangements for generating radiation specially adapted for radiation diagnosis the source being combined with a filter or grating
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/405—Source units specially adapted to modify characteristics of the beam during the data acquisition process
-
- 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/48—Diagnostic techniques
- A61B6/482—Diagnostic techniques involving multiple energy imaging
-
- 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/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
-
- 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/54—Control of apparatus or devices for radiation diagnosis
-
- 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/54—Control of apparatus or devices for radiation diagnosis
- A61B6/542—Control of apparatus or devices for radiation diagnosis involving control of exposure
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4064—Arrangements for generating radiation specially adapted for radiation diagnosis specially adapted for producing a particular type of beam
- A61B6/4085—Cone-beams
-
- 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/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
Definitions
- CT computed tomography
- Automated tube current modulation or dose modulation can adj st the x-r y tube output mAs to the patient's anatomy in order to maintain a desired level of image quality ( See Refs. 10 and 1 1). Greater mAs can be used in. areas with greater x-ray attenuating path-lengths such as the shoulders or the pelvis, and lower mAs can be used in regions containing less attenuating material such as the lungs.
- Iterative reconstruction procedures can solve an optimization problem (See Ref, 13) in which the cost function can he designed to reduce noise or artifacts resulting from reductions in x-ray tube current. Iterative reconstruction procedures can solve an optimization problem (See Ref, 13) in which the cost function can he designed to reduce noise or artifacts resulting from reductions in x-ray tube current. Iterative reconstruction procedures can solve an optimization problem (See Ref, 13) in which the cost function can he designed to reduce noise or artifacts resulting from reductions in x-ray tube current. Iterative reconstruction procedures can solve an optimization problem (See Ref, 13) in which the cost function can he designed to reduce noise or artifacts resulting from reductions in x-ray tube current. Iterative reconstruction procedures can solve an optimization problem (See Ref, 13) in which the cost function can he designed to reduce noise or artifacts resulting from reductions in x-ray tube current. Iterative reconstruction procedures can
- multiplicative dose saving effects See ef. 18
- four 30% dose reduction strategies can in principle be used in combination to effect a 76% dose reduction relative to the starting point.
- existing methods cannot achieve routine subraiHisievert scanning in the major ity' of patients, and sobmillisievert scans remain elusive in almost ail abdominal applications.
- Signature low dose results often involve small patients scanned at low kVp in the chest (e.g., where inherent x-ray attenuation can be low), but doses can be substantively higher in the abdomen and for typical patient sizes,
- CS compressed sensing
- CT presents favorable conditions for the application of CS since (a) medical images are naturally compressible by various transforms which lead to sparse representations, and (b) discarding projections can result in low-value streaking artifacts that add incoherently to the image representation.
- nonlinear reconstruction techniques e.g., iterative reconstruction or nonlinear filtering methods
- CT techniques there are currently no CT techniques that can undersample the x-ray projection data itself, which can be essential for making effective use of techniques such as compressed sensing to further reduce CT radiation dose beyond current limits. Tims, it may be beneficial to provide an exemplary system, method, and computer accessible medium that can modulate x-ray beam intensity for use with CS reconstruction, and which can overcome at least some of the deficiencies described herein above.
- Exemplary system, method and computer accessible medium can facilitate orders of magnitude radiation dose reductions in CT. While current radiation dose reduction techniques hi CI " can rely on techniques to reduce the x-ray tube output throughout the CT scan acquisition, the exemplary system, method and computer accessible medium can utilize nonlinear reconstruction techniques to reconstruct images with, full image quality, information content and diagnostic utility from undersampled dataseis, in which a significant fraction of x-ray exposures have been eliminated.
- an exemplary imaging system can include, for example, a source first arrangement that can provide x-ray radiation, and a collimator second arrangement that can be provided in a path of the x-ray radiation, and can be configured to undersampSe the radiation which can be forwarded to the portion(s) of the anatomical structure.
- a third hard ware arrangement can be configured to receive a further radiation from the portton(s) that can be based on the undersam ied radiation, j0016f S i certain exemplary embodiments of the present disclosure
- the collimator second arrangement can include a multi-hole collimator, which can be configured to move in such a wa that different pathways between the source first arrangement and the third arrangement are blocked or left unobstructed at different times.
- the collimator second arrangement can be moved by a piezoelectric drive, a stepper motor or a pneumatic drive.
- the second arrangement can be further configured undersample the radiation, and provide resultant radiation to the portionis) of an anatomical structure that can be based on the coiHmated radiation, and the third arrangement
- a hardware fourth arrangement can be configured to generate an image(s) of the portion(s) using a compressed sensing reconstruction procedure, or other suitable reconstruction procedure, as a function of the resultant radiation.
- the system can .further include a gantry bore and a gantry table located in the gantry bore, and the source first arrangement can be located on the gantry bore.
- the collimator second arrangement can be configured to reduce a number of the radiations along a detector array fan or cone angle between. the source first arrangement and the third arrangement.
- the second arrangement can be further configured to isrulersampk the radiation, and provide resultant radiation to the portion(s) of tire anatomical structure that can be based on the undersarnpled radiation, and an imaging fourth arrangement can be configured to receive further radiations from the portion(s), and generate an rraage(s) of the portion(s) based on the further radiations, hi certain exemplary embodiments of the present disclosure, the source first arrangement and the third arrangement ca be part of a mulrisHce computed tomography ("CT”) scanner.
- CT computed tomography
- ⁇ 001.8f in a further exemplary embodiment of the present disclosure can be an imaging system, which can nclude, for example, a source first arrangement that can be configured to provide an x-ray radiation, a multi-hole collimator second arrangement, and a hardware third arrangement.
- the second and third arrangements can be provided in a path of the radiation, and the second arrangement can be provided between the first and third arrangements.
- a further exemplary embodiment of the present disclosure can be a system, method and computer-accessible medium fo imaging a portion(s) of a anatomical structure that can include, tor example moving a multi-hole collimator arrangement in and out of a path of an x-ray radiation provided by an x-ray source arrangement; the path can extend to or through a radiation receiving arrangement
- An image of the portion(s) can be generated based on further radiation received from the portion(s) by the radiation receiving arrangement using a compressed sensing reconstruction procedure.
- the further radiation can be associated with the x-ray radiation.
- the image can be generated using an inverse problem solution algorithm incorporating a sparsity transform, which can include wavelets, finite differences, curvetets, or other elements know in the art.
- an apparatus can be provided that can include a multi-hole collimator arrangement configured to undersaraple an x- ray radiation provided from a source arrangement.
- the particular tmdersampled pattern can different at adjacent angles of a gantry bore.
- ⁇ 002 J In a further exemplary embodiment of the present disclosure is an exemplary system.. method, and computer-accessible medium for generating an image(s) of a portion(s) of an anatomical structure including incoherently interrupting an x-ray radiation provided, by an x-ray source arrangement that is directed at the at least one portion to generate a resultant x-ray radiation having a particular ursdersample pattern, receiving, at a radiation receiving
- Figures LA- 1 C are exemplary images generated with three distinct reconstruction methods: A) traditional filtered back protection, B) SAFIRE. a commercial iterative
- Figures 2 A and 2B are exemplary images of ⁇ -on!y undersamplmg and ⁇ - ⁇ undersampl ng according to an exemplary embodiment of the present disclosure
- Figure 3 is a set. of images illustrating reconstruction of phantom data according to an exemplary embodiment of the present disclosure
- Figure 4 is an exemplary illustration of an exemplary compressed sensing system according to an exemplary embodiment of the present disclosure:
- Figures 5 A and SB are exemplary schematics illustrating movement of an exemplary collimator arrangement according to an exemplary embodiment of the present disclosure
- (00291 figure 6 is an exemplary image illustrating further exemplary results provided by the exemplary compressed sensing system according to an exemplar embodiment of the present disclosure
- figure 7 is an exemplary image illustrating further exemplary results provided by the exemplary compressed sensing system according to an exemplary embodiment of the present disclosure
- Figure 8 is an exemplary image illustrating a comparison of conventional imaging techniques and the exemplary compressed sensing technique according to an exemplary embodiment of the present disclosure
- Figure 9 is an image of an exemplary pulsed x-ray tube according to an exemplary embodiment of the present disclosure.
- figure 10 is an flow diagram illustrating a method for generating an image of an anatomical structure according to an exemplary embodiment of the present disclosure: and.
- Figure 1 1 is an illustration of an exemplary block diagram of an exemplary system in accordance with certain exemplary embodiments of the present disclosure.
- the exemplary system, method, and computer-accessible medium according to an exemplary embodiment of the present disclosure can be used for rapid and incoherent interruption of the x-ray source on the CT gantry which, when combined with exemplary sparsiiy-basei! CS reconstruction procedures, can enable reconstruction of high-quality images from a markedly reduced number of projections.
- various methods can be used to interrupt the x-ray source including an exemplary multi-hole collimator and/or an exemplary pulsed x-ray tube.
- Figures 1A-IC show experimental results demonstrating the improved performance of reconstruction with incoherent, undersampling as compared with state-of-the-art iterative reconstruction with reduced tube current using the exemplary system, method, and computer- accessible medium.
- Moving collimators in particular, and potentially pulsed X-ray tubes as well, can he able to be implemented as a modular hardware upgrade which, together with software upgrades to reconstruction procedures, can be deployed on existing CT systems, thereby significantly expanding the impact of the work, and facil itating rapid and broad dissemination of the method.
- Tailored sparsiiy-enforcing compressed-sensing reconstruction procedures can. be utilized to take maximal advantage of the resulting undersampling.
- the exemplary system, method, and computer-accessible medium can enable orders-of-magnitude rather than incremental dose reductions, which can represent a game-changing strategy in CT dose reduction methods.
- the exemplary system, method, and computer-accessible medium can be broadly applicable to all anatomy from head to toe, and can also be used in combination with existing approaches for further dose reductions.
- Exemplary image reconstruction procedures can be used to enforce sparsity subject to data consistency constraints. Given the acquisition model y ⁇ Ax, where y can be the
- A can be the acquisition matrix ot: forward model which can include the scanner geometry, and x can be the image to reconstruct, the reconstruction for the exemplary system, method, and computer-accessible medium can be given by, for example:
- T can be a sparsifying transform under which x can have a sparse representation
- ⁇ can be a regularization parameter that can weight the contribution of the sparsity-term (e.g., right-hand term) relative to the data consistency term (e.g., left-hand term).
- Various procedures can be used to solve Eq. (I) in a very efficient way using, for example, non-linear conjugate gradient, iterative soft -thresholding, and/or gradient-descent methods.
- the most straightforward undersampiing pattern can eliminate data at selected gantry rotation angles ⁇ (e.g., ⁇ undersampiing).
- ⁇ e.g., ⁇ undersampiing
- an improved ' undersampiing pattern can be used, where a different subset of z-axis positions can be irradiated by the incident x-ray beam for each gantry rotation angle, resulting in a ⁇ - ⁇ undersampiing pattern (e.g., Figure 2A).
- the subset of z-slices for each ⁇ angle can be chosen randomly to maximize incoherence.
- This type of undersampiing pattern can be analogous to an exemplary undersampiing pattern for dynamic MR! studies, where a different radial k-space undersampiing pattern can be used for each time point ( See Ref, 28), if the time-dimension is replaced with the slice-dimension.
- Higher-dimensional undersampling patterns can present favorable conditions for compressed sensing since aliasing artifacts can be effectively distributed over a larger space with lower- intensity, and sparsity can be exploited along the extra dimensions ( See Refs, 19 and 23).
- ⁇ - ⁇ unde.rsampl.ing can distribute artifacts incoherently throughout the full three- dimensional ("3D") data volume, thereby improvin performance.
- an 8-fold undersampling factor on numerical phantom data generated can be simulated using sequentially-acquired (e.g., non- helical) cone-beam geometry (e.g.. Figure 2B).
- the fully-sampled data set included. 800 projections and 192 slices.
- the matrix A was given by the forward operator of the FD reconstruction procedure ( See Ret 29).
- Sparse reconstruction was performed using 3 - wavelets (e.g., 4-tap, Daubmüs basis functions) as the sparsifying transform.
- ⁇ - ⁇ undersampling improves performance as compared with ⁇ -only undersampling, due to the increased incoherence, and the ability to exploit extra sparsity along the slice dimension.
- the Q-z undersampling approach can be used, for the exemplary system, method, and computer-accessible medium, which can be implemented using the moving multi-hole collimator described below.
- the exemplary system, method, and computer-accessible medium can be compared against Sinogram Affirmed Iterative Reconstruction ("SAFIRE") , a state-of-the-art reduced- tube-current technique employing iterative reconstruction and provided by Siemens . Medical Solutions on its commercial CT scanners. Data for the comparison were generated from phantom experiments performed on Siemens SOMATO Definition Edge CT scanner using a step and shoot protocol Comparison was performed using a high dose-reduction factor (e.g., 8- foid or about 87.5% reduction in this case, which can be much higher than the usual 20-30% used in typical clinical scans) in order to demonstrate more-than-incremental gains enabled by the exemplary system, method, and computer-accessible medium.
- SAFIRE Sinogram Affirmed Iterative Reconstruction
- S AFIRE data were acquired with tube-current of about 25 mAs and reconstruction was performed online on the scanner.
- images acquired by sampling fully at about 200 mAs were forward projected, and retrospectively undersampled, using the ⁇ - ⁇ undersampiing pattern described above.
- the forward and adjoint acquisition operator were implemented i the EMreeon.
- software package See ef. 30 ⁇ which can. support the use of full cone beam projections.
- Reconstruction for the exemplary system, method, and computer-accessible medium was performed using 3D finite differences as the sparsifying transform, which can be equivalent to minimization of 3D total variation. The reconstruction procedure can enforce data consistency in the original cone-beam space.
- Figure 3 is a set of images illustrating reconstruction of phantom data acquired with full dose (e.g., 200 mAs, no undersampiing) using liitered-back projection ("FBP"), an 8-fold reduced tube-current (e.g., 25 mAs, no undersampiing) using SAFIRE, and an 8-fold retrospective undersampiing (e.g., 200 mAs, 8-fold ⁇ - ⁇ undersampiing) using the exemplary system, method, and computer-accessible medium.
- full dose e.g., 200 mAs, no undersampiing
- SAFIRE liitered-back projection
- 8-fold retrospective undersampiing e.g., 200 mAs, 8-fold ⁇ - ⁇ undersampiing
- the leftmost col m shows full field of view reconstructions and the columns on the right show zoomed panels to emphasize the preservation of small (305) and low-contrast (310) features in the exemplary system, method, and computer- accessible medium compared to the exemplary system, method, and computer-accessible medium.
- Figure 3 illustrates he ' higher reconstruction performance of the exemplary system, method, and computer-accessible medium in comparison, to SAFIRE, particularly in the preservation of high-resolution and low-contrast features, which can be buried in the residual noise of the SAFIRE reconstruction.
- the nonlinear reconstruction .1 the exemplary sysiem.
- the method, and computer-accessible medium can also denoise the .reconstructed images, presenting superior performance even compared with, the full-dose FBP reconstruction, which cm be a linear reconstruction approach. This can be due to the feet that the phantom can consist of piece wise smooth features with truly sparse representations in the finite difference domain, and eonse quently can enable improved compressed sensing reconstructions
- ears mechanically disrupt an x-ray beam close to the .x-ray source, which can be in, or adjacent to, an x-ray collimator arrangement.
- the x-ray beam can be narrowly confined such that ver slight linear motions of a heterogeneously attenuating structure can effecti vely modulate the x-ray beam intensity at a high frequency as compared to the standard CT gantry rotation time.
- the z-axis can be the position of a CT table 405 and a patient 410 (e.g., from head to foot) normal to the gantry bore, and the gantry angle ⁇ can be 0° zero when the x-ray source can be above the CT table 405, and 180° when the x-ray source can be below the CT table 405.
- a moving multi-hole collimator can be utilized. which can expose a different subset of z-axis slices per gantry rotation angle, and thus create an incoherent undersamplmg pattern along 0 and z.
- FIGS. 5A and SB are schematics illustrating movement of the exemplary collimator arrangement (e.g., a moving multi-hole collimator for ⁇ - ⁇ incoherent undersamplmg).
- the muitihole coMiraator can be mounted on the standard fixed collimator housing and needs only- linear motion. As the gantry rotates, for example, from ⁇ 1 ::: 0o ( Figure 5 A) to ⁇ 2 :::: 30" ( Figure 5B), the muitihole collimator can move linearly along the z direction, exposing a different subset of slices at each gantry rotati n angle ⁇ .
- Different slices can be irradiated at adjacent gantry angles, yielding the exemplary Q-z undersampling pattern in Figure. 3.
- a particular pattern of x-ray radiation can be emitted (e.g., the areas denoted by element 505 can facilitate the passing of x-rays while the areas denoted by element 510 can facilitate the blocking of the x-rays).
- This exemplary undersamplmg pattern can result in a different complement of radial protections for each z-axis position, and can enabie a Ml 3D compressed sensing reconstruction, which can have improved performance over 2D compressed sensing due to increased
- Image reconstruction can be performed by utilizing 3D image sparsity using a 3D transform (e.g., wavelets, finite differences, and curvelets).
- a 3D transform e.g., wavelets, finite differences, and curvelets.
- the exemplary system, method and computer accessible medium can interrupt the x-ray output at the x-ray source 515.
- exemplary source collimator designs that, can incoherently reduce the number of x-ray projections along the detector array fan angle ⁇ (See e.g., Figure 4), which can. create an additional imdersarapling dimension by exposing different detector element combinations within each slice.
- the incoherent undersampling in the (e.g., 6, ⁇ , z) space can facilitate orders of magnitude redactions in radiation exposure without any effect on the CT image quality.
- the exemplary system, method and computer accessible medium can move the multi-hole collimator 520 via, for example, a piezoelectric drive, a stepper motor, a pneumatic drive, or other similar drive mechanisms.
- the incident x-ray flux ca be reduced at the source, which can result in a diminished radiation exposure to the patient.
- the exemplary system, method and computer accessible medium can also be fully compatible with and or utilize commonly used- dose- reduction technologies (e.g., including z-axis and in-plane tube current modulation, reduced kVp imaging, adaptive colHmation and/or other hardware and software methods currently in use), facilitating multiplicative dose reductions by capturing the benefits of multiple imaging techniques simultaneously.
- the x-ray radiation can also be interrupted by pulsing the beam itself at its exemplary source, for example, by turning the source on and off in a desired undersarapling pattera (e.g..
- I S using the exemplary x-ray tube of Figure 9 I S using the exemplary x-ray tube of Figure 9).
- a standard x-ray tube electrons can be emitted at the cathode by heating, and then accelerated by a high voltage potential difference (e.g. 120 kV) towards the anode where they can generate x-rays on impact.
- a high voltage potential difference e.g. 120 kV
- the flow of electrons must be interrupted, as opposed to interrupting the flow of x-rays following x-ray generation.
- Various means of pulsing X-ray tubes may be devised, and the undersampied patterns thereby generated can be used for compressed sensing reconstruction, radiation dose reduction, etc.
- the exemplary data undersamp!ing model used to generate preliminary data in Figures 2A-4 assumes full interruption of the X-ray beam in the non-selected detector-rows, whereas partial blockage can be possible at certain angles and collimator positions.
- the exemplary system, method, and computer-accessible medium, according to an exemplary embodiment of the present disclosure, can also incorporate no effects of the collimator holes on beam profile, such as penumbra! effects or focal spot blurring.
- Various exemplary sparse reconstruction procedures can be used For robust operation with the realistic undersatnpHng model.
- beam profile effects can be included in the acquisition matrix.
- A which can be used in the iterations to improve data consistency.
- E ⁇ j. (1) can also be modified to a maximum- likelihood estimation ("MLE") framework in consideration of CT noise statistics (e.g., shot- noise, Poisson distribution), and this exemplary approach can be compared to an. exemplary least-squares franiework,
- MLE maximum- likelihood estimation
- convex optimization procedures can be used with fester convergence, such as, fast iterative soft thresholding ("FiSTA”) ( See Ref. 41), primal-dual approaches, and various other approaches.
- FiSTA fast iterative soft thresholding
- the undersampling model can he applied to fully-sampled acquisitions performed on the American College of Radiology ("ACR") CT accreditation phantom (e.g., Gammex 464) (See Ref. 42).
- ACR American College of Radiology
- the ACR phantom can be composed of 4 separate sections to individually evaluate slice width, low-contrast resolution, noise and uniformity, and high-contrast resolution.
- the image quality assessment can focus on detectabs lity of low-contrast and high-resolution features, which can be most affected by the usual artifacts in compressed sensing reconstructions at. high undersampimg factors.
- Quantitative image quality assessment can be perforated, using fully-sampled images as ground- truth to compute root mean square error ("RMSE") and structural similarity inde (“SSI ”) ( See Ref. 43). Additionally, qualitative impressions of the detectability of low-contrast and high- resolution image features with respect, to the fully-sampled images can be catalogued.
- RMSE root mean square error
- SSI structural similarity inde
- Figure 10 is an exemplary method for generating an image of an anatomical structure.
- x-ray radiation can be generated.
- the x-ra radiation can be incoherently interrupted.
- the incoherently interrupted x-ray radiation is received at an x-ray receiving arrangement.
- the image can be generated, for example, using a compressed sensing procedure.
- Figure 1 1 shows a block diagram of an exemplary embodiment of a system according to the present disclosure.
- exemplary procedures in accordance with the present disclosure described herein can be performed by a processing arrangement and/or a computing arrangement 1 102.
- processing/computing arrangement 1102 can be, for example entirely or a part of, or include, but not limited to, a computer/processor 1 1 4 that can include, for example one or more microprocessors, and use instructions stored on a computer-accessible medium ⁇ e.g., RAM, ROM, hard drive, or other storage device).
- a computer-accessible medium e.g., RAM, ROM, hard drive, or other storage device.
- a computer-accessible medium 1106 e.g., as described herein above, a storage device such as a hard disk, floppy disk, memory stick, CD- ROM, RAM, ROM, etc, or a collection thereof
- the computer-accessible medium 1 106 can contain executable instructions 1 108 thereon, in addition or alternatively, a storage arrangement i i 10 can be provided separatel from the computer-accessible medi um ⁇ 106, which can provide the instructions to the processing arrangement .1.102 so as to configure the processing arrangement to execute certain exemplary procedures, processes, and methods, as described herein above, for example,
- the exemplary processing arrangement 1 102 can be provided with or include an input/ouipui arrangement 1 1 14, which can include, for example a wired network, a wireless network, the internet, an intranet, a data collection probe, a sensor, etc.
- the exemplary processing arrangement 1 102 can be in communication with an exemplary display arrangement 1 1 12, which, according to certain exemplary embodiments of the present disclosure, can be a touch-screen configured for inputting information to the processing arrangement in addition to outputting information from the processing arrangement, for example.
- the exemplary display 1 1 12 and/or a storage arrangement 1 1 10 can be used to display and or store data in a user-accessible format and/or user-readable format.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
An imaging system for imaging a portion(s) of an anatomical structure can be provided. For example, a source first arrangement can provide x-ray radiation, and a multi-hole collimator second arrangement can be provided in a path of the x-ray radiation, and can be configured to undersample the radiation beam which can be forwarded to the portion(s) of the anatomical structure. A third hardware arrangement can be configured to receive a further radiation from the portion(s) that can be based on the undersampled radiation.
Description
SYSTEM, METHOD AND COMPUTER ACCESSIBLE MEDIUM FOR
MODULATING X-RAY BEAM INTENSITY
CROSS-REFERENCE TO RELATED APPLICATIO '
{000 f This application relates to and claims priority from U.S. Patent Application No.
1/804,415, filed on March 22, 2013, the entire disclosure of which is incorporated herein, by reference.
FIELD OF THE DISCLOSURE
{0002} The present disclosure relates generally to computed tomography ("CT"), and more specifically, to exemplary systems, methods and computer-accessible medium for modulating x- ray beam intensity in a C'T imaging system,
BACKGROUND INFORMATION
(0003$ Rapidly increasing utilization of X-ray CT has heightened concerns about the coiiective radiation exposure to the population as a whole, and about potential risks to patients undergoing recurrent imaging for chronic conditions or persistent complaints. These concerns have generated, a great deal of attention to practical radiation dose reduction strategies. For example, a recent addition to the commercially-available dose-reduction arsenal is on iterative reconstruction, which uses reduced tube current, and applies nonlinear procedures to denoise the resulting images. However, such approaches can be limited by the minimum useful tube current (e.g., the photon starvation problem), and in practice, only moderate dose reductions of up to
about 30-40% below the prior state of the art. Recent advances in detector technology can offer similar reductions.
(Θ0Θ4| CT has received a lot of scrutiny because of its relatively high radiation dose per exam. While i t comprises only abou t 17% of all medical imaging procedures, i t produces approximately half of the population's medical radiation exposure, wit nuclear medicine contributing approximately one quarter of the collective dose to the population, and fluoroscopy and conventional x-ray exams accounting for the remainder (See Refs. 5 and 6).
jOGOSJ Various strategies have been proposed to reduce radiation exposure to the population as a whole, and to individual patients (See Refs. 7 and 8). Due to challenges of education and dissemination, the most successful and widely implemented dose-reduction measures can be those that perform with little need for end-use optimization. Examples of currentl used measures are below.
[0000J Exemplary automated tube current modulation: Automated tube current modulation or dose modulation, can adj st the x-r y tube output mAs to the patient's anatomy in order to maintain a desired level of image quality ( See Refs. 10 and 1 1). Greater mAs can be used in. areas with greater x-ray attenuating path-lengths such as the shoulders or the pelvis, and lower mAs can be used in regions containing less attenuating material such as the lungs.
[OO07J Adaptive Collimation: The development of wider detector arrays facilitating rapid anatomic coverage has had the unintended consequence of increasing z-overscatmmg, in which the anatomy above and below the desired scan range can be inadvertently irradiated.
implementation of adaptive collimators that open and close asymmetrical ly as tlie patient can be
7
moved through the gantry can mostly eliminate this z-overscanning, with dose savings on the order of 10-35%, depending on the detector configuration and scan length (See Ref. 12).
{"0008} Iterative Reconstruction; iterative reconstruction procedures can solve an optimization problem (See Ref, 13) in which the cost function can he designed to reduce noise or artifacts resulting from reductions in x-ray tube current. Iterative reconstruction procedures can
transform back and forth between the raw data "projection-space" and the image domain with each iteration, but, as this approach can be computationally expensive, most manufacturers have implemented more rapid shortcut procedures designed to achieve similar ends. Implementation details vary between manufacturers, but generall invol ve a variety of procedures to shift some of the iterative "correction steps" into the raw data or image domains, combined with advanced modeling of the CT acquisition system, and nonlinear image filtering, to reduce noise in homoaeneous resions while attempting to preserve anatomic edae information. The noise reduction effected b iterative reconstruction techniques can enable substantial reductions in radiation exposure. The somewhat different noise texture of the resultant images can need some acclimatization on the part of radiologists, but in practical routine clinical use, dose reductions on die order of 20-40% can be typical (See Ref. 16).
|0009| Reduced kVp imaging: In many applications, reduced kVp can result in substantial dose savings. The greatest benefits can be realized in vascular imaging applications, as iodine attenuates lower energy x-rays more strongly than higher-energy x-rays. Adoptio of low kVp approaches has been hastened through approaches such as CarekV which can direct the CT scanner to select the kVp setting that can result in the lowest dose to the patient, while meeting user-selected image quality constraints and inherent x-ray tube output limits.
|ft01 | Detector technology advances; CT manufacturers have spent substantial effort optimizing the materials (e.g. Siemens UFC detectors) and construction of their detector elements to maximize the conversion of incident x-rays to meaningful signal. Most recently, analog to digital conversion functionality has been moved to chips located on the detector elements themselves {e.g., Siemens Stellar detectors), thus reducing electronic noise. This has been shown to permit dose reductions on the order of about 30%.
[001 !} Each of the dose reduction strategies above lias the potential to reduce a dose incrementally from previous techniques. Some of these techniques may only be appropriate i select circumstances (e.g. dose reductions from tube current modulation for small patients and low kVp for vascular imaging particularly in smaller patients) . In other cases, multiple
techniques can be used in combination, producing multiplicative dose saving effects (See ef. 18) (e.g. four 30% dose reduction strategies can in principle be used in combination to effect a 76% dose reduction relative to the starting point). However, existing methods cannot achieve routine subraiHisievert scanning in the major ity' of patients, and sobmillisievert scans remain elusive in almost ail abdominal applications. Signature low dose results often involve small patients scanned at low kVp in the chest (e.g., where inherent x-ray attenuation can be low), but doses can be substantively higher in the abdomen and for typical patient sizes,
f00i2| An alternative approach to reduce radiation dose in CT is the application of compressed sensing ("CS") techniques (See Refs. 19 and 20). CS can exploit image compressibility to reconstruct an image with full information from a reduced set of incoherent measurements sampled below the Nyquist-Shannon rate. The reconstruction can aiso be iterative, but in contrast to the iterative reconstruction techniques mentioned above, the cost
function can be designed to enforce sparsity in the space where the image can. be compressible. Significant reductions in the number of measurements can be accompiished if the image can be sparse in a known transform domain (e.g., sparsity condition) and if the resulting aliasing artifacts can add incoherently to the sparse image representation (e.g., incoherence condition). CT presents favorable conditions for the application of CS since (a) medical images are naturally compressible by various transforms which lead to sparse representations, and (b) discarding projections can result in low-value streaking artifacts that add incoherently to the image representation.
|<H}13| Some current radiation dose reduction techniques in CT can. rely on 'techniques to reduce the x-ra tube output throughout the CT scan acquisition (e.g., without rapid.
modulations}. In some cases, nonlinear reconstruction techniques (e.g., iterative reconstruction or nonlinear filtering methods) can be used to reduce noise in the resulting images. However, there are currently no CT techniques that can undersample the x-ray projection data itself, which can be essential for making effective use of techniques such as compressed sensing to further reduce CT radiation dose beyond current limits. Tims, it may be beneficial to provide an exemplary system, method, and computer accessible medium that can modulate x-ray beam intensity for use with CS reconstruction, and which can overcome at least some of the deficiencies described herein above.
SUMMARY OF EXEMPLARY EMBODIMENTS
|0O14| Exemplary system, method and computer accessible medium according to an exemplary embodiment of the present disclosure can facilitate orders of magnitude radiation
dose reductions in CT. While current radiation dose reduction techniques hi CI" can rely on techniques to reduce the x-ray tube output throughout the CT scan acquisition, the exemplary system, method and computer accessible medium can utilize nonlinear reconstruction techniques to reconstruct images with, full image quality, information content and diagnostic utility from undersampled dataseis, in which a significant fraction of x-ray exposures have been eliminated. jOOlSf These and other objects of the present disclosure can be achieved an exemplary imaging system that can include, for example, a source first arrangement that can provide x-ray radiation, and a collimator second arrangement that can be provided in a path of the x-ray radiation, and can be configured to undersampSe the radiation which can be forwarded to the portion(s) of the anatomical structure. A third hard ware arrangement can be configured to receive a further radiation from the portton(s) that can be based on the undersam ied radiation, j0016f S i certain exemplary embodiments of the present disclosure, the collimator second arrangement can include a multi-hole collimator, which can be configured to move in such a wa that different pathways between the source first arrangement and the third arrangement are blocked or left unobstructed at different times. The collimator second arrangement can be moved by a piezoelectric drive, a stepper motor or a pneumatic drive. The second arrangement can be further configured undersample the radiation, and provide resultant radiation to the portionis) of an anatomical structure that can be based on the coiHmated radiation, and the third arrangement A hardware fourth arrangement can be configured to generate an image(s) of the portion(s) using a compressed sensing reconstruction procedure, or other suitable reconstruction procedure, as a function of the resultant radiation.
{0017| In certain exemplary embodiments of the present disclosure, the system can .further include a gantry bore and a gantry table located in the gantry bore, and the source first arrangement can be located on the gantry bore. The collimator second arrangement can be configured to reduce a number of the radiations along a detector array fan or cone angle between. the source first arrangement and the third arrangement. The second arrangement can be further configured to isrulersampk the radiation, and provide resultant radiation to the portion(s) of tire anatomical structure that can be based on the undersarnpled radiation, and an imaging fourth arrangement can be configured to receive further radiations from the portion(s), and generate an rraage(s) of the portion(s) based on the further radiations, hi certain exemplary embodiments of the present disclosure, the source first arrangement and the third arrangement ca be part of a mulrisHce computed tomography ("CT") scanner.
{001.8f in a further exemplary embodiment of the present disclosure can be an imaging system, which can nclude, for example, a source first arrangement that can be configured to provide an x-ray radiation, a multi-hole collimator second arrangement, and a hardware third arrangement. The second and third arrangements can be provided in a path of the radiation, and the second arrangement can be provided between the first and third arrangements.
{0019\ In a further exemplary embodiment of the present disclosure can be a system, method and computer-accessible medium fo imaging a portion(s) of a anatomical structure that can include, tor example moving a multi-hole collimator arrangement in and out of a path of an x-ray radiation provided by an x-ray source arrangement; the path can extend to or through a radiation receiving arrangement An image of the portion(s) can be generated based on further radiation received from the portion(s) by the radiation receiving arrangement using a compressed sensing
reconstruction procedure. The further radiation can be associated with the x-ray radiation. The image can be generated using an inverse problem solution algorithm incorporating a sparsity transform, which can include wavelets, finite differences, curvetets, or other elements know in the art.
|0020f in a further exemplary embodiment of the present disclosure, an apparatus can be provided that can include a multi-hole collimator arrangement configured to undersaraple an x- ray radiation provided from a source arrangement. The particular tmdersampled pattern can different at adjacent angles of a gantry bore.
{002 J In a further exemplary embodiment of the present disclosure is an exemplary system.. method, and computer-accessible medium for generating an image(s) of a portion(s) of an anatomical structure including incoherently interrupting an x-ray radiation provided, by an x-ray source arrangement that is directed at the at least one portion to generate a resultant x-ray radiation having a particular ursdersample pattern, receiving, at a radiation receiving
arrangement, the resultant x-ray radiation provided from the portion(s), and generating the image(s) of the portiou(s) based on the resultant x-ray radiation received at the radiation receiving arrangement
|0022| These and other objects, features and advantages of the exemplary embodiments of the present disclosure will become apparent upon reading the following detailed description of the exemplary embodiments of the present disclosure, when taken in conjunction with the appended claims.
BRIEF DESCRIPTION OF THE DRAWINGS
[9023] Fiinher objects, features, and advantages of the present disclosure will become apparent from the following detailed description taken in conjunction with the accompanying Figures showing illustrative embodiments of the present disclosure, in which:
[9024] Figures LA- 1 C are exemplary images generated with three distinct reconstruction methods: A) traditional filtered back protection, B) SAFIRE. a commercial iterative
reconstruction, with tube current reduced uniformly to yield 1/8 the radiation dose of the traditional case, and C) use o the exemplary system, method, and computer-accessible medium according to the present disclosure to reconstruct an image from an undersampled dataset omitting 7/8 of the traditional proj ections and corresponding to 1/8 the radiation dose of the traditional case;
[9025] Figures 2 A and 2B are exemplary images of θ-on!y undersamplmg and θ-ζ undersampl ng according to an exemplary embodiment of the present disclosure;
[9026] Figure 3 is a set. of images illustrating reconstruction of phantom data according to an exemplary embodiment of the present disclosure;
[9027] Figure 4 is an exemplary illustration of an exemplary compressed sensing system according to an exemplary embodiment of the present disclosure:
[0028] Figures 5 A and SB are exemplary schematics illustrating movement of an exemplary collimator arrangement according to an exemplary embodiment of the present disclosure; (00291 figure 6 is an exemplary image illustrating further exemplary results provided by the exemplary compressed sensing system according to an exemplar embodiment of the present disclosure;
|ft030| figure 7 is an exemplary image illustrating further exemplary results provided by the exemplary compressed sensing system according to an exemplary embodiment of the present disclosure;
[9031) Figure 8 is an exemplary image illustrating a comparison of conventional imaging techniques and the exemplary compressed sensing technique according to an exemplary embodiment of the present disclosure;
[9032) Figure 9 is an image of an exemplary pulsed x-ray tube according to an exemplary embodiment of the present disclosure;
[0033) figure 10 is an flow diagram illustrating a method for generating an image of an anatomical structure according to an exemplary embodiment of the present disclosure: and. j0034| Figure 1 1 is an illustration of an exemplary block diagram of an exemplary system in accordance with certain exemplary embodiments of the present disclosure.
[0035) Throughout the drawings, the same reference numerals and characters, unless otherwise stated, are used to denote like features, elements, components, or portions of the illustrated embodiments. Moreover, while the present disclosure will now be described in detail with reference to the figures, it is done so in connection with the illustrati e embodiments and is not limited by the particular embodiments illustrated in the figures or the appended claims.
DETAILED DESCRIPTIO OF EXEMPLARY EMBODIMENTS
[0036) The exemplary system, method, and computer-accessible medium according to an exemplary embodiment of the present disclosure can be used for rapid and incoherent interruption of the x-ray source on the CT gantry which, when combined with exemplary
sparsiiy-basei! CS reconstruction procedures, can enable reconstruction of high-quality images from a markedly reduced number of projections. In particular, various methods can be used to interrupt the x-ray source including an exemplary multi-hole collimator and/or an exemplary pulsed x-ray tube.
[9937] Previous CS techniques suggested the reconstruction of an. image from a reduced number of radial projections, and subsequent work further explored this idea using real CT data (See Ref. 21). Recent work has also demonstrated , in simulations, that reducing the number of projections can achieve a higher radiation dose reduction than reducing the tube current (See Ref. 22).
[9038] la contrast to previous work, which was not effective in reducing radiation dosages. Figures 1A-IC show experimental results demonstrating the improved performance of reconstruction with incoherent, undersampling as compared with state-of-the-art iterative reconstruction with reduced tube current using the exemplary system, method, and computer- accessible medium.
[9039] In Magnetic Resonance Imaging (MRI),. CS using angularly- rtdersam led radial k- space acquisitions, has enabled orders-of-magnitu.de accelerations without degradation of MR image quality ( See Ref, 28). Radial MRI and CT share similar data acquisition schemes, and therefore, similar orders-of-magnifude undersampling factors can be expected in CT, which can lead to the target. sub-mSv acquisitions in all parts of the body. One additional advantage of CS approaches can be that they can be largely complementary to existing dose reduction approaches, and can be used in combination with other exemplary technique for greater improvements.
|ft040| E em lar 1 incoherent undersampling procedure can be a key to unlocking the benefits of CS for CT, Two exemplary procedures of incoherent undersampling can be employed: (a) exemplary moving multi-hole collimators, and (b) exemplary pulsed X-ray tube designs.
Moving collimators in particular, and potentially pulsed X-ray tubes as well, can he able to be implemented as a modular hardware upgrade which, together with software upgrades to reconstruction procedures, can be deployed on existing CT systems, thereby significantly expanding the impact of the work, and facil itating rapid and broad dissemination of the method. Tailored sparsiiy-enforcing compressed-sensing reconstruction procedures can. be utilized to take maximal advantage of the resulting undersampling.
10041 j The exemplary system, method, and computer-accessible medium, according to an exemplary embodiment of the present disclosure, can enable orders-of-magnitude rather than incremental dose reductions, which can represent a game-changing strategy in CT dose reduction methods. The exemplary system, method, and computer-accessible medium can be broadly applicable to all anatomy from head to toe, and can also be used in combination with existing approaches for further dose reductions.
Exemplary Image .Reco truetion Framework
[9042.1 Exemplary image reconstruction procedures can be used to enforce sparsity subject to data consistency constraints. Given the acquisition model y~Ax, where y can be the
undersampted set of projections, A can be the acquisition matrix ot: forward model which can include the scanner geometry, and x can be the image to reconstruct, the reconstruction for the exemplary system, method, and computer-accessible medium can be given by, for example:
where T can be a sparsifying transform under which x can have a sparse representation, and λ can be a regularization parameter that can weight the contribution of the sparsity-term (e.g., right-hand term) relative to the data consistency term (e.g., left-hand term). Various procedures can be used to solve Eq. (I) in a very efficient way using, for example, non-linear conjugate gradient, iterative soft -thresholding, and/or gradient-descent methods.
f0043| Much of the existing work on sparse reconstruction for CT has employed decimated parallel-beam data (e.g., for an undersampiing factor of 5 only every th projection can be used). This exemplary theoretical model can be useful to show proof of concept, but this type of rapid x-ray beam interruption at frequencies on the order of 1-5 kHz can be impossible with current X-ray tube and scanner designs. First, realistic scanner geometry can be included.
Assuming a multisiice CT scanner geometry, the most straightforward undersampiing pattern can eliminate data at selected gantry rotation angles Θ (e.g., Θ undersampiing). Second, as multiple slices along the z-axis {e.g.. cone-beam direction of the detector array) can be acquired for each Θ, an improved 'undersampiing pattern can be used, where a different subset of z-axis positions can be irradiated by the incident x-ray beam for each gantry rotation angle, resulting in a θ-ζ undersampiing pattern (e.g., Figure 2A). The subset of z-slices for each Θ angle can be chosen randomly to maximize incoherence. This type of undersampiing pattern can be analogous to an exemplary undersampiing pattern for dynamic MR! studies, where a different radial k-space undersampiing pattern can be used for each time point ( See Ref, 28), if the time-dimension is replaced with the slice-dimension.
|ft044| Higher-dimensional undersampling patterns can present favorable conditions for compressed sensing since aliasing artifacts can be effectively distributed over a larger space with lower- intensity, and sparsity can be exploited along the extra dimensions ( See Refs, 19 and 23). in CT, θ-ζ unde.rsampl.ing can distribute artifacts incoherently throughout the full three- dimensional ("3D") data volume, thereby improvin performance. In order to compare the performance of θ-only 'undersampling and θ-ζ undersampling, an 8-fold undersampling factor on numerical phantom data generated can be simulated using sequentially-acquired (e.g., non- helical) cone-beam geometry (e.g.. Figure 2B). The fully-sampled data set included. 800 projections and 192 slices. The matrix A was given by the forward operator of the FD reconstruction procedure ( See Ret 29). Sparse reconstruction, was performed using 3 - wavelets (e.g., 4-tap, Daubechtes basis functions) as the sparsifying transform. As can be seen, θ-ζ undersampling improves performance as compared with θ-only undersampling, due to the increased incoherence, and the ability to exploit extra sparsity along the slice dimension. Based on the above, the Q-z undersampling approach can be used, for the exemplary system, method, and computer-accessible medium, which can be implemented using the moving multi-hole collimator described below.
Exemplary Results
|0O4Sf The exemplary system, method, and computer-accessible medium can be compared against Sinogram Affirmed Iterative Reconstruction ("SAFIRE") , a state-of-the-art reduced- tube-current technique employing iterative reconstruction and provided by Siemens. Medical Solutions on its commercial CT scanners. Data for the comparison were generated from
phantom experiments performed on Siemens SOMATO Definition Edge CT scanner using a step and shoot protocol Comparison was performed using a high dose-reduction factor (e.g., 8- foid or about 87.5% reduction in this case, which can be much higher than the usual 20-30% used in typical clinical scans) in order to demonstrate more-than-incremental gains enabled by the exemplary system, method, and computer-accessible medium. S AFIRE data were acquired with tube-current of about 25 mAs and reconstruction was performed online on the scanner. To simulate the exemplary system, method, and computer-accessible medium undersampiing, images acquired by sampling fully at about 200 mAs were forward projected, and retrospectively undersampled, using the θ-ζ undersampiing pattern described above. The forward and adjoint acquisition operator were implemented i the EMreeon. software package ( See ef. 30} which can. support the use of full cone beam projections. Reconstruction for the exemplary system, method, and computer-accessible medium was performed using 3D finite differences as the sparsifying transform, which can be equivalent to minimization of 3D total variation. The reconstruction procedure can enforce data consistency in the original cone-beam space.
[0046} Figure 3 is a set of images illustrating reconstruction of phantom data acquired with full dose (e.g., 200 mAs, no undersampiing) using liitered-back projection ("FBP"), an 8-fold reduced tube-current (e.g., 25 mAs, no undersampiing) using SAFIRE, and an 8-fold retrospective undersampiing (e.g., 200 mAs, 8-fold θ-ζ undersampiing) using the exemplary system, method, and computer-accessible medium. The leftmost col m shows full field of view reconstructions and the columns on the right show zoomed panels to emphasize the preservation of small (305) and low-contrast (310) features in the exemplary system, method, and computer- accessible medium compared to the exemplary system, method, and computer-accessible
medium. Figure 3 illustrates he 'higher reconstruction performance of the exemplary system, method, and computer-accessible medium in comparison, to SAFIRE, particularly in the preservation of high-resolution and low-contrast features, which can be buried in the residual noise of the SAFIRE reconstruction. The nonlinear reconstruction .1» the exemplary sysiem. method, and computer-accessible medium can also denoise the .reconstructed images, presenting superior performance even compared with, the full-dose FBP reconstruction, which cm be a linear reconstruction approach. This can be due to the feet that the phantom can consist of piece wise smooth features with truly sparse representations in the finite difference domain, and eonse quently can enable improved compressed sensing reconstructions
Exemplary Collimator Arrangement
{00471 The exemplary system, method and computer accessible medium, according to an exemplary embodiment of the present disclosure, ears mechanically disrupt an x-ray beam close to the .x-ray source, which can be in, or adjacent to, an x-ray collimator arrangement. For example, near the source, the x-ray beam can be narrowly confined such that ver slight linear motions of a heterogeneously attenuating structure can effecti vely modulate the x-ray beam intensity at a high frequency as compared to the standard CT gantry rotation time.
[9048.1 As shown in Figure 4, the z-axis can be the position of a CT table 405 and a patient 410 (e.g., from head to foot) normal to the gantry bore, and the gantry angle Θ can be 0° zero when the x-ray source can be above the CT table 405, and 180° when the x-ray source can be below the CT table 405. Using a trrultislice CT scanner 400 using which x-ray attenuation data can foe acquired from several slices per rotation, a moving multi-hole collimator can be utilized.
which can expose a different subset of z-axis slices per gantry rotation angle, and thus create an incoherent undersamplmg pattern along 0 and z.
{"0049} Figures 5A and SB are schematics illustrating movement of the exemplary collimator arrangement (e.g., a moving multi-hole collimator for θ-ζ incoherent undersamplmg). The muitihole coMiraator can be mounted on the standard fixed collimator housing and needs only- linear motion. As the gantry rotates, for example, from Θ1 :::0o (Figure 5 A) to Θ2::::30" (Figure 5B), the muitihole collimator can move linearly along the z direction, exposing a different subset of slices at each gantry rotati n angle Θ. Different slices can be irradiated at adjacent gantry angles, yielding the exemplary Q-z undersampling pattern in Figure. 3. For example, in Figures 5A and SB, a particular pattern of x-ray radiation can be emitted (e.g., the areas denoted by element 505 can facilitate the passing of x-rays while the areas denoted by element 510 can facilitate the blocking of the x-rays).
f0050| This exemplary undersamplmg pattern can result in a different complement of radial protections for each z-axis position, and can enabie a Ml 3D compressed sensing reconstruction, which can have improved performance over 2D compressed sensing due to increased
incoherence and sparsity. (See e.g., exemplary results shown in Figures 6 and 7). Image reconstruction can be performed by utilizing 3D image sparsity using a 3D transform (e.g., wavelets, finite differences, and curvelets).
|0O5I| The exemplary system, method and computer accessible medium, according to exemplary embodiments of the present disclosure, can interrupt the x-ray output at the x-ray source 515. including exemplary source collimator designs that, can incoherently reduce the number of x-ray projections along the detector array fan angle Φ (See e.g., Figure 4), which can.
create an additional imdersarapling dimension by exposing different detector element combinations within each slice. Using the exemplary system and method, the incoherent undersampling in the (e.g., 6, Φ, z) space can facilitate orders of magnitude redactions in radiation exposure without any effect on the CT image quality.
|0052| The exemplary system, method and computer accessible medium, according to exemplary embodiments of the present disciosisre, can move the multi-hole collimator 520 via, for example, a piezoelectric drive, a stepper motor, a pneumatic drive, or other similar drive mechanisms. The incident x-ray flux ca be reduced at the source, which can result in a diminished radiation exposure to the patient. The exemplary system, method and computer accessible medium can also be fully compatible with and or utilize commonly used- dose- reduction technologies (e.g., including z-axis and in-plane tube current modulation, reduced kVp imaging, adaptive colHmation and/or other hardware and software methods currently in use), facilitating multiplicative dose reductions by capturing the benefits of multiple imaging techniques simultaneously.
|0053| As shown in Figure S, a comparison is provided of conventional imaging techniques and the exemplary compressed sensing technique according to an exemplary embodiment of the present disclosure, can have an incoherence along both Θ and zs while conventional imaging techniques likely cannot achieve incoherence along both Θ and z. Exemplary Puked X-ray Tubes
{ 054J The x-ray radiation can also be interrupted by pulsing the beam itself at its exemplary source, for example, by turning the source on and off in a desired undersarapling pattera (e.g..
I S
using the exemplary x-ray tube of Figure 9). In a standard x-ray tube, electrons can be emitted at the cathode by heating, and then accelerated by a high voltage potential difference (e.g. 120 kV) towards the anode where they can generate x-rays on impact. In order to pulse an exemplary x- ray tube, the flow of electrons must be interrupted, as opposed to interrupting the flow of x-rays following x-ray generation. Various means of pulsing X-ray tubes may be devised, and the undersampied patterns thereby generated can be used for compressed sensing reconstruction, radiation dose reduction, etc.
Exemplary Realisiic UndersampUng Mo el
[9055] The exemplary data undersamp!ing model used to generate preliminary data in Figures 2A-4 assumes full interruption of the X-ray beam in the non-selected detector-rows, whereas partial blockage can be possible at certain angles and collimator positions. The exemplary system, method, and computer-accessible medium, according to an exemplary embodiment of the present disclosure, can also incorporate no effects of the collimator holes on beam profile, such as penumbra! effects or focal spot blurring.
Exemplary Sparse Reconstruction Procedures
[0056| Various exemplary sparse reconstruction procedures can be used For robust operation with the realistic undersatnpHng model. For example, beam profile effects can be included in the acquisition matrix. A, which can be used in the iterations to improve data consistency. By including the expected profile of partially attenuated projections and differential blurring as a function of time / gantry-angle / collimator-motion explicitly in forward and. adjoint operators, an
improvement can be achieved in the sharpness and the accuracy of the exemplary
reconstructions. The data consistency term in E<j. (1) can also be modified to a maximum- likelihood estimation ("MLE") framework in consideration of CT noise statistics (e.g., shot- noise, Poisson distribution), and this exemplary approach can be compared to an. exemplary least-squares franiework,
j0057| In order to represent realistic anatomy efficiently, several 3D sparsifying transforms, such as orthogonal wavelets, finite differences (e.g., total variation), and geometrically oriented wavelets such as curvelets (see Ref. 37) and shearlets (See Ref. 38), and various combinations of these transforms, can be evaluated based on the compression ratio they can achieve (e.g., as evaluated first in fully-sampled images) and on computational burden. The resulting sparse reconstruction procedures and sparsifying transforms can be implemented using parallel computing techniques on graphical processor units ("GPU") for fast reconstruction. The major limitation on reconstruction speed can be the evaluation of the forward and backward acquisition operator A in each, iteration. Prior work has shown that GPUs can be particularly efficient for these types of operations (see Refs. 39 and 40). Besides parallel implementation of the
reconstruction procedures, convex optimization procedures can be used with fester convergence, such as, fast iterative soft thresholding ("FiSTA") ( See Ref. 41), primal-dual approaches, and various other approaches.
{"0058 j For reconstruction performance evaluation, the undersampling model can he applied to fully-sampled acquisitions performed on the American College of Radiology ("ACR") CT accreditation phantom (e.g., Gammex 464) (See Ref. 42). The ACR phantom can be composed of 4 separate sections to individually evaluate slice width, low-contrast resolution, noise and
uniformity, and high-contrast resolution. Specifically, the image quality assessment can focus on detectabs lity of low-contrast and high-resolution features, which can be most affected by the usual artifacts in compressed sensing reconstructions at. high undersampimg factors.
Quantitative image quality assessment can be perforated, using fully-sampled images as ground- truth to compute root mean square error ("RMSE") and structural similarity inde ("SSI ") ( See Ref. 43). Additionally, qualitative impressions of the detectability of low-contrast and high- resolution image features with respect, to the fully-sampled images can be catalogued.
|¾059] Figure 10 is an exemplary method for generating an image of an anatomical structure. At procedure 1005, x-ray radiation can be generated. At procedure 1010, the x-ra radiation can be incoherently interrupted. At procedure 1015, the incoherently interrupted x-ray radiation is received at an x-ray receiving arrangement. At procedure 1020, the image can be generated, for example, using a compressed sensing procedure.
0060j Figure 1 1 shows a block diagram of an exemplary embodiment of a system according to the present disclosure. For example, exemplary procedures in accordance with the present disclosure described herein can be performed by a processing arrangement and/or a computing arrangement 1 102. Such processing/computing arrangement 1102 can be, for example entirely or a part of, or include, but not limited to, a computer/processor 1 1 4 that can include, for example one or more microprocessors, and use instructions stored on a computer-accessible medium {e.g., RAM, ROM, hard drive, or other storage device).
[00611 As shown in Figure 11 , for example a computer-accessible medium 1106 (e.g., as described herein above, a storage device such as a hard disk, floppy disk, memory stick, CD- ROM, RAM, ROM, etc, or a collection thereof) can be provided (e.g., in communication with
the processing arrangement .102). The computer-accessible medium 1 106 can contain executable instructions 1 108 thereon, in addition or alternatively, a storage arrangement i i 10 can be provided separatel from the computer-accessible medi um ί 106, which can provide the instructions to the processing arrangement .1.102 so as to configure the processing arrangement to execute certain exemplary procedures, processes, and methods, as described herein above, for example,
[9062) Further, the exemplary processing arrangement 1 102 can be provided with or include an input/ouipui arrangement 1 1 14, which can include, for example a wired network, a wireless network, the internet, an intranet, a data collection probe, a sensor, etc. As shown in Figure 1 1, the exemplary processing arrangement 1 102 can be in communication with an exemplary display arrangement 1 1 12, which, according to certain exemplary embodiments of the present disclosure, can be a touch-screen configured for inputting information to the processing arrangement in addition to outputting information from the processing arrangement, for example. Further, the exemplary display 1 1 12 and/or a storage arrangement 1 1 10 can be used to display and or store data in a user-accessible format and/or user-readable format.
(0063$ The foregoing merely illustrates the principles of the disclosure. Various
modifications and alterations to the described embodiments will be apparent to those skilled in the art in view of the teachings herein, it will thus be appreciated that those skilled in the art will be able to devise numerous systems, arrangements, and procedures which, although not explicitly shown or described herein, embody the principles of the disclosure and can be thus within the spirit and scope of the disclosure. Various different exemplary embodiments can be used together with one another, as well as interchangeably therewith, as should be understood by
those ha ving ordinary skill in the art. In addition, certain terms used in the present disclosure, including the specification, drawings and claims thereof, can be used synonymoosly in certain instances, including, but not limited to, e.g., data and information, it should be understood that, while these words, and/or other words that can be syn mous to one another, can he used synonymously herein, that there can be instances when such words can be intended to not be used synonymously. Further, to the extent that the prior art knowledge has not been, explicitly incorporated by reference herein above, it is explicitly incorporated herein in its entirety. All publications referenced are incorporated herein by reference in thei entireties.
EXEMPLARY REFERENCES
The following references are hereby incorporated by reference in their entirety.
L Benchmark report CT. Des Plains; IMV Medical Information Division; 2007,
2. Brenner DJ, Hall EJ. Computed tomography - an increasing source of radiation exposure. N Engl J Med. 20G7;357(22):2277-84.
3. Sodickson A, Baeyens PF, Andriole KP, Prevedello LM, Nawfel RD, Hanson R, Khorasani R. Recurrent CT, cumulative radiation exposure, and associated radiation- induced cancer risks from CT of adults. Radiology. 2009;25l(l): 175-84.
4. Berrington de Gonzalez A. Mahesh , Kim KP, Bhargavan , Lewis R, Mettler F, Land C.
Projected cancer risks from computed tomographic scans performed in the United States in 2007. Arch intern Med. 2009;I69(22):2071-7.
5. Mettler FAJ, Bhargavan M, Faulkner K, Giliey DB, Gray JE, Ibbott GS, Lipoti JA, Mahesh M., McCrohan J'L. Stabin MG, Thomadsen BR, Yoshizunri TT. Radiologic and nuclear medicine studies in the United States and worldwide: frequency, radiation dose, and
comparison with other radiation sources- 1950-2007. Radiology, 2009;253(2):520-31.
6. Faze! R, Krumholz HM, Wang Y, Ross JS, Chen J, Ting H . Exposure to low-dose ionizing radiation from medical imaging procedures, N Engl J Med. 2009;361 (9):849~57.
7. Amis ESJr, Butler PF, Applegate K.E, Birnbaum SB, Brateman LP, Hevezi JM. Mettler FA, Morra RL, Pentecost MJ, Smith GG, Strauss KJ, Zeman RK. American College of Radiology white paper on radiation dose in medicine. J Am Coll Radiol. 2007;4(5);272-84.
8. Sodickson A . Strategies for reducing radiation exposure in multi-detector row CT. Radiol Clin North Am. 2OJ 2;50(1):M4.
9. http:// ^w.healthcare.s½mem
iO0t'''wcm idc sieiuens hwem-hwem ssxa,
root/wcm' dk^groups/publ^
de-to-Sow-dose-2013-00849104.pdf.
10, Kalendar WA, Computed Tomography. 3rd ed. Erlangen, Germany: Pubiicis Publishing;
201 1.
1 1. McCollough CI!, Bruesewitz MR, Kofler JM. CT dose reduction and dose management tools: Overview of available options- Radiographics, 20O6;26(2):5O3~ 12.
12. Deak PD, Langner O, Lell M, Kalender WA. Effects of adaptive section colliraation on
patient radiation dose in multisection spiral CT. Radiology. 2009;252(1): 40-7,
13. Shepp LA, Vardi Y, Maximum likelihood reconstruction for emission tomography. IEEE Trans Med Imag. 1982; 1 (2): 1 1.3-22.
14. Grant , Flohr T. Iterative Reconstruction in Image Space (IRIS). Siemens Healthcare White Paper: A9115-101492. 2010.
15. Grant K, Raupach R. SAFIRE: Sinogram Affirmed Iterative Reconstruction. Siemens
Healthcare White Paper. 2011.
16. Hara A 5 Baden RGS Silva AC, ujak JL, Lawder HJ, Pavlicek W, 'iterative reconstruction technique for reducing body radiation dose at CT: feasibility study. AIR Am J Roentgenol. 2009; 193(3): 764-71.
17. Sodickson A, Weiss M. Effects of patient size on radiation dose reduction and image quality in low-kVp CT pulmonary angiography performed with reduced IV contrast dose. Emerg Radiol 2052: 59{5}:437-45.
18. McCol ugh CH. Chen GF1. Kalender W, Leiig S, Samei E. Taguehi K, Wang G, Yu L5 Pettigre RI, Achieving routine submillisievert CT scanning; report from the summit on management of radiation dose In CT. Radiology. 2012;264(2):567-80,
19. Cande'S E. Romberg J, Tao T. Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information. IEEE Trans Inf Theory. 2006;52:489-509.
20. Dono'ho D. Compressed sensing. IEEE Trans Inf Theory. 2006;S2: 1289-306.
21. Chen GH, Tang J, Leng S. Prior image constrained compressed sensing (PfCCS): a method to accurately reconstruct dynamic CT images from highly undersampled projection dat sets. Med Phys. 2008;35{2):660-3.
22. Brown K, Koeh!er T, Bergner F, Bippus R, Brendei B, Karl W, Singh S, Padole A, Do S,
Sparse sampling for CT dose red«ction20l3: The 12th International Meeting on Fully Three- Dimensional Image Reconstruction in Radiology and Nuclear Medicine, p. 428-31.
23. Otazo R, Kim D, Axel L, Sodickson DK. Combination of compressed sensing and parallel imaging for highly accelerated first-pass cardiac perfusion MRL Magn Reson Med.
2010;64(3):767-76. PubMed Central PMCID: PMC2932824.
24. Feng L, Otazo R, Jung H, Jensen Jl I, Ye JC, Sodickson DK, Kim D. Accelerated cardiac T2 mapping using breath-hold multiecho fast spin-echo pulse sequence with k-t FOCUSS. Magn Reson Med. 201 1 ;65(ό): Ι661 -9. PubMed Central PMCID: PMC3097270.
25. Kim D, Dyvorne HA, Otazo R, Feng L, Sodickson DR. Lee VS. Accelerated phase-contrast cine MR! using k-t SPARSE-SENSE, Magn Reson Med. 2012;67(4): 1054-64. PubM ed
Centra! PMCID: PMC3306497.
26. Feng L. Srichai M. Lim PL, Harrison A, King W, Adluru G, Dibella EV, Sodickson DK, Otazo R, Kim D. Highly accelerated real-time cardiac eioe MR] using k-t SPARSE-SENSE. Magn Reson Med. 2012;7§( l);6 -74.
27. Knoll F, Clason C, Bredies K, Uecker M, Stoll.be.rger R. Parallel imaging with nonlinear reconstruction using variational penalties. Magn Reson Med. 2012;67(! )i34-4 L
28. Chandarana H, Feng L, Block TK, Rosenkraatz AB, Lim RP, Babb JS, Sodickson DK, Otazo R. Free-breathing contrast-enhanced multiphase MRI of the liver using a combinatio of compressed sensing, parallel imaging, and golden-angle radial sampling. Invest Radiol.
2013;48(1 ): 10-6.
29, Feidka p LA, Davis EC, Kress JW. Practical cone-beam algorithm. J Opt Soc Amer.
1984;Α1 (6):612- .
30. Koesters T, Sehaefers KP,, Wuehbeling F. EMRECO'N: An expectation maximization based image reconstruction framework for emission tomography data 201 1 : IEEE Nuclear Science
Symposium, p. 4365-68.
31. Loftus EV, Jr., Schoenfeld P. Sandborn WJ. The epidemiology and natural history of Crohn's disease in population-based patient cohorts from North America: a systematic review.
Aliment Pharmacol Ther. 2002; 16U ):51-60,
32. Wu YW, Tang YH, Hao NX, Tang CY, Miao F. Crohn's disease: CT enterography
manifestations before and after treatment. Eur J Radiol. 2012;81(l):52-9.
33. Paulsen SR, Huprich JE, Hara AK. CT enterography: noninvasive evaluation of Crohn's disease and obscure gastrointestinal bleed. Radiol Clin North Am. 2007;45(2):303-l 5.
34. Hara AK, Swartz PG. CT enterography of Crohn's disease. Abdom imaging. 2009;34(3):289-
95.
35. Hupricii JE, Fletcher JG. CT enterography: principles, technique and utility in Crohn's
disease. Eur J Radiol. 2009;69(3): 393-7.
36, Herrinton Li, Lin L, Lewis JD, Griffin PM. Allison J. Incidence and prevalence of
inflammatory bowel disease in a Northern California managed care organization. 1996-2002. Am j Gastroenterol. 2008; 103(8): 1998-2006.
37. Candes E, Demanet L, Donoho D, Ying LX. Fast discrete curvelet transforms. Multiscie Model Sim. 2006;5(3):86I-99.
38. Guo K, Lahate D. Optimally sparse multidimensional representation using shearlets. SIAM J Math Ami 2007;39:298-318.
39. Knoll F, linger M, Diwoky , Clason C, Pock T, Stollberger R. Fast reductio of
undersampiing artifacts in radial MR angiography with 3D total variation on graphics hardware. MAGMA. 2010;23(2); 10344.
40. Freiberger M, Knoll F, Bredies K, Scharfetter H, Stollberger R. The AGILE library for image reconstruction in biomedical sciences using graphics card hardware acceleration. Computing in Science and Engineerin 2013; 15:34-44.
41. Beck A, Teboulle M. Fast gradient-based algorithms for constrained total variation image denoising and debt luring problems. IEEE Trans Image Process. 2009; 1 §( i ] ):2419-34.
42. McCoilough CH, Bruesewitz MR, McNitt-Gray MF, Bush K, Ruckdeschel T, Payne JT, Brink JA:. Zeman RK. The phantom portion of the American College of Radiology (ACR)
computed tomography (CT) accreditation program; practical tips, artifact examples, and pitfalls to avoid. Med Pliys. 2O04;31(9):2423-42.
Wang Z, Bovik AC, Sheikh HR, Simonce!li EP, Image quality assessment: From error visibility to structural similarity. IBEE Trans Image Proc. 2004; 1.3(4):6(KM 2.
Claims
WHAT IS CLAIMED IS:
1 , A method generating at least one image of at least one portion of an anatomical structure, comprising:
incoherently interrupting an x-ray radiation beam, provided by an x-ray source
arrangement that is directed at the at least one portion to generate a resultant x-ray radiation beam having a particular nndersarapied pattern;
receiving, at a radiation receiving arrangement, the resultant x-ray radiation provided from the at least one portion; and
generating the at least one image of the at least one portion based on the resultant x-ray radiation received at the radiation receiving arrangement.
2. The method of claim .1.. wherein the resultant x-ray radiation is generated using a .multi- hole collimator arrangement. 3. The method of claim 2, further comprising moving the multi-hole collimator arrangement across a path of the x-ray radiation to generate the resultant x-ray radiation.
4. The method of claim 3, further comprising moving the collimator arrangement using at least one of a piezoelectric drive, a stepper motor or a pneumatic drive.
5. The method of claim I , wherein the resultant radiation is generated using at least one pulsed x-ray tube.
6. The method of claim 1, wherein th e x-ray source arrangement and the radiation receiving arrangement are part, of a multislice computed tomography ("CT") scanner,
7. The method of claim 6, further comprising utilizing a gantry bore and a gantry table located in the gantry bore, wherein the x-ray source arrangement is located on the gantry bore.
8. The method of claim 7, wherein the particular undersampied pattern is different at adjacent angles of the gan try bore.
9. The method of claim 1, wherein the image is generated using a compressed sensing reconstruction procedure .
10. The method of c laim 9, w herein the compressed sensing reconstruction procedure includes involves a sparsifymg transform including at least one of wavelets, finite differences, or curve lets.
1 1.. An imaging system, comprising:
a source first arrangement providing an x-ray radiation; and
a second arrangement configured to incoherently interrupt the x-ray radiation provided by the source first arrangement to generate a resultant x-ray radiation having a particular undersampied pattern which is forwarded to at least one portion of an anatomical structure; and
a third hardware arrangement confia redto receive a further radiation from the at least one portion that is based on the resultant x-ray radiation.
12. The system 'of claim 11 , wherein the second arrangement is a multi-hole collimator arrangemen
13, The system of claim 12, wherein the multi-hole collimator arrangement is configured to move across a path, of the x-ray radiation. 14 , The system of claim 13, further comprising at least one of a piezoelectric drive, a stepper motor or a pneumatic drive which is configured to move the multi-hole collimator arrangement.
15. The system of claim 12, wherein the multi-hole collimator arrangement is configured to reduce a number of the radiations along a detector array fan or cone angle between the source first arrangement and the third hardware arrangement.
16. The system of claim 11 , wherein the second, arrangement is ai least one pulsed x-ray tube.
.17. The system of claim 1 1 , wherein the source first arrangement and the third hardware arrangement are part of a mu!tisHce computed tomography ("CT") scanner.
18. The system of claim 17, further comprising a gantry bore and a gantry table located in the gantry bore, wherein the source first arrangement is located on the gantry bore,
1.9. The system of claim 18, wherein the particular undersampied pattern is different at adjacent angles of the gantr bore.
20, The system of claim i 1, wherein the third arrangement includes a hardware fourth arrangement which is configured to generate at least one image of the at least one portion using a compressed sensing reconstruction procedure as a function of the further radiation.
21 , The system of claim 20, wherein the compressed sensing reconstruction procedure includes involves a sparsifying transform including at least one of wavelets, finite differences, or curvelets.
22, A. non-transitory computer-accessible medium having stored thereon computer- executable instructions for generating at least one image of at least one portion of an anatomical structure, wherein, when a computer hardware arrangement executes the instructions, the computer arrangement is configured to perform procedures comprising:
incoherently .interrupting an x-ray radiation beam provided by an x-ray source
arrangement that is directed at the at least one portion to generate a resultant x-ray radiation, beam having a particular undersampied pattern;
receiving, at a radiation receiving arrangement, the resultant x-ray radiation provided from the at least one portion; and
generating the at. least one image of the at least one portion based on the resultant x-ray radiation received at the radiation receiving arrangement.
23. The method of claim 22, wherein the resultant x-ray radiation is generated using a raiiiti- hole collimator arrangement.
24. The method of claim 23, further comprising moving the multi-hole collimator arrangement across a path of the x-ray radiation to generate the resultant x-ray radiation.
25. The method of claim 24, farther comprising moving the collimator arrangement using at least one of a piezoelectric drive, a stepper motor or a pneumatic drive.
26. The method of claim 22, wherein the resultant radiation is generated using a t least one pulsed x-ray tube.
27. The method of claim 22 wherein the x-ray source arrangement and the radiation receiving arrangement are part of a tmtltisHce computed tomography ("CT") scanner.
28. The method of claim 27, further comprising utilizing a gantry bore and a gantry table located in the gantry bore, wherein the x-ray source arrangement is located on the gantry bore.
29, "The method of claim 28, wherein the particular undersampied pattern is different at adjacent angles of the gantry bore. 30. The method of claim 22, wherem the image is generated using a compressed sensing reconstruction procedure.
3.1. The method of claim 30, wherein the compressed sensing reconstruction procedure includes involves a sparstfying transform including at least one of wavelets, finite differences, or curvelets.
32. An apparatus, comprising:
a multi-hole coilimator arrangement configured to incoherently interrupt an x-ray radiation provided from a source arrangement so as to generate a resultant radiation having a particular undersampied pattern.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP14771044.6A EP2976626A4 (en) | 2013-03-22 | 2014-03-24 | Computer accessible medium for modulating x-ray beam intensity |
US14/861,312 US10092253B2 (en) | 2013-03-22 | 2015-09-22 | System, method, and computer accessible medium for modulating X-ray beam intensity |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201361804415P | 2013-03-22 | 2013-03-22 | |
US61/804,415 | 2013-03-22 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US14/861,312 Continuation-In-Part US10092253B2 (en) | 2013-03-22 | 2015-09-22 | System, method, and computer accessible medium for modulating X-ray beam intensity |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2014153561A1 true WO2014153561A1 (en) | 2014-09-25 |
Family
ID=51581538
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2014/031609 WO2014153561A1 (en) | 2013-03-22 | 2014-03-24 | Computer accessible medium for modulating x-ray beam intensity |
Country Status (3)
Country | Link |
---|---|
US (1) | US10092253B2 (en) |
EP (1) | EP2976626A4 (en) |
WO (1) | WO2014153561A1 (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11399780B2 (en) * | 2016-11-07 | 2022-08-02 | Rensselaer Polytechnic Institute | Methods and systems for stationary computed tomography |
US11276208B2 (en) | 2017-06-12 | 2022-03-15 | The Research Foundation For The State University Of New York | System, method and computer-accessible medium for ultralow dose computed tomography image reconstruction |
CN113367718A (en) * | 2021-05-31 | 2021-09-10 | 明峰医疗系统股份有限公司 | CT scanning device and method for sparse angle sampling |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4819259A (en) * | 1985-09-20 | 1989-04-04 | Kabushiki Kaisha Toshiba | Rotary anode type x-ray apparatus |
US5812629A (en) * | 1997-04-30 | 1998-09-22 | Clauser; John F. | Ultrahigh resolution interferometric x-ray imaging |
US20060023842A1 (en) * | 2004-07-28 | 2006-02-02 | General Electric Company | Single-leaf X-ray collimator |
US20110006768A1 (en) * | 2009-07-09 | 2011-01-13 | Lei Ying | Systems and methods for accelerating the acquisition and reconstruction of magnetic resonance images |
US20120265050A1 (en) * | 2011-04-04 | 2012-10-18 | Ge Wang | Omni-Tomographic Imaging for Interior Reconstruction using Simultaneous Data Acquisition from Multiple Imaging Modalities |
Family Cites Families (45)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE2756659A1 (en) * | 1977-12-19 | 1979-06-21 | Philips Patentverwaltung | ARRANGEMENT FOR DETERMINING THE ABSORPTION DISTRIBUTION |
US4506374A (en) * | 1982-04-08 | 1985-03-19 | Technicare Corporation | Hybrid collimator |
FR2526575A1 (en) * | 1982-05-04 | 1983-11-10 | Thomson Csf | RADIOLOGICAL IMAGE PROCESSING METHOD FOR CORRECTING SAID IMAGE OF DEFECTS DUE TO DIFFUSED RADIATION |
US6987831B2 (en) * | 1999-11-18 | 2006-01-17 | University Of Rochester | Apparatus and method for cone beam volume computed tomography breast imaging |
US6480565B1 (en) * | 1999-11-18 | 2002-11-12 | University Of Rochester | Apparatus and method for cone beam volume computed tomography breast imaging |
US6396902B2 (en) * | 2000-07-31 | 2002-05-28 | Analogic Corporation | X-ray collimator |
US7072436B2 (en) * | 2001-08-24 | 2006-07-04 | The Board Of Trustees Of The Leland Stanford Junior University | Volumetric computed tomography (VCT) |
US7145981B2 (en) * | 2001-08-24 | 2006-12-05 | The Board Of Trustees Of The Leland Stanford Junior University | Volumetric computed tomography (VCT) |
DE10222701C1 (en) * | 2002-05-22 | 2003-10-30 | Siemens Ag | X-ray dosage distribution measuring method for computer tomography apparatus using detector elements of X-ray detector for measuring dosage distribution during adjustment of beam stop |
US6876718B2 (en) * | 2003-06-27 | 2005-04-05 | Ge Medical Systems Global Technology Company, Llc | Scatter correction methods and apparatus |
US7382853B2 (en) * | 2004-11-24 | 2008-06-03 | General Electric Company | Method and system of CT data correction |
EP1780676A1 (en) * | 2005-10-25 | 2007-05-02 | GSF-Forschungszentrum für Umwelt und Gesundheit GmbH | Imaging method and device with dual reading scanner |
US7463712B2 (en) * | 2006-05-18 | 2008-12-09 | The Board Of Trustees Of The Leland Stanford Junior University | Scatter correction for x-ray imaging using modulation of primary x-ray spatial spectrum |
WO2008127385A2 (en) * | 2006-10-24 | 2008-10-23 | Thermo Niton Analyzers Llc | Apparatus for inspecting objects using coded beam |
WO2009082531A1 (en) * | 2007-12-20 | 2009-07-02 | Wisconsin Alumni Research Foundation | Method for prior image constrained image reconstruction |
WO2009082419A1 (en) * | 2007-12-20 | 2009-07-02 | Wisconsin Alumni Research Foundation | Method for dynamic prior image constrained image reconstruction |
US8218907B2 (en) * | 2008-01-14 | 2012-07-10 | Wisconsin Alumni Research Foundation | Method for prior image constrained progressive image reconstruction |
EP2349008B1 (en) * | 2008-11-26 | 2015-02-25 | Wisconsin Alumni Research Foundation | Method for prior image constrained image reconstruction in cardiac cone beam computed tomography |
US8310233B2 (en) * | 2009-02-18 | 2012-11-13 | Mayo Foundation For Medical Education And Research | Method for image reconstruction from undersampled medical imaging data |
WO2010121043A2 (en) * | 2009-04-15 | 2010-10-21 | Virginia Tech Intellectual Properties, Inc. | Exact local computed tomography based on compressive sampling |
EP2250967A1 (en) * | 2009-05-13 | 2010-11-17 | Helmholtz Zentrum München Deutsches Forschungszentrum für Gesundheit und Umwelt (GmbH) | Tomography apparatus and tomography method |
US8111893B2 (en) * | 2009-06-09 | 2012-02-07 | Wisconsin Alumni Research Foundation | Method for dynamic prior image constrained image reconstruction |
AU2010262843B2 (en) * | 2009-06-19 | 2015-07-16 | Viewray Technologies, Inc. | System and method for performing tomographic image acquisition and reconstruction |
US8204172B1 (en) * | 2010-03-17 | 2012-06-19 | General Electric Company | System and method of prior image constrained image reconstruction using short scan image data and objective function minimization |
US8483463B2 (en) * | 2010-05-19 | 2013-07-09 | Wisconsin Alumni Research Foundation | Method for radiation dose reduction using prior image constrained image reconstruction |
WO2012009725A1 (en) * | 2010-07-16 | 2012-01-19 | Mayo Foundation For Medical Education And Research | System and method for improved energy series of images using multi-energy ct |
US8189735B2 (en) * | 2010-07-22 | 2012-05-29 | General Electric Company | System and method for reconstruction of X-ray images |
US8781243B2 (en) * | 2011-01-07 | 2014-07-15 | Wisconsin Alumni Research Foundation | Method for constrained reconstruction of high signal-to-noise ratio images |
US9691168B2 (en) * | 2011-03-18 | 2017-06-27 | The Regents Of The University Of California | Image reconstruction using gradient projection for medical imaging applications |
US9268046B2 (en) * | 2011-07-12 | 2016-02-23 | Koninklijke Philips N.V. | Imaging system detector calibration |
EP2732437B1 (en) * | 2011-07-13 | 2017-09-06 | Rayvis GmbH | Transmission image reconstruction and imaging using poissonian detector data |
US20130016805A1 (en) * | 2011-07-15 | 2013-01-17 | Toshiba Medical Systems Corporation | Method and system for acquiring sparse channel data and for image processing utilizing iterative reconstruction algorithms |
AU2012326218B2 (en) * | 2011-10-17 | 2017-03-09 | Butterfly Network, Inc. | Transmissive imaging and related apparatus and methods |
US9121809B2 (en) * | 2011-11-18 | 2015-09-01 | Visuum, Llc | Multi-linear X-ray scanning systems and methods for X-ray scanning |
US8989348B2 (en) * | 2011-11-18 | 2015-03-24 | Visuum, Llc | Multi-linear X-ray scanning systems and methods for X-ray scanning |
US9013691B2 (en) * | 2012-01-29 | 2015-04-21 | Ramot At Tel-Aviv University | Snapshot spectral imaging based on digital cameras |
US9237874B2 (en) * | 2012-04-30 | 2016-01-19 | General Electric Company | Method and system for non-invasive imaging of a target region |
WO2013188957A1 (en) * | 2012-06-18 | 2013-12-27 | University Health Network | Method and system for compressed sensing image reconstruction |
US9460823B2 (en) * | 2012-09-10 | 2016-10-04 | Telesecurity Sciences, Inc. | Dynamic beam aperture control to reduce radiation dose using collimator |
US9036885B2 (en) * | 2012-10-28 | 2015-05-19 | Technion Research & Development Foundation Limited | Image reconstruction in computed tomography |
US9468409B2 (en) * | 2012-11-30 | 2016-10-18 | General Electric Company | Systems and methods for imaging dynamic processes |
US9384566B2 (en) * | 2013-03-14 | 2016-07-05 | Wisconsin Alumni Research Foundation | System and method for simulataneous image artifact reduction and tomographic reconstruction |
US9070218B2 (en) * | 2013-07-01 | 2015-06-30 | Xerox Corporation | Reconstructing an image of a scene captured using a compressed sensing device |
US9240071B2 (en) * | 2013-08-05 | 2016-01-19 | Siemens Aktiengesellschaft | Three-dimensional X-ray imaging |
US9508157B2 (en) * | 2014-12-12 | 2016-11-29 | Siemens Medical Solutions Usa, Inc. | Reconstruction of aneurysm wall motion |
-
2014
- 2014-03-24 EP EP14771044.6A patent/EP2976626A4/en not_active Withdrawn
- 2014-03-24 WO PCT/US2014/031609 patent/WO2014153561A1/en active Application Filing
-
2015
- 2015-09-22 US US14/861,312 patent/US10092253B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4819259A (en) * | 1985-09-20 | 1989-04-04 | Kabushiki Kaisha Toshiba | Rotary anode type x-ray apparatus |
US5812629A (en) * | 1997-04-30 | 1998-09-22 | Clauser; John F. | Ultrahigh resolution interferometric x-ray imaging |
US20060023842A1 (en) * | 2004-07-28 | 2006-02-02 | General Electric Company | Single-leaf X-ray collimator |
US20110006768A1 (en) * | 2009-07-09 | 2011-01-13 | Lei Ying | Systems and methods for accelerating the acquisition and reconstruction of magnetic resonance images |
US20120265050A1 (en) * | 2011-04-04 | 2012-10-18 | Ge Wang | Omni-Tomographic Imaging for Interior Reconstruction using Simultaneous Data Acquisition from Multiple Imaging Modalities |
Non-Patent Citations (1)
Title |
---|
See also references of EP2976626A4 * |
Also Published As
Publication number | Publication date |
---|---|
EP2976626A4 (en) | 2016-09-21 |
US10092253B2 (en) | 2018-10-09 |
EP2976626A1 (en) | 2016-01-27 |
US20160007938A1 (en) | 2016-01-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kennedy et al. | Super-resolution in PET imaging | |
EP2232444B1 (en) | Method for dynamic prior image constrained image reconstruction | |
EP2232446B1 (en) | Method for prior image constrained image reconstruction | |
JP6192542B2 (en) | Truncation correction for iterative cone beam CT reconstruction for SPECT / CT systems | |
US8600137B2 (en) | Iterative CT image reconstruction with a four-dimensional noise filter | |
JP5348855B2 (en) | Object image reconstruction method and apparatus for performing the method | |
US9336612B2 (en) | Contrast-dependent resolution image | |
WO2010127241A2 (en) | System and methods for fast implementation of equally-sloped tomography | |
EP1522045A2 (en) | Motion artifact correction of tomographical images | |
US20120308101A1 (en) | Method and system for reconstruction of tomographic images | |
US9953440B2 (en) | Method for tomographic reconstruction | |
US10092253B2 (en) | System, method, and computer accessible medium for modulating X-ray beam intensity | |
Hu et al. | Sinogram Restoration for Low‐Dosed X‐Ray Computed Tomography Using Fractional‐Order Perona‐Malik Diffusion | |
Chen et al. | Prior image constrained compressed sensing (PICCS) and applications in X-ray computed tomography | |
Muckley et al. | Image reconstruction for interrupted-beam x-ray CT on diagnostic clinical scanners | |
EP2715671A1 (en) | Fast computation of the sensitivity matrix in iterative algorithms | |
Choi et al. | 4D digital tomosynthesis image reconstruction using brute force-based adaptive total variation (BF-ATV) in a prototype LINAC system | |
Kalantari et al. | Attenuation correction in 4D‐PET using a single‐phase attenuation map and rigidity‐adaptive deformable registration | |
Muckley et al. | Reconstruction of reduced-dose sparse CT data acquired with an interrupted-beam prototype on a clinical scanner | |
Baek et al. | Initial results with a multisource inverse-geometry CT system | |
KR101379316B1 (en) | Multiple fan-beam helical computed tomography system and computed tomography method using thereof | |
Gilat et al. | Three-dimensional reconstruction algorithm for a reverse geometry volumetric CT system with a large-array scanned source | |
Wu et al. | A Practical Truncation Correction Method for Digital Breast Tomosynthesis | |
Lebedev | Motion Compensation in Cardiac X-Ray Computed Tomography | |
Beaudry | 4D cone-beam CT image reconstruction of Varian TrueBeam v1. 6 projection images for clinical quality assurance of stereotactic ablative radiotherapy to the lung |
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: 14771044 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2014771044 Country of ref document: EP |