X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Fappletgui%2FAlignViewport.java;h=43831a247ac8504410961eab86aff4710372c95b;hb=cbd8cafe5f1d302cf490d34a96f8231659afd402;hp=18a42069865b4d953020d306ae54a54484cc9afd;hpb=174230b4233d9ce80f94527768d2cd2f76da11ab;p=jalview.git diff --git a/src/jalview/appletgui/AlignViewport.java b/src/jalview/appletgui/AlignViewport.java index 18a4206..43831a2 100755 --- a/src/jalview/appletgui/AlignViewport.java +++ b/src/jalview/appletgui/AlignViewport.java @@ -1,6 +1,6 @@ /* * Jalview - A Sequence Alignment Editor and Viewer - * Copyright (C) 2006 AM Waterhouse, J Procter, G Barton, M Clamp, S Searle + * 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 @@ -36,7 +36,6 @@ public class AlignViewport int startSeq; int endSeq; - boolean cursorMode = false; boolean showJVSuffix = true; @@ -50,6 +49,7 @@ public class AlignViewport boolean showConservation = true; boolean showQuality = true; boolean showConsensus = true; + boolean upperCasebold = false; boolean colourAppliesToAllGroups = true; ColourSchemeI globalColourScheme = null; @@ -78,15 +78,14 @@ public class AlignViewport boolean scaleRightWrapped = true; // The following vector holds the features which are - // currently visible, in the correct order or rendering - Hashtable featuresDisplayed; + // currently visible, in the correct order or rendering + public Hashtable featuresDisplayed; boolean hasHiddenColumns = false; boolean hasHiddenRows = false; boolean showHiddenMarkers = true; - - public Vector vconsensus; + public Hashtable[] hconsensus; AlignmentAnnotation consensus; AlignmentAnnotation conservation; AlignmentAnnotation quality; @@ -95,12 +94,24 @@ public class AlignViewport public int ConsPercGaps = 25; // JBPNote : This should be a scalable property! - private java.beans.PropertyChangeSupport changeSupport = new java.beans.PropertyChangeSupport(this); + private java.beans.PropertyChangeSupport changeSupport = new java.beans. + PropertyChangeSupport(this); boolean ignoreGapsInConsensusCalculation = false; jalview.bin.JalviewLite applet; + Hashtable sequenceColours; + + boolean MAC = false; + + Stack historyList = new Stack(); + Stack redoList = new Stack(); + + String sequenceSetID; + + Hashtable hiddenRepSequences; + public AlignViewport(AlignmentI al, JalviewLite applet) { this.applet = applet; @@ -111,6 +122,11 @@ public class AlignViewport this.endSeq = al.getHeight() - 1; setFont(font); + if (System.getProperty("os.name").startsWith("Mac")) + { + MAC = true; + } + if (applet != null) { String param = applet.getParameter("showFullId"); @@ -142,40 +158,90 @@ public class AlignViewport { showConsensus = Boolean.valueOf(param).booleanValue(); } - } - // We must set conservation and consensus before setting colour, - // as Blosum and Clustal require this to be done - updateConservation(); - updateConsensus(); + param = applet.getParameter("upperCase"); + if (param != null) + { + if (param.equalsIgnoreCase("bold")) + { + upperCasebold = true; + } + } + + } if (applet != null) { String colour = applet.getParameter("defaultColour"); - if(colour == null) + if (colour == null) { colour = applet.getParameter("userDefinedColour"); - if(colour !=null) + if (colour != null) + { colour = "User Defined"; + } } - if(colour != null) + if (colour != null) { globalColourScheme = ColourSchemeProperty.getColour(alignment, colour); if (globalColourScheme != null) { - globalColourScheme.setConsensus(vconsensus); + globalColourScheme.setConsensus(hconsensus); } } - if(applet.getParameter("userDefinedColour")!=null) + if (applet.getParameter("userDefinedColour") != null) { - ((UserColourScheme)globalColourScheme).parseAppletParameter( + ( (UserColourScheme) globalColourScheme).parseAppletParameter( applet.getParameter("userDefinedColour")); } + if (hconsensus == null) + { + if (!alignment.isNucleotide()) + { + conservation = new AlignmentAnnotation("Conservation", + "Conservation of total alignment less than " + + ConsPercGaps + "% gaps", + new Annotation[1], 0f, + 11f, + AlignmentAnnotation.BAR_GRAPH); + conservation.hasText = true; + conservation.autoCalculated = true; + + if (showConservation) + { + alignment.addAnnotation(conservation); + } + if (showQuality) + { + quality = new AlignmentAnnotation("Quality", + "Alignment Quality based on Blosum62 scores", + new Annotation[1], + 0f, + 11f, + AlignmentAnnotation.BAR_GRAPH); + quality.hasText = true; + quality.autoCalculated = true; + + alignment.addAnnotation(quality); + } + } + + consensus = new AlignmentAnnotation("Consensus", "PID", + new Annotation[1], 0f, 100f, + AlignmentAnnotation.BAR_GRAPH); + consensus.hasText = true; + consensus.autoCalculated = true; + + if (showConsensus) + { + alignment.addAnnotation(consensus); + } + } } } @@ -189,168 +255,322 @@ public class AlignViewport return showSequenceFeatures; } - - public void updateConservation() + class ConservationThread + extends Thread { - if(alignment.isNucleotide()) - return; + AlignmentPanel ap; + public ConservationThread(AlignmentPanel ap) + { + this.ap = ap; + } - Conservation cons = new jalview.analysis.Conservation("All", - jalview.schemes.ResidueProperties.propHash, 3, - alignment.getSequences(), 0, - alignment.getWidth() - 1); - cons.calculate(); - cons.verdict(false, ConsPercGaps); - cons.findQuality(); - int alWidth = alignment.getWidth(); - Annotation[] annotations = new Annotation[alWidth]; - Annotation[] qannotations = new Annotation[alWidth]; - String sequence = cons.getConsSequence().getSequence(); - float minR, minG, minB, maxR, maxG, maxB; - minR = 0.3f; - minG = 0.0f; - minB = 0f; - maxR = 1.0f - minR; - maxG = 0.9f - minG; - maxB = 0f - minB; // scalable range for colouring both Conservation and Quality - float min = 0f; - float max = 11f; - float qmin = cons.qualityRange[0].floatValue(); - float qmax = cons.qualityRange[1].floatValue(); - - for (int i = 0; i < alWidth; i++) + public void run() { - float value = 0; try { - value = Integer.parseInt(sequence.charAt(i) + ""); - } - catch (Exception ex) - { - if (sequence.charAt(i) == '*') + updatingConservation = true; + + while (UPDATING_CONSERVATION) { - value = 11; + try + { + if (ap != null) + { + ap.paintAlignment(true); + } + Thread.sleep(200); + } + catch (Exception ex) + { + ex.printStackTrace(); + } } - if (sequence.charAt(i) == '+') + + UPDATING_CONSERVATION = true; + + int alWidth = alignment.getWidth(); + if (alWidth < 0) { - value = 10; + return; } - } - float vprop = value - min; - vprop /= max; - - annotations[i] = new Annotation(sequence.charAt(i) + "", - "", ' ', value, - new Color(minR + maxR * vprop, - minG + maxG * vprop, - minB + maxB * vprop)); - // Quality calc - value = ( (Double) cons.quality.elementAt(i)).floatValue(); - vprop = value - qmin; - vprop /= qmax; - qannotations[i] = new Annotation(" ", - String.valueOf(value), ' ', value, - new - Color(minR + maxR * vprop, - minG + maxG * vprop, - minB + maxB * vprop)); - } - if (conservation == null) - { - conservation = new AlignmentAnnotation("Conservation", - "Conservation of total alignment less than " + - ConsPercGaps + "% gaps", - annotations, - 0f, // cons.qualityRange[0].floatValue(), - 11f, // cons.qualityRange[1].floatValue() - AlignmentAnnotation.BAR_GRAPH); - if (showConservation) + Conservation cons = new jalview.analysis.Conservation("All", + jalview.schemes.ResidueProperties.propHash, 3, + alignment.getSequences(), 0, alWidth - 1); + + cons.calculate(); + cons.verdict(false, ConsPercGaps); + + if (quality != null) + { + cons.findQuality(); + } + + char[] sequence = cons.getConsSequence().getSequence(); + float minR; + float minG; + float minB; + float maxR; + float maxG; + float maxB; + minR = 0.3f; + minG = 0.0f; + minB = 0f; + maxR = 1.0f - minR; + maxG = 0.9f - minG; + maxB = 0f - minB; // scalable range for colouring both Conservation and Quality + + float min = 0f; + float max = 11f; + float qmin = 0f; + float qmax = 0f; + + char c; + + conservation.annotations = new Annotation[alWidth]; + + if (quality != null) + { + quality.graphMax = cons.qualityRange[1].floatValue(); + quality.annotations = new Annotation[alWidth]; + qmin = cons.qualityRange[0].floatValue(); + qmax = cons.qualityRange[1].floatValue(); + } + + for (int i = 0; i < alWidth; i++) + { + float value = 0; + + c = sequence[i]; + + if (Character.isDigit(c)) + { + value = (int) (c - '0'); + } + else if (c == '*') + { + value = 11; + } + else if (c == '+') + { + value = 10; + } + + float vprop = value - min; + vprop /= max; + conservation.annotations[i] = + new Annotation(String.valueOf(c), + String.valueOf(value), ' ', value, + new Color(minR + (maxR * vprop), + minG + (maxG * vprop), + minB + (maxB * vprop))); + + // Quality calc + if (quality != null) + { + value = ( (Double) cons.quality.elementAt(i)).floatValue(); + vprop = value - qmin; + vprop /= qmax; + quality.annotations[i] = new Annotation(" ", String.valueOf(value), + ' ', + value, + new Color(minR + (maxR * vprop), + minG + (maxG * vprop), + minB + (maxB * vprop))); + } + } + } + catch (OutOfMemoryError error) { - alignment.addAnnotation(conservation); + System.out.println("Out of memory calculating conservation!!"); + conservation = null; + quality = null; + System.gc(); } - quality = new AlignmentAnnotation("Quality", - "Alignment Quality based on Blosum62 scores", - qannotations, - cons.qualityRange[0].floatValue(), - cons.qualityRange[1].floatValue(), - AlignmentAnnotation.BAR_GRAPH); - if (showQuality) + + UPDATING_CONSERVATION = false; + updatingConservation = false; + + if (ap != null) { - alignment.addAnnotation(quality); + ap.paintAlignment(true); } + } - else + } + + ConservationThread conservationThread; + + ConsensusThread consensusThread; + + boolean consUpdateNeeded = false; + + static boolean UPDATING_CONSENSUS = false; + + static boolean UPDATING_CONSERVATION = false; + + boolean updatingConsensus = false; + + boolean updatingConservation = false; + + /** + * DOCUMENT ME! + */ + public void updateConservation(final AlignmentPanel ap) + { + if (alignment.isNucleotide() || conservation == null) { - conservation.annotations = annotations; - quality.annotations = qannotations; - quality.graphMax = cons.qualityRange[1].floatValue(); + return; } + conservationThread = new ConservationThread(ap); + conservationThread.start(); } - public void updateConsensus() + /** + * DOCUMENT ME! + */ + public void updateConsensus(final AlignmentPanel ap) { - Annotation[] annotations = new Annotation[alignment.getWidth()]; + consensusThread = new ConsensusThread(ap); + consensusThread.start(); + } - // this routine prevents vconsensus becoming a new object each time - // consenus is calculated. Important for speed of Blosum62 - // and PID colouring of alignment - if (vconsensus == null) + class ConsensusThread + extends Thread + { + AlignmentPanel ap; + public ConsensusThread(AlignmentPanel ap) { - vconsensus = alignment.getAAFrequency(); + this.ap = ap; } - else + + public void run() { - Vector temp = alignment.getAAFrequency(); - vconsensus.removeAllElements(); - Enumeration e = temp.elements(); - while (e.hasMoreElements()) + updatingConsensus = true; + while (UPDATING_CONSENSUS) { - vconsensus.addElement(e.nextElement()); + try + { + if (ap != null) + { + ap.paintAlignment(true); + } + + Thread.sleep(200); + } + catch (Exception ex) + { + ex.printStackTrace(); + } } - } - Hashtable hash = null; - for (int i = 0; i < alignment.getWidth(); i++) - { - hash = (Hashtable) vconsensus.elementAt(i); - float value = 0; - if(ignoreGapsInConsensusCalculation) - value = ((Float)hash.get("pid_nogaps")).floatValue(); - else - value = ((Float)hash.get("pid_gaps")).floatValue(); - String maxRes = hash.get("maxResidue").toString(); - String mouseOver = hash.get("maxResidue") + " "; - if (maxRes.length() > 1) + UPDATING_CONSENSUS = true; + + try { - mouseOver = "[" + maxRes + "] "; - maxRes = "+"; - } + int aWidth = alignment.getWidth(); + if (aWidth < 0) + { + return; + } + + consensus.annotations = null; + consensus.annotations = new Annotation[aWidth]; + hconsensus = new Hashtable[aWidth]; + AAFrequency.calculate(alignment.getSequencesArray(), + 0, + alignment.getWidth(), + hconsensus); - mouseOver += (int) value + "%"; - annotations[i] = new Annotation(maxRes, mouseOver, ' ', value); + for (int i = 0; i < aWidth; i++) + { + float value = 0; + if (ignoreGapsInConsensusCalculation) + { + value = ( (Float) hconsensus[i].get(AAFrequency.PID_NOGAPS)). + floatValue(); + } + else + { + value = ( (Float) hconsensus[i].get(AAFrequency.PID_GAPS)). + floatValue(); + } + + String maxRes = hconsensus[i].get(AAFrequency.MAXRESIDUE).toString(); + String mouseOver = hconsensus[i].get(AAFrequency.MAXRESIDUE) + " "; + + if (maxRes.length() > 1) + { + mouseOver = "[" + maxRes + "] "; + maxRes = "+"; + } + mouseOver += ( (int) value + "%"); + consensus.annotations[i] = new Annotation(maxRes, mouseOver, ' ', + value); + } + + if (globalColourScheme != null) + { + globalColourScheme.setConsensus(hconsensus); + } + + } + catch (OutOfMemoryError error) + { + alignment.deleteAnnotation(consensus); + + consensus = null; + hconsensus = null; + System.out.println("Out of memory calculating consensus!!"); + System.gc(); + } + UPDATING_CONSENSUS = false; + updatingConsensus = false; + + if (ap != null) + { + ap.paintAlignment(true); + } } + } + /** + * get the consensus sequence as displayed under the PID consensus annotation row. + * @return consensus sequence as a new sequence object + */ + /** + * get the consensus sequence as displayed under the PID consensus annotation row. + * @return consensus sequence as a new sequence object + */ + public SequenceI getConsensusSeq() + { if (consensus == null) { - consensus = new AlignmentAnnotation("Consensus", - "PID", annotations, 0f, 100f, AlignmentAnnotation.BAR_GRAPH); - if (showConsensus) - { - alignment.addAnnotation(consensus); - } + return null; } - else + StringBuffer seqs = new StringBuffer(); + for (int i = 0; i < consensus.annotations.length; i++) { - consensus.annotations = annotations; + if (consensus.annotations[i] != null) + { + if (consensus.annotations[i].description.charAt(0) == '[') + { + seqs.append(consensus.annotations[i].description.charAt(1)); + } + else + { + seqs.append(consensus.annotations[i].displayCharacter); + } + } } - - if(globalColourScheme!=null) - globalColourScheme.setConsensus(vconsensus); - + SequenceI sq = new Sequence("Consensus", seqs.toString()); + sq.setDescription("Percentage Identity Consensus " + + ( (ignoreGapsInConsensusCalculation) ? " without gaps" : + "")); + return sq; } public SequenceGroup getSelectionGroup() @@ -454,7 +674,7 @@ public class AlignViewport public void setFont(Font f) { font = f; - if(nullFrame == null) + if (nullFrame == null) { nullFrame = new java.awt.Frame(); nullFrame.addNotify(); @@ -462,7 +682,14 @@ public class AlignViewport java.awt.FontMetrics fm = nullFrame.getGraphics().getFontMetrics(font); setCharHeight(fm.getHeight()); - setCharWidth(fm.charWidth('M')); + charWidth = fm.charWidth('M'); + + if (upperCasebold) + { + Font f2 = new Font(f.getName(), Font.BOLD, f.getSize()); + fm = nullFrame.getGraphics().getFontMetrics(f2); + charWidth = fm.stringWidth("MMMMMMMMMMM") / 10; + } } public Font getFont() @@ -470,11 +697,6 @@ public class AlignViewport return font; } - public void setCharWidth(int w) - { - this.charWidth = w; - } - public int getCharWidth() { return charWidth; @@ -591,8 +813,10 @@ public class AlignViewport public void setHiddenColumns(ColumnSelection colsel) { this.colSel = colsel; - if(colSel.getHiddenColumns()!=null) + if (colSel.getHiddenColumns() != null) + { hasHiddenColumns = true; + } } public ColumnSelection getColumnSelection() @@ -678,11 +902,11 @@ public class AlignViewport public void setIgnoreGapsConsensus(boolean b) { ignoreGapsInConsensusCalculation = b; - updateConsensus(); - if (globalColourScheme!=null) + updateConsensus(null); + if (globalColourScheme != null) { globalColourScheme.setThreshold(globalColourScheme.getThreshold(), - ignoreGapsInConsensusCalculation); + ignoreGapsInConsensusCalculation); } } @@ -695,7 +919,7 @@ public class AlignViewport public void addPropertyChangeListener( java.beans.PropertyChangeListener listener) { - changeSupport.addPropertyChangeListener(listener); + changeSupport.addPropertyChangeListener(listener); } /** @@ -706,7 +930,7 @@ public class AlignViewport public void removePropertyChangeListener( java.beans.PropertyChangeListener listener) { - changeSupport.removePropertyChangeListener(listener); + changeSupport.removePropertyChangeListener(listener); } /** @@ -718,19 +942,20 @@ public class AlignViewport */ public void firePropertyChange(String prop, Object oldvalue, Object newvalue) { - changeSupport.firePropertyChange(prop, oldvalue, newvalue); + changeSupport.firePropertyChange(prop, oldvalue, newvalue); } - - public boolean getIgnoreGapsConsensus() { return ignoreGapsInConsensusCalculation; } + public void hideSelectedColumns() { if (colSel.size() < 1) + { return; + } colSel.hideSelectedColumns(); setSelectionGroup(null); @@ -740,61 +965,106 @@ public class AlignViewport public void invertColumnSelection() { - int column; for (int i = 0; i < alignment.getWidth(); i++) { - column = i; - - if (colSel.contains(column)) - colSel.removeElement(column); + if (colSel.contains(i)) + { + colSel.removeElement(i); + } else - colSel.addElement(column); - + { + if (!hasHiddenColumns || colSel.isVisible(i)) + { + colSel.addElement(i); + } + } } } - public void hideColumns(int start, int end) { - if(start==end) + if (start == end) + { colSel.hideColumns(start); + } else + { colSel.hideColumns(start, end); + } hasHiddenColumns = true; } - public void hideSequence(SequenceI seq) + public void hideRepSequences(SequenceI repSequence, SequenceGroup sg) { - if(seq!=null) + int sSize = sg.getSize(); + if (sSize < 2) { - alignment.getHiddenSequences().hideSequence(seq); - hasHiddenRows = true; - firePropertyChange("alignment", null, alignment.getSequences()); + return; + } + + if (hiddenRepSequences == null) + { + hiddenRepSequences = new Hashtable(); + } + + hiddenRepSequences.put(repSequence, sg); + + //Hide all sequences except the repSequence + SequenceI[] seqs = new SequenceI[sSize - 1]; + int index = 0; + for (int i = 0; i < sSize; i++) + { + if (sg.getSequenceAt(i) != repSequence) + { + if (index == sSize - 1) + { + return; + } + + seqs[index++] = sg.getSequenceAt(i); + } } + + hideSequence(seqs); + } public void hideAllSelectedSeqs() { - if (selectionGroup == null) + if (selectionGroup == null || selectionGroup.getSize()<1) + { return; + } SequenceI[] seqs = selectionGroup.getSequencesInOrder(alignment); - for (int i = 0; i < seqs.length; i++) + hideSequence(seqs); + + setSelectionGroup(null); + } + + public void hideSequence(SequenceI[] seq) + { + if (seq != null) { - alignment.getHiddenSequences().hideSequence(seqs[i]); + for (int i = 0; i < seq.length; i++) + { + alignment.getHiddenSequences().hideSequence(seq[i]); + } + + hasHiddenRows = true; + firePropertyChange("alignment", null, alignment.getSequences()); } - firePropertyChange("alignment", null, alignment.getSequences()); - hasHiddenRows = true; - setSelectionGroup(null); } public void showColumn(int col) { colSel.revealHiddenColumns(col); - if(colSel.getHiddenColumns()==null) + if (colSel.getHiddenColumns() == null) + { hasHiddenColumns = false; + } } public void showAllHiddenColumns() @@ -805,22 +1075,23 @@ public class AlignViewport public void showAllHiddenSeqs() { - if(alignment.getHiddenSequences().getSize()>0) + if (alignment.getHiddenSequences().getSize() > 0) { - if(selectionGroup==null) + if (selectionGroup == null) { selectionGroup = new SequenceGroup(); - selectionGroup.setEndRes(alignment.getWidth()-1); + selectionGroup.setEndRes(alignment.getWidth() - 1); } - Vector tmp = alignment.getHiddenSequences().showAll(); - for(int t=0; t