import sys
import struct
import numpy as np
import time

s = time.time()

r, g, b = 0, 0, 0
r, g, b = 1, 1, 1
w=4192
h=3104

pixels = np.fromfile(sys.argv[1], dtype="<u2").reshape(h,w)
d = time.time()
for rg_row in pixels[::2]:
    r += np.sum(rg_row[::2])
    g += np.sum(rg_row[1::2])
for gb_row in pixels[1::2]:
    g += np.sum(gb_row[::2])
    b += np.sum(gb_row[1::2])

r=r/w/h*4
g=g/w/h*2
b=b/w/h*4
m=(r+g+b)/3
e = time.time()
print('\rR', round(r, 2), round(m/r, 3), round(g/r, 3))
print('G', round(g, 2), round(m/g, 3), '1.000')
print('B', round(b, 2), round(m/b, 3), round(g/b, 3))
print("total time: %.3f seconds. Data read/prepare: %.3f" % ((e - s), (d - s)))