Beispiel #1
0
        public Form2()
        {
            InitializeComponent();

            FrequencyThroughTime ftt = new FrequencyThroughTime()
            {
                SamplesPerSecond = 44010
            };


            Bitmap bmpA = new Bitmap(Image.FromFile("d:\\demo1_filtereda.png", true));
            Bitmap bmpB = new Bitmap(Image.FromFile("d:\\demo1_filteredb.png", true));

            for (int y = 0; y < bmpA.Height; y++)
            {
                FrequencySample fs = new FrequencySample()
                {
                    windowLength = (uint)bmpA.Width, Samples = new List <Complex>()
                };

                for (int x = 0; x < bmpA.Width; x++)
                {
                    Color c  = bmpA.GetPixel(x, y);
                    Color c2 = bmpB.GetPixel(x, y);

                    fs.Samples.Add(new Complex(((c.B) << 8) + c2.B - short.MaxValue, ((c.G) << 8) + c2.G - short.MaxValue));
                }

                ftt.Samples.Add(fs);
            }

            WAV.Write("d:\\sample2.wav", ftt.Samples.SelectMany(x => SimpleDFT.SimpleDFT.IDFT(x.Samples)).Select(x => (short)x).ToList(), ftt.SamplesPerSecond);
        }
Beispiel #2
0
        private void doCalc()
        {
            FrequencyThroughTime ftt = new FrequencyThroughTime();
            SoundSample          ss  = new SoundSample()
            {
                Sound = sound, Start = 0, Length = (uint)sampleLength
            };

            int i = 0;

            while ((ss.Start <= sound.Samples.Count - ss.Length) && i < 5)
            {
                //dft
                FrequencySample s = FrequencySample.FromSoundSample(ss);

                ftt.Samples.Add(s);

                ss.Start += ss.Length;
                i++;
            }

            Func <FrequencyThroughTime, Bitmap[]> bitmaptize = (FrequencyThroughTime _ftt) =>
            {
                Bitmap _bmp  = new Bitmap(sampleLength / 2, ftt.Samples.Count, PixelFormat.Format32bppArgb);
                Bitmap _bmp2 = new Bitmap(sampleLength / 2, ftt.Samples.Count, PixelFormat.Format32bppArgb);

                for (int x = 0; x < ftt.Samples.Count; x++)
                {
                    FrequencySample s = _ftt.Samples[x];

                    for (int y = 0; y < s.Samples.Count; y++)
                    {
                        int im = (int)(ushort.MaxValue / 2 + (int)(s.Samples[y].Imaginary));
                        int re = (int)(ushort.MaxValue / 2 + (int)(s.Samples[y].Real));

                        _bmp.SetPixel(y, x, Color.FromArgb(
                                          0, (im & 0xFF00) >> 8,
                                          (re & 0xFF00) >> 8
                                          ));

                        _bmp2.SetPixel(y, x, Color.FromArgb(
                                           0, im & 0xFF,
                                           re & 0xFF
                                           ));
                    }
                }
                return(new Bitmap[] { _bmp, _bmp2 });
            };

            var reconstructed = ftt.Samples.SelectMany(x => SimpleDFT.SimpleDFT.IDFT(x.Samples)).Select(x => (short)x).ToList();

            /*
             * foreach (var r in sound.Samples.Skip(4200).Take(400))
             *  chart3.Series[0].Points.Add(r);
             *
             * foreach (var r in reconstructed.Skip(4200).Take(400))
             *  chart3.Series[1].Points.Add(r);
             * return;
             */

            var bmpB = bitmaptize(ftt);

            bmpB[0].Save("d:\\demo1a.bmp", ImageFormat.Bmp);
            bmpB[1].Save("d:\\demo1b.bmp", ImageFormat.Bmp);
            WAV.Write("d:\\sample1.wav", reconstructed, ftt.SamplesPerSecond);

            FrequencyThroughTime ftt2 = new FrequencyThroughTime()
            {
                SamplesPerSecond = 22100,
                Samples          = ftt.Samples.Select(fs => new FrequencySample()
                {
                    Samples = fs.Samples.Select(x => x.Magnitude > 20 ? x : Complex.Zero).ToList(), windowLength = fs.windowLength
                }).ToList()
            };

            var bmpA = bitmaptize(ftt2);

            bmpA[0].Save("d:\\demo1_filtereda.png", ImageFormat.Png);
            bmpA[1].Save("d:\\demo1_filteredb.png", ImageFormat.Png);

            WAV.Write("d:\\sample1_filteredx1.wav", ftt2.Samples.SelectMany(x => SimpleDFT.SimpleDFT.IDFT(x.Samples)).Select(x => (short)x).ToList(), ftt2.SamplesPerSecond);

            pictureBox1.SizeMode = PictureBoxSizeMode.Normal;
            pictureBox1.Image    = bmp;
        }