We report on an improved numerical method for hysteresis compensation within a hybrid piezoelectric fiber optic voltage sensor. The previous technique relied on laborious and relatively inaccurate individual curve-fitting procedures to create multiple sets of lookup functions used for the final derivation of the measured voltage. The new technique uses an aggregate 3-D arrangement of the captured hysteresis loops and applies a surface-fitting algorithm based on the Levenberg-Marquardt method to create two 3-D lookup functions, which are then used to derive the instantaneous value of the measured voltage. Furthermore, a more advanced algorithm for selecting top and bottom hysteresis-loop parts has been applied to eliminate errors associated ...