A recent study demonstrated that three-dimensional (3D) continuous displacement fields in transparent soft gels can be constructed from discrete displacement data obtained by optically tracking fluorescent particles embedded in the gels. Strain and stress fields were subsequently determined from gradients of the displacement field. This process was achieved through the moving least-square (MLS) interpolation method. The goal of this study is to evaluate the numerical accuracy of MLS in determining the displacement, strain, and stress fields in soft materials subjected to large deformation. Using an indentation model as the benchmark, we extract displacement at a set of randomly distributed data points from the results of a finite-element model, utilize these data points as the input for MLS, and compare resulting displacement, strain, and stress fields with the corresponding finite-element results. The calculation of strain and stress is based on finite strain kinematics and hyperelasticity theory. We also perform a parametric study in order to understand how parameters of the MLS method affect the accuracy of the interpolated displacement, strain, and stress fields. We further apply the MLS method to two additional cases with highly nonuniform deformation: a plate with a circular cavity subjected to large uniaxial stretch and a plane stress crack under large mode I loading. The results demonstrate the feasibility of using optical particle tracking together with MLS interpolation to map local strain and stress field in highly deformed soft materials.