Обработка аудио

[OmegaBerkut]

Здравствуйте.
У меня есть два аудио-трека.
В первом - фоновая музыка.
Во втором - фоновая музыка с наложенным голосом.

Вопрос: имея только эти два трека, как можно выделить только голос ?

Попытался поиграться с профилями шумов в Adobe Audition: получилось так что очистка не полная (остаются некоторые звуки фоновой музыки), и в удаляемый шум захватывается значительная часть голоса, прям отчётливо слышно весь текст.

Что посоветуете ?

[WorldMaster]

А записано как? Там наверное одна запись чистая а вторая с шумами и голосом.
Может Фурье .. разложить оба сигнала и в кажду момент подавлять те частоты которые выделены в сигнале без голоса

[Desc]

OmegaBerkut, выкладывайте в тему оба файла, при-+-думаем как разделить.

[OmegaBerkut]

Сообщение от WorldMaster

А записано как? Там наверное одна запись чистая а вторая с шумами и голосом

Изначально это один трек, у которого в левом канале только фоновый звук, в правом канале - фоновый звук с наложенным голосом.
То есть - голос только в одном ухе.
Мне поставили задачу: зеркально продублировать правый и левый каналы, в записи голос получился в обоих каналах; с этой задачей я справился обычным наложением (микширование).

Потом мне стало интересно, можно ли вытянуть только голос.
Вытянуть голос из правого канала, отсеять фон; по этому вопросу сначала попробовал что знаю, потом пришёл сюда.

По сути, мне нужно не обычное наложение звука, а какое то хитрое.
Если говорить про картинки - в фотошопе это наложение с вычитанием, и тогда одинаковые пиксели становятся чёрными, а разница с целевой картинкой - инверсия.
Соответственно, обратной инверсией можно получить исходную разницу изображений (если это jpeg - будет минимальный остаточный шум, который нивелируется на этапе обратной инверсии).

Выкладываю ссылку на архив с тремя подписанными файлами. Сразу скажу - на содержимое внимание не обращайте, так как это тема для отдельного обсуждения (за что люди готовы платить).
http://pro-gramm.ru/_download/_Mixer.zip

P. S. Это вопрос исключительно личного интереса, за это никто не платит. У меня есть ещё пара треков, которые я хочу потом так же обработать

[Desc]

Файлы прослушал, теперь не совсем понятна постановка задачи. Одно дело, извлечь из гаммы необходимый спектр. А другое, из исходников создать себе геморрой и думать как от него избавиться.

Сообщение от OmegaBerkut

Мне поставили задачу: зеркально продублировать правый и левый каналы

Здесь вопросов нет, все прозрачно.

Сообщение от OmegaBerkut

Потом мне стало интересно, можно ли вытянуть только голос.

Т. е. из полученного "микса" получить исходный вариант, правильно понимаю постановку задания?

[OmegaBerkut]

Сообщение от Desc

Т. е. из полученного "микса" получить исходный вариант, правильно понимаю постановку задания?

Вот есть у меня фоновая музыка (правый канал), и фон+голос (левый канал). Так получилось ввиду того, что изначально смонтировано криво.
Я их разделил на две отдельные записи. У меня получилась первая запись - "шум", вторая запись - "шум"+голос. Задача: отсеять "шум" на второй записи.
Имея эти вводные, как можно из второй записи (изначально левого канала) вытянуть голос ?
Аналогию с изображениями написал выше.

[Pavia]

Музыка из треков разная.

[OmegaBerkut]

Сообщение от Pavia

Музыка из треков разная.

Может на слух так и кажется, правом канале голос может создавать видимость отличия музыки

[Pavia]

Для разных октав пришлось сделать разный порог чем выше частота тем больше разница. Зависимость взял квадратичную. Но нужно более крутую делать.

procedure TForm1.Button9Click(Sender: TObject);
const
  WindowSize=8*1024; //число отсчётов
var fd:Integer;
  SamplesPerSec:DWord;
  First, Second:TArrayArrayReal;
  First1, Second1, Thri:TArrayReal;
  a,b,a1,b1,c:TArrayReal;
  NilArray:TArrayReal;
  z1,z2:TarrayComplex;
  i,j:Integer;
  M1,M2:Real;
begin

ReadWave('E:\tmp\_Mixer\Микшер\0. Изначальная запись.wav', First, SamplesPerSec);
ReadWave('E:\tmp\_Mixer\Микшер\2. Левый канал.wav', Second, SamplesPerSec);
fd:=SamplesPerSec;
First1:=First[1];
Second1:=Second[0];
SetLength(Thri,Length(First1));


SetLength(NilArray,WindowSize);
SetLength(z1,WindowSize);
SetLength(z2,WindowSize);
SetLength(a,WindowSize);
SetLength(b,WindowSize);
SetLength(c,WindowSize);

SetLength(a1,WindowSize);
SetLength(b1,WindowSize);

 i:=0;
 while i<Length(First1)-WindowSize do
 begin
 Move(First1[i], a[0], (WindowSize)*SizeOf(a[0]));
 Move(Second1[i], b[0], (WindowSize)*SizeOf(b[0]));

 // Вычисляем спектр
 MulLanczosWindow(WindowSize, @a[0]);
 CopyInRe(z1,a);
 CopyInIm(z1,NilArray);
 FFT(z1,False);
 Amp(a1,Z1);
 Do_db(a1,-100);

 // Вычисляем спектр
 MulLanczosWindow(WindowSize, @b[0]);
 CopyInRe(z2,b);
 CopyInIm(z2,NilArray);
 FFT(z2,False);
 Amp(b1,Z2);
 Do_db(b1,-100);

 // Строем огибающие спектров
 for j:=0 to 256 do
  begin
  Blur(a1);
  Blur(b1);
  end;
 // Нормируем громкость
 M1:=Mean(a1);
 M2:=Mean(b1);
 add(b1,M1-M2);
 // Фильтруем близкие спектры
 for j:=0 to WindowSize-1 do
   if abs(a1[j]-b1[j])<5+sqr((8*j) div WindowSize) then
     begin
     Z1[j].Re:=Z1[j].Re*0.001; 
     Z1[j].Im:=Z1[j].Im*0.001; 
     end;

 // Выводим спектры
 if (i div 10 div WindowSize)<>(i+Round(0.2*WindowSize))div WindowSize div 10 then
   begin
   Amp(c,Z1);
   Do_db(c,-100);
   sub(c,20);
   AddArrayInChart(Chart1,0, a1, $80000000+1, fd);
   AddArrayInChart(Chart1,1, b1, $80000000+1, fd);
   AddArrayInChart(Chart1,2, c, $80000000+1, fd);
   Caption:=Format('%.2f', [i/Length(First1)]);
   repaint;        {}
   end;

 //Обратно в отсчеты
 FFT(Z1,true);
 CopyFromRe(a,z1);
 for j:=0 to Length(a)-1 do
   Thri[i+j]:=0.5*Thri[i+j]+0.5*a[j];

 i:=i+Round(0.2*WindowSize);  // Обрабатываем окнами с 80% перекрытием
 end;

for j:=0 to Length(Thri)-1 do
  begin
  First[0][j]:=Thri[j];
  First[1][j]:=Thri[j];
  end;
WriteWave('E:\tmp\_Mixer\Микшер\3. разница.wav', First, fd);

end;

[OmegaBerkut]

Сообщение от Pavia

Для разных октав пришлось сделать разный порог

Сильно =)
Спасибо конечно, гляну. Но вопрос: код скомпилирован ?
А то у меня Delphi не установлен, где-то на антресолях пылится =)