Subsurface scattering is ubiquitous in organic materials and has been widely researched in computer graphics. Inverse rendering of subsurface scattering, however, is often constrained by the planar geometry assumption of traditional analytic Bidirectional Scattering Surface Reflectance Distribution Functions (BSSRDF). To address this issue, a shape-adaptive BSSRDF model has been proposed to render translucent objects on curved geometry with high accuracy. In this paper, we leverage this model to estimate parameters of subsurface scattering for inverse rendering. We compute the finite difference of the rendering equation for subsurface scattering and iteratively update material parameters. We demonstrate the performance of our shape-adaptive inverse rendering model by analyzing the estimation accuracy and comparing to inverse rendering with plane-based BSSRDF models and volumetric methods.