Init locallabgtk3
This commit is contained in:
Desmis
@@ -42,10 +42,14 @@ namespace rtengine
|
||||
{
|
||||
|
||||
static const int maxlevel = 6;
|
||||
static const int maxlevelloc = 5;
|
||||
|
||||
static const float noise = 2000;
|
||||
|
||||
//sequence of scales
|
||||
static const int scales[6] = {1, 2, 4, 8, 16, 32};
|
||||
static const int scalesloc[5] = {1, 2, 4, 8, 16};
|
||||
|
||||
extern const Settings* settings;
|
||||
|
||||
//sequence of scales
|
||||
@@ -257,6 +261,212 @@ SSEFUNCTION void ImProcFunctions :: dirpyr_equalizer(float ** src, float ** dst,
|
||||
|
||||
}
|
||||
|
||||
SSEFUNCTION void ImProcFunctions :: cbdl_local_temp(float ** src, float ** dst, float ** loctemp, int srcwidth, int srcheight, const float * mult, const double dirpyrThreshold, const double skinprot, const bool gamutlab, float b_l, float t_l, float t_r, float b_r, int choice, int scaleprev)
|
||||
{
|
||||
int lastlevel = maxlevelloc;
|
||||
|
||||
if(settings->verbose) {
|
||||
printf("Dirpyr scaleprev=%i\n", scaleprev);
|
||||
}
|
||||
|
||||
float atten123 = (float) settings->level123_cbdl;
|
||||
|
||||
if(atten123 > 50.f) {
|
||||
atten123 = 50.f;
|
||||
}
|
||||
|
||||
if(atten123 < 0.f) {
|
||||
atten123 = 0.f;
|
||||
}
|
||||
|
||||
float atten0 = (float) settings->level0_cbdl;
|
||||
|
||||
if(atten0 > 40.f) {
|
||||
atten123 = 40.f;
|
||||
}
|
||||
|
||||
if(atten0 < 0.f) {
|
||||
atten0 = 0.f;
|
||||
}
|
||||
|
||||
if((t_r - t_l) < 0.55f) {
|
||||
t_l = t_r + 0.55f; //avoid too small range
|
||||
}
|
||||
|
||||
while (lastlevel > 0 && fabs(mult[lastlevel - 1] - 1) < 0.001) {
|
||||
|
||||
lastlevel--;
|
||||
//printf("last level to process %d \n",lastlevel);
|
||||
}
|
||||
|
||||
if (lastlevel == 0) {
|
||||
return;
|
||||
}
|
||||
|
||||
int level;
|
||||
float multi[5] = {1.f, 1.f, 1.f, 1.f, 1.f};
|
||||
float scalefl[5];
|
||||
|
||||
for(int lv = 0; lv < 5; lv++) {
|
||||
scalefl[lv] = ((float) scalesloc[lv]) / (float) scaleprev;
|
||||
|
||||
if(lv >= 1) {
|
||||
if(scalefl[lv] < 1.f) {
|
||||
multi[lv] = (atten123 * ((float) mult[lv] - 1.f) / 100.f) + 1.f; //modulate action if zoom < 100%
|
||||
} else {
|
||||
multi[lv] = (float) mult[lv];
|
||||
}
|
||||
} else {
|
||||
if(scalefl[lv] < 1.f) {
|
||||
multi[lv] = (atten0 * ((float) mult[lv] - 1.f) / 100.f) + 1.f; //modulate action if zoom < 100%
|
||||
} else {
|
||||
multi[lv] = (float) mult[lv];
|
||||
}
|
||||
}
|
||||
|
||||
}
|
||||
|
||||
if(settings->verbose) {
|
||||
printf("CbDL local mult0=%f 1=%f 2=%f 3=%f 4=%f\n", multi[0], multi[1], multi[2], multi[3], multi[4]);
|
||||
}
|
||||
|
||||
multi_array2D<float, maxlevelloc> dirpyrlo (srcwidth, srcheight);
|
||||
|
||||
level = 0;
|
||||
|
||||
//int thresh = 100 * mult[5];
|
||||
int scale = (int)(scalesloc[level]) / scaleprev;
|
||||
|
||||
if(scale < 1) {
|
||||
scale = 1;
|
||||
}
|
||||
|
||||
|
||||
dirpyr_channel(src, dirpyrlo[0], srcwidth, srcheight, 0, scale);
|
||||
|
||||
level = 1;
|
||||
|
||||
while(level < lastlevel) {
|
||||
|
||||
scale = (int)(scalesloc[level]) / scaleprev;
|
||||
|
||||
if(scale < 1) {
|
||||
scale = 1;
|
||||
}
|
||||
|
||||
dirpyr_channel(dirpyrlo[level - 1], dirpyrlo[level], srcwidth, srcheight, level, scale);
|
||||
|
||||
level ++;
|
||||
}
|
||||
|
||||
float **tmpHue, **tmpChr;
|
||||
/*
|
||||
if(skinprot != 0.f) {
|
||||
// precalculate hue and chroma, use SSE, if available
|
||||
// by precalculating these values we can greatly reduce the number of calculations in idirpyr_eq_channel()
|
||||
// but we need two additional buffers for this preprocessing
|
||||
tmpHue = new float*[srcheight];
|
||||
|
||||
for (int i = 0; i < srcheight; i++) {
|
||||
tmpHue[i] = new float[srcwidth];
|
||||
}
|
||||
|
||||
#ifdef __SSE2__
|
||||
#pragma omp parallel for
|
||||
|
||||
for(int i = 0; i < srcheight; i++) {
|
||||
int j;
|
||||
|
||||
for(j = 0; j < srcwidth - 3; j += 4) {
|
||||
_mm_storeu_ps(&tmpHue[i][j], xatan2f(LVFU(l_b[i][j]), LVFU(l_a[i][j])));
|
||||
}
|
||||
|
||||
for(; j < srcwidth; j++) {
|
||||
tmpHue[i][j] = xatan2f(l_b[i][j], l_a[i][j]);
|
||||
}
|
||||
}
|
||||
|
||||
#else
|
||||
#pragma omp parallel for
|
||||
|
||||
for(int i = 0; i < srcheight; i++) {
|
||||
for(int j = 0; j < srcwidth; j++) {
|
||||
tmpHue[i][j] = xatan2f(l_b[i][j], l_a[i][j]);
|
||||
}
|
||||
}
|
||||
|
||||
#endif
|
||||
tmpChr = new float*[srcheight];
|
||||
|
||||
for (int i = 0; i < srcheight; i++) {
|
||||
tmpChr[i] = new float[srcwidth];
|
||||
}
|
||||
|
||||
#ifdef __SSE2__
|
||||
#pragma omp parallel
|
||||
{
|
||||
__m128 div = _mm_set1_ps(327.68f);
|
||||
#pragma omp for
|
||||
|
||||
for(int i = 0; i < srcheight; i++) {
|
||||
int j;
|
||||
|
||||
for(j = 0; j < srcwidth - 3; j += 4) {
|
||||
_mm_storeu_ps(&tmpChr[i][j], _mm_sqrt_ps(SQRV(LVFU(l_b[i][j])) + SQRV(LVFU(l_a[i][j]))) / div);
|
||||
}
|
||||
|
||||
for(; j < srcwidth; j++) {
|
||||
tmpChr[i][j] = sqrtf(SQR((l_b[i][j])) + SQR((l_a[i][j]))) / 327.68f;
|
||||
}
|
||||
}
|
||||
}
|
||||
#else
|
||||
#pragma omp parallel for
|
||||
|
||||
for(int i = 0; i < srcheight; i++) {
|
||||
for(int j = 0; j < srcwidth; j++) {
|
||||
tmpChr[i][j] = sqrtf(SQR((l_b[i][j])) + SQR((l_a[i][j]))) / 327.68f;
|
||||
}
|
||||
}
|
||||
|
||||
#endif
|
||||
}
|
||||
*/
|
||||
// with the current implementation of idirpyr_eq_channel we can safely use the buffer from last level as buffer, saves some memory
|
||||
float ** buffer = dirpyrlo[lastlevel - 1];
|
||||
|
||||
for(int level = lastlevel - 1; level > 0; level--) {
|
||||
idirpyr_eq_channel_loc(dirpyrlo[level], dirpyrlo[level - 1], loctemp, buffer, srcwidth, srcheight, level, multi, dirpyrThreshold, tmpHue, tmpChr, skinprot, gamutlab, b_l, t_l, t_r, b_r, choice );
|
||||
}
|
||||
|
||||
scale = scalesloc[0];
|
||||
|
||||
idirpyr_eq_channel_loc(dirpyrlo[0], dst, loctemp, buffer, srcwidth, srcheight, 0, multi, dirpyrThreshold, tmpHue, tmpChr, skinprot, gamutlab, b_l, t_l, t_r, b_r, choice );
|
||||
/*
|
||||
if(skinprot != 0.f) {
|
||||
for (int i = 0; i < srcheight; i++) {
|
||||
delete [] tmpChr[i];
|
||||
}
|
||||
|
||||
delete [] tmpChr;
|
||||
|
||||
for (int i = 0; i < srcheight; i++) {
|
||||
delete [] tmpHue[i];
|
||||
}
|
||||
|
||||
delete [] tmpHue;
|
||||
}
|
||||
*/
|
||||
//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
||||
#pragma omp parallel for
|
||||
|
||||
for (int i = 0; i < srcheight; i++)
|
||||
for (int j = 0; j < srcwidth; j++) {
|
||||
dst[i][j] = src[i][j];
|
||||
loctemp[i][j] = CLIP(buffer[i][j]); // TODO: Really a clip necessary?
|
||||
}
|
||||
|
||||
}
|
||||
|
||||
|
||||
void ImProcFunctions :: dirpyr_equalizercam (CieImage *ncie, float ** src, float ** dst, int srcwidth, int srcheight, float ** h_p, float ** C_p, const double * mult, const double dirpyrThreshold, const double skinprot, bool execdir, const bool gamutlab, float b_l, float t_l, float t_r, float b_r, int choice, int scaleprev)
|
||||
@@ -624,6 +834,102 @@ SSEFUNCTION void ImProcFunctions::dirpyr_channel(float ** data_fine, float ** da
|
||||
}
|
||||
|
||||
//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
||||
void ImProcFunctions::idirpyr_eq_channel_loc(float ** data_coarse, float ** data_fine, float ** loctemp, float ** buffer, int width, int height, int level, float mult[5], const double dirpyrThreshold, float ** hue, float ** chrom, const double skinprot, const bool gamutlab, float b_l, float t_l, float t_r, float b_r , int choice)
|
||||
{
|
||||
const float skinprotneg = -skinprot;
|
||||
const float factorHard = (1.f - skinprotneg / 100.f);
|
||||
|
||||
float offs;
|
||||
|
||||
if(skinprot == 0.f) {
|
||||
offs = 0.f;
|
||||
} else {
|
||||
offs = -1.f;
|
||||
}
|
||||
|
||||
float multbis[5];
|
||||
|
||||
multbis[level] = mult[level]; //multbis to reduce artifacts for high values mult
|
||||
|
||||
if(level == 4 && mult[level] > 1.f) {
|
||||
multbis[level] = 1.f + 0.65f * (mult[level] - 1.f);
|
||||
}
|
||||
|
||||
// if(level == 5 && mult[level] > 1.f) {
|
||||
// multbis[level] = 1.f + 0.45f * (mult[level] - 1.f);
|
||||
// }
|
||||
|
||||
LUTf irangefn (0x20000);
|
||||
{
|
||||
const float noisehi = 1.33f * noise * dirpyrThreshold / expf(level * log(3.0)), noiselo = 0.66f * noise * dirpyrThreshold / expf(level * log(3.0));
|
||||
//printf("level=%i multlev=%f noisehi=%f noiselo=%f skinprot=%f\n",level,mult[level], noisehi, noiselo, skinprot);
|
||||
|
||||
for (int i = 0; i < 0x20000; i++) {
|
||||
if (abs(i - 0x10000) > noisehi || multbis[level] < 1.0) {
|
||||
irangefn[i] = multbis[level] + offs;
|
||||
} else {
|
||||
if (abs(i - 0x10000) < noiselo) {
|
||||
irangefn[i] = 1.f + offs ;
|
||||
} else {
|
||||
irangefn[i] = 1.f + offs + (multbis[level] - 1.f) * (noisehi - abs(i - 0x10000)) / (noisehi - noiselo + 0.01f) ;
|
||||
}
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
if(skinprot == 0.f)
|
||||
#ifdef _OPENMP
|
||||
#pragma omp parallel for schedule(dynamic,16)
|
||||
#endif
|
||||
for(int i = 0; i < height; i++) {
|
||||
for(int j = 0; j < width; j++) {
|
||||
float hipass = (data_fine[i][j] - data_coarse[i][j]);
|
||||
buffer[i][j] += irangefn[hipass + 0x10000] * hipass;
|
||||
}
|
||||
}
|
||||
|
||||
/*
|
||||
else if(skinprot > 0.f)
|
||||
#ifdef _OPENMP
|
||||
#pragma omp parallel for schedule(dynamic,16)
|
||||
#endif
|
||||
for(int i = 0; i < height; i++) {
|
||||
for(int j = 0; j < width; j++) {
|
||||
float scale = 1.f;
|
||||
float hipass = (data_fine[i][j] - data_coarse[i][j]);
|
||||
// These values are precalculated now
|
||||
float modhue = hue[i][j];
|
||||
float modchro = chrom[i][j];
|
||||
Color::SkinSatCbdl ((data_fine[i][j]) / 327.68f, modhue, modchro, skinprot, scale, true, b_l, t_l, t_r);
|
||||
buffer[i][j] += (1.f + (irangefn[hipass + 0x10000]) * scale) * hipass ;
|
||||
}
|
||||
}
|
||||
else
|
||||
#ifdef _OPENMP
|
||||
#pragma omp parallel for schedule(dynamic,16)
|
||||
#endif
|
||||
for(int i = 0; i < height; i++) {
|
||||
for(int j = 0; j < width; j++) {
|
||||
float scale = 1.f;
|
||||
float hipass = (data_fine[i][j] - data_coarse[i][j]);
|
||||
// These values are precalculated now
|
||||
float modhue = hue[i][j];
|
||||
float modchro = chrom[i][j];
|
||||
Color::SkinSatCbdl ((data_fine[i][j]) / 327.68f, modhue, modchro, skinprotneg, scale, false, b_l, t_l, t_r);
|
||||
float correct = irangefn[hipass + 0x10000];
|
||||
|
||||
if (scale == 1.f) {//image hard
|
||||
buffer[i][j] += (1.f + (correct) * (factorHard)) * hipass ;
|
||||
} else { //image soft with scale < 1 ==> skin
|
||||
buffer[i][j] += (1.f + (correct)) * hipass ;
|
||||
}
|
||||
}
|
||||
}
|
||||
*/
|
||||
}
|
||||
|
||||
|
||||
|
||||
|
||||
void ImProcFunctions::idirpyr_eq_channel(float ** data_coarse, float ** data_fine, float ** buffer, int width, int height, int level, float mult[5], const double dirpyrThreshold, float ** hue, float ** chrom, const double skinprot, const bool gamutlab, float b_l, float t_l, float t_r, float b_r , int choice)
|
||||
{
|
||||
|
||||
Reference in New Issue
Block a user