src/hg/logCrawl/ticksToWig/ticksToWig.c 1.1

Index: src/hg/logCrawl/ticksToWig/ticksToWig.c
===================================================================
RCS file: src/hg/logCrawl/ticksToWig/ticksToWig.c
diff -N src/hg/logCrawl/ticksToWig/ticksToWig.c
--- /dev/null	1 Jan 1970 00:00:00 -0000
+++ src/hg/logCrawl/ticksToWig/ticksToWig.c	9 Sep 2009 19:11:27 -0000	1.1
@@ -0,0 +1,129 @@
+/* ticksToWig - Convert tab-delimited file of Unix time ticks, and possibly also numerical values to wig file(s).. */
+#include "common.h"
+#include "linefile.h"
+#include "hash.h"
+#include "options.h"
+#include "apacheLog.h"
+
+static char const rcsid[] = "$Id$";
+
+int tickCol = 0;
+int valCol = 1;
+
+void usage()
+/* Explain usage and exit. */
+{
+errAbort(
+  "ticksToWig - Convert tab-delimited file of Unix time ticks, and possibly also numerical values to wig file(s).\n"
+  "usage:\n"
+  "   ticksToWig in.tab density.wig [averageVal.wig]\n"
+  "options:\n"
+  "   -tickCol=N - Column for ticks, defaults to 0\n"
+  "   -valCol=N -  Column for value, defaults to 1\n"
+  "   -startTime=day/Mon/Year:HR:MN:SC\n"
+  "   -startTick=tick - start time as a unix tick.\n"
+  "The startTime and startTick values set where the zero time should be.  Otherwise it defaults\n"
+  "to the tick in the first line of in.tab.\n"
+  );
+}
+
+static struct optionSpec options[] = {
+   {"tickCol", OPTION_INT},
+   {"valCol", OPTION_INT},
+   {"startTime", OPTION_STRING},
+   {"startTick", OPTION_INT},
+   {NULL, 0},
+};
+
+void printVarStepHead(FILE *f)
+/* Print variable step header. */
+{
+fprintf(f, "variableStep chrom=chr1 span=1\n");
+}
+
+void ticksToWig(int startTick, char *inTable, char *outDensity, char *outAverage)
+/* ticksToWig - Convert tab-delimited file of Unix time ticks, and possibly also 
+ * numerical values to wig file(s).. */
+{
+struct lineFile *lf = lineFileOpen(inTable, TRUE);
+FILE *densityFile = mustOpen(outDensity, "w");
+printVarStepHead(densityFile);
+FILE *averageFile = NULL;
+if (outAverage != NULL)
+    {
+    averageFile = mustOpen(outAverage, "w");
+    printVarStepHead(averageFile);
+    }
+int colsToParse = 1 + max(tickCol, valCol);
+char *row[colsToParse];
+
+time_t curTick = 0;
+int sameTickCount = 0;
+double tickTotal = 0;
+double val = 0;
+time_t tick;
+while (lineFileNextRow(lf, row, colsToParse))
+    {
+    tick = lineFileNeedNum(lf, row, tickCol);
+    if (averageFile != NULL)
+       val = lineFileNeedDouble(lf, row, valCol);
+    if (curTick != tick)
+        {
+	if (curTick > tick)
+	    errAbort("Input isn't sorted - %ld > %ld line %d of %s\n", 
+	    	(long)curTick, (long)tick, lf->lineIx, lf->fileName);
+	if (startTick == 0)
+	    startTick = tick;
+        if (sameTickCount > 0)
+	    {
+	    fprintf(densityFile, "%ld\t%d\n", curTick - startTick + 1, sameTickCount);
+	    time_t i;
+	    for (i=curTick+1; i<tick; ++i)
+		{
+		fprintf(densityFile, "%ld\t%d\n", i - startTick + 1, 0);
+		}
+	   if (averageFile != NULL)
+	       {
+	       fprintf(averageFile, "%ld\t%f\n", 
+	       		(long)curTick - startTick + 1, tickTotal/sameTickCount);
+	       tickTotal = 0;
+	       }
+	    sameTickCount = 0;
+	    }
+        curTick = tick;
+	}
+    tickTotal += val;
+    sameTickCount += 1;
+    }
+if (sameTickCount > 0)
+   {
+   fprintf(densityFile, "%ld\t%d\n", curTick - startTick + 1, sameTickCount);
+   if (averageFile != NULL)
+       fprintf(averageFile, "%ld\t%f\n", 
+       		(long)curTick - startTick + 1, tickTotal/sameTickCount);
+   }
+carefulClose(&densityFile);
+carefulClose(&averageFile);
+}
+
+int main(int argc, char *argv[])
+/* Process command line. */
+{
+optionInit(&argc, argv, options);
+if (argc != 3 && argc != 4)
+    usage();
+time_t startTick = 0;
+char *startTime = optionVal("startTime", NULL);
+if (startTime != NULL)
+    {
+    startTick = apacheAccessLogTimeToTick(startTime);
+    }
+else if (optionExists("startTick"))
+    {
+    startTick = optionInt("startTick", 0);
+    }
+tickCol = optionInt("tickCol", tickCol);
+valCol = optionInt("valCol", valCol);
+ticksToWig(startTick, argv[1], argv[2], (argc == 4 ? argv[3] : NULL));
+return 0;
+}