Integral Image

import matplotlib.pyplot as plt
from PIL import Image

img = Image.open("Cameraman_noise.bmp").convert("L")

integral_img = img.copy()

value = 0
w, h = img.size[0], img.size[1];
new_image = [[0 for x in range(h)] for y in range(w)]
integral_array = [[0 for x in range(h)] for y in range(w)]

# CALCULATE FOR ALL PIXELS THAT ARE ON THE LEFT
for i in range (0,img.size[0]):
    value = 0
    for j in range (0,img.size[1]):
        value += img.getpixel((i, j))
        new_image[i][j] = value

# CALCULATE FOR ALL PIXELS THAT ARE ABOVE
for i in range (0,img.size[1]):
    value = 0
    for j in range (0,img.size[0]):
        value += new_image[j][i]
        integral_array[j][i] = value

value = 0

for i in range (0,img.size[0]):
    for j in range (0,img.size[1]):
        value = (255-0)*(integral_array[i][j] - integral_array[0][0])//(integral_array[img.size[0]-1][img.size[1]-1] + integral_array[0][0])
        integral_img.putpixel((i,j),(int) (value))



plt.imshow(integral_img)
integral_img.save('Camera_Integ.jpg')