Bézier表面修补表面

Question

使用“BezierFunction”生成表面形状立方Bézier贴片（作为3D点列表）很容易在Mathematica中

示例：

pts = {{{0, 0, 0}, {0, 1, 0}, {0, 2, 0}, {0, 3, 0}}, {{1, 0, 0}, {1, 
     1, 1}, {1, 2, 1}, {1, 3, 0}},
   {{2, 0, 0}, {2, 1, 1}, {2, 2, 1}, {2, 3, 0}},
   {{3, 0, 0}, {3, 1, 0}, {3, 2, 0}, {3, 3, 0}}};

f = BezierFunction[pts]

Show[Graphics3D[{PointSize[Medium], Red, Map[Point, pts]}],
 Graphics3D[{Gray, Line[pts], Line[Transpose[pts]]}], 
 ParametricPlot3D[f[u, v], {u, 0, 1}, {v, 0, 1}, Mesh -> None]]

但现在我有一个相反的问题，如果我们有一个3D表面，我们如何得到一个基于它的立方Bézier补丁（按照一些规则）？