【Cesium】在着色器中计算片元的高程

问题提出

有时候在 Cesium 上做一些效果时,需要知道片元的高程。而渲染是基于笛卡尔坐标的,没有高程这个概念。怎么在着色器中计算片元的高程呢?

基本思路

通常做法是在着色器里对深度纹理采样,将片元的深度值转换为相机坐标。然后把相机坐标转换为世界坐标,世界坐标是片元相对于地球椭球中心的偏移量,下图中假设片元的世界坐标为 (x, y, z)。用世界坐标的模减去片元对应的椭球半径就可以得到高程。但是逐片元去计算它们对应的椭球半径在着色器里不好实现,在CPU中计算相机位置处的椭球半径(下图中的 earthRadius )来近似地作为相机观察范围内所有片元对应的椭球半径从性能和实现难度上都是更优解。

图1 计算片元高程基本思路

具体步骤

1、在CPU中计算相机所在位置对应的地球半径,作为 uniform 变量传入着色器中:相机世界坐标的模 - 相机的高程

javascript 复制代码
const earthRadius = Cartesian3.magnitude(camera.positionWC) - camera.positionCartographic.height

2、在着色器中做如下转换:深度值 -> 相机坐标 -> 世界坐标,最后用片元世界坐标的模 - 相机位置对应的地球半径即可,glsl 代码如表1所示。虽然是用相机位置处的椭球半径近似替代了片元处的地球半径,但是误差很难看得出来,这种方法适用于绝大部分应用场景。

cpp 复制代码
float altitude = 0.0;
float depth = czm_unpackDepth(texture2D(depthTexture, texCoords));
if (depth == 0.0) {
	altitude = czm_infinity;
}
vec4 eyeCoordinate = czm_windowToEyeCoordinates(gl_FragCoord.xy, depth);
float distance = -eyeCoordinate.z / eyeCoordinate.w;
vec4 worldCoordinate4 = czm_inverseView * eyeCoordinate;
vec3 worldCoordinate = worldCoordinate4.xyz / worldCoordinate4.w;
altitude = length(worldCoordinate.xyz) - earthRadius;

表1 glsl 中计算片元的高程

相关推荐
Struggler2816 分钟前
Chrome插件开发
前端
前端 贾公子18 分钟前
Monorepo + vite 怎么热更新
前端
然我1 小时前
不用 Redux 也能全局状态管理?看我用 useReducer+Context 搞个 Todo 应用
前端·javascript·react.js
前端小巷子1 小时前
Web 实时通信:从短轮询到 WebSocket
前端·javascript·面试
神仙别闹1 小时前
基于C#+SQL Server实现(Web)学生选课管理系统
前端·数据库·c#
web前端神器1 小时前
指定阿里镜像原理
前端
枷锁—sha1 小时前
【DVWA系列】——CSRF——Medium详细教程
android·服务器·前端·web安全·网络安全·csrf
枷锁—sha1 小时前
跨站请求伪造漏洞(CSRF)详解
运维·服务器·前端·web安全·网络安全·csrf
群联云防护小杜2 小时前
深度隐匿源IP:高防+群联AI云防护防绕过实战
运维·服务器·前端·网络·人工智能·网络协议·tcp/ip
汉得数字平台2 小时前
【鲲苍提效】全面洞察用户体验,助力打造高性能前端应用
前端·前端监控