MarchineCubes实现思路总结
MarchineCubes定义
是一种基于体素构建三维模型的方式，有些类似《我的世界》中的堆方块，但实际上，建模是以方块之间的交点为中心点，每个cube表示八个象限的相交模型
实现思路
在三维空间中划分网格，每个网格是一个cube，其8个顶点各有两个状态{in、out}，分别表示该顶点是否位于三维模型的内部
根据这8个顶点的状态，可确定当前cube自身的模型，这个排列组合正好是一个8位二进制
这里需要对这8个顶点进行约定，为后续代码提供方便
这样约定的好处是，某些临边、对角位置可用二进制位运算方便获取，临边是相应位异或、对角是按位取反
Cube建模
总共有255种组合，其中有些是重复形状的翻转、和旋转，且0x00、0xff实际上是空模型，那么一个cube就需要254个模型，我没有采用Polygonising中的三角形列表，而是用python在maya中建模
import maya.cmds as cmd
for i in range(1,255):names=[]for j in range(0,8):if 1<1:cmd.polyCBoolOp(names,n="aaa%d"%(i))elif len(names)>0:cmd.rename(names[-1],"aaa%d"%(i))if len(names)>0:cmd.makeIdentity("aaa%d"%(i),a=1,t=1)
这样就建立好1~254所有模型，并重叠摆在原点中心。然后清空历史，选中所有模型，手动将模型上、下、左、右、前、后0.5米之外的部分切割掉
| 切割前 | 切割后 |
|---|---|
 |  |
切割之后，就是一个1x1单位中心对齐的模型，再通过脚本将其一字排开
for i in range(1,255):cmd.move(i*1.5,0,0,"aaa%d"%(i))
在Outliner中可以看到每个mesh被命名为aaa1 ~ aaa254，方便后续引擎中按索引直接取出对应模型
数据结构及伪代码
这里采用自研引擎+lua环境实现
由于lua数组实际是hash结构，所以可以将网格中的三维坐标转为一个整形数值用来做hash key，也就没必要存储为三维数组了
-- id表示为三维空间坐标，由于double有效数为9007199254740992，所以用1024^3足够，哈希值：x*1024^2+y*1024+z
function coordToId(x,y,z)return x*1024^2+y*1024+z
end-- 那么当某个hash反推：x=math.floor(v/1024/1024)、y=math.floor(v/1024%1024)、z=v%1024
function idToCoord(id)return math.floor(id/1024/1024),math.floor(id/1024%1024),id%1024
end
定义grids、cubes数据
function EditorBehavior:__init(name)...self.grids={}self.cubes={}
end
grids表示网格的交点坐标
cubes表示围绕交点周围的模型，也就是说cubes在每一维度都比grids多一个单位
当修改某个grids的状态，in：1、out：nil
function EditorBehavior:AddGrid(x,y,z)self.grids[coordToId(x,y,z)]=1self:UpdateTileRoundTheGrid(x,y,z,true)
endfunction EditorBehavior:RemoveGrid(x,y,z)self.grids[coordToId(x,y,z)]=nilself:UpdateTileRoundTheGrid(x,y,z,false)
end
然后根据这个grids内容更新周围八个象限cubes的状态值，若某个cube更新后数值为0，则清除这个cube
function EditorBehavior:UpdateTileRoundTheGrid(x,y,z,inside)for i=0,1 do for j=0,1 do for k=0,1 dolocal pos=i*4+j*2+klocal ipos=bit.band(bit.bnot(pos),0x7)local mask=bit.lshift(0x1,ipos)local cube=self.cubes[coordToId(x+i,y+j,z+k)]if cube thenif inside thenself.cubes[coordToId(x+i,y+j,z+k)]=bit.bor(cube,mask)elsecube=bit.band(cube,bit.bnot(mask))if 0==cube thenself.cubes[coordToId(x+i,y+j,z+k)]=nilelseself.cubes[coordToId(x+i,y+j,z+k)]=cubeendendelseif inside thenself.cubes[coordToId(x+i,y+j,z+k)]=maskendendend end end
end
在渲染时，将所有cube的状态值，以及模型渲染出来
这里cube的状态值，正好对应之前maya中模型的aaa1 ~ aaa254
参考
https://www.cs.carleton.edu/cs_comps/0405/shape/marching_cubes.html
http://paulbourke.net/geometry/polygonise/
https://www.researchgate.net/publication/202232897_Marching_Cubes_A_High_Resolution_3D_Surface_Construction_Algorithm
https://www.bilibili.com/video/BV1yJ411r73v