src/jalview/controller/AlignViewController.java

   1 /*
   2  * Jalview - A Sequence Alignment Editor and Viewer ($$Version-Rel$$)
   3  * Copyright (C) $$Year-Rel$$ The Jalview Authors
   4  *
   5  * This file is part of Jalview.
   6  *
   7  * Jalview is free software: you can redistribute it and/or
   8  * modify it under the terms of the GNU General Public License
   9  * as published by the Free Software Foundation, either version 3
  10  * of the License, or (at your option) any later version.
  11  *
  12  * Jalview is distributed in the hope that it will be useful, but
  13  * WITHOUT ANY WARRANTY; without even the implied warranty
  14  * of MERCHANTABILITY or FITNESS FOR A PARTICULAR
  15  * PURPOSE.  See the GNU General Public License for more details.
  16  *
  17  * You should have received a copy of the GNU General Public License
  18  * along with Jalview.  If not, see <http://www.gnu.org/licenses/>.
  19  * The Jalview Authors are detailed in the 'AUTHORS' file.
  20  */
  21 package jalview.controller;
  22
  23 import jalview.analysis.AlignmentSorter;
  24 import jalview.api.AlignViewControllerGuiI;
  25 import jalview.api.AlignViewControllerI;
  26 import jalview.api.AlignViewportI;
  27 import jalview.api.AlignmentViewPanel;
  28 import jalview.api.FeatureRenderer;
  29 import jalview.commands.OrderCommand;
  30 import jalview.datamodel.AlignmentI;
  31 import jalview.datamodel.ColumnSelection;
  32 import jalview.datamodel.SequenceCollectionI;
  33 import jalview.datamodel.SequenceFeature;
  34 import jalview.datamodel.SequenceGroup;
  35 import jalview.datamodel.SequenceI;
  36 import jalview.io.DataSourceType;
  37 import jalview.io.FeaturesFile;
  38 import jalview.util.MessageManager;
  39
  40 import java.awt.Color;
  41 import java.util.BitSet;
  42 import java.util.List;
  43
  44 public class AlignViewController implements AlignViewControllerI
  45 {
  46   AlignViewportI viewport = null;
  47
  48   AlignmentViewPanel alignPanel = null;
  49
  50   /**
  51    * the GUI container that is handling interactions with the user
  52    */
  53   private AlignViewControllerGuiI avcg;
  54
  55   @Override
  56   protected void finalize() throws Throwable
  57   {
  58     viewport = null;
  59     alignPanel = null;
  60     avcg = null;
  61   };
  62
  63   public AlignViewController(AlignViewControllerGuiI alignFrame,
  64           AlignViewportI viewport, AlignmentViewPanel alignPanel)
  65   {
  66     this.avcg = alignFrame;
  67     this.viewport = viewport;
  68     this.alignPanel = alignPanel;
  69   }
  70
  71   @Override
  72   public void setViewportAndAlignmentPanel(AlignViewportI viewport,
  73           AlignmentViewPanel alignPanel)
  74   {
  75     this.alignPanel = alignPanel;
  76     this.viewport = viewport;
  77
  78   }
  79
  80   @Override
  81   public boolean makeGroupsFromSelection()
  82   {
  83     SequenceGroup sg = viewport.getSelectionGroup();
  84     ColumnSelection cs = viewport.getColumnSelection();
  85     SequenceGroup[] gps = null;
  86     if (sg != null && (cs == null || cs.isEmpty()))
  87     {
  88       gps = jalview.analysis.Grouping.makeGroupsFrom(
  89               viewport.getSequenceSelection(),
  90               viewport.getAlignmentView(true)
  91                       .getSequenceStrings(viewport.getGapCharacter()),
  92               viewport.getAlignment().getGroups());
  93     }
  94     else
  95     {
  96       if (cs != null)
  97       {
  98         gps = jalview.analysis.Grouping.makeGroupsFromCols(
  99                 (sg == null) ? viewport.getAlignment().getSequencesArray()
 100                         : sg.getSequences().toArray(new SequenceI[0]),
 101                 cs, viewport.getAlignment().getGroups());
 102       }
 103     }
 104     if (gps != null)
 105     {
 106       viewport.getAlignment().deleteAllGroups();
 107       viewport.clearSequenceColours();
 108       viewport.setSelectionGroup(null);
 109       // set view properties for each group
 110       for (int g = 0; g < gps.length; g++)
 111       {
 112         // gps[g].setShowunconserved(viewport.getShowUnconserved());
 113         gps[g].setshowSequenceLogo(viewport.isShowSequenceLogo());
 114         viewport.getAlignment().addGroup(gps[g]);
 115         Color col = new Color((int) (Math.random() * 255),
 116                 (int) (Math.random() * 255), (int) (Math.random() * 255));
 117         col = col.brighter();
 118         for (SequenceI sq : gps[g].getSequences(null))
 119         {
 120           viewport.setSequenceColour(sq, col);
 121         }
 122       }
 123       return true;
 124     }
 125     return false;
 126   }
 127
 128   @Override
 129   public boolean createGroup()
 130   {
 131
 132     SequenceGroup sg = viewport.getSelectionGroup();
 133     if (sg != null)
 134     {
 135       viewport.getAlignment().addGroup(sg);
 136       return true;
 137     }
 138     return false;
 139   }
 140
 141   @Override
 142   public boolean unGroup()
 143   {
 144     SequenceGroup sg = viewport.getSelectionGroup();
 145     if (sg != null)
 146     {
 147       viewport.getAlignment().deleteGroup(sg);
 148       return true;
 149     }
 150     return false;
 151   }
 152
 153   @Override
 154   public boolean deleteGroups()
 155   {
 156     if (viewport.getAlignment().getGroups() != null
 157             && viewport.getAlignment().getGroups().size() > 0)
 158     {
 159       viewport.getAlignment().deleteAllGroups();
 160       viewport.clearSequenceColours();
 161       viewport.setSelectionGroup(null);
 162       return true;
 163     }
 164     return false;
 165   }
 166
 167   @Override
 168   public boolean markColumnsContainingFeatures(boolean invert,
 169           boolean extendCurrent, boolean toggle, String featureType)
 170   {
 171     // JBPNote this routine could also mark rows, not just columns.
 172     // need a decent query structure to allow all types of feature searches
 173     BitSet bs = new BitSet();
 174     SequenceCollectionI sqcol = (viewport.getSelectionGroup() == null
 175             || extendCurrent) ? viewport.getAlignment()
 176                     : viewport.getSelectionGroup();
 177
 178     int nseq = findColumnsWithFeature(featureType, sqcol, bs);
 179
 180     ColumnSelection cs = viewport.getColumnSelection();
 181     if (cs == null)
 182     {
 183       cs = new ColumnSelection();
 184     }
 185
 186     if (bs.cardinality() > 0 || invert)
 187     {
 188       boolean changed = cs.markColumns(bs, sqcol.getStartRes(),
 189               sqcol.getEndRes(), invert, extendCurrent, toggle);
 190       if (changed)
 191       {
 192         viewport.setColumnSelection(cs);
 193         alignPanel.paintAlignment(false, false);
 194         int columnCount = invert
 195                 ? (sqcol.getEndRes() - sqcol.getStartRes() + 1)
 196                         - bs.cardinality()
 197                 : bs.cardinality();
 198         avcg.setStatus(MessageManager.formatMessage(
 199                 "label.view_controller_toggled_marked", new String[]
 200                 { toggle ? MessageManager.getString("label.toggled")
 201                         : MessageManager.getString("label.marked"),
 202                     String.valueOf(columnCount),
 203                     invert ? MessageManager
 204                             .getString("label.not_containing")
 205                             : MessageManager.getString("label.containing"),
 206                     featureType, Integer.valueOf(nseq).toString() }));
 207         return true;
 208       }
 209     }
 210     else
 211     {
 212       avcg.setStatus(MessageManager
 213               .formatMessage("label.no_feature_of_type_found", new String[]
 214               { featureType }));
 215       if (!extendCurrent)
 216       {
 217         cs.clear();
 218         alignPanel.paintAlignment(false, false);
 219       }
 220     }
 221     return false;
 222   }
 223
 224   /**
 225    * Sets a bit in the BitSet for each column (base 0) in the sequence
 226    * collection which includes the specified feature type. Returns the number of
 227    * sequences which have the feature in the selected range.
 228    *
 229    * @param featureType
 230    * @param sqcol
 231    * @param bs
 232    * @return
 233    */
 234   static int findColumnsWithFeature(String featureType,
 235           SequenceCollectionI sqcol, BitSet bs)
 236   {
 237     final int startColumn = sqcol.getStartRes() + 1; // converted to base 1
 238     final int endColumn = sqcol.getEndRes() + 1;
 239     List<SequenceI> seqs = sqcol.getSequences();
 240     int nseq = 0;
 241     for (SequenceI sq : seqs)
 242     {
 243       if (sq != null)
 244       {
 245         // int ist = sq.findPosition(sqcol.getStartRes());
 246         List<SequenceFeature> sfs = sq.findFeatures(startColumn,
 247                 endColumn, featureType);
 248
 249         if (!sfs.isEmpty())
 250         {
 251           nseq++;
 252         }
 253
 254         for (SequenceFeature sf : sfs)
 255         {
 256           int sfStartCol = sq.findIndex(sf.getBegin());
 257           int sfEndCol = sq.findIndex(sf.getEnd());
 258
 259           if (sf.isContactFeature())
 260           {
 261             /*
 262              * 'contact' feature - check for 'start' or 'end'
 263              * position within the selected region
 264              */
 265             if (sfStartCol >= startColumn && sfStartCol <= endColumn)
 266             {
 267               bs.set(sfStartCol - 1);
 268             }
 269             if (sfEndCol >= startColumn && sfEndCol <= endColumn)
 270             {
 271               bs.set(sfEndCol - 1);
 272             }
 273             continue;
 274           }
 275
 276           /*
 277            * contiguous feature - select feature positions (if any)
 278            * within the selected region
 279            */
 280           if (sfStartCol < startColumn)
 281           {
 282             sfStartCol = startColumn;
 283           }
 284           // not sure what the point of this is
 285           // if (sfStartCol < ist)
 286           // {
 287           // sfStartCol = ist;
 288           // }
 289           if (sfEndCol > endColumn)
 290           {
 291             sfEndCol = endColumn;
 292           }
 293           for (; sfStartCol <= sfEndCol; sfStartCol++)
 294           {
 295             bs.set(sfStartCol - 1); // convert to base 0
 296           }
 297         }
 298       }
 299     }
 300     return nseq;
 301   }
 302
 303   @Override
 304   public void sortAlignmentByFeatureDensity(List<String> typ)
 305   {
 306     sortBy(typ, "Sort by Density", AlignmentSorter.FEATURE_DENSITY);
 307   }
 308
 309   protected void sortBy(List<String> typ, String methodText,
 310           final String method)
 311   {
 312     FeatureRenderer fr = alignPanel.getFeatureRenderer();
 313     if (typ == null && fr != null)
 314     {
 315       typ = fr.getDisplayedFeatureTypes();
 316     }
 317     List<String> gps = null;
 318     if (fr != null)
 319     {
 320       gps = fr.getDisplayedFeatureGroups();
 321     }
 322     AlignmentI al = viewport.getAlignment();
 323
 324     int start, stop;
 325     SequenceGroup sg = viewport.getSelectionGroup();
 326     if (sg != null)
 327     {
 328       start = sg.getStartRes();
 329       stop = sg.getEndRes();
 330     }
 331     else
 332     {
 333       start = 0;
 334       stop = al.getWidth();
 335     }
 336     SequenceI[] oldOrder = al.getSequencesArray();
 337     AlignmentSorter.sortByFeature(typ, gps, start, stop, al, method);
 338     avcg.addHistoryItem(new OrderCommand(methodText, oldOrder,
 339             viewport.getAlignment()));
 340     alignPanel.paintAlignment(true, false);
 341
 342   }
 343
 344   @Override
 345   public void sortAlignmentByFeatureScore(List<String> typ)
 346   {
 347     sortBy(typ, "Sort by Feature Score", AlignmentSorter.FEATURE_SCORE);
 348   }
 349
 350   @Override
 351   public boolean parseFeaturesFile(String file, DataSourceType protocol,
 352           boolean relaxedIdMatching)
 353   {
 354     boolean featuresFile = false;
 355     try
 356     {
 357       featuresFile = new FeaturesFile(false, file, protocol).parse(
 358               viewport.getAlignment().getDataset(),
 359               alignPanel.getFeatureRenderer().getFeatureColours(), false,
 360               relaxedIdMatching);
 361     } catch (Exception ex)
 362     {
 363       ex.printStackTrace();
 364     }
 365
 366     if (featuresFile)
 367     {
 368       avcg.refreshFeatureUI(true);
 369       if (alignPanel.getFeatureRenderer() != null)
 370       {
 371         // update the min/max ranges where necessary
 372         alignPanel.getFeatureRenderer().findAllFeatures(true);
 373       }
 374       if (avcg.getFeatureSettingsUI() != null)
 375       {
 376         avcg.getFeatureSettingsUI().discoverAllFeatureData();
 377       }
 378       alignPanel.paintAlignment(true, true);
 379     }
 380
 381     return featuresFile;
 382
 383   }
 384
 385   @Override
 386   public boolean markHighlightedColumns(boolean invert,
 387           boolean extendCurrent, boolean toggle)
 388   {
 389     if (!viewport.hasSearchResults())
 390     {
 391       // do nothing if no selection exists
 392       return false;
 393     }
 394     // JBPNote this routine could also mark rows, not just columns.
 395     BitSet bs = new BitSet();
 396     SequenceCollectionI sqcol = (viewport.getSelectionGroup() == null
 397             || extendCurrent) ? viewport.getAlignment()
 398                     : viewport.getSelectionGroup();
 399
 400     // this could be a lambda... - the remains of the method is boilerplate,
 401     // except for the different messages for reporting selection.
 402     int nseq = viewport.getSearchResults().markColumns(sqcol, bs);
 403
 404     ColumnSelection cs = viewport.getColumnSelection();
 405     if (cs == null)
 406     {
 407       cs = new ColumnSelection();
 408     }
 409
 410     if (bs.cardinality() > 0 || invert)
 411     {
 412       boolean changed = cs.markColumns(bs, sqcol.getStartRes(),
 413               sqcol.getEndRes(), invert, extendCurrent, toggle);
 414       if (changed)
 415       {
 416         viewport.setColumnSelection(cs);
 417         alignPanel.paintAlignment(false, false);
 418         int columnCount = invert
 419                 ? (sqcol.getEndRes() - sqcol.getStartRes() + 1)
 420                         - bs.cardinality()
 421                 : bs.cardinality();
 422         avcg.setStatus(MessageManager.formatMessage(
 423                 "label.view_controller_toggled_marked", new String[]
 424                 { toggle ? MessageManager.getString("label.toggled")
 425                         : MessageManager.getString("label.marked"),
 426                     String.valueOf(columnCount),
 427                     invert ? MessageManager
 428                             .getString("label.not_containing")
 429                             : MessageManager.getString("label.containing"),
 430                     "Highlight", Integer.valueOf(nseq).toString() }));
 431         return true;
 432       }
 433     }
 434     else
 435     {
 436       avcg.setStatus(MessageManager
 437               .formatMessage("No highlighted regions marked"));
 438       if (!extendCurrent)
 439       {
 440         cs.clear();
 441         alignPanel.paintAlignment(false, false);
 442       }
 443     }
 444     return false;
 445   }
 446
 447 }