首页 最新 热门 推荐

  • 首页
  • 最新
  • 热门
  • 推荐

【图像处理】透视变换 Perspective Transformation

  • 23-09-22 20:02
  • 2771
  • 12099
blog.csdn.net

透视变换(Perspective Transformation)是将图片投影到一个新的视平面(Viewing Plane),也称作投影映射(Projective Mapping)。通用的变换公式为:


u,v是原始图片左边,对应得到变换后的图片坐标x,y,其中。
变换矩阵可以拆成4部分,表示线性变换,比如scaling,shearing和ratotion。用于平移,产生透视变换。所以可以理解成仿射等是透视变换的特殊形式。经过透视变换之后的图片通常不是平行四边形(除非映射视平面和原来平面平行的情况)。

重写之前的变换公式可以得到:


所以,已知变换对应的几个点就可以求取变换公式。反之,特定的变换公式也能新的变换后的图片。简单的看一个正方形到四边形的变换:
变换的4组对应点可以表示成:

根据变换公式得到:


定义几个辅助变量:


都为0时变换平面与原来是平行的,可以得到:


不为0时,得到:


求解出的变换矩阵就可以将一个正方形变换到四边形。反之,四边形变换到正方形也是一样的。于是,我们通过两次变换:四边形变换到正方形+正方形变换到四边形就可以将任意一个四边形变换到另一个四边形。


看一段代码:

  1. PerspectiveTransform::PerspectiveTransform(float inA11, float inA21,
  2. float inA31, float inA12,
  3. float inA22, float inA32,
  4. float inA13, float inA23,
  5. float inA33) :
  6. a11(inA11), a12(inA12), a13(inA13), a21(inA21), a22(inA22), a23(inA23),
  7. a31(inA31), a32(inA32), a33(inA33) {}
  8. PerspectiveTransform PerspectiveTransform::quadrilateralToQuadrilateral(float x0, float y0, float x1, float y1,
  9. float x2, float y2, float x3, float y3, float x0p, float y0p, float x1p, float y1p, float x2p, float y2p,
  10. float x3p, float y3p) {
  11. PerspectiveTransform qToS = PerspectiveTransform::quadrilateralToSquare(x0, y0, x1, y1, x2, y2, x3, y3);
  12. PerspectiveTransform sToQ =
  13. PerspectiveTransform::squareToQuadrilateral(x0p, y0p, x1p, y1p, x2p, y2p, x3p, y3p);
  14. return sToQ.times(qToS);
  15. }
  16. PerspectiveTransform PerspectiveTransform::squareToQuadrilateral(float x0, float y0, float x1, float y1, float x2,
  17. float y2, float x3, float y3) {
  18. float dx3 = x0 - x1 + x2 - x3;
  19. float dy3 = y0 - y1 + y2 - y3;
  20. if (dx3 == 0.0f && dy3 == 0.0f) {
  21. PerspectiveTransform result(PerspectiveTransform(x1 - x0, x2 - x1, x0, y1 - y0, y2 - y1, y0, 0.0f,
  22. 0.0f, 1.0f));
  23. return result;
  24. } else {
  25. float dx1 = x1 - x2;
  26. float dx2 = x3 - x2;
  27. float dy1 = y1 - y2;
  28. float dy2 = y3 - y2;
  29. float denominator = dx1 * dy2 - dx2 * dy1;
  30. float a13 = (dx3 * dy2 - dx2 * dy3) / denominator;
  31. float a23 = (dx1 * dy3 - dx3 * dy1) / denominator;
  32. PerspectiveTransform result(PerspectiveTransform(x1 - x0 + a13 * x1, x3 - x0 + a23 * x3, x0, y1 - y0
  33. + a13 * y1, y3 - y0 + a23 * y3, y0, a13, a23, 1.0f));
  34. return result;
  35. }
  36. }
  37. PerspectiveTransform PerspectiveTransform::quadrilateralToSquare(float x0, float y0, float x1, float y1, float x2,
  38. float y2, float x3, float y3) {
  39. // Here, the adjoint serves as the inverse:
  40. return squareToQuadrilateral(x0, y0, x1, y1, x2, y2, x3, y3).buildAdjoint();
  41. }
  42. PerspectiveTransform PerspectiveTransform::buildAdjoint() {
  43. // Adjoint is the transpose of the cofactor matrix:
  44. PerspectiveTransform result(PerspectiveTransform(a22 * a33 - a23 * a32, a23 * a31 - a21 * a33, a21 * a32
  45. - a22 * a31, a13 * a32 - a12 * a33, a11 * a33 - a13 * a31, a12 * a31 - a11 * a32, a12 * a23 - a13 * a22,
  46. a13 * a21 - a11 * a23, a11 * a22 - a12 * a21));
  47. return result;
  48. }
  49. PerspectiveTransform PerspectiveTransform::times(PerspectiveTransform other) {
  50. PerspectiveTransform result(PerspectiveTransform(a11 * other.a11 + a21 * other.a12 + a31 * other.a13,
  51. a11 * other.a21 + a21 * other.a22 + a31 * other.a23, a11 * other.a31 + a21 * other.a32 + a31
  52. * other.a33, a12 * other.a11 + a22 * other.a12 + a32 * other.a13, a12 * other.a21 + a22
  53. * other.a22 + a32 * other.a23, a12 * other.a31 + a22 * other.a32 + a32 * other.a33, a13
  54. * other.a11 + a23 * other.a12 + a33 * other.a13, a13 * other.a21 + a23 * other.a22 + a33
  55. * other.a23, a13 * other.a31 + a23 * other.a32 + a33 * other.a33));
  56. return result;
  57. }
  58. void PerspectiveTransform::transformPoints(vector<float> &points) {
  59. int max = points.size();
  60. for (int i = 0; i < max; i += 2) {
  61. float x = points[i];
  62. float y = points[i + 1];
  63. float denominator = a13 * x + a23 * y + a33;
  64. points[i] = (a11 * x + a21 * y + a31) / denominator;
  65. points[i + 1] = (a12 * x + a22 * y + a32) / denominator;
  66. }
  67. }
对一张透视图片变换回正面图的效果:

  1. int main(){
  2. Mat img=imread("boy.png");
  3. int img_height = img.rows;
  4. int img_width = img.cols;
  5. Mat img_trans = Mat::zeros(img_height,img_width,CV_8UC3);
  6. PerspectiveTransform tansform = PerspectiveTransform::quadrilateralToQuadrilateral(
  7. 0,0,
  8. img_width-1,0,
  9. 0,img_height-1,
  10. img_width-1,img_height-1,
  11. 150,250, // top left
  12. 771,0, // top right
  13. 0,1023,// bottom left
  14. 650,1023
  15. );
  16. vector<float> ponits;
  17. for(int i=0;i
  18. for(int j=0;j
  19. ponits.push_back(j);
  20. ponits.push_back(i);
  21. }
  22. }
  23. tansform.transformPoints(ponits);
  24. for(int i=0;i
  25. uchar* t= img_trans.ptr(i);
  26. for (int j=0;j
  27. int tmp = i*img_width+j;
  28. int x = ponits[tmp*2];
  29. int y = ponits[tmp*2+1];
  30. if(x<0||x>(img_width-1)||y<0||y>(img_height-1))
  31. continue;
  32. uchar* p = img.ptr(y);
  33. t[j*3] = p[x*3];
  34. t[j*3+1] = p[x*3+1];
  35. t[j*3+2] = p[x*3+2];
  36. }
  37. }
  38. imwrite("trans.png",img_trans);
  39. return 0;
  40. }




另外在OpenCV中也实现了基础的透视变换操作,有关函数使用请见下一篇: 【OpenCV】透视变换 Perspective Transformation(续)



(转载请注明作者和出处:http://iyenn.com/index/link?url=http://blog.csdn.net/xiaowei_cqu 未经允许请勿用于商业用途)



注:本文转载自blog.csdn.net的xiaowei_cqu的文章"https://xiaowei.blog.csdn.net/article/details/26471527"。版权归原作者所有,此博客不拥有其著作权,亦不承担相应法律责任。如有侵权,请联系我们删除。
复制链接
复制链接
相关推荐
发表评论
登录后才能发表评论和回复 注册

/ 登录

评论记录:

未查询到任何数据!
回复评论:

分类栏目

后端 (14832) 前端 (14280) 移动开发 (3760) 编程语言 (3851) Java (3904) Python (3298) 人工智能 (10119) AIGC (2810) 大数据 (3499) 数据库 (3945) 数据结构与算法 (3757) 音视频 (2669) 云原生 (3145) 云平台 (2965) 前沿技术 (2993) 开源 (2160) 小程序 (2860) 运维 (2533) 服务器 (2698) 操作系统 (2325) 硬件开发 (2492) 嵌入式 (2955) 微软技术 (2769) 软件工程 (2056) 测试 (2865) 网络空间安全 (2948) 网络与通信 (2797) 用户体验设计 (2592) 学习和成长 (2593) 搜索 (2744) 开发工具 (7108) 游戏 (2829) HarmonyOS (2935) 区块链 (2782) 数学 (3112) 3C硬件 (2759) 资讯 (2909) Android (4709) iOS (1850) 代码人生 (3043) 阅读 (2841)

热门文章

101
推荐
关于我们 隐私政策 免责声明 联系我们
Copyright © 2020-2025 蚁人论坛 (iYenn.com) All Rights Reserved.
Scroll to Top