A kind of joint frequency-DOA spectrum estimation using virtual tapped-delay-line array is proposed, which generalizes the classical subspace algorithm to two-dimensional space-time domain. The proposed method has three advantages: firstly, increasing the freedom degrees of the corresponding array leading to an enhanced ability to identify much more emitters; secondly, almost free from estimation bias due to the divergence of signal wavelength from what the array size expresses; thirdly, estimating not only the DOAs but also the frequency domain characters of emitters. The obtained characteristic of emitters is more integral.