mirrors/JKQtPlotter
1
0
Fork 0
mirror of https://github.com/jkriege2/JKQtPlotter.git synced 2024-12-26 10:31:39 +08:00
Code Issues Actions Packages Projects Releases Wiki Activity
0d71ddc68a
JKQtPlotter/lib/jkqtplottertools/jkqtpimagetools.cpp

2165 lines
91 KiB
C++
Raw Blame History

/*
Copyright (c) 2008-2018 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 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 "jkqtplottertools/jkqtpimagetools.h"
#include <QPainter>
#include <QApplication>
#include <QDir>
#include <QFileInfo>
#ifndef NO_JKQTPLOTTER
#include "jkqtplotter/jkqtpbaseelements.h"
#include "jkqtplotter/jkqtpbaseplotter.h"
#endif
#ifdef QT_XML_LIB
# include <QtXml/QtXml>
#endif
#define JKQTP_PALETTE_ICON_WIDTH 64
#define JKQTP_PALETTE_IMAGEICON_HEIGHT 64
QList<int*> 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(JKQTPimagePlot_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<=JKQTPimagePlot_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,(int)round(255*r), 255), qBound(0,(int)round(255*g), 255), qBound(0,(int)round(255*b), 255)));
} else {
pal<<qMakePair(x, qRgb(qBound(0,(int)round(r), 255), qBound(0,(int)round(g), 255), qBound(0,(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,(int)round(255*r), 255), qBound(0,(int)round(255*g), 255), qBound(0,(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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=l/(double)(JKQTPimagePlot_LUTSIZE);
plut[l]=qRgb(static_cast<int>(255.0*v), 0, 0);
}
plut[JKQTPimagePlot_LUTSIZE+1]=plut[JKQTPimagePlot_LUTSIZE];
}
}
lut_used=lutstore[palette];
}
else if (palette == JKQTPMathImageGREEN)
{
//int*& plut=&lut_green;
if (lutstore[palette]==nullptr) {
lutstore[palette]=(int*)malloc((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=l/(double)(JKQTPimagePlot_LUTSIZE);
plut[l]=qRgb(0, static_cast<int>(255.0*v), 0);
}
plut[JKQTPimagePlot_LUTSIZE+1]=plut[JKQTPimagePlot_LUTSIZE];
}
}
lut_used=lutstore[palette];
}
else if (palette == JKQTPMathImageBLUE)
{
//int*& plut=&lut_blue;
if (lutstore[palette]==nullptr) {
lutstore[palette]=(int*)malloc((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=l/(double)(JKQTPimagePlot_LUTSIZE);
plut[l]=qRgb(0, 0, static_cast<int>(255.0*v));
}
plut[JKQTPimagePlot_LUTSIZE+1]=plut[JKQTPimagePlot_LUTSIZE];
}
}
lut_used=lutstore[palette];
}
else if (palette == JKQTPMathImageGRAY)
{
//int*& plut=&lut_gray;
if (lutstore[palette]==nullptr) {
lutstore[palette]=(int*)malloc((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=l/(double)(JKQTPimagePlot_LUTSIZE);
plut[l]=qRgb(static_cast<int>(255.0*v),
static_cast<int>(255.0*v),
static_cast<int>(255.0*v));
}
plut[JKQTPimagePlot_LUTSIZE+1]=plut[JKQTPimagePlot_LUTSIZE];
}
}
lut_used=lutstore[palette];
}
else if (palette == JKQTPMathImageALPHA)
{
//int*& plut=&lut_gray;
if (lutstore[palette]==nullptr) {
lutstore[palette]=(int*)malloc((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=l/(double)(JKQTPimagePlot_LUTSIZE);
plut[l]=qRgba(255,255,255,
static_cast<int>(255.0*v));
}
plut[JKQTPimagePlot_LUTSIZE+1]=plut[JKQTPimagePlot_LUTSIZE];
}
}
lut_used=lutstore[palette];
}
else if (palette == JKQTPMathImageINVERTED_ALPHA)
{
//int*& plut=&lut_gray;
if (lutstore[palette]==nullptr) {
lutstore[palette]=(int*)malloc((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=l/(double)(JKQTPimagePlot_LUTSIZE);
plut[l]=qRgba(255,255,255,
static_cast<int>(255.0*v));
}
plut[JKQTPimagePlot_LUTSIZE+1]=plut[JKQTPimagePlot_LUTSIZE];
}
}
lut_used=lutstore[palette];
}
else if (palette == JKQTPMathImageINVERTEDRED)
{
//int*& plut=&lut_invred;
if (lutstore[palette]==nullptr) {
lutstore[palette]=(int*)malloc((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=l/(double)(JKQTPimagePlot_LUTSIZE);
plut[l]=qRgb(static_cast<int>(255.0*(1.0-v)), 0, 0);
}
plut[JKQTPimagePlot_LUTSIZE+1]=plut[JKQTPimagePlot_LUTSIZE];
}
}
lut_used=lutstore[palette];
}
else if (palette == JKQTPMathImageINVERTEDGREEN)
{
//int*& plut=&lut_invgreen;
if (lutstore[palette]==nullptr) {
lutstore[palette]=(int*)malloc((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=l/(double)(JKQTPimagePlot_LUTSIZE);
plut[l]=qRgb(0, static_cast<int>(255.0*(1.0-v)), 0);
}
plut[JKQTPimagePlot_LUTSIZE+1]=plut[JKQTPimagePlot_LUTSIZE];
}
}
lut_used=lutstore[palette];
}
else if (palette == JKQTPMathImageINVERTEDBLUE)
{
//int*& plut=&lut_invblue;
if (lutstore[palette]==nullptr) {
lutstore[palette]=(int*)malloc((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=l/(double)(JKQTPimagePlot_LUTSIZE);
plut[l]=qRgb(0, 0, static_cast<int>(255.0*(1.0-v)));
}
plut[JKQTPimagePlot_LUTSIZE+1]=plut[JKQTPimagePlot_LUTSIZE];
}
}
lut_used=lutstore[palette];
}
else if (palette == JKQTPMathImageINVERTEDGRAY)
{
//int*& plut=&lut_invgray;
if (lutstore[palette]==nullptr) {
lutstore[palette]=(int*)malloc((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=1.0-(l/(double)(JKQTPimagePlot_LUTSIZE));
plut[l]=qRgb(static_cast<int>(255.0*v),
static_cast<int>(255.0*v),
static_cast<int>(255.0*v));
}
plut[JKQTPimagePlot_LUTSIZE+1]=plut[JKQTPimagePlot_LUTSIZE];
}
}
lut_used=lutstore[palette];
}
else if (palette == JKQTPMathImageMATLAB)
{
//int*& plut=&lut_matlab;
if (lutstore[palette]==nullptr) {
lutstore[palette]=(int*)malloc((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=l/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=JKQTPimagePlot_LUTSIZE; l>=0; l--) {
double v=(JKQTPimagePlot_LUTSIZE-l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=l/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=JKQTPimagePlot_LUTSIZE; l>=0; l--) {
double v=(JKQTPimagePlot_LUTSIZE-l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=l/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=l/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
//std::cout<<"!!! creating rainbow lut\n";
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=l/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
//std::cout<<"!!! creating rainbow lut\n";
if (plut!=nullptr) {
for (int l=JKQTPimagePlot_LUTSIZE; l>=0; l--) {
double v=(JKQTPimagePlot_LUTSIZE-l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=l/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=JKQTPimagePlot_LUTSIZE; l>=0; l--) {
double v=(JKQTPimagePlot_LUTSIZE-l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=l/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=JKQTPimagePlot_LUTSIZE; l>=0; l--) {
double v=(JKQTPimagePlot_LUTSIZE-l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=l/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=(JKQTPimagePlot_LUTSIZE-l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=l/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=(JKQTPimagePlot_LUTSIZE-l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=l/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=(JKQTPimagePlot_LUTSIZE-l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=l/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=(JKQTPimagePlot_LUTSIZE-l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_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((JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
//std::cout<<"!!! creating rainbow lut\n";
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=(l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
//std::cout<<"!!! creating rainbow lut\n";
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=(l)/(double)(JKQTPimagePlot_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 == JKQTPMathImageBLACKBLUEVIOLETYELLOWWHITE) //gnuplot: 30,31,32
{
//int*& plut=&lut_BLACKBLUEVIOLETYELLOWWHITE;
if (lutstore[palette]==nullptr) {
lutstore[palette]=(int*)malloc((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
//std::cout<<"!!! creating rainbow lut\n";
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=(l)/(double)(JKQTPimagePlot_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 == JKQTPMathImageWHITEYELLOWVIOLETBLUEBLACK) //gnuplot: 30,31,32
{
//int*& plut=&lut_WHITEYELLOWVIOLETBLUEBLACK;
if (lutstore[palette]==nullptr) {
lutstore[palette]=(int*)malloc((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
//std::cout<<"!!! creating rainbow lut\n";
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=(JKQTPimagePlot_LUTSIZE-l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_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((JKQTPimagePlot_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((JKQTPimagePlot_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((JKQTPimagePlot_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((JKQTPimagePlot_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((JKQTPimagePlot_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((JKQTPimagePlot_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((JKQTPimagePlot_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((JKQTPimagePlot_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((JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
//std::cout<<"!!! creating rainbow lut\n";
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=(l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
//std::cout<<"!!! creating rainbow lut\n";
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=(JKQTPimagePlot_LUTSIZE-l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
//std::cout<<"!!! creating rainbow lut\n";
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=(l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
//std::cout<<"!!! creating rainbow lut\n";
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=(JKQTPimagePlot_LUTSIZE-l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
//std::cout<<"!!! creating rainbow lut\n";
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=(l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
//std::cout<<"!!! creating rainbow lut\n";
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=(JKQTPimagePlot_LUTSIZE-l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
//std::cout<<"!!! creating rainbow lut\n";
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=(l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
//std::cout<<"!!! creating rainbow lut\n";
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=(JKQTPimagePlot_LUTSIZE-l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
//std::cout<<"!!! creating rainbow lut\n";
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=(l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_LUTSIZE+2)*sizeof(int));
int* plut=lutstore[palette];
//std::cout<<"!!! creating rainbow lut\n";
if (plut!=nullptr) {
for (int l=0; l<=JKQTPimagePlot_LUTSIZE; l++) {
double v=(JKQTPimagePlot_LUTSIZE-l)/(double)(JKQTPimagePlot_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((JKQTPimagePlot_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((JKQTPimagePlot_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((JKQTPimagePlot_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((JKQTPimagePlot_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((JKQTPimagePlot_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((JKQTPimagePlot_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((JKQTPimagePlot_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((JKQTPimagePlot_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=(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((int)palette, width);
}
QImage JKQTPMathImageGetPaletteImage(JKQTPMathImageColorPalette palette, int width, int height)
{
return JKQTPMathImageGetPaletteImage((int)palette, width, height);
}
QImage JKQTPMathImageGetPaletteImage(int *lut, int lut_size, int width)
{
QImage img;
double* pic=(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, JKQTP_PALETTE_ICON_WIDTH);
QPixmap pix(JKQTP_PALETTE_ICON_WIDTH,8);
QRect r(0,0,JKQTP_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((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=(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;
}
#ifndef NO_JKQTPLOTTER
JKQTPColorPaletteTools::JKQTPColorPaletteTools(JKQtBasePlotter *parent)
{
cbParent=parent;
colorBarRightAxis=new JKQTPverticalIndependentAxis(0, 100, 0, 100, parent);
if (parent) colorBarRightAxis->loadSettings(parent->get_yAxis());
colorBarRightAxis->set_drawMode1(JKQTPCADMline);
colorBarRightAxis->set_drawMode2(JKQTPCADMcomplete);
colorBarRightAxis->set_axisLabel("");
colorBarRightAxis->set_minTicks(3);
colorBarRightAxis->set_showZeroAxis(false);
colorBarRightAxis->set_minorTicks(0);
colorBarRightAxis->set_tickOutsideLength(0);
colorBarRightAxis->set_minorTickOutsideLength(0);
colorBarTopAxis=new JKQTPhorizontalIndependentAxis(0, 100, 0, 100, parent);
if (parent) colorBarTopAxis->loadSettings(parent->get_xAxis());
colorBarTopAxis->set_drawMode1(JKQTPCADMline);
colorBarTopAxis->set_drawMode2(JKQTPCADMcomplete);
colorBarTopAxis->set_axisLabel("");
colorBarTopAxis->set_minTicks(3);
colorBarTopAxis->set_showZeroAxis(false);
colorBarTopAxis->set_minorTicks(0);
colorBarTopAxis->set_tickOutsideLength(0);
colorBarTopAxis->set_minorTickOutsideLength(0);
this->palette=JKQTPMathImageMATLAB;
this->imageNameFontName=parent->get_keyFont();
this->imageNameFontSize=parent->get_keyFontSize();
this->imageName="";
this->showColorBar=true;
this->colorBarWidth=14;
this->colorBarRelativeHeight=0.75;
this->autoImageRange=true;
this->imageMin=0;
this->imageMax=1;
this->colorBarOffset=4;
this->rangeMinFailAction=JKQTPMathImageLastPaletteColor;
this->rangeMaxFailAction=JKQTPMathImageLastPaletteColor;
this->rangeMinFailColor=QColor("black");
this->rangeMaxFailColor=QColor("black");
this->nanColor=QColor("black");
this->infColor=QColor("black");
this->colorBarTopVisible=false;
this->colorBarRightVisible=true;
}
JKQTPColorPaletteTools::~JKQTPColorPaletteTools()
{
}
void JKQTPColorPaletteTools::cbGetOutsideSize(JKQTPEnhancedPainter& painter, int& leftSpace, int& rightSpace, int& topSpace, int& bottomSpace) {
Q_UNUSED(leftSpace);
Q_UNUSED(bottomSpace);
if (showColorBar) {
double internalDataMin=getInternalDataMin();
double internalDataMax=getInternalDataMax();
cbGetDataMinMax(internalDataMin, internalDataMax);
if (colorBarRightVisible) {
rightSpace+=cbParent->pt2px(painter, colorBarWidth+colorBarOffset);
colorBarRightAxis->setRange(internalDataMin, internalDataMax);
colorBarRightAxis->set_axisWidth(colorBarRelativeHeight*cbParent->get_plotHeight());
QSize s=colorBarRightAxis->getSize2(painter);
/*cbParent->get_mathText()->set_fontSize(imageNameFontSize*cbParent->get_fontSizeMultiplier());
cbParent->get_mathText()->set_fontRoman(imageNameFontName);
cbParent->get_mathText()->parse(imageName);
QSizeF names=cbParent->get_mathText()->getSize(painter);*/
QSizeF names=cbParent->getTextSizeSize(imageNameFontName, imageNameFontSize*cbParent->get_fontSizeMultiplier(), imageName, painter);
rightSpace+=cbParent->pt2px(painter, colorBarWidth+colorBarOffset)+qMax((double)s.width(), (double)names.width());
}
if (colorBarTopVisible) {
//topSpace+=cbParent->pt2px(painter, colorBarWidth+colorBarOffset);
colorBarTopAxis->setRange(internalDataMin, internalDataMax);
colorBarTopAxis->set_axisWidth(colorBarRelativeHeight*cbParent->get_plotWidth());
QSize s=colorBarTopAxis->getSize2(painter);
/*cbParent->get_mathText()->set_fontSize(imageNameFontSize*cbParent->get_fontSizeMultiplier());
cbParent->get_mathText()->set_fontRoman(imageNameFontName);
cbParent->get_mathText()->parse(imageName);
QSizeF names=cbParent->get_mathText()->getSize(painter);*/
QSizeF names=cbParent->getTextSizeSize(imageNameFontName, imageNameFontSize*cbParent->get_fontSizeMultiplier(), imageName, painter);
topSpace+=cbParent->pt2px(painter, colorBarWidth+colorBarOffset)+qMax((double)s.height(), (double)names.height());
}
}
}
void JKQTPColorPaletteTools::cbDrawOutside(JKQTPEnhancedPainter& painter, QRect leftSpace, QRect rightSpace, QRect topSpace, QRect bottomSpace) {
Q_UNUSED(leftSpace);
Q_UNUSED(bottomSpace);
if (showColorBar) {
double internalDataMin=getInternalDataMin();
double internalDataMax=getInternalDataMax();
cbGetDataMinMax(internalDataMin, internalDataMax);
if (colorBarRightVisible) {
painter.save();
static const uint8_t h=1;
static const uint8_t dSize = 200*h;
uint8_t d[dSize];//, dd[200*h];
for (int i=0; i<200; i++) {
for (int j=0; j<h; j++) {
d[i*h+j]=i;
//dd[i*h+j]=j;
}
}
QImage b(h, 200, QImage::Format_ARGB32);
JKQTPimagePlot_array2image<uint8_t>(d, h, 200, b, palette, 0, 199);
/*cbParent->get_mathText()->set_fontSize(imageNameFontSize*cbParent->get_fontSizeMultiplier());
cbParent->get_mathText()->set_fontRoman(imageNameFontName);
cbParent->get_mathText()->parse(imageName);
QSizeF names=cbParent->get_mathText()->getSize(painter);*/
QSizeF names=cbParent->getTextSizeSize(imageNameFontName, imageNameFontSize*cbParent->get_fontSizeMultiplier(), imageName, painter);
double icolorBarRelativeHeight=colorBarRelativeHeight;
int barHeight=rightSpace.height()*icolorBarRelativeHeight;
if (barHeight<1) barHeight=1;
// find a height for the bar that allows to show the image name
while ((barHeight>1) && ((rightSpace.height()-barHeight)/2<names.height())) {
barHeight--;
}
//qDebug()<<"r: "<<imageName<<rightSpace;
QRect cb(rightSpace.x()+cbParent->pt2px(painter, colorBarOffset), rightSpace.top()+(rightSpace.height()-barHeight)/2, cbParent->pt2px(painter, colorBarWidth), barHeight);
painter.drawImage(cb, b.mirrored(true, false));
QPen p=painter.pen();
p.setColor(colorBarRightAxis->get_axisColor());
p.setWidthF(qMax(JKQTPLOTTER_ABS_MIN_LINEWIDTH, cbParent->pt2px(painter, colorBarRightAxis->get_lineWidth()*cbParent->get_lineWidthMultiplier())));
painter.setPen(p);
painter.drawRect(cb);
colorBarRightAxis->setRange(internalDataMin, internalDataMax);
colorBarRightAxis->set_axisWidth(cb.height());
colorBarRightAxis->set_axisOffset(cb.top());
colorBarRightAxis->set_otherAxisOffset(cb.left());
colorBarRightAxis->set_otherAxisWidth(cb.width());
colorBarRightAxis->drawAxes(painter);
cbParent->get_mathText()->set_fontSize(imageNameFontSize*cbParent->get_fontSizeMultiplier());
cbParent->get_mathText()->set_fontRoman(imageNameFontName);
cbParent->get_mathText()->parse(imageName);
cbParent->get_mathText()->draw(painter, Qt::AlignHCenter|Qt::AlignVCenter, QRect(rightSpace.x(), rightSpace.y(), rightSpace.width(), (rightSpace.height()-barHeight)/2));
painter.restore();
}
if (colorBarTopVisible) {
painter.save();
/*uint8_t d[200];
for (int i=0; i<200; i++) d[i]=i;
QImage b(200,1, QImage::Format_ARGB32);
JKQTPimagePlot_array2image<uint8_t>(d, 200, 1, b, palette, 0, 199);*/
static const uint8_t h=1;
static const uint8_t dSize = 200*h;
uint8_t d[dSize];//, dd[200*h];
for (int i=0; i<200; i++) {
for (int j=0; j<h; j++) {
d[i*h+j]=i;
//dd[i*h+j]=j;
}
}
QImage b(h,200, QImage::Format_ARGB32);
JKQTPimagePlot_array2image<uint8_t>(d,h,200, b, palette, 0, 199);
/*cbParent->get_mathText()->set_fontSize(imageNameFontSize*cbParent->get_fontSizeMultiplier());
cbParent->get_mathText()->set_fontRoman(imageNameFontName);
cbParent->get_mathText()->parse(imageName);
QSizeF names=cbParent->get_mathText()->getSize(painter);*/
QSizeF names=cbParent->getTextSizeSize(imageNameFontName, imageNameFontSize*cbParent->get_fontSizeMultiplier(), imageName, painter);
double icolorBarRelativeHeight=colorBarRelativeHeight;
int barWidth=topSpace.width()*icolorBarRelativeHeight;
if (barWidth<1) barWidth=1;
// find a height for the bar that allows to show the image name
while ((barWidth>1) && ((topSpace.width()-barWidth)/2<names.width())) {
barWidth--;
}
QRect cb(topSpace.x()+(topSpace.width()-barWidth)/2, topSpace.bottom()-cbParent->pt2px(painter, colorBarOffset+(colorBarWidth)), barWidth, cbParent->pt2px(painter, colorBarWidth));
//qDebug()<<"t: "<<imageName<<topSpace<<topSpace.bottom()<<colorBarOffset<<"\n "<<cb;
QMatrix rm;
rm.rotate(90);
painter.drawImage(cb, b.transformed(rm));
QPen p=painter.pen();
p.setColor(colorBarTopAxis->get_axisColor());
p.setWidthF(qMax(JKQTPLOTTER_ABS_MIN_LINEWIDTH, cbParent->pt2px(painter, colorBarTopAxis->get_lineWidth()*cbParent->get_lineWidthMultiplier())));
painter.setPen(p);
painter.drawRect(cb);
colorBarTopAxis->setRange(internalDataMin, internalDataMax);
colorBarTopAxis->set_axisWidth(cb.width());
colorBarTopAxis->set_axisOffset(cb.left());
colorBarTopAxis->set_otherAxisOffset(cb.top());
colorBarTopAxis->set_otherAxisWidth(cb.height());
colorBarTopAxis->drawAxes(painter);
cbParent->get_mathText()->set_fontSize(imageNameFontSize*cbParent->get_fontSizeMultiplier());
cbParent->get_mathText()->set_fontRoman(imageNameFontName);
cbParent->get_mathText()->parse(imageName);
cbParent->get_mathText()->draw(painter, Qt::AlignHCenter|Qt::AlignVCenter, QRect(topSpace.right()-(topSpace.width()-barWidth)/2, topSpace.y(), (topSpace.width()-barWidth)/2, topSpace.height()));
painter.restore();
}
}
}
QStringList JKQTPColorPaletteTools::getPalettes() {
return JKQTPimagePlot_getPredefinedPalettes();
}
int JKQTPColorPaletteTools::getPalettesCount()
{
return getPalettes().size();
}
QIcon JKQTPColorPaletteTools::getPaletteIcon(int i) {
QImage img=getPaletteImage(i, JKQTP_PALETTE_ICON_WIDTH);
QPixmap pix(JKQTP_PALETTE_ICON_WIDTH,8);
QRect r(0,0,JKQTP_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 JKQTPColorPaletteTools::getPaletteIcon(JKQTPMathImageColorPalette palette) {
return getPaletteIcon((int)palette);
}
QImage JKQTPColorPaletteTools::getPaletteImage(int i, int width)
{
QImage img;
double* pic=(double*)malloc(width*sizeof(double));
for (int j=0; j<width; j++) {
pic[j]=j;
}
JKQTPimagePlot_array2image<double>(pic, width, 1, img, (JKQTPMathImageColorPalette)i, 0, width-1);
free(pic);
return img;
}
QImage JKQTPColorPaletteTools::getPaletteImage(JKQTPMathImageColorPalette palette, int width)
{
return getPaletteImage((int)palette, width);
}
QIcon JKQTPColorPaletteTools::getPaletteKeyIcon(int i)
{
QImage img=getPaletteKeyImage(i, JKQTP_PALETTE_ICON_WIDTH, JKQTP_PALETTE_IMAGEICON_HEIGHT);
QPixmap pix(JKQTP_PALETTE_ICON_WIDTH,JKQTP_PALETTE_IMAGEICON_HEIGHT);
QRect r(0,0,JKQTP_PALETTE_ICON_WIDTH-1,JKQTP_PALETTE_IMAGEICON_HEIGHT-1);
QPainter p(&pix);
p.drawImage(r, img);
p.setPen(QPen(QColor("black")));
p.drawRect(r);
p.end();
return QIcon(pix);
}
QIcon JKQTPColorPaletteTools::getPaletteKeyIcon(JKQTPMathImageColorPalette palette)
{
return getPaletteIcon((int)palette);
}
QImage JKQTPColorPaletteTools::getPaletteKeyImage(int i, int width, int height)
{
QImage img;
const double x01=double(width)/3.0;
const double y01=double(height)/3.0*2.0;
const double w1x=double(width*width)/(5.0*5.0);
const double w1y=double(height*height)/(2.0*2.0);
const double x02=double(width)/3.0*2.0;
const double y02=double(height)/4.0;
const double w2x=double(width*width)/(8.0*8.0);
const double w2y=double(height*height)/(8.0*8.0);
double* pic=(double*)malloc(width*height*sizeof(double));
double mmax=0;
for (int j=0; j<width*height; j++) {
const int x=j%width;
const int y=j/width;
pic[j]=exp(-0.5*(double((x-x01)*double(x-x01))/w1x+double((y-y01)*double(y-y01))/w1y))+0.7*exp(-0.5*(double((x-x02)*double(x-x02))/w2x+double((y-y02)*double(y-y02))/w2y));
if (pic[j]>mmax) mmax=pic[j];
}
JKQTPimagePlot_array2image<double>(pic, width, height, img, (JKQTPMathImageColorPalette)i, 0, mmax);
free(pic);
return img;
}
QImage JKQTPColorPaletteTools::getPaletteKeyImage(JKQTPMathImageColorPalette palette, int width, int height)
{
return getPaletteKeyImage((int)palette, width, height);
}
void JKQTPColorPaletteTools::set_palette(int pal) {
palette=(JKQTPMathImageColorPalette)pal;
}
void JKQTPColorPaletteTools::cbSetParent(JKQtBasePlotter* parent) {
cbParent=parent;
colorBarRightAxis->set_parent(parent);
colorBarTopAxis->set_parent(parent);
}
double JKQTPColorPaletteTools::getInternalDataMin() const
{
return 0;
}
double JKQTPColorPaletteTools::getInternalDataMax() const
{
return 0;
}
#endif
Reference in New Issue View Git Blame Copy Permalink