setmemorymax(512)
setmtmode(5,6)
ImageSource("raw.png",start=0,end=0)
RatioResize(320.0,"adjust2w",kernel="point")
r = ShowRed("YV12")
g = ShowGreen("YV12")
b = ShowBlue("YV12")
interleave(r,g,b)
raw=last
ImageSource("Dogway approach.png",start=0,end=0)
RatioResize(320.0,"adjust2w",kernel="point")
r = ShowRed("YV12")
g = ShowGreen("YV12")
b = ShowBlue("YV12")
interleave(r,g,b)
blurred=last
ImageSource("gdapt.png",start=0,end=0) # using settings as 0:2:1:0
RatioResize(320.0,"adjust2w",kernel="point")
r = ShowRed("YV12")
g = ShowGreen("YV12")
b = ShowBlue("YV12")
interleave(r,g,b)
gdapt=last
setmtmode(2,6)

# get difference of both approaches
mt_lutxy(blurred,gdapt," x y - abs",u=-128,v=-128)
# ensure we always get the highest change regardless of plane
mt_logic(mt_logic(SelectEvery (3, 0), SelectEvery (3, 1),"max"),SelectEvery (3, 2),"max")
# cleaning 1px dots
removegrain(1)
# strengthen the mask
mt_binarize(25,u=-128,v=-128) # sames as mt_lut(yexpr=" x 25 > 255 0 ?")

# mix my original approach with source using the diff of gdapt/blurred as a mask
mt_merge(raw,blurred,last.invert,luma=true)

mergergb(SelectEvery (3, 0), SelectEvery (3, 1), SelectEvery (3, 2))