Emgu 3.0+中图像的DFT

Question

如何在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

Answer 1

也许您已经解决了这个问题，但是这里有一些信息可能会使其他遇到相同问题的人走上正确的道路。我自己曾尝试在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