We consider fast and accurate solution methods for the direct and inverse scattering problems by a few three dimensional piecewise homogeneous dielectric obstacles around the resonance region. The forward problem is reduced to a system of second kind boundary integral equations. For the numerical solution of these coupled integral equations we modify a fast and accurate spectral algorithm, proposed by Ganesh and Hawkins [doi:10.1016/j.jcp.2008.01.016], by transporting these equations onto the unit sphere using the Piola transform of the boundary parametrisations. The computational performances of the forward solver are demonstrated on numerical examples for a variety of three-dimensional smooth and non smooth obstacles. The algorithm, that requires the knowledge of the boundary parametrisation and leads