Return to smoothing.c CVS log

Up to [local] / OpenXM_contrib2 / asir2018 / plot

Annotation of OpenXM_contrib2/asir2018/plot/smoothing.c, Revision 1.1

1.1     ! noro        1: /* $OpenXM$ */
        !             2: #include <stdio.h>
        !             3: #include <time.h>
        !             4: #include <math.h>
        !             5: #include "ca.h"
        !             6: #include "parse.h"
        !             7: #include "ox.h"
        !             8: #include "ifplot.h"
        !             9: #include "cursor.h"
        !            10:
        !            11: #define MAG 1
        !            12: #define PRINT_XOFFSET   100
        !            13: #define PRINT_YOFFSET   100
        !            14: #define DPI 300
        !            15:
        !            16: static FILE *Fp = NULL;
        !            17: static void PSProlog(int, int, int, int);
        !            18: static void PSDashedLine(int, int, int, int);
        !            19: static void polyLine_outputPS_dashed_line(int x0,int y0,int x1,int y1);
        !            20: void PSFromImage(FILE *fp,XImage *image,struct canvas *can);
        !            21:
        !            22: static int Strategy_generate_PS = 0;
        !            23:
        !            24: void PSFromImage(FILE *fp,XImage *image,struct canvas *can) {
        !            25:   int x, y, ir, ig, ib;
        !            26:   float r, g, b;
        !            27:   int srcc, trg,grux, gruy, ix, iy;
        !            28:   float xwd, ywd, t;
        !            29:
        !            30:   Fp = fp;
        !            31:   srcc = 0xffffff;
        !            32:   PSProlog(0,0,can->width,can->height);
        !            33:   for(y=0;y< can->height;y++)
        !            34:     for(x=0;x < can->width; x++){
        !            35:       trg = XGetPixel(image,x,y) & 0xffffff;
        !            36:       if ( trg != 0xffffff ){
        !            37:         if ( trg != srcc ){
        !            38:           srcc = trg;
        !            39:           r = ir = (srcc >>16) & 0xff;
        !            40:           r = r/256;
        !            41:           g = ig = (srcc >> 8) & 0xff;
        !            42:           g = g/256;
        !            43:           b = ib = (srcc) & 0xff;
        !            44:           b = b/256;
        !            45:           if(ir == 0 && ig ==0 && ib == 0)
        !            46:             fprintf(Fp,"0 0 0 C\n");
        !            47:           else
        !            48:             fprintf(Fp,"%f %f %f C\n",r,g,b);
        !            49:         }
        !            50:         fprintf(Fp,"%d %d R\n", x,can->height - y);
        !            51:       }
        !            52:     }
        !            53:   fprintf(Fp,"0 0 0 C\n");
        !            54:   if ( can->noaxis == 0 ){
        !            55:     grux = ceil(can->xmax);
        !            56:     xwd = can->width/(can->xmax - can->xmin);
        !            57:     for(ix=ceil(can->xmin); ix < grux; ix++){
        !            58:       t = ix;
        !            59:       x = (t - can->xmin) * xwd;
        !            60:       PSDashedLine(x, 0, x, can->height);
        !            61:     }
        !            62:     gruy = ceil(can->ymax);
        !            63:     ywd = can->height/(can->ymax - can->ymin);
        !            64:     for(iy=can->ymin; iy < gruy; iy++){
        !            65:       t = iy;
        !            66:       y = (t - can->ymin ) * ywd;
        !            67:       PSDashedLine(0, y, can->width, y);
        !            68:     }
        !            69:   }
        !            70:   fprintf(Fp,"[] 0 setdash 0 0 moveto ");
        !            71:   fprintf(Fp,"%d 0 lineto %d %d lineto 0 %d lineto closepath stroke\n",
        !            72:     can->width, can->width, can->height, can->height);
        !            73:   fprintf(Fp,"showpage grestore\n");
        !            74:   fflush(Fp);
        !            75: }
        !            76:
        !            77: static void PSProlog(int xmin, int ymin,int xmax, int ymax) {
        !            78:   int width, height;
        !            79:   struct tm *systime;
        !            80:   time_t t;
        !            81:
        !            82:   t = time(NULL);
        !            83:   systime = localtime(&t);
        !            84:   width = 72*(xmax - xmin-1)/300;
        !            85:   height = 72*(ymax - ymin-1)/300;
        !            86:   fprintf(Fp,"%%!PS-Adobe-2.0 EPSF-2.0\n");
        !            87:   fprintf(Fp,"%%%%BoundingBox: -3 -2 %d %d \n", width*MAG+3, height*MAG+3);
        !            88:   fprintf(Fp,"%%%%Creator: This is generated by ifplot\n");
        !            89:   fprintf(Fp,"%%%%Title: ifplot\n");
        !            90:   fprintf(Fp,"%%%%CreationDate: %.4d-%.2d-%.2d\n",
        !            91:     systime->tm_year+1900,systime->tm_mon+1,systime->tm_mday);
        !            92:   fprintf(Fp,"%%%%Pages: 0\n");
        !            93:   fprintf(Fp,"%%%%EndComments: \n");
        !            94:   fprintf(Fp,"%%%%EndProlog\n%%%%Page: one 1\n%%!\n");
        !            95:   fprintf(Fp,"/dpi %d def\n", DPI);
        !            96:   fprintf(Fp,"/SCALE 72 dpi div def\n");
        !            97:   fprintf(Fp,"/R {moveto 1 0 rlineto 0 1 rlineto -1 0 rlineto closepath");
        !            98:   fprintf(Fp," fill} bind def\n");
        !            99:   fprintf(Fp,"/C {setrgbcolor} bind def\n");
        !           100:   fprintf(Fp,"gsave newpath\n");
        !           101:   fprintf(Fp,"0 0 moveto SCALE SCALE scale 1 setlinewidth\n");
        !           102:   fflush(Fp);
        !           103: }
        !           104:
        !           105: static void PSDashedLine(int x0,int y0,int x1,int y1) {
        !           106:   extern FILE *Fp;
        !           107:   fprintf(Fp,"0.1 setlinewidth [3 3] 0 setdash\n");
        !           108:   fprintf(Fp," %d %d moveto %d %d lineto stroke \n",x0,y0,x1,y1);
        !           109: }
        !           110: