examples/qt/rastervisual/ReadPixelTool.cpp

Go to the documentation of this file.

/*  Copyright (C) 2008 National Institute For Space Research (INPE) - Brazil.

    This file is part of the TerraLib - a Framework for building GIS enabled applications.

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

    TerraLib 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 for more details.

    You should have received a copy of the GNU Lesser General Public License
    along with TerraLib. See COPYING. If not, write to
    TerraLib Team at <terralib-team@terralib.org>.
 */

/*!
  \file ReadPixel.cpp

  \brief An example of MapDisplay Tool. The only purpose of this tool is to show how you can implement a new tool. Do not consider it as a final application.
*/

// Example
#include "ReadPixelTool.h"

// TerraLib
#include <terralib/raster.h>
#include <terralib/raster/RasterSummary.h>
#include <terralib/raster/RasterSummaryManager.h>
#include <terralib/se.h>
#include <terralib/geometry.h>
#include <terralib/common/STLUtils.h>
//#include <terralib/maptools/RasterLayer.h> // * Under revision *
#include <terralib/maptools/RasterTransform.h>
#include <terralib/maptools/RasterTransformConfigurer.h>
#include <terralib/qt/widgets/canvas/Canvas.h>
#include <terralib/qt/widgets/canvas/MapDisplay.h>
#include <terralib/se/Utils.h>


// Qt
#include <QMouseEvent>
#include <QToolTip>

// STL
#include <cassert>

ReadPixelTool::ReadPixelTool(te::qt::widgets::MapDisplay* display, /*te::map::RasterLayer* layer,*/ QObject* parent) 
  : te::qt::widgets::AbstractTool(display, parent)
{
  // * Under revision *
  //assert(layer);
  //m_layer = dynamic_cast<te::map::RasterLayer*>(layer);
  //assert(m_layer);
}

ReadPixelTool::~ReadPixelTool()
{
}

bool ReadPixelTool::mouseReleaseEvent(QMouseEvent* e)
{
  if(e->button() != Qt::LeftButton)
    return false;

  // Converts clicked point to world coordinates
#if QT_VERSION >= 0x050000
  QPointF qpoint = m_display->transform(e->localPos());
#else
  QPointF qpoint = m_display->transform(e->posF());
#endif

  // * Under revision *
  //te::se::RasterSymbolizer* rs = (te::se::RasterSymbolizer*)m_layer->getRasterSymbolizer()->clone();
  te::se::RasterSymbolizer* rs = 0;

  // * Under revision *
  //te::rst::Raster* raster = m_layer->getRaster();
  te::rst::Raster* raster = 0;

  //create the raster transform
  te::map::RasterTransform rt(raster, 0);

  double rMin, rMax;

  const te::rst::RasterSummary* rsMin = te::rst::RasterSummaryManager::getInstance().get(raster, te::rst::SUMMARY_MIN);
  const te::rst::RasterSummary* rsMax = te::rst::RasterSummaryManager::getInstance().get(raster, te::rst::SUMMARY_MAX);
  const std::complex<double>* cmin = rsMin->at(0).m_minVal;
  const std::complex<double>* cmax = rsMax->at(0).m_maxVal;
  rMin = cmin->real();
  rMax = cmax->real();

  rt.setLinearTransfParameters(0, 255, rMin, rMax);

  //configure the raster transformation using the raster symbolizer
  te::map::RasterTransformConfigurer rtc(rs, &rt);

  rtc.configure();

  te::gm::Coord2D coord = raster->getGrid()->geoToGrid(qpoint.x(), qpoint.y());

  int x = te::rst::Round(coord.x);
  int y = te::rst::Round(coord.y);

  if((x >=0 && x < (int)raster->getNumberOfColumns()) && (y >=0 && y < (int)raster->getNumberOfRows()))
  {
    QString information("<h2>Read Pixel</h2><ul>");

    information += "<h3>Location</h3>";

    QString column, line, xCoord, yCoord;
    column.setNum(x);
    line.setNum(y);
    xCoord.setNum(qpoint.x());
    yCoord.setNum(qpoint.y());

    information += "<li><b>" + QString("Line ") + ":</b> " + line + "</li>";
    information += "<li><b>" + QString("Column ") + ":</b> " + column + "</li>";
    information += "<li><b>" + QString("Coord X ") + ":</b> " + xCoord + "</li>";
    information += "<li><b>" + QString("Coord Y ") + ":</b> " + yCoord + "</li>";

    information += "<h3>Original Values</h3>";

    for(size_t t = 0; t < raster->getNumberOfBands(); ++t )
    {
      double val;
      raster->getValue(x, y, val, t);

      QString band, value;
      band.setNum(t);
      value.setNum(val);

      information += "<li><b>" + QString("Value for band ") + band + ":</b> " + value + "</li>";

    }
    
    te::color::RGBAColor color = rt.apply(x, y);

    information += "<h3>Transformed Values</h3>";

    QString red, green, blue, alpha;
    red.setNum(color.getRed());
    green.setNum(color.getGreen());
    blue.setNum(color.getBlue());
    alpha.setNum(color.getAlpha());

    QString rBand, gBand, bBand;
    rBand.setNum(rt.getRGBMap()[te::map::RasterTransform::RED_CHANNEL]);
    gBand.setNum(rt.getRGBMap()[te::map::RasterTransform::GREEN_CHANNEL]);
    bBand.setNum(rt.getRGBMap()[te::map::RasterTransform::BLUE_CHANNEL]);

    information += "<li><b>" + QString("Red Channel: </b>") + rBand + "<b> - Value:</b> " + red + "</li>";
    information += "<li><b>" + QString("Green Channel: </b>") + gBand + "<b> - Value:</b> " + green + "</li>";
    information += "<li><b>" + QString("Blue Channel: </b>") + bBand + "<b> - Value:</b> " + blue + "</li>";
    information += "<li><b>" + QString("Alpha Channel value") + ":</b> " + alpha + "</li>";


    information += "<h3>Symbolizer</h3>";

    QString type;

    if(rs->getChannelSelection()->getColorCompositionType() == te::se::RGB_COMPOSITION)
    {
      type = "RGB Composition";
    }
    else if(rs->getChannelSelection()->getColorCompositionType() == te::se::GRAY_COMPOSITION)
    {
      type = "Gray Composition";
    }
    else
    {
      type = "Unknown Composition";
    }

    QString opacity = te::se::GetString(rs->getOpacity()).c_str();
    QString gain = te::se::GetString(rs->getGain()).c_str();
    QString offset = te::se::GetString(rs->getOffset()).c_str();

    information += "<li><b>" + QString("Color Composition ") + ":</b> " + type + "</li>";
    information += "<li><b>" + QString("Opacity Value ") + ":</b> " + opacity + "</li>";
    information += "<li><b>" + QString("Gain Value ") + ":</b> " + gain + "</li>";
    information += "<li><b>" + QString("Offset Value ") + ":</b> " + offset + "</li>";

    if(rs->getChannelSelection()->getRedChannel() && 
       rs->getChannelSelection()->getRedChannel()->getContrastEnhancement())
    {
      QString c;
      c.setNum(rs->getChannelSelection()->getRedChannel()->getContrastEnhancement()->getGammaValue());

      information += "<li><b>" + QString("Red Contrast") + ":</b> " + c + "</li>";
    }

    if(rs->getChannelSelection()->getGreenChannel() && 
       rs->getChannelSelection()->getGreenChannel()->getContrastEnhancement())
    {
      QString c;
      c.setNum(rs->getChannelSelection()->getGreenChannel()->getContrastEnhancement()->getGammaValue());

      information += "<li><b>" + QString("Green Contraste") + ":</b> " + c + "</li>";
    }

    if(rs->getChannelSelection()->getBlueChannel() && 
       rs->getChannelSelection()->getBlueChannel()->getContrastEnhancement())
    {
      QString c;
      c.setNum(rs->getChannelSelection()->getBlueChannel()->getContrastEnhancement()->getGammaValue());

      information += "<li><b>" + QString("Blue Contraste") + ":</b> " + c + "</li>";
    }

    if(rs->getChannelSelection()->getGrayChannel() && 
       rs->getChannelSelection()->getGrayChannel()->getContrastEnhancement())
    {
      QString c;
      c.setNum(rs->getChannelSelection()->getGrayChannel()->getContrastEnhancement()->getGammaValue());

      information += "<li><b>" + QString("Gray Contraste") + ":</b> " + c + "</li>";
    }

    information += "</ul>";

    // Show attributes
    QToolTip::showText(QCursor::pos(), information, m_display, m_display->rect());
  }

  return true;
}