Method computeNewVector

src/ETF.cpp:102–122

* Paper's Eq(1) */

 * Paper's Eq(1)
 */
cv::Vec3f ETF::computeNewVector(const int x, const int y, const int kernel)
{
    const cv::Vec3f t_cur_x = flowField.at<cv::Vec3f>(y, x);
    cv::Vec3f t_new         = cv::Vec3f(0, 0, 0);

    for (int r = y - kernel; r <= y + kernel; ++r) {
        for (int c = x - kernel; c <= x + kernel; ++c) {
            if (r < 0 || r >= refinedETF.rows || c < 0 || c >= refinedETF.cols) continue;

            const cv::Vec3f t_cur_y = flowField.at<cv::Vec3f>(r, c);

            const float phi = computePhi(t_cur_x, t_cur_y);
            const float w_s = computeWs(cv::Point2f(x, y), cv::Point2f(c, r), kernel);
            const float w_m =
                computeWm(cv::norm(gradientMag.at<cv::Vec3f>(y, x)), cv::norm(gradientMag.at<float>(r, c)));
            const float w_d = computeWd(t_cur_x, t_cur_y);
            t_new += phi * t_cur_y * w_s * w_m * w_d;
        }
    }
    return cv::normalize(t_new);
}

/*
 * Paper's Eq(5)