The purpose of the present paper is to model a physical polarimetric bidirectional reflectance distribution function (pBRDF), which can character not only the non-Lambertian but also the polarized features in order that the pBRDF can be applied to analyze the relationship between the degree of polarization and the physiological and biochemical parameters of leaves quantitatively later. Firstly, the bidirectional polarized reflectance distributions from several leaves surfaces were measured by the polarized goniometer developed by Changchun Institute of Optics, Fine Mechanics and Physics, Chinese Academy of Sciences. The samples of leaves include two pieces of zea mays L. leaves (young leaf and mature leaf) and a piece of E. palcherrima wild leaf. Non-Lambertian characteristics of directional reflectance from the surfaces of these three leaves are obvious. A Cook-Torrance model was modified by coupling the polarized Fresnel equations to simulate the bidirectional polarized reflectance properties of leaves surfaces. The three parameters in the modified pBRDF model, such as diffuse reflectivity, refractive index and roughness of leaf surface were inversed with genetic algorithm (GA). It was found that the pBRDF model can fit with the measured data well. In addition, these parameters in the model are related with both the physiological and biochemical properties and the polarized characteristics of leaves, therefore it is possible to build the relationships between them later.