在学习《UnityShader入门精要》的一个笔记
一.13.2节中,缺乏对深度纹理采样后的深度值如何得到世界空间下的坐标表示的解释
对应代码如下
//深度纹理采样,本质是返回纹理采样结果的r通道
float d = SAMPLE_DEPTH_TEXTURE(_CameraDepthTexture, i.uv_depth);
//获取投影空间下的坐标(数值范围在-1~1)
float4 H = float4(i.uv.x * 2 - 1, i.uv.y * 2 - 1, d * 2 - 1, 1);
// 通过视角*投影矩阵的逆矩阵,得到视角空间下的坐标
float4 D = mul(_CurrentViewProjectionInverseMatrix, H);
// 除以w分量,将其转为世界坐标
float4 worldPos = D / D.w;
令世界坐标系中有一点P,齐次坐标记为
‘
v
w
o
r
l
d
=
(
x
,
y
,
z
,
w
)
T
‘
` v_{world}=(x,y,z,w)^T `
‘vworld=(x,y,z,w)T‘,其中
‘
v
w
o
r
l
d
.
w
=
1.0
‘
`v_{world}.w=1.0`
‘vworld.w=1.0‘, 则有
‘
v
c
l
i
p
=
M
p
r
o
j
∗
M
v
i
e
w
∗
v
w
o
r
l
d
‘
` v_{clip} = M_{proj}*M_{view}*v_{world}`
‘vclip=Mproj∗Mview∗vworld‘ (1)
其中,
‘
v
c
l
i
p
‘
` v_{clip} `
‘vclip‘表示P在裁剪空间中的坐标,
‘
M
v
i
e
w
‘
`M_{view} `
‘Mview‘表示视角矩阵,
‘
M
p
r
o
j
‘
`M_{proj} `
‘Mproj‘表示投影矩阵,将上述两个矩阵相乘得到当前的观察投影矩阵
‘
M
v
p
‘
`M_{vp} `
‘Mvp‘,则有
‘
v
w
o
r
l
d
=
M
v
p
−
1
∗
v
c
l
i
p
‘
` v_{world} = M_{vp}^{-1}*v_{clip}`
‘vworld=Mvp−1∗vclip‘ (2)
P点的深度纹理采样结果为
‘
d
‘
`d`
‘d‘,P点的NDC坐标记为
‘
v
n
d
c
‘
` v_{ndc} `
‘vndc‘,根据Unity的深度纹理的存储规则,d是NDC映射而来,使用原映射的反函数即可,则有
‘
v
n
d
c
=
d
∗
2
−
1
‘
` v_{ndc} = d*2-1`
‘vndc=d∗2−1‘ (3)
透视除法的过程是裁剪空间坐标的4个分量分别除以w分量,即
‘
v
n
d
c
=
v
c
l
i
p
/
v
c
l
i
p
.
w
‘
` v_{ndc} = v_{clip}/v_{clip}.w`
‘vndc=vclip/vclip.w‘ (4)
由公式(2)和公式(3)得
‘
v
w
o
r
l
d
=
M
v
p
−
1
∗
v
n
d
c
∗
v
c
l
i
p
.
w
‘
` v_{world} = M_{vp}^{-1}*v_{ndc}*v_{clip}.w`
‘vworld=Mvp−1∗vndc∗vclip.w‘ (5)
在公式(5)中,
‘
M
v
p
−
1
‘
`M_{vp}^{-1} `
‘Mvp−1‘由_CurrentViewProjectionInverseMatrix
提供,
‘
v
n
d
c
‘
` v_{ndc} `
‘vndc‘由
‘
v
n
d
c
=
d
∗
2
−
1
‘
` v_{ndc} = d*2-1`
‘vndc=d∗2−1‘得到,因为
‘
v
w
o
r
l
d
.
w
=
1.0
‘
`v_{world}.w=1.0`
‘vworld.w=1.0‘,则有
‘
v
w
o
r
l
d
.
w
=
(
M
v
p
−
1
∗
v
n
d
c
)
.
w
∗
v
c
l
i
p
.
w
=
1.0
‘
` v_{world}.w = (M_{vp}^{-1}*v_{ndc}).w*v_{clip}.w = 1.0`
‘vworld.w=(Mvp−1∗vndc).w∗vclip.w=1.0‘ (6)
由公式(5)可得
‘
v
c
l
i
p
.
w
=
1
,
0
(
M
v
p
−
1
∗
v
n
d
c
)
.
w
‘
` v_{clip}.w = \frac{1,0}{(M_{vp}^{-1}*v_{ndc}).w} `
‘vclip.w=(Mvp−1∗vndc).w1,0‘ (7)
将公式(7)带入公式(5),则有
‘
v
w
o
r
l
d
=
M
v
p
−
1
∗
v
n
d
c
(
M
v
p
−
1
∗
v
n
d
c
)
.
w
‘
` v_{world} = \frac{M_{vp}^{-1}*v_{ndc}}{(M_{vp}^{-1}*v_{ndc}).w}`
‘vworld=(Mvp−1∗vndc).wMvp−1∗vndc‘
且由公式(3),则可知上文H为
‘
v
n
d
c
‘
` v_{ndc} `
‘vndc‘,D为
‘
(
M
v
p
−
1
∗
v
n
d
c
)
‘
`{(M_{vp}^{-1}*v_{ndc})} `
‘(Mvp−1∗vndc)‘,则世界坐标可由D / D.w得到
参考(侵删):osg中NDC坐标转换为世界坐标
Github:TipsDuringReadingUnityShadersBook