SoftimageScript : Fillet っぽいなにか

メモです

ちょっと　パイプが山の用にあって
その曲がりの角ががが・・・　orz

Filletっぽい何かを作ってお茶を濁します。

FilletLinearCurve

リニアカーブのコントロールポイントを選択し、
実行します

カーブの再描画ですので、ご注意を。

# FilletLinerCurve
#
#　リニアカーブの選択したポイントを分割、再描画してくれます
#　正確な円ではありません～　付け焼刃
#
#　v0.1 2013 0417

app = Application; log = app.LogMessage
from win32com.client import constants as c

 
#ログの表示を消す
cmdLog = app.Dictionary.GetObject( "preferences.scripting" )
cmdLog.cmdlog.value = False

#半径と、間に入る頂点数
gRadius   = 1
gDivision = 2


def CreateNewCurve(aPos):
 
 myNewCurve = app.SICreateCurve("",1,0)
 for i in range(len(aPos[0])):
  app.SIAddPointOnCurveAtEnd(myNewCurve,aPos[0][i],aPos[1][i],aPos[2][i],False,0)
 app.SelectObj(myNewCurve)
 return


def l2gArray(aPos,m4):
 v = XSIMath.CreateVector3()
 for i in range(len(aPos[0])):
  v.Set(aPos[0][i],aPos[1][i],aPos[2][i],)
  v.MulByMatrix4InPlace(m4)
  aPos[0][i] = v.X
  aPos[1][i] = v.Y
  aPos[2][i] = v.Z


def alog(aPos):
 log("cnt = " + str(len(aPos)))
 for i in range(len(aPos)):
  log( "V : " + 
  str( aPos[0][i]) + " , " +
  str( aPos[1][i]) + " , " +
  str( aPos[2][i])
  )


def GeiFilletPoint(posA,posB,posC,in_t):

 v1 = XSIMath.CreateVector3()
 v2 = XSIMath.CreateVector3()
 v3 = XSIMath.CreateVector3()
 
 p0 = XSIMath.CreateVector3()
 p1 = XSIMath.CreateVector3()
 p2 = XSIMath.CreateVector3()
 
 v1.Sub(posB,posA)
 v2.Sub(posC,posB)
 
 p1.ScaleAdd(in_t,v1,posA)
 p2.ScaleAdd(in_t,v2,posB)
 
 v3.Sub(p2,p1)
 p0.ScaleAdd(in_t,v3,p1)
 
 return p0
 
 
def addFiletPnts(x,y,z,i,aPos):
 
 #カーブのセグメントを追加 
 #ContorolPoints 0,1,2
 pos00 = XSIMath.CreateVector3(aPos[0][i - 1],aPos[1][i - 1],aPos[2][i - 1])
 pos10 = XSIMath.CreateVector3(aPos[0][i    ],aPos[1][i    ],aPos[2][i    ])
 pos20 = XSIMath.CreateVector3(aPos[0][i + 1],aPos[1][i + 1],aPos[2][i + 1])

 v1 = XSIMath.CreateVector3()
 v2 = XSIMath.CreateVector3()
 
 v1.Sub(pos00,pos10)
 v2.Sub(pos20,pos10)
 
 v1.NormalizeInPlace()
 v2.NormalizeInPlace()

 #new Point 
 pos01 = XSIMath.CreateVector3()
 pos02 = XSIMath.CreateVector3()
 
 #半径分Pos00からの距離
 pos01.ScaleAdd(gRadius,v1,pos10)
 pos02.ScaleAdd(gRadius,v2,pos10)
 
 #分割数分 FilletPositonの追加
 num = gDivision + 2
 
 for i in range(num):
  t = i * 1.0 / (num - 1)
  v = GeiFilletPoint(pos01,pos10,pos02,t)
  x.append( v.X )
  y.append( v.Y )
  z.append( v.Z )


def CreateNewPos(aPos,oSelIDs):

 x = []
 y = []
 z = []
 
 for i in range(len(aPos[0])):
  log(i)
  flg = True
  for j in oSelIDs:
   #最初と最後のポイントは無視（Closeには対応してません
   if i == j and i != 0 and i != len(aPos[0]) - 1:
   
    #選択したポイントであれば
    #log("match")
    flg = False
    
    #カーブのポイントを挿入
    addFiletPnts(x,y,z,i,aPos)
    
  if(flg == True):
   #log("NoMatch")
   x.append( aPos[0][i] )
   y.append( aPos[1][i] )
   z.append( aPos[2][i] )
   
 nPos = [x,y,z]
 
 return nPos
 

def main():
 if app.Selection.Count ==0:
  return
  
 oSel = app.Selection(0)
 if oSel.Type != "pntSubComponent":
  log("Select LinerCurve Control Points")
  return
  
 oSub = oSel.SubComponent
 oCurveList = oSub.Parent3DObject
 if oCurveList.Type != "crvlist":
  log("Select LinerCurve Control Points")
  return
  
 oColl = oSub.ComponentCollection
 oSelIDs = oColl.IndexArray
 oTrans = oCurveList.Kinematics.Global.Transform
 m4 = XSIMath.CreateMatrix4()
 oTrans.GetMatrix4(m4)
 
 #カーブからコントロールポイントのリスト、Posを取得
 oCurve = oCurveList.ActivePrimitive.Geometry.Curves(0)
 aPos = oCurve.ControlPoints.PositionArray
 aPos = [list(aPos[0]),list(aPos[1]),list(aPos[2])]
 
 l2gArray(aPos,m4)
 
 #新しいコントロールポイントのリストを作成
 nPos = CreateNewPos(aPos,oSelIDs)
 
 #新規カーブの作成
 CreateNewCurve(nPos)
 
main()

スミマセン、こうでもして書いておかないと
多忙なときは作ったことすら忘れますのでご容赦を