开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 数值算法/人工智能 > 生物技术 > 查看源码
phylo_tree_rect_cladogram.cpp
资源名称：ncbi_cxx_win--Jul_08_2004.zip　[点击查看]
上传用户：yhdzpy8989
上传日期：2007-06-13
资源大小：13604k
文件大小：11k
源码类别：

生物技术

开发平台：

C/C++

phylo_tree_rect_cladogram.cpp：源码内容
  1. /*
  2.  * ===========================================================================
  3.  * PRODUCTION $Log: phylo_tree_rect_cladogram.cpp,v $
  4.  * PRODUCTION Revision 1000.1  2004/06/01 21:11:48  gouriano
  5.  * PRODUCTION PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.9
  6.  * PRODUCTION
  7.  * ===========================================================================
  8.  */
  10. /*  $Id: phylo_tree_rect_cladogram.cpp,v 1000.1 2004/06/01 21:11:48 gouriano Exp $
  11.  * ===========================================================================
  12.  *
  13.  *                            PUBLIC DOMAIN NOTICE
  14.  *               National Center for Biotechnology Information
  15.  *
  16.  *  This software/database is a "United States Government Work" under the
  17.  *  terms of the United States Copyright Act.  It was written as part of
  18.  *  the author's official duties as a United States Government employee and
  19.  *  thus cannot be copyrighted.  This software/database is freely available
  20.  *  to the public for use. The National Library of Medicine and the U.S.
  21.  *  Government have not placed any restriction on its use or reproduction.
  22.  *
  23.  *  Although all reasonable efforts have been taken to ensure the accuracy
  24.  *  and reliability of the software and data, the NLM and the U.S.
  25.  *  Government do not and cannot warrant the performance or results that
  26.  *  may be obtained by using this software or data. The NLM and the U.S.
  27.  *  Government disclaim all warranties, express or implied, including
  28.  *  warranties of performance, merchantability or fitness for any particular
  29.  *  purpose.
  30.  *
  31.  *  Please cite the author in any work or product based on this material.
  32.  *
  33.  * ===========================================================================
  34.  *
  35.  * Authors:  Vladimir Tereshkov
  36.  *
  37.  * File Description:
  38.  *
  39.  */
  41. #include <ncbi_pch.hpp>
  42. #include <corelib/ncbistl.hpp>
  44. #include <gui/opengl/glhelpers.hpp>
  45. #include <gui/graph/igraph_utils.hpp>
  46. #include <gui/widgets/phylo_tree/phylo_tree_render.hpp>
  47. #include <gui/widgets/phylo_tree/phylo_tree_rect_cladogram.hpp>
  51. #include <FL/Fl.H>
  52. #include <math.h>
  54. BEGIN_NCBI_SCOPE
  57. CPhyloRectCladogram::CPhyloRectCladogram()
  58. {  
  59. }
  61. CPhyloRectCladogram::~CPhyloRectCladogram()
  62. {
  63. }
  66. void  CPhyloRectCladogram::x_Layout(CPhyloTreeDataSource& ds)
  67. {    
  68.     CPhyloTreeNode * tree = ds.GetTree();    
  69.     
  70.     Int4 leafs = ds.GetSize();
  71.     Int4 width = ds.GetWidth();
  72.     
  73.     m_xStep = (m_DimX - m_LeftMargin - m_RightMargin) /  width;
  74.     m_yStep = (m_DimY - m_TopMargin - m_BottomMargin) / (leafs-1); 
  76.     if (ds.GetNormDistance() > 0){
  77.         m_NormDistance = (m_DimX - m_LeftMargin - m_RightMargin) / ds.GetNormDistance();
  78.     }
  80.     x_Calculate(ds.GetTree());    
  81. }
  83. void  CPhyloRectCladogram::x_Render(CGlPane& pane, CPhyloTreeDataSource& ds)
  84. {    
  85.     pane.OpenOrtho();    
  86.     x_DrawTree(ds.GetTree());
  87.     pane.Close();   
  88. }
  90. void CPhyloRectCladogram:: x_Calculate(CPhyloTreeNode * node)
  91. {
  92.     
  94.     if (node->IsLeaf()) {
  95.         
  96.         if (m_bDistMode){
  97.             node->XY().first  = m_LeftMargin + m_NormDistance * node->GetDistFromRoot();               
  98.         }
  99.         else {
  100.             node->XY().first  = (m_DimX - m_RightMargin);                                  // all leafs on the right
  101.         }
  103.         node->XY().second = (m_DimY - m_TopMargin - node->IDX().second * m_yStep);     // y coordinate 
  104.         node->SetAngle(0); // for correct labeling
  105.     }
  106.     else {        
  107.         for(CPhyloTreeNode::TNodeList_I  it = node->SubNodeBegin();  it != node->SubNodeEnd(); it++ )  x_Calculate((*it)->GetValue());         
  108.         
  109.         CPhyloTreeNode::TNodeList_I fst = node->SubNodeBegin();
  110.         CPhyloTreeNode::TNodeList_I lst = node->SubNodeEnd(); 
  111.         lst--;
  112.         
  113.         double yf = (*fst)->GetValue()->XY().second;
  114.         double yl = (*lst)->GetValue()->XY().second;
  115.                
  116.         
  117.         if (m_bDistMode){            
  118.             node->XY().first  = m_LeftMargin + m_NormDistance * node->GetDistFromRoot();               
  119.         }
  120.         else {
  121.             node->XY().first  = (m_LeftMargin + node->IDX().first * m_xStep);
  122.         }
  123.     
  124.         node->XY().second = (yl + (yf - yl)/2);            
  125.     }   
  126. }
  128. void CPhyloRectCladogram:: x_DrawTree(CPhyloTreeNode * node)
  129. {   
  130.     if (node->IsLeaf()) {
  131.          // label
  132.         if (m_pLabelFont) {
  133.             CGlPoint<int> labelPos(node->GetValue()->XY().first, node->GetValue()->XY().second);
  134.             m_Label.Render(*m_pPane, labelPos, node->GetLabel(), node->GetSelected());
  135.             //x_RenderText(node->GetValue()->XY().first, node->GetValue()->XY().second, node->GetLabel());
  136.         }
  137.             
  138.         // drawing leaf node                     
  139.         x_RenderLine(node->GetParent()->GetValue()->XY().first + m_LineWidth, node->GetValue()->XY().second,
  140.             node->GetValue()->XY().first, node->GetValue()->XY().second, m_LineWidth, node->GetSelected()?m_LineSelColor:m_LineColor, node->GetSelected()?m_LineSelHiColor:m_LineHiColor);
  141.        
  142.         x_RenderNode(node->GetValue()->XY().first, node->GetValue()->XY().second, m_NodeSize,  node->GetSelected()?m_NodeSelColor:m_NodeColor);
  143.     }
  144.     else {
  145.         // drawing internal node                   
  146.         CPhyloTreeNode::TNodeList_I fst = node->SubNodeBegin();
  147.         CPhyloTreeNode::TNodeList_I lst = node->SubNodeEnd(); 
  148.         lst--;
  150.         double ratio = m_pPane->GetScaleY()/m_pPane->GetScaleX();
  151.     
  152.         // vertical line
  153.         x_RenderLine(node->GetValue()->XY().first, (*fst)->GetValue()->XY().second - m_LineWidth * ratio,
  154.                      node->GetValue()->XY().first, (*lst)->GetValue()->XY().second + m_LineWidth * ratio, m_LineWidth, node->GetSelected()?m_LineSelColor:m_LineColor, node->GetSelected()?m_LineSelHiColor:m_LineHiColor);
  155.         
  157.         // horizontal line        
  158.         double xleft = node->GetParent()?node->GetParent()->GetValue()->XY().first:m_LeftMargin;
  159.         bool   bSel  = (node->GetSelected() && (node->GetParent() && node->GetParent()->GetValue()->GetSelected()));       
  161.         x_RenderLine(xleft + m_LineWidth, node->GetValue()->XY().second,
  162.                      node->GetValue()->XY().first + m_LineWidth, node->GetValue()->XY().second,m_LineWidth, bSel?m_LineSelColor:m_LineColor, bSel?m_LineSelHiColor:m_LineHiColor);
  164.         
  165.         // drawing top corner        
  166.         x_DrawCorner(node->GetValue()->XY().first, (*fst)->GetValue()->XY().second, eTop,    node->GetSelected()?m_LineSelColor:m_LineColor, node->GetSelected()?m_LineSelHiColor:m_LineHiColor);
  167.         
  168.         // bottom corner
  169.         x_DrawCorner(node->GetValue()->XY().first, (*lst)->GetValue()->XY().second, eBottom, node->GetSelected()?m_LineSelColor:m_LineColor, node->GetSelected()?m_LineSelHiColor:m_LineHiColor);        
  170.         
  171.         // middle on parent
  172.         if (node->GetParent() && (node!=*(node->GetParent()->SubNodeBegin()) && node!=*(--node->GetParent()->SubNodeEnd()))){
  173.             x_DrawCorner(xleft, node->GetValue()->XY().second, eMiddle, node->GetParent()->GetValue()->GetSelected()?m_LineSelColor:m_LineColor, node->GetParent()->GetValue()->GetSelected()?m_LineSelHiColor:m_LineHiColor);        
  174.         }
  175.         
  176.         //  node itself
  177.         x_RenderNode(node->GetValue()->XY().first, node->GetValue()->XY().second, m_NodeSize,   node->GetSelected()?m_NodeSelColor:m_NodeColor);
  179.         // draw childs
  180.          for(CPhyloTreeNode::TNodeList_I  it = node->SubNodeBegin();  it != node->SubNodeEnd(); it++ )  {
  181.             x_DrawTree((*it)->GetValue());                     
  182.         }
  183.        
  184.     }
  185. }
  186.     
  187. void  CPhyloRectCladogram::x_DrawCorner(double x, double y, ECorner corner, CGlColor color1, CGlColor color2)
  188. {
  189.     // reset line width
  190.     glEnable(GL_BLEND);
  191.     glBlendFunc(GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);  
  192.     glLineWidth(1.0);        
  193.     glShadeModel(GL_SMOOTH);
  194.     glEnable(GL_LINE_SMOOTH);
  197.     double ratio = m_pPane->GetScaleY()/m_pPane->GetScaleX();  
  199.     double dx = m_LineWidth;
  200.     double dy = m_LineWidth * ratio * ((corner==eBottom || corner==eMiddle)? 1: -1); 
  201.     
  202.     if (corner==eTop || corner==eBottom){
  204.         glBegin(GL_QUADS);
  205.         glColor4fv(color2.GetGlColor());
  206.         glVertex2d(x, y+dy); 
  207.         glColor4fv(color1.GetGlColor());
  208.         glVertex2d(x+dx,    y+dy);         
  209.         glColor4fv(color2.GetGlColor());
  210.         glVertex2d(x+dx, y); 
  211.         glVertex2d(x, y); 
  214.         glColor4fv(color2.GetGlColor());
  215.         glVertex2d(x, y); 
  216.         glVertex2d(x+dx, y); 
  217.         glColor4fv(color1.GetGlColor());
  218.         glVertex2d(x+dx, y-dy); 
  219.         glVertex2d(x, y-dy); 
  221.         glColor4fv(color2.GetGlColor());
  222.         glVertex2d(x, y+dy); 
  223.         glVertex2d(x, y); 
  224.         glColor4fv(color1.GetGlColor());
  225.         glVertex2d(x-dx, y); 
  226.         glVertex2d(x-dx, y+dy);         
  227.         glEnd();
  230.         glBegin(GL_TRIANGLES);    
  231.         glColor4fv(color2.GetGlColor());        
  232.         glVertex2d(x, y); 
  233.         glColor4fv(color1.GetGlColor());
  234.         glVertex2d(x, y-dy); 
  235.         glVertex2d(x-dx, y);        
  236.         glEnd();
  238.         glBegin(GL_LINES);    
  239.         glColor4fv(color1.GetGlColor());        
  240.         glVertex2d(x+dx, y-dy); glVertex2d(x, y-dy);    
  241.         glVertex2d(x, y-dy);    glVertex2d(x-dx, y);        
  242.         glVertex2d(x-dx, y);    glVertex2d(x-dx, y+dy);
  243.         glEnd();
  244.     }
  245.     else { // middle joint
  246.         glBegin(GL_QUADS);
  247.         glColor4fv(color2.GetGlColor());
  248.         glVertex2d(x+dx, y); 
  249.         glVertex2d(x, y); 
  250.         glVertex2d(x, y+dy); 
  251.         glColor4fv(color1.GetGlColor());
  252.         glVertex2d(x+dx, y+dy); 
  254.         glBegin(GL_QUADS);
  255.         glColor4fv(color2.GetGlColor());
  256.         glVertex2d(x, y-dy); 
  257.         glVertex2d(x, y); 
  258.         glVertex2d(x+dx, y);         
  259.         glColor4fv(color1.GetGlColor());
  260.         glVertex2d(x+dx, y-dy); 
  261.         glEnd();
  262.     }
  264.     glDisable(GL_BLEND);    
  265. }
  268. string CPhyloRectCladogram::GetDescription(void)
  269. {
  270.     return "Rectangle Cladogram";
  271. }
  273. END_NCBI_SCOPE
  275. /*
  276.  * ===========================================================================
  277.  * $Log: phylo_tree_rect_cladogram.cpp,v $
  278.  * Revision 1000.1  2004/06/01 21:11:48  gouriano
  279.  * PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.9
  280.  *
  281.  * Revision 1.9  2004/05/21 22:27:54  gorelenk
  282.  * Added PCH ncbi_pch.hpp
  283.  *
  284.  * Revision 1.8  2004/05/06 19:42:21  tereshko
  285.  * Rendering improvements
  286.  *
  287.  * Revision 1.7  2004/04/28 19:27:10  tereshko
  288.  * Added support for distances rendering
  289.  *
  290.  * Revision 1.6  2004/04/20 21:57:19  tereshko
  291.  * Major rendering, labeling and performance improvements
  292.  *
  293.  * Revision 1.5  2004/04/13 20:28:53  tereshko
  294.  * Numerous renderers improvements
  295.  *
  296.  * Revision 1.4  2004/04/01 21:47:25  tereshko
  297.  * Code clean-up
  298.  *
  299.  * Revision 1.3  2004/02/24 21:43:45  tereshko
  300.  * Added support for rendering joints
  301.  *
  302.  * Revision 1.2  2004/02/23 22:52:27  tereshko
  303.  * Rendering Improvements
  304.  *
  305.  * Revision 1.1  2004/02/13 17:05:08  tereshko
  306.  * Phylogenetic Tree Widget initial revision
  307.  *
  308.  * ===========================================================================
  309.  */