110f5e12a634db49ea9aa1ea23ff4965f1c2befe galt Tue Aug 14 13:48:57 2018 -0700 changing cse to soe in domains, sometimes gi. ref #21876 diff --git src/utils/wigToBigWig/wigToBigWig.c src/utils/wigToBigWig/wigToBigWig.c index 6f3b6d1..177a480 100644 --- src/utils/wigToBigWig/wigToBigWig.c +++ src/utils/wigToBigWig/wigToBigWig.c @@ -1,84 +1,84 @@ /* wigToBigWig - Convert ascii format wig file (in fixedStep, variableStep or bedGraph format) * to binary big wig format.. */ /* Copyright (C) 2014 The Regents of the University of California * See README in this or parent directory for licensing information. */ #include "common.h" #include "obscure.h" #include "linefile.h" #include "hash.h" #include "options.h" #include "bigWig.h" #include "bwgInternal.h" #include "zlibFace.h" static int blockSize = 256; static int itemsPerSlot = 1024; static boolean clipDontDie = FALSE; static boolean doCompress = FALSE; static boolean fixedSummaries = FALSE; static boolean keepAllChromosomes = FALSE; void usage() /* Explain usage and exit. */ { errAbort( "wigToBigWig v %d - Convert ascii format wig file (in fixedStep, variableStep\n" "or bedGraph format) to binary big wig format.\n" "usage:\n" " wigToBigWig in.wig chrom.sizes out.bw\n" "Where in.wig is in one of the ascii wiggle formats, but not including track lines\n" "and chrom.sizes is a two-column file/URL: <chromosome name> <size in bases>\n" "and out.bw is the output indexed big wig file.\n" "If the assembly <db> is hosted by UCSC, chrom.sizes can be a URL like\n" - " http://hgdownload.cse.ucsc.edu/goldenPath/<db>/bigZips/<db>.chrom.sizes\n" + " http://hgdownload.soe.ucsc.edu/goldenPath/<db>/bigZips/<db>.chrom.sizes\n" "or you may use the script fetchChromSizes to download the chrom.sizes file.\n" "If not hosted by UCSC, a chrom.sizes file can be generated by running\n" "twoBitInfo on the assembly .2bit file.\n" "options:\n" " -blockSize=N - Number of items to bundle in r-tree. Default %d\n" " -itemsPerSlot=N - Number of data points bundled at lowest level. Default %d\n" " -clip - If set just issue warning messages rather than dying if wig\n" " file contains items off end of chromosome.\n" " -unc - If set, do not use compression.\n" " -fixedSummaries - If set, use a predefined sequence of summary levels.\n" " -keepAllChromosomes - If set, store all chromosomes in b-tree." , bbiCurrentVersion, blockSize, itemsPerSlot ); } static struct optionSpec options[] = { {"blockSize", OPTION_INT}, {"itemsPerSlot", OPTION_INT}, {"clip", OPTION_BOOLEAN}, {"unc", OPTION_BOOLEAN}, {"fixedSummaries", OPTION_BOOLEAN}, {"keepAllChromosomes", OPTION_BOOLEAN}, {NULL, 0}, }; void wigToBigWig(char *inName, char *chromSizes, char *outName) /* wigToBigWig - Convert ascii format wig file (in fixedStep, variableStep or bedGraph format) * to binary big wig format.. */ { bigWigFileCreateEx(inName, chromSizes, blockSize, itemsPerSlot, clipDontDie, doCompress, keepAllChromosomes, fixedSummaries, outName); } int main(int argc, char *argv[]) /* Process command line. */ { optionInit(&argc, argv, options); blockSize = optionInt("blockSize", blockSize); itemsPerSlot = optionInt("itemsPerSlot", itemsPerSlot); clipDontDie = optionExists("clip"); doCompress = !optionExists("unc"); keepAllChromosomes = optionExists("keepAllChromosomes"); fixedSummaries = optionExists("fixedSummaries"); if (argc != 4) usage(); wigToBigWig(argv[1], argv[2], argv[3]); if (verboseLevel() > 1) printVmPeak(); return 0; }