I have an application that uses bezier triangle patches to define a hexagonal game board. I am finding that in certain situations the lack of local control provided by bezier surfaces is creating unacceptable visual artifacts. I am now hoping to experiment with triangular NURBS surfaces to see if I can solve some of these issues.

I found a paper that seems rather thorough on the subject here, but I'm have trouble understanding some of the math terminology with my lack of formal math training. Here are a list of things from section 2.2 in that paper that I don't understand.

In the first paragraph of section 2.2 it says "Let T = {(delta?)(i) = [r,s,t] | i = (i0,i1,i2) (thing that looks like an E) Z3+". I assume that r,s,t are the triangle vertices and that the index i is a number combo like 300, 030, 003 etc. If anyone can explain some of these symbols I would be grateful. I don't know what the vertical line means either.

After that there is a knot sequence, which I understand is what is called a knot vector elsewhere, and is just a sequence of increasing numbers like {0,0,1,1,2,2} etc.

Then it says we define a convex hull, and the beta symbol seems to be sort of like the index, I don't know what the 'absolute value' lines mean in this context but that symbol seems to have 3 components that add up to the 'order' of the triangle. I am unsure how the index and this beta symbol should be combined to get the desired number sequence.

Then it defines a multivariate simplex spline M(u | Vibeta). The u is a two component coordinate (u,v), but I don't know what the line means.

Then it defines the basis function with that line again, then a lower case d, with the triangle vertices without the indices but with the beta thingies in parenthesis, 'lined' against the M(u | Vibeta). What does that line mean anyway? The lower case d is also unexplained.

Then we have the unweighted B-spline surface. The Sigma operator is something i always have trouble with, but I interpret this as the sum of every index combo multiplied by the point and the basis function and multiplied again with the sum of all the beta combos multiplied by the point and basis function. However, I do not understand what is meant exactly by Pi,|beta|. Pi would be P300, P030, P003 etc, but I don't know what the beta subscript does to that interpretation.

Finally, there are arbitrary weights that can be set, and once I understand the indexing system, I think I understand this part as a weighted mean of the above.

If anyone can help me make sense of this I would be super grateful.

Thanks so much! I will look at that paper and see if that makes more sense. I have a couple of questions.

When you said 'set of triangles', is that triangles within the patch or neighbor triangles?

Also, from what you can tell, how similar are triangular nurbs surfaces to rectangular ones? I can probably find better explanations of the rectangular ones, hopefully with pictures, but I wanted to get right into the thing I need.

These papers would be so much more useful if they had some example code. English is a very ambiguous language for describing mathematical concepts.

I'm looking at the paper linked by LorenzoGatti, Triangular NURBS Surface Modeling of Scattered Data and I'm trying to understand the calculation of the basis functions.

In section 2, they describe how many knots are needed per vertex and how they can be shared among adjacent triangles, then they take 5 of the 9 knots, and arrange them in 6 configurations which seems to be every possibility of a combination of 5 knots. They say nothing of how the extra knots are discarded or if that matters. Each of the 6 scenarios requires a basis function defined from each 5 knot set K as:

B(u|K) = a0B(u|K\v0) + a1B(u|K\v1) + a2B(u|K\v2)

'u' is a point in the 2D domain (cartesian or barycentric?) and 'a' is related? to the barycentric coords. In the context of an arbitrary collection of multiple triangles, I don't know what that means. Then there is a bit about how the last term in the above equation can either be 1 if the 'u' point are in the triangle formed by the three vertex knots and 0 if it is not.

I can't begin from this to understand how to put this in code.

EDIT: by the way, for some context regarding my intentions, my surface patch is a triangle with 10 control points forming 9 sub-triangles like so:

[attachment=18991:np2.gif]

Take my post with a big pinch of salt, I may be entirely wrong here (and welcome anyone to correct me entirely!)

First thing to note: in the paper they are using quadratic splines, not cubic ones. As a result the numbers of their control points and knots won't match the configuration you have above.

I 'think' this is their equivalent of their cox-de-boor (ish) for a triangle. I'm guessing the 3 terms on the right [B(u|K\v0), B(u|K\v1), and B(u|K\v2)] are simply the standard cox-de-boor algorithms?

When you say that the paper is using quadradic splines, you mean that the entirety of the basis function has only to do with the sub triangles, correct? My patch is currently based on a cubic bezier, where moving any point alters the entire patch, but since NURBS offer local control, I can treat the sub triangle basis functions separately correct? What is mean by quadratic and cubic in this case?

Wow...a paper that combines computational geometry, Lagrangian mechanics, and finite element theory. That's impressive. Anyway, it's hard to decipher what's actually going on through the math. I can't see how the basis functions are constructed from the paper.

Have you tried creating a regular square NURBS surface but making all of the control points along one edge the same point in 3D space? The parameter space would be square, but geometrically it will be triangular. That would be a simpler way to do it, IMHO.

Well, I'm getting to the point where I think I just have to start implementing to the best of my understanding, and find out how to do it experimentally. Then once I understand how to get the positions, I have to find the partial derivatives to get the normal. Blegh!