diff options
Diffstat (limited to 'sci_gateway/cpp/lab2xyz.cpp')
| -rw-r--r-- | sci_gateway/cpp/lab2xyz.cpp | 271 |
1 files changed, 271 insertions, 0 deletions
diff --git a/sci_gateway/cpp/lab2xyz.cpp b/sci_gateway/cpp/lab2xyz.cpp new file mode 100644 index 0000000..28ed422 --- /dev/null +++ b/sci_gateway/cpp/lab2xyz.cpp @@ -0,0 +1,271 @@ +/* ==================================================================== */ +/* Author :Priyanka Hiranandani NIT Surat */ +/* ==================================================================== */ +/* overloaded function */ +/* Syntax : return_matrix_3_xyz_values=lab2xyz(matrix_of_3_labvalues a) */ +/* ==================================================================== */ +/* ==================================================================== */ +#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 95.047; // Observer= 2°, Illuminant= D65 + #define REF_Y 100.000; + #define REF_Z 108.883; + + struct Color + { + float R,G,B,X,Y,Z,L,a,b; + }; + + Color lab2xyz1(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; + } + + Color xyz; + xyz.X= x * REF_X; + xyz.Y = y * REF_Y; + xyz.Z = z * REF_Z; + return xyz; + + } + + Color lab2xyz2(float l, float a, float b,char *s) + { + double REF_X1; + double REF_Y1; + double REF_Z1; + //if else condition for checking various input arguments and assign values accordingly + 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; + } + Color xyz; + xyz.X= x * REF_X1; + xyz.Y = y * REF_Y1; + xyz.Z = z * REF_Z1; + return xyz; + } + int lab2xyz(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 R,G,B; + int k=0; + int *piAddrNew = NULL; + double *pstDataR=NULL; + double *pstDataG=NULL; + double *pstDataB=NULL; + double *r,*g,*b; + int m=0; + //checking output argument + CheckOutputArgument(pvApiCtx,1,1); + //retrive address of the list + 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; + } + r=(double *)malloc(sizeof(double)*iRows*iCols); + g=(double *)malloc(sizeof(double)*iRows*iCols); + b=(double *)malloc(sizeof(double)*iRows*iCols); + for(int i=0;i<iRows;i++) + { + for(int j=0;j<iCols;j++) + { + R=(pstDataR[k]); + G=(pstDataG[k]); + B= (pstDataB[k++]); + out=lab2xyz1(R,G,B); + r[m]=(out.X)/100; + g[m]=(out.Y)/100; + b[m++]=(out.Z)/100; + } + } + sciErr = createList(pvApiCtx, nbInputArgument(pvApiCtx) + 1, 3, &piAddrNew); + if(sciErr.iErr) + { + printError(&sciErr, 0); + return 0; + } + //Adding the R value matrix to the list + sciErr = createMatrixOfDoubleInList(pvApiCtx, nbInputArgument(pvApiCtx)+1 , piAddrNew, 1, iRows,iCols, r); + free(r); + 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, g); + free(g); + 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, b); + free(b); + if(sciErr.iErr) + { + printError(&sciErr, 0); + return 0; + } + AssignOutputVariable(pvApiCtx, 1) = nbInputArgument(pvApiCtx) + 1; + ReturnArguments(pvApiCtx); + return 0; + } +} |