38#define _USE_MATH_DEFINES
94 ui32 bit_depth,
bool is_signed)
170 if (num_decompositions > 32)
172 "maximum number of decompositions cannot exceed 32");
181 if (width == 0 || width != (1u << log_width)
182 || height == 0 || height != (1u << log_height)
183 || log_width < 2 || log_height < 2
184 || log_width + log_height > 12)
185 OJPH_ERROR(0x00050011,
"incorrect code block dimensions");
193 if (num_levels == 0 || precinct_size == NULL)
200 size t = precinct_size[i < num_levels ? i : num_levels - 1];
204 if (t.
w == 0 || t.
h == 0)
205 OJPH_ERROR(0x00050021,
"precinct width or height cannot be 0");
206 if (t.
w != (1u<<PPx) || t.
h != (1u<<PPy))
208 "precinct width and height should be a power of 2");
209 if (PPx > 15 || PPy > 15)
210 OJPH_ERROR(0x00050023,
"precinct size is too large");
211 if (i > 0 && (PPx == 0 || PPy == 0))
212 OJPH_ERROR(0x00050024,
"precinct size is too small");
222 size_t len = strlen(name);
236 OJPH_ERROR(0x00050031,
"unknown progression order");
239 OJPH_ERROR(0x00050032,
"improper progression order");
445 bool& is_signed,
ui8& nl_type)
const
462 size_t t = strlen(str);
465 "COM marker string length cannot be larger than 65531");
476 "COM marker string length cannot be larger than 65531");
496 return (
ui16)((t << 8) | (t >> 8));
546 1.4021e+00f, 2.0304e+00f, 2.9012e+00f, 4.1153e+00f, 5.8245e+00f,
547 8.2388e+00f, 1.1652e+01f, 1.6479e+01f, 2.3304e+01f, 3.2957e+01f,
548 4.6609e+01f, 6.5915e+01f, 9.3217e+01f, 1.3183e+02f, 1.8643e+02f,
549 2.6366e+02f, 3.7287e+02f, 5.2732e+02f, 7.4574e+02f, 1.0546e+03f,
550 1.4915e+03f, 2.1093e+03f, 2.9830e+03f, 4.2185e+03f, 5.9659e+03f,
551 8.4371e+03f, 1.1932e+04f, 1.6874e+04f, 2.3864e+04f, 3.3748e+04f,
552 4.7727e+04f, 6.7496e+04f, 9.5454e+04f };
554 1.9669e+00f, 2.8839e+00f, 4.1475e+00f, 5.8946e+00f, 8.3472e+00f,
555 1.1809e+01f, 1.6701e+01f, 2.3620e+01f, 3.3403e+01f, 4.7240e+01f,
556 6.6807e+01f, 9.4479e+01f, 1.3361e+02f, 1.8896e+02f, 2.6723e+02f,
557 3.7792e+02f, 5.3446e+02f, 7.5583e+02f, 1.0689e+03f, 1.5117e+03f,
558 2.1378e+03f, 3.0233e+03f, 4.2756e+03f, 6.0467e+03f, 8.5513e+03f,
559 1.2093e+04f, 1.7103e+04f, 2.4187e+04f, 3.4205e+04f, 4.8373e+04f,
560 6.8410e+04f, 9.6747e+04f, 1.3682e+05f };
562 1.2247e+00f, 1.3229e+00f, 1.5411e+00f, 1.7139e+00f, 1.9605e+00f,
563 2.2044e+00f, 2.5047e+00f, 2.8277e+00f, 3.2049e+00f, 3.6238e+00f,
564 4.1033e+00f, 4.6423e+00f, 5.2548e+00f, 5.9462e+00f, 6.7299e+00f,
565 7.6159e+00f, 8.6193e+00f, 9.7544e+00f, 1.1039e+01f, 1.2493e+01f,
566 1.4139e+01f, 1.6001e+01f, 1.8108e+01f, 2.0493e+01f, 2.3192e+01f,
567 2.6246e+01f, 2.9702e+01f, 3.3614e+01f, 3.8041e+01f, 4.3051e+01f,
568 4.8721e+01f, 5.5138e+01f, 6.2399e+01f };
570 1.3975e+00f, 1.4389e+00f, 1.7287e+00f, 1.8880e+00f, 2.1841e+00f,
571 2.4392e+00f, 2.7830e+00f, 3.1341e+00f, 3.5576e+00f, 4.0188e+00f,
572 4.5532e+00f, 5.1494e+00f, 5.8301e+00f, 6.5963e+00f, 7.4663e+00f,
573 8.4489e+00f, 9.5623e+00f, 1.0821e+01f, 1.2247e+01f, 1.3860e+01f,
574 1.5685e+01f, 1.7751e+01f, 2.0089e+01f, 2.2735e+01f, 2.5729e+01f,
575 2.9117e+01f, 3.2952e+01f, 3.7292e+01f, 4.2203e+01f, 4.7761e+01f,
576 5.4051e+01f, 6.1170e+01f, 6.9226e+01f };
597 1.3328e+00f, 1.3067e+00f, 1.3028e+00f, 1.3001e+00f, 1.2993e+00f,
598 1.2992e+00f, 1.2992e+00f, 1.2992e+00f, 1.2992e+00f, 1.2992e+00f,
599 1.2992e+00f, 1.2992e+00f, 1.2992e+00f, 1.2992e+00f, 1.2992e+00f,
600 1.2992e+00f, 1.2992e+00f, 1.2992e+00f, 1.2992e+00f, 1.2992e+00f,
601 1.2992e+00f, 1.2992e+00f, 1.2992e+00f, 1.2992e+00f, 1.2992e+00f,
602 1.2992e+00f, 1.2992e+00f, 1.2992e+00f, 1.2992e+00f, 1.2992e+00f,
603 1.2992e+00f, 1.2992e+00f };
605 1.2757e+00f, 1.2352e+00f, 1.2312e+00f, 1.2285e+00f, 1.2280e+00f,
606 1.2278e+00f, 1.2278e+00f, 1.2278e+00f, 1.2278e+00f, 1.2278e+00f,
607 1.2278e+00f, 1.2278e+00f, 1.2278e+00f, 1.2278e+00f, 1.2278e+00f,
608 1.2278e+00f, 1.2278e+00f, 1.2278e+00f, 1.2278e+00f, 1.2278e+00f,
609 1.2278e+00f, 1.2278e+00f, 1.2278e+00f, 1.2278e+00f, 1.2278e+00f,
610 1.2278e+00f, 1.2278e+00f, 1.2278e+00f, 1.2278e+00f, 1.2278e+00f,
611 1.2278e+00f, 1.2278e+00f };
613 1.6250e+00f, 1.6875e+00f, 1.6963e+00f, 1.7067e+00f, 1.7116e+00f,
614 1.7129e+00f, 1.7141e+00f, 1.7145e+00f, 1.7151e+00f, 1.7152e+00f,
615 1.7155e+00f, 1.7155e+00f, 1.7156e+00f, 1.7156e+00f, 1.7156e+00f,
616 1.7156e+00f, 1.7156e+00f, 1.7156e+00f, 1.7156e+00f, 1.7156e+00f,
617 1.7156e+00f, 1.7156e+00f, 1.7156e+00f, 1.7156e+00f, 1.7156e+00f,
618 1.7156e+00f, 1.7156e+00f, 1.7156e+00f, 1.7156e+00f, 1.7156e+00f,
619 1.7156e+00f, 1.7156e+00f };
621 2.7500e+00f, 2.8047e+00f, 2.8198e+00f, 2.8410e+00f, 2.8558e+00f,
622 2.8601e+00f, 2.8628e+00f, 2.8656e+00f, 2.8662e+00f, 2.8667e+00f,
623 2.8669e+00f, 2.8670e+00f, 2.8671e+00f, 2.8671e+00f, 2.8671e+00f,
624 2.8671e+00f, 2.8671e+00f, 2.8671e+00f, 2.8671e+00f, 2.8671e+00f,
625 2.8671e+00f, 2.8671e+00f, 2.8671e+00f, 2.8671e+00f, 2.8671e+00f,
626 2.8671e+00f, 2.8671e+00f, 2.8671e+00f, 2.8671e+00f, 2.8671e+00f,
627 2.8671e+00f, 2.8671e+00f };
649 result &= file->
write(&buf, 2) == 2;
651 result &= file->
write(&buf, 2) == 2;
653 result &= file->
write(&buf, 2) == 2;
655 result &= file->
write(&buf, 4) == 4;
657 result &= file->
write(&buf, 4) == 4;
659 result &= file->
write(&buf, 4) == 4;
661 result &= file->
write(&buf, 4) == 4;
663 result &= file->
write(&buf, 4) == 4;
665 result &= file->
write(&buf, 4) == 4;
667 result &= file->
write(&buf, 4) == 4;
669 result &= file->
write(&buf, 4) == 4;
671 result &= file->
write(&buf, 2) == 2;
672 for (
int c = 0; c <
Csiz; ++c)
677 result &= file->
write(&buf, 3) == 3;
687 OJPH_ERROR(0x00050041,
"error reading SIZ marker");
689 int num_comps = (
Lsiz - 38) / 3;
690 if (
Lsiz != 38 + 3 * num_comps)
691 OJPH_ERROR(0x00050042,
"error in SIZ marker length");
693 OJPH_ERROR(0x00050043,
"error reading SIZ marker");
695 if ((
Rsiz & 0x4000) == 0)
697 "Rsiz bit 14 is not set (this is not a JPH file)");
698 if ((
Rsiz & 0x8000) != 0 && (
Rsiz & 0xD5F) != 0)
699 OJPH_WARN(0x00050001,
"Rsiz in SIZ has unimplemented fields");
701 OJPH_ERROR(0x00050045,
"error reading SIZ marker");
704 OJPH_ERROR(0x00050046,
"error reading SIZ marker");
707 OJPH_ERROR(0x00050047,
"error reading SIZ marker");
710 OJPH_ERROR(0x00050048,
"error reading SIZ marker");
713 OJPH_ERROR(0x00050049,
"error reading SIZ marker");
716 OJPH_ERROR(0x0005004A,
"error reading SIZ marker");
719 OJPH_ERROR(0x0005004B,
"error reading SIZ marker");
722 OJPH_ERROR(0x0005004C,
"error reading SIZ marker");
725 OJPH_ERROR(0x0005004D,
"error reading SIZ marker");
727 if (
Csiz != num_comps)
728 OJPH_ERROR(0x0005004E,
"Csiz does not match the SIZ marker size");
736 for (
int c = 0; c <
Csiz; ++c)
739 OJPH_ERROR(0x00050051,
"error reading SIZ marker");
741 OJPH_ERROR(0x00050052,
"error reading SIZ marker");
743 OJPH_ERROR(0x00050053,
"error reading SIZ marker");
798 result &= file->
write(&buf, 2) == 2;
800 result &= file->
write(&buf, 2) == 2;
802 result &= file->
write(&buf, 4) == 4;
805 result &= file->
write(&buf, 2) == 2;
814 OJPH_ERROR(0x00050061,
"error reading CAP marker");
817 OJPH_ERROR(0x00050062,
"error reading CAP marker");
820 if (
Pcap & 0xFFFDFFFF)
822 "error Pcap in CAP has options that are not supported");
823 if ((
Pcap & 0x00020000) == 0)
825 "error Pcap should have its 15th MSB set, Pcap^15. "
826 " This is not a JPH file");
827 for (
ui32 i = 0; i < count; ++i)
829 OJPH_ERROR(0x00050065,
"error reading CAP marker");
830 if (
Lcap != 6 + 2 * count)
831 OJPH_ERROR(0x00050066,
"error in CAP marker length");
867 result &= file->
write(&buf, 2) == 2;
869 result &= file->
write(&buf, 2) == 2;
871 result &= file->
write(&buf, 1) == 1;
873 result &= file->
write(&buf, 1) == 1;
875 result &= file->
write(&buf, 2) == 2;
877 result &= file->
write(&buf, 1) == 1;
882 result &= file->
write(&buf, 4) == 4;
884 result &= file->
write(&buf, 1) == 1;
889 result &= file->
write(&buf, 1) == 1;
916 Lcod = num_comps < 257 ? 9 : 10;
924 result &= file->
write(&buf, 2) == 2;
926 result &= file->
write(&buf, 2) == 2;
930 result &= file->
write(&buf, 1) == 1;
935 result &= file->
write(&buf, 2) == 2;
938 result &= file->
write(&buf, 1) == 1;
943 result &= file->
write(&buf, 4) == 4;
945 result &= file->
write(&buf, 1) == 1;
950 result &= file->
write(&buf, 1) == 1;
962 OJPH_ERROR(0x00050071,
"error reading COD segment");
965 OJPH_ERROR(0x00050072,
"error reading COD segment");
967 OJPH_ERROR(0x00050073,
"error reading COD segment");
969 {
OJPH_ERROR(0x00050074,
"error reading COD segment"); }
973 OJPH_ERROR(0x00050075,
"error reading COD segment");
975 OJPH_ERROR(0x00050076,
"error reading COD segment");
977 OJPH_ERROR(0x00050077,
"error reading COD segment");
979 OJPH_ERROR(0x00050078,
"error reading COD segment");
981 OJPH_ERROR(0x00050079,
"error reading COD segment");
983 OJPH_ERROR(0x0005007A,
"error reading COD segment");
987 OJPH_ERROR(0x0005007B,
"error reading COD segment");
989 OJPH_ERROR(0x0005007C,
"error in COD segment length");
1002 OJPH_ERROR(0x00050121,
"error reading COC segment");
1004 if (num_comps < 257) {
1006 if (file->
read(&t, 1) != 1)
1007 OJPH_ERROR(0x00050122,
"error reading COC segment");
1012 OJPH_ERROR(0x00050123,
"error reading COC segment");
1016 OJPH_ERROR(0x00050124,
"error reading COC segment");
1019 "Unsupported options in Scoc field of the COC segment");
1021 OJPH_ERROR(0x00050125,
"error reading COC segment");
1023 OJPH_ERROR(0x00050126,
"error reading COC segment");
1025 OJPH_ERROR(0x00050127,
"error reading COC segment");
1027 OJPH_ERROR(0x00050128,
"error reading COC segment");
1029 OJPH_ERROR(0x00050129,
"error reading COC segment");
1033 OJPH_ERROR(0x0005012A,
"error reading COC segment");
1035 t += num_comps < 257 ? 0 : 1;
1038 OJPH_ERROR(0x0005012B,
"error in COC segment length");
1046 if (this->atk == NULL)
1047 OJPH_ERROR(0x00050131,
"A COD segment employs the DWT kernel "
1048 "atk = %d, but a corresponding ATK segment cannot be found.",
1055 OJPH_ERROR(0x00050132,
"A COC segment employs the DWT kernel "
1056 "atk = %d, but a corresponding ATK segment cannot be found",
1090 while (p->
next != NULL)
1112 bool all_same =
true;
1113 bool other_comps_exist =
false;
1114 ui32 first_comp = 0xFFFF;
1116 ui32 num_decompositions = 0;
1118 bool is_signed =
false;
1121 for (
ui32 c = 0; c < num_comps; ++c)
1144 other_comps_exist =
true;
1149 ui32 qcd_num_decompositions;
1152 ui32 qcd_wavelet_kern;
1154 ui32 qcd_component = first_comp != 0xFFFF ? first_comp : 0;
1158 qcd_is_signed = siz.
is_signed(qcd_component);
1163 qcd_component < 3 ? employing_color_transform :
false);
1167 this->
base_delta = 1.0f / (float)(1 << qcd_bit_depth);
1178 for (
ui32 c = 0; c < num_comps; ++c)
1197 c < 3 ? employing_color_transform :
false);
1201 qp->
base_delta = 1.0f / (float)(1 << bit_depth);
1208 else if (other_comps_exist)
1211 for (
ui32 c = 0; c < num_comps; ++c)
1222 c < 3 ? employing_color_transform :
false);
1226 qp->
base_delta = 1.0f / (float)(1 << bit_depth);
1247 bool is_employing_color_transform)
1250 B += is_employing_color_transform ? 1 : 0;
1253 ui32 X = (
ui32) ceil(log(bibo_l * bibo_l) / M_LN2);
1256 for (
ui32 d = num_decomps; d > 0; --d)
1260 X = (
ui32) ceil(log(bibo_h * bibo_l) / M_LN2);
1262 max_B_plus_X =
ojph_max(max_B_plus_X, B + X);
1264 max_B_plus_X =
ojph_max(max_B_plus_X, B + X);
1265 X = (
ui32) ceil(log(bibo_h * bibo_h) / M_LN2);
1267 max_B_plus_X =
ojph_max(max_B_plus_X, B + X);
1270 if (max_B_plus_X > 38)
1271 OJPH_ERROR(0x00050151,
"The specified combination of bit_depth, "
1272 "colour transform, and type of wavelet transform requires more than "
1273 "38 bits; it requires %d bits. This is beyond what is allowed in "
1274 "the JPEG2000 image coding format.", max_B_plus_X);
1277 Sqcd = (
ui8)(guard_bits << 5);
1281 for (
ui32 d = num_decomps; d > 0; --d)
1296 Sqcd = (
ui8)((guard_bits<<5)|0x2);
1299 float delta_b =
base_delta / (gain_l * gain_l);
1300 int exp = 0, mantissa;
1301 while (delta_b < 1.0f)
1302 { exp++; delta_b *= 2.0f; }
1305 mantissa = (int)round(delta_b * (
float)(1<<11)) - (1<<11);
1306 mantissa = mantissa < (1<<11) ? mantissa : 0x7FF;
1307 SPqcd.u16[s++] = (
ui16)((exp << 11) | mantissa);
1308 for (
ui32 d = num_decomps; d > 0; --d)
1315 int exp = 0, mantissa;
1316 while (delta_b < 1.0f)
1317 { exp++; delta_b *= 2.0f; }
1318 mantissa = (int)round(delta_b * (
float)(1<<11)) - (1<<11);
1319 mantissa = mantissa < (1<<11) ? mantissa : 0x7FF;
1320 SPqcd.u16[s++] = (
ui16)((exp << 11) | mantissa);
1321 SPqcd.u16[s++] = (
ui16)((exp << 11) | mantissa);
1327 { exp++; delta_b *= 2.0f; }
1328 mantissa = (int)round(delta_b * (
float)(1<<11)) - (1<<11);
1329 mantissa = mantissa < (1<<11) ? mantissa : 0x7FF;
1330 SPqcd.u16[s++] = (
ui16)((exp << 11) | mantissa);
1346 int irrev = p->
Sqcd & 0x1F;
1353 else if (irrev == 2)
1356 ui32 nb = num_decomps - (i ? (i - 1) / 3 : 0);
1371 ui32 num_decompositions,
1374 float arr[] = { 1.0f, 2.0f, 2.0f, 4.0f };
1375 assert((
Sqcd & 0x1F) == 2);
1378 if (dfs != NULL && dfs->
exists())
1383 OJPH_INFO(0x00050101,
"Trying to access quantization step size for "
1384 "subband %d when the QCD/QCC marker segment specifies "
1385 "quantization step sizes for %d subbands only. To continue "
1386 "decoding, we are using the step size for subband %d, which can "
1387 "produce incorrect results",
1391 int eps =
SPqcd.u16[idx] >> 11;
1393 mantissa = (float)((
SPqcd.u16[idx] & 0x7FF) | 0x800) * arr[
subband];
1394 mantissa /= (float)(1 << 11);
1395 mantissa /= (float)(1u << eps);
1406 if (
main->is_employing_color_transform() &&
comp_idx < 3)
1408 for (
ui32 i = 0; i < 3; ++i) {
1421 return precision + 1 + 1;
1435 if (dfs != NULL && dfs->
exists())
1440 OJPH_INFO(0x00050111,
"Trying to access quantization step size for "
1441 "subband %d when the QCD/QCC marker segment specifies "
1442 "quantization step sizes for %d subbands only. To continue "
1443 "decoding, we are using the step size for subband %d, which can "
1444 "produce incorrect results",
1449 int irrev =
Sqcd & 0x1F;
1454 num_bits = num_bits == 0 ? 0 : num_bits - 1;
1456 else if (irrev == 1)
1458 else if (irrev == 2)
1459 num_bits = (
SPqcd.u16[idx] >> 11) - 1;
1469 int irrev =
Sqcd & 0x1F;
1475 num_bits =
ojph_max(num_bits, t == 0 ? 0 : t - 1);
1478 else if (irrev == 1)
1480 else if (irrev == 2)
1496 int irrev =
Sqcd & 0x1F;
1502 else if (irrev == 2)
1512 result &= file->
write(&buf, 2) == 2;
1514 result &= file->
write(&buf, 2) == 2;
1516 result &= file->
write(&buf, 1) == 1;
1522 result &= file->
write(&buf, 1) == 1;
1524 else if (irrev == 2)
1528 result &= file->
write(&buf, 2) == 2;
1554 int irrev =
Sqcd & 0x1F;
1557 Lqcd = (
ui16)(4 + (num_comps < 257 ? 0 : 1));
1560 else if (irrev == 2)
1570 result &= file->
write(&buf, 2) == 2;
1572 result &= file->
write(&buf, 2) == 2;
1573 if (num_comps < 257)
1576 result &= file->
write(&buf, 1) == 1;
1581 result &= file->
write(&buf, 2) == 2;
1584 result &= file->
write(&buf, 1) == 1;
1589 result &= file->
write(&buf, 1) == 1;
1591 else if (irrev == 2)
1595 result &= file->
write(&buf, 2) == 2;
1620 OJPH_ERROR(0x00050081,
"error reading QCD marker");
1623 OJPH_ERROR(0x00050082,
"error reading QCD marker");
1624 if ((
Sqcd & 0x1F) == 0)
1628 OJPH_ERROR(0x00050083,
"wrong Lqcd value in QCD marker");
1631 OJPH_ERROR(0x00050084,
"error reading QCD marker");
1633 else if ((
Sqcd & 0x1F) == 1)
1637 "Scalar derived quantization is not supported yet in QCD marker");
1639 OJPH_ERROR(0x00050085,
"wrong Lqcd value in QCD marker");
1641 else if ((
Sqcd & 0x1F) == 2)
1645 OJPH_ERROR(0x00050086,
"wrong Lqcd value in QCD marker");
1649 OJPH_ERROR(0x00050087,
"error reading QCD marker");
1654 OJPH_ERROR(0x00050088,
"wrong Sqcd value in QCD marker");
1661 OJPH_ERROR(0x000500A1,
"error reading QCC marker");
1663 if (num_comps < 257)
1666 if (file->
read(&v, 1) != 1)
1667 OJPH_ERROR(0x000500A2,
"error reading QCC marker");
1673 OJPH_ERROR(0x000500A3,
"error reading QCC marker");
1677 OJPH_ERROR(0x000500A4,
"error reading QCC marker");
1678 ui32 offset = num_comps < 257 ? 4 : 5;
1679 if ((
Sqcd & 0x1F) == 0)
1683 OJPH_ERROR(0x000500A5,
"wrong Lqcd value in QCC marker");
1686 OJPH_ERROR(0x000500A6,
"error reading QCC marker");
1688 else if ((
Sqcd & 0x1F) == 1)
1692 "Scalar derived quantization is not supported yet in QCC marker");
1694 OJPH_ERROR(0x000500A7,
"wrong Lqcc value in QCC marker");
1696 else if ((
Sqcd & 0x1F) == 2)
1700 OJPH_ERROR(0x000500A8,
"wrong Lqcc value in QCC marker");
1704 OJPH_ERROR(0x000500A9,
"error reading QCC marker");
1709 OJPH_ERROR(0x000500AA,
"wrong Sqcc value in QCC marker");
1740 while (p->
next != NULL)
1760 if (this->
enabled && this->
Tnlt == nonlinearity::OJPH_NLT_NO_NLT)
1764 this->
Tnlt == nonlinearity::OJPH_NLT_BINARY_COMPLEMENT_NLT)
1766 bool all_same =
true;
1774 bool is_signed =
false;
1775 for (
ui32 c = 0; c < num_comps; ++c)
1784 all_same = all_same && (is_signed == siz.
is_signed(c));
1800 if (all_same && bit_depth != 0)
1802 this->
BDnlt = (
ui8)(bit_depth - 1);
1808 for (
ui32 c = 0; c < num_comps; ++c)
1816 p->
Tnlt = nonlinearity::OJPH_NLT_BINARY_COMPLEMENT_NLT;
1826 for (
ui32 c = 0; c < num_comps; ++c)
1848 OJPH_ERROR(0x00050171,
"Nonliearities other than type 0 "
1849 "(No Nonlinearity) or type 3 (Binary Binary Complement to Sign "
1850 "Magnitude Conversion) are not supported yet");
1861 bool& is_signed,
ui8& nl_type)
const
1863 assert(
Cnlt == special_comp_num::ALL_COMPS);
1865 p = (p && p->
enabled) ? p :
this;
1868 bit_depth = (
ui8)((p->
BDnlt & 0x7F) + 1);
1869 bit_depth = bit_depth <= 38 ? bit_depth : 38;
1870 is_signed = (p->
BDnlt & 0x80) == 0x80;
1892 result &= file->
write(&buf, 2) == 2;
1894 result &= file->
write(&buf, 2) == 2;
1896 result &= file->
write(&buf, 2) == 2;
1898 result &= file->
write(&p->
Tnlt, 1) == 1;
1910 if (file->
read(buf, 6) != 6)
1911 OJPH_ERROR(0x00050141,
"error reading NLT marker segment");
1914 if (length != 6 || (buf[5] != 3 && buf[5] != 0))
1915 OJPH_ERROR(0x00050142,
"Unsupported NLT type %d\n", buf[5]);
1940 while (p && p->
Cnlt != comp_num)
1948 assert(comp_num != special_comp_num::ALL_COMPS);
1949 assert(
Cnlt == special_comp_num::ALL_COMPS);
1951 while (p->
next != NULL) {
1952 assert(p->
Cnlt != comp_num);
1967 while (p && p->
enabled ==
false)
1979 OJPH_INFO(0x00050161,
"The NLT marker segment for the "
1980 "non-existing component %d has been removed.", p->
Cnlt);
2001 this->
Psot = payload_len + 14;
2005 result &= file->
write(&buf, 2) == 2;
2007 result &= file->
write(&buf, 2) == 2;
2009 result &= file->
write(&buf, 2) == 2;
2011 result &= file->
write(&buf, 4) == 4;
2027 result &= file->
write(&buf, 2) == 2;
2029 result &= file->
write(&buf, 2) == 2;
2031 result &= file->
write(&buf, 2) == 2;
2033 result &= file->
write(&buf, 4) == 4;
2047 OJPH_INFO(0x00050091,
"error reading SOT marker");
2054 OJPH_INFO(0x00050092,
"error in SOT length");
2060 OJPH_INFO(0x00050093,
"error reading tile index");
2067 OJPH_INFO(0x00050094,
"tile index in SOT marker cannot be 0xFFFF");
2073 OJPH_INFO(0x00050095,
"error reading SOT marker");
2080 OJPH_INFO(0x00050096,
"error reading SOT marker");
2086 OJPH_INFO(0x00050097,
"error reading SOT marker");
2094 OJPH_ERROR(0x00050091,
"error reading SOT marker");
2097 OJPH_ERROR(0x00050092,
"error in SOT length");
2099 OJPH_ERROR(0x00050093,
"error reading SOT tile index");
2102 OJPH_ERROR(0x00050094,
"tile index in SOT marker cannot be 0xFFFF");
2104 OJPH_ERROR(0x00050095,
"error reading SOT marker");
2107 OJPH_ERROR(0x00050096,
"error reading SOT marker");
2109 OJPH_ERROR(0x00050097,
"error reading SOT marker");
2126 OJPH_ERROR(0x000500B1,
"Trying to allocate more than 65535 bytes for "
2127 "a TLM marker; this can be resolved by having more than "
2128 "one TLM marker, but the code does not support this. "
2129 "In any case, this limit means that we have 10922 "
2130 "tileparts or more, which is a huge number.");
2156 result &= file->
write(&buf, 2) == 2;
2158 result &= file->
write(&buf, 2) == 2;
2164 result &= file->
write(&buf, 2) == 2;
2166 result &= file->
write(&buf, 4) == 4;
2183 while (p && p->
Sdfs != index)
2192 ui32 d = decomp_level - 1;
2194 ui32 bits = d & 0x3;
2195 ui32 val = (
Ddfs[idx] >> (6 - 2 * bits)) & 0x3;
2206 ui32 ns[4] = { 0, 3, 1, 1 };
2228 ui32 decomp_level = 1;
2229 while (skipped_resolutions > 0)
2233 { factor.
x *= 2; factor.
y *= 2; }
2240 --skipped_resolutions;
2250 while (p->
next != NULL)
2254 return p->
read(file);
2258 OJPH_ERROR(0x000500D1,
"error reading DFS-Ldfs parameter");
2261 OJPH_ERROR(0x000500D2,
"error reading DFS-Sdfs parameter");
2264 OJPH_ERROR(0x000500D3,
"The DFS-Sdfs parameter is %d, which is "
2265 "larger than the permissible 15",
Sdfs);
2267 if (file->
read(&l_Ids, 1) != 1)
2268 OJPH_ERROR(0x000500D4,
"error reading DFS-Ids parameter");
2269 constexpr int max_Ddfs =
sizeof(
Ddfs) * 4;
2270 if (l_Ids > max_Ddfs)
2271 OJPH_INFO(0x000500D5,
"The DFS-Ids parameter is %d; while this is "
2272 "valid, the number is unnessarily large -- you do not need more "
2273 "than %d. Please contact me regarding this issue.",
2275 Ids = l_Ids < max_Ddfs ? l_Ids : max_Ddfs;
2276 for (
int i = 0; i <
Ids; i += 4)
2277 if (file->
read(&
Ddfs[i / 4], 1) != 1)
2278 OJPH_ERROR(0x000500D6,
"error reading DFS-Ddfs parameters");
2279 for (
int i =
Ids; i < l_Ids; i += 4)
2280 if (file->
read(&t, 1) != 1)
2281 OJPH_ERROR(0x000500D7,
"error reading DFS-Ddfs parameters");
2306 if (coeff_type == 0) {
2308 if (file->
read(&v, 1) != 1)
return false;
2311 else if (coeff_type == 1) {
2313 if (file->
read(&v, 2) != 2)
return false;
2316 else if (coeff_type == 2) {
2321 if (file->
read(&v.i, 4) != 4)
return false;
2325 else if (coeff_type == 3) {
2330 if (file->
read(&v.i, 8) != 8)
return false;
2334 else if (coeff_type == 4) {
2336 if (file->
read(&v, 8) != 8)
return false;
2337 if (file->
read(&v1, 8) != 8)
return false;
2347 si32 e = (
si32)((v >> 48) & 0x7FFF);
2353 s.i |= ((
ui32)(v >> 32) & 0x80000000);
2355 s.i |= (
ui32)((v >> 25) & 0x007FFFFF);
2366 if (coeff_type == 0) {
2368 if (file->
read(&v, 1) != 1)
return false;
2371 else if (coeff_type == 1) {
2373 if (file->
read(&v, 2) != 2)
return false;
2386 while (p->
next != NULL)
2391 return p->
read(file);
2395 OJPH_ERROR(0x000500E1,
"error reading ATK-Latk parameter");
2398 OJPH_ERROR(0x000500E2,
"error reading ATK-Satk parameter");
2401 OJPH_ERROR(0x000500E3,
"ATK-Satk parameter sets m_init to 1, "
2402 "requiring odd-indexed subsequence in first reconstruction step, "
2403 "which is not supported yet.");
2405 OJPH_ERROR(0x000500E4,
"ATK-Satk parameter specified ARB filter, "
2406 "which is not supported yet.");
2408 OJPH_ERROR(0x000500E5,
"ATK-Satk parameter does not make sense. "
2409 "It employs floats with reversible filtering.");
2411 OJPH_ERROR(0x000500E6,
"ATK-Satk parameter requires constant "
2412 "boundary extension, which is not supported yet.");
2415 OJPH_ERROR(0x000500E7,
"error reading ATK-Katk parameter");
2417 OJPH_ERROR(0x000500E8,
"error reading ATK-Natk parameter");
2427 for (
int s = 0; s <
Natk; ++s)
2430 OJPH_ERROR(0x000500E9,
"error reading ATK-Eatk parameter");
2432 OJPH_ERROR(0x000500EA,
"error reading ATK-Batk parameter");
2435 if (file->
read(&LCatk, 1) != 1)
2436 OJPH_ERROR(0x000500EB,
"error reading ATK-LCatk parameter");
2438 OJPH_ERROR(0x000500EC,
"Encountered a ATK-LCatk value of zero; "
2439 "something is wrong.");
2441 OJPH_ERROR(0x000500ED,
"ATK-LCatk value greater than 1; "
2442 "that is, a multitap filter is not supported");
2444 OJPH_ERROR(0x000500EE,
"Error reding ATK-Aatk parameter");
2449 for (
int s = 0; s <
Natk; ++s)
2452 if (file->
read(&LCatk, 1) != 1)
2453 OJPH_ERROR(0x000500EF,
"error reading ATK-LCatk parameter");
2455 OJPH_ERROR(0x000500F0,
"Encountered a ATK-LCatk value of zero; "
2456 "something is wrong.");
2458 OJPH_ERROR(0x000500F1,
"ATK-LCatk value greater than 1; "
2459 "that is, a multitap filter is not supported.");
2461 OJPH_ERROR(0x000500F2,
"Error reding ATK-Aatk parameter");
2472 Katk = (float)1.230174104914001;
2476 d[0].
irv.
Aatk = (float)0.443506852043971;
2477 d[1].
irv.
Aatk = (float)0.882911075530934;
2478 d[2].
irv.
Aatk = (float)-0.052980118572961;
2479 d[3].
irv.
Aatk = (float)-1.586134342059924;
int main(int argc, char *argv[])
virtual size_t read(void *ptr, size_t size)=0
static const float gain_5x3_l[34]
static float get_bibo_gain_l(ui32 num_decomp, bool reversible)
static const float gain_5x3_h[34]
static float get_bibo_gain_h(ui32 num_decomp, bool reversible)
static const float gain_9x7_h[34]
static const float gain_9x7_l[34]
static const float gain_5x3_l[34]
static const float gain_5x3_h[34]
static float get_gain_l(ui32 num_decomp, bool reversible)
static const float gain_9x7_l[34]
static float get_gain_h(ui32 num_decomp, bool reversible)
static const float gain_9x7_h[34]
virtual size_t write(const void *ptr, size_t size)=0
bool is_reversible() const
void set_precinct_size(int num_levels, size *precinct_size)
size get_block_dims() const
void set_reversible(bool reversible)
size get_precinct_size(ui32 level_num) const
ui32 get_num_decompositions() const
size get_log_block_dims() const
size get_log_precinct_size(ui32 level_num) const
void set_num_decomposition(ui32 num_decompositions)
bool get_block_vertical_causality() const
void set_block_dims(ui32 width, ui32 height)
size get_block_dims() const
int get_progression_order() const
bool is_using_color_transform() const
param_coc get_coc(ui32 component_idx)
void set_num_decomposition(ui32 num_decompositions)
ui32 get_num_decompositions() const
size get_log_block_dims() const
bool packets_may_use_sop() const
size get_precinct_size(ui32 level_num) const
const char * get_progression_order_as_string() const
void set_precinct_size(int num_levels, size *precinct_size)
bool packets_use_eph() const
bool is_reversible() const
void set_progression_order(const char *name)
bool get_block_vertical_causality() const
void set_block_dims(ui32 width, ui32 height)
size get_log_precinct_size(ui32 level_num) const
int get_num_layers() const
void set_color_transform(bool color_transform)
void set_reversible(bool reversible)
@ OJPH_NLT_BINARY_COMPLEMENT_NLT
bool get_nonlinear_transform(ui32 comp_num, ui8 &bit_depth, bool &is_signed, ui8 &nl_type) const
get the nonlinearity type associated with comp_num, which should be one from enum nonlinearity
void set_nonlinear_transform(ui32 comp_num, ui8 nl_type)
enables or disables type 3 nonlinearity for a component or the default setting
void set_irrev_quant(float delta)
Set the irreversible quantization base delta.
void set_tile_size(size s)
point get_image_extent() const
void set_component(ui32 comp_num, const point &downsampling, ui32 bit_depth, bool is_signed)
void set_num_components(ui32 num_comps)
ui32 get_bit_depth(ui32 comp_num) const
void set_tile_offset(point offset)
point get_image_offset() const
void set_image_offset(point offset)
size get_tile_size() const
ui32 get_recon_height(ui32 comp_num) const
point get_downsampling(ui32 comp_num) const
void set_image_extent(point extent)
point get_tile_offset() const
ui32 get_recon_width(ui32 comp_num) const
bool is_signed(ui32 comp_num) const
ui32 get_num_components() const
static ui16 swap_byte(ui16 t)
const char OJPH_PO_STRING_PCRL[]
static ui32 population_count(ui32 val)
const char OJPH_PO_STRING_RLCP[]
const char OJPH_PO_STRING_RPCL[]
const char OJPH_PO_STRING_CPRL[]
static ui32 count_leading_zeros(ui32 val)
const char OJPH_PO_STRING_LRCP[]
#define ojph_div_ceil(a, b)
#define OJPH_INFO(t,...)
MACROs to insert file and line number for info, warning, and error.
#define OJPH_ERROR(t,...)
const param_atk * get_atk(int index) const
bool is_reversible() const
bool read(infile_base *file)
int get_coeff_type() const
bool read_coefficient(infile_base *file, float &K)
bool is_whole_sample() const
bool is_using_ws_extension() const
void read(infile_base *file)
bool write(outfile_base *file)
bool write(outfile_base *file)
size get_block_dims() const
const param_cod * get_coc(ui32 comp_idx) const
bool internal_write_coc(outfile_base *file, ui32 num_comps)
size get_log_block_dims() const
ui16 get_dfs_index() const
bool write_coc(outfile_base *file, ui32 num_comps)
void set_reversible(bool reversible)
bool is_employing_color_transform() const
void employ_color_transform(ui8 val)
ui32 get_comp_idx() const
bool get_block_vertical_causality() const
void read(infile_base *file)
size get_log_precinct_size(ui32 res_num) const
void read_coc(infile_base *file, ui32 num_comps, param_cod *top_cod)
bool is_dfs_defined() const
bool packets_use_eph() const
ui8 get_num_decompositions() const
void update_atk(const param_atk *atk)
bool packets_may_use_sop() const
bool is_reversible() const
size get_precinct_size(ui32 res_num) const
param_cod(param_cod *top_cod=NULL, ui16 comp_idx=OJPH_COD_DEFAULT)
ui8 get_wavelet_kern() const
param_cod * add_coc_object(ui32 comp_idx)
bool read(infile_base *file)
dfs_dwt_type get_dwt_type(ui32 decomp_level) const
point get_res_downsamp(ui32 skipped_resolutions) const
ui32 get_subband_idx(ui32 num_decompositions, ui32 resolution, ui32 subband) const
const param_dfs * get_dfs(int index) const
bool write(outfile_base *file) const
param_nlt * add_object(ui32 comp_num)
void trim_non_existing_components(ui32 num_comps)
void read(infile_base *file)
ojph::param_nlt::nonlinearity nonlinearity
bool get_nonlinear_transform(ui32 comp_num, ui8 &bit_depth, bool &is_signed, ui8 &nl_type) const
const param_nlt * get_nlt_object(ui32 comp_num) const
void check_validity(param_siz &siz)
bool is_any_enabled() const
void set_nonlinear_transform(ui32 comp_num, ui8 nl_type)
ui8 encode_SPqcd(ui8 v) const
bool write_qcc(outfile_base *file, ui32 num_comps)
void set_rev_quant(ui32 num_decomps, ui32 bit_depth, bool is_employing_color_transform)
void set_irrev_quant(ui32 num_decomps)
ui32 get_largest_Kmax() const
ui32 get_num_guard_bits() const
float get_irrev_delta(const param_dfs *dfs, ui32 num_decompositions, ui32 resolution, ui32 subband) const
void set_delta(float delta)
void read_qcc(infile_base *file, ui32 num_comps)
void check_validity(const param_siz &siz, const param_cod &cod)
bool write(outfile_base *file)
ui32 propose_precision(const param_cod *cod) const
void read(infile_base *file)
param_qcd * add_qcc_object(ui32 comp_idx)
ui32 get_Kmax(const param_dfs *dfs, ui32 num_decompositions, ui32 resolution, ui32 subband) const
ui8 decode_SPqcd(ui8 v) const
param_qcd * get_qcc(ui32 comp_idx)
param_qcd(param_qcd *top_qcd=NULL, ui16 comp_idx=OJPH_QCD_DEFAULT)
union ojph::local::param_qcd::@6 SPqcd
void trim_non_existing_components(ui32 num_comps)
bool internal_write_qcc(outfile_base *file, ui32 num_comps)
ui32 get_bit_depth(ui32 comp_num) const
bool ws_kern_support_needed
ui32 get_recon_height(ui32 comp_num) const
bool is_signed(ui32 comp_num) const
bool write(outfile_base *file)
ui16 get_num_components() const
point get_recon_downsampling(ui32 comp_num) const
void set_Rsiz_flag(ui16 flag)
point get_recon_size(ui32 comp_num) const
void set_comp_info(ui32 comp_num, const point &downsampling, ui32 bit_depth, bool is_signed)
point get_downsampling(ui32 comp_num) const
void read(infile_base *file)
void set_num_components(ui32 num_comps)
ui32 get_recon_width(ui32 comp_num) const
bool read(infile_base *file, bool resilient)
bool write(outfile_base *file, ui32 payload_len)
void set_next_pair(ui16 Ttlm, ui32 Ptlm)
bool write(outfile_base *file)
void init(ui32 num_pairs, Ttlm_Ptlm_pair *store)