如何计算IplImage的“Arctan”和“Pow”?

时间:2014-11-04 14:50:12

标签: opencv image-processing opencvsharp

我正在尝试使用OpenCvSharp计算Garyscale图像的渐变幅度和方向。问题是“Pow”功能似乎不适合IplImage。 我也想知道如何计算featureImage的tan-1(或arctan)。 谢谢

        using (IplImage cvImage = new IplImage("grayImage.png", LoadMode.AnyDepth | 
                LoadMode.GrayScale))
        using (IplImage dstXImage = new IplImage(cvImage.Size, cvImage.Depth, cvImage.NChannels))
        using (IplImage dstYImage = new IplImage(cvImage.Size, cvImage.Depth, cvImage.NChannels))
        {
            float[] data = {    0, -1, -1, 2    };
            CvMat kernel = new CvMat(2, 2, MatrixType.F32C1, data);

            Cv.Sobel(cvImage, dstXImage, 1, 0, ApertureSize.Size1);
            Cv.Sobel(cvImage, dstYImage, 0, 1, ApertureSize.Size1);

            Cv.Normalize(dstXImage, dstXImage, 1.0, 0, NormType.L1);
            Cv.Filter2D(cvImage, dstXImage, kernel, new CvPoint(0, 0));

            Cv.Normalize(dstYImage, dstYImage, 1.0, 0, NormType.L1);
            Cv.Filter2D(cvImage, dstYImage, kernel, new CvPoint(0, 0));

            // to calculate gradient magnitude, sqrt[(dy)power 2  +  (dx)power 2]
            dstXImage.Mul(dstXImage, dstXImage);
            dstYImage.Mul(dstYImage, dstYImage);
            IplImage dstXYImage = new IplImage(cvImage.Size, cvImage.Depth, cvImage.NChannels);
            dstXImage.Add(dstYImage, dstXYImage);
            dstXYImage.Pow(dstXYImage, 1/2); //this line not working,output image is black page


            // to calculate gradient orientation, arctan(dy/dx)
            IplImage thetaImage = new IplImage(cvImage.Size, cvImage.Depth, cvImage.NChannels);
            dstYImage.Div(dstXImage, thetaImage); //afterwards need help to calculate arctan 


            using (new CvWindow("SrcImage", cvImage))
            using (new CvWindow("DstXImage", dstXImage))
            using (new CvWindow("DstYImage", dstYImage))
            using (new CvWindow("DstXYImage", dstXYImage))
            using (new CvWindow("thetaImage", thetaImage))

            {
                Cv.WaitKey(0);
            }

1 个答案:

答案 0 :(得分:0)

您可以使用“cartToPolar”功能。 此函数计算2D矢量的大小和角度。

幅度(I)= sqrt(x(I)^ 2 + y(I)^ 2), angle(I)= atan2(y(I),x(I))[* 180 / pi]

例如:

IplImage dstXYImage;
IplImage thetaImage;
CartToPolar(dstXImage, dstYImage, dstXYImage, thetaImage, true);