X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Fanalysis%2FSequenceIdMatcher.java;h=3d4cbe7963fb04cb003e543de34943c0441fa923;hb=3d0101179759ef157b088ea135423cd909512d9f;hp=5bde22598ba356163e2735c7fadb9c844bdee162;hpb=f5db0ff259cbb1132f6c061615d912b0668fd728;p=jalview.git diff --git a/src/jalview/analysis/SequenceIdMatcher.java b/src/jalview/analysis/SequenceIdMatcher.java index 5bde225..3d4cbe7 100755 --- a/src/jalview/analysis/SequenceIdMatcher.java +++ b/src/jalview/analysis/SequenceIdMatcher.java @@ -1,92 +1,121 @@ /* - * Jalview - A Sequence Alignment Editor and Viewer (Version 2.6) - * Copyright (C) 2010 J Procter, AM Waterhouse, G Barton, M Clamp, S Searle + * Jalview - A Sequence Alignment Editor and Viewer ($$Version-Rel$$) + * Copyright (C) $$Year-Rel$$ The Jalview Authors * * This file is part of Jalview. * * Jalview is free software: you can redistribute it and/or * modify it under the terms of the GNU General Public License - * as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version. - * + * as published by the Free Software Foundation, either version 3 + * of the License, or (at your option) any later version. + * * Jalview is distributed in the hope that it will be useful, but * WITHOUT ANY WARRANTY; without even the implied warranty * of MERCHANTABILITY or FITNESS FOR A PARTICULAR * PURPOSE. See the GNU General Public License for more details. * - * You should have received a copy of the GNU General Public License along with Jalview. If not, see . + * You should have received a copy of the GNU General Public License + * along with Jalview. If not, see . + * The Jalview Authors are detailed in the 'AUTHORS' file. */ package jalview.analysis; -import java.util.*; +import jalview.datamodel.DBRefEntry; +import jalview.datamodel.SequenceI; -import jalview.datamodel.*; +import java.util.ArrayList; +import java.util.Arrays; +import java.util.HashMap; +import java.util.List; +import java.util.Vector; /** - *

- * Title: - *

- * SequenceIdMatcher - *

- * Description: - *

- * Routine which does approximate Sequence Id resolution by name using string + * Routines for approximate Sequence Id resolution by name using string * containment (on word boundaries) rather than equivalence. It also attempts to * resolve ties where no exact match is available by picking the the id closest * to the query. - *

- * Copyright: Copyright (c) 2004 - *

- * - *

- * Company: Dundee University - *

- * - * @author not attributable - * @version 1.0 */ public class SequenceIdMatcher { - private Hashtable names; + private HashMap names; - public SequenceIdMatcher(SequenceI[] seqs) + public SequenceIdMatcher(List seqs) + { + names = new HashMap(); + addAll(seqs); + } + + /** + * Adds sequences to this matcher + * + * @param seqs + */ + public void addAll(List seqs) { - names = new Hashtable(); - for (int i = 0; i < seqs.length; i++) + for (SequenceI seq : seqs) { - // TODO: deal with ID collisions - SequenceI should be appended to list associated with this key. - names.put(new SeqIdName(seqs[i].getDisplayId(true)), seqs[i]); - // add in any interesting identifiers - if (seqs[i].getDBRef() != null) + add(seq); + } + } + + /** + * Adds one sequence to this matcher + * + * @param seq + */ + public void add(SequenceI seq) + { + // TODO: deal with ID collisions - SequenceI should be appended to list + // associated with this key. + names.put(new SeqIdName(seq.getDisplayId(true)), seq); + SequenceI dbseq = seq; + while (dbseq.getDatasetSequence() != null) + { + dbseq = dbseq.getDatasetSequence(); + } + // add in any interesting identifiers + if (dbseq.getDBRefs() != null) + { + DBRefEntry dbr[] = dbseq.getDBRefs(); + SeqIdName sid = null; + for (int r = 0; r < dbr.length; r++) { - DBRefEntry dbr[] = seqs[i].getDBRef(); - SeqIdName sid = null; - for (int r = 0; r < dbr.length; r++) + sid = new SeqIdName(dbr[r].getAccessionId()); + if (!names.containsKey(sid)) { - sid = new SeqIdName(dbr[r].getAccessionId()); - if (!names.contains(sid)) - { - names.put(sid, seqs[i]); - } + names.put(sid, seq); } } } } /** + * convenience method to make a matcher from concrete array + * + * @param sequences + */ + public SequenceIdMatcher(SequenceI[] sequences) + { + this(Arrays.asList(sequences)); + } + + /** * returns the closest SequenceI in matches to SeqIdName and returns all the * matches to the names hash. * * @param candName * SeqIdName * @param matches - * Vector of SequenceI objects + * List of SequenceI objects * @return SequenceI closest SequenceI to SeqIdName */ - private SequenceI pickbestMatch(SeqIdName candName, Vector matches) + private SequenceI pickbestMatch(SeqIdName candName, + List matches) { - SequenceI[] st= pickbestMatches(candName, matches); - return st==null || st.length==0 ? null : st[0]; + List st = pickbestMatches(candName, matches); + return st == null || st.size() == 0 ? null : st.get(0); } + /** * returns the closest SequenceI in matches to SeqIdName and returns all the * matches to the names hash. @@ -95,18 +124,18 @@ public class SequenceIdMatcher * SeqIdName * @param matches * Vector of SequenceI objects - * @return Object[] { SequenceI closest SequenceI to SeqIdName, SequenceI[] ties } + * @return Object[] { SequenceI closest SequenceI to SeqIdName, SequenceI[] + * ties } */ - private SequenceI[] pickbestMatches(SeqIdName candName, Vector matches) + private List pickbestMatches(SeqIdName candName, + List matches) { - ArrayList best=new ArrayList(); - SequenceI match = null; + ArrayList best = new ArrayList(); if (candName == null || matches == null || matches.size() == 0) { return null; } - match = (SequenceI) matches.elementAt(0); - matches.removeElementAt(0); + SequenceI match = matches.remove(0); best.add(match); names.put(new SeqIdName(match.getName()), match); int matchlen = match.getName().length(); @@ -114,27 +143,30 @@ public class SequenceIdMatcher while (matches.size() > 0) { // look through for a better one. - SequenceI cand = (SequenceI) matches.elementAt(0); - matches.remove(0); + SequenceI cand = matches.remove(0); names.put(new SeqIdName(cand.getName()), cand); - int q,w,candlen = cand.getName().length(); + int q, w, candlen = cand.getName().length(); // keep the one with an id 'closer' to the given seqnam string - if ((q=Math.abs(matchlen - namlen)) > (w=Math.abs(candlen - namlen)) - && candlen > matchlen) + if ((q = Math.abs(matchlen - namlen)) > (w = Math + .abs(candlen - namlen)) && candlen > matchlen) { best.clear(); match = cand; matchlen = candlen; best.add(match); } - if (q==w && candlen==matchlen) + if (q == w && candlen == matchlen) { // record any ties best.add(cand); } } - if (best.size()==0) { return null; }; - return (SequenceI[]) best.toArray(new SequenceI[0]); + if (best.size() == 0) + { + return null; + } + ; + return best; } /** @@ -157,14 +189,22 @@ public class SequenceIdMatcher } /** - * Find all matches for a given sequence name. - * @param seqnam string to query Matcher with. + * Find all matches for a given sequence name. + * + * @param seqnam + * string to query Matcher with. + * @return a new array or (possibly) null */ public SequenceI[] findAllIdMatches(String seqnam) { SeqIdName nam = new SeqIdName(seqnam); - return findAllIdMatches(nam); + List m = findAllIdMatches(nam); + if (m != null) + { + return m.toArray(new SequenceI[m.size()]); + } + return null; } /** @@ -221,6 +261,7 @@ public class SequenceIdMatcher } return pickbestMatch(nam, matches); } + /** * core findIdMatch search method for finding all equivalent matches * @@ -228,20 +269,19 @@ public class SequenceIdMatcher * SeqIdName * @return SequenceI[] */ - private SequenceI[] findAllIdMatches( + private List findAllIdMatches( jalview.analysis.SequenceIdMatcher.SeqIdName nam) { - Vector matches = new Vector(); + ArrayList matches = new ArrayList(); while (names.containsKey(nam)) { - matches.addElement(names.remove(nam)); + matches.add(names.remove(nam)); } - SequenceI[] r=pickbestMatches(nam, matches); + List r = pickbestMatches(nam, matches); return r; } - - private class SeqIdName + class SeqIdName { String id; @@ -249,7 +289,7 @@ public class SequenceIdMatcher { if (s != null) { - id = new String(s); + id = s.toLowerCase(); } else { @@ -257,23 +297,29 @@ public class SequenceIdMatcher } } + @Override public int hashCode() { - return ((id.length() >= 4) ? id.substring(0, 4).hashCode() : id - .hashCode()); + return ((id.length() >= 4) ? id.substring(0, 4).hashCode() + : id.hashCode()); } + @Override public boolean equals(Object s) { + if (s == null) + { + return false; + } if (s instanceof SeqIdName) { - return this.equals((SeqIdName) s); + return this.stringequals(((SeqIdName) s).id); } else { if (s instanceof String) { - return this.equals((String) s); + return this.stringequals(((String) s).toLowerCase()); } } @@ -286,7 +332,8 @@ public class SequenceIdMatcher * arbritrarily extended sequence id's (like portions of an aligned set of * repeats from one sequence) */ - private String WORD_SEP = "~. |#\\/<>!\""+((char)0x00A4)+"$%^*)}[@',?_"; + private String WORD_SEP = "~. |#\\/<>!\"" + ((char) 0x00A4) + + "$%^*)}[@',?_"; /** * matches if one ID properly contains another at a whitespace boundary. @@ -294,36 +341,33 @@ public class SequenceIdMatcher * todo: (JBPNote) Set separator characters appropriately * * @param s - * SeqIdName * @return boolean */ - public boolean equals(SeqIdName s) + private boolean stringequals(String s) { - // TODO: JAL-732 patch for cases when name includes a list of IDs, and the match contains one ID flanked - if (id.length() > s.id.length()) + if (id.length() > s.length()) { - return id.startsWith(s.id) ? (WORD_SEP.indexOf(id.charAt(s.id - .length())) > -1) : false; + return id.startsWith(s) + ? (WORD_SEP.indexOf(id.charAt(s.length())) > -1) + : false; } else { - return s.id.startsWith(id) ? (s.id.equals(id) ? true : (WORD_SEP - .indexOf(s.id.charAt(id.length())) > -1)) : false; + return s.startsWith(id) + ? (s.equals(id) ? true + : (WORD_SEP.indexOf(s.charAt(id.length())) > -1)) + : false; } } - public boolean equals(String s) + /** + * toString method returns the wrapped sequence id. For debugging purposes + * only, behaviour not guaranteed not to change. + */ + @Override + public String toString() { - if (id.length() > s.length()) - { - return id.startsWith(s) ? (WORD_SEP.indexOf(id.charAt(s.length())) > -1) - : false; - } - else - { - return s.startsWith(id) ? (s.equals(id) ? true : (WORD_SEP - .indexOf(s.charAt(id.length())) > -1)) : false; - } + return id; } } }