/*
    Copyright (c) 2008-2019 Jan W. Krieger (<jan@jkrieger.de>)

    

    This software is free software: you can redistribute it and/or modify
    it under the terms of the GNU Lesser General Public License (LGPL) as published by
    the Free Software Foundation, either version 2.1 of the License, or
    (at your option) any later version.

    This program 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 Lesser General Public License (LGPL) for more details.

    You should have received a copy of the GNU Lesser General Public License (LGPL)
    along with this program.  If not, see <http://www.gnu.org/licenses/>.
*/

#include "jkqtcommon/jkqtpbasicimagetools.h"
#include <QPainter>
#include <QApplication>
#include <QDir>
#include <QFileInfo>
#ifdef QT_XML_LIB
#  include <QtXml/QtXml>
#endif

const int JKQTPImageTools::PALETTE_ICON_WIDTH = 64;
const int JKQTPImageTools::PALETTE_IMAGEICON_HEIGHT = 64;
const int JKQTPImageTools::LUTSIZE = 256;

QList<int*> JKQTPImageTools::global_jkqtpimagetools_lutstore = QList<int*>();


bool JKQTPImagePlot_buildDefinedPaletteLessThan(const QPair<double, QRgb> &s1, const QPair<double, QRgb> &s2)
 {
     return s1.first<s2.first;
 }

void JKQTPImagePlot_buildDefinedPaletteLinInterpolate(int* lut, QList<QPair<double, QRgb> > items, int lut_size) {
    qSort(items.begin(), items.end(), JKQTPImagePlot_buildDefinedPaletteLessThan);

    double dmin=items.first().first;
    double dmax=items.last().first;
    double delta=(dmax-dmin)/double(lut_size);
    int i=0;
    int j=0;
    for (double v=dmin; v<=dmax; v=v+delta) {
     double r=0;
     double g=0;
     double b=0;
     if (j+1<items.size()) {
         double del=items[j+1].first-items[j].first;
         double rel=(v-items[j].first)/del;
         r=double(qRed(items[j].second))+rel*(double(qRed(items[j+1].second))-double(qRed(items[j].second)));
         g=double(qGreen(items[j].second))+rel*(double(qGreen(items[j+1].second))-double(qGreen(items[j].second)));
         b=double(qBlue(items[j].second))+rel*(double(qBlue(items[j+1].second))-double(qBlue(items[j].second)));
         while (j+1<items.size() && v>=items[j+1].first) j++;
     } else {
         r=double(qRed(items.last().second));
         g=double(qGreen(items.last().second));
         b=double(qBlue(items.last().second));
     }
     if (i<=lut_size) lut[i]=qRgb(static_cast<int>(r), static_cast<int>(g), static_cast<int>(b));
     i++;
    }
}

void JKQTPImagePlot_buildDefinedPaletteLinInterpolate(int* lut, int N, ...) {
    QList<QPair<double, QRgb> > items;

    va_list args;
    va_start(args, N);

    for (int i=0; i<N; i++) {
     double colval=va_arg(args, double);
     QRgb col=va_arg(args,QRgb);
     items<<qMakePair(colval, col);
    }

    va_end(args);

    JKQTPImagePlot_buildDefinedPaletteLinInterpolate(lut, items);
}


void JKQTPImagePlot_buildDefinedPalette(int* lut, QList<QPair<double, QRgb> > items)  {
    qSort(items.begin(), items.end(), JKQTPImagePlot_buildDefinedPaletteLessThan);

    double dmin=items.first().first;
    double dmax=items.last().first;
    double delta=(dmax-dmin)/double(JKQTPImageTools::LUTSIZE);
    int i=0;
    int j=0;
    for (double v=dmin; v<=dmax; v=v+delta) {
     double r=0;
     double g=0;
     double b=0;
     if (j+1<items.size()) {
         r=qRed(items[j].second);
         g=qGreen(items[j].second);
         b=qBlue(items[j].second);
         while (j+1<items.size() && v>=items[j+1].first) j++;
     }
     if (i<=JKQTPImageTools::LUTSIZE) lut[i]=qRgb(static_cast<int>(r), static_cast<int>(g), static_cast<int>(b));
     i++;
    }
}

void JKQTPImagePlot_buildDefinedPalette(int* lut, int N, ...) {
    QList<QPair<double, QRgb> > items;

    va_list args;
    va_start(args, N);

    for (int i=0; i<N; i++) {
     double colval=va_arg(args, double);
     QRgb col=va_arg(args,QRgb);
     items<<qMakePair(colval, col);
    }

    va_end(args);

    JKQTPImagePlot_buildDefinedPalette(lut, items);
}


double JKQTPImagePlot_QStringToDouble(QString value) {
    QString v=value;
    if (value.contains(',')) {
     v=value.replace(',', '.');
    }
    QLocale loc=QLocale::c();
    loc.setNumberOptions(QLocale::OmitGroupSeparator);
    return loc.toDouble(v) ;
}

static QStringList JKQTPImagePlot_palettesSearchPaths=QStringList();
static QMap<int, QList<QPair<double, QRgb> > > JKQTPImagePlot_lutsFromFiles;

template <class T1, class T2>
bool JKQTPImagePlot_QPairCompareFirst(const QPair<T1, T2> &s1, const QPair<T1, T2> &s2) {
    return s1.first < s2.first;
}

QStringList JKQTPImagePlot_getPredefinedPalettes()  {
    static QStringList sl;

    if (JKQTPImagePlot_palettesSearchPaths.isEmpty()) {
     JKQTPImagePlot_palettesSearchPaths<<":/colorpalettes/";
     JKQTPImagePlot_palettesSearchPaths<<(QApplication::applicationDirPath()+"/colorpalettes/");
     JKQTPImagePlot_palettesSearchPaths<<(QApplication::applicationDirPath()+"/assets/colorpalettes/");
    }

    if (sl.size()<=0) {
     sl<<QObject::tr("gray");
     sl<<QObject::tr("inv. gray");
     sl<<QObject::tr("red");
     sl<<QObject::tr("inv. red");
     sl<<QObject::tr("green");
     sl<<QObject::tr("inv. green");
     sl<<QObject::tr("blue");
     sl<<QObject::tr("inv. blue");
     sl<<QObject::tr("cyan");
     sl<<QObject::tr("inv. cyan");
     sl<<QObject::tr("yellow");
     sl<<QObject::tr("inv. yellow");
     sl<<QObject::tr("magenta");
     sl<<QObject::tr("inv. magenta");

     sl<<QObject::tr("Matlab");
     sl<<QObject::tr("inv. Matlab");
     sl<<QObject::tr("RYGB");
     sl<<QObject::tr("inv. RYGB");
     sl<<QObject::tr("HSV");
     sl<<QObject::tr("inv. HSV");
     sl<<QObject::tr("rainbow");
     sl<<QObject::tr("inv. rainbow");
     sl<<QObject::tr("AFM hot");
     sl<<QObject::tr("inv. AFM hot");
     sl<<QObject::tr("ocean");
     sl<<QObject::tr("inv. ocean");
     sl<<QObject::tr("trafficlight");
     sl<<QObject::tr("inv. trafficl.");
     sl<<QObject::tr("BlMaYe");
     sl<<QObject::tr("YeMaBl");
     sl<<QObject::tr("BlYe");
     sl<<QObject::tr("YeBl");
     sl<<QObject::tr("blue-white-red");
     sl<<QObject::tr("red-white-blue");
     sl<<QObject::tr("BBlRdYe");
     sl<<QObject::tr("GnRdVi");
     sl<<QObject::tr("BW print");
     sl<<QObject::tr("inv. BW print");
     sl<<QObject::tr("BrBG");
     sl<<QObject::tr("PuOr");
     sl<<QObject::tr("green-blue");
     sl<<QObject::tr("blue-green");
     sl<<QObject::tr("YeGnBu");
     sl<<QObject::tr("steps: BrBG");
     sl<<QObject::tr("steps: PuOr");
     sl<<QObject::tr("steps: GnBl");
     sl<<QObject::tr("steps: BlGn");
     sl<<QObject::tr("steps: YeGnBu");

     sl<<QObject::tr("cyan-white");
     sl<<QObject::tr("white-cyan");
     sl<<QObject::tr("yellow-white");
     sl<<QObject::tr("white-yellow");
     sl<<QObject::tr("magenta-white");
     sl<<QObject::tr("white-magenta");

     sl<<QObject::tr("blue-green-red");
     sl<<QObject::tr("red-green-blue");
     sl<<QObject::tr("magenta-yellow");
     sl<<QObject::tr("yellow-magenta");
     sl<<QObject::tr("red-blue");
     sl<<QObject::tr("blue-red");
     //sl<<QObject::tr("");

#ifdef QT_XML_LIB
     int palID=sl.size();
#endif
     for (int i=0; i<JKQTPImagePlot_palettesSearchPaths.size(); i++) {
         QDir d(JKQTPImagePlot_palettesSearchPaths[i]);
         QStringList nameFilters;
         nameFilters<<"*.pal";
         nameFilters<<"*.rgb";
         nameFilters<<"*.csv";
#ifdef QT_XML_LIB
         nameFilters<<"*.xml";
#endif
         QFileInfoList l=d.entryInfoList(nameFilters);
         for (int j=0; j<l.size(); j++) {
          QString ext=l[j].suffix().toLower();
          if (ext=="rgb" || ext=="pal" || ext=="csv") {
              QFile file(l[j].absoluteFilePath());
              if (file.open(QIODevice::ReadOnly|QIODevice::Text)) {
               QString txt=file.readAll();
               QStringList slt=txt.split('\n');
               bool has4=false;
               bool rgb255=false;
               QList<QPair<double, QRgb> > pal;
               QRegExp rx3("\\s*([0-9eE.+-]+)\\s*([,\\t ])\\s*([0-9eE.+-]+)\\s*\\2\\s*([0-9eE.+-]+)\\s*", Qt::CaseInsensitive);
               rx3.setMinimal(false);
               QRegExp rx4("\\s*([0-9eE.+-]+)\\s*([,\\t ])\\s*([0-9eE.+-]+)\\s*\\2\\s*([0-9eE.+-]+)\\s*\\2\\s*([0-9eE.+-]+)\\s*", Qt::CaseInsensitive);
               rx4.setMinimal(false);
               // determine format
               for (int i=slt.size()-1; i>=0; i--) {
                   slt[i]=slt[i].trimmed();
                   if (rx4.indexIn(slt[i])>=0) {
                    has4=true;
                    double r=JKQTPImagePlot_QStringToDouble(rx4.cap(3));
                    double g=JKQTPImagePlot_QStringToDouble(rx4.cap(4));
                    double b=JKQTPImagePlot_QStringToDouble(rx4.cap(5));
                    if (r>1.0 || g>1.0 || b>1.0) {
                        rgb255=true;
                    }
                   } else if (rx3.indexIn(slt[i])>=0) {
                    has4=false;
                    double r=JKQTPImagePlot_QStringToDouble(rx3.cap(1));
                    double g=JKQTPImagePlot_QStringToDouble(rx3.cap(3));
                    double b=JKQTPImagePlot_QStringToDouble(rx3.cap(4));
                    if (r>1.0 || g>1.0 || b>1.0) {
                        rgb255=true;
                    }
                   } else {
                    slt.removeAt(i);
                   }
               }
               // read data
               bool ok=true;
               for (int i=0; i<slt.size(); i++) {
                   double x=0;
                   double r=0, g=0, b=0;
                   if (has4 && rx4.indexIn(slt[i])>=0) {
                    x=JKQTPImagePlot_QStringToDouble(rx4.cap(1));
                    r=JKQTPImagePlot_QStringToDouble(rx4.cap(3));
                    g=JKQTPImagePlot_QStringToDouble(rx4.cap(4));
                    b=JKQTPImagePlot_QStringToDouble(rx4.cap(5));
                   } else if (!has4 && rx3.indexIn(slt[i])>=0) {
                    x=i;
                    r=JKQTPImagePlot_QStringToDouble(rx3.cap(1));
                    g=JKQTPImagePlot_QStringToDouble(rx3.cap(3));
                    b=JKQTPImagePlot_QStringToDouble(rx3.cap(4));
                    //qDebug()<<r<<g<<b;
                   } else {
                    ok=false;
                    break;
                   }
                   if (ok) {
                    if (!rgb255) {
                        pal<<qMakePair(x, qRgb(qBound(0,static_cast<int>(round(255*r)), 255), qBound(0,static_cast<int>(round(255*g)), 255), qBound(0,static_cast<int>(round(255*b)), 255)));
                    } else {
                        pal<<qMakePair(x, qRgb(qBound(0,static_cast<int>(round(r)), 255), qBound(0,static_cast<int>(round(g)), 255), qBound(0,static_cast<int>(round(b)), 255)));
                    }
                   }
               }
               //qDebug()<<"read: "<<l[j].fileName()<<": size="<<pal.size()<<" ok="<<ok<<" has4="<<has4<<" rgb255="<<rgb255;
               if (ok && pal.size()>1) {
                   sl<<l[j].fileName();
                   qStableSort(pal.begin(), pal.end(), JKQTPImagePlot_QPairCompareFirst<double, QRgb>);
                   JKQTPImagePlot_lutsFromFiles[sl.size()-1]=pal;
                   //qDebug()<<"   added as "<<sl.size()-1;
#ifdef QT_XML_LIB
                   palID=l.size();
#endif
               }

              }
#ifdef QT_XML_LIB
          } else if (ext=="xml") {
              QDomDocument doc("palette");
              QFile file(l[j].absoluteFilePath());
              if (file.open(QIODevice::ReadOnly)) {
               if (doc.setContent(&file)) {
                   file.close();
                   QDomElement docElem = doc.documentElement();

                   QDomElement n = docElem.firstChildElement("ColorMap");
                   if (docElem.nodeName()=="ColorMap") n=docElem;
                   while(!n.isNull()) {
                    QString name=n.attribute("name", QString("palette %1").arg(palID+1));
                    QString colorspace=n.attribute("space", "RGB").toLower();
                    QList<QPair<double, QRgb> > pal;
                    if (colorspace=="rgb"){
                        QDomElement e = n.firstChildElement("Point");
                        while(!e.isNull()) {
                         double x=JKQTPImagePlot_QStringToDouble(e.attribute("x", "0"));
                         //double o=JKQTPImagePlot_QStringToDouble(e.attribute("o", "0"));
                         double r=JKQTPImagePlot_QStringToDouble(e.attribute("r", "0"));
                         double g=JKQTPImagePlot_QStringToDouble(e.attribute("g", "0"));
                         double b=JKQTPImagePlot_QStringToDouble(e.attribute("b", "0"));
                         pal<<qMakePair(x, qRgb(qBound(0,static_cast<int>(round(255*r)), 255), qBound(0,static_cast<int>(round(255*g)), 255), qBound(0,static_cast<int>(round(255*b)), 255)));
                         e = e.nextSiblingElement("Point");
                        }
                        if (pal.size()>1) {
                         sl<<name;
                         qStableSort(pal.begin(), pal.end(), JKQTPImagePlot_QPairCompareFirst<double, QRgb>);
                         JKQTPImagePlot_lutsFromFiles[sl.size()-1]=pal;
                         palID=sl.size();
                         //qDebug()<<"read "<<l[j].fileName()<<":   "<<name<<" ("<<colorspace<<") size="<<pal.size();
                        }
                    }

                    n = n.nextSiblingElement("ColorMap");
                   }
               }
              }
#endif
          }
         }
     }
    }


    return sl;
}

void JKQTPImagePlot_freeLUTs(QList<int *> &lutstore) {
    QList<int *> ls=lutstore;
    lutstore.clear();
    for (int i=0; i<ls.size(); i++) {
     if (ls[i]) free(ls[i]);
    }
}

int* JKQTPImagePlot_getCreateLUT(QList<int*>& lutstore, JKQTPMathImageColorPalette palette) {
    int* lut_used=nullptr;

    if (palette == JKQTPMathImageRED)
    {
     //int*& plut=lutstore[palette];
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=l/static_cast<double>(JKQTPImageTools::LUTSIZE);
              plut[l]=qRgb(static_cast<int>(255.0*v), 0, 0);
          }
          plut[JKQTPImageTools::LUTSIZE+1]=plut[JKQTPImageTools::LUTSIZE];
         }
     }
     lut_used=lutstore[palette];

    }
    else if (palette == JKQTPMathImageGREEN)
    {
     //int*& plut=&lut_green;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=l/static_cast<double>(JKQTPImageTools::LUTSIZE);
              plut[l]=qRgb(0, static_cast<int>(255.0*v), 0);
          }
          plut[JKQTPImageTools::LUTSIZE+1]=plut[JKQTPImageTools::LUTSIZE];
         }
     }
     lut_used=lutstore[palette];

    }
    else if (palette == JKQTPMathImageBLUE)
    {
     //int*& plut=&lut_blue;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=l/static_cast<double>(JKQTPImageTools::LUTSIZE);
              plut[l]=qRgb(0, 0, static_cast<int>(255.0*v));
          }
          plut[JKQTPImageTools::LUTSIZE+1]=plut[JKQTPImageTools::LUTSIZE];
         }
     }
     lut_used=lutstore[palette];

    }
    else if (palette == JKQTPMathImageGRAY)
    {
     //int*& plut=&lut_gray;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=l/static_cast<double>(JKQTPImageTools::LUTSIZE);
              plut[l]=qRgb(static_cast<int>(255.0*v),
                        static_cast<int>(255.0*v),
                        static_cast<int>(255.0*v));
          }
          plut[JKQTPImageTools::LUTSIZE+1]=plut[JKQTPImageTools::LUTSIZE];
         }
     }
     lut_used=lutstore[palette];
    }
    else if (palette == JKQTPMathImageALPHA)
    {
     //int*& plut=&lut_gray;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=l/static_cast<double>(JKQTPImageTools::LUTSIZE);
              plut[l]=qRgba(255,255,255,
                        static_cast<int>(255.0*v));
          }
          plut[JKQTPImageTools::LUTSIZE+1]=plut[JKQTPImageTools::LUTSIZE];
         }
     }
     lut_used=lutstore[palette];
    }
    else if (palette == JKQTPMathImageINVERTED_ALPHA)
    {
     //int*& plut=&lut_gray;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=l/static_cast<double>(JKQTPImageTools::LUTSIZE);
              plut[l]=qRgba(255,255,255,
                        static_cast<int>(255.0*v));
          }
          plut[JKQTPImageTools::LUTSIZE+1]=plut[JKQTPImageTools::LUTSIZE];
         }
     }
     lut_used=lutstore[palette];
    }
    else if (palette == JKQTPMathImageINVERTEDRED)
    {
     //int*& plut=&lut_invred;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=l/static_cast<double>(JKQTPImageTools::LUTSIZE);
              plut[l]=qRgb(static_cast<int>(255.0*(1.0-v)), 0, 0);
          }
          plut[JKQTPImageTools::LUTSIZE+1]=plut[JKQTPImageTools::LUTSIZE];
         }
     }
     lut_used=lutstore[palette];

    }
    else if (palette == JKQTPMathImageINVERTEDGREEN)
    {
     //int*& plut=&lut_invgreen;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=l/static_cast<double>(JKQTPImageTools::LUTSIZE);
              plut[l]=qRgb(0, static_cast<int>(255.0*(1.0-v)), 0);
          }
          plut[JKQTPImageTools::LUTSIZE+1]=plut[JKQTPImageTools::LUTSIZE];
         }
     }
     lut_used=lutstore[palette];
    }
    else if (palette == JKQTPMathImageINVERTEDBLUE)
    {
     //int*& plut=&lut_invblue;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=l/static_cast<double>(JKQTPImageTools::LUTSIZE);
              plut[l]=qRgb(0, 0, static_cast<int>(255.0*(1.0-v)));
          }
          plut[JKQTPImageTools::LUTSIZE+1]=plut[JKQTPImageTools::LUTSIZE];
         }
     }
     lut_used=lutstore[palette];
    }
    else if (palette == JKQTPMathImageINVERTEDGRAY)
    {
     //int*& plut=&lut_invgray;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=1.0-(l/static_cast<double>(JKQTPImageTools::LUTSIZE));
              plut[l]=qRgb(static_cast<int>(255.0*v),
                        static_cast<int>(255.0*v),
                        static_cast<int>(255.0*v));
          }
          plut[JKQTPImageTools::LUTSIZE+1]=plut[JKQTPImageTools::LUTSIZE];
         }
     }
     lut_used=lutstore[palette];
    }

    else if (palette == JKQTPMathImageMATLAB)
    {
     //int*& plut=&lut_matlab;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=l/static_cast<double>(JKQTPImageTools::LUTSIZE);
              double r = 382.5 - 1020.0 * std::abs(v - 0.75);
              if (r > 255.0)
               r = 255.0;
              else if (r < 0.0)
               r = 0.0;

              double g = 382.5 - 1020.0 * std::abs(v - 0.5);
              if (g > 255.0)
               g = 255.0;
              else if (g < 0.0)
               g = 0.0;

              double b = 382.5 - 1020.0 * std::abs(v - 0.25);
              if (b > 255.0)
               b = 255.0;
              else if (b < 0.0)
               b = 0.0;

              plut[l]=qRgb(static_cast<int>(r), static_cast<int>(g), static_cast<int>(b));
          }
         }
     }
     lut_used=lutstore[palette];
    }

     else if (palette == JKQTPMathImageINVERTED_MATLAB)
     {
         //int*& plut=&lut_matlab_inverted;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          if (plut!=nullptr) {
              for (int l=JKQTPImageTools::LUTSIZE; l>=0; l--) {
               double v=(JKQTPImageTools::LUTSIZE-l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
               double r = 382.5 - 1020.0 * std::abs(v - 0.75);
               if (r > 255.0)
                   r = 255.0;
               else if (r < 0.0)
                   r = 0.0;

               double g = 382.5 - 1020.0 * std::abs(v - 0.5);
               if (g > 255.0)
                   g = 255.0;
               else if (g < 0.0)
                   g = 0.0;

               double b = 382.5 - 1020.0 * std::abs(v - 0.25);
               if (b > 255.0)
                   b = 255.0;
               else if (b < 0.0)
                   b = 0.0;

               plut[l]=qRgb(static_cast<int>(r), static_cast<int>(g), static_cast<int>(b));
              }
          }
         }
         lut_used=lutstore[palette];

    }
     else if (palette == JKQTPMathImageRYGB) //gnuplot: 30,-13,-23
     {
         //int*& plut=&lut_rygb;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          if (plut!=nullptr) {
              for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
               double v=l/static_cast<double>(JKQTPImageTools::LUTSIZE);
               double r = 796.875*v - 199.21875;
               if (r > 255.0)
                   r = 255.0;
               else if (r < 0.0)
                   r = 0.0;

               double g = 255.0 * std::sin(M_PI*v);

               double b = 255.0 - 765.0 * v;
               if (b < 0.0)
                   b = 0.0;
               plut[l]=qRgb(static_cast<int>(r), static_cast<int>(g), static_cast<int>(b));
              }
          }
         }
         lut_used=lutstore[palette];

    }
     else if (palette == JKQTPMathImageINVERTED_RYGB) //gnuplot: 30,-13,-23
     {
         //int*& plut=&lut_rygb_inverted;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          if (plut!=nullptr) {
              for (int l=JKQTPImageTools::LUTSIZE; l>=0; l--) {
               double v=(JKQTPImageTools::LUTSIZE-l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
               double r = 796.875*v - 199.21875;
               if (r > 255.0)
                   r = 255.0;
               else if (r < 0.0)
                   r = 0.0;

               double g = 255.0 * std::sin(M_PI*v);

               double b = 255.0 - 765.0 * v;
               if (b < 0.0)
                   b = 0.0;
               plut[l]=qRgb(static_cast<int>(r), static_cast<int>(g), static_cast<int>(b));
              }
          }
         }
         lut_used=lutstore[palette];

    }
    else if (palette == JKQTPMathImageHSV)
    {
     //int*& plut=&lut_hsv;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=l/static_cast<double>(JKQTPImageTools::LUTSIZE);
              int h = static_cast<int>(floor(6*v));
              double f = 6*v-double(h);

              switch (h)
              {
               case 0: plut[l]=qRgb(255, static_cast<int>(255.0*f), 0); break;
               case 1: plut[l]=qRgb(static_cast<int>(255.0*(1-f)), 255, 0); break;
               case 2: plut[l]=qRgb(0, 255, static_cast<int>(255.0*f)); break;
               case 3: plut[l]=qRgb(0, static_cast<int>(255.0*(1-f)), 255); break;
               case 4: plut[l]=qRgb(static_cast<int>(255.0*f), 0, 255); break;
               case 5: plut[l]=qRgb(255, 0, static_cast<int>(255.0*(1-f))); break;
               case 6: plut[l]=qRgb(255, static_cast<int>(255.0*f), 0); break;
              }
          }
         }
     }
     lut_used=lutstore[palette];
    }
    else if (palette == JKQTPMathImageINVERTED_HSV)
    {
     //int*& plut=&lut_invhsv;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=l/static_cast<double>(JKQTPImageTools::LUTSIZE);
              int h = static_cast<int>(floor(6.0-6.0*v));
              double f = 6.0-6.0*v-double(h);

              switch (h)
              {
               case 0: plut[l]=qRgb(255, static_cast<int>(255.0*f), 0); break;
               case 1: plut[l]=qRgb(static_cast<int>(255.0*(1-f)), 255, 0); break;
               case 2: plut[l]=qRgb(0, 255, static_cast<int>(255.0*f)); break;
               case 3: plut[l]=qRgb(0, static_cast<int>(255.0*(1-f)), 255); break;
               case 4: plut[l]=qRgb(static_cast<int>(255.0*f), 0, 255); break;
               case 5: plut[l]=qRgb(255, 0, static_cast<int>(255.0*(1-f))); break;
               case 6: plut[l]=qRgb(255, static_cast<int>(255.0*f), 0); break;
              }
          }
         }
     }
     lut_used=lutstore[palette];

    }
    else if (palette == JKQTPMathImageRAINBOW) //gnuplot: 33,13,10
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         //std::cout<<"!!! creating rainbow lut\n";
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=l/static_cast<double>(JKQTPImageTools::LUTSIZE);
              double r = 255.0*std::abs(2.0*v-0.5);
              if (r > 255.0)
               r = 255.0;

              double g = 255.0*sin(M_PI*v);

              double b = 255.0*cos(0.5*M_PI*v);

              plut[l]=qRgb(static_cast<int>(r), static_cast<int>(g), static_cast<int>(b));
          }
         }
     }
     lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageINVERTED_RAINBOW) //gnuplot: 33,13,10
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         //std::cout<<"!!! creating rainbow lut\n";
         if (plut!=nullptr) {
          for (int l=JKQTPImageTools::LUTSIZE; l>=0; l--) {
              double v=(JKQTPImageTools::LUTSIZE-l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
              double r = 255.0*std::abs(2.0*v-0.5);
              if (r > 255.0)
               r = 255.0;

              double g = 255.0*sin(M_PI*v);

              double b = 255.0*cos(0.5*M_PI*v);

              plut[l]=qRgb(static_cast<int>(r), static_cast<int>(g), static_cast<int>(b));
          }
         }
     }
     lut_used=lutstore[palette];
    }
    else if (palette == JKQTPMathImageHOT) //gnuplot: 21,22,23
    {
     //int*& plut=&lut_hot;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=l/static_cast<double>(JKQTPImageTools::LUTSIZE);
              double r = 765.0*v;
              if (r > 255.0)
                   r = 255.0;

              double g = 765.0*v-255.0;
              if (g > 255.0)
                   g = 255.0;
              else if (g < 0.0)
                   g = 0.0;

              double b = 765.0*v-510.0;
              if (b < 0.0)
                   b = 0.0;
              plut[l]=qRgb(static_cast<int>(r), static_cast<int>(g), static_cast<int>(b));
          }
         }
     }
     lut_used=lutstore[palette];
    }
    else if (palette == JKQTPMathImageINVERTED_HOT) //gnuplot: 21,22,23
    {
     //int*& plut=&lut_hot_inverted;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=JKQTPImageTools::LUTSIZE; l>=0; l--) {
              double v=(JKQTPImageTools::LUTSIZE-l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
              double r = 765.0*v;
              if (r > 255.0)
                   r = 255.0;

              double g = 765.0*v-255.0;
              if (g > 255.0)
                   g = 255.0;
              else if (g < 0.0)
                   g = 0.0;

              double b = 765.0*v-510.0;
              if (b < 0.0)
                   b = 0.0;
              plut[l]=qRgb(static_cast<int>(r), static_cast<int>(g), static_cast<int>(b));
          }
         }
     }
     lut_used=lutstore[palette];
    }
    else if (palette == JKQTPMathImageOCEAN) //gnuplot: 23,28,3
    {
     //int*& plut=&lut_ocean;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=l/static_cast<double>(JKQTPImageTools::LUTSIZE);
              double r = 765.0*v-510.0;
              if (r < 0.0)
               r = 0.0;

              double g = std::abs(382.5*v-127.5);

              double b = 255.0*v;
              plut[l]=qRgb(static_cast<int>(r), static_cast<int>(g), static_cast<int>(b));
          }
         }
     }
     lut_used=lutstore[palette];
    }
    else if (palette == JKQTPMathImageINVERTED_OCEAN) //gnuplot: 23,28,3
    {
     //int*& plut=&lut_ocean_inverted;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=JKQTPImageTools::LUTSIZE; l>=0; l--) {
              double v=(JKQTPImageTools::LUTSIZE-l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
              double r = 765.0*v-510.0;
              if (r < 0.0)
               r = 0.0;

              double g = std::abs(382.5*v-127.5);

              double b = 255.0*v;
              plut[l]=qRgb(static_cast<int>(r), static_cast<int>(g), static_cast<int>(b));
          }
         }
     }
     lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageBLUEMAGENTAYELLOW) //gnuplot: 30,31,32
    {
     //int*& plut=&lut_bluemagentayellow;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=l/static_cast<double>(JKQTPImageTools::LUTSIZE);
              double r = (v/0.32-0.78125);
              if (r < 0.0) r = 0.0;
              if (r > 1.0) r = 1.0;

              double g = 2.0*v-0.84;
              if (g < 0.0) g = 0.0;
              if (g > 1.0) g = 1.0;

              double b = 4.0*v;
              if (b>1 || b<0) b = -2.0*v+1.84;
              if (b>1 || b<0) b = v/0.08-11.5;
              if (b>1 || b<0) b=1;

              if (b < 0.0) b = 0.0;
              if (b > 1.0) b = 1.0;
              plut[l]=qRgb(static_cast<int>(255.0*r), static_cast<int>(255.0*g), static_cast<int>(255.0*b));
          }
         }
     }
     lut_used=lutstore[palette];
    }
     else if (palette == JKQTPMathImageINVERTED_BLUEMAGENTAYELLOW) //gnuplot: 30,31,32
       {
        //int*& plut=&lut_bluemagentayellow_inverted;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
            if (plut!=nullptr) {
             for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
                 double v=(JKQTPImageTools::LUTSIZE-l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
                 double r = (v/0.32-0.78125);
                 if (r < 0.0) r = 0.0;
                 if (r > 1.0) r = 1.0;

                 double g = 2.0*v-0.84;
                 if (g < 0.0) g = 0.0;
                 if (g > 1.0) g = 1.0;

                 double b = 4.0*v;
                 if (b>1 || b<0) b = -2.0*v+1.84;
                 if (b>1 || b<0) b = v/0.08-11.5;
                 if (b>1 || b<0) b=1;

                 if (b < 0.0) b = 0.0;
                 if (b > 1.0) b = 1.0;
                 plut[l]=qRgb(static_cast<int>(255.0*r), static_cast<int>(255.0*g), static_cast<int>(255.0*b));
             }
            }
        }
        lut_used=lutstore[palette];
       }
    else if (palette == JKQTPMathImageBLUEYELLOW) //gnuplot: 8,9,10
    {
     //int*& plut=&lut_blueyellow;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=l/static_cast<double>(JKQTPImageTools::LUTSIZE);
              double r = sqrt(sqrt(v));
              if (r < 0.0) r = 0.0;
              if (r > 1.0) r = 1.0;

              double g = sin(M_PI/2.0*v);
              if (g < 0.0) g = 0.0;
              if (g > 1.0) g = 1.0;

              double b = cos(M_PI/2.0*v);
              if (b < 0.0) b = 0.0;
              if (b > 1.0) b = 1.0;

              plut[l]=qRgb(static_cast<int>(255.0*r), static_cast<int>(255.0*g), static_cast<int>(255.0*b));
          }
         }
     }
     lut_used=lutstore[palette];

    } else if (palette == JKQTPMathImageINVERTED_BLUEYELLOW) //gnuplot: 8,9,10
    {
     //int*& plut=&lut_blueyellow_inverted;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=(JKQTPImageTools::LUTSIZE-l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
              double r = sqrt(sqrt(v));
              if (r < 0.0) r = 0.0;
              if (r > 1.0) r = 1.0;

              double g = sin(M_PI/2.0*v);
              if (g < 0.0) g = 0.0;
              if (g > 1.0) g = 1.0;

              double b = cos(M_PI/2.0*v);
              if (b < 0.0) b = 0.0;
              if (b > 1.0) b = 1.0;

              plut[l]=qRgb(static_cast<int>(255.0*r), static_cast<int>(255.0*g), static_cast<int>(255.0*b));
          }
         }
     }
     lut_used=lutstore[palette];

    } else if (palette == JKQTPMathImageCYAN)
    {
     //int*& plut=&lut_cyan;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=l/static_cast<double>(JKQTPImageTools::LUTSIZE);
              double r = v*0.5;
              if (r < 0.0) r = 0.0;
              if (r > 1.0) r = 1.0;

              double g = v;
              if (g < 0.0) g = 0.0;
              if (g > 1.0) g = 1.0;

              double b = v;
              if (b < 0.0) b = 0.0;
              if (b > 1.0) b = 1.0;
              plut[l]=qRgb(static_cast<int>(255.0*r), static_cast<int>(255.0*g), static_cast<int>(255.0*b));
          }
         }
     }
     lut_used=lutstore[palette];

    }
     else if (palette == JKQTPMathImageINVERTED_CYAN)
    {
     //int*& plut=&lut_cyan_inverted;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=(JKQTPImageTools::LUTSIZE-l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
              double r = v*0.5;
              if (r < 0.0) r = 0.0;
              if (r > 1.0) r = 1.0;

              double g = v;
              if (g < 0.0) g = 0.0;
              if (g > 1.0) g = 1.0;

              double b = v;
              if (b < 0.0) b = 0.0;
              if (b > 1.0) b = 1.0;
              plut[l]=qRgb(static_cast<int>(255.0*r), static_cast<int>(255.0*g), static_cast<int>(255.0*b));
          }
         }
     }
     lut_used=lutstore[palette];

    }
     else if (palette == JKQTPMathImageTRAFFICLIGHT)
    {
     //int*& plut=&lut_trafficlight;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=l/static_cast<double>(JKQTPImageTools::LUTSIZE);
              double r = (v < 0.5) ? 128.0*sin(M_PI*(2.0*v-0.5))+128.0 : 255.0;
              if (r > 255.0)
               r = 255.0;

              double g = (v < 0.5) ? 512.0*v+128.0 : 512.0-512.0*v;
              if (g > 255.0)
               g = 255.0;
              plut[l]=qRgb(static_cast<int>(r), static_cast<int>(g), 0);
          }
         }
     }
     lut_used=lutstore[palette];

    }
     else if (palette == JKQTPMathImageINVERTED_TRAFFICLIGHT)
    {
     //int*& plut=&lut_trafficlight_inverted;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         if (plut!=nullptr) {
          for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
              double v=(JKQTPImageTools::LUTSIZE-l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
              double r = (v < 0.5) ? 128.0*sin(M_PI*(2.0*v-0.5))+128.0 : 255.0;
              if (r > 255.0)
               r = 255.0;

              double g = (v < 0.5) ? 512.0*v+128.0 : 512.0-512.0*v;
              if (g > 255.0)
               g = 255.0;
              plut[l]=qRgb(static_cast<int>(r), static_cast<int>(g), 0);
          }
         }
     }
     lut_used=lutstore[palette];

    }
    else if (palette == JKQTPMathImageBLUEWHITERED) //gnuplot: 34,13,-34
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         //std::cout<<"!!! creating rainbow lut\n";
         if (plut!=nullptr) {
          QList<QPair<double, QRgb> > lst;
          lst<<qMakePair<double, QRgb>(8.0, 0xFFB2182B);
          lst<<qMakePair<double, QRgb>(7.0, 0xFFD6604D);
          lst<<qMakePair<double, QRgb>(6.0, 0xFFF4A582);
          lst<<qMakePair<double, QRgb>(5.0, 0xFFFDDBC7);
          lst<<qMakePair<double, QRgb>(4.0, 0xFFD1E5F0);
          lst<<qMakePair<double, QRgb>(3.0, 0xFF92C5DE);
          lst<<qMakePair<double, QRgb>(2.0, 0xFF4393C3);
          lst<<qMakePair<double, QRgb>(1.0, 0xFF2166AC);

          JKQTPImagePlot_buildDefinedPaletteLinInterpolate(plut, lst);
         }
     }
     lut_used=lutstore[palette];
    }
    else if (palette == JKQTPMathImageREDWHITEBLUE) //gnuplot: https://github.com/aschn/gnuplot-colorbrewer/blob/master/diverging/RdBu.plt
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         //std::cout<<"!!! creating rainbow lut\n";
         if (plut!=nullptr) {
          QList<QPair<double, QRgb> > lst;
          lst<<qMakePair<double, QRgb>(0.0, 0xFFB2182B);
          lst<<qMakePair<double, QRgb>(1.0, 0xFFD6604D);
          lst<<qMakePair<double, QRgb>(2.0, 0xFFF4A582);
          lst<<qMakePair<double, QRgb>(3.0, 0xFFFDDBC7);
          lst<<qMakePair<double, QRgb>(4.0, 0xFFD1E5F0);
          lst<<qMakePair<double, QRgb>(5.0, 0xFF92C5DE);
          lst<<qMakePair<double, QRgb>(6.0, 0xFF4393C3);
          lst<<qMakePair<double, QRgb>(7.0, 0xFF2166AC);

          JKQTPImagePlot_buildDefinedPaletteLinInterpolate(plut, lst);
         }
     }
     lut_used=lutstore[palette];



    } else if (palette == JKQTPMathImageBLACKBLUEREDYELLOW) //gnuplot: 7,5,15
     {
//          int*& plut=&lut_BLACKBLUEREDYELLOW;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
               double v=(l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
               double r = 255.0*qBound(0.0,sqrt(v),1.0);
               double g = 255.0*qBound(0.0,v*v*v,1.0);
               double b = 255.0*qBound(0.0,sin(2.0*M_PI*v),1.0);

               plut[l]=qRgb(static_cast<int>(r), static_cast<int>(g), static_cast<int>(b));
              }
          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageGREENREDVIOLET) //gnuplot: 3,11,6
    {
         //int*& plut=&lut_GREENREDVIOLET;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
               double v=(l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
               double r = 255.0*qBound(0.0,v,1.0);
               double g = 255.0*qBound(0.0,fabs(v-0.5),1.0);
               double b = 255.0*qBound(0.0,v*v*v*v,1.0);

               plut[l]=qRgb(static_cast<int>(r), static_cast<int>(g), static_cast<int>(b));
              }
          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageBLACKBLUEWHITEYELLOWWHITE) //gnuplot: 30,31,32
    {
         //int*& plut=&lut_BLACKBLUEVIOLETYELLOWWHITE;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
               double v=(l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
               double r = 255.0*qBound(0.0,v/0.32-0.78125,1.0);
               double g = 255.0*qBound(0.0,v/0.32-0.78125,1.0);
               double b = 255.0*qBound(0.0,(v<0.25)?4*v:(v<0.42)?1.0:(v<0.92)?-2.0*v+1.84:v/0.08-11.5,1.0);

               plut[l]=qRgb(static_cast<int>(r), static_cast<int>(g), static_cast<int>(b));
              }
          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageWHITEYELLOWWHITEBLUEBLACK) //gnuplot: 30,31,32
    {
         //int*& plut=&lut_WHITEYELLOWVIOLETBLUEBLACK;
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
               double v=(JKQTPImageTools::LUTSIZE-l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
               double r = 255.0*qBound(0.0,v/0.32-0.78125,1.0);
               double g = 255.0*qBound(0.0,v/0.32-0.78125,1.0);
               double b = 255.0*qBound(0.0,(v<0.25)?4*v:(v<0.42)?1.0:(v<0.92)?-2.0*v+1.84:v/0.08-11.5,1.0);

               plut[l]=qRgb(static_cast<int>(r), static_cast<int>(g), static_cast<int>(b));
              }
          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageBR_GR) //https://github.com/aschn/gnuplot-colorbrewer/blob/master/diverging/BrBG.plt
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              QList<QPair<double, QRgb> > lst;
              lst<<qMakePair<double, QRgb>(0.0, qRgb(0x8C, 0x51, 0x0A));
              lst<<qMakePair<double, QRgb>(1.0, qRgb(0xBF, 0x81, 0x2D));
              lst<<qMakePair<double, QRgb>(2.0, qRgb(0xDF, 0xC2, 0x7D));
              lst<<qMakePair<double, QRgb>(3.0, 0xFFF6E8C3);
              lst<<qMakePair<double, QRgb>(4.0, 0xFFC7EAE5);
              lst<<qMakePair<double, QRgb>(5.0, 0xFF80CDC1);
              lst<<qMakePair<double, QRgb>(6.0, 0xFF35978F);
              lst<<qMakePair<double, QRgb>(7.0, 0xFF01665E);

              JKQTPImagePlot_buildDefinedPaletteLinInterpolate(plut, lst);
          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageBR_GR_STEP) //https://github.com/aschn/gnuplot-colorbrewer/blob/master/diverging/BrBG.plt
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              QList<QPair<double, QRgb> > lst;
              lst<<qMakePair<double, QRgb>(0.0, qRgb(0x8C, 0x51, 0x0A));
              lst<<qMakePair<double, QRgb>(1.0, qRgb(0xBF, 0x81, 0x2D));
              lst<<qMakePair<double, QRgb>(2.0, qRgb(0xDF, 0xC2, 0x7D));
              lst<<qMakePair<double, QRgb>(3.0, 0xFFF6E8C3);
              lst<<qMakePair<double, QRgb>(4.0, 0xFFC7EAE5);
              lst<<qMakePair<double, QRgb>(5.0, 0xFF80CDC1);
              lst<<qMakePair<double, QRgb>(6.0, 0xFF35978F);
              lst<<qMakePair<double, QRgb>(7.0, 0xFF01665E);

              JKQTPImagePlot_buildDefinedPalette(plut, lst);
          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImagePU_OR) //https://github.com/aschn/gnuplot-colorbrewer/blob/master/diverging/PuOr.plt
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              QList<QPair<double, QRgb> > lst;
              lst<<qMakePair<double, QRgb>(0.0, 0xFFB35806);
              lst<<qMakePair<double, QRgb>(1.0, 0xFFE08214);
              lst<<qMakePair<double, QRgb>(2.0, 0xFFFDB863);
              lst<<qMakePair<double, QRgb>(3.0, 0xFFFEE0B6);
              lst<<qMakePair<double, QRgb>(4.0, 0xFFF7F7F7);
              lst<<qMakePair<double, QRgb>(5.0, 0xFFD8DAEB);
              lst<<qMakePair<double, QRgb>(6.0, 0xFFB2ABD2);
              lst<<qMakePair<double, QRgb>(7.0, 0xFF8073AC);
              lst<<qMakePair<double, QRgb>(8.0, 0xFF542788);

              JKQTPImagePlot_buildDefinedPaletteLinInterpolate(plut, lst);
          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImagePU_OR_STEP) //https://github.com/aschn/gnuplot-colorbrewer/blob/master/diverging/PuOr.plt
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              QList<QPair<double, QRgb> > lst;
              lst<<qMakePair<double, QRgb>(0.0, 0xFFB35806);
              lst<<qMakePair<double, QRgb>(1.0, 0xFFE08214);
              lst<<qMakePair<double, QRgb>(2.0, 0xFFFDB863);
              lst<<qMakePair<double, QRgb>(3.0, 0xFFFEE0B6);
              lst<<qMakePair<double, QRgb>(4.0, 0xFFF7F7F7);
              lst<<qMakePair<double, QRgb>(5.0, 0xFFD8DAEB);
              lst<<qMakePair<double, QRgb>(6.0, 0xFFB2ABD2);
              lst<<qMakePair<double, QRgb>(7.0, 0xFF8073AC);
              lst<<qMakePair<double, QRgb>(8.0, 0xFF542788);

              JKQTPImagePlot_buildDefinedPalette(plut, lst);
          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageYL_GN_BU) //http://colorbrewer2.org/
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              QList<QPair<double, QRgb> > lst;
              lst<<qMakePair<double, QRgb>(0.0, 0xFFFFFFD9);
              lst<<qMakePair<double, QRgb>(1.0, 0xFFEDF8B1);
              lst<<qMakePair<double, QRgb>(2.0, 0xFFC7E9B4);
              lst<<qMakePair<double, QRgb>(3.0, 0xFF7FCDBB);
              lst<<qMakePair<double, QRgb>(4.0, 0xFF41B6C4);
              lst<<qMakePair<double, QRgb>(5.0, 0xFF1D91C0);
              lst<<qMakePair<double, QRgb>(6.0, 0xFF225EA8);
              lst<<qMakePair<double, QRgb>(7.0, 0xFF253494);
              lst<<qMakePair<double, QRgb>(8.0, 0xFF081D58);

              JKQTPImagePlot_buildDefinedPaletteLinInterpolate(plut, lst);
          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageYL_GN_BU_STEP) //http://colorbrewer2.org/
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              QList<QPair<double, QRgb> > lst;
              lst<<qMakePair<double, QRgb>(0.0, 0xFFFFFFD9);
              lst<<qMakePair<double, QRgb>(1.0, 0xFFEDF8B1);
              lst<<qMakePair<double, QRgb>(2.0, 0xFFC7E9B4);
              lst<<qMakePair<double, QRgb>(3.0, 0xFF7FCDBB);
              lst<<qMakePair<double, QRgb>(4.0, 0xFF41B6C4);
              lst<<qMakePair<double, QRgb>(5.0, 0xFF1D91C0);
              lst<<qMakePair<double, QRgb>(6.0, 0xFF225EA8);
              lst<<qMakePair<double, QRgb>(7.0, 0xFF253494);
              lst<<qMakePair<double, QRgb>(8.0, 0xFF081D58);

              JKQTPImagePlot_buildDefinedPalette(plut, lst);
          }
         }
         lut_used=lutstore[palette];

    } else if (palette == JKQTPMathImageGN_BU) //http://colorbrewer2.org/
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              QList<QPair<double, QRgb> > lst;
              lst<<qMakePair<double, QRgb>(0.0, 0xFFF7FCF0);
              lst<<qMakePair<double, QRgb>(1.0, 0xFFE0F3DB);
              lst<<qMakePair<double, QRgb>(2.0, 0xFFCCEBC5);
              lst<<qMakePair<double, QRgb>(3.0, 0xFFA8DDB5);
              lst<<qMakePair<double, QRgb>(4.0, 0xFF7BCCC4);
              lst<<qMakePair<double, QRgb>(5.0, 0xFF4EB3D3);
              lst<<qMakePair<double, QRgb>(6.0, 0xFF2B8CBE);
              lst<<qMakePair<double, QRgb>(7.0, 0xFF0868AC);
              lst<<qMakePair<double, QRgb>(8.0, 0xFF084081);

              JKQTPImagePlot_buildDefinedPaletteLinInterpolate(plut, lst);
          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageGN_BU_STEP) //http://colorbrewer2.org/
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              QList<QPair<double, QRgb> > lst;
              lst<<qMakePair<double, QRgb>(0.0, 0xFFF7FCF0);
              lst<<qMakePair<double, QRgb>(1.0, 0xFFE0F3DB);
              lst<<qMakePair<double, QRgb>(2.0, 0xFFCCEBC5);
              lst<<qMakePair<double, QRgb>(3.0, 0xFFA8DDB5);
              lst<<qMakePair<double, QRgb>(4.0, 0xFF7BCCC4);
              lst<<qMakePair<double, QRgb>(5.0, 0xFF4EB3D3);
              lst<<qMakePair<double, QRgb>(6.0, 0xFF2B8CBE);
              lst<<qMakePair<double, QRgb>(7.0, 0xFF0868AC);
              lst<<qMakePair<double, QRgb>(8.0, 0xFF084081);

              JKQTPImagePlot_buildDefinedPalette(plut, lst);
          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageBU_GN) //http://colorbrewer2.org/
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              QList<QPair<double, QRgb> > lst;
              lst<<qMakePair<double, QRgb>(8.0, 0xFFF7FCF0);
              lst<<qMakePair<double, QRgb>(7.0, 0xFFE0F3DB);
              lst<<qMakePair<double, QRgb>(6.0, 0xFFCCEBC5);
              lst<<qMakePair<double, QRgb>(5.0, 0xFFA8DDB5);
              lst<<qMakePair<double, QRgb>(4.0, 0xFF7BCCC4);
              lst<<qMakePair<double, QRgb>(3.0, 0xFF4EB3D3);
              lst<<qMakePair<double, QRgb>(2.0, 0xFF2B8CBE);
              lst<<qMakePair<double, QRgb>(1.0, 0xFF0868AC);
              lst<<qMakePair<double, QRgb>(0.0, 0xFF084081);

              JKQTPImagePlot_buildDefinedPaletteLinInterpolate(plut, lst);
          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageBU_GN_STEP) //http://colorbrewer2.org/
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              QList<QPair<double, QRgb> > lst;
              lst<<qMakePair<double, QRgb>(8.0, 0xFFF7FCF0);
              lst<<qMakePair<double, QRgb>(7.0, 0xFFE0F3DB);
              lst<<qMakePair<double, QRgb>(6.0, 0xFFCCEBC5);
              lst<<qMakePair<double, QRgb>(5.0, 0xFFA8DDB5);
              lst<<qMakePair<double, QRgb>(4.0, 0xFF7BCCC4);
              lst<<qMakePair<double, QRgb>(3.0, 0xFF4EB3D3);
              lst<<qMakePair<double, QRgb>(2.0, 0xFF2B8CBE);
              lst<<qMakePair<double, QRgb>(1.0, 0xFF0868AC);
              lst<<qMakePair<double, QRgb>(0.0, 0xFF084081);

              JKQTPImagePlot_buildDefinedPalette(plut, lst);
          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageINVERTED_MAGENTA) // from CMYK model
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
               double v=(l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
               plut[l]=QColor::fromCmykF(0,1,0,v).rgba();
              }

          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageMAGENTA) // from CMYK model
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
               double v=(JKQTPImageTools::LUTSIZE-l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
               plut[l]=QColor::fromCmykF(0,1,0,v).rgba();
              }

          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageINVERTED_YELLOW) // from CMYK model
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
               double v=(l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
               plut[l]=QColor::fromCmykF(0,0,1,v).rgba();
              }

          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageYELLOW) // from CMYK model
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
               double v=(JKQTPImageTools::LUTSIZE-l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
               plut[l]=QColor::fromCmykF(0,0,1,v).rgba();
              }

          }
         }
         lut_used=lutstore[palette];

    } else if (palette == JKQTPMathImageINVERTED_MAGENTAWHITE) // from CMYK model
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
               double v=(l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
               plut[l]=QColor::fromCmykF(0,v,0,0).rgba();
              }

          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageMAGENTAWHITE) // from CMYK model
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
               double v=(JKQTPImageTools::LUTSIZE-l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
               plut[l]=QColor::fromCmykF(0,v,0,0).rgba();
              }

          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageINVERTED_YELLOWWHITE) // from CMYK model
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
               double v=(l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
               plut[l]=QColor::fromCmykF(0,0,v,0).rgba();
              }

          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageYELLOWWHITE) // from CMYK model
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
               double v=(JKQTPImageTools::LUTSIZE-l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
               plut[l]=QColor::fromCmykF(0,0,v,0).rgba();
              }

          }
         }
         lut_used=lutstore[palette];

    } else if (palette == JKQTPMathImageINVERTED_CYANWHITE) // from CMYK model
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
               double v=(l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
               plut[l]=QColor::fromCmykF(v,0,0,0).rgba();
              }

          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageCYANWHITE) // from CMYK model
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              for (int l=0; l<=JKQTPImageTools::LUTSIZE; l++) {
               double v=(JKQTPImageTools::LUTSIZE-l)/static_cast<double>(JKQTPImageTools::LUTSIZE);
               plut[l]=QColor::fromCmykF(v,0,0,0).rgba();
              }

          }
         }
         lut_used=lutstore[palette];
    } else if (palette == JKQTPMathImageBlueGreenRed)
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              QList<QPair<double, QRgb> > lst;
              lst<<qMakePair<double, QRgb>(1.0, QColor("blue").rgb());
              lst<<qMakePair<double, QRgb>(2.0, QColor("green").rgb());
              lst<<qMakePair<double, QRgb>(3.0, QColor("red").rgb());

              JKQTPImagePlot_buildDefinedPaletteLinInterpolate(plut, lst);
          }
         }
         lut_used=lutstore[palette];

    } else if (palette == JKQTPMathImageRedGreenBlue)
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              QList<QPair<double, QRgb> > lst;
              lst<<qMakePair<double, QRgb>(1.0, QColor("red").rgb());
              lst<<qMakePair<double, QRgb>(2.0, QColor("green").rgb());
              lst<<qMakePair<double, QRgb>(3.0, QColor("blue").rgb());

              JKQTPImagePlot_buildDefinedPaletteLinInterpolate(plut, lst);
          }
         }
         lut_used=lutstore[palette];


    } else if (palette == JKQTPMathImagePREDEFINED_PALETTES_COUNT)
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              QList<QPair<double, QRgb> > lst;
              lst<<qMakePair<double, QRgb>(1.0, QColor("magenta").rgb());
              lst<<qMakePair<double, QRgb>(2.0, QColor("yellow").rgb());

              JKQTPImagePlot_buildDefinedPaletteLinInterpolate(plut, lst);
          }
         }
         lut_used=lutstore[palette];

    } else if (palette == JKQTPMathImagePREDEFINED_PALETTES_COUNT+1)
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              QList<QPair<double, QRgb> > lst;
              lst<<qMakePair<double, QRgb>(1.0, QColor("yellow").rgb());
              lst<<qMakePair<double, QRgb>(2.0, QColor("magenta").rgb());

              JKQTPImagePlot_buildDefinedPaletteLinInterpolate(plut, lst);
          }
         }
         lut_used=lutstore[palette];

    } else if (palette == JKQTPMathImagePREDEFINED_PALETTES_COUNT+2)
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              QList<QPair<double, QRgb> > lst;
              lst<<qMakePair<double, QRgb>(1.0, QColor("red").rgb());
              lst<<qMakePair<double, QRgb>(2.0, QColor("blue").rgb());

              JKQTPImagePlot_buildDefinedPaletteLinInterpolate(plut, lst);
          }
         }
         lut_used=lutstore[palette];

    } else if (palette == JKQTPMathImagePREDEFINED_PALETTES_COUNT+3)
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         //std::cout<<"!!! creating rainbow lut\n";
         if (plut!=nullptr) {
          QList<QPair<double, QRgb> > lst;
          lst<<qMakePair<double, QRgb>(1.0, QColor("blue").rgb());
          lst<<qMakePair<double, QRgb>(2.0, QColor("red").rgb());

          JKQTPImagePlot_buildDefinedPaletteLinInterpolate(plut, lst);
         }
     }
     lut_used=lutstore[palette];


         /*
    } else if (palette == JKQTPMathImagePU_OR) //https://github.com/aschn/gnuplot-colorbrewer/blob/master/diverging/PuOr.plt
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
          //std::cout<<"!!! creating rainbow lut\n";
          if (plut!=nullptr) {
              QList<QPair<double, QRgb> > lst;
              lst<<qMakePair<double, QRgb>(0.0, 0xFF);
              lst<<qMakePair<double, QRgb>(1.0, 0xFF);
              lst<<qMakePair<double, QRgb>(2.0, 0xFF);
              lst<<qMakePair<double, QRgb>(3.0, 0xFF);
              lst<<qMakePair<double, QRgb>(4.0, 0xFF);
              lst<<qMakePair<double, QRgb>(5.0, 0xFF);
              lst<<qMakePair<double, QRgb>(6.0, 0xFF);
              lst<<qMakePair<double, QRgb>(7.0, 0xFF);
              lst<<qMakePair<double, QRgb>(8.0, 0xFF);

              JKQTPImagePlot_buildDefinedPalette(plut, lst);
          }
         }
         lut_used=lutstore[palette];

  */

    } else if (JKQTPImagePlot_lutsFromFiles.contains(palette)) // LUTs read from files in JKQTPImagePlot_lutsFromFiles
    {
     if (lutstore[palette]==nullptr) {
         lutstore[palette]=(int*)malloc((JKQTPImageTools::LUTSIZE+2)*sizeof(int));
         int* plut=lutstore[palette];
         //std::cout<<"!!! creating rainbow lut\n";
         if (plut!=nullptr) {
          if (JKQTPImagePlot_lutsFromFiles[palette].size()>220) JKQTPImagePlot_buildDefinedPalette(plut, JKQTPImagePlot_lutsFromFiles[palette]);
          else JKQTPImagePlot_buildDefinedPaletteLinInterpolate(plut, JKQTPImagePlot_lutsFromFiles[palette]);
         }
     }
     lut_used=lutstore[palette];
    }
    return lut_used;
}


QImage JKQTPMathImageGetPaletteImage(int i, int width)
{
    return JKQTPMathImageGetPaletteImage(i, width, 1);
}

QImage JKQTPMathImageGetPaletteImage(int i, int width, int height)
{
    QImage img;
    double* pic=static_cast<double*>(malloc(width*height*sizeof(double)));
    for (int j=0; j<width*height; j++) {
     pic[j]=j%width;
    }
    JKQTPImagePlot_array2image<double>(pic, width, height, img, (JKQTPMathImageColorPalette)i, Qt::black, Qt::black);
    free(pic);
    return img;
}

QImage JKQTPMathImageGetPaletteImage(JKQTPMathImageColorPalette palette, int width)
{
    return JKQTPMathImageGetPaletteImage(static_cast<int>(palette), width);
}



QImage JKQTPMathImageGetPaletteImage(JKQTPMathImageColorPalette palette, int width, int height)
{
    return JKQTPMathImageGetPaletteImage(static_cast<int>(palette), width, height);
}

QImage JKQTPMathImageGetPaletteImage(int *lut, int lut_size, int width)
{
    QImage img;
    double* pic=static_cast<double*>(malloc(width*sizeof(double)));
    for (int j=0; j<width; j++) {
     pic[j]=j;
    }
    JKQTPImagePlot_array2image<double>(pic, width, 1, img, lut, lut_size, Qt::black, Qt::black);
    free(pic);
    return img;
}



QIcon JKQTPMathImageGetPaletteIcon(int i)  {
    QImage img=JKQTPMathImageGetPaletteImage(i, JKQTPImageTools::PALETTE_ICON_WIDTH);
    QPixmap pix(JKQTPImageTools::PALETTE_ICON_WIDTH,8);
    QRect r(0,0,JKQTPImageTools::PALETTE_ICON_WIDTH-1,7);
    QPainter p(&pix);
    p.drawImage(r, img);
    p.setPen(QPen(QColor("black")));
    p.drawRect(r);
    p.end();
    return QIcon(pix);

}

QIcon JKQTPMathImageGetPaletteIcon(JKQTPMathImageColorPalette palette)  {
    return JKQTPMathImageGetPaletteIcon(static_cast<int>(palette));
}







QImage JKQTPMathImageGetAlphaPaletteImage(int *lut, int lut_size, int width, int height)
{
    QImage img;
    if (!lut || lut_size<=0 || width<=0 || height<=0) return img;
    double* pic=static_cast<double*>(malloc(width*sizeof(double)));
    for (int j=0; j<width; j++) {
     pic[j]=j;
    }
    QVector<int> l1(lut_size, QColor("white").rgba());
    for (int i=0; i<l1.size(); i++) {
     l1[i]=qRgba(255,255,255, lut[i]);
    }
    QImage img1(8,8,QImage::Format_ARGB32);
    img1.fill(Qt::white);
    for (int y=0; y<img1.height()/2; y++) {
     QRgb* sl1=(QRgb*)img1.scanLine(y);
     QRgb* sl2=(QRgb*)img1.scanLine(y+img1.height()/2);
     for (int x=0; x<img1.width()/2; x++) {
         sl1[x]=QColor(Qt::black).rgba();
         sl2[x+img1.width()/2]=QColor(Qt::black).rgba();
     }
    }
    JKQTPImagePlot_array2image<double>(pic, width, 1, img, l1.data(), lut_size, Qt::black, Qt::black);
    free(pic);
    QImage img2(width, height, QImage::Format_ARGB32);
    QPainter p(&img2);
    QBrush b;
    b.setStyle(Qt::TexturePattern);
    b.setTextureImage(img1);
    p.setBrush(b);
    p.fillRect(QRect(0,0,width,height), b);
    p.drawImage(QRect(0,0,width,height), img);
    p.end();
    return img2;
}





QString JKQTPMathImageColorPalette2String(JKQTPMathImageColorPalette p)
{
    switch (p) {
        case JKQTPMathImageGRAY: return QString("gray");
        case JKQTPMathImageINVERTEDGRAY: return QString("invgray");
        case JKQTPMathImageRED: return QString("red");
        case JKQTPMathImageINVERTEDRED: return QString("invred");
        case JKQTPMathImageGREEN: return QString("green");
        case JKQTPMathImageINVERTEDGREEN: return QString("invgreen");
        case JKQTPMathImageBLUE: return QString("blue");
        case JKQTPMathImageINVERTEDBLUE: return QString("invblue");
        case JKQTPMathImageCYAN: return QString("cyan");
        case JKQTPMathImageINVERTED_CYAN: return QString("invcyan");
        case JKQTPMathImageYELLOW: return QString("yellow");
        case JKQTPMathImageINVERTED_YELLOW: return QString("invyellow");
        case JKQTPMathImageMAGENTA: return QString("magenta");
        case JKQTPMathImageINVERTED_MAGENTA: return QString("invmagenta");


        case JKQTPMathImageMATLAB: return QString("Matlab");
        case JKQTPMathImageINVERTED_MATLAB: return QString("invMatlab");
        case JKQTPMathImageRYGB: return QString("RYGB");
        case JKQTPMathImageINVERTED_RYGB: return QString("invRYGB");
        case JKQTPMathImageHSV: return QString("HSV");
        case JKQTPMathImageINVERTED_HSV: return QString("invHSV");
        case JKQTPMathImageRAINBOW: return QString("rainbow");
        case JKQTPMathImageINVERTED_RAINBOW: return QString("invrainbow");
        case JKQTPMathImageHOT: return QString("AFMhot");
        case JKQTPMathImageINVERTED_HOT: return QString("invAFMhot");
        case JKQTPMathImageOCEAN: return QString("ocean");
        case JKQTPMathImageINVERTED_OCEAN: return QString("invocean");
        case JKQTPMathImageTRAFFICLIGHT: return QString("trafficlight");
        case JKQTPMathImageINVERTED_TRAFFICLIGHT: return QString("invtrafficlight");
        case JKQTPMathImageBLUEMAGENTAYELLOW: return QString("BlMaYe");
        case JKQTPMathImageINVERTED_BLUEMAGENTAYELLOW: return QString("YeMaBl");
        case JKQTPMathImageBLUEYELLOW: return QString("BlYe");
        case JKQTPMathImageINVERTED_BLUEYELLOW: return QString("YeBl");

        case JKQTPMathImageBLUEWHITERED: return QString("bluewhitered");
        case JKQTPMathImageREDWHITEBLUE: return QString("redwhiteblue");

        case JKQTPMathImageBLACKBLUEREDYELLOW: return QString("BBlRdYe");
        case JKQTPMathImageGREENREDVIOLET: return QString("GnRdVi");
        case JKQTPMathImageBLACKBLUEWHITEYELLOWWHITE: return QString("BWprint");
        case JKQTPMathImageWHITEYELLOWWHITEBLUEBLACK: return QString("invBWprint");
        case JKQTPMathImageBR_GR: return QString("BrBG");
        case JKQTPMathImagePU_OR: return QString("PuOr");
        case JKQTPMathImageGN_BU: return QString("greenblue");
        case JKQTPMathImageBU_GN: return QString("bluegreen");
        case JKQTPMathImageYL_GN_BU: return QString("YeGnBu");

        case JKQTPMathImageBR_GR_STEP: return QString("stepsBrBG");
        case JKQTPMathImagePU_OR_STEP: return QString("stepsPuOr");
        case JKQTPMathImageGN_BU_STEP: return QString("stepsGnBl");
        case JKQTPMathImageBU_GN_STEP: return QString("stepsBlGn");
        case JKQTPMathImageYL_GN_BU_STEP: return QString("stepsYeGnBu");


        case JKQTPMathImageCYANWHITE: return QString("cyanwhite");
        case JKQTPMathImageINVERTED_CYANWHITE: return QString("whitecyan");
        case JKQTPMathImageYELLOWWHITE: return QString("yellowwhite");
        case JKQTPMathImageINVERTED_YELLOWWHITE: return QString("whiteyellow");
        case JKQTPMathImageMAGENTAWHITE: return QString("magentawhite");
        case JKQTPMathImageINVERTED_MAGENTAWHITE: return QString("whitemagenta");
        case JKQTPMathImageBlueGreenRed: return QString("bluegreenred");
        case JKQTPMathImageRedGreenBlue: return QString("redgreenblue");
        case JKQTPMathImageALPHA: return QString("alpha");
        case JKQTPMathImageINVERTED_ALPHA: return QString("invAlpha");
        default:
             return QString::number(static_cast<int>(p));
    }
}

JKQTPMathImageColorPalette String2JKQTPMathImageColorPalette(const QString &p)
{
    if (QString::compare(p, "gray", Qt::CaseInsensitive)==0) return JKQTPMathImageGRAY;
    if (QString::compare(p, "invgray", Qt::CaseInsensitive)==0) return JKQTPMathImageINVERTEDGRAY;
    if (QString::compare(p, "red", Qt::CaseInsensitive)==0) return JKQTPMathImageRED;
    if (QString::compare(p, "invred", Qt::CaseInsensitive)==0) return JKQTPMathImageINVERTEDRED;
    if (QString::compare(p, "green", Qt::CaseInsensitive)==0) return JKQTPMathImageGREEN;
    if (QString::compare(p, "invgreen", Qt::CaseInsensitive)==0) return JKQTPMathImageINVERTEDGREEN;
    if (QString::compare(p, "blue", Qt::CaseInsensitive)==0) return JKQTPMathImageBLUE;
    if (QString::compare(p, "invblue", Qt::CaseInsensitive)==0) return JKQTPMathImageINVERTEDBLUE;
    if (QString::compare(p, "cyan", Qt::CaseInsensitive)==0) return JKQTPMathImageCYAN;
    if (QString::compare(p, "invcyan", Qt::CaseInsensitive)==0) return JKQTPMathImageINVERTED_CYAN;
    if (QString::compare(p, "yellow", Qt::CaseInsensitive)==0) return JKQTPMathImageYELLOW;
    if (QString::compare(p, "invyellow", Qt::CaseInsensitive)==0) return JKQTPMathImageINVERTED_YELLOW;
    if (QString::compare(p, "magenta", Qt::CaseInsensitive)==0) return JKQTPMathImageMAGENTA;
    if (QString::compare(p, "invmagenta", Qt::CaseInsensitive)==0) return JKQTPMathImageINVERTED_MAGENTA;


    if (QString::compare(p, "Matlab", Qt::CaseInsensitive)==0) return JKQTPMathImageMATLAB;
    if (QString::compare(p, "invMatlab", Qt::CaseInsensitive)==0) return JKQTPMathImageINVERTED_MATLAB;
    if (QString::compare(p, "RYGB", Qt::CaseInsensitive)==0) return JKQTPMathImageRYGB;
    if (QString::compare(p, "invRYGB", Qt::CaseInsensitive)==0) return JKQTPMathImageINVERTED_RYGB;
    if (QString::compare(p, "HSV", Qt::CaseInsensitive)==0) return JKQTPMathImageHSV;
    if (QString::compare(p, "invHSV", Qt::CaseInsensitive)==0) return JKQTPMathImageINVERTED_HSV;
    if (QString::compare(p, "rainbow", Qt::CaseInsensitive)==0) return JKQTPMathImageRAINBOW;
    if (QString::compare(p, "invrainbow", Qt::CaseInsensitive)==0) return JKQTPMathImageINVERTED_RAINBOW;
    if (QString::compare(p, "AFMhot", Qt::CaseInsensitive)==0) return JKQTPMathImageHOT;
    if (QString::compare(p, "invAFMhot", Qt::CaseInsensitive)==0) return JKQTPMathImageINVERTED_HOT;
    if (QString::compare(p, "ocean", Qt::CaseInsensitive)==0) return JKQTPMathImageOCEAN;
    if (QString::compare(p, "invocean", Qt::CaseInsensitive)==0) return JKQTPMathImageINVERTED_OCEAN;
    if (QString::compare(p, "trafficlight", Qt::CaseInsensitive)==0) return JKQTPMathImageTRAFFICLIGHT;
    if (QString::compare(p, "invtrafficlight", Qt::CaseInsensitive)==0) return JKQTPMathImageINVERTED_TRAFFICLIGHT;
    if (QString::compare(p, "BlMaYe", Qt::CaseInsensitive)==0) return JKQTPMathImageBLUEMAGENTAYELLOW;
    if (QString::compare(p, "YeMaBl", Qt::CaseInsensitive)==0) return JKQTPMathImageINVERTED_BLUEMAGENTAYELLOW;
    if (QString::compare(p, "BlYe", Qt::CaseInsensitive)==0) return JKQTPMathImageBLUEYELLOW;
    if (QString::compare(p, "YeBl", Qt::CaseInsensitive)==0) return JKQTPMathImageINVERTED_BLUEYELLOW;

    if (QString::compare(p, "bluewhitered", Qt::CaseInsensitive)==0) return JKQTPMathImageBLUEWHITERED;
    if (QString::compare(p, "redwhiteblue", Qt::CaseInsensitive)==0) return JKQTPMathImageREDWHITEBLUE;

    if (QString::compare(p, "BBlRdYe", Qt::CaseInsensitive)==0) return JKQTPMathImageBLACKBLUEREDYELLOW;
    if (QString::compare(p, "GnRdVi", Qt::CaseInsensitive)==0) return JKQTPMathImageGREENREDVIOLET;
    if (QString::compare(p, "BWprint", Qt::CaseInsensitive)==0) return JKQTPMathImageBLACKBLUEWHITEYELLOWWHITE;
    if (QString::compare(p, "invBWprint", Qt::CaseInsensitive)==0) return JKQTPMathImageWHITEYELLOWWHITEBLUEBLACK;
    if (QString::compare(p, "BrBG", Qt::CaseInsensitive)==0) return JKQTPMathImageBR_GR;
    if (QString::compare(p, "PuOr", Qt::CaseInsensitive)==0) return JKQTPMathImagePU_OR;
    if (QString::compare(p, "greenblue", Qt::CaseInsensitive)==0) return JKQTPMathImageGN_BU;
    if (QString::compare(p, "bluegreen", Qt::CaseInsensitive)==0) return JKQTPMathImageBU_GN;
    if (QString::compare(p, "YeGnBu", Qt::CaseInsensitive)==0) return JKQTPMathImageYL_GN_BU;

    if (QString::compare(p, "stepsBrBG", Qt::CaseInsensitive)==0) return JKQTPMathImageBR_GR_STEP;
    if (QString::compare(p, "stepsPuOr", Qt::CaseInsensitive)==0) return JKQTPMathImagePU_OR_STEP;
    if (QString::compare(p, "stepsGnBl", Qt::CaseInsensitive)==0) return JKQTPMathImageGN_BU_STEP;
    if (QString::compare(p, "stepsBlGn", Qt::CaseInsensitive)==0) return JKQTPMathImageBU_GN_STEP;
    if (QString::compare(p, "stepsYeGnBu", Qt::CaseInsensitive)==0) return JKQTPMathImageYL_GN_BU_STEP;


    if (QString::compare(p, "cyanwhite", Qt::CaseInsensitive)==0) return JKQTPMathImageCYANWHITE;
    if (QString::compare(p, "whitecyan", Qt::CaseInsensitive)==0) return JKQTPMathImageINVERTED_CYANWHITE;
    if (QString::compare(p, "yellowwhite", Qt::CaseInsensitive)==0) return JKQTPMathImageYELLOWWHITE;
    if (QString::compare(p, "whiteyellow", Qt::CaseInsensitive)==0) return JKQTPMathImageINVERTED_YELLOWWHITE;
    if (QString::compare(p, "magentawhite", Qt::CaseInsensitive)==0) return JKQTPMathImageMAGENTAWHITE;
    if (QString::compare(p, "whitemagenta", Qt::CaseInsensitive)==0) return JKQTPMathImageINVERTED_MAGENTAWHITE;
    if (QString::compare(p, "bluegreenred", Qt::CaseInsensitive)==0) return JKQTPMathImageBlueGreenRed;
    if (QString::compare(p, "redgreenblue", Qt::CaseInsensitive)==0) return JKQTPMathImageRedGreenBlue;
    if (QString::compare(p, "alpha", Qt::CaseInsensitive)==0) return JKQTPMathImageALPHA;
    if (QString::compare(p, "invAlpha", Qt::CaseInsensitive)==0) return JKQTPMathImageINVERTED_ALPHA;

    bool ok=false;
    JKQTPMathImageColorPalette res = static_cast<JKQTPMathImageColorPalette>(p.toInt(&ok));
    if (!ok) return JKQTPMathImageMATLAB;
    else return res;

}