p8.py

#!/usr/bin/env python3
# -*- coding: utf8 -*-

import numpy as np

m = np.fromfile('8.in', dtype=np.uint8, count=25*6*100).reshape((100,6,25)) - 48
f = m.reshape((m.shape[0], m.shape[1]*m.shape[2]))
best_layer = np.argmin(f.shape[1] - np.count_nonzero(f, axis=1))
print("Pt 1:", (f.shape[1] - np.count_nonzero(f[best_layer,:]-1)) * (f.shape[1] - np.count_nonzero(f[best_layer,:]-2)))

print("Pt 2:")
for r in range(m.shape[1]):
    for c in range(m.shape[2]):
        color = m[np.argwhere(m[:,r,c] != 2)[0],r,c]
        print('*' if color == 1 else ' ', end='')
    print(' ')