private void buttonHilbOutputStart_Click(object sender, RoutedEventArgs e)
{
var args = new HilbertWorkerArgs();
if (!Int32.TryParse(textBoxHilbFilterN.Text, out args.firLength) ||
(args.firLength & 1) == 0 || args.firLength <= 0)
{
MessageBox.Show("FIRフィルタ長は正の奇数の数値を半角数字で入力してください。処理中断。");
return;
}
args.inputPath = textBoxHilbInputPath.Text;
args.outputPath = textBoxHilbOutputPath.Text;
if (radioButtonHilbSint16.IsChecked == true)
{
args.outputBitsPerSample = 16;
args.valueRepresentationType = PcmData.ValueRepresentationType.SInt;
}
if (radioButtonHilbSint24.IsChecked == true)
{
args.outputBitsPerSample = 24;
args.valueRepresentationType = PcmData.ValueRepresentationType.SInt;
}
if (radioButtonHilbFloat32.IsChecked == true)
{
args.outputBitsPerSample = 32;
args.valueRepresentationType = PcmData.ValueRepresentationType.SFloat;
}
if (!Double.TryParse(textBoxHilbKaiserAlpha.Text, out args.kaiserAlpha) ||
args.kaiserAlpha <= 4.0 || 9.0 <= args.kaiserAlpha)
{
MessageBox.Show("カイザー窓α値は4.0<α<9.0の範囲の数値を半角数字で入力してください。処理中断。");
return;
}
if (radioButtonHilbBlackman.IsChecked == true)
{
args.windowFunc = WindowFuncType.Blackman;
}
else
{
args.windowFunc = WindowFuncType.Kaiser;
}
args.hilbertFilterType = WWHilbert.HilbertFilterType.HighPass;
if (radioButtonHilbertFilterBandlimited.IsChecked == true)
{
args.hilbertFilterType = WWHilbert.HilbertFilterType.Bandlimited;
}
AddHilbLog(string.Format("開始。{0} → {1}\r\n", args.inputPath, args.outputPath));
buttonHilbOutputStart.IsEnabled = false;
groupBoxHilbInput.IsEnabled = false;
groupBoxHilbSettings.IsEnabled = false;
groupBoxHilbOutput.IsEnabled = false;
m_HilbWorker.RunWorkerAsync(args);
}
private void buttonHilbOutputStart_Click(object sender, RoutedEventArgs e)
{
var args = new HilbertWorkerArgs();
if (!Int32.TryParse(textBoxHilbFilterN.Text, out args.firLength) ||
(args.firLength & 1) == 0 || args.firLength <= 0) {
MessageBox.Show("FIRフィルタ長は正の奇数の数値を半角数字で入力してください。処理中断。");
return;
}
args.inputPath = textBoxHilbInputPath.Text;
args.outputPath = textBoxHilbOutputPath.Text;
if (radioButtonHilbSint16.IsChecked == true) {
args.outputBitsPerSample = 16;
args.valueRepresentationType = PcmData.ValueRepresentationType.SInt;
}
if (radioButtonHilbSint24.IsChecked == true) {
args.outputBitsPerSample = 24;
args.valueRepresentationType = PcmData.ValueRepresentationType.SInt;
}
if (radioButtonHilbFloat32.IsChecked == true) {
args.outputBitsPerSample = 32;
args.valueRepresentationType = PcmData.ValueRepresentationType.SFloat;
}
if (!Double.TryParse(textBoxHilbKaiserAlpha.Text, out args.kaiserAlpha) ||
args.kaiserAlpha <= 4.0 || 9.0 <= args.kaiserAlpha) {
MessageBox.Show("カイザー窓α値は4.0<α<9.0の範囲の数値を半角数字で入力してください。処理中断。");
return;
}
if (radioButtonHilbBlackman.IsChecked == true) {
args.windowFunc = WindowFuncType.Blackman;
} else {
args.windowFunc = WindowFuncType.Kaiser;
}
args.hilbertFilterType = WWHilbert.HilbertFilterType.HighPass;
if (radioButtonHilbertFilterBandlimited.IsChecked == true) {
args.hilbertFilterType = WWHilbert.HilbertFilterType.Bandlimited;
}
AddHilbLog(string.Format("開始。{0} → {1}\r\n", args.inputPath, args.outputPath));
buttonHilbOutputStart.IsEnabled = false;
groupBoxHilbInput.IsEnabled = false;
groupBoxHilbSettings.IsEnabled = false;
groupBoxHilbOutput.IsEnabled = false;
m_HilbWorker.RunWorkerAsync(args);
}
private bool HilbertDo(FirWorkerArgs argsFir, PcmData pcmDataIn, out PcmData pcmDataOutput)
{
HilbertWorkerArgs args = argsFir as HilbertWorkerArgs;
var dft = new WWDirectComputeCS.WWDftCpu();
var hilb = WWHilbert.HilbertFirCoeff(args.hilbertFilterType, args.firLength);
System.Diagnostics.Debug.Assert(hilb.Length == args.firLength);
// 窓関数
double [] window;
if (args.windowFunc == WindowFuncType.Blackman)
{
WWWindowFunc.BlackmanWindow(args.firLength, out window);
}
else
{
WWWindowFunc.KaiserWindow(args.firLength, args.kaiserAlpha, out window);
}
// FIR coeffの個数は、window.Length個。
// ヒルベルト変換パラメータは未来から過去の方向に並んでいるので左右をひっくり返す。
double [] coeff = new double[args.firLength];
for (int i = 0; i < coeff.Length; ++i)
{
int pos = coeff.Length - i - 1;
coeff[i] = hilb[pos] * window[i];
}
/*
* for (int i=0; i < coeff.Length; ++i) {
* System.Console.WriteLine("coeff {0:D2} {1}", i, coeff[i]);
* }
* System.Console.WriteLine("");
*/
pcmDataOutput = new PcmData();
pcmDataOutput.CopyFrom(pcmDataIn);
for (int ch = 0; ch < pcmDataOutput.NumChannels; ++ch)
{
// 全てのチャンネルでループ。
var pcm1ch = new double[pcmDataOutput.NumFrames];
for (long i = 0; i < pcm1ch.Length; ++i)
{
pcm1ch[i] = pcmDataOutput.GetSampleValueInDouble(ch, i);
}
// 少しずつFIRする。
var fir = new WWFirCpu();
fir.Setup(coeff, pcm1ch);
const int FIR_SAMPLE = 65536;
for (int offs = 0; offs < pcm1ch.Length; offs += FIR_SAMPLE)
{
int nSample = FIR_SAMPLE;
if (pcm1ch.Length < offs + nSample)
{
nSample = pcm1ch.Length - offs;
}
var pcmFir = new double[nSample];
fir.Do(offs - window.Length / 2, nSample, pcmFir);
// 結果を出力に書き込む。
for (long i = 0; i < pcmFir.Length; ++i)
{
var re = pcmFir[i];
pcmDataOutput.SetSampleValueInDouble(ch, i + offs, re);
}
// 進捗Update。
int percentage = (int)(
(100L * ch / pcmDataOutput.NumChannels) +
(100L * (offs + 1) / pcm1ch.Length / pcmDataOutput.NumChannels));
m_HilbWorker.ReportProgress(percentage);
}
fir.Unsetup();
}
return(true);
}
