Emgu 3.0+中的图像DFT

时间:2016-08-29 17:46:23

标签: vb.net visual-studio opencv emgucv dft

如何在Visual Basic中使用Emgu 3.0+计算图像的DFT?

我有一些图像我用CvInvoke.Imread()作为Mat加载,我可以显示和更改,我想计算该图像的前向DFT,然后显示幅度。我怎样才能做到这一点?我在将有效的IInputArrays输入CvInvoke.Dft()时遇到了麻烦。我的Emgu版本目前是3.1.0。

这是我目前正在使用的代码:

    Dim imgOrig As Mat
    imgOrig = CvInvoke.Imread(ofd.FileName, LoadImageType.AnyColor)

    EmguImageBox.Image = imgOrig

    Dim imgGrayscale As New Mat()
    CvInvoke.CvtColor(imgOrig, imgGrayscale, ColorConversion.Bgr2Gray)

    EmguImageBox.Image = imgGrayscale

    Dim imgDFT As Mat
    imgDFT = New Mat(imgOrig.Rows, imgOrig.Cols, DepthType.Cv32F, 2)

    CvInvoke.Dft(imgGrayscale, imgDFT, DxtType.Forward, 0)

    EmguImageBox.Image = imgDFT

1 个答案:

答案 0 :(得分:0)

也许您已经解决了这个问题,但是这里的一些信息可能会使其他遇到相同问题的人走上正确的道路。我自己曾尝试在Emgu 3.1.0版中执行相同的操作,但仅发现使用旧版Emgu的示例。终于有了一些工作,也许您可​​以在VB中实现相同的功能。

 #region Perform FFT
                //Get optimal DFT size
                int m = CvInvoke.GetOptimalDFTSize(originalImg.Height);
                int n = CvInvoke.GetOptimalDFTSize(originalImg.Width);

                //Resize image
                Mat resizedImg = new Mat();
                CvInvoke.CopyMakeBorder(originalImg, resizedImg, 0, m - originalImg.Height, 0, n - originalImg.Width, Emgu.CV.CvEnum.BorderType.Default);

                //Create imaginary channel for original image    
                Mat imChannel = new Mat(resizedImg.Size, Emgu.CV.CvEnum.DepthType.Cv32F, 1);
                imChannel.SetTo(new MCvScalar(0));
                //Add both real and imaginary channel to Mat vector
                Emgu.CV.Util.VectorOfMat vectorMat = new Emgu.CV.Util.VectorOfMat(new Mat[2] { resizedImg, imChannel });
                //Create 2 channel Mat from the vector
                Mat FFTimg = new Mat(resizedImg.Size, Emgu.CV.CvEnum.DepthType.Cv32F, 2);
                CvInvoke.Merge(vectorMat, FFTimg);

                //Perform DFT
                CvInvoke.Dft(FFTimg, FFTimg, Emgu.CV.CvEnum.DxtType.Forward, 0);

                //Split channels to separate Mat objects
                CvInvoke.Split(FFTimg, vectorMat);

                //Calculate magnitudes
                Mat fftMag = new Mat();
                CvInvoke.Pow(vectorMat[0], 2, vectorMat[0]);
                CvInvoke.Pow(vectorMat[1], 2, vectorMat[1]);
                CvInvoke.Add(vectorMat[0], vectorMat[1], fftMag);
                CvInvoke.Sqrt(fftMag, fftMag);

                //Swith quadrants 
                Mat q1 = new Mat(fftMag, new System.Drawing.Rectangle(0, (int)(fftMag.Height - fftMag.Height / 2), (int)(fftMag.Width / 4), fftMag.Height - (int)(fftMag.Height - fftMag.Height / 2)));
                Mat q2 = new Mat(fftMag, new System.Drawing.Rectangle(0, 0, (int)(fftMag.Width / 4), (int)(fftMag.Height / 2)));

                Mat tmp = new Mat();
                q1.CopyTo(tmp);
                q2.CopyTo(q1);
                tmp.CopyTo(q2);
                #endregion