To study the polarization reflection characteristics of metal surfaces, a parameter optimization method for the polarization bidirectional reflection distribution function (pBRDF) model of metal surfaces based on the improved strawberry algorithm has been proposed. Firstly, the light scattering characteristics of metal surfaces were analyzed and a multi-parameter pBRDF model was constructed. Then, the working mechanism of the strawberry optimization algorithm was investigated and improved by introducing the chaotic mapping and Levy flight strategy to overcome the shortcomings, such as low convergence rate and easily falling into local optimum. Finally, the method proposed in this paper was validated by simulating open-source data from references and the obtained ones with a self-built experimental platform. The results show that the proposed method outperforms those by nonlinear least squares, particle swarm optimization and the original strawberry algorithm in fitting the detected degree of polarization (DOP) data, indicating the modeling accuracy was significantly improved and better suited to characterize the polarized reflection properties of metal surfaces.