firelab-general/examples/pyrometry/formula_testing.py

104 lines
2.0 KiB
Python
Raw Normal View History

2022-10-12 21:22:39 -07:00
import math
rg_value_sets = {
"firebrand_test.png": [
(219, 7),
(227, 14),
(166, 14),
(197, 10),
(230, 25),
(228, 17),
(218, 17),
(221, 15),
(210, 22),
(229, 17),
],
2022-10-12 21:22:39 -07:00
"streaktest.png": [
(50, 11),
(51, 12),
(52, 10),
(240, 115),
(254, 127),
],
"ember_orange.png": [
(240, 147),
(235, 102),
(223, 103),
(232, 103),
(103, 34),
(128, 47),
(92, 27),
2022-10-19 18:36:45 -07:00
],
"test5.png": [
(164, 49),
(177, 57),
(181, 58),
(156, 46),
(150, 41),
(151, 44),
],
"DSC1022 (Dennis Calibration Image, 1200°C)": [
(252, 84),
(252, 77),
(231, 67),
(215, 61),
(182, 48),
(103, 17),
(97, 15),
(63, 6),
]
2022-10-19 18:36:45 -07:00
}
2022-10-12 21:22:39 -07:00
# Settings
I_Darkcurrent = 7.7
exposure_time = 0.500
f_stop = 2.4
ISO = 100 # basically brightness
def pyrometry_calibration_formula(i_ng, i_nr):
"""
Given the green-red ratio, calculates an approximate temperature
in Celsius.
"""
try:
2022-10-16 15:46:26 -07:00
return (
362.73 * math.log10(i_ng/i_nr) ** 3 +
2186.7 * math.log10(i_ng/i_nr) ** 2 +
4466.5 * math.log10(i_ng / i_nr) +
3753.5
)
2022-10-12 21:22:39 -07:00
except:
return 'dropped'
def tprint(*items):
for item in items:
print(item, end="\t")
print()
def grtemp(px):
r_norm = (px[0] - I_Darkcurrent) * (f_stop ** 2) / (ISO * exposure_time)
g_norm = (px[1] - I_Darkcurrent) * (f_stop ** 2) / (ISO * exposure_time)
res = pyrometry_calibration_formula(g_norm, r_norm)
tprint(
px[0],
px[1],
# round(px[0] / px[1], 2),
2022-10-12 21:22:39 -07:00
round(r_norm, 2),
round(g_norm, 2),
2022-10-16 15:46:26 -07:00
round(g_norm / r_norm, 4),
2022-10-12 21:22:39 -07:00
res,
)
for (key, val) in rg_value_sets.items():
print(f"\n{key}\n")
tprint('RED', 'GREEN', 'RNORM', 'GNORM', 'G_n/R_n', 'RES TEMP')
for rg in val:
grtemp(rg)