diff options
| author | Thomas White <taw@physics.org> | 2014-10-24 17:50:22 +0200 |
|---|---|---|
| committer | Thomas White <taw@physics.org> | 2014-10-24 17:50:22 +0200 |
| commit | 74767e4c66003622bdb3d9f1c84d7a73df470137 (patch) | |
| tree | 35051a569144db49966a61889ca1416353a28ce3 /src/dw-hdfsee.c | |
| parent | 90bf0eee85d7f5687c80464cef9333071f1cd72e (diff) | |
hdfsee: Read all possible reflection formats
Diffstat (limited to 'src/dw-hdfsee.c')
| -rw-r--r-- | src/dw-hdfsee.c | 34 |
1 files changed, 31 insertions, 3 deletions
diff --git a/src/dw-hdfsee.c b/src/dw-hdfsee.c index d22f83c6..f4655d0a 100644 --- a/src/dw-hdfsee.c +++ b/src/dw-hdfsee.c @@ -1028,20 +1028,21 @@ static void load_features_from_file(struct image *image, const char *filename) if ( rval == NULL ) continue; chomp(line); - /* Try long reflection format (from stream) */ r = sscanf(line, "%i %i %i %f %s %f %f %f %f %s", &h, &k, &l, &intensity, phs, &sigma, &cts, &fs, &ss, pn); if ( r == 10 ) { + /* Stream reflection list format 2.3 */ char name[32]; snprintf(name, 31, "%i %i %i", h, k, l); p = find_panel_by_name(image->det, pn); if ( p == NULL ) { - ERROR("Unable to find panel %s\n", pn); + ERROR("Unable to find panel %s " + "(no geometry file given?)\n", pn); } else { @@ -1056,8 +1057,34 @@ static void load_features_from_file(struct image *image, const char *filename) image, 1.0, strdup(name)); continue; + } else if ( r == 9 ) { + + /* Stream reflection list format 2.2 or 2.1 */ + char name[32]; + + snprintf(name, 31, "%i %i %i", h, k, l); + + p = find_orig_panel(image->det, fs, ss); + + if ( p == NULL ) { + ERROR("Unable to find panel for %s " + "(no geometry file given?)\n", name); + } else { + + /* Convert coordinates to match rearranged + * panels in memory */ + fs = fs - p->orig_min_fs + p->min_fs; + ss = ss - p->orig_min_ss + p->min_ss; + + } + + image_add_feature(image->features, fs, ss, + image, 1.0, strdup(name)); + continue; + } + /* Try long peak format from stream */ r = sscanf(line, "%f %f %f %f %s", &fs, &ss, &d, &intensity, pn); @@ -1065,7 +1092,8 @@ static void load_features_from_file(struct image *image, const char *filename) p = find_panel_by_name(image->det, pn); if ( p == NULL ) { - ERROR("Unable to find panel %s\n", pn); + ERROR("Unable to find panel %s " + "(no geometry file given?)\n", pn); } else { /* Convert coordinates to match rearranged panels in |