| OLD | NEW |
| (Empty) |
| 1 // qcms | |
| 2 // Copyright (C) 2009 Mozilla Foundation | |
| 3 // Copyright (C) 1998-2007 Marti Maria | |
| 4 // | |
| 5 // Permission is hereby granted, free of charge, to any person obtaining | |
| 6 // a copy of this software and associated documentation files (the "Software"), | |
| 7 // to deal in the Software without restriction, including without limitation | |
| 8 // the rights to use, copy, modify, merge, publish, distribute, sublicense, | |
| 9 // and/or sell copies of the Software, and to permit persons to whom the Softwar
e | |
| 10 // is furnished to do so, subject to the following conditions: | |
| 11 // | |
| 12 // The above copyright notice and this permission notice shall be included in | |
| 13 // all copies or substantial portions of the Software. | |
| 14 // | |
| 15 // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, | |
| 16 // EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO | |
| 17 // THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND | |
| 18 // NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE | |
| 19 // LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION | |
| 20 // OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION | |
| 21 // WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE. | |
| 22 | |
| 23 #include <math.h> | |
| 24 #include <assert.h> | |
| 25 #include <stdlib.h> | |
| 26 #include <string.h> | |
| 27 #include "qcmsint.h" | |
| 28 | |
| 29 typedef uint32_t be32; | |
| 30 typedef uint16_t be16; | |
| 31 | |
| 32 #if 0 | |
| 33 not used yet | |
| 34 /* __builtin_bswap isn't available in older gccs | |
| 35 * so open code it for now */ | |
| 36 static be32 cpu_to_be32(int32_t v) | |
| 37 { | |
| 38 #ifdef IS_LITTLE_ENDIAN | |
| 39 return ((v & 0xff) << 24) | ((v & 0xff00) << 8) | ((v & 0xff0000) >> 8)
| ((v & 0xff000000) >> 24); | |
| 40 //return __builtin_bswap32(v); | |
| 41 return v; | |
| 42 #endif | |
| 43 } | |
| 44 #endif | |
| 45 | |
| 46 static uint32_t be32_to_cpu(be32 v) | |
| 47 { | |
| 48 #ifdef IS_LITTLE_ENDIAN | |
| 49 return ((v & 0xff) << 24) | ((v & 0xff00) << 8) | ((v & 0xff0000) >> 8)
| ((v & 0xff000000) >> 24); | |
| 50 //return __builtin_bswap32(v); | |
| 51 #else | |
| 52 return v; | |
| 53 #endif | |
| 54 } | |
| 55 | |
| 56 static uint16_t be16_to_cpu(be16 v) | |
| 57 { | |
| 58 #ifdef IS_LITTLE_ENDIAN | |
| 59 return ((v & 0xff) << 8) | ((v & 0xff00) >> 8); | |
| 60 #else | |
| 61 return v; | |
| 62 #endif | |
| 63 } | |
| 64 | |
| 65 /* a wrapper around the memory that we are going to parse | |
| 66 * into a qcms_profile */ | |
| 67 struct mem_source | |
| 68 { | |
| 69 const unsigned char *buf; | |
| 70 size_t size; | |
| 71 qcms_bool valid; | |
| 72 const char *invalid_reason; | |
| 73 }; | |
| 74 | |
| 75 static void invalid_source(struct mem_source *mem, const char *reason) | |
| 76 { | |
| 77 mem->valid = false; | |
| 78 mem->invalid_reason = reason; | |
| 79 } | |
| 80 | |
| 81 static uint32_t read_u32(struct mem_source *mem, size_t offset) | |
| 82 { | |
| 83 /* Subtract from mem->size instead of the more intuitive adding to offse
t. | |
| 84 * This avoids overflowing offset. The subtraction is safe because | |
| 85 * mem->size is guaranteed to be > 4 */ | |
| 86 if (offset > mem->size - 4) { | |
| 87 invalid_source(mem, "Invalid offset"); | |
| 88 return 0; | |
| 89 } else { | |
| 90 be32 k; | |
| 91 memcpy(&k, mem->buf + offset, sizeof(k)); | |
| 92 return be32_to_cpu(k); | |
| 93 } | |
| 94 } | |
| 95 | |
| 96 static uint16_t read_u16(struct mem_source *mem, size_t offset) | |
| 97 { | |
| 98 if (offset > mem->size - 2) { | |
| 99 invalid_source(mem, "Invalid offset"); | |
| 100 return 0; | |
| 101 } else { | |
| 102 be16 k; | |
| 103 memcpy(&k, mem->buf + offset, sizeof(k)); | |
| 104 return be16_to_cpu(k); | |
| 105 } | |
| 106 } | |
| 107 | |
| 108 static uint8_t read_u8(struct mem_source *mem, size_t offset) | |
| 109 { | |
| 110 if (offset > mem->size - 1) { | |
| 111 invalid_source(mem, "Invalid offset"); | |
| 112 return 0; | |
| 113 } else { | |
| 114 return *(uint8_t*)(mem->buf + offset); | |
| 115 } | |
| 116 } | |
| 117 | |
| 118 static s15Fixed16Number read_s15Fixed16Number(struct mem_source *mem, size_t off
set) | |
| 119 { | |
| 120 return read_u32(mem, offset); | |
| 121 } | |
| 122 | |
| 123 #if 0 | |
| 124 static uInt16Number read_uInt16Number(struct mem_source *mem, size_t offset) | |
| 125 { | |
| 126 return read_u16(mem, offset); | |
| 127 } | |
| 128 #endif | |
| 129 | |
| 130 #define BAD_VALUE_PROFILE NULL | |
| 131 #define INVALID_PROFILE NULL | |
| 132 #define NO_MEM_PROFILE NULL | |
| 133 | |
| 134 /* An arbitrary 4MB limit on profile size */ | |
| 135 #define MAX_PROFILE_SIZE 1024*1024*4 | |
| 136 #define MAX_TAG_COUNT 1024 | |
| 137 | |
| 138 static void check_CMM_type_signature(struct mem_source *src) | |
| 139 { | |
| 140 //uint32_t CMM_type_signature = read_u32(src, 4); | |
| 141 //TODO: do the check? | |
| 142 | |
| 143 } | |
| 144 | |
| 145 static void check_profile_version(struct mem_source *src) | |
| 146 { | |
| 147 uint8_t major_revision = read_u8(src, 8 + 0); | |
| 148 uint8_t minor_revision = read_u8(src, 8 + 1); | |
| 149 uint8_t reserved1 = read_u8(src, 8 + 2); | |
| 150 uint8_t reserved2 = read_u8(src, 8 + 3); | |
| 151 if (major_revision != 0x4) { | |
| 152 if (major_revision > 0x2) | |
| 153 invalid_source(src, "Unsupported major revision"); | |
| 154 if (minor_revision > 0x40) | |
| 155 invalid_source(src, "Unsupported minor revision"); | |
| 156 } | |
| 157 if (reserved1 != 0 || reserved2 != 0) | |
| 158 invalid_source(src, "Invalid reserved bytes"); | |
| 159 } | |
| 160 | |
| 161 #define INPUT_DEVICE_PROFILE 0x73636e72 // 'scnr' | |
| 162 #define DISPLAY_DEVICE_PROFILE 0x6d6e7472 // 'mntr' | |
| 163 #define OUTPUT_DEVICE_PROFILE 0x70727472 // 'prtr' | |
| 164 #define DEVICE_LINK_PROFILE 0x6c696e6b // 'link' | |
| 165 #define COLOR_SPACE_PROFILE 0x73706163 // 'spac' | |
| 166 #define ABSTRACT_PROFILE 0x61627374 // 'abst' | |
| 167 #define NAMED_COLOR_PROFILE 0x6e6d636c // 'nmcl' | |
| 168 | |
| 169 static void read_class_signature(qcms_profile *profile, struct mem_source *mem) | |
| 170 { | |
| 171 profile->class = read_u32(mem, 12); | |
| 172 switch (profile->class) { | |
| 173 case DISPLAY_DEVICE_PROFILE: | |
| 174 case INPUT_DEVICE_PROFILE: | |
| 175 break; | |
| 176 case OUTPUT_DEVICE_PROFILE: | |
| 177 default: | |
| 178 invalid_source(mem, "Invalid Profile/Device Class signa
ture"); | |
| 179 } | |
| 180 } | |
| 181 | |
| 182 static void read_color_space(qcms_profile *profile, struct mem_source *mem) | |
| 183 { | |
| 184 profile->color_space = read_u32(mem, 16); | |
| 185 switch (profile->color_space) { | |
| 186 case RGB_SIGNATURE: | |
| 187 case GRAY_SIGNATURE: | |
| 188 break; | |
| 189 default: | |
| 190 invalid_source(mem, "Unsupported colorspace"); | |
| 191 } | |
| 192 } | |
| 193 | |
| 194 struct tag | |
| 195 { | |
| 196 uint32_t signature; | |
| 197 uint32_t offset; | |
| 198 uint32_t size; | |
| 199 }; | |
| 200 | |
| 201 struct tag_index { | |
| 202 uint32_t count; | |
| 203 struct tag *tags; | |
| 204 }; | |
| 205 | |
| 206 static struct tag_index read_tag_table(qcms_profile *profile, struct mem_source
*mem) | |
| 207 { | |
| 208 struct tag_index index = {0, NULL}; | |
| 209 unsigned int i; | |
| 210 | |
| 211 index.count = read_u32(mem, 128); | |
| 212 if (index.count > MAX_TAG_COUNT) { | |
| 213 invalid_source(mem, "max number of tags exceeded"); | |
| 214 return index; | |
| 215 } | |
| 216 | |
| 217 index.tags = malloc(sizeof(struct tag)*index.count); | |
| 218 if (index.tags) { | |
| 219 for (i = 0; i < index.count; i++) { | |
| 220 index.tags[i].signature = read_u32(mem, 128 + 4 + 4*i*3)
; | |
| 221 index.tags[i].offset = read_u32(mem, 128 + 4 + 4*i*3
+ 4); | |
| 222 index.tags[i].size = read_u32(mem, 128 + 4 + 4*i*3
+ 8); | |
| 223 } | |
| 224 } | |
| 225 | |
| 226 return index; | |
| 227 } | |
| 228 | |
| 229 // Checks a profile for obvious inconsistencies and returns | |
| 230 // true if the profile looks bogus and should probably be | |
| 231 // ignored. | |
| 232 qcms_bool qcms_profile_is_bogus(qcms_profile *profile) | |
| 233 { | |
| 234 float sum[3], target[3], tolerance[3]; | |
| 235 float rX, rY, rZ, gX, gY, gZ, bX, bY, bZ; | |
| 236 bool negative; | |
| 237 unsigned i; | |
| 238 | |
| 239 // We currently only check the bogosity of RGB profiles | |
| 240 if (profile->color_space != RGB_SIGNATURE) | |
| 241 return false; | |
| 242 | |
| 243 rX = s15Fixed16Number_to_float(profile->redColorant.X); | |
| 244 rY = s15Fixed16Number_to_float(profile->redColorant.Y); | |
| 245 rZ = s15Fixed16Number_to_float(profile->redColorant.Z); | |
| 246 | |
| 247 gX = s15Fixed16Number_to_float(profile->greenColorant.X); | |
| 248 gY = s15Fixed16Number_to_float(profile->greenColorant.Y); | |
| 249 gZ = s15Fixed16Number_to_float(profile->greenColorant.Z); | |
| 250 | |
| 251 bX = s15Fixed16Number_to_float(profile->blueColorant.X); | |
| 252 bY = s15Fixed16Number_to_float(profile->blueColorant.Y); | |
| 253 bZ = s15Fixed16Number_to_float(profile->blueColorant.Z); | |
| 254 | |
| 255 // Check if any of the XYZ values are negative (see mozilla bug 498245) | |
| 256 // CIEXYZ tristimulus values cannot be negative according to the spec. | |
| 257 negative = | |
| 258 (rX < 0) || (rY < 0) || (rZ < 0) || | |
| 259 (gX < 0) || (gY < 0) || (gZ < 0) || | |
| 260 (bX < 0) || (bY < 0) || (bZ < 0); | |
| 261 | |
| 262 if (negative) | |
| 263 return true; | |
| 264 | |
| 265 | |
| 266 // Sum the values; they should add up to something close to white | |
| 267 sum[0] = rX + gX + bX; | |
| 268 sum[1] = rY + gY + bY; | |
| 269 sum[2] = rZ + gZ + bZ; | |
| 270 | |
| 271 // Build our target vector (see mozilla bug 460629) | |
| 272 target[0] = 0.96420; | |
| 273 target[1] = 1.00000; | |
| 274 target[2] = 0.82491; | |
| 275 | |
| 276 // Our tolerance vector - Recommended by Chris Murphy based on | |
| 277 // conversion from the LAB space criterion of no more than 3 in any one | |
| 278 // channel. This is similar to, but slightly more tolerant than Adobe's | |
| 279 // criterion. | |
| 280 tolerance[0] = 0.02; | |
| 281 tolerance[1] = 0.02; | |
| 282 tolerance[2] = 0.04; | |
| 283 | |
| 284 // Compare with our tolerance | |
| 285 for (i = 0; i < 3; ++i) { | |
| 286 if (!(((sum[i] - tolerance[i]) <= target[i]) && | |
| 287 ((sum[i] + tolerance[i]) >= target[i]))) | |
| 288 return true; | |
| 289 } | |
| 290 | |
| 291 // All Good | |
| 292 return false; | |
| 293 } | |
| 294 | |
| 295 #define TAG_bXYZ 0x6258595a | |
| 296 #define TAG_gXYZ 0x6758595a | |
| 297 #define TAG_rXYZ 0x7258595a | |
| 298 #define TAG_rTRC 0x72545243 | |
| 299 #define TAG_bTRC 0x62545243 | |
| 300 #define TAG_gTRC 0x67545243 | |
| 301 #define TAG_kTRC 0x6b545243 | |
| 302 #define TAG_A2B0 0x41324230 | |
| 303 | |
| 304 static struct tag *find_tag(struct tag_index index, uint32_t tag_id) | |
| 305 { | |
| 306 unsigned int i; | |
| 307 struct tag *tag = NULL; | |
| 308 for (i = 0; i < index.count; i++) { | |
| 309 if (index.tags[i].signature == tag_id) { | |
| 310 return &index.tags[i]; | |
| 311 } | |
| 312 } | |
| 313 return tag; | |
| 314 } | |
| 315 | |
| 316 #define XYZ_TYPE 0x58595a20 // 'XYZ ' | |
| 317 #define CURVE_TYPE 0x63757276 // 'curv' | |
| 318 #define LUT16_TYPE 0x6d667432 // 'mft2' | |
| 319 #define LUT8_TYPE 0x6d667431 // 'mft1' | |
| 320 | |
| 321 static struct XYZNumber read_tag_XYZType(struct mem_source *src, struct tag_inde
x index, uint32_t tag_id) | |
| 322 { | |
| 323 struct XYZNumber num = {0, 0, 0}; | |
| 324 struct tag *tag = find_tag(index, tag_id); | |
| 325 if (tag) { | |
| 326 uint32_t offset = tag->offset; | |
| 327 | |
| 328 uint32_t type = read_u32(src, offset); | |
| 329 if (type != XYZ_TYPE) | |
| 330 invalid_source(src, "unexpected type, expected XYZ"); | |
| 331 num.X = read_s15Fixed16Number(src, offset+8); | |
| 332 num.Y = read_s15Fixed16Number(src, offset+12); | |
| 333 num.Z = read_s15Fixed16Number(src, offset+16); | |
| 334 } else { | |
| 335 invalid_source(src, "missing xyztag"); | |
| 336 } | |
| 337 return num; | |
| 338 } | |
| 339 | |
| 340 static struct curveType *read_tag_curveType(struct mem_source *src, struct tag_i
ndex index, uint32_t tag_id) | |
| 341 { | |
| 342 struct tag *tag = find_tag(index, tag_id); | |
| 343 struct curveType *curve = NULL; | |
| 344 if (tag) { | |
| 345 uint32_t offset = tag->offset; | |
| 346 uint32_t type = read_u32(src, offset); | |
| 347 uint32_t count = read_u32(src, offset+8); | |
| 348 unsigned int i; | |
| 349 | |
| 350 if (type != CURVE_TYPE) { | |
| 351 invalid_source(src, "unexpected type, expected CURV"); | |
| 352 return NULL; | |
| 353 } | |
| 354 | |
| 355 #define MAX_CURVE_ENTRIES 40000 //arbitrary | |
| 356 if (count > MAX_CURVE_ENTRIES) { | |
| 357 invalid_source(src, "curve size too large"); | |
| 358 return NULL; | |
| 359 } | |
| 360 curve = malloc(sizeof(struct curveType) + sizeof(uInt16Number)*c
ount); | |
| 361 if (!curve) | |
| 362 return NULL; | |
| 363 | |
| 364 curve->count = count; | |
| 365 for (i=0; i<count; i++) { | |
| 366 curve->data[i] = read_u16(src, offset + 12 + i *2); | |
| 367 } | |
| 368 } else { | |
| 369 invalid_source(src, "missing curvetag"); | |
| 370 } | |
| 371 | |
| 372 return curve; | |
| 373 } | |
| 374 | |
| 375 /* This function's not done yet */ | |
| 376 static struct lutType *read_tag_lutType(struct mem_source *src, struct tag_index
index, uint32_t tag_id) | |
| 377 { | |
| 378 struct tag *tag = find_tag(index, tag_id); | |
| 379 uint32_t offset = tag->offset; | |
| 380 uint32_t type = read_u32(src, offset); | |
| 381 uint16_t num_input_table_entries; | |
| 382 uint16_t num_output_table_entries; | |
| 383 uint8_t in_chan, grid_points, out_chan; | |
| 384 uint32_t clut_size; | |
| 385 struct lutType *lut; | |
| 386 int i; | |
| 387 | |
| 388 num_input_table_entries = read_u16(src, offset + 48); | |
| 389 num_output_table_entries = read_u16(src, offset + 50); | |
| 390 | |
| 391 in_chan = read_u8(src, offset + 8); | |
| 392 out_chan = read_u8(src, offset + 9); | |
| 393 grid_points = read_u8(src, offset + 10); | |
| 394 | |
| 395 if (!src->valid) | |
| 396 return NULL; | |
| 397 | |
| 398 clut_size = in_chan * grid_points * out_chan; | |
| 399 #define MAX_CLUT_SIZE 10000 // arbitrary | |
| 400 if (clut_size > MAX_CLUT_SIZE) { | |
| 401 return NULL; | |
| 402 } | |
| 403 | |
| 404 if (type != LUT16_TYPE && type != LUT8_TYPE) | |
| 405 return NULL; | |
| 406 | |
| 407 lut = malloc(sizeof(struct lutType) + (clut_size + num_input_table_entri
es + num_output_table_entries)*sizeof(uint8_t)); | |
| 408 if (!lut) | |
| 409 return NULL; | |
| 410 lut->num_input_channels = read_u8(src, offset + 8); | |
| 411 lut->num_output_channels = read_u8(src, offset + 9); | |
| 412 lut->num_clut_grid_points = read_u8(src, offset + 10); | |
| 413 lut->e00 = read_s15Fixed16Number(src, offset+12); | |
| 414 lut->e01 = read_s15Fixed16Number(src, offset+16); | |
| 415 lut->e02 = read_s15Fixed16Number(src, offset+20); | |
| 416 lut->e10 = read_s15Fixed16Number(src, offset+24); | |
| 417 lut->e11 = read_s15Fixed16Number(src, offset+28); | |
| 418 lut->e12 = read_s15Fixed16Number(src, offset+32); | |
| 419 lut->e20 = read_s15Fixed16Number(src, offset+36); | |
| 420 lut->e21 = read_s15Fixed16Number(src, offset+40); | |
| 421 lut->e22 = read_s15Fixed16Number(src, offset+44); | |
| 422 | |
| 423 //TODO: finish up | |
| 424 return lut; | |
| 425 } | |
| 426 | |
| 427 static void read_rendering_intent(qcms_profile *profile, struct mem_source *src) | |
| 428 { | |
| 429 profile->rendering_intent = read_u32(src, 64); | |
| 430 switch (profile->rendering_intent) { | |
| 431 case QCMS_INTENT_PERCEPTUAL: | |
| 432 case QCMS_INTENT_SATURATION: | |
| 433 case QCMS_INTENT_RELATIVE_COLORIMETRIC: | |
| 434 case QCMS_INTENT_ABSOLUTE_COLORIMETRIC: | |
| 435 break; | |
| 436 default: | |
| 437 invalid_source(src, "unknown rendering intent"); | |
| 438 } | |
| 439 } | |
| 440 | |
| 441 qcms_profile *qcms_profile_create(void) | |
| 442 { | |
| 443 return calloc(sizeof(qcms_profile), 1); | |
| 444 } | |
| 445 | |
| 446 /* build sRGB gamma table */ | |
| 447 /* based on cmsBuildParametricGamma() */ | |
| 448 static uint16_t *build_sRGB_gamma_table(int num_entries) | |
| 449 { | |
| 450 int i; | |
| 451 /* taken from lcms: Build_sRGBGamma() */ | |
| 452 double gamma = 2.4; | |
| 453 double a = 1./1.055; | |
| 454 double b = 0.055/1.055; | |
| 455 double c = 1./12.92; | |
| 456 double d = 0.04045; | |
| 457 | |
| 458 uint16_t *table = malloc(sizeof(uint16_t) * num_entries); | |
| 459 if (!table) | |
| 460 return NULL; | |
| 461 | |
| 462 for (i=0; i<num_entries; i++) { | |
| 463 double x = (double)i / (num_entries-1); | |
| 464 double y, output; | |
| 465 // IEC 61966-2.1 (sRGB) | |
| 466 // Y = (aX + b)^Gamma | X >= d | |
| 467 // Y = cX | X < d | |
| 468 if (x >= d) { | |
| 469 double e = (a*x + b); | |
| 470 if (e > 0) | |
| 471 y = pow(e, gamma); | |
| 472 else | |
| 473 y = 0; | |
| 474 } else { | |
| 475 y = c*x; | |
| 476 } | |
| 477 | |
| 478 // Saturate -- this could likely move to a separate function | |
| 479 output = y * 65535. + .5; | |
| 480 if (output > 65535.) | |
| 481 output = 65535; | |
| 482 if (output < 0) | |
| 483 output = 0; | |
| 484 table[i] = (uint16_t)floor(output); | |
| 485 } | |
| 486 return table; | |
| 487 } | |
| 488 | |
| 489 static struct curveType *curve_from_table(uint16_t *table, int num_entries) | |
| 490 { | |
| 491 struct curveType *curve; | |
| 492 int i; | |
| 493 curve = malloc(sizeof(struct curveType) + sizeof(uInt16Number)*num_entri
es); | |
| 494 if (!curve) | |
| 495 return NULL; | |
| 496 curve->count = num_entries; | |
| 497 for (i = 0; i < num_entries; i++) { | |
| 498 curve->data[i] = table[i]; | |
| 499 } | |
| 500 return curve; | |
| 501 } | |
| 502 | |
| 503 static uint16_t float_to_u8Fixed8Number(float a) | |
| 504 { | |
| 505 if (a > (255. + 255./256)) | |
| 506 return 0xffff; | |
| 507 else if (a < 0.) | |
| 508 return 0; | |
| 509 else | |
| 510 return floor(a*256. + .5); | |
| 511 } | |
| 512 | |
| 513 static struct curveType *curve_from_gamma(float gamma) | |
| 514 { | |
| 515 struct curveType *curve; | |
| 516 int num_entries = 1; | |
| 517 curve = malloc(sizeof(struct curveType) + sizeof(uInt16Number)*num_entri
es); | |
| 518 if (!curve) | |
| 519 return NULL; | |
| 520 curve->count = num_entries; | |
| 521 curve->data[0] = float_to_u8Fixed8Number(gamma); | |
| 522 return curve; | |
| 523 } | |
| 524 | |
| 525 static void qcms_profile_fini(qcms_profile *profile) | |
| 526 { | |
| 527 free(profile->redTRC); | |
| 528 free(profile->blueTRC); | |
| 529 free(profile->greenTRC); | |
| 530 free(profile->grayTRC); | |
| 531 free(profile); | |
| 532 } | |
| 533 | |
| 534 //XXX: it would be nice if we had a way of ensuring | |
| 535 // everything in a profile was initialized regardless of how it was created | |
| 536 | |
| 537 //XXX: should this also be taking a black_point? | |
| 538 /* similar to CGColorSpaceCreateCalibratedRGB */ | |
| 539 qcms_profile* qcms_profile_create_rgb_with_gamma( | |
| 540 qcms_CIE_xyY white_point, | |
| 541 qcms_CIE_xyYTRIPLE primaries, | |
| 542 float gamma) | |
| 543 { | |
| 544 qcms_profile* profile = qcms_profile_create(); | |
| 545 if (!profile) | |
| 546 return NO_MEM_PROFILE; | |
| 547 | |
| 548 //XXX: should store the whitepoint | |
| 549 if (!set_rgb_colorants(profile, white_point, primaries)) { | |
| 550 qcms_profile_fini(profile); | |
| 551 return INVALID_PROFILE; | |
| 552 } | |
| 553 | |
| 554 profile->redTRC = curve_from_gamma(gamma); | |
| 555 profile->blueTRC = curve_from_gamma(gamma); | |
| 556 profile->greenTRC = curve_from_gamma(gamma); | |
| 557 | |
| 558 if (!profile->redTRC || !profile->blueTRC || !profile->greenTRC) { | |
| 559 qcms_profile_fini(profile); | |
| 560 return NO_MEM_PROFILE; | |
| 561 } | |
| 562 profile->class = DISPLAY_DEVICE_PROFILE; | |
| 563 profile->rendering_intent = QCMS_INTENT_PERCEPTUAL; | |
| 564 profile->color_space = RGB_SIGNATURE; | |
| 565 return profile; | |
| 566 } | |
| 567 | |
| 568 qcms_profile* qcms_profile_create_rgb_with_table( | |
| 569 qcms_CIE_xyY white_point, | |
| 570 qcms_CIE_xyYTRIPLE primaries, | |
| 571 uint16_t *table, int num_entries) | |
| 572 { | |
| 573 qcms_profile* profile = qcms_profile_create(); | |
| 574 if (!profile) | |
| 575 return NO_MEM_PROFILE; | |
| 576 | |
| 577 //XXX: should store the whitepoint | |
| 578 if (!set_rgb_colorants(profile, white_point, primaries)) { | |
| 579 qcms_profile_fini(profile); | |
| 580 return INVALID_PROFILE; | |
| 581 } | |
| 582 | |
| 583 profile->redTRC = curve_from_table(table, num_entries); | |
| 584 profile->blueTRC = curve_from_table(table, num_entries); | |
| 585 profile->greenTRC = curve_from_table(table, num_entries); | |
| 586 | |
| 587 if (!profile->redTRC || !profile->blueTRC || !profile->greenTRC) { | |
| 588 qcms_profile_fini(profile); | |
| 589 return NO_MEM_PROFILE; | |
| 590 } | |
| 591 profile->class = DISPLAY_DEVICE_PROFILE; | |
| 592 profile->rendering_intent = QCMS_INTENT_PERCEPTUAL; | |
| 593 profile->color_space = RGB_SIGNATURE; | |
| 594 return profile; | |
| 595 } | |
| 596 | |
| 597 /* from lcms: cmsWhitePointFromTemp */ | |
| 598 /* tempK must be >= 4000. and <= 25000. | |
| 599 * similar to argyll: icx_DTEMP2XYZ() */ | |
| 600 static qcms_CIE_xyY white_point_from_temp(int temp_K) | |
| 601 { | |
| 602 qcms_CIE_xyY white_point; | |
| 603 double x, y; | |
| 604 double T, T2, T3; | |
| 605 // double M1, M2; | |
| 606 | |
| 607 // No optimization provided. | |
| 608 T = temp_K; | |
| 609 T2 = T*T; // Square | |
| 610 T3 = T2*T; // Cube | |
| 611 | |
| 612 // For correlated color temperature (T) between 4000K and 7000K: | |
| 613 if (T >= 4000. && T <= 7000.) { | |
| 614 x = -4.6070*(1E9/T3) + 2.9678*(1E6/T2) + 0.09911*(1E3/T) + 0.244
063; | |
| 615 } else { | |
| 616 // or for correlated color temperature (T) between 7000K and 250
00K: | |
| 617 if (T > 7000.0 && T <= 25000.0) { | |
| 618 x = -2.0064*(1E9/T3) + 1.9018*(1E6/T2) + 0.24748*(1E3/T)
+ 0.237040; | |
| 619 } else { | |
| 620 assert(0 && "invalid temp"); | |
| 621 } | |
| 622 } | |
| 623 | |
| 624 // Obtain y(x) | |
| 625 | |
| 626 y = -3.000*(x*x) + 2.870*x - 0.275; | |
| 627 | |
| 628 // wave factors (not used, but here for futures extensions) | |
| 629 | |
| 630 // M1 = (-1.3515 - 1.7703*x + 5.9114 *y)/(0.0241 + 0.2562*x - 0.7341*y); | |
| 631 // M2 = (0.0300 - 31.4424*x + 30.0717*y)/(0.0241 + 0.2562*x - 0.7341*y); | |
| 632 | |
| 633 // Fill white_point struct | |
| 634 white_point.x = x; | |
| 635 white_point.y = y; | |
| 636 white_point.Y = 1.0; | |
| 637 | |
| 638 return white_point; | |
| 639 } | |
| 640 | |
| 641 qcms_profile* qcms_profile_sRGB(void) | |
| 642 { | |
| 643 qcms_profile *profile; | |
| 644 uint16_t *table; | |
| 645 | |
| 646 qcms_CIE_xyYTRIPLE Rec709Primaries = { | |
| 647 {0.6400, 0.3300, 1.0}, | |
| 648 {0.3000, 0.6000, 1.0}, | |
| 649 {0.1500, 0.0600, 1.0} | |
| 650 }; | |
| 651 qcms_CIE_xyY D65; | |
| 652 | |
| 653 D65 = white_point_from_temp(6504); | |
| 654 | |
| 655 table = build_sRGB_gamma_table(1024); | |
| 656 | |
| 657 if (!table) | |
| 658 return NO_MEM_PROFILE; | |
| 659 | |
| 660 profile = qcms_profile_create_rgb_with_table(D65, Rec709Primaries, table
, 1024); | |
| 661 free(table); | |
| 662 return profile; | |
| 663 } | |
| 664 | |
| 665 | |
| 666 /* qcms_profile_from_memory does not hold a reference to the memory passed in */ | |
| 667 qcms_profile* qcms_profile_from_memory(const void *mem, size_t size) | |
| 668 { | |
| 669 uint32_t length; | |
| 670 struct mem_source source; | |
| 671 struct mem_source *src = &source; | |
| 672 struct tag_index index; | |
| 673 qcms_profile *profile; | |
| 674 | |
| 675 source.buf = mem; | |
| 676 source.size = size; | |
| 677 source.valid = true; | |
| 678 | |
| 679 length = read_u32(src, 0); | |
| 680 if (length <= size) { | |
| 681 // shrink the area that we can read if appropriate | |
| 682 source.size = length; | |
| 683 } else { | |
| 684 return INVALID_PROFILE; | |
| 685 } | |
| 686 | |
| 687 /* ensure that the profile size is sane so it's easier to reason about *
/ | |
| 688 if (source.size <= 64 || source.size >= MAX_PROFILE_SIZE) | |
| 689 return INVALID_PROFILE; | |
| 690 | |
| 691 profile = qcms_profile_create(); | |
| 692 if (!profile) | |
| 693 return NO_MEM_PROFILE; | |
| 694 | |
| 695 check_CMM_type_signature(src); | |
| 696 check_profile_version(src); | |
| 697 read_class_signature(profile, src); | |
| 698 read_rendering_intent(profile, src); | |
| 699 read_color_space(profile, src); | |
| 700 //TODO read rest of profile stuff | |
| 701 | |
| 702 if (!src->valid) | |
| 703 goto invalid_profile; | |
| 704 | |
| 705 index = read_tag_table(profile, src); | |
| 706 if (!src->valid || !index.tags) | |
| 707 goto invalid_tag_table; | |
| 708 | |
| 709 if (profile->class == DISPLAY_DEVICE_PROFILE || profile->class == INPUT_
DEVICE_PROFILE) { | |
| 710 if (profile->color_space == RGB_SIGNATURE) { | |
| 711 | |
| 712 profile->redColorant = read_tag_XYZType(src, index, TAG_
rXYZ); | |
| 713 profile->blueColorant = read_tag_XYZType(src, index, TAG
_bXYZ); | |
| 714 profile->greenColorant = read_tag_XYZType(src, index, TA
G_gXYZ); | |
| 715 | |
| 716 if (!src->valid) | |
| 717 goto invalid_tag_table; | |
| 718 | |
| 719 profile->redTRC = read_tag_curveType(src, index, TAG_rTR
C); | |
| 720 profile->blueTRC = read_tag_curveType(src, index, TAG_bT
RC); | |
| 721 profile->greenTRC = read_tag_curveType(src, index, TAG_g
TRC); | |
| 722 | |
| 723 if (!profile->redTRC || !profile->blueTRC || !profile->g
reenTRC) | |
| 724 goto invalid_tag_table; | |
| 725 | |
| 726 } else if (profile->color_space == GRAY_SIGNATURE) { | |
| 727 | |
| 728 profile->grayTRC = read_tag_curveType(src, index, TAG_kT
RC); | |
| 729 if (!profile->grayTRC) | |
| 730 goto invalid_tag_table; | |
| 731 | |
| 732 } else { | |
| 733 goto invalid_tag_table; | |
| 734 } | |
| 735 } else if (0 && profile->class == OUTPUT_DEVICE_PROFILE) { | |
| 736 profile->A2B0 = read_tag_lutType(src, index, TAG_A2B0); | |
| 737 } else { | |
| 738 goto invalid_tag_table; | |
| 739 } | |
| 740 | |
| 741 if (!src->valid) | |
| 742 goto invalid_tag_table; | |
| 743 | |
| 744 free(index.tags); | |
| 745 | |
| 746 return profile; | |
| 747 | |
| 748 invalid_tag_table: | |
| 749 free(index.tags); | |
| 750 invalid_profile: | |
| 751 qcms_profile_fini(profile); | |
| 752 return INVALID_PROFILE; | |
| 753 } | |
| 754 | |
| 755 qcms_intent qcms_profile_get_rendering_intent(qcms_profile *profile) | |
| 756 { | |
| 757 return profile->rendering_intent; | |
| 758 } | |
| 759 | |
| 760 icColorSpaceSignature | |
| 761 qcms_profile_get_color_space(qcms_profile *profile) | |
| 762 { | |
| 763 return profile->color_space; | |
| 764 } | |
| 765 | |
| 766 void qcms_profile_release(qcms_profile *profile) | |
| 767 { | |
| 768 if (profile->output_table_r) | |
| 769 precache_release(profile->output_table_r); | |
| 770 if (profile->output_table_g) | |
| 771 precache_release(profile->output_table_g); | |
| 772 if (profile->output_table_b) | |
| 773 precache_release(profile->output_table_b); | |
| 774 | |
| 775 qcms_profile_fini(profile); | |
| 776 } | |
| 777 | |
| 778 #include <stdio.h> | |
| 779 qcms_profile* qcms_profile_from_file(FILE *file) | |
| 780 { | |
| 781 uint32_t length, remaining_length; | |
| 782 qcms_profile *profile; | |
| 783 size_t read_length; | |
| 784 be32 length_be; | |
| 785 void *data; | |
| 786 | |
| 787 fread(&length_be, sizeof(length), 1, file); | |
| 788 length = be32_to_cpu(length_be); | |
| 789 if (length > MAX_PROFILE_SIZE) | |
| 790 return BAD_VALUE_PROFILE; | |
| 791 | |
| 792 /* allocate room for the entire profile */ | |
| 793 data = malloc(length); | |
| 794 if (!data) | |
| 795 return NO_MEM_PROFILE; | |
| 796 | |
| 797 /* copy in length to the front so that the buffer will contain the entir
e profile */ | |
| 798 *((be32*)data) = length_be; | |
| 799 remaining_length = length - sizeof(length_be); | |
| 800 | |
| 801 /* read the rest profile */ | |
| 802 read_length = fread((unsigned char*)data + sizeof(length_be), 1, remaini
ng_length, file); | |
| 803 if (read_length != remaining_length) { | |
| 804 free(data); | |
| 805 return INVALID_PROFILE; | |
| 806 } | |
| 807 | |
| 808 profile = qcms_profile_from_memory(data, length); | |
| 809 free(data); | |
| 810 return profile; | |
| 811 } | |
| 812 | |
| 813 qcms_profile* qcms_profile_from_path(const char *path) | |
| 814 { | |
| 815 qcms_profile *profile = NULL; | |
| 816 FILE *file = fopen(path, "rb"); | |
| 817 if (file) { | |
| 818 profile = qcms_profile_from_file(file); | |
| 819 fclose(file); | |
| 820 } | |
| 821 return profile; | |
| 822 } | |
| OLD | NEW |