# US 10,499,869 B2 · Apparatus and method for material decomposition of spectrally resolved projection data using singles counts > Agents: this data is available as JSON from the patentgrep API and MCP server — see https://patentgrep.com/llms.txt - Number: US 10,499,869 B2 - Document ID: US10499869B2 - Status: Granted - Kind: B2 - Application number: 15/017,310 - Filed: 2016-02-05 - Published: 2017-08-10 - Granted: 2019-12-10 - Earliest claimed: 2016-02-05 - Assignees: Toshiba Medical Systems Corporation - Applicants: Kabushiki Kaisha Toshiba; Toshiba Medical Systems Corporation - Inventors: Adam Petschke; Gin-Chung Wang; Yu Zou - Web page: https://patentgrep.com/patent/US10499869B2 ## CPC - A61B 6/5211: {involving processing of medical diagnostic data} - A61B 6/032: {Transmission computed tomography [CT]} - A61B 6/4241: {using energy resolving detectors, e.g. photon counting} - A61B 6/482: {involving multiple energy imaging} - A61B 6/5205: {involving processing of raw data to produce diagnostic data} - G01T 1/171: {Compensation of dead-time counting losses} - A61B: DIAGNOSIS; SURGERY; IDENTIFICATION (usefulness limited to only animals A61D) - G01T: MEASUREMENT OF NUCLEAR OR X-RADIATION (radiation analysis of materials, mass spectrometry G01N23/00; tubes for determining the presence, intensity, density or energy of radiation or particles H01J47/00) ## Abstract A method and apparatus is provided to decompose spectral computed tomography (CT) projection data into material components using singles-counts and total-counts projection data. The singles-counts projection data more accurately solves the material decomposition problem, but can produce multiple results only one of which is correct. The total-counts projection data generates a unique result, but is less precise. The total-counts projection data is used to disambiguate the multiple results of the singles-counts projection data providing a unique results that is also precise. The unique and precise material decomposition can be achieved by limiting a search region for the singles-counts result to a neighborhood surrounding the total-counts result, choosing a singles-counts result that is closest to the total-counts result, choosing a singles-counts result that minimizes a total-counts cost function, or using a combined cost function that includes a singles-counts projection data and energy-integrated projection data. ## Claims ### Claim 1 An apparatus, comprising: processing circuitry configured to obtain first projection data having a plurality of energy components, the first projection data representing total counts including single X-ray detection events and multiple X-ray detection events, the total counts representing an intensity of X-ray radiation having been transmitted through an object and then detected at a plurality of energy-resolving detector elements; obtain second projection data having the plurality of energy components, wherein the second projection data represents singles counts of the single X-ray detection events, the singles counts representing the intensity of X-ray radiation having been transmitted through the object and then detected at the plurality of energy-resolving detector elements; and determine a material decomposition by estimating a first pair of projection lengths using a material decomposition of the first projection data, the first pair of projection lengths being a solution of a first material-decomposition method, and estimating a second pair of projection lengths using a material decomposition of the second projection data, the second pair of projection lengths being a solution of a second material-decomposition method. ### Claim 2 (depends on 1) The apparatus of claim 1, wherein the processing circuitry is further configured to determine the material decomposition with the first material-decomposition method being one of a cost-function method and a split step method, and the second material-decomposition method being one of a cost-function method and a split step method. ### Claim 3 (depends on 1) The apparatus of claim 1, wherein the processing circuitry is further configured to determine the material decomposition by determining a search region around the first pair of projection lengths using predefined criteria, and the second material-decomposition method is constrained to maintain the second pair of projection lengths within the search region. ### Claim 4 (depends on 3) The apparatus of claim 3, wherein the processing circuitry is further configured to determine the material decomposition by the second material-decomposition method starting a search to estimate the second pair of projection lengths at the first pair of projection lengths. ### Claim 5 (depends on 1) The apparatus of claim 1, wherein the processing circuitry is further configured to determine the material decomposition by estimating a plurality of pairs of projection lengths that are respectively solutions of the second material-decomposition method, the plurality of pairs of projection lengths including the second pair of projection lengths, determining distances between the first pair of projection lengths and respective pairs of the plurality of pairs of projection lengths, and selecting, as the material-decomposition, a pair of projection lengths of the plurality of pairs of projection lengths corresponding to a shortest distance of the distances. ### Claim 6 (depends on 5) The apparatus of claim 5, wherein the processing circuitry is further configured to determine the material decomposition by determining the distances using a Euclidean distance measure. ### Claim 7 An apparatus, comprising: processing circuitry configured to obtain first projection data having a plurality of energy components, the first projection data representing total counts including single X-ray detection events and multiple X-ray detection events, the total counts representing X-ray radiation detected at a plurality of energy-resolving detector elements; obtain second projection data having a plurality of energy components, wherein the second projection data represents singles counts of the single X-ray detection events, the singles counts representing the X-ray radiation detected at the plurality of energy-resolving detector elements; and determine a material decomposition by estimating a plurality of pairs of projection lengths using a material decomposition of the second projection data, each of the plurality of pairs of projection lengths locally optimizing a material-decomposition method; and generating a plurality of estimates of the first projection data using the plurality of pairs of projection lengths to determine an attenuation of the X-ray radiation, and selecting, as the material-decomposition, a pair of the plurality of pairs of projection lengths corresponding to an estimate of the plurality of estimates of the first projection data that is closest to the first projection data. ### Claim 8 (depends on 7) The apparatus of claim 7, wherein the processing circuitry is further configured to determine the material decomposition by determining which of the plurality of estimates of the first projection data is closest to the first projection data using a cost function that includes at least one distance measure, and the least one distance measure includes at least one of a least squares distance measure, a weighted least squares distance measure, and a Poisson likelihood distance measure. ### Claim 9 (depends on 7) The apparatus of claim 7, wherein the processing circuitry is further configured to determine the material decomposition by estimating the plurality of pairs of projection lengths using the material-decomposition method, which is one of a cost-function method and a split-step method. ### Claim 10 (depends on 7) The apparatus of claim 7, wherein the processing circuitry is further configured to estimate the plurality of pairs of projection lengths by optimizing a first pair of projection lengths of the plurality of pairs of projection lengths using optimization, generating an exclusion region in a region along a search path of the optimizing of the first pair of projection lengths, and optimizing a second pair of projection lengths of the plurality of pairs of projection lengths using optimization that is constrained to exclude the exclusion region, and expanding the exclusion region, when the optimizing of the second pair of projection lengths converges to the first pair of projection lengths, to include all optimization search paths converging to the first pair of projection lengths, and continuing the optimizing of the second pair of projection lengths constrained to exclude the exclusion region, until convergence to a second pair of projection lengths that is different from the first pair of projection lengths. ### Claim 11 An apparatus, comprising: processing circuitry configured to obtain first projection data, the first projection data representing integrated energies of X-ray radiation detected at a plurality of energy integrating detector elements; obtain second projection data having a plurality of energy components, wherein the second projection data represents singles counts, the singles counts representing detection events of single X-ray photons detected at a plurality of energy-resolving detector elements; and determine a pair of projection lengths using a combined material decomposition of the second projection data and the first projection data, wherein the pair of projection lengths optimize a combined cost function that uses a first distance measure between the first projection data and a model of the first projection data and uses a second distance measure between the second projection data and a model of the second projection data. ### Claim 12 (depends on 11) The apparatus of claim 11, wherein the processing circuitry is further configured to determine the pair of projection lengths by optimizing the pair of projection lengths to minimize the combined cost function, wherein the combined cost function includes a linear superposition of the first distance measure between the first projection data and a model of the first projection data measure and the second distance measure between the second projection data and a model of the second projection data, the first distance measure including at least one of a least squares distance measure, a weighted least squares distance measure, and a Poisson-likelihood distance measure, and the second distance measure including at least one of a least squares distance measure, a weighted least squares distance measure, and a Poisson-likelihood distance measure. ### Claim 13 (depends on 11) The apparatus of claim 11, wherein the processing circuitry is further configured to determine the pair of projection lengths by optimizing the pair of projection lengths using the model of the second projection data that accounts for absorption along X-ray trajectories using the pair of projection lengths, and optimizing the pair of projection lengths using the model of the first projection data that accounts for absorption along X-ray trajectories using the pair of projection lengths. ### Claim 14 (depends on 13) The apparatus of claim 13, wherein the processing circuitry is further configured to determine the pair of projection lengths by optimizing the pair of projection lengths using a global optimization method. ### Claim 15 An apparatus, comprising: an X-ray source radiating X-rays; a plurality of detector elements each configured to detect a plurality of energy components of the X-rays that are radiated from the X-ray source, generate first projection data representing total counts including a combination of single X-ray detection events and multiple X-ray detection events at each of the plurality of detector elements, and generate second projection data representing singles counts including the single X-ray detection events at each of the plurality of detector elements; and processing circuitry configured to determine a material decomposition by estimating a first pair of projection lengths using a material decomposition of the first projection data, the first pair of projection lengths being a solution of a first material-decomposition method, and estimating a second pair of projection lengths using a material decomposition of the second projection data, the second pair of projection lengths being a solution of a second material-decomposition method. ### Claim 16 (depends on 15) The apparatus of claim 15, wherein the processing circuitry is further configured to determine the material decomposition by determining a search region around the first pair of projection lengths using predefined criteria, and the second material-decomposition method is constrained to maintain the second pair of projection lengths within the search region. ### Claim 17 (depends on 16) The apparatus of claim 16, wherein the processing circuitry is further configured to determine the material decomposition by estimating a plurality of pairs of projection lengths that locally optimize the second material-decomposition method, the plurality of pairs of projection lengths including the second pair of projection lengths, determining distances between the first pair of projection lengths and respective pairs of the plurality of pairs of projection lengths, and selecting, as the material-decomposition, a pair of projection lengths of the plurality of pairs of projection lengths corresponding to a shortest distance of the distances. ### Claim 18 (depends on 15) The apparatus of claim 15, wherein the processing circuitry is further configured to determine the material decomposition by the first material-decomposition method being one of a cost-function method and a split step method, and the second material-decomposition method being one of a cost-function method and a split step method. ### Claim 19 A method, comprising: obtaining first projection data having a plurality of energy components, the first projection data representing total counts including single X-ray detection events and multiple X-ray detection events, the total counts representing an intensity of X-ray radiation having been transmitted through an object and then detected at a plurality of energy-resolving detector elements; obtaining second projection data having the plurality of energy components, wherein the second projection data represents singles counts of the single X-ray detection events, the singles counts representing the intensity of X-ray radiation having been transmitted through the object and then detected at the plurality of energy-resolving detector elements; and determining a material decomposition by estimating a first pair of projection lengths using a material decomposition of the first projection data, the first pair of projection lengths being a solution of a first material-decomposition method, and estimating a second pair of projection lengths using a material decomposition of the second projection data, the second pair of projection lengths being a solution of a second material-decomposition method. ### Claim 20 (depends on 19) A non-transitory computer readable storage medium including executable instruction, wherein the instructions, when executed by circuitry, cause the circuitry to perform the method according to claim 19.