Conventional seismic imaging methods rely on the single-scattering Born approximation, requiring the removal of multiply scattered events from reflection data prior to imaging. Additionally, many methods use an acoustic approximation, representing the solid Earth as an acoustic (fluid) medium. We propose imaging methods for (solid) elastic media that use primaries and internal multiples, including their P S and SP conversions, thus obviating the need for internal multiple removal and improving handling of internal conversions. The methods rely on the elastic autofocusing method which creates multi-component elastodynamic Green's functions from virtual sources interior to the medium to receivers placed on the surface. They require only surface seismic reflection data and estimates of the direct waves from virtual sources interior to the medium, both of which are commonly available at the imaging step of seismic processing. We demonstrate our methods on a synthetic model with constant P and S velocities and vertical and horizontal density variations, by producing for the first time P P and SS images from elastic autofocusing which are compared to 1 reference seismic images based on conventional methods. Effects of multiples are greatly attenuated in the images, with fewer spurious reflectors than are observed when using Born imaging.2