2 * Jalview - A Sequence Alignment Editor and Viewer ($$Version-Rel$$)
3 * Copyright (C) $$Year-Rel$$ The Jalview Authors
5 * This file is part of Jalview.
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.
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.
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.
21 package jalview.datamodel;
23 import jalview.util.Comparison;
24 import jalview.util.ShiftList;
25 import jalview.viewmodel.annotationfilter.AnnotationFilterParameter;
26 import jalview.viewmodel.annotationfilter.AnnotationFilterParameter.SearchableAnnotationField;
28 import java.util.ArrayList;
29 import java.util.BitSet;
30 import java.util.Collections;
31 import java.util.List;
32 import java.util.Vector;
35 * Data class holding the selected columns and hidden column ranges for a view.
38 public class ColumnSelection
41 * A class to hold an efficient representation of selected columns
46 * list of selected columns (ordered by selection order, not column order)
48 private List<Integer> order;
51 * an unmodifiable view of the selected columns list
53 private List<Integer> _uorder;
56 * bitfield for column selection - allows quick lookup
58 private BitSet selected;
65 order = new ArrayList<Integer>();
66 _uorder = Collections.unmodifiableList(order);
67 selected = new BitSet();
75 IntList(IntList other)
81 for (int i = 0; i < j; i++)
83 add(other.elementAt(i));
89 * adds a new column i to the selection - only if i is not already selected
97 order.add(Integer.valueOf(i));
111 Integer colInt = new Integer(col);
113 if (selected.get(col))
115 // if this ever changes to List.remove(), ensure Integer not int
117 // as List.remove(int i) removes the i'th item which is wrong
118 order.remove(colInt);
123 boolean contains(Integer colInt)
125 return selected.get(colInt);
130 return order.isEmpty();
134 * Returns a read-only view of the selected columns list
138 List<Integer> getList()
149 * gets the column that was selected first, second or i'th
159 protected boolean pruneColumnList(final List<int[]> shifts)
161 int s = 0, t = shifts.size();
162 int[] sr = shifts.get(s++);
163 boolean pruned = false;
164 int i = 0, j = order.size();
165 while (i < j && s <= t)
167 int c = order.get(i++).intValue();
170 if (sr[1] + sr[0] >= c)
171 { // sr[1] -ve means inseriton.
190 * shift every selected column at or above start by change
193 * - leftmost column to be shifted
197 void compensateForEdits(int start, int change)
199 BitSet mask = new BitSet();
200 for (int i = 0; i < order.size(); i++)
202 int temp = order.get(i);
206 // clear shifted bits and update List of selected columns
207 selected.clear(temp);
208 mask.set(temp - change);
209 order.set(i, new Integer(temp - change));
212 // lastly update the bitfield all at once
216 boolean isSelected(int column)
218 return selected.get(column);
223 return selected.length() - 1;
228 return selected.get(0) ? 0 : selected.nextSetBit(0);
232 * @return a series of selection intervals along the range
234 List<int[]> getRanges()
236 List<int[]> rlist = new ArrayList<int[]>();
237 if (selected.isEmpty())
241 int next = selected.nextSetBit(0), clear = -1;
244 clear = selected.nextClearBit(next);
245 rlist.add(new int[] { next, clear - 1 });
246 next = selected.nextSetBit(clear);
252 public int hashCode()
254 // TODO Auto-generated method stub
255 return selected.hashCode();
259 public boolean equals(Object obj)
261 if (obj instanceof IntList)
263 return ((IntList) obj).selected.equals(selected);
269 IntList selection = new IntList();
272 * list of hidden column [start, end] ranges; the list is maintained in
273 * ascending start column order
275 Vector<int[]> hiddenColumns;
278 * Add a column to the selection
283 public void addElement(int col)
289 * clears column selection
297 * Removes value 'col' from the selection (not the col'th item)
300 * index of column to be removed
302 public void removeElement(int col)
304 selection.remove(col);
308 * removes a range of columns from the selection
311 * int - first column in range to be removed
315 public void removeElements(int start, int end)
318 for (int i = start; i < end; i++)
320 colInt = new Integer(i);
321 if (selection.contains(colInt))
323 selection.remove(colInt);
329 * Returns a read-only view of the (possibly empty) list of selected columns
331 * The list contains no duplicates but is not necessarily ordered. It also may
332 * include columns hidden from the current view. To modify (for example sort)
333 * the list, you should first make a copy.
335 * The list is not thread-safe: iterating over it could result in
336 * ConcurrentModificationException if it is modified by another thread.
338 public List<Integer> getSelected()
340 return selection.getList();
344 * @return list of int arrays containing start and end column position for
345 * runs of selected columns ordered from right to left.
347 public List<int[]> getSelectedRanges()
349 return selection.getRanges();
355 * index to search for in column selection
357 * @return true if col is selected
359 public boolean contains(int col)
361 return (col > -1) ? selection.isSelected(col) : false;
365 * Answers true if no columns are selected, else false
367 public boolean isEmpty()
369 return selection == null || selection.isEmpty();
373 * rightmost selected column
375 * @return rightmost column in alignment that is selected
379 if (selection.isEmpty())
383 return selection.getMaxColumn();
387 * Leftmost column in selection
389 * @return column index of leftmost column in selection
393 if (selection.isEmpty())
397 return selection.getMinColumn();
401 * propagate shift in alignment columns to column selection
406 * shift in edit (+ve for removal, or -ve for inserts)
408 public List<int[]> compensateForEdit(int start, int change)
410 List<int[]> deletedHiddenColumns = null;
411 selection.compensateForEdits(start, change);
413 if (hiddenColumns != null)
415 deletedHiddenColumns = new ArrayList<int[]>();
416 int hSize = hiddenColumns.size();
417 for (int i = 0; i < hSize; i++)
419 int[] region = hiddenColumns.elementAt(i);
420 if (region[0] > start && start + change > region[1])
422 deletedHiddenColumns.add(region);
424 hiddenColumns.removeElementAt(i);
430 if (region[0] > start)
443 this.revealHiddenColumns(0);
446 return deletedHiddenColumns;
450 * propagate shift in alignment columns to column selection special version of
451 * compensateForEdit - allowing for edits within hidden regions
456 * shift in edit (+ve for removal, or -ve for inserts)
458 private void compensateForDelEdits(int start, int change)
461 selection.compensateForEdits(start, change);
463 if (hiddenColumns != null)
465 for (int i = 0; i < hiddenColumns.size(); i++)
467 int[] region = hiddenColumns.elementAt(i);
468 if (region[0] >= start)
472 if (region[1] >= start)
476 if (region[1] < region[0])
478 hiddenColumns.removeElementAt(i--);
494 * Adjust hidden column boundaries based on a series of column additions or
495 * deletions in visible regions.
500 public ShiftList compensateForEdits(ShiftList shiftrecord)
502 if (shiftrecord != null)
504 final List<int[]> shifts = shiftrecord.getShifts();
505 if (shifts != null && shifts.size() > 0)
508 for (int i = 0, j = shifts.size(); i < j; i++)
510 int[] sh = shifts.get(i);
511 // compensateForEdit(shifted+sh[0], sh[1]);
512 compensateForDelEdits(shifted + sh[0], sh[1]);
516 return shiftrecord.getInverse();
522 * removes intersection of position,length ranges in deletions from the
523 * start,end regions marked in intervals.
529 private boolean pruneIntervalVector(final List<int[]> shifts,
530 Vector<int[]> intervals)
532 boolean pruned = false;
533 int i = 0, j = intervals.size() - 1, s = 0, t = shifts.size() - 1;
534 int hr[] = intervals.elementAt(i);
535 int sr[] = shifts.get(s);
536 while (i <= j && s <= t)
538 boolean trailinghn = hr[1] >= sr[0];
543 hr = intervals.elementAt(++i);
551 int endshift = sr[0] + sr[1]; // deletion ranges - -ve means an insert
552 if (endshift < hr[0] || endshift < sr[0])
553 { // leadinghc disjoint or not a deletion
556 sr = shifts.get(++s);
564 boolean leadinghn = hr[0] >= sr[0];
565 boolean leadinghc = hr[0] < endshift;
566 boolean trailinghc = hr[1] < endshift;
570 { // deleted hidden region.
571 intervals.removeElementAt(i);
576 hr = intervals.elementAt(i);
582 hr[0] = endshift; // clip c terminal region
583 leadinghn = !leadinghn;
599 // sr contained in hr
602 sr = shifts.get(++s);
612 return pruned; // true if any interval was removed or modified by
617 * remove any hiddenColumns or selected columns and shift remaining based on a
618 * series of position, range deletions.
622 public void pruneDeletions(ShiftList deletions)
624 if (deletions != null)
626 final List<int[]> shifts = deletions.getShifts();
627 if (shifts != null && shifts.size() > 0)
629 // delete any intervals intersecting.
630 if (hiddenColumns != null)
632 pruneIntervalVector(shifts, hiddenColumns);
633 if (hiddenColumns != null && hiddenColumns.size() == 0)
635 hiddenColumns = null;
638 if (selection != null && selection.size() > 0)
640 selection.pruneColumnList(shifts);
641 if (selection != null && selection.size() == 0)
646 // and shift the rest.
647 this.compensateForEdits(deletions);
653 * This Method is used to return all the HiddenColumn regions
655 * @return empty list or List of hidden column intervals
657 public List<int[]> getHiddenColumns()
659 return hiddenColumns == null ? Collections.<int[]> emptyList()
664 * Return absolute column index for a visible column index
667 * int column index in alignment view (count from zero)
668 * @return alignment column index for column
670 public int adjustForHiddenColumns(int column)
673 if (hiddenColumns != null)
675 for (int i = 0; i < hiddenColumns.size(); i++)
677 int[] region = hiddenColumns.elementAt(i);
678 if (result >= region[0])
680 result += region[1] - region[0] + 1;
688 * Use this method to find out where a column will appear in the visible
689 * alignment when hidden columns exist. If the column is not visible, then the
690 * left-most visible column will always be returned.
692 * @param hiddenColumn
693 * the column index in the full alignment including hidden columns
694 * @return the position of the column in the visible alignment
696 public int findColumnPosition(int hiddenColumn)
698 int result = hiddenColumn;
699 if (hiddenColumns != null)
705 region = hiddenColumns.elementAt(index++);
706 if (hiddenColumn > region[1])
708 result -= region[1] + 1 - region[0];
710 } while ((hiddenColumn > region[1]) && (index < hiddenColumns.size()));
712 if (hiddenColumn >= region[0] && hiddenColumn <= region[1])
714 // Here the hidden column is within a region, so
715 // we want to return the position of region[0]-1, adjusted for any
716 // earlier hidden columns.
717 // Calculate the difference between the actual hidden col position
718 // and region[0]-1, and then subtract from result to convert result from
719 // the adjusted hiddenColumn value to the adjusted region[0]-1 value
721 return result - (hiddenColumn - region[0] + 1);
724 return result; // return the shifted position after removing hidden columns.
728 * Use this method to determine where the next hiddenRegion starts
730 * @param hiddenRegion
731 * index of hidden region (counts from 0)
732 * @return column number in visible view
734 public int findHiddenRegionPosition(int hiddenRegion)
737 if (hiddenColumns != null)
743 int[] region = hiddenColumns.elementAt(index);
744 if (hiddenRegion == 0)
749 gaps += region[1] + 1 - region[0];
750 result = region[1] + 1;
752 } while (index <= hiddenRegion);
761 * THis method returns the rightmost limit of a region of an alignment with
762 * hidden columns. In otherwords, the next hidden column.
767 public int getHiddenBoundaryRight(int alPos)
769 if (hiddenColumns != null)
774 int[] region = hiddenColumns.elementAt(index);
775 if (alPos < region[0])
781 } while (index < hiddenColumns.size());
789 * This method returns the leftmost limit of a region of an alignment with
790 * hidden columns. In otherwords, the previous hidden column.
795 public int getHiddenBoundaryLeft(int alPos)
797 if (hiddenColumns != null)
799 int index = hiddenColumns.size() - 1;
802 int[] region = hiddenColumns.elementAt(index);
803 if (alPos > region[1])
809 } while (index > -1);
816 public void hideSelectedColumns()
818 synchronized (selection)
820 for (int[] selregions : selection.getRanges())
822 hideColumns(selregions[0], selregions[1]);
830 * Adds the specified column range to the hidden columns
835 public void hideColumns(int start, int end)
837 if (hiddenColumns == null)
839 hiddenColumns = new Vector<int[]>();
843 * traverse existing hidden ranges and insert / amend / append as
846 for (int i = 0; i < hiddenColumns.size(); i++)
848 int[] region = hiddenColumns.elementAt(i);
850 if (end < region[0] - 1)
853 * insert discontiguous preceding range
855 hiddenColumns.insertElementAt(new int[] { start, end }, i);
859 if (end <= region[1])
862 * new range overlaps existing, or is contiguous preceding it - adjust
865 region[0] = Math.min(region[0], start);
869 if (start <= region[1] + 1)
872 * new range overlaps existing, or is contiguous following it - adjust
873 * start and end columns
875 region[0] = Math.min(region[0], start);
876 region[1] = Math.max(region[1], end);
879 * also update or remove any subsequent ranges
880 * that are overlapped
882 while (i < hiddenColumns.size() - 1)
884 int[] nextRegion = hiddenColumns.get(i + 1);
885 if (nextRegion[0] > end + 1)
888 * gap to next hidden range - no more to update
892 region[1] = Math.max(nextRegion[1], end);
893 hiddenColumns.remove(i + 1);
900 * remaining case is that the new range follows everything else
902 hiddenColumns.addElement(new int[] { start, end });
906 * Hides the specified column and any adjacent selected columns
911 public void hideColumns(int col)
914 * deselect column (whether selected or not!)
919 * find adjacent selected columns
921 int min = col - 1, max = col + 1;
922 while (contains(min))
928 while (contains(max))
935 * min, max are now the closest unselected columns
944 hideColumns(min, max);
948 * Unhides, and adds to the selection list, all hidden columns
950 public void revealAllHiddenColumns()
952 if (hiddenColumns != null)
954 for (int i = 0; i < hiddenColumns.size(); i++)
956 int[] region = hiddenColumns.elementAt(i);
957 for (int j = region[0]; j < region[1] + 1; j++)
964 hiddenColumns = null;
968 * Reveals, and marks as selected, the hidden column range with the given
973 public void revealHiddenColumns(int start)
975 for (int i = 0; i < hiddenColumns.size(); i++)
977 int[] region = hiddenColumns.elementAt(i);
978 if (start == region[0])
980 for (int j = region[0]; j < region[1] + 1; j++)
985 hiddenColumns.removeElement(region);
989 if (hiddenColumns.size() == 0)
991 hiddenColumns = null;
995 public boolean isVisible(int column)
997 if (hiddenColumns != null)
999 for (int[] region : hiddenColumns)
1001 if (column >= region[0] && column <= region[1])
1016 public ColumnSelection(ColumnSelection copy)
1020 selection = new IntList(copy.selection);
1021 if (copy.hiddenColumns != null)
1023 hiddenColumns = new Vector<int[]>(copy.hiddenColumns.size());
1024 for (int i = 0, j = copy.hiddenColumns.size(); i < j; i++)
1027 rh = copy.hiddenColumns.elementAt(i);
1030 cp = new int[rh.length];
1031 System.arraycopy(rh, 0, cp, 0, rh.length);
1032 hiddenColumns.addElement(cp);
1042 public ColumnSelection()
1046 public String[] getVisibleSequenceStrings(int start, int end,
1049 int i, iSize = seqs.length;
1050 String selections[] = new String[iSize];
1051 if (hiddenColumns != null && hiddenColumns.size() > 0)
1053 for (i = 0; i < iSize; i++)
1055 StringBuffer visibleSeq = new StringBuffer();
1056 List<int[]> regions = getHiddenColumns();
1058 int blockStart = start, blockEnd = end;
1060 int hideStart, hideEnd;
1062 for (int j = 0; j < regions.size(); j++)
1064 region = regions.get(j);
1065 hideStart = region[0];
1066 hideEnd = region[1];
1068 if (hideStart < start)
1073 blockStart = Math.min(blockStart, hideEnd + 1);
1074 blockEnd = Math.min(blockEnd, hideStart);
1076 if (blockStart > blockEnd)
1081 visibleSeq.append(seqs[i].getSequence(blockStart, blockEnd));
1083 blockStart = hideEnd + 1;
1087 if (end > blockStart)
1089 visibleSeq.append(seqs[i].getSequence(blockStart, end));
1092 selections[i] = visibleSeq.toString();
1097 for (i = 0; i < iSize; i++)
1099 selections[i] = seqs[i].getSequenceAsString(start, end);
1107 * return all visible segments between the given start and end boundaries
1110 * (first column inclusive from 0)
1112 * (last column - not inclusive)
1113 * @return int[] {i_start, i_end, ..} where intervals lie in
1114 * start<=i_start<=i_end<end
1116 public int[] getVisibleContigs(int start, int end)
1118 if (hiddenColumns != null && hiddenColumns.size() > 0)
1120 List<int[]> visiblecontigs = new ArrayList<int[]>();
1121 List<int[]> regions = getHiddenColumns();
1125 int hideStart, hideEnd;
1127 for (int j = 0; vstart < end && j < regions.size(); j++)
1129 region = regions.get(j);
1130 hideStart = region[0];
1131 hideEnd = region[1];
1133 if (hideEnd < vstart)
1137 if (hideStart > vstart)
1139 visiblecontigs.add(new int[] { vstart, hideStart - 1 });
1141 vstart = hideEnd + 1;
1146 visiblecontigs.add(new int[] { vstart, end - 1 });
1148 int[] vcontigs = new int[visiblecontigs.size() * 2];
1149 for (int i = 0, j = visiblecontigs.size(); i < j; i++)
1151 int[] vc = visiblecontigs.get(i);
1152 visiblecontigs.set(i, null);
1153 vcontigs[i * 2] = vc[0];
1154 vcontigs[i * 2 + 1] = vc[1];
1156 visiblecontigs.clear();
1161 return new int[] { start, end - 1 };
1166 * Locate the first and last position visible for this sequence. if seq isn't
1167 * visible then return the position of the left and right of the hidden
1168 * boundary region, and the corresponding alignment column indices for the
1169 * extent of the sequence
1172 * @return int[] { visible start, visible end, first seqpos, last seqpos,
1173 * alignment index for seq start, alignment index for seq end }
1175 public int[] locateVisibleBoundsOfSequence(SequenceI seq)
1177 int fpos = seq.getStart(), lpos = seq.getEnd();
1180 if (hiddenColumns == null || hiddenColumns.size() == 0)
1182 int ifpos = seq.findIndex(fpos) - 1, ilpos = seq.findIndex(lpos) - 1;
1183 return new int[] { ifpos, ilpos, fpos, lpos, ifpos, ilpos };
1186 // Simply walk along the sequence whilst watching for hidden column
1188 List<int[]> regions = getHiddenColumns();
1189 int spos = fpos, lastvispos = -1, rcount = 0, hideStart = seq
1190 .getLength(), hideEnd = -1;
1191 int visPrev = 0, visNext = 0, firstP = -1, lastP = -1;
1192 boolean foundStart = false;
1193 for (int p = 0, pLen = seq.getLength(); spos <= seq.getEnd()
1196 if (!Comparison.isGap(seq.getCharAt(p)))
1198 // keep track of first/last column
1199 // containing sequence data regardless of visibility
1205 // update hidden region start/end
1206 while (hideEnd < p && rcount < regions.size())
1208 int[] region = regions.get(rcount++);
1210 visNext += region[0] - visPrev;
1211 hideStart = region[0];
1212 hideEnd = region[1];
1216 hideStart = seq.getLength();
1218 // update visible boundary for sequence
1230 // look for next sequence position
1236 return new int[] { findColumnPosition(start),
1237 findColumnPosition(lastvispos), fpos, lpos, firstP, lastP };
1239 // otherwise, sequence was completely hidden
1240 return new int[] { visPrev, visNext, 0, 0, firstP, lastP };
1244 * delete any columns in alignmentAnnotation that are hidden (including
1245 * sequence associated annotation).
1247 * @param alignmentAnnotation
1249 public void makeVisibleAnnotation(AlignmentAnnotation alignmentAnnotation)
1251 makeVisibleAnnotation(-1, -1, alignmentAnnotation);
1255 * delete any columns in alignmentAnnotation that are hidden (including
1256 * sequence associated annotation).
1259 * remove any annotation to the right of this column
1261 * remove any annotation to the left of this column
1262 * @param alignmentAnnotation
1263 * the annotation to operate on
1265 public void makeVisibleAnnotation(int start, int end,
1266 AlignmentAnnotation alignmentAnnotation)
1268 if (alignmentAnnotation.annotations == null)
1272 if (start == end && end == -1)
1275 end = alignmentAnnotation.annotations.length;
1277 if (hiddenColumns != null && hiddenColumns.size() > 0)
1279 // then mangle the alignmentAnnotation annotation array
1280 Vector<Annotation[]> annels = new Vector<Annotation[]>();
1281 Annotation[] els = null;
1282 List<int[]> regions = getHiddenColumns();
1283 int blockStart = start, blockEnd = end;
1285 int hideStart, hideEnd, w = 0;
1287 for (int j = 0; j < regions.size(); j++)
1289 region = regions.get(j);
1290 hideStart = region[0];
1291 hideEnd = region[1];
1293 if (hideStart < start)
1298 blockStart = Math.min(blockStart, hideEnd + 1);
1299 blockEnd = Math.min(blockEnd, hideStart);
1301 if (blockStart > blockEnd)
1306 annels.addElement(els = new Annotation[blockEnd - blockStart]);
1307 System.arraycopy(alignmentAnnotation.annotations, blockStart, els,
1310 blockStart = hideEnd + 1;
1314 if (end > blockStart)
1316 annels.addElement(els = new Annotation[end - blockStart + 1]);
1317 if ((els.length + blockStart) <= alignmentAnnotation.annotations.length)
1319 // copy just the visible segment of the annotation row
1320 System.arraycopy(alignmentAnnotation.annotations, blockStart,
1321 els, 0, els.length);
1325 // copy to the end of the annotation row
1326 System.arraycopy(alignmentAnnotation.annotations, blockStart,
1328 (alignmentAnnotation.annotations.length - blockStart));
1337 alignmentAnnotation.annotations = new Annotation[w];
1340 for (Annotation[] chnk : annels)
1342 System.arraycopy(chnk, 0, alignmentAnnotation.annotations, w,
1349 alignmentAnnotation.restrict(start, end);
1354 * Invert the column selection from first to end-1. leaves hiddenColumns
1355 * untouched (and unselected)
1360 public void invertColumnSelection(int first, int width)
1362 boolean hasHidden = hiddenColumns != null && hiddenColumns.size() > 0;
1363 for (int i = first; i < width; i++)
1371 if (!hasHidden || isVisible(i))
1380 * add in any unselected columns from the given column selection, excluding
1381 * any that are hidden.
1385 public void addElementsFrom(ColumnSelection colsel)
1387 if (colsel != null && !colsel.isEmpty())
1389 for (Integer col : colsel.getSelected())
1391 if (hiddenColumns != null && isVisible(col.intValue()))
1400 * set the selected columns the given column selection, excluding any columns
1405 public void setElementsFrom(ColumnSelection colsel)
1407 selection = new IntList();
1408 if (colsel.selection != null && colsel.selection.size() > 0)
1410 if (hiddenColumns != null && hiddenColumns.size() > 0)
1412 // only select visible columns in this columns selection
1413 addElementsFrom(colsel);
1417 // add everything regardless
1418 for (Integer col : colsel.getSelected())
1427 * Add gaps into the sequences aligned to profileseq under the given
1432 * - alignment to have gaps inserted into it
1434 * - alignment view where sequence corresponding to profileseq is
1436 * @return new Column selection for new alignment view, with insertions into
1437 * profileseq marked as hidden.
1439 public static ColumnSelection propagateInsertions(SequenceI profileseq,
1440 AlignmentI al, AlignmentView input)
1444 // return propagateInsertions(profileseq, al, )
1445 char gc = al.getGapCharacter();
1446 Object[] alandcolsel = input.getAlignmentAndColumnSelection(gc);
1447 ColumnSelection nview = (ColumnSelection) alandcolsel[1];
1448 SequenceI origseq = ((SequenceI[]) alandcolsel[0])[profsqpos];
1449 nview.propagateInsertions(profileseq, al, origseq);
1456 * - sequence in al which corresponds to origseq
1458 * - alignment which is to have gaps inserted into it
1460 * - sequence corresponding to profileseq which defines gap map for
1463 public void propagateInsertions(SequenceI profileseq, AlignmentI al,
1466 char gc = al.getGapCharacter();
1467 // recover mapping between sequence's non-gap positions and positions
1469 pruneDeletions(ShiftList.parseMap(origseq.gapMap()));
1470 int[] viscontigs = getVisibleContigs(0, profileseq.getLength());
1473 // input.pruneDeletions(ShiftList.parseMap(((SequenceI[])
1474 // alandcolsel[0])[0].gapMap()))
1475 // add profile to visible contigs
1476 for (int v = 0; v < viscontigs.length; v += 2)
1478 if (viscontigs[v] > spos)
1480 StringBuffer sb = new StringBuffer();
1481 for (int s = 0, ns = viscontigs[v] - spos; s < ns; s++)
1485 for (int s = 0, ns = al.getHeight(); s < ns; s++)
1487 SequenceI sqobj = al.getSequenceAt(s);
1488 if (sqobj != profileseq)
1490 String sq = al.getSequenceAt(s).getSequenceAsString();
1491 if (sq.length() <= spos + offset)
1494 int diff = spos + offset - sq.length() - 1;
1499 while ((diff = spos + offset - sq.length() - 1) > 0)
1502 // + ((diff >= sb.length()) ? sb.toString() : sb
1503 // .substring(0, diff));
1504 if (diff >= sb.length())
1506 sq += sb.toString();
1510 char[] buf = new char[diff];
1511 sb.getChars(0, diff, buf, 0);
1512 sq += buf.toString();
1516 sq += sb.toString();
1520 al.getSequenceAt(s).setSequence(
1521 sq.substring(0, spos + offset) + sb.toString()
1522 + sq.substring(spos + offset));
1526 // offset+=sb.length();
1528 spos = viscontigs[v + 1] + 1;
1530 if ((offset + spos) < profileseq.getLength())
1532 // pad the final region with gaps.
1533 StringBuffer sb = new StringBuffer();
1534 for (int s = 0, ns = profileseq.getLength() - spos - offset; s < ns; s++)
1538 for (int s = 0, ns = al.getHeight(); s < ns; s++)
1540 SequenceI sqobj = al.getSequenceAt(s);
1541 if (sqobj == profileseq)
1545 String sq = sqobj.getSequenceAsString();
1547 int diff = origseq.getLength() - sq.length();
1551 // + ((diff >= sb.length()) ? sb.toString() : sb
1552 // .substring(0, diff));
1553 if (diff >= sb.length())
1555 sq += sb.toString();
1559 char[] buf = new char[diff];
1560 sb.getChars(0, diff, buf, 0);
1561 sq += buf.toString();
1563 diff = origseq.getLength() - sq.length();
1571 * @return true if there are columns marked
1573 public boolean hasSelectedColumns()
1575 return (selection != null && selection.size() > 0);
1580 * @return true if there are columns hidden
1582 public boolean hasHiddenColumns()
1584 return hiddenColumns != null && hiddenColumns.size() > 0;
1589 * @return true if there are more than one set of columns hidden
1591 public boolean hasManyHiddenColumns()
1593 return hiddenColumns != null && hiddenColumns.size() > 1;
1597 * mark the columns corresponding to gap characters as hidden in the column
1602 public void hideInsertionsFor(SequenceI sr)
1604 List<int[]> inserts = sr.getInsertions();
1605 for (int[] r : inserts)
1607 hideColumns(r[0], r[1]);
1611 public boolean filterAnnotations(Annotation[] annotations,
1612 AnnotationFilterParameter filterParams)
1614 // JBPNote - this method needs to be refactored to become independent of
1615 // viewmodel package
1616 this.revealAllHiddenColumns();
1621 if (annotations[count] != null)
1624 boolean itemMatched = false;
1626 if (filterParams.getThresholdType() == AnnotationFilterParameter.ThresholdType.ABOVE_THRESHOLD
1627 && annotations[count].value >= filterParams
1628 .getThresholdValue())
1632 if (filterParams.getThresholdType() == AnnotationFilterParameter.ThresholdType.BELOW_THRESHOLD
1633 && annotations[count].value <= filterParams
1634 .getThresholdValue())
1639 if (filterParams.isFilterAlphaHelix()
1640 && annotations[count].secondaryStructure == 'H')
1645 if (filterParams.isFilterBetaSheet()
1646 && annotations[count].secondaryStructure == 'E')
1651 if (filterParams.isFilterTurn()
1652 && annotations[count].secondaryStructure == 'S')
1657 String regexSearchString = filterParams.getRegexString();
1658 if (regexSearchString != null
1659 && !filterParams.getRegexSearchFields().isEmpty())
1661 List<SearchableAnnotationField> fields = filterParams
1662 .getRegexSearchFields();
1665 if (fields.contains(SearchableAnnotationField.DISPLAY_STRING)
1666 && annotations[count].displayCharacter
1667 .matches(regexSearchString))
1671 } catch (java.util.regex.PatternSyntaxException pse)
1673 if (annotations[count].displayCharacter
1674 .equals(regexSearchString))
1679 if (fields.contains(SearchableAnnotationField.DESCRIPTION)
1680 && annotations[count].description != null
1681 && annotations[count].description
1682 .matches(regexSearchString))
1690 this.addElement(count);
1694 } while (count < annotations.length);
1699 * Returns a hashCode built from selected columns and hidden column ranges
1702 public int hashCode()
1704 int hashCode = selection.hashCode();
1705 if (hiddenColumns != null)
1707 for (int[] hidden : hiddenColumns)
1709 hashCode = 31 * hashCode + hidden[0];
1710 hashCode = 31 * hashCode + hidden[1];
1717 * Answers true if comparing to a ColumnSelection with the same selected
1718 * columns and hidden columns, else false
1721 public boolean equals(Object obj)
1723 if (!(obj instanceof ColumnSelection))
1727 ColumnSelection that = (ColumnSelection) obj;
1730 * check columns selected are either both null, or match
1732 if (this.selection == null)
1734 if (that.selection != null)
1739 if (!this.selection.equals(that.selection))
1745 * check hidden columns are either both null, or match
1747 if (this.hiddenColumns == null)
1749 return (that.hiddenColumns == null);
1751 if (that.hiddenColumns == null
1752 || that.hiddenColumns.size() != this.hiddenColumns.size())
1757 for (int[] thisRange : hiddenColumns)
1759 int[] thatRange = that.hiddenColumns.get(i++);
1760 if (thisRange[0] != thatRange[0] || thisRange[1] != thatRange[1])
1769 * Updates the column selection depending on the parameters, and returns true
1770 * if any change was made to the selection
1772 * @param markedColumns
1773 * a set identifying marked columns (base 0)
1775 * the first column of the range to operate over (base 0)
1777 * the last column of the range to operate over (base 0)
1779 * if true, deselect marked columns and select unmarked
1780 * @param extendCurrent
1781 * if true, extend rather than replacing the current column selection
1783 * if true, toggle the selection state of marked columns
1787 public boolean markColumns(BitSet markedColumns, int startCol,
1788 int endCol, boolean invert, boolean extendCurrent, boolean toggle)
1790 boolean changed = false;
1791 if (!extendCurrent && !toggle)
1793 changed = !this.isEmpty();
1798 // invert only in the currently selected sequence region
1799 int i = markedColumns.nextClearBit(startCol);
1800 int ibs = markedColumns.nextSetBit(startCol);
1801 while (i >= startCol && i <= endCol)
1803 if (ibs < 0 || i < ibs)
1806 if (toggle && contains(i))
1817 i = markedColumns.nextClearBit(ibs);
1818 ibs = markedColumns.nextSetBit(i);
1824 int i = markedColumns.nextSetBit(startCol);
1825 while (i >= startCol && i <= endCol)
1828 if (toggle && contains(i))
1836 i = markedColumns.nextSetBit(i + 1);
1843 * Adjusts column selections, and the given selection group, to match the
1844 * range of a stretch (e.g. mouse drag) operation
1846 * Method refactored from ScalePanel.mouseDragged
1849 * current column position, adjusted for hidden columns
1851 * current selection group
1853 * start position of the stretch group
1855 * end position of the stretch group
1857 public void stretchGroup(int res, SequenceGroup sg, int min, int max)
1864 if (res > sg.getStartRes())
1866 // expand selection group to the right
1869 if (res < sg.getStartRes())
1871 // expand selection group to the left
1872 sg.setStartRes(res);
1876 * expand or shrink column selection to match the
1877 * range of the drag operation
1879 for (int col = min; col <= max; col++)
1881 if (col < sg.getStartRes() || col > sg.getEndRes())
1883 // shrinking drag - remove from selection
1888 // expanding drag - add to selection