Hi Markus,
I was also wondering if you were seeing a 'rounding error' - after all, those original four values are based on the 'average' ADU of all the pixels on an image-by-image basis.
However, plain ImageIntegration returns the 'average' ADU value of the four pixel values at every location in the image. Only then do you look at the statistics, and expect that the two processes should reurn the same values for you.
Further, you are inspecting the ADU data using rational numbers - it might be interesting to see what the values look like if you view the data as 16-bit integers. For example, how much of a change in value do you actually get when viewing as integers? Is this change still enough to give you concern?
Bear in mind that your observed error is only of the order of 0.001%, or just under 1 ADU count in 65,535 (the full range of a 16-bit, unsigned, integer). So, should you really be concerned?
Hope this helps.