/* * Copyright (c) 1999-2000 Image Power, Inc. and the University of * British Columbia. * Copyright (c) 2001-2003 Michael David Adams. * All rights reserved. */ /* __START_OF_JASPER_LICENSE__ * * JasPer License Version 2.0 * * Copyright (c) 2001-2006 Michael David Adams * Copyright (c) 1999-2000 Image Power, Inc. * Copyright (c) 1999-2000 The University of British Columbia * * All rights reserved. * * Permission is hereby granted, free of charge, to any person (the * "User") obtaining a copy of this software and associated documentation * files (the "Software"), to deal in the Software without restriction, * including without limitation the rights to use, copy, modify, merge, * publish, distribute, and/or sell copies of the Software, and to permit * persons to whom the Software is furnished to do so, subject to the * following conditions: * * 1. The above copyright notices and this permission notice (which * includes the disclaimer below) shall be included in all copies or * substantial portions of the Software. * * 2. The name of a copyright holder shall not be used to endorse or * promote products derived from the Software without specific prior * written permission. * * THIS DISCLAIMER OF WARRANTY CONSTITUTES AN ESSENTIAL PART OF THIS * LICENSE. NO USE OF THE SOFTWARE IS AUTHORIZED HEREUNDER EXCEPT UNDER * THIS DISCLAIMER. THE SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS * "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING * BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A * PARTICULAR PURPOSE AND NONINFRINGEMENT OF THIRD PARTY RIGHTS. IN NO * EVENT SHALL THE COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, OR ANY SPECIAL * INDIRECT OR CONSEQUENTIAL DAMAGES, OR ANY DAMAGES WHATSOEVER RESULTING * FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN ACTION OF CONTRACT, * NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF OR IN CONNECTION * WITH THE USE OR PERFORMANCE OF THIS SOFTWARE. NO ASSURANCES ARE * PROVIDED BY THE COPYRIGHT HOLDERS THAT THE SOFTWARE DOES NOT INFRINGE * THE PATENT OR OTHER INTELLECTUAL PROPERTY RIGHTS OF ANY OTHER ENTITY. * EACH COPYRIGHT HOLDER DISCLAIMS ANY LIABILITY TO THE USER FOR CLAIMS * BROUGHT BY ANY OTHER ENTITY BASED ON INFRINGEMENT OF INTELLECTUAL * PROPERTY RIGHTS OR OTHERWISE. AS A CONDITION TO EXERCISING THE RIGHTS * GRANTED HEREUNDER, EACH USER HEREBY ASSUMES SOLE RESPONSIBILITY TO SECURE * ANY OTHER INTELLECTUAL PROPERTY RIGHTS NEEDED, IF ANY. THE SOFTWARE * IS NOT FAULT-TOLERANT AND IS NOT INTENDED FOR USE IN MISSION-CRITICAL * SYSTEMS, SUCH AS THOSE USED IN THE OPERATION OF NUCLEAR FACILITIES, * AIRCRAFT NAVIGATION OR COMMUNICATION SYSTEMS, AIR TRAFFIC CONTROL * SYSTEMS, DIRECT LIFE SUPPORT MACHINES, OR WEAPONS SYSTEMS, IN WHICH * THE FAILURE OF THE SOFTWARE OR SYSTEM COULD LEAD DIRECTLY TO DEATH, * PERSONAL INJURY, OR SEVERE PHYSICAL OR ENVIRONMENTAL DAMAGE ("HIGH * RISK ACTIVITIES"). THE COPYRIGHT HOLDERS SPECIFICALLY DISCLAIM ANY * EXPRESS OR IMPLIED WARRANTY OF FITNESS FOR HIGH RISK ACTIVITIES. * * __END_OF_JASPER_LICENSE__ */ /* * $Id$ */ /******************************************************************************\* Includes.\******************************************************************************/ #include "jpc_enc.h"#include "jpc_flt.h"#include "jpc_fix.h"#include "jpc_tagtree.h"#include "jpc_cs.h"#include "jpc_mct.h"#include "jpc_tsfb.h"#include "jpc_t1cod.h"#include "jpc_t1enc.h"#include "jpc_t2enc.h"#include "jpc_math.h"#include "jpc_util.h" #include "jasper/jas_types.h"#include "jasper/jas_string.h"#include "jasper/jas_malloc.h"#include "jasper/jas_image.h"#include "jasper/jas_tvp.h"#include "jasper/jas_version.h"#include "jasper/jas_math.h"#include "jasper/jas_debug.h" /******************************************************************************\*\******************************************************************************/ #define JPC_POW2(n) \ (1 << (n)) #define JPC_FLOORTOMULTPOW2(x, n) \ (((n) > 0) ? ((x) & (~((1 << n) - 1))) : (x))/* Round to the nearest multiple of the specified power of two in the direction of negative infinity. */ #define JPC_CEILTOMULTPOW2(x, n) \ (((n) > 0) ? JPC_FLOORTOMULTPOW2(((x) + (1 << (n)) - 1), n) : (x))/* Round to the nearest multiple of the specified power of two in the direction of positive infinity. */ #define JPC_POW2(n) \ (1 << (n)) /******************************************************************************\* Local prototypes.\******************************************************************************/ static jpc_enc_tile_t *jpc_enc_tile_create(jpc_enc_cp_t *cp, jas_image_t *image, int tileno);static void jpc_enc_tile_destroy(jpc_enc_tile_t *tile); static jpc_enc_tcmpt_t *tcmpt_create(jpc_enc_tcmpt_t *tcmpt, jpc_enc_cp_t *cp, jas_image_t *image, jpc_enc_tile_t *tile);static void tcmpt_destroy(jpc_enc_tcmpt_t *tcmpt);static jpc_enc_rlvl_t *rlvl_create(jpc_enc_rlvl_t *rlvl, jpc_enc_cp_t *cp, jpc_enc_tcmpt_t *tcmpt, jpc_tsfb_band_t *bandinfos);static void rlvl_destroy(jpc_enc_rlvl_t *rlvl);static jpc_enc_band_t *band_create(jpc_enc_band_t *band, jpc_enc_cp_t *cp, jpc_enc_rlvl_t *rlvl, jpc_tsfb_band_t *bandinfos);static void band_destroy(jpc_enc_band_t *bands);static jpc_enc_prc_t *prc_create(jpc_enc_prc_t *prc, jpc_enc_band_t *band);static void prc_destroy(jpc_enc_prc_t *prcs);static jpc_enc_cblk_t *cblk_create(jpc_enc_cblk_t *cblk, jpc_enc_prc_t *prc);static void cblk_destroy(jpc_enc_cblk_t *cblks);int ratestrtosize(const char *s, uint_fast32_t rawsize, uint_fast32_t *size);static void pass_destroy(jpc_enc_pass_t *pass);void jpc_enc_dump(jpc_enc_t *enc); void calcrdslopes(jpc_enc_cblk_t *cblk);void dump_layeringinfo(jpc_enc_t *enc);void jpc_quantize(jas_matrix_t *data, jpc_fix_t stepsize);jpc_enc_t *jpc_enc_create(jpc_enc_cp_t *cp, jas_stream_t *out, jas_image_t *image);void jpc_enc_destroy(jpc_enc_t *enc);static int jpc_enc_encodemainhdr(jpc_enc_t *enc);static int jpc_enc_encodemainbody(jpc_enc_t *enc);int jpc_enc_encodetiledata(jpc_enc_t *enc);int rateallocate(jpc_enc_t *enc, unsigned numlyrs, uint_fast32_t *cumlens);static jpc_enc_cp_t *cp_create(const char *optstr, jas_image_t *image);void jpc_enc_cp_destroy(jpc_enc_cp_t *cp);static uint_fast32_t jpc_abstorelstepsize(jpc_fix_t absdelta, int scaleexpn); /******************************************************************************\* Types.\******************************************************************************/ typedef enum { OPT_DEBUG, OPT_IMGAREAOFFX, OPT_IMGAREAOFFY, OPT_TILEGRDOFFX, OPT_TILEGRDOFFY, OPT_TILEWIDTH, OPT_TILEHEIGHT, OPT_PRCWIDTH, OPT_PRCHEIGHT, OPT_CBLKWIDTH, OPT_CBLKHEIGHT, OPT_MODE, OPT_PRG, OPT_NOMCT, OPT_MAXRLVLS, OPT_SOP, OPT_EPH, OPT_LAZY, OPT_TERMALL, OPT_SEGSYM, OPT_VCAUSAL, OPT_RESET, OPT_PTERM, OPT_NUMGBITS, OPT_RATE, OPT_ILYRRATES, OPT_JP2OVERHEAD} optid_t; typedef enum { PO_L = 0, PO_R} poid_t; typedef enum { MODE_INT, MODE_REAL} modeid_t; /******************************************************************************\* Data.\******************************************************************************/ static const jas_taginfo_t encopts[] = { {OPT_DEBUG, "debug"}, {OPT_IMGAREAOFFX, "imgareatlx"}, {OPT_IMGAREAOFFY, "imgareatly"}, {OPT_TILEGRDOFFX, "tilegrdtlx"}, {OPT_TILEGRDOFFY, "tilegrdtly"}, {OPT_TILEWIDTH, "tilewidth"}, {OPT_TILEHEIGHT, "tileheight"}, {OPT_PRCWIDTH, "prcwidth"}, {OPT_PRCHEIGHT, "prcheight"}, {OPT_CBLKWIDTH, "cblkwidth"}, {OPT_CBLKHEIGHT, "cblkheight"}, {OPT_MODE, "mode"}, {OPT_PRG, "prg"}, {OPT_NOMCT, "nomct"}, {OPT_MAXRLVLS, "numrlvls"}, {OPT_SOP, "sop"}, {OPT_EPH, "eph"}, {OPT_LAZY, "lazy"}, {OPT_TERMALL, "termall"}, {OPT_SEGSYM, "segsym"}, {OPT_VCAUSAL, "vcausal"}, {OPT_PTERM, "pterm"}, {OPT_RESET, "resetprob"}, {OPT_NUMGBITS, "numgbits"}, {OPT_RATE, "rate"}, {OPT_ILYRRATES, "ilyrrates"}, {OPT_JP2OVERHEAD, "_jp2overhead"}, {-1, 0}}; static const jas_taginfo_t prgordtab[] = { {JPC_COD_LRCPPRG, "lrcp"}, {JPC_COD_RLCPPRG, "rlcp"}, {JPC_COD_RPCLPRG, "rpcl"}, {JPC_COD_PCRLPRG, "pcrl"}, {JPC_COD_CPRLPRG, "cprl"}, {-1, 0}}; static const jas_taginfo_t modetab[] = { {MODE_INT, "int"}, {MODE_REAL, "real"}, {-1, 0}}; /******************************************************************************\\******************************************************************************/ /** * @return UINT_FAST32_MAX on error */static uint_fast32_t jpc_abstorelstepsize(jpc_fix_t absdelta, int scaleexpn){ int p; uint_fast32_t mant; uint_fast32_t expn; int n; if (absdelta < 0) { return UINT_FAST32_MAX; } p = jpc_fix_firstone(absdelta) - JPC_FIX_FRACBITS; n = 11 - jpc_fix_firstone(absdelta); mant = ((n < 0) ? (absdelta >> (-n)) : (absdelta << n)) & 0x7ff; expn = scaleexpn - p; if (scaleexpn < p) { return UINT_FAST32_MAX; } if (expn >= 0x1f) return UINT_FAST32_MAX; return JPC_QCX_EXPN(expn) | JPC_QCX_MANT(mant);} /******************************************************************************\* The main encoder entry point.\******************************************************************************/ int jpc_encode(jas_image_t *image, jas_stream_t *out, const char *optstr){ jpc_enc_t *enc; jpc_enc_cp_t *cp; enc = 0; jpc_init(); if (!(cp = cp_create(optstr, image))) { jas_logerrorf("invalid JP encoder options\n"); goto error; } if (!(enc = jpc_enc_create(cp, out, image))) { jas_logerrorf("jpc_enc_create failed\n"); goto error; } cp = 0; /* Encode the main header. */ if (jpc_enc_encodemainhdr(enc)) { jas_logerrorf("cannot encode main header\n"); goto error; } /* Encode the main body. This constitutes most of the encoding work. */ if (jpc_enc_encodemainbody(enc)) { jas_logerrorf("cannot encode main body\n"); goto error; } /* Write EOC marker segment. */ if (!(enc->mrk = jpc_ms_create(JPC_MS_EOC))) { jas_logerrorf("cannot create EOC marker\n"); goto error; } if (jpc_putms(enc->out, enc->cstate, enc->mrk)) { jas_logerrorf("cannot write EOC marker\n"); goto error; } jpc_ms_destroy(enc->mrk); enc->mrk = 0; if (jas_stream_flush(enc->out)) { jas_logerrorf("stream flush failed\n"); goto error; } jpc_enc_destroy(enc); return 0; error: if (cp) { jpc_enc_cp_destroy(cp); } if (enc) { jpc_enc_destroy(enc); } jas_logerrorf("jpc_encode failed\n"); return -1;} /******************************************************************************\* Option parsing code.\******************************************************************************/ static jpc_enc_cp_t *cp_create(const char *optstr, jas_image_t *image){ jpc_enc_cp_t *cp; jas_tvparser_t *tvp; int numilyrrates; double *ilyrrates; int i; int tagid; jpc_enc_tcp_t *tcp; jpc_enc_tccp_t *tccp; jpc_enc_ccp_t *ccp; uint_fast16_t rlvlno; uint_fast16_t prcwidthexpn; uint_fast16_t prcheightexpn; bool enablemct; uint_fast32_t jp2overhead; uint_fast16_t lyrno; uint_fast32_t hsteplcm; uint_fast32_t vsteplcm; bool mctvalid; tvp = 0; ilyrrates = 0; numilyrrates = 0; if (!(cp = jas_malloc(sizeof(jpc_enc_cp_t)))) { goto error; } prcwidthexpn = 15; prcheightexpn = 15; enablemct = true; jp2overhead = 0; cp->ccps = 0; cp->debug = 0; cp->imgareatlx = UINT_FAST32_MAX; cp->imgareatly = UINT_FAST32_MAX; cp->refgrdwidth = 0; cp->refgrdheight = 0; cp->tilegrdoffx = UINT_FAST32_MAX; cp->tilegrdoffy = UINT_FAST32_MAX; cp->tilewidth = 0; cp->tileheight = 0; cp->numcmpts = jas_image_numcmpts(image); cp->tcp.ilyrrates = NULL; hsteplcm = 1; vsteplcm = 1; for (unsigned cmptno = 0; cmptno < jas_image_numcmpts(image); ++cmptno) { if (jas_image_cmptbrx(image, cmptno) + jas_image_cmpthstep(image, cmptno) <= jas_image_brx(image) || jas_image_cmptbry(image, cmptno) + jas_image_cmptvstep(image, cmptno) <= jas_image_bry(image)) { jas_logerrorf("unsupported image type\n"); goto error; } /* Note: We ought to be calculating the LCMs here. Fix some day. */ hsteplcm *= jas_image_cmpthstep(image, cmptno); vsteplcm *= jas_image_cmptvstep(image, cmptno); } if (!(cp->ccps = jas_alloc2(cp->numcmpts, sizeof(jpc_enc_ccp_t)))) { goto error; } unsigned cmptno; for (cmptno = 0, ccp = cp->ccps; cmptno < cp->numcmpts; ++cmptno, ++ccp) { ccp->sampgrdstepx = jas_image_cmpthstep(image, cmptno); ccp->sampgrdstepy = jas_image_cmptvstep(image, cmptno); /* XXX - this isn't quite correct for more general image */ ccp->sampgrdsubstepx = 0; ccp->sampgrdsubstepx = 0; ccp->prec = jas_image_cmptprec(image, cmptno); ccp->sgnd = jas_image_cmptsgnd(image, cmptno); ccp->numstepsizes = 0; memset(ccp->stepsizes, 0, sizeof(ccp->stepsizes)); } cp->rawsize = jas_image_rawsize(image); if (cp->rawsize == 0) { /* prevent division by zero in cp_create() */ goto error; } cp->totalsize = UINT_FAST32_MAX; tcp = &cp->tcp; tcp->csty = 0; tcp->intmode = true; tcp->prg = JPC_COD_LRCPPRG; tcp->numlyrs = 1; tcp->ilyrrates = 0; tccp = &cp->tccp; tccp->csty = 0; tccp->maxrlvls = 6; tccp->cblkwidthexpn = 6; tccp->cblkheightexpn = 6; tccp->cblksty = 0; tccp->numgbits = 2; if (!(tvp = jas_tvparser_create(optstr ? optstr : ""))) { goto error; } while (!jas_tvparser_next(tvp)) { switch (jas_taginfo_nonull(jas_taginfos_lookup(encopts, jas_tvparser_gettag(tvp)))->id) { case OPT_DEBUG: cp->debug = atoi(jas_tvparser_getval(tvp)); break; case OPT_IMGAREAOFFX: cp->imgareatlx = atoi(jas_tvparser_getval(tvp)); break; case OPT_IMGAREAOFFY: cp->imgareatly = atoi(jas_tvparser_getval(tvp)); break; case OPT_TILEGRDOFFX: cp->tilegrdoffx = atoi(jas_tvparser_getval(tvp)); break; case OPT_TILEGRDOFFY: cp->tilegrdoffy = atoi(jas_tvparser_getval(tvp)); break; case OPT_TILEWIDTH: cp->tilewidth = atoi(jas_tvparser_getval(tvp)); break; case OPT_TILEHEIGHT: cp->tileheight = atoi(jas_tvparser_getval(tvp)); break; case OPT_PRCWIDTH: prcwidthexpn = jpc_floorlog2(atoi(jas_tvparser_getval(tvp))); break; case OPT_PRCHEIGHT: prcheightexpn = jpc_floorlog2(atoi(jas_tvparser_getval(tvp))); break; case OPT_CBLKWIDTH: tccp->cblkwidthexpn = jpc_floorlog2(atoi(jas_tvparser_getval(tvp))); break; case OPT_CBLKHEIGHT: tccp->cblkheightexpn = jpc_floorlog2(atoi(jas_tvparser_getval(tvp))); break; case OPT_MODE: if ((tagid = jas_taginfo_nonull(jas_taginfos_lookup(modetab, jas_tvparser_getval(tvp)))->id) < 0) { jas_logwarnf("ignoring invalid mode %s\n", jas_tvparser_getval(tvp)); } else { tcp->intmode = (tagid == MODE_INT); } break; case OPT_PRG: if ((tagid = jas_taginfo_nonull(jas_taginfos_lookup(prgordtab,Lines 1–500 of 2665. Load more
/* * Copyright (c) 1999-2000 Image Power, Inc. and the University of * British Columbia. * Copyright (c) 2001-2003 Michael David Adams. * All rights reserved. */ /* __START_OF_JASPER_LICENSE__ * * JasPer License Version 2.0 * * Copyright (c) 2001-2006 Michael David Adams * Copyright (c) 1999-2000 Image Power, Inc. * Copyright (c) 1999-2000 The University of British Columbia * * All rights reserved. * * Permission is hereby granted, free of charge, to any person (the * "User") obtaining a copy of this software and associated documentation * files (the "Software"), to deal in the Software without restriction, * including without limitation the rights to use, copy, modify, merge, * publish, distribute, and/or sell copies of the Software, and to permit * persons to whom the Software is furnished to do so, subject to the * following conditions: * * 1. The above copyright notices and this permission notice (which * includes the disclaimer below) shall be included in all copies or * substantial portions of the Software. * * 2. The name of a copyright holder shall not be used to endorse or * promote products derived from the Software without specific prior * written permission. * * THIS DISCLAIMER OF WARRANTY CONSTITUTES AN ESSENTIAL PART OF THIS * LICENSE. NO USE OF THE SOFTWARE IS AUTHORIZED HEREUNDER EXCEPT UNDER * THIS DISCLAIMER. THE SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS * "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING * BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A * PARTICULAR PURPOSE AND NONINFRINGEMENT OF THIRD PARTY RIGHTS. IN NO * EVENT SHALL THE COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, OR ANY SPECIAL * INDIRECT OR CONSEQUENTIAL DAMAGES, OR ANY DAMAGES WHATSOEVER RESULTING * FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN ACTION OF CONTRACT, * NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF OR IN CONNECTION * WITH THE USE OR PERFORMANCE OF THIS SOFTWARE. NO ASSURANCES ARE * PROVIDED BY THE COPYRIGHT HOLDERS THAT THE SOFTWARE DOES NOT INFRINGE * THE PATENT OR OTHER INTELLECTUAL PROPERTY RIGHTS OF ANY OTHER ENTITY. * EACH COPYRIGHT HOLDER DISCLAIMS ANY LIABILITY TO THE USER FOR CLAIMS * BROUGHT BY ANY OTHER ENTITY BASED ON INFRINGEMENT OF INTELLECTUAL * PROPERTY RIGHTS OR OTHERWISE. AS A CONDITION TO EXERCISING THE RIGHTS * GRANTED HEREUNDER, EACH USER HEREBY ASSUMES SOLE RESPONSIBILITY TO SECURE * ANY OTHER INTELLECTUAL PROPERTY RIGHTS NEEDED, IF ANY. THE SOFTWARE * IS NOT FAULT-TOLERANT AND IS NOT INTENDED FOR USE IN MISSION-CRITICAL * SYSTEMS, SUCH AS THOSE USED IN THE OPERATION OF NUCLEAR FACILITIES, * AIRCRAFT NAVIGATION OR COMMUNICATION SYSTEMS, AIR TRAFFIC CONTROL * SYSTEMS, DIRECT LIFE SUPPORT MACHINES, OR WEAPONS SYSTEMS, IN WHICH * THE FAILURE OF THE SOFTWARE OR SYSTEM COULD LEAD DIRECTLY TO DEATH, * PERSONAL INJURY, OR SEVERE PHYSICAL OR ENVIRONMENTAL DAMAGE ("HIGH * RISK ACTIVITIES"). THE COPYRIGHT HOLDERS SPECIFICALLY DISCLAIM ANY * EXPRESS OR IMPLIED WARRANTY OF FITNESS FOR HIGH RISK ACTIVITIES. * * __END_OF_JASPER_LICENSE__ */ /* * $Id$ */ /******************************************************************************\ * Includes. \******************************************************************************/ #include "jpc_enc.h" #include "jpc_flt.h" #include "jpc_fix.h" #include "jpc_tagtree.h" #include "jpc_cs.h" #include "jpc_mct.h" #include "jpc_tsfb.h" #include "jpc_t1cod.h" #include "jpc_t1enc.h" #include "jpc_t2enc.h" #include "jpc_math.h" #include "jpc_util.h" #include "jasper/jas_types.h" #include "jasper/jas_string.h" #include "jasper/jas_malloc.h" #include "jasper/jas_image.h" #include "jasper/jas_tvp.h" #include "jasper/jas_version.h" #include "jasper/jas_math.h" #include "jasper/jas_debug.h" /******************************************************************************\ * \******************************************************************************/ #define JPC_POW2(n) \ (1 << (n)) #define JPC_FLOORTOMULTPOW2(x, n) \ (((n) > 0) ? ((x) & (~((1 << n) - 1))) : (x)) /* Round to the nearest multiple of the specified power of two in the direction of negative infinity. */ #define JPC_CEILTOMULTPOW2(x, n) \ (((n) > 0) ? JPC_FLOORTOMULTPOW2(((x) + (1 << (n)) - 1), n) : (x)) /* Round to the nearest multiple of the specified power of two in the direction of positive infinity. */ #define JPC_POW2(n) \ (1 << (n)) /******************************************************************************\ * Local prototypes. \******************************************************************************/ static jpc_enc_tile_t *jpc_enc_tile_create(jpc_enc_cp_t *cp, jas_image_t *image, int tileno); static void jpc_enc_tile_destroy(jpc_enc_tile_t *tile); static jpc_enc_tcmpt_t *tcmpt_create(jpc_enc_tcmpt_t *tcmpt, jpc_enc_cp_t *cp, jas_image_t *image, jpc_enc_tile_t *tile); static void tcmpt_destroy(jpc_enc_tcmpt_t *tcmpt); static jpc_enc_rlvl_t *rlvl_create(jpc_enc_rlvl_t *rlvl, jpc_enc_cp_t *cp, jpc_enc_tcmpt_t *tcmpt, jpc_tsfb_band_t *bandinfos); static void rlvl_destroy(jpc_enc_rlvl_t *rlvl); static jpc_enc_band_t *band_create(jpc_enc_band_t *band, jpc_enc_cp_t *cp, jpc_enc_rlvl_t *rlvl, jpc_tsfb_band_t *bandinfos); static void band_destroy(jpc_enc_band_t *bands); static jpc_enc_prc_t *prc_create(jpc_enc_prc_t *prc, jpc_enc_band_t *band); static void prc_destroy(jpc_enc_prc_t *prcs); static jpc_enc_cblk_t *cblk_create(jpc_enc_cblk_t *cblk, jpc_enc_prc_t *prc); static void cblk_destroy(jpc_enc_cblk_t *cblks); int ratestrtosize(const char *s, uint_fast32_t rawsize, uint_fast32_t *size); static void pass_destroy(jpc_enc_pass_t *pass); void jpc_enc_dump(jpc_enc_t *enc); void calcrdslopes(jpc_enc_cblk_t *cblk); void dump_layeringinfo(jpc_enc_t *enc); void jpc_quantize(jas_matrix_t *data, jpc_fix_t stepsize); jpc_enc_t *jpc_enc_create(jpc_enc_cp_t *cp, jas_stream_t *out, jas_image_t *image); void jpc_enc_destroy(jpc_enc_t *enc); static int jpc_enc_encodemainhdr(jpc_enc_t *enc); static int jpc_enc_encodemainbody(jpc_enc_t *enc); int jpc_enc_encodetiledata(jpc_enc_t *enc); int rateallocate(jpc_enc_t *enc, unsigned numlyrs, uint_fast32_t *cumlens); static jpc_enc_cp_t *cp_create(const char *optstr, jas_image_t *image); void jpc_enc_cp_destroy(jpc_enc_cp_t *cp); static uint_fast32_t jpc_abstorelstepsize(jpc_fix_t absdelta, int scaleexpn); /******************************************************************************\ * Types. \******************************************************************************/ typedef enum { OPT_DEBUG, OPT_IMGAREAOFFX, OPT_IMGAREAOFFY, OPT_TILEGRDOFFX, OPT_TILEGRDOFFY, OPT_TILEWIDTH, OPT_TILEHEIGHT, OPT_PRCWIDTH, OPT_PRCHEIGHT, OPT_CBLKWIDTH, OPT_CBLKHEIGHT, OPT_MODE, OPT_PRG, OPT_NOMCT, OPT_MAXRLVLS, OPT_SOP, OPT_EPH, OPT_LAZY, OPT_TERMALL, OPT_SEGSYM, OPT_VCAUSAL, OPT_RESET, OPT_PTERM, OPT_NUMGBITS, OPT_RATE, OPT_ILYRRATES, OPT_JP2OVERHEAD } optid_t; typedef enum { PO_L = 0, PO_R } poid_t; typedef enum { MODE_INT, MODE_REAL } modeid_t; /******************************************************************************\ * Data. \******************************************************************************/ static const jas_taginfo_t encopts[] = { {OPT_DEBUG, "debug"}, {OPT_IMGAREAOFFX, "imgareatlx"}, {OPT_IMGAREAOFFY, "imgareatly"}, {OPT_TILEGRDOFFX, "tilegrdtlx"}, {OPT_TILEGRDOFFY, "tilegrdtly"}, {OPT_TILEWIDTH, "tilewidth"}, {OPT_TILEHEIGHT, "tileheight"}, {OPT_PRCWIDTH, "prcwidth"}, {OPT_PRCHEIGHT, "prcheight"}, {OPT_CBLKWIDTH, "cblkwidth"}, {OPT_CBLKHEIGHT, "cblkheight"}, {OPT_MODE, "mode"}, {OPT_PRG, "prg"}, {OPT_NOMCT, "nomct"}, {OPT_MAXRLVLS, "numrlvls"}, {OPT_SOP, "sop"}, {OPT_EPH, "eph"}, {OPT_LAZY, "lazy"}, {OPT_TERMALL, "termall"}, {OPT_SEGSYM, "segsym"}, {OPT_VCAUSAL, "vcausal"}, {OPT_PTERM, "pterm"}, {OPT_RESET, "resetprob"}, {OPT_NUMGBITS, "numgbits"}, {OPT_RATE, "rate"}, {OPT_ILYRRATES, "ilyrrates"}, {OPT_JP2OVERHEAD, "_jp2overhead"}, {-1, 0} }; static const jas_taginfo_t prgordtab[] = { {JPC_COD_LRCPPRG, "lrcp"}, {JPC_COD_RLCPPRG, "rlcp"}, {JPC_COD_RPCLPRG, "rpcl"}, {JPC_COD_PCRLPRG, "pcrl"}, {JPC_COD_CPRLPRG, "cprl"}, {-1, 0} }; static const jas_taginfo_t modetab[] = { {MODE_INT, "int"}, {MODE_REAL, "real"}, {-1, 0} }; /******************************************************************************\ \******************************************************************************/ /** * @return UINT_FAST32_MAX on error */ static uint_fast32_t jpc_abstorelstepsize(jpc_fix_t absdelta, int scaleexpn) { int p; uint_fast32_t mant; uint_fast32_t expn; int n; if (absdelta < 0) { return UINT_FAST32_MAX; } p = jpc_fix_firstone(absdelta) - JPC_FIX_FRACBITS; n = 11 - jpc_fix_firstone(absdelta); mant = ((n < 0) ? (absdelta >> (-n)) : (absdelta << n)) & 0x7ff; expn = scaleexpn - p; if (scaleexpn < p) { return UINT_FAST32_MAX; } if (expn >= 0x1f) return UINT_FAST32_MAX; return JPC_QCX_EXPN(expn) | JPC_QCX_MANT(mant); } /******************************************************************************\ * The main encoder entry point. \******************************************************************************/ int jpc_encode(jas_image_t *image, jas_stream_t *out, const char *optstr) { jpc_enc_t *enc; jpc_enc_cp_t *cp; enc = 0; jpc_init(); if (!(cp = cp_create(optstr, image))) { jas_logerrorf("invalid JP encoder options\n"); goto error; } if (!(enc = jpc_enc_create(cp, out, image))) { jas_logerrorf("jpc_enc_create failed\n"); goto error; } cp = 0; /* Encode the main header. */ if (jpc_enc_encodemainhdr(enc)) { jas_logerrorf("cannot encode main header\n"); goto error; } /* Encode the main body. This constitutes most of the encoding work. */ if (jpc_enc_encodemainbody(enc)) { jas_logerrorf("cannot encode main body\n"); goto error; } /* Write EOC marker segment. */ if (!(enc->mrk = jpc_ms_create(JPC_MS_EOC))) { jas_logerrorf("cannot create EOC marker\n"); goto error; } if (jpc_putms(enc->out, enc->cstate, enc->mrk)) { jas_logerrorf("cannot write EOC marker\n"); goto error; } jpc_ms_destroy(enc->mrk); enc->mrk = 0; if (jas_stream_flush(enc->out)) { jas_logerrorf("stream flush failed\n"); goto error; } jpc_enc_destroy(enc); return 0; error: if (cp) { jpc_enc_cp_destroy(cp); } if (enc) { jpc_enc_destroy(enc); } jas_logerrorf("jpc_encode failed\n"); return -1; } /******************************************************************************\ * Option parsing code. \******************************************************************************/ static jpc_enc_cp_t *cp_create(const char *optstr, jas_image_t *image) { jpc_enc_cp_t *cp; jas_tvparser_t *tvp; int numilyrrates; double *ilyrrates; int i; int tagid; jpc_enc_tcp_t *tcp; jpc_enc_tccp_t *tccp; jpc_enc_ccp_t *ccp; uint_fast16_t rlvlno; uint_fast16_t prcwidthexpn; uint_fast16_t prcheightexpn; bool enablemct; uint_fast32_t jp2overhead; uint_fast16_t lyrno; uint_fast32_t hsteplcm; uint_fast32_t vsteplcm; bool mctvalid; tvp = 0; ilyrrates = 0; numilyrrates = 0; if (!(cp = jas_malloc(sizeof(jpc_enc_cp_t)))) { goto error; } prcwidthexpn = 15; prcheightexpn = 15; enablemct = true; jp2overhead = 0; cp->ccps = 0; cp->debug = 0; cp->imgareatlx = UINT_FAST32_MAX; cp->imgareatly = UINT_FAST32_MAX; cp->refgrdwidth = 0; cp->refgrdheight = 0; cp->tilegrdoffx = UINT_FAST32_MAX; cp->tilegrdoffy = UINT_FAST32_MAX; cp->tilewidth = 0; cp->tileheight = 0; cp->numcmpts = jas_image_numcmpts(image); cp->tcp.ilyrrates = NULL; hsteplcm = 1; vsteplcm = 1; for (unsigned cmptno = 0; cmptno < jas_image_numcmpts(image); ++cmptno) { if (jas_image_cmptbrx(image, cmptno) + jas_image_cmpthstep(image, cmptno) <= jas_image_brx(image) || jas_image_cmptbry(image, cmptno) + jas_image_cmptvstep(image, cmptno) <= jas_image_bry(image)) { jas_logerrorf("unsupported image type\n"); goto error; } /* Note: We ought to be calculating the LCMs here. Fix some day. */ hsteplcm *= jas_image_cmpthstep(image, cmptno); vsteplcm *= jas_image_cmptvstep(image, cmptno); } if (!(cp->ccps = jas_alloc2(cp->numcmpts, sizeof(jpc_enc_ccp_t)))) { goto error; } unsigned cmptno; for (cmptno = 0, ccp = cp->ccps; cmptno < cp->numcmpts; ++cmptno, ++ccp) { ccp->sampgrdstepx = jas_image_cmpthstep(image, cmptno); ccp->sampgrdstepy = jas_image_cmptvstep(image, cmptno); /* XXX - this isn't quite correct for more general image */ ccp->sampgrdsubstepx = 0; ccp->sampgrdsubstepx = 0; ccp->prec = jas_image_cmptprec(image, cmptno); ccp->sgnd = jas_image_cmptsgnd(image, cmptno); ccp->numstepsizes = 0; memset(ccp->stepsizes, 0, sizeof(ccp->stepsizes)); } cp->rawsize = jas_image_rawsize(image); if (cp->rawsize == 0) { /* prevent division by zero in cp_create() */ goto error; } cp->totalsize = UINT_FAST32_MAX; tcp = &cp->tcp; tcp->csty = 0; tcp->intmode = true; tcp->prg = JPC_COD_LRCPPRG; tcp->numlyrs = 1; tcp->ilyrrates = 0; tccp = &cp->tccp; tccp->csty = 0; tccp->maxrlvls = 6; tccp->cblkwidthexpn = 6; tccp->cblkheightexpn = 6; tccp->cblksty = 0; tccp->numgbits = 2; if (!(tvp = jas_tvparser_create(optstr ? optstr : ""))) { goto error; } while (!jas_tvparser_next(tvp)) { switch (jas_taginfo_nonull(jas_taginfos_lookup(encopts, jas_tvparser_gettag(tvp)))->id) { case OPT_DEBUG: cp->debug = atoi(jas_tvparser_getval(tvp)); break; case OPT_IMGAREAOFFX: cp->imgareatlx = atoi(jas_tvparser_getval(tvp)); break; case OPT_IMGAREAOFFY: cp->imgareatly = atoi(jas_tvparser_getval(tvp)); break; case OPT_TILEGRDOFFX: cp->tilegrdoffx = atoi(jas_tvparser_getval(tvp)); break; case OPT_TILEGRDOFFY: cp->tilegrdoffy = atoi(jas_tvparser_getval(tvp)); break; case OPT_TILEWIDTH: cp->tilewidth = atoi(jas_tvparser_getval(tvp)); break; case OPT_TILEHEIGHT: cp->tileheight = atoi(jas_tvparser_getval(tvp)); break; case OPT_PRCWIDTH: prcwidthexpn = jpc_floorlog2(atoi(jas_tvparser_getval(tvp))); break; case OPT_PRCHEIGHT: prcheightexpn = jpc_floorlog2(atoi(jas_tvparser_getval(tvp))); break; case OPT_CBLKWIDTH: tccp->cblkwidthexpn = jpc_floorlog2(atoi(jas_tvparser_getval(tvp))); break; case OPT_CBLKHEIGHT: tccp->cblkheightexpn = jpc_floorlog2(atoi(jas_tvparser_getval(tvp))); break; case OPT_MODE: if ((tagid = jas_taginfo_nonull(jas_taginfos_lookup(modetab, jas_tvparser_getval(tvp)))->id) < 0) { jas_logwarnf("ignoring invalid mode %s\n", jas_tvparser_getval(tvp)); } else { tcp->intmode = (tagid == MODE_INT); } break; case OPT_PRG: if ((tagid = jas_taginfo_nonull(jas_taginfos_lookup(prgordtab,