Doppler ultrasound systems, which are widely used in cardiovascular disorders detection, have quadrature format outputs. Various types of algorithms were described in literature to process quadrature Doppler signals (QDS), such as phasing filter technique (PFT), fast Fourier transform method, frequency domain Hilbert transform method and complex continuous wavelet transform. However for the discrete wavelet transform (DWT) case, which becomes a common method for processing QDSs, there was not a direct method to recover flow direction from quadrature signals. Traditionally, to process QDSs with DWT, firstly directional signals have to be extracted and later two DWTs must be applied. Although there exists a complex DWT algorithm called dual tree complex discrete wavelet transform (DTCWT), it does not provide directional signal decoding during analysis because of the unwanted energy leaks into its negative frequency bands. Modified DTCWT, which is a combination of PFT and DTCWT, has the capability of extracting directional information while decomposing QDSs into different frequency bands, but it uses an additional Hilbert transform filter and it increases the computational complexity of whole transform. Discrete wavelet packet transform (DWPT), which is a generalization of the ordinary DWT allowing subband analysis without the constraint of dyadic decomposition, can perform an adaptive decomposition of the frequency axis. In this study, a novel complex DWPT, which maps directional information while processing QDSs, is proposed. The success of proposed method will be measured by using simulated quadrature signals.