First CommitHEADmaster

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;                
+    }
+}