X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Fdatamodel%2FSequenceGroup.java;h=0944c3465bad8d23900a57ceac74543f0d1d1649;hb=ab43013b7e357b84b4abade0dba949668dfb2a0e;hp=6bee4167fa15f0f8bba1dbbc0b4e5cfdeaf4b499;hpb=bbfdb203ce3a889600caa52478478b75ab03809f;p=jalview.git diff --git a/src/jalview/datamodel/SequenceGroup.java b/src/jalview/datamodel/SequenceGroup.java index 6bee416..0944c34 100755 --- a/src/jalview/datamodel/SequenceGroup.java +++ b/src/jalview/datamodel/SequenceGroup.java @@ -1,28 +1,37 @@ /* - * Jalview - A Sequence Alignment Editor and Viewer (Version 2.7) - * Copyright (C) 2011 J Procter, AM Waterhouse, J Engelhardt, LM Lui, G Barton, M Clamp, S Searle + * Jalview - A Sequence Alignment Editor and Viewer (Version 2.8.2b1) + * Copyright (C) 2014 The Jalview Authors * * This file is part of Jalview. * * Jalview is free software: you can redistribute it and/or * modify it under the terms of the GNU General Public License - * as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version. - * + * as published by the Free Software Foundation, either version 3 + * of the License, or (at your option) any later version. + * * Jalview is distributed in the hope that it will be useful, but * WITHOUT ANY WARRANTY; without even the implied warranty * of MERCHANTABILITY or FITNESS FOR A PARTICULAR * PURPOSE. See the GNU General Public License for more details. * - * You should have received a copy of the GNU General Public License along with Jalview. If not, see . + * You should have received a copy of the GNU General Public License + * along with Jalview. If not, see . + * The Jalview Authors are detailed in the 'AUTHORS' file. */ package jalview.datamodel; -import java.util.*; - -import java.awt.*; +import jalview.analysis.AAFrequency; +import jalview.analysis.Conservation; +import jalview.schemes.ColourSchemeI; +import jalview.schemes.ResidueProperties; -import jalview.analysis.*; -import jalview.schemes.*; +import java.awt.Color; +import java.util.ArrayList; +import java.util.Enumeration; +import java.util.Hashtable; +import java.util.List; +import java.util.Map; +import java.util.Vector; /** * Collects a set contiguous ranges on a set of sequences @@ -30,7 +39,7 @@ import jalview.schemes.*; * @author $author$ * @version $Revision$ */ -public class SequenceGroup +public class SequenceGroup implements AnnotatedCollectionI { String groupName; @@ -54,7 +63,7 @@ public class SequenceGroup /** * group members */ - private Vector sequences = new Vector(); + private Vector sequences = new Vector(); /** * representative sequence for this group (if any) @@ -91,12 +100,12 @@ public class SequenceGroup * consensus calculation property */ private boolean showSequenceLogo = false; + /** * flag indicating if logo should be rendered normalised */ private boolean normaliseSequenceLogo; - /** * @return the includeAllConsSymbols */ @@ -152,7 +161,7 @@ public class SequenceGroup if (seqsel != null) { sequences = new Vector(); - Enumeration sq = seqsel.sequences.elements(); + Enumeration sq = seqsel.sequences.elements(); while (sq.hasMoreElements()) { sequences.addElement(sq.nextElement()); @@ -169,7 +178,9 @@ public class SequenceGroup endRes = seqsel.endRes; cs = seqsel.cs; if (seqsel.description != null) + { description = new String(seqsel.description); + } hidecols = seqsel.hidecols; hidereps = seqsel.hidereps; idColour = seqsel.idColour; @@ -229,7 +240,9 @@ public class SequenceGroup } } if (!found) + { continue; + } } AlignmentAnnotation newannot = new AlignmentAnnotation( seq.getAnnotation()[a]); @@ -285,7 +298,15 @@ public class SequenceGroup return eres; } - public Vector getSequences(Hashtable hiddenReps) + @Override + public List getSequences() + { + return sequences; + } + + @Override + public List getSequences( + Map hiddenReps) { if (hiddenReps == null) { @@ -294,17 +315,16 @@ public class SequenceGroup else { Vector allSequences = new Vector(); - SequenceI seq, seq2; + SequenceI seq; for (int i = 0; i < sequences.size(); i++) { - seq = (SequenceI) sequences.elementAt(i); + seq = sequences.elementAt(i); allSequences.addElement(seq); if (hiddenReps.containsKey(seq)) { - SequenceGroup hsg = (SequenceGroup) hiddenReps.get(seq); - for (int h = 0; h < hsg.getSize(); h++) + SequenceCollectionI hsg = hiddenReps.get(seq); + for (SequenceI seq2 : hsg.getSequences()) { - seq2 = hsg.getSequenceAt(h); if (seq2 != seq && !allSequences.contains(seq2)) { allSequences.addElement(seq2); @@ -317,20 +337,15 @@ public class SequenceGroup } } - public SequenceI[] getSequencesAsArray(Hashtable hiddenReps) + public SequenceI[] getSequencesAsArray( + Map map) { - Vector tmp = getSequences(hiddenReps); + List tmp = getSequences(map); if (tmp == null) { return null; } - SequenceI[] result = new SequenceI[tmp.size()]; - for (int i = 0; i < result.length; i++) - { - result[i] = (SequenceI) tmp.elementAt(i); - } - - return result; + return tmp.toArray(new SequenceI[tmp.size()]); } /** @@ -468,12 +483,29 @@ public class SequenceGroup } /** - * Max Gaps Threshold for performing a conservation calculation TODO: make - * this a configurable property - or global to an alignment view + * Max Gaps Threshold (percent) for performing a conservation calculation */ private int consPercGaps = 25; /** + * @return Max Gaps Threshold for performing a conservation calculation + */ + public int getConsPercGaps() + { + return consPercGaps; + } + + /** + * set Max Gaps Threshold (percent) for performing a conservation calculation + * + * @param consPercGaps + */ + public void setConsPercGaps(int consPercGaps) + { + this.consPercGaps = consPercGaps; + } + + /** * calculate residue conservation for group - but only if necessary. */ public void recalcConservation() @@ -482,23 +514,17 @@ public class SequenceGroup { return; } - try { Hashtable cnsns[] = AAFrequency.calculate(sequences, startRes, endRes + 1, showSequenceLogo); if (consensus != null) { - _updateConsensusRow(cnsns); + _updateConsensusRow(cnsns, sequences.size()); } if (cs != null) { cs.setConsensus(cnsns); - - if (cs instanceof ClustalxColourScheme) - { - ((ClustalxColourScheme) cs).resetClustalX(sequences, getWidth()); - } } if ((conservation != null) @@ -518,15 +544,13 @@ public class SequenceGroup if (cs.conservationApplied()) { cs.setConservation(c); - - if (cs instanceof ClustalxColourScheme) - { - ((ClustalxColourScheme) cs).resetClustalX(sequences, - getWidth()); - } } } } + if (cs != null) + { + cs.alignmentChanged(context != null ? context : this, null); + } } catch (java.lang.OutOfMemoryError err) { // TODO: catch OOM @@ -557,7 +581,7 @@ public class SequenceGroup public Hashtable[] consensusData = null; - private void _updateConsensusRow(Hashtable[] cnsns) + private void _updateConsensusRow(Hashtable[] cnsns, long nseq) { if (consensus == null) { @@ -574,9 +598,10 @@ public class SequenceGroup consensus.annotations = new Annotation[aWidth]; // should be alignment width AAFrequency.completeConsensus(consensus, cnsns, startRes, endRes + 1, - ignoreGapsInConsensus, showSequenceLogo); // TODO: setting container - // for - // ignoreGapsInConsensusCalculation); + ignoreGapsInConsensus, showSequenceLogo, nseq); // TODO: setting + // container + // for + // ignoreGapsInConsensusCalculation); } /** @@ -617,20 +642,21 @@ public class SequenceGroup } /** - * DOCUMENT ME! * - * @return DOCUMENT ME! + * + * @return the first column selected by this group. Runs from 0<=i 0) { - width = ((SequenceI) sequences.elementAt(0)).getLength(); + width = sequences.elementAt(0).getLength(); } for (int i = 1; i < sequences.size(); i++) { - SequenceI seq = (SequenceI) sequences.elementAt(i); + SequenceI seq = sequences.elementAt(i); if (seq.getLength() > width) { @@ -953,18 +980,19 @@ public class SequenceGroup * * @param alignment * (may not be null) - * @param hashtable + * @param map * (may be null) * @return new group containing sequences common to this group and alignment */ - public SequenceGroup intersect(AlignmentI alignment, Hashtable hashtable) + public SequenceGroup intersect(AlignmentI alignment, + Map map) { SequenceGroup sgroup = new SequenceGroup(this); SequenceI[] insect = getSequencesInOrder(alignment); sgroup.sequences = new Vector(); for (int s = 0; insect != null && s < insect.length; s++) { - if (hashtable == null || hashtable.containsKey(insect[s])) + if (map == null || map.containsKey(insect[s])) { sgroup.sequences.addElement(insect[s]); } @@ -1038,12 +1066,12 @@ public class SequenceGroup { consensus = new AlignmentAnnotation("", "", new Annotation[1], 0f, 100f, AlignmentAnnotation.BAR_GRAPH); + consensus.hasText = true; + consensus.autoCalculated = true; + consensus.groupRef = this; + consensus.label = "Consensus for " + getName(); + consensus.description = "Percent Identity"; } - consensus.hasText = true; - consensus.autoCalculated = true; - consensus.groupRef = this; - consensus.label = "Consensus for " + getName(); - consensus.description = "Percent Identity"; return consensus; } @@ -1175,14 +1203,136 @@ public class SequenceGroup /** * set flag indicating if logo should be normalised when rendered + * * @param norm */ public void setNormaliseSequenceLogo(boolean norm) { - normaliseSequenceLogo=norm; + normaliseSequenceLogo = norm; } + public boolean isNormaliseSequenceLogo() { return normaliseSequenceLogo; } + + @Override + /** + * returns a new array with all annotation involving this group + */ + public AlignmentAnnotation[] getAlignmentAnnotation() + { + // TODO add in other methods like 'getAlignmentAnnotation(String label), + // etc' + ArrayList annot = new ArrayList(); + for (SequenceI seq : sequences) + { + AlignmentAnnotation[] aa = seq.getAnnotation(); + if (aa != null) + { + for (AlignmentAnnotation al : aa) + { + if (al.groupRef == this) + { + annot.add(al); + } + } + } + } + if (consensus != null) + { + annot.add(consensus); + } + if (conservation != null) + { + annot.add(conservation); + } + return annot.toArray(new AlignmentAnnotation[0]); + } + + @Override + public Iterable findAnnotation(String calcId) + { + ArrayList aa = new ArrayList(); + for (AlignmentAnnotation a : getAlignmentAnnotation()) + { + if (a.getCalcId() == calcId) + { + aa.add(a); + } + } + return aa; + } + + /** + * Returns a list of annotations that match the specified sequenceRef, calcId + * and label, ignoring null values. + * + * @return list of AlignmentAnnotation objects + */ + @Override + public Iterable findAnnotations(SequenceI seq, + String calcId, String label) + { + ArrayList aa = new ArrayList(); + for (AlignmentAnnotation ann : getAlignmentAnnotation()) + { + if (ann.getCalcId() != null && ann.getCalcId().equals(calcId) + && ann.sequenceRef != null && ann.sequenceRef == seq + && ann.label != null && ann.label.equals(label)) + { + aa.add(ann); + } + } + return aa; + } + + /** + * Answer true if any annotation matches the calcId passed in (if not null). + * + * @param calcId + * @return + */ + public boolean hasAnnotation(String calcId) + { + if (calcId != null && !"".equals(calcId)) + { + for (AlignmentAnnotation a : getAlignmentAnnotation()) + { + if (a.getCalcId() == calcId) + { + return true; + } + } + } + return false; + } + + public void clear() + { + sequences.clear(); + } + + private AnnotatedCollectionI context; + + /** + * set the alignment or group context for this group + * + * @param context + */ + public void setContext(AnnotatedCollectionI context) + { + this.context = context; + } + + /* + * (non-Javadoc) + * + * @see jalview.datamodel.AnnotatedCollectionI#getContext() + */ + @Override + public AnnotatedCollectionI getContext() + { + return context; + } }