g2Num = np.zeros_like(imgSomma, dtype = np.float64)

for i,e in enumerate(os.scandir(r".\1slit_thermal_npz")):
    tmpImg = np.load(e)['arr_0']
    print(i)

    if e.name in strBuio: continue

    tmpImg = tmpImg - matNoise

    Ix1 = tmpImg[110, 325]
    g2Num += tmpImg * Ix1

g2Num = g2Num/numFrame

#avgx1 = np.sum(imgSommaFilt[110:121, 325:356]) / numFrame
avgx1 = imgSommaFilt[110, 325] / numFrame
avgx2 = imgSommaFilt / numFrame

g2 = g2Num / (avgx1 * avgx2)