CS 284: CAGD 
Lecture #6 -- Th 9/14, 2006.

Preparation:

Analysis: DoFs in an Interpolation Problem

Main Topic: Comparison of Interpolation Methods.

Some Comments on Bezier-stitching Homework

• Not as easy as it first seems -- everybod showed good effort, -- not everybody succeeded completely, -- several of you did a great investigation!
  
• First crucial step is to get the tangent direction right.
• Line parallel to P₀ - P₂ does not lead to great solutions, but it can maintain affine invariance.
• Using a tangent perpendicular to the angle-bisector in the control polygon, gives better results.
• For even better results, some more weight should be given to the shorter segment to define tangent direction;
it is an open-ended question how this should be done for the very best, most robust results ...
(perhaps divide the turning angle between the control polygon segments proportional to the ratio of the segment lengths).
• For G¹-continuity the Bezier control points adjacent to the joint need to lie both on this tangent line.
• They can be at un-equal distances; C¹-continuity does NOT lead to better looking shapes, it just adds one more constraint!
• The distances should be scaled (somewhat) with the length of the underlying control segment.
• The construction should be structurally symmetrical, i.e., use essentially the same formula for both inner control points.

Compare your own solution with the different methods above:

• Trigonometric angle-interpolated circle splines, reparameterized for C²-continuity [Sequin et al]
• Linear angle-interpolated circle splines [Sequin et al]
• Trigonometric positional circle interpolation [Szilvasi-Nagi]
• Linear positional circle interpolation [Wenz]
• Non-uniform Lagrange interpolation.

Bottom Line: Polynomials are not the only possible primitives !

Circle Splines

• SIGGRAPH'03 Presentation
• Kiha's demos: use Explorer!

An interpolation task (see B+B+B, Chapter 3):  Results with different approaches:

• Cubic Hermite interpolation [see B+B+B, Chapter 3]:
• Continuity =?; DoF = ?; Extent of support = ?.
• Non-uniform Lagrange interpolation:
  
• Continuity = ?; DoF = ?; Extent of support = ?.
• Linear positional circle interpolation [Wenz]:
• Continuity = ?; DoF = ?; Extent of support = ?.
• Trigonometric positional circle interpolation [Szilvasi-Nagi]:
• Continuity = ?; DoF = ?; Extent of support = ?.
• Linear angle-interpolated circle splines [Sequin et al]:
• Continuity = ?; DoF = ?; Extent of support = ?.
• Trigonometric angle-interpolated circle splines, reparameterized for C2-continuity [Sequin et al]:
• Continuity = ?; DoF = ?; Extent of support = ?.
• Your stitched-together composite Bezier curve:
• when G¹-continuity is OK: --> DoF = ?; Extent of support = ?. 
• when C¹-continuity is required: --> DoF = ?; Extent of support = ?.

Minimum-Variation Curves

• Choose any suitable representation
• Minimize the arc-length integral over square of derivative of curvature, by using the un-constraint DoF.
• One possible approach uses quintic Hermite splines [Moreton and Sequin],
and a good heuristic to set the initial tangent directions (1st derivatives) and curvatures (2nd derivatives) at the interpolation points.

Next Reading Assignments:

Extended Homework Assignment:  Smooth Sweeps along Complex 3D Curves

   1. Implement smooth curve interpolation in 3D.
     - - Start from the file http://www.cs.berkeley.edu/~sequin/CS284/CODE/pa2.slf
     - - Write the corresponding formulae for the z-components of the inner Bezier control points.
     - - Add the necessary z-expressions in all places where I have not already done so.
     - - I have added a "Z-scale" slider that multiplies all the new z-coordinates with a slider-controlled constant,
          so that you can readily go back and forth between 2D (multiplier=0.0) and 3D (multiplier=1.0).
     - - Lift the data points out of the x-y-plane, and convince yourself that everything still "works smoothly."

   2. Use the SLIDE "sweep" construct to form a small ribbon going along a closed B-spline curve approximating the original 15 control points.
     - - This can be done by sweeping a "star" type cross section, and giving the star two skinny spokes.
     - - I have already put in most of the needed code; if you turn on "drawSweep" in the "slf_swp" slider menu,
          you should see such a ribbon in red, displayed edge-on.
     - - However, the ribbon seems to have 10 pinched-off, twisted points in its loop; WHY is this ? -- Write down your thoughts.
     - - Under what conditions do you get a nice "prismatic" sweep with no funny points ?
     - - How do you have to set parameters such as : "minimizeTorsion" "closed (curve)" "symmetry" "azimuth" ...
          in order to achieve a nice result ?
          (It is advisable to consult the manual pages on SLIDE ...)
     - - Now lift the original control points out of the x-y-plane by increasing the "Z-scale" slider at the bottom of the "slf_bez" slider menu..
     - - How is the behavior of the sweep different now ?
     - - Pay particular attention to the end condition where the loop closes; under what circumstances do you get a smooth, seamless closure ?
     - - Answer the same questions also for a 5-pointed star (set "Spokes" to 5).

ATTENTION: -- the "closed" flag does not seem to work when set dynamically: 
closed {expr $oSweep(closed)}
but will work when set explicitly in the SLIDE file, either as:
closed 0
or:
closed 1

   3. Now try to make such a nice sweep along your own smooth, composite, interpolating curve.
     - - It is somewhat tedious, since you have to specify a separate sweep construct for each of the 15 Bezier segments.
     - - But again, the code is already in place, and you can see such a ribbon in green, if you turn on "drawSweep" in the "slf_bez" menu.
     - - You may now want to turn off the B-spline at this time, so that you can focus on the Bezier segments.
     - - For a 2D curve, as long as the ribbon is perpendicular to the curve's plane, everything should work fairly nicely.
     - - But there may be problems in 3D. -- WHY ? What can you do about them ?
>>> Just discuss the issues you encounter in the 3D case. Find out what "hooks" are or are not available in the SLIDE sweep.
Briefly summarize what would have to be done to make nice prismatic sweeps along a composite Bezier curve in 3D,
and how the necessary information may be computed.
No need to actually do all the programming -- unless you want to.

   4. Now use all the originally given points, as well as your own calculated control points,
      to make a closed polyline with 45 joints, and then sweep a 5-pointed star along that polyline.
     - - The code for this is right next to the B-spline "backbone" code (near line 467 in pa2.slf).
     - - Under what conditions do you get nice prismatic beams with nicely mitered joints ?
     - - How do you have to set parameters such as : "minimize Torsion" "closed (curve)" "symmetry" "azimuth" ...
          in order to achieve a nice result ?
     - - Can the "twist" parameter be of any use in achieving nice, properly mitred joints ? HOW ?
>>> Implement the sweep for the polyline case. The SLIDE machinery should handle all that.

DUE: Thu. 9/21/2006, 2:10pm.
Set the initial values for the sliders to the proper values, so that when your program is run,
one will see a nice "prismatic" sweep of a 3-pointed star along your best composite Bezier curve.
(You may want to change some of the parameters from the values that you gave them in pa1...)
 >>> Without the implementation of (3), this will not work well in 3D.
But the 2D case is much simpler. You should be able to get nice results for the 2D case.

• E-mail me a copy of your final modified pa2.slf SLIDE file.
  

Capture a hardcopy print-out of this sweep.

Hand in at the beginning of class:

• A hardcopy print-out of your sweep;
• and a page or two of text, answering the questions raised above.