plplot绘制3d图

需要用到的plplot主要函数：

// Allocates a block of memory for use as a 2-d grid of PLFLT's.
void Alloc2dGrid( PLFLT ***f, PLINT nx, PLINT ny );
// Frees a block of memory allocated with plAlloc2dGrid().
void Free2dGrid( PLFLT **f, PLINT nx, PLINT ny );

// Plots a mesh representation of the function z[x][y].
    void mesh( const PLFLT *x, const PLFLT *y, const PLFLT * const *z, PLINT nx, PLINT ny, PLINT opt );
// Plots a mesh representation of the function z[x][y] with contour.
    void meshc( const PLFLT *x, const PLFLT *y, const PLFLT * const *z, PLINT nx, PLINT ny, PLINT opt,
                const PLFLT *clevel, PLINT nlevel );

程序源码如下：

    #include "plc++demos.h"

    #ifdef PL_USE_NAMESPACE
    using namespace std;
    #endif
    plstream           *pls;
    const int   x_PointNum   = 35;
    const int   y_PointNum   = 46;
    const int   LEVELS = 10;
    const int   opt[] = { 3, 3 };
    const PLFLT alt[] = { 33.0, 17.0 };
    const PLFLT az[] = { 24.0, 115.0 };

    const char *title[] = {
       "#frPLplot Example 11 - Alt=33, Az=24, Opt=3",
       "#frPLplot Example 11 - Alt=17, Az=115, Opt=3"
    };
    void cmap1_init()
    {
       PLFLT *i   = new PLFLT[2];
       PLFLT *h   = new PLFLT[2];
       PLFLT *l   = new PLFLT[2];
       PLFLT *s   = new PLFLT[2];
       bool *rev = new bool[2];

       i[0] = 0.0;       // left boundary
       i[1] = 1.0;       // right boundary

       h[0] = 240;       // blue -> green -> yellow ->
       h[1] = 0;         // -> red

       l[0] = 0.6;
       l[1] = 0.6;

       s[0] = 0.8;
       s[1] = 0.8;

       rev[0] = false;       // interpolate on front side of colour wheel.
       rev[1] = false;       // interpolate on front side of colour wheel.

       pls->scmap1n( 256 );// Set number of colors in cmap 1
       pls->scmap1l( false, 2, i, h, l, s, rev );

       delete[] i;
       delete[] h;
       delete[] l;
       delete[] s;
       delete[] rev;
    }

    // Does a series of mesh plots for a given data set, with different viewing
    // options in each plot.
    int main( int argc, const char **argv )
    {
       int   i, j, k;

       PLFLT *x = new PLFLT[ x_PointNum ];
       PLFLT *y = new PLFLT[ y_PointNum ];
       PLFLT **z;
       PLFLT zmin = 1E9, zmax = -1E-9;

       PLFLT xx, yy;
       int   nlevel = LEVELS;
       PLFLT *clevel = new PLFLT[LEVELS];
       PLFLT step;
       pls = new plstream();
       // Parse and process command line arguments.
       pls->parseopts( &argc, argv, PL_PARSE_FULL );
       // Initialize plplot.
       pls->sdev("qtwidget");
       //pls->init();
       pls->star(4,2);
       // Allocates a block of memory for use as a 2-d grid of PLFLT's.
       pls->Alloc2dGrid( &z, x_PointNum, y_PointNum );

       for ( i = 0; i < x_PointNum; i++ )
           x[i] = 3. * (PLFLT) ( i - ( x_PointNum / 2 ) ) / (PLFLT) ( x_PointNum / 2 );

       for ( j = 0; j < y_PointNum; j++ )
           y[j] = 3. * (PLFLT) ( j - ( y_PointNum / 2 ) ) / (PLFLT) ( y_PointNum / 2 );

       for ( i = 0; i < x_PointNum; i++ )
       {
           xx = x[i];
           for ( j = 0; j < y_PointNum; j++ )
           {
               yy      = y[j];
               z[i][j] = 3. * ( 1. - xx ) * ( 1. - xx ) * exp( -( xx * xx ) - ( yy + 1. ) * ( yy + 1. ) ) -
                   10. * ( xx / 5. - pow( (double) xx, 3. ) - pow( (double) yy, 5. ) ) * exp( -xx * xx - yy * yy ) -
                   1. / 3. * exp( -( xx + 1 ) * ( xx + 1 ) - ( yy * yy ) );
               if ( false ) // Jungfraujoch/Interlaken
               {
                   if ( z[i][j] < -1. )
                       z[i][j] = -1.;
               }
               if ( zmin > z[i][j] )
                   zmin = z[i][j];
               if ( zmax < z[i][j] )
                   zmax = z[i][j];
           }
       }

       step = ( zmax - zmin ) / ( nlevel + 1 );
       for ( i = 0; i < nlevel; i++ )
           clevel[i] = zmin + step + step * i;

       cmap1_init();
       for ( k = 0; k < 2; k++ )
       {
           for ( i = 0; i < 4; i++ )
           {
               pls->adv( 0 );
               pls->col(Red);
               pls->vpor( 0.0, 1.0, 0.0, 0.9 );
               pls->wind( -1.0, 1.0, -1.0, 1.5 );
               //尺寸范围
               pls->w3d( 1.0, 1.0, 1.2, -3.0, 3.0, -3.0, 3.0, zmin, zmax,
                   alt[k], az[k] );
               //box set
               pls->box3( "bnstu", "x axis", 0.0, 0,
                   "bnstu", "y axis", 0.0, 0,
                   "bcdmnstuv", "z axis", 0.0, 4 );

               pls->col(Yellow);
               // wireframe plot
               if ( i == 0 )
                   pls->mesh( x, y, z, x_PointNum, y_PointNum, opt[k] );//单色显示

               // magnitude colored wireframe plot
               else if ( i == 1 )
                   pls->mesh( x, y, z, x_PointNum, y_PointNum, opt[k] | MAG_COLOR );//等高颜色渐变显示

               // magnitude colored wireframe plot with sides
               else if ( i == 2 )
               {
                   pls->plot3d( x, y, z, x_PointNum, y_PointNum, opt[k] | MAG_COLOR ,true );
               }

               // magnitude colored wireframe plot with base contour
               else if ( i == 3 )
                   pls->meshc( x, y, z, x_PointNum, y_PointNum, opt[k] | MAG_COLOR | BASE_CONT,
                   clevel, nlevel );

               pls->col(Green);
               pls->mtex( "t", 1.0, 0.5, 0.5, title[k] );
           }
       }

       // Frees a block of memory allocated with plAlloc2dGrid().
       pls->Free2dGrid( z, x_PointNum, y_PointNum );

       delete[] x;
       delete[] y;
       delete[] clevel;

       delete pls;
       return 0;
    }

vs2010工程地址：plplot_3d

参考文献：http://plplot.sourceforge.net/examples.php?demo=11&lbind=C

plplot绘制3d图

相关推荐