1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
|
// Copyright (c) the JPEG XL Project Authors. All rights reserved.
//
// Use of this source code is governed by a BSD-style
// license that can be found in the LICENSE file.
#include "lib/jxl/quantizer.h"
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include "lib/jxl/base/compiler_specific.h"
#include "lib/jxl/field_encodings.h"
#include "lib/jxl/fields.h"
#include "lib/jxl/image.h"
#include "lib/jxl/image_ops.h"
#include "lib/jxl/quant_weights.h"
namespace jxl {
static const int kDefaultQuant = 64;
constexpr int Quantizer::kQuantMax;
Quantizer::Quantizer(const DequantMatrices* dequant)
: Quantizer(dequant, kDefaultQuant, kGlobalScaleDenom / kDefaultQuant) {}
Quantizer::Quantizer(const DequantMatrices* dequant, int quant_dc,
int global_scale)
: global_scale_(global_scale), quant_dc_(quant_dc), dequant_(dequant) {
JXL_ASSERT(dequant_ != nullptr);
RecomputeFromGlobalScale();
inv_quant_dc_ = inv_global_scale_ / quant_dc_;
memcpy(zero_bias_, kZeroBiasDefault, sizeof(kZeroBiasDefault));
}
void Quantizer::ComputeGlobalScaleAndQuant(float quant_dc, float quant_median,
float quant_median_absd) {
// Target value for the median value in the quant field.
const float kQuantFieldTarget = 3.80987740592518214386f;
// We reduce the median of the quant field by the median absolute deviation:
// higher resolution on highly varying quant fields.
float scale = kGlobalScaleDenom * (quant_median - quant_median_absd) /
kQuantFieldTarget;
// Ensure that new_global_scale is positive and no more than 1<<15.
if (scale < 1) scale = 1;
if (scale > (1 << 15)) scale = 1 << 15;
int new_global_scale = static_cast<int>(scale);
// Ensure that quant_dc_ will always be at least
// kGlobalScaleDenom/kGlobalScaleNumerator.
const int scaled_quant_dc =
static_cast<int>(quant_dc * kGlobalScaleNumerator);
if (new_global_scale > scaled_quant_dc) {
new_global_scale = scaled_quant_dc;
if (new_global_scale <= 0) new_global_scale = 1;
}
global_scale_ = new_global_scale;
// Code below uses inv_global_scale_.
RecomputeFromGlobalScale();
float fval = quant_dc * inv_global_scale_ + 0.5f;
fval = std::min<float>(1 << 16, fval);
const int new_quant_dc = static_cast<int>(fval);
quant_dc_ = new_quant_dc;
// quant_dc_ was updated, recompute values.
RecomputeFromGlobalScale();
}
void Quantizer::SetQuantFieldRect(const ImageF& qf, const Rect& rect,
ImageI* JXL_RESTRICT raw_quant_field) {
for (size_t y = 0; y < rect.ysize(); ++y) {
const float* JXL_RESTRICT row_qf = rect.ConstRow(qf, y);
int32_t* JXL_RESTRICT row_qi = rect.Row(raw_quant_field, y);
for (size_t x = 0; x < rect.xsize(); ++x) {
int val = ClampVal(row_qf[x] * inv_global_scale_ + 0.5f);
row_qi[x] = val;
}
}
}
void Quantizer::SetQuantField(const float quant_dc, const ImageF& qf,
ImageI* JXL_RESTRICT raw_quant_field) {
JXL_CHECK(SameSize(*raw_quant_field, qf));
std::vector<float> data(qf.xsize() * qf.ysize());
for (size_t y = 0; y < qf.ysize(); ++y) {
const float* JXL_RESTRICT row_qf = qf.Row(y);
for (size_t x = 0; x < qf.xsize(); ++x) {
float quant = row_qf[x];
data[qf.xsize() * y + x] = quant;
}
}
std::nth_element(data.begin(), data.begin() + data.size() / 2, data.end());
const float quant_median = data[data.size() / 2];
std::vector<float> deviations(data.size());
for (size_t i = 0; i < data.size(); i++) {
deviations[i] = fabsf(data[i] - quant_median);
}
std::nth_element(deviations.begin(),
deviations.begin() + deviations.size() / 2,
deviations.end());
const float quant_median_absd = deviations[deviations.size() / 2];
ComputeGlobalScaleAndQuant(quant_dc, quant_median, quant_median_absd);
SetQuantFieldRect(qf, Rect(qf), raw_quant_field);
}
void Quantizer::SetQuant(float quant_dc, float quant_ac,
ImageI* JXL_RESTRICT raw_quant_field) {
ComputeGlobalScaleAndQuant(quant_dc, quant_ac, 0);
int val = ClampVal(quant_ac * inv_global_scale_ + 0.5f);
FillImage(val, raw_quant_field);
}
Status QuantizerParams::VisitFields(Visitor* JXL_RESTRICT visitor) {
JXL_QUIET_RETURN_IF_ERROR(visitor->U32(
BitsOffset(11, 1), BitsOffset(11, 2049), BitsOffset(12, 4097),
BitsOffset(16, 8193), 1, &global_scale));
JXL_QUIET_RETURN_IF_ERROR(visitor->U32(Val(16), BitsOffset(5, 1),
BitsOffset(8, 1), BitsOffset(16, 1), 1,
&quant_dc));
return true;
}
Status Quantizer::Encode(BitWriter* writer, size_t layer,
AuxOut* aux_out) const {
QuantizerParams params;
params.global_scale = global_scale_;
params.quant_dc = quant_dc_;
return Bundle::Write(params, writer, layer, aux_out);
}
Status Quantizer::Decode(BitReader* reader) {
QuantizerParams params;
JXL_RETURN_IF_ERROR(Bundle::Read(reader, ¶ms));
global_scale_ = static_cast<int>(params.global_scale);
quant_dc_ = static_cast<int>(params.quant_dc);
RecomputeFromGlobalScale();
return true;
}
void Quantizer::DumpQuantizationMap(const ImageI& raw_quant_field) const {
printf("Global scale: %d (%.7f)\nDC quant: %d\n", global_scale_,
global_scale_ * 1.0 / kGlobalScaleDenom, quant_dc_);
printf("AC quantization Map:\n");
for (size_t y = 0; y < raw_quant_field.ysize(); ++y) {
for (size_t x = 0; x < raw_quant_field.xsize(); ++x) {
printf(" %3d", raw_quant_field.Row(y)[x]);
}
printf("\n");
}
}
} // namespace jxl
|
