Data visualization with SAS programming

Fun With SAS ODS Graphics, July 4th Edition

Reply
Contributor tc
Contributor
Posts: 67

Fun With SAS ODS Graphics, July 4th Edition

fourthofjuly.gif

 

OK, animated GIF GTL vector plots admittedly are not nearly as fun as the firecrackers of my youth, but at least I'll still have all my fingers at the end of the day!

 

Acknowledgements: If some of the code/concepts look familiar to readers of Sanjay Matange's Graphically Speaking blog posts on some real use cases for vector plots and animated plots, it's no coincidence (bugs/kludges are my own!).

 

SAS CODE

* Fun with GTL vectorplots, 4th of July animated GIF "fireworks";

data fireworks;                * Generate points for fireworks;
pi=constant('pi');
do firework=1 to 7;            * Seven fireworks "explosions";
  r=.2+ranuni(2)*.3;           * Random radius between .2 and .5;
  if firework=7 then r=.5;     * Big finish! (radius=.5);
  xo=r+ranuni(4)*(1-r*2);      * Random origin (x, y);
  yo=r+ranuni(6)*(1-r*2);
  do a=0 to 352.5 by 7.5;      * Line every 7.5 degrees;
    x=r*cos(a*pi/180)+xo;      * Calc x and y coordinates for fully-exploded fireworks;
    y=r*sin(a*pi/180)+yo;
    xi1=.01*cos(a*pi/180)+xo;  * Initial unexploded x, y points (r=.01, .02);
    yi1=.01*sin(a*pi/180)+yo;
    xi2=.02*cos(a*pi/180)+xo;
    yi2=.02*sin(a*pi/180)+yo;
    output;
  end;
end;

data empty;                    * Generate out-of-frame points to create "empty" frame;
x=-1; y=-1; xo=-1; yo=-1; r=1;
run;

ods _all_ close;               * Use generated points to produce animated GIF;
options papersize=('6 in', '6 in') printerpath=gif animation=start 
        nodate nonumber animduration=.1 animloop=YES NOANIMOVERLAY;
ods printer file='/folders/myfolders/fireworks/fourthofjuly.gif';
ods graphics / border=off width=6in height=6in imagefmt=GIF;

%macro fireworks;              * Vector plot (change pattern to "fade" fireworks;
%macro modtemplate(pattern=);
proc template;
define statgraph vectorplot;
begingraph / opaque=true border=false backgroundcolor=black;
  layout overlayequated / equatetype=square opaque=true border=false backgroundcolor=black 
         wallcolor=black WALLDISPLAY=(FILL) commonaxisopts=(viewmin=0 viewmax=1) 
         xaxisopts=(display=none) yaxisopts=(display=none);
    vectorplot y=y x=x xorigin=xo yorigin=yo / 
               arrowheads=false lineattrs=(pattern=&pattern color=&linecolor thickness=1.6pt);
  endlayout;
endgraph;
end;
run;
%mend;

%do f=1 %to 7;                  * Generate frames for animated GIF;
   data _null_;
   call symput("linecolor",scan("hotpink orange dodgerblue yellow lime white red", &f)); 
   run;
   %modtemplate(pattern=solid); * Initial unexploded frames (2 frames);
   proc sgrender data=fireworks(where=(firework=&f) drop=x y rename=(xi1=x yi1=y)) template=vectorplot;
   proc sgrender data=fireworks(where=(firework=&f) drop=x y rename=(xi2=x yi2=y)) template=vectorplot;
   %modtemplate(pattern=solid);
   %do i=1 %to 5;               * Fully exploded frames (5 frames);
     proc sgrender data=fireworks(where=(firework=&f)) template=vectorplot;
   %end;                        * Fading/empty frames (3 frames);
   %modtemplate(pattern=shortdashdot);
   proc sgrender data=fireworks(where=(firework=&f)) template=vectorplot;
   %modtemplate(pattern=dot);
   proc sgrender data=fireworks(where=(firework=&f)) template=vectorplot;
   proc sgrender data=empty template=vectorplot;
%end;
run; 
%mend;

%fireworks;                     * Light 'em up!; 
options printerpath=gif animation=stop;
ods printer close;

 

LINE PATTERNS USED TO SHOW/"FADE" FIREWORKS (SOLID, SHORTDASHDOT, DOT)

Fireworks123.gif

SAS Super FREQ
Posts: 3,233

Re: Fun With SAS ODS Graphics, July 4th Edition

I like it! A fun application of animated plots. For another fun example, see The Game of Life animation.

 

One minor comment on the SAS program. Only the first random number seed is relevant in your DATA step program. Changing the seed for each call doesn't do anything.

Contributor tc
Contributor
Posts: 67

Re: Fun With SAS ODS Graphics, July 4th Edition

Thanks for the "random tip"!

 

Up until now, I think my RANUNI needs have been one-statement-per-DATA-step simple.

 

Like The Game of Life in SAS/IML Studio, too!

SAS Super FREQ
Posts: 1,044

Re: Fun With SAS ODS Graphics, July 4th Edition

Awesome July 4 firework animation, Ted.  Love the way you fade out the light using line patterns.  Made my day.

SAS Employee
Posts: 963

Re: Fun With SAS ODS Graphics, July 4th Edition

That's pretty cool! Smiley Happy

Valued Guide
Posts: 505

Re: Fun With SAS ODS Graphics, July 4th Edition

WOW!!!

 

I had ro comment out  opaque=true   'begingraph /   /*opaque=true*/  border=false backgroundcolor=black; ".

Ran nicely after that

 

Win 7 64 SAS 64 9.4M2

SAS Employee
Posts: 24

Re: Fun With SAS ODS Graphics, July 4th Edition

Breathtaking!

 

May I suggest a tiny tweak? If you were to use mvar or dynamics for the vectorplot statement's pattern= and linecolor=, you can do without the macro modtemplate(). Tweaked code attached.

Attachment
Contributor tc
Contributor
Posts: 67

Take 2: Static B/W image of a GTL vectorplot "fireworks" over downtown Chicago skyline

Chicago4thJuly.gif

And with the help of an adapted version of a Wikimedia Commons image attributed to Eduardo Manchon and a tip from an old SGF Paper of Dan Heath's, here's a simpler version of the code that generates a static B/W image of a GTL vectorplot fireworks show over the downtown Chicago skyline!

 

CODE

* Chicago skyline with GTL vectorplot "fireworks";
 
data fireworks;                * Generate points for fireworks;
pi=constant('pi');
input r xo yo;                 * Read radius, x/y origin for fireworks;
firework+1;                    * Assign group ID for later vector plot;
do a=0 to 352.5 by 7.5;        * Line every 7.5 degrees;
  x=r*cos(a*pi/180)+xo;        * Calc x and y coordinates for fireworks;
  y=r*sin(a*pi/180)+yo;
  output;
end;
datalines;
.15 .13 .85
.23 .62 .775
.07 .95 .925
;
ods listing gpath='/folders/myfolders/fireworks';
ods graphics on / reset=index imagename='Chicago4thJuly' border=off
                  width=5in height=5in imagefmt=GIF;

data chicago;                  * Image of Chicago skyline;
retain function "Image" anchor "bottomleft" x1 0 y1 0 width 100 height 100
       widthunit 'percent' heightunit 'percent' imagescale  'fit'
       drawspace "layoutpercent" layer 'back' 
       id 'chicago' Image "/folders/myfolders/chicago.gif";

proc template;
define statgraph vectorplot;
begingraph / opaque=true border=false drawspace=layoutpercent backgroundcolor=black;
  layout overlayequated / equatetype=square border=false  
         WALLDISPLAY=NONE commonaxisopts=(viewmin=0 viewmax=1) 
         xaxisopts=(display=none) yaxisopts=(display=none);
    vectorplot y=y x=x xorigin=xo yorigin=yo / group=firework  
               arrowheads=false lineattrs=(pattern=solid color=white);
    annotate / id="chicago";           
  endlayout;
endgraph;
end;

proc sgrender data=fireworks template=vectorplot sganno=chicago;

BACKGROUND IMAGE

chicago.gif

 

 

 

 

 

 

Super User
Posts: 1,208

Re: Take 2: Static B/W image of a GTL vectorplot "fireworks" over downtown Chicago skyline

Having a 'blast' with your SAS graphic creativity Ted! Great work. 

Post a Question
Discussion Stats
  • 8 replies
  • 989 views
  • 7 likes
  • 7 in conversation