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

生物技术

开发平台:

C/C++

  1. /*
  2.  * ===========================================================================
  3.  * PRODUCTION $Log: featid_table.cpp,v $
  4.  * PRODUCTION Revision 1000.1  2004/06/01 21:00:37  gouriano
  5.  * PRODUCTION PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.3
  6.  * PRODUCTION
  7.  * ===========================================================================
  8.  */
  9. /*  $Id: featid_table.cpp,v 1000.1 2004/06/01 21:00:37 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:  Mike DiCuccio
  35.  *
  36.  * File Description:
  37.  *
  38.  */
  39. #include <ncbi_pch.hpp>
  40. #include "featid_table.hpp"
  41. #include <objects/general/Dbtag.hpp>
  42. #include <objects/general/Object_id.hpp>
  43. #include <objects/seqloc/Seq_interval.hpp>
  44. #include <objmgr/util/sequence.hpp>
  45. #include <set>
  46. #include <serial/iterator.hpp>
  47. BEGIN_NCBI_SCOPE
  48. USING_SCOPE(objects);
  49. CFeatIdTable::CFeatIdTable(int x, int y, int w, int h, const char* label)
  50.     : CTablePanel< CConstRef<CSeq_loc> >(x, y, w, h, label)
  51. {
  52.     SetColumn(0, "Source",     eString, FL_ALIGN_LEFT, 0.5f);
  53.     SetColumn(1, "Database",   eString, FL_ALIGN_LEFT, 0.5f);
  54.     SetColumn(2, "Identifier", eString, FL_ALIGN_LEFT, 1.0f);
  55. }
  56. void CFeatIdTable::Update(const CSeq_feat& feat)
  57. {
  58.     size_t row = 0;
  59.     // IDs from main location
  60.     {{
  61.         SetCell(row, 0) = "Sequence";
  62.         // use the whole spanned location, not the original location
  63.         const CSeq_id& id = sequence::GetId(feat.GetLocation());
  64.         CRef<CSeq_loc> loc(new CSeq_loc());
  65.         loc->SetInt().SetFrom(feat.GetLocation().GetTotalRange().GetFrom());
  66.         loc->SetInt().SetTo(feat.GetLocation().GetTotalRange().GetTo());
  67.         loc->SetInt().SetStrand(sequence::GetStrand(feat.GetLocation()));
  68.         loc->SetId(id);
  69.         SetData(row, loc);
  70.         SetCell(row, 1, x_GetSource(id));
  71.         SetCell(row, 2, id.AsFastaString());
  72.         ++row;
  73.     }}
  74.     // IDs from the product location, if it exists
  75.     if (feat.IsSetProduct()) {
  76.         SetCell(row, 0, "Product");
  77.         SetData(row, CConstRef<CSeq_loc>(&feat.GetProduct()));
  78.         const CSeq_id& id = sequence::GetId(feat.GetProduct());
  79.         SetCell(row, 1, x_GetSource(id));
  80.         SetCell(row, 2, id.AsFastaString());
  81.         ++row;
  82.     }
  83.     // IDs from DB xrefs
  84.     if (feat.IsSetDbxref()) {
  85.         SetCell(row, 0, "External DB");
  86.         ITERATE (CSeq_feat::TDbxref, iter, feat.GetDbxref()) {
  87.             const CDbtag& tag = **iter;
  88.             CRef<CSeq_loc> loc(new CSeq_loc());
  89.             loc->SetWhole(*x_GetSeqId(tag));
  90.             SetData(row, loc);
  91.             SetCell(row, 1, tag.GetDb());
  92.             if (tag.GetTag().IsStr()) {
  93.                 SetCell(row, 2, tag.GetTag().GetStr());
  94.             } else {
  95.                 SetCell(row, 2, NStr::UIntToString(tag.GetTag().GetId()));
  96.             }
  97.             ++row;
  98.         }
  99.     }
  100. }
  101. // get a string describing the source of a given seq-id
  102. string CFeatIdTable::x_GetSource(const CSeq_id& id) const
  103. {
  104.     switch (id.Which()) {
  105.     case CSeq_id::e_Local:
  106.         return "Local";
  107.     case CSeq_id::e_Gibbsq:
  108.     case CSeq_id::e_Gibbmt:
  109.         return "Geninfo Backbone";
  110.     case CSeq_id::e_Giim:
  111.         return "Geninfo Import";
  112.     case CSeq_id::e_Genbank:
  113.         return "Genbank";
  114.     case CSeq_id::e_Embl:
  115.         return "EMBL";
  116.     case CSeq_id::e_Pir:
  117.         return "PIR";
  118.     case CSeq_id::e_Swissprot:
  119.         return "Swiss-PROT";
  120.     case CSeq_id::e_Patent:
  121.         return "Patent";
  122.     case CSeq_id::e_Other:
  123.         return "Other";
  124.     case CSeq_id::e_General:
  125.         return "General";
  126.     case CSeq_id::e_Gi:
  127.         return "Genbank ID";
  128.     case CSeq_id::e_Ddbj:
  129.         return "DDBJ";
  130.     case CSeq_id::e_Prf:
  131.         return "Prf Sequence";
  132.     case CSeq_id::e_Pdb:
  133.         return "PDB";
  134.     case CSeq_id::e_Tpg:
  135.         return "Third-Party";
  136.     case CSeq_id::e_Tpe:
  137.         return "Third-Party/EMBL";
  138.     case CSeq_id::e_Tpd:
  139.         return "Third-Party/DDBJ";
  140.     default:
  141.         return "Unknown";
  142.         break;
  143.     }
  144. }
  145. // get a seq-id for a given Dbtag
  146. CSeq_id* CFeatIdTable::x_GetSeqId(const CDbtag& dbtag) const
  147. {
  148.     static map<string, string> s_IdDictionary;
  149.     if (s_IdDictionary.size() == 0) {
  150.         s_IdDictionary["ENSEMBL"] = "lcl|";
  151.         s_IdDictionary["GenBank"] = "gb|";
  152.         s_IdDictionary["RefSeq" ] = "ref|";
  153.         s_IdDictionary["EMBL"   ] = "gb|";
  154.         s_IdDictionary["DDBJ"   ] = "gb|";
  155.     }
  156.     string str(dbtag.GetDb() + "|");
  157.     map<string, string>::const_iterator iter =
  158.         s_IdDictionary.find(dbtag.GetDb());
  159.     if (iter != s_IdDictionary.end()) {
  160.         str = iter->second;
  161.     }
  162.     if (dbtag.GetTag().IsStr()) {
  163.         str += dbtag.GetTag().GetStr();
  164.     } else {
  165.         str += NStr::UIntToString(dbtag.GetTag().GetId());
  166.     }
  167.     CRef<CSeq_id> id(new CSeq_id(str));
  168.     if (id->Which() == CSeq_id::e_not_set) {
  169.         id.Reset(new CSeq_id("lcl|" + str));
  170.     }
  171.     return id.Release();
  172. }
  173. END_NCBI_SCOPE
  174. /*
  175.  * ===========================================================================
  176.  * $Log: featid_table.cpp,v $
  177.  * Revision 1000.1  2004/06/01 21:00:37  gouriano
  178.  * PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.3
  179.  *
  180.  * Revision 1.3  2004/05/21 22:27:49  gorelenk
  181.  * Added PCH ncbi_pch.hpp
  182.  *
  183.  * Revision 1.2  2004/01/20 18:17:53  dicuccio
  184.  * Changed to match new API in CTablePanel
  185.  *
  186.  * Revision 1.1  2004/01/13 20:38:46  dicuccio
  187.  * Added new view: feature IDs
  188.  *
  189.  * ===========================================================================
  190.  */