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         for (i = start; i <= end; i++)\r
  65         {\r
  66           res = ( (Integer) ResidueProperties.aaHash.get(seq.charAt(i) + "")).intValue();\r
  67           cons2[i][res]++;\r
  68         }\r
  69         j++;\r
  70       }\r
  71 \r
  72     this.size = seqs.size();\r
  73     makeColours();\r
  74   }\r
  75   public void makeColours() {\r
  76     conses[0] = new Consensus("WLVIMAFCYHP",60);\r
  77     conses[1] = new Consensus("WLVIMAFCYHP",80);\r
  78     conses[2] = new Consensus("ED",50);\r
  79     conses[3] = new Consensus("KR",60);\r
  80     conses[4] = new Consensus("G",50);\r
  81     conses[5] = new Consensus("N",50);\r
  82     conses[6] = new Consensus("QE",50);\r
  83     conses[7] = new Consensus("P",50);\r
  84     conses[8] = new Consensus("TS",50);\r
  85 \r
  86     conses[26] = new Consensus("A",85);\r
  87     conses[27] = new Consensus("C",85);\r
  88     conses[10] = new Consensus("E",85);\r
  89     conses[11] = new Consensus("F",85);\r
  90     conses[12] = new Consensus("G",85);\r
  91     conses[13] = new Consensus("H",85);\r
  92     conses[14] = new Consensus("I",85);\r
  93     conses[15] = new Consensus("L",85);\r
  94     conses[16] = new Consensus("M",85);\r
  95     conses[17] = new Consensus("N",85);\r
  96     conses[18] = new Consensus("P",85);\r
  97     conses[19] = new Consensus("Q",85);\r
  98     conses[20] = new Consensus("R",85);\r
  99     conses[21] = new Consensus("S",85);\r
 100     conses[22] = new Consensus("T",85);\r
 101     conses[23] = new Consensus("V",85);\r
 102     conses[24] = new Consensus("W",85);\r
 103     conses[25] = new Consensus("Y",85);\r
 104     conses[28] = new Consensus("K",85);\r
 105     conses[29] = new Consensus("D",85);\r
 106 \r
 107     conses[30] = new Consensus("G",0);\r
 108     conses[31] = new Consensus("P",0);\r
 109 \r
 110     // We now construct the colours\r
 111     colours = new ConsensusColour[11];\r
 112 \r
 113 \r
 114     Consensus[] tmp8 = new Consensus[1];\r
 115     tmp8[0] = conses[30];//G\r
 116     colours[7] =new ConsensusColour((Color)colhash.get("ORANGE"),tmp8);\r
 117 \r
 118     Consensus[] tmp9 = new Consensus[1];\r
 119     tmp9[0] = conses[31];//P\r
 120     colours[8] =new ConsensusColour((Color)colhash.get("YELLOW"),tmp9);\r
 121 \r
 122     Consensus[] tmp10 = new Consensus[1];\r
 123     tmp10[0] = conses[27];//C\r
 124     colours[9] =new ConsensusColour((Color)colhash.get("PINK"),tmp8);\r
 125 \r
 126     Consensus[] tmp1 = new Consensus[14];\r
 127     tmp1[0] = conses[0]; //%\r
 128     tmp1[1] = conses[1]; //#\r
 129     tmp1[2] = conses[26]; //A\r
 130     tmp1[3] = conses[27]; //C\r
 131     tmp1[4] = conses[11]; //F\r
 132     tmp1[5] = conses[13]; //H\r
 133     tmp1[6] = conses[14]; //I\r
 134     tmp1[7] = conses[15]; //L\r
 135     tmp1[8] = conses[16]; //M\r
 136     tmp1[9] = conses[23]; //V\r
 137     tmp1[10] = conses[24]; //W\r
 138     tmp1[11] = conses[25]; //Y\r
 139     tmp1[12] = conses[18]; //P\r
 140     tmp1[13] = conses[19]; //p\r
 141     colours[0] = new ConsensusColour((Color)colhash.get("BLUE"),tmp1);\r
 142 \r
 143     colours[10] = new ConsensusColour((Color)colhash.get("CYAN"),tmp1);\r
 144 \r
 145     Consensus[] tmp2 = new Consensus[5];\r
 146     tmp2[0] = conses[8];  //t\r
 147     tmp2[1] = conses[21]; //S\r
 148     tmp2[2] = conses[22]; //T\r
 149     tmp2[3] = conses[0];  //%\r
 150     tmp2[4] = conses[1];  //#\r
 151     colours[1] = new ConsensusColour((Color)colhash.get("GREEN"),tmp2);\r
 152 \r
 153     Consensus[] tmp3 = new Consensus[3];\r
 154 \r
 155     tmp3[0] = conses[17]; //N\r
 156     tmp3[1] = conses[29]; //D\r
 157     tmp3[2] = conses[5];  //n\r
 158     colours[2] = new ConsensusColour((Color)colhash.get("GREEN"),tmp3);\r
 159 \r
 160     Consensus[] tmp4 = new Consensus[6];\r
 161     tmp4[0] = conses[6];  // q = QE\r
 162     tmp4[1] = conses[19]; //Q\r
 163     tmp4[2] = conses[22]; //E\r
 164     tmp4[3] = conses[3];  //+\r
 165     tmp4[4] = conses[28]; //K\r
 166     tmp4[5] = conses[20]; //R\r
 167     colours[3] = new ConsensusColour((Color)colhash.get("GREEN"),tmp4);\r
 168     Consensus[] tmp5 = new Consensus[4];\r
 169     tmp5[0] = conses[3];  //+\r
 170     tmp5[1] = conses[28]; //K\r
 171     tmp5[2] = conses[20]; //R\r
 172     tmp5[3] = conses[19]; //Q\r
 173     colours[4] = new ConsensusColour((Color)colhash.get("RED"),tmp5);\r
 174     Consensus[] tmp6 = new Consensus[5];\r
 175     tmp6[0] = conses[3];  //-\r
 176     tmp6[1] = conses[29]; //D\r
 177     tmp6[2] = conses[10]; //E\r
 178     tmp6[3] = conses[6];  //q\r
 179     tmp6[4] = conses[19]; //Q\r
 180     colours[5] = new ConsensusColour((Color)colhash.get("MAGENTA"),tmp6);\r
 181 \r
 182     Consensus[] tmp7 = new Consensus[5];\r
 183     tmp7[0] = conses[3];  //-\r
 184     tmp7[1] = conses[29]; //D\r
 185     tmp7[2] = conses[10]; //E\r
 186     tmp7[3] = conses[17];  //N\r
 187     tmp7[4] = conses[2]; //DE\r
 188     colours[6] = new ConsensusColour((Color)colhash.get("MAGENTA"),tmp7);\r
 189 \r
 190     // Now attach the ConsensusColours to the residue letters\r
 191     ResidueColour = new ConsensusColour[20];\r
 192     ResidueColour[0] = colours[0]; // A\r
 193     ResidueColour[1] = colours[4]; // R\r
 194     ResidueColour[2] = colours[2]; // N\r
 195     ResidueColour[3] = colours[6]; // D\r
 196     ResidueColour[4] = colours[0]; // C\r
 197     ResidueColour[5] = colours[3]; // Q\r
 198     ResidueColour[6] = colours[5]; // E\r
 199     ResidueColour[7] = colours[7]; // G\r
 200     ResidueColour[8] = colours[10]; // H\r
 201     ResidueColour[9] = colours[0]; // I\r
 202     ResidueColour[10] = colours[0]; // L\r
 203     ResidueColour[11] = colours[4]; // K\r
 204     ResidueColour[12] = colours[0]; // M\r
 205     ResidueColour[13] = colours[0]; // F\r
 206     ResidueColour[14] = colours[8]; // P\r
 207     ResidueColour[15] = colours[1]; // S\r
 208     ResidueColour[16] = colours[1]; // T\r
 209     ResidueColour[17] = colours[0]; // W\r
 210     ResidueColour[18] = colours[10]; // Y\r
 211     ResidueColour[19] = colours[0]; // V\r
 212   }\r
 213 \r
 214   public Color findColour(String s)\r
 215   {\r
 216     return Color.pink;\r
 217   }\r
 218 \r
 219  public void setGroupSizeChanged(int i)\r
 220  {\r
 221    size = i;\r
 222  }\r
 223 \r
 224   public Color findColour(String s, int j, Vector seqs)\r
 225    {\r
 226      int i = ((Integer)ResidueProperties.aaHash.get(s)).intValue();\r
 227 \r
 228      Color c = Color.white;\r
 229 \r
 230      if(i>20)\r
 231        return c;\r
 232 \r
 233      for (int k=0; k < ResidueColour[i].conses.length ; k++)\r
 234        if (ResidueColour[i].conses[k].isConserved(cons2,j,size))\r
 235          c =  ResidueColour[i].c;\r
 236 \r
 237      if (i == 4)\r
 238     {\r
 239       if (conses[27].isConserved(cons2, j, size))\r
 240         c = (Color) colhash.get("PINK");\r
 241     }\r
 242 \r
 243 \r
 244      return c;\r
 245 \r
 246    }\r
 247 \r
 248 \r
 249   public boolean canThreshold()\r
 250   {\r
 251     return false;\r
 252   }\r
 253 \r
 254   public boolean isUserDefinable()\r
 255   {\r
 256     return false;\r
 257   }\r
 258 \r
 259 \r
 260 }\r
 261 \r
 262 class ConsensusColour {\r
 263 \r
 264   Consensus[] conses;\r
 265   Color c;\r
 266 \r
 267   public ConsensusColour(Color c,Consensus[] conses) {\r
 268     this.conses = conses;\r
 269     //    this.list = list;\r
 270     this.c = c;\r
 271   }\r
 272 }\r
 273 \r
 274 \r
 275 \r
 276 \r
 277 \r
 278 \r
 279 \r
 280 \r
 281 \r
 282 \r
 283 \r
 284 \r
 285 \r
 286 \r
 287 \r