Two-dimensional imaging of $V_s$ profile becomes more important in Korea because of the large horizontal variation of soil stiffness. To obtain a shear-wave velocity profile in geotechnical practice, various seismic nondestructive investigation methods are being frequently used. In this study, harmonic wavelet analysis of wave (HWAW) method is applied to the determination of $V_s$ profile to overcome some of weaknesses in the existing surface wave methods. HWAW method which is based on time-frequency analysis using harmonic wavelet transform has been developed to determine phase and group velocities of waves. Field testing of this method is relatively simple and fast because one experimental setup which consists of one pair of receivers is needed to determine $V_s$ profile of site. The proposed method uses the signal portion of the maximum local signal/noise ratio to evaluate the phase velocity to minimize the effects of noise, and uses single array inversion which considers receiver locations. Field tests were performed in 2 sites in order to evaluate accuracy of test method and estimate the applicability of 2-D imaging by HWAW method. Through field applications and comparison with other test results, the good accuracy and applicability of the proposed method were verified.