A method includes generating landmarks, for different compositions of materials in scanned structure, wherein a landmark represents a location of a composition of materials, distorted by beam hardening, in a multi-energy Hounsfield Unit space and using the landmarks to perform an image domain material separation and quantification free of beam hardening artifacts.