WebGL的颜色渲染-渲染一张DEM(数字高程模型)
目录
1. 具体实例
通过WebGL，可以渲染生成DEM(数字高程模型)。DEM(数字高程模型)是网格点组成的模型，每个点都有x，y，z值；x，y根据一定的间距组成网格状，同时根据z值的高低来选定每个点的颜色RGB。通过这个例子可以熟悉WebGL颜色渲染的过程。
2. 解决方案
1) DEM数据.XYZ文件
这里使用的DEM文件的数据组织如下，如下图所示。
其中每一行表示一个点，前三个数值表示位置XYZ，后三个数值表示颜色RGB。
2) showDEM.html
<!DOCTYPE html>
<html>
<head>
<meta charset="UTF-8">
<title> 显示地形 </title>
<script src="lib/webgl-utils.js"></script>
<script src="lib/webgl-debug.js"></script>
<script src="lib/cuon-utils.js"></script>
<script src="lib/cuon-matrix.js"></script>
<script src="showDEM.js"></script>
</head>
<body>
<div><input type = 'file' id = 'demFile' ></div>
<!-- <div><textarea id="output" rows="300" cols="200"></textarea></div> -->
<div>
<canvas id ="demCanvas" width="600" height="600">
请使用支持WebGL的浏览器
</canvas>
</div>
</body>
</html>
3) showDEM.js
// Vertex shader program
var VSHADER_SOURCE =
//'precision highp float;\n' +
'attribute vec4 a_Position;\n' +
'attribute vec4 a_Color;\n' +
'uniform mat4 u_MvpMatrix;\n' +
'varying vec4 v_Color;\n' +
'void main() {\n' +
' gl_Position = u_MvpMatrix * a_Position;\n' +
' v_Color = a_Color;\n' +
'}\n';
// Fragment shader program
var FSHADER_SOURCE =
'#ifdef GL_ES\n' +
'precision mediump float;\n' +
'#endif\n' +
'varying vec4 v_Color;\n' +
'void main() {\n' +
' gl_FragColor = v_Color;\n' +
'}\n';
//
var col = 89; //DEM宽
var row = 245; //DEM高
// Current rotation angle ([x-axis, y-axis] degrees)
var currentAngle = [0.0, 0.0];
//当前lookAt()函数初始视点的高度
var eyeHight = 2000.0;
//setPerspective()远截面
var far = 3000;
//
window.onload = function () {
var demFile = document.getElementById('demFile');
if (!demFile) {
console.log("Error!");
return;
}
//demFile.onchange = openFile(event);
demFile.addEventListener("change", function (event) {
//判断浏览器是否支持FileReader接口
if (typeof FileReader == 'undefined') {
console.log("你的浏览器不支持FileReader接口！");
return;
}
//
var reader = new FileReader();
reader.onload = function () {
if (reader.result) {
//
var stringlines = reader.result.split("\n");
verticesColors = new Float32Array(stringlines.length * 6);
//
var pn = 0;
var ci = 0;
for (var i = 0; i < stringlines.length; i++) {
if (!stringlines[i]) {
continue;
}
var subline = stringlines[i].split(',');
if (subline.length != 6) {
console.log("错误的文件格式！");
return;
}
for (var j = 0; j < subline.length; j++) {
verticesColors[ci] = parseFloat(subline[j]);
ci++;
}
pn++;
}
if (ci < 3) {
console.log("错误的文件格式！");
}
//
var minX = verticesColors[0];
var maxX = verticesColors[0];
var minY = verticesColors[1];
var maxY = verticesColors[1];
var minZ = verticesColors[2];
var maxZ = verticesColors[2];
for (var i = 0; i < pn; i++) {
minX = Math.min(minX, verticesColors[i * 6]);
maxX = Math.max(maxX, verticesColors[i * 6]);
minY = Math.min(minY, verticesColors[i * 6 + 1]);
maxY = Math.max(maxY, verticesColors[i * 6 + 1]);
minZ = Math.min(minZ, verticesColors[i * 6 + 2]);
maxZ = Math.max(maxZ, verticesColors[i * 6 + 2]);
}
//包围盒中心
var cx = (minX + maxX) / 2.0;
var cy = (minY + maxY) / 2.0;
var cz = (minZ + maxZ) / 2.0;
//根据视点高度算出setPerspective()函数的合理角度
var fovy = (maxY - minY) / 2.0 / eyeHight;
fovy = 180.0 / Math.PI * Math.atan(fovy) * 2;
startDraw(verticesColors, cx, cy, cz, fovy);
}
};
//
var input = event.target;
reader.readAsText(input.files[0]);
});
}
function startDraw(verticesColors, cx, cy, cz, fovy) {
// Retrieve <canvas> element
var canvas = document.getElementById('demCanvas');
// Get the rendering context for WebGL
var gl = getWebGLContext(canvas);
if (!gl) {
console.log('Failed to get the rendering context for WebGL');
return;
}
// Initialize shaders
if (!initShaders(gl, VSHADER_SOURCE, FSHADER_SOURCE)) {
console.log('Failed to intialize shaders.');
return;
}
// Set the vertex coordinates and color (the blue triangle is in the front)
n = initVertexBuffers(gl, verticesColors); //, verticesColors, n
if (n < 0) {
console.log('Failed to set the vertex information');
return;
}
// Get the storage location of u_MvpMatrix
var u_MvpMatrix = gl.getUniformLocation(gl.program, 'u_MvpMatrix');
if (!u_MvpMatrix) {
console.log('Failed to get the storage location of u_MvpMatrix');
return;
}
// Register the event handler
initEventHandlers(canvas);
// Specify the color for clearing <canvas>
gl.clearColor(0, 0, 0, 1);
gl.enable(gl.DEPTH_TEST);
// Start drawing
var tick = function () {
//setPerspective()宽高比
var aspect = canvas.width / canvas.height;
//
draw(gl, n, aspect, cx, cy, cz, fovy, u_MvpMatrix);
requestAnimationFrame(tick, canvas);
};
tick();
}
//
function initEventHandlers(canvas) {
var dragging = false; // Dragging or not
var lastX = -1, lastY = -1; // Last position of the mouse
// Mouse is pressed
canvas.onmousedown = function (ev) {
var x = ev.clientX;
var y = ev.clientY;
// Start dragging if a moue is in <canvas>
var rect = ev.target.getBoundingClientRect();
if (rect.left <= x && x < rect.right && rect.top <= y && y < rect.bottom) {
lastX = x;
lastY = y;
dragging = true;
}
};
//鼠标离开时
canvas.onmouseleave = function (ev) {
dragging = false;
};
// Mouse is released
canvas.onmouseup = function (ev) {
dragging = false;
};
// Mouse is moved
canvas.onmousemove = function (ev) {
var x = ev.clientX;
var y = ev.clientY;
if (dragging) {
var factor = 100 / canvas.height; // The rotation ratio
var dx = factor * (x - lastX);
var dy = factor * (y - lastY);
// Limit x-axis rotation angle to -90 to 90 degrees
//currentAngle[0] = Math.max(Math.min(currentAngle[0] + dy, 90.0), -90.0);
currentAngle[0] = currentAngle[0] + dy;
currentAngle[1] = currentAngle[1] + dx;
}
lastX = x, lastY = y;
};
//鼠标缩放
canvas.onmousewheel = function (event) {
var lastHeight = eyeHight;
if (event.wheelDelta > 0) {
eyeHight = Math.max(1, eyeHight - 80);
} else {
eyeHight = eyeHight + 80;
}
far = far + eyeHight - lastHeight;
};
}
function draw(gl, n, aspect, cx, cy, cz, fovy, u_MvpMatrix) {
//模型矩阵
var modelMatrix = new Matrix4();
modelMatrix.rotate(currentAngle[0], 1.0, 0.0, 0.0); // Rotation around x-axis
modelMatrix.rotate(currentAngle[1], 0.0, 1.0, 0.0); // Rotation around y-axis
modelMatrix.translate(-cx, -cy, -cz);
//视图矩阵
var viewMatrix = new Matrix4();
viewMatrix.lookAt(0, 0, eyeHight, 0, 0, 0, 0, 1, 0);
//投影矩阵
var projMatrix = new Matrix4();
projMatrix.setPerspective(fovy, aspect, 10, far);
//模型视图投影矩阵
var mvpMatrix = new Matrix4();
mvpMatrix.set(projMatrix).multiply(viewMatrix).multiply(modelMatrix);
// Pass the model view projection matrix to u_MvpMatrix
gl.uniformMatrix4fv(u_MvpMatrix, false, mvpMatrix.elements);
// Clear color and depth buffer
gl.clear(gl.COLOR_BUFFER_BIT | gl.DEPTH_BUFFER_BIT);
// Draw the cube
gl.drawElements(gl.TRIANGLES, n, gl.UNSIGNED_SHORT, 0);
}
function initVertexBuffers(gl, verticesColors) {
//DEM的一个网格是由两个三角形组成的
// 0------1 1
// | |
// | |
// col col------col+1
var indices = new Uint16Array((row - 1) * (col - 1) * 6);
var ci = 0;
for (var yi = 0; yi < row - 1; yi++) {
for (var xi = 0; xi < col - 1; xi++) {
indices[ci * 6] = yi * col + xi;
indices[ci * 6 + 1] = (yi + 1) * col + xi;
indices[ci * 6 + 2] = yi * col + xi + 1;
indices[ci * 6 + 3] = (yi + 1) * col + xi;
indices[ci * 6 + 4] = (yi + 1) * col + xi + 1;
indices[ci * 6 + 5] = yi * col + xi + 1;
ci++;
}
}
//创建缓冲区对象
var vertexColorBuffer = gl.createBuffer();
var indexBuffer = gl.createBuffer();
if (!vertexColorBuffer || !indexBuffer) {
return -1;
}
// 将缓冲区对象绑定到目标
gl.bindBuffer(gl.ARRAY_BUFFER, vertexColorBuffer);
// 向缓冲区对象中写入数据
gl.bufferData(gl.ARRAY_BUFFER, verticesColors, gl.STATIC_DRAW);
//
var FSIZE = verticesColors.BYTES_PER_ELEMENT;
// 向缓冲区对象分配a_Position变量
var a_Position = gl.getAttribLocation(gl.program, 'a_Position');
if (a_Position < 0) {
console.log('Failed to get the storage location of a_Position');
return -1;
}
gl.vertexAttribPointer(a_Position, 3, gl.FLOAT, false, FSIZE * 6, 0);
//开启a_Position变量
gl.enableVertexAttribArray(a_Position);
// 向缓冲区对象分配a_Color变量
var a_Color = gl.getAttribLocation(gl.program, 'a_Color');
if (a_Color < 0) {
console.log('Failed to get the storage location of a_Color');
return -1;
}
gl.vertexAttribPointer(a_Color, 3, gl.FLOAT, false, FSIZE * 6, FSIZE * 3);
//开启a_Color变量
gl.enableVertexAttribArray(a_Color);
// 写入并绑定顶点数组的索引值
gl.bindBuffer(gl.ELEMENT_ARRAY_BUFFER, indexBuffer);
gl.bufferData(gl.ELEMENT_ARRAY_BUFFER, indices, gl.STATIC_DRAW);
return indices.length;
}
4) 运行结果
用chrome打开showDEM.html，选择DEM文件，界面就会显示DEM的渲染效果：
3. 详细讲解
1) 读取文件
程序的第一步是通过JS的FileReader（）函数读取DEM文件，在其回调函数中读取到数组verticesColors中，它包含了位置和颜色信息。读取完成后调用绘制函数startDraw（）。
//
var reader = new FileReader();
reader.onload = function () {
if (reader.result) {
//
var stringlines = reader.result.split("\n");
verticesColors = new Float32Array(stringlines.length * 6);
//
var pn = 0;
var ci = 0;
for (var i = 0; i < stringlines.length; i++) {
if (!stringlines[i]) {
continue;
}
var subline = stringlines[i].split(',');
if (subline.length != 6) {
console.log("错误的文件格式！");
return;
}
for (var j = 0; j < subline.length; j++) {
verticesColors[ci] = parseFloat(subline[j]);
ci++;
}
pn++;
}
if (ci < 3) {
console.log("错误的文件格式！");
}
//
var minX = verticesColors[0];
var maxX = verticesColors[0];
var minY = verticesColors[1];
var maxY = verticesColors[1];
var minZ = verticesColors[2];
var maxZ = verticesColors[2];
for (var i = 0; i < pn; i++) {
minX = Math.min(minX, verticesColors[i * 6]);
maxX = Math.max(maxX, verticesColors[i * 6]);
minY = Math.min(minY, verticesColors[i * 6 + 1]);
maxY = Math.max(maxY, verticesColors[i * 6 + 1]);
minZ = Math.min(minZ, verticesColors[i * 6 + 2]);
maxZ = Math.max(maxZ, verticesColors[i * 6 + 2]);
}
//包围盒中心
var cx = (minX + maxX) / 2.0;
var cy = (minY + maxY) / 2.0;
var cz = (minZ + maxZ) / 2.0;
//根据视点高度算出setPerspective()函数的合理角度
var fovy = (maxY - minY) / 2.0 / eyeHight;
fovy = 180.0 / Math.PI * Math.atan(fovy) * 2;
startDraw(verticesColors, cx, cy, cz, fovy);
}
};
//
var input = event.target;
reader.readAsText(input.files[0]);
2) 绘制函数
绘制DEM跟绘制一个简单三角形的步骤是差不多的：
- 获取WebGL环境。
- 初始化shaders，构建着色器。
- 初始化顶点数组，分配到缓冲对象。
- 绑定鼠标键盘事件，设置模型视图投影变换矩阵。
- 在重绘函数中调用WebGL函数绘制。
其中最关键的步骤是第三步，初始化顶点数组initVertexBuffers()。
function startDraw(verticesColors, cx, cy, cz, fovy) {
// Retrieve <canvas> element
var canvas = document.getElementById('demCanvas');
// Get the rendering context for WebGL
var gl = getWebGLContext(canvas);
if (!gl) {
console.log('Failed to get the rendering context for WebGL');
return;
}
// Initialize shaders
if (!initShaders(gl, VSHADER_SOURCE, FSHADER_SOURCE)) {
console.log('Failed to intialize shaders.');
return;
}
// Set the vertex coordinates and color (the blue triangle is in the front)
n = initVertexBuffers(gl, verticesColors); //, verticesColors, n
if (n < 0) {
console.log('Failed to set the vertex information');
return;
}
// Get the storage location of u_MvpMatrix
var u_MvpMatrix = gl.getUniformLocation(gl.program, 'u_MvpMatrix');
if (!u_MvpMatrix) {
console.log('Failed to get the storage location of u_MvpMatrix');
return;
}
// Register the event handler
initEventHandlers(canvas);
// Specify the color for clearing <canvas>
gl.clearColor(0, 0, 0, 1);
gl.enable(gl.DEPTH_TEST);
// Start drawing
var tick = function () {
//setPerspective()宽高比
var aspect = canvas.width / canvas.height;
//
draw(gl, n, aspect, cx, cy, cz, fovy, u_MvpMatrix);
requestAnimationFrame(tick, canvas);
};
tick();
}
3) 使用缓冲区对象
在函数initVertexBuffers（）中包含了使用缓冲区对象向顶点着色器传入多个顶点数据的过程：
- 创建缓冲区对象（gl.createBuffer()）；
- 绑定缓冲区对象（gl.bindBuffer()）；
- 将数据写入缓冲区对象（gl.bufferData）；
- 将缓冲区对象分配给一个attribute变量(gl.vertexAttribPointer)
- 开启attribute变量(gl.enableVertexAttribArray)；
在本例中，在JS中申请的数组verticesColors分成位置和颜色两部分分配给缓冲区对象，并传入顶点着色器；vertexAttribPointer（）是其关键的函数，需要详细了解其参数的用法。最后，把顶点数据的索引值绑定到缓冲区对象，WebGL可以访问索引来间接访问顶点数据进行绘制。
function initVertexBuffers(gl, verticesColors) {
//DEM的一个网格是由两个三角形组成的
// 0------1 1
// | |
// | |
// col col------col+1
var indices = new Uint16Array((row - 1) * (col - 1) * 6);
var ci = 0;
for (var yi = 0; yi < row - 1; yi++) {
for (var xi = 0; xi < col - 1; xi++) {
indices[ci * 6] = yi * col + xi;
indices[ci * 6 + 1] = (yi + 1) * col + xi;
indices[ci * 6 + 2] = yi * col + xi + 1;
indices[ci * 6 + 3] = (yi + 1) * col + xi;
indices[ci * 6 + 4] = (yi + 1) * col + xi + 1;
indices[ci * 6 + 5] = yi * col + xi + 1;
ci++;
}
}
//创建缓冲区对象
var vertexColorBuffer = gl.createBuffer();
var indexBuffer = gl.createBuffer();
if (!vertexColorBuffer || !indexBuffer) {
return -1;
}
// 将缓冲区对象绑定到目标
gl.bindBuffer(gl.ARRAY_BUFFER, vertexColorBuffer);
// 向缓冲区对象中写入数据
gl.bufferData(gl.ARRAY_BUFFER, verticesColors, gl.STATIC_DRAW);
//
var FSIZE = verticesColors.BYTES_PER_ELEMENT;
// 向缓冲区对象分配a_Position变量
var a_Position = gl.getAttribLocation(gl.program, 'a_Position');
if (a_Position < 0) {
console.log('Failed to get the storage location of a_Position');
return -1;
}
gl.vertexAttribPointer(a_Position, 3, gl.FLOAT, false, FSIZE * 6, 0);
//开启a_Position变量
gl.enableVertexAttribArray(a_Position);
// 向缓冲区对象分配a_Color变量
var a_Color = gl.getAttribLocation(gl.program, 'a_Color');
if (a_Color < 0) {
console.log('Failed to get the storage location of a_Color');
return -1;
}
gl.vertexAttribPointer(a_Color, 3, gl.FLOAT, false, FSIZE * 6, FSIZE * 3);
//开启a_Color变量
gl.enableVertexAttribArray(a_Color);
// 写入并绑定顶点数组的索引值
gl.bindBuffer(gl.ELEMENT_ARRAY_BUFFER, indexBuffer);
gl.bufferData(gl.ELEMENT_ARRAY_BUFFER, indices, gl.STATIC_DRAW);
return indices.length;
}
4. 其他
1.这里用到了几个《WebGL编程指南》书中提供的JS组件。全部源代码(包含DEM数据)地址链接：https://share.weiyun.com/5cvt8PJ ，密码：4aqs8e。
2.如果关心如何设置模型视图投影变换矩阵，以及绑定鼠标键盘事件，可参看这篇文章：WebGL或OpenGL关于模型视图投影变换的设置技巧。
3.渲染的结果如果加入光照，效果会更好。
WebGL的颜色渲染-渲染一张DEM(数字高程模型)的更多相关文章
- WebGIS 利用 WebGL 在 MapboxGL 上渲染 DEM 三维空间数据
毕业两年,一直在地图相关的公司工作,虽然不是 GIS 出身,但是也对地图有些耳濡目染:最近在看 WebGl 的东西,就拿 MapboxGL 做了一个关于 WebGL 的三维数据渲染的 DEMO 练手. ...
- WebGL 纹理颜色原理
本文由云+社区发表 作者:ivweb qcyhust 导语 WebGL绘制图像时,往着色器中传入颜色信息就可以给图形绘制出相应的颜色,现在已经知道顶点着色器和片段着色器一起决定着向颜色缓冲区写入颜色信 ...
- [WebGL入门]四，渲染准备
注意:文章翻译http://wgld.org/,原作者杉本雅広(doxas),文章中假设有我的额外说明,我会加上［lufy:］.另外.鄙人webgl研究还不够深入,一些专业词语,假设翻译有误,欢迎大家 ...
- 基于WebGL的三维地形渲染
1.生成WebMap页面 #!/usr/bin/env python # -*- coding: utf-8 -*- import subprocess from jinja2 import Envi ...
- vue打包空白，图片没加载，背景颜色没有渲染出来-配置秘诀
找到config文件夹下的index.js文件修改一下位置 看清楚是 build(上边还有个dev 是开发环境下的配置,不需要改动)下的 assetsPublicPath :将‘/’改为‘./’ 在c ...
- 用ARCGIS配出一张DEM专题图
专题图是指突出而尽可能完善.详尽地表达制图区内的一种或几种自然或社会经济要素的地图.专题图的制图领域宽广,凡具有空间属性的信息数据都可以用其来表示.由于DEM描述的是地面高程信息,它在测绘.水文.气象 ...
- layui下拉框数据过万渲染渲染问题解决方案
方案一:layui下拉框分页插件 https://fly.layui.com/jie/29002/ 此插件我用了下浏览器缓存有问题,而且当下拉框数据量过万后,会一直渲染不出来,期待后期作者优化 如图下 ...
- 赛门铁克通配符SSL证书，一张通配型证书实现全站加密
赛门铁克通配型SSL证书,验证域名所有权和企业信息,属于企业验证(OV) 级SSL证书,最高支持256位加密.申请通配符SSL证书可以保护相同主域名下无限数量的多个子域名(主机).例如,一个通配符 ...
- 全球数字高程数据（DEM）详解，还有地形晕渲、等高线等干货
1 基本概念 DEM是数字高程模型的英文简称(Digital Elevation Model),是研究分析地形.流域.地物识别的重要原始资料.由于DEM 数据能够反映一定分辨率的局部地形特征,因此通过 ...
随机推荐
- 图片上传是否为空，以及类型的js验证
function check2() { var file = document.getElementsByName("file").value; if(file=="&q ...
- jaspersoft中分组打印
一:前言 使用IReport已经四个月了,最近在做一个保镖,是要按照类型分类,并且这些类型要横着打印,最后还要算这个类型金额的总值,这张报表现是说需要用到子报表,最后和一个同事一起用group来分组做 ...
- GML3示例
GML3示例:https://svn.osgeo.org/geotools/trunk/modules/extension/xsd/xsd-gml3/src/test/resources/org/ge ...
- time,random，os,sys,序列化模块
一.time模块 表示时间的三种方式 在Python中,通常有这三种方式来表示时间:时间戳.元组(struct_time).格式化的时间字符串: (1)时间戳(timestamp) :通常来说,时间戳 ...
- 详解SHOW PROCESSLIST显示哪些线程正在运行列出的状态
SHOW PROCESSLIST显示哪些线程正在运行.您也可以使用mysqladmin processlist语句得到此信息.如果您有SUPER权限,您可以看到所有线程.否则,您只能看到您自己的线程( ...
- python 监控redis的进程与端口
#!/usr/bin/python # -*- coding:utf-8 -*- import glob,psutil import json,os,datetime import collectio ...
- Android之进程通信--Binder
Cilent从ServiceManger哪里获得BnMediaService的BnBinder引用就可以调用BnMediaPlayerService的方法了,BnMediaPlayerService是 ...
- bindingSource具体使用案例
界面如下: using DevExpress.XtraBars.Docking; using DevExpress.XtraEditors; using NewPwrDY.DBEntity; usin ...
- 【C++】嵌套类、友元
黄邦勇帅 里面关于嵌套类的介绍我有疑惑.里面11.9说在创建一个外围类的对象时先执行嵌套类的构造函数然后再执行外围类的构造函数,析构函数则以相反的方式执行. 可是我编程实验了一下,创建外围类对象时并不 ...
- Selenium2+python自动化60-异常后截图（screenshot）【转载】
前言 在执行用例过程中由于是无人值守的,用例运行报错的时候,我们希望能对当前屏幕截图,留下证据. 在写用例的时候,最后一步是断言,可以把截图的动作放在断言这里,那么如何在断言失败后截图呢? 一.截图方 ...