double>(i,j) = imCum.at<double>(i,j+r); } } //imDst(:, r+2:wid-r) = imCum(:, 2*r+2:wid) - imCum(:, 1:wid imCum.at<double>(i,j-r-1); } } //imDst(:, wid-r+1:wid) = repmat(imCum(:, wid), [1, r]) - imCum(:, wid