Горизонтальные полосы почти исчезли, время от времени вылезает непонятная сетчатая рябь вокруг руки.

К своему удивлению, я обнаружил зависимость помех от положения фильтра в пространстве: иногда они полностью отсутствуют, но если повернуть его чуть в сторону, то снова появляется рябь.

print('Les Wright 21 June 2023')
print('https://youtube.com/leslaboratory')
print('A Python program to read, parse and display thermal data from the Topdon TC001, Infiray P2 Pro and other compatible USB Thermal Cameras!')
print('')
print('THIS IS A PATCHED VERSION; patch created by CATCH22MANIA \'07/24 - https://www.youtube.com/@catch22mania')
print('')
print('Key bindings:')
print('')
print('c x: cycle through colormaps - 0 = raw (no hud)')
print('f w: fullscreen, windowed')
print('h : toggle HUD')
print('q : quit')
print('')
print('1 : show highest temperature')
print('2 : show centre crosshairs temperature')
print('3 4: change scale - in fullscreen: zoom in/out font')
print('5 6: +/- contrast')
print('7 8: +/- blur')
# print('9 0: floating high and low temp label threshold')
print('r t: Record and Stop')
# print('p : Snapshot')
print('')
print('check out the command line options > python ~/testing/ptc-patched-c22-014-test1.py -h')
print('')

import cv2
import numpy as np
#print(np.__version__)
import argparse
import time
import io

#256x192 General settings
width = 256 #Sensor; width for output window size (* scale)
height = 192 #sensor; height for output windows size (* scale)
scale = 3 #scale multiplier
alpha = 1.0 # Contrast control (1.0-3.0)
colormap = 0
font=cv2.FONT_HERSHEY_SIMPLEX
dispFullscreen = False
rad = 0 #blur radius
threshold = 2
hud = False
recording = False
showcrosshairs = False
showtemps = False
elapsed = "00:00:00"
snaptime = "None"


#We need to know if we are running on the Pi, because openCV behaves a little oddly on all the builds!
#https://raspberrypi.stackexchange.com/questions/5100/detect-that-a-python-program-is-running-on-the-pi
def is_raspberrypi():
    try:
        with io.open('/sys/firmware/devicetree/base/model', 'r') as m:
            if 'raspberry pi' in m.read().lower(): return True
    except Exception: pass
    return False

isPi = is_raspberrypi()

parser = argparse.ArgumentParser()
parser.add_argument("--device", type=int, default=0, help="video device number; e.g. 0, use v4l2-ctl --list-devices")
parser.add_argument("--width", type=int, default=0, help="video output width; e.g. 720, scale to width (scale factor = 1)")
parser.add_argument("--height", type=int, default=0, help="video output height; e.g. 576, scale to height (scale factor = 1)")
parser.add_argument("--fullscreen", dest="fullscreen",action='store_true', help="start in fullscreen; no argument")
parser.add_argument("--colormap", type=int, default=3, help="colormap; range: 1-22; 0 = no colormap (raw); default: 3")
parser.add_argument("--crosshairs", dest="crosshairs", action='store_true', help="show crosshairs; no argument")
parser.add_argument("--temps", dest="temps", action='store_true', help="show temperatures; no argument")
parser.add_argument("--blur", type=int, default=0, help="blur radius; default: 0")
parser.add_argument("--contrast", type=float, default=1.0, help="contrast; default: 1.0")

args = parser.parse_args()
	
if args.device:
	dev = args.device
else:
	dev = 0
if args.width:
	scale = 1
	width = args.width	
if args.height:
	scale = 1
	height = args.height	
if args.fullscreen:
	dispFullscreen = args.fullscreen
if args.colormap:
	colormap = args.colormap	
if args.blur:
	rad = args.blur	
if args.contrast:
	alpha = args.contrast		
if args.crosshairs:
	showcrosshairs = args.crosshairs
if args.temps:
	showtemps = args.temps

newWidth = width*scale 
newHeight = height*scale


cap = cv2.VideoCapture('/dev/video'+str(dev), cv2.CAP_V4L)
#pull in the video but do NOT automatically convert to RGB, else it breaks the temperature data!
#https://stackoverflow.com/questions/63108721/opencv-setting-videocap-property-to-cap-prop-convert-rgb-generates-weird-boolean
if isPi == True:
	cap.set(cv2.CAP_PROP_CONVERT_RGB, 0.0)
else:
	cap.set(cv2.CAP_PROP_CONVERT_RGB, False)
	

cv2.namedWindow('ThermalWindow',cv2.WND_PROP_FULLSCREEN)
cv2.setWindowProperty('ThermalWindow',cv2.WND_PROP_FULLSCREEN,cv2.WINDOW_FULLSCREEN)
if dispFullscreen == False:
	cv2.resizeWindow('ThermalWindow', newWidth,newHeight)

current_time_str = time.strftime("%Y%m%d--%H%M%S")
fourcc = cv2.VideoWriter_fourcc(*'XVID')
out = cv2.VideoWriter(f"thermal_output-{current_time_str}.avi", fourcc, 25, (newWidth, newHeight))

while(cap.isOpened()):
	# Capture frame-by-frame
	ret, frame = cap.read()
	if ret == True:
		imdata,thdata = np.array_split(frame, 2)
		hi = thdata[96][128][0]
		lo = thdata[96][128][1]
		lo = lo*256
		rawtemp = hi+lo
		temp = (rawtemp/64)-273.15
		temp = round(temp)

		#find the max temperature in the frame
		wat = np.add(thdata[...,0], np.multiply(thdata[...,1], 256))
		maxtemp = wat[...].max()
		posmax = wat[...].argmax()
		mcol, mrow  = divmod(posmax, 256)
		maxtemp = (maxtemp/64) - 273.15
		maxtemp = round(maxtemp)
		
		# hotspot: quick n dirty for testing - scaling width, height
		mcol2 = mcol / (192 / height) 
		mrow2 = mrow / (256 / width)		
		mcol = round (mcol2)
		mrow = round (mrow2)				

		bgr = cv2.cvtColor(imdata,  cv2.COLOR_YUV2BGR_YUYV)

		bgr = cv2.resize(bgr,(newWidth,newHeight),interpolation=cv2.INTER_CUBIC)#Scale up!

		if colormap != 0: # 0 = no colormap, no blur, no contrast
				
			# contrast
			if alpha> 1.0:					
				bgr = cv2.convertScaleAbs(bgr, alpha=alpha)#Contrast
			# blur
			if rad>0:
				bgr = cv2.blur(bgr,(rad,rad))

			# apply colormap
			if colormap == 1:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_AUTUMN)
				cmapText = 'autumn'
			if colormap == 2:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_BONE)
				cmapText = 'bone'
			if colormap == 3:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_JET)
				cmapText = 'jet'			
			if colormap == 4:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_WINTER)
				cmapText = 'winter'			
			if colormap == 5:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_RAINBOW)
				cmapText = 'rainbow'			
			if colormap == 6:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_OCEAN)
				cmapText = 'ocean'			
			if colormap == 7:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_SUMMER)
				cmapText = 'summer'			
			if colormap == 8:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_SPRING)
				cmapText = 'spring'			
			if colormap == 9:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_COOL)
				cmapText = 'cool'			
			if colormap == 10:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_HSV)
				cmapText = 'HSV'			
			if colormap == 11:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_PINK)
				cmapText = 'pink'			
			if colormap == 12:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_HOT)
				cmapText = 'hot'
			if colormap == 13:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_PARULA)
				cmapText = 'parula'			
			if colormap == 14:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_MAGMA)
				cmapText = 'magma'			
			if colormap == 15:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_INFERNO)
				cmapText = 'inferno'			
			if colormap == 16:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_PLASMA)
				cmapText = 'plasma'			
			if colormap == 17:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_VIRIDIS)
				cmapText = 'viridis'			
			if colormap == 18:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_CIVIDIS)
				cmapText = 'cividis'
			if colormap == 19:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_TWILIGHT)
				cmapText = 'twilight'			
			if colormap == 20:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_TWILIGHT_SHIFTED)
				cmapText = 'twilight s.'			
			if colormap == 21:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_TURBO)
				cmapText = 'turbo'			
			if colormap == 22:
				heatmap = cv2.applyColorMap(bgr, cv2.COLORMAP_DEEPGREEN)
				cmapText = 'deepgreen'

			if showcrosshairs == True:
				# draw crosshairs
				cv2.line(heatmap,(int(newWidth/2),int(newHeight/2)+20),\
				(int(newWidth/2),int(newHeight/2)+5),(255,255,255),1) #vline 1
				cv2.line(heatmap,(int(newWidth/2),int(newHeight/2)-5),\
				(int(newWidth/2),int(newHeight/2)-20),(255,255,255),1) #vline 2
							
				cv2.line(heatmap,(int(newWidth/2)-5,int(newHeight/2)),\
				(int(newWidth/2)-20,int(newHeight/2)),(255,255,255),1) #hline 1
				cv2.line(heatmap,(int(newWidth/2)+5,int(newHeight/2)),\
				(int(newWidth/2)+20,int(newHeight/2)),(255,255,255),1) #hline 2			
			
				cv2.putText(heatmap,str(temp)+" C", (int(newWidth/2)+10, int(newHeight/2)-10),\
				cv2.FONT_HERSHEY_SIMPLEX, 0.55,(0, 0, 0), 2, cv2.LINE_AA)
				cv2.putText(heatmap,str(temp)+" C", (int(newWidth/2)+10, int(newHeight/2)-10),\
				cv2.FONT_HERSHEY_SIMPLEX, 0.55,(0, 255, 255), 1, cv2.LINE_AA)

			if hud==True:
				# display black box for our data
				cv2.rectangle(heatmap, (0, 0),(180, 55), (0,0,0), -1)
				# put text in the box
		
				cv2.putText(heatmap,'Colormap: '+str(colormap)+' ('+cmapText+')', (10, 16),\
				cv2.FONT_HERSHEY_SIMPLEX, 0.4,(0, 255, 255), 1, cv2.LINE_AA)

				cv2.putText(heatmap,'Blur: '+str(rad)+' ', (10, 30),\
				cv2.FONT_HERSHEY_SIMPLEX, 0.4,(0, 255, 255), 1, cv2.LINE_AA)

				cv2.putText(heatmap,'Contrast: '+str(alpha)+' ', (10, 44),\
				cv2.FONT_HERSHEY_SIMPLEX, 0.4,(0, 255, 255), 1, cv2.LINE_AA)
					
			if showtemps == True:
				cv2.circle(heatmap, (mrow*scale, mcol*scale), 5, (0,0,0), 1)
				cv2.putText(heatmap,str(maxtemp)+' C', ((mrow*scale)+10, (mcol*scale)-10),\
				cv2.FONT_HERSHEY_SIMPLEX, 0.55,(0,0,0), 2, cv2.LINE_AA)
				cv2.putText(heatmap,str(maxtemp)+' C', ((mrow*scale)+10, (mcol*scale)-10),\
				cv2.FONT_HERSHEY_SIMPLEX, 0.55,(0, 255, 255), 1, cv2.LINE_AA)

				cv2.line(heatmap,(int(mrow*scale),int(mcol*scale)+20),\
				(int(mrow*scale),int(mcol*scale)+5),(255,255,255),1) #vline 1			
				cv2.line(heatmap,(int(mrow*scale),int(mcol*scale)-5),\
				(int(mrow*scale),int(mcol*scale)-20),(255,255,255),1) #vline 2			
				cv2.line(heatmap,(int(mrow*scale)-5,int(mcol*scale)),\
				(int(mrow*scale)-20,int(mcol*scale)),(255,255,255),1) #hline 1
				cv2.line(heatmap,(int(mrow*scale)+5,int(mcol*scale)),\
				(int(mrow*scale)+20,int(mcol*scale)),(255,255,255),1) #hline 2

			#display image
			cv2.imshow('ThermalWindow',heatmap) # rendered image
			out.write(heatmap)
		else:
			cv2.imshow('ThermalWindow',bgr)		# "raw" frame	
			out.write(bgr)

		keyPress = cv2.waitKey(1)
		if keyPress == ord('7'): #Increase blur radius
			rad += 1
		if keyPress == ord('8'): #Decrease blur radius
			rad -= 1
			if rad <= 0:
				rad = 0

		if keyPress == ord('3'): #Increase scale
			scale += 1
			if scale >=5:
				scale = 5
			newWidth = width*scale
			newHeight = height*scale
			if dispFullscreen == False and isPi == False:
				cv2.resizeWindow('ThermalWindow', newWidth,newHeight)
		if keyPress == ord('4'): #Decrease scale
			scale -= 1
			if scale <= 1:
				scale = 1
			newWidth = width*scale
			newHeight = height*scale
			if dispFullscreen == False and isPi == False:
				cv2.resizeWindow('ThermalWindow', newWidth,newHeight)

		if keyPress == ord('f'): #enable fullscreen
			dispFullscreen = True
			cv2.namedWindow('ThermalWindow',cv2.WND_PROP_FULLSCREEN)
			cv2.setWindowProperty('ThermalWindow',cv2.WND_PROP_FULLSCREEN,cv2.WINDOW_FULLSCREEN)
		if keyPress == ord('w'): #disable fullscreen
			dispFullscreen = False
			cv2.namedWindow('ThermalWindow',cv2.WINDOW_GUI_NORMAL)
			cv2.setWindowProperty('ThermalWindow',cv2.WND_PROP_AUTOSIZE,cv2.WINDOW_GUI_NORMAL)
			cv2.resizeWindow('ThermalWindow', newWidth,newHeight)

		if keyPress == ord('5'): #contrast+
			alpha += 0.1
			alpha = round(alpha,1)#fix round error
			if alpha >= 3.0:
				alpha=3.0
		if keyPress == ord('6'): #contrast-
			alpha -= 0.1
			alpha = round(alpha,1)#fix round error
			if alpha<=0:
				alpha = 0.0
 
		if keyPress == ord('h'):
			if hud==True:
				hud=False
			elif hud==False:
				hud=True

		if keyPress == ord('c'): #c to cycle through colormaps
			colormap += 1
			if colormap  > 22:
				colormap = 0
				
		if keyPress == ord('x'): #x to cycle through colormaps
			colormap -= 1
			if colormap < 0:
				colormap = 22				

		if keyPress == ord('q'):
			break
			out.release()
			capture.release()
			cv2.destroyAllWindows()
			exit(0)

		if keyPress == ord('1'):
			if showtemps==True:
				showtemps=False
			elif showtemps==False:
				showtemps=True

		if keyPress == ord('2'):
			if showcrosshairs==True:
				showcrosshairs=False
			elif showcrosshairs==False:
				showcrosshairs=True

Не знаю, как это работает, но за счёт более узкого угла обзора картинка на глаз выглядит резче и чётче. Секрета здесь нет: андроид-приложение растягивает картинку до 1080х810, в то время как "малина" честно умножает 256х192 датчика на 3 (в нашем коде scale = 3), отдавая в итоге 768х576 (что близко к разрешению шлема).

Наблюдаем за людьми и машинами: