Bug bei audioread

WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

22fuzzi

Forum-Anfänger


	Beiträge: 15

	Anmeldedatum: 17.05.16

	Wohnort: Graz

	Version: Octave ;)

Verfasst am: 26.04.2017, 21:13 Titel: Bug bei audioread

Liebes Forum!

Scheinbar habe ich ein Octave Problem entdeckt - welches mir sehr Kopfzerbrechen bereitet:
der Befehl "audioread" mit der DATATYPE 'native' macht Mist...

Code:

[sigINT,fs]=audioread('Motor_Hochlauf.wav','native');

load Data

sig_diff = sig_matlab-sig_octave; % bilde die Differenz der beiden Signale (sollte 0 sein)
sig_diff_fault = find(sig_diff~=0); % Wo ist sie ungleich Null?
figure, stem(sig_diff(sig_diff_fault))

sig_diff2 = sigINT-sig_octave; % bilde die Differenz der beiden Signale (sollte 0 sein)
sig_diff2_fault = find(sig_diff2~=0); % Wo ist sie ungleich Null?
figure, stem(sig_diff2(sig_diff2_fault))

save('-mat', 'Data.mat', 'sig_octave', 'sig_matlab')

Funktion ohne Link?

Bitte die Endung von Motor_Hochlauf in .wav ändern!!! (konnte ich nicht uploaden)

Würde mich über Feedback freuen - jenachdem welches Programm ihr verwendet!
Danke!
j

Motor_Hochlauf.wmv

Beschreibung:

Bitte die Endung in .wav ändern!!!

Download

Dateiname:

Motor_Hochlauf.wmv

Dateigröße:

1.16 MB

Heruntergeladen:

569 mal

Data.mat

Beschreibung:

Download

Dateiname:

Data.mat

Dateigröße:

2.32 MB

Heruntergeladen:

575 mal

AndyydnA

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 27.04.2017, 19:08 Titel:

Wo macht Octave hier Mist? Du schreibst ja nichts dazu aber man könnte meinen du meinst die +/-1 Differenz in sig_diff. Das könnte durch Rundung entstehen und ist, wenn du das SNR durch die Quantisierung betrachtest, vernachlässigbar.

22fuzzi

Themenstarter

Forum-Anfänger


	Beiträge: 15

	Anmeldedatum: 17.05.16

	Wohnort: Graz

	Version: Octave ;)

Verfasst am: 27.04.2017, 19:51 Titel:

Lieber Andy!

Ja genau, es geht um diese +/- 1!

Das Problem an der Sache ist:
im LSB von dem 16 bit Audio Signal ist eine wichtige Information codiert...
dieses Bit gehört eigentlich nicht zum Audio (wird bei der Aufnahme schon so gemacht...).
Kann ich dieses Bit nicht richtig auslesen - kann ich alle Messdaten schmeißen...

Ich habe 600 Mess .wav's und bei ca. 20 gibts da Probleme - und ich weis nicht warum...

Octave braucht da gar nichts runden - ich möchte einfach das LSB einlesen...
... bin etwas verzweifelt... sorry...

hab nen weiteren Thread aufgemacht:
http://www.gomatlab.de/alternative-zu-audioread-t43599.html

LG j

AndyydnA

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 28.04.2017, 15:25 Titel:

Liebster fuzzy,

wenn du in den Code schaust (UTSL) siehst du in audioread.cc Ab Zeile 123 (sic!)

Code:

Matrix audio (items, info.channels);
double *paudio = audio.fortran_vec ();
data += start * info.channels;
for (int i = 0; i < items; i++)
{
for (int channel = 0; channel < info.channels; channel++)
paudio[items*channel+i] = *data++;
}

octave_value ret_audio;
if ((nargin == 2 && args(1).is_string ()) || nargin == 3)
{
std::string type;
if (nargin == 3)
type = args(2).string_value ();
else
type = args(1).string_value ();

if (type == "native")
{
switch (info.format & SF_FORMAT_SUBMASK)
{
case SF_FORMAT_PCM_S8:
ret_audio = int8NDArray (audio * 127);
break;
case SF_FORMAT_PCM_U8:
ret_audio = uint8NDArray (audio * 127 + 127);
break;
case SF_FORMAT_PCM_16:
ret_audio = int16NDArray (audio * 32767);
break;
case SF_FORMAT_PCM_24:
ret_audio = int32NDArray (audio * 8388608);
break;
case SF_FORMAT_PCM_32:
ret_audio = int32NDArray (audio * 2147483648);
break;

Funktion ohne Link?

!

bei "native" wird also das wav als double gelesen und dann als uint16/whatever skaliert. Wie schon von mir gesagt, entstehen deine +/-1 genau durch diese Rundung.

Hättest du gleich zu Beginn gesagt, dass das LSB in der Rohdaten wichtig ist, hätten wir uns etwas Zeit und Nerven sparen können.

Eine mögliche Lösung ist offensichtlich: mach ein kleines.oct und verwende sf_read_int anstatt sf_read_float, kompilieren mit mkoctfile

Viel Erfolg, Gruß Andy

22fuzzi

Themenstarter

Forum-Anfänger


	Beiträge: 15

	Anmeldedatum: 17.05.16

	Wohnort: Graz

	Version: Octave ;)

Verfasst am: 02.05.2017, 15:17 Titel:

Lieber Andy!
Danke für die Aufklärung!
Ich versuche nun ein oct file zu erstellen - kome aber nicht weiter - bitte deshalb wieder um Hilfe!

Ich bekomme beim Kompilieren folgenden Fehler:

Code:

In file included from defun-dld.h:32:0,
from audioreadJF.cc:34:
audioreadJF.cc: In function 'octave_function* Gaudioread__JF(const octave::dynamic_library
&, bool)':
audioreadJF.cc:54:34: error: 'doc' was not declared in this scope
DEFUN_DLD (audioread__JF, args, ,doc: /* -*- texinfo -*-
^
C:\Octave\Octave-4.2.0\include\octave-4.2.0\octave/defun-int.h:123:56: note: in definition
of macro 'DEFINE_FUNX_INSTALLER_FUN'
= octave_dld_function::create (fname, shl, name, doc); \
^
defun-dld.h:47:3: note: in expansion of macro 'DEFINE_FUN_INSTALLER_FUN'
DEFINE_FUN_INSTALLER_FUN (name, doc) \
^
audioreadJF.cc:54:1: note: in expansion of macro 'DEFUN_DLD'
DEFUN_DLD (audioread__JF, args, ,doc: /* -*- texinfo -*-

Funktion ohne Link?

Habe eigentlich alles gelassen nur auf int geändert...
Leider bin ich kein c++ Auskenner...
Bitte um Hilfe / erklärung des Problems! Danke!
mein File:

Code:

#if defined (HAVE_CONFIG_H)
# include "config.h"
#endif

#include <string>
#include <map>

#include "oct-locbuf.h"
#include "unwind-prot.h"

#include "defun-dld.h"
#include "error.h"
#include "errwarn.h"
#include "ovl.h"
#include "ov.h"
#include "ov-struct.h"
#include "pager.h"

#if defined (HAVE_SNDFILE)
# include <sndfile.h>
#endif

#if defined (HAVE_SNDFILE)
static void
safe_close (SNDFILE *file)
{
sf_close (file);
}
#endif

DEFUN_DLD (audioread__JF, args, ,
doc: /* -*- texinfo -*-
@deftypefn {} {[@var{y}, @var{fs}] =} audioread__JF (@var{filename})
@deftypefnx {} {[@var{y}, @var{fs}] =} audioread__JF (@var{filename}, @var{samples})

@deftypefnx {} {[@var{y}, @var{fs}] =} audioread__JF (@var{filename}, @var{datatype})
@deftypefnx {} {[@var{y}, @var{fs}] =} audioread__JF (@var{filename}, @var{samples}, @var{datatype})
Read the audio file @var{filename} and return the audio data @var{y} and
sampling rate @var{fs}.

The audio data is stored as matrix with rows corresponding to audio frames
and columns corresponding to channels.

The optional two-element vector argument @var{samples} specifies starting
and ending frames.

The optional argument @var{datatype} specifies the datatype to return.
If it is @qcode{"native"}, then the type of data depends on how the data
is stored in the audio file.
@end deftypefn */)
{
#if defined (HAVE_SNDFILE)

int nargin = args.length ();

if (nargin < 1 || nargin > 3)
print_usage ();

std::string filename = args(0).xstring_value ("audioread__JF: FILENAME must be a string");

SF_INFO info;
info.format = 0;
SNDFILE *file = sf_open (filename.c_str (), SFM_READ, &info);

if (! file)
error ("audioread__JF: failed to open input file %s", filename.c_str ());

octave::unwind_protect frame;

frame.add_fcn (safe_close, file);

OCTAVE_LOCAL_BUFFER (int, data, info.frames * info.channels);

sf_read_int (file, data, info.frames * info.channels);

sf_count_t start = 0;
sf_count_t end = info.frames;

if ((nargin == 2 && ! args(1).is_string ()) || nargin == 3)
{
RowVector range = args(1).row_vector_value ();

if (range.numel () != 2)
error ("audioread__JF: invalid specification for range of frames");

double dstart = octave::math::isinf (range(0)) ? info.frames : range(0);
double dend = octave::math::isinf (range(1)) ? info.frames : range(1);

if (dstart < 1 || dstart > dend || dend > info.frames
|| octave::math::x_nint (dstart) != dstart || octave::math::x_nint (dend) != dend)
error ("audioread__JF: invalid specification for range of frames");

start = dstart - 1;
end = dend;
}

sf_count_t items = end - start;

Matrix audio (items, info.channels);

int *paudio = audio.fortran_vec ();

data += start * info.channels;

for (int i = 0; i < items; i++)
{
for (int channel = 0; channel < info.channels; channel++)
paudio[items*channel+i] = *data++;
}

octave_value ret_audio;

ret_audio = audio;

return ovl (ret_audio, info.samplerate);

#else

octave_unused_parameter (args);

err_disabled_feature ("audioread__JF",
"reading and writing sound files through libsndfile");

#endif
}

Funktion ohne Link?

AndyydnA

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 04.05.2017, 13:04 Titel:

So, ich habe mal etwas nach https://github.com/octave-de/macgyv.....master/audioread_int16.cc gelegt.

Kompilieren wie im Kommentar beschrieben. Ich habe soweit alles rausgeworfen, was sonst nicht gebraucht wird. Ggf. musst du bissle testen, dass ich die Kanäle nicht vertauscht oder sowas.

Btw, wenn du von einem GPL Programm kopierst, solltest du immer einen Verweis auf das original drin lassen und dein Programm auch unter die GPL stellen

Viel Erfolg

22fuzzi

Themenstarter

Forum-Anfänger


	Beiträge: 15

	Anmeldedatum: 17.05.16

	Wohnort: Graz

	Version: Octave ;)

Verfasst am: 29.05.2017, 14:01 Titel:

Lieber Andy!
Danke für deine Arbeit - leider funktioniert es nicht so wie gehofft!
Die eingelesenen Daten ergeben noch weniger Sinn...

Die beiden Kanäle (L,R) sind nicht in die Spalten 1 und 2 aufgeteilt, sondern in der ersten Spalte befindet sich zuerst der Wert von L und der zweite Wert ist der R Wert.

Da sind die beiden Kanäle ineinander verschachtelt...
Das Problem lässt sich zwar in Matlab recht einfach für eine Seite mit

Code:

signal=signal(1:2:end)

Funktion ohne Link?

beheben, im C-Code bring ichs leider nicht zusammen... bitte um Hilfe.

Danke im Voraus,
johannes

AndyydnA

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 30.05.2017, 10:07 Titel:

Hallo fuzzi,
da kam so lange keine Antwort, dass ich dachte das würde schon alles funktionieren.

Ich habe die Anordnung der Kanäle mit
https://github.com/octave-de/macgyv.....9757bf7ec1de5f6dd6bd2d67d
gefixed.

Ein "git pull && make" sollte ausreichen

HTH, Andy

AndyydnA

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 13.06.2017, 07:10 Titel:

Bitte, gern geschehen....


Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:	Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen

Impressum | Nutzungsbedingungen | Datenschutz | FAQ |

RSS

Hosted by:

Copyright © 2007 - 2024 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.