[转]gluProject 和 gluUnproject 的详解

gluProject 和 gluUnproject 的详解

简介

三维空间中,经常需要将 3D 空间中的点转换到 2D(屏幕坐标),或者将 2D 点转换到 3D 空间中。当你使用 OpenGL 的时候,简单使用 gluProject() 和 gluUnproject() 函数就可以实现这个功能了。但这两个神奇的函数是怎样实现的,一直困扰着我,经过一番仔细研究,将自己的思路写在这里:

gluPorject()

先通过看代码,来一步一步分析它的数学原理吧!(其实代码是次要的,数学原理在这里才是关键所在!)这里的代码据说是赖在 mesa OpenGL 中的!

PS:这里为了更好理解,我修改了一下里面的代码,但没有兼顾效率的!

#include<GL/gl.h>

/*

* Transform a point (column vector) by a 4x4 matrix. I.e. out = m * in

* Input: m - the 4x4 matrix

* in - the 4x1 vector

* Output: out - the resulting 4x1 vector.

*/

static void transform_point(GLdouble out[4],

              const GLdouble m[16],

              const GLdouble in[4])

{

#define M(row,col) m[col* 4+ row]

out[0] = M(0, 0) * in[0] + M(0, 1) * in[1] + M(0, 2) * in[2] + M(0, 3) * in[3];

out[1] = M(1, 0) * in[0] + M(1, 1) * in[1] + M(1, 2) * in[2] + M(1, 3) * in[3];

out[2] = M(2, 0) * in[0] + M(2, 1) * in[1] + M(2, 2) * in[2] + M(2, 3) * in[3];

out[3] = M(3, 0) * in[0] + M(3, 1) * in[1] + M(3, 2) * in[2] + M(3, 3) * in[3];

#undef M

}

GLint gluProject(GLdouble objx,

         GLdouble objy,

         GLdouble objz,

         const GLdouble model[16],

         const GLdouble proj[16],

         const GLint viewport[4],

         GLdouble*winx,

         GLdouble*winy,

         GLdouble*winz) /* transformation matrix */

{

GLdouble objCoor[4];

GLdouble objProj[4], objModel[4];

/* initilise matrix and vector transform */

// 4x4 matrix must be multi to a 4 dimension vector( it a 1 x 4 matrix)

// so we need to put the original vertex to a 4D vector

objCoor[0] = objx;

objCoor[1] = objy;

objCoor[2] = objz;

objCoor[3] = 1.0;

// 由于原来的向量位于标准基向量 (1, 0, 0), (0, 1, 0), (0, 0, 1) 中,所以需要先转换到当前的模型矩阵中

transform_point(objModel, model, objCoor);

// 然后将模型矩阵中的顶点转换到投影矩阵所在坐标系的矩阵中

transform_point(objProj, proj, objModel);

// scale matrix

/* GLdouble scaleMat[4][4] =

{ {0.5, 0, 0, objPr0j[3]}, {0, 0.5, 0,

objProj[3]}, {0, 0, 0.5, objProj[3]}, {1, 1, 1,

1} }; GLdouble objProjTemp[4];

memcpy(objProjTemp, objProj, sizeof(objProjTemp);

transfrom_point(objProj, scaleMat, objProjTemp); */

/* or the result of normalized between -1 and 1 */

if(objProj[3]== 0.0)

return GL_FALSE;

objProj[0] /= objProj[3];

objProj[1] /= objProj[3];

objProj[2] /= objProj[3]; /* in screen coordinates */

// 由于投影矩阵投影在 [-1, 1] 之间,所以需要将转换后的投影坐标放置到 [0, 1] 之间

// 最后再在一个 offset 矩形中转换为屏幕坐标就可以了( viewport[4] 可以简单的认为一个 offset 矩形)

#define SCALE_FROM_0_TO_1(_pt) (((_pt)+ 1) / 2)

objProj[0] = SCALE_FROM_0_TO_1(objProj[0]);

objProj[1] = SCALE_FROM_0_TO_1(objProj[1]);

objProj[2] = SCALE_FROM_0_TO_1(objProj[2]);

#undef SCALE_FROM_0_TO_1

*winx= viewport[0] + objProj[0] * viewport[2];

*winy= viewport[1] + objProj[1] * viewport[3]; /* between 0 and 1 */

*winz= objProj[2];

return GL_TRUE;

}

基本的思路就是:

1、将输入的顶点,通过模型视图矩阵,变换到模型视图矩阵的坐标系中;

2、将模型视图矩阵中的顶点,再变换到投影矩阵中;

3、将顶点缩放到 [0, 1] 的映射区间中;

4、通过视口的位置和大小, 计算出当前 3D 顶点中的屏幕坐标( 2D 坐标)

gluUnproject

其实 gluUnproject 和 gluProject 是非常类似的, 代码我暂时没有去找,但我认为应该是这样的(其实就是 gluPorject 反过来的过程,只是有一些数学运算要注意一下) :

1、首先,需要将输入的顶点,通过视口变换到 [0, 1] 之间;

2、然后将顶点缩放到 [-1, 1] 之间,就是上面代码中的 scaleMat矩阵的逆矩阵

3、然后乘上投影矩阵的逆矩阵;

4、最后就是乘上模型视图矩阵的逆矩阵; gluProject,这里我暂时还没有验证,待我有空的时候再检查一下是否确实如此!呵呵

逆变换和模拟变换:http://book.51cto.com/art/201002/185504.htm

引用:

gluProject/gluUnProject don‘t modify these matrices. In the case of gluProject it applies the transformations (the modelview, projection, and viewport) the same way the standard pipeline does.

In the case of gluUnProject(4) they use the inverses of those matrices.

if v is the vertex, M is the modelview matrix, P is the projection matrix, and V is the viewport matrix:

v given in object space
gluProject returns V*P*M*v

v given in screen space
gluUnProject returns Mi*Pi*Vi*v
(where i represents the matrix inverse)

原文地址:https://www.cnblogs.com/MakeView660/p/10482988.html

时间: 2024-10-08 00:49:14

[转]gluProject 和 gluUnproject 的详解的相关文章

Spring事务管理(详解+实例)

写这篇博客之前我首先读了<Spring in action>,之后在网上看了一些关于Spring事务管理的文章,感觉都没有讲全,这里就将书上的和网上关于事务的知识总结一下,参考的文章如下: Spring事务机制详解 Spring事务配置的五种方式 Spring中的事务管理实例详解 1 初步理解 理解事务之前,先讲一个你日常生活中最常干的事:取钱. 比如你去ATM机取1000块钱,大体有两个步骤:首先输入密码金额,银行卡扣掉1000元钱:然后ATM出1000元钱.这两个步骤必须是要么都执行要么都

转载:DenseNet算法详解

原文连接:http://blog.csdn.net/u014380165/article/details/75142664 参考连接:http://blog.csdn.net/u012938704/article/details/53468483 本文这里仅当学习笔记使用,具体细节建议前往原文细度. 论文:Densely Connected Convolutional Networks 论文链接:https://arxiv.org/pdf/1608.06993.pdf 代码的github链接:h

MariaDB(MySQL)创建、删除、选择及数据类型使用详解

一.MariaDB简介(MySQL简介略过) MariaDB数据库管理系统是MySQL的一个分支,主要由开源社区在维护,采用GPL授权许可 MariaDB的目的是完全兼容MySQL,包括API和命令行,使之能轻松成为MySQL的代替品.在存储引擎方面,使用XtraDB(英语:XtraDB)来代替MySQL的InnoDB. MariaDB由MySQL的创始人Michael Widenius(英语:Michael Widenius)主导开发,他早前曾以10亿美元的价格,将自己创建的公司MySQL A

HttpServletResponse和HttpServletRequest详解

HttpServletResponse,HttpServletRequest详解 1.相关的接口 HttpServletRequest HttpServletRequest接口最常用的方法就是获得请求中的参数,这些参数一般是客户端表单中的数据.同时,HttpServletRequest接口可以获取由客户端传送的名称,也可以获取产生请求并且接收请求的服务器端主机名及IP地址,还可以获取客户端正在使用的通信协议等信息.下表是接口HttpServletRequest的常用方法. 说明:HttpServ

POSIX 线程详解(经典必看)

总共三部分: 第一部分:POSIX 线程详解                                   Daniel Robbins ([email protected]), 总裁/CEO, Gentoo Technologies, Inc.  2000 年 7 月 01 日 第二部分:通用线程:POSIX 线程详解,第 2部分       Daniel Robbins ([email protected]), 总裁/CEO, Gentoo Technologies, Inc.  20

.NET深入解析LINQ框架(五:IQueryable、IQueryProvider接口详解)

阅读目录: 1.环路执行对象模型.碎片化执行模型(假递归式调用) 2.N层对象执行模型(纵横向对比链式扩展方法) 3.LINQ查询表达式和链式查询方法其实都是空壳子 4.详细的对象结构图(对象的执行原理) 5.IQueryable<T>与IQueryProvider一对一的关系能否改成一对多的关系 6.完整的自定义查询 1]. 环路执行对象模型.碎片化执行模型(假递归式调用) 这个主题扯的可能有点远,但是它关系着整个LINQ框架的设计结构,至少在我还没有搞懂LINQ的本意之前,在我脑海里一直频

netstat状态详解

一.生产服务器netstat tcp连接状态................................................................................ 2 1.1生产服务器某个业务LVS负载均衡上连接状态数量............................................... 2 1.2生产服务器某个业务web上连接状态数量...............................................

详解go语言的array和slice 【二】

上一篇  详解go语言的array和slice [一]已经讲解过,array和slice的一些基本用法,使用array和slice时需要注意的地方,特别是slice需要注意的地方比较多.上一篇的最后讲解到创建新的slice时使用第三个索引来限制slice的容量,在操作新slice时,如果新slice的容量大于长度时,添加新元素依然后使源的相应元素改变.这一篇里我会讲解到如何避免这些问题,以及迭代.和做为方法参数方面的知识点. slice的长度和容量设置为同一个值 如果在创建新的slice时我们把

13.Linux键盘按键驱动 (详解)

版权声明:本文为博主原创文章,未经博主允许不得转载. 在上一节分析输入子系统内的intput_handler软件处理部分后,接下来我们开始写input_dev驱动 本节目标: 实现键盘驱动,让开发板的4个按键代表键盘中的L.S.空格键.回车键 1.先来介绍以下几个结构体使用和函数,下面代码中会用到 1)input_dev驱动设备结构体中常用成员如下: struct input_dev { void *private; const char *name; //设备名字 const char *ph