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
|
/*=========================================================================
Program: GDCM (Grassroots DICOM). A DICOM library
Copyright (c) 2006-2011 Mathieu Malaterre
All rights reserved.
See Copyright.txt or http://gdcm.sourceforge.net/Copyright.html for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notice for more information.
=========================================================================*/
#include "gdcmImageChangePhotometricInterpretation.h"
#include <math.h>
#include <string.h>
template <typename T>
double diff(T rgb1[3], T rgb2[2])
{
double sum = 0;
for(int i = 0; i < 3; ++i)
sum += fabs((double)(rgb1[i] - rgb2[i]));
return sum < 1e-3 ? 0 : sum;
//return sum;
}
int TestImageChangePhotometricInterpretation2(int argc, char *argv[])
{
//typedef float Type;
//typedef double Type;
//typedef int Type;
typedef unsigned char Type;
double sdiff = 0;
double max = 0;
int nerrors = 0;
int res = 0;
Type error[3] = {};
Type error2[3] = {};
Type yerror[3] = {};
for(int r = 0; r < 256; ++r)
for(int g = 0; g < 256; ++g)
for(int b = 0; b < 256; ++b)
{
Type rgb[3];
Type ybr[3] = {};
Type rgb2[3] = {};
rgb[0] = r;
rgb[1] = g;
rgb[1] = 128;
rgb[2] = b;
rgb[2] = 128;
// convert rgb 2 ybr:
//gdcm::ImageChangePhotometricInterpretation::RGB2YBR(ybr,rgb);
gdcm::ImageChangePhotometricInterpretation::YBR2RGB(ybr,rgb);
// convert back:
//gdcm::ImageChangePhotometricInterpretation::YBR2RGB(rgb2,ybr);
gdcm::ImageChangePhotometricInterpretation::RGB2YBR(rgb2,ybr);
if( memcmp(rgb,rgb2,3*sizeof(Type)) != 0 )
{
//std::cerr << "Problem with R,G,B=" << r << "," << g << "," << b <<
//" instead of " << (int)rgb2[0] << "," << (int)rgb2[1] << "," << (int)rgb2[2] << std::endl;
//std::cerr << "diff:" << diff(rgb,rgb2) << std::endl;
double d = diff(rgb,rgb2);
sdiff += d;
if( d > max )
{
error2[0] = rgb2[0];
error2[1] = rgb2[1];
error2[2] = rgb2[2];
error[0] = rgb[0];
error[1] = rgb[1];
error[2] = rgb[2];
yerror[0] = ybr[0];
yerror[1] = ybr[1];
yerror[2] = ybr[2];
}
max = std::max(d,max);
res = 1;
++nerrors;
}
}
std::cerr << "nerrors=" << nerrors << std::endl;
std::cerr << "sdiff=" << sdiff<< std::endl;
std::cerr << "max=" << max << std::endl;
std::cerr << "max error=" << (double)error[0] << "," << (double)error[1] << "," << (double)error[2] << std::endl;
std::cerr << "max error2=" << (double)error2[0] << "," << (double)error2[1] << "," << (double)error2[2] << std::endl;
std::cerr << "max yerror=" << (double)yerror[0] << "," << (double)yerror[1] << "," << (double)yerror[2] << std::endl;
return res;
}
|
