place_denoise.py

#!/usr/bin/env python3

# grab data from https://place.thatguyalex.com/

import glob

import numpy
import PIL.Image

import scipy.stats

EXPECTED_PALETTE = None

END = 1649112396

ims = []
qs = []

for n in range(4):
    si = []
    fs = []
    for f in sorted(glob.glob(f"images_quadro/{n}/*-164911*")):
        if END - 3600 <= int(f.split('.')[0].split('-')[-1]) <= END:
            fs.append(f)

    for f in fs:
        i = PIL.Image.open(f)
        if not EXPECTED_PALETTE:
            EXPECTED_PALETTE = i.palette.getdata()
        assert i.palette.getdata() == EXPECTED_PALETTE
        si.append(i)

    print(n, len(fs))

    quadrant = numpy.dstack([numpy.array(i) for i in si])
    mode = scipy.stats.mode(quadrant, axis=2).mode
    mode = mode.transpose()[0].transpose()

    im = PIL.Image.fromarray(mode)
    im.putpalette(EXPECTED_PALETTE[1])
    im.save(f"Q{n}.png")
    qs.append(im)

    ims.extend(fs)


print(len(ims))

out = PIL.Image.new('P', (2000, 2000))
out.putpalette(EXPECTED_PALETTE[1])

for n, q in enumerate(qs):
    out.paste(q, (n % 2 * 1000, n // 2 * 1000))

out.save('averaged.png')