NCF参数化建筑论坛

标题: 马赛克【附代码】 [打印本页]

作者: 夜神 时间: 2010-6-22 16:11
标题: 马赛克【附代码】
[attach]8853[/attach]
代码未经过完全测试，可能个别情况下会出错。
发现bug请指出，谢谢。
'========================================

Option Explicit
run
SUB run
makeMosiac Rhino.GetObject("请选择物体", 0),Rhino.IntegerBox("最短方向的方块数",6),Rhino.RealBox("方块缩小比例",0.9)
END SUB
'--------------------------------------------
' nDensity：最小boundingBox边长的分段数
' rShrink：方块缩小比例
'--------------------------------------------
FUNCTION makeMosiac(obj,nDensity,rShrink)
DIM objs : objs=array(obj)
DIM pBox, boundingBox : pBox=Rhino.BoundingBox(objs)
DIM nX,nY,nZ  '包围盒三边的分段数
DIM minSide  '最短边代码
DIM lSide(2)  '包围盒三边长
DIM objTemp
DIM lCube '马赛克边长
DIM iX,iY,iZ  '三个边的截线数
'计算最短边、马赛克边长、三个方向的马赛克数
lSide(0)=pBox(1)(0)-pBox(0)(0)
lSide(1)=pBox(2)(1)-pBox(1)(1)
lSide(2)=pBox(4)(2)-pBox(0)(2)
IF lSide(0)<=lSide(1) AND lSide(0)<=lSide(2) THEN
  minSide=0
  lCube=lSide(minSide)/nDensity
  iX=nDensity
  iY=int(lSide(1)/lCube+0.5) '四舍五入
  iZ=int(lSide(2)/lCube+0.5)
END IF
IF lSide(1)<=lSide(0) AND lSide(1)<=lSide(2) THEN
  minSide=1
  lCube=lSide(minSide)/nDensity
  iX=int(lSide(0)/lCube+0.5)
  iY=nDensity
  iZ=int(lSide(2)/lCube+0.5)
END IF
IF lSide(2)<=lSide(0) AND lSide(2)<=lSide(1) THEN
  minSide=2
  lCube=lSide(minSide)/nDensity
  iX=int(lSide(0)/lCube+0.5)
  iY=int(lSide(1)/lCube+0.5)
  iZ=nDensity
END IF
'定义马赛克空间，并初始化为零（该空间没有马赛克）
DIM i,j,k,m,s
DIM p() : REDIM p((iX+1)*(iY+1)*(iZ+1))
FOR i=0 TO iX
  FOR j=0 TO iY
FOR k=0 TO iZ
p(i*(iY+1)*(iZ+1)+j*(iZ+1)+k)=0
NEXT
  NEXT
NEXT
DIM iRoundUp '圆整后的坐标序号
DIM pCube(7)
'拆分物体
DIM nObjs : nObjs=UBound(objs)  '选择要马赛克的物体的数目
DIM nSfs  : nSfs=0 '打散后的NURBS曲面总数
DIM sfs()  : REDIM sfs(nSfs) '打散后的NURBS曲面
DIM sfsInEach,sfsInEachPlane,pVertice,pMeshFaces
DIM arrVerticeNum,nMeshFace '顶点序号;Mesh内面片数 '
FOR i=0 TO nObjs '对每一个所选的物体
  IF Rhino.IsSurface(objs(i))=vbTRUE AND Rhino.IsSurfaceTrimmed(objs(i))=vbFALSE THEN  '是单曲面，且不是裁剪曲面
sfs(UBound(sfs))=objs(i)
REDIM PRESERVE sfs(UBound(sfs)+1)
  ELSE  '是复合曲面或裁剪曲面
Rhino.SelectObject objs(i)
Rhino.Command "-Mesh PolygonDensity 12  Enter Enter "  '数字12表示离散后的mesh精度，越大越精细，计算量也越大
objTemp=Rhino.FirstObject(vbTRUE) 'Mesh面
pMeshFaces=Rhino.MeshFaces(objTemp, vbTRUE)  'Mesh内所有顶点的集合
arrVerticeNum=Rhino.MeshFaceVertices(objTemp)  'Mesh内各面角点序号集合
nMeshFace=UBound(arrVerticeNum)  '面片总数
Rhino.DeleteObject objTemp  '删除Mesh
FOR j=0 TO nMeshFace '为每一个Mesh面片建立曲面
objTemp=Rhino.AddSrfPt(Array(pMeshFaces(j*4+0),pMeshFaces(j*4+1),pMeshFaces(j*4+2),pMeshFaces(j*4+3)))
sfs(UBound(sfs))=objTemp
REDIM PRESERVE sfs(UBound(sfs)+1)
NEXT
  END IF
NEXT
Rhino.UnSelectAllObjects
REDIM PRESERVE sfs(UBound(sfs)-1)
nSfs=UBound(sfs)
DIM pNorm,pStart,pEnd,sectLine,pSect,nSectPoints
'第一轮：平行Z轴
dim pTemp() :redim pTemp(0)
FOR i=0 TO iX
  FOR j=0 TO iY
pStart=Array(pBox(0)(0)+i*lCube, pBox(0)(1)+j*lCube, pBox(0)(2))
pEnd=Array(pBox(0)(0)+i*lCube, pBox(0)(1)+j*lCube, pBox(4)(2))
sectLine=Rhino.AddLine(pStart,pEnd)
FOR m=0 TO nSfs
pSect=Rhino.CurveSurfaceIntersection(sectLine, sfs(m))
IF IsArray(pSect)=vbTRUE THEN  '有交叉点，(此句费解，去掉IF后pSect就不是数组了)
   nSectPoints=UBound(pSect,1) '交叉点个数
   FOR k=0 TO nSectPoints
   iRoundUp=int((pSect(k,1)(2)-pStart(2))/lCube+0.5) '把交叉点位置圆整到马赛克中心
   p(i*(iY+1)*(iZ+1)+j*(iZ+1)+iRoundUp)=1 '该位置的马赛克状态为1（有马赛克）
   pTemp(ubound(pTemp))=Rhino.AddPoint(pSect(k,1))
   redim preserve pTemp(ubound(pTemp)+1)
   NEXT
END IF
NEXT
Rhino.DeleteObject sectLine
  NEXT
NEXT
'第二轮：平行X轴
FOR k=0 TO iZ
  FOR j=0 TO iY
pStart=Array(pBox(0)(0), pBox(0)(1)+j*lCube, pBox(0)(2)+k*lCube)
pEnd=Array(pBox(2)(0), pBox(0)(1)+j*lCube, pBox(0)(2)+k*lCube)
sectLine=Rhino.AddLine(pStart,pEnd)
FOR m=0 TO nSfs
pSect=Rhino.CurveSurfaceIntersection(sectLine, sfs(m))
IF IsArray(pSect) THEN '有交叉点
   nSectPoints=UBound(pSect,1) '交叉点个数
   FOR i=0 TO nSectPoints
   iRoundUp=int((pSect(i,1)(0)-pStart(0))/lCube+0.5)
   p(iRoundUp*(iY+1)*(iZ+1)+j*(iZ+1)+k)=1
   pTemp(ubound(pTemp))=Rhino.AddPoint(pSect(i,1))
   redim preserve pTemp(ubound(pTemp)+1)
   NEXT
END IF
NEXT
Rhino.DeleteObject sectLine
  NEXT
NEXT
'第三轮：平行Y轴
FOR k=0 TO iZ
  FOR i=0 TO iX
pStart=Array(pBox(0)(0)+i*lCube, pBox(0)(1), pBox(0)(2)+k*lCube)
pEnd=Array(pBox(0)(0)+i*lCube, pBox(3)(1), pBox(0)(2)+k*lCube)
sectLine=Rhino.AddLine(pStart,pEnd)
FOR m=0 TO nSfs
pSect=Rhino.CurveSurfaceIntersection(sectLine, sfs(m))
IF IsArray(pSect) THEN '有交叉点
   nSectPoints=UBound(pSect,1) '交叉点个数
   FOR j=0 TO nSectPoints
   iRoundUp=int((pSect(j,1)(1)-pStart(1))/lCube+0.5)
   p(i*(iY+1)*(iZ+1)+iRoundUp*(iZ+1)+k)=1
   pTemp(ubound(pTemp))=Rhino.AddPoint(pSect(j,1))
   redim preserve pTemp(ubound(pTemp)+1)
   NEXT
END IF
NEXT
Rhino.DeleteObject sectLine
  NEXT
NEXT
redim preserve pTemp(ubound(pTemp)-1)
'根据记录的方块属性画方块
dim lCube2  '缩后的Cube边长
lCube2=rShrink*lCube
FOR i=0 TO iX
  FOR j=0 TO iY
FOR k=0 TO iZ
IF p(i*(iY+1)*(iZ+1)+j*(iZ+1)+k)=1  THEN
   pCube(0)=Array(pBox(0)(0)+i*lCube-lCube2/2,pBox(0)(1)+j*lCube-lCube2/2,pBox(0)(2)+k*lCube-lCube2/2)
   pCube(1)=Array(pBox(0)(0)+i*lCube+lCube2/2,pBox(0)(1)+j*lCube-lCube2/2,pBox(0)(2)+k*lCube-lCube2/2)
   pCube(2)=Array(pBox(0)(0)+i*lCube+lCube2/2,pBox(0)(1)+j*lCube+lCube2/2,pBox(0)(2)+k*lCube-lCube2/2)
   pCube(3)=Array(pBox(0)(0)+i*lCube-lCube2/2,pBox(0)(1)+j*lCube+lCube2/2,pBox(0)(2)+k*lCube-lCube2/2)
   pCube(4)=Array(pBox(0)(0)+i*lCube-lCube2/2,pBox(0)(1)+j*lCube-lCube2/2,pBox(0)(2)+k*lCube+lCube2/2)
   pCube(5)=Array(pBox(0)(0)+i*lCube+lCube2/2,pBox(0)(1)+j*lCube-lCube2/2,pBox(0)(2)+k*lCube+lCube2/2)
   pCube(6)=Array(pBox(0)(0)+i*lCube+lCube2/2,pBox(0)(1)+j*lCube+lCube2/2,pBox(0)(2)+k*lCube+lCube2/2)
   pCube(7)=Array(pBox(0)(0)+i*lCube-lCube2/2,pBox(0)(1)+j*lCube+lCube2/2,pBox(0)(2)+k*lCube+lCube2/2)
   Rhino.AddBox pCube
END IF
NEXT
  NEXT
NEXT
Rhino.DeleteObjects sfs '删除原始物体
Rhino.DeleteObjects pTemp '删除点
makeMosiac=sfs
END FUNCTION

作者: 夜神 时间: 2010-6-22 16:13
这是个简化版，是为了方便初学者阅读。
还有个效率更高的版本，代码长度有这个三倍长。

作者: yanhui314 时间: 2010-6-22 17:00
唉，我还不能上升到简化的程度，只求能做出来，呵呵

作者: njyqqq 时间: 2010-6-22 17:29
简化版和那个版本效果一样？？差别在那？？

作者: 夜神 时间: 2010-6-22 19:56

简化版和那个版本效果一样？？差别在那？？
njyqqq 发表于 2010-6-22 17:29

登录/注册后可看大图

功能都一样，不打折扣，只是高效版本的速度要快的多。当然，作为代价，程序复杂了许多。

作者: 城市边缘工厂 时间: 2010-6-22 21:30
楼主犀利。进来看看。

作者: panhao1 时间: 2010-6-22 21:43
就单纯做一个类似的效果而言这种算法显然过于冗长
使用Rhino.CurveSurfaceIntersection
这个就很耗时间

为什么不用判定点是否在brep内部那个命令呢？
那个命令就可以判定点是否在brep表面
完全被有必要拆分brep 再一个个面去求

作者: 夜神 时间: 2010-6-22 22:25
本帖最后由夜神于 2010-6-22 22:28 编辑

就单纯做一个类似的效果而言这种算法显然过于冗长
使用Rhino.CurveSurfaceIntersection
这个就很耗时间

为什么不用判定点是否在brep内部那个命令呢？
那个命令就可以判定点是否在brep表面
完全被有必要拆分 ...
panhao1 发表于 2010-6-22 21:43

登录/注册后可看大图

有这个命令吗？晕！
我在另一个程序里写过这个算法，就是为了避开使用Rhino.CurveSurfaceIntersection

作者: stratagem 时间: 2010-6-23 00:05
只能路过看看。。。。

作者: claudemit 时间: 2010-6-23 00:17
爬进来学习 lz大好人啊

作者: uuuf8 时间: 2010-6-24 02:30
本帖最后由 uuuf8 于 2010-6-24 02:36 编辑

确实没必要拆分

ps:虽然做的东西类似但是彼此的逻辑还是有很大不同啊这点很有意思

作者: 李薇薇 时间: 2010-6-24 17:57
好有趣，想学

作者: 董羽天 时间: 2010-6-25 12:40
嗯学习了

作者: divedragon 时间: 2010-6-26 12:42
感谢楼主分享啊。。。。。。。。

作者: qiqiyunyan1986 时间: 2010-6-29 19:05
看着就好复杂

作者: sharkliang 时间: 2010-10-23 17:21
感谢楼主分享啊

作者: taliesin 时间: 2010-10-26 10:12
谢谢楼主分享好东西大家学习

作者: liuchong625 时间: 2010-11-17 15:03
还是不错的说！

作者: Viking 时间: 2010-11-17 17:00
啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊
嗷嗷嗷嗷嗷嗷嗷嗷嗷嗷嗷嗷嗷嗷嗷嗷嗷嗷嗷嗷嗷嗷嗷嗷嗷
甚好甚好甚好甚好甚好甚好甚好甚好甚好甚好甚好甚好甚好
这货是什么?! lz 我爱你!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

作者: Viking 时间: 2010-11-17 17:02
强烈要求斑竹再给他加n分!
喵了个咪真是赞哉! 哦喵~

作者: ccocean 时间: 2010-11-17 21:56
谢谢楼主，崇拜啊。。。。什么我才能写出这样的代码？

作者: huangchang0528 时间: 2010-12-11 23:03
这个····太厉害了

作者: triko 时间: 2010-12-31 17:19
完全是看热闹

作者: ukcolour 时间: 2011-1-30 08:37
很不错的一段

作者: renping88 时间: 2011-3-20 02:05
MAD中钢大厦表皮的模拟

作者: bxsqrym 时间: 2011-4-2 11:44
好复杂哈倒了

作者: gaga 时间: 2011-4-2 22:45
貌似很复杂的样子

作者: gaga 时间: 2011-4-5 11:10
好漂亮啊，哪天我才会啊

作者: bxsqrym 时间: 2011-4-10 19:00
NNNNNNBBBBBB

作者: renping88 时间: 2011-4-13 04:32
ds w

作者: renping88 时间: 2011-4-13 04:32
w ew

作者: wawa 时间: 2011-4-28 12:34
very good!!!!!!!!

作者: kewangzhenxiang 时间: 2011-4-28 13:17
我没照度了，，，送电给我吧。。。

作者: mqcjackson 时间: 2011-10-14 14:03
我觉得做得挺好的……小白~~

作者: anboyh 时间: 2011-11-7 15:23
谢谢楼主分享。

作者: M-chibitora 时间: 2011-11-8 00:46
牛人啊！！！学习了哈哈

作者: hying 时间: 2011-12-10 14:02
顶。。。！

作者: howldyx 时间: 2011-12-14 21:09
感谢楼主分享啊。。。。。。。。

作者: howldyx 时间: 2011-12-14 22:07
收藏之哈哈哈

作者: s.k. 时间: 2011-12-14 23:09
我要学编程啊

作者: jasonbb88 时间: 2011-12-27 11:23
厉害~~~~~~~

作者: 蚊子 时间: 2011-12-27 12:54
什么都看不见啊

作者: bizquit 时间: 2012-1-1 17:12

楼主强大啊

作者: 鹰眼黑客 时间: 2012-2-10 23:48
楼主的分享精神值得赞扬。

作者: 多吃水果蔬菜 时间: 2012-2-28 11:44
做的真好！羡慕你！

作者: liuliu8244 时间: 2012-10-16 21:19
学习一下，lz继续加油

作者: 阿兹海默青年 时间: 2014-1-8 00:37
到现在都还没这么牛逼

作者: 左手◆23年 时间: 2014-2-13 09:43
感谢楼主分享，赞一个

作者: mgltju 时间: 2022-7-11 09:18
论坛高手云集，厉害厉害

欢迎光临 NCF参数化建筑论坛 (http://ncf-china.com/)