X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;ds=sidebyside;f=src%2Fjalview%2Fio%2FAnnotationFile.java;h=00b7ead85bde0ac4d6b5c104d33c3752c971904f;hb=96f41aa46d3e568f8734834ab36fdc9ad8e54ac1;hp=8ec234d12128745ff90f873ed3fa50611d9448c0;hpb=f22d8c3bfc1892fe88b4aad3f6381e5018f175f1;p=jalview.git diff --git a/src/jalview/io/AnnotationFile.java b/src/jalview/io/AnnotationFile.java index 8ec234d..00b7ead 100755 --- a/src/jalview/io/AnnotationFile.java +++ b/src/jalview/io/AnnotationFile.java @@ -1,51 +1,99 @@ /* - * Jalview - A Sequence Alignment Editor and Viewer - * Copyright (C) 2006 AM Waterhouse, J Procter, G Barton, M Clamp, S Searle - * + * Jalview - A Sequence Alignment Editor and Viewer (Development Version 2.4.1) + * Copyright (C) 2009 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 */ - package jalview.io; import java.io.*; -import jalview.datamodel.*; +import java.net.*; import java.util.*; -import java.net.URL; -import jalview.schemes.*; -import jalview.analysis.Conservation; +import jalview.analysis.*; +import jalview.datamodel.*; +import jalview.schemes.*; public class AnnotationFile { - StringBuffer text = new StringBuffer( - "JALVIEW_ANNOTATION\n" - +"# Created: " - +new java.util.Date()+"\n\n"); + StringBuffer text = new StringBuffer("JALVIEW_ANNOTATION\n" + + "# Created: " + new java.util.Date() + "\n\n"); + + /** + * convenience method for pre-2.4 feature files which have no view, hidden + * columns or hidden row keywords. + * + * @param annotations + * @param groups + * @param properties + * @return feature file as a string. + */ + public String printAnnotations(AlignmentAnnotation[] annotations, + Vector groups, Hashtable properties) + { + return printAnnotations(annotations, groups, properties, null); + + } - public String printAnnotations(AlignmentAnnotation [] annotations, - Vector groups) + /** + * hold all the information about a particular view definition read from or + * written out in an annotations file. + */ + public class ViewDef + { + public String viewname; + + public HiddenSequences hidseqs; + + public ColumnSelection hiddencols; + + public Vector visibleGroups; + + public Hashtable hiddenRepSeqs; + + public ViewDef(String viewname, HiddenSequences hidseqs, + ColumnSelection hiddencols, Hashtable hiddenRepSeqs) + { + this.viewname = viewname; + this.hidseqs = hidseqs; + this.hiddencols = hiddencols; + this.hiddenRepSeqs = hiddenRepSeqs; + } + } + /** + * Prepare an annotation file given a set of annotations, groups, alignment properties and views. + * @param annotations + * @param groups + * @param properties + * @param views + * @return annotation file + */ + public String printAnnotations(AlignmentAnnotation[] annotations, + Vector groups, Hashtable properties, ViewDef[] views) { - if(annotations!=null) + // TODO: resolve views issue : annotationFile could contain visible region, or full data + hidden region specifications for a view. + if (annotations != null) { + boolean oneColour = true; AlignmentAnnotation row; String comma; - SequenceI seqref = null; + SequenceI refSeq = null; StringBuffer colours = new StringBuffer(); StringBuffer graphLine = new StringBuffer(); - + StringBuffer rowprops = new StringBuffer(); Hashtable graphGroup = new Hashtable(); java.awt.Color color; @@ -54,100 +102,195 @@ public class AnnotationFile { row = annotations[i]; - if (!row.visible) + if (!row.visible && !row.hasScore()) + { continue; + } color = null; + oneColour = true; if (row.sequenceRef == null) { - if (seqref != null) + if (refSeq != null) + { text.append("\nSEQUENCE_REF\tALIGNMENT\n"); + } - seqref = null; + refSeq = null; } - else if (seqref == null || seqref != row.sequenceRef) + else if (refSeq == null || refSeq != row.sequenceRef) { - seqref = row.sequenceRef; - text.append("\nSEQUENCE_REF\t" + seqref.getName() + "\n"); + refSeq = row.sequenceRef; + text.append("\nSEQUENCE_REF\t" + refSeq.getName() + "\n"); + } + boolean hasGlyphs = false, hasLabels = false, hasValues = false, hasText = false; + // lookahead to check what the annotation row object actually contains. + for (int j = 0; row.annotations != null + && j < row.annotations.length + && (!hasGlyphs || !hasLabels || !hasValues); j++) + { + if (row.annotations[j] != null) + { + hasLabels |= (row.annotations[j].displayCharacter != null + && row.annotations[j].displayCharacter.length() > 0 && !row.annotations[j].displayCharacter + .equals(" ")); + hasGlyphs |= (row.annotations[j].secondaryStructure != 0 && row.annotations[j].secondaryStructure != ' '); + hasValues |= (row.annotations[j].value != Float.NaN); // NaNs can't + // be + // rendered.. + hasText |= (row.annotations[j].description != null && row.annotations[j].description + .length() > 0); + } } if (row.graph == AlignmentAnnotation.NO_GRAPH) { text.append("NO_GRAPH\t"); + hasValues = false; // only secondary structure + // hasLabels = false; // and annotation description string. } else { if (row.graph == AlignmentAnnotation.BAR_GRAPH) + { text.append("BAR_GRAPH\t"); + hasGlyphs = false; // no secondary structure + + } else if (row.graph == AlignmentAnnotation.LINE_GRAPH) + { + hasGlyphs = false; // no secondary structure text.append("LINE_GRAPH\t"); + } if (row.getThreshold() != null) - graphLine.append("GRAPHLINE\t" - + row.label + "\t" - + row.getThreshold().value + "\t" - + row.getThreshold().label + "\t" - + jalview.util.Format.getHexString( - row.getThreshold().colour) + "\n" - ); + { + graphLine + .append("GRAPHLINE\t" + + row.label + + "\t" + + row.getThreshold().value + + "\t" + + row.getThreshold().label + + "\t" + + jalview.util.Format.getHexString(row + .getThreshold().colour) + "\n"); + } if (row.graphGroup > -1) { String key = String.valueOf(row.graphGroup); if (graphGroup.containsKey(key)) - graphGroup.put(key, graphGroup.get(key) - + "\t" + row.label); + { + graphGroup.put(key, graphGroup.get(key) + "\t" + row.label); + } else + { graphGroup.put(key, row.label); + } } } text.append(row.label + "\t"); if (row.description != null) + { text.append(row.description + "\t"); - - for (int j = 0; j < row.annotations.length; j++) + } + for (int j = 0; row.annotations != null + && j < row.annotations.length; j++) { - if (seqref != null && - jalview.util.Comparison.isGap(seqref.getCharAt(j))) + if (refSeq != null + && jalview.util.Comparison.isGap(refSeq.getCharAt(j))) + { continue; + } if (row.annotations[j] != null) { comma = ""; - if (row.annotations[j].secondaryStructure != ' ') + if (hasGlyphs && row.annotations[j].secondaryStructure != ' ') { + text.append(comma + row.annotations[j].secondaryStructure); comma = ","; } - if (row.annotations[j].displayCharacter.length() > 0 - && !row.annotations[j].displayCharacter.equals(" ")) + if (hasValues) { - text.append(comma + row.annotations[j].displayCharacter); + if (row.annotations[j].value != Float.NaN) + { + text.append(comma + row.annotations[j].value); + } + else + { + System.err.println("Skipping NaN - not valid value."); + text.append(comma + 0f);// row.annotations[j].value); + } comma = ","; } + if (hasLabels) + { + // TODO: labels are emitted after values for bar graphs. + if // empty labels are allowed, so + (row.annotations[j].displayCharacter != null + && row.annotations[j].displayCharacter.length() > 0 + && !row.annotations[j].displayCharacter.equals(" ")) + { + text.append(comma + row.annotations[j].displayCharacter); + comma = ","; + } + } + if (hasText) + { + if (row.annotations[j].description != null + && row.annotations[j].description.length() > 0 + && !row.annotations[j].description + .equals(row.annotations[j].displayCharacter)) + { + text.append(comma + row.annotations[j].description); + comma = ","; + } + } + if (color != null && !color.equals(row.annotations[j].colour)) + { + oneColour = false; + } - if (row.annotations[j] != null) + color = row.annotations[j].colour; + + if (row.annotations[j].colour != null + && row.annotations[j].colour != java.awt.Color.black) { - color = row.annotations[j].colour; - if (row.annotations[j].value != 0f) - text.append(comma + row.annotations[j].value); + text.append(comma + + "[" + + jalview.util.Format + .getHexString(row.annotations[j].colour) + + "]"); + comma = ","; } } text.append("|"); } + if (row.hasScore()) + text.append("\t" + row.score); + text.append("\n"); - if (color != null && color != java.awt.Color.black) + if (color != null && color != java.awt.Color.black && oneColour) { - colours.append("COLOUR\t" - + row.label + "\t" - + jalview.util.Format.getHexString(color) + "\n"); + colours.append("COLOUR\t" + row.label + "\t" + + jalview.util.Format.getHexString(color) + "\n"); + } + if (row.scaleColLabel || row.showAllColLabels || row.centreColLabels) + { + rowprops.append("ROWPROPERTIES\t"+row.label); + rowprops.append("\tscaletofit="+row.scaleColLabel); + rowprops.append("\tshowalllabs="+row.showAllColLabels); + rowprops.append("\tcentrelabs="+row.centreColLabels); + rowprops.append("\n"); } - } text.append("\n"); @@ -163,71 +306,127 @@ public class AnnotationFile text.append(en.nextElement() + "\n"); } } + text.append(rowprops.toString()); } - - if(groups!=null) + if (groups != null) { printGroups(groups); } + if (properties != null) + { + text.append("\n\nALIGNMENT"); + Enumeration en = properties.keys(); + while (en.hasMoreElements()) + { + String key = en.nextElement().toString(); + text.append("\t" + key + "=" + properties.get(key)); + } + // TODO: output alignment visualization settings here if required + + } + return text.toString(); } - public void printGroups(Vector sequenceGroups) { SequenceGroup sg; + SequenceI seqrep = null; for (int i = 0; i < sequenceGroups.size(); i++) { sg = (SequenceGroup) sequenceGroups.elementAt(i); - text.append("SEQUENCE_GROUP\t" - + sg.getName() + "\t" - + (sg.getStartRes()+1) + "\t" - + (sg.getEndRes()+1) + "\t" + "-1\t"); - for(int s=0; s -1) { @@ -281,196 +484,325 @@ public class AnnotationFile } } - if(!jvAnnotationFile) + if (!jvAnnotationFile) { in.close(); return false; } - while ( (line = in.readLine()) != null) + while ((line = in.readLine()) != null) { - if(line.indexOf("#")==0 - || line.indexOf("JALVIEW_ANNOTATION")>-1 - || line.length()==0) + if (line.indexOf("#") == 0 + || line.indexOf("JALVIEW_ANNOTATION") > -1 + || line.length() == 0) + { continue; + } st = new StringTokenizer(line, "\t"); token = st.nextToken(); - if(token.equalsIgnoreCase("COLOUR")) + if (token.equalsIgnoreCase("COLOUR")) { + // TODO: use graduated colour def'n here too colourAnnotations(al, st.nextToken(), st.nextToken()); continue; } - else if(token.equalsIgnoreCase("COMBINE") ) + else if (token.equalsIgnoreCase("COMBINE")) { combineAnnotations(al, st); continue; } - + else if (token.equalsIgnoreCase("ROWPROPERTIES")) + { + addRowProperties(al, st); + continue; + } else if (token.equalsIgnoreCase("GRAPHLINE")) { addLine(al, st); continue; } - - else if(token.equalsIgnoreCase("SEQUENCE_REF") ) + else if (token.equalsIgnoreCase("SEQUENCE_REF")) { - refSeq = al.findName(st.nextToken()); - try{ - refSeqIndex = Integer.parseInt(st.nextToken()); - if(refSeqIndex<1) + if (st.hasMoreTokens()) + { + refSeq = al.findName(refSeqId = st.nextToken()); + if (refSeq == null) + { + refSeqId = null; + } + try + { + refSeqIndex = Integer.parseInt(st.nextToken()); + if (refSeqIndex < 1) + { + refSeqIndex = 1; + System.out + .println("WARNING: SEQUENCE_REF index must be > 0 in AnnotationFile"); + } + } catch (Exception ex) { refSeqIndex = 1; - System.out.println("WARNING: SEQUENCE_REF index must be > 0 in AnnotationFile"); } } - catch(Exception ex) + else { - refSeqIndex = 1; + refSeq = null; + refSeqId = null; } - continue; } - else if(token.equalsIgnoreCase("SEQUENCE_GROUP")) + else if (token.equalsIgnoreCase("SEQUENCE_GROUP")) { addGroup(al, st); continue; } - else if(token.equalsIgnoreCase("PROPERTIES")) + else if (token.equalsIgnoreCase("PROPERTIES")) { addProperties(al, st); continue; } + else if (token.equalsIgnoreCase("BELOW_ALIGNMENT")) + { + setBelowAlignment(al, st); + continue; + } + else if (token.equalsIgnoreCase("ALIGNMENT")) + { + addAlignmentDetails(al, st); + continue; + } graphStyle = AlignmentAnnotation.getGraphValueFromString(token); label = st.nextToken(); - if(st.countTokens()>1) - description = st.nextToken(); - else - description = null; - - line = st.nextToken(); - - st = new StringTokenizer(line, "|", true); - annotations = new Annotation[alWidth]; - index = 0; - boolean emptyColumn = true; - + annotations = new Annotation[alWidth]; + description = null; + float score = Float.NaN; - while (st.hasMoreElements() && index -1 && j > -1) + { + UserColourScheme ucs = new UserColourScheme(); + + colour = ucs.getColourFromString(string.substring(i + 1, j)); + if ( i>0 && string.charAt(i-1)==',') + { + // clip the preceding comma as well + i--; + } + string = string.substring(0, i) + string.substring(j + 1); + } + + StringTokenizer st = new StringTokenizer(string, ",", true); String token; - while(st.hasMoreTokens()) + boolean seenContent=false; + while (st.hasMoreTokens()) { token = st.nextToken().trim(); - if(token.length()==0) + if (token.equals(",")) + { + if (!seenContent && parsedValue && !dcset) + { + // allow the value below the bar/line to be empty + dcset = true; + displayChar = " "; + } + seenContent = false; continue; + } else { + seenContent = true; + } - if(!parsedValue) + if (!parsedValue) { - try{ - displayChar = token; + try + { + displayChar = token; + // foo value = new Float(token).floatValue(); parsedValue = true; continue; - }catch(NumberFormatException ex){} + } catch (NumberFormatException ex) + { + } } - if(token.equals("H") || token.equals("E")) + if (hasSymbols && (token.equals("H") || token.equals("E"))) { // Either this character represents a helix or sheet // or an integer which can be displayed ss = token.charAt(0); - if(displayChar.equals(token.substring(0,1))) + if (displayChar.equals(token.substring(0, 1))) + { displayChar = ""; + } } - else if(desc==null) + else if (desc == null) + { desc = token; + } } - - if(desc == null) - desc = value+""; - - if(displayChar.length()>1 && desc.length()==1) +// if (!dcset && string.charAt(string.length() - 1) == ',') +// { +// displayChar = " "; // empty display char symbol. +// } + if (displayChar != null && displayChar.length() > 1 && desc != null + && desc.length() == 1) { String tmp = displayChar; displayChar = desc; desc = tmp; } + /* + * In principle, this code will ensure that the Annotation element generated + * is renderable by any of the applet or application rendering code but + * instead we check for null strings when the display character is rendered. + * if (displayChar==null) { displayChar=""; } + */ + Annotation anot = new Annotation(displayChar, desc, ss, value); + + anot.colour = colour; - return new Annotation(displayChar, desc, ss, value); + return anot; } void colourAnnotations(AlignmentI al, String label, String colour) { UserColourScheme ucs = new UserColourScheme(colour); Annotation[] annotations; - for(int i=0; i 0 && !rng.startsWith("*")) + { + sg.setStartRes(Integer.parseInt(rng) - 1); + } + else + { + sg.setStartRes(0); + } + rng = st.nextToken(); + if (rng.length() > 0 && !rng.startsWith("*")) + { + sg.setEndRes(Integer.parseInt(rng) - 1); + } + else + { + sg.setEndRes(al.getWidth() - 1); + } + } catch (Exception e) + { + System.err + .println("Couldn't parse Group Start or End Field as '*' or a valid column or sequence index: '" + + rng + "' - assuming alignment width for group."); + // assume group is full width + sg.setStartRes(0); + sg.setEndRes(al.getWidth() - 1); + } String index = st.nextToken(); - if(index.equals("-1")) + if (index.equals("-1")) { while (st.hasMoreElements()) { @@ -571,63 +933,111 @@ public class AnnotationFile int end = (Integer.parseInt(st3.nextToken())); if (end > start) + { for (int i = start; i <= end; i++) - sg.addSequence(al.getSequenceAt(i-1), false); + { + sg.addSequence(al.getSequenceAt(i - 1), false); + } + } } else { - sg.addSequence(al.getSequenceAt(Integer.parseInt(tmp)-1 ), false); + sg + .addSequence(al.getSequenceAt(Integer.parseInt(tmp) - 1), + false); } } } - if(refSeq!=null) + if (refSeq != null) { - sg.setStartRes( refSeq.findIndex( sg.getStartRes()+1 )-1 ); - sg.setEndRes( refSeq.findIndex( sg.getEndRes() +1) -1 ); + sg.setStartRes(refSeq.findIndex(sg.getStartRes() + 1) - 1); + sg.setEndRes(refSeq.findIndex(sg.getEndRes() + 1) - 1); + sg.setSeqrep(refSeq); } - if(sg.getSize()>0) + if (sg.getSize() > 0) { al.addGroup(sg); } } + void addRowProperties(AlignmentI al, StringTokenizer st) + { + String label = st.nextToken(),keyValue,key,value; + boolean scaletofit=false,centerlab=false,showalllabs=false; + while (st.hasMoreTokens()) { + keyValue=st.nextToken(); + key = keyValue.substring(0, keyValue.indexOf("=")); + value = keyValue.substring(keyValue.indexOf("=") + 1); + if (key.equalsIgnoreCase("scaletofit")) { + scaletofit = Boolean.valueOf(value).booleanValue(); + } + if (key.equalsIgnoreCase("showalllabs")) { + showalllabs = Boolean.valueOf(value).booleanValue(); + } + if (key.equalsIgnoreCase("centrelabs")) { + centerlab = Boolean.valueOf(value).booleanValue(); + } + AlignmentAnnotation[] alr = al.getAlignmentAnnotation(); + for (int i = 0; i < alr.length; i++) + { + if (alr[i].label.equalsIgnoreCase(label)) + { + alr[i].centreColLabels = centerlab; + alr[i].scaleColLabel = scaletofit; + alr[i].showAllColLabels = showalllabs; + } + } + } + } void addProperties(AlignmentI al, StringTokenizer st) { - //So far we have only added groups to the annotationHash, - //the idea is in the future properties can be added to - //alignments, other annotations etc - if(al.getGroups()==null) + // So far we have only added groups to the annotationHash, + // the idea is in the future properties can be added to + // alignments, other annotations etc + if (al.getGroups() == null) + { return; + } SequenceGroup sg = null; String name = st.nextToken(); Vector groups = al.getGroups(); - for(int i=0; i p) + { + sp.append(atos.substring(p, cp + 1)); + } + else + { + sp.append(atos.substring(p)); + sp.append("\n"); + } + p = cp + 1; + } while (p > 0); + } + return sp.toString(); + } }