alg/glsequence.hpp

   1 #ifndef _GL_SEQUENCE_H_
   2 #define _GL_SEQUENCE_H_
   3
   4 #include <boost/shared_ptr.hpp>
   5
   6 #include "alg/annotation_colors.hpp"
   7 #include "alg/sequence.hpp"
   8 #include "alg/color.hpp"
   9 // isn't it lovely how the header file can live in different places?
  10 #ifdef __APPLE__
  11 #include <OpenGL/gl.h>
  12 #else
  13 #include <GL/gl.h>
  14 #endif
  15
  16 class GlSequence;
  17 typedef boost::shared_ptr<GlSequence> GlSequenceRef;
  18
  19 //! Manage rendering a mussa sequence track
  20 /*! The idea is this will keep track of the location of where the sequence
  21  *  is being rendered, and handle displaying annotations on that track
  22  */
  23 class GlSequence : public Sequence
  24 {
  25 public:
  26   GlSequence(const Sequence & s,
  27              boost::shared_ptr<AnnotationColors> cm);
  28   GlSequence(const GlSequence & s);
  29   GlSequence(const GlSequence *);
  30   //! Make a new GlSequence, using a copy of SeqSpan
  31   GlSequence(const GlSequenceRef);
  32
  33   GlSequence &operator=(const GlSequence &s);
  34
  35   static DrawableRef default_drawable();
  36   static DrawableRef copy_drawable(DrawableRef s);
  37
  38   //! search through sequence and add an annotation highlighting all the non AGCT characters
  39   void add_annotations_for_defined_sequence(Drawable::draw_func_ptr draw=0);
  40   /**!search through our annotations and set the draw function for a particular "type"
  41    * this is really a temporary solution to configuring annotation display styles
  42    */
  43   void update_annotation_draw_function(
  44          std::string type,
  45          Drawable::draw_func_ptr draw,
  46          ColorRef color
  47        );
  48   //! set our starting x (horizontal) coordinate
  49   void setX(float x);
  50   //! get our starting x (horizontal) coordinate
  51   float x() const;
  52   //! set our current y (vertical) position
  53   void setY(float y);
  54   //! get our current y (vertical) position
  55   float y() const;
  56   //! get our current z (depth) position
  57   float z() const;
  58   //! how thick (high) the track we're drawing is
  59   float height() const;
  60
  61   //! set our default draw color
  62   void setColor(ColorRef &c );
  63   //! return our draw color
  64   const ColorRef color() const;
  65   //! return our draw color
  66   ColorRef color();
  67   //! return default gene color
  68   static ColorRef default_gene_color();
  69   //! return default track color
  70   static ColorRef default_track_color();
  71
  72   //! return our annotation color mapper
  73   AnnotationColorsRef colorMapper() { return color_mapper; }
  74
  75   //! draw a from left to right +/- height/2
  76   static
  77   void draw_box(GLfloat world_left, GLfloat world_right,
  78                 GLfloat left, GLfloat right, GLfloat height,
  79                 GLfloat y, GLfloat z,
  80                 GLint primitive=GL_QUADS);
  81   //! draw a track
  82   /*! left and right are the edges of the current viewport
  83    */
  84   void draw(GLfloat left, GLfloat right) const;
  85
  86   //! get our right (horizontal) coordinate (size-x)
  87   float right() const;
  88   //! how long is our sequence track? (computed from the sequence)
  89   GLfloat size() const;
  90
  91   //! return the left (lowest) base index that is fully visible
  92   Sequence::size_type leftbase(GLfloat left) const;
  93   //! return one past rightmost (highest ) base index that is fully visible
  94   //! done mostly so all the iterator logic continues to work correctly.
  95   Sequence::size_type rightbase(GLfloat right) const;
  96
  97   //! provide an iterator to the sequence starting at world coordinate left
  98   Sequence::const_iterator region_begin(GLfloat left, GLfloat right) const;
  99   //! provide an iterator to the sequence ending at world coordinate right
 100   Sequence::const_iterator region_end(GLfloat left, GLfloat right) const;
 101
 102   //! return a subsequence as a GlSequence (instead of a Sequence subsequence)
 103   GlSequence subseq(size_type start, size_type count) const;
 104
 105   //! Return the pixel width of the opengl viewport.
 106   static int get_viewport_width_in_pixels();
 107   //! how big is a pixel in world coordinats
 108   static GLfloat pixelWidth(GLfloat, GLfloat);
 109   //! how big is a pixel in world coordinats (specifying viewport size)
 110   static GLfloat pixelWidth(GLfloat, GLfloat, int);
 111
 112   //! are we close enough that it would make sense to view the base pairs?
 113   static bool is_sequence_renderable(GLfloat left, GLfloat right);
 114   //! are we close enough that it would make sense to view the base pairs?
 115   /*! though we don't actually check to see if there's sequence in our
 116    *  view, just that there's enough pixels to render something if
 117    *  there were.
 118    *  \param[in] left the left edge of the viewable region in world coordinates
 119    *  \param[in] right the right edge of the viewable region in world
 120    *             coordinates
 121    *  \param[in] pixel_width allow setting the current viewport pixel width
 122    */
 123   static bool is_sequence_renderable(GLfloat, GLfloat, int);
 124
 125   friend bool operator==(const GlSequence &left, const GlSequence &right);
 126
 127 protected:
 128   AnnotationColorsRef color_mapper;
 129
 130   //! initalize a SeqSpanRef covering start, count
 131   SeqSpanRef make_drawable_annotation(
 132     Drawable::draw_func_ptr draw,
 133     std::string name,
 134     size_type start,
 135     size_type count,
 136     ColorRef color);
 137
 138   //! draw sequence as a bar
 139   void draw_track(GLfloat, GLfloat) const;
 140   void draw_annotations(GLfloat, GLfloat) const;
 141   //! render a sequence as ATGC
 142   /*! left and right are the current edges of the viewable world
 143    */
 144   void draw_sequence(GLfloat, GLfloat) const;
 145 };
 146
 147 //! draw sequence track when we're far enough out not to see the AGCTs.
 148 void draw_summarized_track(SeqSpanRef ref, float left, float right);
 149 //! draw a 1/4 width centered track
 150 void draw_narrow_track(SeqSpanRef ref, float left, float right);
 151 #endif