Embodiments of the subject matter disclosed herein relate to medical imaging, and more particularly, to computed tomography (CT).
In computed tomography (CT) imaging systems, an x-ray source emits an x-ray beam toward a subject or object, such as a patient. After attenuation by the subject, the x-ray beam impinges upon a detector array. An intensity of the attenuated beam radiation received at the detector array depends on upon attenuation of the x-ray beam by the subject. Each detector element of the detector array produces a separate electrical signal which is transmitted to a data processing system for analysis and generation of a medical image. CT scans may be prone to beam hardening artifacts when imaging a patient.
In one example, a method includes obtaining projection data of a region of interest (ROI) of a subject, the projection data acquired at two different x-ray tube energy levels, generating a first monoenergetic image at a first energy level from the projection data, generating a second monoenergetic image at a second energy level from the projection data, subtracting the first monoenergetic image from the second monoenergetic image to obtain a subtracted image, and outputting the subtracted image for display and/or storage.
It should be understood that the brief description above is provided to introduce in simplified form a selection of concepts that are further described in the detailed description. It is not meant to identify key or essential features of the claimed subject matter, the scope of which is defined uniquely by the claims that follow the detailed description. Furthermore, the claimed subject matter is not limited to implementations that solve any disadvantages noted above or in any part of this disclosure.
The present disclosure will be better understood from reading the following description of non-limiting embodiments, with reference to the attached drawings, wherein below:
The following description relates to correcting metal artifacts present in images generated by a dual-energy spectral computed tomography (CT) system. In particular, systems and methods are provided for obtaining projection data of a region of interest (ROI) of a subject at two different energy levels, generating a first monoenergetic image at a first energy level and a second monoenergetic image at a second energy level from the projection data, the lower energy level and the higher energy level being different than the first energy level and the second energy level, and subtracting the two monoenergetic images to obtain a final image, which is output for display and/or storage.
When a metal component is positioned within a subject, such as a coil positioned in a brain aneurysm of the subject, images generated of the subject may include artifacts as a result of the metal, such as beam hardening artifacts, which may be referred to as metal artifacts. Metal artifacts may hinder accurate visualization of underlying and surrounding anatomical features of the subject proximate to the metal component. Hence, metal artifacts may result in a missed diagnosis, misdiagnosis, or hinder monitoring of a current diagnosis. In particular, metal artifacts due to the coil may hinder accurate visualization of underlying and surrounding vessels around the coil and within the coil itself, which in turn, hinders monitoring of the brain aneurysm.
Visualization of underlying and surrounding anatomical features of the subject proximate to the metal component may be increased by simulating a removal of the metal component and reducing the presence of metal artifacts within the generated images. Removal of the coil may be simulated by performing a dual-energy CT imaging protocol at two different x-ray tube energy levels following injection of contrast agent to obtain projection data and generating a first monoenergetic image at a first energy level from the projection data and a second monoenergetic image at a second, different energy level from the projection data, wherein the first energy level results in metal artifacts being retained in the first monoenergetic image and the second energy level results in reduced metal artifacts in the second monoenergetic image. By subtracting the first monoenergetic image from the second monoenergetic image, a subtracted image may be generated wherein the metal component is not present within the subtracted image. Since the metal component is not included in the subtracted image, the underlying and surrounding anatomical features of the subject may be accurately visualized. In this way, an existing diagnosis of the subject may be monitored and a frequency of missed diagnoses and misdiagnoses of the subject may be reduced.
An example medical imaging system that may be used to generate medical images is shown in
In certain embodiments, the CT imaging system 100 further includes an image processor unit 110 configured to reconstruct images of a target volume of the subject 112 using an iterative or analytic image reconstruction method. For example, the image processor unit 110 may use an analytic image reconstruction approach such as filtered back projection (FBP) to reconstruct images of a target volume of the patient. As another example, the image processor unit 110 may use an iterative image reconstruction approach such as advanced statistical iterative reconstruction (ASIR), conjugate gradient (CG), maximum likelihood expectation maximization (MLEM), model-based iterative reconstruction (MBIR), and so on to reconstruct images of a target volume of the subject 112. In some examples the image processor unit 110 may use an analytic image reconstruction approach such as FBP in addition to an iterative image reconstruction approach. In some embodiments, the image processor unit 110 may use a direct image reconstruction approach, such as using deep-learning trained neural networks.
In some CT imaging system configurations, the x-ray source 104 emits a cone-shaped beam which is collimated to lie within a plane of an X-Y-Z Cartesian coordinate system and generally referred to as an “imaging plane.” The radiation beam passes through an object being imaged, such as the patient or subject 112. The beam, after being attenuated by the object, impinges upon the detector array 108 comprising radiation detectors. The intensity of the attenuated radiation beam received at the detector array 108 is dependent upon the attenuation of the radiation beam by the object. Each detector element of the array produces a separate electrical signal that is a measurement of the beam attenuation of a ray path between the source and the detector element. The attenuation measurements from all the detector elements are acquired separately to produce a transmission profile.
In some CT imaging systems, the radiation source and the detector array are rotated with a gantry within the imaging plane and around the object to be imaged such that an angle at which the radiation beam intersects the object constantly changes. A group of radiation attenuation measurements, e.g., projection data, from the detector array at one gantry angle is referred to as a “view.” A “scan” of the object includes a set of views made at different gantry angles, or view angles, during one revolution of the radiation source and detector. It is contemplated that the benefits of the methods described herein accrue to medical imaging modalities other than CT, so as used herein the term “view” is not limited to the use as described above with respect to projection data from one gantry angle. The term “view” is used to mean one data acquisition whenever there are multiple data acquisitions from different angles, whether from a CT, a positron emission tomography (PET), a single-photon emission CT (SPECT) acquisition, and/or any other modality including modalities yet to be developed as well as combinations thereof in fused or hybrid embodiments.
The projection data is processed to reconstruct an image that corresponds to a two-dimensional slice taken through the object or, in some examples where the projection data includes multiple rotations or scans or two-dimensional (2D) arrays of detectors, a three-dimensional (3D) rendering of the object. One method for reconstructing an image from a set of projection data is referred to in the art as the filtered back projection technique. Transmission and emission tomography reconstruction techniques also include statistical iterative methods, such as maximum likelihood expectation maximization (MLEM) and ordered-subsets expectation-reconstruction techniques, as well as iterative reconstruction techniques. This process may convert the attenuation measurements from a scan into values called “CT numbers” or “Hounsfield units” (HU), which are used to control the brightness of a corresponding pixel on a display device.
To reduce the total scan time, a “helical” scan may be performed. To perform a “helical” scan, the patient is moved while the data for the prescribed number of slices are acquired. The position of the source with respect to the patient in such a system traces a helix. The helix mapped out by the source yields projection data from which images in each prescribed slice may be reconstructed.
As used herein, the phrase “reconstructing an image” is not intended to exclude embodiments of the present invention in which data representing an image are generated but a viewable image is not. Therefore, as used herein, the term “image” broadly refers to both viewable images and data representing a viewable image. However, many embodiments generate (or are configured to generate) at least one viewable image.
In certain embodiments, the imaging system 200 is configured to traverse different angular positions around the subject 112 for acquiring desired projection data. Accordingly, the gantry 102 and the components mounted thereon may be configured to rotate about a center of rotation 206 for acquiring the projection data, for example, at different energy levels. Alternatively, in embodiments where a projection angle relative to the subject 112 varies as a function of time, the mounted components may be configured to move along a general curve rather than along a segment of a circle.
As the x-ray source 104 and the detector array 108 rotate, the detector array 108 collects data of the attenuated x-ray beams. The data collected by the detector array 108 undergoes pre-processing and calibration to condition the data to represent the line integrals of the attenuation coefficients of the scanned subject 112. The processed data are commonly called projections.
In some examples, the individual detectors or detector elements 202 of the detector array 108 may include photon-counting detectors which register the interactions of individual photons into one or more energy bins. It should be appreciated that the methods described herein may also be implemented with energy-integrating detectors.
The acquired sets of projection data may be used for basis material decomposition (BMD). During BMD, the measured projections are converted to a set of material-density projections. The material-density projections may be reconstructed to form a pair or a set of material-density map or image of each respective basis material, such as bone, soft tissue, and/or contrast agent maps. The density maps or images may be, in turn, associated to form a volume rendering of the basis material, for example, bone, soft tissue, and/or contrast agent, in the imaged volume.
Once reconstructed, the basis material image produced by the imaging system 200 reveals internal features of the subject 112, expressed in the densities of two basis materials. The density image may be displayed to show these features. In traditional approaches to diagnosis of medical conditions, such as disease states, and more generally of medical events, a radiologist or physician may consider a hard copy or display of the density image to discern characteristic features of interest. Such features might include lesions, sizes and shapes of particular anatomies or organs, and other features that would be discernable in the image based upon the skill and knowledge of the individual practitioner.
In one embodiment, the imaging system 200 includes a control mechanism 208 to control movement of the components such as rotation of the gantry 102 and the operation of the x-ray source 104. In certain embodiments, the control mechanism 208 further includes an x-ray controller 210 configured to provide power and timing signals to the x-ray source 104. Additionally, the control mechanism 208 includes a gantry motor controller 212 configured to control a rotational speed and/or position of the gantry 102 based on imaging requirements.
In certain embodiments, the control mechanism 208 further includes a data acquisition system (DAS) 214 configured to sample analog data received from the detector elements 202 and convert the analog data to digital signals for subsequent processing. The DAS 214 may be further configured to selectively aggregate analog data from a subset of the detector elements 202. The data sampled and digitized by the DAS 214 is transmitted to a computer or computing device 216. In one example, the computing device 216 stores the data in a mass storage device or storage device 218. The storage device 218, for example, may include a hard disk drive, a floppy disk drive, a compact disk-read/write (CD-R/W) drive, a Digital Versatile Disc (DVD) drive, a flash drive, and/or a solid-state storage drive.
Additionally, the computing device 216 provides commands and parameters to one or more of the DAS 214, the x-ray controller 210, and the gantry motor controller 212 for controlling system operations such as data acquisition and/or processing. In certain embodiments, the computing device 216 controls system operations based on operator input. The computing device 216 receives the operator input, for example, including commands and/or scanning parameters via an operator console 220 operatively coupled to the computing device 216. The operator console 220 may include a user interface (not shown), which may include one or more of a keyboard, a touchscreen, a mouse, a trackpad, and the like to allow the operator to specify the commands and/or scanning parameters.
Although
In one embodiment, the imaging system 200 either includes, or is coupled to, a picture archiving and communications system (PACS) 224. In an exemplary implementation, the PACS 224 is further coupled to a remote system such as a radiology department information system, hospital information system, and/or to an internal or external network (not shown) to allow operators at different locations to supply commands and parameters and/or gain access to the image data.
The computing device 216 uses the operator-supplied and/or system-defined commands and parameters to operate a table motor controller 226, which, in turn, may control a table 114 (see
As previously noted, the DAS 214 samples and digitizes the projection data acquired by the detector elements 202. Subsequently, an image reconstructor 230 uses the sampled and digitized x-ray data to perform high-speed reconstruction. Although
In one embodiment, the image reconstructor 230 stores the images reconstructed in the storage device 218. Alternatively, the image reconstructor 230 may transmit the reconstructed images to the computing device 216 for generating useful patient information for diagnosis and evaluation. In certain embodiments, the computing device 216 may transmit the reconstructed images and/or the patient information to a display or display device 232 communicatively coupled to the computing device 216 and/or the image reconstructor 230. In some embodiments, the reconstructed images may be transmitted from the computing device 216 or the image reconstructor 230 to the storage device 218 for short-term or long-term storage.
As shown in
Referring to
As set forth above, CT imaging system 100 operates as a dual-energy system, thereby providing for correction for the beam hardening and/or other artifacts caused by the metal component. CT imaging system 100 utilizes dual-energy scanning to obtain diagnostic CT images that enhance contrast separation within the image by acquiring data of the same locations at two different tube peak voltage levels (kVp). Using the images obtained during the dual energy CT scan, monochromatic images (also referred to as monoenergetic images) are generated. That is, a monochromatic image can be created at a desired energy level with each image being at a specific x-ray energy (keV). While embodiments of the disclosure described below describe generation of monochromatic images at two levels, it is recognized that monochromatic images can be created at a greater number of energy levels, such as three, four, or five energy levels, for example. Thus, embodiments of the disclosure are not limited to the acquiring and generating monochromatic images at two energy levels (e.g., a pair of “high” and “low” energy images), and discussion of such is not to be construed as such a limitation.
By utilizing dual energy CT, the effect of beam hardening artifact of the metal component can be minimized by looking at the same material at different virtual monochromatic energy levels. That is, high monochromatic energy levels reduce the intensity and corresponding beam hardening artifacts caused by the metal component. By obtaining an image at a high monochromatic energy level, the beam hardening artifacts of the metal component may be reduced when subtracted from an image obtained at a low monochromatic energy level, thus allowing for accurate visualization of the anatomical features surrounding and within the metal component.
An uncorrected image 502 may be generated based on the projection data obtained during the dual-energy CT scan protocol. In some examples, the dual-energy CT scan protocol may be an angiography protocol that includes scanning a head of a subject after administration of contrast agent, wherein the subject includes a coil (e.g., a metal component) that was inserted in a vessel of the subject to treat a condition (e.g., an aneurysm). Hence, the uncorrected image 502 includes beam hardening artifacts due to the coil being scanned, which may render the uncorrected image 502 non-diagnostic and/or prevent the image from being used to monitor a current diagnosis of the subject. In some embodiments, the uncorrected image 502 may be a maximum intensity projection (MIP) at 70 keV.
The dual-energy CT scan protocol may acquire projection data at two different peak energy levels. More specifically, projection data is acquired by nearly simultaneously interleaving the two peak energy levels. However, each energy level includes a spectrum of energies centered at the respective peak. A first projection dataset at a lower energy level and a second dataset at a higher energy level may be extracted from the projection data to form a lower energy sinogram and a high energy sinogram. Material-density (MD) projection maps of basis materials are generated based on the lower energy sinogram and the higher energy sinogram. The MD projection maps may be reconstructed to generate basis material images, such as iodine and water images.
After noise filtering is applied to the iodine and water images, a first monoenergetic image 504 and a second monoenergetic image 506 may be generated based on linear combinations of the reconstructed basis material images at a single energy level between 40 keV and 140 keV. For example, the first monoenergetic image 504 may be generated from a first linear combination of the basis material images and the second monoenergetic image 506 may be generated from a second linear combination of the basis material images. The first monoenergetic image 504 is at a different energy level (e.g., 40 keV) than the second monoenergetic image 506 (e.g., 140 keV). The first monoenergetic image 504 may include significant beam hardening artifacts whereas the second monoenergetic image 506 may include reduced beam hardening artifacts when compared to the first monoenergetic image 504.
The subtraction of the two monoenergetic images may be performed in response to user input received at a user interface 508 of a display device of the computing device. In particular, the user interface 508 may include a combination panel having an operation area 510, an operation saving area 512, and a mode selection area 514. The mode selection area 514 may include a plurality of widgets (e.g., buttons). When user input is received at the plurality of widgets, the user interface performs a pre-determined action. User input may be received at the plurality of widgets when a user input device receives user input (e.g., via clicking a mouse or touchscreen). As one example, a button 522 may be configured to receive user input. In response to receiving user input at the button 522, a subtraction mode of the user interface 508 may be selected. After selecting the subtraction mode, a first image, such as the first monoenergetic image 504, may be selected in the operation area 510 when a first widget 516 receives user input. Additionally, a second image, such as the second monoenergetic image 506, may be selected in the operation area 510 when user input is received at a second widget 518. The second monoenergetic image 506 may be subtracted from the first monoenergetic image 504 when user input is received at a third widget 520 to generate a final image 524 of the subject. In other examples, the first and second images may be selected via input to the first widget 516 and the second widget 518 and the button 522 and/or third widget 520 may be selected to generate the final image. The final image 524 may be saved in memory when user input is received in an operation saving area 512. The final image 524 may be an image of the head of the subject with reduced metal artifacts relative to the uncorrected image 502. In some embodiments, an MIP image may be generated from the final image 524 with a multiplanar reconstruction application integrated in the image processing system.
At 602, the method 600 includes obtaining projection data of a region of interest (ROI) of a subject at two different x-ray tube energy levels. In some embodiments, the ROI may include a metal component inserted within or placed next to an anatomical part. Projection data may be obtained with a dual-energy computed tomography (CT) scan protocol at two different x-ray tube energy levels, which changes the peak and spectrum of energy of the incident photons comprising the emitted x-ray beams. The projection data may be obtained with the CT imaging system of
The scan protocol may include exposing the subject to radiation at a lower x-ray tube energy level and a higher energy x-ray tube energy level. In particular, exposure may be performed by fast switching (or fast kV switching) of the peak tube voltage of the x-ray source of the CT imaging system. In some embodiments, x-ray tube may be controlled to switch between energy levels of 80 kVp and 140 kVp. However, it is recognized that other suitable energy levels may also be selected based on which anatomical parts and features are being scanned.
The x-ray source energy levels may be interleaved as a function of the rotation angle and thus one gantry rotation around the subject is performed to acquire the projection data at the lower x-ray tube energy level and the higher x-ray tube energy level. In this way, projection data may be acquired without having to acquire projection data for each x-ray tube energy level separately (e.g., via the use of two x-ray tubes or two rotations of the gantry). As one example, during one rotation of the gantry, the energy level in the x-ray tube may switch between the lower x-ray tube energy level and the higher x-ray tube energy according to a timing scheme which simulates the subject being exposed to the lower x-ray tube energy level and the higher x-ray tube energy level at around the same time. In other words, the subject is exposed to the lower x-ray tube energy level for a short duration of time followed by the subject being exposed to the higher x-ray tube energy level for another short duration of time, which is repeated until the scan is completed.
When the x-ray tube is at the lower x-ray tube energy level, the subject is exposed to x-ray beams at the lower x-ray tube energy level and the x-ray beams are attenuated by the subject. The attenuated x-ray beams impinge on the detector elements of the detector, such as the detector 400 of
After exposing the subject to radiation at the lower x-ray tube energy level for the short duration of time, the x-ray source may switch from the lower x-ray tube energy to the higher x-ray tube energy level for the short duration of time. When the x-ray tube is at the higher x-ray tube energy level, the subject is exposed to x-ray beams at the higher x-ray tube energy level and the x-ray beams are attenuated by the subject. As described, the attenuated x-ray beams impinging on the detector elements of the detector may be converted to an analog signal by the detector element. The analog signals are then transmitted to the DAS where the analog signals are converted to an electrical signal, which may be transmitted to the computing device for processing and image reconstruction. After exposing the subject to radiation at the higher x-ray tube energy level, the x-ray source switches back to the lower x-ray tube energy level.
The exposure sequence continues until the gantry finishes a single rotation. In this way, projection data may be obtained when the subject is exposed to radiation at the higher x-ray tube energy level and the lower x-ray tuber energy level. Even though the subject is exposed to radiation at two different x-ray tube peak energy levels, projection data corresponding to various energy levels that are different from the two different x-ray tube peak energy levels is collected since the x-ray tube is polychromatic. It may be understood that the entirety of the projection data may be obtained after injection of a contrast agent.
At 604, the method 600 includes generating a first monoenergetic image at a first energy level from the projection data. As described herein, projection data for a spectrum of energy levels is obtained with the dual-energy CT scan protocol, forming an overall projection dataset. Projection data corresponding to a lower energy level, and projection data corresponding to a higher energy level may be extracted from the overall projection dataset to form a first projection dataset and second projection dataset, respectively. The first projection dataset and the second projection dataset may be sinograms.
Basis material decomposition (BMD) may be performed using the first projection dataset and the second projection dataset. During BMD, the measured projections are converted to a set of material-density projections. The material-density projections may be reconstructed to form a pair or a set of material-density map or images of each respective basis material, such as water and iodine. The density maps or images may be, in turn, associated to form a volume rendering of the basis material, for example, water and iodine. A linear combination of the material-density images may be used to generate monoenergetic images at a specific energy level with an image reconstructor, such as image reconstructor 230 of
In some embodiments, the first energy level may be 40 keV. However, in other embodiments, the first energy level may differ based on which anatomical parts and/or features are being scanned during the dual-energy CT scan protocol. Regardless of a value of the first energy level, the first energy level may be lower than a second energy level. Additionally, the first energy level is selected such that the first monoenergetic image generated includes image artifacts that significantly reduce image quality, such as beam hardening artifacts and/or metal artifacts. As such, it may be difficult to visualize anatomical parts or features positioned at or near the ROI in the first monoenergetic image.
At 606, the method 600 includes generating a second monoenergetic image at the second energy level from the projection data. The second monoenergetic image may be generated as explained above, using a second linear combination of the material-density images. In some embodiments, the second energy level may be 140 keV. Similar to the first energy level, the second energy level may differ based on which anatomical parts and/or features are being scanned and imaged during the dual energy CT scan protocol. Further, regardless of the value of the second energy level, the second energy level may be higher than the first energy level. The second energy level is selected such that the second monoenergetic image at the second energy level includes reduced image artifacts, such as beam hardening artifacts or metal artifacts. As such, it may be easier to visualize anatomical parts or features positioned at or near the ROI in the second monoenergetic image when compared with the first monoenergetic image.
At 608, the method 600 includes subtracting the first monoenergetic image from the second monoenergetic image to obtain a subtracted image. The subtraction may be performed in the image domain (e.g., the images themselves are subtracted on a pixel-wise basis). In some examples, the first monoenergetic image may be subtracted from the second monoenergetic image according to a method described below with regards to
At 610, the method 600 includes outputting the subtracted image for display and/or storage. The subtracted image may be displayed on a user interface of a display device. In some embodiments, the subtracted image may be displayed on the user interface alongside an uncorrected image generated from the overall projection dataset. The uncorrected image may be a monoenergetic image at a specific image level. In one example, the uncorrected image may be at 70 keV. In other embodiments, the subtracted image may be displayed alongside a metal artifact reduction (MAR)-corrected image generated from the overall projection dataset and corrected using a conventional metal artifact reduction method. The conventional MAR method may include using information from spectral energy data to segment regions of the detector, such as detector 400 of
An embodiment of the user interface is illustrated in
At 702, the method 700 includes performing an injection of a contrast agent. Prior to performing the CTA protocol to obtain projection data of a subject, the contrast agent may be administered to the subject. In particular, the injection of the contrast agent may be performed with intravenous injection. In other embodiments, the contrast agent may be administered orally and/or administered according to other methods. In some embodiments, the contrast agent may be iodine-based or barium-based. As described previously, the subject may include a metal component, such as a coil, positioned within a head of the subject, for example within an aneurysm of a brain of the subject. By imaging the head of the subject, the aneurysm may be monitored by visualizing the area and/or anatomical features at and near the coil.
At 704, the method 700 includes, at each acquisition of the CTA protocol, controlling an x-ray source of the CT imaging system to emit x-ray radiation at two peak energy levels to obtain an overall projection dataset of the head of the subject. The CTA protocol may include one or more acquisitions that occur at specified timepoints relative to the contrast uptake of the subject, such as a first acquisition at the arterial peak of contrast uptake and optionally one or more additional acquisitions at the venous peak and/or venous return to baseline of the contrast agent. In particular, instructions configured, stored, and executed in memory of the computing device by the one or more processors may cause the processor(s) to, for each acquisition, adjust a peak energy level of the x-ray source of the CT imaging system to a first peak energy for a first pre-determined amount of time. Following the first pre-determined amount of time, the processor(s) may adjust the peak energy level of the x-ray source to a second peak energy for a second pre-determined amount of time.
The peak energy level of the x-ray source may be adjusted between the first peak energy level for the first pre-determined amount of time and the second peak energy level for the second pre-determined amount of time for the duration of each rotation of the gantry. The first pre-determined amount of time and the second pre-determined amount of time may be selected such that the head of the subject is nearly simultaneously exposed to both the first peak energy level and the second peak energy level (e.g., the x-ray source may switch energy levels every millisecond or less). In this way, the head of a subject may be exposed to radiation of the x-ray source at both of the first peak energy level and the second peak energy level at around the same time.
The x-ray source of the CT imaging system may emit radiation that is polychromatic, meaning that although the x-ray source emits radiation at the first peak energy level and the second peak energy level, additional energy levels around the peaks are also emitted, and thus the projection data obtained includes information for various energy levels. The various energy levels are different from the first peak energy level and the second peak energy level. The projection data may be included as part of an overall projection dataset that encompasses projection data of a spectrum of energy levels. It may be understood that the entirety of the projection data is obtained after injection of the contrast agent.
At 706, the method 700 includes generating a first monoenergetic image at a first energy level from the overall projection dataset. The first monoenergetic image may be a contrast image generated from the overall projection dataset according to the method described with respect to
At 708, the method 700 includes generating a second monoenergetic image at a second energy level from the overall projection dataset. The second monoenergetic image may be a contrast image generated from the overall projection dataset according to the method described with respect to
As described previously, the first and second energy levels are selected such that the second monoenergetic image includes reduced metal artifacts relative to the first monoenergetic image. In particular, image artifacts including beam hardening artifacts and/or metal artifacts associated with the high density coil may be present in the second monoenergetic image. However, when compared with the first monoenergetic image, the second monoenergetic image includes fewer metal artifacts. As such, it may be easier to visualize anatomical features of the brain positioned at the coil or near the coil in the second monoenergetic image compared to the first monoenergetic image. However, the extent of visualization of the anatomical features of the brain positioned at the coil or near the coil in the second monoenergetic image may not be suitable for monitoring the current diagnosis of the subject (e.g., the brain aneurysm).
At 710, the method 700 includes subtracting the first monoenergetic image from the second monoenergetic image to form a subtracted image. The subtraction may be performed in an image domain. In some examples, the first monoenergetic image may be subtracted from the second monoenergetic image according to the method described below with regards to
At 712, the method 700 includes outputting the subtracted image, optionally alongside an uncorrected image, to an interface of a display device and/or for storage in memory. Similar to
In one example, the coil in the brain may be the metal component inserted within or placed next to the anatomical brain that induces metal artifacts according to a duel-energy subtraction method described herein with respect to
At 802, the method 800 includes generating a first monoenergetic image in response to a first user input and a second monoenergetic image in response to a second user input. The first monoenergetic image and the second monoenergetic image may be generated according to the methods described above with respect to
In some examples, the user interface may include a combination panel, wherein the user may select images to combine (e.g., add, subtract, etc.). After generating the first and second monoenergetic images, the user may select a user interface element to launch the combination panel (such as the combination panel/user interface 508 of
At 804, the method 800 includes receiving a user input to the user interface selecting the first monoenergetic image at the first energy level. The user interface of the image processing system (and more specifically, the combination panel) may include an operation area that includes a first interface widget. In response to user selection of the first interface widget, a first image may be selected, such as the first monoenergetic image. As explained previously, the first image may be obtained at a lower x-ray tube energy level and an extent of beam hardening artifacts and/or metal artifacts in the first image may be significant. User input may be received at the first interface widget via a user input device (e.g., via a touchscreen or mouse).
At 806, the method 800 includes receiving a user input at the user interface selecting a second monoenergetic image at the second energy level. Referring to the operation area described above, a second interface widget displayed in the operation area may receive user input. In response to reception of user input at the second interface widget, a second image may be selected, such as the second monoenergetic image, with reduced beam hardening artifacts and/or metal artifacts. As explained previously, the second image may be obtained at a higher x-ray tube energy level and an extent of beam hardening artifacts and/or metal artifacts in the second image may be reduced when compared with the first image. User input may be received at the second interface widget via a user input device (e.g., via a touchscreen or mouse).
At 808, the method 800 include receiving a user input selecting a third interface widget of the user interface (e.g., the display panel) to subtract the first monoenergetic image from the second monoenergetic image. The third interface widget may be displayed in the operation area and may receive user input to trigger subtraction of the first image and the second image. As one example, the third interface widget may be a button. User input may be received at the third interface widget via a user input device (e.g., via a touchscreen or mouse). Accordingly, the first monoenergetic image may be subtracted from the second monoenergetic image, and the subtracted image may be displayed on the user interface. The method 800 then ends.
Thus, method 800 described herein provides for generation of a dual-energy subtracted image in response to user input, wherein the user input specifies the generation of each of the first monoenergetic image at the first energy level, the second monoenergetic image at the second energy level, and the subtracted image. However, it may be understood that the subtracted image may be generated in various ways. For example, the first monoenergetic image, the second monoenergetic image, and the subtracted image may be generated automatically as part of a scan protocol, such that the subtracted image is generated and displayed without explicit user input. In another example, the subtracted image may be generated in response to user input requesting generation of the subtracted image, and the first monoenergetic image, the second monoenergetic image, and the subtracted image may be generated in response to the user request to generate the subtracted image. In this way, the first and second monoenergetic images may be generated without explicit user input (e.g., the computing device may know to generate the monoenergetic images at the specific energy levels in order to generate the subtracted image).
The uncorrected image 1106 is generated from projection data obtained with a dual-energy computed tomography angiography (CTA) scan protocol. The projection data includes data for a spectrum of energy levels. The uncorrected image 1106 is an image of a brain of a subject with significant beam hardening artifacts and metal artifacts. A coil is positioned within an aneurysm of the brain of the subject as shown in the uncorrected image 1106.
The corrected image 1108 is generated from a subtraction of a first monoenergetic image at a first energy level and a second monoenergetic image at a second energy level according to the embodiments described herein. The corrected image 1108 is the image of the brain of the subject without significant beam hardening artifacts and metal artifacts. A correction of the uncorrected image 1106 results in the generation of the corrected image 1108. Therefore, the coil is no longer visible in the corrected image 1108, which enables a user to visually identify and monitor vessels positioned at or surrounding the coil. In this way, the user may monitor a condition of an aneurysm in the brain.
The user interface 1100 may additionally include an image output area 1110. The image output area may include information regarding both of the uncorrected image 1106 and the corrected image 1108. In some embodiments, the information may be measurements. For example, the image output area 1110 may include attenuation coefficients for the uncorrected image 1106 and the corrected image.
It may be understood that the user interface 1100 may include more or fewer components and other configurations than that shown in
The technical effect of subtracting two monoenergetic images generated from projection data at two different energy levels is that beam hardening artifacts and metal artifacts due to a metal component may be reduced and/or are not visible in the subtracted image. In this way, anatomical features positioned at the metal component or near the metal component may be accurately visualized, enabling a medical professional to make an accurate diagnosis or monitor a current diagnosis.
The disclosure also provides support for a method for computed tomography (CT) imaging, comprising: obtaining projection data of a region of interest (ROI) of a subject, the projection data acquired at two different x-ray tube energy levels, generating a first monoenergetic image at a first energy level from the projection data, generating a second monoenergetic image at a second energy level from the projection data, subtracting the first monoenergetic image from the second monoenergetic image to obtain a subtracted image, and outputting the subtracted image for display and/or storage. In a first example of the method, the projection data is acquired at 80 kVp and 140 kVp. In a second example of the method, optionally including the first example, an entirety of the projection data is acquired following an injection of contrast agent to the subject. In a third example of the method, optionally including one or both of the first and second examples, the ROI includes a metal component. In a fourth example of the method, optionally including one or more or each of the first through third examples, the metal component is a coil. In a fifth example of the method, optionally including one or more or each of the first through fourth examples, the first energy level is 40 keV. In a sixth example of the method, optionally including one or more or each of the first through fifth examples, the second energy level is 140 keV.
The disclosure also provides support for an image processing system, comprising: one or more processors, and memory storing instructions executable by the one or more processors to: obtain projection data of a region of interest (ROI) of a subject, the projection data acquired at two different x-ray tube energy levels, generate a first monoenergetic image at a first energy level from the projection data, generate a second monoenergetic image at a second energy level from the projection data, subtract the first monoenergetic image from the second monoenergetic image to obtain a subtracted image, and output the subtracted image for display and/or storage. In a first example of the system, subtracting the first monoenergetic image from the second monoenergetic image to obtain the subtracted image comprises: receiving user input selecting the first monoenergetic image at a user interface displayed on a display device coupled to the image processing system, receiving user input selecting the second monoenergetic image at the user interface, and receiving user input selecting an interface widget at the user interface, wherein instructions are executable to subtract the first monoenergetic image from the second monoenergetic image in response to receiving the user input selecting the interface widget. In a second example of the system, optionally including the first example, the subtracted image is output on the user interface of the display device alongside an uncorrected image generated from the projection data.
The disclosure also provides support for a method for a computed tomography (CT) imaging system, comprising: upon injection of a contrast agent to a subject, executing a dual-energy CT angiography scan protocol to acquire projection data of a head of the subject, the head of the subject including a coil positioned therein, generating a first monoenergetic image at a first energy level from the projection data, generating a second monoenergetic image at a second energy level from the projection data, subtracting the first monoenergetic image from the second monoenergetic image to obtain a subtracted image, and outputting the subtracted image for display and/or storage. In a first example of the method, executing the dual-energy CT angiography scan protocol includes controlling an x-ray source of the CT imaging system to emit x-ray radiation at a first peak energy level and at a second peak energy level. In a second example of the method, optionally including the first example, the first peak energy level is 80 kVp and the second peak energy level is 140 kVp. In a third example of the method, optionally including one or both of the first and second examples, the first energy level is 40 keV and the second energy level is 140 keV.
In a fourth example of the method, optionally including one or more or each of the first through third examples, a first projection dataset and a second projection dataset are extracted from the projection data to generate a first material-density image and a second material-density image. In a fifth example of the method, optionally including one or more or each of the first through fourth examples, first monoenergetic image at the first energy level is generated based on a first linear combination of the first material-density image and the second material-density image. In a sixth example of the method, optionally including one or more or each of the first through fifth examples, the second monoenergetic image at the second energy level is generated based on a second linear combination of the first material-density image and the second material-density image. In a seventh example of the method, optionally including one or more or each of the first through sixth examples, the first peak energy level is different than the first energy level. In an eighth example of the method, optionally including one or more or each of the first through seventh examples, metal artifact reduction is not applied to the subtracted image. In a ninth example of the method, optionally including one or more or each of the first through eighth examples, outputting the subtracted image for display comprises outputting the subtracted image for display alongside an uncorrected image generated from the projection data.
As used herein, an element or step recited in the singular and preceded with the word “a” or “an” should be understood as not excluding plural of said elements or steps, unless such exclusion is explicitly stated. Furthermore, references to “one embodiment” of the present invention are not intended to be interpreted as excluding the existence of additional embodiments that also incorporate the recited features. Moreover, unless explicitly stated to the contrary, embodiments “comprising,” “including,” or “having” an element or a plurality of elements having a particular property may include additional such elements not having that property. The terms “including” and “in which” are used as the plain-language equivalents of the respective terms “comprising” and “wherein.” Moreover, the terms “first,” “second,” and “third,” etc. are used merely as labels, and are not intended to impose numerical requirements or a particular positional order on their objects.
This written description uses examples to disclose the invention, including the best mode, and also to enable a person of ordinary skill in the relevant art to practice the invention, including making and using any devices or systems and performing any incorporated methods. The patentable scope of the invention is defined by the claims, and may include other examples that occur to those of ordinary skill in the art. Such other examples are intended to be within the scope of the claims if they have structural elements that do not differ from the literal language of the claims, or if they include equivalent structural elements with insubstantial differences from the literal languages of the claims.