We introduce the Option Interpolation Model (OIM) for accurate approximation of embedded option values in insurance liabilities. Accurate approximation is required for ex-ante risk management applications. The OIM is based on interpolation with radial basis functions, which can interpolate scattered data, and does not suffer from the curse of dimensionality. To reduce computation time we present an inversion method to determine the interpolation function weights. The robustness, accuracy and efficiency of the OIM are investigated in several numerical experiments. We show that the OIM results in highly accurate approximations.