Receiver function (RF) is widely used for imaging crustal and uppermost mantle velocity structure based on horizontal layered assumption. However this traditional technique has limitations on resolving non-layered anomalies, such as the dipping Moho, the subduction zone, and plumes. Thus, we develop a innovative technique to invert RFs based on adjoint waveform tomography method for high-resolution seismic array imaging.