Восстановление прогрессивного потока с NTSC-DVD, транскодированного из PAL

#########################
###
### DoubleFPS()
###
#########################
###
### v.0.3 fixed typo
### changed defaults
### improved accuracy
###
### Function for doubling framerate from original FPS through frame interpolation.
### Useful for better motion smoothness or slow motion. It is prone to artefacting
### so is probably better to fine tune settings apart from function parameters and
### do some craft work.
###
### Didйe function adapted by Dogway (23/02/2011) from following documentation:
### http://forum.doom9.org/showpost.php?p=1409324&postcount=18
### http://www.avsforum.com/avs-vb/showthread.php?p=20045304
###
### It is recommended to supply a clean source to feed this function. It is slow.
###
###
### mode [default: 2]
### -------------------
### Interpolation method:
### 1 = MBlockFps (block interpolation)
### 2 = MFlowFps (pixel interpolation)
### You might want to play with ThSCD.
###
### mopro [default: false]
### -------------------
### motion protection; it duplicates the last frame when there are big
### changes (prone to artifacting), useful for scenechanges, motion blur...
### you may not want it as it might create stutter. [true, false]
###
###
### bias [default: 128]
### -------------------
### Threshold for frame difference detection.
### You wouldn't want to set it higher from
### default value. [0...255]
###
###
### show [default: false]
### -------------------
### Shows a "dot" in the interpolated frames. [true, false]
###
###_________________________________
### |
### Example: |
### |
### DoubleFPS(bias=32,show=true) |
###_________________________________|
###
###
### FEEL FREE TO REMOVE THIS USELESS HELP
##############################################################
function DoubleFPS(clip clp, int "numerator", int "denominator", int "pel",\
int "mode", bool "mopro", int "bias", bool "show")
{
# Defaults
#
numerator = Default(numerator, 2)
denominator = Default(denominator, 1)
pel = Default(pel, 2)
mode = Default(mode, 2)
mopro = Default(mopro, true)
bias = String(Default(bias, 128))
showdot = Default(show, false)
# Here is code of function...
#
v = addborders(clp, 16, 16, 16, 16)
o = v.assumefps(FramerateNumerator(clp))
ox = v.width()
oy = v.height()
S = (showdot == true) ? v.subtitle(".", x=32, y=24).MSuper(pel=pel, hpad=16, vpad=16, rfilter=4) :
\ v.MSuper(pel=pel, hpad=16, vpad=16, rfilter=4)
A = S.MAnalyse(isb=true, blksize=16, overlap=4, plevel=0, search=3, searchparam=8, divide = 2, badrange=(-24), dct=0)
B = S.MAnalyse(isb=false, blksize=16, overlap=4, plevel=0, search=3, searchparam=8, divide = 2, badrange=(-24), dct=0)
C = S.MRecalculate(A, blksize=8, overlap=2, search=3, searchparam=4, thSAD=80)
D = S.MRecalculate(B, blksize=8, overlap=2, search=3, searchparam=4, thSAD=80)
E = S.MRecalculate(C, blksize=4, overlap=0, search=3, searchparam=1, thSAD=100)
F = S.MRecalculate(D, blksize=4, overlap=0, search=3, searchparam=1, thSAD=100)
double = (mode == 1) ? v.MBlockFps(S, E, F, num = FramerateNumerator(v) * numerator, den = FramerateDenominator(v) * denominator, mode = 0, blend = false, ThSCD1 = 220)
\ : (mode == 2) ? v.MFlowFps(S, E, F, num = FramerateNumerator(v) * numerator, den = FramerateDenominator(v) * denominator, mask = 2, ml = 180, blend = false)
\ : NOP()
diff2next = mt_makediff(o, o.selectevery(1, 1)).mt_lut("x "+bias+" - abs 32 / 1 2.0 / ^ 128 *", U=-128, V=-128)
diff2next = mt_lutf(diff2next, diff2next, yexpr="x", mode="average").pointresize(32, 32)
diff2next = interleave(diff2next.selectevery(4, 0).tsg(2), diff2next.selectevery(4, 1).tsg(2),
\ diff2next.selectevery(4, 2).tsg(2), diff2next.selectevery(4, 3).tsg(2))
max = diff2next.mt_logic(diff2next.selectevery(1, -3), "max")
\ .mt_logic(diff2next.selectevery(1, -2), "max")
\ .mt_logic(diff2next.selectevery(1, -1), "max")
\ .mt_logic(diff2next.selectevery(1, 1), "max")
\ .mt_logic(diff2next.selectevery(1, 2), "max")
\ .mt_logic(diff2next.selectevery(1, 3), "max")
ismax = mt_lutxy(diff2next, max, "x y < 255 0 ?", U=-128, V=-128).pointresize(ox, oy) # swap 255 0 to interpolate only on big changes.
themask = interleave(o.mt_lut("0"), ismax)
interleave(o, o).mt_merge(double, themask, luma=true, U=3, V=3)
(mopro == true) ? last.crop(16, 16, -16, -16) : double.crop(16, 16, -16, -16)
return (last)
}
function tsg(clip clp, int t)
{
last = clp
(t < 5) ? last : last.temporalsoften(1, 255, 0, 255, 2).merge(last, 0.25)
(t < 4) ? last : last.temporalsoften(1, 255, 0, 255, 2).merge(last, 0.25)
(t < 3) ? last : last.temporalsoften(1, 255, 0, 255, 2).merge(last, 0.25)
(t < 2) ? last : last.temporalsoften(1, 255, 0, 255, 2).merge(last, 0.25)
(t < 1) ? last : last.temporalsoften(1, 255, 0, 255, 2).merge(last, 0.25)
return (last)
}

# Frame Rate Converter
# Version: 25-Apr-2017
# By Etienne Charland
# Based on Oleg Yushko's YFRC artifact masking,
# johnmeyer's frame interpolation code, and
# raffriff42's "weak mask" and output options.
# Slightly modified interface by manolito
#
# This program is free software; you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation; either version 2 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program; if not, write to the Free Software
# Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA, or visit
# http:#www.gnu.org/copyleft/gpl.html.
#######################################################################################
### mx_fps
### Changes the frame rate with interpolation and strong artifact removal.
##
## YV12/YV24/Y8/YUY2
## Requires: MaskTools2, MvTools2, RemoveGrain and fftw3.dll (in the System32 or SysWOW64 folder)
##
## @ fps - The new framerate
##
## @ Preset - The speed/quality preset [slow|normal|fast|faster]. (default=normal)
##
## @ BlkSize - The horizontal block size (default = Width>=1200||Height>=900 ? 32 : Width>=720||C.Height>=480 ? 16 : 8)
##
## @ BlkSizeV - The vertical block size (default = BlkSize)
##
## @ Output - Output mode [auto|inter|none|mask|over] (default=auto)
## auto=normal artifact masking; flow=interpolation only
## none=ConvertFPS only; mask=mask only; over=mask as cyan overlay for debugging
##
## @ MaskGam - A gamma to be applied to the raw mask, between 0 and 1.
## Higher value means stronger mask. (Default=.6)
##
## @ MaskTrh - The treshold where a block is considered bad, between 0 and 255. 0 to disable artifact masking.
## Higher value means stronger mask. (Default=90)
##
## @ SkipOver - Skip interpolation of frames when artifacts cover more than specified treshold,
## 0 to disable. Does not work in Avisynth+ MT mode. (Default=20)
##
function mx_fps(clip C, float "fps", string "Preset", int "BlkSize")
{
fps = default(fps, 25.000)
NewNum = int(fps * 1000)
NewDen = 1000
Output = "auto"
Preset = Default(Preset, "normal")
BlkSize = Default(BlkSize, C.Width>=1200||C.Height>=900 ? 32 : C.Width>=720||C.Height>=480 ? 16 : 8)
BlkSizeV = BlkSize
blkmin = BlkSize > BlkSizeV ? BlkSizeV : BlkSize
MaskGam = .6
MaskTrh = 90
SkipOver = 20
Assert(Preset == "slow" || Preset == "normal" || Preset == "fast" || Preset == "faster", "mx_fps: Preset must be slow, normal, fast or faster")
Assert(BlkSize == 8 || BlkSize == 16 || BlkSize == 32, "mx_fps: BlkSize must be 8, 16 or 32")
# Performance settings: slow, normal, fast, faster
Recalculate= preset == "slow" || preset == "normal"
Prefilter = preset != "faster"
DCT = preset == "slow" ? 1 : 0
## "B" - Blending, "BHard" - No blending
BHard = C.ChangeFPS(NewNum, NewDen)
B = C.ConvertFPS(NewNum, NewDen)
## jm_fps interpolation
prefiltered = Prefilter ? RemoveGrain(C, 22) : C
superfilt = MSuper(prefiltered, hpad = 16, vpad = 16) # all levels for MAnalyse
super = Prefilter ? MSuper(C, hpad = 16, vpad = 16, levels = 1) : superfilt # one level is enough for MRecalculate
bak = MAnalyse(superfilt, isb = true, blksize = BlkSize, blksizeV = BlkSizeV, overlap = blkmin>8?4:blkmin>4?2:0, search = 3, dct = DCT)
fwd = MAnalyse(superfilt, isb = false, blksize = BlkSize, blksizeV = BlkSizeV, overlap = blkmin>8?4:blkmin>4?2:0, search = 3, dct = DCT)
fwd = Recalculate ? MRecalculate(super, fwd, blksize = BlkSize/2, blksizeV = BlkSizeV/2, overlap = blkmin>8?2:0, thSAD = 100) : fwd
bak = Recalculate ? MRecalculate(super, bak, blksize = BlkSize/2, blksizeV = BlkSizeV/2, overlap = blkmin>8?2:0, thSAD = 100) : bak
Flow = MFlowFps(C, super, bak, fwd, num = NewNum, den = NewDen, blend = false, ml = 200, mask = 2)
## "EM" - error or artifact mask
EM = MaskTrh > 0 ? C.MMask(bak, ml=255, kind=1, ysc=255, gamma=1.0/MaskGam) : BlankClip(C) # kind=SAD
EMfwd = MaskTrh > 0 ? C.MMask(fwd, ml=255, kind=1, gamma=1.0/MaskGam) : EM # kind=temporal blending
EM = MaskTrh > 0 ? EM.Overlay(EMfwd, opacity=0.5, mode="lighten") : EM
SkipEM = EM
BlkSize = C.Width>1200||C.Height>900 ? BlkSize : BlkSize == 32 ? 16 : BlkSize
BlkSizeV = C.Width>1200||C.Height>900 ? BlkSizeV : BlkSizeV == 32 ? 16 : BlkSizeV
EM = EM.BicubicResize(Round(C.Width/BlkSize/4.0)*4, Round(C.Height/BlkSizeV/4.0)*4)
\ .mt_expand(mode= mt_circle(zero=true, radius=1))
\ .mt_binarize(255-MaskTrh)
\ .Blur(.6)
\ .BicubicResize(C.Width, C.Height)
## "Sc" - scene detection / SkipOver
Sc = SkipOver > 0 ? ConditionalFilter(SkipEM, BlankClip(EM, color=color_white), BlankClip(EM),
\ "AverageLuma()", ">", string(SkipOver)) : BlankClip(EM)
## Convert masks to desired frame rate
EM = EM.ChangeFPS(NewNum, NewDen)
Sc = Sc.ChangeFPS(NewNum, NewDen)
## the mx_fps magic happens
M = mt_merge(Flow, B, EM, luma=true)
M = SkipOver > 0 ? mt_merge(M, BHard, Sc, luma=true) : M
R = (StrCmpi(Output, "auto")==0) [** auto: artifact masking *]
\ ? M
\ : (StrCmpi(Output, "flow")==0) [** flow: interpolation only *]
\ ? Flow
\ : (StrCmpi(Output, "none")==0) [** none: ConvertFPS only *]
\ ? B
\ : (StrCmpi(Output, "mask")==0) [** mask: mask only *]
\ ? EM
\ : (StrCmpi(Output, "over")==0) [** over: mask as cyan overlay *]
\ ? mt_merge(
\ Flow.Overlay(MergeRGB(BlankClip(EM), EM, EM), mode="Add", opacity=0.40, pc_range=true),
\ BlankClip(EM, color=color_darkgoldenrod), Sc.mt_lut("x 2 / "), luma=true)
\ : NOP()
return R
}

function ConvertMFlowFPS(clip input, int "numerator", int "denominator")
{
prefiltered = input.blur(0.6)
super = MSuper(input, hpad = 16, vpad = 16, levels = 1) # one level is enough for MRecalculate
superfilt = MSuper(prefiltered, hpad = 16, vpad = 16) # all levels for MAnalyse
A = superfilt.MAnalyse(isb=true, blksize=16, overlap=4, plevel=0, search=3, searchparam=2, badrange=(-24), dct=0)
B = superfilt.MAnalyse(isb=false, blksize=16, overlap=4, plevel=0, search=3, searchparam=2, badrange=(-24), dct=0)
C = superfilt.MRecalculate(A, blksize=8, overlap=2, search=3, searchparam=1, thSAD=100)
D = superfilt.MRecalculate(B, blksize=8, overlap=2, search=3, searchparam=1, thSAD=100)
out = MFlowFps(input, super, C, D, num = numerator, den = denominator, blend = false, ml = 200, mask = 2)
return out
}