imagesc.cpp

Go to the documentation of this file.

//=============================================================================================================
 
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
#include "imagesc.h"
#include "helpers/colormap.h"
 
//=============================================================================================================
// QT INCLUDES
//=============================================================================================================
 
#include <QPainter>
 
//=============================================================================================================
// USED NAMESPACES
//=============================================================================================================
 
using namespace DISPLIB;
using namespace Eigen;
 
//=============================================================================================================
// DEFINE MEMBER METHODS
//=============================================================================================================
 
ImageSc::ImageSc(QWidget* parent)
: Graph(parent)
{
    init();
}
 
//=============================================================================================================
 
ImageSc::ImageSc(MatrixXd& p_dMat, QWidget* parent)
: Graph(parent)
{
    init();
    updateData(p_dMat);
}
 
//=============================================================================================================
 
ImageSc::ImageSc(MatrixXf& p_fMat, QWidget* parent)
: Graph(parent)
{
    init();
    updateData(p_fMat);
}
 
//=============================================================================================================
 
ImageSc::ImageSc(MatrixXi& p_iMat, QWidget* parent)
: Graph(parent)
{
    init();
    updateData(p_iMat);
}
 
//=============================================================================================================
 
ImageSc::~ImageSc() = default;
 
//=============================================================================================================
 
void ImageSc::init()
{
    //Colormap
    pColorMapper = ColorMap::valueToColor;
    m_sColorMap = "Hot";
 
    //Parent init
    Graph::init();
 
    m_iBorderTopBottom = 20;
    m_iBorderLeftRight = 60;
 
    //Colorbar
    m_bColorbar = true;
    m_qFontColorbar.setPixelSize(10);
    m_qPenColorbar = QPen(Qt::black);
    m_iColorbarWidth = 12;
    m_iColorbarSteps = 7; //>= 2!!
    m_iColorbarGradSteps = 200;
}
 
//=============================================================================================================
 
void ImageSc::updateData(MatrixXd& p_dMat)
{
    if (p_dMat.rows() > 0 && p_dMat.cols() > 0) {
        m_dMinValue = p_dMat.minCoeff();
        m_dMaxValue = p_dMat.maxCoeff();
 
        // -- data --
        m_matCentNormData = p_dMat;
        double minValue = m_matCentNormData.minCoeff();
        m_matCentNormData.array() -= minValue;
        double maxValue = m_matCentNormData.maxCoeff();
        // A constant matrix has no range: keep it at 0 instead of dividing 0 by 0
        if (maxValue > 0.0) {
            m_matCentNormData.array() /= maxValue;
        }
 
        updateMaps();
    }
}
 
//=============================================================================================================
 
void ImageSc::updateData(MatrixXf& p_fMat)
{
    MatrixXd t_dMat = p_fMat.cast<double>();
    updateData(t_dMat);
}
 
//=============================================================================================================
 
void ImageSc::updateData(MatrixXi& p_iMat)
{
    MatrixXd t_dMat = p_iMat.cast<double>();
    updateData(t_dMat);
}
 
//=============================================================================================================
 
void ImageSc::updateMaps()
{
    m_pixmapData = QPixmap();
    m_pixmapColorbar = QPixmap();
 
    if (m_matCentNormData.rows() > 0 && m_matCentNormData.cols() > 0) {
        // --Data--
        qint32 x = m_matCentNormData.cols();
        qint32 y = m_matCentNormData.rows();
        QImage t_qImageData(x, y, QImage::Format_RGB32);
        for (qint32 i = 0; i < x; ++i)
            for (qint32 j = 0; j < y; ++j)
                t_qImageData.setPixel(i, j, pColorMapper(m_matCentNormData(j, i), m_sColorMap));
 
        m_pixmapData = QPixmap::fromImage(t_qImageData);
 
        // --Colorbar--
        QImage t_qImageColorbar(1, m_iColorbarGradSteps, QImage::Format_RGB32);
 
        double t_dQuantile = 1.0 / (static_cast<double>(m_iColorbarGradSteps) - 1);
        for (qint32 j = 0; j < m_iColorbarGradSteps; ++j) {
            QRgb t_qRgb = pColorMapper(t_dQuantile * (static_cast<double>(m_iColorbarGradSteps - 1 - j)) * 1.0, m_sColorMap);
            t_qImageColorbar.setPixel(0, j, t_qRgb);
        }
        m_pixmapColorbar = QPixmap::fromImage(t_qImageColorbar);
 
        // --Scale Values--
        m_qVecScaleValues.clear();
        if (m_dMaxValue <= m_dMinValue) {
            update();
            return;
        }
 
        double scale = pow(10, floor(log(m_dMaxValue - m_dMinValue) / log(10.0)));
 
        //Zero Based Scale?
        if (m_dMaxValue > 0 && m_dMinValue < 0) {
            double quantum = floor((((m_dMaxValue - m_dMinValue) / scale) / (m_iColorbarSteps - 1)) * 10.0) * (scale / 10.0);
            double start = 0;
            while (m_dMinValue < (start - quantum))
                start -= quantum;
            //Create Steps
            m_qVecScaleValues.push_back(start);
            for (qint32 i = 1; i < m_iColorbarSteps - 1; ++i)
                m_qVecScaleValues.push_back(m_qVecScaleValues[i - 1] + quantum);
        } else {
            double quantum = floor((((m_dMaxValue - m_dMinValue) / scale) / (m_iColorbarSteps - 1)) * 10.0) * (scale / 10.0);
            double start = floor(((m_dMaxValue - m_dMinValue) / 2.0 + m_dMinValue) / scale) * scale;
            while (m_dMinValue < (start - quantum))
                start -= quantum;
            //Create Steps
            m_qVecScaleValues.push_back(start);
            for (qint32 i = 1; i < m_iColorbarSteps - 1; ++i)
                m_qVecScaleValues.push_back(m_qVecScaleValues[i - 1] + quantum);
        }
        update();
    }
}
 
//=============================================================================================================
 
void ImageSc::setColorMap(const QString& p_sColorMap)
{
    m_sColorMap = p_sColorMap;
 
    updateMaps();
}
 
//=============================================================================================================
 
void ImageSc::paintEvent(QPaintEvent* event)
{
    Q_UNUSED(event);
 
    QPainter painter(this);
    if (!m_pixmapData.isNull()) {
        QPoint t_qPointTopLeft(0, 0);
 
        // -- Data --
        QSize t_qSizePixmapData = m_qSizeWidget;
 
        t_qSizePixmapData.setHeight(t_qSizePixmapData.height() - m_iBorderTopBottom * 2);
        t_qSizePixmapData.setWidth(t_qSizePixmapData.width() - m_iBorderLeftRight * 2);
        // Scale data
        QPixmap t_qPixmapScaledData = m_pixmapData.scaled(t_qSizePixmapData, Qt::IgnoreAspectRatio); //Qt::KeepAspectRatio);
        // Calculate data position
        t_qPointTopLeft.setX((m_qSizeWidget.width() - t_qPixmapScaledData.width()) / 2);
        t_qPointTopLeft.setY((m_qSizeWidget.height() - t_qPixmapScaledData.height()) / 2);
        //Draw data
        painter.drawPixmap(t_qPointTopLeft, t_qPixmapScaledData);
        //Draw border
        painter.drawRect(t_qPointTopLeft.x() - 1, t_qPointTopLeft.y() - 1, t_qPixmapScaledData.width() + 1, t_qPixmapScaledData.height() + 1);
 
        // -- Colorbar --
        if (m_bColorbar && !m_pixmapColorbar.isNull() && m_qVecScaleValues.size() >= 2) {
            QSize t_qSizePixmapColorbar = m_qSizeWidget;
 
            t_qSizePixmapColorbar.setWidth(m_iColorbarWidth);
            t_qSizePixmapColorbar.setHeight(t_qPixmapScaledData.height());
 
            // Scale colorbar
            QPixmap t_qPixmapScaledColorbar = m_pixmapColorbar.scaled(t_qSizePixmapColorbar, Qt::IgnoreAspectRatio);
            // Calculate colorbar position
            t_qPointTopLeft.setY(t_qPointTopLeft.y());
            t_qPointTopLeft.setX(t_qPointTopLeft.x() + t_qPixmapScaledData.width() + m_iBorderLeftRight / 3);
            //Draw colorbar
            painter.drawPixmap(t_qPointTopLeft, t_qPixmapScaledColorbar);
            //Draw border
            painter.drawRect(t_qPointTopLeft.x() - 1, t_qPointTopLeft.y() - 1, m_iColorbarWidth + 1, t_qPixmapScaledData.height() + 1);
 
            // -- Scale --
            painter.setPen(m_qPenColorbar);
            painter.setFont(m_qFontColorbar);
 
            qint32 x = t_qPointTopLeft.x() + m_iColorbarWidth + m_qFontColorbar.pixelSize() / 2;
            qint32 x_markLeft = x - m_iColorbarWidth - m_qFontColorbar.pixelSize() / 2;
            qint32 x_markRight = x_markLeft + m_iColorbarWidth;
 
            // max
            painter.save();
            qint32 y_max = t_qPointTopLeft.y() - m_qFontColorbar.pixelSize() / 2;
            painter.translate(x, y_max - 1);
            painter.drawText(QRect(0, 0, 100, 12), Qt::AlignLeft, QString::number(m_dMaxValue));
            painter.restore();
            //draw max marks
            qint32 y_max_mark = y_max + m_qFontColorbar.pixelSize() / 2;
            painter.drawLine(x_markLeft, y_max_mark, x_markLeft + 2, y_max_mark);
            painter.drawLine(x_markRight - 3, y_max_mark, x_markRight - 1, y_max_mark);
 
            // min
            painter.save();
            qint32 y_min = t_qPointTopLeft.y() + t_qSizePixmapColorbar.height() - 1 - m_qFontColorbar.pixelSize() / 2;
            painter.translate(x, y_min - 1);
            painter.drawText(QRect(0, 0, 100, 12), Qt::AlignLeft, QString::number(m_dMinValue));
            painter.restore();
            //draw min marks
            qint32 y_min_mark = y_min + m_qFontColorbar.pixelSize() / 2;
            painter.drawLine(x_markLeft, y_min_mark, x_markLeft + 2, y_min_mark);
            painter.drawLine(x_markRight - 3, y_min_mark, x_markRight - 1, y_min_mark);
 
            //Scale values
            qint32 y_dist = y_min - y_max;
            double minPercent = (m_qVecScaleValues[0] - m_dMinValue) / (m_dMaxValue - m_dMinValue);
            double distPercent = (m_qVecScaleValues[1] - m_qVecScaleValues[0]) / (m_dMaxValue - m_dMinValue);
            qint32 y_current = y_min - (minPercent * y_dist);
            qint32 y_current_mark;
            //draw scale
            for (qint32 i = 0; i < m_qVecScaleValues.size(); ++i) {
                painter.save();
                painter.translate(x, y_current - 1);
                painter.drawText(QRect(0, 0, 100, 12), Qt::AlignLeft, QString::number(m_qVecScaleValues[i]));
                painter.restore();
                //draw marks
                y_current_mark = y_current + m_qFontColorbar.pixelSize() / 2;
                painter.drawLine(x_markLeft, y_current_mark, x_markLeft + 2, y_current_mark);
                painter.drawLine(x_markRight - 3, y_current_mark, x_markRight - 1, y_current_mark);
                //update y_current
                y_current -= distPercent * y_dist;
            }
        }
 
        //Draw title & axes
        Graph::drawLabels(t_qPixmapScaledData.width(), t_qPixmapScaledData.height());
    }
}