summaryrefslogtreecommitdiff
path: root/sci_gateway/cpp/opencv_lab2rgb.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'sci_gateway/cpp/opencv_lab2rgb.cpp')
-rw-r--r--sci_gateway/cpp/opencv_lab2rgb.cpp331
1 files changed, 331 insertions, 0 deletions
diff --git a/sci_gateway/cpp/opencv_lab2rgb.cpp b/sci_gateway/cpp/opencv_lab2rgb.cpp
new file mode 100644
index 0000000..5de6b2d
--- /dev/null
+++ b/sci_gateway/cpp/opencv_lab2rgb.cpp
@@ -0,0 +1,331 @@
+/* ==================================================================== */
+/* Author :Tess Zacharias */
+/* ==================================================================== */
+/* overloaded function */
+/* Syntax : return_matrix_3_xyz_values=lab2xyz(matrix_of_3_labvalues a) */
+/* Syntax : return_matrix_3_xyz_values=lab2xyz(matrix_of_3_labvalues a,String1 Name,String2 value) */
+/* ==================================================================== */
+/* ==================================================================== */
+#include <numeric>
+#include "opencv2/core/core.hpp"
+#include "opencv2/highgui/highgui.hpp"
+#include "opencv2/opencv.hpp"
+#include <iostream>
+#include<stdio.h>
+#include <sciprint.h>
+#include<string.h>
+using namespace cv;
+using namespace std;
+ extern "C"
+ {
+ #include "api_scilab.h"
+ #include "Scierror.h"
+ #include "BOOL.h"
+ #include <localization.h>
+ // #include "../common.h"
+ #define REF_X .95047; // Observer= 2°, Illuminant= D65
+#define REF_Y 1.00000;
+#define REF_Z 1.08883;
+
+struct Color
+{
+ float R,G,B,X,Y,Z,L,a,b;
+};
+
+Color lab2rgb1(float L, float a, float b)
+{
+
+
+ float y = (L + 16) / 116;
+ float x = a / 500 + y;
+ float z = y - b / 200;
+ if ( pow( y , 3 ) > 0.008856 )
+ {
+ y = pow( y , 3 );
+ }
+ else
+ {
+ y = ( y - 16 / 116 ) / 7.787;
+ }
+ if ( pow( x , 3 ) > 0.008856 )
+ {
+ x = pow( x , 3 );
+ }
+ else
+ {
+ x = ( x - 16 / 116 ) / 7.787;
+ }
+ if ( pow( z , 3 ) > 0.008856 )
+ {
+ z = pow( z , 3 );
+ }
+ else
+ {
+ z = ( z - 16 / 116 ) / 7.787;
+ }
+ float X= x *REF_X;
+ float Y= y *REF_Y;
+ float Z = z *REF_Z;
+ //X = X/100.f; //R 0..1
+ //Y= Y/100.f; //G 0..1
+ //Z= Z/100.f; //B 0..1
+ float R1=3.240479*X-1.537150*Y -0.498535 *Z;
+ float G1=-0.969256*X+1.875992 *Y+0.041556* Z;
+ float B1=0.055648*X-0.204043*Y +1.057311*Z;
+ if( R1 > 0.0031308 )
+ R1 = 1.055 * ( pow( R1, 1./2.4 ) ) - 0.055;
+ else
+ R1 = 12.92 * R1;
+ if( G1 > 0.0031308 )
+ G1 = 1.055 * ( pow( G1, 1./2.4 ) ) - 0.055;
+ else
+ G1 = 12.92 * G1;
+ if( B1 > 0.0031308 )
+ B1 = 1.055 * ( pow( B1, 1./2.4 ) ) - 0.055;
+ else
+ B1 = 12.92 * B1;
+
+
+ Color rgb;
+ rgb.R= R1;
+ rgb.G= G1;
+ rgb.B= B1;
+ return rgb;
+
+}
+ Color lab2rgb2(float L, float a, float b,char *s)
+{
+ double REF_X1;
+ double REF_Y1;
+ double REF_Z1;
+ if(!strcmp(s,"a"))
+ {
+ REF_X1=1.0985;
+ REF_Y1= 1.0000;
+ REF_Z1=0.3558;
+ }
+ else if(!strcmp(s,"c"))
+ {
+ REF_X1=0.9807;
+ REF_Y1= 1.0000;
+ REF_Z1=1.1822;
+ }
+ else if(!strcmp(s,"d50"))
+ {
+ REF_X1=0.9642;
+ REF_Y1= 1.0000;
+ REF_Z1=0.8251;
+ }
+ else if(!strcmp(s,"d65"))
+ {
+ REF_X1=0.9504;
+ REF_Y1= 1.0000;
+ REF_Z1= 1.0888;
+ }
+ else if(!strcmp(s,"icc"))
+ {
+ REF_X1=0.962;
+ REF_Y1=1.000;
+ REF_Z1= 0.8249;
+ }
+ else if(!strcmp(s,"e"))
+ {
+ REF_X1=1.000;
+ REF_Y1=1.000;
+ REF_Z1=1.000;
+ }
+
+ float y = (L + 16) / 116;
+ float x = a / 500 + y;
+ float z = y - b / 200;
+ if ( pow( y , 3 ) > 0.008856 )
+ {
+ y = pow( y , 3 );
+ }
+ else
+ {
+ y = ( y - 16 / 116 ) / 7.787;
+ }
+ if ( pow( x , 3 ) > 0.008856 )
+ {
+ x = pow( x , 3 );
+ }
+ else
+ {
+ x = ( x - 16 / 116 ) / 7.787;
+ }
+ if ( pow( z , 3 ) > 0.008856 )
+ {
+ z = pow( z , 3 );
+ }
+ else
+ {
+ z = ( z - 16 / 116 ) / 7.787;
+ }
+ float X= x *REF_X;
+ float Y= y *REF_Y;
+ float Z = z *REF_Z;
+ X = X/100.f; //R 0..1
+ Y= Y/100.f; //G 0..1
+ Z= Z/100.f; //B 0..1
+ float R1=3.240479*X-1.537150*Y -0.498535 *Z;
+ float G1=-0.969256*X+1.875992 *Y+0.041556* Z;
+ float B1=0.055648*X-0.204043*Y +1.057311*Z;
+ if( R1 > 0.0031308 )
+ R1 = 1.055 * ( pow( R1, 1./2.4 ) ) - 0.055;
+ else
+ R1 = 12.92 * R1;
+ if( G1 > 0.0031308 )
+ G1 = 1.055 * ( pow( G1, 1./2.4 ) ) - 0.055;
+ else
+ G1 = 12.92 * G1;
+ if( B1 > 0.0031308 )
+ B1 = 1.055 * ( pow( B1, 1./2.4 ) ) - 0.055;
+ else
+ B1 = 12.92 * B1;
+
+
+ Color rgb;
+ rgb.R= R1*100;
+ rgb.G= G1*100;
+ rgb.B= B1*100;
+ return rgb;
+
+}
+ int opencv_lab2rgb(char *fname, unsigned long fname_len)
+ {
+ // Error management variable
+ SciErr sciErr;
+ //variable info
+ int iRows = 0;
+ int iCols = 0;
+ int piRows = 0;
+ int piCols = 0;
+ int* piAddr=NULL;
+ int* piAddr1=NULL;
+ int* piAddr3=NULL;
+ double *value=NULL;
+ int* piLen = NULL;
+ char **val;
+ int i,j;
+ int error;
+ Color out;
+ double check;
+ Mat img,dst;
+ float L,a,b;
+ int *piAddrNew = NULL;
+ //checking output argument
+ CheckOutputArgument(pvApiCtx,1,1);
+ //checking if number of input arguments are 3
+
+ int k=0;
+ // SciErr getMatrixOfDoubleInList(void* _pvCtx, int* _piParent, int _iItemPos, int* _piRows, int* _piCols, double** _pdblReal)
+ double *pstDataR = NULL;
+ double *pstDataG = NULL;
+ double *pstDataB = NULL;
+ sciErr = getVarAddressFromPosition(pvApiCtx,1,&piAddr);
+ if(sciErr.iErr)
+ {
+ printError(&sciErr, 0);
+ return 0;
+ }
+ //retrive the matrix of the R values
+ sciErr = getMatrixOfDoubleInList(pvApiCtx, piAddr, 1, &iRows, &iCols, &pstDataR);
+ if(sciErr.iErr)
+ {
+ printError(&sciErr, 0);
+ return 0;
+ }
+ //retrive address of the list
+ sciErr = getVarAddressFromPosition(pvApiCtx,1,&piAddr);
+ if(sciErr.iErr)
+ {
+ printError(&sciErr, 0);
+ return 0;
+ }
+ //retrive the matrix of the G values
+ sciErr = getMatrixOfDoubleInList(pvApiCtx, piAddr, 2, &iRows, &iCols, &pstDataG);
+ if(sciErr.iErr)
+ {
+ printError(&sciErr, 0);
+ return 0;
+ }
+ //retrive address of the list
+ sciErr = getVarAddressFromPosition(pvApiCtx,1,&piAddr);
+ if(sciErr.iErr)
+ {
+ printError(&sciErr, 0);
+ return 0;
+ }
+ //retrive the matrix of the B values
+ sciErr = getMatrixOfDoubleInList(pvApiCtx, piAddr, 3, &iRows, &iCols, &pstDataB);
+ if(sciErr.iErr)
+ {
+ printError(&sciErr, 0);
+ return 0;
+ }
+ int m=0;
+ double *r1,*g1,*b1;
+ r1=(double *)malloc(sizeof(double)*iRows*iCols);
+ g1=(double *)malloc(sizeof(double)*iRows*iCols);
+ b1=(double *)malloc(sizeof(double)*iRows*iCols);
+ for(int i=0;i<iRows;i++)
+ {
+ for(int j=0;j<iCols;j++)
+ {
+ // Vec3d intensity = img.at<Vec3d>(i,j);
+ L= (pstDataR[k]);
+ a= (pstDataG[k]);
+ b= (pstDataB[k++]);
+ out=lab2rgb1(L,a,b);
+ //cout<<R<<"\n";
+ //cout<<G<<"\n";
+ //cout<<B<<"\n";
+ r1[m]=(out.R);
+ g1[m]=(out.G);
+ b1[m++]=(out.B);
+ }
+ }
+
+ sciErr = createList(pvApiCtx, nbInputArgument(pvApiCtx) + 1, 3, &piAddrNew);
+ if(sciErr.iErr)
+ {
+ printError(&sciErr, 0);
+ return 0;
+ }
+
+
+ //Adding the R value matrix to the list
+ //Syntax : createMatrixOfInteger32InList(void* _pvCtx, int _iVar, int* _piParent, int _iItemPos, int _iRows, int _iCols, const int* _piData)
+ sciErr = createMatrixOfDoubleInList(pvApiCtx, nbInputArgument(pvApiCtx)+1 , piAddrNew, 1, iRows,iCols, r1);
+ free(r1);
+ if(sciErr.iErr)
+ {
+ printError(&sciErr, 0);
+ return 0;
+ }
+
+ //Adding the G value matrix to the list
+ sciErr = createMatrixOfDoubleInList(pvApiCtx, nbInputArgument(pvApiCtx)+1 , piAddrNew, 2, iRows, iCols, g1);
+ free(g1);
+ if(sciErr.iErr)
+ {
+ printError(&sciErr, 0);
+ return 0;
+ }
+
+ //Adding the B value matrix to the list
+ sciErr = createMatrixOfDoubleInList(pvApiCtx, nbInputArgument(pvApiCtx)+1 , piAddrNew, 3, iRows, iCols, b1);
+ free(b1);
+ if(sciErr.iErr)
+ {
+ printError(&sciErr, 0);
+ return 0;
+ }
+
+
+ AssignOutputVariable(pvApiCtx, 1) = nbInputArgument(pvApiCtx) + 1;
+ ReturnArguments(pvApiCtx);
+ return 0;
+}
+}
generated by cgit (git 2.34.1) at 2024-12-19 10:07:41 +0000