429 lines
14 KiB
C++
429 lines
14 KiB
C++
/*
|
|
* This file is part of RawTherapee.
|
|
*
|
|
* Copyright (c) 2004-2010 Gabor Horvath <hgabor@rawtherapee.com>
|
|
*
|
|
* RawTherapee is free software: you can redistribute it and/or modify
|
|
* it under the terms of the GNU General Public License as published by
|
|
* the Free Software Foundation, either version 3 of the License, or
|
|
* (at your option) any later version.
|
|
*
|
|
* RawTherapee is distributed in the hope that it will be useful,
|
|
* but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
* GNU General Public License for more details.
|
|
*
|
|
* You should have received a copy of the GNU General Public License
|
|
* along with RawTherapee. If not, see <http://www.gnu.org/licenses/>.
|
|
*/
|
|
#include "rtengine.h"
|
|
#include "improcfun.h"
|
|
#include <glibmm.h>
|
|
#include "iccstore.h"
|
|
#include "iccmatrices.h"
|
|
#include "../rtgui/options.h"
|
|
#include "settings.h"
|
|
#include "curves.h"
|
|
#include "alignedbuffer.h"
|
|
#include "color.h"
|
|
|
|
namespace rtengine
|
|
{
|
|
|
|
extern const Settings* settings;
|
|
|
|
// Used in ImProcCoordinator::updatePreviewImage (rtengine/improccoordinator.cc)
|
|
// Crop::update (rtengine/dcrop.cc)
|
|
// Thumbnail::processImage (rtengine/rtthumbnail.cc)
|
|
//
|
|
// If monitorTransform, divide by 327.68 then apply monitorTransform (which can integrate soft-proofing)
|
|
// otherwise divide by 327.68, convert to xyz and apply the sRGB transform, before converting with gamma2curve
|
|
void ImProcFunctions::lab2monitorRgb (LabImage* lab, Image8* image)
|
|
{
|
|
if (monitorTransform) {
|
|
|
|
int W = lab->W;
|
|
int H = lab->H;
|
|
unsigned char * data = image->data;
|
|
|
|
// cmsDoTransform is relatively expensive
|
|
#ifdef _OPENMP
|
|
#pragma omp parallel firstprivate(lab, data, W, H)
|
|
#endif
|
|
{
|
|
AlignedBuffer<float> pBuf(3 * lab->W);
|
|
float *buffer = pBuf.data;
|
|
|
|
#ifdef _OPENMP
|
|
#pragma omp for schedule(dynamic,16)
|
|
#endif
|
|
|
|
for (int i = 0; i < H; i++) {
|
|
|
|
const int ix = i * 3 * W;
|
|
int iy = 0;
|
|
|
|
float* rL = lab->L[i];
|
|
float* ra = lab->a[i];
|
|
float* rb = lab->b[i];
|
|
|
|
for (int j = 0; j < W; j++) {
|
|
buffer[iy++] = rL[j] / 327.68f;
|
|
buffer[iy++] = ra[j] / 327.68f;
|
|
buffer[iy++] = rb[j] / 327.68f;
|
|
}
|
|
|
|
cmsDoTransform (monitorTransform, buffer, data + ix, W);
|
|
}
|
|
|
|
} // End of parallelization
|
|
|
|
} else {
|
|
|
|
int W = lab->W;
|
|
int H = lab->H;
|
|
unsigned char * data = image->data;
|
|
|
|
#ifdef _OPENMP
|
|
#pragma omp parallel for schedule(dynamic,16) if (multiThread)
|
|
#endif
|
|
|
|
for (int i = 0; i < H; i++) {
|
|
float* rL = lab->L[i];
|
|
float* ra = lab->a[i];
|
|
float* rb = lab->b[i];
|
|
int ix = i * 3 * W;
|
|
|
|
float R, G, B;
|
|
float fy, fx, fz, x_, y_, z_, LL;
|
|
|
|
for (int j = 0; j < W; j++) {
|
|
|
|
//float L1=rL[j],a1=ra[j],b1=rb[j];//for testing
|
|
|
|
fy = (0.00862069 * rL[j]) / 327.68 + 0.137932; // (L+16)/116
|
|
fx = (0.002 * ra[j]) / 327.68 + fy;
|
|
fz = fy - (0.005 * rb[j]) / 327.68;
|
|
LL = rL[j] / 327.68;
|
|
|
|
x_ = 65535.0 * Color::f2xyz(fx) * Color::D50x;
|
|
// y_ = 65535.0 * Color::f2xyz(fy);
|
|
z_ = 65535.0 * Color::f2xyz(fz) * Color::D50z;
|
|
y_ = (LL > (float)Color::epskap) ? 65535.0 * fy * fy * fy : 65535.0 * LL / (float)Color::kappa;
|
|
|
|
Color::xyz2srgb(x_, y_, z_, R, G, B);
|
|
|
|
/* copy RGB */
|
|
//int R1=((int)gamma2curve[(R)])
|
|
data[ix++] = ((int)Color::gamma2curve[R]) >> 8;
|
|
data[ix++] = ((int)Color::gamma2curve[G]) >> 8;
|
|
data[ix++] = ((int)Color::gamma2curve[B]) >> 8;
|
|
}
|
|
}
|
|
}
|
|
}
|
|
|
|
|
|
|
|
// Used in ImProcCoordinator::updatePreviewImage (rtengine/improccoordinator.cc)
|
|
// Crop::update (rtengine/dcrop.cc)
|
|
//
|
|
// Generate an Image8
|
|
//
|
|
// If output profile used, divide by 327.68 then apply the "profile" profile (eventually with a standard gamma)
|
|
// otherwise divide by 327.68, convert to xyz and apply the RGB transform, before converting with gamma2curve
|
|
Image8* ImProcFunctions::lab2rgb (LabImage* lab, int cx, int cy, int cw, int ch, const procparams::ColorManagementParams &icm)
|
|
{
|
|
//gamutmap(lab);
|
|
|
|
if (cx < 0) {
|
|
cx = 0;
|
|
}
|
|
|
|
if (cy < 0) {
|
|
cy = 0;
|
|
}
|
|
|
|
if (cx + cw > lab->W) {
|
|
cw = lab->W - cx;
|
|
}
|
|
|
|
if (cy + ch > lab->H) {
|
|
ch = lab->H - cy;
|
|
}
|
|
|
|
Image8* image = new Image8 (cw, ch);
|
|
Glib::ustring profile;
|
|
cmsHPROFILE oprof = iccStore->getProfile (profile);
|
|
|
|
bool standard_gamma;
|
|
|
|
if(settings->HistogramWorking) {
|
|
profile = icm.working;
|
|
standard_gamma = true;
|
|
} else {
|
|
profile = icm.output;
|
|
if (icm.output.empty() || icm.output == ColorManagementParams::NoICMString) {
|
|
profile = "sRGB";
|
|
}
|
|
standard_gamma = false;
|
|
}
|
|
|
|
if (oprof) {
|
|
cmsHPROFILE oprofG = oprof;
|
|
|
|
if (standard_gamma) {
|
|
oprofG = ICCStore::makeStdGammaProfile(oprof);
|
|
}
|
|
|
|
cmsUInt32Number flags = cmsFLAGS_NOOPTIMIZE | cmsFLAGS_NOCACHE;
|
|
if (icm.outputBPC) {
|
|
flags |= cmsFLAGS_BLACKPOINTCOMPENSATION;
|
|
printf("lab2rgb / bpc=true\n");
|
|
}
|
|
else printf("lab2rgb / bpc=false\n");
|
|
lcmsMutex->lock ();
|
|
cmsHPROFILE LabIProf = cmsCreateLab4Profile(NULL);
|
|
cmsHTRANSFORM hTransform = cmsCreateTransform (LabIProf, TYPE_Lab_DBL, oprofG, TYPE_RGB_8, icm.outputIntent, flags); // NOCACHE is important for thread safety
|
|
cmsCloseProfile(LabIProf);
|
|
lcmsMutex->unlock ();
|
|
|
|
unsigned char *data = image->data;
|
|
|
|
// cmsDoTransform is relatively expensive
|
|
#ifdef _OPENMP
|
|
#pragma omp parallel
|
|
#endif
|
|
{
|
|
AlignedBuffer<double> pBuf(3 * cw);
|
|
double *buffer = pBuf.data;
|
|
int condition = cy + ch;
|
|
|
|
#ifdef _OPENMP
|
|
#pragma omp for firstprivate(lab) schedule(dynamic,16)
|
|
#endif
|
|
|
|
for (int i = cy; i < condition; i++) {
|
|
const int ix = i * 3 * cw;
|
|
int iy = 0;
|
|
float* rL = lab->L[i];
|
|
float* ra = lab->a[i];
|
|
float* rb = lab->b[i];
|
|
|
|
for (int j = cx; j < cx + cw; j++) {
|
|
buffer[iy++] = rL[j] / 327.68f;
|
|
buffer[iy++] = ra[j] / 327.68f;
|
|
buffer[iy++] = rb[j] / 327.68f;
|
|
}
|
|
|
|
cmsDoTransform (hTransform, buffer, data + ix, cw);
|
|
}
|
|
} // End of parallelization
|
|
|
|
cmsDeleteTransform(hTransform);
|
|
|
|
if (oprofG != oprof) {
|
|
cmsCloseProfile(oprofG);
|
|
}
|
|
} else {
|
|
|
|
const auto xyz_rgb = iccStore->workingSpaceInverseMatrix (profile);
|
|
|
|
#ifdef _OPENMP
|
|
#pragma omp parallel for schedule(dynamic,16) if (multiThread)
|
|
#endif
|
|
|
|
for (int i = cy; i < cy + ch; i++) {
|
|
float R, G, B;
|
|
float* rL = lab->L[i];
|
|
float* ra = lab->a[i];
|
|
float* rb = lab->b[i];
|
|
uint8_t* dest = image->r(i - cy) + cx;
|
|
|
|
for (int j = cx; j < cx + cw; j++) {
|
|
|
|
float fy = (0.0086206897f * (*rL)) / 327.68f + 0.1379310345f; // (L+16)/116
|
|
float fx = (0.002f * *(ra++)) / 327.68f + fy;
|
|
float fz = fy - (0.005f * *(rb++)) / 327.68f;
|
|
float LL = *(rL++) / 327.68f;
|
|
|
|
float x_ = 65535.0f * Color::f2xyz(fx) * Color::D50x;
|
|
//float y_ = 65535.0f * Color::f2xyz(fy);
|
|
float z_ = 65535.0f * Color::f2xyz(fz) * Color::D50z;
|
|
float y_ = (LL > (float)Color::epskap) ? 65535.0f * fy * fy * fy : 65535.0f * LL / (float)Color::kappa;
|
|
|
|
Color::xyz2rgb(x_, y_, z_, R, G, B, xyz_rgb);
|
|
|
|
*(dest++) = (int)Color::gamma2curve[CLIP(R)] >> 8;
|
|
*(dest++) = (int)Color::gamma2curve[CLIP(G)] >> 8;
|
|
*(dest++) = (int)Color::gamma2curve[CLIP(B)] >> 8;
|
|
}
|
|
}
|
|
}
|
|
|
|
return image;
|
|
}
|
|
|
|
|
|
/** @brief Convert the final Lab image to the output RGB color space
|
|
*
|
|
* Used in processImage (rtengine/simpleprocess.cc)
|
|
*
|
|
* Provide a pointer to a 7 floats array for "ga" (uninitialized ; this array will be filled with the gamma values) if you want
|
|
* to use the custom gamma scenario. Thoses gamma values will correspond to the ones of the chosen standard output profile
|
|
* (Prophoto if non standard output profile given)
|
|
*
|
|
* If "ga" is NULL, then we're considering standard gamma with the chosen output profile.
|
|
*
|
|
* Generate an Image16
|
|
*
|
|
* If a custom gamma profile can be created, divide by 327.68, convert to xyz and apply the custom gamma transform
|
|
* otherwise divide by 327.68, convert to xyz and apply the sRGB transform, before converting with gamma2curve
|
|
*/
|
|
Image16* ImProcFunctions::lab2rgb16 (LabImage* lab, int cx, int cy, int cw, int ch, const procparams::ColorManagementParams &icm, bool bw, GammaValues *ga)
|
|
{
|
|
|
|
//gamutmap(lab);
|
|
|
|
if (cx < 0) {
|
|
cx = 0;
|
|
}
|
|
|
|
if (cy < 0) {
|
|
cy = 0;
|
|
}
|
|
|
|
if (cx + cw > lab->W) {
|
|
cw = lab->W - cx;
|
|
}
|
|
|
|
if (cy + ch > lab->H) {
|
|
ch = lab->H - cy;
|
|
}
|
|
|
|
Image16* image = new Image16 (cw, ch);
|
|
|
|
cmsHPROFILE oprof = NULL;
|
|
if (ga) {
|
|
iccStore->getGammaArray(icm, *ga);
|
|
oprof = iccStore->createGammaProfile(icm, *ga);
|
|
printf("iccStore->createGammaProfile(icm, *ga);\n");
|
|
} else {
|
|
oprof = iccStore->getProfile (icm.output);
|
|
printf("iccStore->getProfile (%s);\n", icm.output.c_str());
|
|
}
|
|
|
|
if (oprof) {
|
|
|
|
|
|
/*
|
|
|
|
|
|
#ifdef _OPENMP
|
|
#pragma omp parallel for schedule(dynamic,16) if (multiThread)
|
|
#endif
|
|
|
|
for (int i = cy; i < cy + ch; i++) {
|
|
float* rL = lab->L[i];
|
|
float* ra = lab->a[i];
|
|
float* rb = lab->b[i];
|
|
short* xa = (short*)image->r(i - cy) + cx;
|
|
short* ya = (short*)image->g(i - cy) + cx;
|
|
short* za = (short*)image->b(i - cy) + cx;
|
|
|
|
for (int j = 0; j < cw; j++) {
|
|
|
|
float fy = (0.0086206897f * *rL) / 327.68f + 0.1379310345f; // (L+16)/116
|
|
float fx = (0.002f * *(ra++)) / 327.68f + fy;
|
|
float fz = fy - (0.005f * *(rb++)) / 327.68f;
|
|
float LL = *(rL++) / 327.68f;
|
|
|
|
float x_ = 65535.0f * Color::f2xyz(fx) * Color::D50x;
|
|
//float y_ = 65535.0f * Color::f2xyz(fy);
|
|
float z_ = 65535.0f * Color::f2xyz(fz) * Color::D50z;
|
|
float y_ = (LL > (float)Color::epskap) ? 65535.0f * fy * fy * fy : 65535.0f * LL / (float)Color::kappa;
|
|
|
|
*xa = CLIP((int) round(x_)) ;
|
|
*(ya++) = CLIP((int) round(y_));
|
|
*za = CLIP((int) round(z_));
|
|
|
|
if(bw && y_ < 65535.f) { //force Bw value and take highlight into account
|
|
*xa = (int) round(y_ * Color::D50x);
|
|
*za = (int) round(y_ * Color::D50z);
|
|
}
|
|
++xa;
|
|
++za;
|
|
}
|
|
}
|
|
|
|
cmsHPROFILE iprof = iccStore->getXYZProfile ();
|
|
|
|
|
|
*/
|
|
|
|
cmsHPROFILE iprof = cmsCreateLab4Profile(nullptr);
|
|
|
|
|
|
|
|
|
|
// ----------------------------------------------------------------------------
|
|
|
|
|
|
|
|
|
|
cmsUInt32Number flags = cmsFLAGS_NOOPTIMIZE | cmsFLAGS_NOCACHE;
|
|
if (icm.outputBPC) {
|
|
flags |= cmsFLAGS_BLACKPOINTCOMPENSATION;
|
|
printf("lab2rgb16 / icm.outputBPC=true / outputIntent=%d\n", icm.outputIntent);
|
|
}
|
|
else printf("lab2rgb16 / icm.outputBPC=false / outputIntent=%d\n", icm.outputIntent);
|
|
lcmsMutex->lock ();
|
|
//cmsHTRANSFORM hTransform = cmsCreateTransform (iprof, TYPE_RGB_16, oprof, TYPE_RGB_16, icm.outputIntent, flags);
|
|
cmsHTRANSFORM hTransform = cmsCreateTransform (iprof, TYPE_Lab_FLT, oprof, TYPE_RGB_16, icm.outputIntent, flags);
|
|
lcmsMutex->unlock ();
|
|
|
|
//image->ExecCMSTransform(hTransform);
|
|
image->ExecCMSTransform(hTransform, *lab);
|
|
cmsDeleteTransform(hTransform);
|
|
} else {
|
|
//
|
|
#ifdef _OPENMP
|
|
#pragma omp parallel for schedule(dynamic,16) if (multiThread)
|
|
#endif
|
|
for (int i = cy; i < cy + ch; i++) {
|
|
float R, G, B;
|
|
float* rL = lab->L[i];
|
|
float* ra = lab->a[i];
|
|
float* rb = lab->b[i];
|
|
uint16_t* rR = image->r(i - cy) + cx;
|
|
uint16_t* rG = image->g(i - cy) + cx;
|
|
uint16_t* rB = image->b(i - cy) + cx;
|
|
|
|
for (int j = 0; j < cw; j++) {
|
|
|
|
float fy = (0.0086206897f * *rL) / 327.68f + 0.1379310345f; // (L+16)/116
|
|
float fx = (0.002f * *(ra++)) / 327.68f + fy;
|
|
float fz = fy - (0.005f * *(rb++)) / 327.68f;
|
|
float LL = *(rL++) / 327.68f;
|
|
|
|
float x_ = 65535.0f * Color::f2xyz(fx) * Color::D50x;
|
|
//float y_ = 65535.0 * Color::f2xyz(fy);
|
|
float z_ = 65535.0f * Color::f2xyz(fz) * Color::D50z;
|
|
float y_ = (LL > (float)Color::epskap) ? 65535.0f * fy * fy * fy : 65535.0f * LL / (float)Color::kappa;
|
|
|
|
Color::xyz2srgb(x_, y_, z_, R, G, B);
|
|
|
|
*(rR++) = (int)Color::gamma2curve[CLIP(R)];
|
|
*(rG++) = (int)Color::gamma2curve[CLIP(G)];
|
|
*(rB++) = (int)Color::gamma2curve[CLIP(B)];
|
|
}
|
|
}
|
|
}
|
|
|
|
return image;
|
|
}
|
|
|
|
//#include "sRGBgamutbdy.cc"
|
|
|
|
}
|