TweetFollow Us on Twitter

Line Art Rotation
Volume Number:6
Issue Number:5
Column Tag:C Forum

Related Info: Quickdraw

Line Art Rotation

By Jeffrey J. Martin, College Station, TX

Note: Source code files accompanying article are located on MacTech CD-ROM or source code disks.

[ Jeff Martin is a student at Texas A&M University working on his bachelors in computer science. He has been a personal computer technician at the campus computer center, a system operator on the campus mainframes, and now freelances graphic work for various professors. He hopes that one day a motion picture computer animation company will take him away from all of this.]

This being my first stab at an article, I will try to keep it short while leaving in all of the essential vitamins and nutrients. In that spirit my user interface will bring back nostalgic thoughts to those past Apple II and TRS-80 users, and any PC people will feel right at home.

The essence of this program is to show how a seemingly complicated transformation and rotation can be applied to an array of points that form any arbitrary line art.

Of course to form a transformation on the array of points (e.g. offset the points to the left) we simply add some delta x(dx) and/or delta y(dy) to every point:

/* 1 */


Now rotation is a little harder, but to spare you the heartache, it can be shown that for rotation about the origin(fig 1):

So the trick of rotating about some arbitrary point is to first transform that pivot point to be the origin(transforming every other point by the save amount). Second, perform the rotation of all points by the angle theta. Third, transform the pivot back(once again transforming all other points as well).

Now all of this may seem to be a costly maneuver, but the fact is that we can roll all of these into a single matrix multiplication, using homogeneous coordinates:


form one matrix.

Fig. 2 shows the multiplication of a homogeneous coordinate and a translation matrix. Please verify that this results in (X+dx,Y+dy) (if unfamiliar with matrix multiplication see mult procedure in program).

Similarly figure 3 shows multiplication with a rotation matrix - an exact translation of our rotation equations in matix form.

So the translation, rotation, and inverse translation matrices are as shown in figure 4. Which forms one matrix to be multiplied times the vertices.

The following program allows the user to enter in points with the mouse until a key is pressed. At that time the user then uses the mouse to enter a pivot point. The program uses the pivot point to form the translation and inverse translation matrices(from the x and y coordinates). The program then forms a rotation matrix of a constant rotation angle(Π/20) and calculates the new vertices based on the values of the old ones. The program undraws the old lines and redraws the new and calculates again until the object has rotated through a shift of 4Π(2 rotations). press the mouse button again to exit program.

Once again, I point out that the code does not follow the user guidelines, but then it is not exactly meant to be an application in itself. Build your own program around it and see what you can do. One suggestion is to cancel the erasing of the object to achieve spirograph patterns. I think too many of the submissions to MacTutor contain an interface that we all know too well, and for those just interested in the algorithms it can mean a lot of extra work. Have Fun.

/* 2 */

int errno;

void mult();  /*out matrix mult proc*/
/*floating value of points to avoid roundoff*/
typedef struct rec {float h,v;} points;
  int buttondown=0, /*flagg for mouse       */
      n=-1,         /*number of vertices    */
      keypressed=0, /*flagg for key         */
      flip=0,       /*to allow alternating  */
      flop=1,       /*vertices to be drawn  */
      i;            /*array counter         */
  float x,          /*angle counter         */
      T[3][3],      /*translation matrix    */
      Tinv[3][3],   /*translate back        */
      Rz[3][3],     /*rotate matrix         */
      c[3][3],      /*result of T&R         */
      d[3][3];      /*result of c&Tinv      */
  long curtick,     /*for delay loop        */
       lastick;     /*for delay loop        */
  EventRecord nextevent;/*to get mouse&key  */
  Point origin,dummy;   /*pivot and locator */
  points points[2][30];/*vertices(don’t draw Eiffel tower)  */
  WindowPtr scnwdw;    /*window pointer     */
  Rect      scnrect;   /*window rect        */
*  Set things up                     *
scnwdw=NewWindow(0,&scnrect,”\p”,TRUE,dBoxProc, -1,FALSE,0);
*  Get points                        *
      if(nextevent.what==mouseDown) buttondown=1;
      else if(nextevent.what==keyDown) keypressed=1;
    if(buttondown) /*get a point and draw it*/ 
    } /*end of get point*/
  }  /*end of get points*/
*  Get origin                        *
      if(nextevent.what==mouseDown) buttondown=1;
*  Make translation matrix           *
*  Rotate                            *
  x=0.157;  /*rotation angle - about 9 degrees*/
    }  /*end update points*/
    ForeColor(whiteColor);  /*undraw flop*/
    lastick=TickCount(); /*time delay for retace to improve animation*/
    do{curtick=TickCount();} while(lastick+1>curtick);
    for(i=1;i<=n;i++) LineTo((int)points[flop][i].h,(int)points[flop][i].v);
    ForeColor(blackColor);  /*draw flip*/
    do{curtick=TickCount();} while(lastick+1>curtick);
    for(i=1;i<=n;i++) LineTo((int)points[flip][i].h,(int)points[flip][i].v);
  }  /*end rotate*/
*  End everything                    *
      if(nextevent.what==mouseDown) buttondown=1;
}  /*program end*/

void mult(A,B,C)
  float A[][3],B[][3],C[][3];
  int i,j,k;
}  /*end mult*/

3D Modeling & Rotation

The main thrust of this exercise is to extend the line art rotation into 3D object rotation using the same techniques as the 2D, while also implementing parallel projection as our means of 3D modeling.

The first part of the exercise requires that we define an object in a structure that we can easily manipulate. Using a cube for simplicity, we will start by defining the center of the cube and an array of vertices, vertex[2][# of pts] (see GetPoints in program). Referring to fig. 1, each vertex corresponds to a corner of the cube. The second dimension of the array is to provide a destination for transformed vertices. Having both sets will allow us to undraw and immediately redraw the shape - minimizing the hangtime between redrawing allows for smoother animation.

Figure 1.

Next let us construct an array of lines connecting these vertices. Each element of the line array refers to the index of the beginning and ending vertex of that particular line. This array will never change. Think of when you roll a die - the edges still go between the same corners, but the position of the corners has changed.

The next construct is the translation and inverse translation matrixes. As in 2D rotation, we must transform our local center of rotation to the origin, rotate, then translate back.

The idea of homogeneous coordinates was introduced in the last article and is now extended into 3D by adding a fourth term. Fig. 2 shows our homogeneous coordinate as a 1x4 matrix times our translation matrix(4x4). The purpose of this multiplication is to add a dx, dy and dz to every point, in order to center our vertices about the origin. Please verify that the matrix multiplication results in X+dx,Y+dy,Z+dz (if unfamiliar with matrix multiplication see matmult in program).

Figure 2.

Now we once again reach the challenging concept of rotation. Although similar to 2D, we now have the option of rotating around the X and Y as well as the Z-axis.

The simplest, rotation about the z-axis, is just as in our 2D rotations, because none of the z-values change. If this is hard to understand, think about this: if you look straight down a pencil with the point a foot away from you and spin it a half turn, the point is still a foot away, but the writing is now on the other side. The equations for the changes in the X and Y are as follows:

  Xnew=XoldCos(Ø) + YoldSin(Ø)
  Ynew=-XoldSin(Ø) + YoldCos(Ø)

The 3D representation in matrix form with a vertex multiplication is in fig. 3. And the proof of all this is in that dusty old trigonometry book up on your shelf. (once again direct multiplication of fig. 3 will yield the preceding equations).

Figure 3.

Similarly rotation about the X axis changes none of the x-values, and rotation about Y changes none of the y-values. The transformation equations are given as follows:

Rotation about the X:

 Ynew=YoldCos(Ø) + ZoldSin(Ø)

Rotation about the Y:

   Xnew=XoldCos(Ø) - ZoldSin(Ø)
 Znew=XoldSin(Ø) + ZoldCos(Ø)

The corresponding matrices are shown in figures 4 and 5.

Figure 4.

Figure 5.

Once again we will construct a new array of vertices from a single transformation matrix formed from the translation to the origin, rotation about an axis, and translation back. Therefore creating the new vertices:


or after combining T*Rz*Tinv into a single Master Transformation(MT):


Finally the trick of parallel projection when viewing an object from down the Z axis is that all you have to do is draw lines between the x,y components of the points (ignore the z). For those mathematically inclined, you will realize that this is just the projection of those 3D lines on the X-Y plane (see fig. 6).

Figure 6.

The particular stretch of code I’ve included implements this transformation on the cube for rotation along the X and Y axes of the center of the cube using the arrow keys. The successive transformations of the vertices are loaded into the flip of the array (vertex[flip][pnt.#]). Then the flop is undrawn while the flip is drawn as mentioned previously and flip and flop are changed to their corresponding 0 or 1.

After launching, the application immediately draws the cube and then rotates it in response to the arrows. The program exits after a single mouse click.

Once again the code is not intended to match up to the guidelines - but is intended for use with other code or simple instructional purposes. It is concise as possible and should be easy to type in. A quick change to numofpts and numoflines as well as your own vertex and and line definitions would allow you to spin your favorite initial into its most flattering orientation.

The inspiration for this program came from the floating couch problem presented in Dirk Gently’s Holistic Detective Agency, by Douglas Adams. If enough interest is shown, perhaps a future article would include hidden line removal and color rendering techniques. After all, it was a red couch.

One last suggestion for those truly interested is to pull your shape definition in from a 3D cad program that will export in text format, such as Super 3D or AutoCad.

Anyway, on with the show

/* 3 */

/* Following is inline macro for drawing lines */
#define viewpts(s) {for(i=0;i<numoflns;i++)  \
                     { MoveTo((int)vertex[s][line[i].v1].x,  \
                       (int)vertex[s][line[i].v1].y); \
                       LineTo((int)vertex[s][line[i].v2].x, \
                       (int)vertex[s][line[i].v2].y); }}  
#define numofpts 8 /* A cube has eight vertices */
#define numoflns 12    /* lines for every face. */

/* the following are the data structs for vertices and lines*/ typedef 
struct rec1 {float x,y,z;} point3d;
typedef struct rec2 {int v1,v2;} edge;
void mult();/* Matrices multiplication */

  point3d vertex[2][8], /* array of 3D pts   */
          center;/* centroid of cube */
  edge    line[12];/* array of lines */
  int     buttondown=0, /* mousedwn flag(for prog end)*/
          keypressed=0,       /* keydwn flg(for arrows)     */
          flip=0,             /* This is index for vertex so*/
          flop=1,             /* can undraw flip & draw flop*/
          i,                  /* counter           */
          rot=0; /* Flag for direction of rotat*/
  long    low;   /* low word of keydwn message */
  float   a,/* Particular angle of rotat     */
          R[4][4], /* Rotation matrix*/
          c[4][4], /* Product of trans & rot mats*/
          d[4][4], /* Product of c and inv trans */
          T[4][4],Tinv[4][4], /* Translation & inv trans    */
          x=0.087266;/* Algle of rot in rad  */
  EventRecord nextevent;
  KeyMap    thekeys;
  WindowPtr scnwdw;
  Rect      scnrect;
*  Set things up *
*  Get points. Arbitrary cube.*

*  Rotate *
viewpts(flip);   /* This draws first set of pts*/
  while(!buttondown) /* Mini event loop*/
      if(nextevent.what==mouseDown) buttondown=1;
      else if(nextevent.what==keyDown) keypressed=1;
      else if(nextevent.what==autoKey) keypressed=1;
    if(keypressed) /* Find out which one     */
      if(low==126) {rot=1;a=-x;} /* Set dir flag and-*/
      if(low==124) {rot=2;a=-x;} /* angle(pos or neg */
      if(low==125) {rot=3;a=x;}
      if(low==123) {rot=4;a=x;}
        case 1:/* Both of these are rot about the X axis */
        case 3: R[0][0]=1;R[0][1]=0;R[0][2]=0;R[0][3]=0;
        case 2:/* Both of these are rot about the Y axis */
        case 4: 
      }  /*end switch*/
      mult(T,R,c); /* Combine trans & rotation */
      mult(c,Tinv,d);/* Combine that and inv trans */
      flip++;flip=flip%2;flop++;flop=flop%2; /* flip flop   */
      /* The following actually calculates new vert of rotat*/
       viewpts(flop);/* Undraw*/
       viewpts(flip);/* Draw*/
    }  /*end update points*/

*  End everything*
}  /*program end*/

void mult(A,B,C)
  float A[][4],B[][4],C[][4];
  int i,j,k;
}  /*end mult*/


Community Search:
MacTech Search:

Software Updates via MacUpdate

Warhammer 40,000: Freeblade guide - How...
Warhammer 40,000: Freebladejust launched in the App Store and it lets you live your childhood dream of blowing up and slashing a bunch of enemies as a massive, hulking Space Marine. It's not easy being a Space Marine though - and particularly if... | Read more »
Gopogo guide - How to bounce like the be...
Nitrome just launched a new game and, as to be expected, it's a lot of addictive fun. It's called Gopogo, and it challenges you to hoparound a bunch of platforms, avoiding enemies and picking up shiny stuff. It's not easy though - just like the... | Read more »
Sago Mini Superhero (Education)
Sago Mini Superhero 1.0 Device: iOS Universal Category: Education Price: $2.99, Version: 1.0 (iTunes) Description: KAPOW! Jack the rabbit bursts into the sky as the Sago Mini Superhero! Fly with Jack as he lifts impossible weights,... | Read more »
Star Wars: Galaxy of Heroes guide - How...
Star Wars: Galaxy of Heroes is all about collecting heroes, powering them up, and using them together to defeat your foes. It's pretty straightforward stuff for the most part, but increasing your characters' stats can be a bit confusing because it... | Read more »
The best cooking apps (just in time for...
It’s that time of year again, where you’ll be gathering around the dinner table with your family and a huge feast in front of you. [Read more] | Read more »
Square Rave guide - How to grab those te...
Square Rave is an awesome little music-oriented puzzle game that smacks of games like Lumines, but with its own unique sense of gameplay. To help wrap your head around the game, keep the following tips and tricks in mind. [Read more] | Read more »
Snowboard Party 2 (Games)
Snowboard Party 2 1.0 Device: iOS Universal Category: Games Price: $1.99, Version: 1.0 (iTunes) Description: Crowned the best snowboarding game available on the market, Snowboard Party is back to fulfill all your adrenaline needs in... | Read more »
One Button Travel (Games)
One Button Travel 1.0 Device: iOS Universal Category: Games Price: $2.99, Version: 1.0 (iTunes) Description: “To cut a long story short, If you like interactive fiction, just go buy this one.” - “Oozes the polish that... | Read more »
Light Apprentice Volume 1 (Games)
Light Apprentice Volume 1 1.0 Device: iOS Universal Category: Games Price: $4.99, Version: 1.0 (iTunes) Description: Light Apprentice Volume 1 includes Chapters 1 to 4, all gathered in a new exclusive game. When life in the world of... | Read more »
The best games like Animal Crossing on m...
Animal Crossing amiibo Festival is out right now for the Wii U, reminding us of just how much fun that world can be. Or at least to go back and check in on our villages once in a while. [Read more] | Read more »

Price Scanner via

iMobie Releases its Ace iOS Cleaner PhoneClea...
iMobie Inc. has announced the new update of PhoneClean 4, its iOS cleaner designed to reclaim wasted space on iPhone/iPad for use and keep the device fast. Alongside, iMobie hosts a 3-day giveaway of... Read more
U.S. Cellular Offering iPad Pro
U.S. Cellular today announced that it is offering the new iPad Pro with Wi-Fi + Cellular, featuring a 12.9-inch Retina display with 5.6 million pixels — the most ever in an iOS device. U.S. Cellular... Read more
Newegg Canada Unveils Black Friday Deals for...
Newegg Canada is offering more than 1,000 deep discounts to Canadian customers this Black Friday, available now through Cyber Monday, with new deals posted throughout the week. “Black Friday is... Read more
Black Friday: Macs on sale for up to $500 off...
BLACK FRIDAY B&H Photo has all new Macs on sale for up to $500 off MSRP as part of their early Black Friday sale including free shipping plus NY sales tax only: - 15″ 2.2GHz Retina MacBook Pro: $... Read more
Black Friday: Up to $125 off iPad Air 2s at B...
BLACK FRIDAY Walmart has the 16GB iPad Air 2 WiFi on sale for $100 off MSRP on their online store. Choose free shipping or free local store pickup (if available): - 16GB iPad Air 2 WiFi: $399, save $... Read more
Black Friday: iPad mini 4s on sale for $100 o...
BLACK FRIDAY Best Buy has iPad mini 4s on sale for $100 off MSRP on their online store for Black Friday. Choose free shipping or free local store pickup (if available): - 16GB iPad mini 4 WiFi: $299.... Read more
Black Friday: Apple Watch for up to $100 off...
BLACK FRIDAY Apple resellers are offering discounts and bundles with the purchase of an Apple Watch this Black Friday. Below is a roundup of the deals being offered by authorized Watch resellers:... Read more
Black Friday: Target offers 6th Generation iP...
BLACK FRIDAY Save $40 to $60 on a 6th generation iPod touch at Target with free shipping or free local store pickup (if available). Sale prices for online orders only, in-store prices may vary: -... Read more
Black Friday: Walmart and Target offer iPod n...
BLACK FRIDAY Walmart has the 16GB iPod nano (various colors) on sale for $119.20 on their online store for a limited time. That’s $30 off MSRP. Choose free shipping or free local store pickup (if... Read more
Black Friday deals on the Apple Watch and App...
Apple resellers are offering discounts and bundles with the purchase of an Apple Watch this Black Friday weekend. Below is a roundup of the deals being offered by authorized Watch resellers: Apple... Read more

Jobs Board

Storefront Operations Coordinator, *Apple* -...
# Storefront Operations Coordinator, Apple -Latin America Job Number: 43587750 Miami, Florida, United States Posted: Oct. 16, 2015 Weekly Hours: 40.00 **Job Summary** The Read more
*Apple* Enterprise / Government Professional...
# Apple Enterprise / Gove ment Professional Services Engineer Job Number: 42292976 Reston, Virginia, United States Posted: Aug. 18, 2015 Weekly Hours: 40.00 **Job Read more
iOS Wallet & *Apple* Pay Engineer - App...
# iOS Wallet & Apple Pay Engineer Job Number: 40586801 Santa Clara Valley, Califo ia, United States Posted: Nov. 16, 2015 Weekly Hours: 40.00 **Job Summary** The iOS Read more
Software Engineer, *Apple* Watch - Clock Fa...
# Software Engineer, Apple Watch - Clock Face Team Job Number: 44368761 Santa Clara Valley, Califo ia, United States Posted: Nov. 14, 2015 Weekly Hours: 40.00 **Job Read more
Administrative Assistant, *Apple* Online St...
# Administrative Assistant, Apple Online Store Job Number: 43992352 Santa Clara Valley, Califo ia, United States Posted: Nov. 9, 2015 Weekly Hours: 40.00 **Job Summary** Read more
All contents are Copyright 1984-2011 by Xplain Corporation. All rights reserved. Theme designed by Icreon.