In the sections of this chapter the reader will be introduced to the sequence of mathematical processes which, starting from a model interpolation function yield to the corresponding SRE-based paradigm. Particularly, this chapter addresses the development of the SRE-based bivariate interpolation function. The mathematical procedure is consistently iterated in the rest of the book for all the other model functions that the unifying theory embraces. The first step of the procedure is that of the calculation of the intensity-curvature terms and through their ratio the Intensity-Curvature Functional is calculated for the model function. The second step is that of calculating the first order partial derivatives of the Intensity-Curvature Functional. Thirdly, the polynomial consisting of the first order partial derivatives is solved to obtain the Sub-pixel Efficacy Region. At this point, the formula of the unifying theory (equation [21]) sets the stage to obtain the novel re-sampling locations. Worth noting that this formula can be adapted to cover cases of one, two, and three dimensional interpolation functions and it is also consistently employed for linear quadratic cubic and trigonometric (Sinc) models. This shall be manifest throughout the remainder of the book. The remainder of this chapter discusses on the nature of the SRE, also makes a connection with Chapter XX of the book within the context of the relationship existing between resolution and interpolation error, and in the last section, the concept of resilient interpolation is introduced and the relevant math is illustrated for the case of the bivariate linear function.