接着上一篇讲，前面已经展示了纹理把图片贴到立方体上，立方体表面是平面的，纹理坐标不需要设计，下面的两个例子都是基于球面贴图的，把鱼眼和全景图片贴到球面上，即实现曲面贴图。实际上，曲面贴图和平面贴图使用的纹理方面的知识没变化，主要难点在于纹理坐标设计，纹理坐标设计涉及到变换模型设计，需要一定的数学基础，下面例子中的两种贴图方法是比较简单的，网上的资料也非常多，学习本章例子的话，最好能推导例子中用的公式（都是简单的三角函数应用）。
1. 鱼眼球面贴图
一般的鱼眼图是一张正方形图（长方形也无所谓），有效的鱼眼图像都是中心的圆部分，四角部分都是无效的像素区域（一般是黑的），这种鱼眼图片映射到半球面上，即把圆面映射到半球面上，使用这种模型，可以做鱼眼矫正（效果并非最好，实现方便），鱼眼图片一般分前向（对应摄像机壁装）及俯视、仰视（对应摄像机顶装与底装）等几种拍摄方式，下面的例子采用的是仰视方式贴图，另外两种场景算法原理是一致的，会推导一种就能掌握这种贴图方法。
圆映射到半球面，推导方法如下：
1. 假定半球面圆心在原点，半球面在z轴正方形，圆面在x，y平面上（z=0），圆半径R=0.5（直径为1.0，方便对应纹理坐标范围）
2. 半球面上每个点（x，y，z）映射到圆面的坐标为（x，y，0），也就是把z坐标直接丢掉，即圆面是半球面在x，y平面（z=0）上的投影。
3. 所以z坐标可以由x，y坐标表示出来，圆面的（x，y）确定了，半球面的（x，y，z）的坐标也就确定了，反之亦然。
4. 圆面坐标（x，y）与纹理坐标的映射就很简单了，分别加上0.5即可，即把范围有-0.5~0.5映射到0~1.0之间。
绘制的时候，半球面的纹理坐标直接根据顶点坐标计算即可，不需要单独传纹理坐标。
下面的例子使用一张鱼眼图片上下对称映射成一个球体（沿用上篇球的画法），如想只画半球，顶点坐标设计成半球即可，这是从网上找的一张全景图，如下(使用以下源码，需把图片转成bmp格式的)：
源码如下：
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <GL/glew.h>
#include <GL/glut.h>
/* 该头文件请到第6章下载*/
#include "vmath.h"
using namespace vmath;
#pragma pack(2)
#define SPLIT_NUM 10
#define POINT_CNT (2*SPLIT_NUM + 1) /* 每条边对应的点数*/
#define QUAD_CNT (2*SPLIT_NUM) /* 每条边对应的正方形个数*/
#define FISHEYE_IMG "fisheye.bmp"
typedef unsigned char U8;
typedef unsigned short int U16;
typedef unsigned long int U32;
typedef struct
{
U16 bfType; /* windows下该值必需是x4D42,即字符'BM'*/
U32 bfSize; /* bmp文件大小,包含bmp文件头,信息头,调色板大小,数据大小*/
U16 bfReserved1; /* 保留，必须设置为*/
U16 bfReserved2; /* 保留，必须设置为*/
U32 bfOffBits; /* rgb数据相对文件头偏移量*/
} BitmapFileHeader;
typedef struct
{
U32 biSize; /* 信息头大小sizeof(BitmapInfoHeader) */
U32 biWidth; /* 图象的宽度,以象素为单位*/
U32 biHeight; /* 图象的高度，以象素为单位,正值倒向位图,负值正向位图*/
U16 biPlanes; /* 位面数,设为*/
U16 biBitCount; /* 位图位数,可以为,24,16,8,4,1 */
U32 biCompression; /* 说明图象数据压缩的类型,BI_RGB(0) BI_BITFIELDS(3)等*/
U32 biSizeImage; /* 图像数据大小,包括位图信息大小和数据大小*/
U32 biXPelsPerMeter;/* 水平分辨率,一般可填*/
U32 biYPelsPerMeter;/* 垂直分辨率,一般可填*/
U32 biClrUsed; /* 颜色索引使用数,一般填,表示都使用*/
U32 biClrImportant; /* 重要颜色索引数,一般填,表示都重要*/
} BitmapInfoHeader;
typedef struct
{
BitmapFileHeader bfHeader;
BitmapInfoHeader biHeader;
} BitmapInfo;
typedef struct
{
int width;
int height;
char * data;
} BitmapLoad;
typedef struct
{
GLfloat x;
GLfloat y;
GLfloat z;
} POINT_S;
static const GLchar * vertex_source =
"#version 330 core\n"
"uniform mat4 pos_matrix;\n"
"uniform mat4 face_matrix;\n"
"layout (location = 0) in vec3 in_position;\n"
"out vec2 tex_coord;\n"
"void main(void)\n"
"{\n"
" float dist = 2 * sqrt(pow(in_position.x,2)+pow(in_position.y,2)+pow(in_position.z,2));\n"
" vec4 vpos = vec4(in_position.x/dist,in_position.y/dist,in_position.z/dist,1.0);\n"
" vec4 tpos = face_matrix * vpos;\n"
" tex_coord = vec2(tpos.x+0.5,tpos.y+0.5);\n"
" gl_Position = pos_matrix * face_matrix * vpos;\n"
"}\n";
static const GLchar * frag_source =
"#version 330 core\n"
"in vec2 tex_coord;\n"
"out vec4 color;\n"
"uniform sampler2D tex;\n"
"void main(void)\n"
"{\n"
" color = texture(tex, tex_coord);\n"
"}\n";
void loadShader(GLuint program, GLuint type, const GLchar * source)
{
GLint status = 0;
const GLchar * shaderSource[] = {source};
GLuint shader = glCreateShader(type);
glShaderSource(shader, 1, shaderSource, 0);
glCompileShader(shader);
glGetShaderiv(shader, GL_COMPILE_STATUS, &status);
glAttachShader(program, shader);
}
BitmapLoad * LoadBmp(const char * filename, BitmapLoad * texture_image)
{
int pos = 0, dataSize = 0;
FILE * fp = NULL;
BitmapInfo bmpInfo;
fp = fopen(filename, "rb");
if (fp == NULL)
{
perror("fopen");
return NULL;
}
memset(&bmpInfo, 0, sizeof(BitmapInfo));
fread(&bmpInfo, 1, sizeof(BitmapInfo), fp);
printf("width:%lu height:%lu dataSize:%lu\n", bmpInfo.biHeader.biWidth, bmpInfo.biHeader.biHeight, bmpInfo.biHeader.biSizeImage);
if (bmpInfo.biHeader.biSizeImage != bmpInfo.biHeader.biWidth * bmpInfo.biHeader.biHeight * 3)
{
printf("biSizeImage size error! header:%d bmp size:%d\n", sizeof(BitmapInfo), bmpInfo.bfHeader.bfSize);
bmpInfo.biHeader.biSizeImage = bmpInfo.biHeader.biWidth * bmpInfo.biHeader.biHeight * 3;
}
dataSize = bmpInfo.biHeader.biSizeImage;
texture_image->data = (char *)malloc(dataSize);
if (texture_image->data == NULL)
{
fclose(fp);
return NULL;
}
fread(texture_image->data, 1, dataSize, fp);
texture_image->width = bmpInfo.biHeader.biWidth;
texture_image->height = bmpInfo.biHeader.biHeight;
pos = ftell(fp);
printf("bmp pos:%d\n", pos);
fclose(fp);
return texture_image;
}
GLuint vao, vbo, ebo,texture;
mat4 pos_matrix,face_matrix;
GLuint pos_matrix_idx,face_matrix_idx;
void display(void)
{
glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
glBindBuffer(GL_ARRAY_BUFFER, vao);
glUniformMatrix4fv(pos_matrix_idx, 1, GL_FALSE, pos_matrix);
face_matrix = vmath::translate(0.0f, 0.0f, 0.0f);
glUniformMatrix4fv(face_matrix_idx, 1, GL_FALSE, face_matrix);
glDrawElements(GL_QUADS, QUAD_CNT * QUAD_CNT * 4, GL_UNSIGNED_INT, (GLvoid *)(0));
face_matrix = vmath::rotate(90.0f, 1.0f, 0.0f, 0.0f);
glUniformMatrix4fv(face_matrix_idx, 1, GL_FALSE, face_matrix);
glDrawElements(GL_QUADS, QUAD_CNT * QUAD_CNT * 4, GL_UNSIGNED_INT, (GLvoid *)(0));
face_matrix = vmath::rotate(-90.0f, 1.0f, 0.0f, 0.0f);
glUniformMatrix4fv(face_matrix_idx, 1, GL_FALSE, face_matrix);
glDrawElements(GL_QUADS, QUAD_CNT * QUAD_CNT * 4, GL_UNSIGNED_INT, (GLvoid *)(0));
face_matrix = vmath::rotate(90.0f, 0.0f, 1.0f, 0.0f);
glUniformMatrix4fv(face_matrix_idx, 1, GL_FALSE, face_matrix);
glDrawElements(GL_QUADS, QUAD_CNT * QUAD_CNT * 4, GL_UNSIGNED_INT, (GLvoid *)(0));
face_matrix = vmath::rotate(-90.0f, 0.0f, 1.0f, 0.0f);
glUniformMatrix4fv(face_matrix_idx, 1, GL_FALSE, face_matrix);
glDrawElements(GL_QUADS, QUAD_CNT * QUAD_CNT * 4, GL_UNSIGNED_INT, (GLvoid *)(0));
face_matrix = vmath::rotate(180.0f, 1.0f, 0.0f, 0.0f);
glUniformMatrix4fv(face_matrix_idx, 1, GL_FALSE, face_matrix);
glDrawElements(GL_QUADS, QUAD_CNT * QUAD_CNT * 4, GL_UNSIGNED_INT, (GLvoid *)(0));
glutSwapBuffers();
}
/* 只分配一个面的顶点数据,其他面通过旋转绘制*/
POINT_S vertex_list[POINT_CNT][POINT_CNT];
/* 绘制索引,这里使用点的形式划线,注意空间分配*/
GLuint index_list[QUAD_CNT][QUAD_CNT][4];
void init(void)
{
int i = 0, j = 0, n = 0, cnt = 0;
POINT_S * p = NULL;
BitmapLoad textureImage = {0};
memset(vertex_list, 0, sizeof(vertex_list));
memset(index_list, 0, sizeof(index_list));
/* 添加立方体上表面坐标信息*/
for (i = 0; i < POINT_CNT ; i++)
{
for (j = 0; j < POINT_CNT; j++)
{
p = &vertex_list[i][j];
p->x = (-0.5f + j * 0.5f / SPLIT_NUM)*128;
p->y = (0.5f - i * 0.5f / SPLIT_NUM)*128;
p->z = 0.5f * 128;
}
}
/* 添加索引信息,这里按四边形绘制*/
for (i = 0; i < QUAD_CNT; i++)
{
for (j = 0; j < QUAD_CNT; j++)
{
index_list[i][j][0] = POINT_CNT * i + j;
index_list[i][j][1] = POINT_CNT * i + j + 1;
index_list[i][j][2] = POINT_CNT * (i + 1) + j + 1;
index_list[i][j][3] = POINT_CNT * (i + 1) + j;
}
}
GLuint program = glCreateProgram();
loadShader(program, GL_VERTEX_SHADER, vertex_source);
loadShader(program, GL_FRAGMENT_SHADER, frag_source);
glLinkProgram(program);
glUseProgram(program);
glGenBuffers(1, &vao);
glBindBuffer(GL_ARRAY_BUFFER, vao);
glGenBuffers(1, &vbo);
glBindBuffer(GL_ARRAY_BUFFER, vbo);
glBufferData(GL_ARRAY_BUFFER, sizeof(vertex_list), vertex_list, GL_STATIC_DRAW);
glVertexAttribPointer(0, 3, GL_FLOAT, GL_FALSE, 0, (GLvoid *)0);
glEnableVertexAttribArray(0);
glGenBuffers(1, &ebo);
glBindBuffer(GL_ELEMENT_ARRAY_BUFFER, ebo);
glBufferData(GL_ELEMENT_ARRAY_BUFFER, sizeof(index_list), index_list, GL_STATIC_DRAW);
glGenTextures(1, &texture);
if (LoadBmp(FISHEYE_IMG, &textureImage) != NULL)
{
glBindTexture(GL_TEXTURE_2D, texture);
glTexImage2D(GL_TEXTURE_2D, 0, 3, textureImage.width, textureImage.height, 0, GL_BGR, GL_UNSIGNED_BYTE, textureImage.data);
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);//GL_LINEAR GL_NEAREST
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);//
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_S, GL_CLAMP);
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_T, GL_CLAMP);
free(textureImage.data);
}
/* 定义旋转数组,主要存放当前朝向和旋转信息*/
face_matrix_idx = glGetUniformLocation(program, "face_matrix");
pos_matrix_idx = glGetUniformLocation(program, "pos_matrix");
pos_matrix = vmath::translate(0.0f, 0.0f, 0.0f);
glClearColor(0.5f, 0.5f, 1.0f, 1.0f);
glClearDepth(1.0);
glEnable(GL_DEPTH_TEST);
}
void keyboard(unsigned char key, int x, int y)
{
switch (key)
{
case '-':
pos_matrix *= vmath::scale(0.95f);
break;
case '=':
case '+':
pos_matrix *= vmath::scale(1.05f);
break;
case 'i':
pos_matrix = vmath::frustum(-0.5, 0.5, -0.5, 0.5, 0.4, 10.0);
break;
case 'o':
pos_matrix = vmath::translate(0.0f, 0.0f, 0.0f);
break;
default:
break;
}
glutPostRedisplay();
}
void specialKey(GLint key, GLint x, GLint y)
{
float step = 2.0f;
switch (key)
{
case GLUT_KEY_UP:
pos_matrix *= vmath::rotate(step, 1.0f, 0.0f, 0.0f);
break;
case GLUT_KEY_DOWN:
pos_matrix *= vmath::rotate(-1.0f * step, 1.0f, 0.0f, 0.0f);
break;
case GLUT_KEY_LEFT:
pos_matrix *= vmath::rotate(step, 0.0f, 1.0f, 0.0f);
break;
case GLUT_KEY_RIGHT:
pos_matrix *= vmath::rotate(-1.0f * step, 0.0f, 1.0f, 0.0f);
break;
case GLUT_KEY_PAGE_UP :
pos_matrix *= vmath::rotate(step, 0.0f, 0.0f, 1.0f);
break;
case GLUT_KEY_PAGE_DOWN:
pos_matrix *= vmath::rotate(-1.0f * step, 0.0f, 0.0f, 1.0f);
break;
default:
break;
}
glutPostRedisplay();
}
int main(int argc, char * argv[])
{
glutInit(&argc, argv);
glutInitDisplayMode(GLUT_RGBA | GLUT_DOUBLE | GLUT_DEPTH);
glutInitWindowPosition(200, 200);
glutInitWindowSize(400, 400);
glutCreateWindow("FishEye");
glewInit();
init();
glutDisplayFunc(display);
glutKeyboardFunc(keyboard);
glutSpecialFunc(specialKey);
glutMainLoop();
return 0;
}
效果如下：
2. 全景球面贴图
全景图和鱼眼图不一样，可以看作是球面映射的两个相反的方面，鱼眼是半球面投影，而全景图是曲面展开，所以鱼眼图像是可以转成全景图的。这里从网上找了一张全景图（世界地图），宽高比例是2：1，为什么是2：1呢？原因是：
1. 全景图是按球面展开，球的圆周（赤道）展开为宽，经线展开为高。
2. 圆周长度为2πR，经线长为πR，所以是2：1。
注意：这种是等弧长映射，球面经线方向的点和全景图高度方向的点是一一对应的，而不是简单的认为全景图高度等于球体的高度。
纹理坐标映射推导过程：
1. 与上例一样，球面半径为0.5，圆心在原点。
2. 画一个外接圆柱体，半径0.5，总高度为πR，即π/2，把全景图映射到圆柱面上。
3. 横坐标的关系是，宽度=水平弧长，最终总宽度=圆周。
4. 纵坐标的关系是，高度=垂直弧长，总高度=半圆周。（所以叫等弧长映射）
5. 有了上面两个关系量，可以得到球面坐标和纹理坐标的关系（利用好等弧长这个条件即可）。
上面的过程比较难说清楚，如果理解了映射过程（平面->柱面->球面），自己推导比看人家讲解的快（只使用了最基本的三角和反三角函数），先把网上找的图贴一下(使用以下源码，需把图片转成bmp格式的)。
源码如下：
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <GL/glew.h>
#include <GL/glut.h>
/* 该头文件请到第6章下载*/
#include "vmath.h"
using namespace vmath;
#pragma pack(2)
#define SPLIT_NUM 10
#define POINT_CNT (2*SPLIT_NUM + 1) /* 每条边对应的点数*/
#define QUAD_CNT (2*SPLIT_NUM) /* 每条边对应的正方形个数*/
#define FISHEYE_IMG "world2.bmp"
typedef unsigned char U8;
typedef unsigned short int U16;
typedef unsigned long int U32;
typedef struct
{
U16 bfType; /* windows下该值必需是x4D42,即字符'BM'*/
U32 bfSize; /* bmp文件大小,包含bmp文件头,信息头,调色板大小,数据大小*/
U16 bfReserved1; /* 保留，必须设置为*/
U16 bfReserved2; /* 保留，必须设置为*/
U32 bfOffBits; /* rgb数据相对文件头偏移量*/
} BitmapFileHeader;
typedef struct
{
U32 biSize; /* 信息头大小sizeof(BitmapInfoHeader) */
U32 biWidth; /* 图象的宽度,以象素为单位*/
U32 biHeight; /* 图象的高度，以象素为单位,正值倒向位图,负值正向位图*/
U16 biPlanes; /* 位面数,设为*/
U16 biBitCount; /* 位图位数,可以为,24,16,8,4,1 */
U32 biCompression; /* 说明图象数据压缩的类型,BI_RGB(0) BI_BITFIELDS(3)等*/
U32 biSizeImage; /* 图像数据大小,包括位图信息大小和数据大小*/
U32 biXPelsPerMeter;/* 水平分辨率,一般可填*/
U32 biYPelsPerMeter;/* 垂直分辨率,一般可填*/
U32 biClrUsed; /* 颜色索引使用数,一般填,表示都使用*/
U32 biClrImportant; /* 重要颜色索引数,一般填,表示都重要*/
} BitmapInfoHeader;
typedef struct
{
BitmapFileHeader bfHeader;
BitmapInfoHeader biHeader;
} BitmapInfo;
typedef struct
{
int width;
int height;
char * data;
} BitmapLoad;
typedef struct
{
GLfloat x;
GLfloat y;
GLfloat z;
} POINT_S;
static const GLchar * vertex_source =
"#version 330 core\n"
"uniform mat4 pos_matrix;\n"
"uniform mat4 face_matrix;\n"
"layout (location = 0) in vec3 in_position;\n"
"out vec2 tex_coord;\n"
"void main(void)\n"
"{\n"
" float dist = 2 * sqrt(pow(in_position.x,2)+pow(in_position.y,2)+pow(in_position.z,2));\n"
" vec4 vpos = vec4(in_position.x/dist,in_position.y/dist,in_position.z/dist,1.0);\n"
" vec4 tpos = face_matrix * vpos;\n"
" float theta = atan(tpos.y/tpos.x);\n"
" if( tpos.x < 0.0) theta+=3.14; \n"
" else if(tpos.y<0.0) theta+=6.28; \n"
" float ty = asin(2*tpos.z)/3.14;\n"
" tex_coord = vec2(theta/3.14/2,ty * (-1.0)+0.5);\n"
" gl_Position = pos_matrix * face_matrix * vpos;\n"
"}\n";
static const GLchar * frag_source =
"#version 330 core\n"
"in vec2 tex_coord;\n"
"out vec4 color;\n"
"uniform sampler2D tex;\n"
"void main(void)\n"
"{\n"
" color = texture(tex, tex_coord);\n"
"}\n";
void loadShader(GLuint program, GLuint type, const GLchar * source)
{
GLint status = 0;
const GLchar * shaderSource[] = {source};
GLuint shader = glCreateShader(type);
glShaderSource(shader, 1, shaderSource, 0);
glCompileShader(shader);
glGetShaderiv(shader, GL_COMPILE_STATUS, &status);
glAttachShader(program, shader);
}
BitmapLoad * LoadBmp(const char * filename, BitmapLoad * texture_image)
{
int pos = 0, dataSize = 0;
FILE * fp = NULL;
BitmapInfo bmpInfo;
fp = fopen(filename, "rb");
if (fp == NULL)
{
perror("fopen");
return NULL;
}
memset(&bmpInfo, 0, sizeof(BitmapInfo));
fread(&bmpInfo, 1, sizeof(BitmapInfo), fp);
printf("width:%lu height:%lu dataSize:%lu\n", bmpInfo.biHeader.biWidth, bmpInfo.biHeader.biHeight, bmpInfo.biHeader.biSizeImage);
if (bmpInfo.biHeader.biSizeImage != bmpInfo.biHeader.biWidth * bmpInfo.biHeader.biHeight * 3)
{
printf("biSizeImage size error! header:%d bmp size:%d\n", sizeof(BitmapInfo), bmpInfo.bfHeader.bfSize);
bmpInfo.biHeader.biSizeImage = bmpInfo.biHeader.biWidth * bmpInfo.biHeader.biHeight * 3;
}
dataSize = bmpInfo.biHeader.biSizeImage;
texture_image->data = (char *)malloc(dataSize);
if (texture_image->data == NULL)
{
fclose(fp);
return NULL;
}
fread(texture_image->data, 1, dataSize, fp);
texture_image->width = bmpInfo.biHeader.biWidth;
texture_image->height = bmpInfo.biHeader.biHeight;
pos = ftell(fp);
printf("bmp pos:%d\n", pos);
fclose(fp);
return texture_image;
}
GLuint vao, vbo, ebo,texture;
mat4 pos_matrix,face_matrix;
GLuint pos_matrix_idx,face_matrix_idx;
void display(void)
{
glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
glBindBuffer(GL_ARRAY_BUFFER, vao);
glUniformMatrix4fv(pos_matrix_idx, 1, GL_FALSE, pos_matrix);
face_matrix = vmath::translate(0.0f, 0.0f, 0.0f);
glUniformMatrix4fv(face_matrix_idx, 1, GL_FALSE, face_matrix);
glDrawElements(GL_QUADS, QUAD_CNT * QUAD_CNT * 4, GL_UNSIGNED_INT, (GLvoid *)(0));
face_matrix = vmath::rotate(90.0f, 1.0f, 0.0f, 0.0f);
glUniformMatrix4fv(face_matrix_idx, 1, GL_FALSE, face_matrix);
glDrawElements(GL_QUADS, QUAD_CNT * QUAD_CNT * 4, GL_UNSIGNED_INT, (GLvoid *)(0));
face_matrix = vmath::rotate(-90.0f, 1.0f, 0.0f, 0.0f);
glUniformMatrix4fv(face_matrix_idx, 1, GL_FALSE, face_matrix);
glDrawElements(GL_QUADS, QUAD_CNT * QUAD_CNT * 4, GL_UNSIGNED_INT, (GLvoid *)(0));
face_matrix = vmath::rotate(90.0f, 0.0f, 1.0f, 0.0f);
glUniformMatrix4fv(face_matrix_idx, 1, GL_FALSE, face_matrix);
glDrawElements(GL_QUADS, QUAD_CNT * QUAD_CNT * 4, GL_UNSIGNED_INT, (GLvoid *)(0));
face_matrix = vmath::rotate(-90.0f, 0.0f, 1.0f, 0.0f);
glUniformMatrix4fv(face_matrix_idx, 1, GL_FALSE, face_matrix);
glDrawElements(GL_QUADS, QUAD_CNT * QUAD_CNT * 4, GL_UNSIGNED_INT, (GLvoid *)(0));
face_matrix = vmath::rotate(180.0f, 1.0f, 0.0f, 0.0f);
glUniformMatrix4fv(face_matrix_idx, 1, GL_FALSE, face_matrix);
glDrawElements(GL_QUADS, QUAD_CNT * QUAD_CNT * 4, GL_UNSIGNED_INT, (GLvoid *)(0));
glutSwapBuffers();
}
/* 只分配一个面的顶点数据,其他面通过旋转绘制*/
POINT_S vertex_list[POINT_CNT][POINT_CNT];
/* 绘制索引,这里使用点的形式划线,注意空间分配*/
GLuint index_list[QUAD_CNT][QUAD_CNT][4];
void init(void)
{
int i = 0, j = 0, n = 0, cnt = 0;
POINT_S * p = NULL;
BitmapLoad textureImage = {0};
memset(vertex_list, 0, sizeof(vertex_list));
memset(index_list, 0, sizeof(index_list));
/* 添加立方体上表面坐标信息*/
for (i = 0; i < POINT_CNT ; i++)
{
for (j = 0; j < POINT_CNT; j++)
{
p = &vertex_list[i][j];
p->x = (-0.5f + j * 0.5f / SPLIT_NUM)*128;
p->y = (0.5f - i * 0.5f / SPLIT_NUM)*128;
p->z = 0.5f * 128;
}
}
/* 添加索引信息,这里按四边形绘制*/
for (i = 0; i < QUAD_CNT; i++)
{
for (j = 0; j < QUAD_CNT; j++)
{
index_list[i][j][0] = POINT_CNT * i + j;
index_list[i][j][1] = POINT_CNT * i + j + 1;
index_list[i][j][2] = POINT_CNT * (i + 1) + j + 1;
index_list[i][j][3] = POINT_CNT * (i + 1) + j;
}
}
GLuint program = glCreateProgram();
loadShader(program, GL_VERTEX_SHADER, vertex_source);
loadShader(program, GL_FRAGMENT_SHADER, frag_source);
glLinkProgram(program);
glUseProgram(program);
glGenBuffers(1, &vao);
glBindBuffer(GL_ARRAY_BUFFER, vao);
glGenBuffers(1, &vbo);
glBindBuffer(GL_ARRAY_BUFFER, vbo);
glBufferData(GL_ARRAY_BUFFER, sizeof(vertex_list), vertex_list, GL_STATIC_DRAW);
glVertexAttribPointer(0, 3, GL_FLOAT, GL_FALSE, 0, (GLvoid *)0);
glEnableVertexAttribArray(0);
glGenBuffers(1, &ebo);
glBindBuffer(GL_ELEMENT_ARRAY_BUFFER, ebo);
glBufferData(GL_ELEMENT_ARRAY_BUFFER, sizeof(index_list), index_list, GL_STATIC_DRAW);
glGenTextures(1, &texture);
if (LoadBmp(FISHEYE_IMG, &textureImage) != NULL)
{
glBindTexture(GL_TEXTURE_2D, texture);
glTexImage2D(GL_TEXTURE_2D, 0, 3, textureImage.width, textureImage.height, 0, GL_BGR, GL_UNSIGNED_BYTE, textureImage.data);
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);//GL_LINEAR GL_NEAREST
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);//
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_S, GL_CLAMP);
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_T, GL_CLAMP);
free(textureImage.data);
}
/* 定义旋转数组,主要存放当前朝向和旋转信息*/
face_matrix_idx = glGetUniformLocation(program, "face_matrix");
pos_matrix_idx = glGetUniformLocation(program, "pos_matrix");
pos_matrix = vmath::translate(0.0f, 0.0f, 0.0f);
glClearColor(0.5f, 0.5f, 1.0f, 1.0f);
glClearDepth(1.0);
glEnable(GL_DEPTH_TEST);
}
void keyboard(unsigned char key, int x, int y)
{
switch (key)
{
case '-':
pos_matrix *= vmath::scale(0.95f);
break;
case '=':
case '+':
pos_matrix *= vmath::scale(1.05f);
break;
case 'i':
pos_matrix = vmath::frustum(-0.5, 0.5, -0.5, 0.5, 0.4, 10.0);
break;
case 'o':
pos_matrix = vmath::translate(0.0f, 0.0f, 0.0f);
break;
default:
break;
}
glutPostRedisplay();
}
void specialKey(GLint key, GLint x, GLint y)
{
float step = 2.0f;
switch (key)
{
case GLUT_KEY_UP:
pos_matrix *= vmath::rotate(step, 1.0f, 0.0f, 0.0f);
break;
case GLUT_KEY_DOWN:
pos_matrix *= vmath::rotate(-1.0f * step, 1.0f, 0.0f, 0.0f);
break;
case GLUT_KEY_LEFT:
pos_matrix *= vmath::rotate(step, 0.0f, 1.0f, 0.0f);
break;
case GLUT_KEY_RIGHT:
pos_matrix *= vmath::rotate(-1.0f * step, 0.0f, 1.0f, 0.0f);
break;
case GLUT_KEY_PAGE_UP :
pos_matrix *= vmath::rotate(step, 0.0f, 0.0f, 1.0f);
break;
case GLUT_KEY_PAGE_DOWN:
pos_matrix *= vmath::rotate(-1.0f * step, 0.0f, 0.0f, 1.0f);
break;
default:
break;
}
glutPostRedisplay();
}
int main(int argc, char * argv[])
{
glutInit(&argc, argv);
glutInitDisplayMode(GLUT_RGBA | GLUT_DOUBLE | GLUT_DEPTH);
glutInitWindowPosition(200, 200);
glutInitWindowSize(400, 400);
glutCreateWindow("WorldMap");
glewInit();
init();
glutDisplayFunc(display);
glutKeyboardFunc(keyboard);
glutSpecialFunc(specialKey);
glutMainLoop();
return 0;
}
效果如下：
说明：上面的贴图方式并不完善，大家可以看到接缝处有一条黑线，原因是顶点模型选择的问题，立方体逼近的方法，在上例鱼眼贴图中效果会好一些，但本例的全景贴图，球面顶点模型用经纬法的更好，否则没办法处理这种黑线问题，黑线产生的原因是，通过顶点坐标计算出的纹理坐标，在接缝线左（靠近1.0）右（靠近0.0）差不多相差1.0，假如有一个四边形跨了接缝，实际跨度很小，只跨了一条缝，但顶点坐标跨度范围会非常大，如0.1~0.9，这个小四边形差不多把整个纹理都填充在里面，所以表现为一条黑线，实际上我们要的效果是0.0~0.1，0.9~1.0这样的填充，这就必须要求接缝线是直线（意思是经度相同的线，球面上的直线），立方体逼近方法没办法保证接缝线是直线，所以很难去除黑线问题，而经纬法画的球体，刚好有这种特性。大家可以试着把顶点模型改成经纬法的球面，来看效果。
3. 鱼眼矫正
鱼眼矫正的算法有很多种（网上有很多论文），上面例子中的映射方法是鱼眼矫正方法的一种，矫正的效果一般，但胜在原理简单，速度快，有产品使用该模型进行实时视频矫正。
再看一下上面两个例子，是有一些联系的，全景图可以和鱼眼图相互转化，即全景图片<-->柱面映射<-->球面映射<-->鱼眼图片，理解这中间的映射方式，就可以做鱼眼矫正了，鱼眼矫正有转全景、云台、柱面、分割、缩放等多种操作方式，在这种方法中，最关键的就是球面模型建立与处理。鱼眼矫正产品化有一定难度，但原理不复杂，可以用来练手，如把视点定在原点，对球面做投影视图变换，调整角度，就可以实现类似云台的功能。
4. 关于调试
shader程序的调试很麻烦，涉及到模型变换时，稍有不对，得出的结果跟预想的完全不一样。shader调试没有很直接的方法，连打印调试都做不到，这正是并行运算的特点，即使可以打印，打印的数据也是乱序，且数据会太多，网上提到的颜色调试、调试器调试不直观也比较麻烦，个人主要调试方法为：
1. 分析与尝试，屏蔽无关顶点，只保留关注的顶点，方便分析，只需要在顶点着色器中把正常的顶点过滤掉即可（判断顶点坐标为不关注顶点时，main函数空函数）。
2. 逻辑部分代码移到CPU上跑，调试OK后，再移到shader里，由于shader语法与C语法类似，移植还是比较方便的，虽然麻烦但有效。
3. 调纹理的时候使用测试图，调好过后再替换成效果图，测试图带定位信息，更容易看清楚顶点和纹理坐标的对应关系。
下面是调试时使用的测试图（最好看图就能知道对应纹理的位置，格子，环状图都可以），本图做的比较烂，大家可以参考。
