skia库的3D坐标变换其实也是无奈之举。参照SKCamera.cpp：

首先，定义了一个虚拟相机：

void SkCamera3D::reset() {
fLocation.set(0, 0, -SkIntToScalar(576));   // 8 inches backward
fAxis.set(0, 0, SK_Scalar1);                // forward
fZenith.set(0, -SK_Scalar1, 0);             // up
fObserver.set(0, 0, fLocation.fZ);
fNeedToUpdate = true;
}

求出相机的投影矩阵。这个代码里面叫做方位（orientation）

这个方位矩阵完成了从3D空间到2D空间的投影，

这个方法研究了很久，具体分析见注释，不多讲了。

相当于OpenGL的glFrustum函数。代码如下：

void SkCamera3D::doUpdate() const {
SkUnit3D    axis, zenith, cross;
fAxis.normalize(&axis);
{
SkScalar dot = SkUnit3D::Dot(*(const SkUnit3D*)(const void*)&fZenith, axis);
zenith.fX = fZenith.fX - SkUnitScalarMul(dot, axis.fX);
zenith.fY = fZenith.fY - SkUnitScalarMul(dot, axis.fY);
zenith.fZ = fZenith.fZ - SkUnitScalarMul(dot, axis.fZ);
(void)((SkPoint3D*)(void*)&zenith)->normalize(&zenith);
}
/*
[-z, 0, x]   [cross.fX,     cross.fY,       cross.fZ]
[0, -z, y]* [zenith.fx,     zenith.fY,  zenith.fZ]
[0,  0,  1]  [axis.fx,  axis.fY,        axis.fZ]
下面这个矩阵式是投影矩阵:
[-z, 0, x]
[0, -z, y]
[0,  0,  1]
投影到x,y z位置的投影矩阵，其中x, y恒为0。如果不为0，经推算，
这个矩阵是有错误的，用三角形相似，可以算出，应该是下面这个样子
[-z, 0, x-zx]
[0, -z, y-zy]
[0,  0,  1]
*/
SkUnit3D::Cross(axis, zenith, &cross);
{
SkMatrix* orien = &fOrientation;
SkScalar x = fObserver.fX;
SkScalar y = fObserver.fY;
SkScalar z = fObserver.fZ;
orien->set(SkMatrix::kMScaleX, SkUnitScalarMul(x, axis.fX) - SkUnitScalarMul(z, cross.fX));
orien->set(SkMatrix::kMSkewX,  SkUnitScalarMul(x, axis.fY) - SkUnitScalarMul(z, cross.fY));
orien->set(SkMatrix::kMTransX, SkUnitScalarMul(x, axis.fZ) - SkUnitScalarMul(z, cross.fZ));
orien->set(SkMatrix::kMSkewY,  SkUnitScalarMul(y, axis.fX) - SkUnitScalarMul(z, zenith.fX));
orien->set(SkMatrix::kMScaleY, SkUnitScalarMul(y, axis.fY) - SkUnitScalarMul(z, zenith.fY));
orien->set(SkMatrix::kMTransY, SkUnitScalarMul(y, axis.fZ) - SkUnitScalarMul(z, zenith.fZ));
orien->set(SkMatrix::kMPersp0, axis.fX);
orien->set(SkMatrix::kMPersp1, axis.fY);
orien->set(SkMatrix::kMPersp2, axis.fZ);
}
}

然后在patchToMatrix方法中，使用此投影矩阵对当前的变换UVO坐标系进行了变换，代码如下：

void SkCamera3D::patchToMatrix(const SkPatch3D& quilt, SkMatrix* matrix) const {
if (fNeedToUpdate) {
this->doUpdate();
fNeedToUpdate = false;
}
const SkScalar* mapPtr = (const SkScalar*)(const void*)&fOrientation;
const SkScalar* patchPtr;
SkPoint3D       diff;
SkScalar        dot;
diff.fX = quilt.fOrigin.fX - fLocation.fX;
diff.fY = quilt.fOrigin.fY - fLocation.fY;
diff.fZ = quilt.fOrigin.fZ - fLocation.fZ;
dot = SkUnit3D::Dot(*(const SkUnit3D*)(const void*)&diff,
*(const SkUnit3D*)(((const SkScalar*)(const void*)&fOrientation) + 6));
//  patchPtr的结构为{U,V,ORIGIN}其中U，V 代表列向量
// ORIGIN 是坐标原点
patchPtr = (const SkScalar*)&quilt;
/*
其中matrix表示一个3x3的矩阵
第一行代表U的系数，第二行是V的系数，第三行是diff 的系数
matrix 的每一列代表的是一个坐标轴。
*/
//第一列
matrix->set(SkMatrix::kMScaleX, SkScalarDotDiv(3, patchPtr, 1, mapPtr, 1, dot));
matrix->set(SkMatrix::kMSkewY,  SkScalarDotDiv(3, patchPtr, 1, mapPtr+3, 1, dot));
matrix->set(SkMatrix::kMPersp0, SkScalarDotDiv(3, patchPtr, 1, mapPtr+6, 1, dot));
//第二列
patchPtr += 3;
matrix->set(SkMatrix::kMSkewX,  SkScalarDotDiv(3, patchPtr, 1, mapPtr, 1, dot));
matrix->set(SkMatrix::kMScaleY, SkScalarDotDiv(3, patchPtr, 1, mapPtr+3, 1, dot));
matrix->set(SkMatrix::kMPersp1, SkScalarDotDiv(3, patchPtr, 1, mapPtr+6, 1, dot));
//第三列
patchPtr = (const SkScalar*)(const void*)&diff;
matrix->set(SkMatrix::kMTransX, SkScalarDotDiv(3, patchPtr, 1, mapPtr, 1, dot));
matrix->set(SkMatrix::kMTransY, SkScalarDotDiv(3, patchPtr, 1, mapPtr+3, 1, dot));
matrix->set(SkMatrix::kMPersp2, SK_UnitScalar1);
}

到此为止，也就把这个变换的思路给理解清楚了。跟OpenGL的思路类似。

第一。进行模型视图变换，通过一个patch的东西，UVO坐标系。

第二。进行投影变换。定位相机，具体参照如下函数：

 /*
由于x，y，z 与cross的点积，就是在cross方向上的投影。因此，对于x，y，z向量来说。
从世界坐标转换到相机坐标系的矩阵就是
[cross.fX,  cross.fY,       cross.fZ]
[zenith.fx,     zenith.fY,  zenith.fZ]
[axis.fx,     axis.fY,        axis.fZ]
由于是标准正交基，
他的转置或者说，逆矩阵，就是相机矩阵
[cross.fX,     zenith.fx,     axis.fx]
[cross.fY,   zenith.fY,  axis.fY]
[ cross.fZ ,zenith.fZ, axis.fZ]
[-z, 0, x]   [cross.fX,     cross.fY,       cross.fZ]
[0, -z, y]* [zenith.fx,     zenith.fY,  zenith.fZ]
[0,  0,  1]  [axis.fx,  axis.fY,        axis.fZ]
下面这个矩阵式是投影矩阵:
[-z, 0, x]
[0, -z, y]
[0,  0,  1]
*/

增加一些注释

下面增加一些系统介绍

SKCamera3D中关于模拟3D的算法研究

SkCamera是一个可以支持3D变换的组件，通过一些矩阵变换，最终获得了3D的效果。与OpenGL的思路基本相同，分为两个矩阵，一个是投影矩阵，一个是模型视图矩阵。

下面就这两个矩阵，并结合SkCamera的代码进行讲解。

1.       投影矩阵。

投影矩阵是通过SkCamera3D这个类进行生成和管理。SkCamera3D这个的作用还在于将一个SkPatch3D的一个坐标系转换为投影后的坐标系，通过patchToMatrix完成。

总而言之，这个类负责投影的相关工作。

class SkCamera3D {

public:

SkCamera3D();

void reset();

void update();

void patchToMatrix(const SkPatch3D&, SkMatrix* matrix) const;

SkPoint3D   fLocation;

SkPoint3D   fAxis;

SkPoint3D   fZenith;

SkPoint3D   fObserver;

private:

mutable SkMatrix    fOrientation;

mutable bool        fNeedToUpdate;

void doUpdate() const;

};

2.       基本投影算法

世界坐标系如图所示

X

Z

Y

O

相机方向

相机顶的方向

然后在这个世界坐标系中，根据想要的投影方式，初始化了一个相机（类似于OpenGL中的glFrustum函数）。相机顶朝下，相机的摄像头朝着屏幕里面。

对应SkCamera.cpp的代码如下：

void SkCamera3D::reset() {

fLocation.set(0, 0, -SkIntToScalar(576));   // 8 inches backward

fAxis.set(0, 0, SK_Scalar1);                // forward

fZenith.set(0, -SK_Scalar1, 0);             // up

fObserver.set(0, 0, fLocation.fZ);

fNeedToUpdate = true;

}

然后根据相机的位置获得一个投影矩阵。

公式如下：

【相机投影矩阵】*【相机世界矩阵】

根据2维齐次坐标的几何意义，也就是在Z=1平面上的投影。现在变为在Z=-z平面上的投影。具体参照如下示意图：这个图也展示了二维齐次坐标系的几何意义（投影+映射）：

Z=0

P

Z=1

Z=-z

-z， 0， x

0， -z， y

0，    0， 1

具体可以验证这个，假设x’,y’,z’

那么，可得：

X’’ = -x’z + xz’

Y’’ = -y’z + yz’

Z’’ =  z’

进而得到了相机的坐标

X’’ = -x’z/z’ + x

Y’’ = -y’z/z’ + y

Z’’ =  1

另外，从三角形相似原理可以知道：

X’’ = -x’z/z’ 和Y’’ = -y’z/z’刚好把坐标刚好映射到了z平面上。然后进行平移，再加上相机的位置x，y 从而得到如下公式：

X’’ = -x’z/z’ + x

Y’’ = -y’z/z’ + y

这两个公式是相同的，从而验证了映射的正确性。

下面就是投影矩阵的具体实现方法：

void SkCamera3D::doUpdate() const {

SkUnit3D    axis, zenith, cross;

fAxis.normalize(&axis);

{

SkScalar dot = SkUnit3D::Dot(*(const SkUnit3D*)(const void*)&fZenith, axis);

zenith.fX = fZenith.fX - SkUnitScalarMul(dot, axis.fX);

zenith.fY = fZenith.fY - SkUnitScalarMul(dot, axis.fY);

zenith.fZ = fZenith.fZ - SkUnitScalarMul(dot, axis.fZ);

(void)((SkPoint3D*)(void*)&zenith)->normalize(&zenith);

}

SkUnit3D::Cross(axis, zenith, &cross);

{

SkMatrix* orien = &fOrientation;

SkScalar x = fObserver.fX;

SkScalar y = fObserver.fY;

SkScalar z = fObserver.fZ;

orien->set(SkMatrix::kMScaleX, SkUnitScalarMul(x, axis.fX) - SkUnitScalarMul(z, cross.fX));

orien->set(SkMatrix::kMSkewX,  SkUnitScalarMul(x, axis.fY) - SkUnitScalarMul(z, cross.fY));

orien->set(SkMatrix::kMTransX, SkUnitScalarMul(x, axis.fZ) - SkUnitScalarMul(z, cross.fZ));

orien->set(SkMatrix::kMSkewY,  SkUnitScalarMul(y, axis.fX) - SkUnitScalarMul(z, zenith.fX));

orien->set(SkMatrix::kMScaleY, SkUnitScalarMul(y, axis.fY) - SkUnitScalarMul(z, zenith.fY));

orien->set(SkMatrix::kMTransY, SkUnitScalarMul(y, axis.fZ) - SkUnitScalarMul(z, zenith.fZ));

orien->set(SkMatrix::kMPersp0, axis.fX);

orien->set(SkMatrix::kMPersp1, axis.fY);

orien->set(SkMatrix::kMPersp2, axis.fZ);

}

}

3.       模型变换矩阵的生成。

模型变换是通过如下结构体完成的。

struct SkMatrix3D {

SkScalar    fMat[3][4];

void reset();

void setRow(int row, SkScalar a, SkScalar b, SkScalar c, SkScalar d = 0)

{

SkASSERT((unsigned)row < 3);

fMat[row][0] = a;

fMat[row][1] = b;

fMat[row][2] = c;

fMat[row][3] = d;

}

void setRotateX(SkScalar deg);

void setRotateY(SkScalar deg);

void setRotateZ(SkScalar deg);

void setTranslate(SkScalar x, SkScalar y, SkScalar z);

void preRotateX(SkScalar deg);

void preRotateY(SkScalar deg);

void preRotateZ(SkScalar deg);

void preTranslate(SkScalar x, SkScalar y, SkScalar z);

void setConcat(const SkMatrix3D& a, const SkMatrix3D& b);

void mapPoint(const SkPoint3D& src, SkPoint3D* dst) const;

void mapVector(const SkVector3D& src, SkVector3D* dst) const;

void mapPoint(SkPoint3D* v) const

{

this->mapPoint(*v, v);

}

void mapVector(SkVector3D* v) const

{

this->mapVector(*v, v);

}

};

这样得到的矩阵式3D世界的3*4的矩阵。而投影矩阵只能对3*3的矩阵进行变换。而patch构成的坐标系却是3*3的。无法相乘。那么程序如何实现patch从3*4 到 3*3的一个转换呢。

下面是推导过程：

Ux    Vx     Nx    Dx

Uy    Vy     Ny    Dy

Uz    Vz     Nz    Dz

0       0       0       1

这个就是patch坐标系对应的3D世界的变换矩阵。其中U，V，N为基向量。D为坐标原点。

我们称之为D-UVN坐标系。

最重要的信息是，在D-UVN坐标系中，N轴上的值横为0.因为一开始，一个坐标x，y一定是平面的。因此，只需要求出对【x，y， 0， 1】的一个变换就行了。

可以得出：

x‘= Ux*x + Vx*y + Dx

y’ =Uy*x + Vy*y + Dy

z’=Uz*x + Vz*y + Dz

也就是如下矩阵：

Ux    Vx     Dx

Uy    Vy     Dy

Uz    Vz     Dz

代码如下：

void SkCamera3D::patchToMatrix(const SkPatch3D& quilt, SkMatrix* matrix) const {

    if (fNeedToUpdate) {

        this->doUpdate();

        fNeedToUpdate = false;

    }

 

    const SkScalar* mapPtr = (const SkScalar*)(const void*)&fOrientation;

    const SkScalar* patchPtr;

    SkPoint3D       diff;

    SkScalar        dot;

 

    diff.fX = quilt.fOrigin.fX - fLocation.fX;

    diff.fY = quilt.fOrigin.fY - fLocation.fY;

    diff.fZ = quilt.fOrigin.fZ - fLocation.fZ;

 

    dot = SkUnit3D::Dot(*(const SkUnit3D*)(const void*)&diff,

                        *(const SkUnit3D*)(((const SkScalar*)(const void*)&fOrientation) + 6));

 

    patchPtr = (const SkScalar*)&quilt;

    matrix->set(SkMatrix::kMScaleX, SkScalarDotDiv(3, patchPtr, 1, mapPtr, 1, dot));

    matrix->set(SkMatrix::kMSkewY,  SkScalarDotDiv(3, patchPtr, 1, mapPtr+3, 1, dot));

    matrix->set(SkMatrix::kMPersp0, SkScalarDotDiv(3, patchPtr, 1, mapPtr+6, 1, dot));

 

    patchPtr += 3;

    matrix->set(SkMatrix::kMSkewX,  SkScalarDotDiv(3, patchPtr, 1, mapPtr, 1, dot));

    matrix->set(SkMatrix::kMScaleY, SkScalarDotDiv(3, patchPtr, 1, mapPtr+3, 1, dot));

    matrix->set(SkMatrix::kMPersp1, SkScalarDotDiv(3, patchPtr, 1, mapPtr+6, 1, dot));

 

    patchPtr = (const SkScalar*)(const void*)&diff;

    matrix->set(SkMatrix::kMTransX, SkScalarDotDiv(3, patchPtr, 1, mapPtr, 1, dot));

    matrix->set(SkMatrix::kMTransY, SkScalarDotDiv(3, patchPtr, 1, mapPtr+3, 1, dot));

    matrix->set(SkMatrix::kMPersp2, SK_UnitScalar1);

}

 

 

 

 

 

总结：

   Sk3DCamera提供了一种方法，可以对图形进行仿3D变换，这对于想在2D空间进行3D处理的程序来说，是一个不错的选择。