cross_aln_demo_ui.cpp
上传用户:yhdzpy8989
上传日期:2007-06-13
资源大小:13604k
文件大小:6k
源码类别:

生物技术

开发平台:

C/C++

  1. /*
  2.  * ===========================================================================
  3.  * PRODUCTION $Log: cross_aln_demo_ui.cpp,v $
  4.  * PRODUCTION Revision 1000.1  2004/06/01 21:06:44  gouriano
  5.  * PRODUCTION PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.3
  6.  * PRODUCTION
  7.  * ===========================================================================
  8.  */
  9. /*  $Id: cross_aln_demo_ui.cpp,v 1000.1 2004/06/01 21:06:44 gouriano Exp $
  10.  * ===========================================================================
  11.  *
  12.  *                            PUBLIC DOMAIN NOTICE
  13.  *               National Center for Biotechnology Information
  14.  *
  15.  *  This software/database is a "United States Government Work" under the
  16.  *  terms of the United States Copyright Act.  It was written as part of
  17.  *  the author's official duties as a United States Government employee and
  18.  *  thus cannot be copyrighted.  This software/database is freely available
  19.  *  to the public for use. The National Library of Medicine and the U.S.
  20.  *  Government have not placed any restriction on its use or reproduction.
  21.  *
  22.  *  Although all reasonable efforts have been taken to ensure the accuracy
  23.  *  and reliability of the software and data, the NLM and the U.S.
  24.  *  Government do not and cannot warrant the performance or results that
  25.  *  may be obtained by using this software or data. The NLM and the U.S.
  26.  *  Government disclaim all warranties, express or implied, including
  27.  *  warranties of performance, merchantability or fitness for any particular
  28.  *  purpose.
  29.  *
  30.  *  Please cite the author in any work or product based on this material.
  31.  *
  32.  * ===========================================================================
  33.  *
  34.  * Authors:  Vlad Lebedev
  35.  *
  36.  * File Description:
  37.  *
  38.  */
  39. #include <ncbi_pch.hpp>
  40. #include "cross_aln_demo_ui.hpp"
  41. #include <gui/utils/message_box.hpp>
  42. #include <gui/utils/fltk_utils.hpp>
  43. #include <objtools/data_loaders/genbank/gbloader.hpp>
  44. #include <serial/objistr.hpp>
  45. #include <objects/seq/Seq_annot.hpp>
  46. #include <serial/serial.hpp>
  47. BEGIN_NCBI_SCOPE
  48. USING_SCOPE(objects);
  49. #include "cross_aln_demo_ui_.cpp"
  50. CCrossAlnDemoUI::CCrossAlnDemoUI()
  51. {
  52.     m_Window.reset(x_CreateWindow());
  53.     CFltkEvent::SetOSDefaults();
  54. }
  55. CCrossAlnDemoUI::~CCrossAlnDemoUI()
  56. {
  57.     // these must be deleted in proper order
  58.     m_DataSource.Reset();
  59.     m_Window.reset();
  60.     m_Scope.Reset();
  61.     m_ObjMgr.Reset();
  62. }
  63. void CCrossAlnDemoUI::Show(int argc, char** argv)
  64. {
  65.     m_Accession->value("gi|");
  66.     m_Accession->redraw();
  67.     m_InputFile->value("/Users/lebedev/align2.asn");
  68.     m_InputFile->redraw();
  69.     m_Window->show(argc, argv);
  70.     while (m_Window->shown()) {
  71.         Fl::wait();
  72.     }
  73. }
  74. void CCrossAlnDemoUI::x_OnLoadAccession()
  75. {
  76.     CSeq_id id(m_Accession->value());
  77.     if (id.Which() == CSeq_id::e_not_set) {
  78.         string msg("Accession '");
  79.         msg += m_Accession->value();
  80.         msg += "' not recognized as a valid accession";
  81.         NcbiMessageBox(msg, eDialog_Ok, eIcon_Exclamation,
  82.                        "Unhandled Accession");
  83.         return;
  84.     }
  85.     if ( !m_ObjMgr ) {
  86.         m_ObjMgr.Reset(new CObjectManager());
  87.         m_ObjMgr->RegisterDataLoader(*new CGBDataLoader(),
  88.                                      CObjectManager::eDefault);
  89.         m_Scope.Reset(new CScope(*m_ObjMgr));
  90.         m_Scope->AddDefaults();
  91.     }
  92.     // retrieve our sequence
  93.     CBioseq_Handle handle = m_Scope->GetBioseqHandle(id);
  94.     if ( !handle ) {
  95.         string msg("Can't find sequence for accession '");
  96.         msg += m_Accession->value();
  97.         msg += "'";
  98.         NcbiMessageBox(msg, eDialog_Ok, eIcon_Exclamation,
  99.                        "Sequence Not Found");
  100.         return;
  101.     }
  102.     /*CRef<CSeq_align_set> aset(new CSeq_align_set());
  103.     CAlign_CI iter(handle, 0, 0, SAnnotSelector());
  104.     for ( ;  iter;  ++iter) {
  105.         const CSeq_align& align = *iter;
  106.         aset.insert(align);
  107.     }*/
  108.     //m_DataSource.Reset(new CCrossAlnDataSource(aset.GetObject(), m_Scope.GetObject()));
  109. }
  110. void CCrossAlnDemoUI::x_OnLoadFile()
  111. {
  112.     string filename(m_InputFile->value());
  113.     auto_ptr<CObjectIStream>    is(CObjectIStream::Open(eSerial_AsnText, filename));
  114.     //CRef<CSeq_align_set> aset(new CSeq_align_set());
  115.     //*is >> *aset;
  116.     CRef<CSeq_annot> annot(new CSeq_annot());
  117.     *is >> *annot;
  118.     if ( !m_ObjMgr ) {
  119.         m_ObjMgr.Reset(new CObjectManager());
  120.         m_ObjMgr->RegisterDataLoader(*new CGBDataLoader(),
  121.                                      CObjectManager::eDefault);
  122.         m_Scope.Reset(new CScope(*m_ObjMgr));
  123.         m_Scope->AddDefaults();
  124.     }
  125.     //m_DataSource.Reset(new CCrossAlnDataSource(aset.GetObject(), m_Scope.GetObject()));zzzzz
  126.     m_DataSource.Reset(new CCrossAlnDataSource(annot.GetObject(), m_Scope.GetObject()));
  127.     CCrossAlnDataSource::TIdRefCont all_ids =  m_DataSource->GetSeqIDs();
  128.     // choosing sequences to display
  129.     if(all_ids.size() >1 )  {
  130.         CCrossAlnDataSource::TIdRefCont::const_iterator  it = all_ids.begin();
  131.         CConstRef<CSeq_id>   s_id(*it);
  132.         CConstRef<CSeq_id>   q_id(*(++it));
  133.         m_DataSource->SelectIds(s_id, q_id);
  134.         m_AlnWidget->SetDataSource(m_DataSource.GetPointer());    
  135.         x_FillLists();
  136.     }
  137. }
  138. void CCrossAlnDemoUI::x_FillLists()
  139. {
  140.     m_Idx1->clear();
  141.     m_Idx2->clear();
  142.     const CCrossAlnDataSource::TIdRefCont& ids = m_DataSource->GetSeqIDs();
  143.   
  144.     ITERATE(CCrossAlnDataSource::TIdRefCont, iter, ids) {
  145.         CNcbiOstrstream oss;
  146.         (*iter)->WriteAsFasta(oss);
  147.         string tmp = (string)CNcbiOstrstreamToString(oss);
  148.         m_Idx1->add( tmp.c_str() );
  149.         m_Idx2->add( tmp.c_str() );
  150.     }
  151.     m_Idx1->select(1);
  152.     m_Idx2->select(2);
  153. }
  154. void CCrossAlnDemoUI::x_QuitApp()
  155. {
  156.     delete m_MainWindow;
  157. }
  158. END_NCBI_SCOPE
  159. /*
  160.  * ===========================================================================
  161.  * $Log: cross_aln_demo_ui.cpp,v $
  162.  * Revision 1000.1  2004/06/01 21:06:44  gouriano
  163.  * PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.3
  164.  *
  165.  * Revision 1.3  2004/05/21 22:27:52  gorelenk
  166.  * Added PCH ncbi_pch.hpp
  167.  *
  168.  * Revision 1.2  2004/01/05 18:50:29  vasilche
  169.  * Fixed path to include files.
  170.  *
  171.  * Revision 1.1  2003/12/22 13:09:17  lebedev
  172.  * Initial revision
  173.  *
  174.  * ===========================================================================
  175.  */