In this paper, we present a novel method for the identification of the local bending stiffness of a beam. We use shearography to capture measurements of vibrating beams, so the input data for the identification is the modal slope-the differential of the modal shape. The modal slope is fitted by two Fourier-series functions, one of which is derived from a thin-beam model. The local bending stiffness is identified as the one corresponding with the best match between the measured and the two fitted modal slopes. This identification method, which we call simultaneous Fourier-series fitting, is demonstrated on numerically-generated inputs, as well as on experimental measurements. We use a flat, concave and convex beam, as well as beams with locally varying bending stiffness mimicking local damage to verify the method. It is shown that the method gives accurate results and is robust to noise. Additionally, it has advantageous properties that make it useful and practical: using this method, it is possible to perform the identification from only a sub-region of a beam and even without specifying the boundary conditions.