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

生物技术

开发平台:

C/C++

  1. /*
  2.  * ===========================================================================
  3.  * PRODUCTION $Log: sequence_item.cpp,v $
  4.  * PRODUCTION Revision 1000.1  2004/06/01 19:45:28  gouriano
  5.  * PRODUCTION PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.7
  6.  * PRODUCTION
  7.  * ===========================================================================
  8.  */
  9. /*  $Id: sequence_item.cpp,v 1000.1 2004/06/01 19:45:28 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. * Author:  Mati Shomrat, NCBI
  35. *
  36. * File Description:
  37. *   flat-file generator -- origin item implementation
  38. *
  39. */
  40. #include <ncbi_pch.hpp>
  41. #include <corelib/ncbistd.hpp>
  42. #include <objects/seq/Bioseq.hpp>
  43. #include <objects/seqloc/Seq_loc.hpp>
  44. #include <objects/seqloc/Seq_interval.hpp>
  45. #include <objmgr/util/sequence.hpp>
  46. #include <objtools/format/formatter.hpp>
  47. #include <objtools/format/text_ostream.hpp>
  48. #include <objtools/format/items/sequence_item.hpp>
  49. #include <objtools/format/context.hpp>
  50. BEGIN_NCBI_SCOPE
  51. BEGIN_SCOPE(objects)
  52. CSequenceItem::CSequenceItem
  53. (TSeqPos from,
  54.  TSeqPos to,
  55.  bool first, 
  56.  CBioseqContext& ctx) :
  57.     CFlatItem(&ctx),
  58.     m_From(from), m_To(to - 1), m_First(first)
  59. {
  60.     x_GatherInfo(ctx);
  61. }
  62. void CSequenceItem::Format
  63. (IFormatter& formatter,
  64.  IFlatTextOStream& text_os) const
  65. {
  66.     formatter.FormatSequence(*this, text_os);
  67. }
  68. const CSeqVector& CSequenceItem::GetSequence(void) const
  69. {
  70.     return m_Sequence;
  71. }
  72. TSeqPos CSequenceItem::GetFrom(void) const
  73. {
  74.     return m_From + 1;
  75. }
  76. TSeqPos CSequenceItem::GetTo(void) const
  77. {
  78.     return m_To + 1;
  79. }
  80. bool CSequenceItem::IsFirst(void) const
  81. {
  82.     return m_First;
  83. }
  84. /***************************************************************************/
  85. /*                                  PRIVATE                                */
  86. /***************************************************************************/
  87. void CSequenceItem::x_GatherInfo(CBioseqContext& ctx)
  88. {
  89.     x_SetObject(*ctx.GetHandle().GetCompleteBioseq());
  90.     const CSeq_loc& loc = ctx.GetLocation();
  91.     TSeqPos offset = loc.GetStart(kInvalidSeqPos);
  92.     TSeqPos start = m_From + offset;
  93.     TSeqPos end   = min(m_To + offset, ctx.GetHandle().GetInst_Length() - 1);
  94.     
  95.     CSeq_loc region;
  96.     CSeq_loc::TInt& ival = region.SetInt();
  97.     region.SetId(*ctx.GetPrimaryId());
  98.     ival.SetFrom(start);
  99.     ival.SetTo(end);
  100.     
  101.     m_Sequence = 
  102.         ctx.GetHandle().GetSequenceView(region, CBioseq_Handle::eViewConstructed);
  103.     CSeqVector::TCoding code = CSeq_data::e_Iupacna;
  104.     if ( ctx.IsProt() ) {
  105.         code = ctx.Config().IsModeRelease() ?
  106.             CSeq_data::e_Iupacaa : CSeq_data::e_Ncbieaa;
  107.     }
  108.     m_Sequence.SetCoding(code);
  109. }
  110. END_SCOPE(objects)
  111. END_NCBI_SCOPE
  112. /*
  113. * ===========================================================================
  114. *
  115. * $Log: sequence_item.cpp,v $
  116. * Revision 1000.1  2004/06/01 19:45:28  gouriano
  117. * PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.7
  118. *
  119. * Revision 1.7  2004/05/21 21:42:54  gorelenk
  120. * Added PCH ncbi_pch.hpp
  121. *
  122. * Revision 1.6  2004/05/07 15:23:51  shomrat
  123. * Removed unreferenced varaible
  124. *
  125. * Revision 1.5  2004/04/27 15:14:36  shomrat
  126. * Bug fix
  127. *
  128. * Revision 1.4  2004/04/22 15:56:17  shomrat
  129. * Support partial region
  130. *
  131. * Revision 1.3  2004/02/11 22:56:37  shomrat
  132. * use IsModeRelease method
  133. *
  134. * Revision 1.2  2003/12/18 17:43:36  shomrat
  135. * context.hpp moved
  136. *
  137. * Revision 1.1  2003/12/17 20:24:39  shomrat
  138. * Initial Revision (adapted from flat lib)
  139. *
  140. *
  141. * ===========================================================================
  142. */