int main(int argc, char *argv[])
{
RGBPic *p;
PLOT *plot;
p=make_RGBPic(800,600);
RGBPic_clear_area(p,COLOR(60,80,50),0,0,p->width,p->height);
plot=make_plot(640,480);
plot->lower_x=3.1415;
plot->upper_x=12423;
plot->lower_y=2.34E-10;
plot->upper_y=43.5E-8;
draw_grid(p,plot,50,50);
free_plot(plot);
//RGBPic_printf(p,30,30,COLOR(255,240,210),COLOR(1,2,3),"Hello all !\n line 2 %d %p",argc,p);
//RGBPic_printf_v(p,300,300,COLOR(255,240,210),COLOR(1,2,3),"Hello all 2 !\n line 2 %d %p",argc,p);
RGBPic_draw_line(p,COLOR(40,255,80), 40,50,60,70);
RGBPic_draw_line(p,COLOR(250,255,80), 30,50,60,30);
RGBPic_dump_ppm("test1.ppm",p);
free_RGBPic(p);
return 0;
}
int main(int argc, char *argv[])
{
FILE *fin;
int i;
float a;
LOG=stderr;
if(argc<7){
fprintf(stderr, "\nUsage: %s mode resolution lower_cut upper_cut file.dat skymap.png\n", argv[0]);
fprintf(stderr, "\tmodes: plain, log10\n");
fprintf(stderr, "\tcutoff: number or \"auto\"\n\n");
return -1;
}
resolution=atof(argv[2]);
fin=fopen(argv[5], "r");
if(fin==NULL){
perror("Cannot read file:");
return -1;
}
grid=make_sin_theta_grid(resolution);
data=do_alloc(grid->npoints, sizeof(*data));
fread(data, sizeof(*data), grid->npoints, fin);
fprintf(stderr, "%g %g %g %g\n", data[0], data[1], data[2], data[3]);
if(!strcasecmp(argv[1], "log10")){
for(i=0;i<grid->npoints;i++)data[i]=log10(data[i]);
}
if(strcasecmp(argv[3], "auto")){
a=atof(argv[3]);
for(i=0;i<grid->npoints;i++)
if(data[i]<=a)data[i]=a;
}
if(strcasecmp(argv[4], "auto")){
a=atof(argv[4]);
for(i=0;i<grid->npoints;i++)
if(data[i]>=a)data[i]=a;
}
pic=make_RGBPic(grid->max_n_ra+140, grid->max_n_dec);
plot_grid_f(pic, grid, data, 1);
RGBPic_dump_png(argv[6], pic);
return 0;
}
