octave-nkf: scripts/audio/wavread.m annotate

annotate scripts/audio/wavread.m @ 14271:e2a14d1b4eaa

rgbplot.m: Add new Matlab compatible function for plotting colormaps * NEWS: Add section for 3.8 release and list of new functions. * image.txi: Add rgbplot to documentation. * scripts/image/module.mk: Add rgbplot to build system. * rgbplot.m: Added new function.

author	Rik <octave@nomad.inbox5.com>
date	Fri, 27 Jan 2012 19:09:02 -0800
parents	72c96de7a403
children	f3d52523cde1

rev	line source
14138 72c96de7a403 maint: update copyright notices for 2012 John W. Eaton <jwe@octave.org> parents: 13003 diff changeset	1 ## Copyright (C) 2005-2012 Michael Zeising
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	2 ##
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	3 ## This file is part of Octave.
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	4 ##
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	5 ## Octave is free software; you can redistribute it and/or modify it
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	6 ## under the terms of the GNU General Public License as published by
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6549 diff changeset	7 ## the Free Software Foundation; either version 3 of the License, or (at
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6549 diff changeset	8 ## your option) any later version.
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	9 ##
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	10 ## Octave is distributed in the hope that it will be useful, but
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	11 ## WITHOUT ANY WARRANTY; without even the implied warranty of
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	13 ## General Public License for more details.
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	14 ##
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	15 ## You should have received a copy of the GNU General Public License
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6549 diff changeset	16 ## along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6549 diff changeset	17 ## <http://www.gnu.org/licenses/>.
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	18
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	19 ## -- texinfo --
6547 4fb053f24fd6 [project @ 2007-04-19 21:47:40 by jwe] jwe parents: 6304 diff changeset	20 ## @deftypefn {Function File} {@var{y} =} wavread (@var{filename})
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	21 ## Load the RIFF/WAVE sound file @var{filename}, and return the samples
80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	22 ## in vector @var{y}. If the file contains multichannel data, then
80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	23 ## @var{y} is a matrix with the channels represented as columns.
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	24 ##
12701 de3e90a420e3 Overhaul wavwrite, wavread and fix normalization problem (Bug #33420). Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	25 ## @deftypefnx {Function File} {[@var{y}, @var{Fs}, @var{bps}] =} wavread (@var{filename})
11587 c792872f8942 all script files: untabify and strip trailing whitespace John W. Eaton <jwe@octave.org> parents: 11563 diff changeset	26 ## Additionally return the sample rate (@var{fs}) in Hz and the number of bits
12701 de3e90a420e3 Overhaul wavwrite, wavread and fix normalization problem (Bug #33420). Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	27 ## per sample (@var{bps}).
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	28 ##
6547 4fb053f24fd6 [project @ 2007-04-19 21:47:40 by jwe] jwe parents: 6304 diff changeset	29 ## @deftypefnx {Function File} {[@dots{}] =} wavread (@var{filename}, @var{n})
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	30 ## Read only the first @var{n} samples from each channel.
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	31 ##
11563 3c6e8aaa9555 Grammarcheck m-files before 3.4 release. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	32 ## @deftypefnx {Function File} {[@dots{}] =} wavread (@var{filename}, @var{n1} @var{n2})
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	33 ## Read only samples @var{n1} through @var{n2} from each channel.
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	34 ##
6547 4fb053f24fd6 [project @ 2007-04-19 21:47:40 by jwe] jwe parents: 6304 diff changeset	35 ## @deftypefnx {Function File} {[@var{samples}, @var{channels}] =} wavread (@var{filename}, "size")
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	36 ## Return the number of samples (@var{n}) and channels (@var{ch})
80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	37 ## instead of the audio data.
5642 2618a0750ae6 [project @ 2006-03-06 21:26:48 by jwe] jwe parents: 5572 diff changeset	38 ## @seealso{wavwrite}
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	39 ## @end deftypefn
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	40
7117 373c39904877 [project @ 2007-11-07 18:48:00 by jwe] jwe parents: 7017 diff changeset	41 ## Author: Michael Zeising <michael@michaels-website.de>
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	42 ## Created: 06 December 2005
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	43
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	44 function [y, samples_per_sec, bits_per_sample] = wavread (filename, param)
80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	45
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	46 FORMAT_PCM = 0x0001; # PCM (8/16/32 bit)
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	47 FORMAT_IEEE_FLOAT = 0x0003; # IEEE float (32/64 bit)
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	48 BYTEORDER = "ieee-le";
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	49
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	50 if (nargin < 1 \|\| nargin > 2)
6046 34f96dd5441b [project @ 2006-10-10 16:10:25 by jwe] jwe parents: 5937 diff changeset	51 print_usage ();
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	52 endif
80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	53
80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	54 if (! ischar (filename))
10738 a4b8364e04c7 wavread.m: Correctly handle non-word aligned chunks (bug #30309). Rik <octave@nomad.inbox5.com> parents: 10635 diff changeset	55 error ("wavread: FILENAME must be a character string");
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	56 endif
80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	57
13003 74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	58 fid = -1;
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	59
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	60 unwind_protect
8796 9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	61
13003 74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	62 [fid, msg] = fopen (filename, "rb");
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	63
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	64 if (fid < 0)
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	65 error ("wavread: %s", msg);
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	66 endif
8796 9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	67
13003 74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	68 ## Get file size.
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	69 fseek (fid, 0, "eof");
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	70 file_size = ftell (fid);
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	71 fseek (fid, 0, "bof");
8796 9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	72
13003 74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	73 ## Find RIFF chunk.
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	74 riff_size = find_chunk (fid, "RIFF", file_size);
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	75 riff_pos = ftell (fid);
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	76 if (riff_size == -1)
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	77 error ("wavread: file contains no RIFF chunk");
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	78 endif
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	79
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	80 riff_type = char (fread (fid, 4))';
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	81 if (! strcmp (riff_type, "WAVE"))
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	82 error ("wavread: file contains no WAVE signature");
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	83 endif
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	84 riff_pos = riff_pos + 4;
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	85 riff_size = riff_size - 4;
8796 9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	86
13003 74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	87 ## Find format chunk inside the RIFF chunk.
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	88 fseek (fid, riff_pos, "bof");
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	89 fmt_size = find_chunk (fid, "fmt ", riff_size);
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	90 fmt_pos = ftell(fid);
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	91 if (fmt_size == -1)
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	92 error ("wavread: file contains no format chunk");
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	93 endif
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	94
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	95 ## Find data chunk inside the RIFF chunk.
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	96 ## We don't assume that it comes after the format chunk.
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	97 fseek (fid, riff_pos, "bof");
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	98 data_size = find_chunk (fid, "data", riff_size);
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	99 data_pos = ftell (fid);
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	100 if (data_size == -1)
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	101 error ("wavread: file contains no data chunk");
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	102 endif
11587 c792872f8942 all script files: untabify and strip trailing whitespace John W. Eaton <jwe@octave.org> parents: 11563 diff changeset	103
13003 74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	104 ### Read format chunk.
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	105 fseek (fid, fmt_pos, "bof");
11587 c792872f8942 all script files: untabify and strip trailing whitespace John W. Eaton <jwe@octave.org> parents: 11563 diff changeset	106
13003 74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	107 ## Sample format code.
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	108 format_tag = fread (fid, 1, "uint16", 0, BYTEORDER);
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	109 if (format_tag != FORMAT_PCM && format_tag != FORMAT_IEEE_FLOAT)
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	110 error ("wavread: sample format %#x is not supported", format_tag);
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	111 endif
11587 c792872f8942 all script files: untabify and strip trailing whitespace John W. Eaton <jwe@octave.org> parents: 11563 diff changeset	112
13003 74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	113 ## Number of interleaved channels.
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	114 channels = fread (fid, 1, "uint16", 0, BYTEORDER);
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	115
13003 74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	116 ## Sample rate.
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	117 samples_per_sec = fread (fid, 1, "uint32", 0, BYTEORDER);
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	118
13003 74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	119 ## Bits per sample.
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	120 fseek (fid, 6, "cof");
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	121 bits_per_sample = fread (fid, 1, "uint16", 0, BYTEORDER);
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	122
13003 74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	123 ### Read data chunk.
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	124 fseek (fid, data_pos, "bof");
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	125
13003 74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	126 ## Determine sample data type.
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	127 if (format_tag == FORMAT_PCM)
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	128 switch (bits_per_sample)
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	129 case 8
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	130 format = "uint8";
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	131 case 16
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	132 format = "int16";
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	133 case 24
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	134 format = "uint8";
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	135 case 32
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	136 format = "int32";
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	137 otherwise
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	138 error ("wavread: %d bits sample resolution is not supported with PCM",
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	139 bits_per_sample);
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	140 endswitch
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	141 else
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	142 switch (bits_per_sample)
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	143 case 32
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	144 format = "float32";
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	145 case 64
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	146 format = "float64";
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	147 otherwise
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	148 error ("wavread: %d bits sample resolution is not supported with IEEE float",
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	149 bits_per_sample);
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	150 endswitch
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	151 endif
11587 c792872f8942 all script files: untabify and strip trailing whitespace John W. Eaton <jwe@octave.org> parents: 11563 diff changeset	152
13003 74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	153 ## Parse arguments.
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	154 if (nargin == 1)
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	155 length = idivide (8 * data_size, bits_per_sample);
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	156 else
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	157 nparams = numel (param);
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	158 if (nparams == 1)
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	159 ## Number of samples is given.
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	160 length = param * channels;
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	161 elseif (nparams == 2)
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	162 ## Sample range is given.
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	163 if (fseek (fid, (param(1)-1) * channels * (bits_per_sample/8), "cof") < 0)
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	164 warning ("wavread: seeking failed");
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	165 endif
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	166 length = (param(2)-param(1)+1) * channels;
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	167 elseif (nparams == 4 && char (param) == "size")
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	168 ## Size of the file is requested.
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	169 tmp = idivide (8 * data_size, channels * bits_per_sample);
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	170 y = [tmp, channels];
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	171 return;
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	172 else
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	173 error ("wavread: invalid PARAM argument");
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	174 endif
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	175 endif
11587 c792872f8942 all script files: untabify and strip trailing whitespace John W. Eaton <jwe@octave.org> parents: 11563 diff changeset	176
13003 74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	177 ## Read samples and close file.
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	178 if (bits_per_sample == 24)
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	179 length *= 3;
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	180 endif
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	181
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	182 [yi, n] = fread (fid, length, format, 0, BYTEORDER);
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	183
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	184 unwind_protect_cleanup
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	185
74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	186 if (fid >= 0)
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	187 fclose (fid);
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	188 endif
5839 0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	189
13003 74c5fa5cef47 use unwind_protect to ensure wavread closes file. John W. Eaton <jwe@octave.org> parents: 13001 diff changeset	190 end_unwind_protect
5839 0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	191
8506 bc982528de11 comment style fixes John W. Eaton <jwe@octave.org> parents: 7461 diff changeset	192 ## Check data.
6304 78a3254741b2 [project @ 2007-02-14 21:22:59 by jwe] jwe parents: 6046 diff changeset	193 if (mod (numel (yi), channels) != 0)
78a3254741b2 [project @ 2007-02-14 21:22:59 by jwe] jwe parents: 6046 diff changeset	194 error ("wavread: data in %s doesn't match the number of channels",
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9598 diff changeset	195 filename);
6304 78a3254741b2 [project @ 2007-02-14 21:22:59 by jwe] jwe parents: 6046 diff changeset	196 endif
78a3254741b2 [project @ 2007-02-14 21:22:59 by jwe] jwe parents: 6046 diff changeset	197
5839 0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	198 if (bits_per_sample == 24)
0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	199 yi = reshape (yi, 3, rows(yi)/3)';
0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	200 yi(yi(:,3) >= 128, 3) -= 256;
0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	201 yi = yi * [1; 256; 65536];
7151 aeeb646f6538 [project @ 2007-11-09 19:34:17 by jwe] jwe parents: 7117 diff changeset	202 endif
aeeb646f6538 [project @ 2007-11-09 19:34:17 by jwe] jwe parents: 7117 diff changeset	203
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	204 if (format_tag == FORMAT_PCM)
8506 bc982528de11 comment style fixes John W. Eaton <jwe@octave.org> parents: 7461 diff changeset	205 ## Normalize samples.
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	206 switch (bits_per_sample)
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	207 case 8
12701 de3e90a420e3 Overhaul wavwrite, wavread and fix normalization problem (Bug #33420). Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	208 yi = (yi - 128)/128;
5572 c45cf76df06f [project @ 2005-12-12 20:54:06 by jwe] jwe parents: 5567 diff changeset	209 case 16
12701 de3e90a420e3 Overhaul wavwrite, wavread and fix normalization problem (Bug #33420). Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	210 yi /= 32768;
5839 0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	211 case 24
12701 de3e90a420e3 Overhaul wavwrite, wavread and fix normalization problem (Bug #33420). Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	212 yi /= 8388608;
5572 c45cf76df06f [project @ 2005-12-12 20:54:06 by jwe] jwe parents: 5567 diff changeset	213 case 32
12701 de3e90a420e3 Overhaul wavwrite, wavread and fix normalization problem (Bug #33420). Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	214 yi /= 2147483648;
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	215 endswitch
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	216 endif
11587 c792872f8942 all script files: untabify and strip trailing whitespace John W. Eaton <jwe@octave.org> parents: 11563 diff changeset	217
8506 bc982528de11 comment style fixes John W. Eaton <jwe@octave.org> parents: 7461 diff changeset	218 ## Deinterleave.
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	219 nr = numel (yi) / channels;
80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	220 y = reshape (yi, channels, nr)';
11587 c792872f8942 all script files: untabify and strip trailing whitespace John W. Eaton <jwe@octave.org> parents: 11563 diff changeset	221
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	222 endfunction
8796 9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	223
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	224 ## Given a chunk_id, scan through chunks from the current file position
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	225 ## though at most size bytes. Return the size of the found chunk, with
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	226 ## file position pointing to the start of the chunk data. Return -1 for
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	227 ## size if chunk is not found.
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	228
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	229 function chunk_size = find_chunk (fid, chunk_id, size)
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	230 id = "";
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	231 offset = 8;
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	232 chunk_size = 0;
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	233
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	234 while (! strcmp (id, chunk_id) && (offset < size))
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	235 fseek (fid, chunk_size, "cof");
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	236 id = char (fread (fid, 4))';
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	237 chunk_size = fread (fid, 1, "uint32", 0, "ieee-le");
10738 a4b8364e04c7 wavread.m: Correctly handle non-word aligned chunks (bug #30309). Rik <octave@nomad.inbox5.com> parents: 10635 diff changeset	238 ## Chunk sizes must be word-aligned (2 byte)
11587 c792872f8942 all script files: untabify and strip trailing whitespace John W. Eaton <jwe@octave.org> parents: 11563 diff changeset	239 chunk_size += rem (chunk_size, 2);
8796 9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	240 offset = offset + 8 + chunk_size;
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	241 endwhile
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	242 if (! strcmp (id, chunk_id))
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	243 chunk_size = -1;
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	244 endif
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	245 endfunction
12701 de3e90a420e3 Overhaul wavwrite, wavread and fix normalization problem (Bug #33420). Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	246
12799 e0b174e99449 codesprint: fix comment John W. Eaton <jwe@octave.org> parents: 12701 diff changeset	247 ## Mark file as being tested. Tests for wavread/wavwrite pair are in
12819 66af8b914607 codesprint: Add comment for dummy test in surface.m, Fix typo in wavread.m Kai Habel <kai.habel@gmx.de> parents: 12799 diff changeset	248 ## wavwrite.m
12799 e0b174e99449 codesprint: fix comment John W. Eaton <jwe@octave.org> parents: 12701 diff changeset	249 %!assert(1)

Mercurial > hg > octave-nkf

annotate scripts/audio/wavread.m @ 14271:e2a14d1b4eaa