A technique based on the combination of Fourier pseudospectral time-domain method (PSTD) and PML absorbing boundary conditions is widely used to simulate the scattering and propagation problems. However, one of the known disadvantages of this method is that it requires a uniformly distributed spatial grid-set along each orthogonal direction. But its accuracy is lower when applied to the objects involving the curved boundary or thein dimensions are not commensurate with the cell size. The PSTD method involving space transformations (TSNU-PSTD) alleviates this limitation. The CFS-PML is implemented in the PSTD and the TSNU-PSTD algorithm combined with the CFS-PML has been applied to model ground-penetrating radar (GPR) application involving objects in a large-scale lossy half space, and the results are compared with the FDTD. Numerical simulations show that it can achieve a satisfactory result with a grid density of only 3.0 nodes per minimum wavelength, and can simulate the scattering of curved objects in large-scale lossy media efficaciously