Abstract: To avoid the difficulty of finite element mesh generation for complex sites, and thus the difficulty of site seismic response analysis, this study proposes a complex site modeling and seismic response analysis method based on borehole and surface elevation data. Sampling points are selected along boreholes to construct a dataset of scattered spatial coordinates and material properties. A spatial scattered-point interpolation algorithm then converts these data into a digital image representing soil layer distribution at the site. Surface elevation data are integrated to adjust the image pixels, accurately capturing terrain geometry. The adjusted image is processed via quadtree/octree algorithms to generate high-quality meshes. The scaled boundary finite element method is employed to calculate the seismic response of the quadtree/octree mesh of the site. The proposed framework achieves high level of automation and mesh quality, offering an efficient solution for complex site analysis. The accuracy and effectiveness of the proposed method have been verified through numerical examples with known soil stratigraphy, and the method has been applied to the modeling and seismic response analysis of a real-world case study.