liz/rawTherapee
liz/rawTherapee
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Improve action slider chroma mask Retinex

Browse Source
This commit is contained in:
Desmis 2019-11-15 07:25:07 +01:00
parent ce282f5c97
commit abf3ca53dc
1 changed files with 40 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

45
rtengine/ipretinex.cc
View File

@ -811,10 +811,45 @@ void ImProcFunctions::maskforretinex(int sp, int before, float ** luminance, flo
}
}
float chromult1 = 1.f + 0.01f * chro;
float chromult2 = 1.f - 0.01f * chro;
float fab = chromult1 * 20000.f;//value must be good in most cases
if(fab < 50.f) fab = 50.f;
float chromult = 1.f - 0.01f * chro;
//fab
float fab = 50.f;
float meanfab = 0.f;
const int nbfab = W_L * H_L;
double sumab = 0.0;
#ifdef _OPENMP
#pragma omp parallel for reduction(+:sumab)
#endif
for (int y = 0; y < H_L; y++) {
for (int x = 0; x < W_L; x++) {
sumab += fabs(bufreti->a[y][x]);
sumab += fabs(bufreti->b[y][x]);
}
}
meanfab = sumab / (2.f * nbfab);
double som = 0.0;
#ifdef _OPENMP
#pragma omp parallel for reduction(+:som)
#endif
for (int y = 0; y < H_L; y++) {
for (int x = 0; x < W_L; x++) {
som += SQR(fabs(bufreti->a[y][x]) - meanfab) + SQR(fabs(bufreti->b[y][x]) - meanfab);
}
}
const float multsigma = (chro >= 0.f ? 0.035f : 0.018f) * chro + 1.f;
const float stddv = sqrt(som / nbfab);
fab = meanfab + multsigma * stddv;
if (fab <= 0.f) {
fab = 50.f;
}
//end fab
#ifdef _OPENMP
#pragma omp parallel for schedule(dynamic,16)
@ -854,7 +889,7 @@ void ImProcFunctions::maskforretinex(int sp, int before, float ** luminance, flo
float valHH = LIM01(1.f - lochhmasretiCurve[500.f * h]);
if (llretiMask != 4) {
kmaskCH += chromult2 * valHH;
kmaskCH += chromult * valHH;
}
kmaskLexp += 32768.f * valHH;