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
|
// 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 <stdint.h>
#include <stdio.h>
#include <string>
#include <vector>
#include "lib/extras/codec.h"
#include "lib/extras/dec/color_hints.h"
#include "lib/jxl/base/data_parallel.h"
#include "lib/jxl/base/file_io.h"
#include "lib/jxl/base/padded_bytes.h"
#include "lib/jxl/base/printf_macros.h"
#include "lib/jxl/base/status.h"
#include "lib/jxl/base/thread_pool_internal.h"
#include "lib/jxl/butteraugli/butteraugli.h"
#include "lib/jxl/codec_in_out.h"
#include "lib/jxl/color_encoding_internal.h"
#include "lib/jxl/color_management.h"
#include "lib/jxl/enc_butteraugli_comparator.h"
#include "lib/jxl/enc_butteraugli_pnorm.h"
#include "lib/jxl/enc_color_management.h"
#include "lib/jxl/image.h"
#include "lib/jxl/image_bundle.h"
#include "lib/jxl/image_ops.h"
namespace jxl {
namespace {
Status WriteImage(Image3F&& image, const std::string& filename) {
ThreadPoolInternal pool(4);
CodecInOut io;
io.metadata.m.SetUintSamples(8);
io.metadata.m.color_encoding = ColorEncoding::SRGB();
io.SetFromImage(std::move(image), io.metadata.m.color_encoding);
return EncodeToFile(io, filename, &pool);
}
Status RunButteraugli(const char* pathname1, const char* pathname2,
const std::string& distmap_filename,
const std::string& colorspace_hint, double p,
float intensity_target) {
extras::ColorHints color_hints;
if (!colorspace_hint.empty()) {
color_hints.Add("color_space", colorspace_hint);
}
CodecInOut io1;
ThreadPoolInternal pool(4);
if (!SetFromFile(pathname1, color_hints, &io1, &pool)) {
fprintf(stderr, "Failed to read image from %s\n", pathname1);
return false;
}
CodecInOut io2;
if (!SetFromFile(pathname2, color_hints, &io2, &pool)) {
fprintf(stderr, "Failed to read image from %s\n", pathname2);
return false;
}
if (io1.xsize() != io2.xsize()) {
fprintf(stderr, "Width mismatch: %" PRIuS " %" PRIuS "\n", io1.xsize(),
io2.xsize());
return false;
}
if (io1.ysize() != io2.ysize()) {
fprintf(stderr, "Height mismatch: %" PRIuS " %" PRIuS "\n", io1.ysize(),
io2.ysize());
return false;
}
ImageF distmap;
ButteraugliParams ba_params;
ba_params.hf_asymmetry = 0.8f;
ba_params.xmul = 1.0f;
ba_params.intensity_target = intensity_target;
const float distance = ButteraugliDistance(io1.Main(), io2.Main(), ba_params,
GetJxlCms(), &distmap, &pool);
printf("%.10f\n", distance);
double pnorm = ComputeDistanceP(distmap, ba_params, p);
printf("%g-norm: %f\n", p, pnorm);
if (!distmap_filename.empty()) {
float good = ButteraugliFuzzyInverse(1.5);
float bad = ButteraugliFuzzyInverse(0.5);
JXL_CHECK(
WriteImage(CreateHeatMapImage(distmap, good, bad), distmap_filename));
}
return true;
}
} // namespace
} // namespace jxl
int main(int argc, char** argv) {
if (argc < 3) {
fprintf(stderr,
"Usage: %s <reference> <distorted> [--distmap <distmap>] "
"[--intensity_target <intensity_target>]\n"
"[--colorspace <colorspace_hint>]\n"
"NOTE: images get converted to linear sRGB for butteraugli. Images"
" without attached profiles (such as ppm or pfm) are interpreted"
" as nonlinear sRGB. The hint format is RGB_D65_SRG_Rel_Lin for"
" linear sRGB. Intensity target is viewing conditions screen nits"
", defaults to 80.\n",
argv[0]);
return 1;
}
std::string distmap;
std::string colorspace;
double p = 3;
float intensity_target = 80.0; // sRGB intensity target.
for (int i = 3; i < argc; i++) {
if (std::string(argv[i]) == "--distmap" && i + 1 < argc) {
distmap = argv[++i];
} else if (std::string(argv[i]) == "--colorspace" && i + 1 < argc) {
colorspace = argv[++i];
} else if (std::string(argv[i]) == "--intensity_target" && i + 1 < argc) {
intensity_target = std::stof(std::string(argv[i + 1]));
} else if (std::string(argv[i]) == "--pnorm" && i + 1 < argc) {
char* end;
p = strtod(argv[++i], &end);
if (end == argv[i]) {
fprintf(stderr, "Failed to parse pnorm \"%s\".\n", argv[i]);
return 1;
}
} else {
fprintf(stderr, "Unrecognized flag \"%s\".\n", argv[i]);
return 1;
}
}
return jxl::RunButteraugli(argv[1], argv[2], distmap, colorspace, p,
intensity_target)
? 0
: 1;
}
|
