// radio array, 5x5x5, random spacing // // compile with cc -o a3d05 a3d05.c -lgd -lpng -lm // // Uses the libgd library. For documentation see the file: // /usr/share/doc/gd-*/index.html // or the website: http://www.libgd.org/Reference // // You will need truetype fonts. If you don't have them, you can // copy ../fonts/truetype/.. from openoffice.org to /usr/share/ // // Uses swftools to build a swf movie from individual png images // for more information, see http://www.swftools.org/ // // This constructs half the slides, then mirrors them #include "gd.h" #include "math.h" #include #include #include #define RMMKDIR "rm -rf a3d05dir ; mkdir a3d05dir" #define PNG2SWF "png2swf -o a3d05.swf -r 10 a3d05dir/*.png" #define PNGFMT "a3d05dir/a%04d.png" #define YCENT 130 #define WAVELENGTH 10.0 // radio wavelength in pixel units #define SPACE 40.0 // array spacing in pixel units #define NPICS 50 // ( half the ) number of frames #define RSPACE 30.0 // random variation #define TANGX 0.0 // target X angle in degrees #define TANGY 0.0 // target Y angle in degrees #define DANGX 0.0 // display center X angle in degrees #define DANGY 0.0 // display center Y angle in degrees #define DANGW 150.0 // display angle width in degrees #define DANGS 10.0 // display angle step #define ZSCALE 400.0 // depth scaling // #define TILT 0.04 // array tilt around X axis, radians #define TILT 0.00 // array tilt around X axis, radians #define SPOT 10 // size of drawn spot #define ASIZE 5 // number of elements per side #define XLEFT 40 // left side pixel count #define XSIZE 1000 // display window in pixels #define YSIZE 750 // display window in pixels #define FNT "DejaVuMonoSans" #define FSZ 10 // ========================================================================== int main() { gdImagePtr im ; // pixel image map used by gd double space = (double) SPACE ; // nominal spacing between server-sats double xleft = (double) XLEFT ; // left of pixel array double yacent = (double) YCENT ; // center of upper drawing double xsize = (double) XSIZE ; // width of pixels double ysize = (double) YSIZE ; // height of pixels double dangw = (double) DANGW ; // degree width of output pixel array double dangx = (double) DANGX ; // degree center of output pixel array double dangy = (double) DANGY ; // degree center of output pixel array double dangw2 = 0.5*dangw ; // half width of output pixel array double xwide = xsize-xleft ; // pixel width of output pixel array int ybot = 2*YCENT ; // top of output pixel array double dybot = (double) ybot ; // real "" double ytall = ysize-dybot ; // height of output pixel array double pi2 = 8.0*atan(1.0) ; // two times PI double d2r = pi2/360.0 ; // degrees to radians conversion factor int oxc =(XSIZE+XLEFT)/2 ; // center of output pixel array int oyc =(YSIZE+ybot)/2 ; // center of output pixel array double xacent = (double) oxc ; // real "" double aspect =ytall/xwide ; // aspect ration of output pixel array double angH =dangw*aspect ; // angular height of output pixel array double angH2 = 0.5*angH ; // half height of output pixel array double angx ; // pixel X angle radians double angx0=d2r*(dangx-dangw2); // left pixel value, radians double angxs=d2r*dangw/xwide ; // pixel X slope radians per pixel double angxf=angx0-angxs*xleft ; // pixel X offset radians double angy ; // pixel Y angle double angy0=d2r*(dangy+angH2) ; // top pixel value, radians double angys=d2r*angH/ytall ; // pixel Y slope radians double angyf=angy0-angys*ysize ; // pixel Y offset radians double k = pi2/WAVELENGTH ; // wavenumber double rspace= RSPACE ; // random spacing amount double tilt = TILT ; // tilt of array (around Y) in radians FILE *pngout ; // file handle for PNG output frame char dirname[80] ; // frame output directory char framename[80] ; // directory/name of frame char labstring[80] ; // used for labellin char link0[80] ; // directory/name of source char link1[80] ; // directory/name of target hardlink int black, sun1, white, red ; // color map numbers int green, blue, gray ; // color map numbers int dgray, trans ; // color map numbers int gcolor[256] ; // color map numbers double ax[ASIZE][ASIZE][ASIZE] ; // x position of element double ay[ASIZE][ASIZE][ASIZE] ; // y position of element double az[ASIZE][ASIZE][ASIZE] ; // z position of element double rx[ASIZE][ASIZE][ASIZE] ; // x variation of element double ry[ASIZE][ASIZE][ASIZE] ; // y variation of element double rz[ASIZE][ASIZE][ASIZE] ; // z variation of element double ap[ASIZE][ASIZE][ASIZE] ; // phase of element double xx, yy, zz ; // double xr, yr, zr ; // double bx, by, bz ; // coordinates rotated on y axis double cx, cy, cz ; // coordinates rotated on x axis int ix, iy, iz ; // integer counters int ox, oy ; // output pixel position int i ; // general counter double zscale ; // scaling for perspective depth int icolor ; // int cir ; // double stilt = sin(tilt) ; // sin of rotation (around y) of array double ctilt = cos(tilt) ; // cos of rotation (around y) of array int mag ; // pixel magnitude, 0-255 double ss, cc ; // summed sin and cos components double stx = sin(TANGX ) ; // sin of X target angle double ctx = cos(TANGX ) ; // cos of X target angle double sty = sin(TANGY ) ; // sin of Y target angle double cty = cos(TANGY ) ; // cos of Y target angle double spx ; // sin of pixel X angle double cpx ; // cos of pixel X angle double spy ; // sin of pixel Y angle double cpy ; // cos of pixel Y angle double phase ; // phase int xm ; // tickmark max int xs ; // tickmark scale int frame ; // animation frame count int pct=0 ; // percentage of complete frame int pct0=0 ; // previous percentage double rfactor=0.0 ; // random amount added to position int asize = ASIZE ; // nominal spacing of server-sats double offset=0.5*(1.0-asize) ; // offsets counter integer to center double asize3=asize*asize*asize; // double norm = 255.0 /asize3 ; // normalizes summed amplitude #ifdef DEBUG printf( "%12.6f%12.6f%15.9f\n", angx0, angxf, angxs ); printf( "%12.6f%12.6f%15.9f\n", angy0, angyf, angys ); #endif // set up target directory system( RMMKDIR ); // random value generation ---------------------------------------------- for( iz=0; iz 255 ) { icolor = 255 ; } gdImageFilledEllipse(im, ox, oy, cir, cir, gcolor[icolor] ); } } } #ifdef SINGLE printf( "array drawn\n"); #endif // now compute phases for( ox=XLEFT ; ox