SHOW:
|
|
- or go back to the newest paste.
1 | ################################################################################################ | |
2 | ### ### | |
3 | ### Simple MDegrain Mod - SMDegrain() ### | |
4 | ### ### | |
5 | ### Mod by Dogway - Original idea by Caroliano ### | |
6 | ### ### | |
7 | ### Special Thanks: Sagekilla, Didée, cretindesalpes, Gavino and MVtools people ### | |
8 | ### ### | |
9 | ### v3.1.2d (Dogway's mod) - 21 July 2015 ### | |
10 | - | ### mod in 2016-01-24 to fix subpixel=3 and add Prefilter auto 16 (lsb) support ### |
10 | + | ### mod in 2016-06-23 to fix subpixel=3 and add Prefilter auto 16 (lsb) support ### |
11 | ### and yuy2 support in avs 2.6 with masktools 2.6 ### | |
12 | ### and medianblur2 now will be used in Minblur instead of Quantile RemoveGrainHD in avs 2.6 ### | |
13 | ### ### | |
14 | ################################################################################################ | |
15 | ### | |
16 | ### General purpose simple degrain function. Pure temporal denoiser. Basically a wrapper(function)/frontend of mvtools2+mdegrain | |
17 | ### with some added common related options. Goal is accessibility and quality but not targeted to any specific kind of source. | |
18 | ### The reason behind is to keep it simple so aside masktools2 you will only need MVTools2. | |
19 | ### | |
20 | ### Check documentation for deep explanation on settings and defaults. | |
21 | ### VideoHelp thread: (http://forum.videohelp.com/threads/369142) | |
22 | ### | |
23 | ################################################################################################ | |
24 | ||
25 | function SMDegrain (clip input , int "tr" , int "thSAD" , int "thSADC" , bool "RefineMotion", val "contrasharp", clip "CClip" , bool "interlaced" , int "plane" , int "Globals", \ | |
26 | int "pel" , int "subpixel", val "prefilter", clip "mfilter", int "blksize" , int "overlap" , int "search" , bool "truemotion" , int "limit" , int "limitc" , int "thSCD1" , int "thSCD2", \ | |
27 | bool "chroma", int "hpad" , int "vpad" , bool "lsb" , bool "lsb_in" , bool "lsb_out" , int "mode" , bool "slices" , val "Show" , float "Str" , float "Amp"){ | |
28 | ||
29 | ||
30 | ||
31 | # Defaults & Conditionals | |
32 | ||
33 | ||
34 | tr = default( tr, 2 ) | |
35 | thSAD = default( thSAD, 300 ) | |
36 | thSAD2=int(thSAD/2) | |
37 | thSADC = default( thSADC, thSAD2) | |
38 | RefineMotion = default( RefineMotion,false) | |
39 | ||
40 | lsb_in = default( lsb_in, false) | |
41 | lsb_out = default( lsb_out, false) | |
42 | lsb = default( lsb , lsb_in || lsb_out) | |
43 | mode = default( mode, 0 ) | |
44 | mode = mode != 9 ? mode : -1 | |
45 | ampn = mode != 9 ? 0.0 : 1.0 | |
46 | slices = default(slices, true) | |
47 | ||
48 | Globals = default( Globals, 0 ) | |
49 | GlobalR = Globals == 1 ? true : false | |
50 | GlobalO = Globals == 3 ? true : false | |
51 | if1 = Defined(CClip) ? true : false | |
52 | ||
53 | Contrasharp = default( Contrasharp, !GlobalO && if1 ? true : false) | |
54 | Interlaced = default( Interlaced, false) | |
55 | ||
56 | w = input.width () | |
57 | h = input.height() | |
58 | sShow = IsString(show) ? true : false | |
59 | preclip = IsClip(prefilter) ? true : false | |
60 | planar = IsYUY2(input) ? true : false | |
61 | ifC = IsBool(Contrasharp) ? true : false | |
62 | if0 = ifC ? (Contrasharp ? true : false) : (Contrasharp > 0.0 ? true : false) | |
63 | if2 = if0 && lsb && !GlobalO ? true : false | |
64 | if3 = tr > 3 ? true : false | |
65 | if4 = (w > 1099 || h > (lsb_in ? 1199 : 599)) ? true : false | |
66 | if5 = interlaced ? (GetParity(input) ? true : false) : nop() | |
67 | avs26 = VersionNumber() < 2.60 ? false : true | |
68 | yuy25 = !avs26 && planar | |
69 | ||
70 | ||
71 | pel = default( pel, if4 ? 1 : 2 ) | |
72 | subpixel = default( subpixel, 2 ) | |
73 | prefilter = default( prefilter, -1 ) | |
74 | MFilterB = defined(mfilter) ? true : false | |
75 | pelclip = pel>1 && subpixel == 3 ? true : false | |
76 | Str = default( Str, 1.0 ) | |
77 | Amp = default( Amp, 0.0625 ) | |
78 | ||
79 | blksize = default( blksize, if4 ? 16 : 8 ) | |
80 | blk2=blksize/2 | |
81 | overlap = default( overlap, blk2 ) | |
82 | ovl2=overlap/2 | |
83 | search = default( search, 4 ) | |
84 | truemotion = default( truemotion, !if4) | |
85 | thSCD1 = default( thSCD1, int(pow((blksize*2.5),2))) | |
86 | thSCD2 = default( thSCD2, 130 ) | |
87 | ||
88 | ||
89 | plane = default( plane, 4 ) | |
90 | chroma = default( chroma, true) | |
91 | Chr = chroma ? 3 : 1 | |
92 | Chr2 = chroma ? 3 : (prefilter==3 ? 2 : 1) | |
93 | plane0 = plane!=0 ? true : false | |
94 | U = plane0 && plane!=2 ? 3 : 2 | |
95 | V = plane0 && plane!=1 ? 3 : 2 | |
96 | Uin = lsb_in ? 3 : u | |
97 | Vin = lsb_in ? 3 : v | |
98 | ||
99 | ||
100 | hpad = blksize | |
101 | vpad = blksize | |
102 | limit = default( limit, 255 ) | |
103 | limitc = default( limitc,limit) | |
104 | ||
105 | Show = default( Show, false) | |
106 | ||
107 | ||
108 | ||
109 | ||
110 | # Error Report | |
111 | ||
112 | Assert(blksize==4 || blksize==8 || blksize==16 || blksize==32,"MAnalyse: Block's size must be 4x4, 8x4, 8x8, 16x2, 16x8, 16x16, 32x16, 32x32") | |
113 | Assert(overlap <= blk2, "Overlap must be at least half blksize or less") | |
114 | Assert(overlap == overlap-overlap%2, "MAnalyse: overlap must be an even figure") | |
115 | Assert(IsBool(Show) || sShow, "'Show' only accepts bool inputs (true,false), or string modes ('Speed','Memory',Quality')") | |
116 | Assert(preclip || IsInt(prefilter), "'prefilter' only accepts integers and clip inputs") | |
117 | Assert(ifC || IsInt(Contrasharp),"'Contrasharp' only accepts integers and bool inputs") | |
118 | !chroma ? Assert(plane==0, "Denoising chroma with luma only vectors is bugged in mvtools and thus unsupported") : nop() | |
119 | MFilterB ? Assert(IsClip(mfilter), "'mfilter' only accepts clip inputs") : nop() | |
120 | interlaced ? Assert(h%4==0, "Interlaced sources require mod 4 height sizes") : nop() | |
121 | lsb_in ? Assert(!interlaced, "Interlaced 16 bit stacked clip is not a recognized format") : nop() | |
122 | planar && if0 ? Assert(ifC, "LSFmod doesn't support YUY2 sources") : nop() | |
123 | RefineMotion ? Assert(blksize > 7, "For RefineMotion you need a blksize of at least 8") : nop() | |
124 | lsb_in ? Assert(lsb, "lsb_in requires: lsb=true") : nop() | |
125 | lsb_out ? Assert(lsb, "lsb_out requires: lsb=true") : nop() | |
126 | ||
127 | ||
128 | ||
129 | # RefineMotion Variables | |
130 | ||
131 | halfblksize = RefineMotion ? blk2 : nop() # MRecalculate works with half block size | |
132 | halfoverlap = RefineMotion ? (overlap == 2 ? overlap : ovl2+ovl2%2) : nop() # Halve the overlap to suit the halved block size | |
133 | halfthSAD = RefineMotion || if3 ? thSAD2 : nop() # MRecalculate uses a more strict thSAD, which defaults to 150 (half of function's default of 300) | |
134 | halfthSADC = if3 ? int(thSADC/2) : nop() # For MDegrainN() | |
135 | ||
136 | ||
137 | ||
138 | ||
139 | ||
140 | ||
141 | # Input preparation for: LSB_IN, Interlacing, Planar and MSuper optimization when pel=2 | |
142 | ||
143 | inputY = planar ? (lsb_in ? Dither_YUY2toPlanar16(input) : Interleaved2planar(input)) : input | |
144 | ||
145 | inputP = !interlaced ? (pel == 2 ? inputY.AssumeFrameBased() : inputY) : \ | |
146 | (if5 ? inputY.AssumeTFF().SeparateFields() : inputY.AssumeBFF().SeparateFields()) | |
147 | ||
148 | input8h = lsb_in ? inputP. Ditherpost(mode=6, slice=slices) : nop() | |
149 | input8y = planar ? (lsb_in ? input8h.Dither_YUY2toInterleaved() : inputP) : inputP | |
150 | input8 = lsb_in ? (planar ? input8y.Interleaved2planar() : \ | |
151 | input8h) : input8y | |
152 | ||
153 | ||
154 | # Prefilter & Motion Filter | |
155 | ||
156 | Mfilter = MFilterB ? Mfilter : input8 | |
157 | ||
158 | bug_wa = interlaced && planar && chroma && !avs26 ? 2 : Chr # bug: crash prevention workaround (http://forum.doom9.org/showthread.php?p=1546684#post1546684) | |
159 | ||
160 | # The mt_merge() line for prefilter=3 should be swapped with a high bitdepth variant (Dither_merge16_8() ?) due to a 1 point limited range | |
161 | # in both range ends, but then it won't work with planar sources. This isn't as critical since we are only trying to find motion vectors. | |
162 | pref = !GlobalR ? ( preclip ? prefilter : \ | |
163 | ((prefilter==-1) ? inputP : \ | |
164 | (prefilter==0) ? input8.MinBlur(0,Chr,planar) : \ | |
165 | (prefilter==1) ? input8.MinBlur(1,Chr,planar) : \ | |
166 | (prefilter==2) ? input8.MinBlur(2,Chr,planar) : \ | |
167 | (prefilter==3) ? (!planar && lsb ? Dither_merge16_8(inputP.Dfttest(sstring="0.0:4.0 0.2:9.0 1.0:15.0",tbsize=1,u=chroma,v=chroma,lsb=true,lsb_in=lsb_in,quiet=true), lsb_in?inputP:inputP.Dither_convert_8_to_16(), \ | |
168 | lsb_in?inputP.Dither_lut16("x 4096 < 65535 x 19200 > 0 65535 x 4096 - 65535 19200 4096 - / * - ? ?",u=1,v=1).Ditherpost(mode=6, slice=slices, u=Chr, v=Chr) \ | |
169 | :inputP.mt_lut( "x 16 < 255 x 75 > 0 255 x 16 - 255 75 16 - / * - ? ?",u=1,v=1), luma=chroma, u=Chr2, v=Chr2) : \ | |
170 | avs26 && planar && lsb ? Dither_merge16_8(inputP.Dfttest(sstring="0.0:4.0 0.2:9.0 1.0:15.0",tbsize=1,u=chroma,v=chroma,lsb=true,lsb_in=lsb_in,quiet=true), lsb_in?inputP:inputP.Dither_convert_8_to_16(), \ | |
171 | lsb_in?inputP.Dither_lut16("x 4096 < 65535 x 19200 > 0 65535 x 4096 - 65535 19200 4096 - / * - ? ?",u=1,v=1).Ditherpost(mode=6, slice=slices, u=Chr, v=Chr) \ | |
172 | :inputP.Planar2Interleaved().ConvertToYV16().mt_lut( "x 16 < 255 x 75 > 0 255 x 16 - 255 75 16 - / * - ? ?",u=1,v=1), luma=chroma, u=Chr2, v=Chr2).ConvertToYUY2().Interleaved2planar() : \ | |
173 | planar && avs26 ? mt_merge( Dfttest(!lsb_in?inputP.Planar2Interleaved(!chroma):input8y,sstring="0.0:4.0 0.2:9.0 1.0:15.0",tbsize=1,u=chroma,v=chroma,dither=1).ConvertToYV16(),!lsb_in?input8y.Planar2Interleaved().ConvertToYV16():input8y.ConvertToYV16(), \ | |
174 | input8y.Planar2Interleaved().ConvertToYV16().mt_lut( "x 16 < 255 x 75 > 0 255 x 16 - 255 75 16 - / * - ? ?",u=1,v=1), luma=planar?false:chroma, u=bug_wa,v=bug_wa).ConvertToYUY2().Interleaved2planar() : \ | |
175 | mt_merge( planar ? Dfttest(!lsb_in?inputP.Planar2Interleaved(!chroma):input8y,sstring="0.0:4.0 0.2:9.0 1.0:15.0",tbsize=1,u=chroma,v=chroma,dither=1).Interleaved2planar(!chroma) : \ | |
176 | Dfttest( inputP, sstring="0.0:4.0 0.2:9.0 1.0:15.0",tbsize=1,u=chroma,v=chroma,dither=1),!lsb_in?input8y:input8y.Interleaved2planar(!chroma), \ | |
177 | input8y.mt_lut( "x 16 < 255 x 75 > 0 255 x 16 - 255 75 16 - / * - ? ?",u=1,v=1), luma=planar?false:chroma, u=bug_wa,v=bug_wa)) : \ | |
178 | (prefilter==4) ? eval(""" | |
179 | # In a more lucid state I could probably have laid out this block much better... or not... | |
180 | NL_in = !planar ? (lsb_in ? inputP : input8.Dither_convert_8_to_16()) : Interlaced ? (if5 ? input.AssumeTFF().SeparateFields().Dither_convert_8_to_16() : input.AssumeBFF().SeparateFields().Dither_convert_8_to_16()) : (lsb_in ? input : input.Dither_convert_8_to_16()) | |
181 | NL_W = width(NL_in) | |
182 | U8 = chroma || planar ? UToY8(NL_in) : nop() | |
183 | V8 = chroma || planar ? VToY8(NL_in) : nop() | |
184 | NL_lsb = chroma ? (!planar ? StackVertical(StackVertical(Dither_get_msb(u8),Dither_get_msb(v8)),\ | |
185 | StackVertical(Dither_get_lsb(u8),Dither_get_lsb(v8))) :\ | |
186 | StackHorizontal(u8.ConvertToYV12(),v8.ConvertToYV12())) : nop() | |
187 | ||
188 | Y8 = chroma ? (!planar ? StackHorizontal(ConvertToY8(NL_in),NL_lsb) : \ | |
189 | StackHorizontal(StackVertical(Dither_get_msb(ConvertToY8(NL_in)),Dither_get_lsb(ConvertToY8(NL_in))).ConvertToYV12(),NL_lsb)) : \ | |
190 | (!planar ? NL_in : NL_in.converttoY8()) | |
191 | ||
192 | Y8 = Y8.KNLMeansCL(D=1, A=1, h=7.0,lsb_inout=true,device_type="GPU") | |
193 | ||
194 | my = planar ? Y8.Dither_get_msb().crop(0,0,chroma?-NL_W:0,0) : nop() | |
195 | ly = planar ? Y8.Dither_get_lsb().crop(0,0,chroma?-NL_W:0,0) : nop() | |
196 | ||
197 | chroma ? eval (" | |
198 | u_h = u8.height()/2 | |
199 | v_h = v8.height()/2 | |
200 | ||
201 | u16 = !planar ? StackVertical(Dither_get_msb(Y8).crop(NL_W,0,0,-v_h),Dither_get_lsb(Y8).crop(NL_W,0,0,-v_h)) : \ | |
202 | StackVertical(Dither_get_msb(Y8).crop(NL_W,0,-NL_W/2 ,0),Dither_get_lsb(Y8).crop(NL_W,0,-NL_W/2 ,0)) | |
203 | v16 = !planar ? StackVertical(Dither_get_msb(Y8).crop(NL_W,u_h, 0,0),Dither_get_lsb(Y8).crop(NL_W,u_h, 0,0)) : \ | |
204 | StackVertical(Dither_get_msb(Y8).crop(NL_W+NL_W/2,0,0,0),Dither_get_lsb(Y8).crop(NL_W+NL_W/2,0,0,0)) | |
205 | ||
206 | !planar ? YToUV(u16,v16,Y8.crop(0,0,-u8.width(),0)) : avs26 ? YToUV(u16, v16, StackVertical(my, ly)) : YToUV(u16, v16, StackVertical(StackVertical (my,my), StackVertical (ly,ly))) ") : \ | |
207 | ( !planar || avs26 ? Y8 : YToUV (u8, v8, StackVertical (StackVertical (my, my), StackVertical (ly, ly)))) """) : \ | |
208 | Assert(false, "Mode must be between -1~3: "+string(prefilter)))) : \ | |
209 | input8 | |
210 | ||
211 | ||
212 | ||
213 | ||
214 | ||
215 | # Default Auto-Prefilter - Luma expansion TV->PC (up to 16% more values for motion estimation) | |
216 | ||
217 | pref = avs26 && planar && ((lsb_in && prefilter==-1) || prefilter==4) ? pref.Dither_YUY2toInterleaved16().Interleaved2planar() : pref | |
218 | pref16 = (Height(pref)) == (Height(input8)*2) ? true : false | |
219 | pref = !GlobalR ? (!yuy25 && (lsb || (!preclip ? prefilter==4 : false) ) ? ((!preclip ? ((lsb_in && prefilter==-1) || prefilter==3 || prefilter==4 ) : false) ? \ | |
220 | pref.Dither_Luma_Rebuild(S0=Str,c=Amp,slice=slices,lsb_in=true,uv=Chr2) : \ | |
221 | pref.Dither_Luma_Rebuild(S0=Str,c=Amp,slice=slices,lsb =true,lsb_in=pref16,uv=Chr2)) : \ | |
222 | ((!preclip ? ((lsb_in && prefilter==-1) || prefilter==4) : false) ? \ | |
223 | pref.Ditherpost(mode=6,slice=slices,u=Chr,v=Chr).Dither_YUY2toInterleaved().Interleaved2planar(!chroma).Dither_Luma_Rebuild(S0=Str,c=Amp,uv=Chr ) : \ | |
224 | pref. Dither_Luma_Rebuild(S0=Str,lsb_in=pref16,slice=slices,c=Amp,uv=bug_wa))) : \ | |
225 | pref | |
226 | ||
227 | ||
228 | ||
229 | ||
230 | ||
231 | (sShow ? false : !Show) ? eval(""" | |
232 | ||
233 | ||
234 | # Subpixel 3 | |
235 | ||
236 | pclip = pelclip ? pel==4 ? pref.nnedi3_rpow2(rfactor=2,nns=4).nnedi3_rpow2(rfactor=2,nns=4) : pref.nnedi3_rpow2(rfactor=2,nns=4) : nop() | |
237 | pclip2 = pelclip && !GlobalR ? pel==4 ? input8.nnedi3_rpow2(rfactor=2,nns=4).nnedi3_rpow2(rfactor=2,nns=4) : input8.nnedi3_rpow2(rfactor=2,nns=4) : nop() | |
238 | ||
239 | ||
240 | ||
241 | ||
242 | # Motion vectors search | |
243 | ||
244 | ||
245 | super_search = pelclip ? MSuper(pref, pel=pel, chroma=chroma, hpad=hpad, vpad=vpad, pelclip=pclip, rfilter=4, planar=planar) : \ | |
246 | MSuper(pref, pel=pel,sharp=subpixel,chroma=chroma, hpad=hpad, vpad=vpad , rfilter=4, planar=planar) | |
247 | super_render = !GlobalR ? (pelclip ? MSuper(input8, pel=pel, chroma=plane0, hpad=hpad, vpad=vpad,levels=1, pelclip=pclip2 , planar=planar) : \ | |
248 | MSuper(input8, pel=pel,sharp=subpixel,chroma=plane0, hpad=hpad, vpad=vpad,levels=1 , planar=planar)): super_search | |
249 | Recalculate = !GlobalR && RefineMotion ? (pelclip ? MSuper(pref, pel=pel, chroma=chroma, hpad=hpad, vpad=vpad,levels=1, pelclip=pclip , planar=planar) : \ | |
250 | MSuper(pref, pel=pel,sharp=subpixel,chroma=chroma, hpad=hpad, vpad=vpad,levels=1 , planar=planar)): nop() | |
251 | !if3 ? eval(" | |
252 | bv6 = interlaced && (tr > 2) ? (!GlobalR ? super_search.MAnalyse(isb = true, delta = 6, overlap= overlap,blksize= blksize, search=search,chroma=chroma, truemotion=truemotion) : bv6) : nop() | |
253 | bv6 = interlaced && (tr > 2) && RefineMotion && !GlobalR ? MRecalculate (Recalculate, bv6, overlap=halfoverlap,blksize=halfblksize, thSAD=halfthSAD,chroma=chroma, truemotion=truemotion) : bv6 | |
254 | bv4 = interlaced && (tr > 1) ? (!GlobalR ? super_search.MAnalyse(isb = true, delta = 4, overlap= overlap,blksize= blksize, search=search,chroma=chroma, truemotion=truemotion) : bv4) : nop() | |
255 | bv4 = interlaced && (tr > 1) && RefineMotion && !GlobalR ? MRecalculate (Recalculate, bv4, overlap=halfoverlap,blksize=halfblksize, thSAD=halfthSAD,chroma=chroma, truemotion=truemotion) : bv4 | |
256 | bv3 = !interlaced ? ((tr > 2) ? (!GlobalR ? super_search.MAnalyse(isb = true, delta = 3, overlap= overlap,blksize= blksize, search=search,chroma=chroma, truemotion=truemotion) : bv3) : nop()) : nop() | |
257 | bv3 = !interlaced ? ((tr > 2 && RefineMotion && !GlobalR) ? MRecalculate (Recalculate, bv3, overlap=halfoverlap,blksize=halfblksize, thSAD=halfthSAD,chroma=chroma, truemotion=truemotion) : bv3 ) : nop() | |
258 | bv2 = (interlaced || tr > 1) ? (!GlobalR ? super_search.MAnalyse(isb = true, delta = 2, overlap= overlap,blksize= blksize, search=search,chroma=chroma, truemotion=truemotion) : bv2) : nop() | |
259 | bv2 = ((interlaced || tr > 1) && RefineMotion && !GlobalR) ? MRecalculate (Recalculate, bv2, overlap=halfoverlap,blksize=halfblksize, thSAD=halfthSAD,chroma=chroma, truemotion=truemotion) : bv2 | |
260 | bv1 = !interlaced ? (!GlobalR ? super_search.MAnalyse(isb = true, delta = 1, overlap= overlap,blksize= blksize, search=search,chroma=chroma, truemotion=truemotion) : bv1 ) : nop() | |
261 | bv1 = !interlaced ? ( RefineMotion && !GlobalR ? MRecalculate (Recalculate, bv1, overlap=halfoverlap,blksize=halfblksize, thSAD=halfthSAD,chroma=chroma, truemotion=truemotion) : bv1 ) : nop() | |
262 | ||
263 | fv1 = !interlaced ? (!GlobalR ? super_search.MAnalyse(isb = false, delta = 1, overlap= overlap,blksize= blksize, search=search,chroma=chroma, truemotion=truemotion) : fv1 ) : nop() | |
264 | fv1 = !interlaced ? ( RefineMotion && !GlobalR ? MRecalculate (Recalculate, fv1, overlap=halfoverlap,blksize=halfblksize, thSAD=halfthSAD,chroma=chroma, truemotion=truemotion) : fv1 ) : nop() | |
265 | fv2 = (interlaced || tr > 1) ? (!GlobalR ? super_search.MAnalyse(isb = false, delta = 2, overlap= overlap,blksize= blksize, search=search,chroma=chroma, truemotion=truemotion) : fv2) : nop() | |
266 | fv2 = ((interlaced || tr > 1) && RefineMotion && !GlobalR) ? MRecalculate (Recalculate, fv2, overlap=halfoverlap,blksize=halfblksize, thSAD=halfthSAD,chroma=chroma, truemotion=truemotion) : fv2 | |
267 | fv3 = !interlaced ? ((tr > 2) ? (!GlobalR ? super_search.MAnalyse(isb = false, delta = 3, overlap= overlap,blksize= blksize, search=search,chroma=chroma, truemotion=truemotion) : fv3) : nop()) : nop() | |
268 | fv3 = !interlaced ? ((tr > 2 && RefineMotion && !GlobalR) ? MRecalculate (Recalculate, fv3, overlap=halfoverlap,blksize=halfblksize, thSAD=halfthSAD,chroma=chroma, truemotion=truemotion) : fv3 ) : nop() | |
269 | fv4 = interlaced && (tr > 1) ? (!GlobalR ? super_search.MAnalyse(isb = false, delta = 4, overlap= overlap,blksize= blksize, search=search,chroma=chroma, truemotion=truemotion) : fv4) : nop() | |
270 | fv4 = interlaced && (tr > 1) && RefineMotion && !GlobalR ? MRecalculate (Recalculate, fv4, overlap=halfoverlap,blksize=halfblksize, thSAD=halfthSAD,chroma=chroma, truemotion=truemotion) : fv4 | |
271 | fv6 = interlaced && (tr > 2) ? (!GlobalR ? super_search.MAnalyse(isb = false, delta = 6, overlap= overlap,blksize= blksize, search=search,chroma=chroma, truemotion=truemotion) : fv6) : nop() | |
272 | fv6 = interlaced && (tr > 2) && RefineMotion && !GlobalR ? MRecalculate (Recalculate, fv6, overlap=halfoverlap,blksize=halfblksize, thSAD=halfthSAD,chroma=chroma, truemotion=truemotion) : fv6 " ) : \ | |
273 | \ | |
274 | eval(" | |
275 | tr2 = tr*2 | |
276 | vmulti = !GlobalR ? super_search.MAnalyse( multi=true,overlap= overlap,blksize= blksize, search=search,chroma=chroma,truemotion=truemotion, delta = interlaced ? tr2 : tr) : vmulti | |
277 | vmulti = RefineMotion ? (!GlobalR ? MRecalculate (Recalculate, vmulti,overlap=halfoverlap,blksize=halfblksize, thSAD=halfthSAD,chroma=chroma,truemotion=truemotion, tr = interlaced ? tr2 : tr) : vmulti) : vmulti | |
278 | vmulti = !GlobalR ? (interlaced ? vmulti.SelectEvery (4, 2, 3) : vmulti) : vmulti.SelectRangeEvery (Rtr*2, tr2, 0, false)") | |
279 | ||
280 | ||
281 | ||
282 | # Finally, MDegrain | |
283 | ||
284 | ||
285 | output = !GlobalO ? (interlaced ? | |
286 | \ ((lsb ? ( if3 ? mfilter.MDegrainN(super_render, vmulti, tr, thSAD2=halfthSAD, thSADC2=halfthSADC, thSAD=thSAD, thSADC=thSADC, thSCD1=thSCD1, thSCD2=thSCD2, limit=limit, limitc=limitc, plane=plane, planar=planar,lsb=true) : | |
287 | \ (tr == 3) ? mfilter.MDegrain3(super_render, bv2, fv2, bv4, fv4, bv6, fv6, thSAD=thSAD, thSADC=thSADC, thSCD1=thSCD1, thSCD2=thSCD2, limit=limit, limitc=limitc, plane=plane, planar=planar,lsb=true) : | |
288 | \ (tr == 2) ? mfilter.MDegrain2(super_render, bv2, fv2, bv4, fv4, thSAD=thSAD, thSADC=thSADC, thSCD1=thSCD1, thSCD2=thSCD2, limit=limit, limitc=limitc, plane=plane, planar=planar,lsb=true) : | |
289 | \ mfilter.MDegrain1(super_render, bv2, fv2, thSAD=thSAD, thSADC=thSADC, thSCD1=thSCD1, thSCD2=thSCD2, limit=limit, limitc=limitc, plane=plane, planar=planar,lsb=true)): | |
290 | \ ( if3 ? mfilter.MDegrainN(super_render, vmulti, tr, thSAD2=halfthSAD, thSADC2=halfthSADC, thSAD=thSAD, thSADC=thSADC, thSCD1=thSCD1, thSCD2=thSCD2, limit=limit, limitc=limitc, plane=plane, planar=planar,lsb=false): | |
291 | \ (tr == 3) ? mfilter.MDegrain3(super_render, bv2, fv2, bv4, fv4, bv6, fv6, thSAD=thSAD, thSADC=thSADC, thSCD1=thSCD1, thSCD2=thSCD2, limit=limit, limitc=limitc, plane=plane, planar=planar) : | |
292 | \ (tr == 2) ? mfilter.MDegrain2(super_render, bv2, fv2, bv4, fv4, thSAD=thSAD, thSADC=thSADC, thSCD1=thSCD1, thSCD2=thSCD2, limit=limit, limitc=limitc, plane=plane, planar=planar) : | |
293 | \ mfilter.MDegrain1(super_render, bv2, fv2, thSAD=thSAD, thSADC=thSADC, thSCD1=thSCD1, thSCD2=thSCD2, limit=limit, limitc=limitc, plane=plane, planar=planar)))) : | |
294 | \ (lsb ? ( if3 ? mfilter.MDegrainN(super_render, vmulti, tr, thSAD2=halfthSAD, thSADC2=halfthSADC, thSAD=thSAD, thSADC=thSADC, thSCD1=thSCD1, thSCD2=thSCD2, limit=limit, limitc=limitc, plane=plane, planar=planar,lsb=true) : | |
295 | \ (tr == 3) ? mfilter.MDegrain3(super_render, bv1, fv1, bv2, fv2, bv3, fv3, thSAD=thSAD, thSADC=thSADC, thSCD1=thSCD1, thSCD2=thSCD2, limit=limit, limitc=limitc, plane=plane, planar=planar,lsb=true) : | |
296 | \ (tr == 2) ? mfilter.MDegrain2(super_render, bv1, fv1, bv2, fv2, thSAD=thSAD, thSADC=thSADC, thSCD1=thSCD1, thSCD2=thSCD2, limit=limit, limitc=limitc, plane=plane, planar=planar,lsb=true) : | |
297 | \ mfilter.MDegrain1(super_render, bv1, fv1, thSAD=thSAD, thSADC=thSADC, thSCD1=thSCD1, thSCD2=thSCD2, limit=limit, limitc=limitc, plane=plane, planar=planar,lsb=true)): | |
298 | \ ( if3 ? mfilter.MDegrainN(super_render, vmulti, tr, thSAD2=halfthSAD, thSADC2=halfthSADC, thSAD=thSAD, thSADC=thSADC, thSCD1=thSCD1, thSCD2=thSCD2, limit=limit, limitc=limitc, plane=plane, planar=planar,lsb=false): | |
299 | \ (tr == 3) ? mfilter.MDegrain3(super_render, bv1, fv1, bv2, fv2, bv3, fv3, thSAD=thSAD, thSADC=thSADC, thSCD1=thSCD1, thSCD2=thSCD2, limit=limit, limitc=limitc, plane=plane, planar=planar) : | |
300 | \ (tr == 2) ? mfilter.MDegrain2(super_render, bv1, fv1, bv2, fv2, thSAD=thSAD, thSADC=thSADC, thSCD1=thSCD1, thSCD2=thSCD2, limit=limit, limitc=limitc, plane=plane, planar=planar) : | |
301 | \ mfilter.MDegrain1(super_render, bv1, fv1, thSAD=thSAD, thSADC=thSADC, thSCD1=thSCD1, thSCD2=thSCD2, limit=limit, limitc=limitc, plane=plane, planar=planar)))) : nop() | |
302 | ||
303 | ||
304 | # LSB_IN merging | |
305 | ||
306 | outputI= !GlobalO && planar && (lsb || ifC) ? output.Planar2Interleaved() : nop() | |
307 | !GlobalO && lsb ? eval(" | |
308 | blnkclp= lsb_in || (if2 && !planar) ? BlankClip (lsb_in?input8h:input8, pixel_type="+Chr(34)+"YV12"+Chr(34)+", color_yuv=0) : nop() | |
309 | output = lsb && planar ? Dither_YUY2toPlanar16(outputI) : output | |
310 | output = lsb_in ? inputY.Dither_add16 (Dither_sub16(output, planar?inputY:input8.Dither_convert_8_to_16 (), dif=true,u=u,v=v), dif=true,u=u,v=v) : output") : nop() | |
311 | ||
312 | ||
313 | ||
314 | ||
315 | # Contrasharp (only sharpens luma) | |
316 | ||
317 | ifC && if0 && if4 ? eval(" | |
318 | global Super = super_render | |
319 | global cb1 = !if3 ? (interlaced ? bv2 : bv1) : vmulti.SelectEvery (tr2, 0) | |
320 | global cf1 = !if3 ? (interlaced ? fv2 : fv1) : vmulti.SelectEvery (tr2, 1)") : nop() | |
321 | ||
322 | !GlobalO && if0 ? eval(" | |
323 | ||
324 | CCh = if1 ? CClip.height() : nop() | |
325 | ref8 = if1 ? (h == CCh*2) : nop() | |
326 | ref16 = if1 ? (h == CCh/2) : nop() | |
327 | ||
328 | interlaced && if1 ? Assert(!ref16, "+Chr(34)+"Interlaced 16 bit stacked CClip is not a recognized format"+Chr(34)+") : nop() | |
329 | ||
330 | CClip = planar ? (if1 ? (lsb_in ? (ref8 ? CClip.Interleaved2planar(true) : CClip.Dither_YUY2toPlanar16().DitherPost(mode=-1,u=1,v=1).Dither_YUY2toInterleaved().Interleaved2planar(true)) : \ | |
331 | (interlaced ? (if5 ? CClip.AssumeTFF().SeparateFields().Interleaved2planar(true) : CClip.AssumeBFF().SeparateFields().Interleaved2planar(true)) : \ | |
332 | (ref16 ? CClip.Dither_YUY2toPlanar16().DitherPost(mode=-1,u=1,v=1).Dither_YUY2toInterleaved().Interleaved2planar(true) : CClip.Interleaved2planar(true)))) : \ | |
333 | (lsb_in ? inputY.DitherPost(mode=-1,u=1,v=1).Dither_YUY2toInterleaved().Interleaved2planar(true) : \ | |
334 | inputP)) : \ | |
335 | (if1 ? (lsb_in ? (ref8 ? CClip : CClip.DitherPost(mode=-1,u=1,v=1)) : \ | |
336 | (interlaced ? (if5 ? CClip.AssumeTFF().SeparateFields() : CClip.AssumeBFF().SeparateFields()) : (ref16 ? CClip.DitherPost(mode=-1,u=1,v=1) : CClip))) : \ | |
337 | (lsb_in ? input.DitherPost(mode=-1,u=1,v=1) : \ | |
338 | inputP))") : nop() | |
339 | ||
340 | ||
341 | !GlobalO && if2 ? eval(" | |
342 | OutTO8 = output.DitherPost(mode=-1,u=1,v=1) | |
343 | OutTO8P = planar ? OutTO8.Dither_YUY2toInterleaved() : nop() | |
344 | ctr16 = ifC ? ContrasharpeningHD (planar ? OutTO8P.Interleaved2planar(true):OutTO8, CClip, planar=planar, HD=if4) : \ | |
345 | LSFmod(OutTO8,source=CClip,defaults="+Chr(34)+"slow"+Chr(34)+",strength=Contrasharp,soothe=false,Lmode=0,soft=-2,keep=0) | |
346 | ||
347 | ||
348 | ctr16P = planar ? ctr16.Planar2Interleaved(true) : ctr16 | |
349 | ctr16 = StackVertical (ctr16P, planar ? BlankClip (ctr16P, pixel_type="+Chr(34)+"YUY2"+Chr(34)+", color_yuv=0,length=1) : blnkclp ) | |
350 | ctr16 = Dither_add16 (output,mt_makediff(!planar?ctr16P:ctr16P.Dither_YUY2toPlanar(), OutTO8,u=2,v=2).Dither_convert_8_to_16(), dif=true,u=2,v=2) | |
351 | ctr16 = MergeChroma(ctr16, output) ") : nop() | |
352 | ||
353 | ||
354 | ||
355 | # Globals Output | |
356 | ||
357 | GlobalO || Globals == 2 ? (!if3 ? eval(" | |
358 | global Super = super_render | |
359 | global bv6 = bv6 | |
360 | global bv4 = bv4 | |
361 | global bv3 = bv3 | |
362 | global bv2 = bv2 | |
363 | global bv1 = bv1 | |
364 | ||
365 | global fv1 = fv1 | |
366 | global fv2 = fv2 | |
367 | global fv3 = fv3 | |
368 | global fv4 = fv4 | |
369 | global fv6 = fv6 | |
370 | ") : eval(" | |
371 | global bv1 = vmulti.SelectEvery (tr2, 0) | |
372 | global bv2 = interlaced ? bv1 : vmulti.SelectEvery (tr2, 2) | |
373 | global bv3 = vmulti.SelectEvery (tr2, 4) | |
374 | global bv4 = interlaced ? vmulti.SelectEvery (tr2, 2) : nop() | |
375 | global bv6 = interlaced ? bv3 : nop() | |
376 | ||
377 | global fv1 = vmulti.SelectEvery (tr2, 1) | |
378 | global fv2 = interlaced ? fv1 : vmulti.SelectEvery (tr2, 3) | |
379 | global fv3 = vmulti.SelectEvery (tr2, 5) | |
380 | global fv4 = interlaced ? vmulti.SelectEvery (tr2, 3) : nop() | |
381 | global fv6 = interlaced ? fv3 : nop() | |
382 | global Super = super_render | |
383 | global vmulti = vmulti | |
384 | global Rtr = tr")) : nop() | |
385 | ||
386 | ||
387 | ||
388 | # Output | |
389 | ||
390 | !GlobalO ? (lsb_out ? (if0 ? ctr16 : output) : \ | |
391 | (if0 ? (lsb ? (interlaced ? ctr16.weave().DitherPost(u=u,v=v,mode=6,slice=slices,interlaced=true) : ctr16.DitherPost(mode=mode,slice=slices,ampn=ampn,u=Uin,v=Vin)) : \ | |
392 | (interlaced ? (ifC ? ContraSharpeningHD(output, CClip, planar=planar, HD=if4).weave() : LSFmod(output,source=CClip,defaults="slow",strength=Contrasharp,soothe=false,Lmode=0,soft=-2,keep=0).weave()) : \ | |
393 | (ifC ? ContraSharpeningHD(output, CClip, planar=planar, HD=if4) : LSFmod(output,source=CClip,defaults="slow",strength=Contrasharp,soothe=false,Lmode=0,soft=-2,keep=0)))) : \ | |
394 | (lsb ? (interlaced ? output.weave().DitherPost(u=u,v=v,mode=6,slice=slices,interlaced=true) : output.DitherPost(mode=mode,slice=slices,ampn=ampn,u=Uin,v=Vin)) : (interlaced ? output.weave() : output)))) \ | |
395 | : input | |
396 | ||
397 | !GlobalO ? (planar ? (lsb ? Dither_YUY2toInterleaved() : (if0 ? MergeChroma(Planar2Interleaved(true), interlaced ? outputI.weave() : outputI) : Planar2Interleaved())) : last) : last """) : \ | |
398 | \ | |
399 | \ | |
400 | \ | |
401 | \ | |
402 | \ | |
403 | eval(""" | |
404 | ||
405 | ############## | |
406 | # SHOW PANEL # | |
407 | ############## | |
408 | ||
409 | ||
410 | # Common Variables | |
411 | ||
412 | h = lsb_in ? h/2 : h | |
413 | w2 = w*1.0 | |
414 | h2 = h*1.0 | |
415 | hi = interlaced ? round(h/2.0) : h | |
416 | AR = w2/h2 | |
417 | lsp = interlaced?15:30 | |
418 | lsp = AR < 1.9 ? lsp : floor(AR < 1.9 ? lsp : ((AR>2.4) ? 0 : lsp-((AR-1.9)*(lsp/(2.4-1.9))))) | |
419 | ys = round(lsb_in ? h2/6 : (interlaced?h2/15:h2/6)) | |
420 | s = floor(((interlaced?w2/1.6:w2)/720.0)*16.0) | |
421 | type = planar?"YUY2":"YV12" | |
422 | ||
423 | Wadd = w2/40 | |
424 | Wadd2 = Wadd/2 | |
425 | w1 = round(w2+Wadd+Wadd2) | |
426 | w3 = (w2/3)+Wadd/2.1 | |
427 | w2 = round(w1+w3) | |
428 | w3 = round(w1+(w3+Wadd2)*1.96) | |
429 | w3sq = string(round(w3*1.01)) | |
430 | w3s = string(w3) | |
431 | ||
432 | # Panel Separator | |
433 | ||
434 | blkclp = BlankClip(interlaced?assumeframebased(pref):pref,length=1,pixel_type=type,color=$161616) | |
435 | \.crop(0,interlaced?h/2%8:h%8,0,0).separatefields().assumeframebased().separatefields() | |
436 | ||
437 | Interleave(selecteven(blkclp).addborders(2,0,0,0,color=$7f7f7f).crop(0,0,-2,0),selectodd(blkclp)) | |
438 | ||
439 | WOut=weave().assumefieldbased().weave() | |
440 | sep=round((hi-height(WOut))/2.0) | |
441 | ||
442 | StackHorizontal(planar? pref.Planar2Interleaved() : pref,WOut.addborders(0,sep-sep%2,0,hi-(height(WOut)+(sep-sep%2)),color=$161616)) | |
443 | ||
444 | ||
445 | ||
446 | # Color Guide Squares | |
447 | ||
448 | a=interlaced?h*0.37:h*0.74 | |
449 | b=interlaced?h*0.40:h*0.80 | |
450 | c=interlaced?h*0.43:h*0.86 | |
451 | ||
452 | !planar ? mt_merge(mt_lut(yexpr="127",uexpr="183",vexpr="47" ,u=3,v=3) ,mt_lutspa(mode="absolute", yexpr="x "+w3s+" >= x "+w3sq+" <= & y "+string(a)+" >= y "+string(interlaced?h*0.3835:h*0.767)+" <= & & 255 0 ?").trim(0,-1).FreezeFrame(0, FrameCount(pref)-1, 0),luma=true) : nop() | |
453 | !planar ? mt_merge(mt_lut(yexpr="81" ,uexpr="91" ,vexpr="240",u=3,v=3) ,mt_lutspa(mode="absolute", yexpr="x "+w3s+" >= x "+w3sq+" <= & y "+string(b)+" >= y "+string(interlaced?h*0.4135:h*0.827)+" <= & & 255 0 ?").trim(0,-1).FreezeFrame(0, FrameCount(pref)-1, 0),luma=true) : nop() | |
454 | !planar ? mt_merge(mt_lut(yexpr="210",uexpr="16" ,vexpr="146",u=3,v=3) ,mt_lutspa(mode="absolute", yexpr="x "+w3s+" >= x "+w3sq+" <= & y "+string(c)+" >= y "+string(interlaced?h*0.4435:h*0.887)+" <= & & 255 0 ?").trim(0,-1).FreezeFrame(0, FrameCount(pref)-1, 0),luma=true) : nop() | |
455 | ||
456 | ||
457 | # Parameters List | |
458 | ||
459 | blue = sShow ? (show=="Speed" ? true : false) : false | |
460 | red = sShow ? (show=="Memory" ? true : false) : false | |
461 | yellow= sShow ? (show=="Quality" ? true : false) : false | |
462 | ||
463 | Color = sShow ? (blue ? $00aef0 : \ | |
464 | red ? $ff0000 : \ | |
465 | yellow ? $ffff00 : \ | |
466 | Assert(false, "Revise typo for possible color guide modes: 'Speed','Memory','Quality'")) : nop() | |
467 | ||
468 | ||
469 | ||
470 | varn = "\n\n\n\n\n" | |
471 | sp = !GlobalR && str==1.0 ? " " : "" | |
472 | ||
473 | Subtitle("SMDegrain Mod v3.1.2",text_color=red||blue||yellow?Color:$7f7f7f,halo_color=$000000,font="VERDANA",size=s*2,x=w1,y=lsb_in?h/40:(interlaced?h/50:h/20)) | |
474 | ||
475 | Subtitle("BASIC" ,lsp=lsp,text_color=$ffffff,font="COURIER NEW",size=s,x=w1,y=ys) | |
476 | Subtitle( "\n\ntr = " + string(tr) + \ | |
477 | varn+varn+ "chroma = " + string(chroma) ,lsp=lsp,text_color=red||blue ?Color:$ffffff,font="COURIER NEW",size=s,x=w1,y=ys) | |
478 | Subtitle( "\n\n\nthSAD(C) = " + string( GlobalO?"bypass":string(thSAD) + "(" + string(plane0?thSADC:"bypass") + ")")+ \ | |
479 | varn+ "CClip = " + string( GlobalO?"bypass":(defined(CClip)?"defined":"undefined")) + \ | |
480 | "\n\nInterlaced = " + string(interlaced) ,lsp=lsp,text_color= $ffffff,font="COURIER NEW",size=s,x=w1,y=ys) | |
481 | Subtitle(varn+ "RefineMotion = " + string( GlobalO?"bypass":RefineMotion) ,lsp=lsp,text_color=yellow ?Color:$ffffff,font="COURIER NEW",size=s,x=w1,y=ys) | |
482 | Subtitle(varn+ "\n\nContrasharp = " + string( GlobalO?"bypass":contrasharp) ,lsp=lsp,text_color=yellow&&!GlobalO ?Color:$ffffff,font="COURIER NEW",size=s,x=w1,y=ys) | |
483 | Subtitle(varn+varn+ "\nplane = " + string( GlobalO?"bypass":(plane==0?"0 (Luma)":plane==1?"1 (U plane)":plane==2?"2 (V plane)":plane==3?"3 (Chroma)":"4(Luma+Chroma)")),lsp=lsp,text_color=(red||blue)&&!GlobalO?Color:$ffffff,font="COURIER NEW",size=s,x=w1,y=ys) | |
484 | Subtitle("<-- "+string(IsClip(prefilter)||prefilter!=-1?"Prefiltered clip converted\n to PC levels":"Clip converted to PC levels\n ")+" for the motion search",lsp=lsp,text_color=$ffffff,font="COURIER NEW",size=s,x=w1-10,y=(interlaced?h/2:h)-int(ys/1.5)) | |
485 | ||
486 | Subtitle("ADVANCED",lsp=lsp,text_color=$7f7f7f,font="COURIER NEW",size=s,x=w2,y=ys) | |
487 | Subtitle( "\n\npel = " + string(pel) ,lsp=lsp,text_color=red ?Color:$7f7f7f,font="COURIER NEW",size=s,x=w2,y=ys) | |
488 | Subtitle( "\n\n\nsubpixel = " + string( pel==1 ?"bypass":subpixel) ,lsp=lsp,text_color=(yellow||blue)&&pel!=1 ?Color:$7f7f7f,font="COURIER NEW",size=s,x=w2,y=ys) | |
489 | Subtitle(varn+varn+ "\n\n\ntruemotion = " + string( GlobalR?"bypass":truemotion) ,lsp=lsp,text_color=(yellow||blue)&&!GlobalR?Color:$7f7f7f,font="COURIER NEW",size=s,x=w2,y=ys) | |
490 | Subtitle(varn+ "mfilter = " + string( GlobalR?"bypass":MFilterB?"Clip"+sp+" Variable":"OFF") ,lsp=lsp,text_color= $7f7f7f,font="COURIER NEW",size=s,x=w2,y=ys) | |
491 | Subtitle(varn+ "\n\nStr = " + (GlobalR?string("bypass"):string(Str,"%1.2f"))+sp+" Brightness" + \ | |
492 | varn+ "search = " + string( GlobalR?"bypass":search) ,lsp=lsp,text_color=yellow&&!GlobalR ?Color:$7f7f7f,font="COURIER NEW",size=s,x=w2,y=ys) | |
493 | Subtitle(varn+varn+ "blksize = " + string( GlobalR?"bypass":blksize) + \ | |
494 | "\noverlap = " + string( GlobalR?"bypass":overlap) ,lsp=lsp,text_color=(red||blue)&&!GlobalR ?Color:$7f7f7f,font="COURIER NEW",size=s,x=w2,y=ys) | |
495 | Subtitle(varn+ "\nprefilter = " + (GlobalR?string("bypass"):(preclip?"Clip"+sp+" Variable":(prefilter==0?"0"+sp+" MinBlur(0)":prefilter==1?"1"+sp+" MinBlur(1)":prefilter==2?"2"+sp+" MinBlur(2)":prefilter==3?"3"+sp+" Dfttest":prefilter==4?"4"+sp+" KNLmeans":"OFF"+sp+" (-1)"))) + \ | |
496 | "\n\nAmp = " + (GlobalR||str==1.0?string("bypass"):string(Amp,"%1.2f"))+" Amplitude" + \ | |
497 | varn+ "\nlimit(C) = " + (GlobalO?string("bypass"):string(limit) + "(" + string(limitc) + ")") ,lsp=lsp,text_color= $7f7f7f,font="COURIER NEW",size=s,x=w2,y=ys) | |
498 | ||
499 | Subtitle("OTHERS" ,lsp=lsp,text_color=$494949,font="COURIER NEW",size=s,x=w3,y=ys) | |
500 | Subtitle( "\n\n\nthSCD1(2) = " + (GlobalO?string("bypass"):string(thSCD1) + "(" + string(thSCD2) + ")") + \ | |
501 | varn+"\n\nGlobals = " + string(Globals==0 ? "0 (OFF)" : Globals==1 ? "1 (Import)" : Globals==2 ? "2 (Export)" : "3 (E.Only)"),lsp=lsp,text_color= $494949,font="COURIER NEW",size=s,x=w3,y=ys) | |
502 | Subtitle( varn+"\n\nslices = " + string(slices) ,lsp=lsp,text_color=blue ?Color:$494949,font="COURIER NEW",size=s,x=w3,y=ys) | |
503 | Subtitle( "\n\nhpad(v) = " + string(hpad) + "(" + string(vpad) + ")" ,lsp=lsp,text_color=red||blue ?Color:$494949,font="COURIER NEW",size=s,x=w3,y=ys) | |
504 | Subtitle(varn+ "lsb = " + string(lsb) ,lsp=lsp,text_color=yellow ?Color:$494949,font="COURIER NEW",size=s,x=w3,y=ys) | |
505 | Subtitle(varn+ "\nlsb_in/out= " + string(lsb_in) +"/"+ string(lsb_out) ,lsp=lsp,text_color=red&&lsb_in ?Color:$494949,font="COURIER NEW",size=s,x=w3,y=ys) | |
506 | Subtitle(varn+ "\n\n\nmode = " + (GlobalO||(lsb_out||!lsb)?string("bypass"):interlaced?"6 (locked)":string(mode)) ,lsp=lsp,text_color=yellow&&!GlobalO&&lsb ?Color:$494949,font="COURIER NEW",size=s,x=w3,y=ys) | |
507 | ||
508 | Subtitle(" Speed" ,lsp=interlaced?15:30,text_color=planar?$00aef0:$7f7f7f,font="COURIER NEW",size=s,x=w3,y=int(a-1)) | |
509 | Subtitle(" Memory" ,lsp=interlaced?17:34,text_color=planar?$ff0000:$7f7f7f,font="COURIER NEW",size=s,x=w3,y=int(b-1)) | |
510 | Subtitle(" Quality" ,lsp=interlaced?19:38,text_color=planar?$ffff00:$7f7f7f,font="COURIER NEW",size=s,x=w3,y=int(c-1)) | |
511 | ||
512 | lsb_out ? StackVertical (last, BlankClip (last, pixel_type=type, color_yuv=0,length=1)) : last""")} | |
513 | ||
514 | ||
515 | ||
516 | ||
517 | ||
518 | ######################################################## | |
519 | # # | |
520 | # HELPER FUNCTIONS: For Contrasharp and YUY2 handling # | |
521 | # # | |
522 | ######################################################## | |
523 | ||
524 | ######################################## | |
525 | ## Didée's functions: | |
526 | ||
527 | # contra-sharpening: sharpen the denoised clip, but don't add more to any pixel than what was removed previously. | |
528 | # script function from Didée, at the VERY GRAINY thread (http://forum.doom9.org/showthread.php?p=1076491#post1076491) | |
529 | ||
530 | # In final version 2.0d ContraHD() was merged, to allow proper HD sharpening. | |
531 | # In this case global variables of (before denoising) source MSuper and forward and backward | |
532 | # compensated motion vectors are necessary as: Super, cb1 and cf1, if used as a stand alone function. | |
533 | # Don't know who made (mod) it, so I can't give proper credits, sorry. | |
534 | ||
535 | FUNCTION ContraSharpeningHD(clip denoised, clip original, bool "HD", bool "planar", int "overshoot"){ | |
536 | HD = default(HD ,false) | |
537 | planar = default(planar,false) | |
538 | overshoot = default(overshoot,0) | |
539 | avs26 = VersionNumber() < 2.60 ? false : true | |
540 | yuy26 = avs26 && Isyuy2(denoised) && planar | |
541 | ||
542 | HD ? eval(""" | |
543 | cb1=original.MCompensate(Super, cb1, planar=planar) | |
544 | - | pmax = yuy26 ? original.Planar2Interleaved().ConvertToYV16().mt_logic(cb1.Planar2Interleaved().ConvertToYV16(), "max").mt_logic(cf1.Planar2Interleaved().ConvertToYV16(), "max") : original.mt_logic(cb1, "max").mt_logic(cf1, "max").ConvertToYUY2().Interleaved2planar() |
544 | + | |
545 | - | pmin = yuy26 ? original.Planar2Interleaved().ConvertToYV16().mt_logic(cb1.Planar2Interleaved().ConvertToYV16(), "min").mt_logic(cf1.Planar2Interleaved().ConvertToYV16(), "min") : original.mt_logic(cb1, "min").mt_logic(cf1, "min").ConvertToYUY2().Interleaved2planar()""") : nop() |
545 | + | pmax = yuy26 ? original.Planar2Interleaved().ConvertToYV16().mt_logic(cb1.Planar2Interleaved().ConvertToYV16(), "max").mt_logic(cf1.Planar2Interleaved().ConvertToYV16(), "max") : original.mt_logic(cb1, "max").mt_logic(cf1, "max") |
546 | pmin = yuy26 ? original.Planar2Interleaved().ConvertToYV16().mt_logic(cb1.Planar2Interleaved().ConvertToYV16(), "min").mt_logic(cf1.Planar2Interleaved().ConvertToYV16(), "min") : original.mt_logic(cb1, "min").mt_logic(cf1, "min")""") : nop() | |
547 | ||
548 | s = denoised.MinBlur(HD?2:1,1,planar=planar) # Damp down remaining spots of the denoised clip. | |
549 | s = yuy26 ? s.Planar2Interleaved().ConvertToYV16() : s | |
550 | original = yuy26 ? original.Planar2Interleaved().ConvertToYV16() : original | |
551 | denoised = yuy26 ? denoised.Planar2Interleaved().ConvertToYV16() : denoised | |
552 | allD = mt_makediff(original,denoised) # The difference achieved by the denoising. | |
553 | ssD = mt_makediff(s,HD?s.removegrain(20,-1,planar=planar).\ | |
554 | removegrain(20,-1,planar=planar):\ | |
555 | s.removegrain(11,-1,planar=planar)) # The difference of a simple kernel blur. | |
556 | allD = yuy26 ? allD.ConvertToYUY2().Interleaved2planar() : allD | |
557 | ssD = yuy26 ? ssD.ConvertToYUY2().Interleaved2planar() : ssD | |
558 | ssDD = ssD.repair(HD?ssD.repair(allD,1,planar=planar):allD,HD?12:1,planar=planar) # Limit the difference to the max of what the denoising removed locally. | |
559 | ssDD = yuy26 ? ssDD.Planar2Interleaved().ConvertToYV16() : ssDD | |
560 | ssD = yuy26 ? ssD.Planar2Interleaved().ConvertToYV16() : ssD | |
561 | ssDD = SSDD.mt_lutxy(ssD,"x 128 - abs y 128 - abs < x y ?") # abs(diff) after limiting may not be bigger than before. | |
562 | ||
563 | denoised.mt_adddiff(ssDD,U=2,V=2) # Apply the limited difference. (Sharpening is just inverse blurring) | |
564 | HD ? mt_clamp(last,pmax,pmin,overshoot,overshoot,chroma="copy first") : last | |
565 | yuy26 ? last.ConvertToYUY2().Interleaved2planar() : last | |
566 | } | |
567 | ||
568 | # MinBlur by Didée (http://avisynth.org/mediawiki/MinBlur) | |
569 | # Nifty Gauss/Median combination (Modified for Planar input) | |
570 | ||
571 | function MinBlur(clip clp, int "r", int "uv", bool "planar"){ | |
572 | ||
573 | r = default(r,1) | |
574 | uv = default(uv,3) | |
575 | planar = default(planar,false) | |
576 | avs26 = VersionNumber() < 2.60 ? false : true | |
577 | yuy26 = avs26 && Isyuy2(clp) && planar | |
578 | ||
579 | uv2 = (uv==2) ? 1 : uv | |
580 | rg4 = (uv==3) ? 4 : -1 | |
581 | rg11 = (uv==3) ? 11 : -1 | |
582 | rg20 = (uv==3) ? 20 : -1 | |
583 | uvm2 = r==2 ? (uv==3?2:uv==2?0:-1) : nop() | |
584 | uvm3 = r==3 ? (uv==3?3:uv==2?0:-1) : nop() | |
585 | medf = (uv==3) ? 1 : -200 | |
586 | ||
587 | RG11D = (r==0) ? yuy26 ? mt_makediff(clp.Planar2Interleaved().ConvertToYV16(),clp.sbr(planar=planar).Planar2Interleaved().ConvertToYV16(),U=uv2,V=uv2).ConvertToYUY2().Interleaved2planar() : mt_makediff(clp,clp.sbr(planar=planar),U=uv2,V=uv2) | |
588 | \ : (r==1) ? yuy26 ? mt_makediff(clp.Planar2Interleaved().ConvertToYV16(),clp.removegrain(11,rg11,planar=planar).Planar2Interleaved().ConvertToYV16(),U=uv2,V=uv2).ConvertToYUY2().Interleaved2planar() : mt_makediff(clp,clp.removegrain(11,rg11,planar=planar),U=uv2,V=uv2) | |
589 | \ : (r==2) ? yuy26 ? mt_makediff(clp.Planar2Interleaved().ConvertToYV16(),clp.removegrain(11,rg11,planar=planar).removegrain(20,rg20,planar=planar).Planar2Interleaved().ConvertToYV16(),U=uv2,V=uv2).ConvertToYUY2().Interleaved2planar() : mt_makediff(clp,clp.removegrain(11,rg11,planar=planar).removegrain(20,rg20,planar=planar),U=uv2,V=uv2) | |
590 | \ : yuy26 ? mt_makediff(clp.Planar2Interleaved().ConvertToYV16(),clp.removegrain(11,rg11,planar=planar).removegrain(20,rg20,planar=planar).removegrain(20,rg20,planar=planar).Planar2Interleaved().ConvertToYV16(),U=uv2,V=uv2).ConvertToYUY2().Interleaved2planar() : mt_makediff(clp,clp.removegrain(11,rg11,planar=planar).removegrain(20,rg20,planar=planar).removegrain(20,rg20,planar=planar),U=uv2,V=uv2) | |
591 | - | \ : (r==2) ? avs26 ? yuy26 ? mt_makediff(clp.Planar2Interleaved().ConvertToYV16(),clp.Planar2Interleaved().ConvertToYV16().medianblur(2,2*medf,2*medf),U=uv2,V=uv2).ConvertToYUY2().Interleaved2planar() : mt_makediff(clp,clp.medianblur(2,2*medf,2*medf),U=uv2,V=uv2) : mt_makediff(clp,clp.Quantile(radius_y=2, radius_u=uvm2, radius_v=uvm2 ,planar=planar), U=uv2,V=uv2) |
591 | + | |
592 | \ : (r==2) ? avs26 ? yuy26 ? mt_makediff(clp.Planar2Interleaved().ConvertToYV16(),clp.Planar2Interleaved().ConvertToYV16().medianblur(2,2*medf,2*medf),U=uv2,V=uv2).ConvertToYUY2().Interleaved2planar() : mt_makediff(clp,clp.medianblur(2,2*medf,2*medf),U=uv2,V=uv2) : mt_makediff(clp,clp.Quantile(radius_y=2, radius_u=uvm2, radius_v=uvm2, planar=planar), U=uv2,V=uv2) | |
593 | \ : avs26 ? yuy26 ? mt_makediff(clp.Planar2Interleaved().ConvertToYV16(),clp.Planar2Interleaved().ConvertToYV16().medianblur(3,3*medf,3*medf),U=uv2,V=uv2).ConvertToYUY2().Interleaved2planar() : mt_makediff(clp,clp.medianblur(3,3*medf,3*medf),U=uv2,V=uv2) : mt_makediff(clp,clp.Quantile(radius_y=3, radius_u=uvm3, radius_v=uvm3, planar=planar), U=uv2,V=uv2) | |
594 | DD = yuy26 ? mt_lutxy(RG11D.Planar2Interleaved().ConvertToYV16(),RG4D.Planar2Interleaved().ConvertToYV16(),"x 128 - y 128 - * 0 < 128 x 128 - abs y 128 - abs < x y ? ?",U=uv2,V=uv2).ConvertToYUY2().Interleaved2planar() : mt_lutxy(RG11D,RG4D,"x 128 - y 128 - * 0 < 128 x 128 - abs y 128 - abs < x y ? ?",U=uv2,V=uv2) | |
595 | return (yuy26 ? clp.Planar2Interleaved().ConvertToYV16().mt_makediff(DD.Planar2Interleaved().ConvertToYV16(),U=uv,V=uv).ConvertToYUY2().Interleaved2planar() : clp.mt_makediff(DD,U=uv,V=uv)) } | |
596 | ||
597 | ||
598 | # helper function: make a highpass on a blur's difference (well, kind of that) | |
599 | # Highpass of spatial r=1 Gaussian (like blur(1.0)) (Modified for Planar input) | |
600 | ||
601 | function sbr(clip c, bool "planar") { | |
602 | ||
603 | planar = default(planar,false) | |
604 | avs26 = VersionNumber() < 2.60 ? false : true | |
605 | yuy26 = avs26 && Isyuy2(c) && planar | |
606 | ||
607 | rg11 = c.removegrain(11,-1,planar=planar) | |
608 | rg11D = yuy26 ? mt_makediff(c.Planar2Interleaved().ConvertToYV16(),rg11.Planar2Interleaved().ConvertToYV16()).ConvertToYUY2().Interleaved2planar() : mt_makediff(c,rg11) | |
609 | rg11DD = yuy26 ? mt_makediff(rg11D.Planar2Interleaved().ConvertToYV16(),rg11D.removegrain(11,-1,planar=planar).Planar2Interleaved().ConvertToYV16()).mt_lutxy(rg11D.Planar2Interleaved().ConvertToYV16(),"x 128 - y 128 - * 0 < 128 x 128 - abs y 128 - abs < x y ? ?").ConvertToYUY2().Interleaved2planar() : mt_makediff(rg11D,rg11D.removegrain(11,-1,planar=planar)).mt_lutxy(rg11D,"x 128 - y 128 - * 0 < 128 x 128 - abs y 128 - abs < x y ? ?") | |
610 | yuy26 ? c.Planar2Interleaved().ConvertToYV16().mt_makediff(rg11DD.Planar2Interleaved().ConvertToYV16(),U=2,V=2).ConvertToYUY2().Interleaved2planar() : c.mt_makediff(rg11DD,U=2,V=2) } | |
611 | ||
612 | ######################################## | |
613 | ## cretindesalpes' functions: (http://forum.doom9.org/showthread.php?p=1505297#post1505297) | |
614 | ||
615 | # Converts 16 bit interleaved YUY2 to 16bit planar YUY2 suitable for Ditherpost. (Use right before Ditherpost; this format is not suitable for filtering) | |
616 | ||
617 | function Dither_YUY2toPlanar16(c){c | |
618 | yuy26 = !(VersionNumber() < 2.60) && isyuy2(c) | |
619 | my = !yuy26 ? Dither_get_msb().ConvertToYV12() : nop() | |
620 | ly = !yuy26 ? Dither_get_lsb().ConvertToYV12() : nop() | |
621 | u = !yuy26 ? UToY ().ConvertToYV12 () : nop() | |
622 | v = !yuy26 ? VToY ().ConvertToYV12 () : nop() | |
623 | yuy26 ? ConvertToYV16() : YToUV (u, v, StackVertical (StackVertical (my, my), StackVertical (ly, ly)))} | |
624 | ||
625 | ||
626 | # Converts 8 bit interleaved YUY2 to 8 bit planar YUY2 ready for either 8 or 16 bit processing also suitable for Ditherpost | |
627 | ||
628 | function Dither_YUY2toPlanar(c){c | |
629 | yuy26 = !(VersionNumber() < 2.60) && isyuy2(c) | |
630 | y = !yuy26 ? ConvertToYV12 () : nop() | |
631 | u = !yuy26 ? UToY ().ConvertToYV12 () : nop() | |
632 | v = !yuy26 ? VToY ().ConvertToYV12 () : nop() | |
633 | yuy26 ? ConvertToYV16() : YToUV (u, v, StackVertical (y, y))} | |
634 | ||
635 | ||
636 | # Converts back the above special planar formats to normal interleaved 16 bit YUY2 | |
637 | ||
638 | function Dither_YUY2toInterleaved16(c){c | |
639 | yuy26 = !(VersionNumber() < 2.60) | |
640 | y = !yuy26 ? StackVertical (Crop (0, 0, 0, (Height () /(4))), Crop (0, (Height () /(2))+(Height () /(4)), 0, 0)).ConvertToYUY2 () : nop() | |
641 | u = !yuy26 ? UToY ().ConvertToYUY2 () : nop() | |
642 | v = !yuy26 ? VToY ().ConvertToYUY2 () : nop() | |
643 | yuy26 ? ConvertToYUY2() : YToUV (u, v, y)} | |
644 | ||
645 | # Converts back the above special planar formats to normal interleaved YUY2 | |
646 | ||
647 | function Dither_YUY2toInterleaved(c){c | |
648 | yuy26 = !(VersionNumber() < 2.60) | |
649 | y = !yuy26 ? Crop (0, 0, 0, Height () /(2)).ConvertToYUY2 () : nop() | |
650 | u = !yuy26 ? UToY ().ConvertToYUY2 () : nop() | |
651 | v = !yuy26 ? VToY ().ConvertToYUY2 () : nop() | |
652 | yuy26 ? ConvertToYUY2() : YToUV (u, v, y)} | |
653 | ||
654 | ||
655 | ||
656 | # Converts luma (and chroma) to PC levels, and optionally allows tweaking for pumping up the darks. (for the clip to be fed to motion search only) | |
657 | # By courtesy of cretindesalpes. (http://forum.doom9.org/showthread.php?p=1548318#post1548318) | |
658 | ||
659 | function Dither_Luma_Rebuild (clip src, float "s0", float "c",int "uv", bool "lsb", bool "lsb_in", bool "lsb_out", int "mode", float "ampn", bool "slice"){ | |
660 | lsb_in = Default( lsb_in ,false) | |
661 | lsb_out = Default( lsb_out ,false) | |
662 | lsb = Default( lsb ,lsb_in || lsb_out) | |
663 | mode = Default( mode ,6) | |
664 | uv = Default(uv, 3) | |
665 | s0 = Default(s0, 2.0) | |
666 | c = Default(c, 1.0/16) | |
667 | ampn = Default(ampn, 0.0) | |
668 | slice = Default(slice, true) | |
669 | yuy26 = !(VersionNumber() < 2.60) && isyuy2(src) | |
670 | src = yuy26 ? src.Planar2Interleaved().ConvertToYV16() : src | |
671 | ||
672 | k = (s0 - 1) * c | |
673 | t = lsb_in ? "x 4096 - 56064 / 0 1 clip" : "x 16 - 219 / 0 1 clip" | |
674 | e = String(k)+" "+String(1+c)+" "+String((1+c)*c)+" "+t+" "+String(c) | |
675 | \ +" + / - * "+t+" 1 "+String(k)+" - * + "+String(lsb?65536:255)+" *" | |
676 | src | |
677 | lsb ? (lsb_in ? Dither_lut16 (yexpr=e,expr="x 32768 - 32768 * 28672 / 32768 +",y=3, u=uv, v=uv) : \ | |
678 | Dither_lut8 (yexpr=e,expr="x 128 - 32768 * 112 / 32768 +" ,y=3, u=uv, v=uv)) : \ | |
679 | mt_lut (yexpr=e,expr="x 128 - 128 * 112 / 128 +" ,y=3, u=uv, v=uv) | |
680 | lsb_out ? last : (lsb ? Ditherpost(mode=mode,ampn=ampn,slice=slice,staticnoise=ampn!=0.0?true:false,u=uv,v=uv) : last) | |
681 | return yuy26 ? last.ConvertToYUY2().Interleaved2planar() : last | |
682 | } |