/*
* The MIT License (MIT)
*
* Copyright (c) 2007-2015 Broad Institute
*
* Permission is hereby granted, free of charge, to any person obtaining a copy
* of this software and associated documentation files (the "Software"), to deal
* in the Software without restriction, including without limitation the rights
* to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
* copies of the Software, and to permit persons to whom the Software is
* furnished to do so, subject to the following conditions:
*
* The above copyright notice and this permission notice shall be included in
* all copies or substantial portions of the Software.
*
*
* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
* IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
* FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
* AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
* LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
* OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
* THE SOFTWARE.
*/
/*
* To change this template, choose Tools | Templates
* and open the template in the editor.
*/
package org.broad.igv.tools.sort;
import org.apache.log4j.Logger;
import htsjdk.tribble.readers.AsciiLineReader;
import java.io.*;
/**
* @author jrobinso
*/
public class CNSorter extends AsciiSorter {
static private Logger log = Logger.getLogger(AsciiSorter.class);
public CNSorter(File inputFile, File outputFile) {
super(inputFile, outputFile);
}
String writeHeader(AsciiLineReader reader, PrintWriter writer) throws IOException {
String nextLine = reader.readLine();
while (nextLine.startsWith("#") || (nextLine.trim().length() == 0)) {
writer.println(nextLine);
nextLine = reader.readLine();
}
// column headers
writer.println(nextLine);
if(!nextLine.startsWith("SNP")){
log.warn("Expected header line not found");
}
return null;
}
Parser getParser() throws IOException {
String tmp = inputFile.getName();
String fn = tmp.endsWith(".txt") ? tmp.substring(0, tmp.length() - 4) : tmp;
if (fn.endsWith(".cn") || fn.endsWith(".xcn") || fn.endsWith(".snp")) {
return new Parser(1, 2);
} else if (fn.endsWith(".igv")) {
return getIGVParser();
} else {
throw new RuntimeException("Unrecognized copy number extension: " + fn);
}
}
// THe IGV formats allow user overrides of column numbers.
Parser getIGVParser() throws IOException {
int chrColumn = 0;
int startColumn = 1;
BufferedReader reader = null;
try {
reader = new BufferedReader(new FileReader(inputFile));
String nextLine = "";
while ((nextLine = reader.readLine()) != null) {
if (!nextLine.startsWith("#")) {
break;
} else if (nextLine.startsWith("#columns")) {
String[] tokens = nextLine.split("\\s+");
if (tokens.length > 1) {
for (int i = 1; i < tokens.length; i++) {
String[] kv = tokens[i].split("=");
if (kv.length == 2) {
if (kv[0].toLowerCase().equals("chr")) {
int c = Integer.parseInt(kv[1]);
if (c < 1) {
throw new RuntimeException("Error parsing column line: " + nextLine + ". Column numbers must be > 0");
} else {
chrColumn = c - 1;
}
} else if (kv[0].toLowerCase().equals("start")) {
int c = Integer.parseInt(kv[1]);
if (c < 1) {
throw new RuntimeException("Error parsing column line: " + nextLine + ". Column numbers must be > 0");
} else {
startColumn = c - 1;
}
}
}
}
}
break;
}
}
return new Parser(chrColumn, startColumn);
} finally {
if (reader != null) {
reader.close();
}
}
}
}