A method of and a system for computing Z (effective atomic number) images from projection data are provided, wherein the projections are acquired using at least two x-ray spectra for a set of scanned objects, including a set of low energy projections and a set of high energy projections; the method comprises decomposing the low energy projections and high energy projections into photoelectric projections, reconstructing the photoelectric projections into photoelectric images, reconstructing one of the two sets of projections into CT images, and computing Z images from the CT images and the photoelectric images with parameters obtained from a calibration procedure.