• Print |
• E-mail  |
• Bookmark

 /*+-------------------------------------------------------------+
   |          S A S   S A M P L E   L I B R A R Y                |
   |                                                             |
   |    NAME: GTGSISS                                            |
   |   TITLE: GTGSISS-Spline Interpolation and a Smoothed Spline |
   | PRODUCT: GRAPH                                              |
   |  SYSTEM: ALL                                                |
   |    KEYS: GRAPHICS GOPTIONS G3D G3GRID AXIS                  |
   |   PROCS: G3D G3GRID                                         |
   |    DATA: INTERNAL                                           |
   |                                                             |
   | SUPPORT: GRAPHICS STAFF              UPDATE:                |
   |     REF: SAS/GRAPH REFERENCE GUIDE                          |
   +-------------------------------------------------------------+*/

   /* Set the graphics environment */
goptions reset=global gunit=pct border cback=white
         colors=(blue black green red)
         ftext=swiss ftitle=swissb htitle=5 htext=3;

   /* Create data set WORK.NUMS using  */
   /* a set of randomly sampled points */
data nums;
   keep x y z;
   do i=1 to 30;
      x=10*ranuni(33)-5;
      y=10*ranuni(35)-5;
      z=sin(sqrt(x*x+y*y));
      output;
   end;
run;

   /* Define title and footnote for graph */
title 'Surface Plot After Spline Interpolation';
footnote j=r 'GTGSISS(a)';

   /* Process points with PROC G3GRID */
proc g3grid data=nums out=spline;
   grid y*x=z / spline
                axis1=-5 to 5 by .5
                axis2=-5 to 5 by .5;
run;

   /* Generate a surface plot */
proc g3d data=spline;
   plot y*x=z;
run;

   /* Define title and footnote for graph */
title 'Surface Plot After Smoothed Spline Interpolation';
footnote j=r 'GTGSISS(b)';

   /* Process points with PROC G3GRID */
proc g3grid data=nums out=smoothed;
   grid y*x=z / spline
                smooth=.05
                axis1=-5 to 5 by .5
                axis2=-5 to 5 by .5;
run;

   /* Generate a surface plot */
proc g3d data=smoothed;
   plot y*x=z;
run;
quit;

Operating System and Release Information