Double Linear Interpolation
x = s - s0, 0 x 1
y = l - l0, 0 y 1
s1 = s0 + 1, l1 = l0 + 1
DN(s,l0) = DN(s0,l0) + x · (DN(s1,l0)-DN(s0,l0))
DN(s,l1) = DN(s0,l1) + x · (DN(s1,l1)-DN(s0,l1))
DN(s,l) = DN(s,l0) + y · (DN(s,l1) - DN(s,l0))
or
DN(s,l) = DN(s0,l0)