X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Fappletgui%2FAlignViewport.java;h=5eb5e468db783321f31f20ac0e477474511fd3dc;hb=32a72d077d2371f4566688fbd103103d8505b27f;hp=99d0794f3c9adb3cfb933b54b31f7ac131fbd7ed;hpb=10cc0aa63a3c0df54838258157485f14e30cd874;p=jalview.git

diff --git a/src/jalview/appletgui/AlignViewport.java b/src/jalview/appletgui/AlignViewport.java
index 99d0794..5eb5e46 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) 2005 AM Waterhouse, J Procter, G Barton, M Clamp, S Searle
+ * Copyright (C) 2006 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
@@ -50,6 +50,7 @@ public class AlignViewport
   boolean showConservation = true;
   boolean showQuality = true;
   boolean showConsensus = true;
+  boolean upperCasebold = false;
 
   boolean colourAppliesToAllGroups = true;
   ColourSchemeI globalColourScheme = null;
@@ -79,14 +80,14 @@ public class AlignViewport
 
   // The following vector holds the features which are
  // currently visible, in the correct order or rendering
-  Hashtable featuresDisplayed;
+  public Hashtable featuresDisplayed;
 
   boolean hasHiddenColumns = false;
   boolean hasHiddenRows = false;
   boolean showHiddenMarkers = true;
 
 
-  public Vector vconsensus;
+  public Hashtable [] hconsensus;
   AlignmentAnnotation consensus;
   AlignmentAnnotation conservation;
   AlignmentAnnotation quality;
@@ -101,6 +102,15 @@ public class AlignViewport
 
   jalview.bin.JalviewLite applet;
 
+  Hashtable sequenceColours;
+
+  boolean MAC = false;
+
+  Stack historyList = new Stack();
+  Stack redoList = new Stack();
+
+  String sequenceSetID;
+
   public AlignViewport(AlignmentI al, JalviewLite applet)
   {
     this.applet = applet;
@@ -111,6 +121,9 @@ 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,12 +155,15 @@ 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)
     {
@@ -165,7 +181,7 @@ public class AlignViewport
         globalColourScheme = ColourSchemeProperty.getColour(alignment, colour);
         if (globalColourScheme != null)
         {
-          globalColourScheme.setConsensus(vconsensus);
+          globalColourScheme.setConsensus(hconsensus);
         }
       }
 
@@ -175,7 +191,48 @@ public class AlignViewport
             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;
+
+
+          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;
+
+            alignment.addAnnotation(quality);
+          }
+        }
+
+        consensus = new AlignmentAnnotation("Consensus", "PID",
+                                             new Annotation[1], 0f, 100f,
+                                             AlignmentAnnotation.BAR_GRAPH);
+        consensus.hasText = true;
 
+         if (showConsensus)
+         {
+           alignment.addAnnotation(consensus);
+         }
+      }
     }
   }
 
@@ -190,169 +247,296 @@ public class AlignViewport
   }
 
 
-  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)
+        {
+          try
+          {
+            if (ap != null)
+            {
+              ap.repaint();
+            }
+            Thread.sleep(200);
+          }
+          catch (Exception ex)
+          {
+            ex.printStackTrace();
+          }
+        }
+
+        UPDATING_CONSERVATION = true;
+
+
+        int alWidth = alignment.getWidth();
+        if(alWidth<0)
+          return;
+
+        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)
         {
-          value = 11;
+          cons.findQuality();
         }
-        if (sequence.charAt(i) == '+')
+
+        String 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)
         {
-          value = 10;
+          quality.graphMax = cons.qualityRange[1].floatValue();
+          quality.annotations = new Annotation[alWidth];
+          qmin = cons.qualityRange[0].floatValue();
+          qmax = cons.qualityRange[1].floatValue();
         }
-      }
-      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)
+        for (int i = 0; i < alWidth; i++)
+        {
+          float value = 0;
+
+          c = sequence.charAt(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.repaint();
       }
+
     }
-    else
-    {
-      conservation.annotations = annotations;
-      quality.annotations = qannotations;
-      quality.graphMax = cons.qualityRange[1].floatValue();
-    }
+  }
+
+
+  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)
+      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.repaint();
+          }
+
+          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];
 
-      mouseOver += (int) value + "%";
-      annotations[i] = new Annotation(maxRes, mouseOver, ' ', value);
 
-    }
+        hconsensus = new Hashtable[aWidth];
+        AAFrequency.calculate(alignment.getSequencesArray(),
+                              0,
+                              alignment.getWidth(),
+                              hconsensus);
+
+        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);
 
-    if (consensus == null)
-    {
-      consensus = new AlignmentAnnotation("Consensus",
-                                          "PID", annotations, 0f, 100f, AlignmentAnnotation.BAR_GRAPH);
-      if (showConsensus)
-      {
-        alignment.addAnnotation(consensus);
       }
-    }
-    else
-    {
-      consensus.annotations = annotations;
-    }
+      catch (OutOfMemoryError error)
+      {
+        alignment.deleteAnnotation(consensus);
 
-    if(globalColourScheme!=null)
-          globalColourScheme.setConsensus(vconsensus);
+        consensus = null;
+        hconsensus = null;
+        System.out.println("Out of memory calculating consensus!!");
+        System.gc();
+      }
+      UPDATING_CONSENSUS = false;
+      updatingConsensus = false;
 
+      if (ap != null)
+      {
+        ap.repaint();
+      }
+    }
   }
 
+  /**
+   * 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)
+      return null;
+    StringBuffer seqs=new StringBuffer();
+    for (int i=0; i<consensus.annotations.length; i++) {
+      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);
+      }
+    }
+    SequenceI sq = new Sequence("Consensus", seqs.toString());
+    sq.setDescription("Percentage Identity Consensus "+((ignoreGapsInConsensusCalculation) ? " without gaps" : ""));
+    return sq;
+  }
   public SequenceGroup getSelectionGroup()
   {
     return selectionGroup;
@@ -462,7 +646,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 +661,6 @@ public class AlignViewport
     return font;
   }
 
-  public void setCharWidth(int w)
-  {
-    this.charWidth = w;
-  }
-
   public int getCharWidth()
   {
     return charWidth;
@@ -678,7 +864,7 @@ public class AlignViewport
   public void setIgnoreGapsConsensus(boolean b)
   {
     ignoreGapsInConsensusCalculation = b;
-    updateConsensus();
+    updateConsensus(null);
     if (globalColourScheme!=null)
     {
       globalColourScheme.setThreshold(globalColourScheme.getThreshold(),
@@ -760,25 +946,32 @@ public class AlignViewport
       colSel.hideColumns(start);
     else
       colSel.hideColumns(start, end);
-    setSelectionGroup(null);
+
     hasHiddenColumns = true;
   }
 
-  public void hideSequence(SequenceI seq)
+  public void hideAllSelectedSeqs()
   {
-    if(seq!=null)
-    {
-      alignment.getHiddenSequences().hideSequence(seq);
-      hasHiddenRows = true;
-    }
+    if (selectionGroup == null)
+      return;
+
+    SequenceI[] seqs = selectionGroup.getSequencesInOrder(alignment);
+
+    hideSequence(seqs);
+
+    setSelectionGroup(null);
   }
 
-  public void showSequence(int index)
+  public void hideSequence(SequenceI [] seq)
   {
-    alignment.getHiddenSequences().showSequence(index);
+    if(seq!=null)
+    {
+      for (int i = 0; i < seq.length; i++)
+        alignment.getHiddenSequences().hideSequence(seq[i]);
 
-    if(alignment.getHiddenSequences().getSize()<1)
-      hasHiddenRows = false;
+      hasHiddenRows = true;
+      firePropertyChange("alignment", null, alignment.getSequences());
+    }
   }
 
   public void showColumn(int col)
@@ -798,7 +991,19 @@ public class AlignViewport
   {
     if(alignment.getHiddenSequences().getSize()>0)
     {
-      alignment.getHiddenSequences().showAll();
+      if(selectionGroup==null)
+      {
+        selectionGroup = new SequenceGroup();
+        selectionGroup.setEndRes(alignment.getWidth()-1);
+      }
+      Vector tmp = alignment.getHiddenSequences().showAll();
+      for(int t=0; t<tmp.size(); t++)
+      {
+        selectionGroup.addSequence(
+            (SequenceI)tmp.elementAt(t), false
+            );
+      }
+      firePropertyChange("alignment", null, alignment.getSequences());
       hasHiddenRows = false;
     }
   }
@@ -896,9 +1101,9 @@ public class AlignViewport
       }
       // Final match if necessary.
       if (last<end)
-        selection.addOperation(CigarArray.M, end-last);
+        selection.addOperation(CigarArray.M, end-last+1);
     } else {
-      selection.addOperation(CigarArray.M, end-start);
+      selection.addOperation(CigarArray.M, end-start+1);
     }
     return selection;
   }
@@ -914,8 +1119,10 @@ public class AlignViewport
     // object. Refactoring of Cigar and alignment view representation should
     // be done to remove redundancy.
     CigarArray aligview = getViewAsCigars(selectedOnly);
-    if (aligview!=null)
-      return new AlignmentView(aligview);
+    if (aligview!=null) {
+      return new AlignmentView(aligview,
+          (selectedOnly && selectionGroup!=null) ? selectionGroup.getStartRes() : 0);
+    }
     return null;
   }
   /**
@@ -1010,5 +1217,103 @@ public class AlignViewport
     showHiddenMarkers = show;
   }
 
+  public Color getSequenceColour(SequenceI seq)
+  {
+    if (sequenceColours == null || !sequenceColours.containsKey(seq))
+      return Color.white;
+    else
+      return (Color) sequenceColours.get(seq);
+  }
+
+  public void setSequenceColour(SequenceI seq, Color col)
+  {
+    if (sequenceColours == null)
+      sequenceColours = new Hashtable();
+
+    if (col == null)
+      sequenceColours.remove(seq);
+    else
+      sequenceColours.put(seq, col);
+  }
+
+  public String getSequenceSetId()
+  {
+    if (sequenceSetID == null)
+      sequenceSetID = alignment.hashCode() + "";
+
+    return sequenceSetID;
+  }
+
+  public void alignmentChanged(AlignmentPanel ap)
+  {
+    alignment.padGaps();
+
+    if (hconsensus != null && autocalculateConsensus)
+    {
+      updateConsensus(ap);
+      updateConservation(ap);
+    }
+
+    //Reset endRes of groups if beyond alignment width
+    int alWidth = alignment.getWidth();
+    Vector groups = alignment.getGroups();
+    if(groups!=null)
+    {
+      for(int i=0; i<groups.size(); i++)
+      {
+        SequenceGroup sg = (SequenceGroup)groups.elementAt(i);
+        if(sg.getEndRes()>alWidth)
+          sg.setEndRes(alWidth-1);
+      }
+    }
+
+    if(selectionGroup!=null && selectionGroup.getEndRes()>alWidth)
+      selectionGroup.setEndRes(alWidth-1);
+
+    resetAllColourSchemes();
+
+    alignment.adjustSequenceAnnotations();
+  }
+
+  void resetAllColourSchemes()
+  {
+    ColourSchemeI cs = globalColourScheme;
+    if(cs!=null)
+    {
+      if (cs instanceof ClustalxColourScheme)
+      {
+        ( (ClustalxColourScheme) cs).
+            resetClustalX(alignment.getSequences(),
+                          alignment.getWidth());
+      }
+
+      cs.setConsensus(hconsensus);
+      if (cs.conservationApplied())
+      {
+        Alignment al = (Alignment) alignment;
+        Conservation c = new Conservation("All",
+                                          ResidueProperties.propHash, 3,
+                                          al.getSequences(), 0,
+                                          al.getWidth() - 1);
+        c.calculate();
+        c.verdict(false, ConsPercGaps);
+
+        cs.setConservation(c);
+      }
+    }
+
+    int s, sSize = alignment.getGroups().size();
+    for(s=0; s<sSize; s++)
+    {
+      SequenceGroup sg = (SequenceGroup)alignment.getGroups().elementAt(s);
+      if(sg.cs!=null && sg.cs instanceof ClustalxColourScheme)
+      {
+        ((ClustalxColourScheme)sg.cs).resetClustalX(
+            sg.getSequences(true), sg.getWidth());
+      }
+      sg.recalcConservation();
+    }
+  }
+
 
 }