aboutsummaryrefslogtreecommitdiff
path: root/src/powder_plot.c
diff options
context:
space:
mode:
authorThomas White <taw@physics.org>2010-02-22 10:25:45 +0100
committerThomas White <taw@physics.org>2010-02-22 10:25:45 +0100
commit79fcc8f0938f0d8056b7b745eb96b89d133a571e (patch)
tree4d65dff5f1332b72db6c12a7ca58ad8d86560836 /src/powder_plot.c
parent3132451eb2a3c88202078177f84c9c8fc96372fe (diff)
Add powder_plot program
Diffstat (limited to 'src/powder_plot.c')
-rw-r--r--src/powder_plot.c111
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;
+}