X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Fdatamodel%2FAlignmentAnnotation.java;h=0fcf43dc0ca1ff782c1482d516ac8c055c67bd08;hb=a8f483d04205bb8273ee311c12968b7e86d205fa;hp=20a3633499bb6ad34b38881de60969c6162266c6;hpb=86bf99474fd78966716a9d82ff925e4ec7c0663b;p=jalview.git diff --git a/src/jalview/datamodel/AlignmentAnnotation.java b/src/jalview/datamodel/AlignmentAnnotation.java index 20a3633..0fcf43d 100755 --- a/src/jalview/datamodel/AlignmentAnnotation.java +++ b/src/jalview/datamodel/AlignmentAnnotation.java @@ -1,36 +1,46 @@ /* - * Jalview - A Sequence Alignment Editor and Viewer - * Copyright (C) 2007 AM Waterhouse, J Procter, G Barton, M Clamp, S Searle - * - * This program 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 2 - * of the License, or (at your option) any later version. - * - * This program 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 this program; if not, write to the Free Software - * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA + * Jalview - A Sequence Alignment Editor and Viewer (Version 2.8.2) + * 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. + * + * 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 . + * The Jalview Authors are detailed in the 'AUTHORS' file. */ package jalview.datamodel; +import jalview.analysis.Rna; +import jalview.analysis.SecStrConsensus.SimpleBP; + +import jalview.analysis.WUSSParseException; + +import java.util.ArrayList; import java.util.Enumeration; import java.util.Hashtable; +import fr.orsay.lri.varna.models.rna.RNA; + /** * DOCUMENT ME! - * + * * @author $author$ * @version $Revision$ */ public class AlignmentAnnotation { - /** If true, this annotations is calculated every edit, - * eg consensus, quality or conservation graphs */ + /** + * If true, this annotations is calculated every edit, eg consensus, quality + * or conservation graphs + */ public boolean autoCalculated = false; public String annotationId; @@ -45,6 +55,54 @@ public class AlignmentAnnotation /** DOCUMENT ME!! */ public Annotation[] annotations; + + + + public ArrayList bps=null; + /** + * RNA secondary structure contact positions + */ + public SequenceFeature[] _rnasecstr = null; + + /** + * position of annotation resulting in invalid WUSS parsing or -1 + */ + private long invalidrnastruc = -1; + + /** + * Updates the _rnasecstr field Determines the positions that base pair and + * the positions of helices based on secondary structure from a Stockholm file + * + * @param RNAannot + */ + private void _updateRnaSecStr(CharSequence RNAannot) + { + try { + _rnasecstr = Rna.GetBasePairs(RNAannot); + bps = Rna.GetModeleBP(RNAannot); + invalidrnastruc=-1; + } + catch (WUSSParseException px) + { + // DEBUG System.out.println(px); + invalidrnastruc=px.getProblemPos(); + } + if (invalidrnastruc > -1) + { + return; + } + Rna.HelixMap(_rnasecstr); + // setRNAStruc(RNAannot); + + if (_rnasecstr != null && _rnasecstr.length > 0) + { + // show all the RNA secondary structure annotation symbols. + isrna = true; + showAllColLabels = true; + scaleColLabel = true; + } + // System.out.println("featuregroup " + _rnasecstr[0].getFeatureGroup()); + } public java.util.Hashtable sequenceMapping; @@ -57,38 +115,39 @@ public class AlignmentAnnotation /** * Score associated with label and description. */ - public double score= Double.NaN; + public double score = Double.NaN; + /** * flag indicating if annotation has a score. */ - public boolean hasScore=false; + public boolean hasScore = false; public GraphLine threshold; // Graphical hints and tips - /** DOCUMENT ME!! */ + /** Can this row be edited by the user ? */ public boolean editable = false; - /** DOCUMENT ME!! */ + /** Indicates if annotation has a graphical symbol track */ public boolean hasIcons; // - /** DOCUMENT ME!! */ + /** Indicates if annotation has a text character label */ public boolean hasText; - /** DOCUMENT ME!! */ + /** is the row visible */ public boolean visible = true; public int graphGroup = -1; - /** DOCUMENT ME!! */ + /** Displayed height of row in pixels */ public int height = 0; public int graph = 0; public int graphHeight = 40; - public boolean padGaps = true; + public boolean padGaps = false; public static final int NO_GRAPH = 0; @@ -98,6 +157,36 @@ public class AlignmentAnnotation public boolean belowAlignment = true; + public SequenceGroup groupRef = null; + + /** + * display every column label, even if there is a row of identical labels + */ + public boolean showAllColLabels = false; + + /** + * scale the column label to fit within the alignment column. + */ + public boolean scaleColLabel = false; + + /** + * centre the column labels relative to the alignment column + */ + public boolean centreColLabels = false; + + private boolean isrna; + + /* + * (non-Javadoc) + * + * @see java.lang.Object#finalize() + */ + protected void finalize() throws Throwable + { + sequenceRef = null; + groupRef = null; + super.finalize(); + } public static int getGraphValueFromString(String string) { @@ -114,17 +203,23 @@ public class AlignmentAnnotation return NO_GRAPH; } } - + // JBPNote: what does this do ? + public void ConcenStru(CharSequence RNAannot) throws WUSSParseException + { + bps = Rna.GetModeleBP(RNAannot); + } /** * Creates a new AlignmentAnnotation object. - * - * @param label DOCUMENT ME! - * @param description DOCUMENT ME! - * @param annotations DOCUMENT ME!about:blank -Loading... + * + * @param label + * short label shown under sequence labels + * @param description + * text displayed on mouseover + * @param annotations + * set of positional annotation elements */ public AlignmentAnnotation(String label, String description, - Annotation[] annotations) + Annotation[] annotations) { // always editable? editable = true; @@ -132,83 +227,278 @@ Loading... this.description = description; this.annotations = annotations; - validateRangeAndDisplay(); + validateRangeAndDisplay(); } + /** + * Checks if annotation labels represent secondary structures + * + */ void areLabelsSecondaryStructure() { boolean nonSSLabel = false; + isrna = false; + StringBuffer rnastring = new StringBuffer(); + + char firstChar = 0; for (int i = 0; i < annotations.length; i++) { if (annotations[i] == null) { - padGaps = false; continue; } - if (annotations[i].secondaryStructure == 'H' || - annotations[i].secondaryStructure == 'E') + if (annotations[i].secondaryStructure == 'H' + || annotations[i].secondaryStructure == 'E') { - hasIcons = true; + hasIcons |= true; + } + else + // Check for RNA secondary structure + { + //System.out.println(annotations[i].secondaryStructure); + if (annotations[i].secondaryStructure == '(' + || annotations[i].secondaryStructure == '[' + || annotations[i].secondaryStructure == '<' + || annotations[i].secondaryStructure == '{' + || annotations[i].secondaryStructure == 'A' + || annotations[i].secondaryStructure == 'B' + || annotations[i].secondaryStructure == 'C' + || annotations[i].secondaryStructure == 'D' + || annotations[i].secondaryStructure == 'E' + || annotations[i].secondaryStructure == 'F' + || annotations[i].secondaryStructure == 'G' + || annotations[i].secondaryStructure == 'H' + || annotations[i].secondaryStructure == 'I' + || annotations[i].secondaryStructure == 'J' + || annotations[i].secondaryStructure == 'K' + || annotations[i].secondaryStructure == 'L' + || annotations[i].secondaryStructure == 'M' + || annotations[i].secondaryStructure == 'N' + || annotations[i].secondaryStructure == 'O' + || annotations[i].secondaryStructure == 'P' + || annotations[i].secondaryStructure == 'Q' + || annotations[i].secondaryStructure == 'R' + || annotations[i].secondaryStructure == 'S' + || annotations[i].secondaryStructure == 'T' + || annotations[i].secondaryStructure == 'U' + || annotations[i].secondaryStructure == 'V' + || annotations[i].secondaryStructure == 'W' + || annotations[i].secondaryStructure == 'X' + || annotations[i].secondaryStructure == 'Y' + || annotations[i].secondaryStructure == 'Z') + { + hasIcons |= true; + isrna |= true; + } } - if(annotations[i].displayCharacter==null) - continue; - + // System.out.println("displaychar " + annotations[i].displayCharacter); - if (annotations[i].displayCharacter.length() == 1 - && !annotations[i].displayCharacter.equals("H") - && !annotations[i].displayCharacter.equals("E") - && !annotations[i].displayCharacter.equals("-") - && !annotations[i].displayCharacter.equals(".")) + if (annotations[i].displayCharacter == null + || annotations[i].displayCharacter.length() == 0) + { + rnastring.append('.'); + continue; + } + if (annotations[i].displayCharacter.length() == 1) + { + firstChar = annotations[i].displayCharacter.charAt(0); + // check to see if it looks like a sequence or is secondary structure + // labelling. + if (annotations[i].secondaryStructure != ' ' + && !hasIcons + && + // Uncomment to only catch case where + // displayCharacter==secondary + // Structure + // to correctly redisplay SS annotation imported from Stockholm, + // exported to JalviewXML and read back in again. + // && + // annotations[i].displayCharacter.charAt(0)==annotations[i].secondaryStructure + firstChar != ' ' + && firstChar != '$' + && firstChar != '�' // JBPNote should explicitly express as unicode number to avoid source code translation problems + && firstChar != '(' + && firstChar != '[' + && firstChar != '>' + && firstChar != '{' + && firstChar != 'A' + && firstChar != 'B' + && firstChar != 'C' + && firstChar != 'D' + && firstChar != 'E' + && firstChar != 'F' + && firstChar != 'G' + && firstChar != 'H' + && firstChar != 'I' + && firstChar != 'J' + && firstChar != 'K' + && firstChar != 'L' + && firstChar != 'M' + && firstChar != 'N' + && firstChar != 'O' + && firstChar != 'P' + && firstChar != 'Q' + && firstChar != 'R' + && firstChar != 'S' + && firstChar != 'T' + && firstChar != 'U' + && firstChar != 'V' + && firstChar != 'W' + && firstChar != 'X' + && firstChar != 'Y' + && firstChar != 'Z' + && firstChar != '-' + && firstChar < jalview.schemes.ResidueProperties.aaIndex.length) { - if (jalview.schemes.ResidueProperties.aaIndex - [annotations[i].displayCharacter.charAt(0)] < 23) + if (jalview.schemes.ResidueProperties.aaIndex[firstChar] < 23) // TODO: + // parameterise + // to + // gap + // symbol + // number { nonSSLabel = true; } } + } + else + { + rnastring.append(annotations[i].displayCharacter.charAt(1)); + } - if (annotations[i].displayCharacter.length() > 0) - { - hasText = true; - } - else - padGaps = false; + if (annotations[i].displayCharacter.length() > 0) + { + hasText = true; } + } if (nonSSLabel) { hasIcons = false; for (int j = 0; j < annotations.length; j++) { - if (annotations[j] != null && annotations[j].secondaryStructure != ' ') + if (annotations[j] != null + && annotations[j].secondaryStructure != ' ') { - annotations[j].displayCharacter - = String.valueOf(annotations[j].secondaryStructure); + annotations[j].displayCharacter = String + .valueOf(annotations[j].secondaryStructure); annotations[j].secondaryStructure = ' '; } } } + else + { + if (isrna) + { + _updateRnaSecStr(new AnnotCharSequence()); + } + } annotationId = this.hashCode() + ""; } + + /** + * flyweight access to positions in the alignment annotation row for RNA + * processing + * + * @author jimp + * + */ + private class AnnotCharSequence implements CharSequence + { + int offset = 0; + + int max = 0; + + public AnnotCharSequence() + { + this(0, annotations.length); + } + + public AnnotCharSequence(int start, int end) + { + offset = start; + max = end; + } + + @Override + public CharSequence subSequence(int start, int end) + { + return new AnnotCharSequence(offset + start, offset + end); + } + + @Override + public int length() + { + return max - offset; + } + + @Override + public char charAt(int index) + { + String dc; + return ((index + offset < 0) || (index + offset) >= max + || annotations[index + offset] == null || (dc = annotations[index + + offset].displayCharacter.trim()).length() < 1) ? '.' : dc + .charAt(0); + } + + public String toString() + { + char[] string = new char[max - offset]; + int mx = annotations.length; + + for (int i = offset; i < mx; i++) + { + String dc; + string[i] = (annotations[i] == null || (dc = annotations[i].displayCharacter + .trim()).length() < 1) ? '.' : dc.charAt(0); + } + return new String(string); + } + }; + + private long _lastrnaannot = -1; + + public String getRNAStruc() + { + if (isrna) + { + String rnastruc = new AnnotCharSequence().toString(); + if (_lastrnaannot != rnastruc.hashCode()) + { + // ensure rna structure contacts are up to date + _lastrnaannot = rnastruc.hashCode(); + _updateRnaSecStr(rnastruc); + } + return rnastruc; + } + return null; + } + /** * Creates a new AlignmentAnnotation object. - * - * @param label DOCUMENT ME! - * @param description DOCUMENT ME! - * @param annotations DOCUMENT ME! - * @param min DOCUMENT ME! - * @param max DOCUMENT ME! - * @param winLength DOCUMENT ME! + * + * @param label + * DOCUMENT ME! + * @param description + * DOCUMENT ME! + * @param annotations + * DOCUMENT ME! + * @param min + * DOCUMENT ME! + * @param max + * DOCUMENT ME! + * @param winLength + * DOCUMENT ME! */ public AlignmentAnnotation(String label, String description, - Annotation[] annotations, float min, float max, - int graphType) + Annotation[] annotations, float min, float max, int graphType) { // graphs are not editable - editable = graphType==0; + editable = graphType == 0; this.label = label; this.description = description; @@ -218,18 +508,17 @@ Loading... graphMax = max; validateRangeAndDisplay(); } + /** - * checks graphMin and graphMax, - * secondary structure symbols, - * sets graphType appropriately, - * sets null labels to the empty string - * if appropriate. + * checks graphMin and graphMax, secondary structure symbols, sets graphType + * appropriately, sets null labels to the empty string if appropriate. */ - private void validateRangeAndDisplay() { + public void validateRangeAndDisplay() + { - if (annotations==null) + if (annotations == null) { - visible=false; // try to prevent renderer from displaying. + visible = false; // try to prevent renderer from displaying. return; // this is a non-annotation row annotation - ie a sequence score. } @@ -237,7 +526,7 @@ Loading... float min = graphMin; float max = graphMax; boolean drawValues = true; - + _linecolour = null; if (min == max) { min = 999999999; @@ -248,9 +537,8 @@ Loading... continue; } - if (drawValues - && annotations[i].displayCharacter!=null - && annotations[i].displayCharacter.length() > 1) + if (drawValues && annotations[i].displayCharacter != null + && annotations[i].displayCharacter.length() > 1) { drawValues = false; } @@ -264,6 +552,22 @@ Loading... { min = annotations[i].value; } + if (_linecolour == null && annotations[i].colour != null) + { + _linecolour = annotations[i].colour; + } + } + // ensure zero is origin for min/max ranges on only one side of zero + if (min > 0) + { + min = 0; + } + else + { + if (max < 0) + { + max = 0; + } } } @@ -278,15 +582,16 @@ Loading... { if (annotations[i] != null) { - annotations[i].displayCharacter = ""; + annotations[i].displayCharacter = "X"; } } } } /** - * Copy constructor - * creates a new independent annotation row with the same associated sequenceRef + * Copy constructor creates a new independent annotation row with the same + * associated sequenceRef + * * @param annotation */ public AlignmentAnnotation(AlignmentAnnotation annotation) @@ -299,6 +604,7 @@ Loading... this.graph = annotation.graph; this.graphHeight = annotation.graphHeight; this.graphGroup = annotation.graphGroup; + this.groupRef = annotation.groupRef; this.editable = annotation.editable; this.autoCalculated = annotation.autoCalculated; this.hasIcons = annotation.hasIcons; @@ -307,77 +613,115 @@ Loading... this.label = annotation.label; this.padGaps = annotation.padGaps; this.visible = annotation.visible; + this.centreColLabels = annotation.centreColLabels; + this.scaleColLabel = annotation.scaleColLabel; + this.showAllColLabels = annotation.showAllColLabels; + this.calcId = annotation.calcId; if (this.hasScore = annotation.hasScore) { this.score = annotation.score; } - if (threshold!=null) { + if (annotation.threshold != null) + { threshold = new GraphLine(annotation.threshold); } - if (annotation.annotations!=null) { + if (annotation.annotations != null) + { Annotation[] ann = annotation.annotations; this.annotations = new Annotation[ann.length]; - for (int i=0; i= annotations.length) + startRes = annotations.length - 1; + if (endRes >= annotations.length) + endRes = annotations.length - 1; + if (annotations == null) return; - Annotation[] temp = new Annotation[endRes-startRes+1]; - if (startRes=spos && pos.intValue()<=epos) + if (pos.intValue() >= spos && pos.intValue() <= epos) { newmapping.put(pos, sequenceMapping.get(pos)); } @@ -386,33 +730,37 @@ Loading... sequenceMapping = newmapping; } } - annotations=temp; + annotations = temp; } + /** * set the annotation row to be at least length Annotations - * @param length minimum number of columns required in the annotation row + * + * @param length + * minimum number of columns required in the annotation row * @return false if the annotation row is greater than length */ - public boolean padAnnotation(int length) { - if (annotations==null) + public boolean padAnnotation(int length) + { + if (annotations == null) { - annotations = new Annotation[length]; - return true; + return true; // annotation row is correct - null == not visible and + // undefined length } - if (annotations.lengthlength; + return annotations.length > length; } /** * DOCUMENT ME! - * + * * @return DOCUMENT ME! */ public String toString() @@ -439,8 +787,8 @@ Loading... buffer.append(", "); } - - if (label.equals("Consensus")) + // TODO: remove disgusting hack for 'special' treatment of consensus line. + if (label.indexOf("Consensus") == 0) { buffer.append("\n"); @@ -469,22 +817,24 @@ Loading... } /** - * Attach the annotation to seqRef, starting from startRes position. If alreadyMapped is true then the indices of the annotation[] array are sequence positions rather than alignment column positions. + * Attach the annotation to seqRef, starting from startRes position. If + * alreadyMapped is true then the indices of the annotation[] array are + * sequence positions rather than alignment column positions. + * * @param seqRef * @param startRes * @param alreadyMapped */ - public void createSequenceMapping(SequenceI seqRef, - int startRes, - boolean alreadyMapped) + public void createSequenceMapping(SequenceI seqRef, int startRes, + boolean alreadyMapped) { if (seqRef == null) { return; } - sequenceRef=seqRef; - if (annotations==null) + sequenceRef = seqRef; + if (annotations == null) { return; } @@ -513,10 +863,10 @@ Loading... public void adjustForAlignment() { - if (sequenceRef==null) + if (sequenceRef == null) return; - if (annotations==null) + if (annotations == null) { return; } @@ -525,7 +875,7 @@ Loading... if (aSize == 0) { - //Its been deleted + // Its been deleted return; } @@ -546,31 +896,45 @@ Loading... annotations = temp; } + /** - * remove any null entries in annotation row and return the - * number of non-null annotation elements. + * remove any null entries in annotation row and return the number of non-null + * annotation elements. + * * @return */ - private int compactAnnotationArray() { - int j=0; - for (int i=0;i