Abstract A method for calculating the two‐dimensional scattering of incident SH waves by canyons of arbitrary shape is presented. The problem is formulated in terms of a Fredholm integral equation of the first kind with the integration path outside the boundary. Point‐source discretization and a least‐squares scheme are used. Numerical results are compared with the known analytic solution for a semi‐cylindrical canyon. Spatial variations of surface amplitudes are computed for triangular and half‐cycle sinusoidal canyons as well.