([formula   1] I(x,y)=(∑_(x'<x,y'<y))S(x',y')(1))