A method of obtaining a computed tomography image of an object includes determining linear terms and non-linear beam hardening terms in a pair of line integral equations for dual-energy projection data from inserting average and difference from average attenuation terms, obtaining an initial solution of the line integral equation by setting the non-linear beam hardening terms to zero, and iteratively solving the line integral equations to obtain one line integral equations for each basis material. Attenuation by the first basis material corresponds to a photoelectric attenuation process, and attenuation by the second basis material corresponds to a Compton attenuation process. The line integral equations can be inverted by an inverse Radon procedure such as filtered backprojection to give images of each basis material. The images of each basis material can then be optionally combined to give monochromatic images, density and effective atomic number images, or photoelectric and Compton processes images.
