Add powder_plot program

1 files changed, 111 insertions, 0 deletions

diff --git a/src/powder_plot.c b/src/powder_plot.c
new file mode 100644
index 00000000..3f0c1abd
--- /dev/null
+++ b/src/powder_plot.c
@@ -0,0 +1,111 @@
+/*
+ * powder_plot.c
+ *
+ * Plot powder patterns
+ *
+ * (c) 2006-2010 Thomas White <taw@physics.org>
+ *
+ * Part of CrystFEL - crystallography with a FEL
+ *
+ */
+
+
+#ifdef HAVE_CONFIG_H
+#include <config.h>
+#endif
+
+#include <stdarg.h>
+#include <stdlib.h>
+#include <stdio.h>
+#include <string.h>
+#include <unistd.h>
+#include <getopt.h>
+
+#include "utils.h"
+#include "image.h"
+#include "detector.h"
+#include "index.h"
+#include "hdf5-file.h"
+
+
+static void show_help(const char *s)
+{
+	printf("Syntax: %s [options] <file.h5>\n\n", s);
+	printf(
+"Compare intensity lists.\n"
+"\n"
+"  -h, --help                 Display this help message.\n"
+"\n");
+}
+
+
+int main(int argc, char *argv[])
+{
+	int c;
+	struct image image;
+	int x, y;
+	struct hdfile *hdfile;
+	char *filename = NULL;
+
+	/* Long options */
+	const struct option longopts[] = {
+		{"help",               0, NULL,               'h'},
+		{"input",              1, NULL,               'i'},
+		{0, 0, NULL, 0}
+	};
+
+	/* Short options */
+	while ((c = getopt_long(argc, argv, "hi:", longopts, NULL)) != -1) {
+
+		switch (c) {
+		case 'h' : {
+			show_help(argv[0]);
+			return 0;
+		}
+
+		case 0 : {
+			break;
+		}
+
+		case 'i' : {
+			filename = strdup(optarg);
+			break;
+		}
+
+		default : {
+			return 1;
+		}
+		}
+
+	}
+
+	if ( filename == NULL ) {
+		ERROR("You must specify the input filename with -i\n");
+		return 1;
+	}
+
+	#include "geometry-lcls.tmp"
+
+	hdfile = hdfile_open(filename);
+	hdfile_set_image(hdfile, "/data/data");
+	hdf5_read(hdfile, &image);
+
+	for ( x=0; x<image.width; x++ ) {
+	for ( y=0; y<image.height; y++ ) {
+
+		double rx, ry, rz;
+		double q;
+		int intensity;
+
+		map_position(&image, x, y, &rx, &ry, &rz);
+		q = modulus(rx, ry, rz);
+
+		intensity = image.data[x + image.width*y];
+
+		printf("%5i\t%5i\t%7.3f\t%7i\n", x, y, q/1.0e9, intensity);
+
+	}
+	}
+
+	return 0;
+}