We introduce a new method for implicit structural modeling. The method belongs to the Discrete Smooth Interpolation class of methods. We exploit the highly symmetric nature of Cartesian grids to propose new regularization operators that discretize very naturally on the Cartesian grid using finite differences: these operators do not have to be defined on boundary nodes, and their generalization to higher dimensions is straightforward. Numerical examples show that the proposed method is both robust and numerically efficient.