需求如下圖，需要將綠色點(diǎn)圍成的區(qū)域（記做inside）到紅色點(diǎn)圍成的區(qū)域（記做outside）拉伸到藍(lán)色點(diǎn)（圖片中有兩圈藍(lán)點(diǎn)，本文以內(nèi)圈藍(lán)點(diǎn)為準(zhǔn)）圍成的區(qū)域（記做affine），并且綠色區(qū)域內(nèi)部的圖像保持不變，僅拉伸邊緣區(qū)域，原圖中紅色區(qū)域以外的圖像信息忽略。

原理：將紅點(diǎn)與綠點(diǎn)間的區(qū)域切分為多個(gè)三角形，記住srcTriangles， 將藍(lán)點(diǎn)與綠點(diǎn)間的區(qū)域也切分為多個(gè)三角形，記住dstTriangles。使用仿射變換，將srcTriangles變換為dstTriangles即可實(shí)現(xiàn)邊緣拉伸的效果。

注意三角形切分時(shí)，srcTriangles與dstTriangles一一對(duì)應(yīng)

涉及Opencv API：getAffineTransform（仿射變換）

仿射變換的用法：（三點(diǎn)確定仿射矩陣）

		src_tri[0] = Point2f(srcTriangles[i].p1.x, srcTriangles[i].p1.y);
		src_tri[1] = Point2f(srcTriangles[i].p2.x, srcTriangles[i].p2.y);
		src_tri[2] = Point2f(srcTriangles[i].p3.x, srcTriangles[i].p3.y);
		dst_tri[0] = Point2f(dstTriangles[i].p1.x, dstTriangles[i].p1.y);
		dst_tri[1] = Point2f(dstTriangles[i].p2.x, dstTriangles[i].p2.y);
		dst_tri[2] = Point2f(dstTriangles[i].p3.x, dstTriangles[i].p3.y);
		Mat warp_mat_org_inv = cv::getAffineTransform(dst_tri, src_tri);

實(shí)現(xiàn)需求的核心代碼：

Mat deform(std::vector<Triangle> srcTriangles, std::vector<Triangle> dstTriangles, Mat srcImage) {
	int nIsoSize = 512;
	Mat result = Mat::zeros(nIsoSize, nIsoSize, CV_8UC3);
	for (int i = 0; i < srcTriangles.size(); i++)
	{
		cv::Point2f src_tri[3];
		cv::Point2f dst_tri[3];
		src_tri[0] = Point2f(srcTriangles[i].p1.x, srcTriangles[i].p1.y);
		src_tri[1] = Point2f(srcTriangles[i].p2.x, srcTriangles[i].p2.y);
		src_tri[2] = Point2f(srcTriangles[i].p3.x, srcTriangles[i].p3.y);
		dst_tri[0] = Point2f(dstTriangles[i].p1.x, dstTriangles[i].p1.y);
		dst_tri[1] = Point2f(dstTriangles[i].p2.x, dstTriangles[i].p2.y);
		dst_tri[2] = Point2f(dstTriangles[i].p3.x, dstTriangles[i].p3.y);
		Mat warp_mat_org_inv = cv::getAffineTransform(dst_tri, src_tri);
		warp_mat_org_inv.convertTo(warp_mat_org_inv, CV_32FC1);
		for (int x = min(dst_tri[0].x, min(dst_tri[1].x, dst_tri[2].x)); x < max(dst_tri[0].x, max(dst_tri[1].x, dst_tri[2].x)); ++x) {
			for (int y = min(dst_tri[0].y, min(dst_tri[1].y, dst_tri[2].y)); y < max(dst_tri[0].y, max(dst_tri[1].y, dst_tri[2].y)); ++y) {
				if (is_point_in_triangle(cv::Point2f(x, y), dst_tri[0], dst_tri[1], dst_tri[2])) {
					// calculate corresponding position of dst_coord pixel center in image (src)
					const Mat homogenous_dst_coord(Vec3f(x, y, 1.0f));
					const Vec2f src_texel = Mat(warp_mat_org_inv * homogenous_dst_coord);
					if ((round(src_texel[1]) < srcImage.rows) && (round(src_texel[0]) < srcImage.cols) && round(src_texel[0]) > 0 && round(src_texel[1]) > 0)
					{
						int a1 = round(src_texel[1]);
						int a2 = round(src_texel[0]);
						result.at<cv::Vec3b>(y, x)[0] = srcImage.at<cv::Vec3b>(round(src_texel[1]), round(src_texel[0]))[0];
						result.at<cv::Vec3b>(y, x)[1] = srcImage.at<cv::Vec3b>(round(src_texel[1]), round(src_texel[0]))[1];
						result.at<cv::Vec3b>(y, x)[2] = srcImage.at<cv::Vec3b>(round(src_texel[1]), round(src_texel[0]))[2];
					}
				}
			}
		}
	}
	return result;
}