API tools faq
paste
Advertisement
Guest User

SMDegrain

a guest
Jun 23rd, 2016
830
0
Never
Add comment
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Avisynth 55.13 KB | None | 0 0
raw download clone embed print report diff
  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-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. cf1=original.MCompensate(Super, cf1, planar=planar)
  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. RG4D  = (r<=1) ? yuy26 ? mt_makediff(clp.Planar2Interleaved().ConvertToYV16(),clp.removegrain(4,rg4,planar=planar).Planar2Interleaved().ConvertToYV16(),U=uv2,V=uv2).ConvertToYUY2().Interleaved2planar() : mt_makediff(clp,clp.removegrain(4,rg4,planar=planar),U=uv2,V=uv2)
  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. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement
Public Pastes
Advertisement
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.  OK, I Understand
Not a member of Pastebin yet?
Sign Up, it unlocks many cool features!