Coseismic landslides pose significant threats, causing widespread destruction of buildings, roads, pipelines, and leading to numerous casualties. In recent years, the frequency of earthquakes has increased, prompting a growing interest in regional-scale assessment techniques for coseismic landslides. The infinite slope model proposed by Newmark is widely used to evaluate coseismic landslide hazard. However, the infinite slope model fails to reflect the impact of rock mass structure on the stability of slopes. This paper proposes a novel approach for mapping the hazards of coseismic landslides by considering the roughness of the potential slide surfaces in the inner slope. The proposed method is validated using data from the 2013 Lushan earthquake. The datasets, including geological units, peak ground acceleration (PGA), and high-resolution digital elevation models of topography, are rasterized at a grid spacing of 30 meters. They are then combined within an infinite slope model based on Newmark permanent-deformation analysis, enabling the estimation of coseismic landslide displacement in each grid area resulting from the Lushan earthquake. The modeled displacements are compared with the inventory of landslides triggered by the Lushan earthquake, allowing the derivation of a confidence level function that relates predicted displacement to the spatial variation of coseismic landslides. Ultimately, a hazard map of coseismic landslides is generated based on the established confidence level function. This map serves as a valuable tool for predicting the hazard zone of seismic regions and offers essential insights for decision-making related to infrastructure development and post-earthquake construction.