WIP: new image_header struct is used and enable used of cstr_info
[openjpeg.git] / applications / codec / j2k_dump.c
1 /*
2  * Copyright (c) 20010, Mathieu Malaterre, GDCM
3  * All rights reserved.
4  *
5  * Redistribution and use in source and binary forms, with or without
6  * modification, are permitted provided that the following conditions
7  * are met:
8  * 1. Redistributions of source code must retain the above copyright
9  *    notice, this list of conditions and the following disclaimer.
10  * 2. Redistributions in binary form must reproduce the above copyright
11  *    notice, this list of conditions and the following disclaimer in the
12  *    documentation and/or other materials provided with the distribution.
13  *
14  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
15  * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
16  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
17  * ARE DISCLAIMED.  IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
18  * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
19  * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
20  * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
21  * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
22  * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
23  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
24  * POSSIBILITY OF SUCH DAMAGE.
25  */
26
27 #include <stdio.h>
28 #include <string.h>
29 #include <stdlib.h>
30 #include <math.h>
31
32 #ifdef _WIN32
33 #include "windirent.h"
34 #else
35 #include <dirent.h>
36 #endif /* _WIN32 */
37
38 #ifdef _WIN32
39 #include <windows.h>
40 #else
41 #include <strings.h>
42 #define _stricmp strcasecmp
43 #define _strnicmp strncasecmp
44 #endif /* _WIN32 */
45
46 #include "opj_config.h"
47 #include "openjpeg.h"
48 #include "j2k.h"
49 #include "jp2.h"
50 #include "opj_getopt.h"
51 #include "convert.h"
52 #include "index.h"
53
54 #include "format_defs.h"
55
56 typedef struct dircnt{
57         /** Buffer for holding images read from Directory*/
58         char *filename_buf;
59         /** Pointer to the buffer*/
60         char **filename;
61 }dircnt_t;
62
63
64 typedef struct img_folder{
65         /** The directory path of the folder containing input images*/
66         char *imgdirpath;
67         /** Output format*/
68         const char *out_format;
69         /** Enable option*/
70         char set_imgdir;
71         /** Enable Cod Format for output*/
72         char set_out_format;
73
74 }img_fol_t;
75
76 void decode_help_display(void) {
77         fprintf(stdout,"HELP for j2k_dump\n----\n\n");
78         fprintf(stdout,"- the -h option displays this help information on screen\n\n");
79
80 /* UniPG>> */
81         fprintf(stdout,"List of parameters for the JPEG 2000 "
82 #ifdef USE_JPWL
83                 "+ JPWL "
84 #endif /* USE_JPWL */
85                 "decoder:\n");
86 /* <<UniPG */
87         fprintf(stdout,"\n");
88         fprintf(stdout,"\n");
89         fprintf(stdout,"  -ImgDir \n");
90         fprintf(stdout,"        Image file Directory path \n");
91         fprintf(stdout,"  -i <compressed file>\n");
92         fprintf(stdout,"    REQUIRED only if an Input image directory not specified\n");
93         fprintf(stdout,"    Currently accepts J2K-files, JP2-files and JPT-files. The file type\n");
94         fprintf(stdout,"    is identified based on its suffix.\n");
95         fprintf(stdout,"  -o <output file>\n");
96         fprintf(stdout,"    OPTIONAL\n");
97         fprintf(stdout,"    Output file where file info will be dump.\n");
98         fprintf(stdout,"    By default it will be in the stdout.\n");
99         fprintf(stdout,"\n");
100 }
101
102 /* -------------------------------------------------------------------------- */
103 static void j2k_dump_image(FILE *fd, opj_image_header_t * img);
104 static void j2k_dump_cp(FILE *fd, opj_image_t * img, opj_cp_v2_t * cp);
105
106 int get_num_images(char *imgdirpath){
107         DIR *dir;
108         struct dirent* content; 
109         int num_images = 0;
110
111         /*Reading the input images from given input directory*/
112
113         dir= opendir(imgdirpath);
114         if(!dir){
115                 fprintf(stderr,"Could not open Folder %s\n",imgdirpath);
116                 return 0;
117         }
118         
119         while((content=readdir(dir))!=NULL){
120                 if(strcmp(".",content->d_name)==0 || strcmp("..",content->d_name)==0 )
121                         continue;
122                 num_images++;
123         }
124         return num_images;
125 }
126
127 int load_images(dircnt_t *dirptr, char *imgdirpath){
128         DIR *dir;
129         struct dirent* content; 
130         int i = 0;
131
132         /*Reading the input images from given input directory*/
133
134         dir= opendir(imgdirpath);
135         if(!dir){
136                 fprintf(stderr,"Could not open Folder %s\n",imgdirpath);
137                 return 1;
138         }else   {
139                 fprintf(stderr,"Folder opened successfully\n");
140         }
141         
142         while((content=readdir(dir))!=NULL){
143                 if(strcmp(".",content->d_name)==0 || strcmp("..",content->d_name)==0 )
144                         continue;
145
146                 strcpy(dirptr->filename[i],content->d_name);
147                 i++;
148         }
149         return 0;       
150 }
151
152 int get_file_format(char *filename) {
153         unsigned int i;
154         static const char *extension[] = {"pgx", "pnm", "pgm", "ppm", "bmp","tif", "raw", "tga", "png", "j2k", "jp2", "jpt", "j2c", "jpc"  };
155         static const int format[] = { PGX_DFMT, PXM_DFMT, PXM_DFMT, PXM_DFMT, BMP_DFMT, TIF_DFMT, RAW_DFMT, TGA_DFMT, PNG_DFMT, J2K_CFMT, JP2_CFMT, JPT_CFMT, J2K_CFMT, J2K_CFMT };
156         char * ext = strrchr(filename, '.');
157         if (ext == NULL)
158                 return -1;
159         ext++;
160         if(ext) {
161                 for(i = 0; i < sizeof(format)/sizeof(*format); i++) {
162                         if(_strnicmp(ext, extension[i], 3) == 0) {
163                                 return format[i];
164                         }
165                 }
166         }
167
168         return -1;
169 }
170
171 char get_next_file(int imageno,dircnt_t *dirptr,img_fol_t *img_fol, opj_dparameters_t *parameters){
172         char image_filename[OPJ_PATH_LEN], infilename[OPJ_PATH_LEN],outfilename[OPJ_PATH_LEN],temp_ofname[OPJ_PATH_LEN];
173         char *temp_p, temp1[OPJ_PATH_LEN]="";
174
175         strcpy(image_filename,dirptr->filename[imageno]);
176         fprintf(stderr,"File Number %d \"%s\"\n",imageno,image_filename);
177         parameters->decod_format = get_file_format(image_filename);
178         if (parameters->decod_format == -1)
179                 return 1;
180         sprintf(infilename,"%s/%s",img_fol->imgdirpath,image_filename);
181         strncpy(parameters->infile, infilename, sizeof(infilename));
182
183         //Set output file
184         strcpy(temp_ofname,strtok(image_filename,"."));
185         while((temp_p = strtok(NULL,".")) != NULL){
186                 strcat(temp_ofname,temp1);
187                 sprintf(temp1,".%s",temp_p);
188         }
189         if(img_fol->set_out_format==1){
190                 sprintf(outfilename,"%s/%s.%s",img_fol->imgdirpath,temp_ofname,img_fol->out_format);
191                 strncpy(parameters->outfile, outfilename, sizeof(outfilename));
192         }
193         return 0;
194 }
195
196 /* -------------------------------------------------------------------------- */
197 int parse_cmdline_decoder(int argc, char **argv, opj_dparameters_t *parameters,img_fol_t *img_fol) {
198         /* parse the command line */
199         int totlen, c;
200         opj_option_t long_option[]={
201                 {"ImgDir",REQ_ARG, NULL ,'y'},
202         };
203         const char optlist[] = "i:o:h";
204
205         totlen=sizeof(long_option);
206         img_fol->set_out_format = 0;
207         do {
208                 c = opj_getopt_long(argc, argv,optlist,long_option,totlen);
209                 if (c == -1)
210                         break;
211                 switch (c) {
212                         case 'i':                       /* input file */
213                         {
214                                 char *infile = opj_optarg;
215                                 parameters->decod_format = get_file_format(infile);
216                                 switch(parameters->decod_format) {
217                                         case J2K_CFMT:
218                                         case JP2_CFMT:
219                                         case JPT_CFMT:
220                                                 break;
221                                         default:
222                                                 fprintf(stderr, 
223                                                         "!! Unrecognized format for infile : %s [accept only *.j2k, *.jp2, *.jpc or *.jpt] !!\n\n", 
224                                                         infile);
225                                                 return 1;
226                                 }
227                                 strncpy(parameters->infile, infile, sizeof(parameters->infile)-1);
228                         }
229                         break;
230
231                         /* ------------------------------------------------------ */
232
233                         case 'o':     /* output file */
234                         {
235                           char *outfile = opj_optarg;
236                           strncpy(parameters->outfile, outfile, sizeof(parameters->outfile)-1);
237                         }
238                         break;
239                                 
240                                 /* ----------------------------------------------------- */
241
242                         case 'h':                       /* display an help description */
243                                 decode_help_display();
244                                 return 1;                               
245
246                                 /* ------------------------------------------------------ */
247
248                         case 'y':                       /* Image Directory path */
249                                 {
250                                         img_fol->imgdirpath = (char*)malloc(strlen(opj_optarg) + 1);
251                                         strcpy(img_fol->imgdirpath,opj_optarg);
252                                         img_fol->set_imgdir=1;
253                                 }
254                                 break;
255
256                                 /* ----------------------------------------------------- */
257                         
258                         default:
259                                 fprintf(stderr,"WARNING -> this option is not valid \"-%c %s\"\n",c, opj_optarg);
260                                 break;
261                 }
262         }while(c != -1);
263
264         /* check for possible errors */
265         if(img_fol->set_imgdir==1){
266                 if(!(parameters->infile[0]==0)){
267                         fprintf(stderr, "Error: options -ImgDir and -i cannot be used together !!\n");
268                         return 1;
269                 }
270                 if(img_fol->set_out_format == 0){
271                         fprintf(stderr, "Error: When -ImgDir is used, -OutFor <FORMAT> must be used !!\n");
272                         fprintf(stderr, "Only one format allowed! Valid format PGM, PPM, PNM, PGX, BMP, TIF, RAW and TGA!!\n");
273                         return 1;
274                 }
275                 if(!((parameters->outfile[0] == 0))){
276                         fprintf(stderr, "Error: options -ImgDir and -o cannot be used together !!\n");
277                         return 1;
278                 }
279         }else{
280                 if((parameters->infile[0] == 0) ) {
281                         fprintf(stderr, "Example: %s -i image.j2k\n",argv[0]);
282                         fprintf(stderr, "    Try: %s -h\n",argv[0]);
283                         return 1;
284                 }
285         }
286
287         return 0;
288 }
289
290 /* -------------------------------------------------------------------------- */
291
292 /**
293 sample error callback expecting a FILE* client object
294 */
295 void error_callback(const char *msg, void *client_data) {
296         FILE *stream = (FILE*)client_data;
297         fprintf(stream, "[ERROR] %s", msg);
298 }
299 /**
300 sample warning callback expecting a FILE* client object
301 */
302 void warning_callback(const char *msg, void *client_data) {
303         FILE *stream = (FILE*)client_data;
304         fprintf(stream, "[WARNING] %s", msg);
305 }
306 /**
307 sample debug callback expecting no client object
308 */
309 void info_callback(const char *msg, void *client_data) {
310         (void)client_data;
311         fprintf(stdout, "[INFO] %s", msg);
312 }
313
314 /* -------------------------------------------------------------------------- */
315
316 int main(int argc, char *argv[])
317 {
318         int ret;
319         opj_dparameters_t parameters;   /* decompression parameters */
320         img_fol_t img_fol;
321         opj_event_mgr_t event_mgr;              /* event manager */
322         opj_image_header_t* image = NULL;
323         FILE *fsrc = NULL, *fout = NULL;
324         opj_bool bResult;
325         int num_images;
326         int i,imageno;
327         dircnt_t *dirptr = NULL;
328         opj_codec_t* dinfo = NULL;      /* handle to a decompressor */
329         opj_stream_t *cio = NULL;
330         opj_codestream_info_t* cstr_info =NULL;  /* Codestream information structure */
331         //OPJ_INT32 l_tile_x0,l_tile_y0;
332         //OPJ_UINT32 l_tile_width,l_tile_height,l_nb_tiles_x,l_nb_tiles_y;
333
334
335         /* FIXME configure the event callbacks (not required) */
336         memset(&event_mgr, 0, sizeof(opj_event_mgr_t));
337         event_mgr.error_handler = error_callback;
338         event_mgr.warning_handler = warning_callback;
339         event_mgr.info_handler = info_callback;
340         event_mgr.client_data = stderr;
341
342         /* set decoding parameters to default values */
343         opj_set_default_decoder_parameters(&parameters);
344
345         /* Initialize img_fol */
346         memset(&img_fol,0,sizeof(img_fol_t));
347
348         /* parse input and get user encoding parameters */
349         if(parse_cmdline_decoder(argc, argv, &parameters,&img_fol) == 1) {
350                 return EXIT_FAILURE;
351         }
352
353         /* Initialize reading of directory */
354         if(img_fol.set_imgdir==1){      
355                 num_images=get_num_images(img_fol.imgdirpath);
356
357                 dirptr=(dircnt_t*)malloc(sizeof(dircnt_t));
358                 if(dirptr){
359                         dirptr->filename_buf = (char*)malloc(num_images*OPJ_PATH_LEN*sizeof(char));     // Stores at max 10 image file names
360                         dirptr->filename = (char**) malloc(num_images*sizeof(char*));
361
362                         if(!dirptr->filename_buf){
363                                 return EXIT_FAILURE;
364                         }
365
366                         for(i=0;i<num_images;i++){
367                                 dirptr->filename[i] = dirptr->filename_buf + i*OPJ_PATH_LEN;
368                         }
369                 }
370                 if(load_images(dirptr,img_fol.imgdirpath)==1){
371                         return EXIT_FAILURE;
372                 }
373
374                 if (num_images==0){
375                         fprintf(stdout,"Folder is empty\n");
376                         return EXIT_FAILURE;
377                 }
378         }else{
379                 num_images=1;
380         }
381
382         //
383         if (parameters.outfile[0] != 0){
384                 fout = fopen(parameters.outfile,"w");
385                 if (!fout){
386                         fprintf(stderr, "ERROR -> failed to open %s for writing\n", parameters.outfile);
387                         return EXIT_FAILURE;
388                 }
389         }
390         else
391                 fout = stdout;
392
393         /*Read the header of each image one by one*/
394         for(imageno = 0; imageno < num_images ; imageno++){
395                 image = NULL;
396                 fprintf(stderr,"\n");
397
398                 if(img_fol.set_imgdir==1){
399                         if (get_next_file(imageno, dirptr,&img_fol, &parameters)) {
400                                 fprintf(stderr,"skipping file...\n");
401                                 continue;
402                         }
403                 }
404
405                 /*NEW V2 STYLE*/
406                 /* read the input file and put it in memory */
407                 /* ---------------------------------------- */
408                 fsrc = fopen(parameters.infile, "rb");
409                 if (!fsrc) {
410                         fprintf(stderr, "ERROR -> failed to open %s for reading\n", parameters.infile);
411                         return EXIT_FAILURE;
412                 }
413
414                 cio = opj_stream_create_default_file_stream(fsrc,1);
415
416                 /* decode the code-stream */
417                 /* ---------------------- */
418
419                 switch(parameters.decod_format) {
420                         case J2K_CFMT:
421                         {       /* JPEG-2000 codestream */
422
423                                 /* get a decoder handle */
424                                 dinfo = opj_create_decompress_v2(CODEC_J2K);
425                                 break;
426                         }
427                         case JP2_CFMT:
428                         {       /* JPEG 2000 compressed image data */
429
430                                 /* get a decoder handle */
431                                 dinfo = opj_create_decompress_v2(CODEC_JP2);
432                                 break;
433                         }
434                         case JPT_CFMT:
435                         {       /* JPEG 2000, JPIP */
436
437                                 /* get a decoder handle */
438                                 dinfo = opj_create_decompress_v2(CODEC_JPT);
439                                 break;
440                         }
441                         default:
442                                 fprintf(stderr, "skipping file..\n");
443                                 opj_stream_destroy(cio);
444                                 continue;
445                 }
446
447                 /* setup the decoder decoding parameters using user parameters */
448                 opj_setup_decoder_v2(dinfo, &parameters, &event_mgr);
449
450                 /* decode the stream and fill the image structure */
451                 /*              if (*indexfilename)                             // If need to extract codestream information
452                                 image = opj_decode_with_info(dinfo, cio, &cstr_info);
453                         else
454                         */
455                 //opj_codestream_info_t *cstr_info = NULL;
456
457                 /*bResult = opj_read_header(    dinfo,
458                                                                         &image,
459                                                                         &l_tile_x0,
460                                                                         &l_tile_y0,
461                                                                         &l_tile_width,
462                                                                         &l_tile_height,
463                                                                         &l_nb_tiles_x,
464                                                                         &l_nb_tiles_y,
465                                                                         cio);*/
466
467                 bResult = opj_read_header(cio, dinfo, &image, &cstr_info);
468
469                 if(!bResult){
470                         fprintf(stderr, "ERROR -> j2k_to_image: failed to read header\n");
471                         return EXIT_FAILURE;
472                 }
473
474                 if (parameters.decod_format == J2K_CFMT){
475                         /* dump image */
476                         j2k_dump_image(fout, image);
477                         /* dump cp */
478                         //j2k_dump_cp(stdout, image, ((opj_codec_private_t)dinfo)->m_codec);
479                         //j2k_dump_cp(fout, image, ((opj_j2k_t*)dinfo->j2k_handle)->cp);
480                 }
481                 else if (parameters.decod_format == JP2_CFMT){
482                         /* dump image */
483                         if(image->icc_profile_buf){
484                                 free(image->icc_profile_buf);
485                                 image->icc_profile_buf = NULL;
486                         }
487                         j2k_dump_image(fout, image);
488                         /* dump cp */
489                         //j2k_dump_cp(stdout, image, dinfo->m_codec);
490                         //j2k_dump_cp(fout, image, ((opj_jp2_t*)dinfo->jp2_handle)->j2k->cp);
491                 }
492
493                 /* close the byte stream */
494                 opj_stream_destroy(cio);
495                 fclose(fsrc);
496
497                 /* Write the index to disk */
498                 /*if (*indexfilename) {
499                         char bSuccess;
500                         bSuccess = write_index_file(&cstr_info, indexfilename);
501                         if (bSuccess) {
502                                 fprintf(stderr, "Failed to output index file\n");
503                                 ret = EXIT_FAILURE;
504                         }
505                         else
506                                 ret = EXIT_SUCCESS;
507                 }
508                 else*/
509                         ret = EXIT_SUCCESS;
510
511                 /* free remaining structures */
512                 if (dinfo) {
513                         opj_destroy_codec(dinfo);
514                 }
515                 /* free codestream information structure */
516                 /*if (*indexfilename)
517                         FIXME A GARDER opj_destroy_cstr_info(&cstr_info);*/
518                 /* free image data structure */
519                 opj_image_header_destroy(image);
520
521         }
522 /*NEW V2 STYLE*/
523
524         fclose(fout);
525
526   return ret;
527 }
528
529
530 static void j2k_dump_image(FILE *fd, opj_image_header_t * img) {
531         int compno;
532         fprintf(fd, "image {\n");
533         fprintf(fd, "  x0=%d, y0=%d, x1=%d, y1=%d\n", img->x0, img->y0, img->x1, img->y1);
534         fprintf(fd, "  numcomps=%d\n", img->numcomps);
535         for (compno = 0; compno < img->numcomps; compno++) {
536                 opj_image_comp_header_t *comp = &img->comps[compno];
537                 fprintf(fd, "  comp %d {\n", compno);
538                 fprintf(fd, "    dx=%d, dy=%d\n", comp->dx, comp->dy);
539                 fprintf(fd, "    prec=%d\n", comp->prec);
540                 //fprintf(fd, "    bpp=%d\n", comp->bpp);
541                 fprintf(fd, "    sgnd=%d\n", comp->sgnd);
542                 fprintf(fd, "  }\n");
543         }
544         fprintf(fd, "}\n");
545 }
546
547 static void j2k_dump_cp(FILE *fd, opj_image_t * img, opj_cp_v2_t * cp) {
548         int tileno, compno, layno, bandno, resno, numbands;
549         fprintf(fd, "coding parameters {\n");
550         fprintf(fd, "  tx0=%d, ty0=%d\n", cp->tx0, cp->ty0);
551         fprintf(fd, "  tdx=%d, tdy=%d\n", cp->tdx, cp->tdy);
552         fprintf(fd, "  tw=%d, th=%d\n", cp->tw, cp->th);
553         for (tileno = 0; tileno < cp->tw * cp->th; tileno++) {
554                 opj_tcp_v2_t *tcp = &cp->tcps[tileno];
555                 fprintf(fd, "  tile %d {\n", tileno);
556                 fprintf(fd, "    csty=%x\n", tcp->csty);
557                 fprintf(fd, "    prg=%d\n", tcp->prg);
558                 fprintf(fd, "    numlayers=%d\n", tcp->numlayers);
559                 fprintf(fd, "    mct=%d\n", tcp->mct);
560                 fprintf(fd, "    rates=");
561                 for (layno = 0; layno < tcp->numlayers; layno++) {
562                         fprintf(fd, "%.1f ", tcp->rates[layno]);
563                 }
564                 fprintf(fd, "\n");
565                 for (compno = 0; compno < img->numcomps; compno++) {
566                         opj_tccp_t *tccp = &tcp->tccps[compno];
567                         fprintf(fd, "    comp %d {\n", compno);
568                         fprintf(fd, "      csty=%x\n", tccp->csty);
569                         fprintf(fd, "      numresolutions=%d\n", tccp->numresolutions);
570                         fprintf(fd, "      cblkw=%d\n", tccp->cblkw);
571                         fprintf(fd, "      cblkh=%d\n", tccp->cblkh);
572                         fprintf(fd, "      cblksty=%x\n", tccp->cblksty);
573                         fprintf(fd, "      qmfbid=%d\n", tccp->qmfbid);
574                         fprintf(fd, "      qntsty=%d\n", tccp->qntsty);
575                         fprintf(fd, "      numgbits=%d\n", tccp->numgbits);
576                         fprintf(fd, "      roishift=%d\n", tccp->roishift);
577                         fprintf(fd, "      stepsizes=");
578                         numbands = tccp->qntsty == J2K_CCP_QNTSTY_SIQNT ? 1 : tccp->numresolutions * 3 - 2;
579                         for (bandno = 0; bandno < numbands; bandno++) {
580                                 fprintf(fd, "(%d,%d) ", tccp->stepsizes[bandno].mant,
581                                         tccp->stepsizes[bandno].expn);
582                         }
583                         fprintf(fd, "\n");
584                         
585                         if (tccp->csty & J2K_CCP_CSTY_PRT) {
586                                 fprintf(fd, "      prcw=");
587                                 for (resno = 0; resno < tccp->numresolutions; resno++) {
588                                         fprintf(fd, "%d ", tccp->prcw[resno]);
589                                 }
590                                 fprintf(fd, "\n");
591                                 fprintf(fd, "      prch=");
592                                 for (resno = 0; resno < tccp->numresolutions; resno++) {
593                                         fprintf(fd, "%d ", tccp->prch[resno]);
594                                 }
595                                 fprintf(fd, "\n");
596                         }
597                         fprintf(fd, "    }\n");
598                 } /*end of component*/
599                 fprintf(fd, "  }\n");
600         } /*end of tile */
601         fprintf(fd, "}\n");
602 }
603