c++实现 griddata函数

你可以使用以下代码实现 C++ 中的 griddata 函数：

#include <vector>
#include <cmath>

// 定义数据点结构体
struct Point {
    double x;
    double y;
    double z;
};

// 计算两个点之间的距离
double distance(const Point& p1, const Point& p2) {
    double dx = p2.x - p1.x;
    double dy = p2.y - p1.y;
    return std::sqrt(dx * dx + dy * dy);
}

// 使用 griddata 插值算法计算给定点的值
double griddata(const std::vector<Point>& points, const Point& queryPoint) {
    double sumZ = 0.0;
    double sumWeights = 0.0;

    // 遍历所有数据点
    for (const Point& point : points) {
        // 计算当前数据点与查询点的距离
        double d = distance(point, queryPoint);

        // 如果距离为0，则直接返回该数据点的值
        if (d == 0.0) {
            return point.z;
        }

        // 计算权重（距离的倒数）
        double weight = 1.0 / d;

        // 更新加权和
        sumZ += weight * point.z;
        sumWeights += weight;
    }

    // 计算插值结果
    return sumZ / sumWeights;
}

int main() {
    // 创建数据点集合
    std::vector<Point> points = {
        {1.0, 1.0, 10.0},
        {2.0, 2.0, 20.0},
        {3.0, 3.0, 30.0}
    };

    // 创建查询点
    Point queryPoint = {2.5, 2.5, 0.0};

    // 调用 griddata 函数计算查询点的值
    double interpolatedValue = griddata(points, queryPoint);

    // 输出插值结果
    std::cout << "Interpolated value: " << interpolatedValue << std::endl;

    return 0;
}

这段代码实现了一个简单的 griddata 函数，它使用了距离加权插值算法来计算给定点的值。你可以根据自己的需求修改数据点集合和查询点，然后调用 griddata 函数来获得插值结果。请注意，这只是一个简单的示例，实际应用中可能需要考虑更复杂的情况和算法。