An interfacial regularized Stokeslet scheme is presented to predict the motion of solid bodies (e.g. proteins or gel-phase domains) embedded within flowing lipid bilayer membranes. The approach provides a numerical route to calculate velocities and angular velocities in complex flow fields that are not amenable to simple Fax\'en-like approximations. Additionally, when applied to shearing motions, the calculations yield predictions for the effective surface viscosity of dilute rigid body-laden membranes. In the case of cylindrical proteins, effective viscosity calculations are compared to two prior analytical predictions from the literature. Effective viscosity predictions for a dilute suspension of rod-shaped objects in the membrane are also presented.