src/jalview/schemes/ClustalxColourScheme.java

   1 /* Jalview - a java multiple alignment editor\r
   2  * Copyright (C) 1998  Michele Clamp\r
   3  *\r
   4  * This program is free software; you can redistribute it and/or\r
   5  * modify it under the terms of the GNU General Public License\r
   6  * as published by the Free Software Foundation; either version 2\r
   7  * of the License, or (at your option) any later version.\r
   8  *\r
   9  * This program is distributed in the hope that it will be useful,\r
  10  * but WITHOUT ANY WARRANTY; without even the implied warranty of\r
  11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the\r
  12  * GNU General Public License for more details.\r
  13  *\r
  14  * You should have received a copy of the GNU General Public License\r
  15  * along with this program; if not, write to the Free Software\r
  16  * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.\r
  17  */\r
  18 \r
  19 package jalview.schemes;\r
  20 \r
  21 import jalview.datamodel.*;\r
  22 import java.util.*;\r
  23 import java.awt.*;\r
  24 \r
  25 public class ClustalxColourScheme implements ColourSchemeI\r
  26 {\r
  27   Hashtable[] cons;\r
  28   int[][] cons2;\r
  29   ConsensusColour[] colours;\r
  30   ConsensusColour[] ResidueColour;\r
  31   int size;\r
  32     Consensus[] conses = new Consensus[32];\r
  33     public static Hashtable colhash = new Hashtable();\r
  34 \r
  35     {\r
  36         colhash.put("RED",new Color((float)0.9,(float)0.2,(float)0.1));\r
  37         colhash.put("BLUE",new Color((float)0.5,(float)0.7,(float)0.9));\r
  38         colhash.put("GREEN",new Color((float)0.1,(float)0.8,(float)0.1));\r
  39         colhash.put("ORANGE",new Color((float)0.9,(float)0.6,(float)0.3));\r
  40         colhash.put("CYAN",new Color((float)0.1,(float)0.7,(float)0.7));\r
  41         colhash.put("PINK",new Color((float)0.9,(float)0.5,(float)0.5));\r
  42         colhash.put("MAGENTA",new Color((float)0.8,(float)0.3,(float)0.8));\r
  43         colhash.put("YELLOW",new Color((float)0.8,(float)0.8,(float)0.0));\r
  44 \r
  45     }\r
  46   Vector colourTable = new Vector();\r
  47 \r
  48 \r
  49   public ClustalxColourScheme(Vector seqs, int maxWidth)\r
  50   {\r
  51     cons2 = new int[maxWidth][24];\r
  52     int start = 0, end = maxWidth-1;\r
  53 \r
  54     // Initialize the array\r
  55     for (int j=0;j<24;j++)\r
  56       for (int i=0; i < maxWidth;i++)\r
  57         cons2[i][j] = 0;\r
  58 \r
  59       int res,i, j = 0;\r
  60       String seq;\r
  61       while(j < seqs.size())\r
  62       {\r
  63         seq = ( (SequenceI) seqs.get(j)).getSequence();\r
  64         int end_j = seq.length()-1;\r
  65         for (i = start; i <= end_j; i++)\r
  66         {\r
  67           if(seq.length()-1<i)\r
  68             res = 23;\r
  69           else\r
  70             res = ( (Integer) ResidueProperties.aaHash.get(seq.charAt(i) + "")).intValue();\r
  71           cons2[i][res]++;\r
  72         }\r
  73         j++;\r
  74       }\r
  75 \r
  76     this.size = seqs.size();\r
  77     makeColours();\r
  78   }\r
  79   public void makeColours() {\r
  80     conses[0] = new Consensus("WLVIMAFCYHP",60);\r
  81     conses[1] = new Consensus("WLVIMAFCYHP",80);\r
  82     conses[2] = new Consensus("ED",50);\r
  83     conses[3] = new Consensus("KR",60);\r
  84     conses[4] = new Consensus("G",50);\r
  85     conses[5] = new Consensus("N",50);\r
  86     conses[6] = new Consensus("QE",50);\r
  87     conses[7] = new Consensus("P",50);\r
  88     conses[8] = new Consensus("TS",50);\r
  89 \r
  90     conses[26] = new Consensus("A",85);\r
  91     conses[27] = new Consensus("C",85);\r
  92     conses[10] = new Consensus("E",85);\r
  93     conses[11] = new Consensus("F",85);\r
  94     conses[12] = new Consensus("G",85);\r
  95     conses[13] = new Consensus("H",85);\r
  96     conses[14] = new Consensus("I",85);\r
  97     conses[15] = new Consensus("L",85);\r
  98     conses[16] = new Consensus("M",85);\r
  99     conses[17] = new Consensus("N",85);\r
 100     conses[18] = new Consensus("P",85);\r
 101     conses[19] = new Consensus("Q",85);\r
 102     conses[20] = new Consensus("R",85);\r
 103     conses[21] = new Consensus("S",85);\r
 104     conses[22] = new Consensus("T",85);\r
 105     conses[23] = new Consensus("V",85);\r
 106     conses[24] = new Consensus("W",85);\r
 107     conses[25] = new Consensus("Y",85);\r
 108     conses[28] = new Consensus("K",85);\r
 109     conses[29] = new Consensus("D",85);\r
 110 \r
 111     conses[30] = new Consensus("G",0);\r
 112     conses[31] = new Consensus("P",0);\r
 113 \r
 114     // We now construct the colours\r
 115     colours = new ConsensusColour[11];\r
 116 \r
 117 \r
 118     Consensus[] tmp8 = new Consensus[1];\r
 119     tmp8[0] = conses[30];//G\r
 120     colours[7] =new ConsensusColour((Color)colhash.get("ORANGE"),tmp8);\r
 121 \r
 122     Consensus[] tmp9 = new Consensus[1];\r
 123     tmp9[0] = conses[31];//P\r
 124     colours[8] =new ConsensusColour((Color)colhash.get("YELLOW"),tmp9);\r
 125 \r
 126     Consensus[] tmp10 = new Consensus[1];\r
 127     tmp10[0] = conses[27];//C\r
 128     colours[9] =new ConsensusColour((Color)colhash.get("PINK"),tmp8);\r
 129 \r
 130     Consensus[] tmp1 = new Consensus[14];\r
 131     tmp1[0] = conses[0]; //%\r
 132     tmp1[1] = conses[1]; //#\r
 133     tmp1[2] = conses[26]; //A\r
 134     tmp1[3] = conses[27]; //C\r
 135     tmp1[4] = conses[11]; //F\r
 136     tmp1[5] = conses[13]; //H\r
 137     tmp1[6] = conses[14]; //I\r
 138     tmp1[7] = conses[15]; //L\r
 139     tmp1[8] = conses[16]; //M\r
 140     tmp1[9] = conses[23]; //V\r
 141     tmp1[10] = conses[24]; //W\r
 142     tmp1[11] = conses[25]; //Y\r
 143     tmp1[12] = conses[18]; //P\r
 144     tmp1[13] = conses[19]; //p\r
 145     colours[0] = new ConsensusColour((Color)colhash.get("BLUE"),tmp1);\r
 146 \r
 147     colours[10] = new ConsensusColour((Color)colhash.get("CYAN"),tmp1);\r
 148 \r
 149     Consensus[] tmp2 = new Consensus[5];\r
 150     tmp2[0] = conses[8];  //t\r
 151     tmp2[1] = conses[21]; //S\r
 152     tmp2[2] = conses[22]; //T\r
 153     tmp2[3] = conses[0];  //%\r
 154     tmp2[4] = conses[1];  //#\r
 155     colours[1] = new ConsensusColour((Color)colhash.get("GREEN"),tmp2);\r
 156 \r
 157     Consensus[] tmp3 = new Consensus[3];\r
 158 \r
 159     tmp3[0] = conses[17]; //N\r
 160     tmp3[1] = conses[29]; //D\r
 161     tmp3[2] = conses[5];  //n\r
 162     colours[2] = new ConsensusColour((Color)colhash.get("GREEN"),tmp3);\r
 163 \r
 164     Consensus[] tmp4 = new Consensus[6];\r
 165     tmp4[0] = conses[6];  // q = QE\r
 166     tmp4[1] = conses[19]; //Q\r
 167     tmp4[2] = conses[22]; //E\r
 168     tmp4[3] = conses[3];  //+\r
 169     tmp4[4] = conses[28]; //K\r
 170     tmp4[5] = conses[20]; //R\r
 171     colours[3] = new ConsensusColour((Color)colhash.get("GREEN"),tmp4);\r
 172     Consensus[] tmp5 = new Consensus[4];\r
 173     tmp5[0] = conses[3];  //+\r
 174     tmp5[1] = conses[28]; //K\r
 175     tmp5[2] = conses[20]; //R\r
 176     tmp5[3] = conses[19]; //Q\r
 177     colours[4] = new ConsensusColour((Color)colhash.get("RED"),tmp5);\r
 178     Consensus[] tmp6 = new Consensus[5];\r
 179     tmp6[0] = conses[3];  //-\r
 180     tmp6[1] = conses[29]; //D\r
 181     tmp6[2] = conses[10]; //E\r
 182     tmp6[3] = conses[6];  //q\r
 183     tmp6[4] = conses[19]; //Q\r
 184     colours[5] = new ConsensusColour((Color)colhash.get("MAGENTA"),tmp6);\r
 185 \r
 186     Consensus[] tmp7 = new Consensus[5];\r
 187     tmp7[0] = conses[3];  //-\r
 188     tmp7[1] = conses[29]; //D\r
 189     tmp7[2] = conses[10]; //E\r
 190     tmp7[3] = conses[17];  //N\r
 191     tmp7[4] = conses[2]; //DE\r
 192     colours[6] = new ConsensusColour((Color)colhash.get("MAGENTA"),tmp7);\r
 193 \r
 194     // Now attach the ConsensusColours to the residue letters\r
 195     ResidueColour = new ConsensusColour[20];\r
 196     ResidueColour[0] = colours[0]; // A\r
 197     ResidueColour[1] = colours[4]; // R\r
 198     ResidueColour[2] = colours[2]; // N\r
 199     ResidueColour[3] = colours[6]; // D\r
 200     ResidueColour[4] = colours[0]; // C\r
 201     ResidueColour[5] = colours[3]; // Q\r
 202     ResidueColour[6] = colours[5]; // E\r
 203     ResidueColour[7] = colours[7]; // G\r
 204     ResidueColour[8] = colours[10]; // H\r
 205     ResidueColour[9] = colours[0]; // I\r
 206     ResidueColour[10] = colours[0]; // L\r
 207     ResidueColour[11] = colours[4]; // K\r
 208     ResidueColour[12] = colours[0]; // M\r
 209     ResidueColour[13] = colours[0]; // F\r
 210     ResidueColour[14] = colours[8]; // P\r
 211     ResidueColour[15] = colours[1]; // S\r
 212     ResidueColour[16] = colours[1]; // T\r
 213     ResidueColour[17] = colours[0]; // W\r
 214     ResidueColour[18] = colours[10]; // Y\r
 215     ResidueColour[19] = colours[0]; // V\r
 216   }\r
 217 \r
 218   public Color findColour(String s)\r
 219   {\r
 220     return Color.pink;\r
 221   }\r
 222 \r
 223  public void setGroupSizeChanged(int i)\r
 224  {\r
 225    size = i;\r
 226  }\r
 227 \r
 228   public Color findColour(String s, int j)\r
 229    {\r
 230      int i = ((Integer)ResidueProperties.aaHash.get(s)).intValue();\r
 231 \r
 232      Color c = Color.white;\r
 233 \r
 234      if(i>19)\r
 235        return c;\r
 236 \r
 237      for (int k=0; k < ResidueColour[i].conses.length ; k++)\r
 238        if (ResidueColour[i].conses[k].isConserved(cons2,j,size))\r
 239          c =  ResidueColour[i].c;\r
 240 \r
 241      if (i == 4)\r
 242     {\r
 243       if (conses[27].isConserved(cons2, j, size))\r
 244         c = (Color) colhash.get("PINK");\r
 245     }\r
 246 \r
 247 \r
 248      return c;\r
 249 \r
 250    }\r
 251 \r
 252   public void setConsensus(Vector consensus)\r
 253   {\r
 254     // Clustal colouring cannot be coloured by PID\r
 255   }\r
 256 \r
 257   public boolean canThreshold()\r
 258   {\r
 259     return false;\r
 260   }\r
 261 \r
 262   public boolean isUserDefinable()\r
 263   {\r
 264     return false;\r
 265   }\r
 266 \r
 267 \r
 268 }\r
 269 \r
 270 class ConsensusColour {\r
 271 \r
 272   Consensus[] conses;\r
 273   Color c;\r
 274 \r
 275   public ConsensusColour(Color c,Consensus[] conses) {\r
 276     this.conses = conses;\r
 277     //    this.list = list;\r
 278     this.c = c;\r
 279   }\r
 280 }\r
 281 \r
 282 \r
 283 \r
 284 \r
 285 \r
 286 \r
 287 \r
 288 \r
 289 \r
 290 \r
 291 \r
 292 \r
 293 \r
 294 \r
 295 \r