View difference between Paste ID: 0bKmyuV8 and m8X85L1b
SHOW: | | - or go back to the newest paste.
1-
###### srestore v2.79 ######
1+
###### srestore v2.791 ######
2
# ==================================================================================================================================
3
# Change List
4
#
5
# 11/22/2017 AmjadSONY
6-
# - Update v2.78 -> v2.79
6+
# - Update v2.78 -> v2.791
7
#  some edit in speed parameter (in avs26 or avs+ will get different result now)
8
#
9
# 06/15/2017 AmjadSONY
10
# - Update v2.77 -> v2.78
11
#  some fix
12
#
13
# 05/25/2017 AmjadSONY
14
# - Update v2.76 -> v2.77
15
#  fix high bit depth
16
#
17
# 03/31/2017 AmjadSONY
18
# - Update v2.75 -> v2.76
19
#  update to masktools 2.2.5
20
#
21
# 03/09/2017 AmjadSONY
22
# - Update v2.7i -> v2.75
23
#  now work with all high bit depth in avs+
24
#
25
# 02/09/2017 AmjadSONY
26
# - Update v2.7h -> v2.7i
27
#  clean up and add YUY2 and other format
28
#
29
# 03/23/2014 AmjadSONY
30
# - Update v2.7g -> v2.7h
31
#  update to use new Average by tp7
32
#
33
# 11/12/2013 martin53
34
# - Update v2.7f -> v2.7g
35
#  Refactored to use a single line ScriptClip() call like v2.7e - Neccessary to save AviSynth string heap space when RTE is called
36
#
37
# 10/28/2013 martin53
38
# - Update v2.7e -> v2.7f
39
# Callable more than once from the same script with help of unique global variables per instance.
40
# Moderate reformatting to re-establish readability with the new variable names.
41
# Needs GRunT plugin
42
#
43
# ==================================================================================================================================
44
45
function Srestore(clip sourceo, float "frate", "omode", float "speed", float "blocks", int "mode", float "thresh", int "cache", clip "dclip")
46
{
47
avs26=!(VersionNumber() < 2.6)
48
nonyuy2=avs26 && sourceo.isyuy2()
49
source=nonyuy2 ? sourceo.nonyuy2clipin(true) : sourceo
50
ssispmt   = Findstr(VersionString(), "AviSynth+") != 0 && Findstr(VersionString(), "r1576") == 0
51
sislumaonly = ssispmt ? sourceo.isy() : !avs26 ? false : sourceo.isy8()
52
sis444 = ssispmt ? sourceo.is444() : !avs26 ? false : sourceo.isyv24()
53
###### parameters & other necessary vars ######
54
	dm      = default(mode, 2)
55
     sisproslumonl = (!(dm<0) && avs26)
56
    srad  = sisproslumonl || sis444 || sislumaonly ? isfloat(speed) && abs(speed)>=1 ? sqrt(abs(speed)) : 9 : isfloat(speed) && abs(speed)>=1 ? sqrt(abs(speed))*4 : 12
57
    cache = default(cache, -1)
58
    irate = framerate(source)
59
    bsize = isfloat(speed) && speed>0 ? 16 : 32
60
61
	om      = default(omode, 6)
62
	bom     = isstring(om)
63
	thr     = abs(default(thresh, 16))+0.01
64
	frfac   = bom || abs(om-3)<2.5 ? 1 : isfloat(frate) ? frate*5<irate || frate>irate ? 1 :
65
			\ abs(frate)/irate : round(irate*10010)%30000==0 ? 1001./2400. : 480./1001.
66
67
	numr    = abs(frfac*1001-round(frfac*1001))<0.01 ? round(frfac*1001) :
68
			\ abs(1001/frfac-round(1001/frfac))<0.01 ? 1001 : round(frfac*9000)
69
	numr    = isfloat(frate) && abs(irate*numr/float(round(numr/frfac))-frate) >
70
			\ abs(irate*round(frate*100)/float(round(irate*100))-frate) ? round(frate*100) : numr
71
	denm    = round(numr/frfac)
72
73
	###### source preparation & lut ######
74
	out = source
75
	mec = abs(dm)<2 || bom ? nop() : mergeluma(mergechroma(out, out.trim(1, 0), 0.5), out.trim(1, 0), 0.5)
76
	det	= default(dclip, source).killaudio()
77
  	sis420  = ssispmt ? det.is420() : det.isyv12()
78
	det	= dm<0 || sislumaonly ? det : ssispmt ? det.converttoy() : avs26 ? det.converttoy8() : det.converttoyv12()
79
	det	= det.pointresize(sisproslumonl || sis444 || sislumaonly ? int(det.width/(sis420 ? (srad*1.77) : srad)) : srad==4 ? det.width : int(det.width/2/srad+4)*4, sisproslumonl || sis444 || sislumaonly ? int(det.width/(sis420 ? (srad*1.77) : srad)) : srad==4 ? det.height : int(det.height/2/srad+4)*4).trim(2, 0)
80
	det	= !sislumaonly && dm<0 && sis420 ? stackvertical(stackhorizontal(ssispmt ? det.ExtractU() : avs26 ? det.utoy8() : det.utoy(), ssispmt ? det.ExtractV() : avs26 ? det.vtoy8() : det.vtoy()), ssispmt ? det.converttoy() : avs26 ? det.converttoy8() : det) : det
81
	det	= !sislumaonly && dm<0 && !sis420 ? stackhorizontal(stackhorizontal(ssispmt ? det.ExtractU() : avs26 ? det.utoy8() : det.utoy(), ssispmt ? det.ExtractV() : avs26 ? det.vtoy8() : det.vtoy()), ssispmt ? det.converttoy() : avs26 ? det.converttoy8() : det) : det
82
	det = bom ? !avs26 ? det.mt_lut("x 0.5 * 64 +", y=3, u=1, v=1) : det.mt_lut("x 0.5 * 64 scalef +", y=3, u=1, v=1) : det
83
84
		code0   = !avs26 ? "x 128 - y 128 - * 0 > x 128 - abs y 128 - abs < x 128 - 128 x - * y 128 - 128 y - * ? x y + 256 - x y + 256 - * ? 0.25 * 128 +"
85
			  \	 : "x range_half - y range_half - * 0 > x range_half - abs y range_half - abs < x range_half - range_half x - * y range_half - range_half y - * ? x y + range_size - dup * ? 0.25 * range_half +"
86
		code1   = !avs26 ? "x y - 2 ^ 3 * x y + 256 - 2 ^ - 128 +" : "x y - dup * 3 * x y + range_size - 2 ^ - range_half +"
87
		diff 	= mt_makediff(det, det.trim(1, 0), y=3, u=1, v=1)
88
	bclp    = !bom ? mt_lutxy(diff, diff.trim(1, 0), expr=code0, y=3, u=1, v=1).bilinearresize(bsize, bsize) :
89
		       \ mt_lutxy(diff.trim(1, 0), mt_adddiff(diff, diff.trim(2, 0), y=3, u=1, v=1), expr=code1, y=3, u=1, v=1).bilinearresize(bsize, bsize)
90
	dclp    = !avs26 ? diff.trim(1, 0).mt_lut("x 128 - abs 1.1 ^ 1 -", y=3, u=1, v=1).bilinearresize(bsize, bsize)
91
		  \	 : diff.trim(1, 0).mt_lut("x range_half - abs 1.1 ^ 1 scalef -", y=3, u=1, v=1).bilinearresize(bsize, bsize)
92
93
	###### POSTPROCESSING ######
94
    unblend1    = bom ? average(out, 2.0, out.loop(2, 0, 0), -1.0) : nop()
95
    unblend2    = bom ? average(out.trim(2, 0), -1.0, out.trim(1, 0), 2.0) : nop()
96
97
    qmask1  = bom ? mt_makediff(unblend1.removegrain(mode=19, modeU=-1, modeV=-1), unblend1, y=3, u=1, v=1) : nop()
98
    qmask2  = bom ? mt_makediff(unblend2.removegrain(mode=19, modeU=-1, modeV=-1), unblend2, y=3, u=1, v=1) : nop()
99
    diffm   = bom ? mt_makediff(out.loop(2, 0, 0), out, y=3, u=1, v=1).mt_lut("x y - abs", y=3, u=1, v=1).mt_expand() : nop()
100
	# 10/28/2013 note: next UPN expression is buggy, but correct expression not known
101
    bmask	= bom ? !avs26 ? mt_lutxy(qmask1, qmask2, expr="x y - == 128 x 128 - 2 ^ x 128 - 2 ^ y 128 - 2 ^ + / 255 *", y=3, u=1, v=1)
102
			\      : mt_lutxy(qmask1, qmask2, expr="clamp_f_i8 x range_half - dup * dup y range_half - dup * + / range_max *", y=3, u=1, v=1) : nop()
103
    dmask	= bom ? !avs26 ? mt_lutxy(diffm, diffm.trim(2, 0), expr="x 2 * y < x 4 < & 0 y 2 * x < y 4 < & 255 x x y + / 200 * 28 + ? ?", y=3, u=1, v=1)
104
			\      : mt_lutxy(diffm, diffm.trim(2, 0), expr="x 2 * y < x 4 scalef < & 0 y 2 * x < y 4 scalef < & range_max x x y + / 200 scalef * 28 scalef + ? ?", y=3, u=1, v=1) : nop()
105
    pmask	= bom ? !avs26 ? mt_lutxy(dmask, bmask, expr="y 0 > y 255 < & x 0 == x 255 == | & x y ?", y=3, u=1, v=1)
106
			\      : mt_lutxy(dmask, bmask, expr="y 0 > y range_max < & x 0 == x range_max == | & x y ?", y=3, u=1, v=1) : nop()
107
108
    pp0 = bom ? average(out.trim(2, 0), -0.5, out.trim(1, 0), 1.0, out, 1.0, out.loop(2, 0, 0), -0.5) : nop()
109
    pp1 = bom ? mt_merge(unblend1, unblend2, dmask.removegrain(mode=12, modeU=-1, modeV=-1).greyscale(), y=3, u=sislumaonly ? 1 : 3, v=sislumaonly ? 1 : 3) : nop()
110
    pp2 = bom ? mt_merge(unblend1, unblend2, bmask.removegrain(mode=12, modeU=-1, modeV=-1), luma=!sislumaonly) : nop()
111
    pp3 = bom ? mt_merge(unblend1, unblend2, pmask.removegrain(mode=12, modeU=-1, modeV=-1), luma=!sislumaonly).removegrain(mode=0, modeU=12, modeV=12) : nop()
112
113
	fin  = bom ? eval(om) : nop()
114
115
	#######################################
116
    # Unique Identifier Definition
117
    #######################################
118
	Try { global Srestore_InstanceNumber = Srestore_InstanceNumber+1 } catch(err) { global Srestore_InstanceNumber = 1 }
119
120
	# Runtime Environment Preparation & Call
121
	Eval(StrReplace("""	#"
122
        #######################################
123
        # 1) Unique Global Variables Initialization
124
        #######################################
125
		global  Srestore_lfr%%%     = -100
126
		global  Srestore_offs%%%    = 0
127
		global  Srestore_ldet%%%    = -100
128
		global  Srestore_lpos%%%    = 0
129
130
        #######################################
131
        # 2) Unique Runtime Call
132
        #######################################
133
		###### evaluation call & output calculation ######
134
		ARGS = "dm, om, bom, thr, frfac, numr, denm, mec, det, bclp, dclp, fin, ssispmt"
135
		source
136
		GScriptClip("srestore_inside_%%%(last, dm, om, bom, thr, frfac, numr, denm, mec, det, bclp, dclp, fin, ssispmt)", local=true, args=ARGS)
137
138
		###### final decimation & caching ######
139
		nonyuy2 ? nonyuy2clipout(true) : last
140
		cache<0 ? last : last.RequestLinear(8, cache, 5, false, false)
141
142
		temp    = frameratenumerator(det)*float(numr)>2147483600. ? det.assumescaledfps(numr, denm) : nop()
143
		isclip(temp) ? last.changefps(temp, linear=true) : last.changefps(frameratenumerator(source)*numr, frameratedenominator(source)*denm, linear=true)
144
145
		return last
146
147
        #######################################
148
        # 3) Unique Runtime Function Definition
149
        #######################################
150
		function srestore_inside_%%%(clip source, dm, om, bom, thr, frfac, numr, denm, mec, det, bclp, dclp, fin, ssispmt) {
151
			ssIsVF  = ssispmt ? source.IsVideoFloat() : false
152
			shbdbpc = ssispmt ? source.BitsPerComponent() : 8
153
			sIshbd  = shbdbpc > 8
154
			sdoscal = !ssIsVF && sIshbd
155
			sscalva = sdoscal ? shbdbpc == 10 ? 4 : shbdbpc == 12 ? 8 : shbdbpc == 14 ? 64 : 256 : 1
156
				ssIsVF ? eval("
157
			v015625 = .0000612745098039
158
			vp125   = .0004901960784314
159
			vp01    = .0000392156862745
160
			vp45    = .0017647058823529
161
			vp55    = .002156862745098
162
			vp51    = .002
163
			vp35    = .0013725490196078
164
			v0p5    = 0.0019607843137255
165
			v1      = 0.003921568627451
166
			v125    = 0.4901960784313725
167
						") : \
168
							eval("
169
			v015625 = .015625
170
			vp125   = .125
171
			vp01    = .01
172
			vp45    = .45
173
			vp55    = .55
174
			vp51    = .51
175
			vp35    = .35
176
			v0p5    = 0.5
177
			v1      = 1
178
			v125    = sdoscal ? 125*sscalva : 125
179
								")
180
			svhalf  = ssIsVF ? 0.5 : sdoscal ? shbdbpc == 10 ? 512 : shbdbpc == 12 ? 1024 : shbdbpc == 14 ? 8192 : 32768 : 128
181
			thr     = ssIsVF ? thr/255 : sdoscal ? thr*sscalva : thr
182
183
			### preparation ###
184
			cfr = current_frame
185
			jmp = Srestore_lfr%%%+1==cfr
186
			cfo = (((cfr%denm)*numr*2+denm+(versionnumber()>2.575 ? numr : 0))%(2*denm))-denm
187
			bfo = cfo>-numr && cfo<=numr
188
189
			global	Srestore_lfr%%%     = cfr
190
			global	Srestore_offs%%%    = bfo && Srestore_offs%%%<=-4*numr
191
												\? Srestore_offs%%%+2*denm
192
												\: bfo && Srestore_offs%%%>=4*numr
193
													\? Srestore_offs%%%-2*denm
194
													\: Srestore_offs%%%
195
196
			pos = frfac==1 ? 0 : bfo ? -round((cfo+Srestore_offs%%%)/(2.*numr)) : Srestore_lpos%%%
197
			cof = cfo+Srestore_offs%%%+2*numr*pos
198
199
			global	Srestore_ldet%%%    = cfr+pos==Srestore_ldet%%% ? -1 : cfr+pos
200
201
			## diff value shifting ##
202
			d_v = yplanemax(dclp)+v015625
203
			Srestore_d43%%% = jmp ? Srestore_d32%%% : d_v
204
			global 	Srestore_d32%%% = jmp ? Srestore_d21%%% : d_v
205
			global	Srestore_d21%%% = jmp ? Srestore_d10%%% : d_v
206
			global	Srestore_d10%%% = jmp ? Srestore_d01%%% : d_v
207
			global	Srestore_d01%%% = jmp ? Srestore_d12%%% : d_v
208
			global	Srestore_d12%%% = jmp ? Srestore_d23%%% : d_v
209
			global	Srestore_d23%%% = jmp ? Srestore_d34%%% : d_v
210
			global  Srestore_d34%%% = d_v
211
212
			## diff value shifting ##
213
			m_v = isfloat(om) && abs(om)>5 ? lumadifference(det, det.trim(2, 0))+v015625 : v1
214
			Srestore_m53%%% = jmp ? Srestore_m42%%% : m_v
215
			global 	Srestore_m42%%% = jmp ? Srestore_m31%%% : m_v
216
			global  Srestore_m31%%% = jmp ? Srestore_m20%%% : m_v
217
			global	Srestore_m20%%% = jmp ? Srestore_m11%%% : m_v
218
			global	Srestore_m11%%% = jmp ? Srestore_m02%%% : m_v
219
			global	Srestore_m02%%% = jmp ? Srestore_m13%%% : m_v
220
			global	Srestore_m13%%% = jmp ? Srestore_m24%%% : m_v
221
			global  Srestore_m24%%% = m_v
222
223
			## get blend and clear values ##
224
			b_v = svhalf-yplanemin(bclp)
225
			b_v = b_v<v1 ? vp125 : b_v
226
			c_v = yplanemax(bclp)-svhalf
227
			c_v = c_v<v1 ? vp125 : c_v
228
229
			## blend value shifting ##
230
			Srestore_bp3%%% = jmp ? Srestore_bp2%%% : bom ? b_v-c_v : b_v
231
			global 	Srestore_bp2%%% = jmp ? Srestore_bp1%%% : Srestore_bp3%%%
232
			global	Srestore_bp1%%% = jmp ? Srestore_bn0%%% : Srestore_bp3%%%
233
			global	Srestore_bn0%%% = jmp ? Srestore_bn1%%% : Srestore_bp3%%%
234
			global	Srestore_bn1%%% = jmp ? Srestore_bn2%%% : Srestore_bp3%%%
235
			global	Srestore_bn2%%% = jmp ? Srestore_bn3%%% : Srestore_bp3%%%
236
			global	Srestore_bn3%%% = bom ? b_v-c_v : b_v
237
238
			## clear value shifting ##
239
			Srestore_cp3%%% = jmp ? Srestore_cp2%%% : c_v
240
			global	Srestore_cp2%%% = jmp ? Srestore_cp1%%% : c_v
241
			global	Srestore_cp1%%% = jmp ? Srestore_cn0%%% : c_v
242
			global	Srestore_cn0%%% = jmp ? Srestore_cn1%%% : c_v
243
			global	Srestore_cn1%%% = jmp ? Srestore_cn2%%% : c_v
244
			global	Srestore_cn2%%% = jmp ? Srestore_cn3%%% : c_v
245
			global	Srestore_cn3%%% = c_v
246
247
			## used detection values ##
248
			bb  = select(pos+2, Srestore_bp3%%%, Srestore_bp2%%%, Srestore_bp1%%%, Srestore_bn0%%%, Srestore_bn1%%%)
249
			bc  = select(pos+2, Srestore_bp2%%%, Srestore_bp1%%%, Srestore_bn0%%%, Srestore_bn1%%%, Srestore_bn2%%%)
250
			bn  = select(pos+2, Srestore_bp1%%%, Srestore_bn0%%%, Srestore_bn1%%%, Srestore_bn2%%%, Srestore_bn3%%%)
251
252
			cb  = select(pos+2, Srestore_cp3%%%, Srestore_cp2%%%, Srestore_cp1%%%, Srestore_cn0%%%, Srestore_cn1%%%)
253
			cc  = select(pos+2, Srestore_cp2%%%, Srestore_cp1%%%, Srestore_cn0%%%, Srestore_cn1%%%, Srestore_cn2%%%)
254
			cn  = select(pos+2, Srestore_cp1%%%, Srestore_cn0%%%, Srestore_cn1%%%, Srestore_cn2%%%, Srestore_cn3%%%)
255
256
			dbb = select(pos+2, Srestore_d43%%%, Srestore_d32%%%, Srestore_d21%%%, Srestore_d10%%%, Srestore_d01%%%)
257
			dbc = select(pos+2, Srestore_d32%%%, Srestore_d21%%%, Srestore_d10%%%, Srestore_d01%%%, Srestore_d12%%%)
258
			dcn = select(pos+2, Srestore_d21%%%, Srestore_d10%%%, Srestore_d01%%%, Srestore_d12%%%, Srestore_d23%%%)
259
			dnn = select(pos+2, Srestore_d10%%%, Srestore_d01%%%, Srestore_d12%%%, Srestore_d23%%%, Srestore_d34%%%)
260
			dn2 = select(pos+2, Srestore_d01%%%, Srestore_d12%%%, Srestore_d23%%%, Srestore_d34%%%, Srestore_d34%%%)
261
262
			mb1 = select(pos+2, Srestore_m53%%%, Srestore_m42%%%, Srestore_m31%%%, Srestore_m20%%%, Srestore_m11%%%)
263
			mb  = select(pos+2, Srestore_m42%%%, Srestore_m31%%%, Srestore_m20%%%, Srestore_m11%%%, Srestore_m02%%%)
264
			mc  = select(pos+2, Srestore_m31%%%, Srestore_m20%%%, Srestore_m11%%%, Srestore_m02%%%, Srestore_m13%%%)
265
			mn  = select(pos+2, Srestore_m20%%%, Srestore_m11%%%, Srestore_m02%%%, Srestore_m13%%%, Srestore_m24%%%)
266
			mn1 = select(pos+2, Srestore_m11%%%, Srestore_m02%%%, Srestore_m13%%%, Srestore_m24%%%, vp01)
267
268
269
			### basic calculation ###
270
			bbool   = .8*bc*cb>bb*cc && .8*bc*cn>bn*cc && bc*bc>cc
271
			blend   = bbool && bc*5>cc && dbc+dcn>1.5*thr && (dbb<7*dbc || dbb<8*dcn) && (dnn<8*dcn || dnn<7*dbc) &&
272
					\ (mb<mb1 && mb<mc || mn<mn1 && mn<mc || (dbb+dnn)*4<dbc+dcn || (bb*cc*5<bc*cb || mb>thr) && (bn*cc*5<bc*cn || mn>thr) && bc>thr)
273
274
			clear   = dbb+dbc>thr && dcn+dnn>thr && (bc<2*bb || bc<2*bn) && (dbb+dnn)*2>dbc+dcn && (mc<.96*mb && mc<.96*mn && (bb*2>cb || bn*2>cn) &&
275
					\ cc>cb && cc>cn || frfac>vp45 && frfac<vp55 && .8*mc>mb1 && .8*mc>mn1 && mb>.8*mn && mn>.8*mb)
276
277
			highd   = dcn>5*dbc && dcn>5*dnn && dcn>thr && dbc<thr && dnn<thr
278
279
			lowd    = dcn*5<dbc && dcn*5<dnn && dbc>thr && dnn>thr && dcn<thr && frfac>vp35 && (frfac<vp51 || dcn*5<dbb)
280
281
			res     = (Srestore_d43%%%<thr
282
					  \&& Srestore_d32%%%<thr
283
					  \&& Srestore_d21%%%<thr
284
					  \&& Srestore_d10%%%<thr
285
					  \&& Srestore_d01%%%<thr
286
					  \&& Srestore_d12%%%<thr
287
					  \&& Srestore_d23%%%<thr
288
					  \&& Srestore_d34%%%<thr)
289
					  \|| (dbc*4<dbb && dcn*4<dbb && dnn*4<dbb && dn2*4<dbb)
290
					  \|| (dcn*4<dbc && dnn*4<dbc && dn2*4<dbc)
291
292
293
			### offset calculation ###
294
			odm     = blend ? denm : clear ? 0 : highd ? denm-numr : lowd ? 2*denm-numr : cof
295
			odm     = odm+round((cof-odm)/(2.*denm))*2*denm
296
			odr     = (blend ? denm-numr : clear ? numr : highd ? numr : frfac<v0p5 ? 2*numr : 2*(denm-numr)) * .9
297
298
			cof     = Srestore_ldet%%%<0
299
						\? cof
300
						\: cof>odm+odr
301
							\? cof-Srestore_offs%%%-odm-odr>denm && res
302
								\? odm+2*denm-odr
303
								\: odm+odr
304
							\: cof<odm-odr
305
								\? Srestore_offs%%%>denm && res
306
									\? odm-2*denm+odr
307
									\: odm-odr
308
								\: Srestore_offs%%%<-1.15*denm && res
309
									\? cof+2*denm
310
									\: Srestore_offs%%%>1.25*denm && res
311
										\? cof-2*denm
312
										\: cof
313
314
			global	Srestore_offs%%%    = frfac==1 ? 0 : cof-cfo-2*numr*pos
315
			global	Srestore_lpos%%%    = pos
316
			opos    = frfac==1 ? 0 : -round((cfo+Srestore_offs%%%+(bfo && Srestore_offs%%%<=-4*numr ? denm : 0))/(2.*numr))
317
			pos     = opos<-2 ? -2 : opos>2 ? 2 : opos
318
319
			### frame output calculation - resync - dup ###
320
			dbb = select(pos+2, Srestore_d43%%%, Srestore_d32%%%, Srestore_d21%%%, Srestore_d10%%%, Srestore_d01%%%)
321
			dbc = select(pos+2, Srestore_d32%%%, Srestore_d21%%%, Srestore_d10%%%, Srestore_d01%%%, Srestore_d12%%%)
322
			dcn = select(pos+2, Srestore_d21%%%, Srestore_d10%%%, Srestore_d01%%%, Srestore_d12%%%, Srestore_d23%%%)
323
			dnn = select(pos+2, Srestore_d10%%%, Srestore_d01%%%, Srestore_d12%%%, Srestore_d23%%%, Srestore_d34%%%)
324
325
			## dup_hq - merge ##
326
			dup     = opos!=pos || abs(dm)<2 || abs(dm)==3 ? 0 : dcn*5<dbc && dnn*5<dbc && (dcn<1.25*thr || bn<bc && pos==Srestore_lpos%%%) ||
327
					\ (dcn*dcn<dbc || dcn*5<dbc) && bn<bc && pos==Srestore_lpos%%% && dnn<.9*dbc || dnn*9<dbc && dcn*3<dbc ? 1 :
328
					\ (dbc*dbc<dcn || dbc*5<dcn) && bb<bc && pos==Srestore_lpos%%% && dbb<.9*dcn || dbb*9<dcn && dbc*3<dcn ||
329
					\ dbb*5<dcn && dbc*5<dcn && (dbc<1.25*thr || bb<bc && pos==Srestore_lpos%%%) ? -1 : 0
330
			mer     = !bom && opos==pos && dup==0 && abs(dm)>2 && (dbc*8<dcn || dbc*8<dbb || dcn*8<dbc || dcn*8<dnn || dbc*2<thr ||
331
					\ dcn*2<thr || dnn*9<dbc && dcn*3<dbc || dbb*9<dcn && dbc*3<dcn)
332
333
			## deblend - doubleblend removal - postprocessing ##
334
			add     = Srestore_bp1%%%*Srestore_cn2%%%>Srestore_bn2%%%*Srestore_cp1%%%*(v1+thr*.01)
335
					  \&& Srestore_bn0%%%*Srestore_cn2%%%>Srestore_bn2%%%*Srestore_cn0%%%*(v1+thr*.01)
336
					  \&& Srestore_cn2%%%*Srestore_bn1%%%>Srestore_cn1%%%*Srestore_bn2%%%*(v1+thr*.01)
337
			dup     = bom
338
						\? (Srestore_bn0%%%>Srestore_bp2%%%
339
							\&& Srestore_bn0%%%>=Srestore_bp1%%%
340
							\&& Srestore_bn0%%%>Srestore_bn1%%%
341
							\&& Srestore_bn0%%%>Srestore_bn2%%%
342
							\&& Srestore_cn0%%%<v125
343
							\? (Srestore_d12%%%*Srestore_d12%%%<Srestore_d10%%%
344
								\|| Srestore_d12%%%*9<Srestore_d10%%%)
345
								\? 1
346
								\: (Srestore_d10%%%*Srestore_d10%%%<Srestore_d12%%%
347
									\|| Srestore_d10%%%*9<Srestore_d12%%%)
348
									\? 0
349
									\: 4
350
							\: Srestore_bp1%%%>Srestore_bp3%%%
351
								\&& Srestore_bp1%%%>=Srestore_bp2%%%
352
								\&& Srestore_bp1%%%>Srestore_bn0%%%
353
								\&& Srestore_bp1%%%>Srestore_bn1%%%
354
								\? 1
355
								\: 0)
356
						\: dup!=0
357
							\? dup
358
							\: om>0 && om<5
359
								\? (bbool==false
360
									\? 0
361
									\: (om==4 && Srestore_bp1%%%*Srestore_cn1%%%<Srestore_bn1%%%*Srestore_cp1%%%)
362
										\|| (om==3 && Srestore_d10%%%<Srestore_d01%%%)
363
										\|| om==1
364
										\? -1
365
										\: 1)
366
								\: om==5
367
									\? (Srestore_bp1%%%*Srestore_cp2%%%>Srestore_bp2%%%*Srestore_cp1%%%*(v1+thr*.01)
368
										\&& Srestore_bn0%%%*Srestore_cp2%%%>Srestore_bp2%%%*Srestore_cn0%%%*(v1+thr*.01)
369
										\&& Srestore_cp2%%%*Srestore_bn1%%%>Srestore_cn1%%%*Srestore_bp2%%%*(v1+thr*.01)
370
										\&& (add==false || Srestore_cp2%%%*Srestore_bn2%%%>Srestore_cn2%%%*Srestore_bp2%%%)
371
										\? -2
372
										\: add
373
											\? 2
374
											\: Srestore_bn0%%%*Srestore_cp1%%%>Srestore_bp1%%%*Srestore_cn0%%%
375
												\&& (Srestore_bn0%%%*Srestore_cn1%%%<Srestore_bn1%%%*Srestore_cn0%%%
376
												\|| Srestore_cp1%%%*Srestore_bn1%%%>Srestore_cn1%%%*Srestore_bp1%%%)
377
												\? -1
378
												\: Srestore_bn0%%%*Srestore_cn1%%%>Srestore_bn1%%%*Srestore_cn0%%%
379
													\? 1
380
													\: 0)
381
									\: 0
382
383
			### output clip ###
384
			oclp    = mer && dup==0 ? mec : source
385
			opos    = opos + dup - (dup==0 && mer && dbc<dcn ? 1 : 0)
386
			dup==4 ? fin : opos<0 ? oclp.loop(1-opos, 0, 0) : oclp.trim(opos, 0)
387
388
		}
389
	""", "%%%", string(Srestore_InstanceNumber)))	#"
390
}
391
392
393
394
# Original:- By Vampiredom, Gavino, IanB
395
# Repeated, string replacements
396
397
Function StrReplace(string s, string find, string replace)
398
{ i = s.FindStr(find) return (i==0 ? s : (s.LeftStr(i-1)+replace+s.MidStr(Strlen(find)+i)).StrReplace(find, replace)) }