First CommitHEADmaster

1 files changed, 244 insertions, 0 deletions

diff --git a/sci_gateway/cpp/opencv_rgb2lab.cpp b/sci_gateway/cpp/opencv_rgb2lab.cpp
new file mode 100644
index 0000000..c031790
--- /dev/null
+++ b/sci_gateway/cpp/opencv_rgb2lab.cpp
@@ -0,0 +1,244 @@
+/*  function-rgb2lab
+ Author :Tess Zacharias      */
+#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,Ls,as,bs,r,g,b;
+};
+Color rgb2lab1(float R, float G, float B)
+{
+        
+         float r = R/255.f; //R 0..1
+         float g = G/255.f; //G 0..1
+         float b = B/255.f; //B 0..1
+         float eps = 216.f/24389.f;
+         float k = 24389.f/27.f;
+         float  fx,fy,fz;
+         if (r <= 0.04045)
+             r = r/12.92;
+         else
+             r = (float)pow((r+0.055)/1.055,2.4);
+         if (g <= 0.04045)
+             g = g/12.92;
+         else
+             g = (float)pow((g+0.055)/1.055,2.4);
+         if (b <= 0.04045)
+             b = b/12.92;
+         else
+             b = (float)pow((b+0.055)/1.055,2.4);
+         float X =  0.4124*r     + 0.3576*g + 0.1805*b;
+         float Y =  0.2126*r     + 0.7152 *g + 0.0722 *b;
+         float Z =  0.0193*r     + 0.1192*g + 0.9505 *b;
+         float xr = X/REF_X;
+         float yr = Y/REF_Y;
+         float zr = Z/REF_Z;
+         if ( xr > eps )
+            fx =  (float)pow(xr, 1/3.);
+         else
+            fx = (float) ((k * xr + 16.) / 116.);
+         if ( yr > eps )
+            fy =  (float)pow(yr, 1/3.);
+         else
+            fy = (float) ((k * yr + 16.) / 116.);
+         if ( zr > eps )
+            fz =  (float)pow(zr, 1/3.);
+         else
+            fz = (float) ((k * zr + 16.) / 116);
+    Color lab;
+    lab.Ls= ( 116 * fy ) - 16;
+    lab.as= 500*(fx-fy);
+    lab.bs= 200*(fy-fz);
+    return lab;
+
+} 
+ Color rgb2lab2(float R, float G, 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 eps = 216.f/24389.f;
+         float k = 24389.f/27.f;
+         float  fx,fy,fz,X,Y,Z,xr,yr,zr;
+         float r = R/255.f; //R 0..1
+         float g = G/255.f; //G 0..1
+         float b = B/255.f; //B 0..1
+         if (r <= 0.04045)
+             r = r/12.92;
+         else
+             r = (float)pow((r+0.055)/1.055,2.4);
+         if (g <= 0.04045)
+             g = g/12.92;
+         else
+             g = (float)pow((g+0.055)/1.055,2.4);
+         if (b <= 0.04045)
+             b = b/12.92;
+         else
+             b = (float)pow((b+0.055)/1.055,2.4);
+         X =  0.4124*r     + 0.3576*g + 0.1805*b;
+         Y =  0.2126*r     + 0.7152 *g + 0.0722 *b;
+         Z =  0.0193*r     + 0.1192*g + 0.9505 *b;
+         xr = X/REF_X1;
+         yr = Y/REF_Y1;
+         zr = Z/REF_Z1;
+         if ( xr > eps )
+            fx =  (float)pow(xr, 1/3.);
+         else
+            fx = (float) ((k * xr + 16.) / 116.);
+         if ( yr > eps )
+            fy =  (float)pow(yr, 1/3.);
+         else
+            fy = (float) ((k * yr + 16.) / 116.);
+         if ( zr > eps )
+            fz =  (float)pow(zr, 1/3.);
+         else
+            fz = (float) ((k * zr + 16.) / 116);
+        Color lab;
+        lab.Ls= ( 116 * fy ) - 16;
+        lab.as= 500*(fx-fy);
+        lab.bs= 200*(fy-fz);
+        return lab;
+
+}  
+int opencv_rgb2lab(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* piAddr2=NULL;
+        int* piAddr1=NULL;
+        int* piAddr3=NULL;
+        int* piaddr=NULL;
+        char *value=NULL;
+        int* piLen		= NULL;
+        char **val;
+        int i,j;
+        int error;
+        Color out;
+        double check;
+        Mat img;
+        float R,G,B;
+        int *piAddrNew = NULL;                      
+        CheckOutputArgument(pvApiCtx,1,1);
+        int m=0;
+        retrieveImage(img,1);
+        double *r1,*g1,*b1;
+        r1=(double *)malloc(sizeof(double)*img.rows*img.cols);
+        g1=(double *)malloc(sizeof(double)*img.rows*img.cols);
+        b1=(double *)malloc(sizeof(double)*img.rows*img.cols);
+        for(int i=0;i<img.rows;i++)  
+          {
+            for(int j=0;j<img.cols;j++)
+               {
+                  Vec3b intensity = img.at<Vec3b>(i,j);
+                  R= intensity.val[2];
+                  G= intensity.val[1];
+                  B= intensity.val[0];
+                  out=rgb2lab1(R,G,B);
+                  r1[m]=out.Ls;
+                  g1[m]=out.as;
+                  b1[m++]=out.bs;  
+                  }
+          }
+    //creating the list that will be passed to the Scilab enviroment
+    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, img.rows,img.cols, 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, img.rows, img.cols, 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, img.rows, img.cols, b1);
+    free(b1);
+    if(sciErr.iErr)
+    {
+        printError(&sciErr, 0);
+        return 0;
+    }  
+
+           // free(checker);  
+            AssignOutputVariable(pvApiCtx, 1) = nbInputArgument(pvApiCtx) + 1;
+            ReturnArguments(pvApiCtx);            
+            return 0;                
+}
+}