us_resids_bitmap.cpp

Go to the documentation of this file.

#include "us_resids_bitmap.h"
#include "us_settings.h"
#include "us_gui_settings.h"

#include <qwt_legend.h>

// constructor:  residuals bitmap widget
US_ResidsBitmap::US_ResidsBitmap( QVector< QVector< double > >& resids,
   QWidget* wparent )
   : US_WidgetsDialog( wparent, 0 )
{

   setWindowTitle( tr( "Residuals Pixel Map" ) );
   setPalette( US_GuiSettings::frameColor() );

   QGridLayout* main = new QGridLayout( this );
   main->setSpacing        ( 2 );
   main->setContentsMargins( 2, 2, 2, 2 );

   lb_bitmap         = new QLabel();      // will draw to a label

   main->addWidget( lb_bitmap );

   replot( resids );

   setAttribute( Qt::WA_DeleteOnClose, true );
}

void US_ResidsBitmap::replot( QVector< QVector< double > >& resids )
{
   // get standard deviation of residuals
   int    scanCount  = resids.size();
   int    valCount   = resids[ 0 ].size();
   double sigma      = 0.0;
   double kfactor    = 2.0;

   for ( int ii = 0; ii < scanCount; ii++ )
   {  // get sum of squares for sigma (std.dev.)
      for ( int jj = 0; jj < valCount; jj++ )
      {
         sigma    += pow( resids[ ii ][ jj ], 2.0 );
      }
   }

   sigma          = sqrt( sigma / (double)( scanCount * valCount ) );
   sigma         *= kfactor;    // standard deviation times k-factor

   int    red     = 0;
   int    green   = 0;
   int    blue    = 0;
   int    kres;
   double vres;

   QImage res_image( valCount, scanCount, QImage::Format_ARGB32 );

   // build image with pixel values colored by residual relation to stdev

   for ( int ii = 0; ii < scanCount; ii++ )
   {
      for ( int jj = 0; jj < valCount; jj++ )
      {
         vres      = resids[ ii ][ jj ];
         kres      = qRound( 255.0 * vres / sigma );

         if ( kres >= 0  )
         {  // positive: scaled redish or red
            red       = kres < 255 ? kres : 255;
            green     = 0;
         }

         else
         {  // negative: scaled greenish or green
            kres      = -kres;
            red       = 0;
            green     = kres < 255 ? kres : 255;
         }

         // set pixel
         res_image.setPixel( jj, ii, QColor( red, green, blue ).rgb() );
      }
   }

   int    mwid    = valCount  < 600 ? valCount  : 600;
   int    mhgt    = scanCount < 100 ? scanCount : 100;

   // create the label pixmap from the image
   resize( mwid + 4, mhgt + 4 );
   lb_bitmap->resize( mwid, mhgt );
   QPixmap pmap( mwid, mhgt );
   pmap = pmap.fromImage( res_image.scaled( mwid, mhgt ) );

   lb_bitmap->setPixmap( pmap );
}