|
Daniel@0
|
1
|
|
Daniel@0
|
2 <!DOCTYPE html
|
|
Daniel@0
|
3 PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
|
|
Daniel@0
|
4 <html><head>
|
|
Daniel@0
|
5 <meta http-equiv="Content-Type" content="text/html; charset=utf-8">
|
|
Daniel@0
|
6 <!--
|
|
Daniel@0
|
7 This HTML is auto-generated from an M-file.
|
|
Daniel@0
|
8 To make changes, update the M-file and republish this document.
|
|
Daniel@0
|
9 --><title>MP3 reading and writing</title><meta name="generator" content="MATLAB 7.10"><meta name="date" content="2010-04-09"><meta name="m-file" content="demo_mp3readwrite"><style type="text/css">
|
|
Daniel@0
|
10
|
|
Daniel@0
|
11 body {
|
|
Daniel@0
|
12 background-color: white;
|
|
Daniel@0
|
13 margin:10px;
|
|
Daniel@0
|
14 }
|
|
Daniel@0
|
15
|
|
Daniel@0
|
16 h1 {
|
|
Daniel@0
|
17 color: #990000;
|
|
Daniel@0
|
18 font-size: x-large;
|
|
Daniel@0
|
19 }
|
|
Daniel@0
|
20
|
|
Daniel@0
|
21 h2 {
|
|
Daniel@0
|
22 color: #990000;
|
|
Daniel@0
|
23 font-size: medium;
|
|
Daniel@0
|
24 }
|
|
Daniel@0
|
25
|
|
Daniel@0
|
26 /* Make the text shrink to fit narrow windows, but not stretch too far in
|
|
Daniel@0
|
27 wide windows. */
|
|
Daniel@0
|
28 p,h1,h2,div.content div {
|
|
Daniel@0
|
29 max-width: 600px;
|
|
Daniel@0
|
30 /* Hack for IE6 */
|
|
Daniel@0
|
31 width: auto !important; width: 600px;
|
|
Daniel@0
|
32 }
|
|
Daniel@0
|
33
|
|
Daniel@0
|
34 pre.codeinput {
|
|
Daniel@0
|
35 background: #EEEEEE;
|
|
Daniel@0
|
36 padding: 10px;
|
|
Daniel@0
|
37 }
|
|
Daniel@0
|
38 @media print {
|
|
Daniel@0
|
39 pre.codeinput {word-wrap:break-word; width:100%;}
|
|
Daniel@0
|
40 }
|
|
Daniel@0
|
41
|
|
Daniel@0
|
42 span.keyword {color: #0000FF}
|
|
Daniel@0
|
43 span.comment {color: #228B22}
|
|
Daniel@0
|
44 span.string {color: #A020F0}
|
|
Daniel@0
|
45 span.untermstring {color: #B20000}
|
|
Daniel@0
|
46 span.syscmd {color: #B28C00}
|
|
Daniel@0
|
47
|
|
Daniel@0
|
48 pre.codeoutput {
|
|
Daniel@0
|
49 color: #666666;
|
|
Daniel@0
|
50 padding: 10px;
|
|
Daniel@0
|
51 }
|
|
Daniel@0
|
52
|
|
Daniel@0
|
53 pre.error {
|
|
Daniel@0
|
54 color: red;
|
|
Daniel@0
|
55 }
|
|
Daniel@0
|
56
|
|
Daniel@0
|
57 p.footer {
|
|
Daniel@0
|
58 text-align: right;
|
|
Daniel@0
|
59 font-size: xx-small;
|
|
Daniel@0
|
60 font-weight: lighter;
|
|
Daniel@0
|
61 font-style: italic;
|
|
Daniel@0
|
62 color: gray;
|
|
Daniel@0
|
63 }
|
|
Daniel@0
|
64
|
|
Daniel@0
|
65 </style></head><body><div class="content"><h1>MP3 reading and writing</h1><!--introduction--><p>These function, mp3read and mp3write, aim to exactly duplicate the operation of wavread and wavwrite for accessing soundfiles, except the soundfiles are in Mpeg-Audio layer 3 (MP3) compressed format. All the hard work is done by external binaries written by others: mp3info to query the format of existing mp3 files, mpg123 to decode mp3 files, and lame to encode audio files. Binaries for these files are widely available (and may be included in this distribution).</p><p>These functions were originally developed for access to very large mp3 files (i.e. many hours long), and so avoid creating the entire uncompressed audio stream if possible. mp3read allows you to specify the range of frames you want to read (as a second argument), and mp3read will construct an mpg123 command that skips blocks to decode only the part of the file that is required. This can be much quicker (and require less memory/temporary disk) than decoding the whole file.</p><p>mpg123 also provides for "on the fly" downsampling at conversion to mono, which are supported as extra options in mp3read.</p><p>mpg123 can read MP3s across the network. This is supported if the FILE argument is a URL (e.g. beginning 'http://...').</p><p>mp3info sometimes gets the file size wrong (as returned by the mp3read(...'size') syntax). I'm not sure when this happens exactly, but it's probably a result of VBR files. In the worst case, figuring the number of samples in such a file requires scanning through the whole file, and mp3info doesn't usually do this.</p><p>For more information, including advice on handling MP4 files, see <a href="http://labrosa.ee.columbia.edu/matlab/mp3read.html">http://labrosa.ee.columbia.edu/matlab/mp3read.html</a></p><!--/introduction--><h2>Contents</h2><div><ul><li><a href="#1">Example usage</a></li><li><a href="#2">Delay, size, and alignment</a></li><li><a href="#3">External binaries</a></li><li><a href="#4">Installation</a></li></ul></div><h2>Example usage<a name="1"></a></h2><p>Here, we read a wav file in, then write it out as an MP3, then read the resulting MP3 back in, and compare it to the original file.</p><pre class="codeinput"><span class="comment">% Read an audio waveform</span>
|
|
Daniel@0
|
66 [d,sr] = wavread(<span class="string">'piano.wav'</span>);
|
|
Daniel@0
|
67 <span class="comment">% Save to mp3 (default settings)</span>
|
|
Daniel@0
|
68 mp3write(d,sr,<span class="string">'piano.mp3'</span>);
|
|
Daniel@0
|
69 <span class="comment">% Read it back again</span>
|
|
Daniel@0
|
70 [d2,sr] = mp3read(<span class="string">'piano.mp3'</span>);
|
|
Daniel@0
|
71 <span class="comment">% mp3 encoding involves some extra padding at each end; we attempt</span>
|
|
Daniel@0
|
72 <span class="comment">% to cut it off at the start, but can't do that at the end, because</span>
|
|
Daniel@0
|
73 <span class="comment">% mp3read doesn't know how long the original was. But we do, so..</span>
|
|
Daniel@0
|
74 <span class="comment">% Chop it down to be the same length as the original</span>
|
|
Daniel@0
|
75 d2 = d2(1:length(d),:);
|
|
Daniel@0
|
76 <span class="comment">% What is the SNR (distortion)?</span>
|
|
Daniel@0
|
77 ddiff = d - d2;
|
|
Daniel@0
|
78 disp([<span class="string">'SNR is '</span>,num2str(10*log10(sum(d(:).^2)/sum(ddiff(:).^2))),<span class="string">' dB'</span>]);
|
|
Daniel@0
|
79 <span class="comment">% Do they look similar?</span>
|
|
Daniel@0
|
80 subplot(211)
|
|
Daniel@0
|
81 specgram(d(:,1),1024,sr);
|
|
Daniel@0
|
82 subplot(212)
|
|
Daniel@0
|
83 plot(1:5000,d(10000+(1:5000),1),1:5000,d2(10000+(1:5000)));
|
|
Daniel@0
|
84 <span class="comment">% Yes, pretty close</span>
|
|
Daniel@0
|
85 <span class="comment">%</span>
|
|
Daniel@0
|
86 <span class="comment">% NB: lame followed by mpg123 causes a little attenuation; you</span>
|
|
Daniel@0
|
87 <span class="comment">% can get a better match by scaling up the read-back waveform:</span>
|
|
Daniel@0
|
88 ddiff = d - 1.052*d2;
|
|
Daniel@0
|
89 disp([<span class="string">'SNR is '</span>,num2str(10*log10(sum(d(:).^2)/sum(ddiff(:).^2))),<span class="string">' dB'</span>]);
|
|
Daniel@0
|
90 </pre><pre class="codeoutput">Warning: popenw not available, writing temporary file
|
|
Daniel@0
|
91 SNR is 22.632 dB
|
|
Daniel@0
|
92 SNR is 24.8699 dB
|
|
Daniel@0
|
93 </pre><img vspace="5" hspace="5" src="demo_mp3readwrite_01.png" alt=""> <h2>Delay, size, and alignment<a name="2"></a></h2><p>In mid-2006 I noticed that mp3read followed by mp3write followed by mp3read effectively delayed the waveform by 2257 samples (at 44 kHz). So I introduced code to discard the first 2257 samples to ensure that the waveforms remained time aligned. As best I could understand, mpg123 (v 0.5.9) was including the "warm-up" samples from the synthesis filterbank which are more properly discarded.</p><p>Then in late 2009 I noticed that some chord recognition code, which used mp3read to read files which were then segmented on the basis of some hand-marked timings, suddenly started getting much poorer results. It turned out that I had upgraded my version of mpg123 to v 1.9.0, and the warm-up samples had been fixed in this version. So my code was discarding 2257 <b>good</b> samples, and the data was skewed 51ms early relative to the hand labels.</p><p>Hence, the current version of mp3read does not discard any samples by default -- appropriate for the recent versions of mpg123 included here. But if you know you're running an old, v 0.5.9, mpg123, you should edit the mp3read.m source to set the flag MPG123059 = 1.</p><p>Note also that the 'size' function relies on the number of blocks reported by mp3info. However, many mp3 files include additional information about the size of the file in the so-called Xing header, embedded in the first frame, which can specify that a certain number of samples from start and end should additionally be dropped. mp3info doesn't read that, and there's no way for my code to probe it except by running mpg123. Hence, the results of mp3read(fn,'size') may sometimes overestimate the length of the actual vector you'll get if you read the whole file.</p><h2>External binaries<a name="3"></a></h2><p>The m files rely on three external binaries, each of which is available for Linux, Mac OS X, or Windows:</p><p><b>mpg123</b> is a high-performance mp3 decoder. Its home page is <a href="http://www.mpg123.de/">http://www.mpg123.de/</a> .</p><p><b>mp3info</b> is a utility to read technical information on an mp3 file. Its home page is <a href="http://www.ibiblio.org/mp3info/">http://www.ibiblio.org/mp3info/</a> .</p><p><b>lame</b> is an open-source MP3 encoder. Its homepage is <a href="http://lame.sourceforge.net/">http://lame.sourceforge.net/</a> .</p><p>The various authors of these packages are gratefully acknowledged for doing all the hard work to make these Matlab functions possible.</p><h2>Installation<a name="4"></a></h2><p>The two routines, mp3read.m and mp3write.m, will look for their binaries (mpg123 and mp3info for mp3read; lame for mp3write) in the same directory where they are installed. Binaries for different architectures are distinguished by their extension, which is the standard Matlab computer code e.g. ".mac" for Mac PPC OS X, ".glnx86" for i386-linux. The exception is Windows, where the binaries have the extension ".exe".</p><p>Temporary files will be written to (a) a directory taken from the environment variable TMPDIR (b) /tmp if it exists, or (c) the current directory. This can easily be changed by editing the m files.</p><pre class="codeinput"><span class="comment">% Last updated: $Date: 2009/03/15 18:29:58 $</span>
|
|
Daniel@0
|
94 <span class="comment">% Dan Ellis <dpwe@ee.columbia.edu></span>
|
|
Daniel@0
|
95 </pre><p class="footer"><br>
|
|
Daniel@0
|
96 Published with MATLAB® 7.10<br></p></div><!--
|
|
Daniel@0
|
97 ##### SOURCE BEGIN #####
|
|
Daniel@0
|
98 %% MP3 reading and writing
|
|
Daniel@0
|
99 %
|
|
Daniel@0
|
100 % These function, mp3read and mp3write, aim to exactly duplicate
|
|
Daniel@0
|
101 % the operation of wavread and wavwrite for accessing soundfiles,
|
|
Daniel@0
|
102 % except the soundfiles are in Mpeg-Audio layer 3 (MP3) compressed
|
|
Daniel@0
|
103 % format. All the hard work is done by external binaries written
|
|
Daniel@0
|
104 % by others: mp3info to query the format of existing mp3 files,
|
|
Daniel@0
|
105 % mpg123 to decode mp3 files, and lame to encode audio files.
|
|
Daniel@0
|
106 % Binaries for these files are widely available (and may be
|
|
Daniel@0
|
107 % included in this distribution).
|
|
Daniel@0
|
108 %
|
|
Daniel@0
|
109 % These functions were originally developed for access to very
|
|
Daniel@0
|
110 % large mp3 files (i.e. many hours long), and so avoid creating
|
|
Daniel@0
|
111 % the entire uncompressed audio stream if possible. mp3read
|
|
Daniel@0
|
112 % allows you to specify the range of frames you want to read
|
|
Daniel@0
|
113 % (as a second argument), and mp3read will construct an mpg123
|
|
Daniel@0
|
114 % command that skips blocks to decode only the part of the file
|
|
Daniel@0
|
115 % that is required. This can be much quicker (and require less
|
|
Daniel@0
|
116 % memory/temporary disk) than decoding the whole file.
|
|
Daniel@0
|
117 %
|
|
Daniel@0
|
118 % mpg123 also provides for "on the fly" downsampling at conversion
|
|
Daniel@0
|
119 % to mono, which are supported as extra options in mp3read.
|
|
Daniel@0
|
120 %
|
|
Daniel@0
|
121 % mpg123 can read MP3s across the network. This is supported
|
|
Daniel@0
|
122 % if the FILE argument is a URL (e.g. beginning 'http://...').
|
|
Daniel@0
|
123 %
|
|
Daniel@0
|
124 % mp3info sometimes gets the file size wrong (as returned by the
|
|
Daniel@0
|
125 % mp3read(...'size') syntax). I'm not sure when this happens
|
|
Daniel@0
|
126 % exactly, but it's probably a result of VBR files. In the worst
|
|
Daniel@0
|
127 % case, figuring the number of samples in such a file requires
|
|
Daniel@0
|
128 % scanning through the whole file, and mp3info doesn't usually do
|
|
Daniel@0
|
129 % this.
|
|
Daniel@0
|
130 %
|
|
Daniel@0
|
131 % For more information, including advice on handling MP4 files,
|
|
Daniel@0
|
132 % see http://labrosa.ee.columbia.edu/matlab/mp3read.html
|
|
Daniel@0
|
133
|
|
Daniel@0
|
134 %% Example usage
|
|
Daniel@0
|
135 % Here, we read a wav file in, then write it out as an MP3, then
|
|
Daniel@0
|
136 % read the resulting MP3 back in, and compare it to the original
|
|
Daniel@0
|
137 % file.
|
|
Daniel@0
|
138
|
|
Daniel@0
|
139 % Read an audio waveform
|
|
Daniel@0
|
140 [d,sr] = wavread('piano.wav');
|
|
Daniel@0
|
141 % Save to mp3 (default settings)
|
|
Daniel@0
|
142 mp3write(d,sr,'piano.mp3');
|
|
Daniel@0
|
143 % Read it back again
|
|
Daniel@0
|
144 [d2,sr] = mp3read('piano.mp3');
|
|
Daniel@0
|
145 % mp3 encoding involves some extra padding at each end; we attempt
|
|
Daniel@0
|
146 % to cut it off at the start, but can't do that at the end, because
|
|
Daniel@0
|
147 % mp3read doesn't know how long the original was. But we do, so..
|
|
Daniel@0
|
148 % Chop it down to be the same length as the original
|
|
Daniel@0
|
149 d2 = d2(1:length(d),:);
|
|
Daniel@0
|
150 % What is the SNR (distortion)?
|
|
Daniel@0
|
151 ddiff = d - d2;
|
|
Daniel@0
|
152 disp(['SNR is ',num2str(10*log10(sum(d(:).^2)/sum(ddiff(:).^2))),' dB']);
|
|
Daniel@0
|
153 % Do they look similar?
|
|
Daniel@0
|
154 subplot(211)
|
|
Daniel@0
|
155 specgram(d(:,1),1024,sr);
|
|
Daniel@0
|
156 subplot(212)
|
|
Daniel@0
|
157 plot(1:5000,d(10000+(1:5000),1),1:5000,d2(10000+(1:5000)));
|
|
Daniel@0
|
158 % Yes, pretty close
|
|
Daniel@0
|
159 %
|
|
Daniel@0
|
160 % NB: lame followed by mpg123 causes a little attenuation; you
|
|
Daniel@0
|
161 % can get a better match by scaling up the read-back waveform:
|
|
Daniel@0
|
162 ddiff = d - 1.052*d2;
|
|
Daniel@0
|
163 disp(['SNR is ',num2str(10*log10(sum(d(:).^2)/sum(ddiff(:).^2))),' dB']);
|
|
Daniel@0
|
164
|
|
Daniel@0
|
165 %% Delay, size, and alignment
|
|
Daniel@0
|
166 %
|
|
Daniel@0
|
167 % In mid-2006 I noticed that mp3read followed by mp3write followed by
|
|
Daniel@0
|
168 % mp3read effectively delayed the waveform by 2257 samples (at 44
|
|
Daniel@0
|
169 % kHz). So I introduced code to discard the first 2257 samples to ensure
|
|
Daniel@0
|
170 % that the waveforms remained time aligned. As best I could understand,
|
|
Daniel@0
|
171 % mpg123 (v 0.5.9) was including the "warm-up" samples from the
|
|
Daniel@0
|
172 % synthesis filterbank which are more properly discarded.
|
|
Daniel@0
|
173 %
|
|
Daniel@0
|
174 % Then in late 2009 I noticed that some chord recognition code, which
|
|
Daniel@0
|
175 % used mp3read to read files which were then segmented on the basis of
|
|
Daniel@0
|
176 % some hand-marked timings, suddenly started getting much poorer
|
|
Daniel@0
|
177 % results. It turned out that I had upgraded my version of mpg123 to v
|
|
Daniel@0
|
178 % 1.9.0, and the warm-up samples had been fixed in this version. So my
|
|
Daniel@0
|
179 % code was discarding 2257 *good* samples, and the data was skewed 51ms
|
|
Daniel@0
|
180 % early relative to the hand labels.
|
|
Daniel@0
|
181 %
|
|
Daniel@0
|
182 % Hence, the current version of mp3read does not
|
|
Daniel@0
|
183 % discard any samples by default REPLACE_WITH_DASH_DASH appropriate for the recent versions
|
|
Daniel@0
|
184 % of mpg123 included here. But if you know you're running an old, v
|
|
Daniel@0
|
185 % 0.5.9, mpg123, you should edit the mp3read.m source to set the flag
|
|
Daniel@0
|
186 % MPG123059 = 1.
|
|
Daniel@0
|
187 %
|
|
Daniel@0
|
188 % Note also that the 'size' function relies on the number of
|
|
Daniel@0
|
189 % blocks reported by mp3info. However, many mp3 files include
|
|
Daniel@0
|
190 % additional information about the size of the file in the
|
|
Daniel@0
|
191 % so-called Xing header, embedded in the first frame, which can
|
|
Daniel@0
|
192 % specify that a certain number of samples from start and end
|
|
Daniel@0
|
193 % should additionally be dropped. mp3info doesn't read that,
|
|
Daniel@0
|
194 % and there's no way for my code to probe it except by running
|
|
Daniel@0
|
195 % mpg123. Hence, the results of mp3read(fn,'size') may sometimes
|
|
Daniel@0
|
196 % overestimate the length of the actual vector you'll get if
|
|
Daniel@0
|
197 % you read the whole file.
|
|
Daniel@0
|
198
|
|
Daniel@0
|
199 %% External binaries
|
|
Daniel@0
|
200 % The m files rely on three external binaries, each of which is
|
|
Daniel@0
|
201 % available for Linux, Mac OS X, or Windows:
|
|
Daniel@0
|
202 %
|
|
Daniel@0
|
203 % *mpg123* is a high-performance mp3 decoder. Its home page is
|
|
Daniel@0
|
204 % http://www.mpg123.de/ .
|
|
Daniel@0
|
205 %
|
|
Daniel@0
|
206 % *mp3info* is a utility to read technical information on an mp3
|
|
Daniel@0
|
207 % file. Its home page is http://www.ibiblio.org/mp3info/ .
|
|
Daniel@0
|
208 %
|
|
Daniel@0
|
209 % *lame* is an open-source MP3 encoder. Its homepage is
|
|
Daniel@0
|
210 % http://lame.sourceforge.net/ .
|
|
Daniel@0
|
211 %
|
|
Daniel@0
|
212 % The various authors of these packages are gratefully acknowledged
|
|
Daniel@0
|
213 % for doing all the hard work to make these Matlab functions possible.
|
|
Daniel@0
|
214
|
|
Daniel@0
|
215 %% Installation
|
|
Daniel@0
|
216 % The two routines, mp3read.m and mp3write.m, will look for their
|
|
Daniel@0
|
217 % binaries (mpg123 and mp3info for mp3read; lame for mp3write) in
|
|
Daniel@0
|
218 % the same directory where they are installed. Binaries for
|
|
Daniel@0
|
219 % different architectures are distinguished by their extension,
|
|
Daniel@0
|
220 % which is the standard Matlab computer code e.g. ".mac" for Mac
|
|
Daniel@0
|
221 % PPC OS X, ".glnx86" for i386-linux. The exception is Windows,
|
|
Daniel@0
|
222 % where the binaries have the extension ".exe".
|
|
Daniel@0
|
223 %
|
|
Daniel@0
|
224 % Temporary files
|
|
Daniel@0
|
225 % will be written to (a) a directory taken from the environment
|
|
Daniel@0
|
226 % variable TMPDIR (b) /tmp if it exists, or (c) the current
|
|
Daniel@0
|
227 % directory. This can easily be changed by editing the m files.
|
|
Daniel@0
|
228
|
|
Daniel@0
|
229 % Last updated: $Date: 2009/03/15 18:29:58 $
|
|
Daniel@0
|
230 % Dan Ellis <dpwe@ee.columbia.edu>
|
|
Daniel@0
|
231
|
|
Daniel@0
|
232 ##### SOURCE END #####
|
|
Daniel@0
|
233 --></body></html> |
