aboutsummaryrefslogtreecommitdiff
path: root/src
diff options
context:
space:
mode:
authorThomas White <taw@physics.org>2010-01-18 16:16:20 +0100
committerThomas White <taw@physics.org>2010-01-18 16:16:20 +0100
commitfc055095a1321784ded839bdcee42182bc3e5469 (patch)
treede58fde0b77a30aad1837c1ae9d95ad78ab18842 /src
parent094e18eed46ea514ca108f0c3dfc602f0e5eeb0c (diff)
Output more information with --dump-peaks
Diffstat (limited to 'src')
-rw-r--r--src/index.c2
-rw-r--r--src/index.h4
-rw-r--r--src/peaks.c18
3 files changed, 22 insertions, 2 deletions
diff --git a/src/index.c b/src/index.c
index 4191ca0b..62d63d06 100644
--- a/src/index.c
+++ b/src/index.c
@@ -29,7 +29,7 @@
/* x,y in pixels relative to central beam */
-static int map_position(struct image *image, double x, double y,
+int map_position(struct image *image, double x, double y,
double *rx, double *ry, double *rz)
{
/* "Input" space */
diff --git a/src/index.h b/src/index.h
index 3f92b96b..076d9fe2 100644
--- a/src/index.h
+++ b/src/index.h
@@ -21,5 +21,9 @@
extern void index_pattern(struct image *image, int no_index, int dump_peaks,
int use_dirax);
+/* x,y in pixels relative to central beam */
+extern int map_position(struct image *image, double x, double y,
+ double *rx, double *ry, double *rz);
+
#endif /* INDEX_H */
diff --git a/src/peaks.c b/src/peaks.c
index d414c912..9722f125 100644
--- a/src/peaks.c
+++ b/src/peaks.c
@@ -22,6 +22,7 @@
#include "image.h"
#include "utils.h"
+#include "index.h"
#define PEAK_WINDOW_SIZE (10)
@@ -192,6 +193,10 @@ void search_peaks(struct image *image, int dump_peaks)
}
image->features = image_feature_list_new();
+ if ( dump_peaks ) {
+ printf("x/px\ty/px\tq/nm^-1\tPeak I\n");
+ }
+
for ( x=1; x<image->width-1; x++ ) {
for ( y=1; y<image->height-1; y++ ) {
@@ -274,7 +279,18 @@ void search_peaks(struct image *image, int dump_peaks)
/* Map and record reflection */
if ( dump_peaks ) {
- printf("%i %i\n", x, y);
+
+ double q, rx, ry, rz;
+
+ map_position(image, mask_x, mask_y,
+ &rx, &ry, &rz);
+
+
+ q = modulus(rx, ry, rz);
+
+ printf("%i\t%i\t%f\t%i\n", x, y, q/1.0e9,
+ data[mask_x+width*mask_y]);
+
}
image_add_feature(image->features,