JAL-2844 entirely X coord based partitioning
[jalview.git] / src / jalview / ext / archaeopteryx / JalviewBinding.java
index dae01f1..1649b62 100644 (file)
@@ -1,18 +1,31 @@
 package jalview.ext.archaeopteryx;
 
+import jalview.analysis.Conservation;
+import jalview.api.AlignViewportI;
 import jalview.datamodel.ColumnSelection;
 import jalview.datamodel.HiddenColumns;
 import jalview.datamodel.SequenceGroup;
 import jalview.datamodel.SequenceI;
 import jalview.ext.treeviewer.ExternalTreeViewerBindingI;
+import jalview.gui.AlignViewport;
+import jalview.gui.Desktop;
+import jalview.gui.JvOptionPane;
 import jalview.gui.PaintRefresher;
+import jalview.schemes.ColourSchemeI;
+import jalview.schemes.ColourSchemeProperty;
+import jalview.schemes.UserColourScheme;
 import jalview.structure.SelectionSource;
 import jalview.structure.StructureSelectionManager;
+import jalview.util.MappingUtils;
+import jalview.util.MessageManager;
 import jalview.viewmodel.AlignmentViewport;
 
+import java.awt.Color;
+import java.awt.Graphics;
 import java.awt.event.ActionEvent;
 import java.awt.event.InputEvent;
 import java.awt.event.MouseEvent;
+import java.util.ArrayList;
 import java.util.HashSet;
 import java.util.List;
 import java.util.Map;
@@ -24,9 +37,11 @@ import javax.swing.event.ChangeEvent;
 import javax.swing.event.ChangeListener;
 
 import org.forester.archaeopteryx.MainFrame;
+import org.forester.archaeopteryx.TreePanelUtil;
 import org.forester.phylogeny.Phylogeny;
 import org.forester.phylogeny.PhylogenyMethods;
 import org.forester.phylogeny.PhylogenyNode;
+import org.forester.phylogeny.data.BranchColor;
 
 /**
  * Class for binding the Archaeopteryx tree viewer to the Jalview alignment that
@@ -43,7 +58,7 @@ public final class JalviewBinding
 
   private AlignmentViewport parentAvport;
 
-  private JTabbedPane treeTabs;
+  private final JTabbedPane treeTabs;
 
   private final StructureSelectionManager ssm;
 
@@ -51,6 +66,10 @@ public final class JalviewBinding
 
   private Map<PhylogenyNode, SequenceI> nodesBoundToSequences;
 
+  private float rootX;
+
+  private float furthestNodeX;
+
   /**
    * 
    * @param archaeopteryx
@@ -71,6 +90,16 @@ public final class JalviewBinding
           final Map<SequenceI, PhylogenyNode> alignMappedToNodes,
           final Map<PhylogenyNode, SequenceI> nodesMappedToAlign)
   {
+
+    if (archaeopteryx.getMainPanel().getTabbedPane().getTabCount() > 1)
+    {
+      JvOptionPane.showMessageDialog(Desktop.desktop,
+              MessageManager.getString("label.tabs_detected_archaeopteryx"),
+              MessageManager.getString("label.problem_reading_tree_file"),
+              JvOptionPane.WARNING_MESSAGE);
+      ;
+    }
+
     // deal with/prohibit null values here as that will cause problems
     parentAvport = jalviewAlignmentViewport;
     sequencesBoundToNodes = alignMappedToNodes;
@@ -79,8 +108,6 @@ public final class JalviewBinding
     treeView = archaeopteryx.getMainPanel().getCurrentTreePanel();
     treeTabs = archaeopteryx.getMainPanel().getTabbedPane();
     ssm = parentAvport.getStructureSelectionManager();
-
-    // archaeopteryx.getMainPanel().getControlPanel().setColorBranches(true);
     
     ssm.addSelectionListener(this);
     treeView.addMouseListener(this);
@@ -106,6 +133,7 @@ public final class JalviewBinding
            */
           public void run()
           {
+            treeView = archaeopteryx.getMainPanel().getCurrentTreePanel();
             parentAvport.sendSelection();
             // PaintRefresher.Refresh(treeView,
             // parentAvport.getSequenceSetId());
@@ -114,7 +142,7 @@ public final class JalviewBinding
         });
 
       }
-      
+
     });
 
   }
@@ -149,11 +177,13 @@ public final class JalviewBinding
         }
         else
         {
-          partitionTree(e);
-        
-        
+
+          partitionTree(e.getX());
       }
-      
+        PaintRefresher.Refresh(treeView, parentAvport.getSequenceSetId());
+
+
+
       }
     });
 
@@ -198,11 +228,22 @@ public final class JalviewBinding
         if (matchingNode != null)
         {
           treeView.getFoundNodes0().add(matchingNode.getId());
+
+
+          if (!matchingNode.getBranchData().isHasBranchColor())
+          {
+            Color foundNodesColour = treeView.getTreeColorSet()
+                    .getFoundColor0();
+          matchingNode.getBranchData()
+                    .setBranchColor(new BranchColor(foundNodesColour));
+
+          }
+
         }
 
       }
-      treeView.repaint();
 
+      treeView.repaint();
     }
 
 
@@ -211,35 +252,217 @@ public final class JalviewBinding
   /**
    * Partially refactored from TreeCanvas
    */
-  public void partitionTree(final MouseEvent e)
+  public void partitionTree(final int x)
   {
-    int x = e.getX();
-    int lineLength = treeView.getHeight();
-
     Phylogeny tree = treeView.getPhylogeny();
-    double treeHeight = tree.calculateHeight(true);
 
+    if (!tree.isEmpty())
+    {
+      // should be calculated on each partition as the tree can theoretically
+      // change in the meantime
+      PhylogenyNode furthestNode = PhylogenyMethods
+              .calculateNodeWithMaxDistanceToRoot(tree);
+      furthestNodeX = furthestNode.getXcoord();
+      rootX = tree.getRoot().getXcoord();
+
+      if (furthestNodeX != rootX && x < furthestNodeX) // don't bother if 0
+                                                       // distance tree or
+                                                       // clicked x lies outside
+                                                       // of tree
+      {
+        Graphics g = treeView.getGraphics();
+        int panelHeight = treeView.getHeight();
+        g.drawLine(x, 0, x, panelHeight);
 
+        float threshold = (x - rootX) / (furthestNodeX - rootX);
+        List<PhylogenyNode> foundNodes = getNodesAboveThreshold(threshold,
+                tree.getRoot());
 
-    if (treeHeight != 0)
+      }
+    }
+
+
+  }
+
+  public List<PhylogenyNode> getNodesAboveThreshold(double threshold,
+          PhylogenyNode node)
+  {
+
+    List<PhylogenyNode> nodesAboveThreshold = new ArrayList<>();
+
+    parentAvport.setSelectionGroup(null);
+    parentAvport.getAlignment().deleteAllGroups();
+    parentAvport.clearSequenceColours();
+    if (parentAvport.getCodingComplement() != null)
     {
-      int viewWidth = treeView.getWidth();
+      parentAvport.getCodingComplement().setSelectionGroup(null);
+      parentAvport.getCodingComplement().getAlignment().deleteAllGroups();
+      parentAvport.getCodingComplement().clearSequenceColours();
+    }
+
+
+    colourNodesAboveThreshold(nodesAboveThreshold, threshold,
+            node);
+    return nodesAboveThreshold;
+
+  }
+
+  /**
+   * Partially refactored from TreeCanvas colourGroups (can be made nicer).
+   * 
+   * @param nodeList
+   * @param threshold
+   * @param treeLength
+   * @param node
+   * @return
+   */
+  private List<PhylogenyNode> colourNodesAboveThreshold(
+          List<PhylogenyNode> nodeList, double threshold,
+          PhylogenyNode node)
+  {
+    // could also use PhylogenyMethods.getAllDescendants
+    for (PhylogenyNode childNode : node.getDescendants())
+    {
+      childNode.getBranchData()
+              .setBranchColor(new BranchColor(Color.black));
+      float nodeCutoff = (childNode.getXcoord() - rootX)
+              / (furthestNodeX - rootX);
+
+      if (nodeCutoff > threshold)
+      {
+        nodeList.add(childNode);
+
+        Color randomColor = new Color((int) (Math.random() * 255),
+                (int) (Math.random() * 255), (int) (Math.random() * 255));
+        TreePanelUtil.colorizeSubtree(childNode,
+                new BranchColor(randomColor));
+        List<PhylogenyNode> descendantNodes = childNode
+                .getAllExternalDescendants();
+        List<SequenceI> descendantSeqs = new ArrayList<>(); // .forEach instead?
+        for (PhylogenyNode descNode : descendantNodes)
+        {
+          descendantSeqs.add(nodesBoundToSequences.get(descNode));
+        }
 
-      // treeView.validate();
+        SequenceGroup sg = new SequenceGroup(descendantSeqs, null, null,
+                true, true, false, 0,
+                parentAvport.getAlignment().getWidth() - 1);
 
-      // System.out.println("selection");
-      // System.out.println(x);
-      // System.out.println("-------------");
-      // System.out.println("width");
-      // System.out.println(viewWidth);
+        ColourSchemeI cs = null;
+        if (parentAvport.getGlobalColourScheme() != null)
+        {
+          if (parentAvport.getGlobalColourScheme() instanceof UserColourScheme)
+          {
+            cs = new UserColourScheme(
+                    ((UserColourScheme) parentAvport.getGlobalColourScheme())
+                            .getColours());
+          }
+          else
+          {
+            cs = ColourSchemeProperty.getColourScheme(sg, ColourSchemeProperty
+                    .getColourName(parentAvport.getGlobalColourScheme()));
+          }
+        }
+        sg.setColourScheme(cs);
+        sg.getGroupColourScheme().setThreshold(
+                parentAvport.getResidueShading().getThreshold(),
+                parentAvport.isIgnoreGapsConsensus());
+        // sg.recalcConservation();
+        sg.setName("Tree Group:" + sg.hashCode());
+        sg.setIdColour(randomColor);
+
+        if (parentAvport.getGlobalColourScheme() != null
+                && parentAvport.getResidueShading().conservationApplied())
+        {
+          Conservation c = new Conservation("Group", sg.getSequences(null),
+                  sg.getStartRes(), sg.getEndRes());
+          c.calculate();
+          c.verdict(false, parentAvport.getConsPercGaps());
+          sg.cs.setConservation(c);
+        }
+
+        parentAvport.getAlignment().addGroup(new SequenceGroup(sg));
+        // TODO can we push all of the below into AlignViewportI?
+        final AlignViewportI codingComplement = parentAvport
+                .getCodingComplement();
+        if (codingComplement != null)
+        {
+          SequenceGroup mappedGroup = MappingUtils.mapSequenceGroup(sg,
+                  parentAvport,
+                  codingComplement);
+          if (mappedGroup.getSequences().size() > 0)
+          {
+            codingComplement.getAlignment().addGroup(mappedGroup);
+            for (SequenceI seq : mappedGroup.getSequences())
+            {
+              codingComplement.setSequenceColour(seq,
+                      randomColor.brighter());
+            }
+          }
+        }
+
+
+      }
 
+      else
+      {
+        colourNodesAboveThreshold(nodeList, threshold,
+                childNode);
+      }
     }
 
+    // GROSS
+    ((AlignViewport) parentAvport).getAlignPanel().updateAnnotation();
+
+    final AlignViewportI codingComplement = parentAvport
+            .getCodingComplement();
+    if (codingComplement != null)
+    {
+      ((AlignViewport) codingComplement).getAlignPanel().updateAnnotation();
+    }
 
+    return nodeList;
   }
-  
 
 
+
+
+  // public List<PhylogenyNode> groupNodes(float threshold, PhylogenyNode root,
+  // double treeHeight)
+  // {
+  // List<PhylogenyNode> groups = new ArrayList<>();
+  // _groupNodes(groups, root, threshold, treeHeight);
+  // System.out.println(groups);
+  // return groups;
+  // }
+  //
+  // protected void _groupNodes(List<PhylogenyNode> groups, PhylogenyNode nd,
+  // float threshold, double treeHeight)
+  // {
+  // if (nd == null)
+  // {
+  // return;
+  // }
+  //
+  // if ((nd.calculateDistanceToRoot() / treeHeight) > threshold)
+  // {
+  // groups.add(nd);
+  // }
+  // else
+  // {
+  // for (PhylogenyNode childNode : nd.getDescendants())
+  // {
+  // _groupNodes(groups, childNode, threshold, treeHeight);
+  // }
+  // }
+  // }
+  //
+
+
+  /**
+   * may or may not need an extra repaint on the alignment view (check what kira
+   * does)
+   */
   @Override
   public void showNodeSelectionOnAlign(final PhylogenyNode node)
   {
@@ -254,6 +477,7 @@ public final class JalviewBinding
         showMatchingSequence(node);
       }
 
+
     }
 
 
@@ -349,9 +573,9 @@ public final class JalviewBinding
   // for (AlignmentPanel ap : PaintRefresher
   // .getAssociatedPanels(parentAvport.getSequenceSetId()))
   // {
-  // commands.add(sortAlignmentIn(ap.av.getAlignPanel()));
+  // commands.add(sortAlignmentIn(ap.parentAvport.getAlignPanel()));
   // }
-  // av.getAlignPanel().alignFrame.addHistoryItem(new CommandI()
+  // parentAvport.getAlignPanel().alignFrame.addHistoryItem(new CommandI()
   // {
   //
   // @Override
@@ -386,7 +610,7 @@ public final class JalviewBinding
   // }
   // });
   // for (AlignmentPanel ap : PaintRefresher
-  // .getAssociatedPanels(av.getSequenceSetId()))
+  // .getAssociatedPanels(parentAvport.getSequenceSetId()))
   // {
   // // ensure all the alignFrames refresh their GI after adding an undo item
   // ap.alignFrame.updateEditMenuBar();