gl_chart_canvas.cpp

Go to the documentation of this file.

/**************************************************************************
 *   Copyright (C) 2014 by David S. Register                               *
 *   Copyright (C) 2014 Sean D'Epagnier
 *                                                                         *
 *   This program 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 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 General Public License for more details.                          *
 *                                                                         *
 *   You should have received a copy of the GNU General Public License     *
 *   along with this program; if not, see <https://www.gnu.org/licenses/>. *
 ***************************************************************************/
 
// For compilers that support precompilation, includes "wx.h".
#include <wx/wxprec.h>
 
#ifndef WX_PRECOMP
#include <wx/wx.h>
#endif  // precompiled headers
 
#include "dychart.h"
 
#include <algorithm>
#include <stdint.h>
#include <vector>
 
#include <wx/brush.h>
#include <wx/colour.h>
#include <wx/dcmemory.h>
#include <wx/dynarray.h>
#include <wx/event.h>
#include <wx/font.h>
#include <wx/gdicmn.h>
#include <wx/glcanvas.h>
#include <wx/image.h>
#include <wx/jsonval.h>
#include <wx/log.h>
#include <wx/pen.h>
#include <wx/progdlg.h>
#include <wx/stopwatch.h>
#include <wx/string.h>
#include <wx/tokenzr.h>
#include <wx/utils.h>
#include <wx/window.h>
 
#include "model/base_platform.h"
#include "model/config_vars.h"
#include "model/gui_vars.h"
#include "model/own_ship.h"
#include "model/plugin_comm.h"
#include "model/route.h"
#include "model/routeman.h"
#include "model/track.h"
 
#include "ais.h"
#include "chartbase.h"
#include "chart_ctx_factory.h"
#include "chartdb.h"
#include "chartimg.h"
#include "chcanv.h"
#include "ch_info_win.h"
#include "cm93.h"  // for chart outline draw
#include "color_handler.h"
#include "compass.h"
#include "emboss_data.h"
#include "font_mgr.h"
#include "gl_chart_canvas.h"
#include "gl_tex_cache.h"
#include "gshhs.h"
#include "ienc_toolbar.h"
#include "lz4.h"
#include "mbtiles.h"
#include "mipmap/mipmap.h"
#include "mui_bar.h"
#include "navutil.h"
#include "ocpn_platform.h"
#include "piano.h"
#include "pluginmanager.h"
#include "quilt.h"
#include "rollover_win.h"
#include "route_gui.h"
#include "route_point_gui.h"
#include "s52plib.h"
#include "s57chart.h"  // for ArrayOfS57Obj
#include "s57_ocpn_utils.h"
#include "shapefile_basemap.h"
#include "tcmgr.h"
#include "toolbar.h"
#include "TexFont.h"
#include "thumbwin.h"
#include "toolbar.h"
#include "track_gui.h"
#include "viewport.h"
 
#ifdef USE_ANDROID_GLES2
#include <GLES2/gl2.h>
#include "linmath.h"
#include "shaders.h"
#endif
 
#ifdef __ANDROID__
#include "androidUTIL.h"
#elif defined(__WXQT__) || defined(__WXGTK__)
#include <GL/glx.h>
#endif
 
#ifndef GL_ETC1_RGB8_OES
#define GL_ETC1_RGB8_OES 0x8D64
#endif
 
#ifndef GL_DEPTH_STENCIL_ATTACHMENT
#define GL_DEPTH_STENCIL_ATTACHMENT 0x821A
#endif
 
#ifdef __WXMSW__
#define printf printf2
int __cdecl printf2(const char *format, ...);
#endif
 
#if defined(__ANDROID__)
#include "androidUTIL.h"
#elif defined(__WXQT__) || defined(__WXGTK__) || defined(FLATPAK)
#include <GL/glew.h>
#endif
 
#ifdef __ANDROID__
//  arm gcc compiler has a lot of trouble passing doubles as function aruments.
//  We don't really need double precision here, so fix with a (faster) macro.
extern "C" void glOrthof(float left, float right, float bottom, float top,
                         float near, float far);
#define glOrtho(a, b, c, d, e, f) \
  ;                               \
  glOrthof(a, b, c, d, e, f);
 
#endif
 
#ifdef USE_ANDROID_GLES2
#include <GLES2/gl2.h>
#endif
 
#if defined(USE_ANDROID_GLES2) || defined(ocpnUSE_GLSL)
#include "linmath.h"
#include "shaders.h"
#include "model/notification_manager.h"
#endif
 
#if defined(__UNIX__) && !defined(__WXOSX__)
// high resolution stopwatch for profiling
class OCPNStopWatch {
public:
  OCPNStopWatch() { Reset(); }
  void Reset() { clock_gettime(CLOCK_REALTIME, &tp); }
 
  double GetTime() {
    timespec tp_end;
    clock_gettime(CLOCK_REALTIME, &tp_end);
    return (tp_end.tv_sec - tp.tv_sec) * 1.e3 +
           (tp_end.tv_nsec - tp.tv_nsec) / 1.e6;
  }
 
private:
  timespec tp;
};
#endif
 
// extern GLenum g_texture_rectangle_format;
 
extern bool g_running;  
 
static wxColor s_regionColor;
static float g_GLMinCartographicLineWidth;
// MacOS has some missing parts:
#ifndef APIENTRY
#define APIENTRY
#endif
#ifndef APIENTRYP
#define APIENTRYP APIENTRY *
#endif
#ifndef GLAPI
#define GLAPI extern
#endif
 
#ifndef GL_COMPRESSED_RGB_FXT1_3DFX
#define GL_COMPRESSED_RGB_FXT1_3DFX 0x86B0
#endif
 
GLuint g_raster_format = GL_RGB;  
 
// OpenGL/GLES  bindings
PFNGLGENFRAMEBUFFERSEXTPROC s_glGenFramebuffers;
PFNGLGENRENDERBUFFERSEXTPROC s_glGenRenderbuffers;
PFNGLFRAMEBUFFERTEXTURE2DEXTPROC s_glFramebufferTexture2D;
PFNGLBINDFRAMEBUFFEREXTPROC s_glBindFramebuffer;
PFNGLFRAMEBUFFERRENDERBUFFEREXTPROC s_glFramebufferRenderbuffer;
PFNGLRENDERBUFFERSTORAGEEXTPROC s_glRenderbufferStorage;
PFNGLBINDRENDERBUFFEREXTPROC s_glBindRenderbuffer;
PFNGLCHECKFRAMEBUFFERSTATUSEXTPROC s_glCheckFramebufferStatus;
PFNGLDELETEFRAMEBUFFERSEXTPROC s_glDeleteFramebuffers;
PFNGLDELETERENDERBUFFERSEXTPROC s_glDeleteRenderbuffers;
 
PFNGLCOMPRESSEDTEXIMAGE2DPROC s_glCompressedTexImage2D;
PFNGLGETCOMPRESSEDTEXIMAGEPROC s_glGetCompressedTexImage;
 
//      Vertex Buffer Object (VBO) support
PFNGLGENBUFFERSPROC s_glGenBuffers;
PFNGLBINDBUFFERPROC s_glBindBuffer;
PFNGLBUFFERDATAPROC s_glBufferData;
PFNGLDELETEBUFFERSPROC s_glDeleteBuffers;
 
#ifndef USE_ANDROID_GLES2
// #define glDeleteFramebuffers(a, b) (s_glDeleteFramebuffers)(a, b);
// #define glDeleteRenderbuffers(a, b) (s_glDeleteRenderbuffers)(a, b);
#endif
 
typedef void(APIENTRYP PFNGLGETBUFFERPARAMETERIV)(GLenum target, GLenum value,
                                                  GLint *data);
PFNGLGETBUFFERPARAMETERIV s_glGetBufferParameteriv;
 
static bool b_timeGL;
static int panx, pany;
 
bool g_true_zoom;
 
#if defined(USE_ANDROID_GLES2) || defined(ocpnUSE_GLSL)
std::vector<GLfloat> s_tess_vertex_work;
static GLenum s_tess_mode;
static ViewPort s_tessVP;
#endif
 
bool glChartCanvas::s_b_useScissorTest;
bool glChartCanvas::s_b_useStencil;
bool glChartCanvas::s_b_useStencilAP;
bool glChartCanvas::s_b_useFBO;
 
#if 0
/* for debugging */
static void print_region(OCPNRegion &Region)
{
    OCPNRegionIterator upd ( Region );
    while ( upd.HaveRects() )
    {
        wxRect rect = upd.GetRect();
        printf("[(%d, %d) (%d, %d)] ", rect.x, rect.y, rect.width, rect.height);
        upd.NextRect();
    }
}
 
#endif
 
GLboolean QueryExtension(const char *extName) {
  /*
   ** Search for extName in the extensions string. Use of strstr()
   ** is not sufficient because extension names can be prefixes of
   ** other extension names. Could use strtok() but the constant
   ** string returned by glGetString might be in read-only memory.
   */
  char *p;
  char *end;
  int extNameLen;
 
  extNameLen = strlen(extName);
 
  p = (char *)glGetString(GL_EXTENSIONS);
  if (NULL == p) {
    return GL_FALSE;
  }
 
  end = p + strlen(p);
 
  while (p < end) {
    int n = strcspn(p, " ");
    if ((extNameLen == n) && (strncmp(extName, p, n) == 0)) {
      return GL_TRUE;
    }
    p += (n + 1);
  }
  return GL_FALSE;
}
 
int test_attribs[] = {WX_GL_RGBA, WX_GL_DOUBLEBUFFER, WX_GL_DEPTH_SIZE,
                      16,         WX_GL_STENCIL_SIZE, 8,
                      0};
 
glTestCanvas::glTestCanvas(wxWindow *parent)
    : wxGLCanvas(parent, wxID_ANY, test_attribs, wxDefaultPosition,
                 wxSize(2, 2)) {}
 
// This attribute set works OK with vesa software only OpenGL renderer
int attribs[] = {WX_GL_RGBA, WX_GL_DOUBLEBUFFER, WX_GL_DEPTH_SIZE,
                 16,         WX_GL_STENCIL_SIZE, 8,
                 0};
BEGIN_EVENT_TABLE(glChartCanvas, wxGLCanvas)
EVT_PAINT(glChartCanvas::OnPaint)
EVT_ACTIVATE(glChartCanvas::OnActivate)
EVT_SIZE(glChartCanvas::OnSize)
EVT_MOUSE_EVENTS(glChartCanvas::MouseEvent)
END_EVENT_TABLE()
 
glChartCanvas::glChartCanvas(wxWindow *parent, wxGLCanvas *share)
    : wxGLCanvas(parent, wxID_ANY, attribs, wxDefaultPosition, wxSize(256, 256),
                 wxFULL_REPAINT_ON_RESIZE | wxBG_STYLE_CUSTOM, "")
 
{
  m_pParentCanvas = dynamic_cast<ChartCanvas *>(parent);
 
  Init();
}
 
std::unordered_map<wxPenStyle, std::array<wxDash, 2>> glChartCanvas::dash_map =
    {
        {wxPENSTYLE_DOT, {1, 1}},
        {wxPENSTYLE_LONG_DASH, {5, 5}},
        {wxPENSTYLE_SHORT_DASH, {1, 5}},
        {wxPENSTYLE_DOT_DASH, {5, 1}},
};
 
void glChartCanvas::Init() {
  m_bsetup = false;
 
  //    m_pParentCanvas = dynamic_cast<ChartCanvas *>( GetParent() );
 
  SetBackgroundStyle(wxBG_STYLE_CUSTOM);  // on WXMSW, this prevents flashing
 
  m_cache_current_ch = NULL;
 
  m_b_paint_enable = true;
  m_in_glpaint = false;
 
  m_cache_tex[0] = m_cache_tex[1] = 0;
  m_cache_page = 0;
 
  m_b_BuiltFBO = false;
  m_b_DisableFBO = false;
 
  ownship_tex = 0;
  ownship_color = -1;
 
  m_piano_tex = 0;
 
  m_binPinch = false;
  m_binPan = false;
  m_bpinchGuard = false;
  m_binGesture = false;
  m_first_zout = false;
  g_true_zoom = true;
 
  b_timeGL = true;
  m_last_render_time = -1;
 
  m_LRUtime = 0;
 
  m_tideTex = 0;
  m_currentTex = 0;
 
  m_gldc.SetGLCanvas(this);
  m_gldc.SetDPIFactor(g_BasePlatform->GetDisplayDIPMult(GetParent()));
 
  m_displayScale = 1.0;
#if defined(__WXOSX__) || defined(__WXGTK3__)
  // Support scaled HDPI displays.
  m_displayScale = GetContentScaleFactor();
#endif
  m_pParentCanvas->VPoint.SetPixelScale(m_displayScale);
 
#ifdef __ANDROID__
  //  Create/connect a dynamic event handler slot for gesture and some timer
  //  events
  Connect(
      wxEVT_QT_PANGESTURE,
      (wxObjectEventFunction)(wxEventFunction)&glChartCanvas::OnEvtPanGesture,
      NULL, this);
 
  Connect(
      wxEVT_QT_PINCHGESTURE,
      (wxObjectEventFunction)(wxEventFunction)&glChartCanvas::OnEvtPinchGesture,
      NULL, this);
 
  Connect(GESTURE_EVENT_TIMER, wxEVT_TIMER,
          (wxObjectEventFunction)(wxEventFunction)&glChartCanvas::
              onGestureTimerEvent,
          NULL, this);
 
  Connect(GESTURE_FINISH_TIMER, wxEVT_TIMER,
          (wxObjectEventFunction)(wxEventFunction)&glChartCanvas::
              onGestureFinishTimerEvent,
          NULL, this);
 
  Connect(
      ZOOM_TIMER, wxEVT_TIMER,
      (wxObjectEventFunction)(wxEventFunction)&glChartCanvas::onZoomTimerEvent,
      NULL, this);
 
  m_gestureEeventTimer.SetOwner(this, GESTURE_EVENT_TIMER);
  m_gestureFinishTimer.SetOwner(this, GESTURE_FINISH_TIMER);
  zoomTimer.SetOwner(this, ZOOM_TIMER);
 
#ifdef USE_ANDROID_GLES2
//   Connect(
//       TEX_FADE_TIMER, wxEVT_TIMER,
//       (wxObjectEventFunction)(wxEventFunction)&glChartCanvas::onFadeTimerEvent,
//       NULL, this);
//   m_fadeTimer.SetOwner(this, TEX_FADE_TIMER);
#endif
 
#else
#ifdef HAVE_WX_GESTURE_EVENTS
 
  Connect(GESTURE_EVENT_TIMER, wxEVT_TIMER,
          (wxObjectEventFunction)(wxEventFunction)&glChartCanvas::
              onGestureTimerEvent,
          NULL, this);
 
  Connect(GESTURE_FINISH_TIMER, wxEVT_TIMER,
          (wxObjectEventFunction)(wxEventFunction)&glChartCanvas::
              onGestureFinishTimerEvent,
          NULL, this);
 
  Connect(
      ZOOM_TIMER, wxEVT_TIMER,
      (wxObjectEventFunction)(wxEventFunction)&glChartCanvas::onZoomTimerEvent,
      NULL, this);
 
  m_gestureEeventTimer.SetOwner(this, GESTURE_EVENT_TIMER);
  m_gestureFinishTimer.SetOwner(this, GESTURE_FINISH_TIMER);
  zoomTimer.SetOwner(this, ZOOM_TIMER);
  m_zoom_inc = 1.0;
#endif
#endif
 
  m_bgestureGuard = false;
  m_total_zoom_val = 1.0;
  m_step_zoom_val = 1.0;
 
//  Gesture support for platforms other than Android
#ifdef HAVE_WX_GESTURE_EVENTS
  if (!EnableTouchEvents(wxTOUCH_ZOOM_GESTURE | wxTOUCH_PRESS_GESTURES)) {
    wxLogError("Failed to enable touch events");
  }
 
  // Bind(wxEVT_GESTURE_PAN, &glChartCanvas::OnEvtPanGesture, this);
  // Connect(
  //     wxEVT_GESTURE_PAN,
  //     (wxObjectEventFunction)(wxEventFunction)&glChartCanvas::OnEvtPanGesture,
  //     NULL, this);
 
  Bind(wxEVT_GESTURE_ZOOM, &glChartCanvas::OnEvtZoomGesture, this);
 
  Bind(wxEVT_LONG_PRESS, &ChartCanvas::OnLongPress, m_pParentCanvas);
  Bind(wxEVT_PRESS_AND_TAP, &ChartCanvas::OnPressAndTap, m_pParentCanvas);
 
  Bind(wxEVT_RIGHT_UP, &ChartCanvas::OnRightUp, m_pParentCanvas);
  Bind(wxEVT_RIGHT_DOWN, &ChartCanvas::OnRightDown, m_pParentCanvas);
 
  Bind(wxEVT_LEFT_UP, &ChartCanvas::OnLeftUp, m_pParentCanvas);
  Bind(wxEVT_LEFT_DOWN, &ChartCanvas::OnLeftDown, m_pParentCanvas);
 
  Bind(wxEVT_MOUSEWHEEL, &ChartCanvas::OnWheel, m_pParentCanvas);
  Bind(wxEVT_MOTION, &ChartCanvas::OnMotion, m_pParentCanvas);
#endif /* HAVE_WX_GESTURE_EVENTS */
 
  if (!g_glTextureManager) g_glTextureManager = new glTextureManager;
}
 
glChartCanvas::~glChartCanvas() {
#ifdef __ANDROID__
  unloadShaders();
#endif
}
 
int glChartCanvas::GetCanvasIndex() { return m_pParentCanvas->m_canvasIndex; }
 
void glChartCanvas::FlushFBO() {
  if (m_bsetup) BuildFBO();
}
 
void glChartCanvas::OnActivate(wxActivateEvent &event) {
  m_pParentCanvas->OnActivate(event);
}
 
void glChartCanvas::OnSize(wxSizeEvent &event) {
  //  Set the shader viewport transform matrix
  //  Always safe, simple math based on parent VP
  ViewPort *vp = m_pParentCanvas->GetpVP();
  vp->SetVPTransformMatrix();
 
#ifndef __ANDROID__
  if (!m_bsetup) return;
#endif
  if (!IsShown()) return;
 
  SetCurrent(*m_pcontext);
 
  if (!g_bopengl) {
    // Invoked immediately after user has disabled OpenGL.
    SetSize(GetSize().x, GetSize().y);
    event.Skip();
    return;
  }
 
  // this is also necessary to update the context on some platforms
  // OnSize can be called with a different OpenGL context (when a plugin uses a
  // different GL context).
  if (m_pcontext && IsShown()) {
    SetCurrent(*m_pcontext);
  }
 
  wxLogDebug("BuildFBO 3");
  BuildFBO();
}
 
void glChartCanvas::MouseEvent(wxMouseEvent &event) {
  if (m_pParentCanvas->MouseEventOverlayWindows(event)) return;
 
#ifndef __ANDROID__
  if (m_pParentCanvas->MouseEventSetup(event))
    return;  // handled, no further action required
 
  bool obj_proc = m_pParentCanvas->MouseEventProcessObjects(event);
 
  if (!obj_proc && !m_pParentCanvas->singleClickEventIsValid)
    m_pParentCanvas->MouseEventProcessCanvas(event);
 
  if (!g_btouch) m_pParentCanvas->SetCanvasCursor(event);
 
#else
 
  if (m_bgestureGuard) {
    m_pParentCanvas->r_rband.x = 0;  // turn off rubberband temporarily
 
    // Sometimes we get a Gesture Pan start on a simple tap operation.
    // When this happens, we usually get no Gesture Finished event.
    // So, we need to process the next LeftUp event normally, to handle things
    // like Measure and Route Create.
 
    // Allow LeftUp() event through if the pan action is very small
    //  Otherwise, drop the LeftUp() event, since it is not wanted for a Pan
    //  Gesture.
    if (event.LeftUp()) {
      // qDebug() << panx << pany;
      if ((abs(panx) > 2) || (abs(pany) > 2)) {
        return;
      } else {  // Cancel the in=process Gesture state
        m_gestureEeventTimer.Start(10, wxTIMER_ONE_SHOT);  // Short Circuit
      }
    } else
      return;
  }
 
  if (m_pParentCanvas->MouseEventSetup(event, false)) {
    if (!event.LeftDClick()) {
      return;  // handled, no further action required
    }
  }
 
  if (m_binPan && event.RightDown()) {
    qDebug() << "Skip right on pan";
    return;
  } else {
    bool obj_proc = m_pParentCanvas->MouseEventProcessObjects(event);
 
    if (!obj_proc && !m_pParentCanvas->singleClickEventIsValid) {
      if (!m_bgestureGuard)
        m_pParentCanvas->MouseEventProcessCanvas(
            event);  // This is where a physical mouse gets processed, if
                     // detected
    }
  }
 
#endif
}
 
#ifndef GL_MAX_RENDERBUFFER_SIZE
#define GL_MAX_RENDERBUFFER_SIZE 0x84E8
#endif
 
#ifndef USE_ANDROID_GLES2
bool glChartCanvas::buildFBOSize(int fboSize) {
  bool retVal = true;
  if (IsShown()) SetCurrent(*m_pcontext);
 
  if (m_b_BuiltFBO) {
    glDeleteTextures(2, m_cache_tex);
    glDeleteFramebuffers(1, &m_fb0);
    glDeleteRenderbuffers(1, &m_renderbuffer);
    m_b_BuiltFBO = false;
  }
 
  if (m_b_DisableFBO) return false;
 
#ifdef __ANDROID__
  // We use the smallest possible (POT) FBO
  int rb_x = GetSize().x;
  int rb_y = GetSize().y;
  int i = 1;
  while (i < rb_x) i <<= 1;
  rb_x = i;
 
  i = 1;
  while (i < rb_y) i <<= 1;
  rb_y = i;
 
  m_cache_tex_x = wxMax(rb_x, rb_y);
  m_cache_tex_y = wxMax(rb_x, rb_y);
 
#else
  m_cache_tex_x = GetSize().x * m_displayScale;
  m_cache_tex_y = GetSize().y * m_displayScale;
#endif
 
  int err = GL_NO_ERROR;
  GLint params;
  glGetIntegerv(GL_MAX_RENDERBUFFER_SIZE, &params);
 
  err = glGetError();
  if (err == GL_INVALID_ENUM) {
    glGetIntegerv(GL_MAX_TEXTURE_SIZE, &params);
    err = glGetError();
  }
 
  if (err == GL_NO_ERROR) {
    if (fboSize > params) {
      wxLogMessage(
          "    OpenGL-> Requested Framebuffer size exceeds "
          "GL_MAX_RENDERBUFFER_SIZE");
      return false;
    }
  }
 
  glGenFramebuffers(1, &m_fb0);
  err = glGetError();
  if (err) {
    wxString msg;
    msg.Printf("    OpenGL-> Framebuffer GenFramebuffers error:  %08X", err);
    wxLogMessage(msg);
    retVal = false;
  }
 
  glGenRenderbuffers(1, &m_renderbuffer);
  err = glGetError();
  if (err) {
    wxString msg;
    msg.Printf("    OpenGL-> Framebuffer GenRenderbuffers error:  %08X", err);
    wxLogMessage(msg);
    retVal = false;
  }
 
  glBindFramebuffer(GL_FRAMEBUFFER_EXT, m_fb0);
  err = glGetError();
  if (err) {
    wxString msg;
    msg.Printf("    OpenGL-> Framebuffer BindFramebuffers error:  %08X", err);
    wxLogMessage(msg);
    retVal = false;
  }
 
  // initialize color textures
  glGenTextures(2, m_cache_tex);
  for (int i = 0; i < 2; i++) {
    glBindTexture(g_texture_rectangle_format, m_cache_tex[i]);
    glTexParameterf(g_texture_rectangle_format, GL_TEXTURE_MIN_FILTER,
                    GL_NEAREST);
    glTexParameteri(g_texture_rectangle_format, GL_TEXTURE_MAG_FILTER,
                    GL_NEAREST);
    glTexImage2D(g_texture_rectangle_format, 0, GL_RGBA, m_cache_tex_x,
                 m_cache_tex_y, 0, GL_RGBA, GL_UNSIGNED_BYTE, NULL);
  }
 
  glBindRenderbuffer(GL_RENDERBUFFER_EXT, m_renderbuffer);
 
  if (m_b_useFBOStencil) {
    // initialize composite depth/stencil renderbuffer
    glRenderbufferStorage(GL_RENDERBUFFER_EXT, GL_DEPTH24_STENCIL8_EXT,
                          m_cache_tex_x, m_cache_tex_y);
 
    int err = glGetError();
    if (err) {
      wxString msg;
      msg.Printf("    OpenGL-> glRenderbufferStorage error:  %08X", err);
      wxLogMessage(msg);
    }
 
    glFramebufferRenderbuffer(GL_FRAMEBUFFER_EXT, GL_DEPTH_ATTACHMENT_EXT,
                              GL_RENDERBUFFER_EXT, m_renderbuffer);
    err = glGetError();
    if (err) {
      wxString msg;
      msg.Printf("    OpenGL-> glFramebufferRenderbuffer depth error:  %08X",
                 err);
      wxLogMessage(msg);
    }
 
    glFramebufferRenderbuffer(GL_FRAMEBUFFER_EXT, GL_STENCIL_ATTACHMENT_EXT,
                              GL_RENDERBUFFER_EXT, m_renderbuffer);
    err = glGetError();
    if (err) {
      wxString msg;
      msg.Printf("    OpenGL-> glFramebufferRenderbuffer stencil error:  %08X",
                 err);
      wxLogMessage(msg);
    }
 
  } else {
    GLenum depth_format = GL_DEPTH_COMPONENT24;
 
    //      Need to check for availability of 24 bit depth buffer extension on
    //      GLES
#ifdef ocpnUSE_GLES
    if (!QueryExtension("GL_OES_depth24")) depth_format = GL_DEPTH_COMPONENT16;
#endif
 
    // initialize depth renderbuffer
    glRenderbufferStorage(GL_RENDERBUFFER_EXT, depth_format, m_cache_tex_x,
                          m_cache_tex_y);
    int err = glGetError();
    if (err) {
      wxString msg;
      msg.Printf("    OpenGL-> Framebuffer Depth Buffer Storage error:  %08X",
                 err);
      wxLogMessage(msg);
      retVal = false;
    }
 
    glFramebufferRenderbuffer(GL_FRAMEBUFFER_EXT, GL_DEPTH_ATTACHMENT_EXT,
                              GL_RENDERBUFFER_EXT, m_renderbuffer);
 
    err = glGetError();
    if (err) {
      wxString msg;
      msg.Printf("    OpenGL-> Framebuffer Depth Buffer Attach error:  %08X",
                 err);
      wxLogMessage(msg);
      retVal = false;
    }
  }
 
  glBindTexture(GL_TEXTURE_2D, 0);
  glBindFramebuffer(GL_FRAMEBUFFER_EXT, 0);
 
  // Check framebuffer completeness at the end of initialization.
  glBindFramebuffer(GL_FRAMEBUFFER_EXT, m_fb0);
 
  glFramebufferTexture2D(GL_FRAMEBUFFER_EXT, GL_COLOR_ATTACHMENT0_EXT,
                         g_texture_rectangle_format, m_cache_tex[0], 0);
 
  GLenum fb_status = glCheckFramebufferStatus(GL_FRAMEBUFFER_EXT);
 
  glBindFramebuffer(GL_FRAMEBUFFER_EXT, 0);
 
  if (fb_status != GL_FRAMEBUFFER_COMPLETE_EXT) {
    wxString msg;
    msg.Printf("    OpenGL-> buildFBOSize->Framebuffer Incomplete:  %08X",
               fb_status);
    wxLogMessage(msg);
    retVal = false;
  }
 
  return retVal;
}
#endif
 
#ifdef USE_ANDROID_GLES2
bool glChartCanvas::buildFBOSize(int fboSize) {
  bool retVal = true;
 
  // We use the smallest possible (POT) FBO
  int rb_x = GetSize().x;
  int rb_y = GetSize().y;
  int i = 1;
  while (i < rb_x) i <<= 1;
  rb_x = i;
 
  i = 1;
  while (i < rb_y) i <<= 1;
  rb_y = i;
 
  m_cache_tex_x = wxMax(rb_x, rb_y);
  m_cache_tex_y = wxMax(rb_x, rb_y);
 
  // qDebug() << "FBO Size: " << GetSize().x << GetSize().y << m_cache_tex_x;
 
  int err = GL_NO_ERROR;
  GLint params;
  glGetIntegerv(GL_MAX_RENDERBUFFER_SIZE, &params);
 
  err = glGetError();
  if (err == GL_INVALID_ENUM) {
    glGetIntegerv(GL_MAX_TEXTURE_SIZE, &params);
    err = glGetError();
  }
 
  if (err == GL_NO_ERROR) {
    if (fboSize > params) {
      wxLogMessage(
          "    OpenGL-> Requested Framebuffer size exceeds "
          "GL_MAX_RENDERBUFFER_SIZE");
      return false;
    }
  }
 
  glGenFramebuffers(1, &m_fb0);
  err = glGetError();
  if (err) {
    wxString msg;
    msg.Printf("    OpenGL-> Framebuffer GenFramebuffers error:  %08X", err);
    wxLogMessage(msg);
    retVal = false;
  }
 
  glGenRenderbuffers(1, &m_renderbuffer);
  err = glGetError();
  if (err) {
    wxString msg;
    msg.Printf("    OpenGL-> Framebuffer GenRenderbuffers error:  %08X", err);
    wxLogMessage(msg);
    retVal = false;
  }
 
  glBindFramebuffer(GL_FRAMEBUFFER, m_fb0);
  err = glGetError();
  if (err) {
    wxString msg;
    msg.Printf("    OpenGL-> Framebuffer BindFramebuffers error:  %08X", err);
    wxLogMessage(msg);
    retVal = false;
  }
 
  // initialize color textures
  glGenTextures(2, m_cache_tex);
  for (int i = 0; i < 2; i++) {
    glBindTexture(g_texture_rectangle_format, m_cache_tex[i]);
    glTexParameterf(g_texture_rectangle_format, GL_TEXTURE_MIN_FILTER,
                    GL_NEAREST);
    glTexParameteri(g_texture_rectangle_format, GL_TEXTURE_MAG_FILTER,
                    GL_NEAREST);
    glTexImage2D(g_texture_rectangle_format, 0, GL_RGBA, m_cache_tex_x,
                 m_cache_tex_y, 0, GL_RGBA, GL_UNSIGNED_BYTE, NULL);
  }
 
  glBindRenderbuffer(GL_RENDERBUFFER, m_renderbuffer);
 
  // initialize composite depth/stencil renderbuffer
  glRenderbufferStorage(GL_RENDERBUFFER, GL_DEPTH24_STENCIL8_OES, m_cache_tex_x,
                        m_cache_tex_y);
 
  err = glGetError();
  if (err) {
    wxString msg;
    msg.Printf("    OpenGL-> glRenderbufferStorage error:  %08X", err);
    wxLogMessage(msg);
  }
 
  glFramebufferRenderbuffer(GL_FRAMEBUFFER, GL_DEPTH_ATTACHMENT,
                            GL_RENDERBUFFER, m_renderbuffer);
  err = glGetError();
  if (err) {
    wxString msg;
    msg.Printf("    OpenGL-> glFramebufferRenderbuffer depth error:  %08X",
               err);
    wxLogMessage(msg);
  }
 
  glFramebufferRenderbuffer(GL_FRAMEBUFFER, GL_STENCIL_ATTACHMENT,
                            GL_RENDERBUFFER, m_renderbuffer);
  err = glGetError();
  if (err) {
    wxString msg;
    msg.Printf("    OpenGL-> glFramebufferRenderbuffer stencil error:  %08X",
               err);
    wxLogMessage(msg);
  }
 
  glBindTexture(GL_TEXTURE_2D, 0);
  glBindFramebuffer(GL_FRAMEBUFFER, 0);
 
  // Check framebuffer completeness at the end of initialization.
  glBindFramebuffer(GL_FRAMEBUFFER, m_fb0);
 
  glFramebufferTexture2D(GL_FRAMEBUFFER, GL_COLOR_ATTACHMENT0,
                         g_texture_rectangle_format, m_cache_tex[0], 0);
 
  GLenum fb_status = glCheckFramebufferStatus(GL_FRAMEBUFFER_EXT);
 
  glBindFramebuffer(GL_FRAMEBUFFER, 0);
 
  if (fb_status != GL_FRAMEBUFFER_COMPLETE) {
    wxString msg;
    msg.Printf("    OpenGL-> buildFBOSize->Framebuffer Incomplete:  %08X %08X",
               fb_status);
    wxLogMessage(msg);
    retVal = false;
  }
 
  return retVal;
}
#endif
 
void glChartCanvas::BuildFBO() {
  if (g_b_needFinish) glFinish();
 
  if (m_b_BuiltFBO) {
    // return;
    glDeleteTextures(2, m_cache_tex);
    glDeleteFramebuffers(1, &m_fb0);
    glDeleteRenderbuffers(1, &m_renderbuffer);
    m_b_BuiltFBO = false;
  }
 
  if (m_b_DisableFBO) return;
 
  //    int initialSize = 2048;
  int gl_width, gl_height;
  m_pParentCanvas->GetClientSize(&gl_width, &gl_height);
  int initialSize = NextPow2(gl_width * m_displayScale);
 
#ifdef __ANDROID__
  //  Some low mem-spec devices have trouble with 2048 FBO size.
  //  Detect here, and choose 1024 size instead
  wxString info = androidGetDeviceInfo();
 
  if (wxNOT_FOUND != info.Find("GT-S6312")) initialSize = 1024;
#endif
 
  if (!buildFBOSize(initialSize)) {
    glDeleteTextures(2, m_cache_tex);
    glDeleteFramebuffers(1, &m_fb0);
    glDeleteRenderbuffers(1, &m_renderbuffer);
 
    if (!buildFBOSize(1024)) {
      wxLogMessage("BuildFBO C");
 
      m_b_DisableFBO = true;
      wxLogMessage("OpenGL-> FBO Framebuffer unavailable");
      m_b_BuiltFBO = false;
 
      return;
    }
  }
 
  //  All OK
 
  wxString msg;
  msg.Printf("OpenGL-> Framebuffer OK, size = %d", m_cache_tex_x);
  wxLogDebug(msg);
 
  /* invalidate cache */
  Invalidate();
 
  // glClear(GL_COLOR_BUFFER_BIT);
  m_b_BuiltFBO = true;
 
  return;
}
 
void glChartCanvas::SetupOpenGL() {
  if (IsShown()) SetCurrent(*m_pcontext);
 
  char *str = (char *)glGetString(GL_RENDERER);
  if (str == NULL) {
    // perhaps we should edit the config and turn off opengl now
    wxLogMessage("Failed to initialize OpenGL");
    exit(1);
  }
 
  char render_string[80];
  strncpy(render_string, str, 79);
  m_renderer = wxString(render_string, wxConvUTF8);
 
  wxString msg;
  if (g_bSoftwareGL) msg.Printf("OpenGL-> Software OpenGL");
  msg.Printf("OpenGL-> Renderer String: ");
  msg += m_renderer;
  wxLogMessage(msg);
 
  if (ps52plib) ps52plib->SetGLRendererString(m_renderer);
 
  char version_string[80];
  strncpy(version_string, (char *)glGetString(GL_VERSION), 79);
  msg.Printf("OpenGL-> Version reported:  ");
  m_version = wxString(version_string, wxConvUTF8);
  msg += m_version;
  wxLogMessage(msg);
 
  char GLSL_version_string[80];
  strncpy(GLSL_version_string, (char *)glGetString(GL_SHADING_LANGUAGE_VERSION),
          79);
  msg.Printf("OpenGL-> GLSL Version reported:  ");
  m_GLSLversion = wxString(GLSL_version_string, wxConvUTF8);
  msg += m_GLSLversion;
  wxLogMessage(msg);
 
#ifndef __ANDROID__
#ifndef __WXOSX__
  GLenum err = glewInit();
#ifdef GLEW_ERROR_NO_GLX_DISPLAY
  if (GLEW_OK != err && GLEW_ERROR_NO_GLX_DISPLAY != err)
#else
  if (GLEW_OK != err)
#endif
  {
    printf("GLEW init failed: %s\n", glewGetErrorString(err));
    exit(1);
  } else {
    wxLogMessage("GLEW init success!n");
  }
#endif
#endif
 
  const GLubyte *ext_str = glGetString(GL_EXTENSIONS);
  m_extensions = wxString((const char *)ext_str, wxConvUTF8);
 
  //  Set the minimum line width
  GLint parms[2];
#ifndef USE_ANDROID_GLES2
  glGetIntegerv(GL_SMOOTH_LINE_WIDTH_RANGE, &parms[0]);
#else
  glGetIntegerv(GL_ALIASED_LINE_WIDTH_RANGE, &parms[0]);
#endif
  g_GLMinSymbolLineWidth = wxMax(parms[0], 1);
  g_GLMinCartographicLineWidth = wxMax(parms[0], 1);
 
  //    Some GL renderers do a poor job of Anti-aliasing very narrow line
  //    widths. This is most evident on rendered symbols which have horizontal
  //    or vertical line segments Detect this case, and adjust the render
  //    parameters.
 
  if (m_renderer.Upper().Find("MESA") != wxNOT_FOUND) {
    GLfloat parf;
    glGetFloatv(GL_SMOOTH_LINE_WIDTH_GRANULARITY, &parf);
 
    g_GLMinSymbolLineWidth = wxMax(((float)parms[0] + parf), 1);
  }
 
  s_b_useScissorTest = true;
  // the radeon x600 driver has buggy scissor test
  if (GetRendererString().Find("RADEON X600") != wxNOT_FOUND)
    s_b_useScissorTest = false;
 
  if (GetRendererString().Find("GeForce") != wxNOT_FOUND)  // GeForce GTX 1070
    s_b_useScissorTest = false;
 
  bool bad_stencil_code = false;
 
  //      And for the lousy Unichrome drivers, too
  if (GetRendererString().Find("UniChrome") != wxNOT_FOUND)
    bad_stencil_code = true;
 
  //      And for the lousy Mali drivers, too
  if (GetRendererString().Find("Mali") != wxNOT_FOUND) bad_stencil_code = true;
 
  //      Stencil buffer test
  glEnable(GL_STENCIL_TEST);
  GLboolean stencil = glIsEnabled(GL_STENCIL_TEST);
  int sb;
  glGetIntegerv(GL_STENCIL_BITS, &sb);
  //        printf("Stencil Buffer Available: %d\nStencil bits: %d\n", stencil,
  //        sb);
  glDisable(GL_STENCIL_TEST);
 
  s_b_useStencil = false;
  if (stencil && (sb == 8)) s_b_useStencil = true;
 
  if (QueryExtension("GL_ARB_texture_non_power_of_two"))
    g_texture_rectangle_format = GL_TEXTURE_2D;
  else if (QueryExtension("GL_OES_texture_npot"))
    g_texture_rectangle_format = GL_TEXTURE_2D;
  else if (QueryExtension("GL_ARB_texture_rectangle"))
    g_texture_rectangle_format = GL_TEXTURE_RECTANGLE_ARB;
  wxLogMessage(wxString::Format("OpenGL-> Texture rectangle format: %x",
                                g_texture_rectangle_format));
 
#ifdef __ANDROID__
  g_texture_rectangle_format = GL_TEXTURE_2D;
#endif
 
  // VBO??
  g_b_EnableVBO = true;
 
  if (g_b_EnableVBO)
    wxLogMessage("OpenGL-> Using Vertexbuffer Objects");
  else
    wxLogMessage("OpenGL-> Vertexbuffer Objects unavailable");
 
    //      Can we use the stencil buffer in a FBO?
#ifdef ocpnUSE_GLES
  m_b_useFBOStencil = QueryExtension("GL_OES_packed_depth_stencil");
#else
  m_b_useFBOStencil = QueryExtension("GL_EXT_packed_depth_stencil") == GL_TRUE;
#endif
 
#ifndef USE_ANDROID_GLES2
  //  On Intel Graphics platforms, don't use stencil buffer at all
  if (bad_stencil_code) s_b_useStencil = false;
#endif
 
  g_GLOptions.m_bUseCanvasPanning = false;
 
  // TODO
  //  Temporarily disable FBO on Windows, pending implementation of MSAA to
  //  buffers
#ifdef __WXMSW__
  // m_b_DisableFBO = true;
#endif
 
  // Accelerated pan is not used for MacOS Retina display
  // So there is no advantage to using FBO
  if (m_displayScale > 1) m_b_DisableFBO = true;
 
  //      Maybe build FBO(s)
  BuildFBO();
 
#ifndef __ANDROID__
 
  if (m_b_BuiltFBO) {
    // Check framebuffer completeness at the end of initialization.
    glBindFramebuffer(GL_FRAMEBUFFER_EXT, m_fb0);
 
    glFramebufferTexture2D(GL_FRAMEBUFFER_EXT, GL_COLOR_ATTACHMENT0_EXT,
                           g_texture_rectangle_format, m_cache_tex[0], 0);
 
    GLenum fb_status = glCheckFramebufferStatus(GL_FRAMEBUFFER_EXT);
    glBindFramebuffer(GL_FRAMEBUFFER_EXT, 0);
 
    if (fb_status != GL_FRAMEBUFFER_COMPLETE_EXT) {
      wxString msg;
      msg.Printf("    OpenGL-> Framebuffer Incomplete:  %08X", fb_status);
      wxLogMessage(msg);
      m_b_DisableFBO = true;
      BuildFBO();
    }
  }
#endif
 
  if (m_b_BuiltFBO && !m_b_useFBOStencil) s_b_useStencil = false;
 
  //  If stencil seems to be a problem, force use of depth buffer clipping for
  //  Area Patterns
  s_b_useStencilAP = s_b_useStencil & !bad_stencil_code;
 
#ifdef USE_ANDROID_GLES2
  s_b_useStencilAP = s_b_useStencil;  // required for GLES2 platform
#endif
 
  //  Check and determine if GLSL is to be used
  m_bUseGLSL = true;
 
  if (m_b_BuiltFBO) {
    wxLogMessage("OpenGL-> Using Framebuffer Objects");
 
    if (m_b_useFBOStencil)
      wxLogMessage("OpenGL-> Using FBO Stencil buffer");
    else
      wxLogMessage("OpenGL-> FBO Stencil buffer unavailable");
  } else
    wxLogMessage("OpenGL-> Framebuffer Objects unavailable");
 
  if (s_b_useStencil)
    wxLogMessage("OpenGL-> Using Stencil buffer clipping");
  else
    wxLogMessage("OpenGL-> Using Depth buffer clipping");
 
  if (s_b_useScissorTest && s_b_useStencil)
    wxLogMessage("OpenGL-> Using Scissor Clipping");
 
  /* we upload non-aligned memory */
  glPixelStorei(GL_UNPACK_ALIGNMENT, 1);
 
  MipMap_ResolveRoutines();
  SetupCompression();
 
  wxString lwmsg;
  lwmsg.Printf("OpenGL-> Minimum cartographic line width: %4.1f",
               g_GLMinCartographicLineWidth);
  wxLogMessage(lwmsg);
  lwmsg.Printf("OpenGL-> Minimum symbol line width: %4.1f",
               g_GLMinSymbolLineWidth);
  wxLogMessage(lwmsg);
 
  if (!g_bGLexpert)
    g_GLOptions.m_bUseAcceleratedPanning = !m_b_DisableFBO && m_b_BuiltFBO;
 
#ifdef USE_ANDROID_GLES2
  g_GLOptions.m_bUseAcceleratedPanning = true;
#endif
 
  if (1)  // for now upload all levels
  {
    int max_level = 0;
    int tex_dim = g_GLOptions.m_iTextureDimension;
    for (int dim = tex_dim; dim > 0; dim /= 2) max_level++;
    g_mipmap_max_level = max_level - 1;
  }
 
  //  Android, even though using GLES, does not require all levels.
#ifdef __ANDROID__
  g_mipmap_max_level = 4;
#endif
 
  s_b_useFBO = m_b_BuiltFBO;
 
  //  Inform the S52 PLIB of options determined
  if (ps52plib)
    ps52plib->SetGLOptions(
        s_b_useStencil, s_b_useStencilAP, s_b_useScissorTest, s_b_useFBO,
        g_b_EnableVBO, g_texture_rectangle_format, g_GLMinCartographicLineWidth,
        g_GLMinSymbolLineWidth);
 
  m_bsetup = true;
 
  SendJSONConfigMessage();
}
 
void glChartCanvas::SendJSONConfigMessage() {
  if (g_pi_manager) {
    wxJSONValue v;
    v["setupComplete"] = m_bsetup;
    v["useStencil"] = s_b_useStencil;
    v["useStencilAP"] = s_b_useStencilAP;
    v["useScissorTest"] = s_b_useScissorTest;
    v["useFBO"] = s_b_useFBO;
    v["useVBO"] = g_b_EnableVBO;
    v["TextureRectangleFormat"] = g_texture_rectangle_format;
    wxString msg_id("OCPN_OPENGL_CONFIG");
    SendJSONMessageToAllPlugins(msg_id, v);
  }
}
void glChartCanvas::SetupCompression() {
  int dim = g_GLOptions.m_iTextureDimension;
 
#ifdef __WXMSW__
  if (!::IsProcessorFeaturePresent(PF_XMMI64_INSTRUCTIONS_AVAILABLE)) {
    wxLogMessage("OpenGL-> SSE2 Instruction set not available");
    goto no_compression;
  }
#endif
 
  g_uncompressed_tile_size = dim * dim * 4;  // stored as 32bpp in vram
  if (!g_GLOptions.m_bTextureCompression) goto no_compression;
 
  g_raster_format = GL_RGB;
 
  // On GLES, we prefer OES_ETC1 compression, if available
#ifdef ocpnUSE_GLES
  if (QueryExtension("GL_OES_compressed_ETC1_RGB8_texture")) {
    g_raster_format = GL_ETC1_RGB8_OES;
 
    wxLogMessage("OpenGL-> Using oes etc1 compression");
  }
#endif
 
  if (GL_RGB == g_raster_format) {
    /* because s3tc is patented, many foss drivers disable
       support by default, however the extension dxt1 allows
       us to load this texture type which is enough because we
       compress in software using libsquish for superior quality anyway */
 
    if ((QueryExtension("GL_EXT_texture_compression_s3tc") ||
         QueryExtension("GL_EXT_texture_compression_dxt1"))) {
      /* buggy opensource nvidia driver, renders incorrectly,
         workaround is to use format with alpha... */
      if (GetRendererString().Find("Gallium") != wxNOT_FOUND &&
          GetRendererString().Find("NV") != wxNOT_FOUND)
        g_raster_format = GL_COMPRESSED_RGBA_S3TC_DXT1_EXT;
      else
        g_raster_format = GL_COMPRESSED_RGB_S3TC_DXT1_EXT;
 
      wxLogMessage("OpenGL-> Using s3tc dxt1 compression");
    } else if (QueryExtension("GL_3DFX_texture_compression_FXT1")) {
      g_raster_format = GL_COMPRESSED_RGB_FXT1_3DFX;
 
      wxLogMessage("OpenGL-> Using 3dfx fxt1 compression");
    } else {
      wxLogMessage("OpenGL-> No Useable compression format found");
      goto no_compression;
    }
  }
 
#ifdef ocpnUSE_GLES            /* gles doesn't have GetTexLevelParameter */
  g_tile_size = 512 * 512 / 2; /* 4bpp */
#else
  /* determine compressed size of a level 0 single tile */
  GLuint texture;
  glGenTextures(1, &texture);
  glBindTexture(GL_TEXTURE_2D, texture);
  glTexImage2D(GL_TEXTURE_2D, 0, g_raster_format, dim, dim, 0, GL_RGB,
               GL_UNSIGNED_BYTE, NULL);
  glGetTexLevelParameteriv(GL_TEXTURE_2D, 0, GL_TEXTURE_COMPRESSED_IMAGE_SIZE,
                           &g_tile_size);
  glDeleteTextures(1, &texture);
#endif
 
  /* disable texture compression if the tile size is 0 */
  if (g_tile_size == 0) goto no_compression;
 
  wxLogMessage(wxString::Format("OpenGL-> Compressed tile size: %dkb (%d:1)",
                                g_tile_size / 1024,
                                g_uncompressed_tile_size / g_tile_size));
  return;
 
no_compression:
  g_GLOptions.m_bTextureCompression = false;
 
  g_tile_size = g_uncompressed_tile_size;
  wxLogMessage(wxString::Format("OpenGL-> Not Using compression"));
}
 
void glChartCanvas::OnPaint(wxPaintEvent &event) {
  wxPaintDC dc(this);
  if (!m_pcontext) return;
  if (ChartData->IsBusy()) return;
 
  Show(g_bopengl);
  if (!g_bopengl) {
    event.Skip();
    return;
  }
 
  if (IsShown()) SetCurrent(*m_pcontext);
 
  if (!m_bsetup) {
    SetupOpenGL();
 
    if (ps52plib) ps52plib->FlushSymbolCaches(ChartCtxFactory());
 
    m_bsetup = true;
    //        g_bDebugOGL = true;
  }
 
  //  Paint updates may have been externally disabled (temporarily, to avoid
  //  Yield() recursion performance loss)
  if (!m_b_paint_enable) return;
  //      Recursion test, sometimes seen on GTK systems when wxBusyCursor is
  //      activated
  if (m_in_glpaint) return;
 
  //  If necessary, reconfigure the S52 PLIB
  m_pParentCanvas->UpdateCanvasS52PLIBConfig();
 
  //     if( m_pParentCanvas->VPoint.b_quilt ){          // quilted
  //         if( !m_pParentCanvas->m_pQuilt ||
  //         !m_pParentCanvas->m_pQuilt->IsComposed() )
  //             return;  // not ready
  //
  //             if(m_pParentCanvas->m_pQuilt->IsQuiltVector()){
  //                 if(ps52plib->GetStateHash() !=
  //                 m_pParentCanvas->m_s52StateHash){
  //                     m_pParentCanvas->UpdateS52State();
  //                     m_pParentCanvas->m_s52StateHash =
  //                     ps52plib->GetStateHash();
  //                 }
  //             }
  //     }
 
  m_in_glpaint++;
  Render();
  m_in_glpaint--;
}
 
//   These routines allow reusable coordinates
bool glChartCanvas::HasNormalizedViewPort(const ViewPort &vp) {
  return false;
#ifndef USE_ANDROID_GLES2
  return vp.m_projection_type == PROJECTION_MERCATOR ||
         vp.m_projection_type == PROJECTION_POLAR ||
         vp.m_projection_type == PROJECTION_EQUIRECTANGULAR;
#else
  return false;
#endif
}
 
/* adjust the opengl transformation matrix so that
   points plotted using the identity viewport are correct.
   and all rotation translation and scaling is now done in opengl
 
   a central lat and lon of 0, 0 can be used, however objects on the far side of
   the world can be up to 3 meters off because limited floating point precision,
   and if the points cross 180 longitude then two passes will be required to
   render them correctly */
#define NORM_FACTOR 4096.0
void glChartCanvas::MultMatrixViewPort(ViewPort &vp, float lat, float lon) {
#ifndef USE_ANDROID_GLES2
 
  wxPoint2DDouble point;
 
  switch (vp.m_projection_type) {
    case PROJECTION_MERCATOR:
    case PROJECTION_EQUIRECTANGULAR:
    case PROJECTION_WEB_MERCATOR:
      // m_pParentCanvas->GetDoubleCanvasPointPixVP(vp, lat, lon, &point);
      point = vp.GetDoublePixFromLL(lat, lon);
      glTranslated(point.m_x, point.m_y, 0);
      glScaled(vp.view_scale_ppm / NORM_FACTOR, vp.view_scale_ppm / NORM_FACTOR,
               1);
      break;
 
    case PROJECTION_POLAR:
      // m_pParentCanvas->GetDoubleCanvasPointPixVP(vp, vp.clat > 0 ? 90 : -90,
      // vp.clon, &point);
      point = vp.GetDoublePixFromLL(vp.clat > 0 ? 90 : -90, vp.clon);
      glTranslated(point.m_x, point.m_y, 0);
      glRotatef(vp.clon - lon, 0, 0, vp.clat);
      glScalef(vp.view_scale_ppm / NORM_FACTOR, vp.view_scale_ppm / NORM_FACTOR,
               1);
      glTranslatef(-vp.pix_width / 2, -vp.pix_height / 2, 0);
      break;
 
    default:
      printf("ERROR: Unhandled projection\n");
  }
 
  double rotation = vp.rotation;
 
  if (rotation) glRotatef(rotation * 180 / PI, 0, 0, 1);
#endif
}
 
ViewPort glChartCanvas::NormalizedViewPort(const ViewPort &vp, float lat,
                                           float lon) {
  ViewPort cvp = vp;
 
  switch (vp.m_projection_type) {
    case PROJECTION_MERCATOR:
    case PROJECTION_EQUIRECTANGULAR:
    case PROJECTION_WEB_MERCATOR:
      cvp.clat = lat;
      break;
 
    case PROJECTION_POLAR:
      cvp.clat = vp.clat > 0 ? 90 : -90;  // either north or south polar
      break;
 
    default:
      printf("ERROR: Unhandled projection\n");
  }
 
  cvp.clon = lon;
  cvp.view_scale_ppm = NORM_FACTOR;
  cvp.rotation = cvp.skew = 0;
  return cvp;
}
 
bool glChartCanvas::CanClipViewport(const ViewPort &vp) {
  return vp.m_projection_type == PROJECTION_MERCATOR ||
         vp.m_projection_type == PROJECTION_WEB_MERCATOR ||
         vp.m_projection_type == PROJECTION_EQUIRECTANGULAR;
}
 
ViewPort glChartCanvas::ClippedViewport(const ViewPort &vp,
                                        const LLRegion &region) {
  if (!CanClipViewport(vp)) return vp;
 
  ViewPort cvp = vp;
  LLBBox bbox = region.GetBox();
 
  if (!bbox.GetValid()) return vp;
 
  /* region.GetBox() will always try to give coordinates from -180 to 180 but in
     the case where the viewport crosses the IDL, we actually want the clipped
     viewport to use coordinates outside this range to ensure the logic in the
     various rendering routines works the same here (with accelerated panning)
     as it does without, so we can adjust the coordinates here */
 
  if (bbox.GetMaxLon() < cvp.GetBBox().GetMinLon()) {
    bbox.Set(bbox.GetMinLat(), bbox.GetMinLon() + 360, bbox.GetMaxLat(),
             bbox.GetMaxLon() + 360);
    cvp.SetBBoxDirect(bbox);
  } else if (bbox.GetMinLon() > cvp.GetBBox().GetMaxLon()) {
    bbox.Set(bbox.GetMinLat(), bbox.GetMinLon() - 360, bbox.GetMaxLat(),
             bbox.GetMaxLon() - 360);
    cvp.SetBBoxDirect(bbox);
  } else
    cvp.SetBBoxDirect(bbox);
 
  return cvp;
}
 
void glChartCanvas::DrawStaticRoutesTracksAndWaypoints(ViewPort &vp) {
  if (!m_pParentCanvas->m_bShowNavobjects) return;
  ocpnDC dc(*this);
 
  for (Track *pTrackDraw : g_TrackList) {
    /* defer rendering active tracks until later */
    ActiveTrack *pActiveTrack = dynamic_cast<ActiveTrack *>(pTrackDraw);
    if (pActiveTrack && pActiveTrack->IsRunning()) continue;
 
    TrackGui(*pTrackDraw).Draw(m_pParentCanvas, dc, vp, vp.GetBBox());
  }
 
  for (Route *pRouteDraw : *pRouteList) {
    if (!pRouteDraw) continue;
 
    /* defer rendering active routes until later */
    if (pRouteDraw->IsActive() || pRouteDraw->IsSelected()) continue;
 
    /* defer rendering routes being edited until later */
    if (pRouteDraw->m_bIsBeingEdited) continue;
 
    RouteGui(*pRouteDraw).DrawGL(vp, m_pParentCanvas, dc);
    //    pRouteDraw->DrawGL(vp, m_pParentCanvas, dc);
  }
 
  /* Waypoints not drawn as part of routes, and not being edited */
  if (vp.GetBBox().GetValid() && pWayPointMan) {
    for (RoutePoint *pWP : *pWayPointMan->GetWaypointList()) {
      if (pWP && (!pWP->m_bRPIsBeingEdited) && (!pWP->m_bIsInRoute))
        if (vp.GetBBox().ContainsMarge(pWP->m_lat, pWP->m_lon, .5))
          RoutePointGui(*pWP).DrawGL(vp, m_pParentCanvas, dc);
    }
  }
}
 
void glChartCanvas::DrawDynamicRoutesTracksAndWaypoints(ViewPort &vp) {
  ocpnDC dc(*this);
 
  for (Track *pTrackDraw : g_TrackList) {
    ActiveTrack *pActiveTrack = dynamic_cast<ActiveTrack *>(pTrackDraw);
    if (pActiveTrack && pActiveTrack->IsRunning())
      TrackGui(*pTrackDraw).Draw(m_pParentCanvas, dc, vp, vp.GetBBox());
    // We need Track::Draw() to dynamically render last (ownship) point.
  }
 
  for (Route *pRouteDraw : *pRouteList) {
    int drawit = 0;
    if (!pRouteDraw) continue;
 
    /* Active routes */
    if (pRouteDraw->IsActive() || pRouteDraw->IsSelected()) drawit++;
 
    /* Routes being edited */
    if (pRouteDraw->m_bIsBeingEdited) drawit++;
 
    /* Routes Selected */
    if (pRouteDraw->IsSelected()) drawit++;
 
    if (drawit) {
      const LLBBox &vp_box = vp.GetBBox(), &test_box = pRouteDraw->GetBBox();
      if (!vp_box.IntersectOut(test_box))
        RouteGui(*pRouteDraw).DrawGL(vp, m_pParentCanvas, dc);
      //        pRouteDraw->DrawGL(vp, m_pParentCanvas, dc);
    }
  }
 
  /* Waypoints not drawn as part of routes, which are being edited right now */
  if (vp.GetBBox().GetValid() && pWayPointMan) {
    for (RoutePoint *pWP : *pWayPointMan->GetWaypointList()) {
      if (pWP && pWP->m_bRPIsBeingEdited && !pWP->m_bIsInRoute)
        RoutePointGui(*pWP).DrawGL(vp, m_pParentCanvas, dc);
      //        pWP->DrawGL(vp, m_pParentCanvas, dc);
    }
  }
}
 
static void GetLatLonCurveDist(const ViewPort &vp, float &lat_dist,
                               float &lon_dist) {
  // This really could use some more thought, and possibly split at different
  // intervals based on chart skew and other parameters to optimize performance
  switch (vp.m_projection_type) {
    case PROJECTION_TRANSVERSE_MERCATOR:
      lat_dist = 4, lon_dist = 1;
      break;
    case PROJECTION_POLYCONIC:
      lat_dist = 2, lon_dist = 1;
      break;
    case PROJECTION_ORTHOGRAPHIC:
      lat_dist = 2, lon_dist = 2;
      break;
    case PROJECTION_POLAR:
      lat_dist = 180, lon_dist = 1;
      break;
    case PROJECTION_STEREOGRAPHIC:
    case PROJECTION_GNOMONIC:
      lat_dist = 2, lon_dist = 1;
      break;
    case PROJECTION_EQUIRECTANGULAR:
      // this is suboptimal because we don't care unless there is
      // a change in both lat AND lon (skewed chart)
      lat_dist = 2, lon_dist = 360;
      break;
    default:
      lat_dist = 180, lon_dist = 360;
  }
}
 
void glChartCanvas::RenderChartOutline(ocpnDC &dc, int dbIndex, ViewPort &vp) {
  if (ChartData->GetDBChartType(dbIndex) == CHART_TYPE_PLUGIN &&
      !ChartData->IsChartAvailable(dbIndex))
    return;
 
  /* quick bounds check */
  LLBBox box;
  ChartData->GetDBBoundingBox(dbIndex, box);
  if (!box.GetValid()) return;
 
  // Don't draw an outline in the case where the chart covers the entire world
  // */
  if (box.GetLonRange() == 360) return;
 
  LLBBox vpbox = vp.GetBBox();
 
  double lon_bias = 0;
  // chart is outside of viewport lat/lon bounding box
  if (box.IntersectOutGetBias(vp.GetBBox(), lon_bias)) return;
 
  wxColour color;
  if (ChartData->GetDBChartType(dbIndex) == CHART_TYPE_CM93)
    color = GetGlobalColor("YELO1");
  else if (ChartData->GetDBChartFamily(dbIndex) == CHART_FAMILY_VECTOR)
    color = GetGlobalColor("GREEN2");
  else
    color = GetGlobalColor("UINFR");
 
#if !defined(USE_ANDROID_GLES2) && !defined(ocpnUSE_GLSL)
  float plylat, plylon;
 
  if (g_GLOptions.m_GLLineSmoothing) glEnable(GL_LINE_SMOOTH);
 
  glColor3ub(color.Red(), color.Green(), color.Blue());
  glLineWidth(g_GLMinSymbolLineWidth);
 
  float lat_dist, lon_dist;
  GetLatLonCurveDist(vp, lat_dist, lon_dist);
 
  //        Are there any aux ply entries?
  int nAuxPlyEntries = ChartData->GetnAuxPlyEntries(dbIndex), nPly;
  int j = 0;
  do {
    if (nAuxPlyEntries)
      nPly = ChartData->GetDBAuxPlyPoint(dbIndex, 0, j, 0, 0);
    else
      nPly = ChartData->GetDBPlyPoint(dbIndex, 0, &plylat, &plylon);
 
    bool begin = false, sml_valid = false;
    double sml[2];
    float lastplylat = 0.0;
    float lastplylon = 0.0;
    // modulo is undefined for zero (compiler can use a div operation)
    int modulo = (nPly == 0) ? 1 : nPly;
    for (int i = 0; i < nPly + 1; i++) {
      if (nAuxPlyEntries)
        ChartData->GetDBAuxPlyPoint(dbIndex, i % modulo, j, &plylat, &plylon);
      else
        ChartData->GetDBPlyPoint(dbIndex, i % modulo, &plylat, &plylon);
 
      plylon += lon_bias;
 
      if (lastplylon - plylon > 180)
        lastplylon -= 360;
      else if (lastplylon - plylon < -180)
        lastplylon += 360;
 
      int splits;
      if (i == 0)
        splits = 1;
      else {
        int lat_splits = floor(fabs(plylat - lastplylat) / lat_dist);
        int lon_splits = floor(fabs(plylon - lastplylon) / lon_dist);
        splits = wxMax(lat_splits, lon_splits) + 1;
      }
 
      double smj[2];
      if (splits != 1) {
        // must perform border interpolation in mercator space as this is what
        // the charts use
        toSM(plylat, plylon, 0, 0, smj + 0, smj + 1);
        if (!sml_valid) toSM(lastplylat, lastplylon, 0, 0, sml + 0, sml + 1);
      }
 
      for (double c = 0; c < splits; c++) {
        double lat, lon;
        if (c == splits - 1)
          lat = plylat, lon = plylon;
        else {
          double d = (double)(c + 1) / splits;
          fromSM(d * smj[0] + (1 - d) * sml[0], d * smj[1] + (1 - d) * sml[1],
                 0, 0, &lat, &lon);
        }
 
        wxPoint2DDouble s;
        m_pParentCanvas->GetDoubleCanvasPointPix(lat, lon, &s);
        if (!std::isnan(s.m_x)) {
          if (!begin) {
            begin = true;
            glBegin(GL_LINE_STRIP);
          }
          glVertex2f(s.m_x, s.m_y);
        } else if (begin) {
          glEnd();
          begin = false;
        }
      }
      if ((sml_valid = splits != 1)) memcpy(sml, smj, sizeof smj);
      lastplylat = plylat, lastplylon = plylon;
    }
 
    if (begin) glEnd();
 
  } while (++j < nAuxPlyEntries);  // There are no aux Ply Point entries
 
  glDisable(GL_LINE_SMOOTH);
  //    glDisable( GL_BLEND );
 
#else
  double nominal_line_width_pix =
      wxMax(2.0, floor(m_pParentCanvas->GetPixPerMM() / 4));
 
  if (ChartData->GetDBChartType(dbIndex) == CHART_TYPE_CM93)
    dc.SetPen(wxPen(GetGlobalColor("YELO1"), nominal_line_width_pix,
                    wxPENSTYLE_SOLID));
 
  else if (ChartData->GetDBChartFamily(dbIndex) == CHART_FAMILY_VECTOR)
    dc.SetPen(wxPen(GetGlobalColor("UINFG"), nominal_line_width_pix,
                    wxPENSTYLE_SOLID));
 
  else
    dc.SetPen(wxPen(GetGlobalColor("UINFR"), nominal_line_width_pix,
                    wxPENSTYLE_SOLID));
 
  float plylat1, plylon1;
  int pixx1, pixy1;
 
  //        Are there any aux ply entries?
  int nAuxPlyEntries = ChartData->GetnAuxPlyEntries(dbIndex);
  if (0 == nAuxPlyEntries)  // There are no aux Ply Point entries
  {
    wxPoint r, r1;
    std::vector<int> points_vector;
 
    std::vector<float> vec = ChartData->GetReducedPlyPoints(dbIndex);
    int nPly = vec.size() / 2;
 
    if (nPly == 0) return;
 
    for (int i = 0; i < nPly; i++) {
      plylon1 = vec[i * 2];
      plylat1 = vec[i * 2 + 1];
 
      m_pParentCanvas->GetCanvasPointPix(plylat1, plylon1, &r1);
      pixx1 = r1.x;
      pixy1 = r1.y;
 
      points_vector.push_back(pixx1);
      points_vector.push_back(pixy1);
    }
 
    ChartData->GetDBPlyPoint(dbIndex, 0, &plylat1, &plylon1);
    plylon1 += lon_bias;
 
    m_pParentCanvas->GetCanvasPointPix(vec[1], vec[0], &r1);
    pixx1 = r1.x;
    pixy1 = r1.y;
 
    points_vector.push_back(pixx1);
    points_vector.push_back(pixy1);
 
    if (points_vector.size()) {
      std::vector<int>::iterator it = points_vector.begin();
      dc.DrawLines(points_vector.size() / 2, (wxPoint *)&(*it), 0, 0, true);
    }
  }
 
  else  // Use Aux PlyPoints
  {
    wxPoint r, r1;
 
    for (int j = 0; j < nAuxPlyEntries; j++) {
      std::vector<int> points_vector;
 
      std::vector<float> vec = ChartData->GetReducedAuxPlyPoints(dbIndex, j);
      int nAuxPly = vec.size() / 2;
 
      if (nAuxPly == 0) continue;
 
      for (int i = 0; i < nAuxPly; i++) {
        plylon1 = vec[i * 2];
        plylat1 = vec[i * 2 + 1];
 
        m_pParentCanvas->GetCanvasPointPix(plylat1, plylon1, &r1);
        pixx1 = r1.x;
        pixy1 = r1.y;
 
        points_vector.push_back(pixx1);
        points_vector.push_back(pixy1);
      }
 
      m_pParentCanvas->GetCanvasPointPix(vec[1], vec[0], &r1);
      pixx1 = r1.x;
      pixy1 = r1.y;
 
      points_vector.push_back(pixx1);
      points_vector.push_back(pixy1);
 
      if (points_vector.size()) {
        std::vector<int>::iterator it = points_vector.begin();
        dc.DrawLines(points_vector.size() / 2, (wxPoint *)&(*it), 0, 0, true);
      }
    }
  }
 
#endif
}
 
extern void CalcGridSpacing(float WindowDegrees, float &MajorSpacing,
                            float &MinorSpacing);
extern wxString CalcGridText(float latlon, float spacing, bool bPostfix);
void glChartCanvas::GridDraw() {
  if (!m_pParentCanvas->m_bDisplayGrid) return;
 
  ViewPort &vp = m_pParentCanvas->GetVP();
 
  if (!vp.IsValid() || !vp.GetBBox().GetValid()) return;
 
  // TODO: make minor grid work all the time
  bool minorgrid =
      fabs(vp.rotation) < 0.0001 && vp.m_projection_type == PROJECTION_MERCATOR;
 
  double nlat, elon, slat, wlon;
  float lat, lon;
  float gridlatMajor, gridlatMinor, gridlonMajor, gridlonMinor;
  wxCoord w, h;
 
  wxColour GridColor = GetGlobalColor("SNDG1");
 
  if (!m_gridfont.IsBuilt()) {
    double dpi_factor = g_BasePlatform->GetDisplayDIPMult(this);
    wxFont *dFont = FontMgr::Get().GetFont(_("GridText"), 0);
    wxFont font = *dFont;
    // Keep point size unscaled here; TexFont::Build handles DPI and display
    // scale.
    int font_size = wxMax(10, dFont->GetPointSize());
    font.SetPointSize(font_size);
    font.SetWeight(wxFONTWEIGHT_NORMAL);
 
    m_gridfont.SetContentScaleFactor(OCPN_GetDisplayContentScaleFactor());
    m_gridfont.Build(font, m_displayScale, dpi_factor);
  }
  m_gridfont.SetColor(GridColor);
 
  w = vp.pix_width;
  h = vp.pix_height;
 
  LLBBox llbbox = vp.GetBBox();
  nlat = llbbox.GetMaxLat();
  slat = llbbox.GetMinLat();
  elon = llbbox.GetMaxLon();
  wlon = llbbox.GetMinLon();
 
  // calculate distance between latitude grid lines
  CalcGridSpacing(vp.view_scale_ppm, gridlatMajor, gridlatMinor);
  CalcGridSpacing(vp.view_scale_ppm, gridlonMajor, gridlonMinor);
 
  // if it is known the grid has straight lines it's a bit faster
  bool straight_latitudes = vp.m_projection_type == PROJECTION_MERCATOR ||
                            vp.m_projection_type == PROJECTION_WEB_MERCATOR ||
                            vp.m_projection_type == PROJECTION_EQUIRECTANGULAR;
  bool straight_longitudes = vp.m_projection_type == PROJECTION_MERCATOR ||
                             vp.m_projection_type == PROJECTION_WEB_MERCATOR ||
                             vp.m_projection_type == PROJECTION_POLAR ||
                             vp.m_projection_type == PROJECTION_EQUIRECTANGULAR;
 
  double latmargin;
  if (straight_latitudes)
    latmargin = 0;
  else
    latmargin = gridlatMajor / 2;  // don't draw on poles
 
  slat = wxMax(slat, -90 + latmargin);
  nlat = wxMin(nlat, 90 - latmargin);
 
  float startlat = ceil(slat / gridlatMajor) * gridlatMajor;
  float startlon = ceil(wlon / gridlonMajor) * gridlonMajor;
  float curved_step = wxMin(sqrt(5e-3 / vp.view_scale_ppm), 3);
 
  ocpnDC gldc(*this);
  wxPen *pen = wxThePenList->FindOrCreatePen(GridColor, g_GLMinSymbolLineWidth,
                                             wxPENSTYLE_SOLID);
  gldc.SetPen(*pen);
 
  // Draw Major latitude grid lines and text
 
  // calculate position of first major latitude grid line
  float lon_step = elon - wlon;
  if (!straight_latitudes) lon_step /= ceil(lon_step / curved_step);
 
  for (lat = startlat; lat < nlat; lat += gridlatMajor) {
    wxPoint r, s;
    s.x = INVALID_COORD;
    s.y = INVALID_COORD;
    for (lon = wlon; lon < elon + lon_step / 2; lon += lon_step) {
      m_pParentCanvas->GetCanvasPointPix(lat, lon, &r);
      if (s.x != INVALID_COORD && s.y != INVALID_COORD) {
        gldc.DrawLine(s.x, s.y, r.x, r.y, false);
      }
      s = r;
    }
  }
 
  if (minorgrid) {
    // draw minor latitude grid lines
    for (lat = ceil(slat / gridlatMinor) * gridlatMinor; lat < nlat;
         lat += gridlatMinor) {
      wxPoint r;
      m_pParentCanvas->GetCanvasPointPix(lat, (elon + wlon) / 2, &r);
      gldc.DrawLine(0, r.y, 10, r.y, true);
      gldc.DrawLine(w - 10, r.y, w, r.y, false);
    }
  }
 
  // draw major longitude grid lines
  float lat_step = nlat - slat;
  if (!straight_longitudes) lat_step /= ceil(lat_step / curved_step);
 
  for (lon = startlon; lon < elon; lon += gridlonMajor) {
    wxPoint r, s;
    s.x = INVALID_COORD;
    s.y = INVALID_COORD;
    for (lat = slat; lat < nlat + lat_step / 2; lat += lat_step) {
      m_pParentCanvas->GetCanvasPointPix(lat, lon, &r);
      if (s.x != INVALID_COORD && s.y != INVALID_COORD) {
        gldc.DrawLine(s.x, s.y, r.x, r.y, false);
      }
      s = r;
    }
  }
 
  if (minorgrid) {
    // draw minor longitude grid lines
    for (lon = ceil(wlon / gridlonMinor) * gridlonMinor; lon < elon;
         lon += gridlonMinor) {
      wxPoint r;
      m_pParentCanvas->GetCanvasPointPix((nlat + slat) / 2, lon, &r);
      gldc.DrawLine(r.x, 0, r.x, 10, false);
      gldc.DrawLine(r.x, h - 10, r.x, h, false);
    }
  }
 
  // draw text labels
  if (abs(vp.rotation) < .1) {
    glEnable(GL_TEXTURE_2D);
    glEnable(GL_BLEND);
    for (lat = startlat; lat < nlat; lat += gridlatMajor) {
      if (fabs(lat - wxRound(lat)) < 1e-5) lat = wxRound(lat);
 
      wxString st =
          CalcGridText(lat, gridlatMajor, true);  // get text for grid line
      int iy;
      m_gridfont.GetTextExtent(st, 0, &iy);
 
      if (straight_latitudes) {
        wxPoint r, s;
        m_pParentCanvas->GetCanvasPointPix(lat, elon, &r);
        m_pParentCanvas->GetCanvasPointPix(lat, wlon, &s);
 
        float x = 0, y = -1;
        y = (float)(r.y * s.x - s.y * r.x) / (s.x - r.x);
        if (y < 0 || y > h) {
          y = h - iy;
          x = (float)(r.x * s.y - s.x * r.y + (s.x - r.x) * y) / (s.y - r.y);
        }
 
        m_gridfont.RenderString(st, x, y);
      } else {
        // iteratively attempt to find where the latitude line crosses x=0
        wxPoint2DDouble r;
        double y1, y2, lat1, lon1, lat2, lon2;
 
        y1 = 0, y2 = vp.pix_height;
        double error = vp.pix_width, lasterror;
        int maxiters = 10;
        do {
          m_pParentCanvas->GetCanvasPixPoint(0, y1, lat1, lon1);
          m_pParentCanvas->GetCanvasPixPoint(0, y2, lat2, lon2);
 
          double y = y1 + (lat1 - lat) * (y2 - y1) / (lat1 - lat2);
 
          m_pParentCanvas->GetDoubleCanvasPointPix(
              lat, lon1 + (y1 - y) * (lon2 - lon1) / (y1 - y2), &r);
 
          if (fabs(y - y1) < fabs(y - y2))
            y1 = y;
          else
            y2 = y;
 
          lasterror = error;
          error = fabs(r.m_x);
          if (--maxiters == 0) break;
        } while (error > 1 && error < lasterror);
 
        if (error < 1 && r.m_y >= 0 && r.m_y <= vp.pix_height - iy)
          r.m_x = 0;
        else
          // just draw at center longitude
          m_pParentCanvas->GetDoubleCanvasPointPix(lat, vp.clon, &r);
 
        m_gridfont.RenderString(st, r.m_x, r.m_y);
      }
    }
 
    for (lon = startlon; lon < elon; lon += gridlonMajor) {
      if (fabs(lon - wxRound(lon)) < 1e-5) lon = wxRound(lon);
 
      wxPoint r, s;
      m_pParentCanvas->GetCanvasPointPix(nlat, lon, &r);
      m_pParentCanvas->GetCanvasPointPix(slat, lon, &s);
 
      float xlon = lon;
      if (xlon > 180.0)
        xlon -= 360.0;
      else if (xlon <= -180.0)
        xlon += 360.0;
 
      wxString st = CalcGridText(xlon, gridlonMajor, false);
      int ix;
      m_gridfont.GetTextExtent(st, &ix, 0);
 
      if (straight_longitudes) {
        float x = -1, y = 0;
        x = (float)(r.x * s.y - s.x * r.y) / (s.y - r.y);
        if (x < 0 || x > w) {
          x = w - ix;
          y = (float)(r.y * s.x - s.y * r.x + (s.y - r.y) * x) / (s.x - r.x);
        }
 
        m_gridfont.RenderString(st, x, y);
      } else {
        // iteratively attempt to find where the latitude line crosses x=0
        wxPoint2DDouble r;
        double x1, x2, lat1, lon1, lat2, lon2;
 
        x1 = 0, x2 = vp.pix_width;
        double error = vp.pix_height, lasterror;
        do {
          m_pParentCanvas->GetCanvasPixPoint(x1, 0, lat1, lon1);
          m_pParentCanvas->GetCanvasPixPoint(x2, 0, lat2, lon2);
 
          double x = x1 + (lon1 - lon) * (x2 - x1) / (lon1 - lon2);
 
          m_pParentCanvas->GetDoubleCanvasPointPix(
              lat1 + (x1 - x) * (lat2 - lat1) / (x1 - x2), lon, &r);
 
          if (fabs(x - x1) < fabs(x - x2))
            x1 = x;
          else
            x2 = x;
 
          lasterror = error;
          error = fabs(r.m_y);
        } while (error > 1 && error < lasterror);
 
        if (error < 1 && r.m_x >= 0 && r.m_x <= vp.pix_width - ix)
          r.m_y = 0;
        else
          // failure, instead just draw the text at center latitude
          m_pParentCanvas->GetDoubleCanvasPointPix(
              wxMin(wxMax(vp.clat, slat), nlat), lon, &r);
 
        m_gridfont.RenderString(st, r.m_x, r.m_y);
      }
    }
 
    glDisable(GL_TEXTURE_2D);
    glDisable(GL_BLEND);
  }
}
 
void glChartCanvas::DrawEmboss(ocpnDC &dc, emboss_data *emboss) {
  if (!emboss) return;
 
  int w = emboss->width, h = emboss->height;
 
  glEnable(GL_TEXTURE_2D);
 
  // render using opengl and alpha blending
  if (!emboss->gltexind) { /* upload to texture */
 
    emboss->glwidth = NextPow2(emboss->width);
    emboss->glheight = NextPow2(emboss->height);
 
    /* convert to luminance alpha map */
    int size = emboss->glwidth * emboss->glheight;
    char *data = new char[2 * size];
    for (int i = 0; i < h; i++) {
      for (int j = 0; j < emboss->glwidth; j++) {
        if (j < w) {
          data[2 * ((i * emboss->glwidth) + j)] =
              (char)(emboss->pmap[(i * w) + j] > 0 ? 0 : 255);
          data[2 * ((i * emboss->glwidth) + j) + 1] =
              (char)abs((emboss->pmap[(i * w) + j]));
        }
      }
    }
 
    glGenTextures(1, &emboss->gltexind);
    glBindTexture(GL_TEXTURE_2D, emboss->gltexind);
    glTexImage2D(GL_TEXTURE_2D, 0, GL_LUMINANCE_ALPHA, emboss->glwidth,
                 emboss->glheight, 0, GL_LUMINANCE_ALPHA, GL_UNSIGNED_BYTE,
                 data);
    glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);
    glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);
 
    delete[] data;
  }
 
  glBindTexture(GL_TEXTURE_2D, emboss->gltexind);
 
  glEnable(GL_BLEND);
 
  int x = emboss->x, y = emboss->y;
 
  float wp = (float)w / emboss->glwidth;
  float hp = (float)h / emboss->glheight;
 
  float coords[8];
  float uv[8];
 
  // normal uv
  uv[0] = 0;
  uv[1] = 0;
  uv[2] = wp;
  uv[3] = 0;
  uv[4] = wp;
  uv[5] = hp;
  uv[6] = 0;
  uv[7] = hp;
 
  // pixels
  coords[0] = 0;
  coords[1] = 0;
  coords[2] = w;
  coords[3] = 0;
  coords[4] = w;
  coords[5] = h;
  coords[6] = 0;
  coords[7] = h;
 
  // FIXME(dave) Find a way to make this thing a little transparaent
  RenderSingleTexture(dc, coords, uv, m_pParentCanvas->GetpVP(), x, y, 0);
 
  glDisable(GL_BLEND);
  glDisable(GL_TEXTURE_2D);
}
 
void glChartCanvas::ShipDraw(ocpnDC &dc) {
  if (!m_pParentCanvas->GetVP().IsValid()) return;
  wxPoint GPSOffsetPixels(0, 0);
  wxPoint2DDouble lGPSPoint, lShipMidPoint;
 
  //  COG/SOG may be undefined in NMEA data stream
  float pCog = std::isnan(gCog) ? 0 : gCog;
  float pSog = std::isnan(gSog) ? 0 : gSog;
 
  // Here, calculate the ownship location on screen, and make it so.
  // Special case:  No need for such precision on chart drag operations
  double shift_dx = 0;
  double shift_dy = 0;
  bool dynamic = m_pParentCanvas->m_animationActive ||
                 m_pParentCanvas->m_MouseDragging ||
                 m_pParentCanvas->m_chart_drag_inertia_active;
  if (m_pParentCanvas->m_bFollow && !dynamic) {
    lGPSPoint.m_x = m_pParentCanvas->GetVP().pix_width / 2;
    lGPSPoint.m_y = m_pParentCanvas->GetVP().pix_height / 2;
    if (m_pParentCanvas->m_bLookAhead) {
      // In "b_follow" mode, we have a-priori information about the desired
      // screen coordinates of ownship, as a pixel offset from center.
      // Use this information as a performance optimization,
      // at larger viewing scales.
      if (m_pParentCanvas->GetVPChartScale() < 2e6) {
        double angle = m_pParentCanvas->dir_to_shift * PI / 180.;
        angle += m_pParentCanvas->GetVPRotation();
        shift_dx = m_pParentCanvas->meters_to_shift * sin(angle) *
                   m_pParentCanvas->GetVPScale();
        lGPSPoint.m_x -= shift_dx / cos(gLat * PI / 180.);
        shift_dy = m_pParentCanvas->meters_to_shift * cos(angle) *
                   m_pParentCanvas->GetVPScale();
        lGPSPoint.m_y += shift_dy / cos(gLat * PI / 180.);
      } else {
        m_pParentCanvas->GetDoubleCanvasPointPix(gLat, gLon, &lGPSPoint);
      }
    } else {
      lGPSPoint.m_x -= m_pParentCanvas->m_OSoffsetx;
      lGPSPoint.m_y += m_pParentCanvas->m_OSoffsety;
    }
  } else {
    m_pParentCanvas->GetDoubleCanvasPointPixVP(m_pParentCanvas->GetVP(), gLat,
                                               gLon, &lGPSPoint);
  }
 
  lShipMidPoint = lGPSPoint;
 
  //  Draw the icon rotated to the COG
  //  or to the Hdt if available
  float icon_hdt = pCog;
  if (!std::isnan(gHdt)) icon_hdt = gHdt;
 
  //  COG may be undefined in NMEA data stream
  if (std::isnan(icon_hdt)) icon_hdt = 0.0;
 
  //    Calculate the ownship drawing angle icon_rad using an assumed 10 minute
  //    predictor
  double osd_head_lat, osd_head_lon;
  wxPoint2DDouble osd_head_point;
 
  ll_gc_ll(gLat, gLon, icon_hdt, pSog * 10. / 60., &osd_head_lat,
           &osd_head_lon);
 
  m_pParentCanvas->GetDoubleCanvasPointPixVP(
      m_pParentCanvas->GetVP(), osd_head_lat, osd_head_lon, &osd_head_point);
 
  double icon_rad = atan2f((float)(osd_head_point.m_y - lShipMidPoint.m_y),
                           (float)(osd_head_point.m_x - lShipMidPoint.m_x));
  icon_rad += (float)PI;
 
  if (pSog < 0.2)
    icon_rad =
        ((icon_hdt + 90.) * PI / 180.) + m_pParentCanvas->GetVP().rotation;
 
  //    Another draw test ,based on pixels, assuming the ship icon is a fixed
  //    nominal size and is just barely outside the viewport        ....
  BoundingBox bb_screen(0, 0, m_pParentCanvas->GetVP().pix_width,
                        m_pParentCanvas->GetVP().pix_height);
 
  // TODO: fix to include actual size of boat that will be rendered
  int img_height = 0;
 
  if (bb_screen.PointInBox(lShipMidPoint, 20)) {
    if (g_GLOptions.m_GLLineSmoothing) glEnable(GL_LINE_SMOOTH);
    if (g_GLOptions.m_GLPolygonSmoothing) glEnable(GL_POLYGON_SMOOTH);
 
    if (m_pParentCanvas->GetVP().chart_scale >
        300000)  // According to S52, this should be 50,000
    {
      float scale_factor = 1.0;
      // Scale the generic icon to ChartScaleFactor, slightly softened....
      if ((g_ChartScaleFactorExp > 1.0) && (g_OwnShipIconType == 0))
        scale_factor = (log(g_ChartScaleFactorExp) + 1.0) * 1.1;
 
      float nominal_ownship_size_mm = m_pParentCanvas->m_display_size_mm / 44.0;
      nominal_ownship_size_mm = wxMin(nominal_ownship_size_mm, 15.0);
      nominal_ownship_size_mm = wxMax(nominal_ownship_size_mm, 7.0);
 
      scale_factor *= m_pParentCanvas->GetContentScaleFactor();
 
      float nominal_ownship_size_pixels =
          wxMax(20.0, m_pParentCanvas->GetPixPerMM() *
                          nominal_ownship_size_mm);  // nominal length, but not
                                                     // less than 20 pixel
      float v = (nominal_ownship_size_pixels * scale_factor) / 3;
 
      wxPen ppSmallScaleShip;
      if (SHIP_NORMAL == m_pParentCanvas->m_ownship_state)
        ppSmallScaleShip =
            wxPen(GetGlobalColor("URED"), v / 5, wxPENSTYLE_SOLID);
      else
        ppSmallScaleShip =
            wxPen(GetGlobalColor("YELO1"), v / 5, wxPENSTYLE_SOLID);
      dc.SetPen(ppSmallScaleShip);
 
      dc.SetBrush(wxBrush(GetGlobalColor("URED"), wxBRUSHSTYLE_TRANSPARENT));
 
      // start with cross
      dc.DrawLine((-v * 1.2) + lShipMidPoint.m_x, lShipMidPoint.m_y,
                  (v * 1.2) + lShipMidPoint.m_x, lShipMidPoint.m_y);
      dc.DrawLine(lShipMidPoint.m_x, (-v * 1.2) + lShipMidPoint.m_y,
                  lShipMidPoint.m_x, (v * 1.2) + lShipMidPoint.m_y);
 
      //  Two circles
      dc.StrokeCircle(lShipMidPoint.m_x, lShipMidPoint.m_y, v);
      dc.StrokeCircle(lShipMidPoint.m_x, lShipMidPoint.m_y, 0.6 * v);
      img_height = 20;
    } else {
      int draw_color = SHIP_INVALID;
      if (SHIP_NORMAL == m_pParentCanvas->m_ownship_state)
        draw_color = SHIP_NORMAL;
      else if (SHIP_LOWACCURACY == m_pParentCanvas->m_ownship_state)
        draw_color = SHIP_LOWACCURACY;
 
      if (!ownship_tex ||
          (draw_color !=
           ownship_color)) { /* initial run, create texture for ownship,
                                also needed at colorscheme changes (not
                                implemented) */
 
        ownship_color = draw_color;
 
        if (ownship_tex) glDeleteTextures(1, &ownship_tex);
 
        glGenTextures(1, &ownship_tex);
        glBindTexture(GL_TEXTURE_2D, ownship_tex);
 
        glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);
        glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);
 
        wxImage image;
        if (m_pParentCanvas->m_pos_image_user) {
          switch (draw_color) {
            case SHIP_INVALID:
              image = *m_pParentCanvas->m_pos_image_user_grey;
              break;
            case SHIP_NORMAL:
              image = *m_pParentCanvas->m_pos_image_user;
              break;
            case SHIP_LOWACCURACY:
              image = *m_pParentCanvas->m_pos_image_user_yellow;
              break;
          }
        } else {
          switch (draw_color) {
            case SHIP_INVALID:
              image = *m_pParentCanvas->m_pos_image_grey;
              break;
            case SHIP_NORMAL:
              image = *m_pParentCanvas->m_pos_image_red;
              break;
            case SHIP_LOWACCURACY:
              image = *m_pParentCanvas->m_pos_image_yellow;
              break;
          }
        }
 
        int w = image.GetWidth(), h = image.GetHeight();
        int glw = NextPow2(w), glh = NextPow2(h);
        ownship_size = wxSize(w, h);
        ownship_tex_size = wxSize(glw, glh);
 
        unsigned char *d = image.GetData();
        unsigned char *a = image.GetAlpha();
        unsigned char *e = new unsigned char[4 * w * h];
 
        if (d && e && a) {
          for (int p = 0; p < w * h; p++) {
            e[4 * p + 0] = d[3 * p + 0];
            e[4 * p + 1] = d[3 * p + 1];
            e[4 * p + 2] = d[3 * p + 2];
            e[4 * p + 3] = a[p];
          }
        }
        glTexImage2D(GL_TEXTURE_2D, 0, GL_RGBA, glw, glh, 0, GL_RGBA,
                     GL_UNSIGNED_BYTE, 0);
 
        glTexSubImage2D(GL_TEXTURE_2D, 0, 0, 0, w, h, GL_RGBA, GL_UNSIGNED_BYTE,
                        e);
        delete[] e;
      }
 
      /* establish ship color */
#ifndef USE_ANDROID_GLES2
      if (m_pParentCanvas->m_pos_image_user)
        glColor4ub(255, 255, 255, 255);
      else if (SHIP_NORMAL == m_pParentCanvas->m_ownship_state)
        glColor4ub(255, 0, 0, 255);
      else if (SHIP_LOWACCURACY == m_pParentCanvas->m_ownship_state)
        glColor4ub(255, 255, 0, 255);
      else
        glColor4ub(128, 128, 128, 255);
#endif
      float scale_factor_y = 1.0;
      float scale_factor_x = 1.0;
 
      int ownShipWidth = 22;  // Default values from s_ownship_icon
      int ownShipLength = 84;
      lShipMidPoint = lGPSPoint;
 
      /* scaled ship? */
      if (g_OwnShipIconType != 0)
        m_pParentCanvas->ComputeShipScaleFactor(
            icon_hdt, ownShipWidth, ownShipLength, lShipMidPoint,
            GPSOffsetPixels, lGPSPoint, scale_factor_x, scale_factor_y);
 
      glEnable(GL_BLEND);
 
      // Scale the generic icon to ChartScaleFactor, slightly softened....
      if ((g_ShipScaleFactorExp > 1.0) && (g_OwnShipIconType == 0)) {
        scale_factor_x = (log(g_ShipScaleFactorExp) + 1.0) * 1.1;
        scale_factor_y = (log(g_ShipScaleFactorExp) + 1.0) * 1.1;
      }
 
      // Correct for scaled displays, e.g. Retina
      scale_factor_x *= m_pParentCanvas->GetContentScaleFactor();
      scale_factor_y *= m_pParentCanvas->GetContentScaleFactor();
 
      // Set the size of the little circle showing the GPS reference position
      // Set a default early, adjust later based on render type
      float gps_circle_radius = 3.0;
 
      if (g_OwnShipIconType == 0) {  // Default Bitmap
 
        glEnable(GL_TEXTURE_2D);
        glBindTexture(GL_TEXTURE_2D, ownship_tex);
 
        // We choose to render the ownship bitmap at roughly the same size( in
        // pixels ) as the DC mode renderer.
        // For ultra-high definition displays, we clamp the actual on-screen
        // size to be no smaller than 7.0 mm Similarly, for lo-res displays, we
        // limit the actual size to be no larger than 15 mm maximum.
 
        // Get bitmap height in pixels
        int image_height_bitmap = m_pParentCanvas->m_pos_image_red->GetHeight();
        if (m_pParentCanvas->m_pos_image_user)
          image_height_bitmap = m_pParentCanvas->m_pos_image_user->GetHeight();
 
        float nominal_ownship_size_mm =
            image_height_bitmap / m_pParentCanvas->GetPixPerMM();
 
        nominal_ownship_size_mm = wxMin(nominal_ownship_size_mm, 15.0);
        nominal_ownship_size_mm = wxMax(nominal_ownship_size_mm, 7.0);
 
        float nominal_ownship_size_pixels =
            m_pParentCanvas->GetPixPerMM() * nominal_ownship_size_mm;
 
        if (m_pParentCanvas->GetContentScaleFactor() == 1.0) {
          nominal_ownship_size_pixels = wxMax(
              20.0, nominal_ownship_size_pixels);  // not less than 20 pixel
        }
 
        float h = nominal_ownship_size_pixels * scale_factor_y;
        float w = nominal_ownship_size_pixels * scale_factor_x *
                  ownship_size.x / ownship_size.y;
        float glw = ownship_tex_size.x, glh = ownship_tex_size.y;
        float u = ownship_size.x / glw, v = ownship_size.y / glh;
 
        //                printf("%g %g %g %g %g %g %g\n",
        //                nominal_ownship_size_mm, nominal_ownship_size_pixels,
        //                h, w, u, v, m_pParentCanvas->m_display_size_mm);
        // tweak GPS reference point indicator size
        gps_circle_radius = w / 5;
 
        float uv[8], coords[8];
        uv[0] = 0;
        uv[1] = 0;
        uv[2] = u;
        uv[3] = 0;
        uv[4] = u;
        uv[5] = v;
        uv[6] = 0;
        uv[7] = v;
 
        coords[0] = -w / 2;
        coords[1] = -h / 2;
        coords[2] = w / 2;
        coords[3] = -h / 2;
        coords[4] = w / 2;
        coords[5] = h / 2;
        coords[6] = -w / 2;
        coords[7] = h / 2;
 
        RenderSingleTexture(dc, coords, uv, m_pParentCanvas->GetpVP(),
                            lShipMidPoint.m_x, lShipMidPoint.m_y,
                            icon_rad - PI / 2);
        glDisable(GL_TEXTURE_2D);
      } else if (g_OwnShipIconType == 1) {  // Scaled Bitmap
 
        glEnable(GL_TEXTURE_2D);
        glBindTexture(GL_TEXTURE_2D, ownship_tex);
 
        float nominal_ownship_size_pixels_y = 84;
        float nominal_ownship_size_pixels_x = 22;
 
        float h = nominal_ownship_size_pixels_y * scale_factor_y;
        float w = nominal_ownship_size_pixels_x * scale_factor_x;
 
        float u = (float)ownship_size.x / ownship_tex_size.x,
              v = (float)ownship_size.y / ownship_tex_size.y;
 
        // tweak GPS reference point indicator size
        gps_circle_radius = w / 5;
 
        float uv[8], coords[8];
        uv[0] = 0;
        uv[1] = 0;
        uv[2] = u;
        uv[3] = 0;
        uv[4] = u;
        uv[5] = v;
        uv[6] = 0;
        uv[7] = v;
 
        coords[0] = -w / 2;
        coords[1] = -h / 2;
        coords[2] = w / 2;
        coords[3] = -h / 2;
        coords[4] = w / 2;
        coords[5] = h / 2;
        coords[6] = -w / 2;
        coords[7] = h / 2;
 
        RenderSingleTexture(dc, coords, uv, m_pParentCanvas->GetpVP(),
                            lShipMidPoint.m_x, lShipMidPoint.m_y,
                            icon_rad - PI / 2);
 
        glDisable(GL_TEXTURE_2D);
      } else if (g_OwnShipIconType == 2) {  // Scaled Vector
        //                 static const GLint s_ownship_icon[] = { 5, -42, 11,
        //                 -28, 11, 42, -11, 42,
        //                                                         -11, -28, -5,
        //                                                         -42, -11, 0,
        //                                                         11, 0, 0, 42,
        //                                                         0, -42 };
 
        wxPoint shipPoints[6];
 
        wxColour colour = m_pParentCanvas->ShipColor();
        wxPen ppPen(*wxBLACK, 1);
        wxBrush ppBrush(colour);
        dc.SetPen(ppPen);
        dc.SetBrush(ppBrush);
 
        shipPoints[0].x = 0 * scale_factor_x;
        shipPoints[0].y = -28 * scale_factor_y;
        shipPoints[1].x = 11 * scale_factor_x;
        shipPoints[1].y = -28 * scale_factor_y;
        shipPoints[2].x = 11 * scale_factor_x;
        shipPoints[2].y = 42 * scale_factor_y;
        shipPoints[3].x = 0 * scale_factor_x;
        shipPoints[3].y = 42 * scale_factor_y;
        dc.DrawPolygon(4, shipPoints, lShipMidPoint.m_x, lShipMidPoint.m_y, 1,
                       icon_rad - PI / 2);
 
        shipPoints[0].x = 0 * scale_factor_x;
        shipPoints[0].y = -42 * scale_factor_y;
        shipPoints[1].x = 5 * scale_factor_x;
        shipPoints[1].y = -42 * scale_factor_y;
        shipPoints[2].x = 11 * scale_factor_x;
        shipPoints[2].y = -28 * scale_factor_y;
        shipPoints[3].x = 0 * scale_factor_x;
        shipPoints[3].y = -28 * scale_factor_y;
        dc.DrawPolygon(4, shipPoints, lShipMidPoint.m_x, lShipMidPoint.m_y, 1,
                       icon_rad - PI / 2);
 
        shipPoints[0].x = 0 * scale_factor_x;
        shipPoints[0].y = -28 * scale_factor_y;
        shipPoints[1].x = -11 * scale_factor_x;
        shipPoints[1].y = -28 * scale_factor_y;
        shipPoints[2].x = -11 * scale_factor_x;
        shipPoints[2].y = 42 * scale_factor_y;
        shipPoints[3].x = 0 * scale_factor_x;
        shipPoints[3].y = 42 * scale_factor_y;
        dc.DrawPolygon(4, shipPoints, lShipMidPoint.m_x, lShipMidPoint.m_y, 1,
                       icon_rad - PI / 2);
 
        shipPoints[0].x = 0 * scale_factor_x;
        shipPoints[0].y = -42 * scale_factor_y;
        shipPoints[1].x = -5 * scale_factor_x;
        shipPoints[1].y = -42 * scale_factor_y;
        shipPoints[2].x = -11 * scale_factor_x;
        shipPoints[2].y = -28 * scale_factor_y;
        shipPoints[3].x = 0 * scale_factor_x;
        shipPoints[3].y = -28 * scale_factor_y;
        dc.DrawPolygon(4, shipPoints, lShipMidPoint.m_x, lShipMidPoint.m_y, 1,
                       icon_rad - PI / 2);
 
        // draw with cross
        double p1x = -11 * scale_factor_x;
        double p2x = 11 * scale_factor_x;
        double p1y = 0;
        double p2y = 0;
        double p1xr =
            ((p1x)*cos(icon_rad - PI / 2)) - ((p1y)*sin(icon_rad - PI / 2));
        double p2xr =
            ((p2x)*cos(icon_rad - PI / 2)) - ((p2y)*sin(icon_rad - PI / 2));
        double p1yr =
            ((p1y)*cos(icon_rad - PI / 2)) + ((p1x)*sin(icon_rad - PI / 2));
        double p2yr =
            ((p2y)*cos(icon_rad - PI / 2)) + ((p2x)*sin(icon_rad - PI / 2));
        dc.DrawLine(p1xr + lShipMidPoint.m_x, p1yr + lShipMidPoint.m_y,
                    p2xr + lShipMidPoint.m_x, p2yr + lShipMidPoint.m_y);
 
        p1x = 0;
        p2x = 0;
        p1y = -42 * scale_factor_y;
        p2y = 42 * scale_factor_y;
        p1xr = ((p1x)*cos(icon_rad - PI / 2)) - ((p1y)*sin(icon_rad - PI / 2));
        p2xr = ((p2x)*cos(icon_rad - PI / 2)) - ((p2y)*sin(icon_rad - PI / 2));
        p1yr = ((p1y)*cos(icon_rad - PI / 2)) + ((p1x)*sin(icon_rad - PI / 2));
        p2yr = ((p2y)*cos(icon_rad - PI / 2)) + ((p2x)*sin(icon_rad - PI / 2));
        dc.DrawLine(p1xr + lShipMidPoint.m_x, p1yr + lShipMidPoint.m_y,
                    p2xr + lShipMidPoint.m_x, p2yr + lShipMidPoint.m_y);
      }
 
      img_height = ownShipLength * scale_factor_y;
 
      //      Reference point, where the GPS antenna is
      if (m_pParentCanvas->m_pos_image_user) gps_circle_radius = 1;
 
      wxPen ppPen1(GetGlobalColor("UBLCK"), 1, wxPENSTYLE_SOLID);
      dc.SetPen(ppPen1);
      dc.SetBrush(wxBrush(GetGlobalColor("CHWHT")));
 
      dc.StrokeCircle(lGPSPoint.m_x, lGPSPoint.m_y, gps_circle_radius);
    }
 
    //        glDisableClientState(GL_VERTEX_ARRAY);
    glDisable(GL_LINE_SMOOTH);
    glDisable(GL_POLYGON_SMOOTH);
    glDisable(GL_BLEND);
  }
 
  m_pParentCanvas->ShipIndicatorsDraw(dc, img_height, GPSOffsetPixels,
                                      lGPSPoint);
}
 
void glChartCanvas::DrawFloatingOverlayObjects(ocpnDC &dc) {
  ViewPort &vp = m_pParentCanvas->GetVP();
 
  //  Draw any active or selected routes now
  extern Routeman *g_pRouteMan;
  //    extern Track                     *g_pActiveTrack;
  Route *active_route = g_pRouteMan->GetpActiveRoute();
 
  //    if( active_route ) active_route->DrawGL( vp, region );
  //    if( g_pActiveTrack ) g_pActiveTrack->Draw( dc, vp );
  //    if( m_pParentCanvas->m_pSelectedRoute )
  //    m_pParentCanvas->m_pSelectedRoute->DrawGL( vp, region );
 
  GridDraw();
 
  g_overlayCanvas = m_pParentCanvas;
  if (g_pi_manager) {
    g_pi_manager->SendViewPortToRequestingPlugIns(vp);
    g_pi_manager->RenderAllGLCanvasOverlayPlugIns(
        m_pcontext, vp, m_pParentCanvas->m_canvasIndex, OVERLAY_LEGACY);
  }
 
  // all functions called with m_pParentCanvas-> are still slow because they go
  // through ocpndc
  AISDrawAreaNotices(dc, m_pParentCanvas->GetVP(), m_pParentCanvas);
 
  m_pParentCanvas->DrawAnchorWatchPoints(dc);
  AISDraw(dc, m_pParentCanvas->GetVP(), m_pParentCanvas);
  ShipDraw(dc);
  m_pParentCanvas->AlertDraw(dc);
 
  m_pParentCanvas->RenderVisibleSectorLights(dc);
 
  m_pParentCanvas->RenderRouteLegs(dc);
  m_pParentCanvas->RenderShipToActive(dc, true);
  m_pParentCanvas->ScaleBarDraw(dc);
  s57_DrawExtendedLightSectorsGL(dc, m_pParentCanvas->VPoint,
                                 m_pParentCanvas->extendedSectorLegs);
  if (g_pi_manager) {
    g_pi_manager->RenderAllGLCanvasOverlayPlugIns(
        m_pcontext, vp, m_pParentCanvas->m_canvasIndex, OVERLAY_OVER_SHIPS);
  }
}
 
void glChartCanvas::DrawChartBar(ocpnDC &dc) {
  if (m_pParentCanvas->GetPiano()) {
    int canvas_height = GetClientSize().y;
    canvas_height *= m_displayScale;
 
    m_pParentCanvas->GetPiano()->DrawGL(
        canvas_height - m_pParentCanvas->GetPiano()->GetHeight());
  }
}
 
void glChartCanvas::DrawQuiting() {
#ifndef USE_ANDROID_GLES2
  GLubyte pattern[8][8];
  for (int y = 0; y < 8; y++)
    for (int x = 0; x < 8; x++) pattern[y][x] = (y == x) * 255;
 
  glEnable(GL_BLEND);
  glEnable(GL_TEXTURE_2D);
  glBindTexture(GL_TEXTURE_2D, 0);
 
  glTexEnvi(GL_TEXTURE_ENV, GL_TEXTURE_ENV_MODE, GL_REPLACE);
  glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);
  glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);
 
  glTexImage2D(GL_TEXTURE_2D, 0, GL_ALPHA, 8, 8, 0, GL_ALPHA, GL_UNSIGNED_BYTE,
               pattern);
  glColor3f(0, 0, 0);
 
  float x = GetSize().x, y = GetSize().y;
  float u = x / 8, v = y / 8;
 
  glBegin(GL_QUADS);
  glTexCoord2f(0, 0);
  glVertex2f(0, 0);
  glTexCoord2f(0, v);
  glVertex2f(0, y);
  glTexCoord2f(u, v);
  glVertex2f(x, y);
  glTexCoord2f(u, 0);
  glVertex2f(x, 0);
  glEnd();
 
  glDisable(GL_TEXTURE_2D);
  glDisable(GL_BLEND);
#endif
}
 
void glChartCanvas::DrawCloseMessage(wxString msg) {
#ifndef USE_ANDROID_GLES2
 
  if (1) {
    wxFont *pfont = FontMgr::Get().FindOrCreateFont(
        12, wxFONTFAMILY_DEFAULT, wxFONTSTYLE_NORMAL, wxFONTWEIGHT_NORMAL);
 
    TexFont texfont;
 
    texfont.Build(*pfont, 1, 1);
    int w, h;
    texfont.GetTextExtent(msg, &w, &h);
    h += 2;
    int yp = m_pParentCanvas->GetVP().pix_height / 2;
    int xp = (m_pParentCanvas->GetVP().pix_width - w) / 2;
 
    glColor3ub(243, 229, 47);
 
    glBegin(GL_QUADS);
    glVertex2i(xp, yp);
    glVertex2i(xp + w, yp);
    glVertex2i(xp + w, yp + h);
    glVertex2i(xp, yp + h);
    glEnd();
 
    glEnable(GL_BLEND);
 
    glColor3ub(0, 0, 0);
    glEnable(GL_TEXTURE_2D);
    texfont.RenderString(msg, xp, yp);
    glDisable(GL_TEXTURE_2D);
    glDisable(GL_BLEND);
  }
#endif
}
 
GLShaderProgram *pStaticShader;
 
static std::list<double *> combine_work_data;
static void combineCallbackD(GLdouble coords[3], GLdouble *vertex_data[4],
                             GLfloat weight[4], GLdouble **dataOut) {
  double *vertex = new double[3];
  combine_work_data.push_back(vertex);
  memcpy(vertex, coords, 3 * (sizeof *coords));
  *dataOut = vertex;
}
 
#if defined(USE_ANDROID_GLES2) || defined(ocpnUSE_GLSL)
void vertexCallbackD_GLSL(GLvoid *vertex) {
  GLdouble *pointer = (GLdouble *)vertex;
  s_tess_vertex_work.push_back(pointer[0]);
  s_tess_vertex_work.push_back(pointer[1]);
}
 
void beginCallbackD_GLSL(GLenum mode) {
  s_tess_mode = mode;
  s_tess_vertex_work.clear();
}
 
void endCallbackD_GLSL() {
  if (!s_tess_vertex_work.empty() && (s_tess_vertex_work.size() % 2 == 0)) {
    GLShaderProgram *shader = pStaticShader;
    shader->Bind();
 
    shader->SetUniformMatrix4fv("MVMatrix",
                                (GLfloat *)s_tessVP.vp_matrix_transform);
 
    mat4x4 identityMatrix;
    mat4x4_identity(identityMatrix);
    shader->SetUniformMatrix4fv("TransformMatrix", (GLfloat *)identityMatrix);
 
    // Use color stored in static variable.
    float colorv[4];
    colorv[0] = s_regionColor.Red() / float(256);
    colorv[1] = s_regionColor.Green() / float(256);
    colorv[2] = s_regionColor.Blue() / float(256);
    colorv[3] = s_regionColor.Alpha() / float(256);
    shader->SetUniform4fv("color", colorv);
    shader->SetAttributePointerf("position", s_tess_vertex_work.data());
 
    glDrawArrays(s_tess_mode, 0, s_tess_vertex_work.size() / 2);
 
    shader->UnBind();
  } else
    wxLogError(
        "::endCallbackD_GLSL() called with wrong tess vertex vector length "
        "(%d).",
        s_tess_vertex_work.size());
}
#else
void vertexCallbackD(GLvoid *vertex) { glVertex3dv((GLdouble *)vertex); }
 
void beginCallbackD(GLenum mode) { glBegin(mode); }
 
void endCallbackD() { glEnd(); }
 
#endif
 
void glChartCanvas::DrawRegion(ViewPort &vp, const LLRegion &region) {
  float lat_dist, lon_dist;
  GetLatLonCurveDist(vp, lat_dist, lon_dist);
 
  GLUtesselator *tobj = gluNewTess();
  if (!pStaticShader) pStaticShader = GetStaticTriShader();
 
#if defined(USE_ANDROID_GLES2) || defined(ocpnUSE_GLSL)
  gluTessCallback(tobj, GLU_TESS_VERTEX, (_GLUfuncptr)&vertexCallbackD_GLSL);
  gluTessCallback(tobj, GLU_TESS_BEGIN, (_GLUfuncptr)&beginCallbackD_GLSL);
  gluTessCallback(tobj, GLU_TESS_END, (_GLUfuncptr)&endCallbackD_GLSL);
  gluTessCallback(tobj, GLU_TESS_COMBINE, (_GLUfuncptr)&combineCallbackD);
  s_tessVP = vp;
 
#else
  gluTessCallback(tobj, GLU_TESS_VERTEX, (_GLUfuncptr)&vertexCallbackD);
  gluTessCallback(tobj, GLU_TESS_BEGIN, (_GLUfuncptr)&beginCallbackD);
  gluTessCallback(tobj, GLU_TESS_END, (_GLUfuncptr)&endCallbackD);
  gluTessCallback(tobj, GLU_TESS_COMBINE, (_GLUfuncptr)&combineCallbackD);
#endif
 
  gluTessNormal(tobj, 0, 0, 1);
 
  gluTessBeginPolygon(tobj, NULL);
  for (auto i = region.contours.begin(); i != region.contours.end(); i++) {
    gluTessBeginContour(tobj);
    contour_pt l = *i->rbegin();
    double sml[2];
    bool sml_valid = false;
    for (poly_contour::const_iterator j = i->begin(); j != i->end(); j++) {
      int lat_splits = floor(fabs(j->y - l.y) / lat_dist);
      int lon_splits = floor(fabs(j->x - l.x) / lon_dist);
      int splits = wxMax(lat_splits, lon_splits) + 1;
 
      double smj[2];
      if (splits != 1) {
        // must perform border interpolation in mercator space as this is what
        // the charts use
        toSM(j->y, j->x, 0, 0, smj + 0, smj + 1);
        if (!sml_valid) toSM(l.y, l.x, 0, 0, sml + 0, sml + 1);
      }
 
      for (int i = 0; i < splits; i++) {
        double lat, lon;
        if (i == splits - 1)
          lat = j->y, lon = j->x;
        else {
          double d = (double)(i + 1) / splits;
          fromSM(d * smj[0] + (1 - d) * sml[0], d * smj[1] + (1 - d) * sml[1],
                 0, 0, &lat, &lon);
        }
        wxPoint2DDouble q = vp.GetDoublePixFromLL(lat, lon);
        if (std::isnan(q.m_x)) continue;
 
        double *p = new double[6];
 
        // p[0] = q.m_x, p[1] = q.m_y, p[2] = 0;
        //  It is reasonable to use wxRound() here,
        //  since we are working with pixel coordinates at this point
        p[0] = wxRound(q.m_x), p[1] = wxRound(q.m_y), p[2] = 0;
 
        // wxPoint pt = vp.GetPixFromLL(lat, lon);
        // p[0] = pt.x, p[1] = pt.y, p[2] = 0;
 
        gluTessVertex(tobj, p, p);
        combine_work_data.push_back(p);
      }
      l = *j;
 
      if ((sml_valid = splits != 1)) memcpy(sml, smj, sizeof smj);
    }
    gluTessEndContour(tobj);
  }
  gluTessEndPolygon(tobj);
 
  gluDeleteTess(tobj);
 
  for (std::list<double *>::iterator i = combine_work_data.begin();
       i != combine_work_data.end(); i++)
    delete[] *i;
  combine_work_data.clear();
}
 
/* set stencil buffer to clip in this region, and optionally clear using the
 * current color */
void glChartCanvas::SetClipRegion(ViewPort &vp, const LLRegion &region) {
  glColorMask(GL_FALSE, GL_FALSE, GL_FALSE, GL_FALSE);  // disable color buffer
 
  if (s_b_useStencil) {
    //    Create a stencil buffer for clipping to the region
    glEnable(GL_STENCIL_TEST);
    glStencilMask(0x1);  // write only into bit 0 of the stencil buffer
    glClear(GL_STENCIL_BUFFER_BIT);
 
    //    We are going to write "1" into the stencil buffer wherever the region
    //    is valid
    glStencilFunc(GL_ALWAYS, 1, 1);
    glStencilOp(GL_KEEP, GL_KEEP, GL_REPLACE);
  }
// #ifndef USE_ANDROID_GLES2
#if !defined(USE_ANDROID_GLES2) && !defined(ocpnUSE_GLSL)
 
  else  //  Use depth buffer for clipping
  {
    glEnable(GL_DEPTH_TEST);  // to enable writing to the depth buffer
    glDepthFunc(GL_ALWAYS);   // to ensure everything you draw passes
    glDepthMask(GL_TRUE);     // to allow writes to the depth buffer
 
    glClear(GL_DEPTH_BUFFER_BIT);  // for a fresh start
 
    //    Decompose the region into rectangles, and draw as quads
    //    With z = 1
    // dep buffer clear = 1
    // 1 makes 0 in dep buffer, works
    // 0 make .5 in depth buffer
    // -1 makes 1 in dep buffer
 
    //    Depth buffer runs from 0 at z = 1 to 1 at z = -1
    //    Draw the clip geometry at z = 0.5, giving a depth buffer value of 0.25
    //    Subsequent drawing at z=0 (depth = 0.5) will pass if using
    //    glDepthFunc(GL_GREATER);
    glTranslatef(0, 0, .5);
  }
#endif
 
  s_regionColor = wxColor(0, 0, 0, 255);
  DrawRegion(vp, region);
 
  if (s_b_useStencil) {
    //    Now set the stencil ops to subsequently render only where the stencil
    //    bit is "1"
    glStencilFunc(GL_EQUAL, 1, 1);
    glStencilOp(GL_KEEP, GL_KEEP, GL_KEEP);
  }
// #ifndef USE_ANDROID_GLES2
#if !defined(USE_ANDROID_GLES2) && !defined(ocpnUSE_GLSL)
  else {
    glDepthFunc(GL_GREATER);  // Set the test value
    glDepthMask(GL_FALSE);    // disable depth buffer
    glTranslatef(0, 0, -.5);  // reset translation
  }
#endif
  glColorMask(GL_TRUE, GL_TRUE, GL_TRUE, GL_TRUE);  // re-enable color buffer
}
 
void glChartCanvas::SetClipRect(const ViewPort &vp, const wxRect &rect,
                                bool b_clear) {
  /* for some reason this causes an occasional bug in depth mode, I cannot
     seem to solve it yet, so for now: */
  if (s_b_useStencil && s_b_useScissorTest) {
    wxRect vp_rect(0, 0, vp.pix_width, vp.pix_height);
    if (rect != vp_rect) {
      glEnable(GL_SCISSOR_TEST);
      glScissor(rect.x, vp.pix_height - rect.height - rect.y, rect.width,
                rect.height);
    }
#ifndef USE_ANDROID_GLES2
#endif
    return;
  }
}
 
void glChartCanvas::DisableClipRegion() {
  glDisable(GL_SCISSOR_TEST);
  glDisable(GL_STENCIL_TEST);
  glDisable(GL_DEPTH_TEST);
}
 
void glChartCanvas::Invalidate() {
  /* should probably use a different flag for this */
  m_cache_vp.Invalidate();
}
 
void glChartCanvas::RenderRasterChartRegionGL(ChartBase *chart, ViewPort &vp,
                                              LLRegion &region) {
  ChartBaseBSB *pBSBChart = dynamic_cast<ChartBaseBSB *>(chart);
  if (!pBSBChart) return;
 
  if (b_inCompressAllCharts)
    return;  // don't want multiple texfactories to exist
 
  //    Look for the texture factory for this chart
  wxString key = chart->GetHashKey();
 
  glTexFactory *pTexFact;
  ChartPathHashTexfactType &hash = g_glTextureManager->m_chart_texfactory_hash;
  ChartPathHashTexfactType::iterator ittf = hash.find(key);
 
  //    Not Found ?
  if (ittf == hash.end()) {
    hash[key] = new glTexFactory(chart, g_raster_format);
    hash[key]->SetHashKey(key);
  }
 
  pTexFact = hash[key];
  pTexFact->SetLRUTime(++m_LRUtime);
 
  // for small scales, don't use normalized coordinates for accuracy (difference
  // is up to 3 meters error)
  bool use_norm_vp =
      glChartCanvas::HasNormalizedViewPort(vp) && pBSBChart->GetPPM() < 1;
  pTexFact->PrepareTiles(vp, use_norm_vp, pBSBChart);
 
  //    For underzoom cases, we will define the textures as having their base
  //    levels equivalent to a level "n" mipmap, where n is calculated, and is
  //    always binary This way we can avoid loading much texture memory
 
  int base_level;
  if (vp.m_projection_type == PROJECTION_MERCATOR &&
      chart->GetChartProjectionType() == PROJECTION_MERCATOR) {
    double scalefactor = pBSBChart->GetRasterScaleFactor(vp);
    base_level = log(scalefactor) / log(2.0);
 
    if (base_level < 0) /* for overzoom */
      base_level = 0;
    if (base_level > g_mipmap_max_level) base_level = g_mipmap_max_level;
  } else
    base_level = 0;  // base level should be computed per tile, for now load all
 
  /* setup opengl parameters */
  glEnable(GL_TEXTURE_2D);
#if !defined(USE_ANDROID_GLES2) && !defined(ocpnUSE_GLSL)
  glTexEnvi(GL_TEXTURE_ENV, GL_TEXTURE_ENV_MODE, GL_REPLACE);
 
  glEnableClientState(GL_VERTEX_ARRAY);
  glEnableClientState(GL_TEXTURE_COORD_ARRAY);
 
  if (use_norm_vp) {
    glPushMatrix();
    double lat, lon;
    pTexFact->GetCenter(lat, lon);
    MultMatrixViewPort(vp, lat, lon);
  }
#endif
 
  LLBBox box = region.GetBox();
  int numtiles;
  int mem_used = 0;
  if (g_memCacheLimit > 0) {
    // GetMemoryStatus is slow on linux
    platform::GetMemoryStatus(0, &mem_used);
  }
 
  glTexTile **tiles = pTexFact->GetTiles(numtiles);
  for (int i = 0; i < numtiles; i++) {
    glTexTile *tile = tiles[i];
    if (region.IntersectOut(tile->box)) {
      /*   user setting is in MB while we count exact bytes */
      bool bGLMemCrunch =
          g_tex_mem_used > g_GLOptions.m_iTextureMemorySize * 1024 * 1024;
      if (bGLMemCrunch) pTexFact->DeleteTexture(tile->rect);
    } else {
      bool texture = pTexFact->PrepareTexture(base_level, tile->rect,
                                              global_color_scheme, mem_used);
 
      float *coords;
      if (use_norm_vp)
        coords = tile->m_coords;
      else {
        coords = new float[2 * tile->m_ncoords];
        for (int i = 0; i < tile->m_ncoords; i++) {
          wxPoint2DDouble p = vp.GetDoublePixFromLL(tile->m_coords[2 * i + 0],
                                                    tile->m_coords[2 * i + 1]);
          coords[2 * i + 0] = p.m_x;
          coords[2 * i + 1] = p.m_y;
        }
      }
 
#if defined(USE_ANDROID_GLES2) || defined(ocpnUSE_GLSL)
      RenderTextures(m_gldc, coords, tile->m_texcoords, 4,
                     m_pParentCanvas->GetpVP());
#else
      if (!texture) {  // failed to load, draw red
        glDisable(GL_TEXTURE_2D);
        glColor3f(1, 0, 0);
      }
 
      glTexCoordPointer(2, GL_FLOAT, 2 * sizeof(GLfloat), tile->m_texcoords);
      glVertexPointer(2, GL_FLOAT, 2 * sizeof(GLfloat), coords);
      glDrawArrays(GL_QUADS, 0, tile->m_ncoords);
#endif
      if (!texture) glEnable(GL_TEXTURE_2D);
 
      if (!use_norm_vp) delete[] coords;
    }
  }
 
  glDisable(GL_TEXTURE_2D);
 
#if !defined(USE_ANDROID_GLES2) && !defined(ocpnUSE_GLSL)
  if (use_norm_vp) glPopMatrix();
 
  glDisableClientState(GL_TEXTURE_COORD_ARRAY);
  glDisableClientState(GL_VERTEX_ARRAY);
#endif
}
 
void glChartCanvas::RenderQuiltViewGL(ViewPort &vp,
                                      const OCPNRegion &rect_region) {
  if (!m_pParentCanvas->m_pQuilt->GetnCharts() ||
      m_pParentCanvas->m_pQuilt->IsBusy())
    return;
 
  //  render the quilt
  ChartBase *chart = m_pParentCanvas->m_pQuilt->GetFirstChart();
  if (!chart) {
    printf("   Chart NULL\n");
    chart = m_pParentCanvas->m_pQuilt->GetFirstChart();
  }
 
  //  Check the first, smallest scale chart
  if (chart) {
    //            if( ! m_pParentCanvas->IsChartLargeEnoughToRender( chart, vp )
    //            )
    //            chart = NULL;
  }
 
  LLRegion region = vp.GetLLRegion(rect_region);
 
  LLRegion rendered_region;
  while (chart) {
    //  This test does not need to be done for raster charts, since
    //  we can assume that texture binding is acceptably fast regardless of the
    //  render region, and that the quilt zoom methods choose a reasonable
    //  reference chart.
    if (chart->GetChartFamily() != CHART_FAMILY_RASTER) {
      //                if( ! m_pParentCanvas->IsChartLargeEnoughToRender(
      //                chart, vp ) ) {
      //                    chart = m_pParentCanvas->m_pQuilt->GetNextChart();
      //                    continue;
      //                }
    }
 
    QuiltPatch *pqp = m_pParentCanvas->m_pQuilt->GetCurrentPatch();
    if (pqp->b_Valid) {
      LLRegion get_region = pqp->ActiveRegion;
      bool b_rendered = false;
 
      if (!pqp->b_overlay) {
        get_region.Intersect(region);
        if (!get_region.Empty()) {
          if (chart->GetChartFamily() == CHART_FAMILY_RASTER) {
            ChartBaseBSB *Patch_Ch_BSB = dynamic_cast<ChartBaseBSB *>(chart);
            if (Patch_Ch_BSB) {
              SetClipRegion(vp, get_region /*pqp->quilt_region*/);
              RenderRasterChartRegionGL(chart, vp, pqp->ActiveRegion);
              DisableClipRegion();
 
              b_rendered = true;
            } else if (chart->GetChartType() == CHART_TYPE_MBTILES) {
              SetClipRegion(vp, pqp->ActiveRegion /*pqp->quilt_region*/);
              chart->RenderRegionViewOnGL(*m_pcontext, vp, rect_region,
                                          get_region);
              DisableClipRegion();
            }
 
          } else if (chart->GetChartFamily() == CHART_FAMILY_VECTOR) {
            if (chart->GetChartType() == CHART_TYPE_CM93COMP) {
              RenderNoDTA(vp, get_region);
              chart->RenderRegionViewOnGL(*m_pcontext, vp, rect_region,
                                          get_region);
            } else {
              s57chart *Chs57 = dynamic_cast<s57chart *>(chart);
              if (Chs57) {
                if (Chs57->m_RAZBuilt) {
                  RenderNoDTA(vp, get_region);
                  Chs57->RenderRegionViewOnGLNoText(*m_pcontext, vp,
                                                    rect_region, get_region);
                  DisableClipRegion();
                } else {
                  // The SENC is quesed for building, so..
                  // Show GSHHS with compatible color scheme in the meantime.
                  ocpnDC gldc(*this);
                  const LLRegion &oregion = get_region;
                  LLBBox box = oregion.GetBox();
 
                  wxPoint p1 =
                      vp.GetPixFromLL(box.GetMaxLat(), box.GetMinLon());
                  wxPoint p2 =
                      vp.GetPixFromLL(box.GetMaxLat(), box.GetMaxLon());
                  wxPoint p3 =
                      vp.GetPixFromLL(box.GetMinLat(), box.GetMaxLon());
                  wxPoint p4 =
                      vp.GetPixFromLL(box.GetMinLat(), box.GetMinLon());
 
                  wxRect srect(p1.x, p1.y, p3.x - p1.x, p4.y - p2.y);
 
                  bool world = false;
                  ViewPort cvp = ClippedViewport(vp, get_region);
                  if (m_pParentCanvas->GetWorldBackgroundChart()) {
                    SetClipRegion(cvp, get_region);
                    m_pParentCanvas->GetWorldBackgroundChart()->SetColorsDirect(
                        GetGlobalColor("LANDA"), GetGlobalColor("DEPMS"));
                    RenderWorldChart(gldc, cvp, srect, world);
                    m_pParentCanvas->GetWorldBackgroundChart()->SetColorScheme(
                        global_color_scheme);
                    DisableClipRegion();
                  }
                }
              } else {
                ChartPlugInWrapper *ChPI =
                    dynamic_cast<ChartPlugInWrapper *>(chart);
                if (ChPI) {
                  SetClipRegion(vp, get_region);
                  RenderNoDTA(vp, get_region);
                  ChPI->RenderRegionViewOnGLNoText(*m_pcontext, vp, rect_region,
                                                   get_region);
                  DisableClipRegion();
 
                } else {
                  SetClipRegion(vp, get_region);
                  RenderNoDTA(vp, get_region);
                  chart->RenderRegionViewOnGL(*m_pcontext, vp, rect_region,
                                              get_region);
                  DisableClipRegion();
                }
              }
            }
          }
        }
      }
 
      if (b_rendered) {
        //                LLRegion get_region = pqp->ActiveRegion;
        //                    get_region.Intersect( Region );  not technically
        //                    required?
        //                rendered_region.Union(get_region);
      }
    }
 
    chart = m_pParentCanvas->m_pQuilt->GetNextChart();
  }
 
  //    Render any Overlay patches for s57 charts(cells)
  if (m_pParentCanvas->m_pQuilt->HasOverlays()) {
    ChartBase *pch = m_pParentCanvas->m_pQuilt->GetFirstChart();
    while (pch) {
      QuiltPatch *pqp = m_pParentCanvas->m_pQuilt->GetCurrentPatch();
      if (pqp->b_Valid && pqp->b_overlay &&
          pch->GetChartFamily() == CHART_FAMILY_VECTOR) {
        LLRegion get_region = pqp->ActiveRegion;
 
        get_region.Intersect(region);
        if (!get_region.Empty()) {
          s57chart *Chs57 = dynamic_cast<s57chart *>(pch);
          if (Chs57)
            Chs57->RenderOverlayRegionViewOnGL(*m_pcontext, vp, rect_region,
                                               get_region);
          else {
            ChartPlugInWrapper *ChPI = dynamic_cast<ChartPlugInWrapper *>(pch);
            if (ChPI) {
              ChPI->RenderRegionViewOnGL(*m_pcontext, vp, rect_region,
                                         get_region);
            }
          }
        }
      }
 
      pch = m_pParentCanvas->m_pQuilt->GetNextChart();
    }
  }
 
  // Hilite rollover of standard chart key
  ViewPort vph = m_pParentCanvas->GetVP();
  for (auto &index : m_pParentCanvas->m_pQuilt->GetHiLiteIndexArray()) {
    const ChartTableEntry &cte = ChartData->GetChartTableEntry(index);
    LLRegion hiregion =
        m_pParentCanvas->m_pQuilt->GetChartQuiltRegion(cte, vph);
 
    if (!hiregion.Empty()) {
      glEnable(GL_BLEND);
 
      double hitrans;
      switch (global_color_scheme) {
        case GLOBAL_COLOR_SCHEME_DAY:
          hitrans = .4;
          break;
        case GLOBAL_COLOR_SCHEME_DUSK:
          hitrans = .2;
          break;
        case GLOBAL_COLOR_SCHEME_NIGHT:
          hitrans = .1;
          break;
        default:
          hitrans = .4;
          break;
      }
 
#if !defined(USE_ANDROID_GLES2) && !defined(ocpnUSE_GLSL)
 
      glColor4f((float).8, (float).4, (float).4, (float)hitrans);
#else
      s_regionColor = wxColor(204, 102, 102, hitrans * 256);
#endif
 
      DrawRegion(vp, hiregion);
 
      glDisable(GL_BLEND);
    }
  }
 
#if 0
  LLRegion hiregion = m_pParentCanvas->m_pQuilt->GetHiliteRegion();
 
  if (!hiregion.Empty()) {
    glEnable(GL_BLEND);
 
    double hitrans;
    switch (global_color_scheme) {
      case GLOBAL_COLOR_SCHEME_DAY:
        hitrans = .4;
        break;
      case GLOBAL_COLOR_SCHEME_DUSK:
        hitrans = .2;
        break;
      case GLOBAL_COLOR_SCHEME_NIGHT:
        hitrans = .1;
        break;
      default:
        hitrans = .4;
        break;
    }
 
//#ifndef USE_ANDROID_GLES2
#if !defined(USE_ANDROID_GLES2) && !defined(ocpnUSE_GLSL)
 
    glColor4f((float).8, (float).4, (float).4, (float)hitrans);
#else
    s_regionColor = wxColor(204, 102, 102, hitrans * 256);
#endif
 
    DrawRegion(vp, hiregion);
 
    glDisable(GL_BLEND);
  }
#endif
 
  m_pParentCanvas->m_pQuilt->SetRenderedVP(vp);
}
 
void glChartCanvas::RenderQuiltViewGLText(ViewPort &vp,
                                          const OCPNRegion &rect_region) {
  if (!m_pParentCanvas->m_pQuilt->GetnCharts() ||
      m_pParentCanvas->m_pQuilt->IsBusy())
    return;
 
  //  render the quilt
  ChartBase *chart = m_pParentCanvas->m_pQuilt->GetLargestScaleChart();
 
  LLRegion region = vp.GetLLRegion(rect_region);
 
  LLRegion rendered_region;
  while (chart) {
    QuiltPatch *pqp = m_pParentCanvas->m_pQuilt->GetCurrentPatch();
    if (pqp->b_Valid) {
      LLRegion get_region = pqp->ActiveRegion;
 
      if (!pqp->b_overlay) {
        if (chart->GetChartFamily() == CHART_FAMILY_VECTOR) {
          s57chart *Chs57 = dynamic_cast<s57chart *>(chart);
          if (Chs57) {
            Chs57->RenderViewOnGLTextOnly(*m_pcontext, vp);
          } else {
            ChartPlugInWrapper *ChPI =
                dynamic_cast<ChartPlugInWrapper *>(chart);
            if (ChPI) {
              ChPI->RenderRegionViewOnGLTextOnly(*m_pcontext, vp, rect_region);
            }
          }
        }
      }
    }
 
    chart = m_pParentCanvas->m_pQuilt->GetNextSmallerScaleChart();
  }
 
  /*
          //    Render any Overlay patches for s57 charts(cells)
          if( m_pParentCanvas->m_pQuilt->HasOverlays() ) {
              ChartBase *pch = m_pParentCanvas->m_pQuilt->GetFirstChart();
              while( pch ) {
                  QuiltPatch *pqp =
     m_pParentCanvas->m_pQuilt->GetCurrentPatch(); if( pqp->b_Valid &&
     pqp->b_overlay && pch->GetChartFamily() == CHART_FAMILY_VECTOR ) { LLRegion
     get_region = pqp->ActiveRegion;
 
                      get_region.Intersect( region );
                      if( !get_region.Empty()  ) {
                          s57chart *Chs57 = dynamic_cast<s57chart*>( pch );
                          if( Chs57 )
                              Chs57->RenderOverlayRegionViewOnGL( *m_pcontext,
     vp, rect_region, get_region );
                      }
                  }
 
                  pch = m_pParentCanvas->m_pQuilt->GetNextChart();
              }
          }
  */
}
 
void glChartCanvas::RenderCharts(ocpnDC &dc, const OCPNRegion &rect_region) {
  ViewPort &vp = m_pParentCanvas->VPoint;
 
  // Only for cm93 (not quilted), SetVPParms can change the valid region of the
  // chart we need to know this before rendering the chart so we can compute the
  // background region and nodta regions correctly.  I would prefer to just
  // perform this here (or in SetViewPoint) for all vector charts instead of in
  // their render routine, but how to handle quilted cases?
  if (!vp.b_quilt &&
      m_pParentCanvas->m_singleChart->GetChartType() == CHART_TYPE_CM93COMP)
    static_cast<cm93compchart *>(m_pParentCanvas->m_singleChart)
        ->SetVPParms(vp);
 
  LLRegion chart_region;
  if (!vp.b_quilt &&
      (m_pParentCanvas->m_singleChart->GetChartType() == CHART_TYPE_PLUGIN)) {
    if (m_pParentCanvas->m_singleChart->GetChartFamily() ==
        CHART_FAMILY_RASTER) {
      // We do this the hard way, since PlugIn Raster charts do not understand
      // LLRegion yet...
      double ll[8];
      ChartPlugInWrapper *cpw =
          dynamic_cast<ChartPlugInWrapper *>(m_pParentCanvas->m_singleChart);
      if (!cpw) return;
 
      cpw->chartpix_to_latlong(0, 0, ll + 0, ll + 1);
      cpw->chartpix_to_latlong(0, cpw->GetSize_Y(), ll + 2, ll + 3);
      cpw->chartpix_to_latlong(cpw->GetSize_X(), cpw->GetSize_Y(), ll + 4,
                               ll + 5);
      cpw->chartpix_to_latlong(cpw->GetSize_X(), 0, ll + 6, ll + 7);
 
      // for now don't allow raster charts to cross both 0 meridian and IDL
      // (complicated to deal with)
      for (int i = 1; i < 6; i += 2)
        if (fabs(ll[i] - ll[i + 2]) > 180) {
          // we detect crossing idl here, make all longitudes positive
          for (int i = 1; i < 8; i += 2)
            if (ll[i] < 0) ll[i] += 360;
          break;
        }
 
      chart_region = LLRegion(4, ll);
    } else {
      Extent ext;
      m_pParentCanvas->m_singleChart->GetChartExtent(&ext);
 
      double ll[8] = {ext.SLAT, ext.WLON, ext.SLAT, ext.ELON,
                      ext.NLAT, ext.ELON, ext.NLAT, ext.WLON};
      chart_region = LLRegion(4, ll);
    }
  } else
    chart_region = vp.b_quilt
                       ? m_pParentCanvas->m_pQuilt->GetFullQuiltRegion()
                       : m_pParentCanvas->m_singleChart->GetValidRegion();
 
  bool world_view = false;
  for (OCPNRegionIterator upd(rect_region); upd.HaveRects(); upd.NextRect()) {
    wxRect rect = upd.GetRect();
    LLRegion background_region = vp.GetLLRegion(rect);
    //    Remove the valid chart area to find the region NOT covered by the
    //    charts
    background_region.Subtract(chart_region);
 
    if (!background_region.Empty()) {
      ViewPort cvp = ClippedViewport(vp, background_region);
      SetClipRect(cvp, rect, false);
      RenderWorldChart(dc, cvp, rect, world_view);
      DisableClipRegion();
    }
  }
 
  if (vp.b_quilt) {
    RenderQuiltViewGL(vp, rect_region);
    int yyp = 4;
  } else {
    LLRegion region = vp.GetLLRegion(rect_region);
    if (m_pParentCanvas->m_singleChart->GetChartFamily() ==
        CHART_FAMILY_RASTER) {
      if (m_pParentCanvas->m_singleChart->GetChartType() == CHART_TYPE_MBTILES)
        m_pParentCanvas->m_singleChart->RenderRegionViewOnGL(
            *m_pcontext, vp, rect_region, region);
      else
        RenderRasterChartRegionGL(m_pParentCanvas->m_singleChart, vp, region);
    } else if (m_pParentCanvas->m_singleChart->GetChartFamily() ==
               CHART_FAMILY_VECTOR) {
      chart_region.Intersect(region);
      RenderNoDTA(vp, chart_region);
      m_pParentCanvas->m_singleChart->RenderRegionViewOnGL(*m_pcontext, vp,
                                                           rect_region, region);
    }
  }
  glUseProgram(0);
}
 
void glChartCanvas::RenderNoDTA(ViewPort &vp, const LLRegion &region,
                                int transparency) {
  wxColour color = GetGlobalColor("NODTA");
#if !defined(USE_ANDROID_GLES2) && !defined(ocpnUSE_GLSL)
  if (color.IsOk())
    glColor4ub(color.Red(), color.Green(), color.Blue(), transparency);
  else
    glColor4ub(163, 180, 183, transparency);
 
  glEnable(GL_BLEND);
  glBlendFunc(GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
 
#else
  // Store the color for tesselator callback pickup.
  s_regionColor = color;
#endif
 
  DrawRegion(vp, region);
}
 
/* render world chart, but only in this rectangle */
void glChartCanvas::RenderWorldChart(ocpnDC &dc, ViewPort &vp, wxRect &rect,
                                     bool &world_view) {
  // set gl color to water
  wxColour water = m_pParentCanvas->pWorldBackgroundChart->water;
 
  glEnable(GL_SCISSOR_TEST);
  glScissor(rect.x, vp.pix_height - rect.height - rect.y, rect.width,
            rect.height);
 
  // clear background
  if (!world_view) {
    if (!world_view) {
      int x1 = rect.x, y1 = rect.y, x2 = x1 + rect.width, y2 = y1 + rect.height;
#if defined(USE_ANDROID_GLES2) || defined(ocpnUSE_GLSL)
 
      GLShaderProgram *shader = pcolor_tri_shader_program[GetCanvasIndex()];
      shader->Bind();
 
      float colorv[4];
      colorv[0] = water.Red() / float(256);
      colorv[1] = water.Green() / float(256);
      colorv[2] = water.Blue() / float(256);
      colorv[3] = 1.0;
      shader->SetUniform4fv("color", colorv);
 
      float pf[8];
      pf[0] = x2;
      pf[1] = y1;
      pf[2] = x2;
      pf[3] = y2;
      pf[4] = x1;
      pf[5] = y1;
      pf[6] = x1;
      pf[7] = y2;
      shader->SetAttributePointerf("position", pf);
 
      glDrawArrays(GL_TRIANGLE_STRIP, 0, 4);
 
      shader->UnBind();
 
#else
#endif
    }
  }
 
  // m_pParentCanvas->pWorldBackgroundChart->RenderViewOnDC(dc, vp);
  gShapeBasemap.RenderViewOnDC(dc, vp);
 
  glDisable(GL_SCISSOR_TEST);
}
 
/* these are the overlay objects which move with the charts and
   are not frequently updated (not ships etc..)
 
   many overlay objects are fixed to a geographical location or
   grounded as opposed to the floating overlay objects. */
void glChartCanvas::DrawGroundedOverlayObjects(ocpnDC &dc, ViewPort &vp) {
  m_pParentCanvas->RenderAllChartOutlines(dc, vp);
 
  DrawStaticRoutesTracksAndWaypoints(vp);
 
  DisableClipRegion();
}
 
void glChartCanvas::DrawGLTidesInBBox(ocpnDC &dc, LLBBox &BBox) {
  // At small scale, we render the Tide icon as a texture for best performance
  if (m_pParentCanvas->GetVP().chart_scale > 500000) {
    // Prepare the texture if necessary
 
    if (!m_tideTex) {
      wxBitmap bmp = m_pParentCanvas->GetTideBitmap();
      if (!bmp.Ok()) return;
 
      wxImage image = bmp.ConvertToImage();
      int w = image.GetWidth(), h = image.GetHeight();
 
      int tex_w, tex_h;
      if (g_texture_rectangle_format == GL_TEXTURE_2D)
        tex_w = w, tex_h = h;
      else
        tex_w = NextPow2(w), tex_h = NextPow2(h);
 
      m_tideTexWidth = tex_w;
      m_tideTexHeight = tex_h;
 
      unsigned char *d = image.GetData();
      unsigned char *a = image.GetAlpha();
 
      unsigned char mr, mg, mb;
      if (!a) image.GetOrFindMaskColour(&mr, &mg, &mb);
 
      unsigned char *e = new unsigned char[4 * w * h];
      if (e && d) {
        for (int y = 0; y < h; y++)
          for (int x = 0; x < w; x++) {
            unsigned char r, g, b;
            int off = (y * w + x);
            r = d[off * 3 + 0];
            g = d[off * 3 + 1];
            b = d[off * 3 + 2];
 
            e[off * 4 + 0] = r;
            e[off * 4 + 1] = g;
            e[off * 4 + 2] = b;
 
            e[off * 4 + 3] =
                a ? a[off] : ((r == mr) && (g == mg) && (b == mb) ? 0 : 255);
          }
      }
 
      glGenTextures(1, &m_tideTex);
 
      glBindTexture(GL_TEXTURE_2D, m_tideTex);
      glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);
      glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);
 
      if (g_texture_rectangle_format == GL_TEXTURE_2D)
        glTexImage2D(GL_TEXTURE_2D, 0, GL_RGBA, w, h, 0, GL_RGBA,
                     GL_UNSIGNED_BYTE, e);
      else {
        glTexImage2D(GL_TEXTURE_2D, 0, GL_RGBA, tex_w, tex_h, 0, GL_RGBA,
                     GL_UNSIGNED_BYTE, 0);
        glTexSubImage2D(GL_TEXTURE_2D, 0, 0, 0, w, h, GL_RGBA, GL_UNSIGNED_BYTE,
                        e);
      }
 
      delete[] e;
    }
 
    // Texture is ready
 
    glBindTexture(GL_TEXTURE_2D, m_tideTex);
    glEnable(GL_TEXTURE_2D);
    glEnable(GL_BLEND);
 
#if !defined(USE_ANDROID_GLES2) && !defined(ocpnUSE_GLSL)
#else
    for (int i = 1; i < ptcmgr->Get_max_IDX() + 1; i++) {
      const IDX_entry *pIDX = ptcmgr->GetIDX_entry(i);
 
      char type = pIDX->IDX_type;          // Entry "TCtcIUu" identifier
      if ((type == 't') || (type == 'T'))  // only Tides
      {
        double lon = pIDX->IDX_lon;
        double lat = pIDX->IDX_lat;
 
        if (BBox.Contains(lat, lon)) {
          wxPoint r;
          m_pParentCanvas->GetCanvasPointPix(lat, lon, &r);
 
          float xp = r.x;
          float yp = r.y;
 
          double scale = 1.0;
#ifdef __ANDROID__
          scale *= getAndroidDisplayDensity();
#endif
          double width2 = scale * m_tideTexWidth / 2;
          double height2 = scale * m_tideTexHeight / 2;
 
          float coords[8];
          float uv[8];
 
          // normal uv
          uv[0] = 0;
          uv[1] = 0;
          uv[2] = 0;
          uv[3] = 1;
          uv[4] = 1;
          uv[5] = 1;
          uv[6] = 1;
          uv[7] = 0;
 
          coords[0] = xp - width2;
          coords[1] = yp - height2;
          coords[2] = xp - width2;
          coords[3] = yp + height2;
          coords[4] = xp + width2;
          coords[5] = yp + height2;
          coords[6] = xp + width2;
          coords[7] = yp - height2;
 
          RenderTextures(dc, coords, uv, 4, m_pParentCanvas->GetpVP());
        }
      }  // type 'T"
    }  // loop
 
#endif
 
    glDisable(GL_TEXTURE_2D);
    glDisable(GL_BLEND);
    glBindTexture(GL_TEXTURE_2D, 0);
  } else
    m_pParentCanvas->DrawAllTidesInBBox(dc, BBox);
}
 
void glChartCanvas::DrawGLCurrentsInBBox(ocpnDC &dc, LLBBox &BBox) {
  m_pParentCanvas->DrawAllCurrentsInBBox(dc, BBox);
}
 
void glChartCanvas::SetColorScheme(ColorScheme cs) {
  if (!m_bsetup) return;
 
  glDeleteTextures(1, &m_tideTex);
  glDeleteTextures(1, &m_currentTex);
  m_tideTex = 0;
  m_currentTex = 0;
  ownship_color = -1;
}
 
void glChartCanvas::RenderGLAlertMessage() {
  if (!m_pParentCanvas->GetAlertString().IsEmpty()) {
    wxString msg = m_pParentCanvas->GetAlertString();
 
    wxFont *pfont = GetOCPNScaledFont(_("Dialog"));
    m_gldc.SetFont(*pfont);
 
    int w, h;
    wxScreenDC sdc;
    sdc.GetTextExtent(msg, &w, &h, NULL, NULL, pfont);
 
    h += 2;
    w += 4;
    int yp =
        m_pParentCanvas->VPoint.pix_height - GetChartbarHeight() - h - (h / 4);
 
    wxRect sbr = m_pParentCanvas->GetScaleBarRect();
    int xp = sbr.x + sbr.width + 5;
 
    wxPen ppPen1(GetGlobalColor("UBLCK"), 1, wxPENSTYLE_SOLID);
    m_gldc.SetPen(ppPen1);
    m_gldc.SetBrush(wxBrush(GetGlobalColor("YELO1")));
 
    m_gldc.DrawRectangle(xp, yp, w, h);
 
    m_gldc.DrawText(msg, xp, yp);
  }
}
 
unsigned long quiltHash;
int refChartIndex;
extern wxLongLong s_t0;
 
int n_render;
void glChartCanvas::Render() {
  if (!m_bsetup || !m_pParentCanvas->m_pQuilt ||
      (m_pParentCanvas->VPoint.b_quilt && !m_pParentCanvas->m_pQuilt) ||
      (!m_pParentCanvas->VPoint.b_quilt && !m_pParentCanvas->m_singleChart)) {
#ifdef __WXGTK__  // for some reason in gtk, a swap is needed here to get an
                  // initial screen update
    SwapBuffers();
#endif
    if (!g_PrintingInProgress) return;
  }
 
  if (!g_true_zoom && m_binPinch) return;
 
  // if (m_binPinch) printf("    %ld Render Start\n", m_glstopwatch.Time());
  long render_start_time = m_glstopwatch.Time();
 
#if defined(USE_ANDROID_GLES2) || defined(ocpnUSE_GLSL)
  loadShaders(GetCanvasIndex());
  configureShaders(m_pParentCanvas->VPoint);
#endif
 
#ifdef USE_ANDROID_GLES2
 
  OCPNStopWatch sw;
 
  if (m_binPinch) return;
 
  // qDebug() << "Render" << m_pParentCanvas->m_canvasIndex << GetPosition().x
  // << GetSize().x << m_pParentCanvas->GetPosition().x << m_pcontext;
 
  // if(m_pParentCanvas->m_canvasIndex == 0) return;
 
  // Do any setup required...
 
  bool recompose = false;
  if (m_pParentCanvas->VPoint.b_quilt && m_pParentCanvas->m_pQuilt &&
      !m_pParentCanvas->m_pQuilt->IsComposed()) {
    if (m_pParentCanvas->VPoint.IsValid()) {
      m_pParentCanvas->m_pQuilt->Compose(m_pParentCanvas->VPoint);
      m_pParentCanvas->UpdateCanvasControlBar();
      recompose = true;
    } else
      return;
  }
 
  //  Check to see if the Compose() call forced a SENC build.
  //  If so, zoom the canvas just slightly to force a deferred redraw of the
  //  full screen.
  if (sw.GetTime() > 2000) {  //  long enough to detect SENC build.
    m_pParentCanvas->ZoomCanvas(1.0001, false);
  }
 
  // qDebug() << "RenderTime1" << sw.GetTime();
 
  quiltHash = m_pParentCanvas->m_pQuilt->GetXStackHash();
  refChartIndex = m_pParentCanvas->m_pQuilt->GetRefChartdbIndex();
 
#endif
 
#ifdef __WXOSX__
  // Support scaled HDPI displays.
  m_displayScale = GetContentScaleFactor();
#endif
  m_pParentCanvas->VPoint.SetPixelScale(m_displayScale);
 
  m_last_render_time = wxDateTime::Now().GetTicks();
 
  // we don't care about jobs that are now off screen
  // clear out and it will be repopulated during render
  if (g_GLOptions.m_bTextureCompression &&
      !g_GLOptions.m_bTextureCompressionCaching)
    g_glTextureManager->ClearJobList();
 
  ocpnDC gldc(*this);
 
  int gl_width, gl_height;
  gl_width = m_pParentCanvas->VPoint.pix_width;
  gl_height = m_pParentCanvas->VPoint.pix_height;
 
  // Take a copy for use later by DC
  m_glcanvas_width = gl_width;
  m_glcanvas_height = gl_height;
 
  // @todo: If the intention was to work with the same ViewPort object, use a
  // reference instead. Making a copy of VPoint here means that any changes to
  // VPoint will not affect m_pParentCanvas->VPoint. It's not clear if this is
  // the intended behavior.
  ViewPort VPoint = m_pParentCanvas->VPoint;
 
  OCPNRegion screen_region(wxRect(0, 0, gl_width, gl_height));
  glViewport(0, 0, (GLint)gl_width, (GLint)gl_height);
 
// #ifndef USE_ANDROID_GLES2
#if !defined(USE_ANDROID_GLES2)
  glMatrixMode(GL_PROJECTION);
  glLoadIdentity();
 
  glOrtho(0, (GLint)gl_width, (GLint)gl_height, 0, -1, 1);
  glMatrixMode(GL_MODELVIEW);
  glLoadIdentity();
#endif
 
  if (s_b_useStencil) {
    glEnable(GL_STENCIL_TEST);
    glStencilMask(0xff);
    glClear(GL_STENCIL_BUFFER_BIT);
    glDisable(GL_STENCIL_TEST);
  }
 
  // set opengl settings that don't normally change
  // this should be able to go in SetupOpenGL, but it's
  // safer here incase a plugin mangles these
  if (g_GLOptions.m_GLLineSmoothing) glHint(GL_LINE_SMOOTH_HINT, GL_NICEST);
  if (g_GLOptions.m_GLPolygonSmoothing)
    glHint(GL_POLYGON_SMOOTH_HINT, GL_NICEST);
  glBlendFunc(GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
 
  //  Delete any textures known to the GPU that
  //  belong to charts which will not be used in this render
  //  This is done chart-by-chart...later we will scrub for unused textures
  //  that belong to charts which ARE used in this render, if we need to....
 
  g_glTextureManager->TextureCrunch(0.8);
 
  //  If we plan to post process the display, don't use accelerated panning
  double scale_factor = VPoint.ref_scale / VPoint.chart_scale;
 
  bool bpost_hilite = !m_pParentCanvas->m_pQuilt->GetHiliteRegion().Empty();
  bool useFBO = false;
  int sx = gl_width;
  int sy = gl_height;
 
  // Try to use the framebuffer object's cache of the last frame
  // to accelerate drawing this frame (if overlapping)
  if (m_b_BuiltFBO && !bpost_hilite
      //&& VPoint.tilt == 0 // disabling fbo in tilt mode gives better quality
      // but slower
  ) {
    //  Is this viewpoint the same as the previously painted one?
    bool b_newview = true;
    bool b_full = false;
 
    // If the view is the same we do no updates,
    // Just render cached texture to the framebuffer
    if (m_cache_vp.view_scale_ppm == VPoint.view_scale_ppm &&
        m_cache_vp.rotation == VPoint.rotation &&
        m_cache_vp.clat == VPoint.clat && m_cache_vp.clon == VPoint.clon &&
        m_cache_vp.IsValid() && m_cache_vp.pix_height == VPoint.pix_height &&
        m_cache_current_ch == m_pParentCanvas->m_singleChart) {
      b_newview = false;
    }
 
#ifdef USE_ANDROID_GLES2
    if (recompose) b_newview = true;
 
    if (m_bforcefull) {
      b_newview = true;
      b_full = true;
    }
 
    // If no charts are to be rendered, we need to refresh the entire display
    //  This fixes a problem with routes/tracks/marks rendering on pans at very
    //  small scale. It is a workaround, so finding root cause should be
    //  considered a TODO
 
    if (VPoint.b_quilt) {
      ChartBase *chart = m_pParentCanvas->m_pQuilt->GetFirstChart();
      if (!chart) b_full = true;
    }
 
#endif
 
    if (b_newview) {
      float dx = 0;
      float dy = 0;
 
      bool accelerated_pan = false;
      //      if (g_in_inertia)
      //        printf("---     accpan condition %d %d\n",
      //               g_GLOptions.m_bUseAcceleratedPanning,
      //               m_cache_vp.IsValid());
      //      else
      //        printf("|||     accpan condition %d %d\n",
      //               g_GLOptions.m_bUseAcceleratedPanning,
      //               m_cache_vp.IsValid());
 
      if (g_GLOptions.m_bUseAcceleratedPanning && m_cache_vp.IsValid() &&
          (VPoint.m_projection_type == PROJECTION_MERCATOR ||
           VPoint.m_projection_type == PROJECTION_EQUIRECTANGULAR) &&
          m_cache_vp.pix_height == VPoint.pix_height) {
        wxPoint2DDouble c_old =
            VPoint.GetDoublePixFromLL(VPoint.clat, VPoint.clon) *
            m_displayScale;
        wxPoint2DDouble c_new =
            m_cache_vp.GetDoublePixFromLL(VPoint.clat, VPoint.clon) *
            m_displayScale;
 
        dy = wxRound(c_new.m_y - c_old.m_y);
        dx = wxRound(c_new.m_x - c_old.m_x);
 
        //   The math below using the previous frame's texture does not really
        //   work for sub-pixel pans.
        //   TODO is to rethink this.
        //   Meanwhile, require the accelerated pans to be whole pixel multiples
        //   only. This is not as bad as it sounds.  Keyboard and mouse pans are
        //   whole_pixel moves. However, autofollow at large scale is certainly
        //   not.
 
        double deltax = c_new.m_x - c_old.m_x;
        double deltay = c_new.m_y - c_old.m_y;
 
        bool b_whole_pixel = true;
        if ((fabs(deltax - dx) > 1e-2) || (fabs(deltay - dy) > 1e-2))
          b_whole_pixel = false;
 
        accelerated_pan = b_whole_pixel && abs(dx) < m_cache_tex_x &&
                          abs(dy) < m_cache_tex_y &&
                          (abs(dx) > 0 || (abs(dy) > 0));
 
        // if (g_in_inertia && !accelerated_pan)
        // printf("---     accpan %d %d    %g %g\n", accelerated_pan,
        // b_whole_pixel, dx, dy);
      }
 
      //  FBO swapping has trouble with Retina display on MacOS Monterey.
      //  So, disable accelerated pan ops on this case.
      if (m_displayScale > 1) accelerated_pan = false;
 
      // FIXME (dave) There are some display artifact troubles using accPan on
      // rotation.
      //  Especially seen on sparse RNC rendering
      if (fabs(VPoint.rotation) > 0) accelerated_pan = false;
 
        // do we allow accelerated panning?  can we perform it here?
#if !defined(USE_ANDROID_GLES2) && !defined(ocpnUSE_GLSL)
#else  // GLES2
       // enable rendering to texture in framebuffer object
      glBindFramebuffer(GL_FRAMEBUFFER, m_fb0);
 
      if (VPoint.chart_scale < 5000) b_full = true;
 
      if (VPoint.chart_scale > 5e7) b_full = true;
 
      if (b_full) accelerated_pan = false;
 
      if (accelerated_pan) {
        if ((dx != 0) || (dy != 0)) {  // Anything to do?
 
          // calculate the new regions to render
          // add extra pixels to avoid coordindate rounding issues at large
          // scale
          OCPNRegion update_region;
 
          int fluff = 2;
 
          // Avoid rendering artifacts caused by Multi Sampling (MSAA)
          if (VPoint.chart_scale < 10000) fluff = 8;
 
          if (dy > 0 && dy < gl_height)
            update_region.Union(
                wxRect(0, gl_height - (dy + fluff), gl_width, dy + fluff));
          else if (dy < 0)
            update_region.Union(wxRect(0, 0, gl_width, -dy + fluff));
 
          if (dx > 0 && dx < gl_width)
            update_region.Union(
                wxRect(gl_width - (dx + fluff), 0, dx + fluff, gl_height));
          else if (dx < 0)
            update_region.Union(wxRect(0, 0, -dx + fluff, gl_height));
 
          m_cache_page = !m_cache_page; /* page flip */
 
          // Bind the destination (target frame) texture to the frame buffer
          glFramebufferTexture2D(GL_FRAMEBUFFER, GL_COLOR_ATTACHMENT0,
                                 GL_TEXTURE_2D, m_cache_tex[m_cache_page], 0);
 
          // Before rendering anything, clear the color buffers
          //           wxColour color = GetGlobalColor("NODTA");
          //           glClearColor(color.Red() / 256., color.Green() / 256.,
          //                        color.Blue() / 256., 1.0);
          //           glClear(GL_COLOR_BUFFER_BIT);
 
          // First render the new content into the update region
          // if (g_in_inertia)
          // printf("---     R2a %g\n",
          //     (wxGetLocalTimeMillis() - s_t0).ToDouble());
          RenderCharts(m_gldc, update_region);
          // if (g_in_inertia)
          // printf("---     R2b %g\n",
          //      (wxGetLocalTimeMillis() - s_t0).ToDouble());
          glDisable(g_texture_rectangle_format);
          glUseProgram(0);
 
          // Next, render the cached texture as quad to FBO(m_blit_tex) with
          // offsets
          glBindTexture(GL_TEXTURE_2D, m_cache_tex[!m_cache_page]);
          glEnable(GL_TEXTURE_2D);
 
          // Blit the existing content onto the alternate FBO, at the correct
          // location
          float x1, x2, y1, y2;
 
          if (dx > 0)
            x1 = dx, x2 = 0;
          else
            x1 = 0, x2 = -dx;
 
          if (dy > 0)
            y1 = dy, y2 = 0;
          else
            y1 = 0, y2 = -dy;
 
          // normalize to texture coordinates range from 0 to 1
          float tx1, tx2, ty1, ty2;
 
          float xcor = 0;
          float ycor = 0;
 
          tx1 = 0;
          tx2 = sx / (float)m_cache_tex_x;
          ty1 = 0;
          ty2 = sy / (float)m_cache_tex_y;
 
          float coords[8];
          float uv[8];
 
          // normal uv
          uv[0] = tx1;
          uv[1] = ty1;
          uv[2] = tx2;
          uv[3] = ty1;
          uv[4] = tx2;
          uv[5] = ty2;
          uv[6] = tx1;
          uv[7] = ty2;
 
          coords[0] = -dx;
          coords[1] = dy;
          coords[2] = -dx + sx;
          coords[3] = dy;
          coords[4] = -dx + sx;
          coords[5] = dy + sy;
          coords[6] = -dx;
          coords[7] = dy + sy;
 
          GLShaderProgram *shader =
              ptexture_2D_shader_program[GetCanvasIndex()];
          shader->Bind();
 
          // Set up the texture sampler to texture unit 0
          shader->SetUniform1i("uTex", 0);
 
          mat4x4 m, mvp, I;
          mat4x4_identity(m);
          mat4x4_scale_aniso(mvp, m, 2.0 / (float)sx, 2.0 / (float)sy, 1.0);
          mat4x4_translate_in_place(mvp, -(float)sx / 2, -(float)sy / 2, 0);
          shader->SetUniformMatrix4fv("MVMatrix", (GLfloat *)mvp);
          mat4x4_identity(I);
          shader->SetUniformMatrix4fv("TransformMatrix", (GLfloat *)I);
 
          float co1[8];
          co1[0] = coords[0];
          co1[1] = coords[1];
          co1[2] = coords[2];
          co1[3] = coords[3];
          co1[4] = coords[6];
          co1[5] = coords[7];
          co1[6] = coords[4];
          co1[7] = coords[5];
 
          float tco1[8];
          tco1[0] = uv[0];
          tco1[1] = uv[1];
          tco1[2] = uv[2];
          tco1[3] = uv[3];
          tco1[4] = uv[6];
          tco1[5] = uv[7];
          tco1[6] = uv[4];
          tco1[7] = uv[5];
 
          shader->SetAttributePointerf("aPos", co1);
          shader->SetAttributePointerf("aUV", tco1);
 
          glDrawArrays(GL_TRIANGLE_STRIP, 0, 4);
 
          // restore the shader matrix
          shader->SetUniformMatrix4fv("MVMatrix",
                                      (GLfloat *)VPoint.vp_matrix_transform);
 
          shader->UnBind();
          glBindTexture(g_texture_rectangle_format, 0);
 
          glDisable(g_texture_rectangle_format);
          glUseProgram(0);
        }
 
      }  // accelerated pan
 
      else {  // must redraw the entire screen
        // qDebug() << "Fullpage";
        glFramebufferTexture2D(GL_FRAMEBUFFER_EXT, GL_COLOR_ATTACHMENT0,
                               g_texture_rectangle_format,
                               m_cache_tex[!m_cache_page], 0);
 
        m_fbo_offsetx = 0;
        m_fbo_offsety = 0;
        m_fbo_swidth = sx;
        m_fbo_sheight = sy;
 
        // FIXME (dave) test on Android
        //  This can be annoying on Android pinch zoom
 
        // Clear the screen to NODTA color
        wxColour color = GetGlobalColor("NODTA");
        glClearColor(color.Red() / 256., color.Green() / 256.,
                     color.Blue() / 256., 1.0);
        glClear(GL_COLOR_BUFFER_BIT);
 
        OCPNRegion rscreen_region(VPoint.rv_rect);
        // if (g_in_inertia)
        // printf("---     R2c %g\n",
        //      (wxGetLocalTimeMillis() - s_t0).ToDouble());
        RenderCharts(m_gldc, rscreen_region);
        // if (g_in_inertia)
        // printf("---     R2d %g\n",
        //      (wxGetLocalTimeMillis() - s_t0).ToDouble());
 
        m_cache_page = !m_cache_page; /* page flip */
 
      }  // full page render
 
      // Disable Render to FBO
      glBindFramebuffer(GL_FRAMEBUFFER, 0);
 
#endif  // gles2 for accpan
 
    }  // newview
 
    useFBO = true;
  }
 
#ifndef __ANDROID__
  if (VPoint.tilt) {
    glMatrixMode(GL_PROJECTION);
    glLoadIdentity();
 
    gluPerspective(2 * 180 / PI * atan2((double)gl_height, (double)gl_width),
                   (GLfloat)gl_width / (GLfloat)gl_height, 1, gl_width);
 
    glMatrixMode(GL_MODELVIEW);
    glLoadIdentity();
 
    glScalef(1, -1, 1);
    glTranslatef(-gl_width / 2, -gl_height / 2, -gl_width / 2);
    glRotated(VPoint.tilt * 180 / PI, 1, 0, 0);
 
    glGetIntegerv(GL_VIEWPORT, viewport);
    glGetDoublev(GL_PROJECTION_MATRIX, projmatrix);
    glGetDoublev(GL_MODELVIEW_MATRIX, mvmatrix);
  }
#endif
 
  // if (g_in_inertia)
  // printf("---     R3 %g\n", (wxGetLocalTimeMillis() - s_t0).ToDouble());
 
  if (useFBO) {
#if 0  // #ifndef USE_ANDROID_GLES2
    glBindFramebuffer(GL_READ_FRAMEBUFFER, m_fb0);
    glBindFramebuffer(GL_DRAW_FRAMEBUFFER, 0);
    glBlitFramebuffer(0, 0, sx, sy, 0, 0, sx*2, sy*2, GL_COLOR_BUFFER_BIT, GL_LINEAR);
 
    glBindFramebuffer(GL_FRAMEBUFFER, 0);
 
#else
    // Render the cached texture as quad to screen
    glBindTexture(g_texture_rectangle_format, m_cache_tex[m_cache_page]);
    glEnable(g_texture_rectangle_format);
 
    float tx, ty, tx0, ty0, divx, divy;
 
    //  Normalize, or not?
    if (GL_TEXTURE_RECTANGLE_ARB == g_texture_rectangle_format) {
      divx = divy = 1.0f;
    } else {
      divx = m_cache_tex_x;
      divy = m_cache_tex_y;
    }
 
    tx0 = m_fbo_offsetx / divx;
    ty0 = m_fbo_offsety / divy;
    tx = (m_fbo_offsetx + m_fbo_swidth) / divx;
    ty = (m_fbo_offsety + m_fbo_sheight) / divy;
 
    float coords[8];
    float uv[8];
 
    // normal uv
    uv[0] = tx0;
    uv[1] = ty;
    uv[2] = tx;
    uv[3] = ty;
    uv[4] = tx;
    uv[5] = ty0;
    uv[6] = tx0;
    uv[7] = ty0;
 
    // pixels
    coords[0] = 0;
    coords[1] = 0;
    coords[2] = sx;
    coords[3] = 0;
    coords[4] = sx;
    coords[5] = sy;
    coords[6] = 0;
    coords[7] = sy;
 
    wxColour color = GetGlobalColor("NODTA");
    glClearColor(color.Red() / 256., color.Green() / 256., color.Blue() / 256.,
                 1.0);
    glClear(GL_COLOR_BUFFER_BIT);
 
    RenderTextures(gldc, coords, uv, 4, m_pParentCanvas->GetpVP());
#endif
 
    glDisable(g_texture_rectangle_format);
 
    m_cache_vp = VPoint;
    m_cache_vp.Validate();
 
    m_cache_current_ch = m_pParentCanvas->m_singleChart;
 
    if (VPoint.b_quilt) m_pParentCanvas->m_pQuilt->SetRenderedVP(VPoint);
 
  } else  // useFBO
  {
    RenderCharts(m_gldc, screen_region);
  }
 
  // if (m_binPinch)
  //   printf("        Render Charts Done  %ld\n",
  //          m_glstopwatch.Time() - render_start_time);
 
  // Done with base charts.
  // Now the overlays
  RenderS57TextOverlay(VPoint);
  RenderMBTilesOverlay(VPoint);
 
  g_overlayCanvas = m_pParentCanvas;
  if (g_pi_manager) {
    g_pi_manager->SendViewPortToRequestingPlugIns(VPoint);
    g_pi_manager->RenderAllGLCanvasOverlayPlugIns(
        m_pcontext, VPoint, m_pParentCanvas->m_canvasIndex, OVERLAY_CHARTS);
  }
 
  // Render static overlay objects
  for (OCPNRegionIterator upd(screen_region); upd.HaveRects(); upd.NextRect()) {
    wxRect rt = upd.GetRect();
    LLRegion region = VPoint.GetLLRegion(rt);
    ViewPort cvp = ClippedViewport(VPoint, region);
    DrawGroundedOverlayObjects(gldc, cvp);
  }
 
  if (m_pParentCanvas->m_bShowTide || m_pParentCanvas->m_bShowCurrent) {
    LLRegion screenLLRegion = VPoint.GetLLRegion(screen_region);
    LLBBox screenBox = screenLLRegion.GetBox();
    // Enlarge the box a bit
    screenBox.EnLarge(screenBox.GetLonRange() * 0.05);
 
    // update the tide/current select points, if necessary
    if (m_pParentCanvas->m_bShowTide) {
      m_pParentCanvas->RebuildTideSelectList(screenBox);  // full screen
      DrawGLTidesInBBox(gldc, VPoint.GetBBox());
    }
 
    if (m_pParentCanvas->m_bShowCurrent) {
      m_pParentCanvas->RebuildCurrentSelectList(screenBox);
      DrawGLCurrentsInBBox(gldc, VPoint.GetBBox());
    }
  }
 
  // If multi-canvas, indicate which canvas has keyboard focus
  // by drawing a simple blue bar at the top.
  if (m_pParentCanvas->m_show_focus_bar &&
      (g_canvasConfig != 0)) {  // multi-canvas?
    if (m_pParentCanvas == wxWindow::FindFocus()) {
      g_focusCanvas = m_pParentCanvas;
 
      wxColour colour = GetGlobalColor("BLUE4");
      wxPen ppBlue(colour, 1);
      wxBrush ppBrush(colour);
      gldc.SetPen(ppBlue);
      gldc.SetBrush(ppBrush);
      int xw = m_pParentCanvas->GetClientSize().x * m_displayScale;
      float rect_pix = m_pParentCanvas->m_focus_indicator_pix * m_displayScale;
      wxPoint barPoints[4];
      barPoints[0].x = 0;
      barPoints[0].y = 0;
      barPoints[1].x = xw;
      barPoints[1].y = 0;
      barPoints[2].x = xw;
      barPoints[2].y = rect_pix;
      barPoints[3].x = 0;
      barPoints[3].y = rect_pix;
 
      gldc.DrawPolygon(4, barPoints, 0, 0, 1, 0);
    }
  }
 
  DrawDynamicRoutesTracksAndWaypoints(VPoint);
 
  // Now draw all the objects which normally move around and are not
  // cached from the previous frame
  DrawFloatingOverlayObjects(m_gldc);
 
#ifndef USE_ANDROID_GLES2
  // from this point on don't use perspective
  if (VPoint.tilt) {
    glMatrixMode(GL_PROJECTION);
    glLoadIdentity();
 
    glOrtho(0, (GLint)gl_width, (GLint)gl_height, 0, -1, 1);
    glMatrixMode(GL_MODELVIEW);
    glLoadIdentity();
  }
#endif
 
  if (!g_bhide_depth_units)
    DrawEmboss(m_gldc, m_pParentCanvas->EmbossDepthScale());
  if (!g_bhide_overzoom_flag)
    DrawEmboss(m_gldc, m_pParentCanvas->EmbossOverzoomIndicator(gldc));
 
  if (g_pi_manager) {
    ViewPort &vp = m_pParentCanvas->GetVP();
    g_pi_manager->SendViewPortToRequestingPlugIns(vp);
    g_pi_manager->RenderAllGLCanvasOverlayPlugIns(
        m_pcontext, vp, m_pParentCanvas->m_canvasIndex, OVERLAY_OVER_EMBOSS);
  }
  if (!g_PrintingInProgress) {
    if (m_pParentCanvas->m_pTrackRolloverWin)
      m_pParentCanvas->m_pTrackRolloverWin->Draw(gldc);
 
    if (m_pParentCanvas->m_pRouteRolloverWin)
      m_pParentCanvas->m_pRouteRolloverWin->Draw(gldc);
 
    if (m_pParentCanvas->m_pAISRolloverWin)
      m_pParentCanvas->m_pAISRolloverWin->Draw(gldc);
 
    if (m_pParentCanvas->GetMUIBar())
      m_pParentCanvas->GetMUIBar()->DrawGL(gldc, m_displayScale);
 
    if (g_MainToolbar && m_pParentCanvas->IsPrimaryCanvas())
      g_MainToolbar->DrawGL(gldc, m_displayScale);
 
    if (g_iENCToolbar && m_pParentCanvas->IsPrimaryCanvas())
      g_iENCToolbar->DrawGL(gldc, m_displayScale);
  }
 
  //  On some platforms, the opengl context window is always on top of any
  //  standard DC windows, so we need to draw the Chart Info Window
  //  as overlayed bmps.
 
#ifdef __WXOSX__
  if (m_pParentCanvas->m_pCIWin && m_pParentCanvas->m_pCIWin->IsShown()) {
    int x, y, width, height;
    m_pParentCanvas->m_pCIWin->GetClientSize(&width, &height);
    m_pParentCanvas->m_pCIWin->GetPosition(&x, &y);
    wxBitmap bmp(width, height, -1);
    wxMemoryDC dc(bmp);
    if (bmp.IsOk()) {
      dc.SetBackground(wxBrush(GetGlobalColor("UIBCK")));
      dc.Clear();
 
      dc.SetTextBackground(GetGlobalColor("UIBCK"));
      dc.SetTextForeground(GetGlobalColor("UITX1"));
 
      int yt = 0;
      int xt = 0;
      wxString s = m_pParentCanvas->m_pCIWin->GetString();
      int h = m_pParentCanvas->m_pCIWin->GetCharHeight();
 
      wxStringTokenizer tkz(s, "\n");
      wxString token;
 
      while (tkz.HasMoreTokens()) {
        token = tkz.GetNextToken();
        dc.DrawText(token, xt, yt);
        yt += h;
      }
      dc.SelectObject(wxNullBitmap);
 
      m_gldc.DrawBitmap(bmp, x, y, false);
    }
  }
 
#endif
  // render the chart bar
  if (g_bShowChartBar) DrawChartBar(m_gldc);
 
  if (m_pParentCanvas->m_Compass && m_pParentCanvas->m_bShowCompassWin &&
      g_bShowCompassWin)
    m_pParentCanvas->m_Compass->Paint(gldc);
 
  if (m_pParentCanvas->IsPrimaryCanvas() &&
      m_pParentCanvas->m_notification_button) {
    auto &noteman = NotificationManager::GetInstance();
    if (noteman.GetNotificationCount()) {
      m_pParentCanvas->m_notification_button->SetIconSeverity(
          noteman.GetMaxSeverity());
      if (m_pParentCanvas->m_notification_button->UpdateStatus()) Refresh();
      m_pParentCanvas->m_notification_button->Show(true);
      m_pParentCanvas->m_notification_button->Paint(gldc);
    } else {
      m_pParentCanvas->m_notification_button->Show(false);
    }
  }
  RenderGLAlertMessage();
 
  if (g_pi_manager) {
    ViewPort &vp = m_pParentCanvas->GetVP();
    g_pi_manager->SendViewPortToRequestingPlugIns(vp);
    g_pi_manager->RenderAllGLCanvasOverlayPlugIns(
        m_pcontext, vp, m_pParentCanvas->m_canvasIndex, OVERLAY_OVER_UI);
    glActiveTexture(GL_TEXTURE0);
  }
 
  // quiting?
  if (g_bquiting) DrawQuiting();
  if (g_bcompression_wait)
    DrawCloseMessage(_("Waiting for raster chart compression thread exit."));
 
  //  Some older MSW OpenGL drivers are generally very unstable.
  //  This helps...
 
  SwapBuffers();
 
  g_glTextureManager->TextureCrunch(0.8);
  g_glTextureManager->FactoryCrunch(0.6);
 
  m_pParentCanvas->PaintCleanup();
  m_bforcefull = false;
 
  // if (m_binPinch)
  //   printf("    Render Finished:  %ld\n",
  //          m_glstopwatch.Time() - render_start_time);
 
  n_render++;
}
 
void glChartCanvas::RenderS57TextOverlay(ViewPort &VPoint) {
  //  Render the decluttered Text overlay for quilted vector charts, except for
  //  CM93 Composite
  if (VPoint.b_quilt) {
    if (m_pParentCanvas->m_pQuilt->IsQuiltVector() && ps52plib &&
        ps52plib->GetShowS57Text()) {
      ChartBase *chart = m_pParentCanvas->m_pQuilt->GetRefChart();
      if (chart && (chart->GetChartType() != CHART_TYPE_CM93COMP)) {
        //        Clear the text Global declutter list
        if (chart) {
          ChartPlugInWrapper *ChPI = dynamic_cast<ChartPlugInWrapper *>(chart);
          if (ChPI)
            ChPI->ClearPLIBTextList();
          else
            ps52plib->ClearTextList();
        }
 
        // Grow the ViewPort a bit laterally, to minimize "jumping" of text
        // elements at right side of screen
        ViewPort vpx = VPoint;
        vpx.BuildExpandedVP(VPoint.pix_width * 12 / 10, VPoint.pix_height);
 
        OCPNRegion screen_region(
            wxRect(0, 0, VPoint.pix_width, VPoint.pix_height));
        RenderQuiltViewGLText(vpx, screen_region);
      }
    }
  }
}
void glChartCanvas::RenderSingleMBTileOverlay(const int dbIndex, bool bOverlay,
                                              ViewPort &vp,
                                              OCPNRegion &screen_region,
                                              LLRegion &screenLLRegion) {
  ChartBase *chart = ChartData->OpenChartFromDBAndLock(dbIndex, FULL_INIT);
 
  // Chart may have been prevented from initial loading due to size, or some
  // other reason...
  if (chart == NULL) return;
 
  ChartMbTiles *pcmbt = dynamic_cast<ChartMbTiles *>(chart);
  if (!pcmbt) return;
 
  // Is tile an OVERLAY type?
  // Render, or not, depending on passed flag.
  if (bOverlay && pcmbt->GetTileType() != MbTilesType::OVERLAY) return;
 
  wxFileName tileFile(chart->GetFullPath());
  // Size test for 5 GByte
  wxULongLong tileSizeMB = tileFile.GetSize() >> 20;
 
  // Auto-show MBTiles in basemap directories
  bool isBasemap = tileFile.GetPath().Lower().Contains("basemap");
 
  // For basemap MBTiles, stop rendering when zoomed in past their max detail
  // to allow S57 vector charts to show through
  if (isBasemap) {
    double chart_native_ppm = m_pParentCanvas->GetCanvasScaleFactor() /
                              (double)chart->GetNativeScale();
    double zoom_ratio = vp.view_scale_ppm / chart_native_ppm;
    if (zoom_ratio > g_tile_basemap_zoom_factor) return;
  }
 
  if (!isBasemap && (!ChartData->CheckAnyCanvasExclusiveTileGroup() ||
                     (tileSizeMB.GetLo() > 5000))) {
    // Check to see if the tile has been "clicked".
    // If so, do not add to no-show array again.
    if (!m_pParentCanvas->IsTileOverlayIndexInYesShow(dbIndex)) {
      if (!m_pParentCanvas->IsTileOverlayIndexInNoShow(dbIndex)) {
        m_pParentCanvas->m_tile_noshow_index_array.push_back(dbIndex);
      }
    }
  }
 
  // This test catches the case where the chart is added to no_show list
  // when first loaded by OpenChartFromDBAndLock
  if (m_pParentCanvas->IsTileOverlayIndexInNoShow(dbIndex)) {
    return;
  }
 
  pcmbt->RenderRegionViewOnGL(*m_pcontext, vp, screen_region, screenLLRegion);
 
  // Light up the piano key if the chart was rendered
  std::vector<int> piano_active_array_tiles =
      m_pParentCanvas->m_Piano->GetActiveKeyArray();
  bool bfound = false;
 
  if (std::find(piano_active_array_tiles.begin(),
                piano_active_array_tiles.end(),
                dbIndex) != piano_active_array_tiles.end()) {
    bfound = true;
  }
 
  if (!bfound) {
    piano_active_array_tiles.push_back(dbIndex);
    m_pParentCanvas->m_Piano->SetActiveKeyArray(piano_active_array_tiles);
  }
}
 
void glChartCanvas::RenderMBTilesOverlay(ViewPort &VPoint) {
  // Render MBTiles as overlay
  std::vector<int> stackIndexArray =
      m_pParentCanvas->m_pQuilt->GetExtendedStackIndexArray();
  unsigned int im = stackIndexArray.size();
  // XXX should
  // assert(!VPoint.b_quilt && im == 0)
  if (VPoint.b_quilt && im > 0) {
    bool regionVPBuilt = false;
    OCPNRegion screen_region;
    LLRegion screenLLRegion;
    LLBBox screenBox;
    ViewPort vp;
 
    std::vector<int> tiles_to_show;
    for (unsigned int is = 0; is < im; is++) {
      const ChartTableEntry &cte =
          ChartData->GetChartTableEntry(stackIndexArray[is]);
      if (cte.GetChartType() == CHART_TYPE_MBTILES) {
        if (m_pParentCanvas->IsTileOverlayIndexInNoShow(stackIndexArray[is])) {
          // Turn off the piano highlite
          std::vector<int> piano_active_array_tiles =
              m_pParentCanvas->m_Piano->GetActiveKeyArray();
          bool bfound = false;
 
          for (unsigned int i = 0; i < piano_active_array_tiles.size(); i++) {
            if (piano_active_array_tiles[i] == stackIndexArray[is]) {
              piano_active_array_tiles.erase(piano_active_array_tiles.begin() +
                                             i);  // erase it
              bfound = true;
              break;
            }
          }
 
          if (bfound)
            m_pParentCanvas->m_Piano->SetActiveKeyArray(
                piano_active_array_tiles);
 
          continue;
        }
 
        tiles_to_show.push_back(stackIndexArray[is]);
        if (!regionVPBuilt) {
          screen_region =
              OCPNRegion(wxRect(0, 0, VPoint.pix_width, VPoint.pix_height));
          screenLLRegion = VPoint.GetLLRegion(screen_region);
          screenBox = screenLLRegion.GetBox();
 
          vp = VPoint;
          wxPoint p;
          p.x = VPoint.pix_width / 2;
          p.y = VPoint.pix_height / 2;
          VPoint.GetLLFromPix(p, &vp.clat, &vp.clon);
 
          regionVPBuilt = true;
        }
      }
    }
 
    // Render in two passes, to render the OVERLAY types last
 
    // Show the tilesets in reverse order to have the largest scale
    // on top
 
    for (std::vector<int>::reverse_iterator rit = tiles_to_show.rbegin();
         rit != tiles_to_show.rend(); ++rit) {
      RenderSingleMBTileOverlay(*rit, FALSE, vp, screen_region, screenLLRegion);
    }
    for (std::vector<int>::reverse_iterator rit = tiles_to_show.rbegin();
         rit != tiles_to_show.rend(); ++rit) {
      RenderSingleMBTileOverlay(*rit, TRUE, vp, screen_region, screenLLRegion);
    }
 
    // Render the HiLite on piano rollover of mbTile key
    LLRegion hiregion = m_pParentCanvas->m_pQuilt->GetHiliteRegion();
 
    if (!hiregion.Empty()) {
      glEnable(GL_BLEND);
 
      double hitrans;
      switch (global_color_scheme) {
        case GLOBAL_COLOR_SCHEME_DAY:
          hitrans = .4;
          break;
        case GLOBAL_COLOR_SCHEME_DUSK:
          hitrans = .2;
          break;
        case GLOBAL_COLOR_SCHEME_NIGHT:
          hitrans = .1;
          break;
        default:
          hitrans = .4;
          break;
      }
 
#if !defined(USE_ANDROID_GLES2) && !defined(ocpnUSE_GLSL)
      glColor4f((float).8, (float).4, (float).4, (float)hitrans);
#else
      s_regionColor = wxColor(204, 102, 102, hitrans * 256);
#endif
 
      DrawRegion(VPoint, hiregion);
 
      glDisable(GL_BLEND);
    }
  }
}
 
#if 0
void glChartCanvas::RenderCanvasBackingChart(ocpnDC &dc,
                                             OCPNRegion valid_region) {
  //  Fill the FBO with the current gshhs world chart
  int w, h;
  GetClientSize(&w, &h);
 
  glViewport(0, 0, (GLint)m_cache_tex_x, (GLint)m_cache_tex_y);
#if !defined(USE_ANDROID_GLES2) && !defined(ocpnUSE_GLSL)
  glMatrixMode(GL_PROJECTION);
  glLoadIdentity();
 
  glOrtho(0, m_cache_tex_x, m_cache_tex_y, 0, -1, 1);
  glMatrixMode(GL_MODELVIEW);
  glLoadIdentity();
#endif
 
  wxRect rtex(0, 0, m_cache_tex_x, m_cache_tex_y);
  ViewPort cvp =
      m_pParentCanvas->GetVP().BuildExpandedVP(m_cache_tex_x, m_cache_tex_y);
 
  bool world_view = false;
  RenderWorldChart(dc, cvp, rtex, world_view);
  gShapeBasemap.RenderViewOnDC(dc, cvp);
 
  //    dc.SetPen(wxPen(wxColour(254,254,0), 3));
  //    dc.DrawLine( 0, 0, m_cache_tex_x, m_cache_tex_y);
 
  //  Reset matrices
  glViewport(0, 0, (GLint)w, (GLint)h);
#if !defined(USE_ANDROID_GLES2) && !defined(ocpnUSE_GLSL)
  glMatrixMode(GL_PROJECTION);
  glLoadIdentity();
 
  glOrtho(0, (GLint)w, (GLint)h, 0, -1, 1);
  glMatrixMode(GL_MODELVIEW);
  glLoadIdentity();
#endif
}
#endif
 
void glChartCanvas::FastPan(int dx, int dy) {
#if !defined(USE_ANDROID_GLES2) && !defined(ocpnUSE_GLSL)
#endif
}
 
void glChartCanvas::ZoomProject(float offset_x, float offset_y, float swidth,
                                float sheight) {
  if (IsShown()) SetCurrent(*m_pcontext);
  float sx = GetSize().x;
  float sy = GetSize().y;
  glClear(GL_COLOR_BUFFER_BIT);
 
  int w, h;
  GetClientSize(&w, &h);
 
  if (s_b_useStencil) {
    glEnable(GL_STENCIL_TEST);
    glStencilMask(0xff);
    glClear(GL_STENCIL_BUFFER_BIT);
    glDisable(GL_STENCIL_TEST);
  }
 
#ifndef __ANDROID__
  // Render backing texture
  if (1) {
    float tx0 = 0;
    float ty0 = 0;
    float tx = sx;
    float ty = sy;
 
    float vx0 = 0;
    float vy0 = 0;
    float vx = sx;
    float vy = sy;
 
    float sxfactor = sx / swidth;
    float syfactor = sy / sheight;
 
    glViewport(-offset_x * sx / swidth - (sx * sxfactor / 2),
               -offset_y * (sy / sheight) - (sy * syfactor / 2),
               sx * sx / swidth * 2, sy * sy / sheight * 2);
    glBindTexture(g_texture_rectangle_format, m_TouchBackingTexture);
    glEnable(g_texture_rectangle_format);
 
    float uv[8];
    float coords[8];
 
    // normal uv Texture coordinates
    uv[0] = 0;
    uv[1] = 0;
    uv[2] = 1;
    uv[3] = 0;
    uv[4] = 1;
    uv[5] = 1;
    uv[6] = 0;
    uv[7] = 1;
 
    // pixels
    coords[0] = vx0;
    coords[1] = vy0;
    coords[2] = vx;
    coords[3] = vy0;
    coords[4] = vx;
    coords[5] = vy;
    coords[6] = vx0;
    coords[7] = vy;
 
    RenderTextures(m_gldc, coords, uv, 4, &m_texVP);
    glBindTexture(g_texture_rectangle_format, 0);
  }
#endif
 
  // render zoomed canvas section
  if (1) {
    float tx, ty, tx0, ty0;
    tx = sx, ty = sy;
 
    tx0 = ty0 = 0.;
 
    float vx0 = 0;
    float vy0 = 0;
    float vy = sy;
    float vx = sx;
 
    glBindTexture(g_texture_rectangle_format, 0);
 
    // Render the cached texture as quad to screen
    glBindTexture(g_texture_rectangle_format, m_cache_tex[m_cache_page]);
    glEnable(g_texture_rectangle_format);
 
    float uv[8];
    float coords[8];
 
    // normal uv Texture coordinates
    uv[0] = tx0 / m_cache_tex_x;
    uv[1] = ty / m_cache_tex_y;
    uv[2] = tx / m_cache_tex_x;
    uv[3] = ty / m_cache_tex_y;
    uv[4] = tx / m_cache_tex_x;
    uv[5] = ty0 / m_cache_tex_y;
    uv[6] = tx0 / m_cache_tex_x;
    uv[7] = ty0 / m_cache_tex_y;
 
    // pixels
    coords[0] = vx0;
    coords[1] = vy0;
    coords[2] = vx;
    coords[3] = vy0;
    coords[4] = vx;
    coords[5] = vy;
    coords[6] = vx0;
    coords[7] = vy;
 
    glViewport(-offset_x * sx / swidth, -offset_y * (sy / sheight),
               sx * sx / swidth, sy * sy / sheight);
 
    RenderTextures(m_gldc, coords, uv, 4, m_pParentCanvas->GetpVP());
 
    glDisable(g_texture_rectangle_format);
    glBindTexture(g_texture_rectangle_format, 0);
  }
 
#if 0
  // For fun, we prove the coordinates of the blank area outside the chart when
  // zooming out. Bottom stripe
  // wxColour color = GetGlobalColor("YELO1");   //GREY1
  wxColour color = GetGlobalColor("GREY1");  //
  float ht = -offset_y * (sy / sheight);
  wxRect r(0, sy - ht, w, ht);
  RenderColorRect(r, color);
 
  // top stripe
  wxRect rt(0, 0, w, sy - (ht + (sy * sy / sheight)));
  RenderColorRect(rt, color);
 
  // left
  float w1 = -offset_x * sx / swidth;
  wxRect rl(0, 0, w1, sy);
  RenderColorRect(rl, color);
 
  // right
  float px = w1 + sx * sx / swidth;
  wxRect rr(px, 0, sx - px, sy);
  RenderColorRect(rr, color);
#endif
  //  When zooming out, if we go too far, then the frame buffer is repeated
  //  on-screen due to address wrapping in the frame buffer. Detect this case,
  //  and render some simple solid covering quads to avoid a confusing display.
 
  SwapBuffers();
}
 
void glChartCanvas::onZoomTimerEvent(wxTimerEvent &event) {
  // If m_zoomFinal is set, stop the timer.
  if (!m_zoomFinal) {
    if (m_nRun < m_nTotal) {
      m_runoffsetx += m_offsetxStep;
      if (m_offsetxStep > 0)
        m_runoffsetx = wxMin(m_runoffsetx, m_fbo_offsetx);
      else
        m_runoffsetx = wxMax(m_runoffsetx, m_fbo_offsetx);
 
      m_runoffsety += m_offsetyStep;
      if (m_offsetyStep > 0)
        m_runoffsety = wxMin(m_runoffsety, m_fbo_offsety);
      else
        m_runoffsety = wxMax(m_runoffsety, m_fbo_offsety);
 
      m_runswidth += m_swidthStep;
      if (m_swidthStep > 0)
        m_runswidth = wxMin(m_runswidth, m_fbo_swidth);
      else
        m_runswidth = wxMax(m_runswidth, m_fbo_swidth);
 
      m_runsheight += m_sheightStep;
      if (m_sheightStep > 0)
        m_runsheight = wxMin(m_runsheight, m_fbo_sheight);
      else
        m_runsheight = wxMax(m_runsheight, m_fbo_sheight);
 
      m_nRun += m_nStep;
    }
 
    ZoomProject(m_runoffsetx, m_runoffsety, m_runswidth, m_runsheight);
 
  } else {
    zoomTimer.Stop();
    if (m_zoomFinal) {
      m_pParentCanvas->ZoomCanvasSimple(m_zoomFinalZoom);
      if (m_zoomFinaldx || m_zoomFinaldy) {
        m_pParentCanvas->PanCanvas(m_zoomFinaldx, m_zoomFinaldy);
      }
    }
    m_zoomFinal = false;
  }
}
 
void glChartCanvas::FastZoom(float factor, float cp_x, float cp_y, float post_x,
                             float post_y) {
  int sx = GetSize().x;
  int sy = GetSize().y;
 
  m_lastfbo_offsetx = m_fbo_offsetx;
  m_lastfbo_offsety = m_fbo_offsety;
  m_lastfbo_swidth = m_fbo_swidth;
  m_lastfbo_sheight = m_fbo_sheight;
 
  float curr_fbo_offset_x = m_fbo_offsetx;
  float curr_fbo_offset_y = m_fbo_offsety;
  float curr_fbo_swidth = m_fbo_swidth;
  float curr_fbo_sheight = m_fbo_sheight;
 
  float fx = (float)cp_x / sx;
  float fy = 1.0 - (float)cp_y / sy;
 
  float fbo_ctr_x = curr_fbo_offset_x + (curr_fbo_swidth * fx);
  float fbo_ctr_y = curr_fbo_offset_y + (curr_fbo_sheight * fy);
 
  m_fbo_swidth = curr_fbo_swidth / factor;
  m_fbo_sheight = curr_fbo_sheight / factor;
 
  m_fbo_offsetx = fbo_ctr_x - (m_fbo_swidth * fx);
  m_fbo_offsety = fbo_ctr_y - (m_fbo_sheight * fy);
 
  m_fbo_offsetx += post_x;
  m_fbo_offsety += post_y;
 
  {
    m_nStep = 20;
    m_nTotal = 100;
 
    // m_nStep = 10;    // Android?
    // m_nTotal = 40;
 
    m_nRun = 0;
 
    float perStep = m_nStep / m_nTotal;
 
    if (zoomTimer.IsRunning()) {
      m_offsetxStep = (m_fbo_offsetx - m_runoffsetx) * perStep;
      m_offsetyStep = (m_fbo_offsety - m_runoffsety) * perStep;
      m_swidthStep = (m_fbo_swidth - m_runswidth) * perStep;
      m_sheightStep = (m_fbo_sheight - m_runsheight) * perStep;
 
    } else {
      m_offsetxStep = (m_fbo_offsetx - m_lastfbo_offsetx) * perStep;
      m_offsetyStep = (m_fbo_offsety - m_lastfbo_offsety) * perStep;
      m_swidthStep = (m_fbo_swidth - m_lastfbo_swidth) * perStep;
      m_sheightStep = (m_fbo_sheight - m_lastfbo_sheight) * perStep;
 
      m_runoffsetx = m_lastfbo_offsetx;
      m_runoffsety = m_lastfbo_offsety;
      m_runswidth = m_lastfbo_swidth;
      m_runsheight = m_lastfbo_sheight;
    }
 
    if (!zoomTimer.IsRunning()) zoomTimer.Start(m_nStep);
    m_zoomFinal = false;
  }
}
 
#ifdef __ANDROID__
 
void glChartCanvas::OnEvtPanGesture(wxQT_PanGestureEvent &event) {
  // qDebug() << "OnEvtPanGesture" << m_pParentCanvas->m_canvasIndex <<
  // event.cursor_pos.x;
 
  if (m_pParentCanvas->isRouteEditing() || m_pParentCanvas->isMarkEditing())
    return;
 
  if (m_binPinch) return;
  if (m_bpinchGuard) return;
 
  int x = event.GetOffset().x;
  int y = event.GetOffset().y;
 
  int lx = event.GetLastOffset().x;
  int ly = event.GetLastOffset().y;
 
  int dx = lx - x;
  int dy = y - ly;
 
  switch (event.GetState()) {
    case GestureStarted:
      if (m_binPan) break;
 
      panx = pany = 0;
      m_binPan = true;
      m_binGesture = true;
      // qDebug() << "pg1";
      break;
 
    case GestureUpdated:
      if (m_binPan) {
        if (!g_GLOptions.m_bUseCanvasPanning) {
          // qDebug() << "slowpan" << dx << dy;
 
          m_pParentCanvas->FreezePiano();
          m_pParentCanvas->PanCanvas(dx, -dy);
          m_pParentCanvas->ThawPiano();
 
        } else {
          FastPan(dx, dy);
        }
 
        panx -= dx;
        pany -= dy;
      }
      break;
 
    case GestureFinished:
      // qDebug() << "panGestureFinished";
 
      m_pParentCanvas->UpdateCanvasControlBar();
 
      m_binPan = false;
      m_gestureFinishTimer.Start(500, wxTIMER_ONE_SHOT);
 
      break;
 
    case GestureCanceled:
      m_binPan = false;
      m_gestureFinishTimer.Start(500, wxTIMER_ONE_SHOT);
      break;
 
    default:
      break;
  }
 
  m_bgestureGuard = true;
  m_gestureEeventTimer.Start(500, wxTIMER_ONE_SHOT);
  m_bforcefull = false;
 
  // qDebug() << "panGestureDone";
}
 
float zoom_inc = 1.0;
 
void glChartCanvas::OnEvtPinchGesture(wxQT_PinchGestureEvent &event) {
  float zoom_gain = 1.0;
  float zout_gain = 1.0;
 
  float zoom_val;
  float total_zoom_val;
 
  float max_zoom_scale = 1000.;
  float min_zoom_scale = 2e8;
 
  if (event.GetScaleFactor() > 1)
    zoom_val = ((event.GetScaleFactor() - 1.0) * zoom_gain) + 1.0;
  else
    zoom_val = 1.0 - ((1.0 - event.GetScaleFactor()) * zout_gain);
 
  if (event.GetTotalScaleFactor() > 1)
    total_zoom_val = ((event.GetTotalScaleFactor() - 1.0) * zoom_gain) + 1.0;
  else
#if 0
        total_zoom_val = 1.0 - ((1.0 - event.GetTotalScaleFactor()) * zout_gain);
 
    double projected_scale = cc1->GetVP().chart_scale / total_zoom_val;
 
    // Max zoom in is set by scale of quilt reference chart, consistent with chart render limits set elsewhere.
    float max_zoom_scale = 1000.;
    if( cc1->GetVP().b_quilt) {
        int ref_index = cc1->GetQuiltRefChartdbIndex();
//         if((ref_index >= 0) && ChartData){
//             max_zoom_scale = ChartData->GetDBChartScale(ref_index) / 8.0;
//         }
    }
 
 
    float min_zoom_scale = 2e8;
 
#endif
 
    total_zoom_val = 1.0 - ((1.0 - event.GetTotalScaleFactor()) * zoom_gain);
 
  double projected_scale =
      m_pParentCanvas->GetVP().chart_scale / total_zoom_val;
 
  switch (event.GetState()) {
    case GestureStarted:
      m_first_zout = false;
      m_binPinch = true;
      m_binPan = false;  // cancel any tentative pan gesture, in case the "pan
                         // cancel" event was lost
      m_binGesture = true;
      // qDebug() << "pg2";
      m_pinchStart = event.GetCenterPoint();
      m_lpinchPoint = m_pinchStart;
 
      m_pParentCanvas->GetCanvasPixPoint(event.GetCenterPoint().x,
                                         event.GetCenterPoint().y, m_pinchlat,
                                         m_pinchlon);
      //            qDebug() << "center" << event.GetCenterPoint().x <<
      //            event.GetCenterPoint().y;
 
      m_cc_x = m_fbo_offsetx + (m_fbo_swidth / 2);
      m_cc_y = m_fbo_offsety + (m_fbo_sheight / 2);
 
      // Render the full charts with overlay objects onto the frame buffer.
      if (IsShown()) SetCurrent(*m_pcontext);
      RenderScene();
 
      zoom_inc = 1.0;
      break;
 
    case GestureUpdated:
      if (g_GLOptions.m_bUseCanvasPanning) {
        if (projected_scale < min_zoom_scale) {
          wxPoint pinchPoint = event.GetCenterPoint();
 
          float dx = pinchPoint.x - m_lpinchPoint.x;
          float dy = pinchPoint.y - m_lpinchPoint.y;
 
          FastZoom(zoom_val, m_pinchStart.x, m_pinchStart.y,
                   -dx / total_zoom_val, dy / total_zoom_val);
 
          m_lpinchPoint = pinchPoint;
        }
      } else {
        // qDebug() << "update totalzoom" << total_zoom_val << projected_scale;
        if (1 || ((total_zoom_val > 1) && !m_first_zout)) {  // Zoom in
          wxPoint pinchPoint = event.GetCenterPoint();
 
          float dx = pinchPoint.x - m_lpinchPoint.x;
          float dy = pinchPoint.y - m_lpinchPoint.y;
 
          if ((projected_scale > max_zoom_scale) &&
              (projected_scale < min_zoom_scale))
            FastZoom(zoom_val, m_pinchStart.x, m_pinchStart.y,
                     -dx / total_zoom_val, dy / total_zoom_val);
 
          m_lpinchPoint = pinchPoint;
 
        } else {
          m_first_zout = true;
          zoom_inc *= zoom_val;
          if ((zoom_inc < 0.9) || (zoom_inc > 1.1)) {
            m_pParentCanvas->ZoomCanvas(zoom_inc, false);
            zoom_inc = 1.0;
          }
 
          wxPoint pinchPoint = event.GetCenterPoint();
          float dx = pinchPoint.x - m_lpinchPoint.x;
          float dy = pinchPoint.y - m_lpinchPoint.y;
          m_pParentCanvas->PanCanvas(-dx, -dy);
          m_lpinchPoint = pinchPoint;
 
          //                         SetCurrent(*m_pcontext);
          //                         RenderScene();
          //                         g_Piano->DrawGL(cc1->m_canvas_height -
          //                         g_Piano->GetHeight()); SwapBuffers();
        }
      }
 
      break;
 
    case GestureFinished: {
      //            qDebug() << "finish totalzoom" << total_zoom_val <<
      //            projected_scale;
 
      float cc_x = m_fbo_offsetx + (m_fbo_swidth / 2);
      float cc_y = m_fbo_offsety + (m_fbo_sheight / 2);
      float dy = 0;
      float dx = 0;
 
      float tzoom = total_zoom_val;
 
      if (projected_scale >= min_zoom_scale)
        tzoom = m_pParentCanvas->GetVP().chart_scale / min_zoom_scale;
 
      if (projected_scale < max_zoom_scale)
        tzoom = m_pParentCanvas->GetVP().chart_scale / max_zoom_scale;
 
      dx = (cc_x - m_cc_x) * tzoom;
      dy = -(cc_y - m_cc_y) * tzoom;
 
      if (zoomTimer.IsRunning()) {
        //                qDebug() << "Final zoom";
        m_zoomFinal = true;
        m_zoomFinalZoom = tzoom;
        m_zoomFinaldx = dx;
        m_zoomFinaldy = dy;
      }
 
      else {
        double final_projected_scale =
            m_pParentCanvas->GetVP().chart_scale / tzoom;
        // qDebug() << "Final pinchB" << tzoom << final_projected_scale;
 
        if (final_projected_scale < min_zoom_scale) {
          // qDebug() << "zoomit";
          m_pParentCanvas->ZoomCanvas(tzoom, false);
          m_pParentCanvas->PanCanvas(dx, dy);
          m_pParentCanvas->m_pQuilt->Invalidate();
          m_bforcefull = true;
 
        } else {
          double new_scale =
              m_pParentCanvas->GetCanvasScaleFactor() / min_zoom_scale;
          // qDebug() << "clampit";
          m_pParentCanvas->SetVPScale(new_scale);
          m_pParentCanvas->m_pQuilt->Invalidate();
          m_bforcefull = true;
        }
      }
 
      //                if( projected_scale < 3e8)
      //                    m_pParentCanvas->ZoomCanvas( total_zoom_val, false
      //                    );
      //                else
      //                    m_pParentCanvas->ZoomCanvas(m_pParentCanvas->GetVP().chart_scale
      //                    / 3e8, false);
 
      m_binPinch = false;
      m_gestureFinishTimer.Start(500, wxTIMER_ONE_SHOT);
      break;
    }
 
    case GestureCanceled:
      m_binPinch = false;
      m_gestureFinishTimer.Start(500, wxTIMER_ONE_SHOT);
      break;
 
    default:
      break;
  }
 
  m_bgestureGuard = true;
  //    m_bpinchGuard = true;
  m_gestureEeventTimer.Start(500, wxTIMER_ONE_SHOT);
}
 
void glChartCanvas::onGestureTimerEvent(wxTimerEvent &event) {
  //  On some devices, the pan GestureFinished event fails to show up
  //  Watch for this case, and fix it.....
  // qDebug() << "onGestureTimerEvent";
 
  if (m_binPan) {
    m_binPan = false;
    Invalidate();
    Update();
  }
  m_bgestureGuard = false;
  m_bpinchGuard = false;
  m_binGesture = false;
  m_bforcefull = false;
}
 
void glChartCanvas::onGestureFinishTimerEvent(wxTimerEvent &event) {
  // qDebug() << "onGestureFinishTimerEvent";
 
  // signal gesture is finished after a delay
  m_binGesture = false;
  m_bforcefull = false;
}
 
#else
#ifdef HAVE_WX_GESTURE_EVENTS
 
void glChartCanvas::OnEvtPanGesture(wxPanGestureEvent &event) {
  // qDebug() << "OnEvtPanGesture" << m_pParentCanvas->m_canvasIndex <<
  // event.cursor_pos.x;
  if (m_pParentCanvas->isRouteEditing() || m_pParentCanvas->isMarkEditing())
    return;
 
  if (m_binPinch) return;
  if (m_bpinchGuard) return;
 
  int dx = event.GetDelta().x;
  int dy = event.GetDelta().y;
 
  if (event.IsGestureStart()) {
    if (m_binPan) return;
 
    panx = pany = 0;
    m_binPan = true;
    m_binGesture = true;
    // qDebug() << "pg1";
  }
 
  else if (event.IsGestureEnd()) {
    // qDebug() << "panGestureFinished";
    m_pParentCanvas->UpdateCanvasControlBar();
    m_binPan = false;
    m_gestureFinishTimer.Start(500, wxTIMER_ONE_SHOT);
  }
 
  else {
    if (m_binPan) {
      if (!g_GLOptions.m_bUseCanvasPanning) {
        // qDebug() << "slowpan" << dx << dy;
 
        m_pParentCanvas->FreezePiano();
        m_pParentCanvas->PanCanvas(dx, -dy);
        m_pParentCanvas->ThawPiano();
 
      } else {
        FastPan(dx, dy);
      }
 
      panx -= dx;
      pany -= dy;
    }
  }
 
  m_bgestureGuard = true;
  m_gestureEeventTimer.Start(500, wxTIMER_ONE_SHOT);
  m_bforcefull = false;
}
 
//  Generic wxWidgets gesture event processor
void glChartCanvas::OnEvtZoomGesture(wxZoomGestureEvent &event) {
  if (g_true_zoom) {
    float zoom_gain = 1.0;
    float zout_gain = 1.0;
 
    float last_zoom_val = m_step_zoom_val;  // m_total_zoom_val;
 
    float max_zoom_scale = 1000.;
    float min_zoom_scale = 2e8;
 
    if (event.GetZoomFactor() > 1)
      m_total_zoom_val = ((event.GetZoomFactor() - 1.0) * zoom_gain) + 1.0;
    else
      m_total_zoom_val = 1.0 - ((1.0 - event.GetZoomFactor()) * zout_gain);
 
    float inc_zoom_val =
        m_total_zoom_val / last_zoom_val;  // the incremental zoom
 
    double projected_scale = m_cache_vp.chart_scale;
 
    if (event.IsGestureStart()) {
      m_glstopwatch.Start();
      printf("\nStart--------------\n");
      m_binPinch = true;
      m_pParentCanvas->m_inPinch = true;
      m_binPan = false;  // cancel any tentative pan gesture, in case the "pan
                         // cancel" event was lost
      m_binGesture = true;
      m_pinchStart = event.GetPosition();
      m_lpinchPoint = m_pinchStart;
      m_total_zoom_val = 1.0;
      m_final_zoom_val = 1.0;
      m_step_zoom_val = 1.0;
 
      // Capture the lat/lon of the "pinch-point"
      m_pParentCanvas->GetCanvasPixPoint(
          event.GetPosition().x, event.GetPosition().y, m_pinchlat, m_pinchlon);
      m_zoom_inc = 1.0;
    }
 
    if (event.IsGestureEnd()) {
      // Some platforms generate spurious gestureEnd events. Guard for this.
      if (!m_binGesture) return;
      printf("EndZoom--------------\n");
 
      Refresh();
 
      // Let the gesture finish timer (500 msec) signal the pinch is done.
      // this will allow any lingering drag operation to avoid starting the
      // inertial drag process
      m_final_zoom_val = 1.0;
      m_total_zoom_val = 1.0;
      m_step_zoom_val = 1.0;
 
      m_gestureFinishTimer.Start(500, wxTIMER_ONE_SHOT);
    }
 
    else {
      // printf("%ld %6g %6g %6g \n", g_glstopwatch.Time(),
      // event.GetZoomFactor(),
      //        inc_zoom_val, projected_scale);
      float zoom_step = 5;
      float zoom_trigger = 0.05;
      if (projected_scale > 1e5)
        zoom_trigger = .1;
      else if (projected_scale < 3e4)
        zoom_trigger = .02;
 
      if (inc_zoom_val != 1.0) {
        if (inc_zoom_val > 1 + zoom_step) {
          m_step_zoom_val = m_step_zoom_val * (1 + zoom_step);
          printf("   Partial zoom: %6g\n", 1 + zoom_step);
          m_pParentCanvas->ZoomCanvasSimple(1 + zoom_step);
          if (IsShown()) SetCurrent(*m_pcontext);
          Render();
        } else {
          if (fabs(inc_zoom_val - 1.) > zoom_trigger) {
            m_step_zoom_val = m_total_zoom_val;
            printf("   Zoom: %6g\n", inc_zoom_val);
            m_pParentCanvas->ZoomCanvasSimple(inc_zoom_val);
 
            // Consider "zoom to pinch-point"
            // Disable ZTC if lookahead is ON, and currently b_follow is active
            bool b_allow_ztp = true;
            if (m_pParentCanvas->m_bFollow && m_pParentCanvas->m_bLookAhead)
              b_allow_ztp = false;
 
            if (g_bEnableZoomToCursor && b_allow_ztp) {
              // get the new (after zoom) pixel position of the original
              // pinch-point
              wxPoint r;
              m_pParentCanvas->GetCanvasPointPix(m_pinchlat, m_pinchlon, &r);
              // Shift the canvas to follow the pinch-lat/lon
              int dx = r.x - event.GetPosition().x;
              int dy = r.y - event.GetPosition().y;
              m_pParentCanvas->PanCanvas(dx, dy);
            }
 
            if (IsShown()) SetCurrent(*m_pcontext);
            Render();
          }
        }
      }
    }
  }  // true-zoom
  else {
    // Legacy zoom behavior
    float zoom_gain = 1.0;
    float zout_gain = 1.0;
 
    float last_zoom_val = m_total_zoom_val;
 
    float max_zoom_scale = 1000.;
    float min_zoom_scale = 2e8;
 
    if (event.GetZoomFactor() > 1)
      m_total_zoom_val = ((event.GetZoomFactor() - 1.0) * zoom_gain) + 1.0;
    else
      m_total_zoom_val = 1.0 - ((1.0 - event.GetZoomFactor()) * zout_gain);
 
    float inc_zoom_val =
        m_total_zoom_val / last_zoom_val;  // the incremental zoom
 
    double projected_scale =
        m_pParentCanvas->GetVP().chart_scale / m_total_zoom_val;
 
    if (event.IsGestureStart()) {
      // printf("\nStart--------------\n");
      m_first_zout = false;
      m_binPinch = true;
      m_binPan = false;  // cancel any tentative pan gesture, in case the "pan
                         // cancel" event was lost
      m_binGesture = true;
      m_pinchStart = event.GetPosition();
      m_lpinchPoint = m_pinchStart;
      m_total_zoom_val = 1.0;
      m_final_zoom_val = 1.0;
 
      m_pParentCanvas->GetCanvasPixPoint(
          event.GetPosition().x, event.GetPosition().y, m_pinchlat, m_pinchlon);
      //            qDebug() << "center" << event.GetCenterPoint().x <<
      //            event.GetCenterPoint().y;
 
      m_cc_x = m_fbo_offsetx + (m_fbo_swidth / 2);
      m_cc_y = m_fbo_offsety + (m_fbo_sheight / 2);
 
      // Render the full charts with overlay objects onto the frame buffer.
      if (IsShown()) SetCurrent(*m_pcontext);
      RenderScene();
#ifndef __ANDROID__
      ViewPort vpr = m_pParentCanvas->VPoint;
      m_texVP = vpr;
      GetTouchBackingBitmap(vpr);
#endif
      m_zoom_inc = 1.0;
    }
 
    if (event.IsGestureEnd()) {
      // printf("End--------------\n");
      //             qDebug() << "finish totalzoom" << total_zoom_val <<
      //             projected_scale;
 
      // Some ptaforms generate spurious gestureEnd events. Guard for this.
      if (!m_binGesture) return;
 
      float cc_x = m_fbo_offsetx + (m_fbo_swidth / 2);
      float cc_y = m_fbo_offsety + (m_fbo_sheight / 2);
      float dy = 0;
      float dx = 0;
 
      float tzoom = m_final_zoom_val;
 
      dx = (cc_x - m_cc_x) * tzoom;
      dy = -(cc_y - m_cc_y) * tzoom;
 
      if (zoomTimer.IsRunning()) {
        //                qDebug() << "Final zoom";
        m_zoomFinal = true;
        m_zoomFinalZoom = tzoom;
        m_zoomFinaldx = dx;
        m_zoomFinaldy = dy;
      }
 
      else {
        double final_projected_scale =
            m_pParentCanvas->GetVP().chart_scale / tzoom;
        // qDebug() << "Final pinchB" << tzoom << final_projected_scale;
 
        if (final_projected_scale < min_zoom_scale) {
          // qDebug() << "zoomit";
          m_pParentCanvas->ZoomCanvas(tzoom, false);
          m_pParentCanvas->PanCanvas(dx, dy);
          m_pParentCanvas->m_pQuilt->Invalidate();
          m_bforcefull = true;
 
        } else {
          double new_scale =
              m_pParentCanvas->GetCanvasScaleFactor() / min_zoom_scale;
          // qDebug() << "clampit";
          m_pParentCanvas->SetVPScale(new_scale);
          m_pParentCanvas->m_pQuilt->Invalidate();
          m_bforcefull = true;
        }
      }
 
      m_binPinch = false;
      m_final_zoom_val = 1.0;
      m_total_zoom_val = 1.0;
      m_gestureFinishTimer.Start(500, wxTIMER_ONE_SHOT);
    }
 
    else {
      if (1 /* g_GLOptions.m_bUseCanvasPanning*/) {
        if (projected_scale < min_zoom_scale) {
          wxPoint pinchPoint = event.GetPosition();
 
          float dx = pinchPoint.x - m_lpinchPoint.x;
          float dy = pinchPoint.y - m_lpinchPoint.y;
 
          FastZoom(inc_zoom_val, m_pinchStart.x, m_pinchStart.y,
                   -dx / m_total_zoom_val, dy / m_total_zoom_val);
 
          m_lpinchPoint = pinchPoint;
          m_final_zoom_val *= inc_zoom_val;
        }
      } else {
        // qDebug() << "update totalzoom" << total_zoom_val << projected_scale;
        if (1 || ((m_total_zoom_val > 1) && !m_first_zout)) {  // Zoom in
          wxPoint pinchPoint = event.GetPosition();
 
          float dx = pinchPoint.x - m_lpinchPoint.x;
          float dy = pinchPoint.y - m_lpinchPoint.y;
 
          if ((projected_scale > max_zoom_scale) &&
              (projected_scale < min_zoom_scale))
            FastZoom(inc_zoom_val, m_pinchStart.x, m_pinchStart.y,
                     -dx / m_total_zoom_val, dy / m_total_zoom_val);
 
          m_lpinchPoint = pinchPoint;
          m_final_zoom_val *= inc_zoom_val;
 
        } else {
          m_first_zout = true;
          m_zoom_inc *= inc_zoom_val;
          if ((m_zoom_inc < 0.9) || (m_zoom_inc > 1.1)) {
            m_pParentCanvas->ZoomCanvas(m_zoom_inc, false);
            m_zoom_inc = 1.0;
          }
 
          wxPoint pinchPoint = event.GetPosition();
          float dx = pinchPoint.x - m_lpinchPoint.x;
          float dy = pinchPoint.y - m_lpinchPoint.y;
          m_pParentCanvas->PanCanvas(-dx, -dy);
          m_lpinchPoint = pinchPoint;
        }
      }
    }
    m_gestureEeventTimer.Start(500, wxTIMER_ONE_SHOT);
  }
 
  m_bgestureGuard = true;
  m_bpinchGuard = true;
}
 
void glChartCanvas::onGestureTimerEvent(wxTimerEvent &event) {
  //  On some devices, the pan GestureFinished event fails to show up
  //  Watch for this case, and fix it.....
 
  if (m_binPan) {
    m_binPan = false;
    Invalidate();
    Update();
  }
  m_bgestureGuard = false;
  m_bpinchGuard = false;
  m_binGesture = false;
  m_bforcefull = false;
}
 
void glChartCanvas::onGestureFinishTimerEvent(wxTimerEvent &event) {
  // qDebug() << "onGestureFinishTimerEvent";
 
  // signal gesture is finished after a delay
  m_binGesture = false;
  m_bforcefull = false;
  m_binPinch = false;
  m_pParentCanvas->m_inPinch = false;
  printf("******Finish\n");
}
 
#endif
#endif
 
void glChartCanvas::configureShaders(ViewPort &vp) {
#if defined(USE_ANDROID_GLES2) || defined(ocpnUSE_GLSL)
  mat4x4 I;
  mat4x4_identity(I);
 
  ViewPort *pvp = (ViewPort *)&vp;
 
  GLShaderProgram *shader = pcolor_tri_shader_program[GetCanvasIndex()];
  shader->Bind();
  shader->SetUniformMatrix4fv("MVMatrix", (GLfloat *)pvp->vp_matrix_transform);
  shader->SetUniformMatrix4fv("TransformMatrix", (GLfloat *)I);
  shader->UnBind();
 
  //       glUseProgram(color_tri_shader_program);
  //       GLint matloc = glGetUniformLocation(color_tri_shader_program,
  //       "MVMatrix"); glUniformMatrix4fv(matloc, 1, GL_FALSE,
  //                          (const GLfloat *)pvp->vp_transform);
  //       GLint transloc =
  //           glGetUniformLocation(color_tri_shader_program,
  //           "TransformMatrix");
  //       glUniformMatrix4fv(transloc, 1, GL_FALSE, (const GLfloat *)I);
 
  //       glUseProgram(texture_2D_shader_program);
  //       matloc = glGetUniformLocation(texture_2D_shader_program, "MVMatrix");
  //       glUniformMatrix4fv(matloc, 1, GL_FALSE,
  //                          (const GLfloat *)pvp->vp_transform);
  //       transloc =
  //           glGetUniformLocation(texture_2D_shader_program,
  //           "TransformMatrix");
  //       glUniformMatrix4fv(transloc, 1, GL_FALSE, (const GLfloat *)I);
 
  shader = ptexture_2D_shader_program[GetCanvasIndex()];
  shader->Bind();
  shader->SetUniformMatrix4fv("MVMatrix", (GLfloat *)pvp->vp_matrix_transform);
  shader->SetUniformMatrix4fv("TransformMatrix", (GLfloat *)I);
  shader->UnBind();
 
  //       glUseProgram(circle_filled_shader_program);
  //       matloc = glGetUniformLocation(circle_filled_shader_program,
  //       "MVMatrix"); glUniformMatrix4fv(matloc, 1, GL_FALSE,
  //                          (const GLfloat *)pvp->vp_transform);
  //       transloc =
  //           glGetUniformLocation(circle_filled_shader_program,
  //           "TransformMatrix");
  //       glUniformMatrix4fv(transloc, 1, GL_FALSE, (const GLfloat *)I);
 
  shader = pcircle_filled_shader_program[GetCanvasIndex()];
  shader->Bind();
  shader->SetUniformMatrix4fv("MVMatrix", (GLfloat *)pvp->vp_matrix_transform);
  shader->SetUniformMatrix4fv("TransformMatrix", (GLfloat *)I);
  shader->UnBind();
 
  shader = ptexture_2DA_shader_program[GetCanvasIndex()];
  shader->Bind();
  shader->SetUniformMatrix4fv("MVMatrix", (GLfloat *)pvp->vp_matrix_transform);
  shader->SetUniformMatrix4fv("TransformMatrix", (GLfloat *)I);
  shader->UnBind();
 
  // glUseProgram(AALine_shader_program);
  // matloc = glGetUniformLocation(AALine_shader_program, "MVMatrix");
  // glUniformMatrix4fv(matloc, 1, GL_FALSE,
  //                    (const GLfloat *)pvp->vp_transform);
 
  shader = pAALine_shader_program[GetCanvasIndex()];
  shader->Bind();
  shader->SetUniformMatrix4fv("MVMatrix", (GLfloat *)pvp->vp_matrix_transform);
  shader->UnBind();
 
  shader = pring_shader_program[GetCanvasIndex()];
  shader->Bind();
  shader->SetUniformMatrix4fv("MVMatrix", (GLfloat *)pvp->vp_matrix_transform);
  shader->SetUniformMatrix4fv("TransformMatrix", (GLfloat *)I);
  shader->UnBind();
 
  //  Leftover shader required by some older Android plugins
  if (texture_2DA_shader_program) {
    glUseProgram(texture_2DA_shader_program);
    GLint matloc = glGetUniformLocation(texture_2DA_shader_program, "MVMatrix");
    glUniformMatrix4fv(matloc, 1, GL_FALSE,
                       (const GLfloat *)pvp->vp_matrix_transform);
    GLint transloc =
        glGetUniformLocation(texture_2DA_shader_program, "TransformMatrix");
    glUniformMatrix4fv(transloc, 1, GL_FALSE, (const GLfloat *)I);
  }
 
  m_gldc.m_texfont.PrepareShader(vp.pix_width, vp.pix_height, vp.rotation);
 
#endif
}
 
void glChartCanvas::RenderTextures(ocpnDC &dc, float *coords, float *uvCoords,
                                   int nVertex, ViewPort *vp) {
// #ifdef USE_ANDROID_GLES2
#if defined(USE_ANDROID_GLES2) || defined(ocpnUSE_GLSL)
  int nl = nVertex / 4;
  float *lc = coords;
  float *luv = uvCoords;
 
  while (nl) {
    RenderSingleTexture(dc, lc, luv, vp, 0, 0, 0);
 
    lc += 8;
    luv += 8;
    nl--;
  }
 
#else
  glEnableClientState(GL_VERTEX_ARRAY);
  glEnableClientState(GL_TEXTURE_COORD_ARRAY);
 
  glTexCoordPointer(2, GL_FLOAT, 2 * sizeof(GLfloat), uvCoords);
  glVertexPointer(2, GL_FLOAT, 2 * sizeof(GLfloat), coords);
  glDrawArrays(GL_QUADS, 0, 4);
 
#endif
 
  return;
}
 
void glChartCanvas::RenderSingleTexture(ocpnDC &dc, float *coords,
                                        float *uvCoords, ViewPort *vp, float dx,
                                        float dy, float angle_rad) {
#if defined(USE_ANDROID_GLES2) || defined(ocpnUSE_GLSL)
 
  GLShaderProgram *shader = ptexture_2D_shader_program[dc.m_canvasIndex];
  if (!shader) return;
 
  shader->Bind();
 
  // Set up the texture sampler to texture unit 0
  shader->SetUniform1i("uTex", 0);
 
  // Rotate
  mat4x4 I, Q;
  mat4x4_identity(I);
  mat4x4_rotate_Z(Q, I, angle_rad);
 
  // Translate
  Q[3][0] = dx;
  Q[3][1] = dy;
 
  shader->SetUniformMatrix4fv("TransformMatrix", (GLfloat *)Q);
 
  float co1[8];
  float tco1[8];
 
  shader->SetAttributePointerf("aPos", co1);
  shader->SetAttributePointerf("aUV", tco1);
 
// Perform the actual drawing.
 
// For some reason, glDrawElements is busted on Android
// So we do this a hard ugly way, drawing two triangles...
#if 0
    GLushort indices1[] = {0,1,3,2};
    glDrawElements(GL_TRIANGLE_STRIP, 4, GL_UNSIGNED_SHORT, indices1);
#else
 
  co1[0] = coords[0];
  co1[1] = coords[1];
  co1[2] = coords[2];
  co1[3] = coords[3];
  co1[4] = coords[6];
  co1[5] = coords[7];
  co1[6] = coords[4];
  co1[7] = coords[5];
 
  tco1[0] = uvCoords[0];
  tco1[1] = uvCoords[1];
  tco1[2] = uvCoords[2];
  tco1[3] = uvCoords[3];
  tco1[4] = uvCoords[6];
  tco1[5] = uvCoords[7];
  tco1[6] = uvCoords[4];
  tco1[7] = uvCoords[5];
 
  // glVertexAttribPointer(mPosAttrib, 2, GL_FLOAT, GL_FALSE, 0, co1);
  // glVertexAttribPointer(mUvAttrib, 2, GL_FLOAT, GL_FALSE, 0, tco1);
 
  glDrawArrays(GL_TRIANGLE_STRIP, 0, 4);
 
  shader->UnBind();
 
#endif
 
#else
#endif
 
  return;
}
 
void glChartCanvas::RenderColorRect(wxRect r, wxColor &color) {
#if defined(USE_ANDROID_GLES2) || defined(ocpnUSE_GLSL)
 
  GLShaderProgram *shader = pcolor_tri_shader_program[GetCanvasIndex()];
  shader->Bind();
 
  shader->SetUniformMatrix4fv(
      "MVMatrix", (GLfloat *)m_pParentCanvas->GetpVP()->vp_matrix_transform);
 
  float colorv[4];
  colorv[0] = color.Red() / float(256);
  colorv[1] = color.Green() / float(256);
  colorv[2] = color.Blue() / float(256);
  colorv[3] = 1.0;
  shader->SetUniform4fv("color", colorv);
 
  float pf[8];
  pf[0] = r.x + r.width;
  pf[1] = r.y;
  pf[2] = r.x;
  pf[3] = r.y;
  pf[4] = r.x + r.width;
  pf[5] = r.y + r.height;
  pf[6] = r.x;
  pf[7] = r.y + r.height;
  shader->SetAttributePointerf("position", pf);
 
  glDrawArrays(GL_TRIANGLE_STRIP, 0, 4);
 
  shader->UnBind();
 
#else
#endif
}
 
void glChartCanvas::RenderScene(bool bRenderCharts, bool bRenderOverlays) {
#if defined(USE_ANDROID_GLES2) || defined(ocpnUSE_GLSL)
 
  ViewPort VPoint = m_pParentCanvas->VPoint;
  ocpnDC gldc(*this);
 
  int w, h;
  GetClientSize(&w, &h);
  int sx = GetSize().x;
  int sy = GetSize().y;
 
  OCPNRegion screen_region(wxRect(0, 0, VPoint.pix_width, VPoint.pix_height));
 
  glViewport(0, 0, (GLint)w, (GLint)h);
 
  if (s_b_useStencil) {
    glEnable(GL_STENCIL_TEST);
    glStencilMask(0xff);
    glClear(GL_STENCIL_BUFFER_BIT);
    glDisable(GL_STENCIL_TEST);
  }
 
  // Make sure we have a valid quilt composition
  m_pParentCanvas->m_pQuilt->Compose(m_pParentCanvas->VPoint);
 
  // set opengl settings that don't normally change
  // this should be able to go in SetupOpenGL, but it's
  // safer here incase a plugin mangles these
  glHint(GL_LINE_SMOOTH_HINT, GL_NICEST);
  glHint(GL_POLYGON_SMOOTH_HINT, GL_NICEST);
  glBlendFunc(GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
 
  // enable rendering to texture in framebuffer object
  glBindFramebuffer(GL_FRAMEBUFFER_EXT, m_fb0);
 
  glFramebufferTexture2D(GL_FRAMEBUFFER_EXT, GL_COLOR_ATTACHMENT0,
                         g_texture_rectangle_format, m_cache_tex[m_cache_page],
                         0);
 
  m_fbo_offsetx = 0;
  m_fbo_offsety = 0;
  m_fbo_swidth = sx;
  m_fbo_sheight = sy;
 
  if (bRenderCharts) RenderCharts(gldc, screen_region);
 
  if (bRenderOverlays) {
    RenderS57TextOverlay(m_pParentCanvas->VPoint);
    RenderMBTilesOverlay(m_pParentCanvas->VPoint);
    g_overlayCanvas = m_pParentCanvas;
 
    if (g_pi_manager) {
      g_pi_manager->SendViewPortToRequestingPlugIns(VPoint);
      g_pi_manager->RenderAllGLCanvasOverlayPlugIns(
          m_pcontext, VPoint, m_pParentCanvas->m_canvasIndex, OVERLAY_CHARTS);
    }
 
    DrawStaticRoutesTracksAndWaypoints(m_pParentCanvas->VPoint);
    DrawDynamicRoutesTracksAndWaypoints(VPoint);
    DrawFloatingOverlayObjects(m_gldc);
  }
 
  // All done, so disable Render to FBO
  glBindFramebuffer(GL_FRAMEBUFFER, 0);
 
#endif
}
 
wxBitmap &glChartCanvas::GetTouchBackingBitmap(ViewPort &vp) {
  wxBitmap tbm(vp.pix_width, vp.pix_height, -1);
  wxMemoryDC tdc(tbm);
  tdc.SetBackground(wxBrush(GetGlobalColor("BLUEBACK")));
  tdc.Clear();
  ocpnDC dc = ocpnDC(tdc);
  ViewPort tvp = vp;
 
  tvp.view_scale_ppm /= 2;
  tvp.SetBoxes();
 
  gShapeBasemap.SetBasemapLandColor(GetGlobalColor("LANDBACK"));
  dc.SetPen(*wxTRANSPARENT_PEN);
 
  gShapeBasemap.RenderViewOnDC(dc, tvp);
  tdc.SelectObject(wxNullBitmap);
  m_touch_backing_bitmap = tbm;
  CreateBackingTexture();
 
  return m_touch_backing_bitmap;
}
 
void glChartCanvas::CreateBackingTexture() {
  wxImage image = m_touch_backing_bitmap.ConvertToImage();
  unsigned char *imgdata = image.GetData();
  unsigned char *imgalpha = image.GetAlpha();
  m_tex_w = image.GetWidth();
  m_tex_h = image.GetHeight();
  m_image_width = m_tex_w;
  m_image_height = m_tex_h;
 
  GLuint format = GL_RGBA;
  GLuint internalformat = g_texture_rectangle_format;
#ifndef __ANDROID__
  internalformat = GL_RGBA;
#endif
  int stride = 4;
 
  if (imgdata) {
    unsigned char *teximage =
        (unsigned char *)malloc(stride * m_tex_w * m_tex_h);
 
    for (int i = 0; i < m_image_height; i++) {
      for (int j = 0; j < m_image_width; j++) {
        int s = (i * 3 * m_image_width) + (j * 3);
        int d = (i * stride * m_tex_w) + (j * stride);
 
        teximage[d + 0] = imgdata[s + 0];
        teximage[d + 1] = imgdata[s + 1];
        teximage[d + 2] = imgdata[s + 2];
        teximage[d + 3] = 255;
      }
    }
 
    glGenTextures(1, &m_TouchBackingTexture);
    glBindTexture(GL_TEXTURE_2D, m_TouchBackingTexture);
 
    glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_S, GL_CLAMP_TO_EDGE);
    glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_T, GL_CLAMP_TO_EDGE);
    glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER,
                    GL_NEAREST);  // No mipmapping
    glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);
 
    glTexImage2D(GL_TEXTURE_2D, 0, internalformat, m_tex_w, m_tex_h, 0, format,
                 GL_UNSIGNED_BYTE, teximage);
 
    free(teximage);
    glBindTexture(GL_TEXTURE_2D, 0);
  }
}