ebcot_lite_encode_passes.c

关于视频压缩的jpeg2000压缩算法,C编写

                        cp[3]              |= TC_SIG | V_PVE_SIG;
                      }
                  }
              }
            
            ctxt = cp[3];
            if (ctxt & IS_REFINED)
              cp[3] = ctxt & ~IS_REFINED;
            else if (!(ctxt & (SELF_SIG | OUT_OF_BOUNDS)))
              {
                val = dp[3];
                symbol = (val & mask);
                csp = csp_base + (ZC_OFFSET + zc_lut[ctxt & ZC_MASK]);
                dst_emit_symbol(areg,creg,ct,state,symbol,csp);
                if (symbol)
                  { /* New significant value; update contexts & code sign */
new_sig3:
                    sign = val & MIN_IFC_INT; /* Save sign bit. */
                    val &= MAX_IFC_INT;
                    val >>= mse_downshift; val <<= mse_upshift;
                    assert(val == (val & MSE_LUT_MASK));
                    cell->delta_mse += mse_lut[val]; /* David T Cvis mod */
                    symbol = ebcot_sc_lut[(ctxt>>SIGN_POS)&0x00FF];
                    csp = csp_base + SC_OFFSET + (symbol & 0x000F);
                    symbol &= MIN_IFC_INT;
                    symbol ^= sign;
                    dst_emit_symbol(areg,creg,ct,state,symbol,csp);
                    cp[3]                   = ctxt | SELF_SIG;
                    cp[2-4]                |= BR_SIG;
                    cp[2+4]                |= BL_SIG;
                    cp[stripe_gap-4]       |= TR_SIG;
                    cp[stripe_gap+4]       |= TL_SIG;
                    if (sign)
                      { /* Negative sample. */
                        cp[3+4]            |= CL_SIG | H_NVE_SIG;
                        cp[3-4]            |= CR_SIG | H_NVE_SIG;
                        cp[2]              |= BC_SIG | V_NVE_SIG;
                        cp[stripe_gap]     |= TC_SIG | V_NVE_SIG;
                      }
                    else
                      { /* Positive sample. */
                        cp[3+4]            |= CL_SIG | H_PVE_SIG;
                        cp[3-4]            |= CR_SIG | H_PVE_SIG;
                        cp[2]              |= BC_SIG | V_PVE_SIG;
                        cp[stripe_gap]     |= TC_SIG | V_PVE_SIG;
                      }
                  }
              }
          }
      }
  state->mqe.areg = areg; state->mqe.creg = creg; state->mqe.ct = ct;
}

/*****************************************************************************/
/* STATIC                  raw_zero_refinement_pass_func                     */
/*****************************************************************************/

static void /* David T Cvis mod */
  raw_zero_refinement_pass_func(block_master_ptr master)

 /* This is a version of `zero_refinement_pass_func' which outputs raw
    binary digits, bypassing the arithmetic coder altogether.  It may be
    used only with the restart option. */

{
  register std_short *cp;
  register int i;
  register ifc_int *dp;
  register std_short ctxt;
  register dst_arith_state_ptr state;
  ifc_int val, symbol, lsb, shift;
  std_int word;
  std_short current_lsb;
  ifc_int mse_upshift, mse_downshift;
  distortion_cell_ptr cell; /* David T Cvis mod */
  int cell_gap; /* David T Cvis mod */
  std_short *mse_lut, non_causal;
  int eighth_cols, stripe_gap, r, c;

  shift = master->bit_idx - 1;
  lsb = 1 << shift;
  mse_lut = master->initial_mse_lut;
  mse_downshift = shift - (MSE_LUT_BITS-1);
  mse_upshift = (mse_downshift<0)?-mse_downshift:0;
  mse_downshift += mse_upshift;
  state = &(master->coder_state);
  current_lsb = state->current_lsb; word = state->word;
  non_causal = 1 - master->causal;
  eighth_cols = (master->width+7) >> 3;
  stripe_gap = master->interleaved_row_gap;
  cell_gap = master->d_cell_row_gap; /* David T Cvis mod */
  cp = master->interleaved_context_buffer;
  dp = master->interleaved_sample_buffer;
  cell = master->d_cells; /* David T Cvis mod */
  for (r=master->stripes; r > 0; r--,
       dp+=stripe_gap-(eighth_cols<<5), cp+=stripe_gap-(eighth_cols<<5),
       cell += ((r&1)?0:cell_gap)-eighth_cols) /* David T Cvis mod */
    for (c=eighth_cols; c > 0; c--, cell++) /* David T Cvis mod */
      {
        for (i=8; i > 0; i--, dp+=4, cp+=4)
          {
            if (((((std_int *) cp)[0] |
                  ((std_int *) cp)[1]) & 0x00FF00FF) == 0)
              continue;
            if ((((std_int *) cp)[0] & 0x00FF00FF) == 0)
              goto skip2;

            ctxt = cp[0];
            if ((ctxt & 0x00FF) &&
                !(ctxt & (IS_REFINED|SELF_SIG|OUT_OF_BOUNDS)))
              {
                ctxt |= IS_REFINED;
                val = dp[0];
                symbol = (val >> shift) & 1;
                dst_emit_raw_bit(current_lsb,word,symbol,state);
                if (symbol)
                  { /* New significant value; update contexts & code sign */
                    ifc_int sign;

                    sign = (val >> (IMPLEMENTATION_PRECISION-1)) & 1;
                    val &= MAX_IFC_INT;
                    val >>= mse_downshift; val <<= mse_upshift;
                    assert(val == (val & MSE_LUT_MASK));
                    cell->delta_mse += mse_lut[val]; /* David T Cvis mod */
                    dst_emit_raw_bit(current_lsb,word,sign,state);
                    ctxt |= SELF_SIG;
                    cp[(-stripe_gap+3)-4]  |= (non_causal << BR_POS);
                    cp[(-stripe_gap+3)+4]  |= (non_causal << BL_POS);
                    cp[1-4]                |= TR_SIG;
                    cp[1+4]                |= TL_SIG;
                    if (sign)
                      { /* Negative sample. */
                        cp[4]              |= CL_SIG | H_NVE_SIG;
                        cp[-4]             |= CR_SIG | H_NVE_SIG;
                        cp[(-stripe_gap+3)]|= (non_causal << BC_POS) |
                                              (non_causal << V_NVE_POS);
                        cp[1]              |= TC_SIG | V_NVE_SIG;
                      }
                    else
                      { /* Positive sample. */
                        cp[4]              |= CL_SIG | H_PVE_SIG;
                        cp[-4]             |= CR_SIG | H_PVE_SIG;
                        cp[(-stripe_gap+3)]|= (non_causal << BC_POS) |
                                              (non_causal << V_PVE_POS);
                        cp[1]              |= TC_SIG | V_PVE_SIG;
                      }
                  }
                cp[0] = ctxt;
              }
            
            ctxt = cp[1];
            if ((ctxt & 0x00FF) &&
                !(ctxt & (IS_REFINED|SELF_SIG|OUT_OF_BOUNDS)))
              {
                ctxt |= IS_REFINED;
                val = dp[1];
                symbol = (val >> shift) & 1;
                dst_emit_raw_bit(current_lsb,word,symbol,state);
                if (symbol)
                  { /* New significant value; update contexts & code sign */
                    ifc_int sign;

                    sign = (val >> (IMPLEMENTATION_PRECISION-1)) & 1;
                    val &= MAX_IFC_INT;
                    val >>= mse_downshift; val <<= mse_upshift;
                    assert(val == (val & MSE_LUT_MASK));
                    cell->delta_mse += mse_lut[val]; /* David T Cvis mod */
                    dst_emit_raw_bit(current_lsb,word,sign,state);
                    ctxt |= SELF_SIG;
                    cp[0-4]                |= BR_SIG;
                    cp[0+4]                |= BL_SIG;
                    cp[2-4]                |= TR_SIG;
                    cp[2+4]                |= TL_SIG;
                    if (sign)
                      { /* Negative sample. */
                        cp[1+4]            |= CL_SIG | H_NVE_SIG;
                        cp[1-4]            |= CR_SIG | H_NVE_SIG;
                        cp[0]              |= BC_SIG | V_NVE_SIG;
                        cp[2]              |= TC_SIG | V_NVE_SIG;
                      }
                    else
                      { /* Positive sample. */
                        cp[1+4]            |= CL_SIG | H_PVE_SIG;
                        cp[1-4]            |= CR_SIG | H_PVE_SIG;
                        cp[0]              |= BC_SIG | V_PVE_SIG;
                        cp[2]              |= TC_SIG | V_PVE_SIG;
                      }
                  }
                cp[1] = ctxt;
              }
            if ((((std_int *) cp)[1]  & 0x00FF00FF) == 0)
              continue;
skip2:
            ctxt = cp[2];
            if ((ctxt & 0x00FF) &&
                !(ctxt & (IS_REFINED|SELF_SIG|OUT_OF_BOUNDS)))
              {
                ctxt |= IS_REFINED;
                val = dp[2];
                symbol = (val >> shift) & 1;
                dst_emit_raw_bit(current_lsb,word,symbol,state);
                if (symbol)
                  { /* New significant value; update contexts & code sign */
                    ifc_int sign;

                    sign = (val >> (IMPLEMENTATION_PRECISION-1)) & 1;
                    val &= MAX_IFC_INT;
                    val >>= mse_downshift; val <<= mse_upshift;
                    assert(val == (val & MSE_LUT_MASK));
                    cell->delta_mse += mse_lut[val]; /* David T Cvis mod */
                    dst_emit_raw_bit(current_lsb,word,sign,state);
                    ctxt |= SELF_SIG;
                    cp[1-4]                |= BR_SIG;
                    cp[1+4]                |= BL_SIG;
                    cp[3-4]                |= TR_SIG;
                    cp[3+4]                |= TL_SIG;
                    if (sign)
                      { /* Negative sample. */
                        cp[2+4]            |= CL_SIG | H_NVE_SIG;
                        cp[2-4]            |= CR_SIG | H_NVE_SIG;
                        cp[1]              |= BC_SIG | V_NVE_SIG;
                        cp[3]              |= TC_SIG | V_NVE_SIG;
                      }
                    else
                      { /* Positive sample. */
                        cp[2+4]            |= CL_SIG | H_PVE_SIG;
                        cp[2-4]            |= CR_SIG | H_PVE_SIG;
                        cp[1]              |= BC_SIG | V_PVE_SIG;
                        cp[3]              |= TC_SIG | V_PVE_SIG;
                      }
                  }
                cp[2] = ctxt;
              }

            ctxt = cp[3];
            if ((ctxt & 0x00FF) &&
                !(ctxt & (IS_REFINED|SELF_SIG|OUT_OF_BOUNDS)))
              {
                ctxt |= IS_REFINED;
                val = dp[3];
                symbol = (val >> shift) & 1;
                dst_emit_raw_bit(current_lsb,word,symbol,state);
                if (symbol)
                  { /* New significant value; update contexts & code sign */
                    ifc_int sign;

                    sign = (val >> (IMPLEMENTATION_PRECISION-1)) & 1;
                    val &= MAX_IFC_INT;
                    val >>= mse_downshift; val <<= mse_upshift;
                    assert(val == (val & MSE_LUT_MASK));
                    cell->delta_mse += mse_lut[val]; /* David T Cvis mod */
                    dst_emit_raw_bit(current_lsb,word,sign,state);
                    ctxt |= SELF_SIG;
                    cp[2-4]                |= BR_SIG;
                    cp[2+4]                |= BL_SIG;
                    cp[stripe_gap-4]       |= TR_SIG;
                    cp[stripe_gap+4]       |= TL_SIG;
                    if (sign)
                      { /* Negative sample. */
                        cp[3+4]            |= CL_SIG | H_NVE_SIG;
                        cp[3-4]            |= CR_SIG | H_NVE_SIG;
                        cp[2]              |= BC_SIG | V_NVE_SIG;
                        cp[stripe_gap]     |= TC_SIG | V_NVE_SIG;
                      }
                    else
                      { /* Positive sample. */
                        cp[3+4]            |= CL_SIG | H_PVE_SIG;
                        cp[3-4]            |= CR_SIG | H_PVE_SIG;
                        cp[2]              |= BC_SIG | V_PVE_SIG;
                        cp[stripe_gap]     |= TC_SIG | V_PVE_SIG;
                      }
                  }
                cp[3] = ctxt;
              }
          }
      }
  state->current_lsb = current_lsb; state->word = word;
}