波形頻譜

saiyaku

嘛 ASS生成太大了 用處不大
只是結合這個 也比較好理解 頻譜是怎麼一回事
有興趣的可以寫一些 自己的頻譜
思路上建議不用拘泥于傳統頻譜 水波式的發散頻譜 或者是某種圖形 四周的點作為波形點 不規則拉升

1. from tcaxPy import *
   
2. 
   
3. from util.tcAudio import *
   
4. 
   
5. SPECHEIGHT = 127   #頻譜高度
   
6. SPECWIDTH = 300    #頻譜寬度（太寬了很卡 測試可以更小一點 到500 我就crash沒壓出來。。。）
   
7. BANDS = 300        #柱狀條數（波形頻譜和寬度一樣）
   
8. 
   
9. 
   
10. def Spectrum(fft):      # convert fft data to graph
    
11. 
    
12.     bars = []
    
13. 
    
14.     b0 = 0
    
15. 
    
16.     for x in range(BANDS):
    
17. 
    
18.         peak = 0        # peak of a certain bar
    
19. 
    
20.         b1 = pow(2, x * 10.0 / (BANDS - 1))
    
21. 
    
22.         if b1 > 1023:
    
23. 
    
24.             b1 = 1023
    
25. 
    
26.         if b1 <= b0:
    
27. 
    
28.             b1 = b0 + 1    # make sure it uses at least 1 FFT bin
    
29. 
    
30.         while b0 < b1:
    
31. 
    
32.             if peak < fft[1 + b0]:
    
33. 
    
34.                 peak = fft[1 + b0]
    
35. 
    
36.             b0 += 1
    
37. 
    
38.         y = int(sqrt(peak) * 3 * SPECHEIGHT - 4)   # scale it (sqrt to make low values more visible)
    
39. 
    
40.         if y > SPECHEIGHT:
    
41. 
    
42.             y = SPECHEIGHT     # cap it
    
43. 
    
44.         bars.append(y)
    
45. 
    
46.     return bars
    
47. 
    
48. def tcaxPy_User():
    
49. 
    
50. 
    
51.     file_name = GetVal(val_OutFile) + '.ass'
    
52. 
    
53.     ass_header = GetVal(val_AssHeader)
    
54. 
    
55.     ASS_FILE = CreateAssFile(file_name, ass_header)
    
56. 
    
57.     tcAudioInit()
    
58. 
    
59.     channel = tcAudioOpen('h.mp3')      # 此處修改 音頻文件名
    
60. 
    
61.     freq = tcAudioGetFreq(channel)      # usually 44100 Hz
    
62. 
    
63.     duration = tcAudioGetDuration(channel)
    
64. 
    
65.     num = int(duration * freq / 2048)
    
66. 
    
67.     hStr1 = "{\p1}m 0 0 l 0 -70 l 1 0 l 0 70 l 0 0"
    
68.     hStr2 = "{\p1}m 0 0 l 0 -50 l 1 0 l 0 50 l 0 0"
    
69. 
    
70.     for i in range(num):
    
71. 
    
72.         ASS_BUF = []
    
73. 
    
74. 
    
75.         fft = tcAudioGetFFT(channel, 2048)
    
76. 
    
77.         bars = Spectrum(fft)
    
78. 
    
79.         bar_width = int(SPECWIDTH / BANDS)
    
80. 
    
81.         count = len(bars)
    
82. 
    
83.         for j in range(count):
    
84. 
    
85.             start = i * 100 * 2048 / freq
    
86. 
    
87.             end = (i + 1) * 100 * 2048 / freq
    
88. 
    
89.             x = (1280 - SPECWIDTH) / 2 + bar_width * j
    
90. 
    
91.             y = 500
    
92. 
    
93.             FSC = fsc((bar_width) * 100, bars[j] * 2)
    
94. 
    
95.             CLR1 = color1("FFBD00") + alpha1(0)   
    
96. 
    
97.             CLR2 = color1("FFFFFF") + alpha1(119)   
    
98.             ass_main(ASS_BUF, SubL(start, end, 0, Pix_Style), an(7) +blur(0.8) + pos(x, y) + FSC + CLR1, hStr1)
    
99.             ass_main(ASS_BUF, SubL(start, end, 1, Pix_Style), an(7) +blur(0.8) + pos(x, y) + FSC + CLR2, hStr2)
    
100.         WriteAssFile(ASS_FILE, ASS_BUF)     # write the buffer in memory to the file
     
101. 
     
102.         progress(i + 1, num)
     
103. 
     
104.     tcAudioFin()
     
105. 
     
106.     FinAssFile(ASS_FILE)

复制代码