from PIL import Image
import numpy as np

a = 40

datei = "archimedes.jpg"

name = datei.split('.')[0]

bild = Image.open(datei) 
breite, hoehe = bild.size
print(breite, hoehe)

pixel_bild = bild.load()

breite_neu = breite // a
hoehe_neu = hoehe // a
print(breite_neu * a, hoehe_neu * a)
verpixelt = Image.new("RGB", (breite_neu * a, hoehe_neu * a))
pixel_verpixelt = verpixelt.load()

for x in range(breite_neu):
    for y in range(hoehe_neu):
        summe_rgb = np.array([0, 0, 0])
        for i in range(a):
            for j in range(a):
                summe_rgb = summe_rgb + pixel_bild[a*x+i, a*y+j]
        mittelwert = summe_rgb / (a*a)
        mittelwert = mittelwert.astype(int)
        for i in range(a):
            for j in range(a):
                pixel_verpixelt[a*x+i, a*y+j] = tuple(mittelwert)
verpixelt.save(f'{name}_verpixelt-mit-numpy.jpg')