A three step model for high harmonic generation from impurities in solids is developed. The process is found to be similar to high harmonic generation in atomic and molecular gases with the main difference coming from the non-parabolic nature of the bands. This opens a new avenue for strong field atomic and molecular physics in the condensed matter phase. As a first application, our conceptual study demonstrates the feasibility of tomographic measurement of impurity orbitals.