问用对角正弦波实现翘曲图像

EN

void deform()
{
    float alpha = 45 * CV_PI / 180.0; // wave direction
    float ox = cos(alpha);
    float oy = sin(alpha);
    
    cv::Mat src = cv::imread("F:/ImagesForTest/lena.jpg");

    for (int i = 0; i < src.rows; i+=8)
    { 
            cv::line(src, cv::Point(i, 0), cv::Point(i, src.rows),cv::Scalar(255,255,255));
    }

    for (int j = 0; j < src.cols; j += 8)
    { 
        cv::line(src, cv::Point(0,j), cv::Point(src.cols,j), cv::Scalar(255, 255, 255));
    }


    cv::Mat result = src.clone();
    for (int i = 0; i < src.rows; i++)
    { // to y
        for (int j = 0; j < src.cols; j++)
        { // to x
            float t =(i * oy)+ (j * ox); // wave parameter
            for (int ch = 0; ch < 3; ch++)
            { // each colour
                int offset_x =ox* (int)(25.0 * (sin(3.14 * t/ 150)));
                int offset_y =oy* (int)(25.0 * (sin(3.14 * t / 150)));
                if (i + offset_y < src.rows && j + offset_x < src.rows && i + offset_y >=0 && j + offset_x>=0)
                {
                    result.at<cv::Vec3b>(i, j)[ch] = src.at<cv::Vec3b>(i + offset_y, j + offset_x )[ch];
                }
                else
                    result.at<cv::Vec3b>(i, j)[ch] = 0.0;
            }
        }
    }

    cv:: imshow("result", result);
    cv::imwrite("result.jpg", result);
    cv::waitKey();
}