src/utils/bigWigToBedGraph/bigWigToBedGraph.c 1.5
1.5 2009/09/08 19:50:25 kent
Adding udcDir to the command line as an option.
Index: src/utils/bigWigToBedGraph/bigWigToBedGraph.c
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/utils/bigWigToBedGraph/bigWigToBedGraph.c,v
retrieving revision 1.4
retrieving revision 1.5
diff -b -B -U 1000000 -r1.4 -r1.5
--- src/utils/bigWigToBedGraph/bigWigToBedGraph.c 12 Aug 2009 21:44:24 -0000 1.4
+++ src/utils/bigWigToBedGraph/bigWigToBedGraph.c 8 Sep 2009 19:50:25 -0000 1.5
@@ -1,75 +1,79 @@
/* bigWigToBedGraph - Convert from bigWig to bedGraph format.. */
#include "common.h"
#include "linefile.h"
#include "hash.h"
#include "options.h"
#include "localmem.h"
+#include "udc.h"
#include "bigWig.h"
static char const rcsid[] = "$Id$";
char *clChrom = NULL;
int clStart = -1;
int clEnd = -1;
void usage()
/* Explain usage and exit. */
{
errAbort(
"bigWigToBedGraph - Convert from bigWig to bedGraph format.\n"
"usage:\n"
" bigWigToBedGraph in.bigWig out.bedGraph\n"
"options:\n"
" -chrom=chr1 - if set restrict output to given chromosome\n"
" -start=N - if set, restrict output to only that over start\n"
" -end=N - if set, restict output to only that under end\n"
+ " -udcDir=/dir/to/cache - place to put cache for remote bigBed/bigWigs\n"
);
}
static struct optionSpec options[] = {
{"chrom", OPTION_STRING},
{"start", OPTION_INT},
{"end", OPTION_INT},
+ {"udcDir", OPTION_STRING},
{NULL, 0},
};
void bigWigToBedGraph(char *inFile, char *outFile)
/* bigWigToBedGraph - Convert from bigWig to bedGraph format.. */
{
struct bbiFile *bwf = bigWigFileOpen(inFile);
FILE *f = mustOpen(outFile, "w");
struct bbiChromInfo *chrom, *chromList = bbiChromList(bwf);
for (chrom = chromList; chrom != NULL; chrom = chrom->next)
{
if (clChrom != NULL && !sameString(clChrom, chrom->name))
continue;
char *chromName = chrom->name;
struct lm *lm = lmInit(0);
int start = 0, end = chrom->size;
if (clStart > 0)
start = clStart;
if (clEnd > 0)
end = clEnd;
struct bbiInterval *interval, *intervalList = bigWigIntervalQuery(bwf, chromName,
start, end, lm);
for (interval = intervalList; interval != NULL; interval = interval->next)
fprintf(f, "%s\t%u\t%u\t%g\n", chromName, interval->start, interval->end, interval->val);
lmCleanup(&lm);
}
bbiChromInfoFreeList(&chromList);
carefulClose(&f);
bbiFileClose(&bwf);
}
int main(int argc, char *argv[])
/* Process command line. */
{
optionInit(&argc, argv, options);
clChrom = optionVal("chrom", clChrom);
clStart = optionInt("start", clStart);
clEnd = optionInt("end", clEnd);
+udcSetDefaultDir(optionVal("udcDir", udcDefaultDir()));
if (argc != 3)
usage();
bigWigToBedGraph(argv[1], argv[2]);
return 0;
}