hugin1/hugin/GLRenderer.cpp

// -*- c-basic-offset: 4 -*-
/** @file GLRenderer.cpp
 *
 *  @author James Legg
 *  @author Darko Makreshanski
 *
 *  This is free software; you can redistribute it and/or
 *  modify it under the terms of the GNU General Public
 *  License as published by the Free Software Foundation; either
 *  version 2 of the License, or (at your option) any later version.
 *
 *  This software 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 General Public
 *  License along with this software. If not, see
 *  <http://www.gnu.org/licenses/>.
 *
 */

#include <wx/wx.h>
#include <wx/platform.h>

#ifdef __WXMAC__
#include <OpenGL/gl.h>
#include <OpenGL/glu.h>
#else
#ifdef __WXMSW__
#include <vigra/windows.h>
#endif
#include <GL/gl.h>
#include <GL/glu.h>

#endif

#include "hugin_config.h"

#include "panoinc.h"

#include "TextureManager.h"
#include "MeshManager.h"
#include "ViewState.h"
#include "GLRenderer.h"
#include "GLViewer.h"
#include "ToolHelper.h"
#include <panodata/PanoramaOptions.h>

GLRenderer::GLRenderer(const wxColour backgroundColour)
{
    m_background_color = backgroundColour;
};

void GLRenderer::SetPreviewBackgroundColor(const wxColour c)
{
    m_background_color = c;
}

void GLRenderer::SetBackground(unsigned char red, unsigned char green, unsigned char blue)
{
    glClearColor((float) red / 255.0, (float) green / 255.0, (float) blue / 255.0, 1.0);
}

GLRenderer::~GLRenderer()
{
}

GLPreviewRenderer::GLPreviewRenderer(HuginBase::Panorama *pano, TextureManager *tex_man,
                       MeshManager *mesh_man, VisualizationState *visualization_state,
                       PreviewToolHelper *tool_helper,const wxColour backgroundColour) : GLRenderer(backgroundColour)
{
    m_pano = pano;
    m_tex_man = tex_man;
    m_mesh_man = mesh_man;
    m_visualization_state = visualization_state;
    m_tool_helper = tool_helper;
}

GLPanosphereOverviewRenderer::GLPanosphereOverviewRenderer(HuginBase::Panorama *pano, TextureManager *tex_man,
                       MeshManager *mesh_man, PanosphereOverviewVisualizationState *visualization_state,
                       PanosphereOverviewToolHelper *tool_helper, const wxColour backgroundColour) : GLRenderer(backgroundColour)
{
    m_pano = pano;
    m_tex_man = tex_man;
    m_mesh_man = mesh_man;
    m_visualization_state = visualization_state;
    m_tool_helper = tool_helper;
}

GLPlaneOverviewRenderer::GLPlaneOverviewRenderer(HuginBase::Panorama *pano, TextureManager *tex_man,
                       MeshManager *mesh_man, PlaneOverviewVisualizationState *visualization_state,
                       PlaneOverviewToolHelper *tool_helper, const wxColour backgroundColour) : GLRenderer(backgroundColour)
{

    m_pano = pano;
    m_tex_man = tex_man;
    m_mesh_man = mesh_man;
    m_visualization_state = visualization_state;
    m_tool_helper = tool_helper;
}

vigra::Diff2D GLPreviewRenderer::Resize(int in_width, int in_height)
{
  width = in_width;
  height = in_height;
  glViewport(0, 0, width, height);
  // notify viewstate about current canvas size,
  // maybe needed for a changed image center
  m_visualization_state->SetCanvasSize(vigra::Size2D(width, height));
  // we use the view_state rather than the panorama to allow interactivity.
  HuginBase::PanoramaOptions *options = m_visualization_state->getViewState()->GetOptions();
  width_o = options->getWidth();
  height_o = options->getHeight();
  double aspect_screen = double(width) / double (height),
        aspect_pano = width_o / height_o;
  glMatrixMode(GL_PROJECTION);
  glLoadIdentity();
  double scale;
  if (aspect_screen < aspect_pano)
  {
      // the panorama is wider than the screen
      scale = width_o / width;
  } else {
      // the screen is wider than the panorama
      scale = height_o / height;
  }
  hugin_utils::FDiff2D center(0.5, 0.5);
  if (!m_visualization_state->GetMeshManager()->GetLayoutMode())
  {
      // apply scale and translate only when not in layout mode
      scale /= m_visualization_state->GetZoomLevel();
      center = m_visualization_state->GetViewingCenter();
  };
  const hugin_utils::FDiff2D lookAt((center.x - 0.5) * width_o, (center.y - 0.5) * height_o);

  double x_offs = (scale * double(width) - width_o) / 2.0,
         y_offs = (scale * double(height) - height_o) / 2.0;
  // set up the projection, so we can use panorama coordinates.
  glOrtho(-x_offs, width * scale - x_offs,
          height * scale - y_offs, -y_offs,
          -1.0, 1.0);
  glTranslatef(-lookAt.x, -lookAt.y, 0);
  // scissor to the panorama.
  glScissor((x_offs-lookAt.x) / scale, (y_offs+lookAt.y) / scale,
            width_o / scale, height_o / scale);
  glMatrixMode(GL_MODELVIEW);
  // tell the view state the region we are displaying.
  // TODO add support for zooming and panning.
  m_visualization_state->SetVisibleArea(vigra::Rect2D((center.x-0.5) * width_o, (center.y-0.5) * height_o,
      (0.5 + center.x) * width_o, (0.5 + center.y) * height_o));
  m_visualization_state->SetScale(1.0 / scale);

  // return the offset from the top left corner of the viewpoer to the top left
  // corner of the panorama.
  return vigra::Diff2D(int (x_offs / scale), int (y_offs / scale));
}

void GLPreviewRenderer::Redraw()
{
    glClear(GL_COLOR_BUFFER_BIT);
    glEnable(GL_SCISSOR_TEST);
    m_tex_man->DisableTexture();
    // background color of flat pano preview
    glColor3f((float)m_background_color.Red()/255, (float)m_background_color.Green()/255, (float)m_background_color.Blue()/255);
    glBegin(GL_QUADS);
        glVertex2f(0.0, 0.0);
        glVertex2f(width_o, 0.0);
        glVertex2f(width_o, height_o);
        glVertex2f(0.0, height_o);
    glEnd();
    glColor3f(1.0, 1.0, 1.0);
    // draw things under the preview images
    // draw each active image.
    int imgs = m_pano->getNrOfImages();
    // offset by a half a pixel
    glPushMatrix();
    glTranslatef(0.5, 0.5, 0.0);
    glEnable(GL_TEXTURE_2D);
    m_tex_man->Begin();
    m_tool_helper->BeforeDrawImages();
    // The old preview shows the lowest numbered image on top, so do the same:
    for (int img = imgs - 1; img != -1; img--)
    {
        // only draw active images
        if (m_pano->getImage(img).getActive())
        {
            // the tools can cancel drawing of images.
            if (m_tool_helper->BeforeDrawImageNumber(img))
            {
                // the texture manager may need to call the display list
                // multiple times with blending, so we pass it the display list
                // rather than switching to the texture and then calling the
                // list ourselves.
                m_tex_man->DrawImage(img, m_mesh_man->GetDisplayList(img));
                m_tool_helper->AfterDrawImageNumber(img);
            }
        }
    }
    m_tex_man->End();
    // drawn things after the active image.
    m_tool_helper->AfterDrawImages();
    m_tex_man->DisableTexture();
    glPopMatrix();
    // darken the cropped out range, except in layout mode
    if (!m_visualization_state->GetMeshManager()->GetLayoutMode())
    {
        glBlendFunc(GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
        glEnable(GL_BLEND);
        glColor4f(0.0, 0.0, 0.0, 0.5);
        // construct a strip of quads, with each pair being one of the corners.
        const vigra::Rect2D roi = m_visualization_state->getViewState()->GetOptions()->getROI();
        glBegin(GL_QUAD_STRIP);
            glVertex2f(0.0, 0.0);          glVertex2i(roi.left(), roi.top());
            glVertex2f(width_o, 0.0);      glVertex2i(roi.right(), roi.top());
            glVertex2f(width_o, height_o); glVertex2i(roi.right(), roi.bottom());
            glVertex2f(0.0, height_o);     glVertex2i(roi.left(), roi.bottom());
            glVertex2f(0.0, 0.0);          glVertex2i(roi.left(), roi.top());
        glEnd();
        // draw lines around cropped area.
        // we want to invert the color to make it stand out.
        glBlendFunc(GL_ONE_MINUS_DST_COLOR, GL_ZERO);
        glColor3f(1.0, 1.0, 1.0);
        glBegin(GL_LINE_LOOP);
            glVertex2i(roi.left(), roi.top());
            glVertex2i(roi.right(), roi.top());
            glVertex2i(roi.right(), roi.bottom());
            glVertex2i(roi.left(), roi.bottom());
        glEnd();
        glDisable(GL_BLEND);
    };
    glEnable(GL_TEXTURE_2D);

    glDisable(GL_SCISSOR_TEST);
}


void GLPanosphereOverviewRenderer::Redraw()
{
    glClearColor(0,0,0,1);
	glClear(GL_COLOR_BUFFER_BIT);

	glMatrixMode(GL_MODELVIEW);

	glLoadIdentity();

    double R = m_visualization_state->getR();
    double angx = m_visualization_state->getAngX();
    double angy = m_visualization_state->getAngY();

	gluLookAt(R * cos(angy) * cos(angx), R * sin(angy), R * cos(angy) * sin(angx), 0, 0, 0, 0, 1, 0);
    //for look from inside
//	gluLookAt(0,0,0,R * cos(angy) * cos(angx), R * sin(angy), R * cos(angy) * sin(angx), 0, 1, 0);

    // draw things under the preview images
    // draw each active image.
    int imgs = m_pano->getNrOfImages();
    // offset by a half a pixel
    glPushMatrix();

    //draw the rectangle around the sphere
    glColor3f(0.5, 0.5, 0.5);

    double side = 150;
    glBegin(GL_LINE_LOOP);

        glVertex3f(-side,side,0);
        glVertex3f(side,side,0);
        glVertex3f(side,-side,0);
        glVertex3f(-side,-side,0);

    glEnd();

    //draw the axes, to give a sense of orientation
    double axis = 200;
    glBegin(GL_LINES);

        glColor3f(1,0,0);
        glVertex3f(-axis,0,0);
        glVertex3f(axis,0,0);

        glColor3f(0,1,0);
        glVertex3f(0,0,0);
        glVertex3f(0,axis,0);

        glColor3f(0,0,1);
        glVertex3f(0,0,0);
        glVertex3f(0,0,axis);

    glEnd();


    glEnable(GL_TEXTURE_2D);

    //To avoid z-order fight of the images if depth buffer is used, depth buffer is disabled and meshes are drawn twice,
    //first with back faces culled so that the inner face of the sphere is visible and below the outter face,
    //and afterwards the meshes are drawn again with the front faces culled

    glEnable(GL_CULL_FACE);
    glCullFace(GL_BACK);

    //event called only before drawing of the images with front faces culled (the inner face of the panosphere)
    static_cast<PanosphereOverviewToolHelper*>(m_tool_helper)->BeforeDrawImagesBack();
    //generic draw before images are drawn (called twice with front and back faces culled)
    m_tex_man->Begin();
    m_tool_helper->BeforeDrawImages();

    // The old preview shows the lowest numbered image on top, so do the same:
    for (int img = imgs - 1; img != -1; img--)
    {
        // only draw active images
        if (m_pano->getImage(img).getActive())
        {
            // the tools can cancel drawing of images.
            if (m_tool_helper->BeforeDrawImageNumber(img))
            {
                // the texture manager may need to call the display list
                // multiple times with blending, so we pass it the display list
                // rather than switching to the texture and then calling the
                // list ourselves.
                m_tex_man->DrawImage(img, m_mesh_man->GetDisplayList(img));
                m_tool_helper->AfterDrawImageNumber(img);
            }
        }
    }

    m_tex_man->End();
    m_tool_helper->AfterDrawImages();
    m_tex_man->DisableTexture();
    static_cast<PanosphereOverviewToolHelper*>(m_tool_helper)->AfterDrawImagesBack();

//    #ifdef __WXGTK__
////    glCullFace(GL_BACK);
////    glPushMatrix();
////    glRotated(90,1,0,0);
//////    if (imgs > 0) {
//////        glEnable( GL_TEXTURE_2D );
//////        glEnable(GL_BLEND);
//////        glBlendFunc (GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
//////        glColor4f(0.5,0.5,0.5,0.5);
//////        GLUquadric* grid = gluNewQuadric();
//////        gluQuadricTexture(grid, GL_TRUE);
//////        m_tex_man->BindTexture(0);
//////        gluSphere(grid, 101,40,20);
//////        glDisable(GL_BLEND);
//////    } else {
////        glColor4f(0.5,0.5,0.5,0.5);
////        glutWireSphere(101,40,20);
//////    }
////    glPopMatrix();
//    #endif

    glMatrixMode(GL_MODELVIEW);
    glCullFace(GL_FRONT);

    static_cast<PanosphereOverviewToolHelper*>(m_tool_helper)->BeforeDrawImagesFront();

    glEnable(GL_TEXTURE_2D);
    m_tex_man->Begin();
    m_tool_helper->BeforeDrawImages();

    // The old preview shows the lowest numbered image on top, so do the same:
    for (int img = imgs - 1; img != -1; img--)
    {
        // only draw active images
        if (m_pano->getImage(img).getActive())
        {
            // the tools can cancel drawing of images.
            if (m_tool_helper->BeforeDrawImageNumber(img))
            {
                // the texture manager may need to call the display list
                // multiple times with blending, so we pass it the display list
                // rather than switching to the texture and then calling the
                // list ourselves.
                m_tex_man->DrawImage(img, m_mesh_man->GetDisplayList(img));
                m_tool_helper->AfterDrawImageNumber(img);
            }
        }
    }

    m_tex_man->End();
    // drawn things after the active image.
    m_tool_helper->AfterDrawImages();
    m_tex_man->DisableTexture();
    static_cast<PanosphereOverviewToolHelper*>(m_tool_helper)->AfterDrawImagesFront();

    m_tex_man->DisableTexture();

    glDisable(GL_CULL_FACE);

    glPopMatrix();
}

vigra::Diff2D GLPanosphereOverviewRenderer::Resize(int w, int h)
{

    width = w;
    height = h;
    glViewport(0, 0, width, height);
    // we use the view_state rather than the panorama to allow interactivity.
    HuginBase::PanoramaOptions *options = m_visualization_state->GetOptions();
    width_o = options->getWidth();
    height_o = options->getHeight();

    //since gluPerspective needs vertical field of view, depending on the aspect ratio we convert from vertical to horizontal FOV if needed
    double fov = m_visualization_state->getFOV();
    double fovy;
    if (h > w) {
        fovy = 2.0 * atan( tan(fov * M_PI / 360.0) * (float) h / (float) w) / M_PI * 180.0;
    } else {
        fovy = fov;
    }

	float ratio = 1.0* w / h;
//	aspect = ratio;
	glMatrixMode(GL_PROJECTION);
	glLoadIdentity();
	glViewport(0, 0, w, h);
	gluPerspective(fovy,ratio,1,1000000);

    m_visualization_state->SetVisibleArea(vigra::Rect2D(0, 0, options->getWidth(),
                                             options->getHeight()));

    //calculate the scale depending on the section of the panosphere in the center of the screen
    double R = m_visualization_state->getR();
    double radius = m_visualization_state->getSphereRadius();
    //height of the screen in screen pixels over the length of the panosphere in panorama pixels when spread out
    double scrscale = (float) h  / (2 * tan(fovy / 360.0 * M_PI) * (R - radius) / (2 * radius * M_PI) * (options->getWidth()));
    m_visualization_state->SetScale(scrscale);
    m_visualization_state->GetViewer()->MarkToolsDirty();
//    DEBUG_DEBUG("renderer " << scrscale << " " << h << " " << R << " " << fovy);
//    DEBUG_DEBUG("renderer scale " << scrscale);

//    return vigra::Diff2D(w / 2, h / 2);
    return vigra::Diff2D(0,0);

}

void GLPlaneOverviewRenderer::Redraw()
{
    // background color of mosaic plane
    glClearColor((float)m_background_color.Red()/255, (float)m_background_color.Green()/255, (float)m_background_color.Blue()/255,1.0);

	glClear(GL_COLOR_BUFFER_BIT);

	glMatrixMode(GL_MODELVIEW);

	glLoadIdentity();

    double R = m_visualization_state->getR();

	double X = m_visualization_state->getX();
	double Y = m_visualization_state->getY();

	gluLookAt(X,Y,R, X, Y, 0, 0, 1, 0);

    // draw things under the preview images
    m_tool_helper->BeforeDrawImages();
    int imgs = m_pano->getNrOfImages();
    glPushMatrix();

    glColor3f(0.5,0.5,0.5);
    double side = 150;
    glBegin(GL_LINE_LOOP);

        glVertex3f(-side,side,0);
        glVertex3f(side,side,0);
        glVertex3f(side,-side,0);
        glVertex3f(-side,-side,0);

    glEnd();

    double axis = 200;
    glBegin(GL_LINES);

        glColor3f(1,0,0);
        glVertex3f(-axis,0,0);
        glVertex3f(axis,0,0);

        glColor3f(0,1,0);
        glVertex3f(0,0,0);
        glVertex3f(0,axis,0);

        glColor3f(0,0,1);
        glVertex3f(0,0,0);
        glVertex3f(0,0,axis);

    glEnd();


    glEnable(GL_TEXTURE_2D);

    m_tex_man->Begin();
    // The old preview shows the lowest numbered image on top, so do the same:
    for (int img = imgs - 1; img != -1; img--)
    {
        // only draw active images
        if (m_pano->getImage(img).getActive())
        {
            // the tools can cancel drawing of images.
            if (m_tool_helper->BeforeDrawImageNumber(img))
            {
                // the texture manager may need to call the display list
                // multiple times with blending, so we pass it the display list
                // rather than switching to the texture and then calling the
                // list ourselves.
                m_tex_man->DrawImage(img, m_mesh_man->GetDisplayList(img));
                m_tool_helper->AfterDrawImageNumber(img);
            }
        }
    }

    m_tex_man->End();
    m_tool_helper->AfterDrawImages();

    m_tex_man->DisableTexture();

    glMatrixMode(GL_MODELVIEW);

    glEnable(GL_TEXTURE_2D);
    m_tex_man->Begin();

    m_tool_helper->BeforeDrawImages();
    // The old preview shows the lowest numbered image on top, so do the same:
    for (int img = imgs - 1; img != -1; img--)
    {
        // only draw active images
        if (m_pano->getImage(img).getActive())
        {
            // the tools can cancel drawing of images.
            if (m_tool_helper->BeforeDrawImageNumber(img))
            {
                // the texture manager may need to call the display list
                // multiple times with blending, so we pass it the display list
                // rather than switching to the texture and then calling the
                // list ourselves.
                m_tex_man->DrawImage(img, m_mesh_man->GetDisplayList(img));
                m_tool_helper->AfterDrawImageNumber(img);
            }
        }
    }

    m_tex_man->End();
    // drawn things after the active image.
    m_tool_helper->AfterDrawImages();
    m_tex_man->DisableTexture();

    glPopMatrix();
}


vigra::Diff2D GLPlaneOverviewRenderer::Resize(int w, int h)
{
    width = w;
    height = h;
    glViewport(0, 0, width, height);
    // we use the view_state rather than the panorama to allow interactivity.
    HuginBase::PanoramaOptions *options = m_visualization_state->getViewState()->GetOptions();
    width_o = options->getWidth();
    height_o = options->getHeight();

    glMatrixMode(GL_PROJECTION);
    glLoadIdentity();

    double fov = m_visualization_state->getFOV();
    double fovy;
    if (h > w) {
        fovy = 2.0 * atan( tan(fov * M_PI / 360.0) * (float) h / (float) w) / M_PI * 180.0;
    } else {
        fovy = fov;
    }

	float ratio = 1.0* w / h;
//	aspect = ratio;
	glMatrixMode(GL_PROJECTION);
	glLoadIdentity();
	gluPerspective(fovy,ratio,1,1000000);

    m_visualization_state->SetVisibleArea(vigra::Rect2D(0, 0, options->getWidth(),
                                             options->getHeight()));

    double R = m_visualization_state->getR();
    double scrscale = (float) h / (2 * tan(fovy / 360.0 * M_PI) * R   * options->getWidth() / MeshManager::PlaneOverviewMeshInfo::scale);
    m_visualization_state->SetScale(scrscale);
//    m_visualization_state->SetGLScale(gl_scale);
    m_visualization_state->GetViewer()->MarkToolsDirty();

    return vigra::Diff2D(0,0);
}