用opencv 实现halcon算子 integrate_funct_1d 功能

Halcon中的integrate_funct_1d算子用于计算一维函数在指定区间内的定积分值。在OpenCV中，可以使用cv::integral函数来实现类似的功能。以下是一个简单的示例代码，演示了如何使用cv::integral计算一维函数在指定区间内的定积分值。

#include <opencv2/opencv.hpp>
#include <iostream>

using namespace cv;
using namespace std;

double my_func(double x)
{
    return x * x;   // 定义一维函数 f(x) = x^2
}

double integral(double a, double b, int n)
{
    double h = (b - a) / n; // 计算步长
    double sum = 0.0;
    for(int i=0; i<n; i++)
    {
        double x1 = a + i * h;
        double x2 = a + (i+1) * h;
        // 计算矩形面积并累加
        sum += my_func((x1+x2)/2) * h;
    }
    return sum;
}

int main(int argc, char** argv)
{
    // 计算 x^2 在 [0, 1] 区间内的定积分值
    double a = 0.0, b = 1.0;
    int n = 1000;
    double result = integral(a, b, n);
    cout << "Integral value: " << result << endl;

    // 使用cv::integral计算定积分值
    Mat img(1, n, CV_64F);
    for(int i=0; i<n; i++)
    {
        double x = a + i * (b-a) / n;
        img.at<double>(0, i) = my_func(x);
    }
    Mat sum;
    integral(img, sum);
    double cv_result = sum.at<double>(0, n-1);
    cout << "CV Integral value: " << cv_result << endl;

    return 0;
}

在上述代码中，我们先定义了一个一维函数my_func(x)，然后使用自己编写的integral函数计算了该函数在[0, 1]区间内的定积分值。接着，我们使用cv::integral函数计算了相同区间内的定积分值，并将结果与自己编写的函数计算结果进行了比较。