What is going on with these scene cuts?

[LaserBones]

I'm trying to encode this DVD. No matter what i do with the fields I get these weird glitched out frames in many of the scene cuts. I'm not sure what to ask since i don't understand what's going on. I can post other examples or a longer clip if needed.

this is the current script I'm using:
QTGMC(Preset="medium")
QTGMC(Preset="Slower", InputType=1)
srestore(frate=12)
MCTemporalDenoise(settings="low")
DeRainbow(thresh=10, interlaced=false)
MiniDeen(radius=1, thrY=10, thrUV=12, Y=3, U=3, V=3)
DeHalo_alpha_mt(rx=2.0, ry=2.0, darkstr=1.0, brightstr=1.0, lowsens=50, highsens=50, ss=1.5)
Sharpen(.1)

[Selur]

The source probably is not interlaced, but telecined.
Try filtering on the separated fields.

If Vapoursynth is an option and you have a 'newer' NVIDIA card, something like:

# converting interlaced to half-height progressive for filtering (vsBasicVSRPPFilter) (separate fields (join))
clip = core.std.SeparateFields(clip, tff=False)
# adjusting color space from YUV420P8 to RGBH for vsBasicVSRPPFilter
clip = core.resize.Bicubic(clip=clip, format=vs.RGBH, matrix_in_s="470bg", range_in_s="limited", range_s="full")
# Quality enhancement using BasicVSR++
from vsbasicvsrpp import basicvsrpp as BasicVSRPP
clip = BasicVSRPP(clip, model=4)
# converting half-height progressive (join) to interlaced 
clip = core.std.DoubleWeave(clip=clip, tff=False) # resolution 720x480
clip = core.std.SelectEvery(clip=clip, cycle=2, offsets=[0])
# converting interlaced to half-height progressive for filtering (vsMCTemporalDenoise) (separate fields (join))
clip = core.std.SeparateFields(clip, tff=False)
# adjusting color space from RGBH to YUV444P16 for vsMCTemporalDenoise
clip = core.resize.Bicubic(clip=clip, format=vs.YUV444P16, matrix_s="470bg", range_in_s="full", range_s="limited")
# denoising using MCTemporalDenoise
clip = denoise.MCTemporalDenoise(clip, settings="very high", ncpu=10)
# converting half-height progressive (join) to interlaced 
clip = core.std.DoubleWeave(clip=clip, tff=False) # resolution 720x480
clip = core.std.SelectEvery(clip=clip, cycle=2, offsets=[0])
# Deinterlacing using TIVTC
clip = core.tivtc.TFM(clip)
clip = core.tivtc.TDecimate(clip) # new fps: 23.976

could be used. (note that I did not filter odd and even fields separately, BasicVSR++ is what requires Vapoursynth and the NVIDIA card)
I first separated the fields, applied BasicVSR++ and MCTemporalDenoise, then undid the separation and applied TIVTC.

Cu Selur

Ps.: Dehalo&co would still be needed.