1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
|
#! /usr/bin/env python
# encoding: utf-8
"""
Spectral domain noise reduction.
"""
from pyo import *
s = Server(duplex=0).boot()
snd = "../snds/flute.aif"
chnls = sndinfo(snd)[3]
src = SfPlayer(snd, loop=True, mul=0.5)
size = 1024
olaps = 4
# bin-by-bin amplitude threshold
thresh = Sig(0.1)
thresh.ctrl([SLMap(0.0001, 1, "log", "value", 0.1)], title="Threshold")
# attenuation
mult = Sig(0.1)
mult.ctrl([SLMap(0.0001, 1, "log", "value", 0.1)], title="Attenuation")
fin = FFT(src, size=size, overlaps=olaps)
mag = Sqrt(fin["real"] * fin["real"] + fin["imag"] * fin["imag"])
amp = Compare(mag * 50, thresh, ">")
scl = amp * (1 - mult) + mult
re = fin["real"] * scl
im = fin["imag"] * scl
fout = IFFT(re, im, size=size, overlaps=olaps)
ffout = Mix(fout, chnls, mul=0.5).out()
s.gui(locals())
|
