Storage phosphor based imaging detectors are widely used in various fields of X-ray radiography. In this work, a numerical procedure based on the Monte Carlo algorithm is presented to simulate the image readout of storage phosphor plates (SPPs). Several experimentally measurable model parameters, including the optical anisotropy coefficient, the optical absorption coefficient and the mean scattering length, are used to deal with multiple light-scattering processes during the image readout. The density distribution profile of the readout laser and escape probability of the luminescence photons are simulated. Such image metrics as line spread function, photostimulated luminescence (PSL) sensitivity and modulation transfer function (MTF) can be calculated. The simulated MTFs and PSL sensitivities are in good agreement with experimental data at different SPP thickness within 5%.