From 2be20ce7d9996e960d79f3ad6ec439a9895849ed Mon Sep 17 00:00:00 2001 From: Even Rouault Date: Wed, 26 Jul 2017 21:04:01 +0200 Subject: Reformat src/bin/wx/OPJViewer/source/OPJThreads.cpp src/bin/wx/OPJViewer/source/imagjpeg2000.cpp wrapping/java/openjp2/JavaOpenJPEG.c --- wrapping/java/openjp2/JavaOpenJPEG.c | 4115 ++++++++++++++++++---------------- 1 file changed, 2141 insertions(+), 1974 deletions(-) (limited to 'wrapping/java/openjp2/JavaOpenJPEG.c') diff --git a/wrapping/java/openjp2/JavaOpenJPEG.c b/wrapping/java/openjp2/JavaOpenJPEG.c index 368cbfdb..dad2dc66 100644 --- a/wrapping/java/openjp2/JavaOpenJPEG.c +++ b/wrapping/java/openjp2/JavaOpenJPEG.c @@ -1,1974 +1,2141 @@ -/* - * Copyright (c) 2002-2014, Universite catholique de Louvain (UCL), Belgium - * Copyright (c) 2002-2014, Professor Benoit Macq - * Copyright (c) 2001-2003, David Janssens - * Copyright (c) 2002-2003, Yannick Verschueren - * Copyright (c) 2003-2007, Francois-Olivier Devaux - * Copyright (c) 2003-2014, Antonin Descampe - * Copyright (c) 2005, Herve Drolon, FreeImage Team - * Copyright (c) 2006-2007, Parvatha Elangovan - * Copyright (c) 2007, Patrick Piscaglia (Telemis) - * All rights reserved. - * - * Redistribution and use in source and binary forms, with or without - * modification, are permitted provided that the following conditions - * are met: - * 1. Redistributions of source code must retain the above copyright - * notice, this list of conditions and the following disclaimer. - * 2. Redistributions in binary form must reproduce the above copyright - * notice, this list of conditions and the following disclaimer in the - * documentation and/or other materials provided with the distribution. - * - * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS' - * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE - * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE - * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE - * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR - * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF - * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS - * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN - * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) - * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE - * POSSIBILITY OF SUCH DAMAGE. - */ -#include -#include -#include -#include -#include - -#include "openjpeg.h" -#include "opj_includes.h" -#include "opj_getopt.h" -#include "convert.h" -#include "index.h" -#include "dirent.h" -#include "org_openJpeg_OpenJPEGJavaEncoder.h" - -#ifndef _WIN32 -#define stricmp strcasecmp -#define strnicmp strncasecmp -#endif - -#include "format_defs.h" - -#define CINEMA_24_CS 1302083 /*Codestream length for 24fps*/ -#define CINEMA_48_CS 651041 /*Codestream length for 48fps*/ -#define COMP_24_CS 1041666 /*Maximum size per color component for 2K & 4K @ 24fps*/ -#define COMP_48_CS 520833 /*Maximum size per color component for 2K @ 48fps*/ - -extern int get_file_format(char *filename); -extern void error_callback(const char *msg, void *client_data); -extern void warning_callback(const char *msg, void *client_data); -extern void info_callback(const char *msg, void *client_data); - -typedef struct callback_variables { - JNIEnv *env; - /** 'jclass' object used to call a Java method from the C */ - jobject *jobj; - /** 'jclass' object used to call a Java method from the C */ - jmethodID message_mid; - jmethodID error_mid; -} callback_variables_t; - -typedef struct dircnt{ - /** Buffer for holding images read from Directory*/ - char *filename_buf; - /** Pointer to the buffer*/ - char **filename; -}dircnt_t; - -typedef struct img_folder{ - /** The directory path of the folder containing input images*/ - char *imgdirpath; - /** Output format*/ - char *out_format; - /** Enable option*/ - char set_imgdir; - /** Enable Cod Format for output*/ - char set_out_format; - /** User specified rate stored in case of cinema option*/ - float *rates; -}img_fol_t; - -static void encode_help_display() { - fprintf(stdout,"HELP\n----\n\n"); - fprintf(stdout,"- the -h option displays this help information on screen\n\n"); - -/* UniPG>> */ - fprintf(stdout,"List of parameters for the JPEG 2000 " -#ifdef USE_JPWL - "+ JPWL " -#endif /* USE_JPWL */ - "encoder:\n"); -/* <> */ -#ifdef USE_JPWL - fprintf(stdout," * No JPWL protection\n"); -#endif /* USE_JPWL */ -/* < \n"); - fprintf(stdout," Currently accepts PGM, PPM, PNM, PGX, BMP, TIF, RAW and TGA formats\n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-i : source file (-i source.pnm also *.pgm, *.ppm, *.bmp, *.tif, *.raw, *.tga) \n"); - fprintf(stdout," When using this option -o must be used\n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-o : destination file (-o dest.j2k or .jp2) \n"); - fprintf(stdout,"\n"); - fprintf(stdout,"Optional Parameters:\n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-h : display the help information \n "); - fprintf(stdout,"\n"); - fprintf(stdout,"-cinema2K : Digital Cinema 2K profile compliant codestream for 2K resolution.(-cinema2k 24 or 48) \n"); - fprintf(stdout," Need to specify the frames per second for a 2K resolution. Only 24 or 48 fps is allowed\n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-cinema4K : Digital Cinema 4K profile compliant codestream for 4K resolution \n"); - fprintf(stdout," Frames per second not required. Default value is 24fps\n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-r : different compression ratios for successive layers (-r 20,10,5)\n "); - fprintf(stdout," - The rate specified for each quality level is the desired \n"); - fprintf(stdout," compression factor.\n"); - fprintf(stdout," Example: -r 20,10,1 means quality 1: compress 20x, \n"); - fprintf(stdout," quality 2: compress 10x and quality 3: compress lossless\n"); - fprintf(stdout,"\n"); - fprintf(stdout," (options -r and -q cannot be used together)\n "); - fprintf(stdout,"\n"); - - fprintf(stdout,"-q : different psnr for successive layers (-q 30,40,50) \n "); - - fprintf(stdout," (options -r and -q cannot be used together)\n "); - - fprintf(stdout,"\n"); - fprintf(stdout,"-n : number of resolutions (-n 3) \n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-b : size of code block (-b 32,32) \n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-c : size of precinct (-c 128,128) \n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-t : size of tile (-t 512,512) \n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-p : progression order (-p LRCP) [LRCP, RLCP, RPCL, PCRL, CPRL] \n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-s : subsampling factor (-s 2,2) [-s X,Y] \n"); - fprintf(stdout," Remark: subsampling bigger than 2 can produce error\n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-POC : Progression order change (-POC T1=0,0,1,5,3,CPRL/T1=5,0,1,6,3,CPRL) \n"); - fprintf(stdout," Example: T1=0,0,1,5,3,CPRL \n"); - fprintf(stdout," : Ttilenumber=Resolution num start,Component num start,Layer num end,Resolution num end,Component num end,Progression order\n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-SOP : write SOP marker before each packet \n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-EPH : write EPH marker after each header packet \n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-M : mode switch (-M 3) [1=BYPASS(LAZY) 2=RESET 4=RESTART(TERMALL)\n"); - fprintf(stdout," 8=VSC 16=ERTERM(SEGTERM) 32=SEGMARK(SEGSYM)] \n"); - fprintf(stdout," Indicate multiple modes by adding their values. \n"); - fprintf(stdout," ex: RESTART(4) + RESET(2) + SEGMARK(32) = -M 38\n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-TP : divide packets of every tile into tile-parts (-TP R) [R, L, C]\n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-x : create an index file *.Idx (-x index_name.Idx) \n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-ROI : c=%%d,U=%%d : quantization indices upshifted \n"); - fprintf(stdout," for component c=%%d [%%d = 0,1,2]\n"); - fprintf(stdout," with a value of U=%%d [0 <= %%d <= 37] (i.e. -ROI c=0,U=25) \n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-d : offset of the origin of the image (-d 150,300) \n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-T : offset of the origin of the tiles (-T 100,75) \n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-I : use the irreversible DWT 9-7 (-I) \n"); - fprintf(stdout,"\n"); - fprintf(stdout,"-jpip : write jpip codestream index box in JP2 output file\n"); - fprintf(stdout," NOTICE: currently supports only RPCL order\n"); - fprintf(stdout,"\n"); -/* UniPG>> */ -#ifdef USE_JPWL - fprintf(stdout,"-W : adoption of JPWL (Part 11) capabilities (-W params)\n"); - fprintf(stdout," The parameters can be written and repeated in any order:\n"); - fprintf(stdout," [h<=type>,s<=method>,a=,...\n"); - fprintf(stdout," ...,z=,g=,p<=type>]\n"); - fprintf(stdout,"\n"); - fprintf(stdout," h selects the header error protection (EPB): 'type' can be\n"); - fprintf(stdout," [0=none 1,absent=predefined 16=CRC-16 32=CRC-32 37-128=RS]\n"); - fprintf(stdout," if 'tilepart' is absent, it is for main and tile headers\n"); - fprintf(stdout," if 'tilepart' is present, it applies from that tile\n"); - fprintf(stdout," onwards, up to the next h<> spec, or to the last tilepart\n"); - fprintf(stdout," in the codestream (max. %d specs)\n", JPWL_MAX_NO_TILESPECS); - fprintf(stdout,"\n"); - fprintf(stdout," p selects the packet error protection (EEP/UEP with EPBs)\n"); - fprintf(stdout," to be applied to raw data: 'type' can be\n"); - fprintf(stdout," [0=none 1,absent=predefined 16=CRC-16 32=CRC-32 37-128=RS]\n"); - fprintf(stdout," if 'tilepart:pack' is absent, it is from tile 0, packet 0\n"); - fprintf(stdout," if 'tilepart:pack' is present, it applies from that tile\n"); - fprintf(stdout," and that packet onwards, up to the next packet spec\n"); - fprintf(stdout," or to the last packet in the last tilepart in the stream\n"); - fprintf(stdout," (max. %d specs)\n", JPWL_MAX_NO_PACKSPECS); - fprintf(stdout,"\n"); - fprintf(stdout," s enables sensitivity data insertion (ESD): 'method' can be\n"); - fprintf(stdout," [-1=NO ESD 0=RELATIVE ERROR 1=MSE 2=MSE REDUCTION 3=PSNR\n"); - fprintf(stdout," 4=PSNR INCREMENT 5=MAXERR 6=TSE 7=RESERVED]\n"); - fprintf(stdout," if 'tilepart' is absent, it is for main header only\n"); - fprintf(stdout," if 'tilepart' is present, it applies from that tile\n"); - fprintf(stdout," onwards, up to the next s<> spec, or to the last tilepart\n"); - fprintf(stdout," in the codestream (max. %d specs)\n", JPWL_MAX_NO_TILESPECS); - fprintf(stdout,"\n"); - fprintf(stdout," g determines the addressing mode: can be\n"); - fprintf(stdout," [0=PACKET 1=BYTE RANGE 2=PACKET RANGE]\n"); - fprintf(stdout,"\n"); - fprintf(stdout," a determines the size of data addressing: can be\n"); - fprintf(stdout," 2/4 bytes (small/large codestreams). If not set, auto-mode\n"); - fprintf(stdout,"\n"); - fprintf(stdout," z determines the size of sensitivity values: can be\n"); - fprintf(stdout," 1/2 bytes, for the transformed pseudo-floating point value\n"); - fprintf(stdout,"\n"); - fprintf(stdout," ex.:\n"); - fprintf(stdout," h,h0=64,h3=16,h5=32,p0=78,p0:24=56,p1,p3:0=0,p3:20=32,s=0,\n"); - fprintf(stdout," s0=6,s3=-1,a=0,g=1,z=1\n"); - fprintf(stdout," means\n"); - fprintf(stdout," predefined EPB in MH, rs(64,32) from TPH 0 to TPH 2,\n"); - fprintf(stdout," CRC-16 in TPH 3 and TPH 4, CRC-32 in remaining TPHs,\n"); - fprintf(stdout," UEP rs(78,32) for packets 0 to 23 of tile 0,\n"); - fprintf(stdout," UEP rs(56,32) for packs. 24 to the last of tilepart 0,\n"); - fprintf(stdout," UEP rs default for packets of tilepart 1,\n"); - fprintf(stdout," no UEP for packets 0 to 19 of tilepart 3,\n"); - fprintf(stdout," UEP CRC-32 for packs. 20 of tilepart 3 to last tilepart,\n"); - fprintf(stdout," relative sensitivity ESD for MH,\n"); - fprintf(stdout," TSE ESD from TPH 0 to TPH 2, byte range with automatic\n"); - fprintf(stdout," size of addresses and 1 byte for each sensitivity value\n"); - fprintf(stdout,"\n"); - fprintf(stdout," ex.:\n"); - fprintf(stdout," h,s,p\n"); - fprintf(stdout," means\n"); - fprintf(stdout," default protection to headers (MH and TPHs) as well as\n"); - fprintf(stdout," data packets, one ESD in MH\n"); - fprintf(stdout,"\n"); - fprintf(stdout," N.B.: use the following recommendations when specifying\n"); - fprintf(stdout," the JPWL parameters list\n"); - fprintf(stdout," - when you use UEP, always pair the 'p' option with 'h'\n"); - fprintf(stdout," \n"); -#endif /* USE_JPWL */ -/* <tile_size_on = OPJ_FALSE; - parameters->cp_tdx=1; - parameters->cp_tdy=1; - - /*Tile part*/ - parameters->tp_flag = 'C'; - parameters->tp_on = 1; - - /*Tile and Image shall be at (0,0)*/ - parameters->cp_tx0 = 0; - parameters->cp_ty0 = 0; - parameters->image_offset_x0 = 0; - parameters->image_offset_y0 = 0; - - /*Codeblock size= 32*32*/ - parameters->cblockw_init = 32; - parameters->cblockh_init = 32; - parameters->csty |= 0x01; - - /*The progression order shall be CPRL*/ - parameters->prog_order = CPRL; - - /* No ROI */ - parameters->roi_compno = -1; - - parameters->subsampling_dx = 1; parameters->subsampling_dy = 1; - - /* 9-7 transform */ - parameters->irreversible = 1; - -} - -static void cinema_setup_encoder(opj_cparameters_t *parameters,opj_image_t *image, img_fol_t *img_fol){ - int i; - float temp_rate; - opj_poc_t *POC = NULL; - - switch (parameters->cp_cinema){ - case CINEMA2K_24: - case CINEMA2K_48: - if(parameters->numresolution > 6){ - parameters->numresolution = 6; - } - if (!((image->comps[0].w == 2048) | (image->comps[0].h == 1080))){ - fprintf(stdout,"Image coordinates %d x %d is not 2K compliant.\nJPEG Digital Cinema Profile-3" - "(2K profile) compliance requires that at least one of coordinates match 2048 x 1080\n", - image->comps[0].w,image->comps[0].h); - parameters->cp_rsiz = STD_RSIZ; - } - break; - - case CINEMA4K_24: - if(parameters->numresolution < 1){ - parameters->numresolution = 1; - }else if(parameters->numresolution > 7){ - parameters->numresolution = 7; - } - if (!((image->comps[0].w == 4096) | (image->comps[0].h == 2160))){ - fprintf(stdout,"Image coordinates %d x %d is not 4K compliant.\nJPEG Digital Cinema Profile-4" - "(4K profile) compliance requires that at least one of coordinates match 4096 x 2160\n", - image->comps[0].w,image->comps[0].h); - parameters->cp_rsiz = STD_RSIZ; - } - parameters->numpocs = initialise_4K_poc(parameters->POC,parameters->numresolution); - break; - } - - switch (parameters->cp_cinema){ - case CINEMA2K_24: - case CINEMA4K_24: - for(i=0 ; itcp_numlayers ; i++){ - temp_rate = 0 ; - if (img_fol->rates[i]== 0){ - parameters->tcp_rates[0]= ((float) (image->numcomps * image->comps[0].w * image->comps[0].h * image->comps[0].prec))/ - (CINEMA_24_CS * 8 * image->comps[0].dx * image->comps[0].dy); - }else{ - temp_rate =((float) (image->numcomps * image->comps[0].w * image->comps[0].h * image->comps[0].prec))/ - (img_fol->rates[i] * 8 * image->comps[0].dx * image->comps[0].dy); - if (temp_rate > CINEMA_24_CS ){ - parameters->tcp_rates[i]= ((float) (image->numcomps * image->comps[0].w * image->comps[0].h * image->comps[0].prec))/ - (CINEMA_24_CS * 8 * image->comps[0].dx * image->comps[0].dy); - }else{ - parameters->tcp_rates[i]= img_fol->rates[i]; - } - } - } - parameters->max_comp_size = COMP_24_CS; - break; - - case CINEMA2K_48: - for(i=0 ; itcp_numlayers ; i++){ - temp_rate = 0 ; - if (img_fol->rates[i]== 0){ - parameters->tcp_rates[0]= ((float) (image->numcomps * image->comps[0].w * image->comps[0].h * image->comps[0].prec))/ - (CINEMA_48_CS * 8 * image->comps[0].dx * image->comps[0].dy); - }else{ - temp_rate =((float) (image->numcomps * image->comps[0].w * image->comps[0].h * image->comps[0].prec))/ - (img_fol->rates[i] * 8 * image->comps[0].dx * image->comps[0].dy); - if (temp_rate > CINEMA_48_CS ){ - parameters->tcp_rates[0]= ((float) (image->numcomps * image->comps[0].w * image->comps[0].h * image->comps[0].prec))/ - (CINEMA_48_CS * 8 * image->comps[0].dx * image->comps[0].dy); - }else{ - parameters->tcp_rates[i]= img_fol->rates[i]; - } - } - } - parameters->max_comp_size = COMP_48_CS; - break; - } - parameters->cp_disto_alloc = 1; -} - - -/* ------------------------------------------------------------------------------------ */ -static int parse_cmdline_encoder(int argc, char **argv, opj_cparameters_t *parameters, - img_fol_t *img_fol, char *indexfilename) { - int i, j,totlen; - opj_option_t long_option[]={ - {"cinema2K",REQ_ARG, NULL ,'w'}, - {"cinema4K",NO_ARG, NULL ,'y'}, - {"ImgDir",REQ_ARG, NULL ,'z'}, - {"TP",REQ_ARG, NULL ,'u'}, - {"SOP",NO_ARG, NULL ,'S'}, - {"EPH",NO_ARG, NULL ,'E'}, - {"OutFor",REQ_ARG, NULL ,'O'}, - {"POC",REQ_ARG, NULL ,'P'}, - {"ROI",REQ_ARG, NULL ,'R'}, - {"jpip",NO_ARG, NULL, 'J'} - }; - - /* parse the command line */ -/* UniPG>> */ - const char optlist[] = "i:o:hr:q:n:b:c:t:p:s:SEM:x:R:d:T:If:P:C:F:u:" -#ifdef USE_JPWL - "W:" -#endif /* USE_JPWL */ - ; - - /*printf("C: parse_cmdline_encoder:"); - for (i=0; iset_out_format=0; - reset_options_reading(); - - while (1) { - int c = opj_getopt_long(argc, argv, optlist,long_option,totlen); - if (c == -1) - break; - switch (c) { - - /* ----------------------------------------------------- */ - - case 'o': /* output file */ - { - char *outfile = opj_optarg; - parameters->cod_format = get_file_format(outfile); - switch(parameters->cod_format) { - case J2K_CFMT: - case JP2_CFMT: - break; - default: - fprintf(stderr, "Unknown output format image %s [only *.j2k, *.j2c or *.jp2]!! \n", outfile); - return 1; - } - strncpy(parameters->outfile, outfile, sizeof(parameters->outfile)-1); - } - break; - - /* ----------------------------------------------------- */ - case 'O': /* output format */ - { - char outformat[50]; - char *of = opj_optarg; - sprintf(outformat,".%s",of); - img_fol->set_out_format = 1; - parameters->cod_format = get_file_format(outformat); - switch(parameters->cod_format) { - case J2K_CFMT: - case JP2_CFMT: - img_fol->out_format = opj_optarg; - break; - default: - fprintf(stderr, "Unknown output format image [only j2k, j2c, jp2]!! \n"); - return 1; - } - } - break; - - - /* ----------------------------------------------------- */ - - - case 'r': /* rates rates/distorsion */ - { - char *s = opj_optarg; - while (sscanf(s, "%f", ¶meters->tcp_rates[parameters->tcp_numlayers]) == 1) { - parameters->tcp_numlayers++; - while (*s && *s != ',') { - s++; - } - if (!*s) - break; - s++; - } - parameters->cp_disto_alloc = 1; - } - break; - - /* ----------------------------------------------------- */ - - case 'q': /* add fixed_quality */ - { - char *s = opj_optarg; - while (sscanf(s, "%f", ¶meters->tcp_distoratio[parameters->tcp_numlayers]) == 1) { - parameters->tcp_numlayers++; - while (*s && *s != ',') { - s++; - } - if (!*s) - break; - s++; - } - parameters->cp_fixed_quality = 1; - } - break; - - /* dda */ - /* ----------------------------------------------------- */ - - case 'f': /* mod fixed_quality (before : -q) */ - { - int *row = NULL, *col = NULL; - int numlayers = 0, numresolution = 0, matrix_width = 0; - - char *s = opj_optarg; - sscanf(s, "%d", &numlayers); - s++; - if (numlayers > 9) - s++; - - parameters->tcp_numlayers = numlayers; - numresolution = parameters->numresolution; - matrix_width = numresolution * 3; - parameters->cp_matrice = (int *) opj_malloc(numlayers * matrix_width * sizeof(int)); - s = s + 2; - - for (i = 0; i < numlayers; i++) { - row = ¶meters->cp_matrice[i * matrix_width]; - col = row; - parameters->tcp_rates[i] = 1; - sscanf(s, "%d,", &col[0]); - s += 2; - if (col[0] > 9) - s++; - col[1] = 0; - col[2] = 0; - for (j = 1; j < numresolution; j++) { - col += 3; - sscanf(s, "%d,%d,%d", &col[0], &col[1], &col[2]); - s += 6; - if (col[0] > 9) - s++; - if (col[1] > 9) - s++; - if (col[2] > 9) - s++; - } - if (i < numlayers - 1) - s++; - } - parameters->cp_fixed_alloc = 1; - } - break; - - /* ----------------------------------------------------- */ - - case 't': /* tiles */ - { - sscanf(opj_optarg, "%d,%d", ¶meters->cp_tdx, ¶meters->cp_tdy); - parameters->tile_size_on = OPJ_TRUE; - } - break; - - /* ----------------------------------------------------- */ - - case 'n': /* resolution */ - { - sscanf(opj_optarg, "%d", ¶meters->numresolution); - } - break; - - /* ----------------------------------------------------- */ - case 'c': /* precinct dimension */ - { - char sep; - int res_spec = 0; - - char *s = opj_optarg; - do { - sep = 0; - sscanf(s, "[%d,%d]%c", ¶meters->prcw_init[res_spec], - ¶meters->prch_init[res_spec], &sep); - parameters->csty |= 0x01; - res_spec++; - s = strpbrk(s, "]") + 2; - } - while (sep == ','); - parameters->res_spec = res_spec; - } - break; - - /* ----------------------------------------------------- */ - - case 'b': /* code-block dimension */ - { - int cblockw_init = 0, cblockh_init = 0; - sscanf(opj_optarg, "%d,%d", &cblockw_init, &cblockh_init); - if (cblockw_init * cblockh_init > 4096 || cblockw_init > 1024 - || cblockw_init < 4 || cblockh_init > 1024 || cblockh_init < 4) { - fprintf(stderr, - "!! Size of code_block error (option -b) !!\n\nRestriction :\n" - " * width*height<=4096\n * 4<=width,height<= 1024\n\n"); - return 1; - } - parameters->cblockw_init = cblockw_init; - parameters->cblockh_init = cblockh_init; - } - break; - - /* ----------------------------------------------------- */ - - case 'x': /* creation of index file */ - { - char *index = opj_optarg; - strncpy(indexfilename, index, OPJ_PATH_LEN); - } - break; - - /* ----------------------------------------------------- */ - - case 'p': /* progression order */ - { - char progression[4]; - - strncpy(progression, opj_optarg, 4); - parameters->prog_order = give_progression(progression); - if (parameters->prog_order == -1) { - fprintf(stderr, "Unrecognized progression order " - "[LRCP, RLCP, RPCL, PCRL, CPRL] !!\n"); - return 1; - } - } - break; - - /* ----------------------------------------------------- */ - - case 's': /* subsampling factor */ - { - if (sscanf(opj_optarg, "%d,%d", ¶meters->subsampling_dx, - ¶meters->subsampling_dy) != 2) { - fprintf(stderr, "'-s' sub-sampling argument error ! [-s dx,dy]\n"); - return 1; - } - } - break; - - /* ----------------------------------------------------- */ - - case 'd': /* coordonnate of the reference grid */ - { - if (sscanf(opj_optarg, "%d,%d", ¶meters->image_offset_x0, - ¶meters->image_offset_y0) != 2) { - fprintf(stderr, "-d 'coordonnate of the reference grid' argument " - "error !! [-d x0,y0]\n"); - return 1; - } - } - break; - - /* ----------------------------------------------------- */ - - case 'h': /* display an help description */ - encode_help_display(); - return 1; - - /* ----------------------------------------------------- */ - - case 'P': /* POC */ - { - int numpocs = 0; /* number of progression order change (POC) default 0 */ - opj_poc_t *POC = NULL; /* POC : used in case of Progression order change */ - - char *s = opj_optarg; - POC = parameters->POC; - - while (sscanf(s, "T%d=%d,%d,%d,%d,%d,%4s", &POC[numpocs].tile, - &POC[numpocs].resno0, &POC[numpocs].compno0, - &POC[numpocs].layno1, &POC[numpocs].resno1, - &POC[numpocs].compno1, &POC[numpocs].progorder) == 7) { - POC[numpocs].prg1 = give_progression(POC[numpocs].progorder); - numpocs++; - while (*s && *s != '/') { - s++; - } - if (!*s) { - break; - } - s++; - } - parameters->numpocs = numpocs; - } - break; - - /* ------------------------------------------------------ */ - - case 'S': /* SOP marker */ - { - parameters->csty |= 0x02; - } - break; - - /* ------------------------------------------------------ */ - - case 'E': /* EPH marker */ - { - parameters->csty |= 0x04; - } - break; - - /* ------------------------------------------------------ */ - - case 'M': /* Mode switch pas tous au point !! */ - { - int value = 0; - if (sscanf(opj_optarg, "%d", &value) == 1) { - for (i = 0; i <= 5; i++) { - int cache = value & (1 << i); - if (cache) - parameters->mode |= (1 << i); - } - } - } - break; - - /* ------------------------------------------------------ */ - - case 'R': /* ROI */ - { - if (sscanf(opj_optarg, "c=%d,U=%d", ¶meters->roi_compno, - ¶meters->roi_shift) != 2) { - fprintf(stderr, "ROI error !! [-ROI c='compno',U='shift']\n"); - return 1; - } - } - break; - - /* ------------------------------------------------------ */ - - case 'T': /* Tile offset */ - { - if (sscanf(opj_optarg, "%d,%d", ¶meters->cp_tx0, ¶meters->cp_ty0) != 2) { - fprintf(stderr, "-T 'tile offset' argument error !! [-T X0,Y0]"); - return 1; - } - } - break; - - /* ------------------------------------------------------ */ - - case 'C': /* add a comment */ - { - parameters->cp_comment = (char*)opj_malloc(strlen(opj_optarg) + 1); - if(parameters->cp_comment) { - strcpy(parameters->cp_comment, opj_optarg); - } - } - break; - - - /* ------------------------------------------------------ */ - - case 'I': /* reversible or not */ - { - parameters->irreversible = 1; - } - break; - - /* ------------------------------------------------------ */ - - case 'u': /* Tile part generation*/ - { - parameters->tp_flag = opj_optarg[0]; - parameters->tp_on = 1; - } - break; - - /* ------------------------------------------------------ */ - - case 'z': /* Image Directory path */ - { - img_fol->imgdirpath = (char*)opj_malloc(strlen(opj_optarg) + 1); - strcpy(img_fol->imgdirpath,opj_optarg); - img_fol->set_imgdir=1; - } - break; - - /* ------------------------------------------------------ */ - - case 'w': /* Digital Cinema 2K profile compliance*/ - { - int fps=0; - sscanf(opj_optarg,"%d",&fps); - if(fps == 24){ - parameters->cp_cinema = CINEMA2K_24; - }else if(fps == 48 ){ - parameters->cp_cinema = CINEMA2K_48; - }else { - fprintf(stderr,"Incorrect value!! must be 24 or 48\n"); - return 1; - } - fprintf(stdout,"CINEMA 2K compliant codestream\n"); - parameters->cp_rsiz = CINEMA2K; - - } - break; - - /* ------------------------------------------------------ */ - - case 'y': /* Digital Cinema 4K profile compliance*/ - { - parameters->cp_cinema = CINEMA4K_24; - fprintf(stdout,"CINEMA 4K compliant codestream\n"); - parameters->cp_rsiz = CINEMA4K; - } - break; - - /* ------------------------------------------------------ */ - -/* UniPG>> */ -#ifdef USE_JPWL - /* ------------------------------------------------------ */ - - case 'W': /* JPWL capabilities switched on */ - { - char *token = NULL; - int hprot, pprot, sens, addr, size, range; - - /* we need to enable indexing */ - if (!indexfilename) { - strncpy(indexfilename, JPWL_PRIVATEINDEX_NAME, OPJ_PATH_LEN); - } - - /* search for different protection methods */ - - /* break the option in comma points and parse the result */ - token = strtok(opj_optarg, ","); - while(token != NULL) { - - /* search header error protection method */ - if (*token == 'h') { - - static int tile = 0, tilespec = 0, lasttileno = 0; - - hprot = 1; /* predefined method */ - - if(sscanf(token, "h=%d", &hprot) == 1) { - /* Main header, specified */ - if (!((hprot == 0) || (hprot == 1) || (hprot == 16) || (hprot == 32) || - ((hprot >= 37) && (hprot <= 128)))) { - fprintf(stderr, "ERROR -> invalid main header protection method h = %d\n", hprot); - return 1; - } - parameters->jpwl_hprot_MH = hprot; - - } else if(sscanf(token, "h%d=%d", &tile, &hprot) == 2) { - /* Tile part header, specified */ - if (!((hprot == 0) || (hprot == 1) || (hprot == 16) || (hprot == 32) || - ((hprot >= 37) && (hprot <= 128)))) { - fprintf(stderr, "ERROR -> invalid tile part header protection method h = %d\n", hprot); - return 1; - } - if (tile < 0) { - fprintf(stderr, "ERROR -> invalid tile part number on protection method t = %d\n", tile); - return 1; - } - if (tilespec < JPWL_MAX_NO_TILESPECS) { - parameters->jpwl_hprot_TPH_tileno[tilespec] = lasttileno = tile; - parameters->jpwl_hprot_TPH[tilespec++] = hprot; - } - - } else if(sscanf(token, "h%d", &tile) == 1) { - /* Tile part header, unspecified */ - if (tile < 0) { - fprintf(stderr, "ERROR -> invalid tile part number on protection method t = %d\n", tile); - return 1; - } - if (tilespec < JPWL_MAX_NO_TILESPECS) { - parameters->jpwl_hprot_TPH_tileno[tilespec] = lasttileno = tile; - parameters->jpwl_hprot_TPH[tilespec++] = hprot; - } - - - } else if (!strcmp(token, "h")) { - /* Main header, unspecified */ - parameters->jpwl_hprot_MH = hprot; - - } else { - fprintf(stderr, "ERROR -> invalid protection method selection = %s\n", token); - return 1; - }; - - } - - /* search packet error protection method */ - if (*token == 'p') { - - static int pack = 0, tile = 0, packspec = 0, lastpackno = 0; - - pprot = 1; /* predefined method */ - - if (sscanf(token, "p=%d", &pprot) == 1) { - /* Method for all tiles and all packets */ - if (!((pprot == 0) || (pprot == 1) || (pprot == 16) || (pprot == 32) || - ((pprot >= 37) && (pprot <= 128)))) { - fprintf(stderr, "ERROR -> invalid default packet protection method p = %d\n", pprot); - return 1; - } - parameters->jpwl_pprot_tileno[0] = 0; - parameters->jpwl_pprot_packno[0] = 0; - parameters->jpwl_pprot[0] = pprot; - - } else if (sscanf(token, "p%d=%d", &tile, &pprot) == 2) { - /* method specified from that tile on */ - if (!((pprot == 0) || (pprot == 1) || (pprot == 16) || (pprot == 32) || - ((pprot >= 37) && (pprot <= 128)))) { - fprintf(stderr, "ERROR -> invalid packet protection method p = %d\n", pprot); - return 1; - } - if (tile < 0) { - fprintf(stderr, "ERROR -> invalid tile part number on protection method p = %d\n", tile); - return 1; - } - if (packspec < JPWL_MAX_NO_PACKSPECS) { - parameters->jpwl_pprot_tileno[packspec] = tile; - parameters->jpwl_pprot_packno[packspec] = 0; - parameters->jpwl_pprot[packspec++] = pprot; - } - - } else if (sscanf(token, "p%d:%d=%d", &tile, &pack, &pprot) == 3) { - /* method fully specified from that tile and that packet on */ - if (!((pprot == 0) || (pprot == 1) || (pprot == 16) || (pprot == 32) || - ((pprot >= 37) && (pprot <= 128)))) { - fprintf(stderr, "ERROR -> invalid packet protection method p = %d\n", pprot); - return 1; - } - if (tile < 0) { - fprintf(stderr, "ERROR -> invalid tile part number on protection method p = %d\n", tile); - return 1; - } - if (pack < 0) { - fprintf(stderr, "ERROR -> invalid packet number on protection method p = %d\n", pack); - return 1; - } - if (packspec < JPWL_MAX_NO_PACKSPECS) { - parameters->jpwl_pprot_tileno[packspec] = tile; - parameters->jpwl_pprot_packno[packspec] = pack; - parameters->jpwl_pprot[packspec++] = pprot; - } - - } else if (sscanf(token, "p%d:%d", &tile, &pack) == 2) { - /* default method from that tile and that packet on */ - if (!((pprot == 0) || (pprot == 1) || (pprot == 16) || (pprot == 32) || - ((pprot >= 37) && (pprot <= 128)))) { - fprintf(stderr, "ERROR -> invalid packet protection method p = %d\n", pprot); - return 1; - } - if (tile < 0) { - fprintf(stderr, "ERROR -> invalid tile part number on protection method p = %d\n", tile); - return 1; - } - if (pack < 0) { - fprintf(stderr, "ERROR -> invalid packet number on protection method p = %d\n", pack); - return 1; - } - if (packspec < JPWL_MAX_NO_PACKSPECS) { - parameters->jpwl_pprot_tileno[packspec] = tile; - parameters->jpwl_pprot_packno[packspec] = pack; - parameters->jpwl_pprot[packspec++] = pprot; - } - - } else if (sscanf(token, "p%d", &tile) == 1) { - /* default from a tile on */ - if (tile < 0) { - fprintf(stderr, "ERROR -> invalid tile part number on protection method p = %d\n", tile); - return 1; - } - if (packspec < JPWL_MAX_NO_PACKSPECS) { - parameters->jpwl_pprot_tileno[packspec] = tile; - parameters->jpwl_pprot_packno[packspec] = 0; - parameters->jpwl_pprot[packspec++] = pprot; - } - - - } else if (!strcmp(token, "p")) { - /* all default */ - parameters->jpwl_pprot_tileno[0] = 0; - parameters->jpwl_pprot_packno[0] = 0; - parameters->jpwl_pprot[0] = pprot; - - } else { - fprintf(stderr, "ERROR -> invalid protection method selection = %s\n", token); - return 1; - }; - - } - - /* search sensitivity method */ - if (*token == 's') { - - static int tile = 0, tilespec = 0, lasttileno = 0; - - sens = 0; /* predefined: relative error */ - - if(sscanf(token, "s=%d", &sens) == 1) { - /* Main header, specified */ - if ((sens < -1) || (sens > 7)) { - fprintf(stderr, "ERROR -> invalid main header sensitivity method s = %d\n", sens); - return 1; - } - parameters->jpwl_sens_MH = sens; - - } else if(sscanf(token, "s%d=%d", &tile, &sens) == 2) { - /* Tile part header, specified */ - if ((sens < -1) || (sens > 7)) { - fprintf(stderr, "ERROR -> invalid tile part header sensitivity method s = %d\n", sens); - return 1; - } - if (tile < 0) { - fprintf(stderr, "ERROR -> invalid tile part number on sensitivity method t = %d\n", tile); - return 1; - } - if (tilespec < JPWL_MAX_NO_TILESPECS) { - parameters->jpwl_sens_TPH_tileno[tilespec] = lasttileno = tile; - parameters->jpwl_sens_TPH[tilespec++] = sens; - } - - } else if(sscanf(token, "s%d", &tile) == 1) { - /* Tile part header, unspecified */ - if (tile < 0) { - fprintf(stderr, "ERROR -> invalid tile part number on sensitivity method t = %d\n", tile); - return 1; - } - if (tilespec < JPWL_MAX_NO_TILESPECS) { - parameters->jpwl_sens_TPH_tileno[tilespec] = lasttileno = tile; - parameters->jpwl_sens_TPH[tilespec++] = hprot; - } - - } else if (!strcmp(token, "s")) { - /* Main header, unspecified */ - parameters->jpwl_sens_MH = sens; - - } else { - fprintf(stderr, "ERROR -> invalid sensitivity method selection = %s\n", token); - return 1; - }; - - parameters->jpwl_sens_size = 2; /* 2 bytes for default size */ - } - - /* search addressing size */ - if (*token == 'a') { - - static int tile = 0, tilespec = 0, lasttileno = 0; - - addr = 0; /* predefined: auto */ - - if(sscanf(token, "a=%d", &addr) == 1) { - /* Specified */ - if ((addr != 0) && (addr != 2) && (addr != 4)) { - fprintf(stderr, "ERROR -> invalid addressing size a = %d\n", addr); - return 1; - } - parameters->jpwl_sens_addr = addr; - - } else if (!strcmp(token, "a")) { - /* default */ - parameters->jpwl_sens_addr = addr; /* auto for default size */ - - } else { - fprintf(stderr, "ERROR -> invalid addressing selection = %s\n", token); - return 1; - }; - - } - - /* search sensitivity size */ - if (*token == 'z') { - - static int tile = 0, tilespec = 0, lasttileno = 0; - - size = 1; /* predefined: 1 byte */ - - if(sscanf(token, "z=%d", &size) == 1) { - /* Specified */ - if ((size != 0) && (size != 1) && (size != 2)) { - fprintf(stderr, "ERROR -> invalid sensitivity size z = %d\n", size); - return 1; - } - parameters->jpwl_sens_size = size; - - } else if (!strcmp(token, "a")) { - /* default */ - parameters->jpwl_sens_size = size; /* 1 for default size */ - - } else { - fprintf(stderr, "ERROR -> invalid size selection = %s\n", token); - return 1; - }; - - } - - /* search range method */ - if (*token == 'g') { - - static int tile = 0, tilespec = 0, lasttileno = 0; - - range = 0; /* predefined: 0 (packet) */ - - if(sscanf(token, "g=%d", &range) == 1) { - /* Specified */ - if ((range < 0) || (range > 3)) { - fprintf(stderr, "ERROR -> invalid sensitivity range method g = %d\n", range); - return 1; - } - parameters->jpwl_sens_range = range; - - } else if (!strcmp(token, "g")) { - /* default */ - parameters->jpwl_sens_range = range; - - } else { - fprintf(stderr, "ERROR -> invalid range selection = %s\n", token); - return 1; - }; - - } - - /* next token or bust */ - token = strtok(NULL, ","); - }; - - - /* some info */ - fprintf(stdout, "Info: JPWL capabilities enabled\n"); - parameters->jpwl_epc_on = true; - - } - break; -#endif /* USE_JPWL */ -/* < Command line not valid\n"); - return 1; - } - } - - /* check for possible errors */ - if (parameters->cp_cinema){ - if(parameters->tcp_numlayers > 1){ - parameters->cp_rsiz = STD_RSIZ; - fprintf(stdout,"Warning: DC profiles do not allow more than one quality layer. The codestream created will not be compliant with the DC profile\n"); - } - } - - if ((parameters->cp_disto_alloc || parameters->cp_fixed_alloc || parameters->cp_fixed_quality) - && (!(parameters->cp_disto_alloc ^ parameters->cp_fixed_alloc ^ parameters->cp_fixed_quality))) { - fprintf(stderr, "Error: options -r -q and -f cannot be used together !!\n"); - return 1; - } /* mod fixed_quality */ - - /* if no rate entered, lossless by default */ - if (parameters->tcp_numlayers == 0) { - parameters->tcp_rates[0] = 0; /* MOD antonin : losslessbug */ - parameters->tcp_numlayers++; - parameters->cp_disto_alloc = 1; - } - - if((parameters->cp_tx0 > parameters->image_offset_x0) || (parameters->cp_ty0 > parameters->image_offset_y0)) { - fprintf(stderr, - "Error: Tile offset dimension is unnappropriate --> TX0(%d)<=IMG_X0(%d) TYO(%d)<=IMG_Y0(%d) \n", - parameters->cp_tx0, parameters->image_offset_x0, parameters->cp_ty0, parameters->image_offset_y0); - return 1; - } - - for (i = 0; i < parameters->numpocs; i++) { - if (parameters->POC[i].prg == -1) { - fprintf(stderr, - "Unrecognized progression order in option -P (POC n %d) [LRCP, RLCP, RPCL, PCRL, CPRL] !!\n", - i + 1); - } - } - - return 0; -} - - -/** Create the same index as j2k_create_index does, but in an int[] instead of in a file ==> easy to pass it back to Java, to transfer it etc. - @param buffer_size, increased by the length of the compressed index, in number of bytes - @return a pointer to a char[] - Syntax of the index: - one char for the version number (1): one byte because no problem with little endian, big endian etc. - one int for each of the following informations: - Image Width - Image Height - progression order - Tile width - Tile height - Nb tiles in X - Nb tiles in Y - Nb of components - Nb of layers - Nb of resolutions - - for each resolution: - Precinct width - Precinct height - - End main header position - codestream size - - For each tile: - tile number - tile start pos in codestream - tile header end position - tile end position in codestream - - For each LRCP, RLCP etc.: - packet number - tile number - layer number - resolution number - component number - precinct number - start position in the codestream - end position of this packet - */ -static char* create_index_into_byte_array(opj_codestream_info_t *cstr_info, int* buffer_size) { - int tileno, compno, layno, resno, precno, pack_nb, x, y; - char* buffer = NULL; - int buffer_pos = 0; - int prec_max = 0; - - prec_max = 0; - for (tileno = 0; tileno < cstr_info->tw * cstr_info->th; tileno++) { - for (resno = 0; resno < cstr_info->numdecompos[0] + 1; resno++) { - prec_max = int_max(prec_max,cstr_info->tile[tileno].pw[resno] * cstr_info->tile[tileno].ph[resno]); - } - } - - /* Compute the size of the index buffer, in number of bytes*/ - *buffer_size = - 1 /* version */ - + (10 /* image_w until decomposition */ - + (cstr_info->numdecompos[0]+1) * 2 /* pdx size for each tile */ - + 2 /* main_head_end + codestream_size */ - + cstr_info->tw * cstr_info->th * 4 /* tile info, without distorsion info */ - + cstr_info->tw*cstr_info->th * cstr_info->numlayers * (cstr_info->numdecompos[0] + 1) * cstr_info->numcomps * prec_max *8 - ) * sizeof(int); - /*printf("C: index buffer size = %d bytes\n", *buffer_size);*/ - buffer = (char*) opj_malloc(*buffer_size); - - if (!buffer) { - /*opj_event_msg(j2k->cinfo, EVT_ERROR, "failed to allocate index buffer for writing %d int\n", *buffer_size);*/ - fprintf(stderr, "failed to allocate index buffer for writing %d int\n", *buffer_size); - return 0; - } - - buffer[0] = 1; /* Version stored on a byte*/ - buffer++; - /* Remaining informations are stored on a int.*/ - ((int*)buffer)[buffer_pos++] = cstr_info->image_w; - ((int*)buffer)[buffer_pos++] = cstr_info->image_h; - ((int*)buffer)[buffer_pos++] = cstr_info->prog; - ((int*)buffer)[buffer_pos++] = cstr_info->tile_x; - ((int*)buffer)[buffer_pos++] = cstr_info->tile_y; - ((int*)buffer)[buffer_pos++] = cstr_info->tw; - ((int*)buffer)[buffer_pos++] = cstr_info->th; - ((int*)buffer)[buffer_pos++] = cstr_info->numcomps; - ((int*)buffer)[buffer_pos++] = cstr_info->numlayers; - ((int*)buffer)[buffer_pos++] = cstr_info->numdecompos[0]; - - for (resno = cstr_info->numdecompos[0]; resno >= 0; resno--) { - /* based on tile 0 */ - ((int*)buffer)[buffer_pos++] = (1 << cstr_info->tile[0].pdx[resno]); - ((int*)buffer)[buffer_pos++] = (1 << cstr_info->tile[0].pdx[resno]); - } - ((int*)buffer)[buffer_pos++] = cstr_info->main_head_end; - ((int*)buffer)[buffer_pos++] = cstr_info->codestream_size; - - for (tileno = 0; tileno < cstr_info->tw * cstr_info->th; tileno++) { - ((int*)buffer)[buffer_pos++] = cstr_info->tile[tileno].tileno; - ((int*)buffer)[buffer_pos++] = cstr_info->tile[tileno].start_pos; - ((int*)buffer)[buffer_pos++] = cstr_info->tile[tileno].end_header; - ((int*)buffer)[buffer_pos++] = cstr_info->tile[tileno].end_pos; - } - - for (tileno = 0; tileno < cstr_info->tw * cstr_info->th; tileno++) { - int start_pos, end_pos; - int max_numdecompos = 0; - pack_nb = 0; - - for (compno = 0; compno < cstr_info->numcomps; compno++) { - if (max_numdecompos < cstr_info->numdecompos[compno]) - max_numdecompos = cstr_info->numdecompos[compno]; - } - - if (cstr_info->prog == LRCP) { /* LRCP */ - - for (layno = 0; layno < cstr_info->numlayers; layno++) { - for (resno = 0; resno < max_numdecompos + 1; resno++) { - for (compno = 0; compno < cstr_info->numcomps; compno++) { - int prec_max; - if (resno > cstr_info->numdecompos[compno]) - break; - prec_max = cstr_info->tile[tileno].pw[resno] * cstr_info->tile[tileno].ph[resno]; - for (precno = 0; precno < prec_max; precno++) { - start_pos = cstr_info->tile[tileno].packet[pack_nb].start_pos; - end_pos = cstr_info->tile[tileno].packet[pack_nb].end_pos; - ((int*)buffer)[buffer_pos++] = pack_nb; - ((int*)buffer)[buffer_pos++] = tileno; - ((int*)buffer)[buffer_pos++] = layno; - ((int*)buffer)[buffer_pos++] = resno; - ((int*)buffer)[buffer_pos++] = compno; - ((int*)buffer)[buffer_pos++] = precno; - ((int*)buffer)[buffer_pos++] = start_pos; - ((int*)buffer)[buffer_pos++] = end_pos; - pack_nb++; - } - } - } - } - } /* LRCP */ - else if (cstr_info->prog == RLCP) { /* RLCP */ - - for (resno = 0; resno < max_numdecompos + 1; resno++) { - for (layno = 0; layno < cstr_info->numlayers; layno++) { - for (compno = 0; compno < cstr_info->numcomps; compno++) { - int prec_max; - if (resno > cstr_info->numdecompos[compno]) - break; - prec_max = cstr_info->tile[tileno].pw[resno] * cstr_info->tile[tileno].ph[resno]; - for (precno = 0; precno < prec_max; precno++) { - start_pos = cstr_info->tile[tileno].packet[pack_nb].start_pos; - end_pos = cstr_info->tile[tileno].packet[pack_nb].end_pos; - ((int*)buffer)[buffer_pos++] = pack_nb; - ((int*)buffer)[buffer_pos++] = tileno; - ((int*)buffer)[buffer_pos++] = resno; - ((int*)buffer)[buffer_pos++] = layno; - ((int*)buffer)[buffer_pos++] = compno; - ((int*)buffer)[buffer_pos++] = precno; - ((int*)buffer)[buffer_pos++] = start_pos; - ((int*)buffer)[buffer_pos++] = end_pos; - pack_nb++; - } - } - } - } - } /* RLCP */ - else if (cstr_info->prog == RPCL) { /* RPCL */ - - for (resno = 0; resno < max_numdecompos + 1; resno++) { - /* I suppose components have same XRsiz, YRsiz */ - int x0 = cstr_info->tile_Ox + tileno - (int)floor( (float)tileno/(float)cstr_info->tw ) * cstr_info->tw * cstr_info->tile_x; - int y0 = cstr_info->tile_Ox + (int)floor( (float)tileno/(float)cstr_info->tw ) * cstr_info->tile_y; - int x1 = x0 + cstr_info->tile_x; - int y1 = y0 + cstr_info->tile_y; - for (compno = 0; compno < cstr_info->numcomps; compno++) { - int prec_max; - if (resno > cstr_info->numdecompos[compno]) - break; - prec_max = cstr_info->tile[tileno].pw[resno] * cstr_info->tile[tileno].ph[resno]; - for (precno = 0; precno < prec_max; precno++) { - int pcnx = cstr_info->tile[tileno].pw[resno]; - int pcx = (int) pow( 2, cstr_info->tile[tileno].pdx[resno] + cstr_info->numdecompos[compno] - resno ); - int pcy = (int) pow( 2, cstr_info->tile[tileno].pdy[resno] + cstr_info->numdecompos[compno] - resno ); - int precno_x = precno - (int) floor( (float)precno/(float)pcnx ) * pcnx; - int precno_y = (int) floor( (float)precno/(float)pcnx ); - for(y = y0; y < y1; y++) { - if (precno_y*pcy == y ) { - for (x = x0; x < x1; x++) { - if (precno_x*pcx == x ) { - for (layno = 0; layno < cstr_info->numlayers; layno++) { - start_pos = cstr_info->tile[tileno].packet[pack_nb].start_pos; - end_pos = cstr_info->tile[tileno].packet[pack_nb].end_pos; - ((int*)buffer)[buffer_pos++] = pack_nb; - ((int*)buffer)[buffer_pos++] = tileno; - ((int*)buffer)[buffer_pos++] = resno; - ((int*)buffer)[buffer_pos++] = precno; - ((int*)buffer)[buffer_pos++] = compno; - ((int*)buffer)[buffer_pos++] = layno; - ((int*)buffer)[buffer_pos++] = start_pos; - ((int*)buffer)[buffer_pos++] = end_pos; - pack_nb++; - } - } - }/* x = x0..x1 */ - } - } /* y = y0..y1 */ - } /* precno */ - } /* compno */ - } /* resno */ - } /* RPCL */ - else if (cstr_info->prog == PCRL) { /* PCRL */ - /* I suppose components have same XRsiz, YRsiz */ - int x0 = cstr_info->tile_Ox + tileno - (int)floor( (float)tileno/(float)cstr_info->tw ) * cstr_info->tw * cstr_info->tile_x; - int y0 = cstr_info->tile_Ox + (int)floor( (float)tileno/(float)cstr_info->tw ) * cstr_info->tile_y; - int x1 = x0 + cstr_info->tile_x; - int y1 = y0 + cstr_info->tile_y; - - for (compno = 0; compno < cstr_info->numcomps; compno++) { - for (resno = 0; resno < max_numdecompos + 1; resno++) { - int prec_max; - if (resno > cstr_info->numdecompos[compno]) - break; - prec_max = cstr_info->tile[tileno].pw[resno] * cstr_info->tile[tileno].ph[resno]; - for (precno = 0; precno < prec_max; precno++) { - int pcnx = cstr_info->tile[tileno].pw[resno]; - int pcx = (int) pow( 2, cstr_info->tile[tileno].pdx[resno] + cstr_info->numdecompos[compno] - resno ); - int pcy = (int) pow( 2, cstr_info->tile[tileno].pdy[resno] + cstr_info->numdecompos[compno] - resno ); - int precno_x = precno - (int) floor( (float)precno/(float)pcnx ) * pcnx; - int precno_y = (int) floor( (float)precno/(float)pcnx ); - for(y = y0; y < y1; y++) { - if (precno_y*pcy == y ) { - for (x = x0; x < x1; x++) { - if (precno_x*pcx == x ) { - for (layno = 0; layno < cstr_info->numlayers; layno++) { - start_pos = cstr_info->tile[tileno].packet[pack_nb].start_pos; - end_pos = cstr_info->tile[tileno].packet[pack_nb].end_pos; - ((int*)buffer)[buffer_pos++] = pack_nb; - ((int*)buffer)[buffer_pos++] = tileno; - ((int*)buffer)[buffer_pos++] = precno; - ((int*)buffer)[buffer_pos++] = compno; - ((int*)buffer)[buffer_pos++] = resno; - ((int*)buffer)[buffer_pos++] = layno; - ((int*)buffer)[buffer_pos++] = start_pos; - ((int*)buffer)[buffer_pos++] = end_pos; - pack_nb++; - } - } - }/* x = x0..x1 */ - } - } /* y = y0..y1 */ - } /* precno */ - } /* resno */ - } /* compno */ - } /* PCRL */ - else { /* CPRL */ - - for (compno = 0; compno < cstr_info->numcomps; compno++) { - /* I suppose components have same XRsiz, YRsiz */ - int x0 = cstr_info->tile_Ox + tileno - (int)floor( (float)tileno/(float)cstr_info->tw ) * cstr_info->tw * cstr_info->tile_x; - int y0 = cstr_info->tile_Ox + (int)floor( (float)tileno/(float)cstr_info->tw ) * cstr_info->tile_y; - int x1 = x0 + cstr_info->tile_x; - int y1 = y0 + cstr_info->tile_y; - - for (resno = 0; resno < max_numdecompos + 1; resno++) { - int prec_max; - if (resno > cstr_info->numdecompos[compno]) - break; - prec_max = cstr_info->tile[tileno].pw[resno] * cstr_info->tile[tileno].ph[resno]; - for (precno = 0; precno < prec_max; precno++) { - int pcnx = cstr_info->tile[tileno].pw[resno]; - int pcx = (int) pow( 2, cstr_info->tile[tileno].pdx[resno] + cstr_info->numdecompos[compno] - resno ); - int pcy = (int) pow( 2, cstr_info->tile[tileno].pdy[resno] + cstr_info->numdecompos[compno] - resno ); - int precno_x = precno - (int) floor( (float)precno/(float)pcnx ) * pcnx; - int precno_y = (int) floor( (float)precno/(float)pcnx ); - for(y = y0; y < y1; y++) { - if (precno_y*pcy == y ) { - for (x = x0; x < x1; x++) { - if (precno_x*pcx == x ) { - for (layno = 0; layno < cstr_info->numlayers; layno++) { - start_pos = cstr_info->tile[tileno].packet[pack_nb].start_pos; - end_pos = cstr_info->tile[tileno].packet[pack_nb].end_pos; - ((int*)buffer)[buffer_pos++] = pack_nb; - ((int*)buffer)[buffer_pos++] = tileno; - ((int*)buffer)[buffer_pos++] = compno; - ((int*)buffer)[buffer_pos++] = precno; - ((int*)buffer)[buffer_pos++] = resno; - ((int*)buffer)[buffer_pos++] = layno; - ((int*)buffer)[buffer_pos++] = start_pos; - ((int*)buffer)[buffer_pos++] = end_pos; - pack_nb++; - } - } - }/* x = x0..x1 */ - } - } /* y = y0..y1 */ - } /* precno */ - } /* resno */ - } /* compno */ - } /* CPRL */ - } /* tileno */ - - if (buffer_pos > *buffer_size) { - /*opj_event_msg(j2k->cinfo, EVT_ERROR, "index creation: buffer_pos (%d) > buffer_size (%d)!\n", buffer_pos, *buffer_size);*/ - fprintf(stderr, "index creation: buffer_pos (%d) > buffer_size (%d)!\n", buffer_pos, *buffer_size); - return 0; - } - - return --buffer; -} - - - - -/* -------------------------------------------------------------------------- - ------------ Get the image byte[] from the Java object -------------------*/ - -static opj_image_t* loadImage(opj_cparameters_t *parameters, JNIEnv *env, jobject obj, jclass cls) { - int i,max,shift,w,h,depth; - opj_image_t * img = NULL; - int compno, numcomps; - opj_image_t * image = NULL; - opj_image_comp_t *comp; - opj_image_cmptparm_t cmptparm[3]; /* maximum of 3 components */ - OPJ_COLOR_SPACE color_space; - jfieldID fid; - jint ji; - jbyteArray jba; - jshortArray jsa; - jintArray jia; - int len; - jbyte *jbBody; - jshort *jsBody; - jint *jiBody; - jboolean isCopy; - - /* Image width, height and depth*/ - fid = (*env)->GetFieldID(env, cls,"width", "I"); - ji = (*env)->GetIntField(env, obj, fid); - w = ji; - - fid = (*env)->GetFieldID(env, cls,"height", "I"); - ji = (*env)->GetIntField(env, obj, fid); - h = ji; - - fid = (*env)->GetFieldID(env, cls,"depth", "I"); - ji = (*env)->GetIntField(env, obj, fid); - depth = ji; - - /* Read the image*/ - if (depth <=16) { - numcomps = 1; - color_space = CLRSPC_GRAY; - } else { - numcomps = 3; - color_space = CLRSPC_SRGB; - } - memset(&cmptparm[0], 0, numcomps * sizeof(opj_image_cmptparm_t)); - - if (numcomps == 1) { - cmptparm[0].x0 = parameters->image_offset_x0; - cmptparm[0].y0 = parameters->image_offset_y0; - cmptparm[0].w = !cmptparm[0].x0 ? (w - 1) * parameters->subsampling_dx + 1 : cmptparm[0].x0 + (w - 1) * parameters->subsampling_dx + 1; - cmptparm[0].h = !cmptparm[0].y0 ? (h - 1) * parameters->subsampling_dy + 1 : cmptparm[0].y0 + (h - 1) * parameters->subsampling_dy + 1; - /* Java types are always signed but we use them as unsigned types (shift of the negative part of - the pixels of the images in Telemis before entering the encoder).*/ - cmptparm[0].sgnd = 0; - if (depth<=16) - cmptparm[0].prec=depth; - else - cmptparm[0].prec = 8; - cmptparm[0].bpp = cmptparm[0].prec; - cmptparm[0].dx = parameters->subsampling_dx; - cmptparm[0].dy = parameters->subsampling_dy; - /*printf("C: component 0 initialised: x0=%d, y0=%d, w=%d, h=%d, sgnd=%d, bpp=%d, dx=%d, dy=%d, color_space=%d\n", cmptparm[0].x0, cmptparm[0].y0, cmptparm[0].w, - cmptparm[0].h, cmptparm[0].sgnd, cmptparm[0].bpp, cmptparm[0].dx, cmptparm[0].dy, color_space);*/ - } else { - for(i = 0; i < numcomps; i++) { - cmptparm[i].prec = 8; - cmptparm[i].bpp = 8; - cmptparm[i].sgnd = 0; - cmptparm[i].dx = parameters->subsampling_dx; - cmptparm[i].dy = parameters->subsampling_dy; - cmptparm[i].w = w; - cmptparm[i].h = h; - } - } - - /* create the image */ - image = opj_image_create(numcomps, &cmptparm[0], color_space); - - if (!image) - return NULL; - - if (depth <=16) { - image->numcomps=1; - } else { - image->numcomps = 3; - } - - /* set image offset and reference grid */ - image->x0 = cmptparm[0].x0; - image->y0 = cmptparm[0].x0; - image->x1 = cmptparm[0].w; - image->y1 = cmptparm[0].h; - - /* set image data */ - for (compno=0; compnocomps[compno]; - max = -100000; - if (depth == 8) { - fid = (*env)->GetFieldID(env, cls,"image8", "[B"); /* byteArray []*/ - jba = (*env)->GetObjectField(env, obj, fid); - len = (*env)->GetArrayLength(env, jba); - - jbBody = (*env)->GetPrimitiveArrayCritical(env, jba, &isCopy); - /*printf("C: before transferring 8 bpp image\n");*/ - if (comp->sgnd) { - for(i=0; i< len;i++) { - comp->data[i] = (char) jbBody[i]; - if (comp->data[i] > max) max = comp->data[i]; - } - } else { - for(i=0; i< len;i++) { - comp->data[i] = (unsigned char) jbBody[i]; - if (comp->data[i] > max) max = comp->data[i]; - } - } - (*env)->ReleasePrimitiveArrayCritical(env, jba, jbBody, 0); - } else if(depth == 16) { - fid = (*env)->GetFieldID(env, cls,"image16", "[S"); /* shortArray []*/ - jsa = (*env)->GetObjectField(env, obj, fid); - len = (*env)->GetArrayLength(env, jsa); - - jsBody = (*env)->GetPrimitiveArrayCritical(env, jsa, &isCopy); - /*printf("C: before transferring 16 bpp image\n");*/ - if (comp->sgnd) { /* Special behaviour to deal with signed elements ??*/ - comp->data[i] = (short) jsBody[i]; - for(i=0; i< len;i++) { - if (comp->data[i] > max) max = comp->data[i]; - } - } else { - for(i=0; i< len;i++) { - comp->data[i] = (unsigned short) jsBody[i]; - if (comp->data[i] > max) max = comp->data[i]; - } - } - (*env)->ReleasePrimitiveArrayCritical(env, jsa, jsBody, 0); - } else if (depth == 24) { - fid = (*env)->GetFieldID(env, cls,"image24", "[I"); /* intArray []*/ - jia = (*env)->GetObjectField(env, obj, fid); - len = (*env)->GetArrayLength(env, jia); - shift = compno*8; - - jiBody = (*env)->GetPrimitiveArrayCritical(env, jia, &isCopy); - /*printf("C: before transferring 24 bpp image (component %d, signed = %d)\n", compno, comp->sgnd);*/ - if (comp->sgnd) { /* Special behaviour to deal with signed elements ?? XXXXX*/ - for(i=0; i< len;i++) { - comp->data[i] = ( ((int) jiBody[i]) & (0xFF << shift) ) >> shift; - if (comp->data[i] > max) max = comp->data[i]; - } - } else { - for(i=0; i< len;i++) { - comp->data[i] = ( ((unsigned int) jiBody[i]) & (0xFF << shift) ) >> shift; - if (comp->data[i] > max) max = comp->data[i]; - } - } - (*env)->ReleasePrimitiveArrayCritical(env, jia, jiBody, 0); - } - comp->bpp = int_floorlog2(max)+1; - comp->prec = comp->bpp; - /*printf("C: component %d: max %d, real bpp = %d\n", compno, max, comp->bpp);*/ - } - return image; -} - - -/* -------------------------------------------------------------------------- - -------------------- MAIN METHOD, CALLED BY JAVA -----------------------*/ -JNIEXPORT jlong JNICALL Java_org_openJpeg_OpenJPEGJavaEncoder_internalEncodeImageToJ2K(JNIEnv *env, jobject obj, jobjectArray javaParameters) { - int argc; /* To simulate the command line parameters (taken from the javaParameters variable) and be able to re-use the */ - char **argv; /* 'parse_cmdline_decoder' method taken from the j2k_to_image project */ - opj_bool bSuccess; - opj_cparameters_t parameters; /* compression parameters */ - img_fol_t img_fol; - opj_event_mgr_t event_mgr; /* event manager */ - opj_image_t *image = NULL; - int i,j,num_images; - int imageno; - opj_codestream_info_t cstr_info; /* Codestream information structure */ - char indexfilename[OPJ_PATH_LEN]; /* index file name */ - - char* compressed_index = NULL; - int compressed_index_size=-1; - /* ==> Access variables to the Java member variables*/ - jsize arraySize; - jclass cls; - jobject object; - jboolean isCopy; - jfieldID fid; - jbyteArray jba; - jbyte *jbBody; - callback_variables_t msgErrorCallback_vars; - /* <== access variable to the Java member variables.*/ - - /* For the encoding and storage into the file*/ - opj_cinfo_t* cinfo; - int codestream_length; - opj_cio_t *cio = NULL; - FILE *f = NULL; - - /* JNI reference to the calling class*/ - cls = (*env)->GetObjectClass(env, obj); - - /* Pointers to be able to call a Java method for all the info and error messages*/ - msgErrorCallback_vars.env = env; - msgErrorCallback_vars.jobj = &obj; - msgErrorCallback_vars.message_mid = (*env)->GetMethodID(env, cls, "logMessage", "(Ljava/lang/String;)V"); - msgErrorCallback_vars.error_mid = (*env)->GetMethodID(env, cls, "logError", "(Ljava/lang/String;)V"); - - arraySize = (*env)->GetArrayLength(env, javaParameters); - argc = (int) arraySize +1; - argv = opj_malloc(argc*sizeof(char*)); - argv[0] = "ProgramName.exe"; /* The program name: useless*/ - j=0; - for (i=1; iGetObjectArrayElement(env, javaParameters, i-1); - argv[i] = (char*)(*env)->GetStringUTFChars(env, object, &isCopy); - } - - /*printf("C: "); - for (i=0; iReleaseStringUTFChars(env, (*env)->GetObjectArrayElement(env, javaParameters, i-1), argv[i]); - return -1; - } - - /* Release the Java arguments array*/ - for (i=1; iReleaseStringUTFChars(env, (*env)->GetObjectArrayElement(env, javaParameters, i-1), argv[i]); - - if (parameters.cp_cinema){ - cinema_parameters(¶meters); - } - - - /* Create comment for codestream */ - if(parameters.cp_comment == NULL) { - const char comment[] = "Created by JavaOpenJPEG version "; - const size_t clen = strlen(comment); - const char *version = opj_version(); -/* UniPG>> */ -#ifdef USE_JPWL - parameters.cp_comment = (char*)opj_malloc(clen+strlen(version)+11); - sprintf(parameters.cp_comment,"%s%s with JPWL", comment, version); -#else - parameters.cp_comment = (char*)opj_malloc(clen+strlen(version)+1); - sprintf(parameters.cp_comment,"%s%s", comment, version); -#endif -/* <numcomps == 3 ? 1 : 0; - - if(parameters.cp_cinema){ - cinema_setup_encoder(¶meters,image,&img_fol); - } - - /* encode the destination image */ - /* ---------------------------- */ - /* get a J2K compressor handle */ - if (parameters.cod_format == J2K_CFMT) { /* J2K format output */ - cinfo = opj_create_compress(CODEC_J2K); - } else { /* JP2 format output */ - cinfo = opj_create_compress(CODEC_JP2); - } - /* catch events using our callbacks and give a local context */ - opj_set_event_mgr((opj_common_ptr)cinfo, &event_mgr, &msgErrorCallback_vars); - - /* setup the encoder parameters using the current image and user parameters */ - opj_setup_encoder(cinfo, ¶meters, image); - - /* open a byte stream for writing */ - /* allocate memory for all tiles */ - cio = opj_cio_open((opj_common_ptr)cinfo, NULL, 0); - - /* encode the image */ - bSuccess = opj_encode_with_info(cinfo, cio, image, &cstr_info); - printf("C: after opj_encode_with_info\n"); - if (!bSuccess) { - opj_cio_close(cio); - fprintf(stderr, "failed to encode image\n"); - return -1; - } - codestream_length = cio_tell(cio); - - /* write the index on disk, if needed (-x 'filename') */ - if (*indexfilename) { - bSuccess = write_index_file(&cstr_info, indexfilename); - if (bSuccess) { - fprintf(stderr, "Failed to output index file into [%s]\n", indexfilename); - } - } - - compressed_index = create_index_into_byte_array(&cstr_info, &compressed_index_size); - /* Allocates the Java compressedIndex byte[] and sends this index into the Java object */ - fid = (*env)->GetFieldID(env, cls,"compressedIndex", "[B"); - jba = (*env)->NewByteArray(env, compressed_index_size+1); - jbBody = (*env)->GetPrimitiveArrayCritical(env, jba, 0); - memcpy(jbBody, compressed_index, compressed_index_size); - (*env)->ReleasePrimitiveArrayCritical(env, jba, jbBody, 0); - (*env)->SetObjectField(env, obj, fid, jba); - opj_free(compressed_index); - - /* write the generated codestream to disk ? */ - if (parameters.outfile[0]!='\0') { - f = fopen(parameters.outfile, "wb"); - if (!f) { - fprintf(stderr, "failed to open [%s] for writing\n", parameters.outfile); - return -1; - } - fwrite(cio->buffer, 1, codestream_length, f); - fclose(f); - fprintf(stdout,"Generated outfile [%s]\n",parameters.outfile); - } - - /* Write the generated codestream to the Java pre-allocated compressedStream byte[] */ - fid = (*env)->GetFieldID(env, cls,"compressedStream", "[B"); - jba = (*env)->GetObjectField(env, obj, fid); - jbBody = (*env)->GetPrimitiveArrayCritical(env, jba, 0); - memcpy(jbBody, cio->buffer, codestream_length); - (*env)->ReleasePrimitiveArrayCritical(env, jba, jbBody, 0); - - /* close and free the byte stream */ - opj_cio_close(cio); - - /* free remaining compression structures */ - opj_destroy_compress(cinfo); - opj_destroy_cstr_info(&cstr_info); - - /* free image data */ - opj_image_destroy(image); - } - - /* free user parameters structure */ - if(parameters.cp_comment) opj_free(parameters.cp_comment); - if(parameters.cp_matrice) opj_free(parameters.cp_matrice); - - return codestream_length; -} +/* + * Copyright (c) 2002-2014, Universite catholique de Louvain (UCL), Belgium + * Copyright (c) 2002-2014, Professor Benoit Macq + * Copyright (c) 2001-2003, David Janssens + * Copyright (c) 2002-2003, Yannick Verschueren + * Copyright (c) 2003-2007, Francois-Olivier Devaux + * Copyright (c) 2003-2014, Antonin Descampe + * Copyright (c) 2005, Herve Drolon, FreeImage Team + * Copyright (c) 2006-2007, Parvatha Elangovan + * Copyright (c) 2007, Patrick Piscaglia (Telemis) + * All rights reserved. + * + * Redistribution and use in source and binary forms, with or without + * modification, are permitted provided that the following conditions + * are met: + * 1. Redistributions of source code must retain the above copyright + * notice, this list of conditions and the following disclaimer. + * 2. Redistributions in binary form must reproduce the above copyright + * notice, this list of conditions and the following disclaimer in the + * documentation and/or other materials provided with the distribution. + * + * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS' + * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE + * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE + * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE + * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR + * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF + * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS + * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN + * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) + * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE + * POSSIBILITY OF SUCH DAMAGE. + */ +#include +#include +#include +#include +#include + +#include "openjpeg.h" +#include "opj_includes.h" +#include "opj_getopt.h" +#include "convert.h" +#include "index.h" +#include "dirent.h" +#include "org_openJpeg_OpenJPEGJavaEncoder.h" + +#ifndef _WIN32 +#define stricmp strcasecmp +#define strnicmp strncasecmp +#endif + +#include "format_defs.h" + +#define CINEMA_24_CS 1302083 /*Codestream length for 24fps*/ +#define CINEMA_48_CS 651041 /*Codestream length for 48fps*/ +#define COMP_24_CS 1041666 /*Maximum size per color component for 2K & 4K @ 24fps*/ +#define COMP_48_CS 520833 /*Maximum size per color component for 2K @ 48fps*/ + +extern int get_file_format(char *filename); +extern void error_callback(const char *msg, void *client_data); +extern void warning_callback(const char *msg, void *client_data); +extern void info_callback(const char *msg, void *client_data); + +typedef struct callback_variables { + JNIEnv *env; + /** 'jclass' object used to call a Java method from the C */ + jobject *jobj; + /** 'jclass' object used to call a Java method from the C */ + jmethodID message_mid; + jmethodID error_mid; +} callback_variables_t; + +typedef struct dircnt { + /** Buffer for holding images read from Directory*/ + char *filename_buf; + /** Pointer to the buffer*/ + char **filename; +} dircnt_t; + +typedef struct img_folder { + /** The directory path of the folder containing input images*/ + char *imgdirpath; + /** Output format*/ + char *out_format; + /** Enable option*/ + char set_imgdir; + /** Enable Cod Format for output*/ + char set_out_format; + /** User specified rate stored in case of cinema option*/ + float *rates; +} img_fol_t; + +static void encode_help_display() +{ + fprintf(stdout, "HELP\n----\n\n"); + fprintf(stdout, "- the -h option displays this help information on screen\n\n"); + + /* UniPG>> */ + fprintf(stdout, "List of parameters for the JPEG 2000 " +#ifdef USE_JPWL + "+ JPWL " +#endif /* USE_JPWL */ + "encoder:\n"); + /* <> */ +#ifdef USE_JPWL + fprintf(stdout, " * No JPWL protection\n"); +#endif /* USE_JPWL */ + /* < \n"); + fprintf(stdout, + " Currently accepts PGM, PPM, PNM, PGX, BMP, TIF, RAW and TGA formats\n"); + fprintf(stdout, "\n"); + fprintf(stdout, + "-i : source file (-i source.pnm also *.pgm, *.ppm, *.bmp, *.tif, *.raw, *.tga) \n"); + fprintf(stdout, " When using this option -o must be used\n"); + fprintf(stdout, "\n"); + fprintf(stdout, "-o : destination file (-o dest.j2k or .jp2) \n"); + fprintf(stdout, "\n"); + fprintf(stdout, "Optional Parameters:\n"); + fprintf(stdout, "\n"); + fprintf(stdout, "-h : display the help information \n "); + fprintf(stdout, "\n"); + fprintf(stdout, + "-cinema2K : Digital Cinema 2K profile compliant codestream for 2K resolution.(-cinema2k 24 or 48) \n"); + fprintf(stdout, + " Need to specify the frames per second for a 2K resolution. Only 24 or 48 fps is allowed\n"); + fprintf(stdout, "\n"); + fprintf(stdout, + "-cinema4K : Digital Cinema 4K profile compliant codestream for 4K resolution \n"); + fprintf(stdout, " Frames per second not required. Default value is 24fps\n"); + fprintf(stdout, "\n"); + fprintf(stdout, + "-r : different compression ratios for successive layers (-r 20,10,5)\n "); + fprintf(stdout, + " - The rate specified for each quality level is the desired \n"); + fprintf(stdout, " compression factor.\n"); + fprintf(stdout, " Example: -r 20,10,1 means quality 1: compress 20x, \n"); + fprintf(stdout, + " quality 2: compress 10x and quality 3: compress lossless\n"); + fprintf(stdout, "\n"); + fprintf(stdout, + " (options -r and -q cannot be used together)\n "); + fprintf(stdout, "\n"); + + fprintf(stdout, + "-q : different psnr for successive layers (-q 30,40,50) \n "); + + fprintf(stdout, + " (options -r and -q cannot be used together)\n "); + + fprintf(stdout, "\n"); + fprintf(stdout, "-n : number of resolutions (-n 3) \n"); + fprintf(stdout, "\n"); + fprintf(stdout, "-b : size of code block (-b 32,32) \n"); + fprintf(stdout, "\n"); + fprintf(stdout, "-c : size of precinct (-c 128,128) \n"); + fprintf(stdout, "\n"); + fprintf(stdout, "-t : size of tile (-t 512,512) \n"); + fprintf(stdout, "\n"); + fprintf(stdout, + "-p : progression order (-p LRCP) [LRCP, RLCP, RPCL, PCRL, CPRL] \n"); + fprintf(stdout, "\n"); + fprintf(stdout, "-s : subsampling factor (-s 2,2) [-s X,Y] \n"); + fprintf(stdout, " Remark: subsampling bigger than 2 can produce error\n"); + fprintf(stdout, "\n"); + fprintf(stdout, + "-POC : Progression order change (-POC T1=0,0,1,5,3,CPRL/T1=5,0,1,6,3,CPRL) \n"); + fprintf(stdout, " Example: T1=0,0,1,5,3,CPRL \n"); + fprintf(stdout, + " : Ttilenumber=Resolution num start,Component num start,Layer num end,Resolution num end,Component num end,Progression order\n"); + fprintf(stdout, "\n"); + fprintf(stdout, "-SOP : write SOP marker before each packet \n"); + fprintf(stdout, "\n"); + fprintf(stdout, "-EPH : write EPH marker after each header packet \n"); + fprintf(stdout, "\n"); + fprintf(stdout, + "-M : mode switch (-M 3) [1=BYPASS(LAZY) 2=RESET 4=RESTART(TERMALL)\n"); + fprintf(stdout, + " 8=VSC 16=ERTERM(SEGTERM) 32=SEGMARK(SEGSYM)] \n"); + fprintf(stdout, + " Indicate multiple modes by adding their values. \n"); + fprintf(stdout, + " ex: RESTART(4) + RESET(2) + SEGMARK(32) = -M 38\n"); + fprintf(stdout, "\n"); + fprintf(stdout, + "-TP : divide packets of every tile into tile-parts (-TP R) [R, L, C]\n"); + fprintf(stdout, "\n"); + fprintf(stdout, + "-x : create an index file *.Idx (-x index_name.Idx) \n"); + fprintf(stdout, "\n"); + fprintf(stdout, + "-ROI : c=%%d,U=%%d : quantization indices upshifted \n"); + fprintf(stdout, " for component c=%%d [%%d = 0,1,2]\n"); + fprintf(stdout, + " with a value of U=%%d [0 <= %%d <= 37] (i.e. -ROI c=0,U=25) \n"); + fprintf(stdout, "\n"); + fprintf(stdout, + "-d : offset of the origin of the image (-d 150,300) \n"); + fprintf(stdout, "\n"); + fprintf(stdout, + "-T : offset of the origin of the tiles (-T 100,75) \n"); + fprintf(stdout, "\n"); + fprintf(stdout, "-I : use the irreversible DWT 9-7 (-I) \n"); + fprintf(stdout, "\n"); + fprintf(stdout, + "-jpip : write jpip codestream index box in JP2 output file\n"); + fprintf(stdout, " NOTICE: currently supports only RPCL order\n"); + fprintf(stdout, "\n"); + /* UniPG>> */ +#ifdef USE_JPWL + fprintf(stdout, + "-W : adoption of JPWL (Part 11) capabilities (-W params)\n"); + fprintf(stdout, + " The parameters can be written and repeated in any order:\n"); + fprintf(stdout, + " [h<=type>,s<=method>,a=,...\n"); + fprintf(stdout, + " ...,z=,g=,p<=type>]\n"); + fprintf(stdout, "\n"); + fprintf(stdout, + " h selects the header error protection (EPB): 'type' can be\n"); + fprintf(stdout, + " [0=none 1,absent=predefined 16=CRC-16 32=CRC-32 37-128=RS]\n"); + fprintf(stdout, + " if 'tilepart' is absent, it is for main and tile headers\n"); + fprintf(stdout, + " if 'tilepart' is present, it applies from that tile\n"); + fprintf(stdout, + " onwards, up to the next h<> spec, or to the last tilepart\n"); + fprintf(stdout, " in the codestream (max. %d specs)\n", + JPWL_MAX_NO_TILESPECS); + fprintf(stdout, "\n"); + fprintf(stdout, + " p selects the packet error protection (EEP/UEP with EPBs)\n"); + fprintf(stdout, " to be applied to raw data: 'type' can be\n"); + fprintf(stdout, + " [0=none 1,absent=predefined 16=CRC-16 32=CRC-32 37-128=RS]\n"); + fprintf(stdout, + " if 'tilepart:pack' is absent, it is from tile 0, packet 0\n"); + fprintf(stdout, + " if 'tilepart:pack' is present, it applies from that tile\n"); + fprintf(stdout, + " and that packet onwards, up to the next packet spec\n"); + fprintf(stdout, + " or to the last packet in the last tilepart in the stream\n"); + fprintf(stdout, " (max. %d specs)\n", + JPWL_MAX_NO_PACKSPECS); + fprintf(stdout, "\n"); + fprintf(stdout, + " s enables sensitivity data insertion (ESD): 'method' can be\n"); + fprintf(stdout, + " [-1=NO ESD 0=RELATIVE ERROR 1=MSE 2=MSE REDUCTION 3=PSNR\n"); + fprintf(stdout, + " 4=PSNR INCREMENT 5=MAXERR 6=TSE 7=RESERVED]\n"); + fprintf(stdout, + " if 'tilepart' is absent, it is for main header only\n"); + fprintf(stdout, + " if 'tilepart' is present, it applies from that tile\n"); + fprintf(stdout, + " onwards, up to the next s<> spec, or to the last tilepart\n"); + fprintf(stdout, " in the codestream (max. %d specs)\n", + JPWL_MAX_NO_TILESPECS); + fprintf(stdout, "\n"); + fprintf(stdout, + " g determines the addressing mode: can be\n"); + fprintf(stdout, " [0=PACKET 1=BYTE RANGE 2=PACKET RANGE]\n"); + fprintf(stdout, "\n"); + fprintf(stdout, + " a determines the size of data addressing: can be\n"); + fprintf(stdout, + " 2/4 bytes (small/large codestreams). If not set, auto-mode\n"); + fprintf(stdout, "\n"); + fprintf(stdout, + " z determines the size of sensitivity values: can be\n"); + fprintf(stdout, + " 1/2 bytes, for the transformed pseudo-floating point value\n"); + fprintf(stdout, "\n"); + fprintf(stdout, " ex.:\n"); + fprintf(stdout, + " h,h0=64,h3=16,h5=32,p0=78,p0:24=56,p1,p3:0=0,p3:20=32,s=0,\n"); + fprintf(stdout, " s0=6,s3=-1,a=0,g=1,z=1\n"); + fprintf(stdout, " means\n"); + fprintf(stdout, + " predefined EPB in MH, rs(64,32) from TPH 0 to TPH 2,\n"); + fprintf(stdout, + " CRC-16 in TPH 3 and TPH 4, CRC-32 in remaining TPHs,\n"); + fprintf(stdout, + " UEP rs(78,32) for packets 0 to 23 of tile 0,\n"); + fprintf(stdout, + " UEP rs(56,32) for packs. 24 to the last of tilepart 0,\n"); + fprintf(stdout, + " UEP rs default for packets of tilepart 1,\n"); + fprintf(stdout, + " no UEP for packets 0 to 19 of tilepart 3,\n"); + fprintf(stdout, + " UEP CRC-32 for packs. 20 of tilepart 3 to last tilepart,\n"); + fprintf(stdout, " relative sensitivity ESD for MH,\n"); + fprintf(stdout, + " TSE ESD from TPH 0 to TPH 2, byte range with automatic\n"); + fprintf(stdout, + " size of addresses and 1 byte for each sensitivity value\n"); + fprintf(stdout, "\n"); + fprintf(stdout, " ex.:\n"); + fprintf(stdout, " h,s,p\n"); + fprintf(stdout, " means\n"); + fprintf(stdout, + " default protection to headers (MH and TPHs) as well as\n"); + fprintf(stdout, " data packets, one ESD in MH\n"); + fprintf(stdout, "\n"); + fprintf(stdout, + " N.B.: use the following recommendations when specifying\n"); + fprintf(stdout, " the JPWL parameters list\n"); + fprintf(stdout, + " - when you use UEP, always pair the 'p' option with 'h'\n"); + fprintf(stdout, " \n"); +#endif /* USE_JPWL */ + /* <tile_size_on = OPJ_FALSE; + parameters->cp_tdx = 1; + parameters->cp_tdy = 1; + + /*Tile part*/ + parameters->tp_flag = 'C'; + parameters->tp_on = 1; + + /*Tile and Image shall be at (0,0)*/ + parameters->cp_tx0 = 0; + parameters->cp_ty0 = 0; + parameters->image_offset_x0 = 0; + parameters->image_offset_y0 = 0; + + /*Codeblock size= 32*32*/ + parameters->cblockw_init = 32; + parameters->cblockh_init = 32; + parameters->csty |= 0x01; + + /*The progression order shall be CPRL*/ + parameters->prog_order = CPRL; + + /* No ROI */ + parameters->roi_compno = -1; + + parameters->subsampling_dx = 1; + parameters->subsampling_dy = 1; + + /* 9-7 transform */ + parameters->irreversible = 1; + +} + +static void cinema_setup_encoder(opj_cparameters_t *parameters, + opj_image_t *image, img_fol_t *img_fol) +{ + int i; + float temp_rate; + opj_poc_t *POC = NULL; + + switch (parameters->cp_cinema) { + case CINEMA2K_24: + case CINEMA2K_48: + if (parameters->numresolution > 6) { + parameters->numresolution = 6; + } + if (!((image->comps[0].w == 2048) | (image->comps[0].h == 1080))) { + fprintf(stdout, + "Image coordinates %d x %d is not 2K compliant.\nJPEG Digital Cinema Profile-3" + "(2K profile) compliance requires that at least one of coordinates match 2048 x 1080\n", + image->comps[0].w, image->comps[0].h); + parameters->cp_rsiz = STD_RSIZ; + } + break; + + case CINEMA4K_24: + if (parameters->numresolution < 1) { + parameters->numresolution = 1; + } else if (parameters->numresolution > 7) { + parameters->numresolution = 7; + } + if (!((image->comps[0].w == 4096) | (image->comps[0].h == 2160))) { + fprintf(stdout, + "Image coordinates %d x %d is not 4K compliant.\nJPEG Digital Cinema Profile-4" + "(4K profile) compliance requires that at least one of coordinates match 4096 x 2160\n", + image->comps[0].w, image->comps[0].h); + parameters->cp_rsiz = STD_RSIZ; + } + parameters->numpocs = initialise_4K_poc(parameters->POC, + parameters->numresolution); + break; + } + + switch (parameters->cp_cinema) { + case CINEMA2K_24: + case CINEMA4K_24: + for (i = 0 ; i < parameters->tcp_numlayers ; i++) { + temp_rate = 0 ; + if (img_fol->rates[i] == 0) { + parameters->tcp_rates[0] = ((float)(image->numcomps * image->comps[0].w * + image->comps[0].h * image->comps[0].prec)) / + (CINEMA_24_CS * 8 * image->comps[0].dx * image->comps[0].dy); + } else { + temp_rate = ((float)(image->numcomps * image->comps[0].w * image->comps[0].h * + image->comps[0].prec)) / + (img_fol->rates[i] * 8 * image->comps[0].dx * image->comps[0].dy); + if (temp_rate > CINEMA_24_CS) { + parameters->tcp_rates[i] = ((float)(image->numcomps * image->comps[0].w * + image->comps[0].h * image->comps[0].prec)) / + (CINEMA_24_CS * 8 * image->comps[0].dx * image->comps[0].dy); + } else { + parameters->tcp_rates[i] = img_fol->rates[i]; + } + } + } + parameters->max_comp_size = COMP_24_CS; + break; + + case CINEMA2K_48: + for (i = 0 ; i < parameters->tcp_numlayers ; i++) { + temp_rate = 0 ; + if (img_fol->rates[i] == 0) { + parameters->tcp_rates[0] = ((float)(image->numcomps * image->comps[0].w * + image->comps[0].h * image->comps[0].prec)) / + (CINEMA_48_CS * 8 * image->comps[0].dx * image->comps[0].dy); + } else { + temp_rate = ((float)(image->numcomps * image->comps[0].w * image->comps[0].h * + image->comps[0].prec)) / + (img_fol->rates[i] * 8 * image->comps[0].dx * image->comps[0].dy); + if (temp_rate > CINEMA_48_CS) { + parameters->tcp_rates[0] = ((float)(image->numcomps * image->comps[0].w * + image->comps[0].h * image->comps[0].prec)) / + (CINEMA_48_CS * 8 * image->comps[0].dx * image->comps[0].dy); + } else { + parameters->tcp_rates[i] = img_fol->rates[i]; + } + } + } + parameters->max_comp_size = COMP_48_CS; + break; + } + parameters->cp_disto_alloc = 1; +} + + +/* ------------------------------------------------------------------------------------ */ +static int parse_cmdline_encoder(int argc, char **argv, + opj_cparameters_t *parameters, + img_fol_t *img_fol, char *indexfilename) +{ + int i, j, totlen; + opj_option_t long_option[] = { + {"cinema2K", REQ_ARG, NULL, 'w'}, + {"cinema4K", NO_ARG, NULL, 'y'}, + {"ImgDir", REQ_ARG, NULL, 'z'}, + {"TP", REQ_ARG, NULL, 'u'}, + {"SOP", NO_ARG, NULL, 'S'}, + {"EPH", NO_ARG, NULL, 'E'}, + {"OutFor", REQ_ARG, NULL, 'O'}, + {"POC", REQ_ARG, NULL, 'P'}, + {"ROI", REQ_ARG, NULL, 'R'}, + {"jpip", NO_ARG, NULL, 'J'} + }; + + /* parse the command line */ + /* UniPG>> */ + const char optlist[] = "i:o:hr:q:n:b:c:t:p:s:SEM:x:R:d:T:If:P:C:F:u:" +#ifdef USE_JPWL + "W:" +#endif /* USE_JPWL */ + ; + + /*printf("C: parse_cmdline_encoder:"); + for (i=0; iset_out_format = 0; + reset_options_reading(); + + while (1) { + int c = opj_getopt_long(argc, argv, optlist, long_option, totlen); + if (c == -1) { + break; + } + switch (c) { + + /* ----------------------------------------------------- */ + + case 'o': { /* output file */ + char *outfile = opj_optarg; + parameters->cod_format = get_file_format(outfile); + switch (parameters->cod_format) { + case J2K_CFMT: + case JP2_CFMT: + break; + default: + fprintf(stderr, + "Unknown output format image %s [only *.j2k, *.j2c or *.jp2]!! \n", outfile); + return 1; + } + strncpy(parameters->outfile, outfile, sizeof(parameters->outfile) - 1); + } + break; + + /* ----------------------------------------------------- */ + case 'O': { /* output format */ + char outformat[50]; + char *of = opj_optarg; + sprintf(outformat, ".%s", of); + img_fol->set_out_format = 1; + parameters->cod_format = get_file_format(outformat); + switch (parameters->cod_format) { + case J2K_CFMT: + case JP2_CFMT: + img_fol->out_format = opj_optarg; + break; + default: + fprintf(stderr, "Unknown output format image [only j2k, j2c, jp2]!! \n"); + return 1; + } + } + break; + + + /* ----------------------------------------------------- */ + + + case 'r': { /* rates rates/distorsion */ + char *s = opj_optarg; + while (sscanf(s, "%f", ¶meters->tcp_rates[parameters->tcp_numlayers]) == + 1) { + parameters->tcp_numlayers++; + while (*s && *s != ',') { + s++; + } + if (!*s) { + break; + } + s++; + } + parameters->cp_disto_alloc = 1; + } + break; + + /* ----------------------------------------------------- */ + + case 'q': { /* add fixed_quality */ + char *s = opj_optarg; + while (sscanf(s, "%f", ¶meters->tcp_distoratio[parameters->tcp_numlayers]) + == 1) { + parameters->tcp_numlayers++; + while (*s && *s != ',') { + s++; + } + if (!*s) { + break; + } + s++; + } + parameters->cp_fixed_quality = 1; + } + break; + + /* dda */ + /* ----------------------------------------------------- */ + + case 'f': { /* mod fixed_quality (before : -q) */ + int *row = NULL, *col = NULL; + int numlayers = 0, numresolution = 0, matrix_width = 0; + + char *s = opj_optarg; + sscanf(s, "%d", &numlayers); + s++; + if (numlayers > 9) { + s++; + } + + parameters->tcp_numlayers = numlayers; + numresolution = parameters->numresolution; + matrix_width = numresolution * 3; + parameters->cp_matrice = (int *) opj_malloc(numlayers * matrix_width * sizeof( + int)); + s = s + 2; + + for (i = 0; i < numlayers; i++) { + row = ¶meters->cp_matrice[i * matrix_width]; + col = row; + parameters->tcp_rates[i] = 1; + sscanf(s, "%d,", &col[0]); + s += 2; + if (col[0] > 9) { + s++; + } + col[1] = 0; + col[2] = 0; + for (j = 1; j < numresolution; j++) { + col += 3; + sscanf(s, "%d,%d,%d", &col[0], &col[1], &col[2]); + s += 6; + if (col[0] > 9) { + s++; + } + if (col[1] > 9) { + s++; + } + if (col[2] > 9) { + s++; + } + } + if (i < numlayers - 1) { + s++; + } + } + parameters->cp_fixed_alloc = 1; + } + break; + + /* ----------------------------------------------------- */ + + case 't': { /* tiles */ + sscanf(opj_optarg, "%d,%d", ¶meters->cp_tdx, ¶meters->cp_tdy); + parameters->tile_size_on = OPJ_TRUE; + } + break; + + /* ----------------------------------------------------- */ + + case 'n': { /* resolution */ + sscanf(opj_optarg, "%d", ¶meters->numresolution); + } + break; + + /* ----------------------------------------------------- */ + case 'c': { /* precinct dimension */ + char sep; + int res_spec = 0; + + char *s = opj_optarg; + do { + sep = 0; + sscanf(s, "[%d,%d]%c", ¶meters->prcw_init[res_spec], + ¶meters->prch_init[res_spec], &sep); + parameters->csty |= 0x01; + res_spec++; + s = strpbrk(s, "]") + 2; + } while (sep == ','); + parameters->res_spec = res_spec; + } + break; + + /* ----------------------------------------------------- */ + + case 'b': { /* code-block dimension */ + int cblockw_init = 0, cblockh_init = 0; + sscanf(opj_optarg, "%d,%d", &cblockw_init, &cblockh_init); + if (cblockw_init * cblockh_init > 4096 || cblockw_init > 1024 + || cblockw_init < 4 || cblockh_init > 1024 || cblockh_init < 4) { + fprintf(stderr, + "!! Size of code_block error (option -b) !!\n\nRestriction :\n" + " * width*height<=4096\n * 4<=width,height<= 1024\n\n"); + return 1; + } + parameters->cblockw_init = cblockw_init; + parameters->cblockh_init = cblockh_init; + } + break; + + /* ----------------------------------------------------- */ + + case 'x': { /* creation of index file */ + char *index = opj_optarg; + strncpy(indexfilename, index, OPJ_PATH_LEN); + } + break; + + /* ----------------------------------------------------- */ + + case 'p': { /* progression order */ + char progression[4]; + + strncpy(progression, opj_optarg, 4); + parameters->prog_order = give_progression(progression); + if (parameters->prog_order == -1) { + fprintf(stderr, "Unrecognized progression order " + "[LRCP, RLCP, RPCL, PCRL, CPRL] !!\n"); + return 1; + } + } + break; + + /* ----------------------------------------------------- */ + + case 's': { /* subsampling factor */ + if (sscanf(opj_optarg, "%d,%d", ¶meters->subsampling_dx, + ¶meters->subsampling_dy) != 2) { + fprintf(stderr, "'-s' sub-sampling argument error ! [-s dx,dy]\n"); + return 1; + } + } + break; + + /* ----------------------------------------------------- */ + + case 'd': { /* coordonnate of the reference grid */ + if (sscanf(opj_optarg, "%d,%d", ¶meters->image_offset_x0, + ¶meters->image_offset_y0) != 2) { + fprintf(stderr, "-d 'coordonnate of the reference grid' argument " + "error !! [-d x0,y0]\n"); + return 1; + } + } + break; + + /* ----------------------------------------------------- */ + + case 'h': /* display an help description */ + encode_help_display(); + return 1; + + /* ----------------------------------------------------- */ + + case 'P': { /* POC */ + int numpocs = 0; /* number of progression order change (POC) default 0 */ + opj_poc_t *POC = NULL; /* POC : used in case of Progression order change */ + + char *s = opj_optarg; + POC = parameters->POC; + + while (sscanf(s, "T%d=%d,%d,%d,%d,%d,%4s", &POC[numpocs].tile, + &POC[numpocs].resno0, &POC[numpocs].compno0, + &POC[numpocs].layno1, &POC[numpocs].resno1, + &POC[numpocs].compno1, &POC[numpocs].progorder) == 7) { + POC[numpocs].prg1 = give_progression(POC[numpocs].progorder); + numpocs++; + while (*s && *s != '/') { + s++; + } + if (!*s) { + break; + } + s++; + } + parameters->numpocs = numpocs; + } + break; + + /* ------------------------------------------------------ */ + + case 'S': { /* SOP marker */ + parameters->csty |= 0x02; + } + break; + + /* ------------------------------------------------------ */ + + case 'E': { /* EPH marker */ + parameters->csty |= 0x04; + } + break; + + /* ------------------------------------------------------ */ + + case 'M': { /* Mode switch pas tous au point !! */ + int value = 0; + if (sscanf(opj_optarg, "%d", &value) == 1) { + for (i = 0; i <= 5; i++) { + int cache = value & (1 << i); + if (cache) { + parameters->mode |= (1 << i); + } + } + } + } + break; + + /* ------------------------------------------------------ */ + + case 'R': { /* ROI */ + if (sscanf(opj_optarg, "c=%d,U=%d", ¶meters->roi_compno, + ¶meters->roi_shift) != 2) { + fprintf(stderr, "ROI error !! [-ROI c='compno',U='shift']\n"); + return 1; + } + } + break; + + /* ------------------------------------------------------ */ + + case 'T': { /* Tile offset */ + if (sscanf(opj_optarg, "%d,%d", ¶meters->cp_tx0, + ¶meters->cp_ty0) != 2) { + fprintf(stderr, "-T 'tile offset' argument error !! [-T X0,Y0]"); + return 1; + } + } + break; + + /* ------------------------------------------------------ */ + + case 'C': { /* add a comment */ + parameters->cp_comment = (char*)opj_malloc(strlen(opj_optarg) + 1); + if (parameters->cp_comment) { + strcpy(parameters->cp_comment, opj_optarg); + } + } + break; + + + /* ------------------------------------------------------ */ + + case 'I': { /* reversible or not */ + parameters->irreversible = 1; + } + break; + + /* ------------------------------------------------------ */ + + case 'u': { /* Tile part generation*/ + parameters->tp_flag = opj_optarg[0]; + parameters->tp_on = 1; + } + break; + + /* ------------------------------------------------------ */ + + case 'z': { /* Image Directory path */ + img_fol->imgdirpath = (char*)opj_malloc(strlen(opj_optarg) + 1); + strcpy(img_fol->imgdirpath, opj_optarg); + img_fol->set_imgdir = 1; + } + break; + + /* ------------------------------------------------------ */ + + case 'w': { /* Digital Cinema 2K profile compliance*/ + int fps = 0; + sscanf(opj_optarg, "%d", &fps); + if (fps == 24) { + parameters->cp_cinema = CINEMA2K_24; + } else if (fps == 48) { + parameters->cp_cinema = CINEMA2K_48; + } else { + fprintf(stderr, "Incorrect value!! must be 24 or 48\n"); + return 1; + } + fprintf(stdout, "CINEMA 2K compliant codestream\n"); + parameters->cp_rsiz = CINEMA2K; + + } + break; + + /* ------------------------------------------------------ */ + + case 'y': { /* Digital Cinema 4K profile compliance*/ + parameters->cp_cinema = CINEMA4K_24; + fprintf(stdout, "CINEMA 4K compliant codestream\n"); + parameters->cp_rsiz = CINEMA4K; + } + break; + + /* ------------------------------------------------------ */ + + /* UniPG>> */ +#ifdef USE_JPWL + /* ------------------------------------------------------ */ + + case 'W': { /* JPWL capabilities switched on */ + char *token = NULL; + int hprot, pprot, sens, addr, size, range; + + /* we need to enable indexing */ + if (!indexfilename) { + strncpy(indexfilename, JPWL_PRIVATEINDEX_NAME, OPJ_PATH_LEN); + } + + /* search for different protection methods */ + + /* break the option in comma points and parse the result */ + token = strtok(opj_optarg, ","); + while (token != NULL) { + + /* search header error protection method */ + if (*token == 'h') { + + static int tile = 0, tilespec = 0, lasttileno = 0; + + hprot = 1; /* predefined method */ + + if (sscanf(token, "h=%d", &hprot) == 1) { + /* Main header, specified */ + if (!((hprot == 0) || (hprot == 1) || (hprot == 16) || (hprot == 32) || + ((hprot >= 37) && (hprot <= 128)))) { + fprintf(stderr, "ERROR -> invalid main header protection method h = %d\n", + hprot); + return 1; + } + parameters->jpwl_hprot_MH = hprot; + + } else if (sscanf(token, "h%d=%d", &tile, &hprot) == 2) { + /* Tile part header, specified */ + if (!((hprot == 0) || (hprot == 1) || (hprot == 16) || (hprot == 32) || + ((hprot >= 37) && (hprot <= 128)))) { + fprintf(stderr, "ERROR -> invalid tile part header protection method h = %d\n", + hprot); + return 1; + } + if (tile < 0) { + fprintf(stderr, + "ERROR -> invalid tile part number on protection method t = %d\n", tile); + return 1; + } + if (tilespec < JPWL_MAX_NO_TILESPECS) { + parameters->jpwl_hprot_TPH_tileno[tilespec] = lasttileno = tile; + parameters->jpwl_hprot_TPH[tilespec++] = hprot; + } + + } else if (sscanf(token, "h%d", &tile) == 1) { + /* Tile part header, unspecified */ + if (tile < 0) { + fprintf(stderr, + "ERROR -> invalid tile part number on protection method t = %d\n", tile); + return 1; + } + if (tilespec < JPWL_MAX_NO_TILESPECS) { + parameters->jpwl_hprot_TPH_tileno[tilespec] = lasttileno = tile; + parameters->jpwl_hprot_TPH[tilespec++] = hprot; + } + + + } else if (!strcmp(token, "h")) { + /* Main header, unspecified */ + parameters->jpwl_hprot_MH = hprot; + + } else { + fprintf(stderr, "ERROR -> invalid protection method selection = %s\n", token); + return 1; + }; + + } + + /* search packet error protection method */ + if (*token == 'p') { + + static int pack = 0, tile = 0, packspec = 0, lastpackno = 0; + + pprot = 1; /* predefined method */ + + if (sscanf(token, "p=%d", &pprot) == 1) { + /* Method for all tiles and all packets */ + if (!((pprot == 0) || (pprot == 1) || (pprot == 16) || (pprot == 32) || + ((pprot >= 37) && (pprot <= 128)))) { + fprintf(stderr, "ERROR -> invalid default packet protection method p = %d\n", + pprot); + return 1; + } + parameters->jpwl_pprot_tileno[0] = 0; + parameters->jpwl_pprot_packno[0] = 0; + parameters->jpwl_pprot[0] = pprot; + + } else if (sscanf(token, "p%d=%d", &tile, &pprot) == 2) { + /* method specified from that tile on */ + if (!((pprot == 0) || (pprot == 1) || (pprot == 16) || (pprot == 32) || + ((pprot >= 37) && (pprot <= 128)))) { + fprintf(stderr, "ERROR -> invalid packet protection method p = %d\n", pprot); + return 1; + } + if (tile < 0) { + fprintf(stderr, + "ERROR -> invalid tile part number on protection method p = %d\n", tile); + return 1; + } + if (packspec < JPWL_MAX_NO_PACKSPECS) { + parameters->jpwl_pprot_tileno[packspec] = tile; + parameters->jpwl_pprot_packno[packspec] = 0; + parameters->jpwl_pprot[packspec++] = pprot; + } + + } else if (sscanf(token, "p%d:%d=%d", &tile, &pack, &pprot) == 3) { + /* method fully specified from that tile and that packet on */ + if (!((pprot == 0) || (pprot == 1) || (pprot == 16) || (pprot == 32) || + ((pprot >= 37) && (pprot <= 128)))) { + fprintf(stderr, "ERROR -> invalid packet protection method p = %d\n", pprot); + return 1; + } + if (tile < 0) { + fprintf(stderr, + "ERROR -> invalid tile part number on protection method p = %d\n", tile); + return 1; + } + if (pack < 0) { + fprintf(stderr, "ERROR -> invalid packet number on protection method p = %d\n", + pack); + return 1; + } + if (packspec < JPWL_MAX_NO_PACKSPECS) { + parameters->jpwl_pprot_tileno[packspec] = tile; + parameters->jpwl_pprot_packno[packspec] = pack; + parameters->jpwl_pprot[packspec++] = pprot; + } + + } else if (sscanf(token, "p%d:%d", &tile, &pack) == 2) { + /* default method from that tile and that packet on */ + if (!((pprot == 0) || (pprot == 1) || (pprot == 16) || (pprot == 32) || + ((pprot >= 37) && (pprot <= 128)))) { + fprintf(stderr, "ERROR -> invalid packet protection method p = %d\n", pprot); + return 1; + } + if (tile < 0) { + fprintf(stderr, + "ERROR -> invalid tile part number on protection method p = %d\n", tile); + return 1; + } + if (pack < 0) { + fprintf(stderr, "ERROR -> invalid packet number on protection method p = %d\n", + pack); + return 1; + } + if (packspec < JPWL_MAX_NO_PACKSPECS) { + parameters->jpwl_pprot_tileno[packspec] = tile; + parameters->jpwl_pprot_packno[packspec] = pack; + parameters->jpwl_pprot[packspec++] = pprot; + } + + } else if (sscanf(token, "p%d", &tile) == 1) { + /* default from a tile on */ + if (tile < 0) { + fprintf(stderr, + "ERROR -> invalid tile part number on protection method p = %d\n", tile); + return 1; + } + if (packspec < JPWL_MAX_NO_PACKSPECS) { + parameters->jpwl_pprot_tileno[packspec] = tile; + parameters->jpwl_pprot_packno[packspec] = 0; + parameters->jpwl_pprot[packspec++] = pprot; + } + + + } else if (!strcmp(token, "p")) { + /* all default */ + parameters->jpwl_pprot_tileno[0] = 0; + parameters->jpwl_pprot_packno[0] = 0; + parameters->jpwl_pprot[0] = pprot; + + } else { + fprintf(stderr, "ERROR -> invalid protection method selection = %s\n", token); + return 1; + }; + + } + + /* search sensitivity method */ + if (*token == 's') { + + static int tile = 0, tilespec = 0, lasttileno = 0; + + sens = 0; /* predefined: relative error */ + + if (sscanf(token, "s=%d", &sens) == 1) { + /* Main header, specified */ + if ((sens < -1) || (sens > 7)) { + fprintf(stderr, "ERROR -> invalid main header sensitivity method s = %d\n", + sens); + return 1; + } + parameters->jpwl_sens_MH = sens; + + } else if (sscanf(token, "s%d=%d", &tile, &sens) == 2) { + /* Tile part header, specified */ + if ((sens < -1) || (sens > 7)) { + fprintf(stderr, "ERROR -> invalid tile part header sensitivity method s = %d\n", + sens); + return 1; + } + if (tile < 0) { + fprintf(stderr, + "ERROR -> invalid tile part number on sensitivity method t = %d\n", tile); + return 1; + } + if (tilespec < JPWL_MAX_NO_TILESPECS) { + parameters->jpwl_sens_TPH_tileno[tilespec] = lasttileno = tile; + parameters->jpwl_sens_TPH[tilespec++] = sens; + } + + } else if (sscanf(token, "s%d", &tile) == 1) { + /* Tile part header, unspecified */ + if (tile < 0) { + fprintf(stderr, + "ERROR -> invalid tile part number on sensitivity method t = %d\n", tile); + return 1; + } + if (tilespec < JPWL_MAX_NO_TILESPECS) { + parameters->jpwl_sens_TPH_tileno[tilespec] = lasttileno = tile; + parameters->jpwl_sens_TPH[tilespec++] = hprot; + } + + } else if (!strcmp(token, "s")) { + /* Main header, unspecified */ + parameters->jpwl_sens_MH = sens; + + } else { + fprintf(stderr, "ERROR -> invalid sensitivity method selection = %s\n", token); + return 1; + }; + + parameters->jpwl_sens_size = 2; /* 2 bytes for default size */ + } + + /* search addressing size */ + if (*token == 'a') { + + static int tile = 0, tilespec = 0, lasttileno = 0; + + addr = 0; /* predefined: auto */ + + if (sscanf(token, "a=%d", &addr) == 1) { + /* Specified */ + if ((addr != 0) && (addr != 2) && (addr != 4)) { + fprintf(stderr, "ERROR -> invalid addressing size a = %d\n", addr); + return 1; + } + parameters->jpwl_sens_addr = addr; + + } else if (!strcmp(token, "a")) { + /* default */ + parameters->jpwl_sens_addr = addr; /* auto for default size */ + + } else { + fprintf(stderr, "ERROR -> invalid addressing selection = %s\n", token); + return 1; + }; + + } + + /* search sensitivity size */ + if (*token == 'z') { + + static int tile = 0, tilespec = 0, lasttileno = 0; + + size = 1; /* predefined: 1 byte */ + + if (sscanf(token, "z=%d", &size) == 1) { + /* Specified */ + if ((size != 0) && (size != 1) && (size != 2)) { + fprintf(stderr, "ERROR -> invalid sensitivity size z = %d\n", size); + return 1; + } + parameters->jpwl_sens_size = size; + + } else if (!strcmp(token, "a")) { + /* default */ + parameters->jpwl_sens_size = size; /* 1 for default size */ + + } else { + fprintf(stderr, "ERROR -> invalid size selection = %s\n", token); + return 1; + }; + + } + + /* search range method */ + if (*token == 'g') { + + static int tile = 0, tilespec = 0, lasttileno = 0; + + range = 0; /* predefined: 0 (packet) */ + + if (sscanf(token, "g=%d", &range) == 1) { + /* Specified */ + if ((range < 0) || (range > 3)) { + fprintf(stderr, "ERROR -> invalid sensitivity range method g = %d\n", range); + return 1; + } + parameters->jpwl_sens_range = range; + + } else if (!strcmp(token, "g")) { + /* default */ + parameters->jpwl_sens_range = range; + + } else { + fprintf(stderr, "ERROR -> invalid range selection = %s\n", token); + return 1; + }; + + } + + /* next token or bust */ + token = strtok(NULL, ","); + }; + + + /* some info */ + fprintf(stdout, "Info: JPWL capabilities enabled\n"); + parameters->jpwl_epc_on = true; + + } + break; +#endif /* USE_JPWL */ + /* < Command line not valid\n"); + return 1; + } + } + + /* check for possible errors */ + if (parameters->cp_cinema) { + if (parameters->tcp_numlayers > 1) { + parameters->cp_rsiz = STD_RSIZ; + fprintf(stdout, + "Warning: DC profiles do not allow more than one quality layer. The codestream created will not be compliant with the DC profile\n"); + } + } + + if ((parameters->cp_disto_alloc || parameters->cp_fixed_alloc || + parameters->cp_fixed_quality) + && (!(parameters->cp_disto_alloc ^ parameters->cp_fixed_alloc ^ + parameters->cp_fixed_quality))) { + fprintf(stderr, "Error: options -r -q and -f cannot be used together !!\n"); + return 1; + } /* mod fixed_quality */ + + /* if no rate entered, lossless by default */ + if (parameters->tcp_numlayers == 0) { + parameters->tcp_rates[0] = 0; /* MOD antonin : losslessbug */ + parameters->tcp_numlayers++; + parameters->cp_disto_alloc = 1; + } + + if ((parameters->cp_tx0 > parameters->image_offset_x0) || + (parameters->cp_ty0 > parameters->image_offset_y0)) { + fprintf(stderr, + "Error: Tile offset dimension is unnappropriate --> TX0(%d)<=IMG_X0(%d) TYO(%d)<=IMG_Y0(%d) \n", + parameters->cp_tx0, parameters->image_offset_x0, parameters->cp_ty0, + parameters->image_offset_y0); + return 1; + } + + for (i = 0; i < parameters->numpocs; i++) { + if (parameters->POC[i].prg == -1) { + fprintf(stderr, + "Unrecognized progression order in option -P (POC n %d) [LRCP, RLCP, RPCL, PCRL, CPRL] !!\n", + i + 1); + } + } + + return 0; +} + + +/** Create the same index as j2k_create_index does, but in an int[] instead of in a file ==> easy to pass it back to Java, to transfer it etc. + @param buffer_size, increased by the length of the compressed index, in number of bytes + @return a pointer to a char[] + Syntax of the index: + one char for the version number (1): one byte because no problem with little endian, big endian etc. + one int for each of the following informations: + Image Width + Image Height + progression order + Tile width + Tile height + Nb tiles in X + Nb tiles in Y + Nb of components + Nb of layers + Nb of resolutions + + for each resolution: + Precinct width + Precinct height + + End main header position + codestream size + + For each tile: + tile number + tile start pos in codestream + tile header end position + tile end position in codestream + + For each LRCP, RLCP etc.: + packet number + tile number + layer number + resolution number + component number + precinct number + start position in the codestream + end position of this packet + */ +static char* create_index_into_byte_array(opj_codestream_info_t *cstr_info, + int* buffer_size) +{ + int tileno, compno, layno, resno, precno, pack_nb, x, y; + char* buffer = NULL; + int buffer_pos = 0; + int prec_max = 0; + + prec_max = 0; + for (tileno = 0; tileno < cstr_info->tw * cstr_info->th; tileno++) { + for (resno = 0; resno < cstr_info->numdecompos[0] + 1; resno++) { + prec_max = int_max(prec_max, + cstr_info->tile[tileno].pw[resno] * cstr_info->tile[tileno].ph[resno]); + } + } + + /* Compute the size of the index buffer, in number of bytes*/ + *buffer_size = + 1 /* version */ + + (10 /* image_w until decomposition */ + + (cstr_info->numdecompos[0] + 1) * 2 /* pdx size for each tile */ + + 2 /* main_head_end + codestream_size */ + + cstr_info->tw * cstr_info->th * 4 /* tile info, without distorsion info */ + + cstr_info->tw * cstr_info->th * cstr_info->numlayers * + (cstr_info->numdecompos[0] + 1) * cstr_info->numcomps * prec_max * 8 + ) * sizeof(int); + /*printf("C: index buffer size = %d bytes\n", *buffer_size);*/ + buffer = (char*) opj_malloc(*buffer_size); + + if (!buffer) { + /*opj_event_msg(j2k->cinfo, EVT_ERROR, "failed to allocate index buffer for writing %d int\n", *buffer_size);*/ + fprintf(stderr, "failed to allocate index buffer for writing %d int\n", + *buffer_size); + return 0; + } + + buffer[0] = 1; /* Version stored on a byte*/ + buffer++; + /* Remaining informations are stored on a int.*/ + ((int*)buffer)[buffer_pos++] = cstr_info->image_w; + ((int*)buffer)[buffer_pos++] = cstr_info->image_h; + ((int*)buffer)[buffer_pos++] = cstr_info->prog; + ((int*)buffer)[buffer_pos++] = cstr_info->tile_x; + ((int*)buffer)[buffer_pos++] = cstr_info->tile_y; + ((int*)buffer)[buffer_pos++] = cstr_info->tw; + ((int*)buffer)[buffer_pos++] = cstr_info->th; + ((int*)buffer)[buffer_pos++] = cstr_info->numcomps; + ((int*)buffer)[buffer_pos++] = cstr_info->numlayers; + ((int*)buffer)[buffer_pos++] = cstr_info->numdecompos[0]; + + for (resno = cstr_info->numdecompos[0]; resno >= 0; resno--) { + /* based on tile 0 */ + ((int*)buffer)[buffer_pos++] = (1 << cstr_info->tile[0].pdx[resno]); + ((int*)buffer)[buffer_pos++] = (1 << cstr_info->tile[0].pdx[resno]); + } + ((int*)buffer)[buffer_pos++] = cstr_info->main_head_end; + ((int*)buffer)[buffer_pos++] = cstr_info->codestream_size; + + for (tileno = 0; tileno < cstr_info->tw * cstr_info->th; tileno++) { + ((int*)buffer)[buffer_pos++] = cstr_info->tile[tileno].tileno; + ((int*)buffer)[buffer_pos++] = cstr_info->tile[tileno].start_pos; + ((int*)buffer)[buffer_pos++] = cstr_info->tile[tileno].end_header; + ((int*)buffer)[buffer_pos++] = cstr_info->tile[tileno].end_pos; + } + + for (tileno = 0; tileno < cstr_info->tw * cstr_info->th; tileno++) { + int start_pos, end_pos; + int max_numdecompos = 0; + pack_nb = 0; + + for (compno = 0; compno < cstr_info->numcomps; compno++) { + if (max_numdecompos < cstr_info->numdecompos[compno]) { + max_numdecompos = cstr_info->numdecompos[compno]; + } + } + + if (cstr_info->prog == LRCP) { /* LRCP */ + + for (layno = 0; layno < cstr_info->numlayers; layno++) { + for (resno = 0; resno < max_numdecompos + 1; resno++) { + for (compno = 0; compno < cstr_info->numcomps; compno++) { + int prec_max; + if (resno > cstr_info->numdecompos[compno]) { + break; + } + prec_max = cstr_info->tile[tileno].pw[resno] * + cstr_info->tile[tileno].ph[resno]; + for (precno = 0; precno < prec_max; precno++) { + start_pos = cstr_info->tile[tileno].packet[pack_nb].start_pos; + end_pos = cstr_info->tile[tileno].packet[pack_nb].end_pos; + ((int*)buffer)[buffer_pos++] = pack_nb; + ((int*)buffer)[buffer_pos++] = tileno; + ((int*)buffer)[buffer_pos++] = layno; + ((int*)buffer)[buffer_pos++] = resno; + ((int*)buffer)[buffer_pos++] = compno; + ((int*)buffer)[buffer_pos++] = precno; + ((int*)buffer)[buffer_pos++] = start_pos; + ((int*)buffer)[buffer_pos++] = end_pos; + pack_nb++; + } + } + } + } + } /* LRCP */ + else if (cstr_info->prog == RLCP) { /* RLCP */ + + for (resno = 0; resno < max_numdecompos + 1; resno++) { + for (layno = 0; layno < cstr_info->numlayers; layno++) { + for (compno = 0; compno < cstr_info->numcomps; compno++) { + int prec_max; + if (resno > cstr_info->numdecompos[compno]) { + break; + } + prec_max = cstr_info->tile[tileno].pw[resno] * + cstr_info->tile[tileno].ph[resno]; + for (precno = 0; precno < prec_max; precno++) { + start_pos = cstr_info->tile[tileno].packet[pack_nb].start_pos; + end_pos = cstr_info->tile[tileno].packet[pack_nb].end_pos; + ((int*)buffer)[buffer_pos++] = pack_nb; + ((int*)buffer)[buffer_pos++] = tileno; + ((int*)buffer)[buffer_pos++] = resno; + ((int*)buffer)[buffer_pos++] = layno; + ((int*)buffer)[buffer_pos++] = compno; + ((int*)buffer)[buffer_pos++] = precno; + ((int*)buffer)[buffer_pos++] = start_pos; + ((int*)buffer)[buffer_pos++] = end_pos; + pack_nb++; + } + } + } + } + } /* RLCP */ + else if (cstr_info->prog == RPCL) { /* RPCL */ + + for (resno = 0; resno < max_numdecompos + 1; resno++) { + /* I suppose components have same XRsiz, YRsiz */ + int x0 = cstr_info->tile_Ox + tileno - (int)floor((float)tileno / + (float)cstr_info->tw) * cstr_info->tw * cstr_info->tile_x; + int y0 = cstr_info->tile_Ox + (int)floor((float)tileno / + (float)cstr_info->tw) * cstr_info->tile_y; + int x1 = x0 + cstr_info->tile_x; + int y1 = y0 + cstr_info->tile_y; + for (compno = 0; compno < cstr_info->numcomps; compno++) { + int prec_max; + if (resno > cstr_info->numdecompos[compno]) { + break; + } + prec_max = cstr_info->tile[tileno].pw[resno] * + cstr_info->tile[tileno].ph[resno]; + for (precno = 0; precno < prec_max; precno++) { + int pcnx = cstr_info->tile[tileno].pw[resno]; + int pcx = (int) pow(2, cstr_info->tile[tileno].pdx[resno] + + cstr_info->numdecompos[compno] - resno); + int pcy = (int) pow(2, cstr_info->tile[tileno].pdy[resno] + + cstr_info->numdecompos[compno] - resno); + int precno_x = precno - (int) floor((float)precno / (float)pcnx) * pcnx; + int precno_y = (int) floor((float)precno / (float)pcnx); + for (y = y0; y < y1; y++) { + if (precno_y * pcy == y) { + for (x = x0; x < x1; x++) { + if (precno_x * pcx == x) { + for (layno = 0; layno < cstr_info->numlayers; layno++) { + start_pos = cstr_info->tile[tileno].packet[pack_nb].start_pos; + end_pos = cstr_info->tile[tileno].packet[pack_nb].end_pos; + ((int*)buffer)[buffer_pos++] = pack_nb; + ((int*)buffer)[buffer_pos++] = tileno; + ((int*)buffer)[buffer_pos++] = resno; + ((int*)buffer)[buffer_pos++] = precno; + ((int*)buffer)[buffer_pos++] = compno; + ((int*)buffer)[buffer_pos++] = layno; + ((int*)buffer)[buffer_pos++] = start_pos; + ((int*)buffer)[buffer_pos++] = end_pos; + pack_nb++; + } + } + }/* x = x0..x1 */ + } + } /* y = y0..y1 */ + } /* precno */ + } /* compno */ + } /* resno */ + } /* RPCL */ + else if (cstr_info->prog == PCRL) { /* PCRL */ + /* I suppose components have same XRsiz, YRsiz */ + int x0 = cstr_info->tile_Ox + tileno - (int)floor((float)tileno / + (float)cstr_info->tw) * cstr_info->tw * cstr_info->tile_x; + int y0 = cstr_info->tile_Ox + (int)floor((float)tileno / + (float)cstr_info->tw) * cstr_info->tile_y; + int x1 = x0 + cstr_info->tile_x; + int y1 = y0 + cstr_info->tile_y; + + for (compno = 0; compno < cstr_info->numcomps; compno++) { + for (resno = 0; resno < max_numdecompos + 1; resno++) { + int prec_max; + if (resno > cstr_info->numdecompos[compno]) { + break; + } + prec_max = cstr_info->tile[tileno].pw[resno] * + cstr_info->tile[tileno].ph[resno]; + for (precno = 0; precno < prec_max; precno++) { + int pcnx = cstr_info->tile[tileno].pw[resno]; + int pcx = (int) pow(2, cstr_info->tile[tileno].pdx[resno] + + cstr_info->numdecompos[compno] - resno); + int pcy = (int) pow(2, cstr_info->tile[tileno].pdy[resno] + + cstr_info->numdecompos[compno] - resno); + int precno_x = precno - (int) floor((float)precno / (float)pcnx) * pcnx; + int precno_y = (int) floor((float)precno / (float)pcnx); + for (y = y0; y < y1; y++) { + if (precno_y * pcy == y) { + for (x = x0; x < x1; x++) { + if (precno_x * pcx == x) { + for (layno = 0; layno < cstr_info->numlayers; layno++) { + start_pos = cstr_info->tile[tileno].packet[pack_nb].start_pos; + end_pos = cstr_info->tile[tileno].packet[pack_nb].end_pos; + ((int*)buffer)[buffer_pos++] = pack_nb; + ((int*)buffer)[buffer_pos++] = tileno; + ((int*)buffer)[buffer_pos++] = precno; + ((int*)buffer)[buffer_pos++] = compno; + ((int*)buffer)[buffer_pos++] = resno; + ((int*)buffer)[buffer_pos++] = layno; + ((int*)buffer)[buffer_pos++] = start_pos; + ((int*)buffer)[buffer_pos++] = end_pos; + pack_nb++; + } + } + }/* x = x0..x1 */ + } + } /* y = y0..y1 */ + } /* precno */ + } /* resno */ + } /* compno */ + } /* PCRL */ + else { /* CPRL */ + + for (compno = 0; compno < cstr_info->numcomps; compno++) { + /* I suppose components have same XRsiz, YRsiz */ + int x0 = cstr_info->tile_Ox + tileno - (int)floor((float)tileno / + (float)cstr_info->tw) * cstr_info->tw * cstr_info->tile_x; + int y0 = cstr_info->tile_Ox + (int)floor((float)tileno / + (float)cstr_info->tw) * cstr_info->tile_y; + int x1 = x0 + cstr_info->tile_x; + int y1 = y0 + cstr_info->tile_y; + + for (resno = 0; resno < max_numdecompos + 1; resno++) { + int prec_max; + if (resno > cstr_info->numdecompos[compno]) { + break; + } + prec_max = cstr_info->tile[tileno].pw[resno] * + cstr_info->tile[tileno].ph[resno]; + for (precno = 0; precno < prec_max; precno++) { + int pcnx = cstr_info->tile[tileno].pw[resno]; + int pcx = (int) pow(2, cstr_info->tile[tileno].pdx[resno] + + cstr_info->numdecompos[compno] - resno); + int pcy = (int) pow(2, cstr_info->tile[tileno].pdy[resno] + + cstr_info->numdecompos[compno] - resno); + int precno_x = precno - (int) floor((float)precno / (float)pcnx) * pcnx; + int precno_y = (int) floor((float)precno / (float)pcnx); + for (y = y0; y < y1; y++) { + if (precno_y * pcy == y) { + for (x = x0; x < x1; x++) { + if (precno_x * pcx == x) { + for (layno = 0; layno < cstr_info->numlayers; layno++) { + start_pos = cstr_info->tile[tileno].packet[pack_nb].start_pos; + end_pos = cstr_info->tile[tileno].packet[pack_nb].end_pos; + ((int*)buffer)[buffer_pos++] = pack_nb; + ((int*)buffer)[buffer_pos++] = tileno; + ((int*)buffer)[buffer_pos++] = compno; + ((int*)buffer)[buffer_pos++] = precno; + ((int*)buffer)[buffer_pos++] = resno; + ((int*)buffer)[buffer_pos++] = layno; + ((int*)buffer)[buffer_pos++] = start_pos; + ((int*)buffer)[buffer_pos++] = end_pos; + pack_nb++; + } + } + }/* x = x0..x1 */ + } + } /* y = y0..y1 */ + } /* precno */ + } /* resno */ + } /* compno */ + } /* CPRL */ + } /* tileno */ + + if (buffer_pos > *buffer_size) { + /*opj_event_msg(j2k->cinfo, EVT_ERROR, "index creation: buffer_pos (%d) > buffer_size (%d)!\n", buffer_pos, *buffer_size);*/ + fprintf(stderr, "index creation: buffer_pos (%d) > buffer_size (%d)!\n", + buffer_pos, *buffer_size); + return 0; + } + + return --buffer; +} + + + + +/* -------------------------------------------------------------------------- + ------------ Get the image byte[] from the Java object -------------------*/ + +static opj_image_t* loadImage(opj_cparameters_t *parameters, JNIEnv *env, + jobject obj, jclass cls) +{ + int i, max, shift, w, h, depth; + opj_image_t * img = NULL; + int compno, numcomps; + opj_image_t * image = NULL; + opj_image_comp_t *comp; + opj_image_cmptparm_t cmptparm[3]; /* maximum of 3 components */ + OPJ_COLOR_SPACE color_space; + jfieldID fid; + jint ji; + jbyteArray jba; + jshortArray jsa; + jintArray jia; + int len; + jbyte *jbBody; + jshort *jsBody; + jint *jiBody; + jboolean isCopy; + + /* Image width, height and depth*/ + fid = (*env)->GetFieldID(env, cls, "width", "I"); + ji = (*env)->GetIntField(env, obj, fid); + w = ji; + + fid = (*env)->GetFieldID(env, cls, "height", "I"); + ji = (*env)->GetIntField(env, obj, fid); + h = ji; + + fid = (*env)->GetFieldID(env, cls, "depth", "I"); + ji = (*env)->GetIntField(env, obj, fid); + depth = ji; + + /* Read the image*/ + if (depth <= 16) { + numcomps = 1; + color_space = CLRSPC_GRAY; + } else { + numcomps = 3; + color_space = CLRSPC_SRGB; + } + memset(&cmptparm[0], 0, numcomps * sizeof(opj_image_cmptparm_t)); + + if (numcomps == 1) { + cmptparm[0].x0 = parameters->image_offset_x0; + cmptparm[0].y0 = parameters->image_offset_y0; + cmptparm[0].w = !cmptparm[0].x0 ? (w - 1) * parameters->subsampling_dx + 1 : + cmptparm[0].x0 + (w - 1) * parameters->subsampling_dx + 1; + cmptparm[0].h = !cmptparm[0].y0 ? (h - 1) * parameters->subsampling_dy + 1 : + cmptparm[0].y0 + (h - 1) * parameters->subsampling_dy + 1; + /* Java types are always signed but we use them as unsigned types (shift of the negative part of + the pixels of the images in Telemis before entering the encoder).*/ + cmptparm[0].sgnd = 0; + if (depth <= 16) { + cmptparm[0].prec = depth; + } else { + cmptparm[0].prec = 8; + } + cmptparm[0].bpp = cmptparm[0].prec; + cmptparm[0].dx = parameters->subsampling_dx; + cmptparm[0].dy = parameters->subsampling_dy; + /*printf("C: component 0 initialised: x0=%d, y0=%d, w=%d, h=%d, sgnd=%d, bpp=%d, dx=%d, dy=%d, color_space=%d\n", cmptparm[0].x0, cmptparm[0].y0, cmptparm[0].w, + cmptparm[0].h, cmptparm[0].sgnd, cmptparm[0].bpp, cmptparm[0].dx, cmptparm[0].dy, color_space);*/ + } else { + for (i = 0; i < numcomps; i++) { + cmptparm[i].prec = 8; + cmptparm[i].bpp = 8; + cmptparm[i].sgnd = 0; + cmptparm[i].dx = parameters->subsampling_dx; + cmptparm[i].dy = parameters->subsampling_dy; + cmptparm[i].w = w; + cmptparm[i].h = h; + } + } + + /* create the image */ + image = opj_image_create(numcomps, &cmptparm[0], color_space); + + if (!image) { + return NULL; + } + + if (depth <= 16) { + image->numcomps = 1; + } else { + image->numcomps = 3; + } + + /* set image offset and reference grid */ + image->x0 = cmptparm[0].x0; + image->y0 = cmptparm[0].x0; + image->x1 = cmptparm[0].w; + image->y1 = cmptparm[0].h; + + /* set image data */ + for (compno = 0; compno < numcomps; compno++) { + comp = &image->comps[compno]; + max = -100000; + if (depth == 8) { + fid = (*env)->GetFieldID(env, cls, "image8", "[B"); /* byteArray []*/ + jba = (*env)->GetObjectField(env, obj, fid); + len = (*env)->GetArrayLength(env, jba); + + jbBody = (*env)->GetPrimitiveArrayCritical(env, jba, &isCopy); + /*printf("C: before transferring 8 bpp image\n");*/ + if (comp->sgnd) { + for (i = 0; i < len; i++) { + comp->data[i] = (char) jbBody[i]; + if (comp->data[i] > max) { + max = comp->data[i]; + } + } + } else { + for (i = 0; i < len; i++) { + comp->data[i] = (unsigned char) jbBody[i]; + if (comp->data[i] > max) { + max = comp->data[i]; + } + } + } + (*env)->ReleasePrimitiveArrayCritical(env, jba, jbBody, 0); + } else if (depth == 16) { + fid = (*env)->GetFieldID(env, cls, "image16", "[S"); /* shortArray []*/ + jsa = (*env)->GetObjectField(env, obj, fid); + len = (*env)->GetArrayLength(env, jsa); + + jsBody = (*env)->GetPrimitiveArrayCritical(env, jsa, &isCopy); + /*printf("C: before transferring 16 bpp image\n");*/ + if (comp->sgnd) { /* Special behaviour to deal with signed elements ??*/ + comp->data[i] = (short) jsBody[i]; + for (i = 0; i < len; i++) { + if (comp->data[i] > max) { + max = comp->data[i]; + } + } + } else { + for (i = 0; i < len; i++) { + comp->data[i] = (unsigned short) jsBody[i]; + if (comp->data[i] > max) { + max = comp->data[i]; + } + } + } + (*env)->ReleasePrimitiveArrayCritical(env, jsa, jsBody, 0); + } else if (depth == 24) { + fid = (*env)->GetFieldID(env, cls, "image24", "[I"); /* intArray []*/ + jia = (*env)->GetObjectField(env, obj, fid); + len = (*env)->GetArrayLength(env, jia); + shift = compno * 8; + + jiBody = (*env)->GetPrimitiveArrayCritical(env, jia, &isCopy); + /*printf("C: before transferring 24 bpp image (component %d, signed = %d)\n", compno, comp->sgnd);*/ + if (comp->sgnd) { /* Special behaviour to deal with signed elements ?? XXXXX*/ + for (i = 0; i < len; i++) { + comp->data[i] = (((int) jiBody[i]) & (0xFF << shift)) >> shift; + if (comp->data[i] > max) { + max = comp->data[i]; + } + } + } else { + for (i = 0; i < len; i++) { + comp->data[i] = (((unsigned int) jiBody[i]) & (0xFF << shift)) >> shift; + if (comp->data[i] > max) { + max = comp->data[i]; + } + } + } + (*env)->ReleasePrimitiveArrayCritical(env, jia, jiBody, 0); + } + comp->bpp = int_floorlog2(max) + 1; + comp->prec = comp->bpp; + /*printf("C: component %d: max %d, real bpp = %d\n", compno, max, comp->bpp);*/ + } + return image; +} + + +/* -------------------------------------------------------------------------- + -------------------- MAIN METHOD, CALLED BY JAVA -----------------------*/ +JNIEXPORT jlong JNICALL +Java_org_openJpeg_OpenJPEGJavaEncoder_internalEncodeImageToJ2K(JNIEnv *env, + jobject obj, jobjectArray javaParameters) +{ + int argc; /* To simulate the command line parameters (taken from the javaParameters variable) and be able to re-use the */ + char **argv; /* 'parse_cmdline_decoder' method taken from the j2k_to_image project */ + opj_bool bSuccess; + opj_cparameters_t parameters; /* compression parameters */ + img_fol_t img_fol; + opj_event_mgr_t event_mgr; /* event manager */ + opj_image_t *image = NULL; + int i, j, num_images; + int imageno; + opj_codestream_info_t cstr_info; /* Codestream information structure */ + char indexfilename[OPJ_PATH_LEN]; /* index file name */ + + char* compressed_index = NULL; + int compressed_index_size = -1; + /* ==> Access variables to the Java member variables*/ + jsize arraySize; + jclass cls; + jobject object; + jboolean isCopy; + jfieldID fid; + jbyteArray jba; + jbyte *jbBody; + callback_variables_t msgErrorCallback_vars; + /* <== access variable to the Java member variables.*/ + + /* For the encoding and storage into the file*/ + opj_cinfo_t* cinfo; + int codestream_length; + opj_cio_t *cio = NULL; + FILE *f = NULL; + + /* JNI reference to the calling class*/ + cls = (*env)->GetObjectClass(env, obj); + + /* Pointers to be able to call a Java method for all the info and error messages*/ + msgErrorCallback_vars.env = env; + msgErrorCallback_vars.jobj = &obj; + msgErrorCallback_vars.message_mid = (*env)->GetMethodID(env, cls, "logMessage", + "(Ljava/lang/String;)V"); + msgErrorCallback_vars.error_mid = (*env)->GetMethodID(env, cls, "logError", + "(Ljava/lang/String;)V"); + + arraySize = (*env)->GetArrayLength(env, javaParameters); + argc = (int) arraySize + 1; + argv = opj_malloc(argc * sizeof(char*)); + argv[0] = "ProgramName.exe"; /* The program name: useless*/ + j = 0; + for (i = 1; i < argc; i++) { + object = (*env)->GetObjectArrayElement(env, javaParameters, i - 1); + argv[i] = (char*)(*env)->GetStringUTFChars(env, object, &isCopy); + } + + /*printf("C: "); + for (i=0; iReleaseStringUTFChars(env, (*env)->GetObjectArrayElement(env, + javaParameters, i - 1), argv[i]); + } + return -1; + } + + /* Release the Java arguments array*/ + for (i = 1; i < argc; i++) { + (*env)->ReleaseStringUTFChars(env, (*env)->GetObjectArrayElement(env, + javaParameters, i - 1), argv[i]); + } + + if (parameters.cp_cinema) { + cinema_parameters(¶meters); + } + + + /* Create comment for codestream */ + if (parameters.cp_comment == NULL) { + const char comment[] = "Created by JavaOpenJPEG version "; + const size_t clen = strlen(comment); + const char *version = opj_version(); + /* UniPG>> */ +#ifdef USE_JPWL + parameters.cp_comment = (char*)opj_malloc(clen + strlen(version) + 11); + sprintf(parameters.cp_comment, "%s%s with JPWL", comment, version); +#else + parameters.cp_comment = (char*)opj_malloc(clen + strlen(version) + 1); + sprintf(parameters.cp_comment, "%s%s", comment, version); +#endif + /* <numcomps == 3 ? 1 : 0; + + if (parameters.cp_cinema) { + cinema_setup_encoder(¶meters, image, &img_fol); + } + + /* encode the destination image */ + /* ---------------------------- */ + /* get a J2K compressor handle */ + if (parameters.cod_format == J2K_CFMT) { /* J2K format output */ + cinfo = opj_create_compress(CODEC_J2K); + } else { /* JP2 format output */ + cinfo = opj_create_compress(CODEC_JP2); + } + /* catch events using our callbacks and give a local context */ + opj_set_event_mgr((opj_common_ptr)cinfo, &event_mgr, &msgErrorCallback_vars); + + /* setup the encoder parameters using the current image and user parameters */ + opj_setup_encoder(cinfo, ¶meters, image); + + /* open a byte stream for writing */ + /* allocate memory for all tiles */ + cio = opj_cio_open((opj_common_ptr)cinfo, NULL, 0); + + /* encode the image */ + bSuccess = opj_encode_with_info(cinfo, cio, image, &cstr_info); + printf("C: after opj_encode_with_info\n"); + if (!bSuccess) { + opj_cio_close(cio); + fprintf(stderr, "failed to encode image\n"); + return -1; + } + codestream_length = cio_tell(cio); + + /* write the index on disk, if needed (-x 'filename') */ + if (*indexfilename) { + bSuccess = write_index_file(&cstr_info, indexfilename); + if (bSuccess) { + fprintf(stderr, "Failed to output index file into [%s]\n", indexfilename); + } + } + + compressed_index = create_index_into_byte_array(&cstr_info, + &compressed_index_size); + /* Allocates the Java compressedIndex byte[] and sends this index into the Java object */ + fid = (*env)->GetFieldID(env, cls, "compressedIndex", "[B"); + jba = (*env)->NewByteArray(env, compressed_index_size + 1); + jbBody = (*env)->GetPrimitiveArrayCritical(env, jba, 0); + memcpy(jbBody, compressed_index, compressed_index_size); + (*env)->ReleasePrimitiveArrayCritical(env, jba, jbBody, 0); + (*env)->SetObjectField(env, obj, fid, jba); + opj_free(compressed_index); + + /* write the generated codestream to disk ? */ + if (parameters.outfile[0] != '\0') { + f = fopen(parameters.outfile, "wb"); + if (!f) { + fprintf(stderr, "failed to open [%s] for writing\n", parameters.outfile); + return -1; + } + fwrite(cio->buffer, 1, codestream_length, f); + fclose(f); + fprintf(stdout, "Generated outfile [%s]\n", parameters.outfile); + } + + /* Write the generated codestream to the Java pre-allocated compressedStream byte[] */ + fid = (*env)->GetFieldID(env, cls, "compressedStream", "[B"); + jba = (*env)->GetObjectField(env, obj, fid); + jbBody = (*env)->GetPrimitiveArrayCritical(env, jba, 0); + memcpy(jbBody, cio->buffer, codestream_length); + (*env)->ReleasePrimitiveArrayCritical(env, jba, jbBody, 0); + + /* close and free the byte stream */ + opj_cio_close(cio); + + /* free remaining compression structures */ + opj_destroy_compress(cinfo); + opj_destroy_cstr_info(&cstr_info); + + /* free image data */ + opj_image_destroy(image); + } + + /* free user parameters structure */ + if (parameters.cp_comment) { + opj_free(parameters.cp_comment); + } + if (parameters.cp_matrice) { + opj_free(parameters.cp_matrice); + } + + return codestream_length; +} + -- cgit v1.2.3