TP9 submission deadline is Thursday, 13/04/2017 at 23:59. The TP on Friday 14/04 will start at 13:30 (instead of 15:15).

TP6 : Bézier Surfaces

17 March 2017 in


Tensor-product Bézier surfaces

Tensor product Bézier surfaces are a direct generalization of Bézier curves. Instead of a control polygon $\mathbf b_i$, we consider a control net $\mathbf b_{ij}$. A Bézier surface is the described as

where $m,n$ are degrees in $u,v$, respectively; $B_{l}^{k}$ are Bernstein polynomials.


To evaluate a point $S(u,v)$ on a surface, let’s fix $j$ and vary $i$.

\begin{equation} \displaystyle B(u,v) = \sum_{j=0}^{n} B_{j}^{n}(v) \underset{= \mathbf b_{j}(u) } { \underbrace{ \left[ \sum_{i=0}^{m} \mathbf b_{ij} B_{i}^{m}(u) \right] } } \end{equation}

The $\mathbf b_{j}(u)$ defines a Bézier curve in $u$ and can be evaluated using the De Casteljau.

\begin{equation} \displaystyle B(u,v) = \sum_{j=0}^{n} \mathbf b_{j}(u)B_{j}^{n}(v) \label{final} \end{equation}

This equation defines a Bézier curve in $v$ with control points $\mathbf b_{j}(u)$ depending on $u$. At the end of the day, we have

  • $n+1$ evaluations of the De Casteljau for the degree $m$;
  • $1$ evaluation of the De Casteljau for the degree $n$.

Alternatively, we can fix $i$ and vary $j$.


Bézier patch evaluation scheme. (image by Pierre-Luc Manteaux)

Coordinate matrices

In the code, the control points $\mathbf b_{ij}$ are actually stored in three coordinate matrices $\texttt{Mx}, \texttt{My}, \texttt{Mz} \; $ so that

This means the code is more comprehensible as the structure of the matrices directly represents the grid topology of the patch. On the other hand, it also means the computation needs to be done for each coordinate individually.

Bézier surface vs. piecewise Bézier surface

In practice, a Bézier surface often consists of multiple surface patches, each having its own control polygon. Therefore, it is sometimes called a piecewise Bézier surface.

Some surfaces in the data/ folder consist of more than one patch, ranging from 2 (heart) to 32 (teapot). They are saved in the BPT format; some datafiles are taken from the website of Ryan Holmes where you’ll also find more details about this format in case you’re interested.

teapot isophotes teapot isophotes

The discontinuity in isophotes shows the piecewise Bézier Utah teapot is not $\mathcal C^1$.


  1. Implement the evaluation of Bézier surfaces for $(u,v) \in [0,1]^2$. Use simple and wave for first tests (these contain only one patch).
  2. When you’re sure the implementation works for the simple cases, test your algorithm on datasets with multiple patches: heart (2), sphere (8), teapot (32), teacup (26), teaspoon (16). Don’t set the density parameter too high, always start with smaller values (5 or 10) as the number of computed points is density².