forked from GreyElaina/SingingPiano
-
Notifications
You must be signed in to change notification settings - Fork 2
/
etuconvert.py
55 lines (54 loc) · 1.95 KB
/
etuconvert.py
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
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
import libs
import wave
from numpy import array
import struct
from tqdm import trange
def decodingDFTData(DFTData):
"""version=0.1"""
DFTList = []
if not len(DFTData)%112 == 0:
print("Data corrupted .")
DDList = [DFTData[k*112:(k+1)*112] for k in range(int(len(DFTData)/112))]
for i in DDList:
L_i=[]
for j in range(128):
if j<16:#filtering low-band noise (f<20hz), used for spectrogram data obtained by STFT with ticklength=0.05 and windowlength=0.1
L_i.append(0)
continue
pos=(j+1)*7//8
offset=j%8
#n__1=N_i[j-1]
n__1 =i[pos-1]
n_0 =i[pos] if len(i) > pos+1 else 0
L_i.append(((n__1 >> (8-offset)) + (n_0 << offset))%128)
DFTList.append(L_i)
return DFTList
if __name__ == "__main__":
import sys
_name=sys.argv[1]#input("Drag etu file into this window:")
if _name[0]=='"' and _name[-1] == '"':
_name=_name[1:-1]
name=_name[:-4]
ticklength=0.05
with open(name+".etu","rb") as fetu:
X=decodingDFTData(fetu.read())
for i in range(len(X)):#trange(len(NFTData),dynamic_ncols=True,ascii=True,smoothing=1,mininterval=0.25,unit="ticks",unit_scale=False,unit_divisor=1024):
for k in range(len(X[0])):
this_n=(X[i][k]*32)**2
X[i][k]=this_n
#imi.append(this_n)
#if this_n >= 128:notesOverFlow+=1
T=len(X)*ticklength
fs = 44100
x=libs.myalgs.distft(array(X), fs, T, ticklength*2, ticklength)
w= array([round(j/64) for j in x],dtype="int16")
print(max(w),min(w),max(x),min(x))
#B=b""
#or i in trange(len(w),dynamic_ncols=True,ascii=True,smoothing=1,mininterval=0.25,unit="ticks",unit_scale=False,unit_divisor=1024):
#B+= struct.pack("h",int(w[i]))
W=wave.open(name+".etu.wav",mode='wb')
W.setframerate(fs)
W.setnchannels(1)
W.setsampwidth(2)
W.writeframes(w)
input("Complete.")