123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727728729730731732733734735736737738739740741742743744745746747748749750751752753754755756757758759760761762763764765766767768769770771772773774775776777778779780781782783784785786787788789790791792793794795 |
- // This file is part of HDRip.
- //
- // HDRip 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.
- //
- // HDRip 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 HDRip. If not, see <https://www.gnu.org/licenses/>.
- //
- // HDRip project
- // Author : Rémi Synave
- // Contact : remi.synave@univ-littoral.fr
- #include "pch.h"
- #include <iostream>
- #include <thread>
- #include "Conversion.hpp"
- #include "MT_linear.hpp"
- #include "MT_channel.hpp"
- /**************************************/
- /******** LINEAR_TO_NON_LINEAR ********/
- /**************************************/
- float Conversion::linear_to_non_linear(float data)
- {
- if (data <= 0.0031308f)
- return data * 12.92f;
- return 1.055f * powf(data, 0.4166666667f) - 0.055f;
- }
- #ifdef _MT_
- void* linear_to_non_linear_MT(void* arg)
- {
- MT_linear* a = (MT_linear*)arg;
- const float* data = a->data;
- float* result = a->result;
- for (unsigned int i = 0; i < a->length; i++)
- result[i] = Conversion::linear_to_non_linear(data[i]);
- return arg;
- }
- float* Conversion::linear_to_non_linear(const float* data, unsigned int length)
- {
- float* non_linear = new float[length];
- std::thread tab_t[_MT_];
- MT_linear tab_a[_MT_];
- unsigned int id;
- unsigned int block_size = length / _MT_;
- for (id = 0; id < _MT_; id++) {
- tab_a[id].data = data + (id * block_size);
- tab_a[id].length = block_size;
- tab_a[id].result = non_linear + (id * block_size);
- if (id == (_MT_ - 1))
- tab_a[id].length = length - ((_MT_ - 1) * block_size);
- tab_t[id] = std::thread(linear_to_non_linear_MT, (void*)(tab_a + id));
- }
- for (id = 0; id < _MT_; id++) {
- tab_t[id].join();
- }
- return non_linear;
- }
- #else
- float* Conversion::linear_to_non_linear(const float* data, unsigned int length)
- {
- float* non_linear = new float[length];
- for (unsigned int i = 0; i < length; i++)
- non_linear[i] = linear_to_non_linear(data[i]);
- return non_linear;
- }
- #endif
- /**************************************/
- /******** NON_LINEAR_TO_LINEAR ********/
- /**************************************/
- float Conversion::non_linear_to_linear(float data)
- {
- if (data <= 0.040449936f)
- return data / 12.92f;
- return powf((data + 0.055f) / 1.055f, 2.4f);
- }
- #ifdef _MT_
- void* non_linear_to_linear_MT(void* arg)
- {
- MT_linear* a = (MT_linear*)arg;
- const float* data = a->data;
- float* result = a->result;
- for (unsigned int i = 0; i < a->length; i++)
- result[i] = Conversion::non_linear_to_linear(data[i]);
- return arg;
- }
- float* Conversion::non_linear_to_linear(const float* data, unsigned int length)
- {
- float* linear = new float[length];
- std::thread tab_t[_MT_];
- MT_linear tab_a[_MT_];
- unsigned int id;
- unsigned int block_size = length / _MT_;
- for (id = 0; id < _MT_; id++) {
- tab_a[id].data = data + (id * block_size);
- tab_a[id].length = block_size;
- tab_a[id].result = linear + (id * block_size);
- if (id == (_MT_ - 1))
- tab_a[id].length = length - ((_MT_ - 1) * block_size);
- tab_t[id] = std::thread(non_linear_to_linear_MT, (void*)(tab_a + id));
- }
- for (id = 0; id < _MT_; id++) {
- tab_t[id].join();
- }
- return linear;
- }
- #else
- float* Conversion::non_linear_to_linear(const float* data, unsigned int length)
- {
- float* linear = new float[length];
- for (unsigned int i = 0; i < length; i++)
- linear[i] = non_linear_to_linear(data[i]);
- return linear;
- }
- #endif
- /*************************************/
- /************ sRGB_TO_XYZ ************/
- /*************************************/
- std::tuple<float, float, float> Conversion::sRGB_to_XYZ(float r, float g, float b)
- {
- float x = r * Conversion::sRGB_to_XYZ_m[0][0] + g * Conversion::sRGB_to_XYZ_m[0][1] + b * Conversion::sRGB_to_XYZ_m[0][2];
- float y = r * Conversion::sRGB_to_XYZ_m[1][0] + g * Conversion::sRGB_to_XYZ_m[1][1] + b * Conversion::sRGB_to_XYZ_m[1][2];
- float z = r * Conversion::sRGB_to_XYZ_m[2][0] + g * Conversion::sRGB_to_XYZ_m[2][1] + b * Conversion::sRGB_to_XYZ_m[2][2];
- return std::make_tuple(x, y, z);
- }
- #ifdef _MT_
- void* sRGB_to_XYZ_MT(void* arg)
- {
- MT_channel* a = (MT_channel*)arg;
- const float* data = a->data;
- float* result = a->channel;
- for (unsigned int i = 0; i < a->length; i++)
- {
- std::tuple<float, float, float> v = Conversion::sRGB_to_XYZ(data[i * 3], data[i * 3 + 1], data[i * 3 + 2]);
- result[i * 3] = std::get<0>(v);
- result[i * 3 + 1] = std::get<1>(v);
- result[i * 3 + 2] = std::get<2>(v);
- }
- return arg;
- }
- float* Conversion::sRGB_to_XYZ(const float* data, const unsigned int length)
- {
- float* xyz = new float[length * 3];
- std::thread tab_t[_MT_];
- MT_channel tab_a[_MT_];
- unsigned int id;
- unsigned int block_size = length / _MT_;
- for (id = 0; id < _MT_; id++) {
- tab_a[id].data = data + (id * block_size * 3);
- tab_a[id].length = block_size;
- tab_a[id].channel = xyz + (id * block_size * 3);
- if (id == (_MT_ - 1))
- tab_a[id].length = length - ((_MT_ - 1) * block_size);
- tab_t[id] = std::thread(sRGB_to_XYZ_MT, (void*)(tab_a + id));
- }
- for (id = 0; id < _MT_; id++) {
- tab_t[id].join();
- }
- return xyz;
- }
- #else
- float* Conversion::sRGB_to_XYZ(const float* data, const unsigned int length)
- {
- float* xyz = new float[length * 3];
- for (unsigned int i = 0; i < length; i++)
- {
- std::tuple<float, float, float> conv = sRGB_to_XYZ(data[i * 3], data[i * 3 + 1], data[i * 3 + 2]);
- xyz[i * 3] = std::get<0>(conv);
- xyz[i * 3 + 1] = std::get<1>(conv);
- xyz[i * 3 + 2] = std::get<2>(conv);
- }
- return xyz;
- }
- #endif
- float Conversion::sRGB_to_Y_of_XYZ(float r, float g, float b)
- {
- return (r * Conversion::sRGB_to_XYZ_m[1][0] + g * Conversion::sRGB_to_XYZ_m[1][1] + b * Conversion::sRGB_to_XYZ_m[1][2]);
- }
- #ifdef _MT_
- void* sRGB_to_Y_of_XYZ_MT(void* arg)
- {
- MT_channel* a = (MT_channel*)arg;
- const float* data = a->data;
- float* channel = a->channel;
- for (unsigned int i = 0; i < a->length; i++)
- channel[i] = Conversion::sRGB_to_Y_of_XYZ(data[i * 3], data[i * 3 + 1], data[i * 3 + 2]);
- return arg;
- }
- float* Conversion::sRGB_to_Y_of_XYZ(const float* data, const unsigned int length)
- {
- float* channelY = new float[length];
- std::thread tab_t[_MT_];
- MT_channel tab_a[_MT_];
- unsigned int id;
- unsigned int block_size = length / _MT_;
- for (id = 0; id < _MT_; id++) {
- tab_a[id].data = data + (id * block_size * 3);
- tab_a[id].length = block_size;
- tab_a[id].channel = channelY + (id * block_size);
- if (id == (_MT_ - 1))
- tab_a[id].length = length - ((_MT_ - 1) * block_size);
- tab_t[id] = std::thread(sRGB_to_Y_of_XYZ_MT, (void*)(tab_a + id));
- }
- for (id = 0; id < _MT_; id++) {
- tab_t[id].join();
- }
- return channelY;
- }
- #else
- float* Conversion::sRGB_to_Y_of_XYZ(const float* data, const unsigned int length)
- {
- float* y = new float[length];
- for (unsigned int i = 0; i < length; i++)
- y[i] = sRGB_to_Y_of_XYZ(data[i * 3], data[i * 3 + 1], data[i * 3 + 2]);
- return y;
- }
- #endif
- /*************************************/
- /************* XYZ_TO_LAB ************/
- /*************************************/
- std::tuple<float, float, float> Conversion::XYZ_to_Lab(float x, float y, float z)
- {
- float xNorm = x / 0.950455927f;
- float yNorm = y;
- float zNorm = z / 1.08905775f;
- float coeff = 16.0f / 116.0f;
- float fx = 7.787f * xNorm + coeff;
- float fy = 7.787f * yNorm + coeff;
- float fz = 7.787f * zNorm + coeff;
- if (xNorm > 0.008856f)
- fx = powf(xNorm, 0.3333333333f);
- if (yNorm > 0.008856f)
- fy = powf(yNorm, 0.3333333333f);
- if (zNorm > 0.008856f)
- fz = powf(zNorm, 0.3333333333f);
- return std::make_tuple(116.0f * fy - 16.0f, 500.0f * (fx - fy), 200.0f * (fy - fz));
- }
- #ifdef _MT_
- void* XYZ_to_Lab_MT(void* arg)
- {
- MT_channel* a = (MT_channel*)arg;
- const float* data = a->data;
- float* result = a->channel;
- for (unsigned int i = 0; i < a->length; i++)
- {
- std::tuple<float, float, float> v = Conversion::XYZ_to_Lab(data[i * 3], data[i * 3 + 1], data[i * 3 + 2]);
- result[i * 3] = std::get<0>(v);
- result[i * 3 + 1] = std::get<1>(v);
- result[i * 3 + 2] = std::get<2>(v);
- }
- return arg;
- }
- float* Conversion::XYZ_to_Lab(const float* data, const unsigned int length)
- {
- float* channelLab = new float[length * 3];
- std::thread tab_t[_MT_];
- MT_channel tab_a[_MT_];
- unsigned int id;
- unsigned int block_size = length / _MT_;
- for (id = 0; id < _MT_; id++) {
- tab_a[id].data = data + (id * block_size * 3);
- tab_a[id].length = block_size;
- tab_a[id].channel = channelLab + (id * block_size * 3);
- if (id == (_MT_ - 1))
- tab_a[id].length = length - ((_MT_ - 1) * block_size);
- tab_t[id] = std::thread(XYZ_to_Lab_MT, (void*)(tab_a + id));
- }
- for (id = 0; id < _MT_; id++) {
- tab_t[id].join();
- }
- return channelLab;
- }
- #else
- float* Conversion::XYZ_to_Lab(const float* data, const unsigned int length)
- {
- float* channelLab = new float[length*3];
- for (unsigned int i = 0; i < length; i++)
- {
- std::tuple<float, float, float> lab = XYZ_to_Lab(data[i * 3], data[i * 3 + 1], data[i * 3 + 2]);
- channelLab[i*3] = std::get<0>(lab);
- channelLab[i*3+1] = std::get<1>(lab);
- channelLab[i*3+2] = std::get<2>(lab);
- }
- return channelLab;
- }
- #endif
- /*************************************/
- /************ SRGB_TO_LAB ************/
- /*************************************/
- std::tuple<float, float, float> Conversion::sRGB_to_Lab(float r, float g, float b)
- {
- std::tuple<float, float, float> xyz = sRGB_to_XYZ(r, g, b);
- return XYZ_to_Lab(std::get<0>(xyz),std::get<1>(xyz),std::get<2>(xyz));
- }
- float* Conversion::sRGB_to_Lab(const float* data, const unsigned int length)
- {
- float* rgb_to_xyz = sRGB_to_XYZ(data, length);
- float* lab = XYZ_to_Lab(rgb_to_xyz,length);
- delete[](rgb_to_xyz);
- return lab;
- }
- float Conversion::sRGB_to_L_of_Lab(float r, float g, float b)
- {
- std::tuple<float, float, float> xyz = sRGB_to_XYZ(r, g, b);
- float fy = 7.787f * std::get<1>(xyz) + (16.0f / 116.0f);
- if (std::get<1>(xyz) > 0.008856f)
- fy = powf(std::get<1>(xyz), 1.0f / 3.0f);
- return (116.0f * fy - 16.0f);
- }
- /************************************/
- /************ LAB_TO_LCH ************/
- /************************************/
- std::tuple<float, float, float> Conversion::Lab_to_LCH(float L, float a, float b)
- {
- float C = sqrtf(a * a + b * b);
- float theta = atan2(b, a);
- while (theta < 0)
- theta += (float)(2.0f * M_PI);
- while (theta > (float)(2.0f * M_PI))
- theta -= (float)(2.0f * M_PI);
- float H = theta / ((float)M_PI) * 180.0f;
- return std::make_tuple(L, C, H);
- }
- #ifdef _MT_
- void* Lab_to_LCH_MT(void* arg)
- {
- MT_channel* a = (MT_channel*)arg;
- const float* data = a->data;
- float* result = a->channel;
- for (unsigned int i = 0; i < a->length; i++)
- {
- std::tuple<float, float, float> v = Conversion::Lab_to_LCH(data[i * 3], data[i * 3 + 1], data[i * 3 + 2]);
- result[i * 3] = std::get<0>(v);
- result[i * 3 + 1] = std::get<1>(v);
- result[i * 3 + 2] = std::get<2>(v);
- }
- return arg;
- }
- float* Conversion::Lab_to_LCH(const float* data, const unsigned int length)
- {
- float* channelLCH = new float[length * 3];
- std::thread tab_t[_MT_];
- MT_channel tab_a[_MT_];
- unsigned int id;
- unsigned int block_size = length / _MT_;
- for (id = 0; id < _MT_; id++) {
- tab_a[id].data = data + (id * block_size * 3);
- tab_a[id].length = block_size;
- tab_a[id].channel = channelLCH + (id * block_size * 3);
- if (id == (_MT_ - 1))
- tab_a[id].length = length - ((_MT_ - 1) * block_size);
- tab_t[id] = std::thread(Lab_to_LCH_MT, (void*)(tab_a + id));
- }
- for (id = 0; id < _MT_; id++) {
- tab_t[id].join();
- }
- return channelLCH;
- }
- #else
- float* Conversion::Lab_to_LCH(const float* data, const unsigned int length)
- {
- float* LCH = new float[length * 3];
- for (unsigned int i = 0; i < length; i++)
- {
- std::tuple<float, float, float> conv = Lab_to_LCH(data[i * 3], data[i * 3 + 1], data[i * 3 + 2]);
- LCH[i * 3] = std::get<0>(conv);
- LCH[i * 3 + 1] = std::get<1>(conv);
- LCH[i * 3 + 2] = std::get<2>(conv);
- }
- return LCH;
- }
- #endif
- float Conversion::Lab_to_C_of_LCH(float a, float b)
- {
- return sqrtf(a * a + b * b);
- }
- float Conversion::Lab_to_H_of_LCH(float a, float b)
- {
- float theta = atan2(b, a);
- while (theta < 0)
- theta += (float)(2.0f * M_PI);
- while (theta > (float)(2.0f * M_PI))
- theta -= (float)(2.0f * M_PI);
- return theta / ((float)M_PI) * 180.0f;
- }
- /************************************/
- /************ LCH_TO_LAB ************/
- /************************************/
- std::tuple<float, float, float> Conversion::LCH_to_Lab(float L, float C, float H)
- {
- float rho = C;
- float phi = (H/180.0f)*((float)M_PI);
- float a = rho*cos(phi);
- float b = rho*sin(phi);
- return std::make_tuple(L,a,b);
- }
- #ifdef _MT_
- void* LCH_to_Lab_MT(void* arg)
- {
- MT_channel* a = (MT_channel*)arg;
- const float* data = a->data;
- float* result = a->channel;
- for (unsigned int i = 0; i < a->length; i++)
- {
- std::tuple<float, float, float> v = Conversion::LCH_to_Lab(data[i * 3], data[i * 3 + 1], data[i * 3 + 2]);
- result[i * 3] = std::get<0>(v);
- result[i * 3 + 1] = std::get<1>(v);
- result[i * 3 + 2] = std::get<2>(v);
- }
- return arg;
- }
- float* Conversion::LCH_to_Lab(const float* data, const unsigned int length)
- {
- float* channelLab = new float[length * 3];
- std::thread tab_t[_MT_];
- MT_channel tab_a[_MT_];
- unsigned int id;
- unsigned int block_size = length / _MT_;
- for (id = 0; id < _MT_; id++) {
- tab_a[id].data = data + (id * block_size * 3);
- tab_a[id].length = block_size;
- tab_a[id].channel = channelLab + (id * block_size * 3);
- if (id == (_MT_ - 1))
- tab_a[id].length = length - ((_MT_ - 1) * block_size);
- tab_t[id] = std::thread(LCH_to_Lab_MT, (void*)(tab_a + id));
- }
- for (id = 0; id < _MT_; id++) {
- tab_t[id].join();
- }
- return channelLab;
- }
- #else
- float* Conversion::LCH_to_Lab(const float* data, const unsigned int length)
- {
- float* Lab = new float[length * 3];
- for (unsigned int i = 0; i < length; i++)
- {
- std::tuple<float, float, float> conv = LCH_to_Lab(data[i * 3], data[i * 3 + 1], data[i * 3 + 2]);
- Lab[i * 3] = std::get<0>(conv);
- Lab[i * 3 + 1] = std::get<1>(conv);
- Lab[i * 3 + 2] = std::get<2>(conv);
- }
- return Lab;
- }
- #endif
- /************************************/
- /************ LAB_TO_XYZ ************/
- /************************************/
- std::tuple<float, float, float> Conversion::Lab_to_XYZ(float L, float a, float b)
- {
- float fy = (L+16.0f)/116.0f;
- float fx = a/500.0f+fy;
- float fz = fy - b/200.0f;
- float xNorm = 0.950455927f*(fx - 0.137931034f)*0.128418549f;
- float yNorm = (fy - 0.137931034f)*0.128418549f;
- float zNorm = 1.08905775f*(fz - 0.137931034f)*0.128418549f;
- if(fx>0.206896552f)
- xNorm = 0.950455927f*fx*fx*fx;
- if(fy>0.206896552f)
- yNorm = fy*fy*fy;
- if(fz>0.206896552f)
- zNorm = 1.08905775f*fz*fz*fz;
-
- return std::make_tuple(xNorm,yNorm,zNorm);
- }
- #ifdef _MT_
- void* Lab_to_XYZ_MT(void* arg)
- {
- MT_channel* a = (MT_channel*)arg;
- const float* data = a->data;
- float* result = a->channel;
- for (unsigned int i = 0; i < a->length; i++)
- {
- std::tuple<float, float, float> v = Conversion::Lab_to_XYZ(data[i * 3], data[i * 3 + 1], data[i * 3 + 2]);
- result[i * 3] = std::get<0>(v);
- result[i * 3 + 1] = std::get<1>(v);
- result[i * 3 + 2] = std::get<2>(v);
- }
- return arg;
- }
- float* Conversion::Lab_to_XYZ(const float* data, const unsigned int length)
- {
- float* channelXYZ = new float[length * 3];
- std::thread tab_t[_MT_];
- MT_channel tab_a[_MT_];
- unsigned int id;
- unsigned int block_size = length / _MT_;
- for (id = 0; id < _MT_; id++) {
- tab_a[id].data = data + (id * block_size * 3);
- tab_a[id].length = block_size;
- tab_a[id].channel = channelXYZ + (id * block_size * 3);
- if (id == (_MT_ - 1))
- tab_a[id].length = length - ((_MT_ - 1) * block_size);
- tab_t[id] = std::thread(Lab_to_XYZ_MT, (void*)(tab_a + id));
- }
- for (id = 0; id < _MT_; id++) {
- tab_t[id].join();
- }
- return channelXYZ;
- }
- #else
- float* Conversion::Lab_to_XYZ(const float* data, const unsigned int length)
- {
- float* Lab = new float[length * 3];
- for (unsigned int i = 0; i < length; i++)
- {
- std::tuple<float, float, float> conv = Lab_to_XYZ(data[i * 3], data[i * 3 + 1], data[i * 3 + 2]);
- Lab[i * 3] = std::get<0>(conv);
- Lab[i * 3 + 1] = std::get<1>(conv);
- Lab[i * 3 + 2] = std::get<2>(conv);
- }
- return Lab;
- }
- #endif
- /*************************************/
- /************ XYZ_TO_sRGB ************/
- /*************************************/
- std::tuple<float, float, float> Conversion::XYZ_to_sRGB(float x, float y, float z)
- {
- float r = x * Conversion::XYZ_to_sRGB_m[0][0] + y * Conversion::XYZ_to_sRGB_m[0][1] + z * Conversion::XYZ_to_sRGB_m[0][2];
- float g = x * Conversion::XYZ_to_sRGB_m[1][0] + y * Conversion::XYZ_to_sRGB_m[1][1] + z * Conversion::XYZ_to_sRGB_m[1][2];
- float b = x * Conversion::XYZ_to_sRGB_m[2][0] + y * Conversion::XYZ_to_sRGB_m[2][1] + z * Conversion::XYZ_to_sRGB_m[2][2];
- return std::make_tuple(r, g, b);
- }
- #ifdef _MT_
- void* XYZ_to_sRGB_MT(void* arg)
- {
- MT_channel* a = (MT_channel*)arg;
- const float* data = a->data;
- float* result = a->channel;
- for (unsigned int i = 0; i < a->length; i++)
- {
- std::tuple<float, float, float> v = Conversion::XYZ_to_sRGB(data[i * 3], data[i * 3 + 1], data[i * 3 + 2]);
- result[i * 3] = std::get<0>(v);
- result[i * 3 + 1] = std::get<1>(v);
- result[i * 3 + 2] = std::get<2>(v);
- }
- return arg;
- }
- float* Conversion::XYZ_to_sRGB(const float* data, const unsigned int length)
- {
- float* channelRGB = new float[length * 3];
- std::thread tab_t[_MT_];
- MT_channel tab_a[_MT_];
- unsigned int id;
- unsigned int block_size = length / _MT_;
- for (id = 0; id < _MT_; id++) {
- tab_a[id].data = data + (id * block_size * 3);
- tab_a[id].length = block_size;
- tab_a[id].channel = channelRGB + (id * block_size * 3);
- if (id == (_MT_ - 1))
- tab_a[id].length = length - ((_MT_ - 1) * block_size);
- tab_t[id] = std::thread(XYZ_to_sRGB_MT, (void*)(tab_a + id));
- }
- for (id = 0; id < _MT_; id++) {
- tab_t[id].join();
- }
- return channelRGB;
- }
- #else
- float* Conversion::XYZ_to_sRGB(const float* data, const unsigned int length)
- {
- float* rgb = new float[length * 3];
- for (unsigned int i = 0; i < length; i++)
- {
- std::tuple<float, float, float> conv = XYZ_to_sRGB(data[i * 3], data[i * 3 + 1], data[i * 3 + 2]);
- rgb[i * 3] = std::get<0>(conv);
- rgb[i * 3 + 1] = std::get<1>(conv);
- rgb[i * 3 + 2] = std::get<2>(conv);
- }
- return rgb;
- }
- #endif
- /*************************************/
- /************ LCH_TO_sRGB ************/
- /*************************************/
- float* Conversion::LCH_to_sRGB(const float* data, const unsigned int length)
- {
- float* Lab = LCH_to_Lab(data,length);
- float* xyz = Lab_to_XYZ(Lab,length);
- float* rgb= XYZ_to_sRGB(xyz,length);
- delete[](Lab);
- delete[](xyz);
- return rgb;
- }
|