def spl_face(px, py, pz):
nx, ny = px.shape
pnt_2d = TColgp_Array2OfPnt(1, nx, 1, ny)
for row in range(pnt_2d.LowerRow(), pnt_2d.UpperRow() + 1):
for col in range(pnt_2d.LowerCol(), pnt_2d.UpperCol() + 1):
i, j = row - 1, col - 1
pnt = gp_Pnt(px[i, j], py[i, j], pz[i, j])
pnt_2d.SetValue(row, col, pnt)
#print (i, j, px[i, j], py[i, j], pz[i, j])
api = GeomAPI_PointsToBSplineSurface(pnt_2d, 3, 8, GeomAbs_G2, 0.001)
api.Interpolate(pnt_2d)
#surface = BRepBuilderAPI_MakeFace(curve, 1e-6)
# return surface.Face()
return BRepBuilderAPI_MakeFace(api.Surface(), 1e-6).Face()
def spl_face(px, py, pz, axs=gp_Ax3()):
nx, ny = px.shape
pnt_2d = TColgp_Array2OfPnt(1, nx, 1, ny)
for row in range(pnt_2d.LowerRow(), pnt_2d.UpperRow() + 1):
for col in range(pnt_2d.LowerCol(), pnt_2d.UpperCol() + 1):
i, j = row - 1, col - 1
pnt = gp_Pnt(px[i, j], py[i, j], pz[i, j])
pnt_2d.SetValue(row, col, pnt)
#print (i, j, px[i, j], py[i, j], pz[i, j])
api = GeomAPI_PointsToBSplineSurface(pnt_2d, 0, 0, GeomAbs_C0, 0.001)
api.Interpolate(pnt_2d)
#surface = BRepBuilderAPI_MakeFace(curve, 1e-6)
# return surface.Face()
surf = api.Surface()
surf.Transform(set_trf(gp_Ax3(), axs))
face = BRepBuilderAPI_MakeFace(surf, 1e-6).Face()
shel = BRepBuilderAPI_MakeShell(surf).Shell()
return surf, face, shel
