firelab-general/ratio_pyrometry.py

import math
import cv2 as cv
import numpy as np
from numba import jit
import json

# camera settings
file = '01-0001.png'
I_Darkcurrent = 150.5
exposure_time = 0.500
f_stop = 2.4
ISO = 64 # basically brightness

# pyrometry config
MAX_GR_RATIO = 1200
MIN_GR_RATIO = 0
# original range from paper
# MAX_GR_RATIO = 1200
# MIN_GR_RATIO = 600

# Cropping config
x1 = 420
x2 = 1200
y1 = 400
y2 = -1

# post-processing
smoothing_radius = 2

@jit(nopython=True)
def rg_ratio_normalize(imgarr):
    tmin = MAX_GR_RATIO
    tmax = 0
    imgnew = imgarr
    for i in range(len(imgarr)):
        for j in range(len(imgarr[i])):
            px = imgarr[i][j]
            r_norm = normalization_func(px[0])
            g_norm = normalization_func(px[1])

            # apply camera calibration func
            temp_C = pyrometry_calibration_formula(g_norm, r_norm)

            # remove pixels outside calibration range
            if MAX_GR_RATIO != None and temp_C > MAX_GR_RATIO or MIN_GR_RATIO != None and temp_C < MIN_GR_RATIO:
                temp_C = 0

            # update min & max
            if temp_C < tmin and temp_C >= 0:
                tmin = temp_C
            if temp_C > tmax:
                tmax = temp_C

            imgnew[i][j] = [temp_C, temp_C, temp_C]
    return imgnew, tmin, tmax


@jit(nopython=True)
def normalization_func(i):
    """
    does something to the pixels that i don't understand lol
    """
    return (i - I_Darkcurrent) * (f_stop ** 2) / (ISO * exposure_time)


@jit(nopython=True)
def pyrometry_calibration_formula(i_ng, i_nr):
    """
    Given the green-red ratio, calculates an approximate temperature 
    in Celsius.
    """
    return 362.73 * math.log10(
        (i_ng/i_nr) ** 3
    ) + 2186.7 * math.log10(
        (i_ng/i_nr) ** 3
    ) + 4466.5 * math.log10(
        (i_ng / i_nr) ** 3
    ) + 3753.5


# read image & crop
file_name = file.split(".")[0]
file_ext = file.split(".")[1]
img = cv.imread(file)
img = img[y1:y2, x1:x2]
cv.imwrite(f'{file_name}-cropped.{file_ext}', img)

# img = cv.imread('ember_test.png')

img, tmin, tmax = rg_ratio_normalize(img)

print(f"min: {tmin}°C")
print(f"max: {tmax}°C")

# build & apply smoothing conv kernel
k = []
for i in range(smoothing_radius):
    k.append([1/(smoothing_radius**2) for i in range(smoothing_radius)])
    # for j in range(smoothing_radius):
kernel = np.array(k)

img = cv.filter2D(src=img, ddepth=-1, kernel=kernel)

# apply jet color map
img = cv.applyColorMap(img, cv.COLORMAP_JET)

cv.imwrite(f'{file_name}-cropped-transformed-ratio.{file_ext}', img)