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

annotate scripts/audio/wavread.m @ 12097:7d396910b6c5 release-3-2-x

wavread.m: rename data_size from obsolete ck_size

author	Christophe Tournery <christophe.tournery@illusonic.com>
date	Fri, 20 Nov 2009 10:44:02 +0100
parents	eb63fbe60fab
children	95c3e38098bf

rev	line source
8920 eb63fbe60fab update copyright notices John W. Eaton <jwe@octave.org> parents: 8796 diff changeset	1 ## Copyright (C) 2005, 2006, 2007, 2008, 2009 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 ##
6547 4fb053f24fd6 [project @ 2007-04-19 21:47:40 by jwe] jwe parents: 6304 diff changeset	25 ## @deftypefnx {Function File} {[@var{y}, @var{Fs}, @var{bits}] =} wavread (@var{filename})
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	26 ## Additionally return the sample rate (@var{fs}) in Hz and the number of bits
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	27 ## per sample (@var{bits}).
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 ##
6549 5a5a09d7deb8 [project @ 2007-04-20 06:55:29 by jwe] jwe parents: 6547 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))
80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	55 error ("wavwrite: expecting filename to be a character string");
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
8796 9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	58 # Open file for binary reading.
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	59 [fid, msg] = fopen (filename, "rb");
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	60 if (fid < 0)
5572 c45cf76df06f [project @ 2005-12-12 20:54:06 by jwe] jwe parents: 5567 diff changeset	61 error ("wavread: %s", msg);
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	62 endif
8796 9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	63
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	64 ## Get file size.
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	65 fseek (fid, 0, "eof");
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	66 file_size = ftell (fid);
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	67 fseek (fid, 0, "bof");
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	68
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	69 ## Find RIFF chunk.
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	70 riff_size = find_chunk (fid, "RIFF", file_size);
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	71 riff_pos = ftell (fid);
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	72 if (riff_size == -1)
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	73 fclose (fid);
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	74 error ("wavread: file contains no RIFF chunk");
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	75 endif
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	76
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	77 riff_type = char (fread (fid, 4))';
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	78 if(! strcmp (riff_type, "WAVE"))
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	79 fclose (fid);
8796 9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	80 error ("wavread: file contains no WAVE signature");
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	81 endif
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	82 riff_pos = riff_pos + 4;
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	83 riff_size = riff_size - 4;
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	84
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	85 ## Find format chunk inside the RIFF chunk.
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	86 fseek (fid, riff_pos, "bof");
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	87 fmt_size = find_chunk (fid, "fmt ", riff_size);
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	88 fmt_pos = ftell(fid);
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	89 if (fmt_size == -1)
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	90 fclose (fid);
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	91 error ("wavread: file contains no format chunk");
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	92 endif
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	93
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	94 ## Find data chunk inside the RIFF chunk.
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	95 ## We don't assume that it comes after the format chunk.
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	96 fseek (fid, riff_pos, "bof");
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	97 data_size = find_chunk (fid, "data", riff_size);
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	98 data_pos = ftell (fid);
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	99 if (data_size == -1)
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	100 fclose (fid);
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	101 error ("wavread: file contains no data chunk");
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	102 endif
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	103
8796 9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	104 ### Read format chunk.
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	105 fseek (fid, fmt_pos, "bof");
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	106
8506 bc982528de11 comment style fixes John W. Eaton <jwe@octave.org> parents: 7461 diff changeset	107 ## Sample format code.
6304 78a3254741b2 [project @ 2007-02-14 21:22:59 by jwe] jwe parents: 6046 diff changeset	108 format_tag = fread (fid, 1, "uint16", 0, BYTEORDER);
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	109 if (format_tag != FORMAT_PCM && format_tag != FORMAT_IEEE_FLOAT)
80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	110 fclose (fid);
80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	111 error ("wavread: sample format %#x is not supported", format_tag);
80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	112 endif
80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	113
8506 bc982528de11 comment style fixes John W. Eaton <jwe@octave.org> parents: 7461 diff changeset	114 ## Number of interleaved channels.
6304 78a3254741b2 [project @ 2007-02-14 21:22:59 by jwe] jwe parents: 6046 diff changeset	115 channels = fread (fid, 1, "uint16", 0, BYTEORDER);
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	116
8506 bc982528de11 comment style fixes John W. Eaton <jwe@octave.org> parents: 7461 diff changeset	117 ## Sample rate.
6304 78a3254741b2 [project @ 2007-02-14 21:22:59 by jwe] jwe parents: 6046 diff changeset	118 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	119
8506 bc982528de11 comment style fixes John W. Eaton <jwe@octave.org> parents: 7461 diff changeset	120 ## Bits per sample.
8796 9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	121 fseek (fid, 6, "cof");
6304 78a3254741b2 [project @ 2007-02-14 21:22:59 by jwe] jwe parents: 6046 diff changeset	122 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	123
8796 9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	124 ### Read data chunk.
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	125 fseek (fid, data_pos, "bof");
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	126
8506 bc982528de11 comment style fixes John W. Eaton <jwe@octave.org> parents: 7461 diff changeset	127 ## Determine sample data type.
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	128 if (format_tag == FORMAT_PCM)
5839 0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	129 switch (bits_per_sample)
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	130 case 8
5572 c45cf76df06f [project @ 2005-12-12 20:54:06 by jwe] jwe parents: 5567 diff changeset	131 format = "uint8";
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	132 case 16
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	133 format = "int16";
5839 0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	134 case 24
0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	135 format = "uint8";
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	136 case 32
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	137 format = "int32";
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	138 otherwise
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	139 fclose (fid);
6304 78a3254741b2 [project @ 2007-02-14 21:22:59 by jwe] jwe parents: 6046 diff changeset	140 error ("wavread: %d bits sample resolution is not supported with PCM",
78a3254741b2 [project @ 2007-02-14 21:22:59 by jwe] jwe parents: 6046 diff changeset	141 bits_per_sample);
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	142 endswitch
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	143 else
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	144 switch (bits_per_sample)
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	145 case 32
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	146 format = "float32";
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	147 case 64
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	148 format = "float64";
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	149 otherwise
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	150 fclose (fid);
6304 78a3254741b2 [project @ 2007-02-14 21:22:59 by jwe] jwe parents: 6046 diff changeset	151 error ("wavread: %d bits sample resolution is not supported with IEEE float",
78a3254741b2 [project @ 2007-02-14 21:22:59 by jwe] jwe parents: 6046 diff changeset	152 bits_per_sample);
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	153 endswitch
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	154 endif
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	155
8506 bc982528de11 comment style fixes John W. Eaton <jwe@octave.org> parents: 7461 diff changeset	156 ## Parse arguments.
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	157 if (nargin == 1)
8796 9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	158 length = 8 * data_size / bits_per_sample;
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	159 else
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	160 if (size (param, 2) == 1)
8506 bc982528de11 comment style fixes John W. Eaton <jwe@octave.org> parents: 7461 diff changeset	161 ## Number of samples is given.
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	162 length = param * channels;
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	163 elseif (size (param, 2) == 2)
8506 bc982528de11 comment style fixes John W. Eaton <jwe@octave.org> parents: 7461 diff changeset	164 ## Sample range is given.
8796 9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	165 if (fseek (fid, (param(1)-1) * channels * (bits_per_sample/8), "cof") < 0)
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	166 warning ("wavread: seeking failed");
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	167 endif
5839 0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	168 length = (param(2)-param(1)+1) * channels;
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	169 elseif (size (param, 2) == 4 && char (param) == "size")
8506 bc982528de11 comment style fixes John W. Eaton <jwe@octave.org> parents: 7461 diff changeset	170 ## Size of the file is requested.
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	171 fclose (fid);
12097 7d396910b6c5 wavread.m: rename data_size from obsolete ck_size Christophe Tournery <christophe.tournery@illusonic.com> parents: 8920 diff changeset	172 y = [data_size/channels/(bits_per_sample/8), channels];
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	173 return
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	174 else
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	175 fclose (fid);
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	176 error ("wavread: invalid argument 2");
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	177 endif
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	178 endif
5839 0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	179
8506 bc982528de11 comment style fixes John W. Eaton <jwe@octave.org> parents: 7461 diff changeset	180 ## Read samples and close file.
5839 0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	181 if (bits_per_sample == 24)
0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	182 length *= 3;
0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	183 endif
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	184 [yi, n] = fread (fid, length, format, 0, BYTEORDER);
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	185 fclose (fid);
5839 0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	186
8506 bc982528de11 comment style fixes John W. Eaton <jwe@octave.org> parents: 7461 diff changeset	187 ## Check data.
6304 78a3254741b2 [project @ 2007-02-14 21:22:59 by jwe] jwe parents: 6046 diff changeset	188 if (mod (numel (yi), channels) != 0)
78a3254741b2 [project @ 2007-02-14 21:22:59 by jwe] jwe parents: 6046 diff changeset	189 error ("wavread: data in %s doesn't match the number of channels",
78a3254741b2 [project @ 2007-02-14 21:22:59 by jwe] jwe parents: 6046 diff changeset	190 filename);
78a3254741b2 [project @ 2007-02-14 21:22:59 by jwe] jwe parents: 6046 diff changeset	191 endif
78a3254741b2 [project @ 2007-02-14 21:22:59 by jwe] jwe parents: 6046 diff changeset	192
5839 0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	193 if (bits_per_sample == 24)
0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	194 yi = reshape (yi, 3, rows(yi)/3)';
0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	195 yi(yi(:,3) >= 128, 3) -= 256;
0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	196 yi = yi * [1; 256; 65536];
7151 aeeb646f6538 [project @ 2007-11-09 19:34:17 by jwe] jwe parents: 7117 diff changeset	197 endif
aeeb646f6538 [project @ 2007-11-09 19:34:17 by jwe] jwe parents: 7117 diff changeset	198
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	199 if (format_tag == FORMAT_PCM)
8506 bc982528de11 comment style fixes John W. Eaton <jwe@octave.org> parents: 7461 diff changeset	200 ## Normalize samples.
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	201 switch (bits_per_sample)
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	202 case 8
7117 373c39904877 [project @ 2007-11-07 18:48:00 by jwe] jwe parents: 7017 diff changeset	203 yi = (yi - 128)/127;
5572 c45cf76df06f [project @ 2005-12-12 20:54:06 by jwe] jwe parents: 5567 diff changeset	204 case 16
7117 373c39904877 [project @ 2007-11-07 18:48:00 by jwe] jwe parents: 7017 diff changeset	205 yi /= 32767;
5839 0bc78d9ddcd4 [project @ 2006-06-01 21:21:11 by jwe] jwe parents: 5642 diff changeset	206 case 24
7117 373c39904877 [project @ 2007-11-07 18:48:00 by jwe] jwe parents: 7017 diff changeset	207 yi /= 8388607;
5572 c45cf76df06f [project @ 2005-12-12 20:54:06 by jwe] jwe parents: 5567 diff changeset	208 case 32
7117 373c39904877 [project @ 2007-11-07 18:48:00 by jwe] jwe parents: 7017 diff changeset	209 yi /= 2147483647;
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	210 endswitch
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	211 endif
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	212
8506 bc982528de11 comment style fixes John W. Eaton <jwe@octave.org> parents: 7461 diff changeset	213 ## Deinterleave.
5567 80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	214 nr = numel (yi) / channels;
80e629357483 [project @ 2005-12-07 06:31:28 by jwe] jwe parents: 5566 diff changeset	215 y = reshape (yi, channels, nr)';
5565 2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	216
2eeed655e801 [project @ 2005-12-06 20:38:20 by jwe] jwe parents: diff changeset	217 endfunction
8796 9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	218
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	219 ## 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	220 ## 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	221 ## 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	222 ## 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	223
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	224 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	225 id = "";
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	226 offset = 8;
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	227 chunk_size = 0;
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 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	230 fseek (fid, chunk_size, "cof");
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	231 id = char (fread (fid, 4))';
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	232 chunk_size = fread (fid, 1, "uint32", 0, "ieee-le");
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	233 offset = offset + 8 + chunk_size;
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	234 endwhile
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	235 if (! strcmp (id, chunk_id))
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	236 chunk_size = -1;
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	237 endif
9662dfb26652 wavread.m: improve search for data chunks Frederick Umminger <Frederick_Umminger@playstation.sony.com> parents: 8506 diff changeset	238 endfunction

Mercurial > octave-nkf

annotate scripts/audio/wavread.m @ 12097:7d396910b6c5 release-3-2-x