Recent work [J.S. Lundeen et al. Nature, 474, 188 (2011)] directly measured the wavefunction by weakly measuring a variable followed by a normal (i.e. 'strong') measurement of the complementary variable. We generalize this method to mixed states by considering the weak measurement of various products of these observables, thereby providing the density matrix an operational definition in terms of a procedure for its direct measurement. The method only requires measurements in two bases and can be performed 'in situ', determining the quantum state without destroying it.PACS numbers: 03.65. Ta, 03.65.Wj, 42.50.Dv, The wavefunction Ψ is at the heart of quantum mechanics, yet its nature has been debated since its inception. It is typically relegated to being a calculational device for predicting measurement outcomes. Recently, Lundeen et al. proposed a simple and general operational definition of the wavefunction based on a method for its direct measurement: "it is the average result of a weak measurement of a variable followed by a strong measurement of the complementary variable [1,2]." By 'direct' it is meant that a value proportional to the wavefunction appears straight on the measurement apparatus itself without further complicated calculations or fitting. The 'wavefunction' referred to here is a special case of a general quantum state, known as a 'pure state.' The general case is represented by the density operator ρ, which can describe both pure and 'mixed' states. The latter incorporates both the effects of classical randomness (e.g., noise) and entanglement with other systems (e.g., decoherence). The density operator plays an important role in quantum statistics, quantum information, and the study of decoherence. Because of its generality and because it follows naturally from classical concepts of probability and measures, some consider ρ, rather than Ψ, to be the fundamental quantum state description. In this letter, we propose two methods to directly measure general quantum states, one of which directly gives the matrix elements of ρ.The standard method for experimentally determining the density operator is Quantum State Tomography [3]. In it, one makes a diverse set of measurements on an ensemble of identical systems and then determines the quantum state that is most compatible with the measurement results. An alternative is our direct measurement method, which may have advantages over tomography, such as simplicity, versatility, and directness. A quantitative comparison of measures such as the signal to noise ratio, resolution, and fidelity, has not been undertaken but some limitations of the direct method have been identified in [4]. As compared to tomography, which works with mixed states, the most significant limitation of the direct measurement of the wavefunction is that it has only been shown to work with pure states.Previous works have developed direct methods to measure quasi-probability distributions, such as the Wigner function [5], Husimi Q-function [6], and the GlauberSudarshan P-function [7]...