/* mat2cps.c */
#include <stdio.h>
#include <ctype.h>
#include <math.h>
#define BOOLEAN char
#define FALSE 0
#define VERS 1.3
main(argc,argv)
int argc;
char *argv[];
{
 FILE *fd_inp;
 FILE *fd_out;
 char name[20];
 int type, imagf;
 double *inpr, *inpi; /* matrix real and imaginary: imaginary not used */
/****** double sqrt(); UNDEFINED _sqrt... ******/
 int M, N, MN; /* size of original matrix */
 int m, n, mn;
 int msd, lsd; /* least most significant digits to output in hex */
 unsigned char inp_file[256]; /* char *filename gives 4 characters/byte */
 unsigned char out_file[256];
 /* HELP MESSAGE */
 if (argc < 2) { /* give help when only program name given (no args) */
 fprintf(stderr,"%s: convert matlab matrix to eps\n",
 argv[0]);
 fprintf(stderr,"version %f Steve Mann; 1990\n",VERS);
 fprintf(stderr,"---------------------------------------------------");
 fprintf(stderr,"----------------------------\n");
 fprintf(stderr,"%s deleteme.mat chirplet_growler.eps ",argv[0]);
 fprintf(stderr,"converts deleteme.mat to eps file\n");
 fprintf(stderr,"%s deleteme.mat &#124; lpr -Pps ",argv[0]);
 fprintf(stderr,"stdout (without creating a file)\n");
 fprintf(stderr,"---------------------------------------------------");
 fprintf(stderr,"----------------------------\n");
 exit(0);
 }
 strcpy(inp_file, argv[1]); /* for now input file is always first */
 fd_inp = fopen(inp_file,"rb");
 if (fd_inp==NULL) {perror(argv[0]); exit(-1);}
 /* read in input matrix */
 if (loadmat(fd_inp, &type, name, &M, &N, &imagf, &inpr, &inpi)) {
		printf("\nRead error when calling loadmat\n");
 }
 fclose(fd_inp);
 fprintf(stderr,"%s: input matrix has been read in\n",argv[0]);
 /* output stream may be stdout or file */
 if (argc < 3) { /* only one extra command line arg given */
 fd_out = stdout;
 fprintf(stderr,"%s : output is stdout\n",argv[0]);
 }
 else { /* output is a file */
 strcpy(out_file, argv[argc-1]); /* last argument is outfile */
 fprintf(stderr,"%s ; output file = %s\n", argv[0], out_file);
 fd_out = fopen(out_file, "w");
 if (fd_out==NULL) {perror(argv[0]); exit(-1);}
 }
 /* write out PostScript header */
 fprintf(fd_out, "%%!PS-Adobe-2.0 EPSF-1.2\n");
 fprintf(fd_out, "%%%%Copyright (c) 1991 by SwimWhere SoftWear\n");
 fprintf(fd_out, "%%%%Creator: mat2eps_tf.c version 1.3, Steve Mann\n");
 fprintf(fd_out, "%%%%Title: time-freq image\n");
 fprintf(fd_out, "%%%%Pages: 1\n");
 fprintf(fd_out, "%%%%BoundingBox: 100 90 495 504\n");
 fprintf(fd_out, "%%%%EndComments\n");
 fprintf(fd_out, "\n");
 fprintf(fd_out, "%% tf_eps dictionary\n");
 fprintf(fd_out, "/tfeps 20 dict begin\n");
 fprintf(fd_out, "%% definition operators\n");
 fprintf(fd_out, "/bdef {bind def} bind def\n");
 fprintf(fd_out, "/xdef {exch def} bdef\n");
 fprintf(fd_out, "%% page state control\n");
 fprintf(fd_out, "/pgsv () def\n");
 fprintf(fd_out, "/bpage {/pgsv save def} bdef\n");
 fprintf(fd_out, "/epage {pgsv restore} bdef\n");
 fprintf(fd_out, "/picstr %d string def\n",M);
 fprintf(fd_out, "/imagefile\n");
 if (M!=N)
 fprintf(stderr,"%s ; still not sure which is M and which is N\n",argv[0]);
 fprintf(fd_out, "{ %d %d 8 [%d 0 0 %d 0 0]\n",M,N,M,N);
 fprintf(fd_out, " {currentfile picstr readhexstring pop} image} def\n");
 fprintf(fd_out, "/xrightshow %% stk: string\n");
 fprintf(fd_out, "{dup stringwidth pop\n");
 fprintf(fd_out, "neg 0 rmoveto %% move left\n");
 fprintf(fd_out, "show} def\n");
 fprintf(fd_out, "/xcentershow %% stk: string\n");
 fprintf(fd_out, "{dup stringwidth pop\n");
 fprintf(fd_out, "2 div %% half width\n");
 fprintf(fd_out, "neg 0 rmoveto %% move left\n");
 fprintf(fd_out, "show} def\n");
 fprintf(fd_out, "currentdict end def %% dictionary\n");
 fprintf(fd_out, "%%%%EndProlog\n");
 fprintf(fd_out, "\n");
 fprintf(fd_out, "%%%%Page: 1 1\n");
 fprintf(fd_out, "tfeps begin\n");
 fprintf(fd_out, "bpage\n");
 fprintf(fd_out, "gsave\n");
 fprintf(fd_out, "360 360 scale\n");
 fprintf(fd_out, "0.35 0.35 translate\n");
 fprintf(fd_out, "/Times-Roman findfont\n");
 fprintf(fd_out, "0.05 scalefont\n");
 fprintf(fd_out, "setfont\n");
 fprintf(fd_out, "0.5 1.1 moveto\n");
 fprintf(fd_out, "(delete this title Title TITLE or change 504 to 534) xcentershow\n");
 fprintf(fd_out, "%% image data\n");
 fprintf(fd_out, "imagefile\n");
 fprintf(stderr,"%s ; written out first part of header\n", argv[0]);
 
 MN = M*N;
 for (mn=0; mn<mn; mn++) { /* put into one long column */
 msd = (unsigned char)inpr[mn]/16;
 lsd = (unsigned char)inpr[mn]%16;
 fprintf(fd_out, "%x%x",msd,lsd);
 if (mn%32 == 31) fprintf(fd_out, "\n"); /* new line for each 32 bytes */
 }
 fprintf(fd_out, "\n"); /* new line at the very end in case not mult of 32 */
 /* tested that it is OK to have extra blank lines */
 fprintf(stderr,"%s ; written out image data\n", argv[0]);
 fprintf(fd_out, "0.003 setlinewidth\n");
 fprintf(fd_out, "%% box around image\n");
 fprintf(fd_out, "0.0 0.0 moveto\n");
 fprintf(fd_out, "0.0 1.0 lineto\n");
 fprintf(fd_out, "1.0 1.0 lineto\n");
 fprintf(fd_out, "1.0 0.0 lineto\n");
 fprintf(fd_out, "0.0 0.0 lineto\n");
 fprintf(fd_out, "0.0 0.0 moveto\n");
 fprintf(fd_out, "%% up and down tick marks for axis to the right\n");
 fprintf(fd_out, "0.0 0.0 moveto\n");
 fprintf(fd_out, "0.0 -0.02 lineto\n");
 fprintf(fd_out, "1.0 0.0 moveto\n");
 fprintf(fd_out, "1.0 -0.02 lineto\n");
 fprintf(fd_out, "%% left to right ticks for axis up\n");
 fprintf(fd_out, "0.0 0.0 moveto\n");
 fprintf(fd_out, "-0.02 0.0 lineto\n");
 fprintf(fd_out, "0.0 0.5 moveto\n");
 fprintf(fd_out, "-0.02 0.5 lineto\n");
 fprintf(fd_out, "0.0 1.0 moveto\n");
 fprintf(fd_out, "-0.02 1.0 lineto\n");
 fprintf(fd_out, "%% tick marks to right and on top, going up\n");
 fprintf(fd_out, "1 .5 moveto\n");
 fprintf(fd_out, "1.02 .5 lineto\n");
 fprintf(fd_out, "1 0 moveto\n");
 fprintf(fd_out, "1.02 0 lineto\n");
 fprintf(fd_out, "1 1 moveto\n");
 fprintf(fd_out, "1.02 1 lineto\n");
 fprintf(fd_out, "0 1 moveto\n");
 fprintf(fd_out, "0 1.02 lineto\n");
 fprintf(fd_out, "1 1 moveto\n");
 fprintf(fd_out, "1 1.02 lineto\n");
 fprintf(fd_out, "%% centerline for X axis y=0\n");
 fprintf(fd_out, "0 .5 moveto\n");
 fprintf(fd_out, "1 .5 lineto\n");
 fprintf(fd_out, "%% centerline for Y axis x=0\n");
 fprintf(fd_out, ".5 0 moveto\n");
 fprintf(fd_out, ".5 1 lineto\n");
 fprintf(fd_out, "stroke\n");
 fprintf(fd_out, "/Times-Roman findfont\n");
 fprintf(fd_out, "0.04 scalefont\n");
 fprintf(fd_out, "setfont\n");
 fprintf(fd_out, "0.5 1.14 moveto\n");
 fprintf(fd_out, "0.5 1.09 moveto\n");
 fprintf(fd_out, "0.5 -0.06 moveto\n");
 fprintf(fd_out, "(physical ``X'' dimension) xcentershow\n");
 fprintf(fd_out, "gsave\n");
 fprintf(fd_out, " 90 rotate\n");
 fprintf(fd_out, " 0.5 .03 moveto\n");
 fprintf(fd_out, " (physical ``Y'' dimension) xcentershow\n");
 fprintf(fd_out, "grestore\n");
 fprintf(fd_out, "/Times-Roman findfont\n");
 fprintf(fd_out, "0.03 scalefont\n");
 fprintf(fd_out, "setfont\n");
 fprintf(fd_out, "gsave\n");
 fprintf(fd_out, " 90 rotate\n");
 fprintf(fd_out, " 0.0 0.03 moveto\n");
 fprintf(fd_out, " (-) show\n");
 fprintf(fd_out, " 1.0 0.03 moveto\n");
 fprintf(fd_out, " (+) xrightshow\n");
 fprintf(fd_out, "grestore\n");
 fprintf(fd_out, "0.0 -0.06 moveto\n");
 fprintf(fd_out, "(-) show\n");
 fprintf(fd_out, "1.0 -0.06 moveto\n");
 fprintf(fd_out, "(+) xrightshow\n");
 fprintf(fd_out, "grestore\n");
 fprintf(fd_out, "epage\n");
 fprintf(fd_out, "end\n");
 fprintf(fd_out, "\n");
 fprintf(fd_out, "%%%%Trailer\n");
 fprintf(fd_out, "\n");
 fprintf(fd_out, "showpage\n");
 fprintf(stderr,"%s done: converted %d by %d (%d long) mat to %d/32 lines of hex text\n"
 ,argv[0], M, N, M*N, M*N);
}
</div><div class="naked_ctrl">
<form action="/index.cgi/contrast" method="get" name="gate">
<p><a href="http://altstyle.alfasado.net">AltStyle</a> によって変換されたページ <a href="http://wearcam.org/orbits/matlab/mat2eps_xy.c">(-&gt;オリジナル)</a>
/ <label>アドレス: <input type="text" name="naked_post_url" value="http://wearcam.org/orbits/matlab/mat2eps_xy.c" size="22" /></label> <label>モード: <select name="naked_post_mode">
<option value="default">デフォルト</option>
<option value="speech">音声ブラウザ</option>
<option value="ruby">ルビ付き</option>
<option value="contrast" selected="selected">配色反転</option>
<option value="larger-text">文字拡大</option>
<option value="mobile">モバイル</option>
</select>
<input type="submit" value="表示" />
</p>
</form>
</div>