src/jalview/io/IdentifyFile.java

   1 /*
   2  * Jalview - A Sequence Alignment Editor and Viewer (Version 2.8.2b1)
   3  * Copyright (C) 2014 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.io;
  22
  23 import java.io.IOException;
  24
  25 /**
  26  * DOCUMENT ME!
  27  *
  28  * @author $author$
  29  * @version $Revision$
  30  */
  31 public class IdentifyFile
  32 {
  33   /**
  34    * Identify a datasource's file content.
  35    *
  36    * @note Do not use this method for stream sources - create a FileParse object
  37    *       instead.
  38    *
  39    * @param file
  40    *          DOCUMENT ME!
  41    * @param protocol
  42    *          DOCUMENT ME!
  43    * @return ID String
  44    */
  45   public String Identify(String file, String protocol)
  46   {
  47     String emessage = "UNIDENTIFIED FILE PARSING ERROR";
  48     FileParse parser = null;
  49     try
  50     {
  51       parser = new FileParse(file, protocol);
  52       if (parser.isValid())
  53       {
  54         return Identify(parser);
  55       }
  56     } catch (Exception e)
  57     {
  58       System.err.println("Error whilst identifying");
  59       e.printStackTrace(System.err);
  60       emessage = e.getMessage();
  61     }
  62     if (parser != null)
  63     {
  64       return parser.errormessage;
  65     }
  66     return emessage;
  67   }
  68
  69   public String Identify(FileParse source)
  70   {
  71     return Identify(source, true); // preserves original behaviour prior to
  72     // version 2.3
  73   }
  74
  75   /**
  76    * Identify contents of source, closing it or resetting source to start
  77    * afterwards.
  78    *
  79    * @param source
  80    * @param closeSource
  81    * @return filetype string
  82    */
  83   public String Identify(FileParse source, boolean closeSource)
  84   {
  85     String reply = "PFAM";
  86     String data;
  87     int length = 0;
  88     boolean lineswereskipped = false;
  89     boolean isBinary = false; // true if length is non-zero and non-printable
  90     // characters are encountered
  91     try
  92     {
  93       if (!closeSource)
  94       {
  95         source.mark();
  96       }
  97       while ((data = source.nextLine()) != null)
  98       {
  99         length += data.length();
 100         if (!lineswereskipped)
 101         {
 102           for (int i = 0; !isBinary && i < data.length(); i++)
 103           {
 104             char c = data.charAt(i);
 105             isBinary = (c < 32 && c != '\t' && c != '\n' && c != '\r'
 106                     && c != 5 && c != 27); // nominal binary character filter
 107             // excluding CR, LF, tab,DEL and ^E
 108             // for certain blast ids
 109           }
 110         }
 111         if (isBinary)
 112         {
 113           // jar files are special - since they contain all sorts of random
 114           // characters.
 115           if (source.inFile != null)
 116           {
 117             String fileStr = source.inFile.getName();
 118             // possibly a Jalview archive.
 119             if (fileStr.lastIndexOf(".jar") > -1
 120                     || fileStr.lastIndexOf(".zip") > -1)
 121             {
 122               reply = "Jalview";
 123             }
 124           }
 125           if (!lineswereskipped && data.startsWith("PK"))
 126           {
 127             reply = "Jalview"; // archive.
 128             break;
 129           }
 130         }
 131         data = data.toUpperCase();
 132
 133         if ((data.indexOf("# STOCKHOLM") > -1))
 134         {
 135           reply = "STH";
 136
 137           break;
 138         }
 139
 140         if ((data.indexOf("<") > -1))
 141         {
 142           reply = "RNAML";
 143
 144           break;
 145         }
 146
 147         if ((data.length() < 1) || (data.indexOf("#") == 0))
 148         {
 149           lineswereskipped = true;
 150           continue;
 151         }
 152
 153         if (data.indexOf("PILEUP") > -1)
 154         {
 155           reply = "PileUp";
 156
 157           break;
 158         }
 159
 160         if ((data.indexOf("//") == 0)
 161                 || ((data.indexOf("!!") > -1) && (data.indexOf("!!") < data
 162                         .indexOf("_MULTIPLE_ALIGNMENT "))))
 163         {
 164           reply = "MSF";
 165
 166           break;
 167         }
 168         else if (data.indexOf("CLUSTAL") > -1)
 169         {
 170           reply = "CLUSTAL";
 171
 172           break;
 173         }
 174
 175         else if (data.indexOf(">") > -1)
 176         {
 177           // FASTA, PIR file or BLC file
 178           boolean checkPIR = false, starterm = false;
 179           if ((data.indexOf(">P1;") > -1) || (data.indexOf(">DL;") > -1))
 180           {
 181             // watch for PIR file attributes
 182             checkPIR = true;
 183             reply = "PIR";
 184           }
 185           // could also be BLC file, read next line to confirm
 186           data = source.nextLine();
 187
 188           if (data.indexOf(">") > -1)
 189           {
 190             reply = "BLC";
 191           }
 192           else
 193           {
 194             // Is this a single line BLC file?
 195             String data1 = source.nextLine();
 196             String data2 = source.nextLine();
 197             int c1;
 198             if (checkPIR)
 199             {
 200               starterm = (data1 != null && data1.indexOf("*") > -1)
 201                       || (data2 != null && data2.indexOf("*") > -1);
 202             }
 203             if (data2 != null && (c1 = data.indexOf("*")) > -1)
 204             {
 205               if (c1 == 0 && c1 == data2.indexOf("*"))
 206               {
 207                 reply = "BLC";
 208               }
 209               else
 210               {
 211                 reply = "FASTA"; // possibly a bad choice - may be recognised as
 212                 // PIR
 213               }
 214               // otherwise can still possibly be a PIR file
 215             }
 216             else
 217             {
 218               reply = "FASTA";
 219               // TODO : AMSA File is indicated if there is annotation in the
 220               // FASTA file - but FASTA will automatically generate this at the
 221               // mo.
 222               if (!checkPIR)
 223               {
 224                 break;
 225               }
 226             }
 227           }
 228           // final check for PIR content. require
 229           // >P1;title\n<blah>\nterminated sequence to occur at least once.
 230
 231           // TODO the PIR/fasta ambiguity may be the use case that is needed to
 232           // have
 233           // a 'Parse as type XXX' parameter for the applet/application.
 234           if (checkPIR)
 235           {
 236             String dta = null;
 237             if (!starterm)
 238             {
 239               do
 240               {
 241                 try
 242                 {
 243                   dta = source.nextLine();
 244                 } catch (IOException ex)
 245                 {
 246                 }
 247                 ;
 248                 if (dta != null && dta.indexOf("*") > -1)
 249                 {
 250                   starterm = true;
 251                 }
 252               } while (dta != null && !starterm);
 253             }
 254             if (starterm)
 255             {
 256               reply = "PIR";
 257               break;
 258             }
 259             else
 260             {
 261               reply = "FASTA"; // probably a bad choice!
 262             }
 263           }
 264           // read as a FASTA (probably)
 265           break;
 266         }
 267         else if (data.indexOf("HEADER") == 0 || data.indexOf("ATOM") == 0)
 268         {
 269           reply = "PDB";
 270           break;
 271         }
 272         else if (data.matches("\\s*\\d+\\s+\\d+\\s*"))
 273         {
 274           reply = PhylipFile.FILE_DESC;
 275           break;
 276         }
 277
 278         /*
 279          * // TODO comment out SimpleBLAST identification for Jalview 2.4.1 else
 280          * if (!lineswereskipped && data.indexOf("BLAST")<4) { reply =
 281          * "SimpleBLAST"; break;
 282          *
 283          * } // end comments for Jalview 2.4.1
 284          */
 285         else if (!lineswereskipped && data.charAt(0) != '*'
 286                 && data.charAt(0) != ' '
 287                 && data.indexOf(":") < data.indexOf(",")) // &&
 288           // data.indexOf(",")<data.indexOf(",",
 289           // data.indexOf(",")))
 290         {
 291           // file looks like a concise JNet file
 292           reply = "JnetFile";
 293           break;
 294         }
 295
 296         lineswereskipped = true; // this means there was some junk before any
 297         // key file signature
 298       }
 299       if (closeSource)
 300       {
 301         source.close();
 302       }
 303       else
 304       {
 305         source.reset(); // so the file can be parsed from the beginning again.
 306       }
 307     } catch (Exception ex)
 308     {
 309       System.err.println("File Identification failed!\n" + ex);
 310       return source.errormessage;
 311     }
 312     if (length == 0)
 313     {
 314       System.err
 315       .println("File Identification failed! - Empty file was read.");
 316       return "EMPTY DATA FILE";
 317     }
 318     return reply;
 319   }
 320
 321   public static void main(String[] args)
 322   {
 323     for (int i = 0; args != null && i < args.length; i++)
 324     {
 325       IdentifyFile ider = new IdentifyFile();
 326       String type = ider.Identify(args[i], AppletFormatAdapter.FILE);
 327       System.out.println("Type of " + args[i] + " is " + type);
 328     }
 329     if (args == null || args.length == 0)
 330     {
 331       System.err.println("Usage: <Filename> [<Filename> ...]");
 332     }
 333   }
 334 }