如何在opencvandroidjava中实现matlab的imgradient()函数

Question

我想在使用opencv.的android应用程序中使用matlab的imgradient()函数。我如何做到这一点？opencv的哪个函数与Matlab imgradient()函数等效。

我正在使用下面提到的函数，对吗？

public Mat imgradient(Mat grayScaleImage)
    {
        Mat grad_x=new Mat();
        Mat grad_y = new Mat();
        Mat abs_grad_x=new Mat();
        Mat abs_grad_y=new Mat();            
        Mat gradientImag = new Mat(grayScaleImage.rows(),grayScaleImage.cols(),CvType.CV_8UC1);
         Imgproc.Sobel(grayScaleImage, grad_x, CvType.CV_16S, 1, 0,3,1,0,Imgproc.BORDER_DEFAULT );
         Core.convertScaleAbs( grad_x, abs_grad_x );             
         Imgproc.Sobel( grayScaleImage, grad_y, CvType.CV_16S, 0, 1, 3, 1,0,Imgproc.BORDER_DEFAULT );
         Core.convertScaleAbs( grad_y, abs_grad_y );                
         double[] buff_grad = new double[1];
         for(int i = 0; i < abs_grad_y.cols(); i++)
            {
                for(int j =0 ; j<abs_grad_y.rows() ; j++)
                {
                    double[] buff_x = abs_grad_x.get(j, i);
                    double[] buff_y = abs_grad_y.get(j, i);
                    double x =  buff_x[0];
                    double y =  buff_y[0];
                    double ans=0;
                    try
                    {
                         ans = Math.sqrt(Math.pow(x,2)+Math.pow(y,2));
                    }catch(NullPointerException e)
                    {
                        ans = 0;
                    }
                    buff_grad[0] =  ans;                        
                    gradientImag.put(j, i, buff_grad);   
                }
            }           
        return gradientImag;
    }

Answer 1

您是否尝试过使用类似sobel或canny的运算符？

Answer 2

当matlab imgradient（）为坐标为x，y的每个像素返回梯度"幅度"（即sqrt（dx（x，y）²+dy（x，y）²）时，您可能需要执行以下操作：

// 1) Get the horizontal gradient
Mat kH = (cv::Mat_<double>(1,3) << -1,0,1); // differential kernel in x
Mat Dx;
filter2D(image, Dx, -1, kH, cv::Point(-1,-1), 0);
// 2) Get the vertical gradient
Mat kV = (cv::Mat_<double>(3,1) << -1,0,1); // differential kernel in y
Mat Dy;
filter2D(image, Dy, -1, kV, cv::Point(-1,-1), 0);
// 3) Get sqrt(dx²+dy²) in each point
for(int i=0; i<Dx.rows; i++)
    for(int j=0; j<Dx.cols; j++)
        Dmag.at<double>(i,j) = sqrt(pow(Dx.at<double>(i,j),2)+pow(Dy.at<double>(i,j),2));

它应该能让你得到你想要的。通过访问梯度数据，而不是对每个像素使用.at（i，j），可以获得更好的性能。

希望它能有所帮助！