Parametric imaging using the Patlak graphical method has been widely used to analyze dynamic PET data. Conventionally a Patlak parametric image is generated by reconstructing a sequence of dynamic images first and then performing Patlak graphical analysis on the time-activity curves pixel-by-pixel. However, because it is rather difficult to model the noise distribution in reconstructed images, the spatially variant noise correlation is simply ignored in the Patlak analysis, which leads to sub-optimal results. In this paper we present a Bayesian method for reconstructing Patlak parametric images directly from raw sinogram data by incorporating the Patlak plot model into the image reconstruction procedure. A preconditioned conjugate gradient algorithm is used to find the maximum a posteriori solution. The proposed direct method is statistically more efficient than the conventional indirect approach because the Poisson noise distribution in PET data can be accurately modeled in the direct reconstruction. The computation cost of the direct method is similar to reconstruction time of two dynamic frames. Therefore, when more than two dynamic frames are used in the Patlak analysis, the direct method is faster than the conventional indirect approach. We conduct computer simulations to validate the proposed direct method. Comparisons with the conventional indirect approach show that the proposed method results in a more accurate estimate of the parametric image. The proposed method has been applied to dynamic fully 3D PET data from a microPET scanner.